Paleohydrologic records provide a valuable perspective on the variability of streamflow and hydroclimate that is critical for water resource planning and placing present day and future conditions into a long-term context. Until now, key insights gained from streamflow reconstructions in the other river basins across the Western U.S. been lacking in the Upper Missouri River Basin due to a lack of extended streamflow records. Here we utilize a new database of naturalized streamflow records for the Upper Missouri and an expanded network of tree-ring records from the region to reconstruct streamflow at 31 gaging locations across the major Mountain Headwaters of the United States’ largest river basin. The database also includes an Upper Missouri Basin Basin composite record of streamflow that is not specific to any streamgage location, but rather summarizes streamflow variability across all the major gaging locations in the Upper Missouri River. The reconstructions explain an average of 68% of the variability in the observed streamflow records and extend records of streamflow to C.E. 886 on average. The network of streamflow reconstructions presented here fills a major geographical void in paleohydrologic understanding and provides important data resources to water managers balancing increasing water demands for hydropower, irrigation, navigation, and ecological resources with increasing flood risk in the basin.

Tree-ring reconstructions of water-year (Oct 1 through Sep 30th) flow for 31 gaging sites in Missouri River basin, with complete data for 1685 through 1977 (n = 293 water years). The complete 105 tree-ring reconstructions of streamflow used in the Nature Scientific Reports paper were obtained from various sources; 74 flow reconstructions were obtained from the web resource, TreeFlow (http://www.treeflow.info/), and an additional 31 flow reconstructions (provided here) for the Missouri River basin were obtained from the U.S. Geological Survey in Bozeman, Montana.

Tree-ring reconstructions of water-year (Oct 1 through Sep 30th) flow for 31 gaging sites in Missouri River basin, with complete data for 1685 through 1977 (n = 293 water years). The complete 105 tree-ring reconstructions of streamflow used in the Nature Scientific Reports paper were obtained from various sources; 74 flow reconstructions were obtained from the web resource, TreeFlow (http://www.treeflow.info/), and an additional 31 flow reconstructions (provided here) for the Missouri River basin were obtained from the U.S. Geological Survey in Bozeman, Montana.

Paleohydrologic records provide a valuable perspective on the variability of streamflow and hydroclimate that is critical for water resource planning and placing present day and future conditions into a long-term context. Until now, key insights gained from streamflow reconstructions in the other river basins across the Western U.S. have been lacking in the Milk and St. Mary River Basin. Here we utilize a new database of naturalized streamflow records for the Milk and St. Mary Rivers and an expanded network of tree-ring records from the region to reconstruct streamflow at eight gaging locations located in the mountains, foothills, and plains reaches of the basins. The network of streamflow reconstructions presented here were generated for use by the Bureau of Reclamation and Montana Department of Natural Resources and Conservation in the Basin Study update for the Milk and St. Mary Rivers, and provides important data resources to water managers balancing increasing water demands for hydropower, irrigation, and ecological resources with increasing drought and flood risk in the basin.

Paleohydrologic records provide a valuable perspective on the variability of streamflow and hydroclimate that is critical for water resource planning and placing present day and future conditions into a long-term context. Until now, key insights gained from streamflow reconstructions in the other river basins across the Western U.S. have been lacking in the Milk and St. Mary River Basin. Here we utilize a new database of naturalized streamflow records for the Milk and St. Mary Rivers and an expanded network of tree-ring records from the region to reconstruct streamflow at eight gaging locations located in the mountains, foothills, and plains reaches of the basins. The network of streamflow reconstructions presented here were generated for use by the Bureau of Reclamation and Montana Department of Natural Resources and Conservation in the Basin Study update for the Milk and St. Mary Rivers, and provides important data resources to water managers balancing increasing water demands for hydropower, irrigation, and ecological resources with increasing drought and flood risk in the basin.

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 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 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 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.