SEARCH · Geology Reports

Results for “Environmental Data Science”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,045 records · Page 58Linked to original sources

USGS2021

Capacity assessment for Earth Monitoring, Analysis, and Prediction (EarthMAP) and future integrated monitoring and predictive science at the U.S. Geological Survey

Executive Summary Managers of our Nation’s resources face unprecedented challenges driven by the convergence of increasing, competing societal demands and a changing climate that affects the stability, vulnerability, and predictability of those resources. To help meet these challenges, the scientific community must take advantage of all available technologies, data, and integrative Earth systems modeling capacity to better inform resource and risk management decisions. This is the overarching goal of the U.S. Geological Survey (USGS) Earth Monitoring, Analysis, and Prediction (EarthMAP) vision: “By 2030, the USGS will deliver well integrated observations and predictions of the future state of natural systems—water, ecosystems, energy, minerals, hazards—at regional and national scales, working primarily with federal, state, and academic partners to develop and operate the capability” (U.S. Geological Survey, 2021). Providing more integrated Earth systems science and actionable information to decision makers, stakeholders, and the public requires a better understanding of the depth and distribution of existing capacity (capabilities, tools, and techniques) across the Bureau. Identifying existing capacity is also a critical first step toward gap analysis and targeted investments to increase capacity over time. The USGS formed a Capacity Assessment Team (CAT) and charged it with (1) conducting a Request for Information (RFI) to identify existing USGS expertise and activities supportive of integrated and predictive science to inform decision making, (2) developing a strategy and proof-of-concept for a continuously updated capacity assessment capability, and (3) identifying lessons learned to inform development of best practices for future capacity assessment efforts. The RFI took the form of a survey, with content guided by the science and technology needs identified in a USGS report titled “Grand Challenges for Integrated U.S. Geological Survey Science—A Workshop Report” (Jenni and others, 2017). The 44-question survey provided respondents the ability to rate their level of experience with a suite of priority disciplines, analysis and modeling approaches, technologies, and stakeholder engagement strategies and to enter optional narrative text for supporting context. An introductory portion focused on general science capacity assessment, followed by three sections targeting capabilities related to the foundational components of EarthMAP: (1) data and information integration, (2) integrated predictive science, and (3) actionable information. The survey results provided a high-level snapshot of USGS capacity in the targeted areas. Respondents (1,035 individuals) represented approximately 13 percent of the USGS across all mission areas and regions. Seventy-four percent of the respondents held a science-focused position title and the remainder had position titles in information technology, computer science, management, administrative, or other (contractors, volunteers, emeritus, and unknown). To provide greater insight into respondent capabilities and activities, information from the U.S. Department of the Interior and USGS enterprise information systems were used to further characterize topical expertise and organizational associations of survey respondents. To address the ongoing need to assess the Bureau’s capacity to address integrated predictive science priorities, the CAT developed a software-based proof-of-concept called the Integrated Science Assessment Information Database (iSAID) for assembling various information sources together toward making the full extent of USGS capabilities and scientific assets available for routine capacity assessment. This proof-of-concept is intended to serve as a catalyst for further development. The process of implementing the EarthMAP capacity assessment survey, analyzing survey responses, and developing the proof-of-concept resulted in lessons learned, findings, and recommendations. Example scenarios throughout the report demonstrate how capacity assessment data can inform science planning. Three overarching findings and recommendations are: (1) Finding: Capacity is limited in some critical disciplines, skills, and technology applications, but “sufficient” depends on the question and the need relative to availability at a given point in time. Recommendation: Develop an on-demand capacity assessment framework that enables rapid identification and evaluation of existing and available expertise to support decision needs as they arise. (2) Finding: Institutional barriers and lack of awareness constrain the ability of USGS staff to adopt new technologies, collaborate across administrative boundaries, and deliver actionable information to stakeholders in a timely manner. However, these barriers are not universally experienced. Recommendation: Pursue more targeted inquiries to clarify which institutional barriers are obstructing the adoption of new technologies and approaches or the sharing of expertise and equipment across organizational and regional boundaries. These inquiries should inform USGS leadership, mission areas, and regions whether policies can be revised or whether a lack of understanding is creating perceived obstacles. Highlight cases when staff have successfully adopted new technologies and approaches to advance EarthMAP priorities and provide actionable information in a timely manner to spread awareness of how perceived obstacles can be navigated and overcome when appropriate. (3) Finding: Examples of people and projects integrating across disciplines and scales and applying advanced approaches to meet complex stakeholder needs exist. Such examples provide transfer value across the spectrum from approach to decision making. Many projects, already underway, appear to meet elements of the EarthMAP vision, and the USGS has people who can provide leadership in multiple types of specific integrated science efforts. Recommendation: Use these findings as a starting point for near-term strategic planning for integrated science. Highlight, incentivize, and build on existing interdisciplinary predictive science and information delivery activities across the USGS to advance toward further realization of an EarthMAP capacity. The CAT efforts to develop and assess existing USGS capacity to advance the EarthMAP vision revealed a fundamental challenge for not only this effort but any effort to assess existing capacity: A considerable amount of thought, time, and effort is required to survey and assess capabilities and tools available to support a given need, yet best results are still likely to provide an incomplete assessment. To better meet the frequent need to assess capabilities, tools, products, and projects that address an expressed strategic priority, the CAT proposes the concept of an on-demand capacity assessment framework supported by a software package that dynamically pulls and integrates information from existing USGS information systems and public domain registries. Although existing USGS enterprise information systems currently lack the structure, cross-system consistency, interoperability, and stability to support a continuously updated capacity assessment capability, we identify reasonable near-term steps to improve the utility of information gathered on expertise and project capacity and to improve the consistency and completeness of information and the ability of USGS systems to share that information. The ability to search and characterize this information will make future assessments of capacity faster, more complete, more efficient, and more targeted. This approach would grow the Bureau’s capacity knowledge over time, iteratively improving the ability to access, leverage, and synthesize existing capabilities and assets as well as identify and fill critical gaps. The greatest promise for developing integrated science could lie in linking across existing projects and expertise to create a multi-project capacity for addressing large, complex environmental issues.

Open-File Report↗

USGS2025

When the wild things are: Defining mammalian diel activity and plasticity

Circadian rhythms are a mechanism by which species adapt to environmental variability and fundamental to understanding species behavior. However, we lack data and a standardized framework to accurately assess and compare temporal activity for species during rapid ecological change. Through a global network representing 38 countries, we leveraged 8.9 million mammalian observations to create a library of 14,587 standardized diel activity estimates for 445 species. We found that less than half the species’ estimates were in agreement with diel classifications from the reference literature and that species commonly used more than one diel classification. Species diel activity was highly plastic when exposed to anthropogenic change. Furthermore, body size and distributional extent were strongly associated with whether a species is diurnal or nocturnal. Our findings provide essential knowledge of species behavior in an era of rapid global change and suggest the need for a new, quantitative framework that defines diel activity logically and consistently while capturing species plasticity.

Science Advances↗

USGS2026

Small creatures, big challenge: A review of current knowledge and assessment of research needs to inform water withdrawal management for the Duck River's freshwater mussel fauna

The Duck River is a remarkable national resource, boasting outstanding biodiversity (Ahlstedt et al. 2017). It contains 151 fish species, over 60 mussel species, and 22 snail species, many of which are federally listed (Etnier and Starnes 1993; Johnson et al. 2013; Ahlstedt et al. 2017; Womble and Rosenberger 2025a). The Duck River watershed is at a pivotal time in its management and ecological history due to significant growth in human population and economic activity. The consequent increase in surface water demand presents challenges for sustainable water resource management. The Tennessee Department of Environment and Conservation (TDEC) has permitted or proposed to permit eight water utilities to withdraw over 77 million gallons per day (mgd) from the Duck River, representing > 40% increase from previously permitted amounts (USFWS 2025). On November 20, 2024, Tennessee Governor Bill Lee signed Executive Order No 108 to ensure long-term stability of the Duck River watershed through balancing economic growth, water resource management, and environmental conservation. Some primary components of this plan are the creation of a Comprehensive Regional Drought Management Plan, the development of a Habitat Conservation Plan, and the establishment of the Duck River Watershed Planning Partnership (DRWPP). Although the planning partnership is not focused on a sole aspect of the Duck River, but the watershed as a whole, the mussel assemblages of the Duck River have been a particular focus due to their importance in providing ecosystem services, sensitivity to water quality, conservation status, and sensitivity to periodic drought conditions in this system. With multiple state and federal agencies involved in managing the Duck River’s unique water resources, a focused review on how future water management could affect the freshwater mussel communities in the Duck River will contribute to this process. Our objectives of this review are to: 1. Review the effects of changing flow regimes on mussels, with emphasis on the effect of droughts/low flows, highlighting key processes/factors affecting mussels. 2. Highlight critical data gaps that could allow future, evidence-based water management decisions in this system. 3. Propose research topics to address critical data gaps and enable evidence-based water management decisions that consider the needs of freshwater mussels. Freshwater mussels are a critical component of riverine ecosystems, providing ecosystem services including water filtration, nutrient cycling, and habitat provision (Vaughn et al. 2004; Haag 2012) that benefit human communities. Environmental flows, also referred to as ecological or e-flows, refer to the quantity, timing, and quality of water flows necessary to maintain the ecological health of rivers, wetlands, and other aquatic ecosystems (Acreman and Dunbar 2004; Acreman 2016). For freshwater mussels, environmental flows include maintaining a flow regime and associated water quality regimes that are favorable for both juveniles and adults while also seasonally supporting host fish and the survival of settling glochidia (i.e., larval mussels) in suitable environments (Gates et al. 2015). Freshwater mussels are adapted to the natural flow regime of their native streams, including seasonal variability in flows. As sedentary filter feeders, mussels require adequate water flow to provide oxygen and food (Allen and Vaughn 2009); however, they also require periods of low flow to increase access to fish hosts and to ensure juveniles are not displaced (Gates et al 2015). Mussels are sensitive to water quality conditions, including pollutants, sedimentation, and artificial alterations of temperature and pH (Strayer and Malcom 2012). Mussel community characteristics, including presence, abundance, and species richness, generally increase with stream size (e.g., stream order); however, responses to changes in flow conditions are often species-specific, and the magnitude and timing of extreme flow conditions have more predictive power for mussel community response than general changes in flow (Lopez and Vaughn 2021). Complex hydraulic variables, such as relative shear stress and Reynolds number, can be strong predictors of mussel community responses because they combine simple hydraulic variables (e.g., velocity, depth) with stream channel characteristics, such as substrate composition and dynamics. Droughts directly and indirectly affect mussels by influencing their physiology, behavior, reproduction, and ecosystem roles. Mass mussel die-offs during droughts and low water conditions are widely documented (Vaughn et al. 2015; Sousa et al. 2018; Dubose et al. 2019). Droughts can also lead to a shift in mussel community composition, with drought-tolerant species becoming more dominant while more sensitive species experience declines (Galbraith et al, 2010; Lopez et al. 2022; Vaughn and Atkinson 2025). Droughts affect mussels through direct mortality from stranding as well as physiological stress from altered water temperature and quality, reduced habitat and food availability and quality, increased predation and competition, disruption of reproductive cycles, and increased susceptibility to disease and pathogens. Although the effects of droughts, water withdrawals, and low flow conditions on mussels are an increasing area of study, there is a general paucity of data-driven consensus regarding optimal mitigation approaches and strategies (Cushway et al. 2024). Long-term management decisions regarding water withdrawals, particularly during extended or multi-year drought conditions, may consider factors beyond direct mortality such as the sublethal effects of drought and low flow conditions on freshwater mussels for a comprehensive understanding of the effects of prolonged low flows on freshwater mussels and aquatic communities.

