Search USGSSearch

SEARCH · Search USGS

Results for “Scientific Investigations Report”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,279 records · Page 71Linked to original sources

Anthropogenic organic compounds in source water of selected community water systems that use groundwater, 2002-05

Source water, defined as groundwater collected from a community water system well prior to water treatment, was sampled from 221 wells during October 2002 to July 2005 and analyzed for 258 anthropogenic organic compounds. Most of these compounds are unregulated in drinking water and include pesticides and pesticide degradates, gasoline hydrocarbons, personal-care and domestic-use products, and solvents. The laboratory analytical methods used in the study have detection levels that commonly are 100 to 1,000 times lower than State and Federal standards and guidelines for protecting water quality. Detections of anthropogenic organic compounds do not necessarily indicate a concern to human health but rather help to identify emerging issues and track changes in occurrence and concentrations over time. Less than one-half (120) of the 258 compounds were detected in at least one source-water sample. Chloroform, in 36 percent of samples, was the most commonly detected of the 12 compounds that were in about 10 percent or more of source-water samples. The herbicides atrazine, metolachlor, prometon, and simazine also were among the commonly detected compounds. The commonly detected degradates of atrazine - deethylatrazine and deisopropylatrazine - as well as degradates of acetochlor and alachlor, generally were detected at concentrations similar to or greater than concentrations of the parent herbicide. The compounds perchloroethene, trichloroethene, 1,1,1-trichloroethane, methyl tert-butyl ether, and cis-1,2-dichloroethene also were detected commonly. The most commonly detected compounds in source-water samples generally were among those detected commonly across the country and reported in previous studies by the U.S. Geological Survey's National Water-Quality Assessment Program. Relatively few compounds were detected at concentrations greater than human-health benchmarks, and 84 percent of the concentrations were two or more orders of magnitude less than benchmarks. Five compounds (perchloroethene, trichloroethene, 1,2-dibromoethane, acrylonitrile, and dieldrin) were detected at concentrations greater than their human-health benchmark. The human-health benchmarks used for comparison were U.S. Environmental Protection Agency Maximum Contaminant Levels (MCLs) for regulated compounds and Health-Based Screening Levels developed by the U.S. Geological Survey in collaboration with the U.S. Environmental Protection Agency and other agencies for unregulated compounds. About one-half of all detected compounds do not have human-health benchmarks or adequate toxicity information to evaluate results in a human-health context. Ninety-four source-water and finished-water (water that has passed through all the treatment processes but prior to distribution) sites were sampled at selected community water systems during June 2004 to September 2005. Most of the samples were analyzed for compounds that were detected commonly or at relatively high concentrations during the initial source-water sampling. The majority of the finished-water samples represented water blended with water from one or more other wells. Thirty-four samples were from water systems that did not blend water from sampled wells with water from other wells prior to distribution. The comparison of source- and finished-water samples represents an initial assessment of whether compounds present in source water also are present in finished water and is not intended as an evaluation of water-treatment efficacy. The treatment used at the majority of the community water systems sampled is disinfection, which, in general, is not designed to remove the compounds monitored in this study. Concentrations of all compounds detected in finished water were less than their human-health benchmarks. Two detections of perchloroethene and one detection of trichloroethene in finished water had concentrations within an order of magnitude of the MCL. Concentrations of disinfection by-products were

Scientific Investigations Report

Sources of mercury in sediments, water, and fish of the lakes of Whatcom County, Washington

Concerns about mercury (Hg) contamination in Lake Whatcom, Washington, were raised in the late 1990s after a watershed protection survey reported elevated concentrations of Hg in smallmouth bass. The U.S. Geological Survey, the Whatcom County Health Department, and the Washington State Department of Ecology (Ecology) cooperated to develop a study to review existing data and collect new data that would lead to a better understanding of Hg deposition to Lake Whatcom and other lakes in Whatcom County, Washington. A simple atmospheric deposition model was developed that allowed comparisons of the deposition of Hg to the surfaces of each lake. Estimates of Hg deposition derived from the model indicated that the most significant deposition of Hg would have occurred to the lakes north of the City of Bellingham. These lakes were in the primary wind pattern of two municipal waste incinerators. Of all the lakes examined, basin 1 of Lake Whatcom would have been most affected by the Hg emissions from the chlor-alkali plant and the municipal sewage-sludge incinerator in the City of Bellingham. The length-adjusted concentrations of Hg in largemouth and smallmouth bass were not related to estimated deposition rates of Hg to the lakes from local atmospheric sources. Total Hg concentrations in the surface sediments of Lake Whatcom are affected by the sedimentation of fine-grained particles, whereas organic carbon regulates the concentration of methyl-Hg in the surface sediments of the lake. Hg concentrations in dated sediment core samples indicate that increases in Hg sedimentation were largest during the first half of the 20th century. Increases in Hg sedimentation were smaller after the chlor-alkali plant and the incinerators began operating between 1964 and 1984. Analysis of sediments recently deposited in basin 1 of Lake Whatcom, Lake Terrell, and Lake Samish indicates a decrease in Hg sedimentation. Concentrations of Hg in Seattle precipitation and in tributary waters were used to calculate current (2002-03) loadings of Hg to Lake Whatcom. Hg in tributaries contributed 59 percent of the total Hg, whereas non-local atmospheric deposition was estimated to have contributed 41 percent of the 303 grams of Hg entering Lake Whatcom each year. However, these inputs cannot be verified without a better understanding of the sources of sediment to Lake Whatcom.

Scientific Investigations Report

Surface-Water and Ground-Water Resources of Kendall County, Illinois

