USGS ScienceSearch

Geology topics

Aldo V. Vecchia

Publications and source records attributed to Aldo V. Vecchia.

At least 19 recordsLinked to original sources

Water-quality trends for selected sites and constituents in the international Red River of the North Basin, Minnesota and North Dakota, United States, and Manitoba, Canada, 1970–2017

A comprehensive study to evaluate water-quality trends, while considering natural hydroclimatic variability, in the Red River of the North Basin and assess water-quality conditions for the Red River of the North crossing the international boundary near Emerson, Manitoba, Canada (the binational site), was completed by the U.S. Geological Survey in cooperation with the International Joint Commission, North Dakota Department of Environmental Quality, and Minnesota Pollution Control Agency and in collaboration with Manitoba Sustainable Development and Environment and Climate Change Canada. The international Red River of the North Basin encompasses 3 U.S. States (South Dakota, North Dakota, and Minnesota) and 1 Canadian Province (Manitoba). Water quality in the Red River of the North Basin is of concern for both Federal governments and State and Provincial governments. Water-quality objectives have been previously established for selected dissolved ions and recently (2019) proposed for selected nutrients for the binational site. In the current (2020) study, water-quality data from State, Provincial, and Federal agencies in the United States and Canada for sites in the Red River of the North Basin from 1970 to 2017 were compiled and used for trend analysis. Trend analysis using a water-quality dataset from multiple agencies that collect water-quality data for various objectives presented multiple challenges. The trend-analysis approach was able to accommodate differences in water-quality data caused by field-collection and laboratory-analytical method differences, disparities in sampling frequencies, and spatial and temporal gaps in data. Most of these challenges were overcome by the statistical tool, R–QWTREND, which identifies trends in concentration unrelated to variability in streamflow. The integrated basin approach used in the current study, combined with comparing current data trends with historical trends, provided valuable insights into understanding how water quality is changing spatially (34 sites analyzed for a recent period, 2000–15) and temporally (5 sites analyzed for a 45-year historical period, 1970–2015) within the Red River of the North Basin. One of the most consistent spatial and temporal changes observed in the current study was increasing concentrations of sulfate among tributary and main-stem sites since 2000. For some sites, increases were detected starting as early as 1985. Total dissolved solids and chloride concentrations had spatial and temporal patterns like sulfate. Although R–QWTREND removes the variability in constituent concentration caused by natural streamflow variability, all variability in sulfate caused by hydroclimatic variability may not be captured because of changes in hydrologic pathways and changes in the contributions of sulfate from various natural sources. Nutrient concentrations demonstrated less consistent spatial and temporal changes than sulfate, and changes in nutrient concentrations were assumed to be more closely tied to human-induced rather than natural changes. Nitrate-plus-nitrite concentrations were mostly increasing in the upper Red River of the North subbasin, and for nitrate plus nitrite and total nitrogen, the Sheyenne River subbasin had consistent decreasing concentrations. Since 2000, total phosphorus has decreased in the upper Red River of the North subbasin, but total phosphorus concentration has increased for sites in the lower Red River of the North subbasin, and for some main-stem sites, concentrations have been increasing since 1985. Unlike sulfate, the pattern in historical trends for total phosphorus for the main-stem sites differed from tributary sites, indicating that human-induced changes affected tributaries and main-stem sites differently. The more detailed evaluation of flow-averaged water-quality conditions for the binational site provided an understanding of how loads have changed over time and what proportion of the year and season concentrations are expected to exceed water-quality objectives. In a basin with highly variable streamflow like the Red River of the North, the trend in flow-averaged load (assuming streamflow conditions are the same year after year) provided a robust measure of change over time. Increasing concentrations of sulfate, chloride, total dissolved solids, and total phosphorus since 1985 for the binational site resulted in longer periods of exceedance of water-quality objectives per year occurring over time. For total nitrogen, decreasing concentrations resulted in shorter periods of exceedance per year during 1980 to 2015, but concentrations were still expected to exceed the water-quality objective about half the year. Periods of when exceedances were likely to occur during the year were affected by the source and transport mechanisms of the constituent. Trend results from this effort identified how water quality has changed across the basin, and further investigation would help to identify causes for the trends observed here. Information from the current study provides a basis for future trend attribution studies, evaluation of water-quality objectives, and development of comprehensive strategies for reducing nutrients to desired targets and establishes a baseline for tracking future progress in the Red River of the North Basin.

Minnesota, North Dakota, South Dakota, Manitoba

Time-series model, statistical methods, and software documentation for R–QWTREND—An R package for analyzing trends in stream-water quality

