SEARCH · Geology Reports

Results for “Nature Scientific Data”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,153 records · Page 64Linked to original sources

USGS2012

Lateral and vertical channel movement and potential for bed-material movement on the Madison River downstream from Earthquake Lake, Montana

The 1959 Hebgen Lake earthquake caused a massive landslide (Madison Slide) that dammed the Madison River and formed Earthquake Lake. The U.S. Army Corps of Engineers excavated a spillway through the Madison Slide to permit outflow from Earthquake Lake. In June 1970, high streamflows on the Madison River severely eroded the spillway channel and damaged the roadway embankment along U.S. Highway 287 downstream from the Madison Slide. Investigations undertaken following the 1970 flood events concluded that substantial erosion through and downstream from the spillway could be expected for streamflows greater than 3,500 cubic feet per second (ft 3 /s). Accordingly, the owners of Hebgen Dam, upstream from Earthquake Lake, have tried to manage releases from Hebgen Lake to prevent streamflows from exceeding 3,500 ft 3 /s measured at the U.S. Geological Survey (USGS) gaging station 0638800 Madison River at Kirby Ranch, near Cameron, Montana. Management of flow releases from Hebgen Lake to avoid exceeding the threshold streamflow at USGS gaging station 06038800 is difficult, and has been questioned for two reasons. First, no road damage was reported downstream from the Earthquake Lake outlet in 1993, 1996, and 1997 when streamflows exceeded the 3,500-ft 3 /s threshold. Second, the 3,500-ft 3 /s threshold generally precludes releases of higher flows that could be beneficial to the blue-ribbon trout fishery downstream in the Madison River. In response to concerns about minimizing streamflow downstream from Earthquake Lake and the possible armoring of the spillway, the USGS, in cooperation with the Madison River Fisheries Technical Advisory Committee (MADTAC; Bureau of Land Management; Montana Department of Environmental Quality; Montana Fish, Wildlife and Parks; PPL-Montana; U.S. Department of Agriculture Forest Service - Gallatin National Forest; and U.S. Fish and Wildlife Service), conducted a study to determine movement of the Madison River channel downstream from Earthquake Lake and to investigate the potential for bed material movement along the same reach. The purpose of this report is to present information about the lateral and vertical movement of the Madison River from 1970 to 2006 for a 1-mile reach downstream from Earthquake Lake and for Raynolds Pass Bridge, and to provide an analysis of the potential for bed-material movement so that MADTAC can evaluate the applicability of the previously determined threshold streamflow for initiation of damaging erosion. As part of this study channel cross sections originally surveyed by the USGS in 1971 were resurveyed in 2006. Incremental channel-movement distances were determined by comparing the stream centerlines from 14 aerial photographs taken between 1970 and 2006. Depths of channel incision and aggregation were determined by comparing the 2006 and 1971 cross-section and water-surface data. Particle sizes of bed and bank materials were measured in 2006 and 2008 using the pebble-count method and sieve analyses. A one-dimensional hydraulic-flow model (HEC-RAS) was used to calculate mean boundary-shear stresses for various streamflows; these calculated boundary-shear stresses were compared to calculated critical-shear stresses for the bed materials to determine the potential for bed-material movement. A comparison of lateral channel movement distances with annual peak streamflows shows that streamflows higher than the 3,500-ft 3 /s threshold were followed by lateral channel movement except from 1991 to 1992 and possibly from 1996 to 1997. However, it was not possible to discern whether the channel moved gradually or suddenly, or in response to one peak flow, to several peak flows, or to sustained flows. The channel moved between 2002 and 2005 even when streamflows were less than the threshold streamflow of 3,500 ft 3 /s. Comparisons of cross sections and aerial photographs show that the channel has moved laterally and incised and aggraded to varying degrees. The channel has developed meander bends and has incised as much as 5–12 feet (ft) through the upstream part of the Madison Slide (cross sections 1400–800). Near cross section 800, the stream has eroded into the steep right bank between the stream and the road where fill was mechanically placed after 1970. Channel movement also was noted downstream from the Madison Slide. Near Raynolds Pass Bridge, about 3 miles (mi) downstream from Earthquake Lake, elevations across the channel have changed by -1.4 ft to +1.9 ft, but these changes were local in nature and could represent a few rocks or depressions in the bed. Overall, it does not appear that the materials eroded from the Madison Slide are causing aggradation in the subreach near the Raynolds Pass Bridge. Comparisons of critical shear stresses to mean boundary-shear stresses indicate that the D50 particle sizes (median size) along the right side of the bed between cross sections 400 and 500 and along the right side of the bed between cross sections 1300 and 1400 could move at the threshold streamflow. In contrast, most of the D84 particle sizes at those two locations probably will not move at the threshold streamflow. This lack of movement for the larger particles at the threshold streamflow could lead to further armoring of the bed as the D50 and smaller-sized particles are removed from the bed and transported downstream. The Shields parameter values from 0.04 to 0.08 that were used to calculate critical shear stresses could be conservative for a high-gradient stream such as the Madison. A higher, less conservative, Shields parameter would result in higher critical shear stresses, meaning that higher streamflows would be required to move material than those reported herein. In addition, because materials in the channel thalweg are exposed to higher boundary-shear stresses than the materials along the sides of the channel, larger, more erosion-resistant materials likely exist in the deeper parts of the channel where high-flow depths and velocities prevented sediment sampling. Movement of these materials might require higher critical shear stresses than estimated in this report. Characterization of sediment sizes in the center of the stream and observation of bed-material movement for a range of streamflows could provide information to help refine the Shields parameter and critical-shear stress estimates for bed materials in the Madison River downstream from Earthquake Lake. Furthermore, resurveying cross sections and water-surface elevations more frequently (either annually or after high streamflows) could better define the relation between streamflow and lateral and vertical channel movement.

Montana↗

USGS2016

Sediment loads and transport at constructed chutes along the Missouri River - Upper Hamburg Chute near Nebraska City, Nebraska, and Kansas Chute near Peru, Nebraska

The U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, monitored suspended sediment within constructed Missouri River chutes during March through October 2012. Chutes were constructed at selected river bends by the U.S. Army Corps of Engineers to help mitigate aquatic habitat lost through the creation and maintenance of the navigation channel on the Missouri River. The restoration and development of chutes is one method for creating shallow-water habitat within the Missouri River to meet requirements established by the amended 2000 Biological Opinion. Understanding geomorphic channel-evolution processes and sediment transport is important for the design of chutes, monitoring and maintenance of existing chutes, and characterizing the habitat that the chutes provide. This report describes the methods used to monitor suspended sediment at two Missouri River chutes and presents the results of the data analysis to help understand the suspended-sediment characteristics of each chute and the effect the chutes have on the Missouri River. Upper Hamburg chute, near Nebraska City, Nebraska, and Kansas chute, near Peru, Nebraska, were selected for monitoring. At each study site, monthly discrete samples were collected from April through October in the Missouri River main-channel transects upstream from the chute inlet, downstream from the chute outlet, at the outlet (downstream transect) of both chutes, and at the inlet (upstream transect) of Kansas chute. In addition, grab samples from all chute sampling locations were collected using autosamplers. Suspended-sediment concentration (SSC) and grain-size metrics were determined for all samples (discrete and grab). Continuous water-quality monitors recorded turbidity and water temperature at 15-minute intervals at the three chute sampling locations. Two acoustic Doppler velocimeters, one within each chute, measured water depth and current velocities continuously. The depth and velocity data were used to estimate streamflow within each chute. The sampling design was developed to understand the suspended-sediment differences within each chute and between the chute and the Missouri River main channel during discrete sampling. The sampling design also allowed for site-specific surrogate relations between SSC and turbidity to be developed, which could be used to compute real-time estimates of SSC and sediment loads within the chutes. Real-time estimates of SSC and sediment loads enable a better understanding of sediment transport within the chutes during times when physical samples are not collected, including periods of high flow. High flows during the summer of 2011 resulted in substantial alterations to both studied chutes; therefore, the U.S. Army Corps of Engineers repaired and modified both chutes during 2012. These unforeseen repairs and modifications within the chutes added uncertainty to the analysis because concentrations were altered by construction equipment and flow alteration. Daily suspended-sediment and suspended-silt loads were estimated based on surrogate relations with turbidity. A linear regression was used to estimate equal-width increment (EWI)-equivalent SSC from autosampler SSC before using the model-calibration dataset to determine the best-fit model for prediction of SSC from the turbidity and, in some cases, discharge. Correlation between suspended-sand concentration ( SSandC ) in EWI samples and concurrent samples collected by an autosampler was low; therefore, SSandC was excluded from development of surrogate relations because a large part of the calibration dataset was from autosamples. Instead, SSandC was estimated as SSC minus suspended-silt-clay concentration ( SSiltC ). At all sites, the best-fit models included the base-10 logarithm of concentration and turbidity, and at Kansas chute upstream, the base-10 logarithm of streamflow was also included in the best-fit models. These surrogate models were used to estimate continuous time series of SSC and SSiltC . Estimated concentrations of suspended sediment were used to estimate instantaneous and daily loads for total suspended sediment, suspended silt-clay, and suspended sand. Estimated daily suspended-sediment loads were not significantly different between upstream and downstream transects within the Kansas chute, and most individual daily loads within the chute were not significantly different between upstream and downstream transects when evaluated using overlap in daily 95-percent confidence intervals. The comparison of daily load values for upstream and downstream chute transects, as estimated from turbidity-based surrogate models for Kansas chute, documents the daily dynamic nature of sediment transport within the chute with a temporal resolution that is not practical with discrete suspended-sediment sampling alone. Comparisons of concentrations and loads from EWI samples collected from different transects within a study site resulted in few significant differences, but comparisons are limited by small sample sizes and large within-transect variability. When comparing the Missouri River upstream transect to the chute inlet transect, similar results were determined in 2012 as were determined in 2008—the chute inlet affected the amount of sediment entering the chute from the main channel. In addition, the Kansas chute is potentially affecting the sediment concentration within the Missouri River main channel, but small sample size and construction activities within the chute limit the ability to fully understand either the effect of the chute in 2012 or the effect of the chute on the main channel during a year without construction. Finally, some differences in SSC were detected between the Missouri River upstream transects and the chute downstream transects; however, the effect of the chutes on the Missouri River main-channel sediment transport was difficult to isolate because of construction activities and sampling variability.