Tennessee↗

USGS2022

Pharmaceutical pollution of the world’s rivers

Environmental exposure to active pharmaceutical ingredients (APIs) can have negative effects on the health of ecosystems and humans. While numerous studies have monitored APIs in rivers, these employ different analytical methods, measure different APIs, and have ignored many of the countries of the world. This makes it difficult to quantify the scale of the problem from a global perspective. Furthermore, comparison of the existing data, generated for different studies/regions/continents, is challenging due to the vast differences between the analytical methodologies employed. Here, we present a global-scale study of API pollution in 258 of the world’s rivers, representing the environmental influence of 471.4 million people across 137 geographic regions. Samples were obtained from 1,052 locations in 104 countries (representing all continents and 36 countries not previously studied for API contamination) and analyzed for 61 APIs. Highest cumulative API concentrations were observed in sub-Saharan Africa, south Asia, and South America. The most contaminated sites were in low- to middle-income countries and were associated with areas with poor wastewater and waste management infrastructure and pharmaceutical manufacturing. The most frequently detected APIs were carbamazepine, metformin, and caffeine (a compound also arising from lifestyle use), which were detected at over half of the sites monitored. Concentrations of at least one API at 25.7% of the sampling sites were greater than concentrations considered safe for aquatic organisms, or which are of concern in terms of selection for antimicrobial resistance. Therefore, pharmaceutical pollution poses a global threat to environmental and human health, as well as to delivery of the United Nations Sustainable Development Goals.

Proceedings of the National Academy of Sciences↗

USGS2013

Borehole-explosion and air-gun data acquired in the 2011 Salton Seismic Imaging Project (SSIP), southern California: description of the survey

The Imperial and Coachella Valleys are being formed by active plate-tectonic processes. From the Imperial Valley southward into the Gulf of California, plate motions are rifting the continent apart. In the Coachella Valley, the plates are sliding past one another along the San Andreas and related faults (fig. 1). These processes build the stunning landscapes of the region, but also produce damaging earthquakes. Rupture of the southern section of the San Andreas Fault (SAF), from the Coachella Valley to the Mojave Desert, is believed to be the greatest natural hazard that California will face in the near future. With an estimated magnitude between 7.2 and 8.1, such an event would result in violent shaking, loss of life, and disruption of infrastructure (freeways, aqueducts, power, petroleum, and communication lines) that might bring much of southern California to a standstill. As part of the nation’s efforts to avert a catastrophe of this magnitude, a number of projects have been undertaken to more fully understand and mitigate the effects of such an event. The Salton Seismic Imaging Project (SSIP), funded jointly by the National Science Foundation (NSF) and the U.S. Geological Survey (USGS), seeks to understand, through seismic imaging, the structure of the Earth surrounding the SAF, including the sedimentary basins on which cities are built. The principal investigators (PIs) of this collaborative project represent the USGS, Virginia Polytechnic Institute and State University (Virginia Tech), California Institute of Technology (Caltech), Scripps Institution of Oceanography (Scripps), University of Nevada, Reno (UNR), and Stanford University. SSIP will create images of underground structure and sediments in the Imperial and Coachella Valleys and adjacent mountain ranges to investigate the earthquake hazards posed to cities in this area. Importantly, the images will help determine the underground geometry of the SAF, how deep the sediments are, and how fast earthquake energy can travel through the sediments. All of these factors determine how hard the earth will shake during a major earthquake. If we can improve on our understanding of how and where earthquakes will occur, and how strong their resultant shaking will be, then buildings can be designed or retrofitted accordingly in order to resist damage and collapse, and emergency plans can be adequately prepared. In addition, SSIP will investigate the processes of rifting and magmatism in the Salton Trough in order to better understand this important plate-boundary region. The Salton Trough is a unique rift in that subsidence is accompanied by huge influxes of infilling sediment from the Colorado River. Volcanism that accompanies the subsidence here is muted by these influxes of sediment. The Salton Trough, in the central part of the Imperial Valley, is apparently made up of entirely new crust: young sediment in the upper crust and basaltic intrusive rocks in the mid-to-lower crust (Fuis and others, 1984). Similar to the ultrasound and computed tomography (CT) scans performed by the medical industry, seismic imaging is a collection of techniques that enable scientists to obtain a picture of what is underground. The petroleum industry routinely uses these techniques to search for oil and gas at relatively shallow depths; however, the scope of this project demanded that we image as much as 30 km into the Earth’s crust. This project generated and recorded seismic waves, similar to sound waves, which move downward into the Earth and are bent (refracted) or echoed (reflected) back to the surface. SSIP acquired data in a series of intersecting lines that cover key areas of the Salton Trough. The sources of sound waves were detonations (shots) in deep boreholes, designed to create energy equivalent to magnitude 1–2 earthquakes. The study region routinely experiences earthquakes of these magnitudes, but earthquakes are not located in such a way as to permit us to create the detailed images we need for earthquake hazard assessment. Air gun bursts, generated in the Salton Sea along extensions of our onshore seismic lines, also were utilized as sound-wave sources. Temporary deployments of portable land seismometers, as well as ocean-bottom seismometers (OBSs) on the floor of the Salton Sea, recorded the energy from the land shots and air gun bursts. SSIP is similar to the Los Angeles Regional Seismic Experiments of 1994 and 1999 (LARSE I and II, respectively; Murphy and others, 1996; Fuis and others, 2001). The LARSE surveys demonstrated that the USGS and collaborators can safely and effectively conduct seismic imaging surveys in urban and nonurban areas, on lands owned and/or managed by many different types of agencies and entities. Information was produced that could not have been obtained any other way, and this information was key to changing the leading ideas about earthquake hazards at that time in the Los Angeles region. These surveys produced no significant environmental impact or damage to structures, and they did not trigger earthquakes.

California↗

USGS2025

Late Quaternary environmental change in eastern Beringia

Eastern Beringia (Alaska and western Yukon) is an extensive, high-latitude region of North America that remained largely unglaciated throughout the Quaternary. Consequently, its sedimentary deposits preserve long-term environmental records that have intrigued scientists for nearly a century. Recent advances in palaeoecological proxies and dating methods have proved critical in addressing long-standing questions about regional late Quaternary environmental change. At the same time, they have led to new and sometimes controversial hypotheses. This review covers recent discoveries and unresolved questions focused on the period 57,000–10,000 calendar years before C.E. 1950 (cal yr BP). The middle Wisconsin interstadial (57,000–30,000 cal yr BP) was a period of relative warmth in eastern Beringia, compared with the late Wisconsin (30,000–14,000 cal yr BP). Early in the interstadial occasional Picea woodland was present amongst widespread shrub tundra. Palaeoecological, sedimentary and isotopic data indicate that climate was cooler and drier than the Holocene, with high rates of aeolian activity. Megafauna typically associated with the ‘mammoth steppe’ ecosystem (woolly mammoth [ Mammuthus primigenius ], horse [ Equus ] and steppe-bison [ Bison priscus ]) were present in some abundance. The transition towards late Wisconsin cold-stage conditions (35,000–30,000 cal yr BP) coincided with the establishment of the Bering Land Bridge and featured expansion of spatially varied, herbaceous vegetation, sometimes associated with deep active layers. Sedimentary DNA ( seda DNA) and macrofossil evidence show vegetation was not a prairie-like grassland, and the term “steppe-tundra” is a better descriptor. Permafrost pore-ice isotopic (δ 18 O) records suggest a step change in one or more climate drivers ca. 30,000 cal yr BP, by which time steppe-tundra was established across eastern Beringia. It remains uncertain whether Picea survived cold-stage conditions within isolated refugia, or whether it recolonized from south of the Laurentide-Cordilleran ice sheets. Genetic data suggest that Picea probably survived in situ ; however, there is no definitive fossil evidence to support this. The end-Pleistocene transition from steppe-tundra to shrub tundra began ca. 15,000 cal yr BP and took place within decades at local scales. The expansion of woody taxa coincided with rising sea levels, reduced sea-ice extent and an abrupt shift in atmospheric circulation that enhanced precipitation. During this time, Earth's orbital configuration caused high early-summer temperatures and strong seasonality, creating growing conditions very different from today. The vegetation consisted of Salix and Betula shrub tundra with open areas of herbs and graminoids. During the deglacial warming trend, the Younger Dryas oscillation (12,800–11,700 cal yr BP) was variably expressed. It is generally evident in records affected by adjacent oceans but can be absent at sites in continental areas. These past conditions and paleoenvironmental changes have implications for contemporary issues: hypotheses about Pleistocene mammalian extinction; sensitivity of eastern Beringia to major oceanic reorganizations and high-frequency climate variability; the nature of woody plant expansion with climate warming; grazing, hydroclimate and fire as controls over ecosystems; the efficacy of “Pleistocene rewilding” for carbon capture.

