Geology topics

Steven K. Sando

Publications and source records attributed to Steven K. Sando.

At least 19 recordsLinked to original sources

USGS2026

Streamflow, water quality, and energy development characteristics and water-quality trends for selected sites in the Green River and Muddy Creek drainage basins, southwestern Wyoming, water years 2006–16

Energy resource development in southwestern Wyoming has been rapidly expanding since the early 2000s to include innovative technologies for extracting oil and gas and enhancing renewable energy. In response to concerns about the potential negative effects of energy development on natural resources, the Wyoming Landscape Conservation Initiative was created to implement a long-term, science-based program of assessing, conserving, and enhancing fish and wildlife habitats while facilitating responsible energy development through local collaborations and partnerships. The U.S. Geological Survey, in cooperation with the Bureau of Land Management, began monitoring streamflow and water quality as part of this initiative in 2006 at four U.S. Geological Survey streamgages within the Green River and Muddy Creek drainage basins to characterize the effects of energy development on streamflow and water quality in these areas. To quantify changes in water quality metrics, trend analyses for the 2006–16 water years were calculated and related to drainage basin characteristics including energy development (wells per square mile). Effects of energy development on water quality varied by watershed and the variability was driven primarily by the runoff efficiency and the geochemistry of the drainage basin. The New Fork River water chemistry was a calcium-bicarbonate type with low ionic strength and showed no statistically significant trends in water quality despite increased energy development. The Green River experienced moderate energy development, with significant upward trends in magnesium and alkalinity, indicating low sensitivity to development. Muddy Creek near Dad had minimal energy development and complex hydrogeochemical characteristics, showing mixed trends in water quality constituents. Muddy Creek near Baggs experienced moderate energy development, with significant upward trends in several water quality indicators indicating a response to the land disturbance associated with development. Overall, geochemistry and runoff efficiency influenced the vulnerability of water quality to energy development activities, highlighting the importance of drainage basin physical characteristics on vulnerability to development.

Wyoming↗

USGS2025

Hydroclimatic and land-use factors affecting peak streamflow in Illinois, Iowa, Michigan, Minnesota, Missouri, Montana, North Dakota, South Dakota, and Wisconsin

Flood-frequency analysis provides the basis for flood risk estimates used by water-resource managers in land-use planning, and it informs the design of essential infrastructure such as bridges and culverts. Federal guidelines for flood-frequency analysis do not offer guidance on addressing changing climate and land-use conditions when estimating floods. However, failing to consider climatic and land-use changes that cause abrupt or gradual changes in flood regimes can result in a poor representation of the true flood risk. In response to concerns about changing flood regimes, the U.S. Geological Survey, in cooperation with nine State agencies (Illinois Department of Transportation, Iowa Department of Transportation, Michigan Department of Transportation, Minnesota Department of Transportation, Missouri Department of Transportation, Montana Department of Natural Resources and Conservation, North Dakota Department of Water Resources, South Dakota Department of Transportation, and Wisconsin Department of Transportation) began a study to examine variability and change in hydrology and climate and the effects of urbanization and tile drainage on flooding. The analyses of patterns and changes in hydrology and climate were reported in a multichapter Scientific Investigations Report, the findings of which are summarized in this U.S. Geological Survey Circular. Additional analyses documenting changes in seasonality of flooding and the effects of urbanization and tile drainage were completed and published as separate studies and are also summarized in this Circular. These studies provide extensive exploratory analysis of peak streamflow, daily streamflow, and climate data, setting the stage for advancements in flood-frequency analysis.

Illinois, Iowa, Michigan, Minnesota, Missouri, Mon↗

USGS2025

Peak streamflow trends in Montana and northern Wyoming and their relation to changes in climate, water years 1921–2020

Frequency analysis on annual peak streamflow (hereinafter, peak flow) is essential to water-resources management applications, including critical structure design (for example, bridges and culverts) and floodplain mapping. Nonstationarity is a statistical property of a peak-flow series such that the distributional properties (the mean, variance, or skew) change either gradually (monotonic trend) or abruptly (shift, step change or change point) through time. Not incorporating or accounting for observed nonstationarity into peak-flow frequency analysis might result in a poor representation of the true probability of large floods and thus misrepresent the actual flood risks to life and property. This report summarizes how hydroclimatic variability might affect the temporal and spatial distributions of peak-flow data in the State of Montana (and northern Wyoming) and is part of a larger study to document peak-flow nonstationarity and hydroclimatic changes across a nine-State region consisting of Illinois, Iowa, Michigan, Minnesota, Missouri, Montana, North Dakota, South Dakota, and Wisconsin. 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 flows. This study was completed in cooperation with the Montana Department of Natural Resources and Conservation. The purpose of this report is to characterize temporal and spatial patterns of nonstationarity in peak flows and hydroclimatology in Montana and northern Wyoming. In this evaluation, peak-flow, daily streamflow, and model-simulated gridded climatic data were examined for monotonic trends, change points, and other statistical properties that might indicate changing climatic and environmental conditions. This report includes background information on the study area, the history of U.S. Geological Survey peak-flow data collection and frequency analysis in Montana, and the review of research relating to hydroclimatic variability and change in Montana. This study might help provide a framework for addressing potential nonstationarity issues in peak-flow frequency updates that commonly are completed by the U.S. Geological Survey in cooperation with other agencies throughout the Nation. The analytical structure of this study includes analyses of monotonic trends and change points in numerous hydroclimatic variables in assigned 30-, 50-, 75-, and 100-year analysis periods. For Montana and part of Wyoming, the 30-, 50-, 75, and 100-year analyses included 157, 70, 48, and 12 streamgages, respectively. For those streamgages, nonstationarities were analyzed in the following variables: (1) climatic variables, including annual and seasonal (winter, spring, summer, and fall) temperature and precipitation; (2) daily streamflow variables, including the annual center of volume duration, annual center of volume median, and peaks over threshold with a mean of four events per year; and (3) annual peak-flow variables, including peak-flow timing and magnitude. A likelihood approach was used to express statistical confidence and assign the nonstationarity results as likely upward or downward (highest statistical confidence), somewhat likely upward or downward (less statistical confidence), or about as likely as not (little statistical confidence; hereinafter, neutral). For the nonstationarity analyses of the climatic, daily streamflow, and peak-flow variables, the results are presented in detail and discussed with respect to statewide patterns and geographic variability. For each of the 30-, 50-, and 75-year analyses, peak-flow change-point and monotonic trend analyses were compiled for streamgages classified with likely downward or likely upward trends. For those streamgages, the associated basin characteristics and nonstationarity results for peak-flow timing, daily streamflow, and climatic variables were investigated and statistically compared to discern associations among other variables that might contribute to the peak-flow nonstationarity results. The 50- and 75-year peak-flow nonstationarities identified in this study are mostly downward, in association with mostly upward temperature and potential evapotranspiration:precipitation monotonic trends. For the 50-, 75-, and 100-year analyses, the peak-flow change points are predominantly downward and are concentrated in the 1970s and 1980s, which indicates general consistency among the longer trend periods. These findings are in association with substantial research documenting globally rising temperature and atmospheric greenhouse gas concentrations that might be largely attributed to anthropogenic activities. Anthropogenic effects might represent long-term (on the order of several decades to more than a century) climate changes that might happen within highly variable natural climate fluctuations. Several paleo studies in the north-central United States have indicated that hydroclimatic extremes (that is, low- and high-streamflow conditions) before European settlement have been outside of extremes since the 1900s. Depending on the interactions of anthropogenic effects and natural climate variability, extreme high-streamflow conditions might occur in the future, even in the presence of long-term downward peak-flow trends.