Nebraska↗

USGS2006

Simulation of nutrient and sediment concentrations and loads in the Delaware inland bays watershed: Extension of the hydrologic and water-quality model to ungaged segments

Rapid population increases, agriculture, and industrial practices have been identified as important sources of excessive nutrients and sediments in the Delaware Inland Bays watershed. The amount and effect of excessive nutrients and sediments in the Inland Bays watershed have been well documented by the Delaware Geological Survey, the Delaware Department of Natural Resources and Environmental Control, the U.S. Environmental Protection Agency’s National Estuary Program, the Delaware Center for Inland Bays, the University of Delaware, and other agencies. This documentation and data previously were used to develop a hydrologic and water-quality model of the Delaware Inland Bays watershed to simulate nutrients and sediment concentrations and loads, and to calibrate the model by comparing concentrations and streamflow data at six stations in the watershed over a limited period of time (October 1998 through April 2000). Although the model predictions of nutrient and sediment concentrations for the calibrated segments were fairly accurate, the predictions for the 28 ungaged segments located near tidal areas, where stream data were not available, were above the range of values measured in the area. The cooperative study established in 2000 by the Delaware Department of Natural Resources and Environmental Control, the Delaware Geological Survey, and the U.S. Geological Survey was extended to evaluate the model predictions in ungaged segments and to ensure that the model, developed as a planning and management tool, could accurately predict nutrient and sediment concentrations within the measured range of values in the area. The evaluation of the predictions was limited to the period of calibration (1999) of the 2003 model. To develop estimates on ungaged watersheds, parameter values from calibrated segments are transferred to the ungaged segments; however, accurate predictions are unlikely where parameter transference is subject to error. The unexpected nutrient and sediment concentrations simulated with the 2003 model were likely the result of inappropriate criteria for the transference of parameter values. From a model-simulation perspective, it is a common practice to transfer parameter values based on the similarity of soils or the similarity of land-use proportions between segments. For the Inland Bays model, the similarity of soils between segments was used as the basis to transfer parameter values. An alternative approach, which is documented in this report, is based on the similarity of the spatial distribution of the land use between segments and the similarity of land-use proportions, as these can be important factors for the transference of parameter values in lumped models. Previous work determined that the difference in the variation of runoff due to various spatial distributions of land use within a watershed can cause substantialloss of accuracy in the model predictions. The incorporation of the spatial distribution of land use to transfer parameter values from calibrated to uncalibrated segments provided more consistent and rational predictions of flow, especially during the summer, and consequently, predictions of lower nutrient concentrations during the same period. For the segments where the similarity of spatial distribution of land use was not clearly established with a calibrated segment, the similarity of the location of the most impervious areas was also used as a criterion for the transference of parameter values. The model predictions from the 28 ungaged segments were verified through comparison with measured in-stream concentrations from local and nearby streams provided by the Delaware Department of Natural Resources and Environmental Control. Model results indicated that the predicted edge-of-stream total suspended solids loads in the Inland Bays watershed were low in comparison to loads reported for the Eastern Shore of Maryland from the Chesapeake Bay watershed model. The flatness of the terrain and the low annual surface runoff are important factors in determining the amount of detached sediment from the land that is delivered to streams. The highest predicted total suspended solids loads were found in the southern part of the watershed, where the values are associated with high total streamflow and a high surface-runoff component, and related to soil and aquifer permeability and land use. Nutrient loads from model segments in the southern part of the Inland Bays watershed were also higher than those measured in the northern part of the basin, due to relatively high runoff and the substantial amount of available organic fertilizer (animal waste) that results in over-application of organic fertilizer to crops. Time series of simulated hourly concentrations indicated a seasonal pattern in the simulated base flow for total nitrogen, with the lowest values occurring during the summer and the highest values during the winter months. Total phosphorus and total-suspended-solids concentrations were less seasonal and were more storm-dependent; in general, base-flow concentrations of total phosphorus and total suspended solids were low. During storm events, the total nitrogen concentrations tended to be diluted and total phosphorus concentrations tended to rise sharply. Nitrogen was transported mainly in the aqueous phase and largely through ground water, whereas phosphorus was strongly associated with sediment, which washes off during rainfall events.

Delaware↗

USGS2007

Status of Water Levels and Selected Water-Quality Conditions in the Sparta-Memphis Aquifer in Arkansas and the Status of Water Levels in the Sparta Aquifer in Louisiana, Spring 2005

The U.S. Geological Survey in cooperation with the Arkansas Natural Resources Commission, the Arkansas Geological Commission, and the Louisiana Department of Transportation and Development has monitored water levels in the Sparta Sand of Claiborne Group and Memphis Sand of Claiborne Group since the 1920's. Ground-water withdrawals have increased while water levels have declined since monitoring was initiated. This report has been produced to describe ground-water levels in the aquifers in the Sparta Sand and Memphis Sand and provide information for the management of this valuable resource. The 2005 potentiometric-surface map of the aquifers in the Sparta Sand and Memphis Sand was constructed using water-level data collected in 333 wells in Arkansas and 120 wells in Louisiana during the spring of 2005. The highest water-level altitude measured in Arkansas was 327 feet above National Geodetic Vertical Datum of 1929 located in Grant County in the outcrop at the western boundary of the study area; the lowest water-level altitude was 189 feet below National Geodetic Vertical Datum of 1929 in Union County. The highest water-level altitude measured in Louisiana was 246 feet above National Geodetic Vertical Datum of 1929 located in Bossier Parish in the outcrop area near the western boundary of the study area; the lowest water-level altitude was 226 feet below National Geodetic Vertical Datum of 1929 in central Ouachita Parish. Three large depressions centered in Columbia, Jefferson, and Union Counties in Arkansas are the result of large withdrawals for industrial and public supplies. In Louisiana, three major pumping centers are in Ouachita, Jackson, and Lincoln Parishes. Water withdrawals from these major pumping centers primarily is used for industrial and public-supply purposes. Withdrawals from Ouachita and Lincoln Parishes and Union County, Arkansas, primarily for industrial purposes, have caused the resulting cones of depression to coalesce so that the -40 foot potentiometric contour encircles the three pumping centers. Seven smaller depressions are evident on the 2005 Sparta-Memphis potentiometric-surface map located in Webster and Winn Parishes, Louisiana, and Calhoun, Cleveland, western Columbia, Desha, and Lafayette Counties, Arkansas. The depression in Calhoun County initially was shown in the 1996-1997 potentiometric surface. The depression in Desha County initially was shown in the 1999 potentiometric surface. The depressions in Webster and Winn Parishes were shown as early as 1975. The depressions in Cleveland, western Columbia, and Lafayette Counties initially were shown in the 2003 potentiometric surface. A map of differences in water-level measurements between 2001 and 2005 was constructed using the difference between water-level measurements from 294 wells in Arkansas and 29 wells in Louisiana. The difference in water levels between 2001 and 2005 ranged from -30.1 to 44.6 feet. The largest rise of 44.6 feet in water level measured was in Union County in Arkansas. The largest decline of 30.1 feet in water level measured was in Columbia County in Arkansas. Areas with a general rise in water levels in Arkansas are shown in Arkansas, Columbia, Craighead, Jefferson, Prairie, and the western half of Union Counties. The area around west-central Union County had rises as much as 44.6 feet, with seven wells showing a rise of 20 feet or greater, which is an annual rise of 5 feet or greater. Areas in Arkansas with a general decline in water level are shown in western Bradley, eastern Calhoun, Cleveland, Cross, Desha, Drew, Lafayette, Lee, Lincoln, Lonoke, Poinsett, and the eastern half of Union Counties. In Louisiana, the water-level difference map showed a general rise in water levels in northern Claiborne, northern Webster, and northwestern Union Parishes mainly because of a decrease in industrial withdrawals in southern Arkansas, particularly Union County. Another rise in water level was indicated in western

