Search USGS⌕ Search

Geology topics

Karen R. Ryberg

Publications and source records attributed to Karen R. Ryberg.

At least 19 recordsLinked to original sources

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

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

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

A framework for understanding the effects of subsurface agricultural drainage on downstream flows

Understanding controls on streamflow volume and magnitude is important to water resource management applications, such as critical water and transportation structure design and floodplain mapping. Changes in land use and agricultural practices, such as subsurface agricultural drainage, may be contributing to changes in streamflow characteristics. Subsurface agricultural drainage, also known as tile drainage, is the practice of installing drains in the subsurface of agricultural fields to improve productivity. Because of the complex interactions between subsurface drainage systems, precipitation, local soil conditions, and land management practices, it is difficult to determine how subsurface agricultural drainage affects downstream flow. Previously developed subsurface agricultural drainage conceptual models under dry, saturated, and winter conditions are summarized, and current literature on the effects of subsurface agricultural drainage on downstream flows, focusing on peak flow, non-event flow, and total flow to develop frameworks for discussing these systems is compiled. The effects that subsurface drainage has on hydrologic systems are expected to vary by site and are seasonally based on system design, soil type, moisture conditions, precipitation characteristics, and land conditions. Subsurface drainage can affect the magnitude of peak flow by converting surface runoff from a storm event to subsurface runoff. By increasing hydrologic connectivity of a catchment, subsurface drainage can increase non-event flow or the flow between two storm events, typically dependent on lateral flow through the subsurface and groundwater. Theoretically, by diverting water from groundwater recharge or by reducing water available for evapotranspiration, subsurface drainage may increase the total volume of flow. Precipitation changes may increase infiltration, excess overland flow, and flood risk regardless of the presence or absence of subsurface drainage.

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

Evaluating drought risk of the Red River of the North Basin using historical and stochastic streamflow upstream from Emerson, Manitoba