As part of a U.S. Geological Survey water-quality study started in 2018, in cooperation with the International Joint Commission, North Dakota Department of Environmental Quality, and Minnesota Pollution Control Agency, a publicly available software package called R–QWTREND was developed for analyzing trends in stream-water quality. The R–QWTREND package is a collection of functions written in R, an open source language and a general environment for statistical computing and graphics. The package uses a parametric time-series model to express logarithmically transformed concentration in terms of flow-related variability, trend, and serially correlated model errors. Flow-related variability captures natural variability in concentration on the basis of concurrent and antecedent streamflow. The trend identifies systematic changes in concentration in terms of potential step trends, piecewise monotonic trends, or user-specified trends. Maximum likelihood estimation is used to estimate model parameters and determine the best-fit trend model. This report describes the time-series model and statistical methodology behind R–QWTREND and provides formal documentation for installing and using the package.

Open-File Report

Landscape drivers of dynamic change in water quality of US rivers

Water security is a top concern for social well-being and dramatic changes in the availability of freshwater have occurred as a result of human uses and landscape management. Elevated nutrient loading and perturbations to major ion composition have resulted from human activities and have degraded freshwater resources. This study addresses the emerging nature of stream water quality in the 21st century through analysis of concentrations and trends in a wide variety of constituents in streams and rivers of the U.S. Concentrations of fifteen separate water quality parameters including nutrients, major ions, sediment, and specific conductance were analyzed over the period 1982-2012 and a targeted trend analysis was performed from 1992-2012. Although environmental policy is geared toward addressing the long-standing problem of nutrient overenrichment, these efforts have had uneven success, with decreasing nutrient concentrations at urbanized sites and little to no change at agricultural sites. However, freshwaters are being salinized rapidly in all human-dominated land use types. Increasing salinity negatively affects biodiversity, mobilizes sediment-bound contaminants, and increases lead contamination of drinking water but the effects are poorly quantified. Therefore, while efforts to control nutrients are ongoing, rapid salinity increases are ushering in a new set of poorly-defined issues.

Environmental Science & Technology

Water-balance modeling of selected lakes for evaluating viability as long-term fisheries in Kidder, Logan, and Stutsman Counties, North Dakota

Water levels in lakes and wetlands in the central North Dakota Missouri Coteau region that were either dry or only sporadically held water since before the 1930s have been rising since the early 1990s in response to an extended wet period. The lakes have remained full since the mid-1990s, which has provided benefits to migratory waterfowl, fisheries, and wildlife. A small shift in climate conditions, either to drier or wetter conditions, can have a large effect on the lake levels of these water bodies. The North Dakota Game and Fish Department identified five lakes as candidates for sustaining long-term fisheries. The lakes are in Kidder, Stutsman, and Logan Counties, and some lakes might receive inflow from mostly freshwater aquifers, such as the Central Dakota and Streeter aquifers, and were mostly dry during the early 1990s. After about 1995, the lakes had filled up and were deep enough to sustain populations of game fish such as walleye, perch, and northern pike. Before investing in development of permanent fisheries and associated infrastructure, such as campgrounds and boat ramps, fisheries biologists needed to know if the lake levels are likely to remain high in coming decades. The U.S. Geological Survey, in cooperation with the North Dakota Game and Fish Department, developed a water-balance model to determine the effects of precipitation, evapotranspiration, and groundwater interaction on lake volumes. The model was developed using climate input data and lake volumes for the calibration period 1992 through 2016, during which historical lake volumes could be estimated using land surface elevation data and Landsat images. Long-term (1940–2018) climate input data were used with the water-balance model to reconstruct historical lake volumes prior to the calibration period, and block-bootstrapping was used to simulate potential future climate input data and lake volumes for 2017 through 2067. The simulated future lake volumes were used to estimate the likelihood of annual lake volumes remaining consistent, increasing, or decreasing through the year 2067. Of the five lakes, Sibley Lake was the most likely to sustain a long-term fishery for a period longer than 50 years. The simulated lake volumes for Alkaline Lake, Big Mallard Marsh, and Remmick Lake indicated the lakes have a 50-percent chance to fall below 75 percent of their 2016 volume by about 2030, 2067, and 2025, respectively. Simulation results for Marvin Miller Lake were substantially different compared to the other four lakes and indicated the lake has a 50-percent chance to fall below 75 percent of its 2016 volume prior to 2025.

North Dakota

Stochastic model for simulating Souris River Basin regulated streamflow upstream from Minot, North Dakota