Water-supply needs in Kendall County, in northern Illinois, are met exclusively from ground water derived from glacial drift aquifers and bedrock aquifers open to Silurian, Ordovician, and Cambrian System units. As a result of population growth in Kendall County and the surrounding area, water use has increased from about 1.2 million gallons per day in 1957 to more than 5 million gallons per day in 2000. The purpose of this report is to characterize the surface-water and ground-water resources of Kendall County. The report presents a compilation of available information on geology, surface-water and ground-water hydrology, water quality, and water use. The Fox River is the primary surface-water body in Kendall County and is used for both wastewater disposal and as a drinking-water supply upstream of the county. Water from the Fox River requires pretreatment for use as drinking water, but the river is a potentially viable additional source of water for the county. Glacial drift aquifers capable of yielding sufficient water for municipal supply are expected to be present in northern Kendall County, along the Fox River, and in the Newark Valley and its tributaries. Glacial drift aquifers capable of yielding sufficient water for residential supply are present in most of the county, with the exception of the southeastern portion. Volatile organic compounds and select trace metals and pesticides have been detected at low concentrations in glacial drift aquifers near waste-disposal sites. Agricultural-related constituents have been detected infrequently in glacial drift aquifers near agricultural areas. However, on the basis of the available data, widespread, consistent problems with water quality are not apparent in these aquifers. These aquifers are a viable source for additional water supply, but would require further characterization prior to full development. The shallow bedrock aquifer is composed of the sandstone units of the Ancell Group, the Prairie du Chien Group, the Galena-Platteville dolomite, the Maquoketa Group, and the Silurian dolomite where these units are at the bedrock surface. The availability of water from the shallow bedrock aquifer depends primarily on the geologic unit utilized. The Silurian dolomite, Galena-Platteville dolomite, and Ancell Group can yield sufficient water for residential and municipal supply in at least some parts of the county. The Cambrian-Ordovician aquifer system is composed of the most widespread, productive aquifers in northern Illinois and is used for water supply by a number of municipalities and industrial facilities. Water levels in the aquifer system have declined by as much as 600 feet in Kendall County and the aquifer frequently contains concentrations of radium above established health guidelines.

Scientific Investigations Report

Nutrient loads in the Lost River and Klamath River Basins, south-central Oregon and northern California, March 2012–March 2015

The U.S. Geological Survey and Bureau of Reclamation collected water-quality data from March 2012 to March 2015 at locations in the Lost River and Klamath River Basins, Oregon, in an effort to characterize water quality and compute a nutrient budget for the Bureau of Reclamation Klamath Reclamation Project. The study described in this report resulted in the following significant findings: Total phosphorus (TP), total nitrogen (TN), 5-day biochemical oxygen demand (BOD5), and 5-day carbonaceous biochemical oxygen demand (CBOD5) loads, calculated using the U.S. Geological Survey LOADEST software package at the upper and lower boundaries of the Klamath Reclamation Project, indicated higher loads at the upper boundary on the southern end of Upper Klamath Lake upstream of the Bureau of Reclamation A Canal diversion compared to the lower boundary on the Klamath River downstream of Keno Dam. Accounting for the diversion of loads down A Canal, BOD5 and CBOD5 loads decreased between these two sites during irrigation season, indicating that the Klamath Reclamation Project is not a large source of oxygen-demanding material and that much of the oxygen demand at study site FMT, the northern boundary of the study area, has been expressed by the time the same water passes through site KRK, the southern boundary of the study area. An evaluation of the nutrient balance along the Klamath River flowpath from sites FMT to KRK indicated that, during irrigation season in the 3 years of the study period (March 2012–March 2015), more loads of TP, TN, BOD5, and CBOD5 were being diverted from the Klamath River than were being added to the Klamath River from the combination of Klamath Straits Drain, regulated point sources along the Klamath River, and internal loading from the bottom sediments in the river. By contrast, during non-irrigation seasons, more loads were added to the Klamath River than were diverted through Ady and North Canals, and this difference primarily was due to additional loads to the river from the Lost River Diversion Channel. At the Lost River Diversion Channel, BOD5 loads were higher during irrigation season than non-irrigation season in all three study years owing to the high concentrations of oxygen-demanding cyanobacterial biomass from the seasonal blooms of Aphanizomenon flos-aquae in the Klamath River and Upper Klamath Lake. The difference between the two seasons was particularly large in years 2 and 3, when the low flows of these two drought years resulted in smaller nonirrigation period loads than in year 1. CBOD5 loads also were higher during irrigation season in years 2 and 3 than during non-irrigation season, indicating that the largest oxygen demand was coming from senescence of Aphanizomenon flos-aquae cells that are present in the Klamath River during the summer. However, during irrigation season in year 1, CBOD5 loads were lower than in the non-irrigation season, which may indicate that at times high concentrations of ammonia or cellular organic nitrogen leaving Upper Klamath Lake contribute a large nitrogenous oxygen demand as well. The smallest loads were computed for the farthest upstream sites in the Lost River Basin, suggesting that the upper Lost River Basin does not contribute substantial loads of TP, TN, BOD5, and CBOD5 to the Klamath Reclamation Project. Median concentrations of BOD5 and CBOD5 were lowest among the upper Lost River Basin sites and highest at site PPD (however, this comparison is based on only four samples collected at site PPD over the 3-year study). Median concentrations of BOD5 and CBOD5 also were elevated at sites KSDH (6.60 and 4.70 milligrams per liter [mg/L], respectively) and KSD97 (4.47 and 3.45 mg/L, respectively). The highest maximum BOD5 and CBOD5 concentrations were reported at the Lost River Diversion Channel (39.0 and 26.5 mg/L, respectively) when water was flowing from the Klamath River toward the Klamath Reclamation Project, and site FMT (25.0 and 23.9 mg/L, respectively), the study site at the southern end of Upper Klamath Lake. Carbonaceous oxygen demand, as represented by CBOD5, typically dominated the composition of the samples at all sites. The highest concentrations of dissolved organic carbon were present at sites KSDH (the headworks of Klamath Straits Drain) and KSD97 (Klamath Straits drain before it enters the Klamath River), and PPD (outlet of Tule Lake). Median concentrations of TN and TP at the upper Lost River Basin sites in years 1 and 2 were variable, but site MCRV showed a smaller range of values in those years compared to the other upper Lost River Basins sites, and an overall lower median concentration during irrigation seasons in years 1 and 2, suggesting that Gerber Reservoir does not contribute high concentrations of nutrients to the Lost River during irrigation season. Total Maximum Daily Load (TMDL) load allocations for TP and TN in Klamath Straits Drain were exceeded in all three study years. BOD5 load allocations were exceeded in years 1 and 2, but not year 3. TMDL load allocations for TP were exceeded in the Lost River Diversion Channel for all 3 years. Load allocations for TN were exceeded in year 1, but not in years 2 and 3. BOD5 loads were less than the TMDL load allocation for all three study years. The dearth of samples collected at the Klamath Straits Drain just downstream of the Lower Klamath National Wildlife Refuge did not allow for direct assessment of the Klamath Straits Drain acting as a nutrient source or sink. TP, TN, BOD5, and CBOD5 loads estimated during the study period likely were smaller than long-term average conditions because of persistent drought conditions in the Upper Klamath Basin. The study results, therefore, fail to characterize loads from the Klamath Reclamation Project to the Klamath River that could be present in typical years, and suggest the need for load assessments during average or aboveaverage streamflow years.

