Geology ReportsSearch

SEARCH · Geology Reports

Results for “Trends”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3Linked to original sources

Water-quality characteristics and trends for selected wells possibly influenced by wastewater disposal at the Idaho National Laboratory, Idaho, 1981-2012

The U.S. Geological Survey, in cooperation with the U.S. Department of Energy, analyzed water-quality data collected from 64 aquifer wells and 35 perched groundwater wells at the Idaho National Laboratory (INL) from 1981 through 2012. The wells selected for the study were wells that possibly were affected by wastewater disposal at the INL. The data analyzed included tritium, strontium-90, major cations, anions, nutrients, trace elements, total organic carbon, and volatile organic compounds. The analyses were performed to examine water-quality trends that might influence future management decisions about the number of wells to sample at the INL and the type of constituents to monitor. The data were processed using custom computer scripts developed in the R programming language. Summary statistics were calculated for the datasets. Water-quality trends were determined using a parametric survival regression model to fit the observed data, including left-censored, interval-censored, and uncensored data. The null hypothesis of the trend test was that no relation existed between time and concentration; the alternate hypothesis was that time and concentration were related through the regression equation. A significance level of 0.05 was selected to determine if the trend was statistically significant. Trend test results for tritium and strontium-90 concentrations in aquifer wells indicated that nearly all wells had decreasing or no trends. Similarly, trends in perched groundwater wells were mostly decreasing or no trends; trends were increasing in two perched groundwater wells near the Advanced Test Reactor Complex. Decreasing trends generally are attributed to lack of recent wastewater disposal and radioactive decay. Trend test results for chloride, sodium, sulfate, nitrite plus nitrate (as nitrogen), chromium, trace elements, and total organic carbon concentrations in aquifer wells indicated that most wells had either decreasing or no trends. The decreasing trends in these constituents are attributed to decrease in disposal of these constituents, as well as discontinued use of the old percolation ponds south of the Idaho Nuclear Technology and Engineering Center (INTEC) and redirection of wastewater to the new percolation ponds 2 miles southwest of the INTEC in 2002. Chloride (along with sodium, sulfate, and some nitrate) concentrations in wells south of the INTEC may be influenced by episodic recharge from the Big Lost River. These constituent concentrations decrease during wetter periods when there is probably more recharge from the Big Lost River and increase during dry periods, when there is less recharge. Some wells downgradient of the Central Facilities Area and near the southern boundary of the INL showed increasing trends in sodium concentration, whereas there was no trend in chloride. The increasing trend for sodium could be due to the long term influence of wastewater disposal from upgradient facilities and the lack of trend for chloride could be because chloride is more mobile than sodium and more dispersed in the aquifer system. Volatile organic compound concentration trends were analyzed for nine aquifer wells. Trend test results indicated an increasing trend for carbon tetrachloride for the Radioactive Waste Management Complex Production Well for the period 1987–2012; however, trend analyses of data collected since 2005 show no statistically significant trend indicating that engineering practices designed to reduce movement of volatile organic compounds to the aquifer may be having a positive effect on the aquifer.

Idaho

Application of a Weighted Regression Model for Reporting Nutrient and Sediment Concentrations, Fluxes, and Trends in Concentration and Flux for the Chesapeake Bay Nontidal Water-Quality Monitoring Network, Results Through Water Year 2012

In the Chesapeake Bay watershed, estimated fluxes of nutrients and sediment from the bay’s nontidal tributaries into the estuary are the foundation of decision making to meet reductions prescribed by the Chesapeake Bay Total Maximum Daily Load (TMDL) and are often the basis for refining scientific understanding of the watershed-scale processes that influence the delivery of these constituents to the bay. Two regression-based flux and trend estimation models, ESTIMATOR and Weighted Regressions on Time, Discharge, and Season (WRTDS), were compared using data from 80 watersheds in the Chesapeake Bay Nontidal Water-Quality Monitoring Network (CBNTN). The watersheds range in size from 62 to 70,189 square kilometers and record lengths range from 6 to 28 years. ESTIMATOR is a constant-parameter model that estimates trends only in concentration; WRTDS uses variable parameters estimated with weighted regression, and estimates trends in both concentration and flux. WRTDS had greater explanatory power than ESTIMATOR, with the greatest degree of improvement evident for records longer than 25 years (30 stations; improvement in median model R 2 = 0.06 for total nitrogen, 0.08 for total phosphorus, and 0.05 for sediment) and the least degree of improvement for records of less than 10 years, for which the two models performed nearly equally. Flux bias statistics were comparable or lower (more favorable) for WRTDS for any record length; for 30 stations with records longer than 25 years, the greatest degree of improvement was evident for sediment (decrease of 0.17 in median statistic) and total phosphorus (decrease of 0.05). The overall between-station pattern in concentration trend direction and magnitude for all constituents was roughly similar for both models. A detailed case study revealed that trends in concentration estimated by WRTDS can operationally be viewed as a less-constrained equivalent to trends in concentration estimated by ESTIMATOR. Estimates of annual mean flow-adjusted (ESTIMATOR) and flow-normalized (WRTDS) concentration for years initially constituting the end of a water-quality record showed a similar degree of variability as data for additional years were incrementally added and the initial estimates “aged.” On the basis of the results of this broad comparison of the two models, the U.S. Geological Survey is adopting WRTDS as the primary model for estimating constituent fluxes and trends throughout the CBNTN. Nutrient and sediment flux and trend estimates, based on WRTDS, are summarized narratively and tabulated in appendixes for all stations for which fluxes or trends were reported through water year 2012. WRTDS also was used to explore the sensitivity of flux and trend estimates to three data-quality issues common in many large-scale monitoring networks and evident in some of the CBNTN records. The potential effects of inconsistency in annual sampling effort and inconsistency in storm sampling effort were explored by way of a subsampling experiment using eight of the most densely sampled long-term (1985–2012) stations in the CBNTN as baseline datasets. From each dataset, a set of 10 “design guideline” subsamples was selected, consisting of 12 monthly samples and 8 targeted storm samples per year. The selection was conducted in a manner that preserved the overall intensity of storm sampling in the baseline data. These 10 subsamples were further manipulated to create “heterogeneous” subsamples by removing storm samples prior to 2003. The maximum relative difference between flow-normalized flux estimated in a single year from any of the 10 design guideline subsamples and values estimated in the corresponding year from baseline data was smallest for dissolved inorganic nitrogen (median of 8 stations = 6 percent of baseline estimate), but more appreciable for total phosphorus and sediment (medians of 22 and 32 percent, respectively). The maximum relative difference between flow-normalized flux estimated from from the 10 heterogeneous subsamples and values estimated in the corresponding year from baseline data was more pronounced, with medians for 8 stations of 15, 30, and 53 percent of the corresponding baseline estimates for dissolved inorganic nitrogen, total phosphorus, and sediment, respectively. The worst-case maximum relative differences between flow-normalize flux estimated in a single year from the 10 heterogeneous subsamples and values estimated in the corresponding year from baseline data were 25 percent for dissolved inorganic nitrogen, 37 percent for total phosphorus, and 250 percent for sediment. The results for the heterogeneous subsamples indicate that changes in storm sampling frequency can result in appreciable distortion of estimated trends in flow-normalized flux, especially for total phosphorus and sediment. Trend lines estimated from heterogeneous subsamples tended to converge with the trend lines estimated from baseline data after 2003. In contrast, 2003–12 trends based on subsamples truncated by discarding all data prior to the induced heterogeneity in 2003 showed appreciable biases and differences in slope, relative to the corresponding 2003–12 segment of the trend computed from the design guideline subsamples. Overall, the results indicate that for particulate constituents, load and trend estimates computed using long-term records recently converted to CBNTN design guideline sampling protocols will be most reliable if the trend is computed using the entire record, but reported only for the period that design guideline sampling protocols were followed. Inconsistencies related to changing laboratory methods were also examined via two manipulative experiments. In the first experiment, increasing and decreasing “stair-step” patterns of changes in censoring level, overall representing a factor-of-five change in the laboratory reporting limit, were artificially imposed on a 27-year record with no censoring and a period-of-record concentration trend of –68.4 percent. Trends estimated on the basis of the manipulated records were broadly similar to the original trend (–63.6 percent for decreasing censoring levels and –70.3 percent for increasing censoring levels), lending a degree of confidence that the survival regression routines upon which WRTDS is based are generally robust to data censoring. The second experiment considered an abrupt disappearance of low-concentration observations of total phosphorus, associated with a laboratory method change and not reflected through censoring, near the middle of a 28-year record. By process of elimination, an upward shift in the estimated flow-normalize concentration trend line around the same time was identified as a likely artifact resulting from the laboratory method change, although a contemporaneous change in watershed processes cannot be ruled out. Decisions as to how to treat records with potential sampling protocol or laboratory methods-related artifacts should be made on a case-by-case basis, and trend results should be appropriately qualified.

