USGS Science⌕ Search

USGS · sir20215053

Analysis of Escherichia coli, total recoverable iron, and dissolved selenium concentrations, loading, and identifying data gaps for selected 303(d) listed streams, Grand Valley, western Colorado, 1980–2018

Abstract

Tributaries to the Colorado River in the Grand Valley in western Colorado (segment COLCLC13b) have been placed on the State of Colorado 303(d) list as impaired for Escherichia coli (E. coli) , total recoverable iron, and dissolved selenium. The Colorado Department of Public Health and Environment Water Quality Control Division is required to develop total maximum daily loads for these constituents in these tributaries. The U.S. Geological Survey, in cooperation with the Grand Valley Drainage District and Colorado Water Conservation Board, conducted a study to (1) characterize concentrations, loads, and load reductions for E. coli , total recoverable iron, and dissolved selenium using existing data and (2) identify water-quality data gaps to inform future monitoring strategies. This study analyzed water-quality and streamflow data for 3 main-stem sites (2 sites along the Colorado River and 1 site along the Gunnison River) and 29 selected sites on tributaries to the Colorado River. Sample data were available at five sites along Adobe Creek and at six sites along Leach Creek, the two tributaries in the study area that are impaired for E. coli . All geometric mean E. coli concentrations at sites along Adobe Creek and Leach Creek exceeded the State recreational use standard of 126 colony forming units per 100 milliliters (CFU/100 mL). In Adobe Creek, E. coli concentrations in samples ranged from 45.7 to more than 2,420 CFU/100 mL (method upper reporting limit for undiluted samples), and geometric mean concentrations at sites ranged from 301 to 1,180 CFU/100 mL. The E. coli concentrations generally increased in the downstream direction in Adobe Creek; however, increases were not seen between all sites. The largest downstream increase in E. coli concentration was measured between the two most upstream sites. In Leach Creek, concentrations of E. coli in samples ranged from 25.9 to more than 2,420 CFU/100 mL, and geometric mean concentrations at sites ranged from 160 to 259 CFU/100 mL. The E. coli concentrations showed no consistent downgradient increase in Leach Creek. In fact, some of the highest E. coli concentrations were measured at the most upstream site, Leach Creek at Summer Hill Drive. Total recoverable iron concentrations and loads were evaluated at 15 tributary sites for samples collected from August 1993 to February 2018. Median total recoverable iron concentrations ranged from 211 to 4,670 micrograms per liter (µg/L). The chronic aquatic-life water-quality standard (1,000 µg/L) was exceeded in most irrigation season (April through October) samples but was rarely exceeded in nonirrigation season (November through March) samples. Concentrations were often an order of magnitude higher in samples collected during irrigation season than in samples collected during nonirrigation season. None of the sites had enough concurrent total recoverable iron and streamflow data to compute annual loads. As with E. coli , the lack of concurrent total recoverable iron and streamflow information represents a data gap, which needs to be addressed to compute annual loads. Dissolved selenium concentrations and loads were evaluated at 20 tributary sites using discrete water-quality data collected 1991–2018. Dissolved selenium concentrations were higher during nonirrigation season than during irrigation season at tributary sites. However, irrigation season dissolved selenium loads were generally higher than nonirrigation selenium loads, because streamflows were higher during irrigation season. Regression analysis was used to estimate daily dissolved selenium concentrations and loads at three main-stem sites for water years (WYs) 1980–2018 (Gunnison River near Grand Junction and Colorado River near Colorado-Utah State Line) and WYs 2002–18 (Colorado River near Cameo). A trend analysis of dissolved selenium concentrations and loads was completed for these sites from the same respective starting dates but ending in 2017. A continuing downward trend in dissolved selenium concentration was observed at all sites and across all seasonal designations of the analysis. The dissolved selenium concentration decreased by 0.12 µg/L from WY 2002 to 2017 at Colorado River near Cameo, representing an 18-percent decrease during the time period. The dissolved selenium concentration at Gunnison River near Grand Junction decreased by 4.2 µg/L from WY 1980 to 2017, representing a 56-percent decrease overall. During the same time period, dissolved selenium concentration at Colorado River near Colorado-Utah State Line decreased by 3.8 µg/L, representing a 56-percent decrease overall. A downward trend in dissolved selenium load was also observed at all sites and across all seasonal designations of the analysis. The relative contribution of dissolved selenium from the Grand Valley near Grand Junction was estimated by comparing loads at main-stem sites bracketing the study area. The two upstream sites, Colorado River near Cameo and Gunnison River near Grand Junction, contributed 60,300 cumulative pounds and 251,000 cumulative pounds, respectively, during WYs 2002–18. At the furthest downstream site, Colorado River near Colorado-Utah State Line, 490,000 cumulative pounds were estimated during the same time period, indicating that the region between Whitewater and State line contributed approximately 179,000 cumulative pounds or a mean annual load of 10,500 lb/yr. Grand Valley dissolved selenium contributions appear to be stable during WYs 2002–18.

Explore related subjects

90° N90° S · 180° W ← longitude → 180° E
Source-reported bounding extent: 38.736946065676° to 39.470125122358176° latitude; -109.083251953125° to -107.99560546875° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Lisa D. Miller, Rachel G. Gidley, Natalie K. Day, Judith C. Thomas. 2021-06-09. Analysis of Escherichia coli, total recoverable iron, and dissolved selenium concentrations, loading, and identifying data gaps for selected 303(d) listed streams, Grand Valley, western Colorado, 1980–2018. https://doi.org/10.3133/sir20215053

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Karst hydrologic processes in carbonate fensters of the western Great Smoky Mountains, Tennessee, 2017–18