Drought and its effect on streamflow are important to understand because of the potential to adversely affect water supply, agricultural production, and ecological conditions. The Red River of the North Basin in north-central United States and central Canada is susceptible to dry conditions. During an extended drought, streamflow conditions in the Red River of the North may become inadequate to support existing water supply needs in the basin for agriculture, industry, human use, and aquatic life. To understand potential future low-streamflow conditions in the Red River of the North Basin, the U.S. Geological Survey, in cooperation with the International Joint Commission, North Dakota Department of Water Resources, Red River Joint Water Resource District, and Red River Watershed Management Board, developed a water-balance model of the Red River of the North Basin upstream from Emerson, Manitoba, Canada, and coupled the model with stochastic weather inputs to simulate possible future low-streamflow conditions. Historical changes in low-streamflow conditions were characterized across the Red River of the North Basin using multiple change-point analysis for 12 streamgages. Across these stations, significant change-point years in 1943 and 1994 marked increases in the magnitude of low-streamflow conditions. During 1920–2015, conversion of primary land (not affected by human use) to agricultural and secondary land was followed by a conversion from smalls grains to corn and soybeans as the dominant crop type. From land-use analysis, 1940–2000 was determined to have relatively stable land use and therefore was used as the calibration period for the water-balance model. A deterministic water-balance model was developed for the Red River of the North Basin upstream from Emerson, Manitoba. The water-balance model was calibrated with data from 37 U.S. Geological Survey streamgages for 1940–2000 and verified using data for 2001–15. The calibrated water-balance model simulated streamflow distributions that mirrored the seasonal patterns of the observed mean monthly streamflow and the standard deviation of the monthly streamflow data, especially during the fall and winter months when streamflow was lowest. For the verification period, during the low-streamflow months of December through January, the difference between simulated and observed data was similar to the calibration comparison and successfully reproduced seasonal trends in the distribution of streamflow, even when using weather data that were outside the calibration period. To determine the future risk of low-streamflow conditions in the Red River of the North Basin, a block-bootstrap method was used to generate multiple possible future climates. These stochastically generated weather time series were then input to a water-balance model to simulate a distribution of possible streamflows. Three sets of experiments were performed, with each experiment containing a set of scenarios. The first set of experiments from the stochastic streamflow model were designed to investigate how changes in reservoir management would affect the distribution of low streamflow. Relative to scenario 1 (present-day [2023] reservoir operation), scenario 2 (no reservoir operation) shifted the low-streamflow frequency curves downward, reducing the annual minimum monthly streamflow for the Emerson subbasin. Subbasins were defined by the contributing area upstream from a selected streamgage station. Relative to scenario 1, scenario 3 (regulated streamflow with an increased reservoir capacity of 10 percent) shifted the low-streamflow frequency curves upward for the Emerson subbasin. The magnitude of this upward shift, caused by increased reservoir capacity, was lower than the magnitude of the shift caused by the absence of the reservoirs, which indicates that the streamflow was most affected when the reservoirs were first constructed. The second set of experiments from the stochastic streamflow model included two scenarios that were performed to better understand how the Red River of the North Basin responds to long periods of low or high precipitation. The results indicate that the model consistently overestimated streamflow, but the relative change between a wet and dry climate state of simulated streamflow distribution reasonably matched the relative change of historical streamflow. Across the subbasins, the model was most accurate for low-streamflow conditions associated with nonexceedance probabilities between 20 and 40 percent. The third set of experiments from the stochastic streamflow model were done to investigate low-streamflow response across the basin to several drought events. Low-end streamflow was reduced when the basin was exposed to a drought, and the magnitude of the reduction increased with longer or more intense droughts. Compared to the low-intensity drought scenarios, the range of percent reductions (as indicated by the interquartile range) was larger for the high-intensity drought scenarios for all subbasins, and the subbasins of Grand Forks and Emerson had a smaller range of reductions compared to the other three subbasins. The larger drainage area—combined with the large contribution of the Red Lake River and several other Minnesota tributaries that generally experience wetter climate conditions—upstream from the Emerson and Grand Forks subbasins may contribute to the smaller range in reductions under the high intensity scenarios. Comparison of the percent reduction in low-end streamflow among subbasins also indicated that the effects of drought duration and intensity could be cumulative. Combining factors of time and intensity produced a larger reduction in streamflow than when each effect was isolated. The array of drought scenarios can be used to determine how a subbasin would respond to multiple possible future conditions. Based on climate predictions, the drought scenario that best matches a future anticipated drought scenario can be used to estimate a low streamflow response for a given subbasin.

Manitoba, Minnesota, North Dakota↗

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

Standardized guidelines for completing flood-flow frequency analyses are presented in a U.S. Geological Survey Techniques and Methods report known as Bulletin 17C, https://doi.org/10.3133/tm4B5 . In recent decades (since about 2000), a better understanding of long-term climatic persistence (periods of clustered floods or droughts, or wet or dry periods) and concerns about potential climate change and land-use change have caused a reexamination of the stationarity assumptions underlying methods in Bulletin 17C. Bulletin 17C does not offer guidance on incorporating nonstationarities and further identifies a need for flood-frequency studies that incorporate changing climate or basin characteristics. As part of that reexamination, a study of annual peak streamflow (peak flow) has begun in the Midwest. This chapter of the study summarizes how hydroclimatic variability affects peak flows in North Dakota. In this analysis of peak flow, daily streamflow, and climate metrics, four periods were selected: (1) a 100-year period, 1921–2020; (2) a 75-year period, 1946–2020; (3) a 50-year period, 1971–2020; and (4) a 30-year period, 1991–2020. Output from a monthly water-balance model was used for the climate data. Statistical analysis of peak flow consisted of evaluations of autocorrelation, trends, and change points and was augmented with analyses of seasonality and daily streamflow. The long-term pattern of decreasing peak flow in the west and increasing peak flow in the east is a pattern of opposing signals on either side of the 100th meridian. Analyses indicate that a key factor in changing hydroclimatology is the increase in fall precipitation. The trends in soil moisture closely match the trends in annual precipitation. Nonstationary flood-frequency analysis necessitates detailed exploratory data analysis and additional data and information about climate, land use, and other factors. This study provides extensive exploratory analysis for peak flow, daily streamflow, and climate data for North Dakota, setting the stage for informed nonstationary flood-frequency analysis.

North Dakota↗