Maryland, New York, Pennsylvania, Virginia

Trends in the quality of water in New Jersey streams, water years 1998-2007

Trends were determined in flow-adjusted values of selected water-quality characteristics measured year-round during water years 1998-2007 (October 1, 1997, through September 30, 2007) at 70 stations on New Jersey streams. Water-quality characteristics included in the analysis are dissolved oxygen, pH, total dissolved solids, total phosphorus, total organic nitrogen plus ammonia, and dissolved nitrate plus nitrite. In addition, trend tests also were conducted on measurements of dissolved oxygen made only during the growing season, April to September. Nearly all the water-quality data analyzed were collected by the New Jersey Department of Environmental Protection and the U.S. Geological Survey as part of the New Jersey Department of Environmental Protection Ambient Surface-Water Quality Monitoring Network. Monotonic trends in flow-adjusted values of water quality were determined by use of procedures in the ESTREND computer program. A 0.05 level of significance was selected to indicate a trend. Results of tests were not reported if there were an insufficient number of measurements or insufficient number of detected concentrations, or if the results of the tests were affected by a change in data-collection methods. Trends in values of dissolved oxygen, pH, and total dissolved solids were identified using the Seasonal Kendall test. Trends or no trends in year-round concentrations of dissolved oxygen were determined for 66 stations; decreases at 4 stations and increases at 0 stations were identified. Trends or no trends in growing-season concentrations of dissolved oxygen were determined for 65 stations; decreases at 4 stations and increases at 4 stations were identified. Tests of pH values determined trends or no trends at 26 stations; decreases at 2 stations and increases at 3 stations were identified. Trends or no trends in total dissolved solids were reported for all 70 stations; decreases at 0 stations and increases at 24 stations were identified. Trends in total phosphorus, total organic nitrogen plus ammonia, and dissolved nitrate plus nitrite were identified by use of Tobit regression. Two sets of trend tests were conducted-one set with all measurements and a second set with all measurements except the most extreme outlier if one could be identified. The result of the test with all measurements is reported if the results of the two tests are equivalent. The result of the test without the outlier is reported if the results of the two tests are not equivalent. Trends or no trends in total phosphorus were determined for 69 stations. Decreases at 12 stations and increases at 5 stations were identified. Of the five stations on the Delaware River included in this study, decreases in concentration were identified at four. Trends or no trends in total organic nitrogen plus ammonia were determined for 69 stations. Decreases and increases in concentrations were identified at six and nine stations, respectively. Trends or no trends in dissolved nitrate plus nitrite were determined for 66 stations. Decreases and increases in concentration were identified at 4 and 19 stations, respectively.

Scientific Investigations Report

Trends in Streamflow and Nutrient and Suspended-Sediment Concentrations and Loads in the Upper Mississippi, Ohio, Red, and Great Lakes River Basins, 1975-2004

