USGS Science⌕ Search

SEARCH · USGS Science

Results for “Environmental Modeling & Assessment”

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

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

At least 1,495 records · Page 83Linked to original sources

Assessment of macroinvertebrate communities in adjacent urban stream basins, Kansas City, Missouri, metropolitan area, 2007 through 2011

Macroinvertebrates were collected as part of two separate urban water-quality studies from adjacent basins, the Blue River Basin (Kansas City, Missouri), the Little Blue River and Rock Creek Basins (Independence, Missouri), and their tributaries. Consistent collection and processing procedures between the studies allowed for statistical comparisons. Seven Blue River Basin sites, nine Little Blue River Basin sites, including Rock Creek, and two rural sites representative of Missouri ecological drainage units and the area’s ecoregions were used in the analysis. Different factors or levels of urban intensity may affect the basins and macroinvertebrate community metrics differently, even though both basins are substantially developed above their downstream streamgages (Blue River, 65 percent; Little Blue River, 52 percent). The Blue River has no flood control reservoirs and receives wastewater effluent and stormflow from a combined sewer system. The Little Blue River has flood control reservoirs, receives no wastewater effluent, and has a separate stormwater sewer system. Analysis of macroinvertebrate community structure with pollution-tolerance metrics and water-quality parameters indicated differences between the Blue River Basin and the Little Blue River Basin. A four-metric score (total taxa richness, Ephemeroptera plus Plecoptera plus Trichoptera taxa richness, Macroinvertebrate Biotic Index, and Shannon Diversity Index) for richest-targeted habitat was used to calculate a Stream Condition Index (SCI) in order to evaluate the aquatic-life status of the streams. About 80 percent of all samples combined were determined to be less than fully biologically supporting, and about 11 percent of spring samples were fully biologically supporting. No sites within the Blue River Basin had a fully supporting score. The aquatic-life status scores for the Little Blue River and its tributaries were higher (indicating more optimal conditions) than for the Blue River and its tributaries. Fall samples scored higher than spring samples. However, fall samples were collected at the Little Blue River Basin and rural sites only. The Little Blue River sites scored higher for fall samples than spring samples; about 39 percent fully biologically supporting and 61 percent partially biologically supporting; more similar to the rural comparison sites, 40 percent fully biologically supporting and 60 percent partially biologically supporting. The SCI was compared to other multimetric indices with more or other component metrics to determine if the SCI effectively described differences among sites. Environmental variables (streamflow, water quality, land use, impervious cover, and population density) were used in statistical analyses to evaluate relations to macroinvertebrate metrics. Multimetric indices (MMIs) were modeled using step regression with a simple urban intensity index (SUII) based on percentage of impervious cover, population density, and forest cover in a 30-meter stream-buffer zone, and two were selected for further analysis. Three other multimetric indices composed of metrics common to local and national studies show results similar to the two modeled MMIs. A common Benthic Index of Biotic Integrity (R 2 equals 0.71) developed for a national study had the highest correlation with urban intensity as measured with the SUII, followed by a modeled 6-metric index (R 2 equals 0.61). The other MMIs and the SCI explained less than a half of the variability in macroinvertebrate communities in relation to the SUII. Wastewater-treatment plant discharges during base flow, which elevated specific conductance and nutrient concentrations, combined sewer overflows, and nonpoint sources likely contributed to water-quality impairment and lower aquatic-life status at the Blue River Basin sites. Releases from upstream reservoirs to the Little Blue River likely decreased specific conductance, suspended-sediment, and dissolved constituent concentrations and may have benefitted water quality and aquatic life of main-stem sites. Chloride concentrations in base-flow samples, attributable to winter road salt application, had the highest correlation with the SUII (Spearman’s ρ equals 0.87), were negatively correlated with the SCI (Spearman’s ρ equals -0.53) and several pollution sensitive Ephemeroptera plus Plecoptera plus Trichoptera abundance and percent richness metrics, and were positively correlated with pollution tolerant Oligochaeta abundance and percent richness metrics. Study results show that the easily calculated SUII and the selected modeled multimetric indices are effective for comparing urban basins and for evaluation of water quality in the Kansas City metropolitan area.

Missouri↗

Groundwater quality, age, and susceptibility and vulnerability to nitrate contamination with linkages to land use and groundwater flow, Upper Black Squirrel Creek Basin, Colorado, 2013

