USGS Science⌕ Search

SEARCH · USGS Science

Results for “Decision Analysis”

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,027 records · Page 57Linked to original sources
USGS2007

A comparison of macroinvertebrate and habitat methods of data collection in the Little Colorado River Watershed, Arizona 2007

The Arizona Department of Environmental Quality (ADEQ) and the U.S. Environmental Protection Agency (USEPA) Ecological Monitoring and Assessment Program (EMAP), use different field methods for collecting macroinvertebrate samples and habitat data for bioassessment purposes. Arizona’s Biocriteria index was developed using a riffle habitat sampling methodology, whereas the EMAP method employs a multi-habitat sampling protocol. There was a need to demonstrate comparability of these different bioassessment methodologies to allow use of the EMAP multi-habitat protocol for both statewide probabilistic assessments for integration of the EMAP data into the national (305b) assessment and for targeted in-state bioassessments for 303d determinations of standards violations and impaired aquatic life conditions. The purpose of this study was to evaluate whether the two methods yield similar bioassessment results, such that the data could be used interchangeably in water quality assessments. In this Regional EMAP grant funded project, a probabilistic survey of 30 sites in the Little Colorado River basin was conducted in the spring of 2007. Macroinvertebrate and habitat data were collected using both ADEQ and EMAP sampling methods, from adjacent reaches within these stream channels. All analyses indicated that the two macroinvertebrate sampling methods were significantly correlated. ADEQ and EMAP samples were classified into the same scoring categories (meeting, inconclusive, violating the biocriteria standard) 82% of the time. When the ADEQ-IBI was applied to both the ADEQ and EMAP taxa lists, the resulting IBI scores were significantly correlated (r=0.91), even though only 4 of the 7 metrics in the IBI were significantly correlated. The IBI scores from both methods were significantly correlated to the percent of riffle habitat, even though the average percent riffle habitat was only 30% of the stream reach. Multivariate analyses found that the percent riffle was an important attribute for both datasets in classifying IBI scores into assessment categories. Habitat measurements generated from EMAP and ADEQ methods were also significantly correlated; 13 of 16 habitat measures were significantly correlated (p<0.01). The visual-based percentage estimates of percent riffle and pool habitats, vegetative cover and percent canopy cover, and substrate measurements of percent fine substrate and embeddedness were all remarkably similar, given the different field methods used. A multivariate analysis identified substrate and flow conditions, as well as canopy cover as important combinations of habitat attributes affecting both IBI scores. These results indicate that similar habitat measures can be obtained using two different field sampling protocols. In addition, similar combinations of these habitat parameters were important to macroinvertebrate community condition in multivariate analyses of both ADEQ and EMAP datasets. These results indicate the two sampling methods for macroinvertebrates and habitat data were very similar in terms of bioassessment results and stressors. While the bioassessment category was not identical for all sites, overall the assessments were significantly correlated, providing similar bioassessment results for the cold water streams used in this study. The findings of this study indicate that ADEQ can utilize either a riffle-based sampling methodology or a multi-habitat sampling approach in cold water streams as both yield similar results relative to the macroinvertebrate assemblage. These results will allow for use of either macroinvertebrate dataset to determine water quality standards compliance with the ADEQ Indexes of Biological Integrity, for which threshold values were just recently placed into the Arizona Surface Water Quality Standards. While this survey did not include warm water desert streams of Arizona, we would predict that EMAP and ADEQ sampling methodologies would provide similar bioassessment results and would not be significantly different, as we have found that the percent riffle habitat in cold and warm water perennial, wadeable streams is not significantly different. However, a comparison study of sampling methodologies in warm water streams should be conducted to confirm the predicted similarity of bioassessment results. ADEQ will continue to implement a monitoring strategy that includes probabilistic monitoring for a statewide ecological assessment of stream conditions. Conclusions from this study will guide decisions regarding the most appropriate sampling methods for future probabilistic monitoring sample plans.

Arizona↗
USGS1993

Economic analysis of critical habitat designation for the desert tortoise (Mojave population)

The U.S. Fish and Wildlife Service emergency 1isted the Mojave population of the desert tortoise as endangered on August 4, 1989. The Mojave population formally was listed as threatened on April 2, 1990. The Endangered Species Act of 1973, as amended, requires that the economic benefits and costs and other relevant effects of critical habitat designation be considered. The Secretary of the Interior may exclude from designation areas where the costs of designation are greater than the benefits, unless the exclusion would result in extinction of the species. Desert tortoises are threatened by an accumulation of human-and disease-related mortality accompanied by habitat destruction, degradation, and fragmentation. Many desert tortoises are illegally collected for pets, food, and commercial trade. Others are accidentally struck and killed by vehicles on roads and highways or are killed by gunshot or vehicles traveling off-highway. Raven predation on hatchling desert tortoises has increased as raven populations in the desert have risen. An upper respiratory tract disease is suspected to be a major cause of mortality in the western Mojave Desert. This presumably incurable affliction presumably is thought to be spread through the release of infected tortoises into the desert. The Service has proposed designating critical habitat in nine counties within four states. The 12 critical habitat units encompass 6.4 million acres of land, more than 80% federally owned. This region is economically and demographically diverse. Most of the land is sparsely settled and characterized as a hot desert ecosystem. Major industries in the region include entertainment and lodging (primarily in Las Vegas), property development to accommodate the rapid population growth, and services. Millions of rural acres in the region are leased by the federal government for livestock grazing and used for mining. Overall economic benefits to the affected states derived from cattle and sheep grazing in the hot desert areas are minimal and, according to a recent U.S. General Accounting Office study (1991), local economies do not depend on the grazing of public lands for economic survival. The economic analysis describes the economy in 1990, prior to designation, and estimates the effects of designation. The report estimates those incremental effects attributable to critical habitat designation. Impacts attributable to listing the species were not considered in this analysis. Although critical habitat units have been designated in nine counties, two counties are omitted from the economic analysis because of the small proportion of critical habitat acreage they include. Three key activities (cattle grazing, mineral extraction, and off-road vehicle use) were studied in detail. Even if livestock grazing and commercial off-road racing events were eliminated in the proposed critical habitat units, the potential incremental regional economic impacts would be extremely small. The findings in the report do not include the assumption that mining would be eliminated from critical habitat units, but rather that consultation may result in added mitigation and/or relocation of features. Studies show that society will realize benefits from preservation of species and ecosystems. Survey-based studies confirm that benefits exist and are substantial in size, although these benefits often are not measured in traditional economic markets. The total benefit to society of desert tortoise preservation includes several components. Biodiversity in the Mojave and Colorado Deserts will be improved, some recreation values may increase, and gains in intrinsic value will be realized. Critical habitat designation should result in the loss of fewer than 425 total jobs in the seven counties. These include 340 direct ranching jobs and 85 indirect jobs in other industries. This estimated employment loss will not be permanent for most laborers, it is anticipated that over 85% will be reemployed within two years. The economic consequences of designating critical habitat includes reduced ranch profits in the seven counties of $4,470,000 [the estimated permanent decrease in ranch profits, capitalized at 10% for a 50-year period, in accordance with the methodology of River, et al. (1978)]. Second, the federal government will compensate allottees with a one-time payment estimated at $376,000 for the loss of permanent improvements to grazing lands (pending BLM administrative decisions of partially affected allotments). Finally, discontinuing grazing fees will result in an annual reduction of $170,000 in collected grazing fees that are divided among range improvements, the Federal Treasury, and local governments.