Many actions have been taken to reduce nutrient and suspended-sediment concentrations and the amount of nutrients and sediment transported in streams as a result of the Clean Water Act and subsequent regulations. This report assesses how nutrient and suspended-sediment concentrations and loads in selected streams have changed during recent years to determine if these actions have been successful. Flow-adjusted and overall trends in concentrations and trends in loads from 1993 to 2004 were computed for total nitrogen, dissolved ammonia, total organic nitrogen plus ammonia, dissolved nitrite plus nitrate, total phosphorus, dissolved phosphorus, total suspended material (total suspended solids or suspended sediment), and total suspended sediment for 49 sites in the Upper Mississippi, Ohio, Red, and Great Lakes Basins. Changes in total nitrogen, total phosphorus, and total suspended-material loads were examined from 1975 to 2003 at six sites to provide a longer term context for the data examined from 1993 to 2004. Flow-adjusted trends in total nitrogen concentrations at 19 of 24 sites showed tendency toward increasing concentrations, and overall trends in total nitrogen concentrations at 16 of the 24 sites showed a general tendency toward increasing concentrations. The trends in these flow-adjusted total nitrogen concentrations are related to the changes in fertilizer nitrogen applications. Flow-adjusted trends in dissolved ammonia concentrations from 1993 to 2004 showed a widespread tendency toward decreasing concentrations. The widespread, downward trends in dissolved ammonia concentrations indicate that some of the ammonia reduction goals of the Clean Water Act are being met. Flow-adjusted and overall trends in total organic plus ammonia nitrogen concentrations from 1993 to 2004 did not show a distinct spatial pattern. Flow-adjusted and overall trends in dissolved nitrite plus nitrate concentrations from 1993 to 2004 also did not show a distinct spatial pattern. Flow-adjusted trends in total phosphorus concentrations were upward at 24 of 40 sites. Overall trends in total phosphorus concentrations were mixed and showed no spatial pattern. Flow-adjusted and overall trends in dissolved phosphorus concentrations were consistently downward at all of the sites in the eastern part of the basins studied. The reduction in phosphorus fertilizer use and manure production east of the Mississippi River could explain most of the observed trends in dissolved phosphorus. Flow-adjusted trends in total suspended-material concentrations showed distinct spatial patterns of increasing tendencies throughout the western part of the basins studied and in Illinois and decreasing concentrations throughout most of Wisconsin, Iowa, and in the eastern part of the basins studied. Flow-adjusted trends in total phosphorus were strongly related to the flow-adjusted trends in suspended materials. The trends in the flow-adjusted suspended-sediment concentrations from 1993 to 2004 resembled those for suspended materials. The long-term, nonmonotonic trends in total nitrogen, total phosphorus, and suspended-material loads for 1975 to 2003 were described by local regression, LOESS, smoothing for six sites. The statistical significance of those trends cannot be determined; however, the long-term changes found for annual streamflow and load data indicate that the monotonic trends from 1993 to 2004 should not be extrapolated backward in time.

Scientific Investigations Report

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

General water-quality conditions, long-term trends, and network analysis at selected sites within the Ambient Water-Quality Monitoring Network in Missouri, water years 1993–2017

The U.S. Geological Survey, in cooperation with the Missouri Department of Natural Resources, collects data pertaining to the surface-water resources of Missouri. Established in 1964, the Ambient Water-Quality Monitoring Network (AWQMN) consisted of 69 sites in 2017. Two additional sites from the National Water-Quality Program are included with the AWQMN sites for the analyses in this report. The sites are sampled typically from 2 to 12 times per year for physical properties, total suspended solids, nutrients, fecal indicator bacteria, and trace elements. The period of analysis for this study was from 1993 through 2017 and data analysis included 71 sites and 15 water-quality constituents plus discharge. Data analysis involved retrieving the data, conditioning the data for analysis, analyzing the data for trends, and analyzing the monitoring network to determine if potential data gaps or data redundancies exist in the network. Results from these analyses can be used to help manage the monitoring network into the future. Water-quality data were analyzed using several software packages to provide graphical and statistical information for interpretation of trends in the data at selected sites. Discharge data at selected sites were analyzed to determine the general trends during the analysis period and how the water-quality samples represented the range of daily mean discharges at each site. Water-quality data also were analyzed at selected sites to determine the relative sensitivity of selected sites and constituents to changes in data collection frequency. Trend analysis at selected sites using a simulated reduction in sampling frequency was completed to compare to trends obtained using monthly data to determine the potential degradation in the ability of determining trends from a reduced sampling frequency. The viability of using estimated discharge to evaluate long-term trends for sites with no continuous discharge was investigated. Data from sites were statistically compared in groups to determine the relative similarity (or difference) between sites for each water-quality constituent to identify potentially redundant sites in the monitoring network. Discharge-weighted long-term trends during 1993 through 2017 were analyzed for 15 water-quality constituents at 58 sites and results indicated there were significant single- or two-period trends in about 17 percent of the analyses. Some trends indicated improvement and some trends indicated deterioration of the general water quality at some sites in the AWQMN. No trend was indicated in about 31 percent of the analyses. The constituents pH, specific conductance, and total phosphorus showed the most frequent significant trends, and each of the 15 constituents examined had a significant trend at one or more sites. A total of 42 sites indicated at least 1 constituent with a significant single- or two-period trend, and 10 sites indicated 6 or more significant trends. Potential data gaps identified for computing discharge-weighted long-term trends in the monitoring network included the lack of collection of continuous discharge at 23 sites, insufficient sampling frequency for some constituents (dissolved chloride and total and dissolved lead and zinc) at some sites, insufficient temporal sample distribution (lack of at least one sample in each season per year) at some sites, and insufficient sampling frequency for some highly censored constituents (nutrients and total and dissolved lead and zinc) at some sites. Potential data gaps based on site spatial distribution were identified in 7 basins greater than 800 square miles. Potential site redundancies were identified in 4 basins that had an area greater than 500 square miles with a site density greater than 2 sites per 1,000 square miles. Potential site redundancies also were identified for nine site pairs by observing statistical similarities in the constituent data distributions. Sampling frequency was investigated to determine if reducing the sampling frequency of select constituents could provide a statistically similar data distribution. At 28 of 71 sites, 11 constituents had sufficient data collection frequency (approximately monthly) to allow for the creation of simulated datasets of various reduced data collection frequency. For the selected monitoring network sites analyzed, the data distribution of a simulated sampling frequency of four times per year or greater, roughly evenly distributed over the year, was not significantly different than the data distribution of the original monthly sampling frequency. Sites analyzed using varying simulated sampling frequencies tended to be more sensitive to sampling frequency changes if they were in basins classified as large or very large size and tended to be least sensitive in basins classified as small and medium size in the Ozark Plateaus Province. Simulated reduced frequency sampling analysis indicated that the constituents and measurements most sensitive to changes in sampling frequencies were water temperature, dissolved oxygen, discharge, and dissolved nitrate, and least sensitive were pH, total suspended solids, dissolved phosphorus, and total phosphorus. Discharge-weighted long-term trend analysis was repeated at 22 sites for 11 constituents using a simulated quarterly sampling frequency, and matched about 46 percent of the significant single-period trends identified using monthly data and about 65 percent of the analyses that indicated no trend using the monthly data.

Missouri

Analysis of nitrate (NO3-N) concentration trends in 25 ground-water-quality management areas, Idaho, 1961-2001