The Upper Black Squirrel Creek Basin is located about 25 kilometers east of Colorado Springs, Colorado. The primary aquifer is a productive section of unconsolidated deposits that overlies bedrock units of the Denver Basin and is a critical resource for local water needs, including irrigation, domestic, and commercial use. The primary aquifer also serves an important regional role by the export of water to nearby communities in the Colorado Springs area. Changes in land use and development over the last decade, which includes substantial growth of subdivisions in the Upper Black Squirrel Creek Basin, have led to uncertainty regarding the potential effects to water quality throughout the basin. In response, the U.S. Geological Survey, in cooperation with Cherokee Metropolitan District, El Paso County, Meridian Service Metropolitan District, Mountain View Electric Association, Upper Black Squirrel Creek Groundwater Management District, Woodmen Hills Metropolitan District, Colorado State Land Board, and Colorado Water Conservation Board, and the stakeholders represented in the Groundwater Quality Study Committee of El Paso County conducted an assessment of groundwater quality and groundwater age with an emphasis on characterizing nitrate in the groundwater. Groundwater-quality samples were collected from 50 randomly selected wells between May and June 2013. The samples were analyzed for major ions, nutrients, dissolved gases, tritium ( 3 H), chlorofluorocarbons (CFC-11, CFC-12, and CFC-113), and fuel products (such as benzene, toluene, ethylbenzene, and xylenes). None of the groundwater samples exceeded the U.S. Environmental Protection Agency (EPA) National Primary Drinking Water Regulations for primary maximum contaminant levels (MCL) for major ions. Secondary maximum contaminant levels, which are not health concerns and affect mainly taste, color, or odor of the water, were observed in rare instances for pH (2 samples), chloride (1 sample), iron (3 samples), and manganese (8 samples). The secondary maximum contaminant level for total dissolved solids was also exceeded for two samples. Nitrate (nitrite plus nitrate as nitrogen in groundwater) was elevated above the estimated background concentration of natural recharge waters of 1 milligram per liter (mg/L) in 44 of the 50 wells sampled and showed a median concentration of 5.4 mg/L. Nitrate concentrations were above the MCL of 10 mg/L in 5 of the 50 wells sampled and above half of the EPA MCL (5 mg/L) in 27 of the 50 wells sampled, which included samples above the MCL. Dissolved-oxygen concentrations exceeded 0.5 mg/L in 95 percent of reported values (40 of 42 samples) and exceeded 2.0 mg/L in 90 percent of reported values (38 of 42 samples). The oxidized conditions observed in most areas indicate that nitrate from fertilizers and animal or human waste was geochemically stable and could persist in the groundwater for decades or perhaps longer. A historical analysis of median nitrate concentrations over nearly three decades showed an increase in nitrate of approximately 1 mg/L from 4.3 to 5.4 mg/L, although the increase was not determined to be significantly different using nonparametric statistical methods. Major-ion data indicate that groundwater representative of the primary aquifer was classified as calcium-sodium bicarbonate type water. Other water samples from wells located mainly along the periphery of the primary aquifer had cation-anion compositions consistent with distinct water sources, including groundwater contributions from the underlying bedrock aquifers. The areas with differentiable water sources were located mainly where alluvial deposits were thin and geologic contacts to the underlying bedrock aquifers were relatively shallow. Nitrate concentrations in the groundwater were evaluated for relations to land use. An agricultural region was defined using a sequence of land satellite imagery. Groundwater flow directions interpreted from median water-table elevations measured from 2000 to 2013 were used in conjunction with cropland locations to define the agricultural region boundaries by encompassing potential pathways of nitrate transport in the groundwater from nitrogen-based fertilizers. A statistically significant higher median nitrate concentration was observed for areas inside the agricultural region (6.7 mg/L) compared to areas outside the agricultural region (2.3 mg/L), although median concentrations in both areas were below the MCL (10 mg/L). Median nitrate concentration was also significantly greater in land parcels with septic use (4.9 mg/L) compared to nonseptic parcels (1.7 mg/L). In general, agriculture or septic use was identified as the primary source of nitrate, depending on location, while commercial, county, grazing, and residential land uses were generally secondary sources of nitrate. Apparent groundwater ages were estimated from chlorofluorocarbons (CFC-11, CFC-12, and CFC-113) and tritium ( 3 H) data using models that assumed piston flow and binary mixing (dilution of a young component with old, tracer-free water). The mean and median groundwater ages were about 30 years and the standard deviation was 6 years, indicating that most groundwater in the primary aquifer was “young” water that had recharged to the aquifer over the last few decades (post-1950s). The median fraction of young water was about 71 percent, and the standard deviation was 29 percent. The remaining water predated the 1950s, which may have originated from deeper geologic formations or may represent slow moving groundwater within the primary aquifer. Some of the oldest groundwater ages (older than 30 years) were observed in the upper reaches of the aquifer to the northwest where the primary aquifer is thin and intersects bedrock, supporting the hypothesis of geochemically distinct groundwater entering the primary aquifer from below. Groundwater that had reached the central part of the aquifer from upgradient areas of the basin was variable in age because of differences in flow paths and travel velocities. The groundwater age analysis showed that current (2013) land-use practices could affect water quality over decades to come, and that responses to remedial actions could be slow, especially for constituents, such as nitrate, that are stable under oxidized conditions. Fuel products (including acetone, benzene, diisopropyl ether, ethylbenzene, methyl acetate, methyl tertiary butyl ether (MTBE), methyl tert-pentyl ether, m- + p-xylene, o-xylene, tert-amyl alcohol, tert-butyl alcohol, tert-butyl ethyl ether, and toluene) were analyzed in groundwater from 49 of the 50 wells. Water from seven sites had detections for fuel compounds; all concentrations were below MCL. The results provided assurance of water quality and a valuable baseline to evaluate future trends of fuel constituents as the region is further developed. Probability maps were developed from logistic regression models to examine the likelihood that nitrate concentrations in groundwater exceeded specified levels. Susceptibility analysis examined relations between mid-level (5.0 mg/L) nitrate concentrations and climatic, hydrologic, and geologic variables; the significant variables were identified as depth to groundwater, soil organic matter, and soil water storage to 25-centimeter (cm) depth. The vulnerability assessments included natural factors driving susceptibility but also human factors related to land use and septic use. Vulnerability to low-level (2.5 mg/L) nitrate was related to depth to groundwater, septic zoning, and soil organic matter. The results highlighted that septic zoning affected low-level nitrate concentrations. Vulnerability to mid-level (5.0 mg/L) nitrate was examined using all 50 samples and also with two data outliers removed, which showed relatively high nitrate concentrations but also anomalous water chemistry or were located beyond the primary study area. Vulnerability to mid-level (5.0 mg/L) nitrate using all 50 samples was related to depth to groundwater, land use, septic use within a 500-meter (m) radius, soil water storage to a 25-cm depth, soil organic matter, and whether a location was within the agricultural region. The mid-level (5.0 mg/L) vulnerability model using 48 samples (two outliers removed) produced the best overall fit and was related to the same variables as when using all samples except septic use. The results for mid-level vulnerability provided additional support that septic use was associated with low levels of nitrate in the groundwater. Soil properties and land use were identified as the main drivers of moderate nitrate concentrations. Probabilities of exceeding low-level nitrate concentrations were high in most areas with the lowest probabilities usually to the northwest along thin geologic deposits in the upper part of the basin. The results of this investigation offer the foundational information needed for developing best management practices to mitigate nitrate contamination, basic concepts on water quality to aid public education, and information to guide regulatory measures if policy makers determine this is warranted. Science-based decision making will require continued monitoring and analysis of water quality in the future.

Colorado↗

Assessing the topographic distribution of legacy soil phosphorus in agricultural fields of the Delmarva Peninsula, Mid-Atlantic Coastal Plain, USA

Phosphorus (P) management remains a challenge in agricultural watersheds. The Choptank River Conservation Effects Assessment Project watershed, located in Maryland and Delaware and draining to the Chesapeake Bay, contains legacy soil P from historical dairy and poultry manure applications. These practices elevated soil P beyond crop needs, contributing to persistent P export to aquatic ecosystems. We assessed spatial P distribution and analyzed GIS (Geographic Information Systems)-derived landscape features driving legacy P movement on a farm (47 ha). We hypothesized that P accumulates in drained lowlands and depressional areas due to gravity-driven processes that accelerate P-enriched water to receiving waters via overland flow. In collaboration with the US Department of Agriculture Legacy P project, we collected 105 soil samples (0- to 5-cm and 5- to 15-cm depths) and 14 ditch sediment samples across five topographic openness classes from a farm with >100 years of dairy manure application. Average Mehlich-III P concentrations were 218 and 179 mg kg −1 at 0- to 5-cm and 5- to 15-cm depths, respectively, with legacy areas defined by P content > 100 mg kg −1 . Soil P and clay particle size were positively correlated ( r = 0.42, p < 0.05), increased as landscape openness decreased, and were negatively correlated with topographic openness (ranging from −0.2 to −0.4, p < 0.05), indicating accumulation of P and clay in low-lying areas. These patterns suggest that historical field-level managements have primarily shaped P distribution, while hydrologic and landscape properties further influence its redistribution via transport pathways and drainage. These findings support the development of landscape models to map critical source areas in low-relief watersheds and guide targeted mitigation in high-risk P export zones.

