USGS Science⌕ Search

SEARCH · USGS Science

Results for “Environmental Monitoring and Analysis”

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

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

717 records · Page 40Linked to original sources

Status of water quality in groundwater resources used for drinking-water supply in the southeastern San Joaquin Valley, 2013–15—California GAMA Priority Basin Project

The California Groundwater Ambient Monitoring and Assessment Program Priority Basin Project (GAMA-PBP) investigated water quality of groundwater resources used for drinking-water supplies in the Madera-Chowchilla, Kings, Kaweah, Tule, and Tulare Lake groundwater subbasins of the southeastern San Joaquin Valley during 2013–15. The study focused primarily on groundwater resources used for domestic-supply wells in the southeastern San Joaquin Valley (SESJV-D), which correspond mostly to shallower parts of aquifer systems, compared to the groundwater resources used for public-supply wells in the southeastern San Joaquin Valley (SESJV-P). The investigation had three components: (1) characterization of the status of water quality in the SESJV-D, (2) comparison between water quality in the SESJV-D and SESJV-P, and (3) identification of natural and anthropogenic factors that potentially could affect water quality in these resources. The characterization of water quality in the SESJV-D was based on data collected from 198 domestic wells sampled during 2013–15 by the U.S. Geological Survey (USGS); characterization of water quality in the SESJV-P was based on data collected from 124 wells sampled by the USGS during 2005–18 and an additional 1,577 wells with publicly available data reported to the California State Water Resources Control Board Division of Drinking Water (SWRCB-DDW). Measured concentrations were compared to regulatory and non-regulatory drinking-water quality benchmarks. A grid-based method was used to estimate the areal proportions of each study area and the whole southeastern San Joaquin Valley with high (greater than benchmark concentration), moderate (greater than half of the benchmark for inorganic and one-tenth of the benchmark for organic), and low concentrations relative to those benchmarks. Natural and anthropogenic factors that could affect groundwater quality for the SESJV-D were identified in the context of the hydrogeologic setting of the southeastern San Joaquin Valley. The considered factors represented hydrologic conditions and position in the groundwater flow system (well depth, lateral position, presence of hydric soils, percentage of coarse-grained sediment, and aridity index), land-use characteristics (percentages of agricultural, urban, and natural land use, percentage of orchard or vineyard land use, and densities of septic tanks and underground storage tanks near the wells), and geochemical conditions (groundwater age class, oxidation-reduction class, pH, and dissolved oxygen and bicarbonate concentrations). Factors are compared between SESJV-D and SESJV-P at the scale of the five study areas. One or more inorganic constituents with U.S. Environmental Protection Agency (EPA) or California maximum contaminant levels (MCLs) were detected at high concentrations in 47 percent of the SESJV-D and in 32 percent of the SESJV-P. The inorganic constituents most commonly present at high concentrations in the SESJV-D were nitrate, uranium, and arsenic. Within the SESJV-D, the proportion of the study area with high concentrations of inorganic constituents ranged from 19 percent in Madera-Chowchilla to 60 percent in Kings and Tulare Lake. One or more inorganic constituents with California State Water Resources Control Board Division of Drinking Water secondary maximum contaminant levels (SMCL-CAs) were detected at high concentrations in 14 percent of the SESJV-D and in 19 percent of the SESJV-P. The constituents most commonly present at high concentrations were manganese, iron, and total dissolved solids (TDS). Although the proportion of SESJV-D and SESJV-P with high concentrations of TDS greater than the upper SMCL were similar at 4 percent, the proportion of the SESJV-D with moderate concentrations (between the recommended and upper SMCL-CA), 30 percent, was greater than the proportion of the SESJV-P with moderate concentrations, 12 percent. One or more organic constituents with MCLs were present at high concentrations in 19 percent of the SESJV-D and in 12 percent of the SESJV-P. All the constituents detected at high concentrations in the SESJV-D were fumigants, primarily 1,2,3-trichloropropane (1,2,3-TCP) and 1,2-dibromo-3-chloropropane (DBCP). Fumigants also were the constituents most commonly detected at high concentrations in the SESJV-P, although high concentrations of solvents also were detected. The SESJV-D dataset included analysis of many organic constituents without MCL benchmarks and with detection levels far below drinking water benchmark concentrations; detections at these low concentrations can be used as tracers of anthropogenic influence on groundwater. Pesticides and degradates of pesticides were detected in 60 percent of the SESJV-D; the most frequently detected pesticides were the herbicides simazine, didealkylatrazine (CAAT, a degradate of simazine and atrazine), diuron, and bromacil.

California↗

Environmental presence and persistence of pharmaceuticals: An overview

Emerging contaminants (ECs) in the environment – that is, chemicals with domestic, municipal, industrial, or agricultural sources that are not commonly monitored but may have the potential for adverse environmental effects – is a rapidly growing field of research. The use of “emerging” is not intended to infer that the presence of these compounds in the environment is new. These chemicals have been released into the environment as long as they have been in production or, in the case of hormones and other endogenous compounds, since the rise of animal life. What is emerging is the interest by the scientific and lay communities in the presence of these chemicals in the environment, the analytical capabilities required for detection, and the subtle effects that very small concentrations of these chemicals appear to have on aquatic biota. In December 2006, Environmental Science & Technology devoted an entire special issue (volume 40, number 23) to the topic of ECs, illustrating the increased interest in the subject. Within the EGs, one particular class that has seen a substantial increase in research over the past 10 years is pharmaceuticals and personal-care products (PPCPs). This increased research interest can be demonstrated by several means, including requests for proposals from funding agencies, but the clearest indication of a focused effort to understand the introduction, transformation, and potential health and environmental effects of PPCPs and ECs, in general, is the number of published reports. This increase can be shown by examining six environmental journals that regularly publish PPCP-related papers – Chemosphere , Environmental Science & Technology , Environmental Toxicology and Chemistry , Science of the Total Environment , Water Research , and Water Science and Technology . In 1998 there were 22 papers published on pharmaceuticals, antibiotics, or drugs in these 6 journals; by 2006, this number increased sixfold to 132 papers (Figure 1.1). This growth can be attributed to a number of factors. The presence of pharmaceuticals in surface-water samples from Europe and the United States was documented in several sentinel papers. These ground-breaking works encouraged other scientists to examine the rivers, streams, lakes, and reservoirs in their regions for such chemicals. In addition, the intense public attention paid to news reports on the environmental detections of these chemicals and possible effects of aquatic life has made this issue visible to the wastewater-treatment, drinking-water treatment, and regulatory communities. This has driven the funding bodies associated with these communities to fund studies or request proposals that address the presence, fate, and effects of PPCPs in aquatic systems. The release of the first comprehensive reconnaissance of pharmaceuticals and other wastewater contaminants in the United States provides an example of the intense media interest in this topic. Within 6 days on online publication of this study, 72 newspapers across the United States had published articles describing the results, either locally written or based on international media syndicate reports. There also was substantial concurrent coverage by local and national radio and television outlets, including the Cable News Network, ABC World News Tonight , and National Public Radio. A substantial fraction of these news stories may be attributable to press releases and media briefings prior to publication. However, the interest by television and print journalists in reporting the results of a peer-reviewed journal article to the general public was motivated by the recognition that describing the presence of PPCPs in water supplies would be of interest to the public. To better convey the results of the study published by Koplin et al. to the public, a separate general-interest fact sheet was published to summarize the important points of the study. Because PPCPs are commonly and widely used by individuals, there is likely a preexisting, personal identification with these compounds that does not occur for the wide range of other organic and inorganic contaminants whose presence in the environment has previously been described. This greater public “name recognition” makes itself known through the media to the regulatory and technical community and has prompted interest in sponsoring research that defines the composition and concentrations of PPCPs in potential sources and their fate and effects following relase into the environment. Independent of the drivers that potentially fuel the interest in studies of PPCPs, it is clear that PPCP research has grown beyond surface-water studies to examine issues such as: • Presence in other matrices, such as groundwater, landfill leachates, sediments, and biosolids. • Environmental transport and fate in surface water, groundwater, and soils amended with reclaimed water or biosolids. • PPCP source elucidation, such as wastewater treatment plant (WWTP) effluents, confined animal feeding operations (CAFOs), and aquaculture. • Removal during wastewater and drinking-water treatment. • Effects on aquatic ecosystems, terrestrial ecosystems, and human health. The chapters in this book provide an extensive examination of current environmental pharmaceutical research and are divided into three sections: “Occurrence and Analysis of Pharmaceuticals in the Environment,” “Environment Fate and Transformations of Veterinary Pharmaceuticals,” and “treatment of Pharmaceuticals in Drinking Water and Wastewater.” The purpose of this introductory overview chapter is to outline current (2004-2006) knowledge about the presence and concentration of PPCPs as described in the published literature. Previous reviews should be consulted for discussions on pre-2004 publications. Those reviews will provide the reader with a comprehensive introduction to the topic of PPCPs in the environment. This chapter describes the sources of PPCPs and other organic contaminants often associated with human wastewater into the environment, the range of concentrations present in various environmental compartments, and the potential routes of removal/sequestration. An overview of the sources and fate of veterinary pharmaceuticals will be discussed in Chapter 5, “Fate and Transport of Veterinary Medicines in the Soil Environment.”

Book chapter↗

A spatially referenced regression model (SPARROW) for suspended sediment in streams of the Conterminous U.S.