In Idaho, drinking-water supplies are pumped from relatively shallow ground-water zones where water quality has great potential for degradation by land- and water-use activities. One indicator of water quality, and one of the most widespread contaminants in Idaho ground water related to land and water uses, is dissolved nitrate. In December 2000, the U.S. Geological Survey, in cooperation with the Idaho Department of Environmental Quality, began a study to compile and assess nitrate data for ground water in 25 ground-water-quality management areas in Idaho. The primary objective of the study was to determine whether statistically significant trends in ground-water concentrations were discernible. Data included ground-water analyses but not analyses of water from springs, drains, or thermal water sources. A total of 8,465 nitrate analyses were compiled from 2,931 wells in the 25 priority areas; analyses dates ranged from June 1961 to February 2001. A time-period comparison was used to assess general trends in nitrate concentrations within individual priority areas. A time-series comparison was used to assess trends in nitrate concentrations from selected wells in each priority area. In time-period assessments, general trends in an area were evaluated by compiling nitrate data for selected time intervals to determine whether the population distributions were significantly different between one period and another. Data within each priority area were sorted into decades— 1970s, 1980s, and 1990s—for long-term trend assessment. The 1990s data also were divided into sets of selected years corresponding to Statewide Ambient Ground-Water Quality Program sampling cycles—1991 through 1994, 1995 through 1998, and a partial cycle, 1999 through 2000—for shortterm trend assessment. Data were analyzed by summary statistics, boxplots, and the Mann- Whitney statistical test. Long-term increasing trends in nitrate concentrations were evident for 6 of 25 priority areas, and long-term decreasing trends were evident for 4 of 25 areas. Short-term increasing trends were evident for 7 of 25 areas, and a short-term decreasing trend was evident for only 1 area. No long-term nitrate trends were evident for 7 of 25 areas, and no shortterm trends were evident for 15 of 25 areas. Data were insufficient for long-term trend assessment in 8 areas and for short-term assessment in 2 areas. Time-series trend analyses were conducted on data from wells with 7 or more nitrate analyses and longest periods of record, at least 10 years between oldest and most recent analyses. Because long-term records were available for so few wells and well construction data were not available for several of these wells, time-series analyses were not helpful to the nitrate trend assessment study. Trend results may be strongly affected by well construction, hydrogeologic environments, and changes in density and areal distribution of wells and analyses. The utility of nitrate trend assessments in priority areas would be improved by more consistent and specific well location descriptions between agencies; well construction and major water-yielding zone information to accompany the water-quality data for each well; and addition of historical nitrate data to data bases, particularly analyses prior to about 1990. Investigations of the possible effects of changing priority area boundaries and time periods on both long-term and short-term trend assessments are needed. Addition of these kinds of information would allow assessment of trends associated with hydrology and geology of each area and would provide a much stronger basis for trend assessment than currently possible.

Idaho

Water-quality trend analysis and sampling design for streams in North Dakota, 1971-2000

This report presents the results of a study conducted by the U.S. Geological Survey, in cooperation with the North Dakota Department of Health, to analyze historical water-quality trends in selected dissolved major ions, nutrients, and dissolved trace metals for 10 streams in southwestern and eastern North Dakota and to develop an efficient sampling design to monitor future water-quality trends. A time-series model for daily streamflow and constituent concentration was used to identify significant concentration trends, separate natural hydroclimatic variability in concentration from variability that could have resulted from anthropogenic causes, and evaluate various sampling designs to monitor future water-quality trends. The interannual variability in concentration as a result of variability in streamflow, referred to as the annual concentration anomaly, generally was high for all constituents and streams used in the trend analysis and was particularly sensitive to the severe drought that occurred in the late 1980's and the very wet period that began in 1993 and has persisted to the present (2002). Although climatic conditions were similar across North Dakota during the trend-analysis period (1971-2000), significant differences occurred in the annual concentration anomalies from constituent to constituent and location to location, especially during the drought and the wet period. Numerous trends were detected in the historical constituent concentrations after the annual concentration anomalies were removed. The trends within each of the constituent groups (major ions, nutrients, and trace metals) showed general agreement among the streams. For most locations, the largest dissolved major-ion concentrations occurred during the late 1970's and concentrations in the mid- to late 1990's were smaller than concentrations during the late 1970's. However, the largest concentrations for three of the Missouri River tributaries and one of the Red River of the North tributaries occurred during the mid- to late 1990's. Concentration trends for total ammonia plus organic nitrogen showed close agreement among the streams for which that constituent was evaluated. The largest concentrations occurred during the early 1980's, and the smallest concentrations occurred during the early 1990's. Nutrient data were not available for the early 1970's or late 1990's. Although a detailed analysis of the causes of the trends was beyond the scope of this report, a preliminary analysis of cropland, livestock-inventory, and oil-production data for 1971-2000 indicated the concentration trends may be related to the livestock-inventory and oil-production activities in the basins. Dissolved iron and manganese concentrations for the southwestern North Dakota streams generally remained stable during 1971-2000. However, many of the recorded concentrations for those streams were less than the detection limit, and trends that were masked by censoring may have occurred. Several significant trends were detected in dissolved iron and manganese concentrations for the eastern North Dakota streams. Concentrations for those streams either remained stable or increased during most of the 1970's and then decreased rapidly for about 2 years beginning in the late 1970's. The concentrations were relatively stable from the early 1980's to 2000 except at two locations where dissolved iron concentrations increased during the early 1990's. The most efficient overall sampling designs for the detection of annual trends (that is, trends that occur uniformly during the entire year) consisted of balanced designs in which the sampling dates and the number of samples collected remained fixed from year to year and in which the samples were collected throughout the year rather than in a short timespan. The best overall design for the detection of annual trends consisted of three samples per year, with samples collected near the beginning of December, April, and August. That design had acceptable sensitivity for the detection of trends in most constituents at all locations. Little improvement in sensitivity was achieved by collecting more than three samples per year. The sampling designs that were first evaluated for annual trends also were evaluated with regard to their sensitivity to detect seasonal trends that occurred during three seasons--April through August, August through December, and December through April. Design results indicated that an average of one extra sample per station per year resulted in an efficient design for detecting seasonal trends. However, allocation of the extra samples varied depending on the station, month, and constituent group (major ions, nutrients, and trace metals).

North Dakota

Water-quality trend analysis and sampling design for the Devils Lake Basin, North Dakota, January 1965 through September 2003