Beringia↗

USGS2021

Synoptic analysis of a decade of daily measurements of SO2 emission in the troposphere from volcanoes of the global ground-based Network for Observation of Volcanic and Atmospheric Change

Volcanic plumes are common and far-reaching manifestations of volcanic activity during and between eruptions. Observations of the rate of emission and composition of volcanic plumes are essential to recognize and, in some cases, predict the state of volcanic activity. Measurements of the size and location of the plumes are important to assess the impact of the emission from sporadic or localized events to persistent or widespread processes of climatic and environmental importance. These observations provide information on volatile budgets on Earth, chemical evolution of magmas, and atmospheric circulation and dynamics. Space-based observations during the last decades have given us a global view of Earth's volcanic emission, particularly of sulfur dioxide ( SO 2 ). Although none of the satellite missions were intended to be used for measurement of volcanic gas emission, specially adapted algorithms have produced time-averaged global emission budgets. These have confirmed that tropospheric plumes, produced from persistent degassing of weak sources, dominate the total emission of volcanic SO 2 . Although space-based observations have provided this global insight into some aspects of Earth's volcanism, it still has important limitations. The magnitude and short-term variability of lower-atmosphere emissions, historically less accessible from space, remain largely uncertain. Operational monitoring of volcanic plumes, at scales relevant for adequate surveillance, has been facilitated through the use of ground-based scanning differential optical absorption spectrometer (ScanDOAS) instruments since the beginning of this century, largely due to the coordinated effort of the Network for Observation of Volcanic and Atmospheric Change (NOVAC). In this study, we present a compilation of results of homogenized post-analysis of measurements of SO 2 flux and plume parameters obtained during the period March 2005 to January 2017 of 32 volcanoes in NOVAC. This inventory opens a window into the short-term emission patterns of a diverse set of volcanoes in terms of magma composition, geographical location, magnitude of emission, and style of eruptive activity. We find that passive volcanic degassing is by no means a stationary process in time and that large sub-daily variability is observed in the flux of volcanic gases, which has implications for emission budgets produced using short-term, sporadic observations. The use of a standard evaluation method allows for intercomparison between different volcanoes and between ground- and space-based measurements of the same volcanoes. The emission of several weakly degassing volcanoes, undetected by satellites, is presented for the first time. We also compare our results with those reported in the literature, providing ranges of variability in emission not accessible in the past. The open-access data repository introduced in this article will enable further exploitation of this unique dataset, with a focus on volcanological research, risk assessment, satellite-sensor validation, and improved quantification of the prevalent tropospheric component of global volcanic emission.

Earth System Science Data↗

USGS2019

Case studies in groundwater contaminant fate and transport

A case study of groundwater contamination is a detailed study of a single site contaminated with a chemical or mixture that is known to be a problem at many sites. The goal of case studies is to provide insights into the physical, chemical, and biological processes controlling migration, natural attenuation, or remediation of common groundwater contaminants. Ideally, processes occurring at a case study site are representative of other sites so that knowledge gained from these intensive studies can be applied at thousands of sites where fewer data are available. Several characteristics of case studies contribute to their value. First, they may have tens to hundreds of monitoring wells, compared to fewer than ten wells at some contaminated sites. Second, some case studies continue for many years or even decades, providing insights into temporal progression of slow processes. Third, analytical methods prohibitively expensive for routine use or under development may be tested at case study sites. Finally, the ongoing characterization typical of case study sites builds a foundation of knowledge that facilitates sophisticated experimental design and testing of new methods. This article is divided into sections based on the contaminant type because the chemical and biological processes required for remediation vary for each contaminant. Most importantly, some contaminants can be biodegraded whereas metals and radionuclides cannot be destroyed but can be immobilized or rendered less toxic. The emphasis is on case studies of natural processes that control the fate and transport of contaminants in groundwater rather than on active remediation methods. The principles learned from these studies may form the basis for design of remedial strategies. The organic contaminants are divided into: petroleum hydrocarbons, fuel oxygenates, coal tar and wastes from manufactured gas plants, and chlorinated solvents. The inorganic contaminants covered are metals and radionuclides, arsenic, and nitrate. Case studies of mixed waste plumes from landfills are also described. Experimental sites where contaminants have been introduced into an aquifer as an emplaced source or a controlled release may not meet the above definition of case studies, but some are included because the overall goal is to impart lessons learned from detailed field studies. It is impossible to cover all case studies in this short format. Conversely, focusing on one or two does not convey the breadth of research results in entire range of case studies. Instead, the strategy is to describe the evolution of knowledge for each contaminant class while providing citations of relevant case studies. Much of the progress in understanding of the fate of contaminants in groundwater is based on laboratory studies; thus whenever possible, papers that included both field and laboratory results have been included among the citations. Two topics of growing importance have not been covered. These are the fate of pharmaceuticals in groundwater and discharge of contaminant plumes to surface water. These topics merit coverage in the future as knowledge grows and case studies increase in number.

Environmental Science↗

USGS2020

Population and habitat analyses for greater sage-grouse (Centrocercus urophasianus) in the bi-state distinct population segment—2018 update

Executive Summary The Bi-State Distinct Population Segment (Bi-State DPS) of greater sage-grouse ( Centrocercus urophasianus , hereinafter “sage-grouse”) represents a genetically distinct and geographically isolated population that straddles the border between Nevada and California. The primary threat to these sage-grouse populations is the expansion of single-leaf pinyon ( Pinus monophylla ) and Utah juniper ( Juniperus osteosperma ) into sagebrush ecosystems, which fragments and reduces population connectivity and survival. Other important threats include low water availability during brood-rearing, particularly during drought, and increased predation by common ravens ( Corvus corax ), a generalist predator often associated with anthropogenic resource subsidies. Although the Bi-State DPS occurs at high elevations relative to sage-grouse range-wide, changes in historical wildfire cycles and the conversion of native shrubs to invasive annual grasslands still threaten these populations. The Bi-State DPS has undergone multiple federal status assessments and associated litigation. For example, in October of 2013, the Bi-State DPS was proposed for listing as threatened under the Endangered Species Act of 1973 by the U.S. Fish and Wildlife Service (USFWS), then withdrawn in April 2015. The withdrawal decision was challenged, and in May 2018, a Federal district court ordered the withdrawal decision to be vacated, and USFWS was required to re-open the October 2013 listing evaluation. In response, the U.S. Geological Survey (USGS), with State and Federal collaborators, embarked on a multipronged analysis to provide current and best available science regarding population status of sage-grouse within the Bi-State DPS. Using data from a long-term monitoring program, we carried out four analytical study objectives, and here, we provide preliminary results of these analyses. First, we used integrated population modeling (IPM) to predict annual population abundance and annual finite rate of population change for the Bi-State DPS, as a whole, and for each subpopulation between 1995 and 2018. Because sage-grouse exhibit population cycles (periodic increases and decreases in abundance across approximately 6- to 10-year wavelengths), we estimated trends across three nested temporal scales that represent one (11 years), two (18 years), and three (24 years) complete population cycles. These estimates of relatively long-term averaged population change account for temporal (that is, interannual) variation. Our model predicted population abundance for the Bi-State DPS during 2018 at 3,305 individuals (2,247–4,683), with the majority occupying Bodie Hills and Long Valley. The model also predicted cyclic dynamics in abundance through time with evidence of 24-year population growth and slight trends of decline over the past 18 years. Specifically, across the Bi-State DPS as a whole, we estimated annual average at 0.99, 0.99, and 1.02 over the one, two, and three population cycles, which equated to a 10.5 percent, 16.6 percent decrease, and 60.0 percent increase in abundance over the 11-, 18-, and 24-year cycles. Estimated abundance in 2018 had not reached numbers lower than those predicted during 1995. However, we found spatial variation in population trends across the three cycles. Bodie Hills subpopulation comprised the greatest (1,521) and exhibited average annual greater than 1.0 across all periods resulting in average annual increases of 7 percent. This relatively large subpopulation has grown 5 times larger than what was predicted in 1995 while experiencing cyclical dynamics within that period. Conversely, other smaller subpopulations within the Bi-State DPS exhibited average annual equal to or less than 1.0 resulting in estimated 10-year risks of extirpation ranging from 2.0 to 76.1 percent. In general, evidence of decline among smaller subpopulations was greatest for the most recent period (2008–18) compared to a period that encompassed three full population cycles (24-year). This difference coincides with an intense period of drought that began in 2012. For comparative purposes as part of this first objective, we conducted a similar analysis for populations of sage-grouse within Nevada and California but outside the Bi-State DPS. We developed a region-wide and distance-weighted IPM using lek count from Nevada Department of Wildlife (NDOW) and California Department of Fish and Wildlife (CDFW) databases and with telemetry data collected by USGS across 12 sage-grouse subpopulations. Our models predicted similar patterns in population cycling outside the Bi-State DPS but with much stronger evidence of long-term declines across 24 years. Specifically, median averaged across each year of the 11-, 18-, and 24-year periods resulted in average annual values of 0.94, 0.97, and 0.99, respectively. These values equate to 41.0 percent, 38.5 percent, and 21.3 percent declines over the corresponding periods. Second, we used lek count data in a state-space modeling framework to compare trends in population abundance across different spatial scales (that is, leks versus Bi-State DPS). This hierarchical framework allowed us to disentangle declines associated with climate conditions as opposed to other local level factors that might signal the need for management intervention. Specifically, we identified 7 leks that were both declining and recently decoupled from larger spatial scale trends, typically governed by climatic conditions (referred to as soft or hard signals). The goal of this analysis was to provide an early warning system that might have implications for conservation actions at local scales. Third, we developed phenological (spring, summer–fall, and winter) and reproductive life stage (nesting, early brood-rearing, and late-brood rearing) based resource selection functions using various environmental covariates. We report rankings of variable importance for each season and life stage, developed maps of habitat selection indices (HSI), binned categories representing low, moderate, and high classes of quality (where any category greater than or equal to low indicated selected habitat) for each phenological season and life stage, and produced composite maps by selected phenological and reproductive stage to estimate annual habitat. Fourth, we used for each lek within the Bi-State DPS to carry out a spatial analysis that quantified substantial changes in the distribution of occupied habitat across long- (24-year) and short- (11-year) term periods. Owing to differences among available datasets, the long-term analysis primarily reflected spatial shifts among subpopulations comprising the majority of the Bi-State DPS (that is, Bodie Hills and Long Valley) while the short-term analysis also quantified changes among subpopulations along the periphery. Over long and short-term periods, the overall distribution of occupied habitat (as measured by 99 percent utilization distributions intersecting any quantified habitat) was reduced by 20,573 ha and 55,492 ha, respectively. Occupied core areas (as measured by 50 percent utilization distributions intersecting any quantified habitat) over long-term periods were solely located in Bodie Hills and Long Valley. Although nearly all subpopulations experienced contractions in occupied overall and core distribution, Bodie Hills experienced spatial expansion that occurred with concomitant spatial contraction at Long Valley over both periods. Subpopulations at the northern (Pine Nuts), central (Sagehen) and southern (White-Mountains) extents of the Bi-State DPS also experienced spatial contraction over the short-term period. These findings, coupled with those of population trends, indicate long-term patterns in redistribution of sage-grouse from Long Valley and peripheral subpopulations to Bodie Hills. That is, sage-grouse subpopulations at the periphery are declining while the largest population at the core is increasing, which could have meaningful impacts on overall metapopulation persistence. We provide evidence for loss of occupied habitat (reduced distribution) given local extirpation of subpopulations. Fifth, we calculated percentages of selected phenological, life stage, and annual habitat that each subpopulation contributed to the Bi-State DPS. We then intersected these maps with a composite estimate of occupied habitat from the fourth objective and calculated percentages of selected habitat likely occupied by sage-grouse that each subpopulation contributed to the Bi-State DPS. These values provide evidence for loss of occupied habitat and subsequent reductions in spatial distribution given reductions in abundance and, in some cases, extirpation of leks within subpopulations. Lastly, we carried out an initial in-depth analysis of selection for irrigated pastures and wet meadows during the brood-rearing stage for the Long Valley subpopulation. We chose this subpopulation because it represents a population core, representing 26.5 percent of total sage-grouse within the Bi-State DPS, and has exhibited long-term declines in abundance and distribution. This subpopulation is highly sensitive to precipitation and other factors that influence water availability. Models predicted higher use of the interior portions of irrigated pastures and wet meadows during late brood-rearing period, which represented a potentially risky use of habitat that was exacerbated during periods of low moisture (for example, drought, reduced water delivery, or both). Sage-grouse typically used edges of riparian areas and pastures, largely because the interior of these mesic areas consisted of considerably less overhead concealment cover (for example, shrubs) that likely resulted in a higher risk of mortality. We found that a lack of water delivery to pastures in the form of overwinter precipitation or diversion ditches increased the movements of sage-grouse to the interior of pastures. Although further investigation of water delivery impacts on chick survival are warrented, our initial findings regarding resource selection may explain recent declines in population growth at Long Valley.