Scientific Investigations Report↗

USGS2024

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

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

South Dakota↗

USGS2018

Land-cover changes associated with oil and natural-gas production and concentrations of selected constituents in surface-water and streambed-sediment samples collected upstream from and within an area of oil and natural-gas production, south Texas, 2008–17

The extensive development of oil and natural-gas resources in south Texas during the past 10 years has led to questions regarding possible environmental effects of processes associated with oil and natural-gas production, in particular the process of hydraulic fracturing, on water and other natural resources. Part of the lower San Antonio River watershed intersects an area of oil and natural-gas production from the sedimentary rocks that compose the Eagle Ford Group. The rapid expansion of infrastructure associated with oil and natural-gas production increases potential pathways for inorganic and organic contaminants to enter surface-water systems. The U.S. Geological Survey, in cooperation with the San Antonio River Authority, analyzed geospatial data from different years (2008 and 2015) to evaluate changes in land cover associated with oil and natural-gas production activities in the lower San Antonio River watershed. Impervious surface in this study is defined as land cover consisting of well pads, oil- and gas-related features, or roads. The areal coverage associated with impervious surface increased from 201 acres to 5,390 acres (net increase of 5,189 acres) between 2008 and 2015. The total percentage of the study area accounted for by impervious surface resulting from oil and natural-gas production activities increased from 0.034 percent to 0.912 percent, which is an increase of approximately 27-fold. Collectively, 0.878 percent of the study area was converted to new impervious surface between 2008 and 2015. If the area associated with new storage ponds (0.066 percent) is added to the estimate of total land-cover changes as a result of oil and natural-gas production, then 0.944 percent of the study area was altered. During 2015–17, surface-water samples collected from 5 sites and streambed-sediment samples collected from 17 sites in the lower San Antonio River watershed were analyzed for a broad range of constituents that might be associated with oil and natural-gas production. All major elements, trace elements, semivolatile organic compounds (SVOCs), and volatile organic compounds (VOCs) measured in surface-water samples were detected at concentrations less than any of the U.S. Environmental Protection Agency’s water-quality standards. In general, the greatest SVOC and VOC concentrations were observed in samples collected from sites upstream from the area of active oil and natural-gas production and just downstream from urban areas. The lack of benzene, toluene, ethylbenzene, and all isomers of xylene (hereinafter referred to as BTEX) for most sites within the area of active oil and natural-gas production indicates that little, if any, local runoff associated with the area of active oil and natural-gas production has contaminated the surface water with BTEX compounds. Glycols, which are commonly used in hydraulic fracturing fluids as scale inhibitors, were detected in one surface-water sample from Ecleto Creek within the area of oil and natural-gas production; however, the presence of glycols does not necessarily indicate contamination from hydraulic fracturing fluid. The glycols detected also have other potential sources including the use of diethylene and ethylene glycols in antifreeze used in vehicles and the use of triethylene glycol in antibacterial air sanitizers. The concentrations of select constituents in the streambed-sediment samples were compared to sediment quality guidelines (SQGs). The SQGs evaluate the potential toxicity of bed sediments to sediment-dwelling organisms. Two SQG concentration levels are used: (1) a lower level, called the threshold effect concentration (TEC), below which harmful effects to benthic biota are not expected, and (2) a higher level, the probable effect concentration (PEC), above which harmful effects are expected to occur frequently. The PEC for arsenic was exceeded in a sample collected from one site on Ecleto Creek. The origin of the elevated arsenic concentration is unknown; the contamination likely is not related to oil and natural-gas production because the site of the sample collection is located upstream from the area of active oil and natural-gas production. Streambed-sediment samples were analyzed for selected polycyclic aromatic hydrocarbons (PAHs) because PAHs can be used as indicators of petroleum hydrocarbons associated with produced waters. Each streambed-sediment sample was analyzed for two size fractions of PAHs: less than (<) 63 micrometers (μm) and < 2 millimeters (mm). Total PAH concentrations in all samples, regardless of size fraction, were less than the TEC for total PAHs of 1,610 micrograms per kilogram. Total PAH concentrations generally were greater in the <63-μm size-fraction samples than in the <2-mm size-fraction samples, indicating that PAHs could potentially sorb more readily to the exclusively silt- and clay-sized particles that compose <63-μm size-fraction samples than to the mixture of silt and clay and larger sized particles that compose the <2-mm size-fraction samples. Total PAH concentrations typically were greater in the samples collected from the sites upstream from the area of active oil and natural-gas production compared to those collected from sites within the area in both the <2-mm and <63-μm size-fraction samples. The smaller PAH concentrations measured in samples collected from within the area of active oil and natural-gas production in comparison to the upstream urbanized areas indicate relatively minor additional local contributions of PAHs of uncertain origin to the watershed.

Texas↗

USGS2014

Status and understanding of groundwater quality in the Sierra Nevada Regional study unit, 2008: California GAMA Priority Basin Project

Groundwater quality in the Sierra Nevada Regional (SNR) study unit was investigated as part of the California State Water Resources Control Board’s Groundwater Ambient Monitoring and Assessment Program Priority Basin Project. The study was designed to provide statistically unbiased assessments of the quality of untreated groundwater within the primary aquifer system of the Sierra Nevada. The primary aquifer system for the SNR study unit was delineated by the depth intervals over which wells in the State of California’s database of public drinking-water supply wells are open or screened. Two types of assessments were made: (1) a status assessment that described the current quality of the groundwater resource, and (2) an evaluation of relations between groundwater quality and potential explanatory factors that represent characteristics of the primary aquifer system. The assessments characterize untreated groundwater quality, rather than the quality of treated drinking water delivered to consumers by water distributors. The status assessment was based on water-quality data collected by the U.S. Geological Survey from 83 wells in the SNR study unit in 2008 and from 117 wells in 3 small study units within the SNR study unit in 2006–07 and on water-quality data compiled in the State’s database for 1,066 wells sampled in 2006–08. To provide some context for the results, water-quality data were converted to relative-concentrations (RCs), which are the sample concentrations divided by the concentrations of Federal or California regulatory and non-regulatory benchmarks for drinking-water quality. RCs for inorganic constituents (major ions, trace elements, nutrients, and radioactive constituents) were classified as “high” (RC > 1.0, indicating that concentration is above the benchmark), “moderate” (1.0 ≥ RC > 0.5), or “low” (RC ≤ 0.5). For organic constituents (volatile organic compounds and pesticides) and special-interest constituents (perchlorate and N -nitrosodimethylamine [NDMA]), the boundary between moderate and low RCs was set at 0.1. All benchmarks used for organic constituents were health-based, whereas health-based and aesthetic-based benchmarks were used for inorganic constituents. The primary metric used for quantifying regional-scale groundwater quality was “aquifer-scale proportion.” Aquifer-scale proportions were calculated as the areal percentages of the primary aquifer system having high, moderate, and low RCs for a given constituent or class of constituents. The SNR study unit area was classified into four aquifer lithologic types—granitic rocks, metamorphic rocks, sedimentary deposits, and volcanic rocks—and aquifer-scale proportions were calculated on an area-weighted basis for each of the four aquifer lithologies and for the study unit as a whole (aggregated system). The results of the status assessment indicated that inorganic constituents were present at high and moderate RCs in greater proportions in the SNR study unit aggregated primary aquifer system than were organic constituents and that there were significant differences (p < 0.05) between the four aquifer lithologies. One or more inorganic constituents with health-based benchmarks were present at high RCs in 16 percent of the aggregated primary aquifer system and at moderate RCs in 21 percent. Arsenic (9.7 percent), uranium (2.9 percent), boron (2.0 percent), fluoride (1.8 percent), and nitrate (1.4 percent) were the constituents most commonly present at high RCs. For inorganic constituents with aesthetic-based benchmarks, 18 percent of the aggregated primary aquifer system had high RCs of one or more constituent, and 6.8 percent had moderate RCs. Iron (15.8 percent), manganese (15.1 percent), and total dissolved solids (1.3 percent) were the constituents most commonly present at high RCs. Organic constituents were not detected in 72 percent of the primary aquifer system. One or more organic constituents had high RCs in 0.1 percent of the primary aquifer system, moderate RCs in 3.0 percent, and low RCs in 25 percent. Proportions of the four lithologic primary aquifer systems with high or moderate concentrations of organic constituents were not significantly different. Three organic constituents had area-weighted detection frequencies greater than 10 percent in the primary aquifer system as a whole or at least one of the four lithologic primary aquifer systems: the gasoline oxygenate methyl tert -butyl ether, the trihalomethane chloroform, and the herbicide simazine. The special-interest constituent perchlorate was detected at high RCs in 0.01 percent of the primary aquifer system and at moderate RCs in 1.0 percent, and detection frequencies could be accounted for by the distribution of perchlorate under natural conditions. Statistical tests were used to evaluate relations between constituent concentrations and potential explanatory factors descriptive of land use, geography, depth, geochemical conditions, and groundwater age. Higher concentrations of trace elements, radioactive constituents, and constituents with aesthetic-based benchmarks generally were associated with anoxic conditions, higher pH, and location within a particular compositional band in the Sierra Nevada batholith corresponding to the southwestern part of the study unit. High concentrations of organic constituents generally were associated with greater proportions of urban land use. No significant relations were observed between the concentrations of organic constituents and measures of well depth or groundwater age, perhaps because of the high proportions of springs and modern groundwater in the dataset.