California, Oregon

Estimating metal concentrations with regression analysis and water-quality surrogates at nine sites on the Animas and San Juan Rivers, Colorado, New Mexico, and Utah

The purpose of this report is to evaluate the use of site-specific regression models to estimate metal concentrations at nine U.S. Geological Survey streamflow-gaging stations on the Animas and San Juan Rivers in Colorado, New Mexico, and Utah. Downstream users could use these regression models to determine if metal concentrations are elevated and pose a risk to water supplies, agriculture, and recreation. Multiple linear-regression models were developed by relating metal concentrations in discrete water-quality samples to continuously monitored streamflow and surrogate parameters (specific conductance, pH, turbidity, and water temperature) collected at the U.S. Geological Survey stations. Models were developed for dissolved and total concentrations of aluminum, arsenic, cadmium, copper, iron, lead, manganese, and zinc using water-quality samples collected from 2005 to 2017 by several Federal, State, Tribal, and local agencies using different collection methods and analytical laboratories. Model performance varied but, in general, models for dissolved metals did not perform as well as those for total metals. Dissolved metals generally were correlated to specific conductance or streamflow and total metals generally were better correlated with turbidity. Explanatory variables in the models reflected hydrologic and geochemical processes within the basin. A larger number of regression models were statistically significant for the most upstream sites, where metal concentrations were elevated by drainage from abandoned mines and mineralized bedrock. Models generally did not perform as well at downstream sites, especially for dissolved metals, which occurred at lower concentrations than at the upstream sites. In the lower reaches of the rivers, the input of more alkaline water from tributaries and groundwater reduced metal solubility and diluted metal concentrations. The number and distribution of samples in the calibration datasets also may have been a factor in model development. At some sites on the San Juan River, calibration datasets were more limited and did not represent the full range of observed hydrologic and water-quality conditions, especially during storm events in summer and fall. Recommendations for model use are given based on estimates of model precision, biases, and adequacy of the calibration datasets in terms of the number of samples and representativeness of the observed range of streamflow and water-quality conditions.

Colorado, New Mexico, Utah

A review of surface energy balance models for estimating actual evapotranspiration with remote sensing at high spatiotemporal resolution over large extents

Many approaches have been developed for measuring or estimating actual evapotranspiration ( ET a ), and research over many years has led to the development of remote sensing methods that are reliably reproducible and effective in estimating ET a . Several remote sensing methods can be used to estimate ET a at the high spatial resolution of agricultural fields and the large extent of river basins. More complex remote sensing methods apply an analytical approach to ET a estimation using physically based models of varied complexity that require a combination of ground-based and remote sensing data, and are grounded in the theory behind the surface energy balance model. This report, funded through cooperation with the International Joint Commission, provides an overview of selected remote sensing methods used for estimating water consumed through ET a and focuses on Mapping Evapotranspiration at High Resolution with Internalized Calibration (METRIC) and Operational Simplified Surface Energy Balance (SSEBop), two energy balance models for estimating ET a that are currently applied successfully in the United States. The METRIC model can produce maps of ET a at high spatial resolution (30 meters using Landsat data) for specific areas smaller than several hundred square kilometers in extent, an improvement in practice over methods used more generally at larger scales. Many studies validating METRIC estimates of ET a against measurements from lysimeters have shown model accuracies on daily to seasonal time scales ranging from 85 to 95 percent. The METRIC model is accurate, but the greater complexity of METRIC results in greater data requirements, and the internalized calibration of METRIC leads to greater skill required for implementation. In contrast, SSEBop is a simpler model, having reduced data requirements and greater ease of implementation without a substantial loss of accuracy in estimating ET a . The SSEBop model has been used to produce maps of ET a over very large extents (the conterminous United States) using lower spatial resolution (1 kilometer) Moderate Resolution Imaging Spectroradiometer (MODIS) data. Model accuracies ranging from 80 to 95 percent on daily to annual time scales have been shown in numerous studies that validated ET a estimates from SSEBop against eddy covariance measurements. The METRIC and SSEBop models can incorporate low and high spatial resolution data from MODIS and Landsat, but the high spatiotemporal resolution of ET a estimates using Landsat data over large extents takes immense computing power. Cloud computing is providing an opportunity for processing an increasing amount of geospatial “big data” in a decreasing period of time. For example, Google Earth Engine TM has been used to implement METRIC with automated calibration for regional-scale estimates of ET a using Landsat data. The U.S. Geological Survey also is using Google Earth Engine TM to implement SSEBop for estimating ET a in the United States at a continental scale using Landsat data.

Scientific Investigations Report

Peak-flow and low-flow magnitude estimates at defined frequencies and durations for nontidal streams in Delaware