California, Nevada↗

USGS2008

Hydrogeology and water quality of the Leetown area, West Virginia

The U.S. Geological Survey’s Leetown Science Center and the co-located U.S. Department of Agriculture’s National Center for Cool and Cold Water Aquaculture both depend on large volumes of cold clean ground water to support research operations at their facilities. Currently, ground-water demands are provided by three springs and two standby production wells used to augment supplies during periods of low spring flow. Future expansion of research operations at the Leetown Science Center is dependent on assessing the availability and quality of water to the facilities and in locating prospective sites for additional wells to augment existing water supplies. The hydrogeology of the Leetown area, West Virginia, is a structurally complex karst aquifer. Although the aquifer is a karst system, it is not typical of most highly cavernous karst systems, but is dominated by broad areas of fractured rock drained by a relatively small number of solution conduits. Characterization of the aquifer by use of fluorometric tracer tests, a common approach in most karst terranes, therefore only partly defines the hydrogeologic setting of the area. In order to fully assess the hydrogeology and water quality in the vicinity of Leetown, a multi-disciplinary approach that included both fractured rock and karst research components was needed. The U.S. Geological Survey developed this multi-disciplinary research effort to include geologic, hydrologic, geophysical, geographic, water-quality, and microbiological investigations in order to fully characterize the hydrogeology and water quality of the Leetown area, West Virginia. Detailed geologic and karst mapping provided the framework on which hydrologic investigations were based. Fracture trace and lineament analysis helped locate potential water-bearing fractures and guided installation of monitoring wells. Monitoring wells were drilled for borehole geophysical surveys, water-quality sampling, water-level measurements, and aquifer tests to characterize the quality of water and the hydraulic properties of the aquifer. Surface geophysical surveys provided a 3-dimensional view of bedrock resistivity in order to assess geologic and lithologic controls on ground-water flow. Borehole geophysical surveys were conducted in monitoring wells to assess the storage and movement of water in subsurface fractures. Numerous single-well, multi-well, and straddle packer aquifer tests and step-drawdown tests were conducted to define the hydraulic properties of the aquifer and to assess the role of bedrock fractures and solution conduits in the flow of ground water. Water samples collected from wells and springs were analyzed to assess the current quality of ground water and provide a baseline for future assessment. Microbiological sampling of wells for indicator bacteria and human and animal DNA provided an analysis of agricultural and suburban development impacts on ground-water quality. Light detection and ranging (LiDAR) data were analyzed to develop digital elevation models (DEMs) for assessing sinkhole distribution, to provide elevation data for development of a ground-water flow model, and to assess the distribution of major fractures and faults in the Leetown area. The flow of ground water in the study area is controlled by lithology and geologic structure. Bedrock, especially low permeability units such as the shale Martinsburg Formation and the Conococheague Limestone, act as barriers to water flowing down gradient and across bedding. This retardation of cross-strike flow is especially pronounced in the Leetown area, where bedding typically dips at steep angles. Highly permeable fault and fracture zones that disrupt the rocks in cross-strike directions provide avenues through which ground water can flow laterally across or through strata of low primary permeability. Significant strike parallel thrust faults and cross-strike faults typically coincide with larger solution conduits and act as drains for the more pervasive network of interconnected diffuse fractures. Results of borehole geophysical surveys indicate that although numerous fractures may intersect a borehole, only one or two of the fractures typically transmit most of the water to a well. The diffuse-flow dominated network of fractures that provides the majority of storage occupies only a small proportion of the total aquifer volume but constitutes the majority of porosity within the aquifer. Solution conduits, while occupying a relatively small volume of the overall aquifer, are especially important because they serve as primary drains for the ground-water flow system. Surface resistivity maps and cross-sectionsshow anomalous areas of low resistivities coincident with the prevailing geologic strike at N. 20º E., with major cross-strike faults, and with major springs in the region. Transmissivity derived from straddle packer tests was highly variable, and ranged over three orders of magnitude (1.8 x 10 -6 to 5.9 x 10 -3 ft 2 /d) in diffuse-flow fractures. A similar large variability in transmissivity was documented by single- and multi-well aquifer tests conducted in conduit-flow dominated portions of the aquifer (2.0 x 10 3 to 1.4 x 10 4 ft 2 /d) in lowland areas immediately adjacent to the Leetown Science Center. A stream-gaging station installed on Hopewell Run near the point where the stream exits the Leetown watershed indicates average daily streamflow for the Hopewell Run of approximately 11.2 ft 3 /s, and ranged from a minimum of 1.80 ft 3 /s on September 28, 2005, to a maximum of 73.0 ft 3 /s on December 11, 2003. Base-flow (ground-water) discharge surveys identified numerous small seeps adjacent to streams in the area. Hydrographs of the stage of Balch Spring show rapid response to individual storms. Strong correlation of the flow of Hopewell Run and Balch Spring indicates the nearby losing stream reach is partly responsible for higher fluctuations in the stage of Balch Spring. A water budget for the study period (2003-2005), based on measured precipitation and hydrograph analyses, is expressed as Precipitation (38.60 in/yr) = Surface Runoff (1.36 in/yr) + Ground-Water Discharge (17.73 in/yr) + Evapotranspiration (24.23 in/yr) – Change in storage (4.72 in/yr). Flow of ground water through the epikarst, a shallow zone of intensely weathered rock and regolith, can be rapid (on the order of days or weeks) as flow is concentrated in solution conduits. Flow within the intermediate and deeper zones is typically much slower. Eight dye-tracer tests conducted in the Leetown area found ground-water flow patterns to be divergent, with velocities ranging from about 12.5 to 610 ft/day and a median velocity of 50 ft/day. Estimates of ground-water age in carbonate rocks in the region are on the order of 15 years in the shallower portions of the aquifer to 50 years or older for deeper portions of the aquifer. Shallow springs can have a significant component of fairly young water (< 5 years in age). Ground-water samples collected from 16 sites (12 wells and 4 springs) in the Leetown area were analyzed for more than 340 constituents. Only turbidity, indicator bacteria, and radon were typically present in concentrations exceeding U.S. Environmental Protection Agency (USEPA) drinking-water or aquatic life standards.

