The hydrologic regime of rivers and streams is a major determinant of habitat quality for fish and aquatic invertebrates. Long-term streamflow data were compiled and multidecadal streamflow trends and ecological flow (EFlow) statistics were calculated in support of the United States Geological Survey (USGS) Chesapeake Bay Science Initiative toward understanding fish habitat and health in the Chesapeake Bay Watershed (CBWS). A dataset comprising all streamgages (n = 409) reporting daily means of streamflow within the CBWS and remaining active as of September 30, 2018 (the end of Water Year [WY] 2018), independent of streamgage installation date, was retrieved from the USGS National Water Information System (NWIS). This dataset was then subset to include only those streamgages with a contiguous timeseries of streamflow data from a start date no earlier than April 1, 1939 (Climate Year [CY] 1940) and no later than October 1, 1999 (WY 2000). The R packages “EGRET” and "Eflowstats" were utilized together to determine streamflow trends and EFlow statistics from the subset (n = 243). Trends and EFlows were computed for the ranges 1940-1969 (n = 90), 1970-1999 (n = 167), and 2000-2018 (n = 243). Streamflow trends were computed for eight annual metrics (1-, 7- and 30-day minima [CY] and maxima [WY], mean and median [WYs]). These streamflow trends provide context for the 178 EFlow statistics (WY) which have been designated to characterize the magnitude, frequency, and duration of extreme high and low flows, the timing of seasonal flows, and the consistency of the historic regime. Files herein include the following Child Items: (1) a table summarizing streamflow trends for three time periods at a minimum of 90 and maximum of 243 streamgages and 500 time-series plots graphically representing those trends; (2) a table summarizing EFlow statistics and the change between each statistic for three time periods at a minimum of 90 and maximum of 243 streamgages; and (3) a GIS shapefile of the original 409 USGS streamgage locations, complete with NWIS attributes, active within the CBWS through September 30, 2018.

This dataset contains mean-annual metrics quantifying intermittent streamflow, climate, topography, land cover and geology for 540 USGS GAGES-2 watersheds experiencing non-perennial flow.

This data release contains trend results computed on the basis of modeled and observed daily streamflows at 502 reference gages across the conterminous U.S. from October 1, 1983 through September 30, 2016. Modeled daily streamflows were computed using the deterministic Precipitation Runoff Modeling System (PRMS), and five statistical techniques: Nearest-Neighbor Drainage Area Ratio (NNDAR), Map-Correlation Drainage Area Ratio (MCDAR), Ordinary Kriging of the logarithms of discharge per unit area (OKDAR), Nearest-Neighbor nonlinear spatial interpolation using flow duration curves (NNQPPQ), and Map-Correlation nonlinear spatial interpolation using flow duration curves (MCQPPQ). Observed daily streamflow data for the study gages were retrieved from the National Water Information System (NWIS). Study gages were selected from among Hydro-Climatic Data Network 2009 (HCDN-2009) gages in the GAGES-II dataset considered to be minimally affected by regulation, diversion, mining, or other anthropogenic activities. Results include trends in annual and monthly means, annual percentiles (1, 5, 10, 25, 50, 75, 90, 95, 99), annual 1-day high, 3-day high, and 7-day low, and annual snowmelt-related runoff timing for a subset of snowmelt dominated basins. Bias and volumetric efficiency statistics between observed and modeled streamflows also are provided.

This data set contains U.S. Geological Survey (USGS) streamgage identification numbers, begin and end years of the periods of streamflow record tested, Sen slope trends in the annual minimum 7-day streamflow for the period of record tested, the p-values (significance) of the trends, and the trend Sen slopes standardized by the standard deviations of the residual errors defined as the difference between observations and the Sen slope lines, for 174 USGS streamgages with 56 to 75 years of record in the Chesapeake Bay watershed, Mid-Atlantic U.S.

Ecological flow (EFlow) statistics have been designated to characterize the magnitude, frequency, and duration of extreme high- and low-flows, the timing of seasonal flows, and the consistency of the historic regime. This Child Item contains a table of 178 EFlows for the time periods 1940-1969, 1970-1999, and 2000-2018, with absolute and percent change between periods, where applicable. Statistics were computed by Water Year (WY) for all 178 metrics and absolute and percent change were calculated by comparing metrics between combinations of two of the three time periods (1940-1969 and 1970-1999; 1940-1969 and 2000-2018; 1970-1999 and 2000-2018). Streamgages from the original dataset (n = 409) were excluded from one or more time periods of analysis because of extensive data gaps that would yield incomplete EFlows; therefore, stations were indexed into the earliest possible time period relative to their installation date (for example, a streamgage with an operating start year of 1958 would be included in the analysis for the time periods 1970-1999 and 2000-2018), which resulted in different sample sizes for each period: 1940-1969 (n = 90), 1970-1999 (n = 167), and 2000-2018 (n = 243). Similarly, multiple stations were wholly excluded because of frequent discontinuities in the daily mean streamflow through all three time periods. Finally, a streamgage must have fallen within at least two time periods to have a change value. As such, not all stations are represented in the change analysis (change between 1940-1969 and 1970-1999 [n = 90]; change between 1940-1969 and 2000-2018 [n = 90]; change between 1970-1999 and 2000-2018 [n = 167]).