Maryland↗

Assessment of water quality trends in the Minnesota River using non-parametric and parametric methods

Excessive loading of sediment and nutrients to rivers is a major problem in many parts of the United States. In this study, we tested the non-parametric Seasonal Kendall (SEAKEN) trend model and the parametric USGS Quality of Water trend program (QWTREND) to quantify trends in water quality of the Minnesota River at Fort Snelling from 1976 to 2003. Both methods indicated decreasing trends in flow-adjusted concentrations of total suspended solids (TSS), total phosphorus (TP), and orthophosphorus (OP) and a generally increasing trend in flow-adjusted nitrate plus nitrite-nitrogen (NO 3 -N) concentration. The SEAKEN results were strongly influenced by the length of the record as well as extreme years (dry or wet) earlier in the record. The QWTREND results, though influenced somewhat by the same factors, were more stable. The magnitudes of trends between the two methods were somewhat different and appeared to be associated with conceptual differences between the flow-adjustment processes used and with data processing methods. The decreasing trends in TSS, TP, and OP concentrations are likely related to conservation measures implemented in the basin. However, dilution effects from wet climate or additional tile drainage cannot be ruled out. The increasing trend in NO 3 -N concentrations was likely due to increased drainage in the basin. Since the Minnesota River is the main source of sediments to the Mississippi River, this study also addressed the rapid filling of Lake Pepin on the Mississippi River and found the likely cause to be increased flow due to recent wet climate in the region. Copyright ?? 2009 by the American Society of Agronomy, Crop Science Society of America, and Soil Science Society of America. All rights reserved.

Journal of Environmental Quality↗

Local water use and climate variability drive water stress and alter ecological flows over the conterminous United States

Consistent, large-scale estimates of water availability are needed to identify and avoid potential conflicts among human and ecosystem uses of water. We present an assessment of water limitation, defined as the monthly balance (difference) between water supply ( ws ) and human consumptive water use ( wc ), for the conterminous United States (CONUS) during water years 2010–2020. We estimate that 26.7 million Americans, 8% of CONUS population, live in areas with chronic high or severe water limitation. Although ws greatly exceeds wc at the CONUS scale, water is limited locally or regionally due to spatial and temporal patterns in climate and wc . Our water limitation metric, the monthly supply and use index (SUI), peaked in 2012 during a widespread drought when 38% of the CONUS land area experienced elevated water stress. The central and Southwestern U.S. experienced the highest SUI due to the combination of low ws and high wc , especially for irrigation. Spatial overlays of SUI and fish habitat ranges, including those of conservation concern, revealed that several species had notable proportions of their habitat exposed to high or severe water limitation during spawning season over the modeled time period, especially the Arkansas River shiner. ws was calculated from two CONUS, physically-based, hydrologic models while wc was calculated from three CONUS models of water use for crop irrigation, thermoelectric power generation, and public supply. The ws and wc values were routed through a stream network and compared to calculate water limitation and SUI for human populations and fish species at the scale of 12-digit hydrologic unit codes. Evaluation of water availability at higher spatial and temporal resolution promotes more comprehensive analyses of the drivers of water availability and can be combined with complementary studies of water quality and water limiting thresholds to better understand the limitations on water availability.

conterminous United States↗

Hydropower assessment of Bolivia—A multisource satellite data and hydrologic modeling approach

This study produced a geospatial database for use in a decision support system by the Bolivian authorities to investigate further development and investment potentials in sustainable hydropower in Bolivia. The study assessed theoretical hydropower of all 1-kilometer (km) stream segments in the country using multisource satellite data and a hydrologic modeling approach. With the assessment covering the 2 million square kilometer (km 2 ) region influencing Bolivia’s drainage network, the potential hydropower figures are based on theoretical yield assuming that the systems generating the power are 100 percent efficient. There are several factors to consider when determining the real-world or technical power potential of a hydropower system, and these factors can vary depending on local conditions. Since this assessment covers a large area, it was necessary to reduce these variables to the two that can be modeled consistently throughout the region, streamflow or discharge, and elevation drop or head. First, the Shuttle Radar Topography Mission high-resolution 30-meter (m) digital elevation model was used to identify stream segments with greater than 10 km 2 of upstream drainage. We applied several preconditioning processes to the 30-m digital elevation model to reduce errors and improve the accuracy of stream delineation and head height estimation. A total of 316,500 1-km stream segments were identified and used in this study to assess the total theoretical hydropower potential of Bolivia. Precipitation observations from a total of 463 stations obtained from the Bolivian Servicio Nacional de Meteorología e Hidrología (Bolivian National Meteorology and Hydrology Service) and the Brazilian Agência Nacional de Águas (Brazilian National Water Agency) were used to validate six different gridded precipitation estimates for Bolivia obtained from various sources. Validation results indicated that gridded precipitation estimates from the Tropical Rainfall Measuring Mission (TRMM) reanalysis product (3B43) had the highest accuracies. The coarse-resolution (25-km) TRMM data were disaggregated to 5-km pixels using climatology information obtained from the Climate Hazards Group Infrared Precipitation with Stations dataset. About a 17-percent bias was observed in the disaggregated TRMM estimates, which was corrected using the station observations. The bias-corrected, disaggregated TRMM precipitation estimate was used to compute stream discharge using a regionalization approach. In regionalization approach, required homogeneous regions for Bolivia were derived from precipitation patterns and topographic characteristics using a k -means clustering approach. Using the discharge and head height estimates for each 1-km stream segment, we computed hydropower potential for 316,490 stream segments within Bolivia and that share borders with Bolivia. The total theoretical hydropower potential (TTHP) of these stream segments was found to be 212 gigawatts (GW). Out of this total, 77.4 GW was within protected areas where hydropower projects cannot be developed; hence, the remaining total theoretical hydropower in Bolivia (outside the protected areas) was estimated as 135 GW. Nearly 1,000 1-km stream segments, however, were within the boundaries of existing hydropower projects. The TTHP of these stream segments was nearly 1.4 GW, so the residual TTHP of the streams in Bolivia was estimated as 133 GW. Care should be exercised to understand and interpret the TTHP identified in this study because all the stream segments identified and assessed in this study cannot be harnessed to their full capacity; furthermore, factors such as required environmental flows, efficiency, economics, and feasibility need to be considered to better identify a more real-world hydropower potential. If environmental flow requirements of 20–40 percent are considered, the total theoretical power available reduces by 60–80 percent. In addition, a 0.72 efficiency factor further reduces the estimation by another 28 percent. This study provides the base theoretical hydropower potential for Bolivia, the next step is to identify optimal hydropower plant locations and factor in the principles to appraise a real-world power potential in Bolivia.