California, Nevada↗
USGS2003

Vulnerability of ground water to contamination, northern Bexar County, Texas

The Trinity aquifer, composed of Lower Cretaceous carbonate rocks, largely controls the ground-water hydrology in the study area of northern Bexar County, Texas. Discharge from the Trinity aquifer recharges the downgradient, hydraulically connected Edwards aquifer one of the most permeable and productive aquifers in the Nation and the sole source of water for more than a million people in south-central Texas. The unconfined, karstic outcrop of the Edwards aquifer makes it particularly vulnerable to contamination resulting from urbanization that is spreading rapidly northward across an "environmentally sensitive" recharge zone of the Edwards aquifer and its upgradient "catchment area," composed mostly of the less permeable Trinity aquifer. A better understanding of the Trinity aquifer is needed to evaluate water-management decisions affecting the quality of water in both the Trinity and Edwards aquifers. A study was made, therefore, in cooperation with the San Antonio Water System to assess northern Bexar County's vulnerability to ground-water contamination. The vulnerability of ground water to contamination in this area varies with the effects of five categories of natural features (hydrogeologic units, faults, caves and (or) sinkholes, slopes, and soils) that occur on the outcrop and in the shallow subcrop of the Glen Rose Limestone. Where faults affect the rates of recharge or discharge or the patterns of ground-water flow in the Glen Rose Limestone, they likewise affect the risk of water-quality degradation. Caves and sinkholes generally increase the vulnerability of ground water to contamination, especially where their occurrences are concentrated. The slope of land surface can affect the vulnerability of ground water by controlling where and how long a potential contaminant remains on the surface. Disregarding the exception of steep slopes which are assumed to have no soil cover the greater the slope, the less the risk of ground-water contamination. Because most soils in the study area are uniformly thin, they have only minimal effect on the vulnerability of ground water to contamination. The results of hydrogeologic mapping during the present study divide the outcrop of the Glen Rose Limestone into five mappable intervals, labeled (youngest to oldest) A through E. Of these intervals, only the middle (C) and the lowermost (E) generally provide appreciable permeability. The vulnerability assessment provided herein was determined by combining the presumed effects of selected natural features (with individual vulnerability ratings ranging from 0 through 35) using a grid-based, multilayer system of digital datasets and geographic information system analysis. The resulting vulnerability map comprises composite vulnerability ratings that range from 26 through 104. The relatively less vulnerable areas those containing no faults, sinkholes, or caves occupy about 92 percent of the study area. The most vulnerable areas are those containing both a fault and one or more caves. The distribution of the most vulnerable areas which trend from southwest to northeast, roughly parallel to the Balcones fault zone occur mainly where faults intersect caves.

Texas↗
USGS2007

Effects of wastewater disinfection on waterborne bacteria and viruses

Wastewater disinfection is practiced with the goal of reducing risks of human exposure to pathogenic microorganisms. In most circumstances, the efficacy of a wastewater disinfection process is regulated and monitored based on measurements of the responses of indicator bacteria. However, inactivation of indicator bacteria does not guarantee an acceptable degree of inactivation among other waterborne microorganisms (e.g., microbial pathogens). Undisinfected effluent samples from several municipal wastewater treatment facilities were collected for analysis. Facilities were selected to provide a broad spectrum of effluent quality, particularly as related to nitrogenous compounds. Samples were subjected to bench-scale chlorination and dechlorination and UV irradiation under conditions that allowed compliance with relevant discharge regulations and such that disinfectant exposures could be accurately quantified. Disinfected samples were subjected to a battery of assays to assess the immediate and long-term effects of wastewater disinfection on waterborne bacteria and viruses. In general, (viable) bacterial populations showed an immediate decline as a result of disinfectant exposure; however, incubation of disinfected samples under conditions that were designed to mimic the conditions in a receiving stream resulted in substantial recovery of the total bacterial community. The bacterial groups that are commonly used as indicators do not provide an accurate representation of the response of the bacterial community to disinfectant exposure and subsequent recovery in the environment. UV irradiation and chlorination/dechlorination both accomplished measurable inactivation of indigenous phage; however, the extent of inactivation was fairly modest under the conditions of disinfection used in this study. UV irradiation was consistently more effective as a virucide than chlorination/dechlorination under the conditions of application, based on measurements of virus (phage) diversity and concentration. Taken together, and when considered in conjunction with previously published research, the results of these experiments illustrate several important limitations of common disinfection processes as applied in the treatment of municipal wastewaters. In general, it is not clear that conventional disinfection processes, as commonly implemented, are effective for control of the risks of disease transmission, particularly those associated with viral pathogens. Microbial quality in receiving streams may not be substantially improved by the application of these disinfection processes; under some circumstances, an argument can be made that disinfection may actually yield a decrease in effluent and receiving water quality. Decisions regarding the need for effluent disinfection must account for site-specific characteristics, but it is not clear that disinfection of municipal wastewater effluents is necessary or beneficial for all facilities. When direct human contact or ingestion of municipal wastewater effluents is likely, disinfection may be necessary. Under these circumstances, UV irradiation appears to be superior to chlorination in terms of microbial quality and chemistry and toxicology. This advantage is particularly evident in effluents that contain appreciable quantities of ammonia-nitrogen or organic nitrogen.

Water Environment Research↗
USGS2023

Potentiometric surface map of the Southern High Plains aquifer in the Cannon Air Force Base area, Curry County, New Mexico, 2020