This report presents the results of a study conducted by the U.S. Geological Survey, in cooperation with the North Dakota State Water Commission, the Devils Lake Basin Joint Water Resource Board, and the Red River Joint Water Resource District, to analyze historical water-quality trends in three dissolved major ions, three nutrients, and one dissolved trace element for eight stations in the Devils Lake Basin in North Dakota and to develop an efficient sampling design to monitor the future trends. A multiple-regression model was used to detect and remove streamflow-related variability in constituent concentrations. To separate the natural variability in concentration as a result of variability in streamflow from the variability in concentration as a result of other factors, the base-10 logarithm of daily streamflow was divided into four components-a 5-year streamflow anomaly, an annual streamflow anomaly, a seasonal streamflow anomaly, and a daily streamflow anomaly. The constituent concentrations then were adjusted for streamflow-related variability by removing the 5-year, annual, seasonal, and daily variability. Constituents used for the water-quality trend analysis were evaluated for a step trend to examine the effect of Channel A on water quality in the basin and a linear trend to detect gradual changes with time from January 1980 through September 2003. The fitted upward linear trends for dissolved calcium concentrations during 1980-2003 for two stations were significant. The fitted step trends for dissolved sulfate concentrations for three stations were positive and similar in magnitude. Of the three upward trends, one was significant. The fitted step trends for dissolved chloride concentrations were positive but insignificant. The fitted linear trends for the upstream stations were small and insignificant, but three of the downward trends that occurred during 1980-2003 for the remaining stations were significant. The fitted upward linear trends for dissolved nitrite plus nitrate as nitrogen concentrations during 1987-2003 for two stations were significant. However, concentrations during recent years appear to be lower than those for the 1970s and early 1980s but higher than those for the late 1980s and early 1990s. The fitted downward linear trend for dissolved ammonia concentrations for one station was significant. The fitted linear trends for total phosphorus concentrations for two stations were significant. Upward trends for total phosphorus concentrations occurred from the late 1980s to 2003 for most stations, but a small and insignificant downward trend occurred for one station. Continued monitoring will be needed to determine if the recent trend toward higher dissolved nitrite plus nitrate as nitrogen and total phosphorus concentrations continues in the future. For continued monitoring of water-quality trends in the upper Devils Lake Basin, an efficient sampling design consists of five major-ion, nutrient, and trace-element samples per year at three existing stream stations and at three existing lake stations. This sampling design requires the collection of 15 stream samples and 15 lake samples per year rather than 16 stream samples and 20 lake samples per year as in the 1992-2003 program. Thus, the design would result in a program that is less costly and more efficient than the 1992-2003 program but that still would provide the data needed to monitor water-quality trends in the Devils Lake Basin.

North Dakota

Trends in peak flows of selected streams in Kansas

The possibility of a systematic change in flood potential led to an investigation of trends in the magnitude of annual peak flows in Kansas. Efficient design of highway bridges and other flood-plain structures depends on accurate understanding of flood characteristics. The Kendall's tau test was used to identify trends at 40 stream-gaging stations during the 40-year period 1958–97. Records from 13 (32 percent) of the stations showed significant trends at the 95-percent confidence level. Only three of the records (8 percent) analyzed had increasing trends, whereas 10 records (25 percent) had decreasing trends, all of which were for stations located in the western one-half of the State. An analysis of flow volume using mean annual discharge at 29 stations in Kansas resulted in 6 stations (21 percent) with significant trends in flow volumes. All six trends were decreasing and occurred in the western one-half of the State. The Kendall's tau test also was used to identify peak-flow trends over the entire period of record for 54 stream-gaging stations in Kansas. Of the 23 records (43 percent) showing significant trends, 16 (30 percent) were decreasing, and 7 (13 percent) were increasing. The trend test then was applied to 30-year periods moving in 5-year increments to identify time periods within each station record when trends were occurring. Systematic changes in precipitation patterns and long-term declines in ground-water levels in some stream basins may be contributing to peak-flow trends. To help explain the cause of the streamflow trends, the Kendall's tau test was applied to total annual precipitation and ground-water levels in Kansas. In western Kansas, the lack of precipitation and presence of decreasing trends in ground-water levels indicated that declining water tables are contributing to decreasing trends in peak streamflow. Declining water tables are caused by ground-water withdrawals and other factors such as construction of ponds and terraces. Peak-flow records containing trends introduce statistical error into flood-frequency analysis. To examine the effect of trends on flood-frequency analysis, statistically significant trends were added systematically to four nontrending station records. Flood magnitudes estimated on the basis of each data series were compared. The added trends resulted in changes in the 100-year flood magnitudes of as much as 70 percent.

Kansas

Nutrient concentrations and loads in the northeastern United States - Status and trends, 1975-2003