Montana, Wyoming↗

USGS2025

Changes in streamflow seasonality associated with hydroclimatic variability in the north-central United States among three discrete temporal periods, 1946–2020

Study region North-central United States Study focus This study uses circular statistics to characterize the seasonal properties of annual maximum (AMS) and peaks-over-threshold (POT) streamflow time series for 841 and 623 selected U.S. Geological Survey (USGS) streamgages, respectively, without regulation or substantial diversion among common 75-, 50-, and 30-year trend periods through water year 2020 (the period from October 1, 2019, through September 30, 2020). A subset of AMS time series with detected change points (abrupt changes) in the median and (or) scale are analyzed on either side of the change point to evaluate changes in their circular statistics. New hydrologic insights for the region In the 50-year trend period, five regions share common mean flood timing in the AMS and POT partial duration series. Changes from asymmetric distributions to reflective symmetric distributions are detected particularly among the 50- and 30-year trend periods in the northernmost States of Minnesota, North Dakota, and Wisconsin. For the subset of streamgages with abrupt change points in the AMS, regional patterns of changes in seasonality are detected between the period of records before and after the change point. These findings can inform decisions related to the AMS used for flood frequency and potential mixed population analyses and flood control operations that may be affected by changes in when seasonal events occur, how long seasonal events last, and the long-term variability in the intensity and frequency of seasonal events.

Illinois, Iowa, Michigan, Minnesota, Missouri, Mon↗

USGS2024

Peak streamflow trends in South Dakota and their relation to changes in climate, water years 1921–2020

Peak-flow (flood) frequency analysis is essential to water-resources management applications, including the design of critical infrastructure such as bridges and culverts, and floodplain mapping. Federal guidelines for performing peak-flow 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, which documents the guidelines for determining annual peak streamflow frequency, is that, for basins without major hydrologic alterations (for example, regulation, diversion, and urbanization), statistical properties of the distribution of annual peak streamflows are stationary; that is, the mean, variance, and skew are constant through time. Nonstationarity is a statistical property of a peak-flow series such that the long-term (on the order of decades) distributional properties change one or more times either gradually or abruptly through time. Individual nonstationarities may be attributed to one source such as flow regulation, land-use change, or climate but are often the result of a combination of sources, making detection and attribution of nonstationarities challenging. In response to a growing concern regarding nonstationarity in peak streamflows in the region, 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, assessed the potential nonstationarity in peak streamflows in the north-central United States. This chapter characterizes the effects of natural hydroclimatic shifts and potential climate change on annual peak streamflows in the State of South Dakota. Annual peak and daily streamflow as well as model-simulated gridded climatic data were examined for temporal monotonic trends, change points, and other statistical properties indicative of changing climatic and environmental conditions. Changes in annual peak and daily flows were evaluated among 13, 35, and 81 qualifying U.S. Geological Survey streamgages for the 75-, 50-, and 30-year trend periods through water year 2020 (the period from October 1, 2019, to September 30, 2020) in South Dakota, respectively. No qualifying streamgages were in the 100-year trend period in the State. Statistical tests for autocorrelation (independent and identically distributed assumption), monotonic trends, and change points in the median and scale are analyzed to evaluate potential stationarity violations (nonstationarity) for performing at-site peak-flow flood-frequency analysis. The trends are reported using a likelihood approach as an alternative to simply reporting significant trends with an arbitrary p -value cutoff point. A distinct east-west spatial pattern of likely upward and downward monotonic trends and change points, respectively, was detected in 75- and 50-year trend periods, but an inconsistent spatial pattern was detected in the 30-year trend period. Additionally, change points in the median annual peak streamflows were detected in the late 1970s and early 1980s in the western part of the State, but in the east, the change point was more commonly detected in 1992–93. A similar east-west spatial pattern of likely upward and downward trends was detected in the annual peak-flow timing, the day of the year of the annal peak streamflow. In the western part of the State, the annual peak streamflows are arriving earlier, but in the east, the annual peak streamflows are arriving later. A peaks-over-threshold (POT) analysis where, on average, there are two events per year (POT2) and four events per year (POT4) was also used to evaluate changes in the frequency (count) of daily streamflows exceeding the threshold. Similar to detected changes in the annual peak streamflow, an east-west likely upward or downward change corresponding to an increase or decrease, respectively, in the frequency of daily streamflow greater than a POT2 and POT4 threshold was detected. A monthly water-balance model was used to evaluate hydroclimatic variation in annual and seasonal precipitation, snowfall, potential evapotranspiration, and soil moisture storage for all qualifying streamgages in the 75-, 50-, and 30-year trend periods. Detected trends in the annual hydroclimatic metrics for the 75- and 50-year trend periods indicate a spatially consistent statewide increase in precipitation, decrease in snowfall, increase in potential evapotranspiration, and increase in soil moisture storage. Furthermore, detected trends in seasonal precipitation in the 75- and 50-year trend periods highlight a pronounced change in precipitation in winter and later into the summer season, especially in the 50-year trend period in the eastern part of the State. Statewide increases in seasonal soil moisture storage were also detected, highlighting year-round increasing flood magnitudes, particularly in the eastern part of the State. Based on the results of these stationarity tests for the qualifying streamgages in South Dakota among the 75-, 50-, and 30-year trend periods, consistent temporal and spatial patterns of nonstationarity were detected among the 75- and 50-year trend periods. Furthermore, when nonstationarity is detected in daily streamflow, increased streamflow and volume (increasing frequency in POT), as well as potentially bridge scour, may have implications on culvert and highway design in the eastern part of South Dakota. Thus, when performing at-site peak-flow flood-frequency analyses in South Dakota, potential nonstationarities and alternative approaches are important considerations.

South Dakota↗

USGS2024

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↗

USGS2023