Suspended sediment has long been recognized as an important contaminant affecting water resources. Besides its direct role in determining water clarity, bridge scour and reservoir storage, sediment serves as a vehicle for the transport of many binding contaminants, including nutrients, trace metals, semi-volatile organic compounds, a nd numerous pesticides (U.S. Environmental Protection Agency, 2000a). Recent efforts to addr ess water-quality concerns through the Total Maximum Daily Load (TMDL) process have iden tified sediment as the single most prevalent cause of impairment in the Nation’s streams a nd rivers (U.S. Environmental Protection Agency, 2000b). Moreover, sediment has been identified as a medium for the tran sport and sequestration of organic carbon, playing a potentia lly important role in understa nding sources and sinks in the global carbon budget (Stallard, 1998). A comprehensive understanding of sediment fate a nd transport is considered essential to the design and implementation of effective plans for sediment management (Osterkamp and others, 1998, U.S. General Accounting Office, 1990). An exte nsive literature addr essing the problem of quantifying sediment transport has produced a nu mber of methods for estimating its flux (see Cohn, 1995, and Robertson and Roerish, 1999, for us eful surveys). The accuracy of these methods is compromised by uncertainty in the concentration measurements and by the highly episodic nature of sediment movement, particul arly when the methods are applied to smaller basins. However, for annual or decadal flux es timates, the methods are generally reliable if calibrated with extended periods of data (Robertson and Roerish, 1999). A substantial literature also supports the Universal Soil Loss Equation (U SLE) (Soil Conservation Service, 1983), an engineering method for estimating sheet and rill erosion, although the empirical credentials of the USLE have recently been questioned (Tri mble and Crosson, 2000). Conversely, relatively little direct evidence is available concerning the fate of sediment. The common practice of quantifying sediment fate with a sediment deliv ery ratio, estimated from a simple empirical relation with upstream basin area, does not artic ulate the relative importance of individual storage sites within a basin (Wolman, 1977). Rates of sediment deposition in reservoirs and flood plains can be determined from empirical measurement s , but only a limited number of sites have been monitored, and net rates of deposition or loss from other potential sinks and sources is largely unknown (Stallard, 1998). In particular, little is known about how much sediment loss from fields ultimately makes its way to stream channels, and how much sediment is subsequently stored in or lost from th e streambed (Meade and Parker, 1985, Trimble and Crosson, 2000). This paper reports on recent progress made to a ddress empirically the question of sediment fate and transport on a national scale. The model pres ented here is based on the SPAtially Referenced Regression On Watershed attr ibutes (SPARROW) methodology, fi rst used to estimate the distribution of nutrients in str eams and rivers of the United Stat es, and subsequently shown to describe land and stream processes affecting the delivery of nutrients (Smith and others, 1997, Alexander and others, 2000, Preston and Brakeb ill, 1999). The model makes use of numerous spatial datasets, available at the national level, to explain long-term sediment water-quality conditions in major streams and rivers throughou t the United States. Sediment sources are identified using sediment erosion rates from the National Resources I nventory (NRI) (Natural Resources Conservation Service, 2000) and apportioned over the landscape according to 30- meter resolution land-use information from th e National Land Cover Data set (NLCD) (U.S. Geological Survey, 2000a). More than 76,000 reservoirs from the National Inventory of Dams (NID) (U.S. Army Corps of Engin eers, 1996) are identified as pot ential sediment sinks. Other, non-anthropogenic sources and sinks are identified using soil in formation from the State Soil Survey Geographic (STATSGO) data base (Schwarz and Alexander, 1995) and spatial coverages representing surficial rock t ype and vegetative cover. The SPA RROW model empirically relates these diverse spatial datasets to estimates of long-term, mean annual sediment flux computed from concentration and flow measurements co llected over the period 1985 -95 from more than 400 monitoring stations maintained by the Na tional Stream Quality Accounting Network (Alexander and others, 1998), the National Wa ter Quality Assessment Program, and U.S. Geological Survey District offices (Turcios and Gray, in press). Th e calibrated model is used to estimate sediment flux for over 60,000 stream segments included in the River Reach File 1 (RF1) stream network (Alexander and others, 1999). SPARROW uses statis tical methods to calibrate a simple, structural model of riverine water quality, one that imposes mass ba lance in accounting for changes in contaminant flux. As applied here, the mass-balance approach facilitates the interpretation of model results in terms of physical processes affecting sediment transport, and makes possible the estimation of various rates of sediment generation and loss associated with stream channels and features of the landscape. The statistical approach provides a basi s for assessing the error of these inferred rates and of the error in extrapolated estimates of sediment flux made for streams in the RF1 network. An important implication of the holistic modeling approach adopted in this analysis is that estimates of sediment production and loss ar e based on, and therefore consistent with, measurements of in-stream flux. Other ancillary information, such as direct measurements of long-term sediment storage and release from rese rvoirs (Steffen, 1996), is incorporated into the analysis by specifying additional equations expl aining these ancillary variables. The imposition of cross-equation constraints affords this info rmation a statistically consistent weight in explaining in-stream sediment flux. Thus, the me thodology described here represents a general framework for synthesizing a wide spectrum of available information relevant to the understanding of sediment fate and transport.

Conterminous United States↗

Historical methane hydrate project review

In 1995, U.S. Geological Survey made the first systematic assessment of the volume of natural gas stored in the hydrate accumulations of the United States. That study, along with numerous other studies, has shown that the amount of gas stored as methane hydrates in the world greatly exceeds the volume of known conventional gas resources. However, gas hydrates represent both a scientific and technical challenge and much remains to be learned about their characteristics and occurrence in nature. Methane hydrate research in recent years has mostly focused on: (1) documenting the geologic parameters that control the occurrence and stability of gas hydrates in nature, (2) assessing the volume of natural gas stored within various gas hydrate accumulations, (3) analyzing the production response and characteristics of methane hydrates, (4) identifying and predicting natural and induced environmental and climate impacts of natural gas hydrates, and (5) analyzing the effects of methane hydrate on drilling safety. Methane hydrates are naturally occurring crystalline substances composed of water and gas, in which a solid water-­‐lattice holds gas molecules in a cage-­‐like structure. The gas and water becomes a solid under specific temperature and pressure conditions within the Earth, called the hydrate stability zone. Other factors that control the presence of methane hydrate in nature include the source of the gas included within the hydrates, the physical and chemical controls on the migration of gas with a sedimentary basin containing methane hydrates, the availability of the water also included in the hydrate structure, and the presence of a suitable host sediment or “reservoir”. The geologic controls on the occurrence of gas hydrates have become collectively known as the “methane hydrate petroleum system”, which has become the focus of numerous hydrate research programs. Recognizing the importance of methane hydrate research and the need for a coordinated effort, the U.S. Congress enacted Public Law 106-­‐193, the Methane Hydrate Research and Development Act of 2000. This Act called for the Secretary of Energy to begin a methane hydrate research and development program in consultation with other U.S. federal agencies. At the same time a new methane hydrate research program had been launched in Japan by the Ministry of International Trade and Industry to develop plans for a methane hydrate exploratory drilling project in the Nankai Trough. Since this early start we have seen other countries including India, China, Canada, and the Republic of Korea establish large gas hydrate research and development programs. These national led efforts have also included the investment in a long list of important scientific research drilling expeditions and production test studies that have provided a wealth of information on the occurrence of methane hydrate in nature. The most notable expeditions and projects have including the following: -­‐Ocean Drilling Program Leg 164 (1995) -­‐Japan Nankai Trough Project (1999-­‐2000) -­‐Ocean Drilling Program Leg 204 (2004) -­‐Japan Tokai-­‐oki to Kumano-­‐nada Project (2004) -­‐Gulf of Mexico JIP Leg I (2005) -­‐Integrated Ocean Drilling Program Expedition 311 (2005) -­‐Malaysia Gumusut-­‐Kakap Project (2006) -­‐India NGHP Expedition 01 (2006) -­‐China GMGS Expedition 01 (2007) -­‐Republic of Korea UBGH Expedition 01 (2007) -­‐Gulf of Mexico JIP Leg II (2009) -­‐Republic of Korea UBGH Expedition 02 (2010) -­‐MH-­‐21 Nankai Trough Pre-­‐Production Expedition (2012-­‐2013) -­‐Mallik Gas Hydrate Testing Projects (1998/2002/2007-­‐2008) -­‐Alaska Mount Elbert Stratigraphic Test Well (2007) -­‐Alaska Iġnik Sikumi Methane Hydrate Production Test Well (2011-­‐2012) Research coring and seismic programs carried out by the Ocean Drilling Program (ODP) and Integrated Ocean Drilling Program (IODP), starting with the ODP Leg 164 drilling of the Blake Ridge in the Atlantic Ocean in 1995, have also contributed greatly to our understanding of the geologic controls on the formation, occurrence, and stability of gas hydrates in marine environments. For the most part methane hydrate research expeditions carried out by the ODP and IODP provided the foundation for our scientific understanding of gas hydrates. The methane hydrate research efforts under ODP-­‐IODP have mostly dealt with the assessment of the geologic controls on the occurrence of gas hydrate, with a specific goal to study the role methane hydrates may play in the global carbon cycle. Over the last 10 years, national led methane hydrate research programs, along with industry interest have led to the development and execution of major methane hydrate production field test programs. Two of the most important production field testing programs have been conducted at the Mallik site in the Mackenzie River Delta of Canada and in the Eileen methane hydrate accumulation on the North Slope of Alaska. Most recently we have also seen the completion of the world’s first marine methane hydrate production test in the Nankai Trough in the offshore of Japan. Industry interest in gas hydrates has also included important projects that have dealt with the assessment of geologic hazards associated with the presence of hydrates. The scientific drilling and associated coring, logging, and borehole monitoring technologies developed in the long list of methane hydrate related field studies are one of the most important developments and contributions associated with methane hydrate research and development activities. Methane hydrate drilling has been conducted from advanced scientific drilling platforms like the JOIDES Resolution and the D/V Chikyu, which feature highly advanced integrated core laboratories and borehole logging capabilities. Hydrate research drilling has also included the use of a wide array of industry, geotechnical and multi-­‐service ships. All of which have been effectively used to collect invaluable geologic and engineering data on the occurrence of methane hydrates throughout the world. Technologies designed specifically for the collection and analysis of undisturbed methane hydrate samples have included the development of a host of pressure core systems and associated specialty laboratory apparatus. The study and use of both wireline conveyed and logging-­‐while-­‐drilling technologies have also contributed greatly to our understanding of the in-­‐situ nature of hydrate-­‐bearing sediments. Recent developments in borehole instrumentation specifically designed to monitor changes associated with hydrates in nature through time or to evaluate the response of hydrate accumulations to production have also contributed greatly to our understanding of the complex nature and evolution of methane hydrate systems. Our understanding of how methane hydrates occur and behave in nature is still growing and evolving – we do not yet know if methane hydrates can be economically produced, nor do we know fully the role of hydrates as an agent of climate change or as a geologic hazard. But it is known for certain that scientific drilling has contributed greatly to our understanding of hydrates in nature and will continue to be a critical source of the information to advance our understanding of methane hydrates.

Report↗

Groundwater-quality data in seven GAMA study units: results from initial sampling, 2004-2005, and resampling, 2007-2008, of wells: California GAMA Program Priority Basin Project