The Souris River Basin is a 24,000 square-mile basin in the Provinces of Saskatchewan and Manitoba in Canada, and the State of North Dakota in the United States. Above-average snowpack during the winter of 2010–11, along with record-setting rains in May and June of 2011, led to record flooding that caused extensive damage to Minot, North Dakota, and numerous smaller communities in Saskatchewan, Manitoba, and North Dakota. As a result, the International Souris River Board created the Souris River Flood Task Force to evaluate potential reservoir operation changes and flood control measures to manage future floods and droughts. Part of this evaluation involved identifying a need for a stochastic streamflow model to estimate the likelihood of future flooding or drought. A stochastic natural (unregulated) streamflow simulation model described in a previous report was built upon in this report to include the effects of regulation of four reservoirs (Rafferty, Alameda, and Boundary Reservoirs and Lake Darling) and their operation guidelines. First, a regulated reservoir storage/streamflow routing model was developed and calibrated from when all four reservoirs were in operation until the end of the reconstructed natural streamflow dataset provided by the U.S. Army Corps of Engineers (1992–2011). The regulated reservoir storage/streamflow routing model then was combined with the stochastic natural (unregulated) streamflow model to provide a stochastic regulated streamflow simulation model for the Souris River Basin upstream from Minot, North Dakota. The stochastic regulated streamflow simulation model was used to estimate regulated flood frequency curves, which are useful for feasibility and design of critical structures such as levees or bridges. Three potential future climatic conditions were considered in this analysis: condition A (wet equilibrium), representing wet (similar to 1970–2017) climatic conditions; condition B (transition), representing transition from wet to dry (similar to 1912–69) climatic conditions; and condition C (dry equilibrium), representing dry climatic conditions. Comparison of the estimated flood frequency curves for regulated flow among the three climatic conditions indicated large differences in flood magnitudes for the more extreme (1-percent or less) annual exceedance probabilities. The estimated 0.2-percent annual exceedance probability flood magnitude for the Souris River upstream from Minot, N. Dak., was 29,300 cubic feet per second for condition A (wet equilibrium), compared to 14,800 cubic feet per second for condition C (dry equilibrium). For comparison, the recorded peak streamflow for 2011 for the Souris River upstream from Minot, N. Dak., was 26,900 cubic feet per second. Although it is not possible to predict how long the current (1970–2017) wet climatic conditions may persist, flood risk for at least the next 25 years, or until about 2040, may be represented best by climatic condition A.

North Dakota

Interactive tool to estimate groundwater elevations in central and eastern North Dakota

This report describes an interactive tool (NDakGWtool) in which a statistical model is developed using locally weighted regression to estimate monthly mean groundwater elevations for a specified latitude and longitude, referred to as the “user-specified location.” For each user-specified location, seven models are developed for each month from April through October. Localized, high spatial-resolution maps of estimated monthly mean groundwater surface elevations are produced from the models. The tool was evaluated for glacial drift aquifers of the 32-county study area in central and eastern North Dakota. Although groundwater elevations from 1960 to 2017 were available to develop the tool, groundwater elevations from 1995 to 2015 were used for model testing and development of the model domain. There are 413 grid cells of 0.1-degree latitude by 0.1-degree longitude size in the model domain, and the tool produces maps of estimated monthly mean groundwater surface elevations for the cell containing the user-specified location. Additionally, the NDakGWtool produces maps of estimated groundwater depth below land surface and ArcGIS files of estimated groundwater surface elevations and groundwater depth below land surface. The tool is composed of four main components: data input, statistical model, output, and user-interactive process.

North Dakota

Model methodology for estimating pesticide concentration extremes based on sparse monitoring data

This report describes a new methodology for using sparse (weekly or less frequent observations) and potentially highly censored pesticide monitoring data to simulate daily pesticide concentrations and associated quantities used for acute and chronic exposure assessments, such as the annual maximum daily concentration. The new methodology is based on a statistical model that expresses log-transformed daily pesticide concentration in terms of a seasonal wave, flow-related variability, long-term trend, and serially correlated errors. Methods are described for estimating the model parameters, generating conditional simulations of daily pesticide concentration given sparse (weekly or less frequent) and potentially highly censored observations, and estimating concentration extremes based on the conditional simulations. The model can be applied to datasets with as few as 3 years of record, as few as 30 total observations, and as few as 10 uncensored observations. The model was applied to atrazine, carbaryl, chlorpyrifos, and fipronil data for U.S. Geological Survey pesticide sampling sites with sufficient data for applying the model. A total of 112 sites were analyzed for atrazine, 38 for carbaryl, 34 for chlorpyrifos, and 33 for fipronil. The results are summarized in this report; and, R functions, described in this report and provided in an accompanying model archive, can be used to fit the model parameters and generate conditional simulations of daily concentrations for use in investigations involving pesticide exposure risk and uncertainty.

Scientific Investigations Report

Water-quality trends in the nation’s rivers and streams, 1972–2012—Data preparation, statistical methods, and trend results