Application of surrogate technology to predict real-time metallic-contaminant concentrations and loads in the Clark Fork near Grant-Kohrs Ranch National Historic Site, Montana, water years 2019–20

Grant-Kohrs Ranch National Historic Site (GRKO) in southwestern Montana commemorates the frontier cattle era and its formative role in shaping the culture and history of the Western United States. The ranch was designated a national historic landmark in 1960 and a unit of the National Park Service (NPS) by Congress in 1972. The GRKO is unique because of its proximity to large-scale extraction, milling, and smelting of gold, silver, copper, and lead ore from the 1860s to the 1980s in the Butte mining district. During this time, mining and milling wastes were discarded in the upper Clark Fork Basin, resulting in the deposition of large amounts of waste materials (tailings) enriched with metallic contaminants (including cadmium, copper, iron, lead, manganese, zinc, and the metalloid trace element arsenic) in soils and in nearby streams and floodplains. Denuded vegetation and fish kills attributed to large concentrations of heavy metals caused the U.S. Environmental Protection Agency to designate a 120-mile section of the Clark Fork River (hereafter referred to as the “Clark Fork”), including GRKO, to be included on the National Priority List for Superfund cleanup in 1989. In 2018, with oversight from the Montana Department of Environmental Quality, the NPS began remediation of 2.6 miles of the Clark Fork as it flows through GRKO property. In 2019, the U.S. Geological Survey (USGS), in collaboration with the NPS, conducted a study using time-series data from backscatter signals from fixed-point turbidity and acoustic sensors with the intent to provide a high-resolution monitoring tool to estimate metallic-contaminant concentrations (MCCs) and loads during NPS remediation of the Clark Fork. Two monitoring sites at USGS streamgages on the Clark Fork on either side of GRKO property were instrumented with turbidity and acoustic sensors and surrogate relations were developed among time-series data and MCCs. The application of high-resolution surrogate data was used to infer contaminant source and fate and evaluate MCC values relative to aquatic-life standards. Using high-resolution surrogate data, it was determined that during spring runoff and storm-related runoff events, MCCs peaked at their highest values at streamflows markedly lower and prior to peak streamflow. Because MCCs peaked prior to streamflow peaks, it could be inferred that the source of MCCs originated from channel bed sediments in close spatial proximity to the monitoring site or from nearby streambanks and floodplains. High-resolution surrogate data revealed that copper concentrations in the Clark Fork exceeded chronic aquatic-life standards 90 percent of the time when streamflow exceeded 200 cubic feet per second (ft 3 /s) and exceeded acute aquatic-life standards 85 percent of the time when streamflow exceeded 260 ft 3 /s. These data helped support NPS management goals for evaluating variation in water quality during remediation of GRKO property, evaluating MCC values relative to aquatic-life standards, and quantifying benefits from Superfund remediation activities.

Montana↗

USGS2022

Potential effects of energy development on environmental resources of the Williston Basin in Montana, North Dakota, and South Dakota—Water resources