This data release contains daily time series estimates of natural streamflow at 5,439 GAGES-II non-reference streamgages in 19 study regions across the conterminous United States from October 1, 1980 through September 30, 2017, using five statistical techniques: nearest-neighbor drainage area ratio (NNDAR), map-correlation drainage area ratio (MCDAR), nearest-neighbor nonlinear spatial interpolation using flow duration curves (NNQPPQ), map-correlation nonlinear spatial interpolation using flow duration curves (MCQPPQ), and ordinary kriging of the logarithms of discharge per unit area (OKDAR). NNDAR, MCDAR, NNQPPQ, and MCQPPQ estimates were computed following methods described in Farmer and others (2014), with updates to the flow-duration curve modeling which is described in Over and others (2018). OKDAR estimates were computed using pooled variograms for each study region following methods described in Farmer (2016). Daily streamflow estimation was conducted by study region (hydrologic unit code level-2 regions as defined in Falcone, 2011) by building statistical models using 1,385 GAGES-II reference streamgages from mostly undisturbed watersheds as index gages (Russell and others, 2020). Estimates were then made at GAGES-II non-reference streamgages. Location information and basin characteristics for study gages were obtained from the GAGES-II dataset (Falcone, 2011). Observed daily streamflow data were retrieved from the National Water Information System (USGS, 2019). This data release contains 19 separate zip files; one for each study region. Each zip file contains an individual tab-delimited text file for each non-reference streamgage in the study region. A text file summarizing period of record information for each non-reference streamgage is provided (non-reference_gages_summary.csv). This data release also contains a text file (Model_info.csv) of regional regression equations for 27 flow quantiles that were developed in each study region in order to implement the QPPQ methods and a text file (BC_transformations.csv) describing transformations made to the GAGES-II derived basin characteristics prior to use in the regression equations. The five sets of streamflow estimates represent expected natural streamflow conditions with minimal disturbance by human activities, in other words, without the effects of regulation, diversion, land development, or other anthropogenic activities. The observed streamflow records at the non-reference streamgages were compared to the five simulated streamflow records. These performance metrics are provided at each gage for all five statistical methods (NonRef_PMs_byStation.csv) and as summaries by region (NonRef_PM_summaries_byRegion.csv). References cited: Falcone, J.A., 2011, GAGES-II: Geospatial Attributes of Gages for Evaluating Streamflow [digital spatial dataset]: U.S. Geological Survey Water Resources NSDI Node web page, https://water.usgs.gov/lookup/getspatial?gagesII_Sept2011. Farmer, W.H., Archfield, S.A., Over, T.M., Hay, L.E., LaFontaine, J.H., and Kiang, J.E., 2014, A comparison of methods to predict historical daily streamflow time series in the southeastern United States: U.S. Geological Survey Scientific Investigations Report 2014–5231, 34 p., http://dx.doi.org/10.3133/sir20145231. Farmer, W. H., 2016, Ordinary kriging as a tool to estimate historical daily streamflow records, Hydrology and Earth System Sciences, 20, 2721-2735, https://doi.org/10.5194/hess-20-2721-2016. Over, T.M., Farmer, W.H., and Russell, A.M., 2018, Refinement of a regression-based method for prediction of flow-duration curves of daily streamflow in the conterminous United States: U.S. Geological Survey Scientific Investigations Report 2018–5072, 34 p., https://doi.org/10.3133/sir20185072. Russell, A.M., Over, T.M., and Farmer, W.H., 2020, Cross-validation results for five statistical methods of daily streamflow estimation at 1,385 reference streamgages in the conterminous United States, Water Years

In 1991, the U.S. Geological Survey (USGS) began a study of more than 50 major river basins across the Nation as part of the National Water-Quality Assessment (NAWQA) project of the National Water-Quality Program. One of the major goals of the NAWQA project is to determine how water-quality conditions change over time. To support that goal, long-term consistent and comparable monitoring has been conducted on streams and rivers throughout the Nation. Outside of the NAWQA project, the USGS and other Federal, State, and local agencies also have collected long-term water-quality data to support their own assessments of changing water-quality conditions. Data from these multiple sources have been combined to support one of the most comprehensive assessment conducted to date of water-quality trends in the United States. Ultimately, these data will provide insight into how natural features and human activities have contributed to water-quality changes over time in Nation’s streams and rivers. This USGS data release contains all of the input and output files necessary to reproduce the results from the Weighted Regressions on Time, Discharge, and Season (WRTDS) models described in the associated U.S. Geological Survey Scientific Investigations Report. Data preparation for input to the models is also fully described in the above-mentioned report.