The Priority Basin Project (PBP) of the Groundwater Ambient Monitoring and Assessment (GAMA) Program was developed in response to the Groundwater Quality Monitoring Act of 2001 and is being conducted by the U.S. Geological Survey (USGS) in cooperation with the California State Water Resources Control Board (SWRCB). The GAMA-PBP began sampling, primarily public supply wells in May 2004. By the end of February 2006, seven (of what would eventually be 35) study units had been sampled over a wide area of the State. Selected wells in these first seven study units were resampled for water quality from August 2007 to November 2008 as part of an assessment of temporal trends in water quality by the GAMA-PBP. The initial sampling was designed to provide a spatially unbiased assessment of the quality of raw groundwater used for public water supplies within the seven study units. In the 7 study units, 462 wells were selected by using a spatially distributed, randomized grid-based method to provide statistical representation of the study area. Wells selected this way are referred to as grid wells or status wells. Approximately 3 years after the initial sampling, 55 of these previously sampled status wells (approximately 10 percent in each study unit) were randomly selected for resampling. The seven resampled study units, the total number of status wells sampled for each study unit, and the number of these wells resampled for trends are as follows, in chronological order of sampling: San Diego Drainages (53 status wells, 7 trend wells), North San Francisco Bay (84, 10), Northern San Joaquin Basin (51, 5), Southern Sacramento Valley (67, 7), San Fernando–San Gabriel (35, 6), Monterey Bay and Salinas Valley Basins (91, 11), and Southeast San Joaquin Valley (83, 9). The groundwater samples were analyzed for a large number of synthetic organic constituents (volatile organic compounds [VOCs], pesticides, and pesticide degradates), constituents of special interest (perchlorate, N-nitrosodimethylamine [NDMA], and 1,2,3-trichloropropane [1,2,3-TCP]), and naturally-occurring inorganic constituents (nutrients, major and minor ions, and trace elements). Naturally-occurring isotopes (tritium, carbon-14, and stable isotopes of hydrogen and oxygen in water) also were measured to help identify processes affecting groundwater quality and the sources and ages of the sampled groundwater. Nearly 300 constituents and water-quality indicators were investigated. Quality-control samples (blanks, replicates, and samples for matrix spikes) were collected at 24 percent of the 55 status wells resampled for trends, and the results for these samples were used to evaluate the quality of the data for the groundwater samples. Field blanks rarely contained detectable concentrations of any constituent, suggesting that contamination was not a noticeable source of bias in the data for the groundwater samples. Differences between replicate samples were mostly within acceptable ranges, indicating acceptably low variability in analytical results. Matrix-spike recoveries were within the acceptable range (70 to 130 percent) for 75 percent of the compounds for which matrix spikes were collected. This study did not attempt to evaluate the quality of water delivered to consumers. After withdrawal, groundwater typically is treated, disinfected, and blended with other waters to maintain acceptable water quality. The benchmarks used in this report apply to treated water that is served to the consumer, not to untreated groundwater. To provide some context for the results, however, concentrations of constituents measured in these groundwater samples were compared with benchmarks established by the U.S. Environmental Protection Agency (USEPA) and California Department of Public Health (CDPH). Comparisons between data collected for this study and benchmarks for drinking water are for illustrative purposes only and are not indicative of compliance or non-compliance with those benchmarks. Most constituents that were detected in groundwater samples from the trend wells were found at concentrations less than drinking-water benchmarks. Four VOCs—trichloroethene, tetrachloroethene, 1,2-dibromo-3-chloropropane, and methyl tert-butyl ether—were detected in one or more wells at concentrations greater than their health-based benchmarks, and six VOCs were detected in at least 10 percent of the samples during initial sampling or resampling of the trend wells. No pesticides were detected at concentrations near or greater than their health-based benchmarks. Three pesticide constituents—atrazine, deethylatrazine, and simazine—were detected in more than 10 percent of the trend-well samples during both sampling periods. Perchlorate, a constituent of special interest, was detected more frequently, and at greater concentrations during resampling than during initial sampling, but this may be due to a change in analytical method between the sampling periods, rather than to a change in groundwater quality. Another constituent of special interest, 1,2,3-TCP, was also detected more frequently during resampling than during initial sampling, but this pattern also may not reflect a change in groundwater quality. Samples from several of the wells where 1,2,3-TCP was detected by low-concentration-level analysis during resampling were not analyzed for 1,2,3-TCP using a low-level method during initial sampling. Most detections of nutrients and trace elements in samples from trend wells were less than health-based benchmarks during both sampling periods. Exceptions include nitrate, arsenic, boron, and vanadium, all detected at concentrations greater than their health-based benchmarks in at least one well during both sampling periods, and molybdenum, detected at concentrations greater than its health-based benchmark during resampling only. The isotopic ratios of oxygen and hydrogen in water and tritium and carbon-14 activities generally changed little between sampling periods, suggesting that the predominant sources and ages of groundwater in most trend wells were consistent between the sampling periods.

California↗

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↗

Genetic integrity, population status, and long-term viability of isolated populations of shoal bass in the upper Chattahoochee River basin, Georgia

This report characterizes the status of multiple isolated Shoal Bass ( Micropterus cataractae ) populations in the upper Chattahoochee River basin (UCRB), Georgia. The Shoal Bass, a sport fish endemic to the Apalachicola-Chattahoochee-Flint River (ACF) basin, is a fluvial-specialist species considered vulnerable to local extirpations and extinction due to habitat fragmentation and introgression with non-native congeners. Perhaps one of the most isolated populations of Shoal Bass exists in a 2-km reach of Big Creek, a tributary of the Chattahoochee River located near Roswell, Georgia. Big Creek is partially contained within the Chattahoochee River National Recreation Area, although the Big Creek watershed is riddled with urban land cover. Roswell Mill Dam limits the upstream extent of the Shoal Bass population at Big Creek, and the downstream extent is presumably limited to the confluence of Big Creek and the Chattahoochee River. This reach of the Chattahoochee River is thermally depressed because of coldwater releases from Lake Lanier, and is considered unsuitable for Shoal Bass. Herein, we examine the genetic integrity, population status, and long-term viability of the Shoal Bass population in Big Creek. We also examine two additional Shoal Bass populations that occur in the UCRB, specifically the Chestatee River and the upper Chattahoochee River, both of which are impounded at Lake Lanier. Together, the Shoal Bass inhabiting these three stream systems comprise a distinct genetic stock of Shoal Bass (Taylor 2017), underscoring the importance of conserving these populations towards maintaining the overall diversity and adaptive potential of the species. We assessed genetic diversity and estimated effective population sizes within these three rivers by genotyping fish with 16 microsatellite DNA markers. Results demonstrated that the Shoal Bass population in Big Creek has experienced high rates of introgression with non-native Smallmouth Bass (M. dolomieu), purportedly introduced into the Chattahoochee River in the past 10-15 years. Alarmingly, only 24% (15 of 62) of putative Shoal Bass collected from Big Creek were genetically pure Shoal Bass, whereas the majority of fish were first-filial (F1) generation hybrids and unidirectional backcrosses towards Shoal Bass. Fleeting opportunity may remain to conserve the native genome of the Shoal Bass population in Big Creek. High hybridization rates prevented genetic diversity analysis for the Big Creek population. Shoal Bass populations in the Chestatee and Chattahoochee rivers displayed levels of genetic diversity similar to populations that persist in other rivers in the ACF basin, namely the Flint and Chipola rivers. Effective population sizes of 93.8– 197.4 for the Chestatee and Chattahoochee rivers (combined) suggest that the conservation status of these populations is stable for the short-term, but may be at risk of losing genetic diversity and adaptive potential in the long-term. To estimate age and mortality of the three populations, we used fish scales and capture-markrecapture (CMR) as complementary, non-lethal methods for age estimation. Estimated ages of phenotypic Shoal Bass ranged from 1-12 years in all three populations, demonstrating increased longevity compared to populations elsewhere within the native range. Catch-curve estimates of annual mortality ranged from 18.4-23.7%, which are markedly lower than those observed in other Shoal Bass populations in the ACF basin. These differences in life-history characteristics underscore the need for the development of population-specific management and conservation strategies for Shoal Bass in the UCRB. The lowest recruitment variability (i.e., the variation in year-class strength) was observed in the Chestatee River, a forested watershed, whereas the highest variability was observed in Big Creek, an urbanized watershed. Recruitment strength in Big Creek was negatively influenced by discharge variability in the summer months, suggesting that flashy, sediment-laden flows hinder survival of recently hatched young. Other statistically significant models from Big Creek and the Chattahoochee River indicated that over-winter survival could be an important pinch-point for recruitment in UCRB populations. A multi-agency sampling effort was conducted from May 2013-May 2016 to estimate the population size of Shoal Bass occupying the 1-km of wadeable shoal habitats in Big Creek. Using CMR models, we estimated that approximately 219-348 Shoal Bass (≥ 70 mm total length) occupied the area throughout the duration of our study. These estimates largely reflect abundance of individuals aged 0-2 years, as only 9% (36 of 408) tagged fish were aged ≥ 7 years. Local abundance appeared similar to that reported for a population that inhabited Little Uchee Creek, a similar-sized tributary of the Chattahoochee River, prior to its recent functional extirpation. The low abundance of large, adult Shoal Bass further suggests the long-term viability of the Big Creek population may be in jeopardy. Perhaps most importantly, CMR estimates reflect abundance of phenotypic Shoal Bass – genetic analyses suggest the abundance of pure Shoal Bass could be an order of magnitude smaller. To evaluate the potential for adult Shoal Bass to emigrate from Big Creek into the mainstem Chattahoochee River, we tagged eight adults with acoustic telemetry tags and assessed their seasonal residency at two stationary receiver locations located in increasing proximity to the confluence with the Chattahoochee River. Fish took up residency near the confluence during the fall and winter months, during which time water temperatures in Big Creek were periodically colder than the Chattahoochee River. Although we were unable to document emigration, we conclude that the potential for emigration is highest during the winter months when the Chattahoochee River may be warmer than Big Creek. Two of the tagged fish were caught by anglers near the confluence, suggesting that angling pressure at Big Creek may be higher than previously suspected. Overall, this study observed unique life-history characteristics and characterized the population status of multiple Shoal Bass populations in the UCRB. Populations in the Chestatee and Chattahoochee rivers appear stable at present and likely represent the last remaining strongholds for pure Shoal Bass in the UCRB. Efforts to preserve forested watershed conditions, natural hydrology, and shoal habitats would contribute to the long-term persistence of Shoal Bass populations in these two rivers. Additionally, the detection of non-native Alabama Bass and their associated hybrids in both rivers is cause for concern. Diligent monitoring of hybridization dynamics between Alabama Bass and Shoal Bass is warranted, along with an assessment of Alabama Bass invasion extent upstream of Lake Lanier. The Shoal Bass population in Big Creek is threatened by elevated levels of introgression with nonnative Smallmouth Bass, recruitment variability, low abundance of adults, and isolation from other populations. Conservation intervention is urgently needed to restore and preserve this genetically distinct population, which would contribute to preservation of range wide genetic diversity and adaptability of the species. Additionally, an urban sport fishery for Shoal Bass at Big Creek has the potential to serve as a tool for increasing public awareness, engagement, and support of Shoal Bass conservation efforts in the UCRB. We suggest strategies for conservation of the remnant shoal habitats and Shoal Bass population in Big Creek, including potential development of a supplemental stocking program, selective removal of non-native congeners, and delivery of environmental education programs that could bolster awareness and appreciation.