Introduction and methods of analysis for peak streamflow trends and their relation to changes in climate in Illinois, Iowa, Michigan, Minnesota, Missouri, Montana, North Dakota, South Dakota, and Wisconsin

Flood-frequency analysis, also called peak-flow frequency or flood-flow frequency analysis, is essential to water resources management applications including critical structure design and floodplain mapping. Federal guidelines for doing flood-frequency analyses are presented in a U.S. Geological Survey Techniques and Methods Report known as Bulletin 17C. A basic assumption within Bulletin 17C is that for drainage basins without major hydrologic alterations, statistical properties of the distribution of annual peak streamflows (peak flows) are stationary; that is, the mean, variance, and skew are constant. The stationarity assumption has been widely accepted within the flood-frequency community; however, a better understanding of long-term climatic persistence and concerns about potential climate change and land-use change has caused a reexamination of the stationarity assumption. Flood-frequency analyses that do not incorporate observed trends and abrupt changes may result in a poor representation of the true flood risk. Bulletin 17C does not offer guidance for incorporating nonstationarities when estimating floods, and it describes a need for studies that incorporate changing climate or basin characteristics. In response to this need and a history of concern regarding nonstationarity peak flows in the region, this study was done by the U.S. Geological Survey, in cooperation with the Departments of Transportation of Illinois, Iowa, Michigan, Minnesota, Missouri, South Dakota, and Wisconsin; the Montana Department of Natural Resources and Conservation; and the North Dakota Department of Water Resources, to assess potential nonstationarity in peak flows in the north-central United States. This chapter summarizes the methods used to detect hydroclimatic changes in peak-flow data in the study region. A wide range of analyses and statistical approaches are applied to document the primary mechanisms controlling floods and characterize temporal changes in hydroclimatic variables and peak flow. Four periods were selected for analysis of peak flow, daily streamflow, and climate data. The periods are (1) a 100-year period, 1921–2020; (2) a 75-year period, 1946–2020; (3) a 50-year period, 1971–2020; and (4) a 30-year period, 1991–2020. The climate data consist of monthly time series estimates of temperature, precipitation, potential evapotranspiration, actual evapotranspiration, snowfall, soil moisture storage, snow water equivalent, and runoff on a 3.1-mile by 3.1-mile grid for the conterminous United States. Statistical and graphical analyses were used to investigate potential changes in hydrology and climate. The starting point for these analyses was the initial data analysis of peak flow described in Bulletin 17C, which includes plotting the peak flow and checking for autocorrelation, monotonic trends, and changes points. Analyses were added to examine additional features in the data. To examine potential causal drivers of changes, the climate data were analyzed graphically and statistically. Results are provided in a U.S. Geological Survey data release. The study limitations are documented for users of the results.

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

Assessing the added value of antecedent streamflow alteration information in modeling stream biological condition

In stream systems, disentangling relationships between biology and flow and subsequent prediction of these relationships to unsampled streams is a common objective of large-scale ecological modeling. Often, streamflow metrics are derived from aggregating continuous streamflow records available at a subset of stream gages into long-term flow regime descriptors. Despite demonstrated value, shortcomings of these long-term approaches include spatial restriction to locations with long-term continuous flow records (commonly, biased toward larger systems) and omission of potentially ecologically important short-term (i.e., ≤1 year) antecedent streamflow information. We used long-term flow regime and short-term antecedent streamflow alteration information to evaluate relative performance in modeling stream fish biological condition. We compared results to understand whether short-term antecedent streamflow information improved models of fish biological condition. Results indicated that models incorporating short-term antecedent data performed better than those relying solely on long-term flow regime data (kappa statistic = 0.29 and 0.23, respectively) and improved prediction accuracy among stream sizes and in six of nine ecoregions. Additionally, models relying solely on short-term streamflow information performed similarly to those with only long-term streamflow information (kappa = 0.23). Incorporating short-term antecedent streamflow metrics may provide added ecological information not fully captured by long-term flow regime summaries in macroscale modeling efforts or perform similarly to long-term streamflow data when long-term data are not available.

Science of the Total Environment↗

Method for identification of reservoir regulation within U.S. Geological Survey streamgage basins in the Central United States using a decadal dam impact metric