The Williston Basin has been a leading oil and gas producing area for more than 50 years. While oil production initially peaked within the Williston Basin in the mid-1980s, production rapidly increased in the mid-2000s, largely because of improved horizontal (directional) drilling and hydraulic fracturing methods. In 2012, energy development associated with the Bakken Formation was identified as a priority requiring collaboration toward improved timeliness of issuing permits for new wells combined with reasonable measures to maintain environmental quality. Shortly thereafter, the Bakken Federal Executive Group was created to address common challenges associated with energy development. The Bakken Federal Executive Group partner agencies identified a gap in current understanding of the cumulative environmental challenges attributed to energy development throughout the area, resulting in an effort to aggregate scientific data and identify additional research and information needs related to natural resources within areas of energy development in the Williston Basin. As part of this effort, water resources in the area (including groundwater; streams and rivers; and lakes, reservoirs, and wetlands) were characterized and described in terms of physical occurrence, flow characteristics, recharge, water quality, and water use. Similarly, waters produced during energy-development activities also were characterized even though these waters are not considered usable resources within the area. Groundwater resources were characterized by the major hydrogeologic units, or aquifers, identifying the units that supply most groundwater used for domestic, stock, agricultural, and industrial purposes. The groundwater characterization included other deeper hydrogeologic units in the Williston Basin that may be a useable source of water with treatment, have utility as a reservoir for reinjection of produced waters, or be a source of minerals and energy resources. A generalized groundwater budget and flow system identifying the sources of recharge (stream infiltration, precipitation, and movement [leakage] from other aquifers) and the general groundwater flow direction is included for each of the major hydrogeologic units. Rivers and streams within the Williston Basin with 10 or more years of continuous streamflow data were identified. For a subset of these sites, streamflow characteristics, including the monthly and annual mean flow, were generated to identify seasonal and interannual changes in streamflow and thus provide information on the drivers and reliability of streamflow at the seasonal or multiyear scale. Daily streamflow and annual extreme flows (peak and low flow) also were estimated for the subset of sites. The daily streamflow and annual extreme flow values provide information on short-term or extreme events that are relevant to infrastructure design and evaluating spills, leaks, or accidental discharges of water or petroleum products. Surface-water features (lakes, ponds, and wetlands) were classified using the Cowardin system and identified on the National Wetlands Inventory maps generated by the U.S. Fish and Wildlife Service. The spatial distribution of the surface-water features was analyzed by State, county, and specifically in comparison to the Prairie Pothole Region. The proximity of the surface-water features to energy development infrastructure (specifically oil or gas well pads) was evaluated. It was determined that, although oil or gas wells are often near a surface-water feature, most surface-water features do not have wells nearby, with the exception of wells in the Prairie Pothole Region. Water-quality data were aggregated from two data sources: (1) the Water-Quality Portal, sponsored by the U.S. Geological Survey (USGS), U.S. Environmental Protection Agency (EPA), and National Water Quality Monitoring Council; and (2) a data compilation completed as part of the USGS National Water-Quality Assessment project. The Water-Quality Portal integrates publicly available water-quality data from databases maintained by the USGS, EPA, and U.S. Department of Agriculture, including water-quality data from Tribal, State, and local databases. Water-quality data for 15 commonly measured water-quality constituents were aggregated for groundwater, rivers and streams, and lakes and reservoirs. For each aggregated dataset (groundwater, rivers and streams, and lakes and reservoirs), analyses of the water-quality data included summary statistics, maps of spatial distribution of constituent values, boxplots of constituent values by timeframe or hydrogeologic unit, spatial comparisons of site locations and constituent values to petroleum well density, and comparisons of the constituent values measured to EPA drinking-water standards/guidelines. Produced water includes all fluids brought to the surface along with the targeted hydrocarbons as part of the oil and gas exploration and extraction processes. These fluids may include formation water (waters that co-exist with rock/oil/gas), hydraulic fracturing fluids, and other combinations of water and chemicals used during oil and gas well drilling, development, treatments, recompletions, and workovers. Produced water datasets were aggregated from two sources: the USGS National Produced Waters Geochemical database (ver. 2.1) and a series of projects focused specifically on sampling produced water in the Williston Basin from 2010 to 2014. The National Produced Waters Geochemical database was useful for a general understanding of produced-water chemistry. Produced waters are characterized by extreme salinity and contain elevated concentrations of other constituents (including arsenic, barium, cadmium, lead, zinc, radium-226/radium-228, and ammonium) that could negatively affect water and aquatic resources if released. Produced waters also have a generally unique chemical (isotopic) signature that may be useful in tracking water from different geologic units; for example, the oxygen/deuterium and strontium ratio values measured in brine waters from the Bakken Formation are distinct from brines collected from other geologic units in the Williston Basin. Water-use information related to energy production in the area also was aggregated and summarized. The summary of water use is not limited to oil and gas production but includes water used to produce all types of energy resources in the Williston Basin, including coal/lignite, thermoelectric power, oil and gas, hydropower, biomass and biofuels, wind, geothermal, and solar. Each State has its own methods for regulating and reporting water usage within its jurisdiction. These methods can introduce problems when examining water use from sources, such as the Missouri River or Fox Hills aquifer, that are shared across political boundaries. Without the one-to-one match for usage types and amounts used from a water source, it is difficult to develop a comprehensive water budget for the water source being evaluated. A large amount of freshwater is required to prepare a well for oil and gas well production; in some cases, 3 to 7 million gallons of water are needed per well. The EPA estimates that hydraulic fracturing in the Williston Basin uses between 70 to 140 billion gallons per year. Water also is used for myriad other purposes related to ancillary oil and gas extraction. In addition to water used for immediate energy development, the expanded human workforce migrating into the area and other support staff who have moved into the area during the development also use water. Research and information needs were identified that could be relevant in the evaluation of the effects of energy development on water resources. Information needs related to the evaluation of groundwater resources include the following: improved potentiometric-surface maps for glacial units; availability of a uniform stream network digital geographic coverage that spans the international boundary with Canada; enhanced surface-water use information with regards to the gain and loss of streamflow to shallow groundwater, which would increase understanding groundwater and surface-water interactions; and expanded geophysical assessments. Gaps in the availability of streamflow data include the lack of information on ice-jam flooding despite potential for effects to infrastructure (pipelines, roads, and facilities) and an understanding of the cumulative effects of largely undocumented stock and diversion dams. Although this study resulted in the aggregation of a large quantity of water-quality data, the availability of consistently collected, systematically processed and reported data over large parts of the Williston Basin is sparse. Few samples have been analyzed for constituents that may indicate the effect of energy development on water resources. Constituents that could be considered include boron, chloride, bromide, iodine, fluoride, manganese, lithium, radium, strontium isotopes, volatile organic compounds, and isotopes of inorganic ions (such as hydrogen and carbon). Collaboration between Tribal, Federal, State, and local entities to identify a common study design, common monitoring constituents, and consistent sampling locations would generate datasets with broad utility and would likely result in overall cost savings for monitoring over time. Similarly, there is a need for standardized sample collection, processing, laboratory analytical methods, and the collection of ancillary data for produced waters sampling. Additional characterization of the range of chemical, microbial, and isotopic compositions and quantities of “end-member” produced waters, and the collection of time-series datasets to document the changes in produced waters during and after well development also were needs identified during this study. Water-use estimates would be improved through the implementation of comprehensive studies of water use from groundwater and surface-water sources using consistent methodologies across the Williston Basin. The submission of chemical and water data related to hydraulic fracturing collected by the oil and gas industry would add to the quantity of available data. Consistent implementation of regulations and monitoring controls across political boundaries (State, county, and international) would further improve the consistency of data available for the estimates of water use.

Montana, North Dakota, South Dakota↗

USGS2021

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↗

USGS2018

Methods for peak-flow frequency analysis and reporting for streamgages in or near Montana based on data through water year 2015

This report documents the methods for peak-flow frequency (hereinafter “frequency”) analysis and reporting for streamgages in and near Montana following implementation of the Bulletin 17C guidelines. The methods are used to provide estimates of peak-flow quantiles for 50-, 42.9-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities for selected streamgages operated by the U.S. Geological Survey Wyoming-Montana Water Science Center (WY–MT WSC). These annual exceedance probabilities correspond to 2-, 2.33-, 5-, 10-, 25-, 50-, 100-, 200-, and 500-year recurrence intervals, respectively. Standard procedures specific to the WY–MT WSC for implementing the Bulletin 17C guidelines include (1) the use of the Expected Moments Algorithm analysis for fitting the log-Pearson Type III distribution, incorporating historical information where applicable; (2) the use of weighted skew coefficients (based on weighting at-site station skew coefficients with generalized skew coefficients from the Bulletin 17B national skew map); and (3) the use of the Multiple Grubbs-Beck Test for identifying potentially influential low flows. For some streamgages, the peak-flow records are not well represented by the standard procedures and require user-specified adjustments informed by hydrologic judgement. The specific characteristics of peak-flow records addressed by the informed-user adjustments include (1) regulated peak-flow records, (2) atypical upper-tail peak-flow records, and (3) atypical lower-tail peak-flow records. In all cases, the informed-user adjustments use the Expected Moments Algorithm fit of the log-Pearson Type III distribution using the at-site station skew coefficient, a manual potentially influential low flow threshold, or both. Appropriate methods can be applied to at-site frequency estimates to provide improved representation of long-term hydroclimatic conditions. The methods for improving at-site frequency estimates by weighting with regional regression equations and by Maintenance of Variance Extension Type III record extension are described. Frequency analyses were conducted for 99 example streamgages to indicate various aspects of the frequency- analysis methods described in this report. The frequency analyses and results for the example streamgages are presented in a separate data release associated with this report consisting of tables and graphical plots that are structured to include information concerning the interpretive decisions involved in the frequency analyses. Further, the separate data release includes the input files to the PeakFQ program, version 7.1, including the peak-flow data file and the analysis specification file that were used in the peak-flow frequency analyses. Peak-flow frequencies are also reported in separate data releases for selected streamgages in the Beaverhead River and Clark Fork Basins and also for selected streamgages in the Ruby, Jefferson, and Madison River Basins.

Montana↗

USGS2016

Water-quality trends and constituent-transport analysis for selected sampling sites in the Milltown Reservoir/Clark Fork River Superfund Site in the upper Clark Fork Basin, Montana, water years 1996–2015