Open-File Report↗

Methods for estimating adsorbed uranium(VI) and distribution coefficients of contaminated sediments

Assessing the quantity of U(VI) that participates in sorption/desorption processes in a contaminated aquifer is an important task when investigating U migration behavior. U-contaminated aquifer sediments were obtained from 16 different locations at a former U mill tailings site at Naturita, CO (U.S.A.) and were extracted with an artificial groundwater, a high pH sodium bicarbonate solution, hydroxylamine hydrochloride solution, and concentrated nitric acid. With an isotopic exchange method, both a KD value for the specific experimental conditions as well as the total exchangeable mass of U(VI) was determined. Except for one sample, KD values determined by isotopic exchange with U-contaminated sediments that were in equilibrium with atmospheric CO2 agreed within a factor of 2 with KD values predicted from a nonelectrostatic surface complexation model (NEM) developed from U(VI) adsorption experiments with uncontaminated sediments. The labile fraction of U(VI) and U extracted by the bicarbonate solution were highly correlated (r2 = 0.997), with a slope of 0.96 ?? 0.01. The proximity of the slope to one suggests that both methods likely access the same reservoir of U(VI) associated with the sediments. The results indicate that the bicarbonate extraction method is useful for estimating the mass of labile U(VI) in sediments that do not contain U(IV). In-situ KD values calculated from the measured labile U(VI) and the dissolved U(VI) in the Naturita alluvial aquifer agreed within a factor of 3 with in-situ K D values predicted with the NEM and groundwater chemistry at each well.

Environmental Science & Technology↗

A review of models and micrometeorological methods used to estimate wetland evapotranspiration

Within the past decade or so, the accuracy of evapotranspiration (ET) estimates has improved due to new and increasingly sophisticated methods. Yet despite a plethora of choices concerning methods, estimation of wetland ET remains insufficiently characterized due to the complexity of surface characteristics and the diversity of wetland types. In this review, we present models and micrometeorological methods that have been used to estimate wetland ET and discuss their suitability for particular wetland types. Hydrological, soil monitoring and lysimetric methods to determine ET are not discussed. Our review shows that, due to the variability and complexity of wetlands, there is no single approach that is the best for estimating wetland ET. Furthermore, there is no single foolproof method to obtain an accurate, independent measure of wetland ET. Because all of the methods reviewed, with the exception of eddy covariance and LIDAR, require measurements of net radiation (Rn) and soil heat flux (G), highly accurate measurements of these energy components are key to improving measurements of wetland ET. Many of the major methods used to determine ET can be applied successfully to wetlands of uniform vegetation and adequate fetch, however, certain caveats apply. For example, with accurate Rn and G data and small Bowen ratio (??) values, the Bowen ratio energy balance method can give accurate estimates of wetland ET. However, large errors in latent heat flux density can occur near sunrise and sunset when the Bowen ratio ?? ??? - 1??0. The eddy covariance method provides a direct measurement of latent heat flux density (??E) and sensible heat flux density (II), yet this method requires considerable expertise and expensive instrumentation to implement. A clear advantage of using the eddy covariance method is that ??E can be compared with Rn-G H, thereby allowing for an independent test of accuracy. The surface renewal method is inexpensive to replicate and, therefore, shows particular promise for characterizing variability in ET as a result of spatial heterogeneity. LIDAR is another method that has special utility in a heterogeneous wetland environment, because it provides an integrated value for ET from a surface. The main drawback of LIDAR is the high cost of equipment and the need for an independent ET measure to assess accuracy. If Rn and G are measured accurately, the Priestley-Taylor equation can be used successfully with site-specific calibration factors to estimate wetland ET. The 'crop' cover coefficient (Kc) method can provide accurate wetland ET estimates if calibrated for the environmental and climatic characteristics of a particular area. More complicated equations such as the Penman and Penman-Monteith equations also can be used to estimate wetland ET, but surface variability and lack of information on aerodynamic and surface resistances make use of such equations somewhat questionable. ?? 2004 John Wiley and Sons, Ltd.

Hydrological Processes↗

Ground-truthing electrical resistivity methods in support of submarine groundwater discharge studies: Examples from Hawaii, Washington, and California

Submarine groundwater discharge (SGD) is an important conduit that links terrestrial and marine environments. SGD conveys both water and water-borne constituents into coastal waters, where these inflows may impact near-shore ecosystem health and sustainability. Multichannel electrical resistivity techniques have proven to be a powerful tool to examine scales and dynamics of SGD and SGD forcings. However, there are uncertainties both in data aquisition and data processing that must be addressed to maximize the effectiveness of this tool in estuarine or marine environments. These issues most often relate to discerning subtle nuances in the flow of electricity through variably saturated media that can also be highly conductive ( i.e., seawater). Three contrasting field sites were examined for this study to assess the effectiveness of electrical resistivity techniques in varying coastal settings by comparing resistivity data to direct salinity and resistivity observations, quantifying changes in lithology and beach geomorphology, and fine-tuning inversion protocols. The three study sites all have substantial (up to 85 cm day &minus;1 ) submarine groundwater discharge rates, but the hydrologic, oceanographic, and geologic characteristics of the sites are all very different. At a site in Pelekane Bay on the Big Island of Hawaii, seasonal flooding introduces very high concentrations of fine to coarse sediment into the bay. Near-shore circulation is limited in Pelekane Bay, so this newly introduced sediment can become deposited in the bay where it accumulates over time. At a site in Hood Canal, a fjord within Puget Sound, Washington, SGD rates can be high because of the large tidal range, abundant recharge, and steep hydrologic gradients. At Younger Lagoon in northern California, the flow of groundwater towards the coast is much more parsimonious, but here marine processes, including recirculated seawater, are important in controlling the flow of material towards the coast. Rigorous ground-truthing at each field site showed that multi-channel electrcial resistivity techniques can reproduce the scales and dynamics of a seepage field when such data are correctly collected, and when the model inversions are tuned to field site characteristics. Such information can provide a unique perspective on the scales and dynamics of exchange processes within a coastal aquifer&mdash;information essential to scientists and resource managers alike.

California, Hawaii, Washington↗

Status and trends of total nitrogen and total phosphorus concentrations, loads, and yields in streams of Mississippi, water years 2008–18