The U.S. Geological Survey (USGS) National Water-Quality Assessment Program (NAWQA) began regional studies in 2003 to synthesize information on nutrient concentrations, trends, stream loads, and sources. In the northeastern United States, a study area that extends from Maine to central Virginia, nutrient data were evaluated for 130 USGS water-quality monitoring stations. Nutrient data were analyzed for trends in flow-adjusted concentrations, modeled instream (non-flow-adjusted) concentrations, and stream loads for 32 stations with 22 to 29 years of water-quality and daily mean streamflow record during 1975-2003 (termed the long-term period), and for 46 stations during 1993-2003 (termed the recent period), by using a coupled statistical model of streamflow and water quality developed by the USGS. Recent trends in flow-adjusted concentrations of one or more nutrients also were analyzed for 90 stations by using Tobit regression. Annual stream nutrient loads were estimated, and annual nutrient yields were calculated, for 47 stations for the long-term and recent periods, and for 37 additional stations that did not have a complete streamflow and water-quality record for 1993-2003. Nutrient yield information was incorporated for 9 drainage basins evaluated in a national NAWQA study, for a total of 93 stations evaluated for nutrient yields. Long-term downward trends in flow-adjusted concentrations of total nitrogen and total phosphorus (18 and 19 of 32 stations, respectively) indicate regional improvements in nutrient-related water-quality conditions. Most of the recent trends detected for total phosphorus were upward (17 of 83 stations), indicating possible reversals to the long-term improvements. Concentrations of nutrients in many streams persist at levels that are likely to affect aquatic habitat adversely and promote freshwater or coastal eutrophication. Recent trends for modeled instream concentrations, and modeled reference concentrations, were evaluated relative to ecoregion-based nutrient criteria proposed by the U.S. Environmental Protection Agency. Instream concentrations of total nitrogen and total phosphorus persist at levels higher than proposed criteria at more than one-third and about one-half, respectively, of the 46 stations analyzed. Long-term trends in nutrient loads were primarily downward, with downward trends in total nitrogen and total phosphorus loads detected at 12 and 17 of 32 stations, respectively. Upward trends were rare, with one upward trend for total nitrogen loads and none for total phosphorus. Trends in loads of nitrite-plus-nitrate nitrogen included 7 upward and 8 downward trends among 32 stations. Downward trends in loads of ammonia nitrogen and total Kjeldahl nitrogen were detected at all six stations evaluated. Long-term downward trends detected in four of the five largest drainage basins evaluated include: total nitrogen loads for the Connecticut, Delaware, and James Rivers; total Kjeldahl nitrogen and ammonia nitrogen loads for the Susquehanna River; ammonia nitrogen and nitrite-plus-nitrate nitrogen loads for the James River; and total phosphorus loads for the Connecticut and Delaware Rivers. No trends in load were detected for the Potomac River. Nutrient yields were evaluated relative to the extent of land development in 93 drainage basins. The undeveloped land-use category included forested drainage basins with undeveloped land ranging from 75 to 100 percent of basin area. Median total nitrogen yields for the 27 undeveloped drainage basins evaluated, including 9 basins evaluated in a national NAWQA study, ranged from 290 to 4,800 pounds per square mile per year (lb/mi 2 /yr). Total nitrogen yields even in the most pristine drainage basins may be elevated relative to natural conditions, because of high rates of atmospheric deposition of nitrogen in parts of the northeastern United States. Median total phosphorus yields ranged from 12 to 330 lb/mi 2 /yr for the 26 undeveloped basins evaluated. The undeveloped category includes some large drainage basins with point-source discharges and small percentages of developed land; in these basins, streamflow from undeveloped headwater areas dilutes streamflow in more urbanized reaches, and dampens but does not eliminate the point-source "signal" of higher nutrient loads. Median total nitrogen yields generally do not exceed 1,700 lb/mi 2 /yr, and median total phosphorus yields generally do not exceed 100 lb/mi 2 /yr, in the drainage basins that are least affected by human land-use and waste-disposal practices. Agricultural and urban land use has increased nutrient yields substantially relative to undeveloped drainage basins. Median total nitrogen yields for 24 agricultural basins ranged from 1,700 to 26,000 lb/mi 2 /yr, and median total phosphorus yields ranged from 94 to 1,000 lb/mi 2 /yr. The maximum estimated total nitrogen and total phosphorus yields, 32,000 and 16,000 lb/mi 2 /yr, respectively, for all stations in the region were in small (less than 50 square miles (mi 2 )) agricultural drainage basins. Median total nitrogen yields ranged from 1,400 to 17,000 lb/mi 2 /yr in 26 urbanized drainage basins, and median total phosphorus yields ranged from 43 to 1,900 lb/mi 2 /yr. Urbanized drainage basins with the highest nutrient yields are generally small (less than 300 mi2) and are drained by streams that receive major point-source discharges. Instream nutrient loads were evaluated relative to loads from point-source discharges in four drainage basins: the Quinebaug River Basin in Connecticut, Massachusetts, and Rhode Island; the Raritan River Basin in New Jersey; the Patuxent River Basin in Maryland; and the James River Basin in Virginia. Long-term downward trends in nutrient loads, coupled with similar trends in flow-adjusted nutrient concentrations, indicate long-term reductions in the delivery of most nutrients to these streams. However, the absence of recent downward trends in load for most nutrients, coupled with instream concentrations that exceed proposed nutrient criteria in several of these waste-receiving streams, indicates that challenges remain in reducing delivery of nutrients to streams from point sources. During dry years, the total nutrient load from point sources in some of the drainage basins approached or equaled the nutrient load transported by the stream.

Connecticut;Delaware;Maine;Maryl;Massachusetts;New

Comparison of two regression-based approaches for determining nutrient and sediment fluxes and trends in the Chesapeake Bay watershed