Since passage of the Clean Water Act in 1972, Federal, State, and local governments have invested billions of dollars to reduce pollution entering rivers and streams. To understand the return on these investments and to effectively manage and protect the Nation’s water resources in the future, we need to know how and why water quality has been changing over time. As part of the National Water-Quality Assessment Project, of the U.S. Geological Survey’s National Water-Quality Program, data from the U.S. Geological Survey, along with multiple other Federal, State, Tribal, regional, and local agencies, have been used to support the most comprehensive assessment conducted to date of surface-water-quality trends in the United States. This report documents the methods used to determine trends in water quality and ecology because these methods are vital to ensuring the quality of the results. Specific objectives are to document (1) the data compilation and processing steps used to identify river and stream sites throughout the Nation suitable for water-quality, pesticide, and ecology trend analysis, (2) the statistical methods used to determine trends in target parameters, (3) considerations for water-quality, pesticide, and ecology data and streamflow data when modeling trends, (4) sensitivity analyses for selecting data and interpreting trend results with the Weighted Regressions on Time, Discharge, and Season method, and (5) the final trend results at each site. The scope of this study includes trends in water-quality concentrations and loads (nutrient, sediment, major ion, salinity, and carbon), pesticide concentrations and loads, and metrics for aquatic ecology (fish, invertebrates, and algae) for four time periods: (1) 1972–2012, (2) 1982–2012, (3) 1992–2012, and (4) 2002–12. In total, nearly 12,000 trends in concentration, load, and ecology metrics were evaluated in this study; there were 11,893 combinations of sites, parameters, and trend periods. The final trend results are presented with examples of how to interpret the results from each trend model. Interpretation of the trend results, such as causal analysis, is not included.

Scientific Investigations Report

2011 Souris River flood—Will it happen again?

The Souris River Basin is a 61,000 square kilometer basin in the provinces of Saskatchewan and Manitoba and the state of North Dakota. Record setting rains in May and June of 2011 led to record flooding with peak annual streamflow values (762 cubic meters per second [m 3 /s]) more than twice that of any previously recorded peak streamflow and more than five times the estimated 100 year postregulation streamflow (142 m 3 /s) at the U.S. Geological Survey (USGS) streamflow-gaging station above Minot, North Dakota. Upstream from Minot, N. Dak., the Souris River is regulated by three reservoirs in Saskatchewan (Rafferty, Boundary, and Alameda) and Lake Darling in North Dakota. During the 2011 flood, the city of Minot, N. Dak., experienced devastating damages with more than 4,000 homes flooded and 11,000 evacuated. As a result, the Souris River Basin Task Force recommended the U.S. Geological Survey (in cooperation with the North Dakota State Water Commission) develop a model for estimating the probabilities of future flooding and drought. The model that was developed took on four parts: (1) looking at past climate, (2) predicting future climate, (3) developing a streamflow model in response to certain climatic variables, and (4) combining future climate estimates with the streamflow model to predict future streamflow events. By taking into consideration historical climate record and trends in basin response to various climatic conditions, it was determined flood risk will remain high in the Souris River Basin until the wet climate state ends.

Souris River

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

Tree-ring-based estimates of long-term seasonal precipitation in the Souris River Region of Saskatchewan, North Dakota and Manitoba

Historically unprecedented flooding occurred in the Souris River Basin of Saskatchewan, North Dakota and Manitoba in 2011, during a longer term period of wet conditions in the basin. In order to develop a model of future flows, there is a need to evaluate effects of past multidecadal climate variability and/or possible climate change on precipitation. In this study, tree-ring chronologies and historical precipitation data in a four-degree buffer around the Souris River Basin were analyzed to develop regression models that can be used for predicting long-term variations of precipitation. To focus on longer term variability, 12-year moving average precipitation was modeled in five subregions (determined through cluster analysis of measures of precipitation) of the study area over three seasons (November–February, March–June and July–October). The models used multiresolution decomposition (an additive decomposition based on powers of two using a discrete wavelet transform) of tree-ring chronologies from Canada and the US and seasonal 12-year moving average precipitation based on Adjusted and Homogenized Canadian Climate Data and US Historical Climatology Network data. Results show that precipitation varies on long-term (multidecadal) time scales of 16, 32 and 64 years. Past extended pluvial and drought events, which can vary greatly with season and subregion, were highlighted by the models. Results suggest that the recent wet period may be a part of natural variability on a very long time scale.

Manitoba, North Dakota, Saskatchewan

Stochastic model for simulating Souris River Basin precipitation, evapotranspiration, and natural streamflow