Researchers routinely study streamflow data to understand the effects of natural climate variability and anthropogenic climate change, and to develop methods for estimating streamflow at ungaged locations. These studies require streamflow data that are not modified or largely altered by other anthropogenic activities, such as reservoirs or diversions. This report discusses a method for identifying basins with reservoir regulation using a decadal impact metric that characterizes the degree of regulation of a given river reach. The method is applied to U.S. Geological Survey streamgage basins from eight States in the Central United States. Using this metric, 140 streamgages with known regulation effects (annual peak streamflow values qualified with a code 6) were evaluated for their impact metric values in decades with annual peak streamflow values qualified with code 6. Based on the distribution of median impact metric values at these regulated basins, a threshold value of 0.1 was identified as the value that when exceeded was the most characteristic of the regulated streamgage basins in the study area. Streamgage basins from nine States with peak streamflow values that were not qualified with code 6 were evaluated for impact metric values equal to or greater than the established threshold. About 13 percent of streamgages (136 of 1,017) had an impact metric equal to or greater than the identified regulated threshold at some point in their periods of record. The method discussed in this report, which has limitations owing to characteristics of the data underlying the dam impact metric, provides a regionally consistent approach to identifying regulated U.S. Geological Survey streamgage basins.

Iowa, Illinois, Michigan, Missouri, Montana, North↗

Flood-frequency analysis in the Midwest: Addressing potential nonstationarity of annual peak-flow records

Flood-frequency analysis is essential in numerous water-resource management applications, including critical structure design and flood-plain mapping. A basic assumption within Bulletin 17C [1], the standardized guidelines for conducting flood-frequency analysis, is that basins without major hydrologic alterations, such as regulation or urbanization, exhibit stationary statistical properties of the distribution of annual peak streamflow. That is, the mean, variance, and skew are constant over time and the peak-flow record is a representative sample of the population of future floods [1]. In recent decades, better understanding of long-term climatic persistence and concerns about climate and land-use change have caused the assumption of stationarity in peak-flow records to be reexamined [2, 3, 4, 5]. Under nonstationary conditions, the long-term distributional properties (mean, variance, and/or skew) of peak-flow series change one or more times, either gradually or abruptly. Nonstationarities may be attributed to one source, but are often a result of a mixture of drivers, making detection and attribution of nonstationarities challenging [6, 7, 8]. Failure to incorporate observed trends and abrupt changes into flood-frequency analysis may result in a poor representation of the true flood risk. Bulletin 17C currently offers no guidance on how to account for nonstationarities when estimating floods and acknowledges the benefit additional flood frequency studies that incorporate changing climate or basin characteristics into the analysis would provide[1].

Midwest↗

Projecting flood frequency curves under near-term climate change

Flood-frequency curves, critical for water infrastructure design, are typically developed based on a stationary climate assumption. However, climate changes are expected to violate this assumption. Here, we propose a new, climate-informed methodology for estimating flood-frequency curves under non-stationary future climate conditions. The methodology develops an asynchronous, semiparametric local-likelihood regression (ASLLR) model that relates moments of annual maximum flood to climate variables using the generalized linear model. We estimate the first two marginal moments (MM) – the mean and variance – of the underlying log-Pearson Type-3 distribution from the ASLLR with the monthly rainfall and temperature as predictors. The proposed methodology, ASLLR-MM, is applied to 40 U.S. Geological Survey streamgages covering 18 water resources regions across the conterminous United States. A correction based on the aridity index was applied on the estimated variance, after which the ASLLR-MM approach was evaluated with both historical (1951–2005) and projected (2006–2035, under RCP4.5 and RCP8.5) monthly precipitation and temperature from eight Global Circulation Models (GCMs) consisting of 39 ensemble members. The estimated flood-frequency quantiles resulting from the ASLLR-MM and GCM members compare well with the flood-frequency quantiles estimated using the historical period of observed climate and flood information for humid basins, whereas the uncertainty in model estimates is higher in arid basins. Considering additional atmospheric and land-surface conditions and a multi-level model structure that includes other basins in a region could further improve the model performance in arid basins.