Nutrient and sediment fluxes and changes in fluxes over time are key indicators that water resource managers can use to assess the progress being made in improving the structure and function of the Chesapeake Bay ecosystem. The U.S. Geological Survey collects annual nutrient (nitrogen and phosphorus) and sediment flux data and computes trends that describe the extent to which water-quality conditions are changing within the major Chesapeake Bay tributaries. Two regression-based approaches were compared for estimating annual nutrient and sediment fluxes and for characterizing how these annual fluxes are changing over time. The two regression models compared are the traditionally used ESTIMATOR and the newly developed Weighted Regression on Time, Discharge, and Season (WRTDS). The model comparison focused on answering three questions: (1) What are the differences between the functional form and construction of each model? (2) Which model produces estimates of flux with the greatest accuracy and least amount of bias? (3) How different would the historical estimates of annual flux be if WRTDS had been used instead of ESTIMATOR? One additional point of comparison between the two models is how each model determines trends in annual flux once the year-to-year variations in discharge have been determined. All comparisons were made using total nitrogen, nitrate, total phosphorus, orthophosphorus, and suspended-sediment concentration data collected at the nine U.S. Geological Survey River Input Monitoring stations located on the Susquehanna, Potomac, James, Rappahannock, Appomattox, Pamunkey, Mattaponi, Patuxent, and Choptank Rivers in the Chesapeake Bay watershed. Two model characteristics that uniquely distinguish ESTIMATOR and WRTDS are the fundamental model form and the determination of model coefficients. ESTIMATOR and WRTDS both predict water-quality constituent concentration by developing a linear relation between the natural logarithm of observed constituent concentration and three explanatory variables—the natural log of discharge, time, and season. ESTIMATOR uses two additional explanatory variables—the square of the log of discharge and time-squared. Both models determine coefficients for variables for a series of estimation windows. ESTIMATOR establishes variable coefficients for a series of 9-year moving windows; all observed constituent concentration data within the 9-year window are used to establish each coefficient. Conversely, WRTDS establishes variable coefficients for each combination of discharge and time using only observed concentration data that are similar in time, season, and discharge to the day being estimated. As a result of these distinguishing characteristics, ESTIMATOR reproduces concentration-discharge relations that are closely approximated by a quadratic or linear function with respect to both the log of discharge and time. Conversely, the linear model form of WRTDS coupled with extensive model windowing for each combination of discharge and time allows WRTDS to reproduce observed concentration-discharge relations that are more sinuous in form. Another distinction between ESTIMATOR and WRTDS is the reporting of uncertainty associated with the model estimates of flux and trend. ESTIMATOR quantifies the standard error of prediction associated with the determination of flux and trends. The standard error of prediction enables the determination of the 95-percent confidence intervals for flux and trend as well as the ability to test whether the reported trend is significantly different from zero (where zero equals no trend). Conversely, WRTDS is unable to propagate error through the many (over 5,000) models for unique combinations of flow and time to determine a total standard error. As a result, WRTDS flux estimates are not reported with confidence intervals and a level of significance is not determined for flow-normalized fluxes. The differences between ESTIMATOR and WRTDS, with regard to model form and determination of model coefficients, have an influence on the determination of nutrient and sediment fluxes and associated changes in flux over time as a result of management activities. The comparison between the model estimates of flux and trend was made for combinations of five water-quality constituents at nine River Input Monitoring stations. The major findings with regard to nutrient and sediment fluxes are as follows: (1)WRTDS produced estimates of flux for all combinations that were more accurate, based on reduction in root mean squared error, than flux estimates from ESTIMATOR; (2) for 67 percent of the combinations, WRTDS and ESTIMATOR both produced estimates of flux that were minimally biased compared to observed fluxes(flux bias = tendency to over or underpredict flux observations); however, for 33 percent of the combinations, WRTDS produced estimates of flux that were considerably less biased (by at least 10 percent) than flux estimates from ESTIMATOR; (3) the average percent difference in annual fluxes generated by ESTIMATOR and WRTDS was less than 10 percent at 80 percent of the combinations; and (4) the greatest differences related to flux bias and annual fluxes all occurred for combinations where the pattern in observed concentration-discharge relation was sinuous (two points of inflection) rather than linear or quadratic (zero or one point of inflection). The major findings with regard to trends are as follows: (1) both models produce water-quality trends that have factored in the year-to-year variations in flow; (2) trends in water-quality condition are represented by ESTIMATOR as a trend in flow-adjusted concentration and by WRTDS as a flow normalized flux; (3) for 67 percent of the combinations with trend estimates, the WRTDS trends in flow-normalized flux are in the same direction and magnitude to the ESTIMATOR trends in flow-adjusted concentration, and at the remaining 33 percent the differences in trend magnitude and direction are related to fundamental differences between concentration and flux; and (4) the majority (85 percent) of the total nitrogen, nitrate, and orthophosphorus combinations exhibited long-term (1985 to 2010) trends in WRTDS flow-normalized flux that indicate improvement or reduction in associated flux and the majority (83 percent) of the total phosphorus (from 1985 to 2010) and suspended sediment (from 2001 to 2010) combinations exhibited trends in WRTDS flow-normalized flux that indicate degradation or increases in the flux delivered.

Delaware, Maryland, New York, Pennsylvania, Virgin

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

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

Trends in precipitation and streamflow and changes in stream morphology in the Fountain Creek watershed, Colorado, 1939-99

The Fountain Creek watershed, located in and along the eastern slope of the Front Range section of the southern Rocky Mountains, drains approximately 930 square miles of parts of Teller, El Paso, and Pueblo Counties in eastern Colorado. Streamflow in the watershed is dominated by spring snowmelt runoff and storm runoff during the summer monsoon season. Flooding during the 1990?s has resulted in increased streambank erosion. Property loss and damage associated with flooding and bank erosion has cost area residents, businesses, utilities, municipalities, and State and Federal agencies millions of dollars. Precipitation (4 stations) and streamflow (6 stations) data, aerial photographs, and channel reconnaissance were used to evaluate trends in precipitation and streamflow and changes in channel morphology. Trends were evaluated for pre-1977, post-1976, and period-of-record time periods. Analysis revealed the lack of trend in total annual and seasonal precipitation during the pre-1977 time period. In general, the analysis also revealed the lack of trend in seasonal precipitation for all except the spring season during the post-1976 time period. Trend analysis revealed a significant upward trend in long-term (period of record) total annual and spring precipitation data, apparently due to a change in total annual precipitation throughout the Fountain Creek watershed. During the pre-1977 time period, precipitation was generally below average; during the post- 1976 time period, total annual precipitation was generally above average. During the post- 1976 time period, an upward trend in total annual and spring precipitation was indicated at two stations. Because two of four stations evaluated had upward trends for the post-1976 period and storms that produce the most precipitation are isolated convection storms, it is plausible that other parts of the watershed had upward precipitation trends that could affect trends in streamflow. Also, because of the isolated nature of convection storms that hit some areas of the watershed and not others, it is difficult to draw strong conclusions on relations between streamflow and precipitation. Trends in annual instantaneous peak streamflow, 70th percentile, 90th percentile, maximum daily-mean streamflow (100th percentile), 7-, 14-, and 30-day high daily-mean stream- flow duration, minimum daily-mean streamflow (0th percentile), 10th percentile, 30th percentile, and 7-, 14-, 30-day low daily-mean streamflow duration were evaluated. In general, instantaneous peak streamflow has not changed significantly at most of the stations evaluated. Trend analysis revealed the lack of a significant upward trend in streamflow at all stations for the pre-1977 time period. Trend tests indicated a significant upward trend in high and low daily-mean streamflow statistics for the post-1976 period. Upward trends in high daily-mean streamflow statistics may be an indication that changes in land use within the watershed have increased the rate and magnitude of runoff. Upward trends in low daily-mean 2 Trends in Precipitation and Streamflow and Changes in Stream Morphology in the Fountain Creek Watershed, Colorado, 1939-99 streamflow statistics may be related to changes in water use and management. An analysis of the relation between streamflow and precipitation indicated that changes in water management have had a marked effect on streamflow. Observable change in channel morphology and changes in distribution and density of vegetation varied with magnitude, duration, and frequency of large streamflow events, and increases in the magnitude and duration of low streamflows. Although more subtle, low stream- flows were an important component of day-to-day channel erosion. Substantial changes in channel morphology were most often associated with infrequent large or catastrophic streamflow events that erode streambed and banks, alter stream course, and deposit large amounts of sediment in the flood plain.

Colorado

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

Evaluation of long-term trends in hydrologic and water-quality conditions, and estimation of water budgets through 2013, Chester County, Pennsylvania