The Souris River Basin is a 61,000-square-kilometer basin in the Provinces of Saskatchewan and Manitoba and the State of North Dakota. In May and June of 2011, record-setting rains were seen in the headwater areas of the basin. Emergency spillways of major reservoirs were discharging at full or nearly full capacity, and extensive flooding was seen in numerous downstream communities. To determine the probability of future extreme floods and droughts, the U.S. Geological Survey, in cooperation with the North Dakota State Water Commission, developed a stochastic model for simulating Souris River Basin precipitation, evapotranspiration, and natural (unregulated) streamflow. Simulations from the model can be used in future studies to simulate regulated streamflow, design levees, and other structures; and to complete economic cost/benefit analyses. Long-term climatic variability was analyzed using tree-ring chronologies to hindcast precipitation to the early 1700s and compare recent wet and dry conditions to earlier extreme conditions. The extended precipitation record was consistent with findings from the Devils Lake and Red River of the North Basins (southeast of the Souris River Basin), supporting the idea that regional climatic patterns for many centuries have consisted of alternating wet and dry climate states. A stochastic climate simulation model for precipitation, temperature, and potential evapotranspiration for the Souris River Basin was developed using recorded meteorological data and extended precipitation records provided through tree-ring analysis. A significant climate transition was seen around1970, with 1912–69 representing a dry climate state and 1970–2011 representing a wet climate state. Although there were some distinct subpatterns within the basin, the predominant differences between the two states were higher spring through early fall precipitation and higher spring potential evapotranspiration for the wet compared to the dry state. A water-balance model was developed for simulating monthly natural (unregulated) mean streamflow based on precipitation, temperature, and potential evapotranspiration at select streamflow-gaging stations. The model was calibrated using streamflow data from the U.S. Geological Survey and Environment Canada, along with natural (unregulated) streamflow data from the U.S. Army Corps of Engineers. Correlation coefficients between simulated and natural (unregulated) flows generally were high (greater than 0.8), and the seasonal means and standard deviations of the simulated flows closely matched the means and standard deviations of the natural (unregulated) flows. After calibrating the model for a monthly time step, monthly streamflow for each subbasin was disaggregated into three values per month, or an approximately 10-day time step, and a separate routing model was developed for simulating 10-day streamflow for downstream gages. The stochastic climate simulation model for precipitation, temperature, and potential evapotranspiration was combined with the water-balance model to simulate potential future sequences of 10-day mean streamflow for each of the streamflow-gaging station locations. Flood risk, as determined by equilibrium flow-frequency distributions for the dry (1912–69) and wet (1970–2011) climate states, was considerably higher for the wet state compared to the dry state. Future flood risk will remain high until the wet climate state ends, and for several years after that, because there may be a long lag-time between the return of drier conditions and the onset of a lower soil-moisture storage equilibrium.

Manitoba, North Dakota, Saskatchewan

An evaluation of methods for estimating decadal stream loads

Effective management of water resources requires accurate information on the mass, or load of water-quality constituents transported from upstream watersheds to downstream receiving waters. Despite this need, no single method has been shown to consistently provide accurate load estimates among different water-quality constituents, sampling sites, and sampling regimes. We evaluate the accuracy of several load estimation methods across a broad range of sampling and environmental conditions. This analysis uses random sub-samples drawn from temporally-dense data sets of total nitrogen, total phosphorus, nitrate, and suspended-sediment concentration, and includes measurements of specific conductance which was used as a surrogate for dissolved solids concentration. Methods considered include linear interpolation and ratio estimators, regression-based methods historically employed by the U.S. Geological Survey, and newer flexible techniques including Weighted Regressions on Time, Season, and Discharge (WRTDS) and a generalized non-linear additive model. No single method is identified to have the greatest accuracy across all constituents, sites, and sampling scenarios. Most methods provide accurate estimates of specific conductance (used as a surrogate for total dissolved solids or specific major ions) and total nitrogen – lower accuracy is observed for the estimation of nitrate, total phosphorus and suspended sediment loads. Methods that allow for flexibility in the relation between concentration and flow conditions, specifically Beale’s ratio estimator and WRTDS, exhibit greater estimation accuracy and lower bias. Evaluation of methods across simulated sampling scenarios indicate that (1) high-flow sampling is necessary to produce accurate load estimates, (2) extrapolation of sample data through time or across more extreme flow conditions reduces load estimate accuracy, and (3) WRTDS and methods that use a Kalman filter or smoothing to correct for departures between individual modeled and observed values benefit most from more frequent water-quality sampling.

Journal of Hydrology

Impact of climate variability on runoff in the north-central United States