To assess the status and trends of conditions of surface waters throughout Mississippi, the U.S. Geological Survey, in cooperation with the Mississippi Department of Environmental Quality (MDEQ), summarized concentrations and estimated loads, yields, trends, and spatial and temporal patterns of total nitrogen (TN) and total phosphorus (TP) at 20 stream sites in MDEQ’s ambient water-quality monitoring network and 2 stream sites in the U.S. Geological Survey’s National Water-Quality Assessment Project’s monitoring network. Comparison of streamflow at the time of water-quality sample collection to flow-duration curves for each site showed that samples were relatively evenly spread over a wide range of flows, indicating that load estimations were representative of a wide range of flows. Relation of streamflow to concentrations of TN and TP varied among sites and land use. Sites with high agriculture land use in the drainage basin tended to have a positive correlation between streamflow and concentration, suggesting influence of event-driven nonpoint-source runoff. Sites near urban (developed) areas tended to have a negative correlation between streamflow and concentration, suggesting chronic point-source influences during low-flow conditions. Sites with high forest land use and lower agriculture and urban (developed) land use showed little to no association between streamflow and concentration. Seasonal distributions of concentrations of TN and TP also corresponded closely with variations in land use. Sites near urban (developed) land had the highest concentrations in late summer and fall, sites with a high percentage of agricultural land had the highest concentrations in the spring, and sites that were primarily forested or with little developed land did not exhibit substantial changes in concentration across seasons. Eight sites had statistical likelihoods for upward trends of TN loads, and seven sites had statistical likelihoods for downward trends. Trends in TN loads at six sites were considered “about as likely as not,” meaning that a site has an equal chance of having an upward or downward trend. Trend results of mean annual flow-normalized loads of TP for the period of analysis (2008–18) showed that 16 sites had upward trends, 3 sites had downward trends, and 2 sites were considered “about as likely as not.” Results from our study were compared to results from existing regional models to assess accuracy of predictions at a local scale. Comparisons of yields predicted from 2012 regional-scale SPAtially Referenced Regressions on Watershed attributes (SPARROW) to results from this study showed the 2012 SPARROW-predicted estimates varied in consistency with results from this study. The 2012 SPARROW-prediction model underestimated TN yields, more often and by a slightly larger degree, more than it overestimated TN yields. The 2012 SPARROW-predicted model tended to underestimate yields at study sites with higher yields. All four sites in the predominantly agricultural area of northwest Mississippi, locally known as the Mississippi Delta, were underestimated by 2012 SPARROW. For TP, yield comparisons at sites with lower yields were consistent, yields at sites with midrange yields tended to be overestimated by SPARROW, and yields at sites with high yields tended to be underestimated by SPARROW. TP yields at four sites in the Mississippi Delta were underestimated by the 2012 SPARROW-predicted model. Results of select sites from our study were also compared to other published load estimates from an earlier time period to evaluate possible trends. Comparison of TN yields at four sites and TP yields at three sites from the study-derived estimates to estimates made from data spanning 1993–2004 showed decreasing TN yields at all four sites and decreasing TP yields at two of three sites, with increasing yields of TP at the Yazoo River lower site. Also, a third comparison of the TN and TP yields of the Yazoo River lower site of this study to estimates made from data spanning 1996–97 showed decreasing TN yields but similar TP yields. This suggests that TN yields may have decreased over the last 20–30 years, but TP yields remain constant or are increasing.

Mississippi↗

Analysis of data on nutrients and organic compounds in ground water in the upper Snake River basin, Idaho and western Wyoming, 1980-91

Nutrient and organic compound data from the U.S. Geological Survey and the U.S. Environmental Protection Agency STORET data bases provided information for development of a preliminary conceptual model of spatial and temporal ground-water quality in the upper Snake River Basin. Nitrite plus nitrate (as nitrogen; hereafter referred to as nitrate) concentrations exceeded the Federal drinking-water regulation of 10 milligrams per liter in three areas in Idaho" the Idaho National Engineering Laboratory, the area north of Pocatello (Fort Hall area), and the area surrounding Burley. Water from many wells in the Twin Falls area also contained elevated (greater than two milligrams per liter) nitrate concentrations. Water from domestic wells contained the highest median nitrate concentrations; water from industrial and public supply wells contained the lowest. Nitrate concentrations decreased with increasing well depth, increasing depth to water (unsaturated thickness), and increasing depth below water table (saturated thickness). Kjeldahl nitrogen concentrations decreased with increasing well depth and depth below water table. The relation between kjeldahl nitrogen concentrations and depth to water was poor. Nitrate and total phosphorus concentrations in water from wells were correlated among three hydrogeomorphic regions in the upper Snake River Basin, Concentrations of nitrate were statistically higher in the eastern Snake River Plain and local aquifers than in the tributary valleys. There was no statistical difference in total phosphorus concentrations among the three hydrogeomorphic regions. Nitrate and total phosphorus concentrations were correlated with land-use classifications developed using the Geographic Information Retrieval and Analysis System. Concentrations of nitrate were statistically higher in area of agricultural land than in areas of rangeland. There was no statistical difference in concentrations between rangeland and urban land and between urban land and agricultural land. There was no statistical difference in total phosphorus concentrations among any of the land-use classifications. Nitrate and total phosphorus concentrations also were correlated with land-use classifications developed by the Idaho Department of Water Resources for the Idaho part of the upper Snake River Basin. Nitrate concentrations were statistically higher in areas of irrigated agriculture than in areas of dryland agriculture and rangeland. There was no statistical difference in total phosphorus concentrations among any of the Idaho Department of Water Resources land-use classifications. Data were sufficient to assess long-term trends of nitrate concentrations in water from only eight wells: four wells north of Burley and four wells northwest of Pocatello. The trend in nitrate concentrations in water from all wells in upward. The following organic compounds were detected in ground water in the upper Snake River Basin: cyanazine, 2,4-D DDT, dacthal, diazinon, dichloropropane, dieldrin, malathion, and metribuzin. Of 211 wells sampled for organic compounds, water from 17 contained detectable concentrations.

Idaho;Montana;Nevada;Utah;Wyoming↗

Coproduction and modeling spatial contact networks prevent bias about infectious hematopoietic necrosis virus transmission for Snake River Basin salmonids