During the extended history of mining in the upper Clark Fork Basin in Montana, large amounts of waste materials enriched with metallic contaminants (cadmium, copper, lead, and zinc) and the metalloid trace element arsenic were generated from mining operations near Butte and milling and smelting operations near Anaconda. Extensive deposition of mining wastes in the Silver Bow Creek and Clark Fork channels and flood plains had substantial effects on water quality. Federal Superfund remediation activities in the upper Clark Fork Basin began in 1983 and have included substantial remediation near Butte and removal of the former Milltown Dam near Missoula. To aid in evaluating the effects of remediation activities on water quality, the U.S. Geological Survey began collecting streamflow and water-quality data in the upper Clark Fork Basin in the 1980s. Trend analysis was done on specific conductance, selected trace elements (arsenic, copper, and zinc), and suspended sediment for seven sampling sites in the Milltown Reservoir/Clark Fork River Superfund Site for water years 1996–2015. The most upstream site included in trend analysis is Silver Bow Creek at Warm Springs, Montana (sampling site 8), and the most downstream site is Clark Fork above Missoula, Montana (sampling site 22), which is just downstream from the former Milltown Dam. Water year is the 12-month period from October 1 through September 30 and is designated by the year in which it ends. Trend analysis was done by using a joint time-series model for concentration and streamflow. To provide temporal resolution of changes in water quality, trend analysis was conducted for four sequential 5-year periods: period 1 (water years 1996–2000), period 2 (water years 2001–5), period 3 (water years 2006–10), and period 4 (water years 2011–15). Because of the substantial effect of the intentional breach of Milltown Dam on March 28, 2008, period 3 was subdivided into period 3A (October 1, 2005–March 27, 2008) and period 3B (March 28, 2008–September 30, 2010) for the Clark Fork above Missoula (sampling site 22). Trend results were considered statistically significant when the statistical probability level was less than 0.01. In conjunction with the trend analysis, estimated normalized constituent loads (hereinafter referred to as “loads”) were calculated and presented within the framework of a constituent-transport analysis to assess the temporal trends in flow-adjusted concentrations (FACs) in the context of sources and transport. The transport analysis allows assessment of temporal changes in relative contributions from upstream source areas to loads transported past each reach outflow. Trend results indicate that FACs of unfiltered-recoverable copper decreased at the sampling sites from the start of period 1 through the end of period 4; the decreases ranged from large for one sampling site (Silver Bow Creek at Warm Springs [sampling site 8]) to moderate for two sampling sites (Clark Fork near Galen, Montana [sampling site 11] and Clark Fork above Missoula [sampling site 22]) to small for four sampling sites (Clark Fork at Deer Lodge, Montana [sampling site 14], Clark Fork at Goldcreek, Montana [sampling site 16], Clark Fork near Drummond, Montana [sampling site 18], and Clark Fork at Turah Bridge near Bonner, Montana [sampling site 20]). For period 4 (water years 2011–15), the most notable changes indicated for the Milltown Reservoir/Clark Fork River Superfund Site were statistically significant decreases in FACs and loads of unfiltered-recoverable copper for sampling sites 8 and 22. The period 4 changes in FACs of unfiltered-recoverable copper for all other sampling sites were not statistically significant. Trend results indicate that FACs of unfiltered-recoverable arsenic decreased at the sampling sites from period 1 through period 4 (water years 1996–2015); the decreases ranged from minor (sampling sites 8–20) to small (sampling site 22). For period 4 (water years 2011–15), the most notable changes indicated for the Milltown Reservoir/Clark Fork River Superfund Site were statistically significant decreases in FACs and loads of unfiltered-recoverable arsenic for sampling site 8 and near statistically significant decreases for sampling site 22. The period 4 changes in FACs of unfiltered-recoverable arsenic for all other sampling sites were not statistically significant. Trend results indicate that FACs of suspended sediment decreased at the sampling sites from period 1 through period 4 (water years 1996–2015); the decreases ranged from moderate (sampling site 8) to small (sampling sites 11–22). For period 4 (water years 2011–15), the changes in FACs of suspended sediment were not statistically significant for any sampling sites. The reach of the Clark Fork from Galen to Deer Lodge is a large source of metallic contaminants and suspended sediment, which strongly affects downstream transport of those constituents. Mobilization of copper and suspended sediment from flood-plain tailings and the streambed of the Clark Fork and its tributaries within the reach results in a contribution of those constituents that is proportionally much larger than the contribution of streamflow from within the reach. Within the reach from Galen to Deer Lodge, unfiltered-recoverable copper loads increased by a factor of about 4 and suspended-sediment loads increased by a factor of about 5, whereas streamflow increased by a factor of slightly less than 2. For period 4 (water years 2011–15), unfiltered-recoverable copper and suspended-sediment loads sourced from within the reach accounted for about 41 and 14 percent, respectively, of the loads at Clark Fork above Missoula (sampling site 22), whereas streamflow sourced from within the reach accounted for about 4 percent of the streamflow at sampling site 22. During water years 1996–2015, decreases in FACs and loads of unfiltered-recoverable copper and suspended sediment for the reach generally were proportionally smaller than for most other reaches. Unfiltered-recoverable copper loads sourced within the reaches of the Clark Fork between Deer Lodge and Turah Bridge near Bonner (just upstream from the former Milltown Dam) were proportionally smaller than contributions of streamflow sourced from within the reaches; these reaches contributed proportionally much less to copper loading in the Clark Fork than the reach between Galen and Deer Lodge. Although substantial decreases in FACs and loads of unfiltered-recoverable copper and suspended sediment were indicated for Silver Bow Creek at Warm Springs (sampling site 8), those substantial decreases were not translated to downstream reaches between Deer Lodge and Turah Bridge near Bonner. The effect of the reach of the Clark Fork from Galen to Deer Lodge as a large source of copper and suspended sediment, in combination with little temporal change in those constituents for the reach, contributes to this pattern. With the removal of the former Milltown Dam in 2008, substantial amounts of contaminated sediments that remained in the Clark Fork channel and flood plain in reach 9 (downstream from Turah Bridge near Bonner) became more available for mobilization and transport than before the dam removal. After the removal of the former Milltown Dam, the Clark Fork above Missoula (sampling site 22) had statistically significant decreases in FACs of unfiltered-recoverable copper in period 3B (March 28, 2008, through water year 2010) that continued in period 4 (water years 2011–15). Also, decreases in FACs of unfiltered-recoverable arsenic and suspended sediment were indicated for period 4 at this site. The decrease in FACs of unfiltered-recoverable copper for sampling site 22 during period 4 was proportionally much larger than the decrease for the Clark Fork at Turah Bridge near Bonner (sampling site 20). Net mobilization of unfiltered-recoverable copper and arsenic from sources within reach 9 are smaller for period 4 than for period 1 when the former Milltown Dam was in place, providing evidence that contaminant source materials have been substantially reduced in reach 9.