Reliable estimates of the magnitude of peak flows in streams are required for the economical and safe design of transportation and water conveyance structures. In addition, reliable estimates of the magnitude of low flows at defined frequencies and durations are needed for meeting regulatory requirements, quantifying base flows in streams and rivers, and evaluating time of travel and dilution of toxic spills. This report, in cooperation with the Delaware Department of Transportation and the Delaware Geological Survey, presents methods for estimating the magnitude of peak flows and low flows at defined frequencies and durations on nontidal streams in Delaware, at locations both monitored by streamflow-gage sites and ungaged. Methods are presented for estimating (1) the magnitude of peak flows for return periods ranging from 2 to 500 years (50-percent to 0.2-percent annual-exceedance probability), and (2) the magnitude of low flows as applied to 7-, 14-, and 30-consecutive day low-flow periods with recurrence intervals of 2, 10, and 20 years (50-, 10-, and 5-percent annual non-exceedance probabilities). These methods are applicable to watersheds that exhibit a full range of development conditions in Delaware. The report also describes StreamStats, a web application that allows users to easily obtain peak-flow and low-flow magnitude estimates for user-selected locations in Delaware. Peak-flow and low-flow magnitude estimates for ungaged sites are obtained using statistical regression analysis through a process known as regionalization, where information from a group of streamflow-gage sites within a region forms the basis for estimates for ungaged sites within the same region. Ninety-four streamflow-gage sites in and near Delaware with at least 10 years of nonregulated annual peak-flow data were used for the peak-flow regression analysis, a subset of the 121 sites for which peak-flow estimates were computed. These sites included both continuous-record streamflow-gage sites as well as partial record sites. Forty-five streamflow-gage sites with at least 10 years of nonregulated low-flow data available were used for the low-flow regression analyses, a subset of the 68 sites for which low-flow estimates were computed. Estimates for gaged sites are obtained by combining (1) the station peak-flow statistics (mean, standard deviation, and skew) and peak-flow estimates using the recent Bulletin 17C guidelines that incorporate the Expected Moments Algorithm with (2) regional estimates of peak-flow magnitude derived from regional regression equations and regional skew derived from sites with records greater than or equal to 35 years. Example peak-flow estimate calculations using the methods presented in the report are given for (1) ungaged sites, (2) gaged sites, (3) sites upstream or downstream from a gaged location, and (4) sites between gaged locations. Estimates for low-flow gaged sites are obtained by combining (1) the station low-flow statistics (mean, standard deviation, and skew) and low-flow estimates with (2) regional estimates of low-flow magnitude derived from regional regression equations. Example low-flow estimate calculations using the methods presented in the report are given for (1) ungaged sites, (2) gaged sites, (3) sites upstream or downstream from a gaged location, and (4) sites between gaged locations. A total of 54 sites in the Coastal Plain region were used to develop peak-flow regressions for the region and 40 sites were used for the Piedmont region. Similarly, 24 sites were used for low-flow regression equation development in the Coastal Plain, with 21 in the Piedmont. Peak and low-flow site inclusion in the Coastal Plain tended to be more restricted with tidal influence and ranges of basin characteristics, including drainage area, limiting regression equation development and application. Regional regression equations for peak flows and low flows, as applicable to ungaged sites in the Piedmont and Coastal Plain Physiographic Provinces in Delaware, are presented. Peak-flow regression equations used variables that quantified drainage area, basin slope, percent area with well-drained soils, percent area with poorly drained soils, impervious area, and percent area of surface water storage in estimating peak-flow estimates, whereas low-flow regression equations used only drainage area and percent poorly drained soils in the estimation of low flows. Average standard errors for peak-flow regressions tended to be lower than those for low- flow regressions, with lower errors in the Piedmont region for both peak- and low-flow regressions. For peak-flow estimates, a sensitivity analysis of Piedmont regression equation estimates to changes in impervious area is also presented. Additional topics associated with the analyses performed during the study are discussed, including (1) the availability and description of 32 basin and climatic characteristics considered during the development of the regional regression equations; (2) the treatment of increasing trends in the annual peak-flow series identified at 18 gaged sites and inclusion in or exclusion from the regional analysis; (3) regional skew analysis and determination of regression regions; (4) sample adjustments and removal of sites owing to regulation and redundancy; and (5) a brief comparison of peak- and low-flow estimates at gages used in previous studies.

Delaware

Use of sediment rating curves and optical backscatter data to characterize sediment transport in the Upper Yuba River watershed, California, 2001-03

Sediment transport in the upper Yuba River watershed, California, was evaluated from October 2001 through September 2003. This report presents results of a three-year study by the U.S. Geological Survey, in cooperation with the California Ecosystem Restoration Program of the California Bay-Delta Authority and the California Resources Agency. Streamflow and suspended-sediment concentration (SSC) samples were collected at four gaging stations; however, this report focuses on sediment transport at the Middle Yuba River (11410000) and the South Yuba River (11417500) gaging stations. Seasonal suspended-sediment rating curves were developed using a group-average method and non-linear least-squares regression. Bed-load transport relations were used to develop bed-load rating curves, and bed-load measurements were collected to assess the accuracy of these curves. Annual suspended-sediment loads estimated using seasonal SSC rating curves were compared with previously published annual loads estimated using the Graphical Constituent Loading Analysis System (GCLAS). The percent difference ranged from -85 percent to +54 percent and averaged -7.5 percent. During water year 2003 optical backscatter sensors (OBS) were installed to assess event-based suspended-sediment transport. Event-based suspended-sediment loads calculated using seasonal SSC rating curves were compared with loads calculated using calibrated OBS output. The percent difference ranged from +50 percent to -369 percent and averaged -79 percent. The estimated average annual sediment yield at the Middle Yuba River (11410000) gage (5 tons/mi2) was significantly lower than that estimated at the South Yuba River (11417500) gage (14 tons/mi2). In both rivers, bed load represented 1 percent or less of the total annual load throughout the project period. Suspended sediment at the Middle Yuba River (11410000) and South Yuba River (11417500) gages was typically greater than 85 percent silt and clay during water year 2003, and sand concentrations at the South Yuba River (11417500) gage were typically higher than those at the Middle Yuba River (11410000) gage for a given streamflow throughout the three year project period. Factors contributing to differences in sediment loads and grain-size distributions at the Middle Yuba River (11410000) and South Yuba River (11417500) gages include contributing drainage area, flow diversions, and deposition of bed-material-sized sediment in reservoirs upstream of the Middle Yuba River (11410000) gage. Owing to its larger drainage area, higher flows, and absence of man-made structures that restrict sediment movement in the lower basin, the South Yuba River transports a greater and coarser sediment load.

Scientific Investigations Report

Peak-flow variability, peak-flow informational needs, and consideration of regional regression analyses in managing the crest-stage gage network in Montana