Nonstationary streamflow due to environmental and human-induced causes can affect water quality over time, yet these effects are poorly accounted for in water-quality trend models. This data release provides instream water-quality trends and estimates of two components of change, for sites across the Nation previously presented in Oelsner et al. (2017). We used previously calibrated Weighted Regressions on Time, Discharge, and Season (WRTDS) models published in De Cicco et al. (2017) to estimate instream water-quality trends and associated uncertainties with the generalized flow normalization procedure available in EGRET version 3.0 (Hirsch et al., 2018a) and EGRETci version 2.0 (Hirsch et al., 2018b). The procedure allows for nonstationarity in the flow regime, whereas previous versions of EGRET assumed streamflow stationarity. Water-quality trends of annual mean concentrations and loads (also referred to as fluxes) are provided as an annual series and the change between the start and end year for four trend periods (1972-2012, 1982-2012, 1992-2012, and 2002-2012). Information about the sites, including the collecting agency and associated streamflow gage, and information about site selection and the data screening process can be found in Oelsner et al. (2017). This data release includes results for 19 water-quality parameters including nutrients (ammonia, nitrate, filtered and unfiltered orthophosphate, total nitrogen, total phosphorus), major ions (calcium, chloride, magnesium, potassium, sodium, sulfate), salinity indicators (specific conductance, total dissolved solids), carbon (alkalinity, dissolved organic carbon, total organic carbon), and sediment (total suspended solids, suspended-sediment concentration) at over 1,200 sites. Note, the number of parameters with data varies by site with most sites having data for 1-4 parameters. Each water-quality trend was parsed into two components of change: (1) the streamflow trend component (QTC) and (2) the watershed management trend component (MTC). The QTC is an indicator of the amount of change in the water-quality trend attributed to changes in the streamflow regime, and the MTC is an indicator of the amount of change in the water-quality trend that may be attributed to human actions and changes in point and non-point sources in a watershed. Note, the MTC is referred to as the concentration-discharge trend component (CQTC) in the EGRET version 3.0 software. For our work, we chose to refer to this trend component as the MTC because it provides a more conceptual description (Murphy and Sprague, 2019). The trend results presented here expand upon the results in De Cicco et al. (2017) and Oelsner et al. (2017), which were analyzed using flow-normalization under the stationary streamflow assumption. The results presented in this data release are intended to complement these previously published results and support investigations into natural and human effects on water-quality trends across the United States. Data preparation information and WRTDS model specifications are described in Oelsner et al. (2017) and Murphy and Sprague (2019). This work was completed as part of the National Water-Quality Assessment (NAWQA) project of the National Water-Quality Program. De Cicco, L.A., Sprague, L.A., Murphy, J.C., Riskin, M.L., Falcone, J.A., Stets, E.G., Oelsner, G.P., and Johnson, H.M., 2017, Water-quality and streamflow datasets used in the Weighted Regressions on Time, Discharge, and Season (WRTDS) models to determine trends in the Nation’s rivers and streams, 1972-2012 (ver. 1.1 July 7, 2017): U.S. Geological Survey data release, https://doi.org/10.5066/F7KW5D4H. Hirsch, R., De Cicco, L., Watkins, D., Carr, L., and Murphy, J., 2018a, EGRET: Exploration and Graphics for RivEr Trends, version 3.0, https://CRAN.R-project.org/package=EGRET. Hirsch, R., De Cicco, L., and Murphy, J., 2018b, EGRETci: Exploration and Graphics for RivEr Trends (EGRET) Confidence Intervals, version 2.0.

In 1991, the U.S. Geological Survey (USGS) began a study of more than 50 major river basins across the Nation as part of the National Water-Quality Assessment (NAWQA) project of the National Water-Quality Program. One of the major goals of the NAWQA project is to determine how water-quality conditions change over time. To support that goal, long-term consistent and comparable monitoring has been conducted on streams and rivers throughout the Nation. Outside of the NAWQA project, the USGS and other Federal, State, and local agencies also have collected long-term water-quality data to support their own assessments of changing water-quality conditions. Data from these multiple sources have been combined to support one of the most comprehensive assessment conducted to date of water-quality trends in the United States. Ultimately, these data will provide insight into how natural features and human activities have contributed to water-quality changes over time in Nation’s streams and rivers. This USGS data release contains all of the input and output files necessary to reproduce the results from the Weighted Regressions on Time, Discharge, and Season (WRTDS) models described in the associated U.S. Geological Survey Scientific Investigations Report. Data preparation for input to the models is also fully described in the above-mentioned report.