West Virginia↗

USGS2022

Fishway Entrance Palisade

This technical report summarizes the work that was conducted by the University of Massachusetts Amherst and the United States Geological Survey (USGS), along with other project partners, on the Fishway Entrance Palisade (EP), a projected funded through the Department of Energy’s (DOE) funding opportunity titled ‘Innovative Solutions for Fish Passage at Hydropower Dams’ (DE‐FOA‐0001662). The period of performance ranged from September 1, 2018 through September 30, 2021. The EP is a novel fish passage engineering technology designed to provide more favorable entry conditions for fish and to reduce costs relative to conventional fishway auxiliary water systems (AWS). The EP project has four primary components. First, the Northeast United States Auxiliary Water Systems Database was created (Northeast Fishway Auxiliary Water Systems Database Section). The database, developed with material provided by the U.S. Fish and Wildlife Service, contains information on fishway type (e.g., lift, Denil, pool and weir) and Auxiliary Water System (AWS) details (e.g., water conveyance method, diffuser type) for 60 hydroelectric sites in the region. Findings indicate that nearly 4 out of every 10 fishway in the region is a fish lift and approximately 1 out of every 4 is a Denil ladder. The remainder are a mix of vertical slot fishways, pool and weirs, and Ice Harbor fishways. Furthermore, over half of all AWS systems use floor diffusers to discharge the auxiliary (or attraction) water into the entrance of a fishway, whereas only 14% use wall diffusers. Second, limited experiments on a conventional AWS with live, actively migrating fish were conducted at the USGS Easter Ecological Science Center (EESC) S.O. Conte Research Laboratory (Conventional Auxiliary Water System Experiments Section). This study determined how water velocity through a wall diffuser, without turning vanes or timber baffles to distribute the flow, affects the behavior and passage of adult American shad, a conservative surrogate species for migratory fish on the East Coast. Two gross diffuser velocity treatments were examined, 0.5 ft/s and 1.0 ft/s. These wall diffuser velocities represented current (0.5 ft/s) and past (1.0 ft/s) design criteria guidelines set forth by the USFWS North Atlantic-Appalachian Region (Rojas 2020; USFWS 2019). Six trials with a total of 151 American Shad were conducted in June of 2019 for the two treatments. No differences in American shad passage efficiency were discovered between the two treatments, while approximately 3 in every 4 attempts were successful at passing the diffuser. While these results may appear to indicate that the generally accepted gross wall diffuser velocity criteria for American shad of 0.5 ft/s could be safely increased to 1.0 ft/s, further analysis is warranted. Furthermore, it is unknown how other migratory and resident fish species that traverse these structures would be impacted by such a change. Studying the wall diffuser hydraulics led to an important AWS observation. Without turning vanes or timber baffles in this study, doubling the diffuser area was insufficient at producing the type of flow field change one may expect by halving the gross diffuser velocity. Instead, the flow fields throughout each treatments study area were similar, which led to similar results in shad performance. This not only highlights the importance of installing flow guidance devices like turning vanes, but also to the importance of properly maintaining them, which can be costly. Third, more expansive experiments on the novel EP were conducted in the spring of 2019 and 2021 (Fishway Entrance Palisade Experiments). The goal of this study was to determine how adult American shad responded to a variety of conditions at a full-scale EP. A total of six treatments were examined by changing the average auxiliary channel velocity between 1.0 and 5.0 ft/s in intervals of 1.0 ft/s and by inserting/removing an entrance gate at the opening of the fishway. Thirty trials with a total of 1,273 shad were conducted over the two years. In all treatments, at least ~7 out of every 10 fish successfully passed the EP diffuser and swam into the entrance channel within the 3.5-hour long trial, highlighting the general effectiveness of the novel AWS technology. In both study years, lower velocities through the EP diffuser led to increased shad performance, though performance peaked for the 2 ft/s velocity treatment. This treatment condition represents an approximate six-fold increase in gross diffuser velocity relative to conventional auxiliary water systems, which in turn presents opportunities for cost savings (e.g., reduction in diffuser size). Shad performance, in general, was worse in 2019 than in 2021, potentially due to the different run timing when our trials were conducted (2019 trials occurred near the end of the migration season, unlike in 2021). Treatments in 2019 had approximately a 20% reduction in entrance efficiency by the trial end, including a 16.7% drop for the 3 ft/s velocity treatment in 2019 relative to 2021 (the only carryover treatment between years). Lastly, adding an entrance gate caused a significant delay to entry. The time to 25% entry raised ~20 minutes from the near instantaneous 25% entry that was reported for the other treatments conducted in the same year (2021). Though by the end of the 3.5-hour trial, the overall entrance efficiency nearly matched those of the other 2021 treatments. The fourth and final component of the EP project was an economic analysis that focused on the cost of attraction and environmental flows (Modeling Power Generation Losses Due to Environmental and Fish Passage Attraction Flows at a Run-Of-River Hydroelectric Operation in the Northeast). The study assessed the economic impact of meeting environmental flow requirements at a representative hydroelectric facility and fish lift in the Northeast. An initial finding of the study was that there is a paucity of published data on the costs of meeting attraction and environmental flows. This is due, in part, to the proprietary nature of this data. To explore the costs associated with these flows, three types of environmental flows were assessed: upstream fishway attraction flows, downstream fishway attraction flows, and habitat maintenance minimum flows. A physics-based model was developed and calibrated with three years of hourly generation and flow data as inputs. Gage flow inputs were adjusted and used to calculate power generated. To address hydrologic variability, the model was executed to simulate 30 years of historical flows. Results indicate that both interannual and seasonal climatic factors impact the costs of meeting environmental flow requirements. Generation potential is most strongly curtailed during dry years in terms of maximizing the capacity factor (the percent of time a plant generates at capacity). Dry years, and especially dry summers, have the most significant costs associated with mitigation flows. Of the three types of flows, habitat flows are most costly in terms of power production, followed by upstream attraction flows. Downstream attraction flows are least costly. This finding is the likely result of differences in both flow rates and duration of the seasonal requirement for each flow. Overall, environmental flows represented a 2-12% loss in annual generation, but losses during a dry summer can reach over 20%.

Final Technical Report↗

USGS2006

Environmental contaminants and biomarker responses in fish from the Columbia River and its tributaries: spatial and temporal trends

Fish were collected from 16 sites on rivers in the Columbia River Basin (CRB) from September 1997 to April 1998 to document temporal and spatial trends in the concentrations of accumulative contaminants and to assess contaminant effects on the fish. Sites were located on the mainstem of the Columbia River and on the Snake, Willamette, Yakima, Salmon, and Flathead Rivers. Common carp (Cyprinus carpio), black bass (Micropterus sp.), and largescale sucker (Catostomus macrocheilus) were the targeted species. Fish were field-examined for external and internal lesions, selected organs were weighed to compute somatic indices, and tissue and fluid samples were preserved for fish health and reproductive biomarker analyses. Composite samples of whole fish, grouped by species and gender, from each site were analyzed for organochlorine and elemental contaminants using instrumental methods and for 2,3,7,8-tetrachloro dibenzo-p-dioxin-like activity (TCDD-EQ) using the H4IIE rat hepatoma cell bioassay. Overall, pesticide concentrations were greatest in fish from lower CRB sites and elemental concentrations were greatest in fish from upper CRB sites. These patterns reflected land uses. Lead (Pb) concentrations in fish from the Columbia River at Northport and Grand Coulee, Washington (WA) exceeded fish and wildlife toxicity thresholds (> 0.4 ??g/g). Selenium (Se) concentrations in fish from the Salmon River at Riggins, Idaho (ID), the Columbia River at Vernita Bridge, WA, and the Yakima River at Granger, WA exceeded toxicity thresholds for piscivorous wildlife (> 0.6 ??g/g). Mercury (Hg) concentrations in fish were elevated throughout the basin but were greatest (> 0.4 ??g/g) in predatory fish from the Salmon River at Riggins, ID, the Yakima River at Granger, WA, and the Columbia River at Warrendale, Oregon (OR). Residues of p,p???-DDE were greatest (> 0.8 ??g/g) in fish from agricultural areas of the Snake, Yakima, and Columbia River basins but were not detected in upper CRB fish. Other organochlorine pesticides did not exceed toxicity thresholds in fish or were detected infrequently. Total polychlorinated biphenyls (PCBs; > 0.11 ??g/g) and TCDD-EQs (> 5 pg/g) exceeded wildlife guidelines in fish from the middle and lower CRB, and ethoxyresorufin O-deethylase (EROD) activity was also elevated at many of the same sites. Temporal trend analysis indicated decreasing or stable concentrations of Pb, Se, Hg, p,p???-DDE, and PCBs at most sites where historical data were available. Altered biomarkers were noted in fish throughout the CRB. Fish from some stations had responded to chronic contaminant exposure as indicated by fish health and reproductive biomarker results. Although most fish from some sites had grossly visible external or internal lesions, histopathological analysis determined these to be inflammatory responses associated with helminth or myxosporidian parasites. Many largescale sucker from the Columbia River at Northport and Grand Coulee, WA had external lesions and enlarged spleens, which were likely associated with infections. Intersex male smallmouth bass (Micropterus dolomieu) were found in the Snake River at Lewiston, ID and the Columbia River at Warrendale, OR. Male bass, carp, and largescale sucker containing low concentrations of vitellogenin were common in the CRB, and comparatively high concentrations (> 0.3 mg/mL) were measured in male fish from the Flathead River at Creston, Montana, the Snake River at Ice Harbor Dam, WA, and the Columbia River at Vernita Bridge, WA and Warrendale, OR. Results from our study and other investigations indicate that continued monitoring in the CRB is warranted to identify consistently degraded sites and those with emerging problems. ?? 2005 Elsevier B.V. All rights reserved.