The U.S. Geological Survey (USGS), in cooperation with the Montana Department of Transportation (MDT), has operated a crest-stage gage (CSG) network in Montana to collect peak-flow data since 1955. The CSG network is vital to collecting peak-flow data on small drainage basins that typically are not addressed by continuous streamflow operations. Discussions between USGS and MDT identified a need for evaluating the CSG network to allow for better decision making in the management of the network. The purpose of this report is to (1) generally describe peak-flow variability in Montana, (2) assess peak-flow informational needs relevant to MDT activities, and (3) consider the characteristics of the active CSG network in relation to addressing the informational needs. The evaluation of the CSG network is intended to assist in prioritization for discontinuation of CSGs and other activities involving changes to the CSG network. Peak-flow variability was investigated by analysis of selected peak-flow characteristics of 659 unregulated streamgages in or near Montana. A generalized peak-flow variability index ( PFVI ) was developed to provide large-scale representation of peak-flow variability in Montana. For unregulated Montana streamgages, PFVI generally monotonically decreases with increasing drainage area, although there is somewhat large (but generally consistent) variability about the locally weighted scatterplot smooth line. Presumably, highly variable small-scale hydroclimatic processes are integrated with increasing drainage area such that variability in many hydrologic characteristics is reduced. PFVI also decreases with increasing mean basin elevation and mean annual precipitation. Presumably, higher elevation and wetter hydroclimatic settings in Montana contribute to reduced variability in hydrologic characteristics. Intuitively, PFVI might be expected to generally decrease with increasing years of record because the standard deviation might typically be expected to decrease with increasing sample size. However, relations among PFVI and years of record are more complex and variable than drainage area, elevation, and precipitation. PFVI variably increases from 10 to about 40 years of record and then generally monotonically decreases from about 40 to about 105 years of record. Relations among PFVI and the years of record might be confounded by effects of drainage area because streamgages with long periods of record (greater than about 60 years) generally have large drainage areas (greater than about 100 square miles). The relations between PFVI and drainage area, mean basin elevation, mean annual precipitation, and years of record substantially differ among the eight hydrologic regions in Montana. As such, the PFVI relations were further investigated within each hydrologic region. A major use of peak-flow information by MDT is for design of road and highway infrastructure, including bridges, culverts, and roadside drainage ditches. As such, basin characteristics (including drainage area, mean basin elevation, and mean annual precipitation) of the Montana streamgage network (735 regulated and unregulated streamgages) were statistically investigated in relation to basin characteristics of 12,639 road and stream intersections in Montana. Both regulated and unregulated streamgages were investigated because the road and stream intersections are on both regulated and unregulated streams. Exploratory analyses indicated that the various relations substantially differ among the hydrologic regions. As such, the relations between the Montana streamgage network and the road and stream intersections were further investigated within each hydrologic region. An important objective of the CSG network is to provide data for developing regional regression equations (RREs) for estimating frequencies at ungaged sites in Montana. Various characteristics of the RREs substantially differ among the eight hydrologic regions in Montana. As such, the RRE characteristics were further investigated within each hydrologic region. For each of the eight hydrologic regions, various characteristics of peak-flow variability, peak-flow informational needs, and regional regression analyses were investigated in detail. Possible shortcomings of the streamgage network in each hydrologic region are identified and possible future improvements to the CSG network are presented.

Montana

Trends in the quality of water in New Jersey streams, water years 1971–2011

In a study conducted by the U.S. Geological Survey in cooperation with the New Jersey Department of Environmental Protection and the Delaware River Basin Commission, trend tests were conducted on selected water-quality characteristics measured at stations on streams in New Jersey during selected periods over water years 1971‒2011. Tests were conducted on 3 nutrients (total nitrogen, filtered nitrate plus nitrite, and total phosphorus) at 28 water-quality stations. At 4 of these stations, tests were also conducted on 3 measures of major ions (specific conductance, filtered chloride, and total dissolved solids). Two methods were used to identify trends—Weighted Regressions on Time, Discharge, and Season (WRTDS) models and seasonal rank-sum tests. For this report, the use of WRTDS models included the use of the WRTDS Bootstrap Test (WBT). WRTDS models identified trends in flow-normalized annual concentrations and flow-normalized annual fluxes over water years 1980‒2011 and 2000‒11 for each nutrient, filtered chloride, and total dissolved solids. WRTDS models were developed for each nutrient at the 20 or 21 stations at which streamflow was measured or estimated. Trends in nutrient concentration were reported for these stations; trends in nutrient fluxes were reported only for 15–17 of these stations. The results of WRTDS models for water years 1980‒2011 identified more stations with downward trends in concentrations of either total nitrogen or total phosphorus than upward trends. For total nitrogen, there were downward trends at 9 stations and an upward trend at 1 station. For total phosphorus, there were downward trends at 8 stations and an upward trend at 1 station. For filtered nitrate plus nitrite, there were downward trends at 6 stations and upward trends at 6 stations. The result of the trend test in flux for a selected nutrient at a selected station (downward trend, no trend, or upward trend) usually matched the trend result in concentration. Seasonal rank-sum tests, the second method used, identified step trends in water-quality measured in different decades—1970s, 1980s, 1990s, and 2000s. Tests were conducted on all nutrients at 28 stations and on all measures of major ions at the 4 selected stations. Results of seasonal rank-sum tests between the 1980s and the 2000s identified more stations with downward trends in concentrations of total nitrogen (14) than stations with upward trends (2) and more stations with downward trends in concentrations of total phosphorus (18) than stations with upward trends (1). A combined dataset of trend results for concentrations over water years 1980‒2011 was created from the results of the two tests for the period. Results of WRTDS models were included in this combined dataset, if available. Otherwise, the results of the seasonal rank-sum tests between water-quality characteristics measured in the 1980s and 2000s were included. Trend results over water years 1980‒2011 in the combined dataset show that few of the 28 stations had upward trends in concentrations of either total nitrogen or total phosphorus. There were only 2 stations with upward trends in total nitrogen concentration and 1 station with an upward trend in total phosphorus concentration. Results for filtered nitrate plus nitrite show about the same number of stations with upward trends (9) as stations with downward trends (7). Results for all measures of major ions show upward trends at the four stations tested.

New Jersey

Natural resource mitigation, adaptation and research needs related to climate change in the Great Basin and Mojave Desert