The geology of the Great Smoky Mountains in Tennessee and North Carolina is dominated by metamorphic strata. Along the western extent of the range, however, several carbonate fensters expose the Ordovician Jonesboro Limestone along the western edge of the Great Smoky Fault. The fensters (Cades Cove, Tuckaleechee Cove, and Wear Cove) create opportunities for allogenic recharge at the contact between the surrounding, relatively insoluble, strata and the underlying soluble formations, resulting in the development of karst systems. Karst environments have a high degree of connectivity between surface streams and groundwater systems; thus, karst aquifers are inherently vulnerable to contamination from surface activities. Karst environments also provide a specific environment for sensitive biota. Protection of the unique resources of karst environments requires an understanding of the behavior of karst aquifers, knowledge of karst recharge areas, and the mapping of groundwater flow paths. A study began in 2017 to characterize the hydrologic behavior of the karst systems in Great Smoky Mountains National Park through hydrologic monitoring, dye tracing, and seepage investigations. Hydrologic monitoring equipment was installed at several karst hydrogeologic features: two in-cave sites, two major springs, Abrams Creek in Cades Cove, and a karst depression wetland. Dye tracing was conducted to determine resurgences for the cave systems, to delineate spring recharge areas, and to estimate karst groundwater velocities. Seepage investigations along Abrams Creek helped determine if streams were sinking into the Jonesboro Limestone, either at the contact with insoluble strata or as streams flowed across the carbonate fenster. Study results indicate considerable differences in groundwater behavior between the conduit-dominated flow systems in Tuckaleechee Cove compared to the lower-gradient, colluvium-covered karst systems in Cades Cove. Study results also indicate that contaminants in the area of Bull sink and Whiteoak Sink could be rapidly transported through karst groundwater towards springs in Tuckaleechee Cove. Karst resources in Cades Cove are most likely to be affected by local events or activities taking place within the cove itself or in upper Abrams Creek upstream of the Great Smoky Fault. This study provides information on the vulnerability and uniqueness of the karst systems in the western Great Smoky Mountains, which are characterized by complex geology, high-gradient subterranean streams, deep caves, and carbonate fensters.

Tennessee↗

Analysis of groundwater geochemistry in the Lansing, Michigan, area

Groundwater resources in the Lansing, Michigan, area are protected by a wellhead protection plan, but previous studies have indicated that groundwater in this area may be susceptible to surface contamination. The U.S. Geological Survey, in cooperation with the Lansing Board of Water & Light, analyzed geochemical and aquifer characteristics to assess the potential for contamination of groundwater resources. Recharge age categories were determined from tritium concentrations for some wells, and by using ternary diagrams and principal component analysis for other wells. These age classifications, combined with hydrogeologic and geochemical data, were analyzed by using statistical methods to identify trends that may indicate areas of increased groundwater susceptibility. Wells with modern recharge tended to have higher sodium, chloride, and chloride:bromide ratios than did wells with mixed recharge, suggesting greater surface connectivity and potential influence from anthropogenic sources such as road salt. In contrast, wells in the mixed age category exhibited more stable geochemistry and lower constituent concentrations. Classification and regression tree models identified aquifer geometry—particularly the thickness of the uppermost bedrock unit and presence of confining layers—as key predictors of recharge age and geochemical indicators. The trends identified from these analyses could help to reveal locations in the Lansing area that are among the most susceptible to contamination.

Michigan↗

Evaluating the number of realizations (train-test splits) necessary to ensure robust estimates of conventional hydrothermal resource favorability using Monte Carlo cross-validation, the Great Basin, Western USA

Supervised machine learning approaches used to produce favorability maps in resource assessments commonly partition an input dataset into training and testing subsets (that is, through train-test splits). However, if a single split is not representative of underlying resource conditions, the resulting favorability map may be biased, particularly with sparse or heterogeneous datasets. To mitigate the potential for unrepresentative favorability maps, recent works by U.S. Geological Survey researchers have used Monte Carlo cross-validation to implement a large number ( n) of randomized train-test splits to produce n realizations. The realizations are then aggregated using cell-wise median predictions to produce the resulting final favorability map, recognizing that the median will remove extreme values that might have been estimated from a comparatively few unrepresentative realizations. Theoretically, the median prediction at every cell converges to a single estimate as n approaches infinity, and the question addressed herein is: How many realizations are necessary to get a robust estimate (that is, an estimate unaffected by additional realizations)? Therefore, convergence metrics are described to identify the number of train-test splits and realizations needed to make a reliable map using the data and supervised machine learning workflow for conventional hydrothermal resource assessment of the Great Basin of the Western United States. For this study, robust prediction occurs with n approximately equal to 600 realizations, and smaller ensembles (for example, n is approximately equal to 200 or even 100) may be suitable during preliminary or exploratory analyses. Although the resulting threshold numbers will vary depending on the data and modeling approach used, the convergence metrics proposed and applied herein are applicable to other supervised machine learning studies that use Monte Carlo cross-validation.

California, Idaho, Nevada, Oregon, Utah↗