California↗

USGS2014

Methods for estimating magnitude and frequency of floods in Arizona, developed with unregulated and rural peak-flow data through water year 2010

Flooding is among the worst natural disasters responsible for loss of life and property in Arizona, underscoring the importance of accurate estimation of flood magnitude for proper structural design and floodplain mapping. Twenty-four years of additional peak-flow data have been recorded since the last comprehensive regional flood frequency analysis conducted in Arizona. Periodically, flood frequency estimates and regional regression equations must be revised to maintain the accurate estimation of flood frequency and magnitude. Annual peak-flow data collected through water year 2010 were compiled from 448 unregulated streamflow-gaging stations, hereafter referred to as streamgages, in Arizona having a minimum of 10 years of record. Flood frequency estimates were first computed with station (or at-site) skew using the Expected Moments Algorithm with a multiple Grubbs-Beck test to identify multiple potentially influential low flows to fit a Pearson Type III distribution. Next, a multiple step Bayesian least-squares-regression approach was used to determine a new statewide regional skew of −0.09. No basin characteristics analyzed were statistically significant in explaining the variation in skew and as a result, the constant model was chosen as the best regional skew model for the Arizona study area. The mean square error used in Bulletin 17B (B17B) of the Interagency Advisory Committee on Water Data is used to describe the precision of the regional skew. The constant model had a mean square error equal to 0.08, which corresponds to an effective record length of 85 years. This is a marked improvement over a previous Arizona regional skew analysis, with a reported mean square error of 0.31, for a corresponding effective record length of around 17 years. Thus the new regional model had almost five times the information content (as measured by effective record length) of that calculated in USGS Water Supply Paper 2433, published in 1997, or the value of 0.302 reported in the B17B generalized skew map. The flood frequency estimates were recalculated using a weighted skew of the station and regional skew. Station flood frequency estimates for each streamgage are presented for the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities. Geographical information systems were used to compute basin characteristic information for each streamgage for the purpose of developing regional equations to estimate flood statistics at ungaged basins. Five hydrologic flood regions in Arizona were defined in a multivariate regionalization process based on mean basin elevation, mean annual precipitation, and soil permeability. A regional generalized least-squares-regression analysis was used to develop five sets of equations from 344 nonredundant streamgages, corresponding to five regions, for estimating the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probabilities at ungaged basins in Arizona. The regression equations developed for these five regions were based on one or more of the statistically significant explanatory variables: drainage area, mean basin elevation, and mean annual precipitation. Average standard errors of prediction for the regression regions for the five regions ranged from 27 to 122 percent and the pseudo-coefficients of determination (pseudo-R 2 ), a measure of the proportion of peak-flow variation that is explained by the basin characteristics, ranged from 68 to 98 percent. Regression equations for Central Highlands (region 4) had the lowest model error and the greatest pseudo-R 2 metrics. The equations for Colorado Plateau (region 2) regression equations generally had greater model error and lower pseudo-R 2 metrics. The improvement of regional regression equation model error and pseudo-R 2 metrics was related to higher numbers of streamgages, longer period of record, and even spatial coverage within a region. The regional regression equations were integrated into the U.S. Geological Survey’s StreamStats program. The StreamStats program is a national map-based web application that allows the public to easily access published flood frequency and basin characteristic statistics. The interactive web application allows a user to select a point within a watershed (gaged or ungaged) and retrieve flood-frequency estimates derived from the current regional regression equations and geographic information system data within the selected basin. StreamStats provides users with an efficient and accurate means for retrieving the most up to date flood frequency and basin characteristic data. StreamStats is intended to provide consistent statistics, minimize user error, and reduce the need for large datasets and costly geographic information system software.

Arizona↗

USGS2021

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↗

USGS2007

Benthic habitats and offshore geological resources of Kaloko-Honokohau National Historical Park, Hawai‘i

Kaloko-Honokōhau National Historical Park (KAHO) is one of three National Park lands located along the western coast of the Island of Hawai‘i and the only one to include submerged lands and marine resources within its official boundaries. The park was established in 1978 and is 1,160 acres in size, including 596 acres of marine area. The submerged lands are currently managed by the State of Hawaii, Department of Land and Natural Resources, Division of Aquatic Resources (DLNR-DAR). Marine resources located within KAHO include coral reef and habitat for many marine animals, including the green sea turtle and a variety of fish and invertebrates. In addition, many archeological, cultural, and recreational resources are located within the marine realm of the park. Potential threats and stressors to the modern marine environment include ground-water and surface-water contamination, invasive plants and algae, fishing pressure and use of monofilament gill nets (which can ensnare marine life or become tangled on reefs and be left behind as fishing debris), and visitor use impacts, such as scuba diving and snorkeling. Illegal dumping, oil releases, boat groundings, and other physical damage to reef resources are potential threats from users of the nearby harbor. Special issues of concern for the park include establishing baseline conditions of the offshore resources before the development of adjacent coastal lands. Until this study, only a general knowledge of the distribution of benthic habitats and the characteristics of the offshore region of Kaloko-Honokōhau National Historical Park was available. In 2003, a collaborative project between the U.S. Geological Survey (USGS) Coastal and Marine Geology Program and the National Park Service (NPS) was initiated to develop detailed benthic-habitat classification maps for the marine lands within and adjacent to the park. The intent of this project is to provide baseline maps and a Geographic Information System (GIS) database and description of the biological and geological resources of these marine lands in order to facilitate the management, interpretation, and understanding of park resources. A benthic-habitat classification map was created for the park using existing color aerial photography, Scanning Hydrographic Operational Airborne Lidar Survey (SHOALS) bathymetric data, georeferenced underwater video, and still photography. Individual habitat polygons were classified using five basic attributes: (1) major structure or substrate, (2) dominant structure, (3) major biologic cover on the substrate, (4) percentage of major biological cover, and (5) geographic zone. Additional information regarding geology, morphology, and coral species were also noted. Marine resources located within KAHO include coral reef and habitat for many marine animals, including the green sea turtle and a variety of fish and invertebrates. In addition, many archeological, cultural, and recreational resources are located within the marine realm of the park. Potential threats and stressors to the modern marine environment include ground-water and surface-water contamination, invasive plants and algae, fishing pressure and use of monofilament gill nets (which can ensnare marine life or become tangled on reefs and be left behind as fishing debris), and visitor use impacts, such as scuba diving and snorkeling. Illegal dumping, oil releases, boat groundings, and other physical damage to reef resources are potential threats from users of the nearby harbor. Special issues of concern for the park include establishing baseline conditions of the offshore resources before the development of adjacent coastal lands. Until this study, only a general knowledge of the distribution of benthic habitats and the characteristics of the offshore region of Kaloko-Honokōhau National Historical Park was available. In 2003, a collaborative project between the U.S. Geological Survey (USGS) Coastal and Marine Geology Program and the National Park Service (NPS) was initiated to develop detailed benthic-habitat classification maps for the marine lands within and adjacent to the park. The intent of this project is to provide baseline maps and a Geographic Information System (GIS) database and description of the biological and geological resources of these marine lands in order to facilitate the management, interpretation, and understanding of park resources. A benthic-habitat classification map was created for the park using existing color aerial photography, Scanning Hydrographic Operational Airborne Lidar Survey (SHOALS) bathymetric data, georeferenced underwater video, and still photography. Individual habitat polygons were classified using five basic attributes: (1) major structure or substrate, (2) dominant structure, (3) major biologic cover on the substrate, (4) percentage of major biological cover, and (5) geographic zone. Additional information regarding geology, morphology, and coral species were also noted.

Hawai'i↗

USGS2023

Evaluation of natural and anthropogenic (human-made) hexavalent chromium

