This dataset provides numerical and categorical descriptions of 48 basin characteristics for 956 basins with observed streamflow information at U.S. Geological Survey (USGS) streamflow-gaging stations. Characteristics are indexed by National Hydrography Dataset (NHD) version 2 COMID (integer that uniquely identifies each feature in the NHD) and USGS station number for streamflow-gaging station. The variables represent mutable and immutable basin characteristics and are organized by characteristic type: physical (5), hydrologic (6), categorical (12), climate (6), landscape alteration (7), and land cover (12). Mutable characteristics such as climate, land cover, and landscape alteration variables are reported in decadal increments (for example, average percent forest for the decade 1950-1959, 1960-1969, etc). The majority of basin characteristics in this dataset were calculated using divergence-routing methods and are often referred to as “network-accumulated”. This method uses a modified routing database to navigate the NHDPlus reach network to aggregate (accumulate) the values derived from the reach catchment scale (Schwarz, G.E., and Wieczorek, M.E., 2018, Database of modified routing for NHDPlus version 2.1 flowlines: ENHDPlusV2_us: U.S. Geological Survey data release, https://doi.org/10.5066/P9PA63SM ). In four instances, values are also provided for the entire catchment above a site and area designated using the “CAT_” prefix.

This dataset provides numerical and categorical descriptions of 48 basin characteristics for 9,314 ungaged basins coinciding with 12-digit hydrologic unit code (HUC12) pour points that drain to the Gulf of Mexico. Characteristics are indexed by National Hydrography Dataset (NHD) version 2 COMID (integer that uniquely identifies each feature in the NHD) and HUC12 identifying number. The variables represent mutable and immutable basin characteristics and are organized by characteristic type: physical (5), hydrologic (6), categorical (12), climate (6), landscape alteration (7), and land cover (12). Mutable characteristics such as climate, land cover, and landscape alteration variables are reported in decadal increments (for example, average percent forest for the decade 1950-1959, 1960-1969, etc). The majority of basin characteristics in this dataset were calculated using divergence-routing methods and are often referred to as “network-accumulated”. This method uses a modified routing database to navigate the NHDPlus reach network to aggregate (accumulate) the values derived from the reach catchment scale (Schwarz, G.E., and Wieczorek, M.E., 2018, Database of modified routing for NHDPlus version 2.1 flowlines: ENHDPlusV2_us: U.S. Geological Survey data release, https://doi.org/10.5066/P9PA63SM ). In four instances, values are also provided for the entire catchment above a site and area designated using the “CAT_” prefix.

This geospatial dataset includes a one-point feature-class shapefile, one-polygon feature-class shapefile, and associated FGDC-compliant metadata to define 193 streamflow and 299 basin characteristics at 1,320 U.S. Geological Survey streamflow gaging stations. Sites included in the dataset either (1) drain to the Gulf of Mexico or (2) are adjacent to watersheds that flow to the Gulf of Mexico and are considered both physiographically similar and valuable for analysis. Drainage area to the sites varies from less than 1 to approximately 67,500 square miles. Data presented describe the streamflow regime (Rossman, 1990; Thompson and Archfield, 2014), climate (Daly and others, 2008), land use and land-use change (Sohl and others, 2014; Sohl and others, 2016), and anthropogenic features. Basins were identified following Hirsch and DiCicco (2015), and daily value streamflow data were retrieved from the USGS National Water Information System (U.S. Geological Survey, 2017). Daily value streamflow data were available beginning in 1892 through the 2016 water year (a 12-month period beginning October 1, for any given year through September 30 of the following year). All characteristics based on time series (streamflow, climate, land use for example) were summarized in terms of period of record and 10 water year increments (for example, 1930 – 1939). Data presented provide a numerical foundation supporting the: (1) development of statistical models of streamflow characteristics; (2) evaluation of spatial and temporal trends in streamflow characteristics; and (3) development of network optimization analysis. Basin characteristics will be used as independent variables to estimate streamflow characteristics (measures of the magnitude, duration, frequency, timing, and rate of change of the annual hydrograph) in a manner similar to Knight and others (2012). Daly, C., Halbleib, M., Smith, J.I., Gibson, W.P., Doggett, M.K., Taylor, G.H., Curtis, J., and Pasteris, P.P., 2008, Physiographically sensitive mapping of climatological temperature and precipitation across the conterminous United States: International Journal of Climatology, v. 28, no. 15, p. 2031–2064. Dunne, T., and Black, R., 1970. “An experimental investigation of runoff production in permeable soils.” Water Resour. Res., 6(2), 478–490 ESRI 2011. ArcGIS Desktop: Release 10.4.1 Redlands, CA: Environmental Systems Research Institute. Falcone, J.A., Carlisle, D.M., Wolock, D.M., and Meador, M.R., 2010b. GAGES: A stream gage database for evaluating natural and altered flow conditions in the conterminous United States, Ecology, 91 (2), p 621; Data Paper in Ecological Archives E091-045-D1; available online at: http://esapubs.org/Archive/ecol/E091/045/metadata.htm. Hamon, W.R., 1961. Estimating Potential Evaporation. Journal of the Hydraulics Division, Proceedings of American Society of Civil Engineers 87:107-120. Horton, Robert E. (1933) "The role of infiltration in the hydrologic cycle" Transactions of the American Geophysics Union, 14th Annual Meeting, pp. 446–460. Hirsch, R.M., and DiCicco, L.A., 2015, User guide to Exploration and Graphics for RivEr Trends (EGRET) and dataRetrieval: R packages for hydrologic data (version 2.0, February 2015):, accessed at https://pubs.usgs.gov/tm/04/a10/. Juracek, K.E., 1999, Estimation of potential runoff contributing areas in the Kansas-Lower Republican River Basin, Kansas: U.S. Geological Survey Water Resources Investigations Report 99-4089, 24 p Kjelstrom, L.C., 1998, Methods for estimating selected flow-duration and flood-frequency characteristics at ungaged sites in central Idaho: U.S. Geological Survey Water-Resources Investigations Report 94-4120, 10 p Knight, R.R., Gain, W.S., and Wolfe, W.J., 2012, Modelling ecological flow regime: an example from the Tennessee and Cumberland River basins: Ecohydrology, v. 5, no. 5, p. 613–627. NAWQA- U.S. Department of the Interior, U.S. Geological Survey. National Water-Quality Assessment (NAWQA) Program.

This dataset contains statistical descriptions of observed daily-mean streamflow for 956 sites in the southeast United States. For each site, statistical descriptions are provided according to decade for up to six decades, beginning in 1950 (1950-59 calendar years) and ending with 2000 (2000 - 2009 calendar years) with no more than 7 missing values per year in total (continuous or noncontinuous). There are 40 statistical descriptions including 28 flow-duration curve values, 8 L-moments, and 4 describing the minimum, maximum, median flow for days not equal to zero, and number of zero-flow days. Site information is provided by decade - the number of rows per site varies from 1 to 6 depending on the number of decades with observed record available. This information was used as the response variable(s) for statistical models for estimating the same characteristics at nearly 10,000 ungaged locations throughout the southeast United States. The dataset has been provided as a shapefile and a comma-delimited file. The comma-delimited file is an exact copy of the attribute table of the shapefile.

Censored and uncensored generalized additive models (GAMs) are developed from 955 U.S. Geological Survey streamflow-gaging stations (streamgages) to predict decadal statistics of streamflow. The streamgages are located on streams draining to the Gulf of Mexico. Decadal statistics include no-flow fractions and selected L-moments of nonzero streamflow for six decades (1950s—2000s). These statistics represent metrics of decadal flow-duration curves (dFDCs) derived from about 10 million daily mean streamflows. The L-moments include the mean, coefficient of L-variation, and the third through fifth L-moment ratios. The models are fit using watershed properties such as basin area and slope, decadal precipitation and temperature and decadal values of flood storage and urban development percentages. The GAMs then estimate decadal statistics for 8,988 prediction locations (stream reaches) coincident with outlets of level-12 hydrologic unit codes. Both the entire dataset (whole model) and leave-one-watershed-out model results are reported. No-flow fractions are censored data and Tobit extensions to GAMs are effective in estimation of ephemeral streamflow conditions. Uncensored GAMs conversely were used for estimation of the L-moments. The R language was used to pull and process the streamflow data, and the scripts can be found online at https://code.usgs.gov/water/restore/fdclmrpplo.

Censored and uncensored generalized additive models (GAMs) are developed from 955 U.S. Geological Survey streamflow-gaging stations (streamgages) to predict decadal statistics of streamflow. The streamgages are located on streams draining to the Gulf of Mexico. Decadal statistics include no-flow fractions and selected L-moments of nonzero streamflow for six decades (1950s—2000s). These statistics represent metrics of decadal flow-duration curves (dFDCs) derived from about 10 million daily mean streamflows. The L-moments include the mean, coefficient of L-variation, and the third through fifth L-moment ratios. The models are fit using watershed properties such as basin area and slope, decadal precipitation and temperature and decadal values of flood storage and urban development percentages. The GAMs then estimate decadal statistics for 8,988 prediction locations (stream reaches) coincident with outlets of level-12 hydrologic unit codes. Both the entire dataset (whole model) and leave-one-watershed-out model results are reported. No-flow fractions are censored data and Tobit extensions to GAMs are effective in estimation of ephemeral streamflow conditions. Uncensored GAMs conversely were used for estimation of the L-moments. The R language was used to pull and process the streamflow data, and the scripts can be found online at https://code.usgs.gov/water/restore/fdclmrpplo.

Censored and uncensored generalized additive models (GAMs) are developed from 955 U.S. Geological Survey streamflow-gaging stations (streamgages) to predict decadal statistics of streamflow. The streamgages are located on streams draining to the Gulf of Mexico. Decadal statistics include no-flow fractions and selected L-moments of nonzero streamflow for six decades (1950s—2000s). These statistics represent metrics of decadal flow-duration curves (dFDCs) derived from about 10 million daily mean streamflows. The L-moments include the mean, coefficient of L-variation, and the third through fifth L-moment ratios. The models are fit using watershed properties such as basin area and slope, decadal precipitation and temperature and decadal values of flood storage and urban development percentages. The GAMs then estimate decadal statistics for 8,988 prediction locations (stream reaches) coincident with outlets of level-12 hydrologic unit codes. Both the entire dataset (whole model) and leave-one-watershed-out model results are reported. No-flow fractions are censored data and Tobit extensions to GAMs are effective in estimation of ephemeral streamflow conditions. Uncensored GAMs conversely were used for estimation of the L-moments. The R language was used to pull and process the streamflow data, and the scripts can be found online at https://code.usgs.gov/water/restore/fdclmrpplo.

Two methods of calculating hydrologic alteration were applied to modeled daily streamflow data for 9,201 12-digit hydrologic unit code (HUC12) pour points draining to the Gulf of Mexico (Robinson and others, 2020). The first method is a new modified method of calculating ecosurplus and ecodeficit called hydro change. For this project, ecosurplus and ecodeficit have been combined to assess overall hydrologic regime change. The second method is the confidence interval hypothesis test (Kroll and others, 2015). The first method is a means of quantifying hydrologic alteration while the second is a hypothesis test to simply determine if statistically significant alteration has occurred. Both methods are employed to determine which is best at analyzing alteration of the hydrologic regime in the Gulf Coast Ecosystem Restoration Council (RESTORE) study area. Statistical analysis was done in RStudio (2020). The data release includes four attached files: (1) metadata .xml file, (2) csv with the p-values for each HUC12, (3) csv with results from the hydrologic change analysis, and (4) the shapefile of the pour point locations for the HUC12s used in the analyses.

This tabular data set represents base flow index (BFI) compiled for two spatial components of the NHDPlus version 2 data suite (NHDPlusv2) for the conterminous United States; 1) individual reach catchments and 2) reach catchments accumulated upstream through the river network. This dataset can be linked to the NHDPlus version 2 data suite by the unique identifier COMID. The source data is the "Base-flow index grid for the conterminous United States" produced by the United States Geological Survey (Wolock, 2003). Units are percent. Base flow is the component of streamflow that can be attributed to ground-water discharge into streams. The BFI is the ratio of base flow to total flow, expressed as a percentage. Reach catchment information characterizes data at the local scale. Reach catchments accumulated upstream through the river network characterizes cumulative upstream conditions. Network-accumulated values are computed using two methods, 1) divergence-routed and 2) total cumulative drainage area. Both approaches use a modified routing database to navigate the NHDPlus reach network to aggregate (accumulate) the metrics derived from the reach catchment scale. (Schwarz and Wieczorek, 2018).

This tabular data set represents estimated the 30 year (1951-1980) mean annual natural groundwater recharge compiled for two spatial components of the NHDPlus version 2 data suite (NHDPlusv2) for the conterminous United States; 1) individual reach catchments and 2) reach catchments accumulated upstream through the river network. This dataset can be linked to the NHDPlus version 2 data suite by the unique identifier COMID. The source data for estimated mean annual natural groundwater recharge was produced by the United States Geological Survey (Wolock, 2003). Units are millimeters per year. Reach catchment information characterizes data at the local scale. Reach catchments accumulated upstream through the river network characterizes cumulative upstream conditions. Network-accumulated values are computed using two methods, 1) divergence-routed and 2) total cumulative drainage area. Both approaches use a modified routing database to navigate the NHDPlus reach network to aggregate (accumulate) the metrics derived from the reach catchment scale. (Schwarz and Wieczorek, 2018).