Much remains unknown about variation in pathogen transmission across the geographic range of a free-ranging fish or animal species and about the influence of movement (associated with husbandry practices or animal behavior) on pathogen transmission. Salmonid hatcheries are an ideal system in which to study these processes. Salmonid hatcheries are managed for endangered species recovery, supplementation of threatened or at-risk fish stocks, support of fisheries, and ecosystem stability. Infectious hematopoietic necrosis virus (IHNV) is a rhabdovirus of significant concern to salmon aquaculture. Landscape IHNV transmission dynamics previously had been estimated only for salmonid hatcheries in the Lower Columbia River Basin (LCRB). The objectives of this study were to estimate IHNV transmission dynamics in a unique geographic region, the Snake River Basin (SRB), and to quantitatively estimate the effect of model coproduction on inference because previous assessments of coproduction have been qualitative. In contrast to the LCRB, the SRB has hatchery complexes consisting of a main hatchery and ≥1 satellite facility. Knowledge about hatchery complexes was held by a subset of project researchers but would not have been available to project modelers without coproduction. Project modelers generated and tested multiple versions of Bayesian susceptible-exposedinfected models to realistically represent the SRB and estimate the effect of coproduction. Models estimated the frequency of transmission routes, route-specific infection probabilities, and infection probabilities for combinations of salmonid hosts and IHNV lineages. Model results indicated that in the SRB, avoiding exposure to IHNV-positive adult salmonids is the most important action to prevent juvenile infections. Migrating adult salmonids exposed juvenile cohort-sites most frequently, and the infection probability was greatest following exposure to migrating adults. Without coproduction, the frequency of exposure by migrating adults would have been overestimated by 70 cohort-sites, and the infection probability following exposure to migrating adults would have been underestimated by∼0.09. The coproduced model had less uncertainty in the infection probability if no transmission route could be identified (Bayesian credible interval (BCI) width = 0.12) compared to the model without coproduction (BCI width = 0.34). Evidence for virus lineage MD specialization on steelhead and rainbow trout (both Oncorhynchus mykiss) was apparent without model coproduction. In the SRB, we found a greater probability of virus lineage UC infection in Chinook salmon (Oncorhynchus tshawytscha) compared to in O. mykiss, whereas in the LCRB, UC more clearly exhibited a generalist approach. Coproduction influenced estimates that depended on transmission routes, which operated differently at main hatcheries and satellite sites within hatchery complexes. Hatchery complexes are found outside of the SRB and are not specific to salmonid hatcheries alone. There is great potential for coproduction and modeling spatial contact networks to advance understanding about infectious disease transmission in complex production systems and surrounding free-ranging animal populations.

Idaho, Oregon, Washington↗

Assessment of water-quality conditions in the J.B. Converse Lake watershed, Mobile County, Alabama, 1990-98

J.B. Converse (Converse) Lake is a 3,600-acre, tributary-storage reservoir in Mobile County, southwestern Alabama. The lake serves as the primary drinking-water supply for the city of Mobile. The Converse Lake watershed lies within the Coastal Plain Physiographic Province. Semiconsolidated to unconsolidated sediments of sand, silt, gravel, and clay underlie the watershed, and are covered by acidic soils. Land use in the watershed is mainly forest (64 percent) and agriculture (31 percent). Residential and commercial development account for only 1 percent of the total land use in the watershed. Converse Lake receives inflow from seven major tributaries. The greatest inflows are from Big Creek, Crooked Creek, and Hamilton Creek that had mean annual streamflows of 72.2, 19.4, and 25.0 cubic feet per second, respectively, for the period 1990 to 1998, which represents about 72 percent of the total annual streamflow to the lake. The total mean annual inflow to the lake is estimated to be about 163 cubic feet per second. In general, water quality in Converse Lake and its tributaries meets the criteria established by the Alabama Department of Environmental Management (ADEM) for drinking-water supplies, whole-body contact, and aquatic life. The exceptions include acidic pH levels, iron and manganese levels above secondary or aesthetic criteria, and fecal bacterial levels in some tributaries above whole-body contact (swimmable) criteria. The pH levels throughout the watershed were commonly below the criteria level of 6.0, but this appears to have been a naturally occurring phenomenon caused by poorly buffered soil types, resistant sediments, and forested land use. Median iron and manganese levels were above aesthetic criteria levels of 300 and 50 micrograms per liter, respectively, in some tributaries. All tributary sites in the Converse Lake watershed had median and minimum dissolved-oxygen concentrations above the ADEM criteria level of 5 milligrams per liter except for Boggy Branch, which had a minimum dissolved-oxygen concentration of 3.7 milligrams per liter. The degree to which nutrient contributions from tributaries were causing nutrient enrichment and eutrophication in Converse Lake was assessed. Trend analysis detected little or no change in nutrient concentrations at the tributary and lake sites in the Converse Lake watershed from the 1991 to 1998 water years. Nutrient concentrations at most tributary sites exhibited a significant, positive relation with streamflow that indicated the dominant source of nutrient input to the watershed is from nonpoint contributions. From 1990 to 1998, computed mean annual loads of 75,400 kilograms of total nitrogen, 36,950 kilograms of total Kjeldahl nitrogen, 28,870 kilograms of total inorganic nitrogen, and 3,480 kilograms of total phosphorus were contributed to the lake by Big Creek, Hamilton Creek, and Crooked Creek combined. These mean annual loads of nutrients corresponded to borderline eutrophic/mesotrophic conditions in the lake. Of the combined loads, 62 percent of the total nitrogen, 70 percent of the total Kjeldahl nitrogen, 54 percent of the total inorganic nitrogen, and 47 percent of the total phosphorus originated from the forested subbasin of Big Creek. The more residential and agricultural subbasins of Crooked Creek and Hamilton Creek, however, yielded over twice the total phosphorus load per hectare of land use. Crooked and Hamilton Creek subbasins also had higher yields of the more bioavailable total inorganic nitrogen. A simplistic empirical model could not explain the relation between year-to-year nutrient contributions to Converse Lake from the tributaries and the lake's ability to assimilate those contributions. The potential presence of pathogens in the lake and its tributaries was assessed based on fecal bacterial concentrations. Fecal bacterial concentrations at some tributary sites were above existing criteria for swimmable uses. Contributions of fecal bacte

Water-Resources Investigations Report↗

Statistical inference from capture data on closed animal populations