Hexavalent chromium, Cr(VI), was released between 1952 and 1964 from the Pacific Gas and Electric Company (PG&E) Hinkley compressor station, in the Mojave Desert about 80 miles northeast of Los Angeles, California. Geologic, geochemical, and hydrologic data from more than 100 wells collected between March 2015 and November 2017 were interpreted using a summative-scale analysis to define the extent of anthropogenic (human-made) Cr(VI) in groundwater. The summative scale consisted of eight questions requiring binary (yes or no) answers for each sampled well. The questions were intended to (1) provide a transparent framework for data interpretation in which all stakeholders participated; (2) provide unbiased interpretation of data traceable to numerical measurements; (3) provide a framework that enabled geologic, geochemical, and hydrologic data to be considered collectively; and (4) consolidate different types of data into a simple, easy-to-understand interpretation. When data from each well are scored using questions and metrics within the summative scale, all stakeholders would score each well the same way and would draw the same summative-scale Cr(VI) plume extent. The areal extent of the summative-scale Cr(VI) plume was 5.5 square miles (mi 2 ); this is larger than the 2.2-mi 2 extent of the October–December 2015 (Q4 2015) regulatory Cr(VI) plume but smaller than the 8.3-mi2 maximum mapped extent of Cr(VI) greater than the interim regulatory Cr(VI) background value of 3.1 micrograms per liter (μg/L). The summative-scale Cr(VI) plume is within the area covered by the PG&E monitoring well network and lies within “Mojave-type” deposits composed of low-chromium stream and near-shore lake deposits sourced from the Mojave River. The summative-scale Cr(VI) plume included all shallow wells within the footprint of the Q4 2015 regulatory Cr(VI) plume, but summative-scale scores indicate that anthropogenic Cr(VI) was not present in several wells within the footprint of the regulatory Cr(VI) plume that were screened within the deep zone of the upper aquifer. The summative-scale Cr(VI) plume extent was consistent with mineralogic and geochemical data collected as part of this study that were not used within the summative-scale analysis. Data from wells outside the summative-scale Cr(VI) plume collected for regulatory purposes from April 2017 through March 2018 were used to estimate Cr(VI) background concentrations as the upper 95-percent tolerance limit (UTL 95 ) in different parts of Hinkley and Water Valleys. The UTL 95 values were calculated using the computer program ProUCL 5.1 and are suitable for use by regulatory agencies in support of (1) updating the regulatory Cr(VI) plume extent and management of Cr(VI) near the plume margins, (2) establishing cleanup goals for Cr(VI) within the updated regulatory Cr(VI) plume, and (3) identifying unusual Cr(VI) concentrations outside the regulatory Cr(VI) plume. The nonparametric UTL 95 values for wells screened in Mojave-type deposits in the eastern, western, and northern subareas of Hinkley Valley were 3.7, 3.9, and 4.0 μg/L, respectively. The normal UTL 95 values for wells screened in undifferentiated, unconsolidated deposits in the eastern and western subareas and the northern subarea upgradient from the Mount General fault were 2.8, 3.8, and 4.8 μg/L, respectively. An overall normal UTL 95 value of 3.8 μg/L was calculated for undifferentiated, unconsolidated deposits in these areas. This value is similar to the overall nonparametric UTL 95 value of 3.9 μg/L calculated for Mojave-type deposits and similar to the maximum Cr(VI) concentration of older groundwater in contact with Mojave-type deposits of 3.6 μg/L. The provenance of most PG&E monitoring wells is not precisely known, and the UTL 95 values for wells screened in undifferentiated, unconsolidated deposits in the different subareas may be more widely applicable for regulatory purposes than the UTL 95 values for Mojave-type deposits. The UTL 95 value of 2.8 μg/L for wells screened in undifferentiated, unconsolidated deposits in the eastern subarea is important for plume management because most of the summative-scale Cr(VI) plume is within the eastern subarea. A UTL 95 value of 5.8 μg/L was calculated for older (pre-1952) groundwater associated with mudflat/playa deposits in the eastern subarea near Mount General. A UTL 95 value of 2.3 μg/L was calculated for Mojave-type deposits within the Cr(VI) plume downgradient from the Hinkley compressor station after regulatory updates. This lower value is consistent with neutral to slightly alkaline, younger (post-1952) groundwater within coarse-textured, low-chromium Mojave-type deposits in this area and may be a suitable metric for Cr(VI) cleanup goals. The UTL 95 value of 4.8 μg/L for wells screened in undifferentiated, unconsolidated deposits in the northern subarea upgradient from the Mount General fault provides for possible increases in Cr(VI) concentrations if water levels continue to decline. Downgradient from the Q4 2015 regulatory Cr(VI) plume and the summative-scale Cr(VI) plume, UTL 95 values of 9.0 and 6.4 μg/L were calculated for wells screened in undifferentiated, unconsolidated deposits in the northern subarea downgradient from the Mount General fault and for Water Valley, respectively, consistent with different geologic and geochemical conditions in these areas. The UTL 95 values calculated as part of this study provide scientifically defensible estimates of background Cr(VI) concentrations that differ with local geologic, geochemical, and hydrologic conditions in Hinkley and Water Valleys. The regulatory Cr(VI) plume extent can be updated on the basis of these values. The summative-scale Cr(VI) plume extent may contain wells having anthropogenic Cr(VI) concentrations less than the UTL 95 values for their respective subareas that may not require regulatory attention, and an updated regulatory Cr(VI) plume extent may be less than the summative-scale Cr(VI) plume extent. The UTL 95 values are not background Cr(VI) concentrations for regulatory purposes, and the authority to establish regulatory values resides solely with the Lahontan Regional Water Quality Control Board.

California↗

USGS2020

Hydrology and geomorphology of the Taiya River near the West Creek Tributary, southeast Alaska

The Taiya River flows through the Chilkoot Trail Unit of Klondike Gold Rush National Historical Park in southeast Alaska, which was founded to preserve cultural and historical resources and further understanding of natural processes active in the surrounding coastal-to-subarctic basin. Riverine processes exert an important influence on ecologically important boreal toad ( Anaxryus boreas boreas ), salmon [chum salmon ( Oncorhynchus keta ), pink salmon ( O. gorbushca ), and coho salmon ( O. kisutch )], and eulachon ( Thaleichthys pacificus ) habitats, erosion of the historic ghost town of Dyea and other cultural and historical artifacts, and recreational opportunities in the lower 7.5 kilometers (km) of the Taiya River valley bottom. Recurrent consideration of hydroelectric development in West Creek upstream of the park since the 1980s has included proposals for damming and diverting West Creek, which could alter the delivery of water and sediment to this section of the Taiya River. To improve understanding of the hydrologic dependence of park resources for the purposes of guiding effective monitoring and conservation, this study, conducted by the U.S. Geological Survey in cooperation with the National Park Service, used a review of hydrologic data, collection of discrete suspended sediment data, geomorphic mapping, and analysis of historical aerial and ground photographs in a reconnaissance of formative geomorphic processes and hydrologic conditions in the lower 7.5 km of the Taiya River valley bottom. Streamflow and suspended sediment data collected at the U.S. Geological Survey streamgages on the Taiya River and West Creek, combined with historical data, document conditions consistent with streams draining strongly glacierized basins in Alaska. Suspended sediment concentrations from samples collected concurrently over six varying flow levels during 2017–18 ranged from 6 to 284 milligrams per liter (mg/L) for the Taiya River and 13 to 162 mg/L for West Creek, which are similar to or slightly higher than historical values. For the common period of record (1970–77), correlation of daily mean discharge between the two streams was strongest (Pearson’s r = 0.97) during the prolonged May–October high-flow season and weakest ( r = 0.90) during the November–April low-flow season, when West Creek daily mean discharge was proportionally higher. For the Taiya River, streamflow data compared between the available periods of record (1970–77 and 2004–17) showed no decadal-scale patterns in mean annual discharge but did show a shift toward an earlier spring snowmelt pulse. Notable flooding in the Taiya River Basin includes glacial lake outburst floods from the Nourse River valley prior to and during the 1897–98 Gold Rush, a 2002 glacial lake outburst flood from the West Creek valley, and a 1967 rainfall-generated flood. Geomorphic mapping identified four categories of surfaces in the valley bottom—active main stem, abandoned main stem, alluvial fans, and emergent tidal surfaces. Using the maps, main-stem surfaces were subdivided into age categories to identify channel migration patterns from prior to 1940s to 2018. The valley bottom is dominated by active or abandoned channels of the Taiya River except at the extensive low-angle West Creek fan. The active main stem presently supports a mostly single-thread channel with bars and a few sloughs, but the channel actively moved and sometimes was braided within multiple, wider unvegetated corridors in 1894 and earlier. An inventory of 29 off-main-stem channels identified for the study indicates that abandoned main stem channels provide local topographic lows that can intercept groundwater or sustain tributary flow, facilitating the formation of most nonestuarine wetlands in the valley and sustaining important boreal toad breeding habitat. Within the active main stem corridor, the channel has episodically built and reworked meanders and bars, eroding more than one-half of the historic Dyea townsite, in response to glacially controlled delivery of water and sediment, flooding, inputs from West Creek, local features including large woody debris and beaver dams, and rapid uplift from isostatic rebound. West Creek has constructed a large, persistent fan, provoked kilometer-scale Taiya River channel change near the confluence, constructively added to high-season streamflow that affects Taiya River channel migration capacity, disproportionately contributed early-season streamflow, and possibly contributed to groundwater levels in the valley bottom. The progressive narrowing and stability of the main stem corridor, possibly a result of reduction in the magnitude or frequency of glacial lake outburst floods or glacial sediment delivery to streams, indicates less active future reworking of abandoned main-stem surfaces or regeneration of wetland features. The fluvial history of the Taiya River valley bottom collectively indicates continued channel change within a limited corridor, relative stability in wetland locations but uncertainty in stability of groundwater supply to them, and channel incision and extension in response to uplift.