Water Resources Research↗

Assessment of streamflow trends in the eastern Dakotas, water years 1960–2019

Hydrologic extremes, whether periods of drought or flooding, are occurring more frequently with greater severity and can have substantial economic impacts. Along with flooding, the timing and volume of streamflow also is changing across the United States. The focus of this report is to characterize a unique trend in mean annual streamflow occurring in eastern North and South Dakota, hereafter referred to as the eastern Dakotas, that is not being observed anywhere else in the conterminous United States. Streamflow records for 1,853 U.S. Geological Survey streamgages obtained from the U.S. Geological Survey National Water Information System database with a continuous record of mean annual streamflow during water years 1960–2019 were included in this study. Using a Kendall tau statistical test ( p -value less than or equal to 0.10), 573 streamgages had a statistically significant upward trend in mean annual streamflow and are primarily located in the Midwest and northeastern United States. Of the streamgages, 182 had a statistically significant downward trend and are located primarily in the western and southeastern States. Several sites had increases in streamflow between 100 and 500 percent. Most of the streamgages with the highest increases in mean annual streamflow are along the same rivers in the eastern Dakotas, regardless of basin size. A comparison of mean annual streamflow of the last decade (2010–19) to the first decade (1960–69) of the study period shows that the largest increases in annual streamflow volumes in the United States also are in the eastern Dakotas. Among all 1,853 streamgages in the United States, the Sheyenne River near Warwick, North Dakota (U.S. Geological Survey station 05056000), has the greatest percent change, with an increase of 486 percent. Several factors may be contributing to increasing trends in streamflow in the eastern Dakotas and may include, in part, precipitation changes owing to climatic variation within the region, geologic makeup of the subsurface, and land-use changes. A better understanding of these research areas will help producers, resource managers, and infrastructure engineers to make more informed environmental and economic decisions.

North Dakota, South Dakota↗

Simulation experiments comparing nonstationary design-flood adjustments based on observed annual peak flows in the conterminous United States

While nonstationary flood frequency analysis (NSFFA) methods have proliferated, few studies have rigorously compared them for modeling changes in both the central tendency and variability of annual peak-flow series, also known as the annual maximum series (AMS), in hydrologically diverse areas. Through Monte Carlo experiments, we appraise five methods for updating estimates of 10- and 100-year floods at gauged sites using synthetic records based on sample moments and change trajectories of observed AMS in the conterminous United States (CONUS). We compare two methods that consider changes in both central tendency and variability - a Gamma generalized linear model estimated with weighted least squares and the Generalized Additive Model for Location, Scale, Shape (GAMLSS) - with a distribution-free approach (quantile regression), and baseline cases assuming stationarity or only changes in central tendency. ‘Trend-space’ plots identify realistic AMS changes for which modeling trends in both central tendency and variability were warranted based on fractional root mean squared errors (fRMSE). They also reveal statistical properties of AMS under which NSFFA models perform especially well or poorly. For instance, quantile regression performed especially well (poorly) under strong negative (positive) skewness. Although the nonstationary LP3 distribution accommodates most AMS with trends well, the sensitivity of NSFFA model performance to different sample moments and trends suggests the need for more flexibility in prescribing design-flood adjustments in the CONUS. A follow-up comparison of regional NSFFA models pooling at-site AMS would further illuminate NSFFA guidance, especially for AMS with properties less conducive to NSFFA modeling, such as positive skewness and increasing variability.

Journal of Hydrology X↗

Climate extremes as drivers of surface-water-quality trends in the United States