Georgia↗

Response curves for phosphorus plume lengths from reactive-solute-transport simulations of onland disposal of wastewater in noncarbonate sand and gravel aquifers

Surface-water resources in Massachusetts often are affected by eutrophication, excessive plant growth, which has resulted in impaired use for a majority of the freshwater ponds and lakes and a substantial number of river-miles in the State. Because supply of phosphorus usually is limiting to plant growth in freshwater systems, control of phosphorus input to surface waters is critical to solving the impairment problem. Wastewater is a substantial source of phosphorus for surface water, and removal of phosphorus before disposal may be necessary. Wastewater disposed onland by infiltration loses phosphorus from the dissolved phase during transport through the subsurface and may be an effective disposal method; quantification of the phosphorus loss can be simulated to determine disposal feasibility. In 2003, the U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection, initiated a project to simulate distance of phosphorus transport in the subsurface for plausible conditions of onland wastewater disposal and subsurface properties. A coupled one-dimensional unsaturated-zone and three-dimensional saturated-zone reactive-solute-transport model (PHAST) was used to simulate lengths of phosphorus plumes. Knowledge of phosphorus plume length could facilitate estimates of setback distances for wastewater-infiltration sites from surface water that would be sufficient to protect the surface water from eutrophication caused by phosphorus transport through the subsurface and ultimate discharge to surface water. The reactive-solute-transport model PHAST was used to simulate ground-water flow, solute transport, equilibrium chemistry for dissolved and sorbed species, and kinetic regulation of organic carbon decomposition and phosphate mineral formation. The phosphorus plume length was defined for the simulations as the maximum extent of the contour for the 0.015 milligram-per-liter concentration of dissolved phosphorus downgradient from the infiltration bed after disposal cessation. Duration of disposal before cessation was assumed to be 50 years into an infiltration bed of 20,000 square feet at the rate of 3 gallons per square foot per day. Time for the maximum extent of the phosphorus plume to develop is on the order of 100 years after disposal cessation. Simulations indicated that phosphorus transport beyond the extent of the 0.015 milligram-per-liter concentration contour was never more than 0.18 kilogram per year, an amount that would likely not alter the ecology of most surface water. Simulations of phosphorus plume lengths were summarized in a series of response curves. Simulated plume lengths ranged from 200 feet for low phosphorus-concentration effluents (0.25 milligram per liter) and thick (50 feet) unsaturated zones to 3,400 feet for high phosphorus-concentration effluents (14 milligrams per liter) discharged directly into the aquifer (unsaturated-zone thickness of 0 feet). Plume length was nearly independent of unsaturated-zone thickness at phosphorus concentrations in the wastewater that were less than 2 milligrams per liter because little or no phosphorus mineral formed at low phosphorus concentrations. For effluents of high phosphorus concentration, plume length varied from 3,400 feet for unsaturated-zone thickness of 0 to 2,550 feet for unsaturated-zone thickness of 50 feet. Model treatments of flow and equilibrium-controlled chemistry likely were more accurate than rates of kinetically controlled reactions, notably precipitation of iron-phosphate minerals; the kinetics of such reactions are less well known and thus less well defined in the model. Sensitivity analysis indicated that many chemical and physical aquifer properties, such as hydraulic gradient and model width, did not affect the simulated plume length appreciably, but duration of discharge, size of infiltration bed, amount of dispersion, and number of sorption sites on the aquifer sediments did affect plume length appreciably. Because simulation of plume length in carbonate-mineral sediments indicated that the plume would be substantially longer than in noncarbonate-mineral sediments, the application of the response curves in locations with carbonate-mineral sediments would be inappropriate. The effect of carbonate minerals in sediments is to increase pH, which causes decreased sorption of phosphorus on aquifer sediments. Phosphorus removal from solution by precipitation onto aquifer sediments is more efficient at high concentrations of disposed phosphorus than at low concentrations. At very low phosphorus concentrations, the solubility product of phosphorus minerals is not exceeded and no phosphorus mineral forms. An important consequence is that removal of dissolved phosphorus from the plume by processes in the subsurface is decreased the more that removal efforts are applied in treatment before wastewater is disposed. Model simulations indicate that removal of phosphorus from wastewater disposed through septic systems would have the advantage of efficient phosphorus removal in the subsurface because phosphorus concentrations are high in septic-system effluent. Short plume lengths result from wastewater disposal through septic systems because of the efficient phosphorus removal and because of the low volume of wastewater involved. The simulation results for small-volume systems are not quantitative, however, because wastewater-infiltration rates are much lower than those of the higher-volume system that was used to calibrate the model and to create the plume-length response curves. The response curves for phosphorus plume lengths, as defined by the maximum extent of the 0.015 milligram-per-liter concentration contour, is clearly defined in the model simulations, although the relation between simulated plume length and protective setback distance is subject to interpretation. Phosphorus does move beyond the point at which the simulated 0.015 milligram-per-liter concentration contour has stopped, so that a determination of protective setback distance must include a consideration of whether that continued flux, or some other flux amount, is appropriate. Also, simulations indicate that phosphorus plumes do not reach their full extent until 50 to 200 years after disposal cessation, depending on concentration of phosphorus disposed. No phosphorus plume has been monitored for that long after cessation, so there is no way to verify the long-term simulation results.

Massachusetts↗

Water quality and quantity and simulated surface-water and groundwater flow in the Laurel Hill Creek Basin, southwestern Pennsylvania, 1991–2007