Columbia River Basin↗

USGS2019

Using environmental DNA (eDNA) to assess the presence of cavefish and cave crayfish populations in caves of the Ozark Highlands

Many cavefishes and cave crayfishes are considered of conservation concern; however, sampling these species is inherently difficult given their occupied environments. The goal of our project was to verify the presence of select karst organisms while developing the foundation for sampling approaches that might be useful to conservation and management agencies. Our project objectives were to develop assays to amplify deoxyribonucleic acid (DNA) from several species of Ozark cavefishes and cave crayfishes and complete an initial surveillance of locations across the Ozark Highlands using environmental DNA (eDNA). Using DNA either provided by agency cooperators or that we extracted from tissue samples, we PCR amplified and then sequenced the Cytochrome Oxidase 1 (CO1) gene for cave crayfishes and the NADH Dehydrogenase Subunit 2 (ND2) gene for cavefishes. We developed species-specific primers and probes for five cave crayfishes and two cavefishes. From February 2017 to May 2017, we sampled 1–5 sampling units from 42 caves, wells, and springs (i.e., sites) using eDNA and traditional visual surveys. We measured physicochemical parameters at each sampling unit to estimate detection probability associated with both techniques. We also calculated two occupancy covariates for each site using geospatial data. We successfully amplified Troglichthys rosae DNA from the environment and detected DNA representing this species at 24 of 40 sites. At 16 of the sites where we detected T. rosae DNA, we did not visually observe the species. Although our assay for Typlichthys eigenmanni successfully amplified the target DNA from the environment, it also resulted in false absences where the species was visually confirmed. Using eDNA to detect cave crayfishes was much more difficult. The assay for Cambarus subterraneus did not work for eDNA samples and we were unable to pick up DNA from the environment, even at locations where it was visually confirmed. Alternatively, the eDNA surveys worked well for C. tartarus and we were able to amplify DNA at every site where it was visually observed. Our assay for C. aculabrum was based on a single sample obtained from GenBank, and did not amplify eDNA from field samples. Lastly, our eDNA results from samples in the known range of Orconectes stygocaneyi suggested the species may be found at an additional cave. Detection using eDNA based on our O. stygocaneyi assay was likely low because it was designed from a pseudogene; however, positive eDNA samples were sequenced to confirm species-specific DNA. Detection probability of both cavefishes and cave crayfishes varied by survey technique and was influenced by water volume, water clarity, water velocity, and substrate. Detection of cavefishes and cave crayfishes via visual surveys decreased when water volume increased, whereas detection using eDNA increased with greater water volume. Detection between taxa using either sample method was highest in habitats classified by fine substrates, except for eDNA detection of crayfishes which was greatest in coarse substrates. Detection of cavefishes increased with water clarity, but detection of cave crayfishes increased with turbidity. Detection probability of both cavefishes and crayfishes using eDNA increased slightly with water velocity, but decreased with visual surveys as water velocity increased. Occupancy by both taxa was positively related to particular geologic series. Crayfish occupancy was negatively related to fine-scale anthropogenic disturbance (i.e., 500-m buffer around the site), whereas crayfish showed no relationship with disturbance. Our results suggest possible range extensions, provide insights to factors driving detection using both sample techniques, and suggest areas where recharge zones may be shared among caves. Future efforts focused on a comprehensive evaluation of genetic diversity among cave crayfishes to improve assay design could improve detection and the applicability of eDNA as a supplemental and non-invasive sampling approach.

Cooperator Science Series↗

USGS2013

Geomorphic and vegetation processes of the Willamette River floodplain, Oregon: current understanding and unanswered science questions

This report summarizes the current understanding of floodplain processes and landforms for the Willamette River and its major tributaries. The area of focus encompasses the main stem Willamette River above Newberg and the portions of the Coast Fork Willamette, Middle Fork Willamette, McKenzie, and North, South and main stem Santiam Rivers downstream of U.S. Army Corps of Engineers dams. These reaches constitute a large portion of the alluvial, salmon-bearing rivers in the Willamette Basin. The geomorphic, or historical, floodplain of these rivers has two zones - the active channel where coarse sediment is mobilized and transported during annual flooding and overbank areas where fine sediment is deposited during higher magnitude floods. Historically, characteristics of the rivers and geomorphic floodplain (including longitudinal patterns in channel complexity and the abundance of side channels, islands and gravel bars) were controlled by the interactions between floods and the transport of coarse sediment and large wood. Local channel responses to these interactions were then shaped by geologic features like bedrock outcrops and variations in channel slope. Over the last 150 years, floods and the transport of coarse sediment and large wood have been substantially reduced in the basin. With dam regulation, nearly all peak flows are now confined to the main channels. Large floods (greater than 10-year recurrence interval prior to basinwide flow regulation) have been largely eliminated. Also, the magnitude and frequency of small floods (events that formerly recurred every 2–10 years) have decreased substantially. The large dams trap an estimated 50–60 percent of bed-material sediment—the building block of active channel habitats—that historically entered the Willamette River. They also trap more than 80 percent of the estimated bed material in the lower South Santiam River and Middle and Coast Forks of the Willamette River. Downstream, revetments further decrease bed-material supply by an unknown amount because they limit bank erosion and entrainment of stored sediment. The rivers, geomorphic floodplain, and vegetation within the study area have changed noticeably in response to the alterations in floods and coarse sediment and wood transport. Widespread decreases have occurred in the rates of meander migration and avulsions and the number and diversity of landforms such as gravel bars, islands, and side channels. Dynamic and, in some cases, multi-thread river segments have become stable, single-thread channels. Preliminary observations suggest that forest area has increased within the active channel, further reducing the area of unvegetated gravel bars. Alterations to floods and sediment transport and ongoing channel, floodplain, and vegetation responses result in a modern Willamette River Basin. Here, the floodplain influenced by the modern flow and sediment regimes, or the functional floodplain, is narrower and inset with the broader and older geomorphic floodplain. The functional floodplain is flanked by higher elevation relict floodplain features that are no longer inundated by modern floods. The corridor of present- day active channel surfaces is narrower, enabling riparian vegetation to establish on formerly active gravel bar surfaces. The modern Willamette River Basin with its fundamental changes in the flood, sediment transport, and large wood regimes has implications for future habitat conditions. System-wide future trends probably include narrower floodplains and a lower diversity of landforms and habitats along the Willamette River and its major tributaries compared to historical patterns and today. Furthermore, specific conditions and future trends will probably vary between geologically stable, anthropogenically stable, and dynamic reaches. The middle and lower segments of the Willamette River are geologically stable, whereas the South Santiam and Middle Fork Willamette Rivers were historically dynamic, but are now largely stable in response to flow regulation and revetment construction. The upper Willamette and North Santiam Rivers retain some dynamic characteristics, and provide the greatest diversity of aquatic and riparian habitats under the current flow and sediment regime. The McKenzie River has some areas that are more dynamic, whereas other sections are stable due to geology or revetments. Historical reductions in channel dynamism also have implications for ongoing and future recruitment and succession of floodplain forests. For instance, the succession of native plants like black cottonwood is currently limited by (1) fewer low-elevation gravel bars for stand initiation; (2) altered streamflow during seed release, germination, and stand initiation; (3) competition from introduced plant species; and (4) frequent erosion of young vegetation in some locations because scouring flows are concentrated within a narrow channel corridor. Despite past alterations, the Willamette River Basin has many of the physical and ecological building blocks necessary for highly functioning rivers. Management strategies, including environmental flow programs, river and floodplain restoration, revetment modifications, and reclamation of gravel mines, are underway to mitigate some historical changes. However, there are some substantial gaps in the scientific understanding of the modern Willamette basin that is needed to efficiently integrate these blocks and to establish realistic objectives for future conditions. Unanswered questions include: 1. What is the distribution and diversity of landforms and habitats along the Willamette River and its tributaries? 2. What is the extent of today’s functional floodplain—the part of the river corridor actively formed and modified by fluvial processes? 3. How are landforms and habitats in the Willamette River Basin created and sustained by present-day flow and sediment conditions? 4. How is the succession of native floodplain vegetation shaped by present-day flow and sediment conditions? Answering these questions will produce baseline data on the current distributions of landforms and habitats (question 1), the extent of the functional floodplain (question 2), and the effects of modern flow and sediment regimes on future floodplain landforms, habitats, and vegetation succession (questions 3 and 4). Addressing questions 1 and 2 is a logical next step because they underlie questions 3 and 4. Addressing these four questions would better characterize the modern Willamette Basin and help in implementing and setting realistic targets for ongoing management strategies, demonstrating their effectiveness at the site and basin scales, and anticipating future trends and conditions.

Oregon↗

USGS2005

Environmental contaminants and biomarker responses in fish from the Rio Grande and its U.S. tributaries: Spatial and temporal trends