Declining water levels and the potential impact on water resources on and around Cannon Air Force Base (AFB), New Mexico, has necessitated an up-to-date review of the potentiometric surface to evaluate the availability of water resources for future use. Analysis of groundwater-flow directions and hydraulic gradients can provide an understanding of depletion by heavy groundwater pumping and recharge through playa lakes, as well as the relationship between the groundwater levels and underlying geology. The objectives of this study, conducted by the U.S. Geological Survey in cooperation with the U.S. Air Force Civil Engineer Center, are to assist Cannon AFB in understanding and interpreting current and local hydrologic conditions and to evaluate groundwater-level change from 2015 to 2020 using new and historical data. A groundwater potentiometric surface contour map was constructed to better understand the Southern High Plains aquifer around Cannon AFB and to show the altitude of the water-table surface and groundwater-flow patterns. Four hydrographs were created from periodic measurements of groundwater levels in wells on and around Cannon AFB to provide information about historical groundwater-level changes and visualize trends from the water-level records. The long-term trend present in all four hydrographs is a steady decline in groundwater levels, with some areas declining faster than others. The groundwater-level change map presented in this study provides a visual representation of the change in groundwater level from the winter 2015 to winter 2020 measuring events. Results show that among corresponding wells measured in 2015 and 2020, 50.7 percent indicated a decline in water levels, 29.9 percent indicated neutral water levels, and 19.4 percent indicated a rise in water levels. The region to the north of the groundwater trough on Cannon AFB contained most of the groundwater-level rises, whereas the regions located near the trough and just west of Clovis, N. Mex., contained most of the declines. These results suggest that continued monitoring of declining groundwater levels in the area would provide valuable decision-support information for assessing the sustainability of this water resource.

New Mexico↗
USGS2013

Selenium in ecosystems within the mountaintop coal mining and valley-fill region of southern West Virginia-assessment and ecosystem-scale modeling

Coal and associated waste rock are among environmental selenium (Se) sources that have the potential to affect reproduction in fish and aquatic birds. Ecosystems of southern West Virginia that are affected by drainage from mountaintop coal mines and valleys filled with waste rock in the Coal, Gauley, and Lower Guyandotte watersheds were assessed during 2010 and 2011. Sampling data from earlier studies in these watersheds (for example, Upper Mud River Reservoir) and other mining-affected watersheds also are included to assess additional hydrologic settings and food webs for comparison. Basin schematics give a comprehensive view of sampled species and Se concentration data specific to location and date. Food-web diagrams document the progression of Se trophic transfer across suspended particulate material, invertebrates, and fish for each site to serve as the basis for developing an ecosystem-scale model to predict Se exposure within the hydrologic conditions and food webs of southern West Virginia. This approach integrates a site-specific predator’s dietary exposure pathway into modeling to ensure an adequate link to Se toxicity and, thus, to species vulnerability. Site-specific fish abundance and richness data in streams documented various species of chub, shiner, dace, darters, bass, minnow, sunfish, sucker, catfish, and central stoneroller (Campostoma anomalum), mottled sculpin (Cottus bairdii), and least brook lamprey (Lampetra aepyptera). However, Se assessment species for streams, and hence, model species for streams, were limited to creek chub (Semotilus atromaculatus) and central stoneroller. Both of these species of fish are generally considered to have a high tolerance for environmental stress based on traditional comparative fish community assessment, with creek chub being present at all sites. Aquatic insects (mayfly, caddisfly, stonefly, dobsonfly, chironomid) were the main invertebrates sampled in streams. Collection of suspended particulate material acted as an integrator of organic-rich, fine-grained biomass present in streams. The base-case food web modeled for streams was suspended particulate material to aquatic insect to creek chub, with comparative modeling of a direct particulate-to-stoneroller food web. Model species for a reservoir setting were based on an earlier study of bluegill sunfish (Lepomis macrochirus), green sunfish (Lepomis cyanellus), and largemouth bass (Micropterus salmoides). Several reservoir food webs were considered based on a variety of invertebrates (insect, snail, clam). For stream and reservoir settings, predicted Se concentrations in exposure scenarios showed a high degree of correlation (r 2 = 0.91 for invertebrates and 0.75 for fish) with field observations of Se concentrations when modeling was initiated from suspended-particulate-material Se concentrations and model transfer parameters defined previously in the literature were used. These strong correlations validate the derived site-specific model and establish sufficient confidence that the predictions from the developed model can be quantitatively applied to the ecosystems in southern West Virginia. An application of modeling used a metric describing the partitioning of Se between particulate material and dissolved phases (K d ) to allow determination of a dissolved Se concentration that would be necessary to attain a site-specific Se fish body burden. The operationally defined K d quantifies the complex process of transformation at the base of a food web on a site-specific basis. The magnitude of this metric is known to vary with such factors as Se speciation, particulate-material type, and hydrology. This application (1) ties dissolved Se concentrations to fish tissue concentrations; (2) allows consideration of different choices for intervening site-specific exposure steps that set Se bioaccumulation, partitioning, and bioavailability; and (3) generates implications for management decisions that define protection through different regulatory pathways and guidelines. The range of model outcomes accounts for critical sources of variability and establishes whether site and food-web characterization were adequate to represent the dynamics of the system with certainty. This is especially true in terms of particulate-material phases at the base of the food web and utilization of K d in different hydrologic settings. For streams, a range of field-derived K d ds were applied to food-web exposure scenarios within a framework of locational and hydrologic variables (area of stream basin; stream gradient and discharge) that may affect the magnitude of K d . Overlaying even a coarse temporal scale that acknowledges variability in stream dissolved Se and Se speciation, such as through seasonal derivation of K d , can substantially narrow model uncertainty. Modeling that constrains the place and time of greatest ecosystem Se sensitivity within a specified food web gives insight into Se risk and identifies controlling management alternatives within a watershed or stream basin. If there is a range of hydrologic settings, specificity is needed to establish a hierarchy of in-stream and off-stream habitats for a watershed approach that takes into account Se-enriched water moving through different K d and food web environments. If there is a range of predator vulnerabilities (measured as a combination of food-web Se biodynamics and response in Se toxicity tests) within the site-specific community of fish species to be protected, then choice of fish species is critical to protection because it determines the food web and, hence, the magnitude of biotransfer through which Se is modeled. Whether creek chub is representative of the vulnerability to Se of all fish species encountered within the study-site ecosystems will require additional species-specific data and analysis. A range of site-specific scenarios illustrated here set model outcomes, but the final quantitative evaluation of alternatives and their implications will be those generated through choices and guidance formulated by state and other agencies in their decisionmaking processes. Proposed additions and refinements to the ecosystem-scale site-specific approach developed here include consideration of: measurement of temporally matched pairs of dissolved and suspended-particulate-material Se concentrations across a broader range of stream sites to expand the stream K d database and to test the representativeness of a suspended-particulate-material sample within a stream; characterization of different phases of particulate material across seasons to better define the base of the food web and connect to invertebrate feeding; refinement of model assumptions concerning dietary preferences and composition for fish to develop additional trophic transfer factors (TTFs) (for example, calculation of TTFinvertebrate composite for mixed diets); expansion of modeling of fish species and their food webs to include Se-vulnerable species; temporal characterization of a predator’s life cycle and habitat use as additional model layers to integrate with Se biodynamics in streams; investigation of the effect of stream gradient on K d based on a finer scale than presented here in terms of such variables as residence time, watershed dilution, and physical habitat attributes (for example, amount of ponding versus run or riffle within a stream); and linkage to discharge through use of stream gaging to record variability and enable model organization within water-year types and discharge seasons. Investigating the presence and variability of prey and predator species in demographically open systems such as streams also is key to model outcomes given the overall environmental stressors (for example, general landscape change, food-web disruption, recolonization potential) imposed on the composition of biological communities in coal mining and valley-fill affected watersheds