Large changes in runoff in the north-central United States have occurred during the past century, with larger floods and increases in runoff tending to occur from the 1970s to the present. The attribution of these changes is a subject of much interest. Long-term precipitation, temperature, and streamflow records were used to compare changes in precipitation and potential evapotranspiration (PET) to changes in runoff within 25 stream basins. The basins studied were organized into four groups, each one representing basins similar in topography, climate, and historic patterns of runoff. Precipitation, PET, and runoff data were adjusted for near-decadal scale variability to examine longer-term changes. A nonlinear water-balance analysis shows that changes in precipitation and PET explain the majority of multidecadal spatial/temporal variability of runoff and flood magnitudes, with precipitation being the dominant driver. Historical changes in climate and runoff in the region appear to be more consistent with complex transient shifts in seasonal climatic conditions than with gradual climate change. A portion of the unexplained variability likely stems from land-use change.

Iowa, Kansas, Minnesota, Missouri, Nebraska, North

Trends in major-ion constituents and properties for selected sampling sites in the Tongue and Powder River watersheds, Montana and Wyoming, based on data collected during water years 1980-2010

The primary purpose of this report is to present information relating to flow-adjusted temporal trends in major-ion constituents and properties for 16 sampling sites in the Tongue and Powder River watersheds based on data collected during 1980–2010. In association with this primary purpose, the report presents background information on major-ion characteristics (including specific conductance, calcium, magnesium, potassium, sodium adsorption ratio, sodium, alkalinity, chloride, fluoride, dissolved sulfate, and dissolved solids) of the sampling sites and coal-bed methane (CBM) produced water (groundwater pumped from coal seams) in the site watersheds, trend analysis methods, streamflow conditions, and factors that affect trend results. The Tongue and Powder River watersheds overlie the Powder River structural basin (PRB) in northeastern Wyoming and southeastern Montana. Limited extraction of coal-bed methane (CBM) from the PRB began in the early 1990’s, and increased dramatically during the late 1990’s and early 2000’s. CBM-extraction activities produce discharges of water with high concentrations of dissolved solids (particularly sodium and bicarbonate ions) relative to most stream water in the Tongue and Powder River watersheds. Water-quality of CBM produced water is of concern because of potential effects of sodium on agricultural soils and potential effects of bicarbonate on aquatic biota. Two parametric trend-analysis methods were used in this study: the time-series model (TSM) and ordinary least squares regression (OLS) on time, streamflow, and season. The TSM was used to analyze trends for 11 of the 16 study sites. For five sites, data requirements of the TSM were not met and OLS was used to analyze trends. Two primary 10-year trend-analysis periods were selected. Trend-analysis period 1 (water years 1986–95; hereinafter referred to as period 1) was selected to represent variability in major-ion concentrations in the Tongue and Powder River watersheds before potential effects of CBM-extraction activities. Trend analysis period 2 (water years 2001–10; hereinafter referred to as period 2) was selected because it encompassed substantial CBM-extraction activities and therefore might indicate potential effects of CBM-extraction activities on water quality of receiving streams in the Tongue and Powder River watersheds. For sites that did not satisfy data requirements for the TSM, OLS was used to analyze trends for period 2 (if complete data were available) or a 6-year period (2005–10). Flow-rate characteristics of CBM-produced water were estimated to allow general comparisons with streamflow characteristics of the sampling sites. The information on flow-rate characteristics of CBM-produced water in relation to streamflow does not account for effects of disposal, treatment, or other remediation activities on the potential quantitative effects of CBM-produced water on receiving streams. In many places, CBM-produced water is discharged into impoundments or channels in upper reaches of tributary watersheds where water infiltrates and does not directly contribute to streamflow. For Tongue River at State line (site 4) mean annual pumping rate of CBM-produced water during water years 2001–10 (hereinafter referred to as mean CBM pumping rate) was 6 percent of the mean of annual median streamflows during water years 2001–10 (hereinafter referred to as 2001–10 median streamflow). For main-stem Tongue River sites 5, 7, and 10, mean CBM pumping rate was 8–12 percent of 2001–10 median streamflow. For main-stem Powder River sites (sites 12, 13, and 16), mean CBM pumping rates were 26, 28, and 34 percent of 2001–10 median streamflows, respectively. For main-stem Tongue River sites analyzed by using the TSM and downstream from substantial CBM-extraction activities [Tongue River at State line (site 4), Tongue River at Tongue River Dam (site 5), Tongue River at Birney Day School (site 7), and Tongue River at Miles City (site 10)], generally small significant or nonsignificant decreases in most constituents are indicated for period 1. For period 2 for these sites, the TSM trend results do not allow confident conclusions concerning detection of effects of CBM-extraction activities on stream water quality. Detection of significant trends in major-ion constituents and properties for period 2 generally was infrequent, and direction, magnitudes, and significance of fitted trends were not strongly consistent with relative differences in water quality between stream water and CBM-produced water. The TSM indicated significant or generally large magnitude increases in median values of sodium adsorption ratio (SAR), sodium, and alkalinity for period 2 for sites 5 and 7, which might indicate potential effects of CBM-extraction activities on stream water. However, other factors, including operations of Tongue River Reservoir, irrigation activities, contributions of saline groundwater, and operations of the Decker coal mine, confound confident determination of causes of detected significant trends for sites 5 and 7. For all mainstem Tongue River sites, trends for period 2 generally are within ranges of those for period 1 before substantial CBM-extraction activities. For main-stem Powder River sites analyzed by using the TSM [Powder River at Sussex (site 11), Powder River at Arvada (site 12), Powder River at Moorhead (site 13), and Powder River near Locate (site 16)], significant or generally large magnitude decreases in median values of SAR, sodium, estimated alkalinity, chloride, fluoride, specific conductance, and dissolved solids are indicated for period 1. Patterns in trend results for period 1 for main-stem Powder River sites are consistent with effects of Salt Creek oil-brine reinjection that started in 1990. Trend results for all main-stem Powder River sites downstream from substantial CBM-extraction activities (sites 12, 13, and 16) indicate evidence of potential effects of CBM-extraction activities on stream water quality, although evidence is stronger for sites 12 and 13 than for site 16. Evidence in support of potential CBM effects includes significant increases in median values of SAR, sodium, and estimated alkalinity for period 2 for sites 12, 13, and 16 that are consistent with relative differences between stream water and CBM-produced water. Significant increases in median values of these constituents for period 2 are not indicated for Powder River at Sussex (site 11) upstream from substantial CBM-extraction activities. In interpreting the trend results, it is notable that the fitted trends evaluate changes in median concentrations and also notable that changes in median concentrations that might be attributed to CBM-extraction activities probably are more strongly evident during low to median streamflow conditions than during mean to high streamflow conditions. This observation is relevant in assessing trend results in relation to specific water-quality concerns, including effects of water-quality changes on irrigators and effects on stream biota and ecology.