We collected, examined, and analyzed 368 fish of seven species from 10 sites on rivers of the Rio Grande Basin (RGB) during late 1997 and early 1998 to document temporal and geographic trends in the concentrations of accumulative contaminants and to assess contaminant effects on the fish. Sites were located on the mainstem of the Rio Grande and on the Arroyo Colorado and Pecos River in Texas (TX), New Mexico (NM), and Colorado. Common carp (Cyprinus carpio) and largemouth bass (Micropterus salmoides) were the targeted species. Fish were examined in the field for internal and external visible gross lesions, selected organs were weighed to compute ponderal and organosomatic indices, and samples of tissues and fluids were obtained and preserved for analysis of fish health and reproductive biomarkers. Whole fish from each station were composited by species and gender and analyzed for organochlorine chemical residues and elemental contaminants using instrumental methods, and for 2,3,7,8-tetrachloro dibenzo-p-dioxin-like activity (TCDD-EQ) using the H4IIE rat hepatoma cell bioassay. Overall, fish from lower RGB stations contained greater concentrations of organochlorine pesticide residues and appeared to be less healthy than those from sites in the central and upper parts of the basin, as indicated by a general gradient of residue concentrations and biomarker responses. A minimal number of altered biomarkers and few or no elevated contaminant concentrations were noted in fish from the upper RGB. The exception was elevated concentrations [up to 0.46 ??g/g wet-weight (ww)] of total mercury (Hg) in predatory species from the Rio Grande at Elephant Butte Reservoir, NM, a condition documented in previous studies. Arsenic (As) and selenium (Se) concentrations were greatest in fish from sites in the central RGB; Se concentrations in fish from the Pecos River at Red Bluff Lake, TX and from the Rio Grande at Langtry, TX and Amistad International Reservoir, TX exceeded published fish and wildlife toxicity thresholds. In the lower RGB, residues of p,p???-DDT metabolites (???1.69 ??g/g ww), chlordane-related compounds (???0.21 ??g/g ww), dieldrin (???0.0.05 ??g/g ww), and toxaphene (???2.4 ??g/g ww) were detected in fish from most sites; maximum concentrations were in channel catfish (Ictalurus punctatus) from the Arroyo Colorado at Harlingen, TX. Concentrations of one or more residues exceeded toxicity thresholds for fish and wildlife in fish from this site and from the Rio Grande at Mission, TX and Brownsville, TX; however, concentrations were lower than those reported by previous studies. In addition, the proportional concentrations of p,p???-DDT at all sites were low, indicating weathered DDT rather than the influx of new material. Concentrations of total PCBs (<0.05 ??g/g ww) and TCDD-EQ (???6 pg/g ww) were comparatively low in all samples. Hepatic ethoxyresorufin O-deethylase (EROD) activity in some fish was elevated relative to reference rates at most sites, but was generally lower than previously reported activity in fish from heavily contaminated locations. The comparatively low PCB and TCDD-EQ concentrations together with elevated EROD activity may reflect exposure to polycyclic aromatic hydrocarbons. Reproductive biomarkers were consistent with chronic contaminant exposure at lower RGB sites; comparatively large percentages of intersex male largemouth bass, relatively low gonadosomatic indices, and elevated plasma vitellogenin concentrations in male fish were noted at three of the four stations. Large percentages of atretic eggs were also observed in the ovaries of female common carp from the Rio Grande at Brownsville, TX. Although many of the conditions noted may have other causes in addition to contaminant exposure, the biomarker results for the lower RGB sites are consistent with subtle responses of fish to contaminants, an interpretation supported by the chemical data of this and other investigations.

Science of the Total Environment↗

USGS1993

1993 Annual Report: San Francisco estuary regional monitoring program for trace substances

This first annual report of the San Francisco Estuary Regional Monitoring Program contains the results of monitoring measurements made in 1993. Measurements of conventional water quality parameters and trace contaminant concentrations were made at 16 stations throughout the Estuary three times during the year: the wet period (March), during declining Delta outflow (May), and during the dry period (September). Water toxicity tests were conducted at 8 of those stations. Measurements of sediment quality and contaminant concentrations were made at the same 16 stations during the wet and dry sampling periods. Sediment toxicity was measured at 8 of those stations. Transplanted, bagged bivalve bioaccumulation and condition was measured at 11 stations during the wet and dry sampling periods. Water Monitoring. Total or near-total (dissolved + particulate, see text) arsenic, cadmium, selenium, and dissolved (0.45 µm filtered) arsenic, cadmium, copper, nickel, silver, and zinc in water were highest in the South Bay. In general, dissolved metals in water were usually lowest in the Central Bay due to ocean influences. Near-total nickel and total mercury in water were highest in the northern estuary (San Pablo and Suisun Bays). Dissolved chromium and lead were highest at the Sacramento and San Joaquin River confluence stations. Six of the ten dissolved trace metals were highest in March during high runoff. Dissolved and total arsenic, selenium, and near-total cadmium were highest in September. Concentrations of trace organic contaminants are reported for the March sampling period. Total PAHs and PCBs were highest in the South Bay, but PCBs were also high in the Napa River. Dissolved PAHs were highest in the Central Bay, and dissolved PCBs were highest in the Napa River. Total and dissolved pesticides were highest in the Sacramento River and in the Extreme South Bay. Concentrations of trace elements in water (except selenium) were usually closely related with other environmental parameters. Total or near-total metals concentrations in water were most often associated with the amount of particulate material (TSS) in the water. Dissolved concentrations were usually associated with salinity or dissolved organic carbon (DOC) content. Dissolved PAHs were well correlated with TSS, but dissolved and total trace organic contaminants were poorly correlated with other water parameters. Based on deviations from conservative mixing of fresh and salt water, three different patterns of possible sources of metals were identified in 1993. For dissolved chromium and lead, rivers and local runoff appeared to be important sources. For dissolved arsenic, cadmium, copper, and nickel year-round inputs from the South Bay appeared to be important sources. Dissolved mercury, selenium, and zinc were associated with local runoff in the South Bay during the wet period. Dissolved silver did not fit any of these patterns. Although most contaminant concentrations were below water quality objectives, several trace contaminants were above the objectives at some stations. Comparisons to water quality objectives are used as a guide for evaluation of contaminant concentrations, but there are some differences in the way the RMP data are measured and that prescribed for regulatory purposes (see text). Concentrations of 5 metals in water were above EPA or Regional Basin Plan water quality objectives at six stations (see Table 30). Most of these elevated levels occurred at the northern estuary stations. Total PCB concentrations were above EPA human health objectives at all RMP stations. The pesticides chlordane, dieldrin, and DDTs were above the EPA objectives at several RMP stations, particularly at the northern-most, and river confluence stations. Although some of the contaminant concentrations were above water quality objectives, water toxicity tests (96 hour algal growth and 48 hour bivalve larval development tests) did not indicate toxicity (sometimes inconclusive) associated with the water samples collected at any of the RMP stations in 1993. Exposure to Bay San Francisco Estuary Regional Monitoring Program Regional Monitoring Program 1993 Report ii water actually enhanced algal growth at most stations. In addition to the Estuary-wide sampling, the Sacramento and San Joaquin Rivers were sampled upstream from their confluence. Stations in each river were sampled six times over a 6 week period of high flows. In the Sacramento River, seven of the ten dissolved metals measured had concentrations lower than those measured at the river confluence stations. Some metals concentrations in the San Joaquin River were higher, and some were lower than concentrations from the river confluence station. Metals concentrations in the Sacramento River were poorly related to river flow because the station at Rio Vista is under considerable tidal influence. In the San Joaquin River, flows were inversely related to 7 of 10 total metals concentrations. Sediment Monitoring. Concentrations of silver, mercury, and lead in sediment were highest in the South Bay. However, concentrations of most trace metals in sediments were highest in the northern estuary at stations with the finest (silt, clay) sediments. The northern estuary stations with the coarsest (sand, shell) sediments generally had the lowest metals concentrations. There were differences in concentrations of cadmium, lead, and selenium in sediments between the sampling periods, but no consistent trend as to which sampling period had higher values. In September, PAHs and PCBs in sediments were highest in the Central Bay, but pesticides in sediments were highest in the northern estuary and Extreme South Bay. NOAA’s Median Effects Ranges (ERM) for sediments were used as a guide for evaluation of sediment contaminant concentrations. Nickel was the only trace contaminant in sediment above the ERM guidelines, and it was high at all RMP stations. These high levels are probably due to natural, geologic sources. Although sediment contaminant concentrations were below ERMs, sediment toxicity tests (10 day amphipod mortality, and 48 hour bivalve larval development in elutriates) indicated toxicity at all stations tested. Sediment factors that could have caused the toxicity were not investigated. Bivalve Bioaccumulation. Mussels, oysters, and freshwater clams were transplanted to the RMP stations to evaluate bioaccumulation of trace substances. Trace metals were bioaccumulated at nearly all RMP stations. However, arsenic, lead, and mercury did not appear to bioaccumulate. There was generally more bioaccumulation during the dry season than during the wet season. In September, PAHs, PCBs, and pesticides accumulated in all samples. Bioaccumulation of PAHs and pesticides was generally highest at the river confluence stations, and the Napa River. PCBs accumulated most at Redwood Creek. There were substantial differences in the degree of bioaccumulation among the species. Oysters appeared to accumulate higher concentrations of trace metals than the other species, especially copper, which may be a natural phenomenon. There are no established tissue contaminant standards for trace metal and organic contaminants. Therefore, comparisons to Median International Standards (MIS) for human consumption, or U.S. Food and Drug Administration (USFDA) action levels for trace organics are used to evaluate the bioaccumulation results. Concentrations of selenium were higher than MIS guidelines at all stations during the wet season. Other trace metal concentrations were higher than MIS guidelines at various stations during one or the other sampling period. However, none of the bivalves contained concentrations above the USFDA or National Academy of Sciences (NAS) guidelines for trace organic contaminants. The transplanted bivalves survived well at all stations except in the Napa River where less than 35% survived during both sampling seasons. Measures of bivalve condition (dry weight, shell volume) indicated that bivalves deployed in the Central Bay grew significantly, but those at most other stations actually lost weight. Whether these differences were due to natural causes such as salinity or food supply, or to contamination, was not determined. Pilot Studies. Two pilot monitoring studies were conducted in 1993. A pilot study of Estuary hydrography and phytoplankton was conducted by scientists from the U.S. Summary Geological Survey in Menlo Park and U.C. Davis. Water column profiles at up to 37 stations were monitored along a transect of the Estuary run monthly between the South Bay and the Delta. The primary objective of this study was to define physical (salinity, temperature, suspended particulate matter, and light penetration), chemical (dissolved oxygen) and biological (chlorophyll a) characteristics of Estuary water that may influence other chemical and biological reactions. A second objective was to investigate planktonic indicators of ecosystem structure and function. The data collected in 1993 showed the extent and duration of the spring phytoplankton bloom in the South Bay, other localized blooms in the northern estuary, the stratification and mixing associated with the entrapment zone in the northern estuary, and mixing in the Estuary resulting from the high rainfall in 1993. Knowledge of the duration and extent of these natural features of the Estuary provide context for interpretation of the RMP contaminant data collected only 3 times per year. Another pilot study of suspended sediment transport processes was conducted by the USGS in Sacramento. This study used continuous recording sensors at Point San Pablo and the Bay Bridge to measure the amount of suspended sediment in the water at mid-depth and near the bottom, as well as tide height. The objectives of this study were to estimate which factors determine suspended solids concentrations in the Central Bay and to collect time series of suspended solids that are appropriate for continuous monitoring of suspended solids and for calibration and validation of numerical models. The investigators determined that spring tides accounted for most of the variation in suspended solids concentrations at the stations monitored, not runoff from the Sacramento or San Joaquin Rivers, or semidiurnal and diurnal tides. Comparisons were also made between measurements made by the continuous recordings and the RMP samples collected during the regular monitoring cruises. The different ways of measuring TSS were generally comparable, however only 3 measurements per year as made by the RMP could not provide the information of TSS variation actually occurring in the Estuary. This information is important because as shown by the RMP data, total contaminant concentrations in Estuary water is largely dependent on the TSS in the water. This implies that the RMP measurements alone cannot determine accurately the range of contaminant concentrations without better characterizing the dynamics of TSS. The RMP Pilot Studies are important to the developing RMP because they will help put RMP measurements into the perspective of Estuary processes and mechanisms at other time scales. The studies can relate those processes to the RMP measurements and will facilitate revision of sampling design and interpretation. Summaries of other monitoring activities pertinent to regional monitoring are also included in the Report: a description of the Regional Board’s Bay Protection Studies, the Sacramento Coordinated Monitoring Program, and a wetlands monitoring plan are included.