An evaluation of trends in hydrologic and water quality conditions and estimation of water budgets through 2013 was done by the U.S. Geological Survey in cooperation with the Chester County Water Resources Authority. Long-term hydrologic, meteorologic, and biologic data collected in Chester County, Pennsylvania, which included streamflow, groundwater levels, surface-water quality, biotic integrity, precipitation, and air temperature were analyzed to determine possible trends or changes in hydrologic conditions. Statistically significant trends were determined by applying the Kendall rank correlation test; the magnitudes of the trends were determined using the Sen slope estimator. Water budgets for eight selected watersheds were updated and a new water budget was developed for the Marsh Creek watershed. An average water budget for Chester County was developed using the eight selected watersheds and the new Marsh Creek water budget. Annual and monthly mean streamflow, base flow, and runoff were analyzed for trends at 10 streamgages. The periods of record at the 10 streamgages ranged from 1961‒2013 to 1988‒2013. The only statistically significant trend for annual mean streamflow was for West Branch Brandywine Creek near Honey Brook, Pa. (01480300) where annual mean streamflow increased 1.6 cubic feet per second (ft 3 /s) per decade. The greatest increase in monthly mean streamflow was for Brandywine Creek at Chadds Ford, Pa. (01481000) for December; the increase was 47 ft 3 /s per decade. No statistically significant trends in annual mean base flow or runoff were determined for the 10 streamgages. The greatest increase in monthly mean base flow was for Brandywine Creek at Chadds Ford, Pa. (01481000) for December; the increase was 26 ft 3 /s per decade. The magnitude of peaks greater than a base streamflow was analyzed for trends for 12 streamgages. The period of record at the 12 stream gages ranged from 1912‒2012 to 2004–11. Fifty percent of the streamgages showed a small statistically significant increase in peaks greater than the base streamflow. The greatest increase was for Brandywine Creek at Chadds Ford, Pa. (01481000) during 1962‒2012; the increase was 1.8 ft 3 /s per decade. There were no statistically significant trends in the number of floods equal to or greater than the 2-year recurrence interval flood flow. Twenty‒one monitoring wells were evaluated for statistically significant trends in annual mean water level, minimum annual water level, maximum annual water level, and annual range in water-level fluctuations. For four wells, a small statistically significant increase in annual mean water level was determined that ranged from 0.16 to 0.7 feet per decade. There was poor or no correlation between annual mean groundwater levels and annual mean streamflow and base flow. No correlation was determined between annual mean groundwater level and annual precipitation. Despite rapid population growth and land-use change since 1950, there appears to have been little or no detrimental effects on groundwater levels in 21 monitoring wells. Long-term precipitation and temperature data were available from the West Chester (1893‒2013) and Phoenixville, Pa. (1915‒2013) National Oceanic and Atmospheric Administration (NOAA) weather stations. No statistically significant trends in annual mean precipitation or annual mean temperature were determined for either station. Both weather stations had a significant decrease in the number of days per year with precipitation greater than or equal to 0.1 inch. Annual mean minimum and maximum temperatures from the NOAA Southeastern Piedmont Climate Division increased 0.2 degrees Fahrenheit (F) per decade between 1896 and 2014. The number of days with a maximum temperature equal to or greater than 90 degrees F increased at West Chester and decreased at Phoenixville. No statistically significant trend was determined for annual snowfall amounts. Data from 1974 to 2013 for three stream water-quality monitors in the Brandywine Creek watershed were evaluated. The monitors are on the West Branch Brandywine Creek at Modena, Pa. (01480617), East Branch Brandywine Creek below Downingtown, Pa. (01480870), and Brandywine Creek at Chadds Ford, Pa. (01481000). Statistically significant upward trends were determined for annual mean specific conductance at all three stations, indicating the total dissolved solids load has been increasing. If the current trend continues, the annual mean specific conductance could almost double from 1974 to 2050. The increase in specific conductance likely is due to increases in chloride concentrations, which have been increasing steadily over time at all three stations. No correlation was found between monthly mean specific conductance and monthly mean streamflow or base flow. Statistically significant upward trends in pH were determined for all three stations. Statistically significant upward trends in stream temperature were determined for East Branch Brandywine Creek below Downingtown, Pa. (01480870) and Brandywine Creek at Chadds Ford, Pa. (01481000). The stream water-quality data indicate substantial increases in the minimum daily dissolved oxygen concentrations in the Brandywine Creek over time. The Chester County Index of Biotic Integrity (CC-IBI) determined for 1998‒2013 was evaluated for the five biological sampling sites collocated with streamgages. CC-IBI scores are based on a 0‒100 scale with higher scores indicating better stream quality. Statistically significant upward trends in the CC-IBI were determined for West Branch Brandywine Creek at Modena, Pa. (01480617) and East Branch Brandywine Creek below Downingtown, Pa. (01480870). No correlation was found between the CC-IBI and streamflow, precipitation, or stream specific conductance, pH, temperature, or dissolved oxygen concentration. A Chester County average water budget was developed using the nine estimated watershed water budgets. Average precipitation was 48.4 inches, and average streamflow was 21.4 inches. Average runoff and base flow were 8.3 and 13.1 inches, respectively, and average evapotranspiration and estimation of errors was 27.2 inches."

Pennsylvania

Consequences of ignoring spatial variation in population trend when conducting a power analysis

Long-term, large-scale monitoring programs are becoming increasingly common to document status and trends of wild populations. A successful program for monitoring population trend hinges on the ability to detect the trend of interest. Power analyses are useful for quantifying the sample size needed for trend detection, given expected variation in the population. Four components of variation (within-year variation at a given site, interannual variation within a site, variation among sites in the interannual variation, and variation among sites in mean abundance or density) are commonly considered in power analyses for population trend, but a fifth is not: variation among sites in the local trend. Spatial variation in trend is expected to reduce statistical power, but the magnitude of this reduction has not been fully explored. We used computer simulations to evaluate the consequences of ignoring spatial variation in trend under a variety of sampling designs and wide ranges of other components of variation. The effect of spatial variation in trend on power was minor when other input parameters took extreme values that made the trend either very difficult or very easy to detect. However, at moderate values of the other parameters, spatial variation in trend had a strong effect, reducing statistical power by up to 60%. In some cases, ignoring spatial variation in trend resulted in an 80% probability of a Type I error (falsely detecting a trend in a stable population). Spatial variation in trend is therefore an important consideration when designing a long-term monitoring program for many species, especially those affected by local conditions at sites that are repeatedly surveyed. If variation in trend is ignored, as in most previous power analyses, the recommended sampling design will likely be insufficient to detect the trend of interest and lead to potentially false conclusions of a stable population.

Ecography