Montana;Wyoming

Water-quality trends for selected sampling sites in the upper Clark Fork Basin, Montana, water years 1996-2010

A large-scale trend analysis was done on specific conductance, selected trace elements (arsenic, cadmium, copper, iron, lead, manganese, and zinc), and suspended-sediment data for 22 sites in the upper Clark Fork Basin for water years 1996–2010. Trend analysis was conducted by using two parametric methods: a time-series model (TSM) and multiple linear regression on time, streamflow, and season (MLR). Trend results for 1996–2010 indicate moderate to large decreases in flow-adjusted concentrations (FACs) and loads of copper (and other metallic elements) and suspended sediment in Silver Bow Creek upstream from Warm Springs. Deposition of metallic elements and suspended sediment within Warm Springs Ponds substantially reduces the downstream transport of those constituents. However, mobilization of copper and suspended sediment from floodplain tailings and stream banks in the Clark Fork reach from Galen to Deer Lodge is a large source of metallic elements and suspended sediment, which also affects downstream transport of those constituents. Copper and suspended-sediment loads mobilized from within this reach accounted for about 40 and 20 percent, respectively, of the loads for Clark Fork at Turah Bridge (site 20); whereas, streamflow contributed from within this reach only accounted for about 8 percent of the streamflow at Turah Bridge. Minor changes in FACs and loads of copper and suspended sediment are indicated for this reach during 1996–2010. Clark Fork reaches downstream from Deer Lodge are relatively smaller sources of metallic elements than the reach from Galen to Deer Lodge. In general, small decreases in loads and FACs of copper and suspended sediment are indicated for Clark Fork sites downstream from Deer Lodge during 1996–2010. Thus, although large decreases in FACs and loads of copper and suspended sediment are indicated for Silver Bow Creek upstream from Warm Springs, those large decreases are not translated to the more downstream reaches largely because of temporal stationarity in constituent transport relations in the Clark Fork reach from Galen to Deer Lodge. Unlike metallic elements, arsenic (a metalloid element) in streams in the upper Clark Fork Basin typically is mostly in dissolved phase, has less variability in concentrations, and has weaker direct relations with suspended-sediment concentrations and streamflow. Arsenic trend results for 1996–2010 indicate generally moderate decreases in FACs and loads in Silver Bow Creek upstream from Opportunity. In general, small temporal changes in loads and FACs of arsenic are indicated for Silver Bow Creek and Clark Fork reaches downstream from Opportunity during 1996–2010. Contribution of arsenic (from Warm Springs Ponds, the Mill-Willow bypass, and groundwater sources) in the Silver Bow Creek reach from Opportunity to Warm Springs is a relatively large source of arsenic. Arsenic loads originating from within this reach accounted for about 11 percent of the load for Clark Fork at Turah Bridge; whereas, streamflow contributed from within this reach only accounted for about 2 percent of the streamflow at Turah Bridge.

Montana