California↗

USGS2013

Sources and characteristics of organic matter in the Clackamas River, Oregon, related to the formation of disinfection by-products in treated drinking water

This study characterized the amount and quality of organic matter in the Clackamas River, Oregon, to gain an understanding of sources that contribute to the formation of chlorinated and brominated disinfection by-products (DBPs), focusing on regulated DBPs in treated drinking water from two direct-filtration treatment plants that together serve approximately 100,000 customers. The central hypothesis guiding this study was that natural organic matter leaching out of the forested watershed, in-stream growth of benthic algae, and phytoplankton blooms in the reservoirs contribute different and varying proportions of organic carbon to the river. Differences in the amount and composition of carbon derived from each source affects the types and concentrations of DBP precursors entering the treatment plants and, as a result, yield varying DBP concentrations and species in finished water. The two classes of DBPs analyzed in this study-trihalomethanes (THMs) and haloacetic acids (HAAs)-form from precursors within the dissolved and particulate pools of organic matter present in source water. The five principal objectives of the study were to (1) describe the seasonal quantity and character of organic matter in the Clackamas River; (2) relate the amount and composition of organic matter to the formation of DBPs; (3) evaluate sources of DBP precursors in the watershed; (4) assess the use of optical measurements, including in-situ fluorescence, for estimating dissolved organic carbon (DOC) concentrations and DBP formation; and (5) assess the removal of DBP precursors during treatment by conducting treatability "jar-test" experiments at one of the treatment plants. Data collection consisted of (1) monthly sampling of source and finished water at two drinking-water treatment plants; (2) event-based sampling in the mainstem, tributaries, and North Fork Reservoir; and (3) in-situ continuous monitoring of fluorescent dissolved organic matter (FDOM), turbidity, chlorophyll- a , and other constituents to continuously track source-water conditions in near real-time. Treatability tests were conducted during the four event-based surveys to determine the effectiveness of coagulant and powdered activated carbon (PAC) on the removal of DBP precursors. Sample analyses included DOC, total particulate carbon (TPC), total and dissolved nutrients, absorbance and fluorescence spectroscopy, and, for regulated DBPs, concentrations of THMs and HAAs in finished water and laboratory-based THM and HAA formation potentials (THMFP and HAAFP, respectively) for source water and selected locations throughout the watershed. The results of this study may not be typical given the record and near record amounts of precipitation that occurred during spring that produced streamflow much higher than average in 2010-11. Although there were algal blooms, lower concentrations of chlorophyll- a were observed in the water column during the study period compared to historical data. Concentrations of DBPs in finished (treated) water averaged 0.024 milligrams per liter (mg/L) for THMs and 0.022 mg/L for HAAs; maximum values were about 0.040 mg/L for both classes of DBPs. Although DBP concentrations were somewhat higher within the distribution system, none of the samples collected for this study or for the quarterly compliance monitoring by the water utilities exceeded levels permissible under existing U.S. Environmental Protection Agency (USEPA) regulations: 0.080 mg/L for THMs and 0.060 mg/L for HAAs. DOC concentrations were generally low in the Clackamas River, typically about 1.0-1.5 mg/L. Concentrations in the mainstem occasionally increased to nearly 2.5 mg/L during storms; DOC concentrations in tributaries were sometimes much higher (up to 7.8 mg/L). The continuous in-situ FDOM measurements indicated sharp rises in DOC concentrations in the mainstem following rainfall events; concentrations were relatively stable during summer base flow. Even though the first autumn storm mobilized appreciable quantities of carbon, higher concentrations of DBPs in finished water were observed 3-weeks later, after the ground was saturated from additional rainfall. The majority of the DOC in the lower Clackamas River appears to originate from the upper basin, suggesting terrestrial carbon was commonly the dominant source. Lower-basin tributaries typically contained the highest concentrations of DOC and DBP precursors and contributed substantially to the overall loads in the mainstem during storms. During low-flow periods, tributaries were not major sources of DOC or DBP precursors to the Clackamas River. Although the dissolved fraction of organic carbon contributed the majority of DBP precursors, at times the particulate fraction (inorganic sediment and organic particles including detritus and algal material) contributed a substantial fraction of DBP precursors. Considering just the main-stem sites, on average, 10 percent of THMFP and 32 percent of HAAFP were attributed to particulate carbon. This finding suggests water-treatment methods that remove particles prior to chlorination would reduce finished-water DBP concentrations to some degree. Overall, concentrations of THM and HAA precursors were closely linked to DOC concentrations; laboratory DBP formation potentials (DBPFPs) clearly showed that THMFP and HAAFP were greatest in the downstream tributaries that contained elevated carbon concentrations. However, carbon-normalized "specific" formation potentials for THMs and HAAs (STHMFP and SHAAFP, respectively) revealed changes in carbon character over time that affected the two types of DBP classes differently. HAA precursors were elevated in waters containing aromatic-rich soil-derived material arising from forested areas. In contrast, THM precursors were associated with carbon having a lower aromatic content; highest STHMFP occurred in autumn 2011 in the mainstem from North Fork Reservoir downstream to LO DWTP. This pattern suggests the potential for a link between THM precursors and algal-derived carbon. The highest STHMFP value was measured within North Fork Reservoir, indicating reservoir derived carbon may be important for this class of DBPs. Weak correlations between STHMFP and SHAAFP emphasize that precursor sources for these types of DBPs may be different. This highlights not only that different locations within the watershed produce carbon with different reactivity (specific DBPFP), but also that different management approaches for each class of DBP precursors could be required for control. Treatability tests conducted on source water during four basin-wide surveys demonstrated that an average of about 40 percent of DOC can be removed by coagulation. While the decrease in THMFP following coagulation was similar to DOC, the decrease in HAAFP was much greater (approximately 70 percent), indicating coagulation is particularly effective at removing HAA precursors'likely because of the aromatic nature of the carbon associated with HAA precursors. Several findings from this study have direct implications for managing drinking-water resources and for providing useful information that may help improve treatment-plant operations. For example, the use of in-situ fluorometers that measure FDOM provided an excellent proxy for DOC concentration in this system and revealed short-term, rapid changes in DOC concentration during storm events. In addition, the strong correlation between FDOM values measured in-situ and HAA5 concentrations in finished water may permit estimation of continuous HAA concentrations, as was done here. As part of this study, multiple in-situ FDOM sensors were deployed continuously and in real-time to characterize the composition of dissolved organic matter. Although the initial results were promising, additional research and engineering developments will be needed to demonstrate the full utility of these sensors for this purpose. In conclusion, although DBPFPs were strongly correlated to DOC concentration, some DBPs formed from particulate carbon, including terrestrial leaf material and algal material such as planktonic species of blue-green algae and sloughed filaments, stalks, and cells of benthic algae. Different precursor sources in the watershed were evident from the data, suggesting specific actions may be available to address some of these sources. In-situ measurements of FDOM proved to be an excellent proxy for DOC concentration as well as HAA formation during treatment, which suggests further development and refinement of these sensors have the potential to provide real-time information about complex watershed processes to operators at the drinking-water treatment plants. Follow-up studies could examine the relative roles that terrestrial and algal sources have on the DBP precursor pool to better understand how watershed-management activities may be affecting the transport of these compounds to Clackamas River drinking-water intakes. Given the low concentrations of algae in the water column during this study, additional surveys during more typical river conditions could provide a more complete understanding of how algae contribute DBP precursors. Further development of FDOM-sensor technology can improve our understanding of carbon dynamics in the river and how concentrations may be trending over time. This study was conducted in collaboration with Clackamas River Water and the City of Lake Oswego water utilities. Other research partners included Oregon Health and Science University in Hillsboro, Oregon, Alexin Laboratory in Tigard, Oregon, U.S. Geological Survey National Research Program Laboratory in Denver, Colorado, and the U.S. Geological Survey Water Science Centers in Portland, Oregon, and Sacramento, California. This project was supported with funding from Clackamas River Water, City of Lake Oswego, the U.S. Geological Survey, and the Water Research Foundation.

Oregon↗

USGS2016

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↗