Laurel Hill Creek is considered one of the most pristine waterways in southwestern Pennsylvania and has high recreational value as a high-quality cold-water fishery; however, the upper parts of the basin have documented water-quality impairments. Groundwater and surface water are withdrawn for public water supply and the basin has been identified as a Critical Water Planning Area (CWPA) under the State Water Plan. The U.S. Geological Survey, in cooperation with the Somerset County Conservation District, collected data and developed modeling tools to support the assessment of water-quality and water-quantity issues for a basin designated as a CWPA. Streams, springs, and groundwater wells were sampled for water quality in 2007. Streamflows were measured concurrent with water-quality sampling at main-stem sites on Laurel Hill Creek and tributaries in 2007. Stream temperatures were monitored continuously at five main-stem sites from 2007 to 2010. Water usage in the basin was summarized for 2003 and 2009 and a Water-Analysis Screening Tool (WAST) developed for the Pennsylvania State Water Plan was implemented to determine whether the water use in the basin exceeded the “safe yield” or “ the amount of water that can be withdrawn from a water resource over a period of time without impairing the long-term utility of a water resource .” A groundwater and surface-water flow (GSFLOW) model was developed for Laurel Hill Creek and calibrated to the measured daily streamflow from 1991 to 2007 for the streamflow-gaging station near the outlet of the basin at Ursina, Pa. The CWPA designation requires an assessment of current and future water use. The calibrated GSFLOW model can be used to assess the hydrologic effects of future changes in water use and land use in the basin. Analyses of samples collected for surface-water quality during base-flow conditions indicate that the highest nutrient concentrations in the main stem of Laurel Hill Creek were at sites in the northeastern part of the basin where agricultural activity is prominent. All of the total nitrogen (N) and a majority of the total phosphorus (P) concentrations in the main stem exceeded regional nutrient criteria levels of 0.31 and 0.01 milligrams per liter (mg/L), respectively. The highest total N and total P concentrations in the main stem were 1.42 and 0.06 mg/L, respectively. Tributary sites with the highest nutrient concentrations are in subbasins where treated wastewater is discharged, such as Kooser Run and Lost Creek. The highest total N and total P concentrations in subbasins were 3.45 and 0.11 mg/L, respectively. Dissolved chloride and sodium concentrations were highest in the upper part of the basin downstream from Interstate 76 because of road deicing salts. The mean base-flow concentrations of dissolved chloride and sodium were 117 and 77 mg/L, respectively, in samples from the main stem just below Interstate 76, and the mean concentrations in Clear Run were 210 and 118 mg/L, compared to concentrations less than 15 mg/L in tributaries that were not affected by highway runoff. Water quality in forested tributary subbasins underlain by the Allegheny and Pottsville Formations was influenced by acidic precipitation and, to a lesser extent, the underlying geology as indicated by pH values less than 5.0 and corresponding specific conductance ranging from 26 to 288 microsiemens per centimeter at 25 degrees Celsius for some samples; in contrast, pH values for main stem sites ranged from 6.6 to 8.5. Manganese (Mn) was the only dissolved constituent in the surface-water samples that exceeded the secondary maximum contaminant level (SMCL). More than one-half the samples from the main stem had Mn concentrations exceeding the SMCL level of 50 micrograms per liter (μg/L), whereas only 19 percent of samples from tributaries exceeded the SMCL for Mn. Stream temperatures along the main stem of Laurel Hill Creek became higher moving downstream. During the summer months of June through August, the daily mean temperatures at the five sites exceeded the limit of 18.9 degrees Celsius (°C) for a cold-water fishery. The maximum instantaneous values for each site ranged from 27.2 to 32.8 °C. Water-quality samples collected at groundwater sites (wells and springs) indicate that wells developed within the Mauch Chunk Formation had the best water quality, whereas wells developed within the Allegheny and Pottsville Formations yielded the poorest water quality. Waters from the Mauch Chunk Formation had the highest median pH (7.6) and alkalinity (80 mg/L calcium carbonate) values. The lowest pH and alkalinity median values were in waters from the Allegheny and Pottsville Formations. Groundwater samples collected from wells in the Allegheny and Pottsville Formations also had the highest concentrations of dissolved iron (Fe) and dissolved Mn. Seventy-eight percent of the groundwater samples collected from the Allegheny Formation exceeded the SMCL of 300 μg/L for Fe and 50 μg/L for Mn. Forty-three and 62 percent of the groundwater samples collected from the Pottsville Formation exceeded the SMCL for iron and Mn, respectively. The highest Fe and Mn concentrations for surface waters were measured for tributaries draining the Pottsville Formation. The highest median Fe concentration for tributaries was in samples from streams draining the Allegheny Formation. During base-flow conditions, the streamflow per unit area along the main stem of Laurel Hill Creek was lowest in the upper parts of the basin [farthest upstream site 0.07 cubic foot per second per square mile (ft 3 /s/mi 2 )] and highest (two sites averaging about 0.20 (ft 3 /s/mi 2 ) immediately downstream from Laurel Hill Lake in the center of the basin. Tributaries with the highest streamflow per unit area were those subbasins that drain the western ridge of the Laurel Hill Creek Basin. The mean streamflow per unit area for tributaries draining areas that extend into the western ridge and draining eastern or central sections was 0.24 and 0.05 ft 3 /s/mi 2 , respectively. In general, as the drainage area increased for tributary basins, the streamflow per unit area increased. Criteria established by the Pennsylvania Department of Environmental Protection indicate that the safe yield of water withdrawals from the Laurel Hill Creek Basin is 1.43 million gallons per day (Mgal/d). Water-use data for 2009 indicate that net (water withdrawals subtracted by water discharges) water withdrawals from groundwater and surface-water sources in the basin were approximately 1.93 Mgal/d. Water withdrawals were concentrated in the upper part of the basin with approximately 80 percent of the withdrawals occurring in the upper 36 mi 2 of the basin. Three subbasins—Allen Creek, Kooser Run, and Shafer Run— in the upper part were affected the most by water withdrawals such that safe yields were exceeded by more than 1,000 percent in the first two and more than 500 percent in the other. In the subbasin of Shafer Run, intermittent streamflow characterizes sections that historically have been perennial. The GSFLOW model of the Laurel Hill Creek Basin is a simple one-layer representation of the groundwater flow system. The GSFLOW model was primarily calibrated to reduce the error term associated with base-flow periods. The total amount of observed streamflow at the Laurel Hill Creek at Ursina, Pa. streamflow-gaging station and the simulated streamflow were within 0.1 percent over the entire modeled period; however, annual differences between simulated and observed streamflow showed a range of -27 to 24 percent from 1992 to 2007 with nine of the years having less than a 10-percent difference. The primary source of simulated streamflow in the GSFLOW model was the subsurface (interflow; 62 percent), followed by groundwater (25 percent) and surface runoff (13 percent). Most of the simulated subsurface flow that reached the stream was in the form of slow flow as opposed to preferential (fast) interflow.

Pennsylvania↗

Evaluation of the Bushy Park Reservoir three-dimensional hydrodynamic and water-quality model, South Carolina, 2012–15

The Bushy Park Reservoir is a relatively shallow impoundment in southeastern South Carolina. The reservoir, located under a semi-tropical climate, is the principal water supply for the city of Charleston, South Carolina, and the surrounding areas including the Bushy Park Industrial Complex. Although there was an adequate supply of freshwater in the reservoir in 2022, water-quality concerns are present over taste-and-odor and saltwater-intrusion issues. From 2013 to 2015, the U.S. Geological Survey (USGS), in cooperation with the Charleston Water System, engaged in a multi-year study of the hydrology and hydrodynamics of Bushy Park Reservoir to better understand factors affecting water-quality conditions in the reservoir. As part of this study, Charleston Water System worked with Tetra Tech, Inc., a consulting and engineering firm, to develop a Bushy Park Reservoir hydrodynamic and water-quality modeling framework, built upon earlier efforts by both Tetra Tech and the U.S. Army Corps of Engineers. At the completion of the new modeling framework, the USGS was requested to evaluate the calibrated hydrodynamic and water-quality model. The Bushy Park Reservoir Environmental Fluid Dynamics Code (EFDC) model was calibrated for the time period from January 1, 2012, to December 31, 2015. The general modeling approach for the newly revised modeling framework, as briefly detailed in this report, was developed with EFDC. The EFDC is a grid-based modeling package that can simulate three-dimensional flow, transport, and water quality in surface-water systems. This report evaluated the capacity of Tetra Tech’s Bushy Park EFDC model to simulate water discharge, water circulation, surface elevations, temperature, salinity, and other water-quality parameters. The USGS model review focused specifically on the following criteria: (1) determine if the model, with additional effort, could be developed into an adequate planning tool for Bushy Park Reservoir; (2) assess the capacity of the model to specifically address water-quality issues in the reservoir related to taste-and-odor and saltwater intrusions; and, (3) evaluate three preliminary water-management scenarios related to reduced water withdrawals in the reservoir and the effect on saltwater intrusion. Overall, the model was able to simulate discharge, flow velocity, and water-surface elevations with generally good agreement between the simulated and measured values. Specifically, the model was able to demonstrate good agreement for discharge at two USGS continuous discharge locations (USGS station 02172002; USGS station 02172040), with Wilmott index of agreements of 0.86 and 0.75, respectively. A total of seven USGS streamgages, located on the West Branch of the Cooper River, Durham Canal, and the Cooper River, were available for water-surface elevations, with index of agreements ranging from 0.74 to 0.99. However, model-simulated water-surface elevation ranges were appreciably high (compared to measured ranges) for two locations near Pinopolis Dam, farthest upstream on the West Branch of the Cooper River. This result may indicate that too much simulated tidal energy propagated through the model domain. For water temperature, 16 calibration stations were available for at least part of the 4-year simulation. The index of agreement range for temperature comparisons was from 0.95 to 1.00, indicating excellent agreement between the measured and simulated results. One of the primary future applications for the Bushy Park Reservoir EFDC model is to determine the extent of saltwater intrusions. A wide range in the salinity prediction quality was simulated with the model. The prediction quality ranged from an index of agreement of 0.15 at Cooper River approximately 2.75 miles southeast of the Tee, South Carolina, to 0.92 at West Branch Cooper River near Moncks Corner, South Carolina. Although the model did not accurately simulate some of the larger salinity deviations resulting during individual hydrologic events, the seasonal salinity trends were adequately simulated with the model during the study period (2012–15). Therefore, it may be difficult to simulate extreme hydrologic events, such as during large storms, where high salinity water is exchanged with Bushy Park Reservoir. There was agreement in model simulation with the measured data either on the quantitative index of agreement values or qualitative agreement in the seasonal salinity data trends. For water quality, the index of agreement values were generally low for total nitrogen, ammonia, nitrate, total Kjeldahl nitrogen, total phosphorus, and orthophosphate. Although general trends were adequately simulated at specific stations, particularly for Bushy Park Reservoir, the model-simulated fit was low across all the constituents described above with index of agreements usually below 0.50. A limitation for simulating nutrient concentrations across the model domain was the lack of characterization for the constituents directly entering Bushy Park Reservoir, or the lack of data directly attributed to the boundary condition (for example, the Cooper River). The other two calibrated water-quality constituents (besides the nutrients mentioned above) were dissolved oxygen and chlorophyll a . Dissolved oxygen varied from index of agreement values from 0.58 to 0.94 for 11 stations, generally indicating agreement with the available measured data. Chlorophyll a , calibrated for seven stations, had a wider range from 0.11 to 0.74 for the index of agreement. With the current modeling framework, taste-and-odor events, related to cyanobacterial blooms, cannot be directly simulated. However, indirect estimates of cyanobacteria concentrations may be obtained by using the chlorophyll a model outputs, which represent total phytoplankton biomass, and the phytoplankton biovolume data by group (diatoms, green algae, cyanobacteria and others) collected from 2012 to 2015. For the Bushy Park Reservoir modeling framework to be used directly for taste-and-odor issues, cyanobacteria must be simulated and calibrated based on observations of cyanobacteria biomass concentrations. In addition to the cyanobacteria sampling conducted within the reservoir between 2012 and 2015, the new model calibration would also require new algae biomass data-collection efforts to characterize the external sources of cyanobacteria entering the Bushy Park Reservoir from tributaries, as well as the internal cycling, production, and decay of cyanobacteria in the hydrologic system. Further improvements to the EFDC model would include expanding the collection of boundary condition datasets, such as water-quality monitoring to determine improved nutrient loads into the model domain. Along with improved water-quality monitoring for the major boundary conditions, continuous discharge, for both Foster Creek and the Back River, would further constrain the flow balance and the loads into Bushy Park Reservoir. In addition to better boundary-condition characterization, it is important to better characterize possible shortcomings specifically to the model domain, such as the grid resolution, bathymetry, and numerical hydrodynamic errors. Further consideration of the model may involve a sensitivity analysis to determine if errors in the simulation outputs, such as discharge, water-surface elevations, and salinity, were more likely caused by poor boundary condition characterization or, specifically, the model setup. Three model scenarios were run with the revised Bushy Park Reservoir model: (1) reduced withdrawals from one of the large intake-discharge locations for Bushy Park Reservoir, the Williams Station; (2) elevated (above background levels) ocean water level causing saltwater intrusion from the ocean through Durham Canal into Bushy Park Reservoir; and (3) overtopping of the Back River Dam at the southernmost end of Bushy Park Reservoir. For the reduced withdrawals scenarios, the largest shift in flow resulted near the Williams Station intake, with the next largest flow change at the southern end of Bushy Park Reservoir, and a net increase in flow out of the Bushy Park Reservoir to the Cooper River by way of the Durham Canal. The effect resulting from scenario 3 on water quality and salinity was small, with larger increases for dissolved oxygen than other constituents at several monitoring stations. For the two scenarios related to saltwater intrusion (including dam overtopping), the changes in salinity generally were found to dissipate in the following 2 weeks and generally back to baseline salinity conditions within 3 months. This result did vary depending on the severity of the storm or length of the dam overtopping event.

South Carolina↗

An update of hydrologic conditions and distribution of selected constituents in water, eastern Snake River aquifer and perched groundwater zones, Idaho National Laboratory, Idaho, emphasis 2019–21

Since 1952, wastewater discharged to infiltration ponds (also called “percolation ponds”) and disposal wells at the Idaho National Laboratory (INL) has affected water quality in the eastern Snake River Plain (ESRP) aquifer and perched groundwater zones underlying the INL. The U.S. Geological Survey (USGS), in cooperation with the U.S. Department of Energy (DOE), maintains groundwater-monitoring networks at the INL to determine hydrologic trends and to delineate the movement of radiochemical and chemical wastes in both the aquifer and perched groundwater zones. This report presents an analysis of water-level and water-quality data collected from the ESRP aquifer and perched groundwater wells from the USGS groundwater monitoring networks during 2019–21. From March–May 2018 to March–May 2021, water levels in wells completed in the ESRP aquifer increased in the northern part of the INL and decreased in the southwestern part. Water-level increases ranged from 0.02 to 1.04 feet in the northern part and decreases ranged from 0.03 to 2.94 feet in the southwestern part of the INL. Detectable concentrations of radiochemical constituents in water samples from wells in the ESRP aquifer at the INL generally decreased or remained constant during 2019–21. Decreases in concentrations were attributed to radioactive decay, changes in waste-disposal methods, and dilution from recharge and underflow. In 2021, tritium was detected above reporting levels in water samples collected from 46 of 105 aquifer wells and ranged from 150±50 to 4,280±150 picocuries per liter (pCi/L). Tritium concentrations from eight wells completed in deep perched groundwater near the Advanced Test Reactor Complex (ATRC) generally were greater than or equal to the reporting level during at least one sampling event during 2019–21, and concentrations ranged from 160±50 to 2,097±107 pCi/L. Concentrations of strontium-90 in water from 12 of 45 aquifer wells sampled in 2021 exceeded the reporting level, and concentrations ranged from 2.5±0.7 to 299±6 pCi/L. During 2021, concentrations of strontium-90 from five wells completed in deep perched groundwater at the ATRC equaled or exceeded the reporting levels, and concentrations ranged from 3±0.9 pCi/L to 27.8±1.3 pCi/L. Concentrations of cesium-137 were less than the reporting level in all but one aquifer well, and concentrations of plutonium-238, plutonium-239, -240 (undivided), and americium-241 were less than the reporting level in water samples from all aquifer wells sampled during this study period. Dissolved chromium concentrations in water samples from 64 ESRP aquifer wells ranged from less than (<) 0.5 to 76.4 micrograms per liter (μg/L). During 2019–21, dissolved chromium was detected in water from wells completed in deep perched groundwater above the ESRP aquifer at the ATRC, and concentrations ranged from <1 to 82.1 μg/L. In 2021, concentrations of dissolved sodium in water from most ESRP aquifer wells in the southern part of the INL were greater than the western tributary groundwater background concentration of 8.3 milligrams per liter (mg/L). During 2021, dissolved sodium concentrations in water from 15 wells completed in deep perched groundwater ranged from 11.7 to 122.5 mg/L. Variations in sodium concentrations in aquifer wells and perched groundwater zones are attributed to either migration of remnant water from the former chemical-waste ponds or disposal volume and composition variability in percolation ponds installed in 2008. In 2021, concentrations of chloride in most water samples from ESRP aquifer wells south of the Idaho Nuclear Technology and Engineering Center (INTEC) and at the Central Facilities Area (CFA) exceeded background concentrations. Chloride concentrations in water from wells south of the INTEC have generally decreased because of discontinued chloride disposal to the legacy percolation ponds since 2002 when the discharge of wastewater was discontinued. During 2019–21, dissolved chloride concentrations in deep perched groundwater above the ESRP aquifer from 18 wells at the ATRC ranged from 8.15 to 231 mg/L. In 2021, sulfate concentrations in water samples from ESRP aquifer wells in the south-central part of the INL that exceeded the background concentration of sulfate, ranged from 21 to 141 mg/L. The greater-than-background concentrations in water from these wells are attributed to sulfate disposal at the ATRC infiltration ponds or the legacy INTEC percolation ponds. In 2021, sulfate concentrations in water samples from aquifer wells near the Radioactive Waste Management Complex (RWMC) were mostly greater than background concentrations. The maximum dissolved sulfate concentration in shallow perched groundwater near the ATRC was 575 mg/L in 2021. During 2021, dissolved sulfate concentrations in water from wells completed in deep perched groundwater near the cold waste ponds at the ATRC ranged from 22.3 to 519 mg/L. In 2021, concentrations of nitrate in water from most ESRP aquifer wells at and near the INTEC exceeded the western tributary groundwater background concentration of 0.655 mg/L. Concentrations of nitrate in aquifer wells southwest of INTEC and farther away from the influence of disposal areas and the Big Lost River, in intermittent source of surface water recharge to the aquifer, show a general decrease in nitrate concentration over time. Two aquifer wells south of INTEC show increasing trends that could result from wastewater beneath the INTEC tank farm being mobilized to the aquifer. During 2019–21, water samples from several ESRP aquifer wells were collected and analyzed for volatile organic compounds (VOCs). Twelve VOCs were detected, and 1–4 VOCs were detected in water samples from 10 wells. The most frequently detected VOCs include carbon tetrachloride (tetrachloromethane), trichloromethane, tetrachloroethene, 1,1,1-trichloroethane, and trichloroethene. In 2019–21, concentrations for all VOCs were less than their respective maximum contaminant levels (MCLs) for drinking water, except carbon tetrachloride in one well, trichloroethene in two wells, and vinyl chloride in one well. During 2019–21, variability and bias were evaluated from 34 replicate and 14 blank quality-assurance samples. Results from replicate analyses were investigated to evaluate sample variability. Constituents with acceptable reproducibility were major ions, trace elements, nutrients, and VOCs. All radiochemical constituents including gross alpha- and beta- radioactivity, strontium-90, cesium-137, and tritium, had acceptable reproducibility. Bias from sample contamination was evaluated from equipment, field, and source-solution blanks. Chloride and sulfate were detected slightly above their respective method detection limits in equipment and field blanks, but at concentrations well below the co-collected sample for that well. These chloride and sulfate detections in the field and equipment blanks were inconsequential because they weren’t detected above the analysis-specific variability for those constituents as determined by replicate sample result evaluation. None of the detections of nutrients and trace inorganic constituents were high enough to indicate environmental sample or analytical procedure bias.

Idaho↗

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↗

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 &ldquo;young&rdquo; 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↗

Evaluating management alternatives for Wyoming elk feedgrounds in consideration of chronic wasting disease

Executive Summary The authors used decision and modeling analyses to evaluate management alternatives for a decision on whether to permit Cervus canadensis (elk) feeding on two sites on Bridger-Teton National Forest, Dell Creek and Forest Park. Supplemental feeding of elk could increase the transmission of chronic wasting disease (CWD) locally and disease spread regionally, potentially impacting elk populations over time with wider implications for Odocoileus hemionus (mule deer) and Odocoileus virginianus (white-tailed deer) populations and hunting, tourism, and regional revenue. Supplemental feeding is thought to improve overwinter elk survival and reduce the commingling of elk with cattle during months when brucellosis transmission risk is highest. We worked with the U.S. Department of Agriculture Forest Service to identify their fundamental objectives and associated performance metrics related to this feedground decision. We then developed disease and habitat selection models to quantify the effect of four management alternatives on select performance metrics. The four alternatives were to continue to permit feeding, phaseout permits to feed in three years, permit feeding on an emergency basis, or stop permitting feeding. In this report, we present methods and summarized results on disease and habitat selection models and summaries of other performance metrics analyzed by BIO-WEST, Inc. and Cirrus Ecological Solutions as part of an Environmental Impact Statement. Data from Wyoming Game and Fish Department (WGFD) supported the assumption that supplemental elk feeding allows for larger elk populations in a region. We documented that herd units (HU) without feedgrounds had 23 percent lower densities of elk per area of winter range when compared against HUs with feedgrounds, after accounting for differences in sightability of elk during counts on and off feedgrounds. Thus, throughout our analyses, we assumed feedground closures would reduce elk carrying capacity resulting in an average decline of previously fed elk population segments by 23 percent (5th and 95th percentiles = [11 percent, 35 percent]) by year 20. Most of that decline occurred within the first few years after a feedground ceases to operate. We used a panel of CWD experts to help estimate CWD trans-mission in fed and unfed elk population segments. In aggregate, the expert panel estimated that median values of direct and indirect transmission of CWD are expected to be 1.9 and 4 times higher, respectively, in fed elk populations compared to unfed elk. We used these disease transmission estimates in combination with local elk demographic rates and carrying capacity estimates to project disease and population dynamics. In year 20, we predicted CWD prevalence would increase to 42 percent (5th and 95th percentiles = [29 percent, 55 percent]), and 13 percent (5th and 95th percentiles = [4 percent, 26 percent]) on average for fed and unfed elk population segments, respectively, given a starting prevalence of 1.6 percent. The prevalence estimates for the unfed elk population segments are in the range of previous observations of CWD in elk in the western United States. The average CWD prevalence from 2016 to 2018 in the unfed elk population of Wind Cave National Park in South Dakota was 18 percent overall but up to 30 percent in some regions (Sargeant and others, 2021). Meanwhile, CWD prevalence in the Iron Mountain and Laramie Peak elk herds in Wyoming from 2016 to 2018 was 14 percent and 7 percent, respectively, despite being present since at least 2002 (Wyoming Game and Fish Department, 2020b). From 2016 to 2020, elk that were fed at Dell Creek and Forest Park constituted on average 12–20 percent of the total elk on their respective HUs. As a result, the differences between management alternatives are modest when considering the closure of only one feedground on a HU. The no feeding alternative for Forest Park resulted in a CWD prevalence of 17 percent (SD = 7 percent) in the Afton HU compared to 20 percent (SD = 7 percent) with continued feeding by year 20. In the Upper Green River HU, no feeding on Dell Creek resulted in a CWD prevalence of 27 percent (SD = 6 percent) compared to 30 percent (SD = 5 percent) with continued feeding. In terms of disease-associated mortality, we predicted the closure of Forest Park and Dell Creek feedgrounds would reduce the total number of CWD mortalities by 9 percent in the Upper Green River HU and 26 percent in the Afton HU during the 20-year timespan. Our spatial analyses predicted that management alternative effects vary by HU as a function of private property and other wildlife winter ranges proximity relative to feedground location. The predicted number of elk abortions on private land, as a proxy for brucellosis risk to cattle, may increase by 8–21 percent in the absence of feeding at Dell Creek and Forest Park. Eight feedgrounds are located on Bridger-Teton National Forest, all of which have permits that have expired or will expire prior to 2028. In addition, WGFD could change their management of feedgrounds given new information; therefore, we also assessed the cumulative effects of continued feeding, phaseout, and no feeding management alternatives across five HUs south of Jackson, Wyoming (Afton HU, Fall Creek HU, Piney HU, Pinedale HU, and Upper Green River HU). These five HUs ranged from about 41 to 85 percent of the elk herd using feedgrounds, which corresponded to a CWD prevalence at year 20 of 23–34 percent if all feedgrounds in those five HUs remained open relative to 12 to 14 percent if all feedgrounds were closed. We predicted feedground closures may result in immediate reductions in population size relative to alternatives that continue feeding (for example, continued feeding and emergency feeding alternatives); however, over longer periods of time, CWD-associated mortality leads to larger population reductions. The no feeding alternative resulted in higher elk population sizes compared to the continued feeding alternative after about 10 years of implementation. Delayed action under a phaseout alternative resulted in increasing the CWD prevalence to 20 percent relative to 12 to 14 percent, on average, without feeding on HUs with a large population of fed elk such as the Upper Green River HU. Summarizing our cumulative results across all five of the analyzed HUs, we predicted continued feeding will lead to fewer elk by year 20 (mean = 8,300, standard deviation [SD] = 740) compared to no feeding at U.S. Department of Agri-culture Forest Service sites (10,700, SD = 890). The closure of all feedgrounds was projected to result in the largest elk populations at year 20 (12,500, SD = 980). No feeding at all sites also resulted in the largest cumulative harvest of 57,700 (SD = 2,600) compared to 51,100 (SD = 3,800) for continued feeding at all current feedground sites on the five HUs. Continued feeding also resulted in the lowest brucellosis costs to producers ($194,600, SD = $11,500) compared to no feeding on all feedgrounds ($243,000, SD = $13,700). Assuming moderate reductions in hunter interest because of increasing CWD prevalence in elk, we predicted that no feeding resulted in regional revenues generated by hunting activities of $190 million (SD = $10 million) compared to $173 million (SD = $10 million) for continued feeding over the 20-year timeframe. Recent CWD detections in mule deer and elk in Grand Teton National Park has elevated the importance of the cur-rent decision on whether, and how, to permit elk feeding on Dell Creek and Forest Park and the management of the other feedgrounds. Aggressive male harvest has slowed, but not stopped, the increasing prevalence of CWD in mule deer (Conner and others, 2021). It is unclear whether harvest management can be an effective tool to slow the spread of CWD in elk. There are also no effective treatments or vaccines for CWD, and it is unlikely that any will be developed that can be easily deployed in the near future. Thus, reducing artificial aggregations is one of the few management approaches suggested by the Western Association of Fish and Wildlife Agencies (Almberg and others, 2017). Future surveillance and monitoring can be designed to resolve uncertainties that can improve future decision-making. If feedgrounds close, research could quantify elk population reductions in the absence of feeding, the redistribution of fed elk to other places, or the consequences of elk movement on private property. If feedgrounds remain open, research could assess how rapidly CWD spreads in artificial aggregations of elk; however, surveillance programs would need to be designed with sufficient power to detect initial changes of CWD prevalence. Delaying action on feedground management was projected to be costly. Results of the phaseout alternative relative to the no feeding alternative suggested a 3-year delay was enough for substantial long-term changes in CWD prevalence. The long-term persistence of infectious CWD prions in the environment suggests that feedground management decisions may have long-lasting consequences. Our results indicated tradeoffs in the ability of a management agency to achieve all their objectives, and all management alternatives resulted in significant reductions in elk population size. This report contains the foundational elements for formal decision analysis methods, which can be implemented to help decision makers transparently evaluate the consequences of decision alternatives and identify the set of actions that best achieve agency and stakeholder priorities.

Wyoming↗

Spatial and temporal variation of stream chemistry associated with contrasting geology and land-use patterns in the Chesapeake Bay watershed—Summary of results from Smith Creek, Virginia; Upper Chester River, Maryland; Conewago Creek, Pennsylvania; and Difficult Run, Virginia, 2010–2013

Despite widespread and ongoing implementation of conservation practices throughout the Chesapeake Bay watershed, water quality continues to be degraded by excess sediment and nutrient inputs. While the Chesapeake Bay Program has developed and maintains a large-scale and long-term monitoring network to detect improvements in water quality throughout the watershed, fewer resources have been allocated for monitoring smaller watersheds, even though water-quality improvements that may result from the implementation of conservation practices are likely to be first detected at smaller watershed scales. In 2010, the U.S. Geological Survey partnered with the U.S. Environmental Protection Agency and the U.S. Department of Agriculture to initiate water-quality monitoring in four selected small watersheds that were targeted for increased implementation of conservation practices. Smith Creek watershed is an agricultural watershed in the Shenandoah Valley of Virginia that is dominated by cattle and poultry production, and the Upper Chester River watershed is an agricultural watershed on the Eastern Shore of Maryland that is dominated by row-cropping activities. The Conewago Creek watershed is an agricultural watershed in southeastern Pennsylvania that is characterized by mixed agricultural activities. The fourth watershed, Difficult Run, is a suburban watershed in northern Virginia that is dominated by medium density residential development. The objective of this study was to investigate spatial and temporal variations in water chemistry and suspended sediment in these four relatively small watersheds that represent a range of land-use patterns and underlying geology to (1) characterize current water-quality conditions in these watersheds, and (2) identify the dominant sources, sinks, and transport processes in each watershed. The general study design involved two components. The first included intensive routine water-quality monitoring at an existing streamgage within each study area (including continuous water-quality monitoring as well as discrete water-quality sampling) to develop a detailed understanding of the temporal and hydrologic variability in stream chemistry and sediment transport in each watershed. The second component involved extensive water-quality monitoring at various sites throughout each watershed to develop a detailed understanding of spatial patterns. Both components were used to improve understanding of sources and transport processes affecting stream chemistry, including nutrients and suspended sediments, and their implications for detecting long-term trends related to best management practices. This report summarizes the results of monitoring that was performed from April 2010 through September 2013. Individual Small Watershed Summaries Summaries for each of the four small watersheds are presented below. Each watershed has a more descriptive and detailed section in the report, but these summaries may be particularly useful for some watershed managers and stakeholders desiring slightly less technical detail. Smith Creek Smith Creek is a 105.39-mi 2 watershed within the Shenandoah Valley that drains to the North Fork Shenandoah River. The long-term Smith Creek base-flow index is 72.3 percent, indicating that on average, approximately 72 percent of Smith Creek flow was base flow, which suggests that Smith Creek streamflow is dominated by groundwater discharge rather than stormwater runoff. A series of cluster and principal components analyses demonstrated that the majority of the variability in Smith Creek water quality could be attributed to hydrologic and seasonal variability. Statistically significant positive correlations with flow were observed for turbidity, suspended sediments, total nitrogen, ammonium, orthophosphate, iron, total phosphorus, and the ratio of calcium to magnesium. Statistically significant inverse correlations with flow were observed for specific conductance, magnesium, δ 15 N of nitrate, pH, bicarbonate, calcium, and δ 18 O of nitrate. Of particular note, flow and nitrate were not statistically significantly correlated, likely because of the relatively complex concentration-discharge relationship observed in continuous and discrete datasets. Statistically significant seasonal patterns were observed for numerous water-quality constituents: water temperature, turbidity, orthophosphate, total phosphorus, suspended-sediment concentration, and silica were higher during the warm season, but pH, dissolved oxygen, and sulfate were higher during the cool season. Surrogate regression models were developed to compute sediment and nutrient loads in Smith Creek using the continuous water-quality monitors. The mean Smith Creek in-stream sediment load was approximately 6,900 tons per year, with nearly 90 percent of the sediment load over the 3-year study period contributed during the eight largest storm events during that period. The Smith Creek total phosphorus load was approximately 21,000 pounds of phosphorus per year, with the majority of the load contributed during stormflow periods, although a substantial phosphorus load still occurs during base-flow conditions. The Smith Creek total nitrogen load was approximately 400,000 pounds per year, with total nitrogen accumulation less dominated by stormflow contributions (as was the case for sediment and total phosphorus) and strongly affected by base-flow export of nitrogen from the basin. Extensive water-quality monitoring throughout the Smith Creek watershed revealed how the complex geology and hydrology interacted to result in variable water chemistry. During relatively dry and low base-flow periods, much of the discharge in Smith Creek was contributed by a single dominant spring—Lacey Spring. During wetter base-flow periods, the flows in Smith Creek were largely generated by a mixture of headwater springs and forested mountain tributaries with very different geochemical composition. The headwater springs generally issued from limestone bedrock and were characterized as having relatively high nitrate, specific conductance, calcium, and magnesium, as well as relatively low concentrations of phosphorus, ammonium, iron, and manganese. The undeveloped, high-gradient, forested mountain sites were generally characterized by low ionic strength waters with low nutrient concentrations. Nitrate isotope data from the limestone springs generally were consistent with manure-derived nitrogen sources (such as cattle and poultry), although the possibility of other mixed sources cannot be excluded. Nitrate isotope data from the undeveloped, high-gradient forested mountain sites were more consistent with nitrogen from undisturbed soils, atmospheric deposition, or nitrogen fixation. Regardless of the nitrogen source, oxygen isotope data indicate that the nitrate was largely a result of nitrification. Land-use data indicate that manure sources of nitrogen dominated watershed nitrogen inputs. Phosphorus sources were less well studied. The presence of a single point-source discharge near the town of New Market contributed the majority of the phosphorus to Smith Creek under base-flow conditions, but nonpoint sources of phosphorus dominated the loading to Smith Creek during stormflow periods. Implementation of conservation practices increased in the Smith Creek watershed during the study period, and even though a broad range of practice types was implemented, the most common practices included stream fencing (for cattle exclusion), the development of nutrient management plans, conservation crop rotation, and the planting of cover crops. While the implementation of these conservation practices is encouraging, results indicate small increases in nitrate concentrations at the streamgage over the last 29 years, concurrent with small decreases in nitrate fluxes. It will likely be years before the cumulative effect of these practices can be detected in the Smith Creek water quality, and the magnitude of the effect of these conservation practices detected in Smith Creek will depend largely on whether nutrient loading (of manure and commercial fertilizer) is reduced over time. Upper Chester River The Upper Chester River watershed includes the 36-square-mile (mi 2 ) watershed area around several nontidal tributaries that drain into the tidal Chester River. The streamgage is on Chesterville Branch, the largest nontidal tributary (approximately 6.12 mi 2 ) and is the site for continuous water-quality monitoring for this project. The base-flow index at Chesterville Branch is about 72 percent and indicates that, as in most of the Coastal Plain, groundwater is the greatest contributor to streamflow. As such, more than 90 percent of the nitrogen in the stream is in the form of nitrate from groundwater. Continuous and discrete data collected at Chesterville Branch show the effects of streamflow and season on water quality. Significantly positive correlations with flow were observed for ammonium, dissolved and total phosphorus, sediment, and turbidity as runoff carried these constituents from the land surface into Chesterville Branch. Other constituents that increased significantly with flow include potassium, sulfate, iron, and manganese, which are likely contributed from near-stream areas and ponds with high organic-matter content. Total nitrogen, pH, and specific conductance, along with chemical constituents associated with groundwater inputs including nitrate, calcium, ratio of calcium to magnesium, silica, bicarbonate, and sodium, were negatively correlated with flow because concentrations of these constituents were diluted by runoff. Seasonal differences in water chemistry, which are most likely related to increased biologic effects on the uptake and release of chemicals in the stream and near-stream areas, also were observed. Water temperature, orthophosphate, δ 15 N of nitrate, bicarbonate, sodium, and the ratio of sodium to chloride were higher during the warm season, and dissolved oxygen, total nitrogen, nitrate, magnesium, sulfate, and manganese were higher during the cool season. Surrogate-regression models developed by using continuous water-quality data showed that the annual sediment load for the 2013 water year was about 2,600 tons, with more than 90 percent of this sediment contributed during two storms. The total phosphorus load in 2013 was about 13,000 pounds with more than 90 percent contributed during the same two storms as sediment. The load of total nitrogen, 140,000 pounds, accumulated steadily throughout the 2013 water year as nitrate in groundwater continuously discharged into the stream. The same two large storms that contributed 90 percent of the suspended-sediment and total phosphorus load only contributed about 20 percent of the annual total nitrogen load. Extensive water-quality monitoring of stream base flow throughout the Upper Chester River watershed identified how differences in land use and hydrogeology affected water chemistry. In parts of the watershed with well-drained soil and thick sandy aquifer sediments, concentrations of nitrate and other chemicals associated with fertilizer and lime application increased in streams as agricultural land use increased. More than 90 percent of the nitrogen in streams from these areas was in the form of nitrate, and concentrations ranged from about 5 milligrams per liter (mg/L) to 8 mg/L as nitrogen in the two largest tributaries. Stream nitrate concentrations were about 1 mg/L as nitrogen where soils were more poorly drained, the surficial aquifer sediments were thinner, and forests and wetlands were more widespread than agriculture. Nitrate isotope data were consistent with inorganic fertilizers ± atmospheric deposition and N 2 fixation as sources of nitrogen, and with nitrification as the dominant nitrate-forming process. Nitrate reduction was indicated by elevated δ 15 N and δ 18 O values in some samples from streams draining watersheds with poorly drained soils. An analysis of land-use data and SPARROW modeling input data attributed almost 90 percent of the nitrogen sources in the Upper Chester River watershed to inorganic fertilizer and fixation of atmospheric nitrogen by legumes, which is in agreement with the isotopic characteristics of nitrate in this watershed. Local sources of manure are limited in this area. Total phosphorus concentrations during base flow ranged from below detection to about 0.2 mg/L. Stream phosphorus concentrations during base flow were generally lower than those measured during storms because most phosphorus transport likely occurs as phosphorus attached to sediment particles during runoff. Because manure is not widely used in this area, the major source of phosphorus is likely fertilizer. The implementation of conservation practices in the Upper Chester River watershed increased substantially during the study period, with a total implementation of 1,194 U.S. Department of Agriculture-compliant practices. The most frequently used practices were oriented towards nutrient and sediment control, including cover crops, nutrient management planning, conservation crop rotation, conservation tillage, and irrigation management. The current Chesapeake Bay model for this area predicts that implementation of best management practices should result in a 13-percent decrease in overall delivery of nitrogen to the Upper Chester River. Because most nitrogen travels through the groundwater system for years to decades before being discharged to streams, the time period of monitoring was not sufficient to see the effects of these practices on water quality. The magnitude of the effect that may eventually be detected will depend on the degree to which nitrate leaching into the groundwater system is reduced over time. Loadings of phosphorus and sediment are primarily transported during large runoff events and are difficult to control and analyze for trends because of their timing and episodic nature. Conewago Creek Conewago Creek has two primary monitoring locations—one near the middle of the 47-mi 2 watershed and the other near the outlet just upstream of the Susquehanna River. The base-flow index was 47.3 percent for 2012–2013, indicating that on average, approximately 53 percent of the streamflow in Conewago Creek exited the watershed as surface flow, which suggests that the stormwater runoff was somewhat greater than groundwater discharge (base flow). A series of cluster and principal components analyses demonstrated that the majority of the variability in the Conewago Creek water quality could be attributed to hydrologic and seasonal variability. Statistically significant positive correlations with flow were observed at both monitoring sites for ammonium, total phosphorus, orthophosphate, iron, and manganese; additionally, at the upstream monitoring station, total nitrogen demonstrated a statistically significant positive correlation with flow. Statistically significant inverse correlations with flow were observed at both sites for water temperature, specific conductance (at the downstream site only), sulfate, chloride, calcium, and magnesium. Statistically significant seasonal patterns were observed for several water-quality constituents. Water temperature, phosphorus (upstream site only), and orthophosphate were higher during the warm season, and nitrate and total nitrogen (upstream site only) were higher during the cool season. Surrogate regression models were developed to compute sediment and nutrient load in Conewago Creek by using the continuous water-quality monitors and water-quality samples. Conewago Creek sediment load was approximately 9,900 tons in 2012 and approximately 18,900 tons in 2013, with nearly 80 percent of the sediment load in 2013 contributed by the three largest storm events. Annual total nitrogen loads could not be estimated due to poor model performance. The addition of continued monitoring or a continuously recording nitrate sensor could improve estimates of total nitrogen loads. During 2012 and 2013, phosphorus loads in Conewago Creek were approximately 50,000 pounds in each year. Combining data from one high-flow synoptic sampling with the data from routine sampling revealed how the geology and hydrology interact to result in variable water chemistry throughout the Conewago Creek watershed. The areas above the upstream gage in the headwaters are generally underlain by forested non-carbonate bedrock and are characterized by relatively low nitrate, specific conductance, calcium, and magnesium, as well as relatively low concentrations of phosphorus, ammonium, iron, and manganese. The more developed, agricultural areas below the upstream site were generally characterized by higher ionic strength waters with higher nutrient and metal concentrations. An analysis of land-use data and SPAtially Referenced Regressions On Watershed (SPARROW) modeling data indicates that manure sources of nitrogen dominate the input of nitrogen to the watershed. Implementation of conservation practices increased in the Conewago Creek watershed during the study period, and while a broad range of practice types were implemented, the most common practices included residue and tillage management, cover crops, nutrient management, terracing, and stream fencing (for animal exclusion or bank restoration). While the implementation of these conservation practices is encouraging, the cumulative effects of these practices probably will not be detected in Conewago Creek water quality for several years. The magnitude of the effects of these conservation practices on water quality in Conewago Creek will depend largely on the extent to which nutrient loading (septic, manure, and commercial fertilizer) and sediment-producing activities are reduced over time. Difficult Run The Difficult Run watershed is a 57.82-mi 2 watershed that drains to the Potomac River. The long-term Difficult Run base-flow index (from 1936 to 2010) was 57.9, indicating that approximately 58 percent of streamflow exited the watershed as base flow and 42 percent as stormflow; however, with continued development and urbanization of the watershed, the base-flow index has decreased to 50 percent during the last 20 years. This base-flow index was less than those of the other watersheds evaluated in this study, likely because the Difficult Run watershed largely is underlain by crystalline piedmont metamorphic rocks and has a greater proportion of impervious urban land cover. A series of cluster and principal components analyses indicated that most of the variability in Difficult Run water quality could be attributed to hydrologic variability and seasonality. Statistically significant positive correlations with flow were observed for turbidity, dissolved oxygen, suspended sediments, ammonium, orthophosphate, iron, and total phosphorus. Statistically significant inverse correlations with flow were observed for water temperature, pH, specific conductance, bicarbonate, calcium, magnesium, nitrate, δ 15 N of nitrate, and silica. Statistically significant seasonal patterns were observed for numerous water-quality constituents: water temperature, ammonium, orthophosphate, and δ 15 N of nitrate were higher during the warm season, and dissolved oxygen, nitrate, and manganese were higher during the cool season. Surrogate regression models were developed to compute sediment and nutrient loading rates. The Difficult Run sediment load was approximately 8,000 tons per year, with greater than 95 percent of the sediment load in the 2013 water year contributed by the seven largest storm events. The total phosphorus load in Difficult Run was approximately 14,000 pounds of phosphorus per year, with the majority of the load contributed during stormflow periods. The total nitrogen load in Difficult Run is estimated to have been approximately 140,000 pounds per year, with total nitrogen accumulation less dominated by stormflow contributions than that of phosphorus and strongly affected by base-flow export of nitrogen from the basin. Extensive water-quality monitoring throughout the Difficult Run watershed revealed relatively uniform generation of flow per unit of watershed area, as well as spatial variation in water quality that is strongly related to land-use activities. Elevated nitrate concentrations were observed in a subset of monitoring sites that are inversely correlated with population density and positively correlated to the septic system density within each subwatershed. The majority of the elevated nitrate concentrations for these sites are hypothesized to be caused by nitrate leaching from septic systems, more so than homeowner fertilizer usage among these subwatersheds that have lower population densities than other parts of the watershed. Nitrate isotope data, temporal patterns in the water-quality data, mass-balance computations, and a separate land-use analysis all generally indicate that leachate from septic systems was the likely source of the elevated nitrate. Another group of water-quality sites have relatively low nitrogen concentrations, are located in areas that are served by city sewer lines, and have experienced stream restoration activities. A final group of sites drained the areas with the highest imperviousness and had strongly elevated specific conductance, chloride, and sodium, which were likely caused by a combination of road salting and other anthropogenic sources draining these urbanized areas in the watershed. A fourth group of sites represents a mixture of water sources and had water quality similar to that at the Difficult Run streamgage. Analysis of the nitrate isotope data generally indicates a broad range of composition indicative of mixed natural and anthropogenic nitrogen sources. Implementation of conservation practices increased in the Difficult Run watershed during the study period, and while a broad range of practice types was implemented, the most common practices included stream restoration. While the implementation of these conservation practices is encouraging, the cumulative effect of these practices probably will not be detected in Difficult Run water quality for several years.

Maryland, Pennsylvania, Virginia↗