West Virginia↗
USGS2025

Estimated annual abundance of migratory Peale's Peregrine Falcons in coastal Washington, USA

Following the recovery of Peregrine Falcons ( Falco peregrinus ), the US Fish and Wildlife Service began a process to allow “take” (capture) of wild peregrines for falconry in the United States. Recently, that effort involved generating updated estimates of the collective abundance of the three North American peregrine subspecies: F. p. anatum , F. p. tundrius , and F. p. pealei (Peale's Peregrine Falcon). Because of the more limited distribution of F. p. pealei , we conducted an analysis specific to its geographic range. We analyzed data from a long-term banding and resighting program on three beaches on the southern coast of Washington, USA, to estimate the annual abundance of migrating and overwintering F. p. pealei , using the capture histories of 250 Peregrine Falcons, nearly all of which were captured during 1277 vehicle surveys between 1995 and 2024. Because we studied an open population of migratory individuals, we used a zero-inflated Poisson log-normal mark-resight model to estimate annual abundance. For the analyses, we partitioned our survey data into sighting periods, each of which extended from 1 September of one year to 31 May of the next. We anticipated that first-year F. p. pealei would be identified for falconry take, and our annual abundance estimates for first-year birds of this subspecies ranged from a high of 24.8 ± 6.1 (SE) individuals in the 2014–2015 sighting period to a low of 1.9 ± 1.4 individuals in the 2023–2024 sighting period. Peregrine Falcon abundance varied annually and appeared to decline during the last two sighting periods. Our sighting rate of marked peregrines was negatively associated with Bald Eagle ( Haliaeetus leucocephalus ) encounter rate. There was a lesser relationship to human activity, and we suspect the change in sighting rate was a behavioral response by Peregrine Falcons to the threat of kleptoparasitism by Bald Eagles. We currently lack comprehensive information about the natal origin of the individual peregrines in our study area, which prevented us from assessing the degree to which falconry take from the pool of falcons migrating to or through Washington might potentially impact local or regional abundances. Although a better understanding of natal origins is needed, our data add clarity to the migration and overwinter abundance of F. p. pealei on the Washington coast and may inform decisions about the take of this subspecies for falconry.

Washington↗
USGS2003

Effectiveness of the New Hampshire stream-gaging network in providing regional streamflow information

The stream-gaging network in New Hampshire was analyzed for its effectiveness in providing regional information on peak-flood flow, mean-flow, and low-flow frequency. The data available for analysis were from stream-gaging stations in New Hampshire and selected stations in adjacent States. The principles of generalized-least-squares regression analysis were applied to develop regional regression equations that relate streamflow-frequency characteristics to watershed characteristics. Regression equations were developed for (1) the instantaneous peak flow with a 100-year recurrence interval, (2) the mean-annual flow, and (3) the 7-day, 10-year low flow. Active and discontinued stream-gaging stations with 10 or more years of flow data were used to develop the regression equations. Each stream-gaging station in the network was evaluated and ranked on the basis of how much the data from that station contributed to the cost-weighted sampling-error component of the regression equation. The potential effect of data from proposed and new stream-gaging stations on the sampling error also was evaluated. The stream-gaging network was evaluated for conditions in water year 2000 and for estimated conditions under various network strategies if an additional 5 years and 20 years of streamflow data were collected. The effectiveness of the stream-gaging network in providing regional streamflow information could be improved for all three flow characteristics with the collection of additional flow data, both temporally and spatially. With additional years of data collection, the greatest reduction in the average sampling error of the regional regression equations was found for the peak- and low-flow characteristics. In general, additional data collection at stream-gaging stations with unregulated flow, relatively short-term record (less than 20 years), and drainage areas smaller than 45 square miles contributed the largest cost-weighted reduction to the average sampling error of the regional estimating equations. The results of the network analyses can be used to prioritize the continued operation of active stations, the reactivation of discontinued stations, or the activation of new stations to maximize the regional information content provided by the stream-gaging network. Final decisions regarding altering the New Hampshire stream-gaging network would require the consideration of the many uses of the streamflow data serving local, State, and Federal interests.

Water-Resources Investigations Report↗
USGS2024

A protocol for assessing bias and robustness of social network metrics using GPS based radio-telemetry data

Background Social network analysis of animal societies allows scientists to test hypotheses about social evolution, behaviour, and dynamic processes. However, the accuracy of estimated metrics depends on data characteristics like sample proportion, sample size, and frequency. A protocol is urgently needed to assess for bias and robustness of social network metrics estimated for the animal populations especially when a limited number of individuals are monitored. Methods We used GPS telemetry datasets of five ungulate species to combine known social network approaches with novel ones into a comprehensive five-step protocol. To quantify the bias and uncertainty in the network metrics obtained from a partial population, we presented novel statistical methods which are particularly suited for autocorrelated data, such as telemetry relocations. The protocol was validated using a sixth species, the fallow deer, with a known population size where ⇠ 85% of the individuals have been directly monitored. Results Through the protocol, we demonstrated how pre-network data permu tations allow researchers to assess non-random aspects of interactions within a population. The protocol assesses bias in global network metrics, obtains confidence intervals, and quantifies uncertainty of global and node-level network metrics based on the number of nodes in the network. We found that global network metrics like density remained robust even with a lowered sample size, while local network metrics like eigenvector centrality were unreliable for four of the species. The fallow deer network showed low uncertainty and bias even at lower sampling proportions, indicating the importance of a thoroughly sampled population while demonstrating the accuracy of our evaluation methods for smaller samples. Conclusions The protocol allows researchers to analyse GPS-based radio telemetry or other data to determine the reliability of social network metrics. The estimates enable the statistical comparison of networks under di↵erent conditions, such as analysing daily and seasonal changes in the density of a network. The methods can also guide methodological decisions in animal social network research, such as sampling design and allow more accurate ecological inferences from the available data. The R package aniSNA enables researchers to implement this workflow on their dataset, generating reliable inferences and guiding methodological decisions

Article↗
USGS2010

Methods for estimating flow-duration and annual mean-flow statistics for ungaged streams in Oklahoma

Flow statistics can be used to provide decision makers with surface-water information needed for activities such as water-supply permitting, flow regulation, and other water rights issues. Flow statistics could be needed at any location along a stream. Most often, streamflow statistics are needed at ungaged sites, where no flow data are available to compute the statistics. Methods are presented in this report for estimating flow-duration and annual mean-flow statistics for ungaged streams in Oklahoma. Flow statistics included the (1) annual (period of record), (2) seasonal (summer-autumn and winter-spring), and (3) 12 monthly duration statistics, including the 20th, 50th, 80th, 90th, and 95th percentile flow exceedances, and the annual mean-flow (mean of daily flows for the period of record). Flow statistics were calculated from daily streamflow information collected from 235 streamflow-gaging stations throughout Oklahoma and areas in adjacent states. A drainage-area ratio method is the preferred method for estimating flow statistics at an ungaged location that is on a stream near a gage. The method generally is reliable only if the drainage-area ratio of the two sites is between 0.5 and 1.5. Regression equations that relate flow statistics to drainage-basin characteristics were developed for the purpose of estimating selected flow-duration and annual mean-flow statistics for ungaged streams that are not near gaging stations on the same stream. Regression equations were developed from flow statistics and drainage-basin characteristics for 113 unregulated gaging stations. Separate regression equations were developed by using U.S. Geological Survey streamflow-gaging stations in regions with similar drainage-basin characteristics. These equations can increase the accuracy of regression equations used for estimating flow-duration and annual mean-flow statistics at ungaged stream locations in Oklahoma. Streamflow-gaging stations were grouped by selected drainage-basin characteristics by using a k-means cluster analysis. Three regions were identified for Oklahoma on the basis of the clustering of gaging stations and a manual delineation of distinguishable hydrologic and geologic boundaries: Region 1 (western Oklahoma excluding the Oklahoma and Texas Panhandles), Region 2 (north- and south-central Oklahoma), and Region 3 (eastern and central Oklahoma). A total of 228 regression equations (225 flow-duration regressions and three annual mean-flow regressions) were developed using ordinary least-squares and left-censored (Tobit) multiple-regression techniques. These equations can be used to estimate 75 flow-duration statistics and annual mean-flow for ungaged streams in the three regions. Drainage-basin characteristics that were statistically significant independent variables in the regression analyses were (1) contributing drainage area; (2) station elevation; (3) mean drainage-basin elevation; (4) channel slope; (5) percentage of forested canopy; (6) mean drainage-basin hillslope; (7) soil permeability; and (8) mean annual, seasonal, and monthly precipitation. The accuracy of flow-duration regression equations generally decreased from high-flow exceedance (low-exceedance probability) to low-flow exceedance (high-exceedance probability) . This decrease may have happened because a greater uncertainty exists for low-flow estimates and low-flow is largely affected by localized geology that was not quantified by the drainage-basin characteristics selected. The standard errors of estimate of regression equations for Region 1 (western Oklahoma) were substantially larger than those standard errors for other regions, especially for low-flow exceedances. These errors may be a result of greater variability in low flow because of increased irrigation activities in this region. Regression equations may not be reliable for sites where the drainage-basin characteristics are outside the range of values of independent vari

Scientific Investigations Report↗
USGS1974

Environmental impact analysis; the example of the proposed Trans-Alaska Pipeline

The environmental impact analysis made as required by the National Environmental Policy Act of 1969 for the proposed trans-Alaska pipeline included consideration of the (1) technologically complex and geographically extensive proposed project, (2) extremely different physical environments across Alaska along the proposed route and elsewhere in Alaska and in Canada along alternative routes, (3) socioeconomic environment of the State of Alaska, and (4) a wide variety of alternatives. The analysis was designed specifically to fit the project and environment that would be affected. The environment was divided into two general parts--natural physical systems and superposed socioeconomic systems--and those parts were further divided into discipline-oriented systems or components that were studied and analyzed by scientists of the appropriate discipline. Particular attention was given to potential feedback loops in the impact network and to linkages between the project's impacting effects and the environment. The results of the analysis as reported in the final environmental impact statement were that both unavoidable and threatened environmental impacts would result from construction, operation, and maintenance of the proposed pipeline system and the developments related to it. The principal unavoidable effects would be (1) disturbances of terrain, fish and wildlife habitat, and human environs, (2) the results of the discharge of effluent from the tanker-ballast-treatment facility into Port Valdez and of some indeterminate amount of oil released into the ocean from tank-cleaning operations at sea, and (3) the results associated with increased human pressures of all kinds on the environment. Other unavoidable effects would be those related to increase of State and Native Corporation revenues, accelerated cultural change of the Native population, and extraction of the oil and gas resource. The main threatened environmental effects would all be related to unintentional oil loss from the pipeline, from tankers, or in the oil field. Oil losses from the pipeline could be caused by direct or indirect effects of earthquakes, destructive sea waves, slope failure caused by natural or artificial processes, thaw-plug instability (in permafrost), differential settlement of permafrost terrain, and bed scour and bank erosion at stream crossings. Oil loss from tankers could be caused by accidents during transfer operations at Valdez and at destination ports and by casualties involving tankers and other ships. Comparison of alternative routes and transportation systems and of their environmental impacts provided information which indicates to the author that one corridor containing both oil and gas pipelines would have less environmental impact than would separate corridors. Considering also the threat to the marine environment that any tanker system would impose and the threat that zones of high earthquake frequency and magnitude would impose on pipelines, it is apparent to the author that environmental impact and cost would be least for a single-corridor on-land route that avoided earthquake zones. The alternative trans-Alaska-Canada routes would meet these criteria. The decisions of the U.S. Department of the Interior, the U.S. Congress, and the President of the United States in favor of the proposed trans-Alaska pipeline system indicate the relative weight given by the decision makers in balancing the importance of potential environmental consequences against the advantages to be derived from rapid resource development.

Circular↗
USGS2016

Sediment load and distribution in the lower Skagit River, Skagit County, Washington

The Skagit River delivers about 40 percent of all fluvial sediment that enters Puget Sound, influencing flood hazards in the Skagit lowlands, critically important estuarine habitat in the delta, and some of the most diverse and productive agriculture in western Washington. A total of 175 measurements of suspended-sediment load, made routinely from 1974 to 1993, and sporadically from 2006 to 2009, were used to develop and evaluate regression models of sediment transport (also known as “sediment-rating curves”) for estimating suspended-sediment load as a function of river discharge. Using a flow-range model and 75 years of daily discharge record (acquired from 1941 to 2015), the mean annual suspended-sediment load for the Skagit River near Mount Vernon, Washington, was estimated to be 2.5 teragrams (Tg, where 1 Tg = 1 million metric tons). The seasonal model indicates that 74 percent of the total annual suspended‑sediment load is delivered to Puget Sound during the winter storm season (from October through March), but also indicates that discharge is a poor surrogate for suspended‑sediment concentration (SSC) during the summer low-flow season. Sediment-rating curves developed for different time periods revealed that the regression model slope of the SSC-discharge relation increased 66 percent between the periods of 1974–76 and 2006–09 when suspended-sediment samples were collected, implying that changes in sediment supply, channel hydraulics, and (or) basin hydrology occurred between the two time intervals. In the relatively wet water year 2007 (October 1, 2006, through September 30, 2007), an automated sampler was used to collect daily samples of suspended sediment from which an annual load of 4.5 Tg was calculated, dominated by a single large flood event that contributed 1.8 Tg, or 40 percent of the total. In comparison, the annual load calculated for water year 2007 using the preferred flow-range model was 4.8 Tg (+6.7 percent), in close agreement with the measured value. Particle size affects sediment transport, fate and distribution across watersheds, and therefore is important for predicting how coastal environments, particularly deltas and beaches, will respond to changes in climate and sea-level. Particle-size analysis of winter storm samples indicated that about one-half of the suspended-sediment load consisted of fines (that is, silt- and clay-sized particles smaller than 0.0625 mm in diameter), and the remainder consisted of mostly fine- to medium-sized sand (0.0625–0.5 mm), whereas bedload during winter storm flows (about 1–3 percent of total sediment load) was predominantly composed of medium to coarse sand (0.25–1 mm). A continuous turbidity record from the Anacortes Water Treatment Plant (water years 1999–2013), used as a surrogate for the concentration of fines ( R 2 = 0.93, p = 4.2E-10, n = 17), confirms that about one-half of the mean annual suspended-sediment load is composed of fines. The distribution of flow through the delta distributaries (that is, the channels into which the main stem splits as it approaches the delta) is dynamic, with twice as much flow through the North Fork of the Skagit River relative to the South Fork during low-flow conditions, and close to equal flows in the two channels during high-flow conditions. Turbidity, monitored at several locations in the lower river in spring 2009, was essentially uniform among sites, indicating that fines are well mixed in the lower Skagit River system (defined as the Skagit River and all its distributaries downstream of the Mount Vernon streamgage). A strong relation ( R 2 = 0.95, p = 3.2E-14, n = 21; linear regression) between the concentration of fines and turbidity measured at various locations in summer 2009 indicates that turbidity is an effective surrogate for the concentration of fines, independent of location in the river, under naturally well-mixed fluvial conditions. This relation is especially useful for monitoring suspended sediment in western Washington rivers that are seasonally dominated by glacier meltwater because glacial melting typically produces suspended-sediment concentrations that are not well correlated with discharge. These results provide a comprehensive set of tools to estimate sediment delivery and delta responses of interest to scientists and resource managers including decision-makers examining options for flood hazard mitigation, estuary restoration, and climate change adaptation.

Washington↗
USGS2013

Vascular plant and vertebrate species richness in national parks of the eastern United States

Given the estimates that species diversity is diminishing at 50-100 times the normal rate, it is critical that we be able to evaluate changes in species richness in order to make informed decisions for conserving species diversity. In this study, we examined the potential of vascular plant species richness to be used as a surrogate for vertebrate species richness in the classes of amphibians, reptiles, birds, and mammals. Vascular plants, as primary producers, represent the biotic starting point for ecological community structure and are the logical place to start for understanding vertebrate species associations. We used data collected by the United States (US) National Park Service (NPS) on species presence within parks in the eastern US to estimate simple linear regressions between plant species richness and vertebrate richness. Because environmental factors may also influence species diversity, we performed simple linear regressions of species richness versus natural logarithm of park area, park latitude, mean annual precipitation, mean annual temperature, and human population density surrounding the parks. We then combined plant species richness and environmental variables in multiple regressions to determine the variables that remained as significant predictors of vertebrate species richness. As expected, we detected significant relationships between plant species richness and amphibian, bird, and mammal species richness. In some cases, plant species richness was predicted by park area alone. Species richness of mammals was only related to plant species richness. Reptile species richness, on the other hand, was related to plant species richness, park latitude and annual precipitation, while amphibian species richness was related to park latitude, park area, and plant species richness. Thus, plant species richness predicted species richness of different vertebrate groups to varying degrees and should not be used exclusively as a surrogate for vertebrate species richness. Plant species richness should be included with other variables such as area and climate when considering strategies to manage and conserve species in US National Parks. It is not always appropriate to draw conclusions about analyses of taxonomic surrogates from one area to another. Two patterns evident from the linear regressions were the increase in species richness with the increase of park area and with increase of vascular plant species richness. To test whether there were differences in these patterns among networks, we used analysis of covariance (ANCOVA). Differences among networks were detected only in bird species richness versus plant species richness and for all taxa except mammals for vertebrate species richness versus park area. Some of these results may be due to small sample size among networks, and therefore, low statistical power. Other factors that could have contributed to these results were differences in average park area and habitat heterogeneity among networks, latitudinal gradients, low variation in mean annual precipitation, and different use of vegetation by migratory species. Based on these results we recommend that management of biodiversity be approached from local and site specific criteria rather than applying management directives derived from other regions of the US. It is also recommended that analyses similar to those presented here be conducted for all national parks, once data become available for all networks in the US, to gain a better understanding of how vascular plant species richness, area, and vertebrate species richness are related in the US.

Natural Resource Technical Report NPS/NCR/NCRO/NRT↗
USGS2014

Simulation of the effects of rainfall and groundwater use on historical lake water levels, groundwater levels, and spring flows in central Florida

The urbanization of central Florida has progressed substantially in recent decades, and the total population in Lake, Orange, Osceola, Polk, and Seminole Counties more than quadrupled from 1960 to 2010. The Floridan aquifer system is the primary source of water for potable, industrial, and agricultural purposes in central Florida. Despite increases in groundwater withdrawals to meet the demand of population growth, recharge derived by infiltration of rainfall in the well-drained karst terrain of central Florida is the largest component of the long-term water balance of the Floridan aquifer system. To complement existing physics-based groundwater flow models, artificial neural networks and other data-mining techniques were used to simulate historical lake water level, groundwater level, and spring flow at sites throughout the area. Historical data were examined using descriptive statistics, cluster analysis, and other exploratory analysis techniques to assess their suitability for more intensive data-mining analysis. Linear trend analyses of meteorological data collected by the National Oceanic and Atmospheric Administration at 21 sites indicate 67 percent of sites exhibited upward trends in air temperature over at least a 45-year period of record, whereas 76 percent exhibited downward trends in rainfall over at least a 95-year period of record. Likewise, linear trend analyses of hydrologic response data, which have varied periods of record ranging in length from 10 to 79 years, indicate that water levels in lakes (307 sites) were about evenly split between upward and downward trends, whereas water levels in 69 percent of wells (out of 455 sites) and flows in 68 percent of springs (out of 19 sites) exhibited downward trends. Total groundwater use in the study area increased from about 250 million gallons per day (Mgal/d) in 1958 to about 590 Mgal/d in 1980 and remained relatively stable from 1981 to 2008, with a minimum of 559 Mgal/d in 1994 and a maximum of 773 Mgal/d in 2000. The change in groundwater-use trend in the early 1980s and the following period of relatively slight trend is attributable to the concomitant effects of increasing public-supply withdrawals and decreasing use of water by the phosphate industry and agriculture. On the basis of available historical data and exploratory analyses, empirical lake water-level, groundwater-level, and spring-flow models were developed for 22 lakes, 23 wells, and 6 springs. Input time series consisting of various frequencies and frequency-band components of daily rainfall (1942 to 2008) and monthly total groundwater use (1957 to 2008) resulted in hybrid signal-decomposition artificial neural network models. The final models explained much of the variability in observed hydrologic data, with 43 of the 51 sites having coefficients of determination exceeding 0.6, and the models matched the magnitude of the observed data reasonably well, such that models for 32 of the 51 sites had root-mean-square errors less than 10 percent of the measured range of the data. The Central Florida Artificial Neural Network Decision Support System was developed to integrate historical databases and the 102 site-specific artificial neural network models, model controls, and model output into a spreadsheet application with a graphical user interface that allows the user to simulate scenarios of interest. Overall, the data-mining analyses indicate that the Floridan aquifer system in central Florida is a highly conductive, dynamic, open system that is strongly influenced by external forcing. The most important external forcing appears to be rainfall, which explains much of the multiyear cyclic variability and long-term downward trends observed in lake water levels, groundwater levels, and spring flows. For most sites, groundwater use explains less of the observed variability in water levels and flows than rainfall. Relative groundwater-use impacts are greater during droughts, however, and long-term trends in water levels and flows were identified that are consistent with historical groundwater-use patterns. The sensitivity of the hydrologic system to rainfall is expected, owing to the well-drained karst terrain and relatively thin confinement of the Floridan aquifer system in much of central Florida. These characteristics facilitate the relatively rapid transmission of infiltrating water from rainfall to the water table and contribute to downward leakage of water to the Floridan aquifer system. The areally distributed nature of rainfall, as opposed to the site-specific nature of groundwater use, and the generally high transmissivity and low storativity properties of the semiconfined Floridan aquifer system contribute to the prevalence of water-level and flow patterns that mimic rainfall patterns. In general, the data-mining analyses demonstrate that the hydrologic system in central Florida is affected by groundwater use differently during wet periods, when little or no system storage is available (high water levels), compared to dry periods, when there is excess system storage (low water levels). Thus, by driving the overall behavior of the system, rainfall indirectly influences the degree to which groundwater use will effect persistent trends in water levels and flows, with groundwater-use impacts more prevalent during periods of low water levels and spring flows caused by low rainfall and less prevalent during periods of high water levels and spring flows caused by high rainfall. Differences in the magnitudes of rainfall and groundwater use during wet and dry periods also are important determinants of hydrologic response. An important implication of the data-mining analyses is that rainfall variability at subannual to multidecadal timescales must be considered in combination with groundwater use to provide robust system-response predictions that enhance sustainable resource management in an open karst aquifer system. The data-driven approach was limited, however, by the confounding effects of correlation between rainfall and groundwater use, the quality and completeness of the historical databases, and the spatial variations in groundwater use. The data-mining analyses indicate that available historical data when used alone do not contain sufficient information to definitively quantify the related individual effects of rainfall and groundwater use on hydrologic response. The knowledge gained from data-driven modeling and the results from physics-based modeling, when compared and used in combination, can yield a more comprehensive assessment and a more robust understanding of the hydrologic system than either of the approaches used separately.

Florida↗
USGS2001

Characterization of surface-water quality based on real-time monitoring and regression analysis, Quivira National Wildlife Refuge, south-central Kansas, December 1998 through June 2001

Because of the considerable wildlife benefits offered by the Quivira National Wildlife Refuge in south-central Kansas, there is a desire to ensure suitable water quality. To assess the quality of water flowing from Rattlesnake Creek into the refuge, the U.S. Geological Survey collected periodic water samples from December 1998 through June 2001 and analyzed the samples for physical properties, dissolved solids, total suspended solids, suspended sediment, major ions, nutrients, metals, pesticides, and indicator bacteria. Concentrations of 10 of the 125 chemicals analyzed did not meet water-quality criteria to protect aquatic life and drinking water in a least one sample. These were pH, turbidity, dissolved oxygen, dissolved solids, sodium, chloride, phosphorus, total coliform bacteria, E. coli bacteria, and fecal coliform bacteria. No metal or pesticide concentrations exceeded water-quality criteria. Twenty-two of the 43 metals analyzed were not detected, and 36 of the 46 pesticides analyzed were not detected. Because dissolved solids, sodium, chloride, fecal coliform bacteria, and other chemicals that are a concern for the health and habitat of fish and wildlife at the refuge cannot be measured continuously, regression equations were developed from a comparison of the analytical results of periodic samples and in-stream monitor measurements of specific conductance, pH, water temperature, turbidity, and dissolved oxygen. A continuous record of estimated chemical concentrations was developed from continuously recorded in-stream measurements. Annual variation in water quality was evaluated by comparing 1999 and 2000 sample data- the 2 years for which complete data sets were available. Median concentrations of alkalinity, fluoride, nitrate, and fecal coliform bacteria were smaller or did not change from 1999 to 2000. Dissolved solids, total suspended solids, sodium, chloride, sulfate, total organic nitrogen, and total phosphorus had increases in median concentrations from 1999 to 2000. Increases in the median concentrations of the major ions were expected due to decreased rainfall in 2000 and very low streamflow late in the year. Increases for solids and nutrients may have been due to the unusually high streamflow in the early spring of 2000. This was the time of year when fields were tilled, exposing solids and nutrients that were transported with runoff to Rattlesnake Creek. Load estimates indicate the chemical mass transported into the refuge and can be used in the development of total maximum daily loads (as specified by the U.S. Environmental Protection Agency) for water-quality contaminants in Rattlesnake Creek. Load estimates also were used to evaluate seasonal variation in water quality. Seasonal variation was most pronounced in the estimates of nutrient loads, and most of the nutrient load transported to the refuge occurred during just a few periods of surface runoff in the spring and summer. This information may be used by resource managers to determine when water-diversion strategies would be most beneficial. Load estimates also were used to calculate yields, which are useful for site comparisons. The continuous and real-time nature of the record of estimated concentrations, loads, and yields may be important for resource managers, recreationalists, or others for evaluating water-diversion strategies, making water-use decisions, or assessing the environmental effects of chemicals in time to prevent adverse effects on fish or other aquatic life at the refuge.

Kansas↗
USGS2020

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↗
USGS2018

Assessment of dissolved-selenium concentrations and loads in the lower Gunnison River Basin, Colorado, as part of the Selenium Management Program, from 2011 to 2016

The Gunnison Basin Selenium Management Program implemented a water-quality monitoring network in 2011 in the lower Gunnison River Basin in Colorado. Selenium is a trace element that bioaccumulates in aquatic food chains and can cause reproductive failure, deformities, and other harmful effects. This report presents the percentile values of selenium because regulatory agencies in Colorado make decisions based on the U.S. Environmental Protection Agency (EPA) Clean Water Act Section 303(d) that uses percentile values of concentration. Also presented are dissolved-selenium loads at 18 sites in the lower Gunnison River Basin for water years (WYs) 2011–2016 (October 1, 2010, through September 30, 2016). Annual dissolved-selenium loads were calculated for five sites with continuous U.S. Geological Survey (USGS) streamflow-gaging stations. Annual dissolved-selenium loads for WY 2011 through WY 2016 ranged from 179 and 391 pounds (lb) at Uncompahgre River at Colona to 11,100 and 17,300 lb at Gunnison River near Grand Junction (herein called Whitewater), respectively. Instantaneous loads were calculated for five sites with continuous U.S. Geological Survey (USGS) streamflow-gaging stations and 13 ancillary sites where discrete water-quality sampling also took place, using discrete water-quality samples and the associated discharge measurements collected during the period. Median instantaneous loads ranged from 0.01 pound per day (lb/d) at Smith Fork near Lazear to 33.0 lb/d at Whitewater. Mean instantaneous loads ranged from 0.06 lb/d at Smith Fork near Lazear to 36.2 lb/d at Whitewater. Most tributary sites in the basin had a median instantaneous dissolved-selenium load of less than 20.0 lb/day. In general, dissolved-selenium loads at Gunnison River main-stem sites showed an increase from upstream to downstream. The State of Colorado water-quality standard for dissolved selenium of 4.6 micrograms per liter (µg/L) was compared to the 85th percentiles for dissolved selenium at selected sites. Annual 85th percentiles for dissolved selenium were calculated for the five core sites having USGS streamflow-gaging stations using estimated dissolved-selenium concentrations from linear regression models. The 85th-percentile concentrations for WYs 2011–2016 based on this method ranged from 0.62 µg/L and 1.1µg/L at Uncompahgre River at Colona to 12.1 µg/L and 18.7 µg/L at Uncompahgre River at Delta. The 85th percentiles for dissolved selenium also were calculated for sites with sufficient data using water-quality samples collected during WYs 2011–2016. The annual 85th-percentile concentrations based on the discrete samples ranged from 0.16 µg/L and 0.17 µg/L at Gunnison River below Gunnison Tunnel to 62.2 µg/L and 170 µg/L at Loutzenhizer Arroyo at North River Road. A trend analysis was completed for Whitewater to determine if dissolved-selenium loads are increasing or decreasing. The trend analysis indicates a decrease of 9,100 lb from WY 1986 to WY 2016, a 40.8 percent reduction during the time period. The trend analysis for the annual dissolved-selenium load for WY 1994 to WY 2016 indicates a decrease of 6,300 lb per year, or 33.3 percent.

Colorado↗
USGS2026

Guidelines for producing integrated 210Pb and 14C age-models

Accurate reconstructions of past environmental changes are crucial in paleoecological research and require reliable chronologies of sedimentary archives. Establishing robust age-models and obtaining the most appropriate proxies for analysis is a complex scientific endeavor, requiring extensive resources and collaboration among specialists, including radiochronologists. Radiometric dating methods, such as 210 Pb and radiocarbon ( 14 C), are frequently employed to establish chronologies in aquatic sedimentary deposits and peat bogs. In this study, we review key aspects of sampling, analysis, and the principles underlying 210 Pb and 14 C age-models, focusing on methods for developing robust joint chronologies for paleoenvironmental research. Drawing largely from the authors' experiences and group discussions during and after a scientific workshop in 2022, we discuss important considerations for site selection, sampling strategies, and radiometric dating to construct integrated 210 Pb and 14 C age-models. Using expert consensus, this group – called Paleostats – aims to provide a set of best practices for other geochronologists with this methods paper. Among our conclusions, we emphasize the importance of accounting for site-specific factors such as prior information on sedimentation rates to establish appropriate sampling and analytical strategies. The use of appropriate coring devices can minimize disturbance to sediments and ensure the core surface remains intact and preserved until sectioning. Where excess 210 Pb is expected, sectioning at intervals of ≤1 cm provides an adequate sampling resolution for 210 Pb dating. Exceptions are possible, allowing for ~2–3 cm sections in areas with confirmed high sedimentation rates (e.g., > 1 cm yr −1 ). Recovering deeper core sections for 14 C dating with sufficient overlap allows for accounting errors in depth estimates made in the field. Special attention is advised during time intervals where validation proxies, such as the human-made radionuclides 137 Cs or post-bomb 14 C, are expected, and to determine the depth of secular equilibrium between 210 Pb and 226 Ra. Radiocarbon analyses are commonly performed by accelerator mass spectrometry, and age models are constructed mainly using Bayesian statistics with Markov Chain Monte Carlo techniques (e.g., Bacon ). A Bayesian approach ( Plum ) is now available for producing 210 Pb age-models, which infers the 210 Pb ex flux, eliminates the need for selecting an equilibrium depth, and allows dating cores with incomplete 210 Pb ex inventory. Plum offers improved chronologies by integrating raw 210 Pb and 14 C data, and these age-models can be enriched with other dating methodologies, such as identifying tephras and other well-recorded historical events. Harmonized reporting would contribute to making radiometric age-models reproducible, which would benefit from an international effort. Using 210 Pb and 14 C to produce integrated age-models may yield better insights into the interplay between natural and recent anthropogenic forcings on ecosystems. This can enhance our understanding of environmental processes and their impacts on climate change, ultimately supporting science-based assessments and decisions.

Earth-Science Reviews↗