The estimation of animal abundance is an important problem in both the theoretical and applied biological sciences. Serious work to develop estimation methods began during the 1950s, with a few attempts before that time. The literature on estimation methods has increased tremendously during the past 25 years (Cormack 1968, Seber 1973). However, in large part, the problem remains unsolved. Past efforts toward comprehensive and systematic estimation of density (D) or population size (N) have been inadequate, in general. While more than 200 papers have been published on the subject, one is generally left without a unified approach to the estimation of abundance of an animal population This situation is unfortunate because a number of pressing research problems require such information. In addition, a wide array of environmental assessment studies and biological inventory programs require the estimation of animal abundance. These needs have been further emphasized by the requirement for the preparation of Environmental Impact Statements imposed by the National Environmental Protection Act in 1970. This publication treats inference procedures for certain types of capture data on closed animal populations. This includes multiple capture-recapture studies (variously called capture-mark-recapture, mark-recapture, or tag-recapture studies) involving livetrapping techniques and removal studies involving kill traps or at least temporary removal of captured individuals during the study. Animals do not necessarily need to be physically trapped; visual sightings of marked animals and electrofishing studies also produce data suitable for the methods described in this monograph. To provide a frame of reference for what follows, we give an exampled of a capture-recapture experiment to estimate population size of small animals using live traps. The general field experiment is similar for all capture-recapture studies (a removal study is, of course, slightly different). A typical field experiment is the following: a number of traps are positioned in the area to be studied, say 144 traps in a 12 X 12 grid, 7 m apart. At the beginning of the study (j=1) a sample size of n 1 is taken from the population, the animals are tagged and marked for future identification, and then returned to the population, usually at the same point where they were trapped. After allowing time of the marked and unmarked animals to mix, a second sample (j=2, often the following day) or n 2 animals is then taken.the second sample normally contains both marked and unmarked animals. The unmarked animals are marked and all captured animals are released back into the population. This procedure continues for t periods where t ≥ 2. The animals should be marked in such a way that the capture-recapture history of each animal caught during the study is known. In practice, toes are often clipped to uniquely identify individual animals (Taber and Cowan 1969) or serially numbered tags are sometimes used on larger animals. Such capture studies are classified by 2 schemes that are directly related to what class of models are appropriate and what parameters can be estimated. The first classification addresses the subject of closure. Closure usually means the size of the population is constant over the priod of investigation, i.e., no recruitment (birth or immigration) or losses (death or emigration). This is a strong assumption and, of course, never completely true in a natural biological population. For greater generality, we define closure to mean there are no unknown changes to the initial population. In practice, this means known losses (trap death), or deliberate removals) do not violate our definition of closure. If the study is properly designed, closure can be met at least approximately. Open or nonclosed populations explicitly allow for one or more types of recruitment or losses to operate during the course of the experiment (Jolly 1965, Seber 1965, Robson 1969, Pollock 1975). Only closed populations will be considered in this monograph. The second classification depends on the type of data collected with 2 possibilities occurring (Pollock 1974, unpublished doctoral dissertation, Cornell University, Ithaca, New York): (1) only information on the recovery of marked animals is available for each sampling occasion, j, j=1, 2, ... t. (2) information on both marked and unmarked animals is available for each sampling occasion, j, j=1, 2, ... t. In case (1), population size (N) is not identifiable, however, other parameters can be estimated (Brownie et al. 1978). In case (2), N can be estimated using a wide variety of approaches depending upon what we wish to assume. Only case (2) will be dealt with here.

Wildlife Monographs↗

Quantitative microbial risk assessment for ingestion of antibiotic resistance genes from private wells contaminated by human and livestock fecal sources

We used quantitative microbial risk assessment to estimate ingestion risk for intI1 , erm (B), sul1 , tet (A), tet (W), and tet (X) in private wells contaminated by human and/or livestock feces. Genes were quantified with five human-specific and six bovine-specific microbial source-tracking (MST) markers in 138 well-water samples from a rural Wisconsin county. Daily ingestion risk (probability of swallowing ≥1 gene) was based on daily water consumption and a Poisson exposure model. Calculations were stratified by MST source and soil depth over the aquifer where wells were drilled. Relative ingestion risk was estimated using wells with no MST detections and >6.1 m soil depth as a referent category. Daily ingestion risk varied from 0 to 8.8 × 10 −1 by gene and fecal source (i.e., human or bovine). The estimated number of residents ingesting target genes from private wells varied from 910 ( tet (A)) to 1,500 ( intI1 and tet (X)) per day out of 12,000 total. Relative risk of tet (A) ingestion was significantly higher in wells with MST markers detected, including wells with ≤6.1 m soil depth contaminated by bovine markers (2.2 [90% CI: 1.1–4.7]), wells with >6.1 m soil depth contaminated by bovine markers (1.8 [1.002–3.9]), and wells with ≤6.1 m soil depth contaminated by bovine and human markers simultaneously (3.1 [1.7–6.5]). Antibiotic resistance genes (ARGs) were not necessarily present in viable microorganisms, and ingestion is not directly associated with infection. However, results illustrate relative contributions of human and livestock fecal sources to ARG exposure and highlight rural groundwater as a significant point of exposure.

Wisconsin↗

The quality of our nation's waters: Nutrients in the nation's streams and groundwater, 1992-2004

National Findings and Their Implications Although the use of artificial fertilizer has supported increasing food production to meet the needs of a growing population, increases in nutrient loadings from agricultural and, to a lesser extent, urban sources have resulted in nutrient concentrations in many streams and parts of aquifers that exceed standards for protection of human health and (or) aquatic life, often by large margins. Do NAWQA findings substantiate national concerns for aquatic and human health? National Water-Quality Assessment (NAWQA) findings indicate that nutrient concentrations in streams and groundwater in basins with significant agricultural or urban development are substantially greater than naturally occurring or “background” levels. For example, median concentrations of total nitrogen and phosphorus in agricultural streams are about 6 times greater than background levels. Findings also indicate that concentrations in streams routinely were 2 to 10 times greater than regional nutrient criteria recommended by the U.S. Environmental Protection Agency (USEPA) to protect aquatic life. Such large differences in magnitude suggest that significant reductions in sources of nutrients, as well as greater use of land management strategies to reduce the transport of nutrients to streams, are needed to meet recommended criteria for streams draining areas with significant agricultural and urban development. Nitrate concentrations above the Federal drinking-water standard—or Maximum Contaminant Level (MCL)—of 10 milligrams per liter (mg/L, as nitrogen) are relatively uncommon in samples from streams used for drinking water or from relatively deep aquifers; the MCL is exceeded, however, in more than 20 percent of shallow (less than 100 feet below the water table) domestic wells in agricultural areas. This finding raises concerns for human health in rural agricultural areas where shallow groundwater is used for domestic supply and may warn of future contamination of deeper groundwater pumped from public‑supply wells. Are levels of nutrients in water increasing or decreasing? A decadal assessment of trends in concentrations of nitrogen and phosphorus from about 1993 to 2003 shows minimal changes in those concentrations in the majority of studied streams across the Nation, and more upward than downward trends in concentrations at sites with changes. These findings underscore the need for reductions in nutrient inputs or management strategies that would reduce transport of nutrients to streams. Upward trends were evident among all land uses, including those only minimally affected by agricultural and (or) urban development, which suggests that additional protection of some of our Nation’s most pristine streams warrants consideration. The median of nitrate concentrations in groundwater from 495 wells also increased significantly from 3.2 to 3.4 mg/L (6 percent) during about the same period, and the proportion of wells with concentrations of nitrate greater than the MCL increased from 16 to 21 percent. Nitrate concentrations in water in deep aquifers are likely to increase during the next decade as shallow groundwater with elevated concentrations moves downward. The potential for future contamination of the deep aquifers requires attention because these aquifers commonly are used for public water supply, and because restoration of groundwater is costly and difficult. Long-term and consistent monitoring of nutrients, improved accounting of nutrient sources, and improved tracking and modeling of climatic and landscape changes will be essential for distinguishing trends in nutrient concentrations, understanding the causes of those trends, and accurately tracking the effectiveness of strategies implemented to manage nutrients.

Circular↗