Surface-water quality can change in response to climate perturbations, such as changes in the frequency of heavy precipitation or droughts, through direct effects, such as dilution and concentration, and through physical processes, such as bank scour. Water quality might also change through indirect mechanisms, such as changing water demand or changes in runoff interaction with organic matter on the landscape. Many studies predict future changes in water-quality related to climate changes; however, fewer studies specifically document changes in water quality related to changes in climate, and they are usually limited in geographic scope. Recently, the U.S. Geological Survey's National Water-Quality Program reported nearly 12,000 trends in concentration and load for numerous water-quality constituents, including nutrients, sediment, major ions, and carbon. The results provide an unprecedented opportunity to examine sites across the conterminous United States for changes in water quality related to climate changes. We used published water-quality trends, modeled using the method of Weighted Regressions on Time, Season and Discharge, and calculated trends in climate extremes indices, using a modified Mann-Kendall trend method. The water-quality and the climate extremes trends were combined to identify areas in the conterminous United States where changes in climate extremes may have changed water quality. We investigated the water-quality trends in these areas to determine whether the trends related to changes in climate. We found that it was important to go beyond spatial correlation and examine trends on a watershed scale to investigate key drivers of trends. We found successful management practices in Iowa to reduce chloride concentrations, despite increases in icing days. For sediment, it appeared that management practices were having a larger effect than climate changes. For nutrients, complex forces affecting water quality make it difficult to unequivocally attribute water-quality change to climate change.

Science of the Total Environment↗

Historical and paleoflood analyses for probabilistic flood-hazard assessments—Approaches and review guidelines

Paleoflood studies are an effective means of providing specific information on the recurrence and magnitude of rare and large floods. Such information can be combined with systematic flood measurements to better assess the frequency of large floods. Paleoflood data also provide valuable information about the linkages among climate, land use, flood-hazard assessments, and channel morphology. This document summarizes methods and techniques for the preparation, gathering, evaluation, and interpretation of paleoflood information, including uncertainties, especially with respect to new statistical approaches available to efficiently use such data. We summarize best practices and strategies for assessing and mitigating uncertainties and provide guidelines on appropriate technical review of paleoflood analyses based on project goals and requirements.

Techniques and Methods↗

Flood-frequency estimation for very low annual exceedance probabilities using historical, paleoflood, and regional information with consideration of nonstationarity

Streamflow estimates for floods with an annual exceedance probability of 0.001 or lower are needed to accurately portray risks to critical infrastructure, such as nuclear powerplants and large dams. However, extrapolating flood-frequency curves developed from at-site systematic streamflow records to very low annual exceedance probabilities (less than 0.001) results in large uncertainties in the streamflow estimates. Traditionally, methods for statistically estimating flood frequency have relied on the systematic streamflow record, which provides a time series of annual maximum flood peaks, often including some historical peaks. However, most peak-flow records are less than 100 years, and uncertainties are large when trying to extrapolate magnitudes of very low annual exceedance probability events. Other data may be available that extend the record beyond the systematic dataset. Historical data are defined as data from outside the period of systematic records but within the period of human records. Examples of historical information include flood estimates from other agencies and newspaper accounts that can be translated to flood magnitude point estimates, interval estimates, or perception thresholds (such as a statement that an 1880 flood was the largest since 1869). Paleoflood data, which may also extend the dataset, include a broad range of information about flood occurrence or magnitude from sources like sediment deposits or tree rings. Several assumptions are made in flood-frequency analysis, and an understanding of whether the data conform to these assumptions is desired. A particularly difficult assumption to evaluate for flood-frequency analysis is the underlying assumption that the flood series is stationary—the assumption that a time series of peak flow varies around a constant mean within a particular range of values (constant variance). As the hydrologic community’s understanding of natural systems and anthropogenic effects on streamflows has evolved, the community has come to understand that many surface-water systems exhibit one or more forms of nonstationarity, and thus the stationarity assumption is often violated to some degree. However, there is currently (2020) no consensus among hydrologists regarding the most appropriate flood-frequency-analysis methods for nonstationary systems, and this topic remains an active area of research. A literature review was completed to summarize the state of the science of flood frequency. The literature review highlights tools available to detect nonstationarities and identifies approaches that include external information to inform flood-frequency analysis. To demonstrate methods for initial data analysis and for incorporating historical and paleoflood information in flood-frequency analysis, five sites were selected: the Red River of the North at James Avenue Pumping Station, Winnipeg, Manitoba, Canada; lower reach, Rapid Creek, South Dakota; Spring Creek, South Dakota; Cherry Creek near Melvin, Colorado; and Escalante River near Escalante, Utah. The sites were chosen for the availability of published historical and paleoflood data and for their geographic diversity and unique characteristics, which highlighted issues such as autocorrelation, change points, trends, outlier peaks, or short periods of record. An initial data analysis that involved examining records for autocorrelation, change points, and trends was completed for all sites. The flood-frequency analysis completed for this study used version 7.2 of the U.S. Geological Survey PeakFQ program. Multiple analyses were done on each site documenting the change in the flood-frequency curve when additional historical or paleoflood data were added. When other flood-frequency studies were available, their results were compared to the results here. The comparisons in some cases simply show the effect of additional years of data, whereas other comparisons show results from probability distributions or fitting methods other than those used in PeakFQ. For the Red River of the North, flood-frequency analysis shows that paleoflood data appear necessary to reasonably estimate very low annual exceedance probabilities. For the analysis of the lower reach of Rapid Creek and Spring Creek, paleoflood information helped put a high outlier from the systematic period in context; however, very low annual exceedance probabilities at these sites still had extraordinarily large confidence bounds. These sites also showed that paleoflood information might be transferred from one site to another, with the caveat that this is a case where we had existing paleoflood data to test the transfer of paleoflood information—this is not the case at many sites, and transferring paleoflood information requires assumptions about the comparability of floods at the sites. The Cherry Creek analysis affirmed the result of an earlier study that showed that the generalized Pareto distribution was not a good distribution for estimating very low annual exceedance probabilities. The Escalante River analysis showed that adding paleoflood information might increase uncertainty for very low annual exceedance probabilities, compared to analysis with the systematic period of record information only, when the paleoflood peaks are of much larger magnitudes than the systematic record.