Montana↗

USGS2016

Montana StreamStats—A method for retrieving basin and streamflow characteristics in Montana: Chapter A in Montana StreamStats

The U.S. Geological Survey (USGS) provides streamflow characteristics and other related information needed by water-resource managers to protect people and property from floods, plan and manage water-resource activities, and protect water quality. Streamflow characteristics provided by the USGS, such as peak-flow and low-flow frequencies for streamflow-gaging stations, are frequently used by engineers, flood forecasters, land managers, biologists, and others to guide their everyday decisions. In addition to providing streamflow characteristics at streamflow-gaging stations, the USGS also develops regional regression equations and drainage area-adjustment methods for estimating streamflow characteristics at locations on ungaged streams. Regional regression equations can be complex and often require users to determine several basin characteristics, which are physical and climatic characteristics of the stream and its drainage basin. Obtaining these basin characteristics for streamflow-gaging stations and ungaged sites traditionally has been time consuming and subjective, and led to inconsistent results. StreamStats is a Web-based geographic information system application that was created by the USGS to provide users with access to an assortment of analytical tools that are useful for water-resource planning and management. StreamStats allows users to easily obtain streamflow and basin characteristics for USGS streamflow-gaging stations and user-selected locations on ungaged streams. The USGS, in cooperation with Montana Department of Transportation, Montana Department of Environmental Quality, and Montana Department of Natural Resources and Conservation, completed a study to develop a StreamStats application for Montana, compute streamflow characteristics at streamflow-gaging stations, and develop regional regression equations to estimate streamflow characteristics at ungaged sites. Chapter A of this Scientific Investigations Report describes the Montana StreamStats application and the datasets, streamflow-gaging stations, streamflow characteristics, and regression equations (as described fully in Chapters B through G of this report) that are used for development of the StreamStats application for Montana.

Montana↗

USGS2016

Temporal trends and stationarity in annual peak flow and peak-flow timing for selected long-term streamflow-gaging stations in or near Montana through water year 2011: Chapter B in Montana StreamStats

A large-scale study by the U.S. Geological Survey, in cooperation with the Montana Department of Transportation and the Montana Department of Natural Resources and Conservation, was done to investigate general patterns in peak-flow temporal trends and stationarity through water year 2011 for 24 long-term streamflow-gaging stations (hereinafter referred to as gaging stations) in Montana. Hereinafter, all years refer to water years; a water year is the 12-month period from October 1 through September 30 and is designated by the year in which it ends. The primary focus of the study was to identify general patterns in peak-flow temporal trends and stationarity that are relevant to application of peak-flow frequency analyses within a statewide gaging-station network. Temporal trends were analyzed for two hydrologic variables: annual peak flow and peak-flow timing. Annual peak flow is the maximum instantaneous discharge, in cubic feet per second, recorded each year a gaging station was operated. Peak-flow timing is the day of the annual peak flow (hereinafter referred to as day of peak), recorded each year a gaging station was operated. Study results provide evidence that annual peak flow for most of the long-term gaging stations can be reasonably considered as stationary for application of peak-flow frequency analyses within a statewide gaging station network. Upward trends in annual peak flow during 1930–76 generally were stronger than downward trends during 1967–2011 for most long-term gaging stations. Statistical distributions of annual peak flow generally were similar among three summary time periods (1930–78, 1979–2011, and the entire period of record). However, for two low-elevation gaging stations in eastern Montana (Poplar River at international boundary [gaging station 06178000] and Powder River at Moorhead, Montana [gaging station 06324500]), substantial downward trends in annual peak flow during 1967–2011 were of similar or stronger magnitude than the upward trends during 1930–76, and the annual-peak-flow medians for 1979–2011 were substantially lower than the medians for the entire period of record. For peak-flow timing for most long-term gaging stations, differences in trends between 1930–76 and 1967–2011 are variable and not particularly strong. Statistical distributions generally are similar among the summary time periods. However, for two high-elevation gaging stations on a Missouri River headwater tributary (Gallatin River near Gallatin Gateway, Montana [gaging station 06043500] and Gallatin River at Logan, Montana [gaging station 06052500]) and for five high-elevation gaging stations in the Yellowstone River Basin (Yellowstone River at Corwin Springs, Montana [gaging station 06191500], Yellowstone River near Livingston, Montana [gaging station 06192500], Clarks Fork Yellowstone River near Belfry, Montana [gaging station 06207500], Clarks Fork Yellowstone River at Edgar, Montana [gaging station 06208500], and Yellowstone River at Billings, Montana [gaging station 06214500]) downward trends in peak-flow timing during 1967–2011 generally were stronger than upward trends during 1930–1976, and day-of-peak medians for 1979–2011 were considerably less than medians for 1930–78. The downward trends in peak-flow timing for 1967–2011 indicate that the timing of annual peak flows changed from later in the year to earlier in the year. For the seven high-elevation gaging stations, the mean change during 1967–2011 was about 13 days (range of 8 to 22 days). For most of the high-elevation gaging stations in the Missouri River headwaters, Yellowstone River Basin, and Columbia River Basin, there was general correspondence between trend patterns for annual peak flow and trend patterns for peak-flow timing; that is, during periods when there were upward trends in annual peak flow, there generally also were upward trends in peak-flow timing. Conversely, during periods when there were downward trends in annual peak flow, there generally also were downward trends in peak-flow timing. The two low-elevation gaging stations in eastern Montana (Poplar River at international boundary [gaging station 06178000] and Powder River at Moorhead, Montana [gaging station 06324500]) had considerable changes in annual-peakflow characteristics after the mid-1970s, which might provide evidence of potential nonstationarity in the peak-flow records. The two low-elevation gaging stations that have potential nonstationarity are located in drainage basins that are strongly affected by agricultural activities that potentially affect the hydrologic regimes. Primary agricultural activities that might alter natural hydrologic conditions include construction of small impoundments (primarily for stock-watering purposes) and irrigation diversions. Temporal variability in these activities might contribute to the potential nonstationarity issues. Changes in climatic characteristics after the mid-1970s also possibly contribute to the potential nonstationarity issues. Lack of considerable indication of potential nonstationarity in annual peak flow for the other long-term gaging stations in this study might indicate that climatic changes have been more pronounced with respect to effects on peak flows in low elevation areas in eastern Montana than in areas represented by the other long-term gaging stations. Another possibility is that climatic changes after the mid-1970s are exacerbated in low-elevation areas where small-impoundment development and potential effects of irrigation diversions might be more extensive.