This report synthesizes the knowledge, opinions, and concerns of many Federal and State land managers, scientists, stakeholders, and partners from a workshop, held at the University of Nevada, Las Vegas, on April 20-22, 2010. Land managers, research scientists, and resource specialists identified common concerns regarding the potential effects of climate change on public lands and natural resources in the Great Basin and Mojave Desert and developed recommendations for mitigation, adaptation, and research needs. Water and, conversely, the effects of drought emerged as a common theme in all breakout sessions on terrestrial and aquatic species at risk, managing across boundaries, monitoring, and ecosystem services. Climate change models for the southwestern deserts predict general warming and drying with increasing precipitation variability year to year. Scientists noted that under these changing conditions the past may no longer be a guide to the future in which managers envision increasing conflicts between human water uses and sustaining ecosystems. Increasing environmental stress also is expected as a consequence of shifting ecosystem boundaries and species distributions, expansion of non-native species, and decoupling of biotic mutualisms, leading to increasingly unstable biologic communities. Managers uniformly expressed a desire to work across management and agency boundaries at a landscape scale but conceded that conflicting agency missions and budgetary constraints often impede collaboration. More and better science is needed to cope with the effects of climate change but, perhaps even more important is the application of science to management issues using the methods of adaptive management based on long-term monitoring to assess the merits of management actions. Access to data is essential for science-based land management. Basic inventories, spatial databases, baseline condition assessments, data quality assurance, and data sharing were identified as top information priorities by all participants at this workshop. Optimizing the utility of ecosystem monitoring data will require standardizing monitoring protocols across agencies. Better communication among researchers and managers and cooperation through partnerships to manage resources across boundaries were emphasized as necessary for adapting to changing climatic conditions. However, even these strategies may be insufficient unless policy mandates, agency missions, and funding are coordinated at a high level.

Scientific Investigations Report

Stream classification and gravel-bar inventory for Buffalo National River and Ozark National Scenic Riverways

The data summarized in this report provide a baseline characterization of the physical attributes of the riverine ecosystems in two landscapes managed by the National Park Service—Ozark National Scenic Riverways, Missouri, and Buffalo National River, Arkansas—to inform understanding and management of aquatic habitat. The study utilized a basin-scale approach and consisted of two components: a basin-scale channel classification and a longitudinal inventory of gravel bars. We evaluated the Jacks Fork and Current River in Ozark National Scenic Riverways and the main stem Buffalo River in Buffalo National River. The primary objective of the study was to characterize geomorphic patterns that affect channel stability and rates of geomorphic change in both national park units. Findings may be used to inform understanding of the distribution and availability of aquatic habitat. For the basin-scale channel classification, we performed exploratory statistical analyses using nine geomorphic variables (channel width, standard deviation in channel width, valley width, distance to valley wall, confinement, bar area, bluff area, braid index, and sinuosity). Each metric was quantified along the length of the river system at 200-meter intervals. We then performed a cluster analysis for each river using a subset of variables, resulting in 2 to 5 distinct geomorphic classes depending on criteria used for determining number of clusters. Longitudinal patterns in clusters vary for each river system but reflect a combination of landscape factors including valley width, influence of tributaries, and lithology, which affect channel stability and aquatic habitat. We developed a longitudinal inventory of gravel bars by quantifying the area of gravel bars from a series of imagery in each park. In Ozark National Scenic Riverways we analyzed five time periods, with the earliest being 1992 and most recent being 2014. In Buffalo National River, we also analyzed five series of aerial imagery, ranging from 1982 to 2013. The analysis indicated a general decrease in gravel storage in upstream reaches of each river evaluated, accompanied by a general increase in storage farther downstream. Local patterns in gravel-bar area are associated with longitudinal patterns in geomorphic setting, such as valley geometry and channel width, that affect depositional patterns and sediment storage at the reach scale.

Arkansas, Missouri

Stream sediment sources in Medicine Creek, northern Missouri and southern Iowa

This report presents the results of a cooperative study by the U.S. Geological Survey and Missouri Department of Natural Resources to quantify sediment transport source contributions in the Medicine Creek drainage basin. Understanding relative source contributions provides valuable information for selecting the conservation practices that may be most effective in reducing sediment and sediment-associated nutrient transport in the Medicine Creek drainage basin and similar areas of the Lower Grand River drainage basin. Sediment samples were collected from potential contributing areas (source samples) and from fluvial-transported samples (target samples). Source sample types included streambanks, row crop fields, and a combined pastures and forests category. Samples were analyzed for particle size and quantity of carbon, nitrogen, stable isotopes of carbon and nitrogen, and 49 mineral elements as potential tracers. Results for the carbon stable isotope ratio of carbon-13/carbon-12 (δ 13 C) and concentrations of total carbon, total nitrogen, calcium, potassium, and copper were selected by discriminant function analysis as the best combination of multiple tracers to differentiate each source type. The discriminant function analysis poorly differentiated pastures and forests, so these source types were combined. The sources defined by the discriminant function analysis were then used in an unmixing model to apportion sources for each target sample. In the study area, transported sediment was predominantly bank sediment, with an overall average of 86.9 percent of suspended-sediment samples and depositional streambed samples attributed to bank material. Suspended-sediment samples from the mainstem of Medicine Creek were dominated by bank sediments (average of 95.8 percent), and depositional streambed samples from throughout the drainage basin had more variable source contributions with an average of 71.1 percent attributed to bank material. The relative importance of upland sources (row crop fields and the combined pastures and forests category) varied seasonally and with streamflow but was not related to land use or drainage basin size. Relative contributions from upland sources were greater in the summer through winter rather than spring and during lower streamflow, though this may be driven by the seasonality of streamflow. These results indicate management practices that reduce bank erosion could be effective strategies for managing the dominant source of sediment and sediment-associated phosphorus.

Iowa, Missouri

Hydrogeologic framework and salinity distribution of the Floridan aquifer system of Broward County, Florida