Modeled sulfate concentrations in North Dakota streams, 1993-2008, based on spatial basin characteristics

Sulfate concentration data collected from North Dakota streams during recent (1993–2008) years indicates generally higher sulfate concentrations across much of the State compared to concentrations during earlier years. The higher sulfate concentrations have been attributed in other studies to wetter climatic conditions, associated increases in contributing drainage areas, and rising water tables. The State’s current (2013) stream classification system, which includes a standard for 30-day average sulfate concentration, is based on earlier data and thus may not reflect natural conditions for more recent years. The U.S. Geological Survey, in cooperation with the North Dakota Department of Health and the North Dakota State Water Commission, completed a study to evaluate the relation of maximum seasonal (30-day moving average) sulfate concentrations during 1993–2008 to characteristics of the contributing basins to model expected naturally-occurring sulfate concentrations in North Dakota streams. Sulfate concentration data for 75 stream sampling sites in North Dakota were analyzed for this study. A spatial analysis was conducted with digital data using a Geographic Information System to obtain selected basin characteristics, which were in turn used as explanatory variables in a regression analysis to model the maximum seasonal (30-day moving average) sulfate concentration. Characteristics used in the regression analysis included mean annual precipitation, mean percent soil clay content, and mean percent saturation overland flow. Modeled sulfate concentrations generally were highest (greater than 750 milligrams per liter) in basins in western North Dakota and lowest (less than 250 milligrams per liter) in basins in the upper Sheyenne River and upper James River. Area-weighted means for the basin characteristics also were computed for 10-digit and 8-digit hydrologic units for streams in North Dakota and modeled sulfate concentrations were computed from the characteristics. The resulting distribution of modeled sulfate concentrations was similar to the distribution of estimates for the 12-digit hydrologic units, but less variable because the basin characteristics were averaged over larger areas.

North Dakota

Pesticide trends in major rivers of the United States, 1992-2010

This report is part of a series of pesticide trend assessments led by the National Water-Quality Assessment Program of the U.S. Geological Survey. This assessment focuses on major rivers of various sizes throughout the United States that have large watersheds with a range of land uses, changes in pesticide use, changes in management practices, and natural influences typical of the regions being drained. Trends were assessed at 59 sites for 40 pesticides and pesticide degradates during each of three overlapping periods: 1992–2001, 1997–2006, and 2001–10. In addition to trends in concentration, trends in agricultural-use intensity (agricultural use) were also assessed at 57 of the sites for 35 parent compounds with agricultural uses during the same three periods. The SEAWAVE-Q model was used to analyze trends in concentration, and parametric survival regression for interval-censored data was used to assess trends in agricultural use. All trends are provided in downloadable electronic files. A subset of 39 sites was chosen to represent non-nested, generally independent basins for a national analysis of pesticide and agricultural-use trends for the most prevalent pesticides (15 pesticides and 2 degradation products). Graphical and numerical results are presented to provide a national overview of concentration and use trends. As another perspective on understanding pesticide concentration trends in large rivers in relation to multiple tributary watersheds, this report also presents a detailed assessment of concentration and use trends for simazine, metolachlor, atrazine, deethylatrazine, and diazinon for a set of 17 nested sites in the Mississippi River Basin (including the Ohio and Missouri River Basins), for the second and third trend periods. Pesticides strongly dominated by agricultural use—cyanazine, metolachlor, atrazine, and alachlor—had widespread agreement between concentration trends and agricultural-use trends. Pesticides with substantial use in agricultural and urban applications—simazine, tebuthiuron, Dacthal, pendimethalin, chlorpyrifos, malathion, diazinon, fipronil, carbofuran, and carbaryl—had concentration trends that were mostly explained by a combination of agricultural-use trends and concentration trends in urban streams that were evaluated in a separate companion study. The importance of the urban stream trends for explaining concentration trends in major rivers indicates the significance of nonagricultural uses of some pesticides to concentrations in major rivers despite the much smaller area of urban land use compared to agriculture. Deethylatrazine, a degradate of atrazine, was the only pesticide compound assessed that had frequent occurrences during 1997–2006 and 2001–10 of concentration trends in the opposite direction of use trends (atrazine use). The nested analysis for the Mississippi River indicates that most trends observed in the largest rivers—multiple Mississippi River sites, the Ohio River, and the Missouri River—are consistent with streamflow contributions and concentration trends observed at tributary sites. Streamflow (incorporated into the trend model and shown in the nested basin analysis), trends in agricultural use of pesticides (quantified in this report), and urban use of pesticides (represented by concentration trends in a companion study of urban streams) are all important influences on pesticide concentrations in streams and rivers. Consideration of these influences is vital to understanding trends in pesticide concentrations.

Scientific Investigations Report