Montana↗

USGS2016

Peak-flow frequency analyses and results based on data through water year 2011 for selected streamflow-gaging stations in or near Montana: Chapter C in Montana StreamStats

Chapter C of this Scientific Investigations Report documents results from a study by the U.S. Geological Survey, in cooperation with the Montana Department of Transportation and the Montana Department of Natural Resources, to provide an update of statewide peak-flow frequency analyses and results for Montana. The purpose of this report chapter is to present peak-flow frequency analyses and results for 725 streamflow-gaging stations in or near Montana based on data through water year 2011. The 725 streamflow-gaging stations included in this study represent nearly all streamflowgaging stations in Montana (plus some from adjacent states or Canadian Provinces) that have at least 10 years of peak-flow records through water year 2011. For 29 of the 725 streamflow-gaging stations, peak-flow frequency analyses and results are reported for both unregulated and regulated conditions. Thus, peak-flow frequency analyses and results are reported for a total of 754 analyses. Estimates of peak-flow magnitudes for 66.7-, 50-, 42.9-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities are reported. These annual exceedance probabilities correspond to 1.5-, 2-, 2.33-, 5-, 10-, 25-, 50-, 100-, 200-, and 500-year recurrence intervals.

Montana↗

USGS2016

Adjusted peak-flow frequency estimates for selected streamflow-gaging stations in or near Montana based on data through water year 2011: Chapter D in Montana StreamStats

The climatic conditions of the specific time period during which peak-flow data were collected at a given streamflow-gaging station (hereinafter referred to as gaging station) can substantially affect how well the peak-flow frequency (hereinafter referred to as frequency) results represent long-term hydrologic conditions. Differences in the timing of the periods of record can result in substantial inconsistencies in frequency estimates for hydrologically similar gaging stations. Potential for inconsistency increases with decreasing peak-flow record length. The representativeness of the frequency estimates for a short-term gaging station can be adjusted by various methods including weighting the at-site results in association with frequency estimates from regional regression equations (RREs) by using the Weighted Independent Estimates (WIE) program. Also, for gaging stations that cannot be adjusted by using the WIE program because of regulation or drainage areas too large for application of RREs, frequency estimates might be improved by using record extension procedures, including a mixed-station analysis using the maintenance of variance type I (MOVE.1) procedure. The U.S. Geological Survey, in cooperation with the Montana Department of Transportation and the Montana Department of Natural Resources and Conservation, completed a study to provide adjusted frequency estimates for selected gaging stations through water year 2011. The purpose of Chapter D of this Scientific Investigations Report is to present adjusted frequency estimates for 504 selected streamflow-gaging stations in or near Montana based on data through water year 2011. Estimates of peak-flow magnitudes for the 66.7-, 50-, 42.9-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities are reported. These annual exceedance probabilities correspond to the 1.5-, 2-, 2.33-, 5-, 10-, 25-, 50-, 100-, 200-, and 500-year recurrence intervals, respectively. The at-site frequency estimates were adjusted by weighting with frequency estimates from RREs using the WIE program for 438 selected gaging stations in Montana. These 438 selected gaging stations (1) had periods of record less than or equal to 40 years, (2) represented unregulated or minor regulation conditions, and (3) had drainage areas less than about 2,750 square miles. The weighted-average frequency estimates obtained by weighting with RREs generally are considered to provide improved frequency estimates. In some cases, there are substantial differences among the at-site frequency estimates, the regression-equation frequency estimates, and the weighted-average frequency estimates. In these cases, thoughtful consideration should be applied when selecting the appropriate frequency estimate. Some factors that might be considered when selecting the appropriate frequency estimate include (1) whether the specific gaging station has peak-flow characteristics that distinguish it from most other gaging stations used in developing the RREs for the hydrologic region; and (2) the length of the peak-flow record and the general climatic characteristics during the period when the peak-flow data were collected. For critical structure-design applications, a conservative approach would be to select the higher of the at-site frequency estimate and the weighted-average frequency estimate. The mixed-station MOVE.1 procedure generally was applied in cases where three or more gaging stations were located on the same large river and some of the gaging stations could not be adjusted using the weighted-average method because of regulation or drainage areas too large for application of RREs. The mixed-station MOVE.1 procedure was applied to 66 selected gaging stations on 19 large rivers. The general approach for using mixed-station record extension procedures to adjust at-site frequencies involved (1) determining appropriate base periods for the gaging stations on the large rivers, (2) synthesizing peak-flow data for the gaging stations with incomplete peak-flow records during the base periods by using the mixed-station MOVE.1 procedure, and (3) conducting frequency analysis on the combined recorded and synthesized peak-flow data for each gaging station. Frequency estimates for the combined recorded and synthesized datasets for 66 gaging stations with incomplete peak-flow records during the base periods are presented. The uncertainties in the mixed-station record extension results are difficult to directly quantify; thus, it is important to understand the intended use of the estimated frequencies based on analysis of the combined recorded and synthesized datasets. The estimated frequencies are considered general estimates of frequency relations among gaging stations on the same stream channel that might be expected if the gaging stations had been gaged during the same long-term base period. However, because the mixed-station record extension procedures involve secondary statistical analysis with accompanying errors, the uncertainty of the frequency estimates is larger than would be obtained by collecting systematic records for the same number of years in the base period.

Montana↗

USGS2016

Methods for estimating peak-flow frequencies at ungaged sites in Montana based on data through water year 2011: Chapter F in Montana StreamStats