Alaska↗

USGS2021

Historical changes in bed elevation and water depth within the Nehalem Bay, Oregon, 1891–2019

Estuaries, at the nexus of rivers and the ocean, are depositional areas that respond to changes in streamflow, tides, sea level, and inputs of sediment from marine and watershed sources. Understanding changes in bed elevations, deposited and eroded sediment, and water depth throughout estuaries is relevant for understanding their present-day status and long-term evolution, identifying potential hazards to human communities, and informing estuarine conservation. In response to observations of sedimentation in the Nehalem Bay, northwestern Oregon, by the Port of Nehalem, the magnitudes and patterns of bathymetric change in the Bay were documented and described by two approaches. The first approach compared changes in bed elevation with estimated volumes of erosion and deposition from overlapping survey data acquired in 1957 and 2019 for the area of the Nehalem Bay from upstream of the Highway 101 bridge to downstream of Fishery Point. The second approach examined changes in water depth for seven zones from the confluence of the North Fork and Nehalem Rivers to the mouth of the Nehalem River using nautical charts (1891, 1947, 1970, 1990, and 2004). These two approaches were used because the bathymetric surveys from 1957 and 2019 could be tied to a common vertical datum, allowing for a direct comparison of changes in bed elevations, whereas the nautical charts could not be tied to a common vertical datum, which limited the analyses to a comparison of changes in water depths over a broader time frame. Bed elevation changes from 1957 to 2019 were assessed from upstream of the Highway 101 bridge to downstream of Fishery Point where the two surveys overlapped (2 square kilometers) using thalweg longitudinal profiles, channel cross sections, and digital elevation models (DEMs) showing the elevation differences between the two surveys (or DEMs of difference). The most prominent change between 1957 and 2019 was the migration of the thalweg (or deepest part of the channel) between the downstream end of Lazarus Island and downstream of Fishery Point; this migration resulted in sediment deposition in the former thalweg and sediment erosion in formerly shallow areas to form the new thalweg. Bed elevation changes in the thalweg also varied longitudinally between 1957 and 2019. The bed elevation of the thalweg in both surveys, however, was generally less than 1 meter (m). The thalweg in the area of overlapping surveys shortened from about 7.0 to 6.7 kilometers in length over that same period. The bed elevation changes between the DEMs showed that maximum erosion and deposition was 4.3 and 4.5 m, respectively. In this same time period, the net change in sediment volume was 230,000 cubic meters (m 3 ), indicating net deposition. However, the error estimated for the 95 percent confidence interval analyses is ±315,000 m 3 , and therefore does not preclude the possibility that net erosion may have occurred. Historical changes in water depth from soundings depicted on nautical charts from 1891, 1947, 1970, 1990, and 2004 were evaluated by assessing spatial and temporal changes for seven zones of the Nehalem Bay. Across all years and zones, water depths ranged from about 0.2 to 9.4 m, whereas median water depths ranged from 0.3 to 6.4 m. Median depths and the range of water depths did not systematically increase or decrease throughout all zones during the same periods. In all nautical charts, the zone at the mouth of the Nehalem River consistently had the deepest soundings (7.9 to 9.4 m) and the greatest range of water depths (7.3 to 8.8 m). Qualitative evaluation of the nautical charts showed minimal changes in the overall shape of the Nehalem Bay. The exception to this observation was at the mouth of the Bay, where two historical outlets to the Pacific Ocean depicted in the 1891 nautical chart were reduced to one outlet following the construction of jetties (1916 and 1918). The results of this study emphasize that bed elevations and water depths within the Nehalem Bay have varied between 1891 and 2019, as illustrated by the lateral and vertical changes in the thalweg and changes in water depths over time. Changes in thalweg position and related patterns of sediment erosion and deposition are expected in the future as the Nehalem Bay continues to respond to changes in tides, sea level, streamflow, and sediment inputs from watershed and marine sources. The results of this study and the surveys from 1957 and 2019 provide a foundation for documenting and evaluating future changes in the Nehalem Bay and prioritizing actions to manage and protect natural resources and recreational access to the Nehalem Bay.

Oregon↗

USGS2004

Integrated monitoring of hydrogeomorphic, vegetative, and edaphic conditions in riparian ecosystems of Great Basin National Park, Nevada

In semiarid regions such as the Great Basin, riparian areas function as oases of cooler and more stable microclimates, greater relative humidity, greater structural complexity, and a steady flow of water and nutrients relative to upland areas. These qualities make riparian areaʼs attractive not only to resident and migratory wildlife, but also to visitors in recreation areas such as Great Basin National Park in the Snake Range, east-central Nevada. To expand upon the system of ten permanent plots sampled in 1992 (Smith et al. 1994) and 2001 (Beever et al. in press ), we established a collection of 31 cross-sectional transects of 50-m width across the mainstems of Strawberry, Lehman, Baker, and Snake creeks. Our aims in this research were threefold: a) map riparian vegetative communities in greater detail than had been done by past efforts; b) provide a monitoring baseline of hydrogeomorphology; structure, composition, and function of upland- and riparianassociated vegetation; and edaphic properties potentially sensitive to management; and c) test whether instream conditions or physiographic variables predicted vegetation patterns across the four target streams. In each of the four watersheds, we performed walking transects from the lower-elevation boundary of the park along creek mainstems to a point well above the point at which vehicle access stopped. In these transects, we ranked, by cover, the riparian and upland woody species on each side of the creek, in 0.32-km segments. These walking transects also facilitated selection of a suite of cross-sectional transects that might serve as an early-warning signal of change for natural (e.g., aggradative) and anthropogenic changes (e.g., due to visitor impacts or climate change). At each cross-sectional transect, we used several methods: a) measurement of the number, approximate volume, and total length of instream logs greater than 10 cm in diameter that were within 5 m up- or downstream of the transect; b) counts of pebbles by size class, following Wolman (1954); c) line-point intercepts, which provided various measures of percent cover; d) gap-intercept transects, following Herrick et al. (in press), to measure susceptibility of uplands to erosion by wind or water; e) 1-m 2 quadrats, to obtain frequency of woody species; f) nested-frequency plots, to measure frequency of all plant species in quadrats of varying size; g) a field-based soil aggregate stability test following Herrick et al. (2001); and h) an impact penetrometer, to measure penetration resistance of soil horizons. We used species-accumulation curves to assess the ability of our methods to detect the majority of plant species at sites, using the most species-rich and species-poor sites as illustrations. We compared characteristics of hydrogeomorphic valley types (designated by Frissell and Liss 1993), vegetation types, and creeks individually and, using multivariate analyses for the first two ʻtypes,ʼ simultaneously. For the latter, using both the nested-frequency and 1-m 2 frequency data, we first used nonmetric multidimensional scaling (NMS) to assess relationships of plant communities among sites. Secondly, we used multi-response permutation procedures (MRPP) to test whether plant-community differences existed among either hydrogeomophic valley types or vegetation types. To increase the value of these comparisons for management, we used indicator species analyses to quantify the indicator value of each individual plant species for separating groups. In contrast to the more incised riparian channels of central Nevada, we observed knickzones, downcutting, and incision only rarely and usually with limited extent in the walking surveys. Downcutting occurred most frequently and extensively in Strawberry and Snake creeks, due in part to their more erodible soils. According to a hydrogeomorphologist with extensive experience in Great Basin riparian systems, the sediment-delivery and hydrologic systems appeared relatively undisturbed in most reaches, with respect to grazing animals and other types of anthropogenic alteration. Site elevation of the 31 transects ranged from 1,950-2,987 m, and stream slope (i.e., gradient) was relatively steep (mean = 9.3%, range 3-16%). Strawberry Creek averaged the lowest maximum water depth, and correspondingly had greatest width/depth ratios. Baker Creek sites averaged the smallest amount of tree-canopy gaps, whereas Snake Creek sites on average had the largest proportion of gaps in understory vegetation. Sites in terrace-bound valley types averaged the lowest slope in the channel as well as the least cover of trees, litter, and vegetation overall, whereas alluviated, boulder-bed canyon sites averaged the greatest widths of the active channel. Sites in Lehman Creek averaged nearly twice as much coarse woody debris as sites from any other creek, whereas Baker Creek sites averaged greatest tree cover (mean = 67%, range 40 – 96%) and species richness (mean = 17.3 species). Multivariate ordinations suggested that sites in leveed outwash valleys and alluvial-fan-influenced valleys had the greatest inter-site heterogeneity in plant composition, whereas sites in incised moraine-filled valleys appeared most homogeneous. Differences among homogeneity of sites within vegetation types were less pronounced, but sites dominated by either aspen and Woodsʼ rose or narrow-leaved cottonwood had the most similar plant communities among sites of the same vegetation type. A number of species were faithful indicators of various valley and vegetation types, using either set of plant-frequency data. We estimate that all 31 sites could be subsequently re-sampled in 14-18 field days by individuals possessing familiarity of the riparian flora of the southern Snake Range. As with any research, monitoring-focused investigations must balance the concerns for number of ecosystem attributes measured, extensiveness in time and space of sampling periods and locations, and the time and cost of sampling.