Scientific Investigations Report↗

seawaveQ—An R package providing a model and utilities for analyzing trends in chemical concentrations in streams with a seasonal wave (seawave) and adjustment for streamflow (Q) and other ancillary variables, version 2.0.0

The seawaveQ R package provides functionality and help to fit a parametric regression model, SEAWAVE-Q, to pesticide concentration data from stream-water samples to assess trends. The model incorporates the strong seasonality and high degree of censoring common in pesticide data, and users can incorporate numerous ancillary variables such as streamflow anomalies. The model is fitted to pesticide data using maximum likelihood methods for censored data and is robust in terms of pesticide, stream location, and degree of censoring of the concentration data. This R package standardizes this methodology for trend analysis, documents the code, and provides help and tutorial information. In previous investigations, the SEAWAVE-Q model assumed a linear trend across the period analyzed. For short trend periods, this assumption of a linear trend is adequate. However, as the period of record analyzed becomes longer, the assumption of linearity is problematic because of changes in pesticide regulation and use, some of which can be abrupt. In this update to the model, a restricted cubic spline option was added for long trend periods. This option allows for more flexibility in the time component of the model. Bootstrap functionality is included to determine statistical significance. Model results with the new restricted cubic spline option are compared to the linear trend option for two pesticide-site combinations.

Open-File Report↗

Statistical methods in water resources

This text began as a collection of class notes for a course on applied statistical methods for hydrologists taught at the U.S. Geological Survey (USGS) National Training Center. Course material was formalized and organized into a textbook, first published in 1992 by Elsevier as part of their Studies in Environmental Science series. In 2002, the work was made available online as a USGS report. The text has now been updated as a USGS Techniques and Methods Report. It is intended to be a text in applied statistics for hydrology, environmental science, environmental engineering, geology, or biology that addresses distinctive features of environmental data. For example, water resources data tend to have many variables with a lower bound of zero, tend to be more skewed than data from many other disciplines, commonly contain censored data (less than values), and assumptions that the data are normally distributed are not appropriate. Computer-intensive methods (bootstrapping and permutation tests) now improve upon and replace the dependence on t-intervals, t-tests, and analysis of variance. A new chapter on sampling design addresses questions such as “How many observations do I need?” The chapter also presents distribution-free methods to help plan sampling efforts. The trends chapter has been updated to include the WRTDS (Weighted Regressions on Time, Discharge, and Season) method for analysis of water-quality data. This new version contains updated graphics and updated guidance on the use of statistical techniques. The text utilizes R, a programming language and open-source software environment, for all exercises and most graphics, and the R code used to generate figures and examples is provided for download.

Techniques and Methods↗