Concerns about water-level decline and seawater intrusion in the surficial Biscayne aquifer, currently the principal source of water supply to Broward County, prompted a study to refine the hydrogeologic framework of the underlying Floridan aquifer system to evaluate its potential as an alternative source of supply. This report presents cross sections that illustrate the stratigraphy and hydrogeology in eastern Broward County; maps of the upper surfaces and thicknesses of several geologic formations or units within the Floridan aquifer system; and maps of two of the potentially productive water-bearing zones within the system, the Upper Floridan aquifer and the Avon Park permeable zone. An analysis of data on rock depositional textures, associated pore networks, and flow zones in the Floridan aquifer system shows that groundwater moves through the system in two ways. These data support a conceptual, dual-porosity model of the system wherein groundwater moves either as concentrated flow in discrete, thin bedding-plane vugs or zones of vuggy megaporosity, or as diffuse flow through rocks with primarily interparticle and moldic-particle porosity. Because considerable exchange of groundwater may occur between the zones of vuggy and matrix-dominated porosity, understanding the distribution of that porosity and flow zone types is important to evaluating the suitability of the several units within the Floridan aquifer system for managing the water through practices such as aquifer storage and recovery (ASR). The salinity of the water in the Floridan aquifer system is highest in the central part of the study area, and lower toward the north and south. Although salinity generally increases with depth, in the western part of the study area a zone of relatively high saline water is perched above water of lower salinity in the underlying Avon Park permeable zone. Overall, the areas of highest salinity in the aquifer system coincide with those with the lowest estimated transmissivity, so that the occurrence of perched saline water in the system may be the consequence of incompletely flushed connate water or intruded seawater. A seismic reflection profile along the Hillsboro Canal, at the northern edge of the study area, shows seven seismic-sag structures that are interpreted as downward deformation of overlying strata into collapsed deep cave systems. These structures may compromise the integrity of the confinement created by the underlying strata by allowing upconing of saline water from depth, which has implications for successful application of ASR and use of the Floridan aquifer system as an alternative water supply.

Florida

Groundwater conditions in Georgia, 2015–16

The U.S. Geological Survey collects groundwater data and conducts studies to monitor hydrologic conditions, define groundwater resources, and address problems related to water supply, water use, and water quality. In Georgia, water levels were monitored continuously at 157 wells during calendar years 2015 and 2016. Because of missing data or short periods of record (less than 5 years) for several of these wells, data for 147 wells are presented in this report. These wells include 15 in the surficial aquifer system, 18 in the Brunswick aquifer system and equivalent sediments, 59 in the Upper Floridan aquifer, 13 in the Lower Floridan aquifer and underlying units, 9 in the Claiborne aquifer, 1 in the Gordon aquifer, 8 in the Clayton aquifer, 16 in the Cretaceous aquifer system, 2 in Paleozoic-rock aquifers, and 6 in crystalline-rock aquifers. Data from the well network indicate that water levels generally rose during the 10-year period from 2007 through 2016, with water levels rising in 105 wells and declining in 31 wells; insufficient data prevented determination of a 10-year trend in 11 wells. Water levels declined over the long-term period of record at 80 wells, increased at 62 wells, and remained relatively constant at 5 wells. In addition to continuous water-level data, periodic water-level data were collected and used to construct potentiometric-surface maps for the Upper Floridan aquifer in the Brunswick–Glynn County area during October 2015 and October 2016 and in the Albany–Dougherty County area during December 2015 and November and December 2016. Periodic water-level measurements were also collected and used to construct potentiometric-surface maps for the Cretaceous aquifer system in the Augusta–Richmond County area during July 2015 and June 2016. In general, water levels in the Upper Floridan aquifer were higher during 2015 than during 2016 in the Brunswick–Glynn County and Albany–Dougherty County areas due to higher precipitation during 2015. Water levels were lower, however, during 2015 than during 2016 in the Cretaceous aquifer system in the Augusta–Richmond County area. In the Brunswick area, maps showing the chloride concentration of water in the Upper Floridan aquifer constructed using data collected from 33 wells during October 2015 and from 30 wells during October 2016 indicate that chloride concentrations remained above the U.S. Environmental Protection Agency’s secondary drinking-water standard in an approximately 2-square-mile area. During calendar years 2015–16, chloride concentrations generally were similar to those measured during 2012–14; however, some wells did show an increase in chloride concentration, likely due to increases in pumping.

Georgia

Comparison of surrogate models to estimate pesticide concentrations at six U.S. Geological Survey National Water Quality Network sites during water years 2013–18

During water years 2013–18, the U.S. Geological Survey National Water-Quality Assessment Project sampled the National Water Quality Network for Rivers and Streams year-round and reported on 221 pesticides at 72 sites across the United States. Pesticides are difficult to measure, their concentrations often represent discrete snapshots in time, and capturing peak concentrations is expensive. Three types of regression models were developed to estimate concentrations for two selected pesticides at each of six National Water Quality Network for Rivers and Streams sites. The regression models used continuously measured streamflow and water-quality properties (differing combinations of pH, specific conductance, turbidity, and water temperature); discrete water-quality samples analyzed for atrazine, azoxystrobin, bentazon, bromacil, imidacloprid, simazine, and triclopyr; and time as an additional explanatory variable for seasonality. The modeling approaches included (1) a standard regression that included surrogates (differing combinations of pH, specific conductance, turbidity, and water temperature) and periodic functions (sine-cosine) of pesticide application use as predictor variables; (2) the seasonal wave with flow adjustment model that included a seasonal component and flow anomalies but excluded surrogates; and (3) the seasonal wave with flow adjustment model that included a seasonal component, flow anomalies, and surrogates. Models were evaluated using three measures of model performance: generalized coefficient of determination (generalized R 2 ), Akaike’s Information Criteria, and scale (the estimated standard deviation of the tobit regression error term). Because of low observation numbers, results from this study can be considered a pilot effort with the possibility that some models are overfit. In all cases, estimated pesticide concentrations modeled with base SEAWAVE-Q were better than the standard surrogate regression models; all 39 generalized R 2 values increased by 3–56 percent (median of 25 percent) when compared to the standard surrogate regression models, and all Akaike’s Information Criteria and scale values decreased. The addition of surrogate variables such as pH, specific conductance, turbidity, and water temperature to the base SEAWAVE-Q model to improve estimates of pesticide concentrations resulted in only modest improvements; generalized R 2 values increased by only 0–10 percent (median of 3 percent). In some instances, combinations of the surrogates produced more appreciative improvements in model results, but in those instances, we hypothesize that the surrogates correlated with some unknown measure that directly relates to pesticide transport.