The U.S. Geological Survey (USGS), in cooperation with the Montana Department of Natural Resources and Conservation, completed a study to update methods for estimating peak-flow frequencies at ungaged sites in Montana based on peak-flow data at streamflow-gaging stations through water year 2011. The methods allow estimation of peak-flow frequencies (that is, peak-flow magnitudes, in cubic feet per second, associated with annual exceedance probabilities of 66.7, 50, 42.9, 20, 10, 4, 2, 1, 0.5, and 0.2 percent) at ungaged sites. The annual exceedance probabilities correspond to 1.5-, 2-, 2.33-, 5-, 10-, 25-, 50-, 100-, 200-, and 500-year recurrence intervals, respectively. Regional regression analysis is a primary focus of Chapter F of this Scientific Investigations Report, and regression equations for estimating peak-flow frequencies at ungaged sites in eight hydrologic regions in Montana are presented. The regression equations are based on analysis of peak-flow frequencies and basin characteristics at 537 streamflow-gaging stations in or near Montana and were developed using generalized least squares regression or weighted least squares regression. All of the data used in calculating basin characteristics that were included as explanatory variables in the regression equations were developed for and are available through the USGS StreamStats application (http://water.usgs.gov/osw/streamstats/) for Montana. StreamStats is a Web-based geographic information system application that was created by the USGS to provide users with access to an assortment of analytical tools that are useful for water-resource planning and management. The primary purpose of the Montana StreamStats application is to provide estimates of basin characteristics and streamflow characteristics for user-selected ungaged sites on Montana streams. The regional regression equations presented in this report chapter can be conveniently solved using the Montana StreamStats application. Selected results from this study were compared with results of previous studies. For most hydrologic regions, the regression equations reported for this study had lower mean standard errors of prediction (in percent) than the previously reported regression equations for Montana. The equations presented for this study are considered to be an improvement on the previously reported equations primarily because this study (1) included 13 more years of peak-flow data; (2) included 35 more streamflow-gaging stations than previous studies; (3) used a detailed geographic information system (GIS)-based definition of the regulation status of streamflow-gaging stations, which allowed better determination of the unregulated peak-flow records that are appropriate for use in the regional regression analysis; (4) included advancements in GIS and remote-sensing technologies, which allowed more convenient calculation of basin characteristics and investigation of many more candidate basin characteristics; and (5) included advancements in computational and analytical methods, which allowed more thorough and consistent data analysis. This report chapter also presents other methods for estimating peak-flow frequencies at ungaged sites. Two methods for estimating peak-flow frequencies at ungaged sites located on the same streams as streamflow-gaging stations are described. Additionally, envelope curves relating maximum recorded annual peak flows to contributing drainage area for each of the eight hydrologic regions in Montana are presented and compared to a national envelope curve. In addition to providing general information on characteristics of large peak flows, the regional envelope curves can be used to assess the reasonableness of peak-flow frequency estimates determined using the regression equations.

Montana↗

USGS2016

Methods for estimating streamflow characteristics at ungaged sites in western Montana based on data through water year 2009: Chapter G in Montana StreamStats

The U.S. Geological Survey, in cooperation with the Montana Department of Environmental Quality and the Montana Department of Natural Resources and Conservation, developed regional regression equations based on basin and streamflow characteristics for streamflow-gaging stations through water year 2009 that can be used to estimate streamflow characteristics for ungaged sites in western Montana. The regression equations allow estimation of low-flow frequencies; mean annual and mean monthly streamflows; and the 20-, 50-, and 80-percent durations for annual and monthly duration streamflows for ungaged sites in western Montana that are unaffected by regulation. Regression equations were developed for 54 streamflow characteristics in each of the 4 hydrologic regions generally located in western Montana. Exploratory data analyses determined that of the 40 basin characteristics used as explanatory variables, only 3 were significant in the final equations. These three basin characteristics are contributing drainage area, mean annual precipitation of the drainage basin, and percentage of basin with slopes greater than 50 percent. Low-flow frequency data and daily mean streamflow data from 152 streamflow-gaging stations were used to develop the regression equations. These 152 streamflow-gaging stations were considered to be unregulated (or only the unregulated period of record for regulated streamflow-gaging stations were used). For this study, a streamflow-gaging station is considered to be unregulated if cumulative drainage area upstream from all upstream dams is less than 20 percent of the drainage basin for the streamflow-gaging station, and if no large diversion canals are upstream from the streamflow-gaging station. In addition to the regional regression equations, two methods for estimating streamflow characteristics based on drainage areas of a nearby streamflow-gaging station are provided. Also, spreadsheet tools are provided for calculating streamflow characteristics based on regression equations or the drainage-area ratio method. All of the data used to calculate basin characteristics were derived from publicly available data sources and are available through the U.S. Geological Survey Streamstats program ( http://water.usgs.gov/osw/streamstats/ ) for Montana. The primary purpose of the Montana StreamStats application is to provide estimates of basin characteristics and streamflow characteristics for user-selected ungaged sites on Montana streams. The regional regression equations presented in this report have been loaded to the Montana StreamStats application and can be used to derive streamflow characteristics for ungaged sites.

Montana↗

USGS2015

Water-quality trends for selected sites in the Boulder River and Tenmile Creek watersheds, Montana, based on data collected during water years 1997-2013

In the Boulder River and Tenmile Creek watersheds in southwestern Montana, there was intensive mining during a 40-year period after the discovery of gold in the early 1860s. Potential effects from the historic mining activities include acid-mine drainage and elevated concentrations of potentially toxic trace elements from mining remnants such as waste rock and tailing piles. In support of remediation efforts, water-quality monitoring by the U.S. Geological Survey began in 1997 in the Boulder River and Tenmile Creek watersheds and has continued to present (2014). The U.S. Geological Survey, in cooperation with the U.S. Forest Service, investigated temporal trends in water quality at 13 sites, including 2 adit (or mine entrance) sites and 11 stream sites. The primary purpose of this report is to present results of trend analysis of specific conductance, selected trace-elements (cadmium, copper, lead, zinc, and arsenic), and suspended sediment for the 13 sites. Trend results for most stream sites in the Boulder River watershed for water years 2000–13 (water year is the 12-month period from October 1 through September 30 and is designated by the year in which it ends) indicate decreasing trends in flow-adjusted specific conductance, in flow-adjusted concentrations (FACs) for most filtered and unfiltered-recoverable trace elements, and in suspended sediment. Overall, magnitudes of the decreasing trends in FACs of metallic contaminants are largest for Bullion Mine tributary at mouth (site 3), Jack Creek at mouth (site 4), and Cataract Creek at Basin (site 8). For sites 3, 4, and 8, magnitudes of decreasing trends generally ranged from about -5 to -10 percent per year. Notably, the watersheds upstream from sites 3, 4, and 8 have been targeted by substantial remediation activities. Consideration of trend patterns among all stream sites in the Boulder River watershed provides strong evidence that remediation activities are the primary cause of decreasing trends in metallic contaminants. Trend results for sites in the Tenmile Creek watershed generally are more variable and difficult to interpret than for sites in the Boulder River watershed. Trend results for Tenmile Creek above City Diversion (site 11) and Minnehaha Creek near Rimini (site 12) for water years 2000–13 indicate decreasing trends in FACs of cadmium, copper, and zinc. The magnitudes of the decreasing trends in FACs of copper generally are moderate and statistically significant for sites 11 and 12. The magnitudes of the decreasing trends in FACs of cadmium and zinc for site 11 are minor to small and not statistically significant; however, the magnitudes for site 12 are moderate and statistically significant. In general, patterns in FACs for Tenmile Creek near Rimini (site 13) are not well represented by fitted trends within the short data collection period, which might indicate that the trend-analysis structure of the study is not appropriate for describing trends in FACs for site 13. The large decreasing trend in FACs of suspended sediment is the strongest indication of change in water quality during the short period of record for site 13; however, this trend is not statistically significant.

Montana↗