Nevada↗

USGS2020

Spatially referenced models of streamflow and nitrogen, phosphorus, and suspended-sediment loads in streams of the midwestern United States

In this report, SPAtially Referenced Regression On Watershed attributes (SPARROW) models developed to describe long-term (2000–14) mean-annual streamflow, total nitrogen (TN), total phosphorus (TP), and suspended-sediment (SS) transport in streams of the Midwestern part of the United States (the Mississippi River, Great Lakes, and Red River of the North Basins) are described. The nutrient and suspended-sediment models have a base year of 2012, which means they were developed based on source inputs and management practices similar to those existing during or near 2012 and average hydrological conditions detrended to 2012 (2000–14), whereas the streamflow model has base years of 2000–14, which means it was developed based on the average input precipitation minus actual evapotranspiration from 2000 to 2014. In developing the models, several updates and improvements were made to the data inputs and statistical approaches used to calibrate/develop the models from those used in the previous 2002 SPARROW models. The 2012 SPARROW models were constructed using a higher resolution stream network, which resulted in a mean catchment size of 2.7 square kilometers compared to 480 square kilometers in the 2002 models; more detailed and updated wastewater treatment plant contribution estimates; inputs from background phosphorus sources that were not included in the 2002 model; and more accurate loads for calibration that were computed using a modified Beale ratio-estimator technique whenever no trend in load was determined. Statistical approaches were added to compensate for the unequal effect of each monitoring site during the calibration process by adjusting for the fraction of the basin included in other upstream monitored sites (nested share) and thinning the calibration sites if a negative statistical correlation between nearby sites was determined. Results from 2012 SPARROW models describe how much of each water, TN, TP, and SS source was delivered to the stream network, and the major landscape factors that affected their delivery. Atmospheric deposition and natural (background) sources of TN and TP, respectively, were the dominant sources in anthropogenically unaffected areas (especially in the Rocky Mountains and north-central areas of the Midwest), whereas fertilizers, manure, and fixation were dominant sources in agricultural areas, especially in the Corn Belt and near the Mississippi River. Urban sources of TN and TP were typically localized, but they were still important for some large areas, especially the Lake Erie Basin. All of the land-to-water delivery variables in the nutrient and sediment SPARROW models, such as runoff, soil erodibility, basin slope, and the amount of tile drains, are commonly included in process-driven models. In the SPARROW TN and TP models, best management practices (BMPs) reduced the delivery of these nutrients to streams. Long-term mean-annual flows and nutrient and sediment loads were simulated in streams throughout the Midwest. The simulated flows from the SPARROW flow model were used in the SPARROW TN, TP, and SS models to help describe nutrient and sediment transport from the watershed and through the stream network. Outputs from the TN, TP, and SS models describe loads and yields of these constituents throughout the Midwest, and from major drainage basins throughout the Midwest. Highest TN, TP, and SS yields and delivered yields were from the Lake Erie, Ohio River, Upper Mississippi River, and Lower Mississippi River Basins, whereas lowest yields were spread over most other areas. Losses during downstream delivery resulted in part of the TN, TP, and SS that reach the stream network not reaching the downstream receiving bodies: 14, 15, and 28 percent of the TN, TP, and SS, respectively, are lost during delivery to the Great Lakes and 19, 23, and 52 percent of the TN, TP, and SS, respectively, are lost during delivery to the Gulf of Mexico. The largest losses of nutrients and sediments during transport were in the Missouri and Arkansas River Basins. Information from these SPARROW models can help guide nutrient and sediment reduction strategies throughout the Midwest. Model results provide information on what may be the most appropriate general type of actions to reduce total loading by describing the relative importance of each source, and where to most efficiently place the efforts to reduce loading by describing the distribution of nutrient and sediment loading. By implementing management efforts addressing the major sources of the loads in areas contributing the highest loads, it may be possible to reduce nutrient loading throughout the Mississippi River Basin and thus reduce the size of the hypoxic zone in the Gulf of Mexico; reduce nutrient loading into lakes, and thus reduce the occurrence of harmful algal blooms; and reduce sediment losses, and thus improve the benthic habitat in streams and rivers throughout the Midwest.

Scientific Investigations Report↗

USGS2024

Evaluation of the characteristics, discharge, and water quality of selected springs at Fort Irwin National Training Center, San Bernardino County, California

Eight springs and seeps at Fort Irwin National Training Center were described and categorized by their general characteristics, discharge, geophysical properties, and water quality between 2015 and 2017. The data collected establish a modern (2017) baseline of hydrologic conditions at the springs. Two types of springs were identified: (1) precipitation-fed upland springs (Cave, Desert King, Devouge, No Name, and Panther Springs) and (2) groundwater discharge-fed basin springs (Garlic, Bitter, and Jack Springs). Comparison of electrical resistivity tomography data collected at groundwater basin springs from 2015 to 2017 indicated that spring discharge and connection to the underlying groundwater system is highly focused, although the springs themselves appear diffuse and are spread out over a large area. Spring discharge was consistently less than reported by Thompson (1929), except at Garlic Spring where discharges and vegetation have increased in recent years. Multiple discrete flume and seepage meter measurements taken between October 2015 and April 2016 indicated that discharge changed predictably on diurnal and seasonal timescales in response to evapotranspiration. These preliminary results and the lush vegetation noted at some of the springs, particularly at Bitter, Garlic, and Jack Springs, indicated plant evapotranspiration accounts for a substantial part of the discharge from these springs. The quality of water ranges from fresh in precipitation-fed upland springs (Cave, Desert King, Devouge, and Panther Springs) to slightly saline (Garlic and Jack Springs) and moderately saline (Bitter Spring) in groundwater-fed discharge springs. Nitrate concentrations from water at most of the springs were less than 3 milligrams per liter, except for samples from Devouge and Desert King Springs and one sample from Jack Spring. An analysis of delta nitrogen-15 in nitrate (δ 15 N-NO 3 ) and delta oxygen-18 in nitrate (δ 18 O-NO 3 ) indicates high nitrate concentrations in excess of the U.S. Environmental Protection Agency maximum contaminant level at Jack Spring and Desert King Spring resulting from the dissolution of nitrate-bearing caliche deposits; nitrate concentrations at Devouge Spring are a result of algal growth within the spring, and the source of nitrate concentrations in Garlic Spring are consistent with a treated wastewater origin from Langford Valley-Irwin subbasin upgradient. The source of water in upland springs, indicated by values of delta oxygen-18 (δ 18 O) and delta deuterium (δD) are consistent with recharge from winter precipitation. In groundwater basin springs, values of δ 18 O and δD are consistent with groundwater sampled from nearby wells. Summer monsoonal precipitation appears to contribute little water to spring flow. Most springs contain low levels of tritium and appear to be primarily older (pre-1950s) groundwater. Groundwater basin springs with detectable tritium may result from occasional streamflow in nearby washes. These springs could be susceptible to decreases in flow during extended dry periods when the localized recharge may be reduced due to the loss of focused recharge through nearby washes. Groundwater samples from Garlic and Bitter Springs contained arsenic concentrations above the U.S. Environmental Protection Agency maximum contaminant level. Groundwater samples from all springs, except Cave, Desert King, and Devouge Springs, exceeded the State of California maximum contaminant level for fluoride. Garlic Spring was the only sampled spring that contained vanadium concentrations that exceeded the State of California notification level. Only a single water sample from Jack Spring contained uranium at a concentration that exceeded the U.S. Environmental Protection Agency maximum contaminant level. Many other constituents of concern were analyzed, including those from anthropogenic sources that may be a result of military activities. Most of these constituents were not detected above their respective reporting levels in spring water; only 15 were detected in spring waters. Diesel and gasoline degradants, many of which also occur naturally, were the most commonly detected compounds. Several other organic compounds, primarily solvents or their degradants, were detected in groundwater basin springs. These constituents, in order of decreasing detection frequency, were carbon disulfide; perchlorate; mercury; acetone; methylnaphthalene; toluene; methyl ethyl ketone; cyanide; and styrene; 4-iso-propyl-toluene; isopropylbenzene; methyl salicylate; and phenol. Except for Garlic Spring, which is affected by discharges of treated wastewater, the quality of water from most springs appears to be relatively unaffected by activities at the Fort Irwin National Training Center.

California↗

USGS2015

Hexavalent and total chromium at low reporting concentrations in source-water aquifers and surface waters used for public supply in Illinois, 2013