Indiana, Kansas, Missouri, North Dakota, Oregon

Introduction and methods of analysis for peak streamflow trends and their relation to changes in climate in Illinois, Iowa, Michigan, Minnesota, Missouri, Montana, North Dakota, South Dakota, and Wisconsin

Flood-frequency analysis, also called peak-flow frequency or flood-flow frequency analysis, is essential to water resources management applications including critical structure design and floodplain mapping. Federal guidelines for doing flood-frequency analyses are presented in a U.S. Geological Survey Techniques and Methods Report known as Bulletin 17C. A basic assumption within Bulletin 17C is that for drainage basins without major hydrologic alterations, statistical properties of the distribution of annual peak streamflows (peak flows) are stationary; that is, the mean, variance, and skew are constant. The stationarity assumption has been widely accepted within the flood-frequency community; however, a better understanding of long-term climatic persistence and concerns about potential climate change and land-use change has caused a reexamination of the stationarity assumption. Flood-frequency analyses that do not incorporate observed trends and abrupt changes may result in a poor representation of the true flood risk. Bulletin 17C does not offer guidance for incorporating nonstationarities when estimating floods, and it describes a need for studies that incorporate changing climate or basin characteristics. In response to this need and a history of concern regarding nonstationarity peak flows in the region, this study was done by the U.S. Geological Survey, in cooperation with the Departments of Transportation of Illinois, Iowa, Michigan, Minnesota, Missouri, South Dakota, and Wisconsin; the Montana Department of Natural Resources and Conservation; and the North Dakota Department of Water Resources, to assess potential nonstationarity in peak flows in the north-central United States. This chapter summarizes the methods used to detect hydroclimatic changes in peak-flow data in the study region. A wide range of analyses and statistical approaches are applied to document the primary mechanisms controlling floods and characterize temporal changes in hydroclimatic variables and peak flow. Four periods were selected for analysis of peak flow, daily streamflow, and climate data. The periods are (1) a 100-year period, 1921–2020; (2) a 75-year period, 1946–2020; (3) a 50-year period, 1971–2020; and (4) a 30-year period, 1991–2020. The climate data consist of monthly time series estimates of temperature, precipitation, potential evapotranspiration, actual evapotranspiration, snowfall, soil moisture storage, snow water equivalent, and runoff on a 3.1-mile by 3.1-mile grid for the conterminous United States. Statistical and graphical analyses were used to investigate potential changes in hydrology and climate. The starting point for these analyses was the initial data analysis of peak flow described in Bulletin 17C, which includes plotting the peak flow and checking for autocorrelation, monotonic trends, and changes points. Analyses were added to examine additional features in the data. To examine potential causal drivers of changes, the climate data were analyzed graphically and statistically. Results are provided in a U.S. Geological Survey data release. The study limitations are documented for users of the results.

Iowa, Illinois, Michigan, Missouri, Minnesota, Mon

Magnitude of flood flows for selected annual exceedance probabilities in Rhode Island through 2010

Heavy persistent rains from late February through March 2010 caused severe widespread flooding in Rhode Island that set or nearly set record flows and water levels at many long-term streamgages in the State. In response, the U.S. Geological Survey, in partnership with the Federal Emergency Management Agency, conducted a study to update estimates of flood magnitudes at streamgages and regional equations for estimating flood flows at ungaged locations. This report provides information needed for flood plain management, transportation infrastructure design, flood insurance studies, and other purposes that can help minimize future flood damages and risks. The magnitudes of floods were determined from the annual peak flows at 43 streamgages in Rhode Island (20 sites), Connecticut (14 sites), and Massachusetts (9 sites) using the standard Bulletin 17B log-Pearson type III method and a modification of this method called the expected moments algorithm (EMA) for 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probability (AEP) floods. Annual-peak flows were analyzed for the period of record through the 2010 water year; however, records were extended at 23 streamgages using the maintenance of variance extension (MOVE) procedure to best represent the longest period possible for determining the generalized skew and flood magnitudes. Generalized least square regression equations were developed from the flood quantiles computed at 41 streamgages (2 streamgages in Rhode Island with reported flood quantiles were not used in the regional regression because of regulation or redundancy) and their respective basin characteristics to estimate magnitude of floods at ungaged sites. Of 55 basin characteristics evaluated as potential explanatory variables, 3 were statistically significant—drainage area, stream density, and basin storage. The pseudo-coefficient of determination (pseudo- R2 ) indicates these three explanatory variables explain 95 to 96 percent of the variance in the flood magnitudes from 20- to 0.2-percent AEPs. Estimates of uncertainty of the at-site and regression flood magnitudes are provided and were combined with their respective estimated flood quantiles to improve estimates of flood flows at streamgages. This region has a long history of urban development, which is considered to have an important effect on flood flows. This study includes basins that have an impervious area ranging from 0.5 to 37 percent. Although imperviousness provided some explanatory power in the regression, it was not statistically significant at the 95-percent confidence level for any of the AEPs examined. Influence of urbanization on flood flows indicates a complex interaction with other characteristics that confounds a statistical explanation of its effects. Standard methods for calculating magnitude of floods for given AEP are based on the assumption of stationarity, that is, the annual peak flows exhibit no significant trend over time. A subset of 16 streamgages with 70 or more years of unregulated systematic record indicates all but 4 streamgages have a statistically significant positive trend at the 95-percent confidence level; three of these are statistically significant at about the 90-percent confidence level or above. If the trend continues linearly in time, the estimated magnitude of floods for any AEP, on average, will increase by 6, 13, and 21 percent in 10, 20, and 30 years' time, respectively. In 2010, new peaks of record were set at 18 of the 21 active streamgages in Rhode Island. The updated flood frequency analysis indicates the peaks at these streamgages ranged from 2- to 0.2-percent AEP. Many streamgages in the State peaked at a 0.5- and 0.2-percent AEP, except for streamgages in the Blackstone River Basin, which peaked from a 4- to 2-percent AEP.

Rhode Island