Avian predation of juvenile Lost River and Shortnose Suckers in Upper Klamath Lake: An assessment of Sucker assisted rearing program releases during 2018–2020

To bolster recruitment in Endangered Species Act (ESA) listed Lost River Suckers (Deltistes luxatus) and Shortnose Suckers (Chasmistes brevirostris) in the Upper Klamath Basin (UKB), the U.S. Fish and Wildlife Service (USFWS) and its partners have implemented the Sucker Assisted Rearing Program (SARP). As part of this program, juvenile suckers were reared in captivity, implanted with passive integrated transponder (PIT) tags (n= 8,857), and released into the Upper Klamath Lake or its tributaries during 2018–2020. Previous research suggests that predation by American White Pelicans (Pelecanus erythrorhynchos), Double-crested Cormorants (Nannopterum auritum), and Caspian Terns (Hydroprogne caspia) may negatively influence sucker survival, particularly predation on juvenile suckers. Estimates of predation impacts from past studies, however, represented minimum estimates of sucker mortality because analyses did not account for the proportion of consumed tags that were deposited by birds on their breeding colony where PIT tag recovery efforts took place. To estimate and account for deposition probabilities, we conducted a field study in which we fed pelicans PIT-tagged juvenile suckers (n = 401). We accounted for deposition probabilities of cormorants and terns by using previously published estimates. Sucker PIT tags were recovered from pelican, cormorant, and tern nesting sites in the UKB following each breeding season and a hierarchical Bayesian model was used to estimate predation rates (percentage of available tagged fish consumed) on SARP releases as well as naturally-reared or wild juvenile suckers and adult suckers that were PIT-tagged in Upper Klamath Lake and Clear Lake Reservoir. Pelican deposition probabilities were estimated at 0.47 (95% credible interval = 0.36–0.60), indicating that for every 100 PIT tags consumed, on average, 47 were deposited by pelicans on breeding colonies. Estimates of predation rates that incorporate corrections for deposition on SARP releases ranged annually from 4.4% (95% credible interval = 2.9–6.8%) to 8.8% (6.2–13.3%) during 2018–2020. Results suggest that colonial waterbird predation impacts on SARP releases likely constituted a small, but unknown, component of total mortality for suckers released into the Upper Klamath Lake system. Predation impacts on SARP juvenile suckers and wild juvenile suckers, which were estimated annually at 4.7% (1.0–13.9%) to 14.9% (7.6–29.3%), were consistently higher than those observed on adult suckers, with predation on adult suckers typically less than 4.0% of available fish annually. Future predation studies may consider models that integrate both live and dead detections of PIT-tagged suckers to generate more accurate and precise estimates of survival following release, as well as models that consider environmental factors that influence sucker susceptibility to colonial waterbird predation. Such models would provide a more holistic understanding of the degree to which avian predation limits the survival of ESA-listed suckers in the UKB.

Oregon↗

Assessment of dissolved-selenium concentrations and loads in the Lower Gunnison River Basin, Colorado, as part of the Selenium Management Program, 2011–17

The Gunnison Basin Selenium Management Program implemented a water-quality monitoring network in 2011 to measure concentrations of selenium in the lower Gunnison River Basin in Colorado. Selenium is a trace element that bioaccumulates in aquatic food chains. Selenium is essential for life, but elevated amounts can cause reproductive failure, deformities, and other harmful effects. The primary goal of the Selenium Management Program is to meet the State of Colorado water-quality standard of 4.6 micrograms per liter (µg/L) for dissolved selenium at the U.S. Geological Survey (USGS) streamflow-gaging station number 09152500—Gunnison River near Grand Junction, Colorado—herein referred to as “Whitewater.” The U.S. Geological Survey, in cooperation with the Bureau of Reclamation, has completed a review of dissolved-selenium data collected from the Selenium Management Program network during Water Year (WY) 2017 (October 1, 2016 through September 30, 2017) to further the understanding of the status and trends of selenium in the basin. This report presents the percentile values for selenium because regulatory agencies in Colorado make decisions based on the U.S. Environmental Protection Agency’s Clean Water Act section 303(d), which uses percentile values for concentrations. Also presented are dissolved-selenium loads at 14 sites in the lower Gunnison River Basin for WYs 2011–17. Annual dissolved-selenium loads were calculated for six sites with continuous U.S. Geological Survey streamflow-gaging stations. These six sites are referred to as “core” sites in this report. The remaining sites, which do not have streamflow-gaging stations, are referred to as “ancillary” sites in this report. During WY 2017, the loads calculated at the six core sites ranged from 306 pounds (lb) at Uncompahgre River at Colona to 12,600 lb at Whitewater, respectively. By using discrete water-quality samples and the associated discharge measurements, instantaneous loads were calculated for 14 sites in WYs 2011–17 where discrete water-quality sampling took place. Median instantaneous loads ranged from 0.52 pounds per day (lb/d) at Uncompahgre River at Colona to 35.7 lb/d at Whitewater. Mean instantaneous loads ranged from 0.63 lb/d at Cummings Gulch at mouth to 35.5 lb/d at Whitewater. Most tributary sites in the basin had a median instantaneous dissolved-selenium load of less than 20.0 lb/d. In general, dissolved-selenium loads at Gunnison River main-stem sites showed an increase from upstream to downstream. The State of Colorado’s water-quality standard for dissolved selenium of 4.6 µg/L was compared to the 85th percentiles for dissolved selenium at selected sites. Annual 85th percentiles for dissolved selenium were calculated by using estimated dissolved-selenium concentrations from linear regression models for the six core sites with U.S. Geological Survey streamflow-gaging stations. The 85th-percentile concentrations for WY 2017 based on this method ranged from 0.68 µg/L at Uncompahgre River at Colona to 140 µg/L at Loutzenhizer Arroyo at North River Road. The 85th percentiles for concentrations of dissolved selenium also were calculated from water-quality samples collected during WY 2017 from sites with sufficient data. The annual 85th-percentile concentrations based on the discrete samples ranged from 0.75 µg/L at Uncompahgre River at Colona to 106 µg/L at Loutzenhizer Arroyo at North River Road. An analysis was completed for Whitewater to determine if an upward or downward trend exists for dissolved-selenium loads during two time periods. The first time period included all data at Whitewater, whereas the second time period focused on more recent data. The trend analysis indicates a decrease from 22,200 to 12,600 lb, which is a 43.1 percent (9,600 lb) reduction during the time period WY 1986 through WY 2017. The trend analysis for the annual dissolved-selenium load for WY 1995 through WY 2017 indicates a decrease of 6,600 lb per year, or 35.5 percent. An evaluation of laboratory bias was completed for selenium data which was used in the trend analysis. Findings indicated a potential positive bias of approximately 12 percent may exist in the data from October 2005 through August 2015.

Colorado↗