On the basis of their recent review of the human health effects of hexavalent chromium [Cr(VI)] in public drinking water, the U.S. Environmental Protection Agency is considering the need for Federal regulation of Cr(VI). Presently, only total chromium is regulated, at a Maximum Contaminant Level (MCL) of 100 micrograms per liter (µg/L). The occurrence of Cr(VI) in groundwater and surface waters generally is attributed to industrial sources, but can be of natural origin. California’s recently established MCL for Cr(VI) of 10 µg/L illustrates the drinking-water concerns associated with Cr(VI). To improve understanding of the possible impact of a Cr(VI)-specific standard that approximates the California level on the management of Illinois’ public drinking water, the U.S. Geological Survey, in cooperation with the Illinois Environmental Protection Agency, assessed the occurrence and distribution of Cr(VI) in the State’s public-water supplies. During 2013, untreated water samples were collected to be analyzed for Cr(VI) and total chromium [Cr(T)] at 119 water-supply wells and 32 surface-water intakes; also, 32 treated surface-water samples were collected near the point of treatment and 32 near the furthest point of distribution. Public-supply sample sites were selected by a stratified random method. Samples typically were analyzed within 24 hours of collection at reporting limits of 0.02 µg/L for Cr(VI) and 0.1 µg/L for Cr(T). The occurrence of Cr(VI) was compared with selected geophysical, physical, and sampling factors that might more fully explain its distribution and magnitude of concentrations. The maximum concentration of Cr(VI) in groundwater was 2.1 µg/L. Maximum concentrations in untreated and treated surface water were 0.29 µg/L and 2.4 µg/L, respectively. All sample concentrations were below the California MCL; only 35 percent were below that State’s non-enforceable public health goal of 0.02 of µg/L. Cr(VI) was undetected in 43 percent of untreated groundwater samples, with a median of 0.06 µg/L when detected. All but two (94 percent) of untreated surface-water samples had detections. In untreated surface water, the median concentration was 0.09 µg/L, whereas in treated (tap and distributed) water the median was 0.20 µg/L. Surface waters treated with lime for softening typically had the greatest Cr(VI) concentrations (maximum, 2.4 µg/L; median, 1.2 µg/L). The maximum concentration of Cr(T) in groundwater was 1.8 µg/L. Maximum concentrations in untreated and treated surface water were 1.8 µg/L and 2.5 µg/L, respectively. All sample concentrations were below the Federal MCL. Total chromium was detected in 65 percent of untreated groundwater samples, with a median of 0.40 µg/L, when detected. All but one (97 percent) of untreated surface-water samples had detections. In untreated surface water, the median concentration was 0.40 µg/L, whereas in treated (tap and distributed) water the median was 0.30 µg/L. As with Cr(VI), surface waters treated with lime typically had the greatest Cr(T) concentrations. Examination of factors that might account for or be associated with the occurrence of Cr(VI) in public-supply source waters found few clearly evident factors. Associations in frequencies of occurrence and range of concentrations indicate that surface waters and groundwaters of shallow, unconsolidated, unconfined aquifers, particularly alluvial aquifers, are possibly most commonly affected by anthropogenic sources of Cr(VI). Groundwaters of deep (greater than 500 feet) bedrock aquifers, particularly the Cambrian-Ordovician aquifer system, are possibly most commonly affected by geologic sources of Cr(VI). Additional study, with supporting geologic and geochemical data that were not collected in this study, would be necessary to verify these associations. There was a weak positive relation (ρ = 0.23) between concentrations of Cr(VI) and Cr(T) in untreated water samples, with a much stronger positive relation (ρ = 0.86 and ρ = 0.90, respectively) in samples collected soon after treatment and near the endpoint of distribution. The stronger relation and greater similarity between Cr(VI) and Cr(T) concentrations in treated water samples indicate that Cr(VI) represents a greater proportion of the measured concentrations of Cr(T) in treated waters than in untreated waters. The analysis of spikes and other quality-assurance samples indicate uncertainties associated with obtaining or confirming consistently accurate analytical results for Cr(VI) at near the applied reporting limit of 0.02 µg/L.

Illinois↗

USGS2016

Application of dimensionless sediment rating curves to predict suspended-sediment concentrations, bedload, and annual sediment loads for rivers in Minnesota

Consistent and reliable sediment data are needed by Federal, State, and local government agencies responsible for monitoring water quality, planning river restoration, quantifying sediment budgets, and evaluating the effectiveness of sediment reduction strategies. Heightened concerns about excessive sediment in rivers and the challenge to reduce costs and eliminate data gaps has guided Federal and State interests in pursuing alternative methods for measuring suspended and bedload sediment. Simple and dependable data collection and estimation techniques are needed to generate hydraulic and water-quality information for areas where data are unavailable or difficult to collect. The U.S. Geological Survey, in cooperation with the Minnesota Pollution Control Agency and the Minnesota Department of Natural Resources, completed a study to evaluate the use of dimensionless sediment rating curves (DSRCs) to accurately predict suspended-sediment concentrations (SSCs), bedload, and annual sediment loads for selected rivers and streams in Minnesota based on data collected during 2007 through 2013. This study included the application of DSRC models developed for a small group of streams located in the San Juan River Basin near Pagosa Springs in southwestern Colorado to rivers in Minnesota. Regionally based DSRC models for Minnesota also were developed and compared to DSRC models from Pagosa Springs, Colorado, to evaluate which model provided more accurate predictions of SSCs and bedload in Minnesota. Multiple measures of goodness-of-fit were developed to assess the effectiveness of DSRC models in predicting SSC and bedload for rivers in Minnesota. More than 600 dimensionless ratio values of SSC, bedload, and streamflow were evaluated and delineated according to Pfankuch stream stability categories of “good/fair” and “poor” to develop four Minnesota-based DSRC models. The basis for Pagosa Springs and Minnesota DSRC model effectiveness was founded on measures of goodness-of-fit that included proximity of the model(s) fitted line to the 95-percent confidence intervals of the site-specific model, Nash-Sutcliffe Efficiency values, model biases, and deviation of annual sediment loads from each model to the annual sediment loads calculated from measured data. Composite plots comparing Pagosa Springs DSRCs, Minnesota DSRCs, site-specific regression models, and measured data indicated that regionally developed DSRCs (Minnesota DSRC models) more closely approximated measured data for nearly every site. Pagosa Springs DSRC models had markedly larger exponents (slopes) when compared to the Minnesota DSRC models and the site-specific regression models, and over-represent SSC and bedload at streamflows exceeding bankfull. The Nash-Sutcliffe Efficiency values for the Minnesota DSRC model for suspended-sediment concentrations closely matched Nash-Sutcliffe Efficiency values of the site-specific regression models for 12 out of 16 sites. Nash-Sutcliffe Efficiency values associated with Minnesota DSRCs were greater than those associated with Pagosa Springs DSRCs for every site except the Whitewater River near Beaver, Minnesota site. Pagosa Springs DSRC models were less accurate than the mean of the measured data at predicting SSC values for one-half of the sites for good/fair stability sites and one-half of the sites for poor stability sites. Relative model biases were calculated and determined to be substantial (greater than 5 percent) for Pagosa Springs and Minnesota models, with Minnesota models having a lower mean model bias. For predicted annual suspended-sediment loads (SSL), the Minnesota DSRC models for good/fair and poor stream stability sites more closely approximated the annual SSLs calculated from the measured data as compared to the Pagosa Springs DSRC model. Results of data analyses indicate that DSRC models developed using data collected in Minnesota were more effective at compensating for differences in individual stream characteristics across a variety of basin sizes and flow regimes than DSRC models developed using data collected for Pagosa Springs, Colorado. Minnesota DSRC models retained a substantial portion of the unique sediment signatures for most rivers, although deviations were observed for streams with limited sediment supply and for rivers in southeastern Minnesota, which had markedly larger regression exponents. Compared to Pagosa Springs DSRC models, Minnesota DSRC models had regression slopes that more closely matched the slopes of site-specific regression models, had greater Nash-Sutcliffe Efficiency values, had lower model biases, and approximated measured annual sediment loads more closely. The results presented in this report indicate that regionally based DSRCs can be used to estimate reasonably accurate values of SSC and bedload. Practitioners are cautioned that DSRC reliability is dependent on representative measures of bankfull streamflow, SSC, and bedload. It is, therefore, important that samples of SSC and bedload, which will be used for estimating SSC and bedload at the bankfull streamflow, are collected over a range of conditions that includes the ascending and descending limbs of the event hydrograph. The use of DSRC models may have substantial limitations for certain conditions. For example, DSRC models should not be used to predict SSC and sediment loads for extreme streamflows, such as those that exceed twice the bankfull streamflow value because this constitutes conditions beyond the realm of current (2016) empirical modeling capability. Also, if relations between SSC and streamflow and between bedload and streamflow are not statistically significant, DSRC models should not be used to predict SSC or bedload, as this could result in large errors. For streams that do not violate these conditions, DSRC estimates of SSC and bedload can be used for stream restoration planning and design, and for estimating annual sediment loads for streams where little or no sediment data are available.

Minnesota↗