SEARCH · USGS Science

Results for “Environmental Modeling & Assessment”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,549 records · Page 86Linked to original sources

USGS2001

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↗

USGS2023

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↗

USGS2010

Effects of low-impact-development (LID) practices on streamflow, runoff quantity, and runoff quality in the Ipswich River Basin, Massachusetts: A summary of field and modeling studies

Low-impact-development (LID) approaches are intended to create, retain, or restore natural hydrologic and water-quality conditions that may be affected by human alterations. Wide-scale implementation of LID techniques may offer the possibility of improving conditions in river basins, such as the Ipswich River Basin in Massachusetts, that have run dry during the summer because of groundwater withdrawals and drought. From 2005 to 2008, the U.S. Geological Survey, in a cooperative funding agreement with the Massachusetts Department of Conservation and Recreation, monitored small-scale installations of LID enhancements designed to diminish the effects of storm runoff on the quantity and quality of surface water and groundwater. Funding for the studies also was contributed by the U.S. Environmental Protection Agency’s Targeted Watersheds Grant Program through a financial assistance agreement with Massachusetts Department of Conservation and Recreation. The monitoring studies examined the effects of replacing an impervious parking-lot surface with a porous surface on groundwater quality, installing rain gardens and porous pavement in a neighborhood of 3 acres on the quantity and quality of stormwater runoff, and installing a 3,000-ft2 (square-foot) green roof on the quantity and quality of rainfall-generated roof runoff. In addition to these small-scale installations, the U.S. Geological Survey’s Ipswich River Basin model was used to simulate the basin-wide effects on streamflow of several changes: broad-scale implementation of LID techniques, reduced water-supply withdrawals, and water-conservation measures. Water-supply and conservation scenarios for application in model simulations were developed with the assistance of two technical advisory committees that included representatives of State agencies responsible for water resources, the U.S. Environmental Protection Agency, the U.S. Geological Survey, water suppliers, and non-governmental organizations. From June 2005 to June 2007, groundwater quality was monitored at the Silver Lake town beach parking lot in Wilmington, Massachusetts, prior to and following the replacement of the conventional, impervious-asphalt surface with a porous surface consisting primarily of porous asphalt and porous pavers designed to enhance rainfall infiltration into the groundwater and to minimize runoff to Silver Lake. Concentrations of phosphorus, nitrogen, cadmium, chromium, copper, lead, nickel, zinc, and total petroleum hydrocarbons in groundwater were monitored. Enhancing infiltration of precipitation did not result in discernible increases in concentrations of these potential groundwater contaminants. Concentrations of dissolved oxygen increased slightly in groundwater profiles following the removal of the impervious asphalt parking-lot surface. In Wilmington, Massachusetts, in a 3-acre neighborhood, stormwater runoff volume and quality were monitored to determine the ability of selected LID enhancements (rain gardens and porous paving stones) to reduce flows and loads of the selected constituents to Silver Lake. Water-quality samples were analyzed for nutrients, metals, total petroleum hydrocarbons, and total-coliform and E. coli bacteria. A decrease in runoff quantity was observed for storms of 0.25 inch or less of precipitation. Water-quality-monitoring results were inconclusive; there were no statistically significant differences in concentrations or loads when the pre- and post-installation-period samples were compared. In a third field study, the characteristics of runoff from a vegetated "green" roof and a conventional, rubber-membrane roof were compared. The two primary factors affecting the green roof’s water-storage capacity were the amount of precipitation and antecedent dry period. Although concentrations of many of the chemicals in roof runoff were higher from the green roof than from the conventional roof, the ability of the green roof to retain water generally resulted in decreased differences between the total amounts (loads) of the chemicals that ran off the roofs. Land-use and water-management changes associated with LID implementation were investigated at multiple spatial scales, using the U.S. Geological Survey’s Ipswich River Basin model, to evaluate the effects of updated water-supply withdrawals for the towns of Reading and Wilmington (representing new baseline conditions for all simulations), potential land-use changes at buildout (potential future development), widespread implementation of retrofitting LID techniques, basin-scale water withdrawal reductions based on water-conservation pilot programs conducted by the Massachusetts Department of Conservation and Recreation, and land-use change and LID applications at a local scale. The new baseline simulation indicated that reduced water-supply withdrawals for the towns of Reading and Wilmington led to substantially higher medium and low flows in most of the reaches upstream from the South Middleton streamgage in the upper Ipswich River basin. Overall, simulations pointed to the importance of spatial scale in determining the effects of land-use change and LID practices on streamflow. Potential land-use changes at buildout had modest effects on streamflow in most subbasins (percent differences of less than 20 percent) because relatively little land in the basin was available for development. Results of the simulations conducted to evaluate widespread effective-impervious-area reductions upstream from the South Middleton streamgage indicated that the percentages of urban land use and associated effective impervious area were too small for even a 50-percent reduction of effective impervious area to appreciably affect streamflow in most subbasins. In contrast, the results of the hypothetical local-scale simulations indicated that for smaller streams, with high percentages of urban land use and associated effective impervious area, land-use change, development patterns, and LID practices may have substantial effects on streamflow. Modeling studies concurred with the results of fieldwork in the assessment that LID enhancements would likely have the greatest effect on decreasing stormwater runoff when broadly applied to highly impervious urban areas.

Massachusetts↗

USGS2014

Evaluation of toxicity to the amphipod, Hyalella azteca , and to the midge, Chironomus dilutus ; and bioaccumulation by the oligochaete, Lumbriculus variegatus , with exposure to PCB-contaminated sediments from Anniston, Alabama

The U.S. Environmental Protection Agency (USEPA) requested that as part of the remedial investigation for the Anniston, Alabama Polychlorinated Biphenyl (PCB) Site (Anniston PCB Site), that Pharmacia Corporation and Solutia Inc. (P/S) perform long-term reproduction toxicity tests with the amphipod, Hyalella azteca , and the midge, Chironomus dilutus , and bioaccumulation tests with the oligochaete, Lumbriculus variegatus , using sediment samples collected from reference locations and from Operable Unit 4 of the Anniston PCB Site. The sediment toxicity testing and sediment bioaccumulation results will be used by ARCADIS U.S., Inc. (ARCADIS) as part of a weight-of-evidence assessment to evaluate risks and establish sediment remediation goals for contaminants to sediment-dwelling organisms inhabiting the Anniston PCB Site. The goal of this study was to characterize relations between sediment chemistry and sediment toxicity and relations between sediment chemistry and sediment bioaccumulation in samples of sediments collected from the Anniston PCB Site. A total of 32 samples were evaluated from six test sites and one reference site to provide a wide range in concentrations of chemicals of potential concern (COPCs) including PCBs in samples of whole sediment. The goal of this study was not to determine the extent of sediment contamination across the Anniston PCB Site. Hence, the test sites or samples collected from within a test site were not selected to represent the spatial extent of sediment contamination across the Anniston PCB Site. Sediment chemistry, pore-water chemistry, and sediment toxicity data were generated for 26 sediment samples from the Anniston PCB Site. All of the samples were evaluated to determine if they qualified as reference sediment samples. Those samples that met the chemical selection criteria and biological selection criteria were identified as reference samples and used to develop the reference envelope for each toxicity test endpoint. Physical characterization of samples of whole sediment included analyses of grain size, TOC, and nutrients. Organic chemical characterization of samples of whole sediment included PCB homologs and select (13) PCB congeners, parent and alkylated polycyclic aromatic hydrocarbons (PAHs), organochlorine pesticides, and polychlorinated dibenzo-p-dioxins; and dibenzofurans. The PCB aroclors analyzed included 1016, 1221, 1232, 1242, 1248, 1254, 1260, 1262 and 1268. Analyses of whole sediment also included total metals, simultaneously extracted metals, and acid volatile sulfide. Chemical characterization of samples of pore water isolated from samples of whole sediment at the start of the sediment toxicity exposures or at the start of the sediment bioaccumulation exposures included metals, major cations, major anions, dissolved organic carbon, and additional water-quality characteristics. Concentrations of metals or PCBs in pore water during the sediment toxicity exposures or during sediment bioaccumulation exposures also were determined using peeper samples (for metals) or solid-phase microextraction (SPME) samplers (for PCBs). The bioavailability and bioaccumulation of PCBs in 14 sediment samples were investigated using SPME passive samplers and the 28-d L. variegatus whole-sediment bioaccumulation exposures In general the accumulation of PCBs consistently was predicted through the use of organic carbon normalization and equilibrium partitioning. In these sediments, PCB homologs were accumulated differently based on bioavailability and potential to accumulate in oligochaetes. As part of this assessment homolog specific biota sediment accumulation factor values were developed that could be applied across the larger site to predict tissue levels of PCBs. The whole-sediment toxicity tests done with H. azteca and C. dilutus met the established ASTM and USEPA test acceptability criteria. The most responsive H. azteca endpoints were day 42 survival normalized young per female and day 28 biomass and that the most responsive C. dilutus endpoints were adult biomass and percent adult emergence. Overall, between the two species, the most responsive endpoint assessed for these two species was H. azteca survival-normalized young per female (67 percent of the samples classified as toxic). Concentration-response models (CRMs) and site-specific sediment toxicity thresholds (TTs) were generated with matching sediment chemistry and sediment toxicity data. Sediment chemistry, pore-water chemistry, and sediment toxicity data were evaluated for as many as 26 sediment samples from the Anniston PCB Site. The reference-envelope approach was used to identify the sediment samples that were toxic to benthic invertebrates. This procedure involved identification of reference sediment samples, normalizing the toxicity data to reflect control responses, developing a reference envelope for each toxicity test endpoint, and designating each sediment sample as toxic or not toxic for each toxicity test endpoint, for each species, and for all species combined. These results demonstrated percent emergence of adult C. dilutus , biomass of adult C. dilutus , and reproduction of H. azteca normalized to percent survival were among the most responsive endpoints that were evaluated. Therefore, these endpoints were selected for CRM development. The site-specific TTs for whole sediment provide a reliable basis for identifying toxic and not toxic sediment samples in the Anniston PCB Site (that is, for correctly classifying the sediment samples used to derive the TTs as toxic or not toxic, for the endpoint used to derive the TTs). Among the 69 TTs for sediment, the TT LRs for total PCB homologs [499 to 1,870 micrograms per kilogram dry weight (μg/kg DW)] and for lead [(9.48 to 10.3 milligrams per kilogram (mg/kg) DW] based on reproduction of H. azteca or based on emergence or biomass of adult C. dilutus , were the most reliable. Such TTs had low rates of false negative errors (that is, only 0 to 11 percent of the samples below the TT were toxic to benthic invertebrates), low rates of false positive errors (only 0 to 6 percent of the samples greater than the TT were not toxic to benthic invertebrates), and high rates of correct classification (that is, 92 to 96 percent). The site-specific TTs for PCBs and other COPCs derived in this study also were compared to empirically based sediment quality guidelines (SQGs), to equilibrium-partitioning based SQGs, and to the results of spiked-sediment toxicity tests. The results of this evaluation indicated that the site-specific sediment TTs for PCBs were comparable to the consensus-based SQGs that were derived for PCBs. In addition, the site-specific sediment TTs for PCBs are well within the range of SQGs derived using the equilibrium partitioning approach. The site-specific sediment TTs for PCBs also are consistent with the results of chronic TTs that have been estimated for benthic invertebrates using the results of spiked-sediment toxicity tests. As the site-specific sediment TTs for PCBs are consistent with empirically based SQGs, equilibrium-partitioning based SQGs, and results of sediment-spiking studies, these site- specific sediment TTs likely represent the concentrations of PCBs that are sufficient to cause toxicity to benthic invertebrates (as opposed to simply being correlated with adverse effects on the survival, weight, or reproduction of benthic invertebrates). Importantly, such site-specific sediment TTs have been demonstrated to accurately classify sediment samples as toxic or not toxic to benthic invertebrates at the Anniston PCB Site. In contrast, the TTs for metals, PAHs, and organochlorine pesticides were generally lower than consensus-based SQGs (that is, probable effect concentrations), and LC 50s (median lethal effect concentrations) generated in spiked-sediment toxicity tests, indicating that these COPCs are likely not the main contributors to the observed toxicity of the site sediments evaluated in this study. The reproduction endpoint for H. azteca provided lower TTs compared to the day 28 biomass endpoint for H. azteca and the emergence or biomass endpoints for adult C. dilutus provided lower TTs compared to the day 13 biomass endpoint for C. dilutus .

Alabama↗

USGS2015

Chloride concentrations, loads, and yields in four watersheds along Interstate 95, southeastern Connecticut, 2008-11: factors that affect peak chloride concentrations during winter storms

Chloride (Cl - ) concentrations and loads and other water chemistry characteristics were assessed to evaluate potential effects of road-deicer applications on streamwater quality in four watersheds along Interstate 95 (I–95) in southeastern Connecticut from November 1, 2008, through September 30, 2011. Streamflow and water quality were studied in the Four Mile River, Oil Mill Brook, Stony Brook, and Jordan Brook watersheds, where developed land ranged from 9 to 32 percent. Water-quality samples were collected and specific conductance was measured continuously at paired water-quality monitoring sites, upstream and downstream from I–95. Specific conductance values were related to Cl - concentrations to assist in determining the effects of road-deicing operations on the levels of Cl - in the streams. Streamflow and water-quality data were compared with weather data and with the timing, amount, and composition of deicers applied to State highways. Grab samples were collected during winter stormwater-runoff events, such as winter storms or periods of rain or warm temperatures in which melting takes place. Grab samples were also collected periodically during the spring and summer and during base-flow conditions. The estimated Cl - concentrations at the eight water-quality monitoring sites during winter storms peaked as high as 270 milligrams per liter (mg/L) and were well below the U.S. Environmental Protection Agency (EPA) recommended acute chloride toxicity criterion of 860 mg/L and the chronic 4-day average toxicity criterion of 230 mg/L. Cl - concentrations in streams, particularly at sites downstream from I–95, peaked during increased streamflow in the winter and early spring as a result of deicers applied to roads and washed off by stormwater or meltwater. Cl - concentrations during most of the nonwinter seasons decreased during increases in streamflow because storm runoff was more dilute than base flow. However, peaks in specific conductance and estimated chloride concentrations at streams in more urbanized watersheds corresponded to peaks in streamflow well after winter snow or ice events; these delayed peaks in Cl - concentration likely resulted from deicer residue that remained in melting snow piles and on roadsides and (or) that were flushed from soils and shallow groundwater, then discharged downstream. Estimated peak Cl - concentrations varied with the type of winter storm event and were highest during or after winter storms of frozen precipitation and rain, in which the rain or meltwater effectively washed off the deicers. Estimated peak Cl - concentrations correlated positively with the duration of deicer application but generally not with streamflow. Estimated peak Cl - concentrations during the winter season were highest during low streamflow at most sites. Chloride concentrations varied considerably in shallow groundwater as a result of land-use differences. Cl - concentrations were very high (as high as 800 mg/L) in shallow groundwater downstream from I–95 at the Four Mile River site. Chloride/bromide mass concentration ratios and the proximity of a former landfill and sewage lagoon upstream indicate a likely source of Cl - is landfill leachate and possibly sewage leachate. Cl - loads in streams generally were highest in the winter and early spring. The estimated daily Cl - yield for the four monitoring sites downstream from I–95 ranged from 0.0004 ton per day per square mile for one of the least developed watersheds to 0.052 ton per day per square mile for the watershed with the highest percentage of urban development and impervious area. The estimated median contribution of Cl - load from atmospheric deposition was small and ranged from 0.07 percent of Cl - load at the Jordan Brook watershed to 0.57 percent at the Oil Mill Brook watershed. The Cl - loads in streams (outputs) were compared with Cl - load inputs, which include atmospheric deposition and deicer applications; Cl - load inputs were slightly larger than the Cl - load outputs at most of the sites during most years but do not account for the Cl - load in groundwater leaving the watersheds. A multiple linear regression model was developed to describe the variability of the natural log of peak specific conductance, as well as estimated Cl - concentrations. Five significant variables best explained the variability in the natural log of the peak specific conductance after deicing events: (1) snow on ground before deicing event; (2) winter precipitation with rain; (3) specific conductance in base flow; (4) State-operated road lane miles divided by watershed area; and (5) amount of Cl - from deicers applied to State-operated roads per lane mile. In this report, winter precipitation is defined as any type of precipitation, including frozen precipitation and rain, that occurs during the active deicing season, typically November through March. Frozen precipitation is defined here as snow, sleet, freezing rain, or any winter precipitation except rain. The addition of a lane mile in both directions on I–95 would result in an estimate of approximately 2 to 11 percent increase in Cl - input from deicers applied to I–95 and other roads maintained by Connecticut Department of Transportation. The largest estimated increase in Cl - load was in the watersheds with the greatest number miles of I–95 corridor relative to the total lane miles maintained by Connecticut Department of Transportation. On the basis of these estimates and the estimated peak Cl - concentrations during the study period, it is unlikely that the increased use of deicers on the additional lanes would lead to Cl - concentrations that exceed the aquatic habitat criteria.

Connecticut↗

USGS2016

Hydrologic analyses in support of the Navajo Generating Station–Kayenta Mine Complex environmental impact statement

Introduction The U.S. Department of Interior’s Bureau of Reclamation, Lower Colorado Region (Reclamation) is preparing an environmental impact statement (EIS) for the Navajo Generating Station-Kayenta Mine Complex Project (NGS-KMC Project). The proposed project involves various Federal approvals that would facilitate continued operation of the Navajo Generating Station (NGS) from December 23, 2019 through 2044, and continued operation of the Kayenta Mine and support facilities (collectively called the Kayenta Mine Complex, or KMC) to supply coal to the NGS for this operational period. The EIS will consider several project alternatives that are likely to produce different effects on the Navajo (N) aquifer; the N aquifer is the principal water resource in the Black Mesa area used by the Navajo Nation, Hopi Tribe, and Peabody Western Coal Company (PWCC). The N aquifer is composed of three hydraulically connected formations—the Navajo Sandstone, the Kayenta Formation, and the Lukachukai Member of the Wingate Sandstone—that function as a single aquifer. The N aquifer is confined under most of Black Mesa, and the overlying stratigraphy limits recharge to this part of the aquifer. The N aquifer is unconfined in areas surrounding Black Mesa, and most recharge occurs where the Navajo Sandstone is exposed in the area near Shonto, Arizona. Overlying the N aquifer is the D aquifer, which includes the Dakota Sandstone, Morrison Formation, Entrada Sandstone, and Carmel Formation. The aquifer is named for the Dakota Sandstone, which is the primary water-bearing unit. The NGS is located near Page, Arizona on the Navajo Nation. The KMC, which delivers coal to NGS by way of a dedicated electric railroad, is located approximately 83 miles southeast of NGS (about 125 miles northeast of Flagstaff, Arizona). The Kayenta Mine permit area is located on about 44,073 acres of land leased within the boundaries of the Hopi and Navajo Indian Reservations. KMC has been conducting mining and reclamation operations within the Kayenta Mine permit boundary since 1973. The KMC part of the proposed project requires approval by the Office of Surface Mining (OSM) of a significant revision of the mine’s permit to operate in accordance with the Surface Mine Control and Reclamation Act (Public Law 95-87, 91 Stat. 445 [30 U.S.C. 1201 et seq. ]) . The revision will identify coal resource areas that may be used to continue extracting coal at the present rate of approximately 8.2 million tons per year. The Kayenta Mine Complex uses water pumped from the D and N aquifers beneath PWCC’s leasehold to support mining and reclamation activities. Prior to 2006, water from the PWCC well field also was used to transport coal by way of a coal-slurry pipeline to the now-closed Mohave Generating Station. Water usage at the leasehold was approximately 4,100 acre-feet per year (acre-ft/yr) during the period the pipeline was in use, and declined to an average 1,255 acre-ft/yr from 2006 to 2011. The Probable Hydrologic Consequences (PHC) section of the mining and reclamation permit must be modified to project the consequences of extended water use by the mine for the duration of the KMC part of the project, including a post-mining reclamation period. Since 1971, the U.S. Geological Survey (USGS) has conducted the Black Mesa Monitoring Program, which consists of monitoring water levels and water quality in the N aquifer, compiling information on water use by PWCC and tribal communities, maintaining several stream-gaging stations, measuring discharge at selected springs, conducting special studies, and reporting findings. These data are useful in evaluating the effects on the N aquifer from PWCC and community pumping, and the effects of variable precipitation. The EIS will assess the impacts of continued pumping on the N aquifer, including changes in storage, water quality, and effects on spring and baseflow discharge, by proposed mining through 2044, and during the reclamation process to 2057. Several groundwater models exist for the area and Reclamation concluded it would conduct a peer review of the groundwater flow model that will be used to assess the direct, reasonably foreseeable indirect, and cumulative effects of future groundwater withdrawals on the D and N aquifers in the Black Mesa area. Reclamation made this determination because of the level of controversy around the effects of continued water use and the comments received from the 2014 draft EIS scoping meetings. Reclamation requested assistance from the USGS in evaluating existing groundwater flow models of the Black Mesa Basin that can be used to predict the effects of different project alternatives on the D and N aquifers.

Arizona↗

USGS2022

Use case development for earth monitoring, analysis, and prediction (EarthMAP)—A road map for future integrated predictive science at the U.S. Geological Survey

Executive Summary The U.S. Geological Survey (USGS) 21st-century science strategy 2020–30 promotes a bureau-wide strategy to develop and deliver an integrated, predictive science capability that works at the scales and timelines needed to inform societally relevant resource management and protection and public safety and environmental health decisions (U.S. Geological Survey, 2021). This is the overarching goal of the USGS Earth Monitoring, Analysis, and Prediction (EarthMAP) vision, which consists of three components: (1) integrated data and information, (2) integrated predictive science, and (3) actionable information—all designed and delivered to respond to user needs. To launch this vision and help shape the design and implementation of integrated predictive science, the USGS Regional Offices each developed a set of use cases (hereafter Use Cases)—short descriptions of potential science applications that could clearly address high priority decision-making needs of our stakeholders and that align with an integrated science focus. Use Cases are not actionable science planning documents, nor stand-alone scholarly works, but should be considered as innovative, next-generation science ideas that can be considered as potential components of science plans still under development. The goal of Use Case development was to (1) identify and characterize existing USGS scientific capacities and expertise that can support science goals and products, (2) identify opportunities to leverage current capacities for next-generation science, and (3) foster engagement across the entire Bureau to further refine the USGS strategy for EarthMAP and integrated predictive science. The Use Case development effort documented in this report was coordinated by the Use Case Development Team (UCDT), consisting of representatives from each region. The UCDT undertook five tasks: (1) develop a unified approach to engage bureau scientists consistently across all regions in aspirational thinking about what can be accomplished; (2) work with the regions and their Science Centers to generate an initial set of Use Cases, authored directly by scientists; (3) characterize, summarize, and document the initial set of Use Case submissions from authors to illuminate bureau-level demand for integrated science; (4) compare existing and needed capacities from the Use Case descriptions with preliminary results of the EarthMAP Capacity Assessment (Keisman and others, 2021); and (5) describe lessons learned from the Use Case development process and provide recommendations to inform future efforts to generate integrated science activities. This report outlines the approach the UCDT developed to solicit Use Cases from the regions and summarizes the high-level qualitative findings from this first-round effort. The UCDT received 36 Use Cases from the regions and identified potential points of convergence and commonalities considered useful in making connections among the participating scientists. The Southwest (SW) Region and the Rocky Mountain (RM) Region asked scientists to give special consideration to Use Cases with applicability to the Colorado River Basin, and seven of the Use Cases specifically named that geographic area as a focus. Coastal hazards and coastal resilience were identified in Use Cases from the Alaska (AK), Northeast (NE), and Southeast (SE) Regions. Aspects of wildfire and post-wildfire response were part of Uses Cases from AK, RM, and SW Regions. The greatest convergence of Use Case themes was related to conservation of public lands and waters, which is a powerful linkage lending strength to future collaborative efforts. The most common type of stakeholder decisions that would be informed by the Use Case science applications were related to adaptation, mitigation, and response (for example, how to increase the resilience of coastal communities to climate-related stressors and how to prevent or respond to harmful algal blooms). Other common types of decisions included water and land management decisions (including operational water management decisions such as reservoir operations and land use planning in the sagebrush biome), decisions about how to manage and conserve habitats and species, and risk management decisions (such as managing the post-wildfire flood risks). These decision types are not exclusive because many Use Cases cross categories. Use Case authors identified existing and needed science and technology capabilities required for Use Case implementation, which were then aligned to capabilities assessed in the EarthMAP Capacity Assessment (Keisman and others, 2021). Strong alignment was found for data and information integration approaches, modeling and prediction approaches, and capabilities related to delivery of actionable information. A majority of Use Cases indicated insufficient current capacity for needed data collection methods, data integration, and modeling and prediction approaches, whereas only 25 percent indicated insufficient capacity for actionable information delivery. Overall, many Use Case capacity demand gaps could potentially be met by existing bureau-wide capacity. In addition, nearly half of the Use Cases could potentially be implemented within 3 years if funding, capabilities, and personnel impediments were removed and science priorities were realigned. Several challenges emerged during the Use Case development process. The first challenge was developing an approach that was flexible enough to accommodate regional differences in planning and implementation, while also ensuring enough guidance to promote meaningful summary analyses. The UCDT encountered a strong demand for continuous communication and education to improve overall understanding of the integrated predictive science strategy. Another challenge was managing expectations about EarthMAP activities as a design effort that was not aligned to an immediate funding opportunity. Connecting the Use Cases to stakeholder needs without the opportunity for direct stakeholder engagement was also challenging. The last notable challenge was in obtaining consistent interpretation and characterization of the qualitative data housed in the narrative descriptions of Use Cases, written in different styles. Overall, the 36 Use Cases can serve as components of a road map for advancing integrated monitoring and predictive science throughout the USGS by revealing opportunities to (1) encourage cross-region initiatives that address shared interests in common themes by integrating similar Use Cases and through direct involvement of stakeholders in identifying needs and designing effective responses, (2) leverage the Use Cases to target investments that are aligned with the Bureau and Department of the Interior (DOI) priorities, (3) connect Use Cases and the results of the companion EarthMAP Capacity Assessment (Keisman and others, 2021) to identify potential priorities for capacity building investments, and (4) raise awareness of common integrated and interdisciplinary science interests within and across the regions through Use Case and Capacity Assessment summary outreach activities.

Open-File Report↗

USGS2004

Selenium loading through the Blackfoot River watershed--linking sources to ecosystem

The upper Blackfoot River watershed in southeast Idaho receives drainage from 11 of 16 phosphate mines that have extracted ore from the Phosphoria Formation, three of which are presently active. Toxic effects from selenium (Se), including death of livestock and deformity in aquatic birds, were documented locally in areas where phosphatic shales are exposed ( Piper et al., 2000 ; Presser et al., Chapter 11). Current drainage conditions are leading to Se bioaccumulation at concentrations that pose a risk to fish in the Blackfoot River and its tributaries (Hamilton et al., Chapter 18). A gaging station on the Blackfoot River was re-activated in April 2001 to assess hydrologic conditions and concentration, load, and speciation for Se discharges on a watershed scale. The gaging-station data are considered to represent regional drainage conditions in the upper Blackfoot River water- shed because of its location near the outlet of the watershed and directly upstream of the Blackfoot Reservoir. Watershed discharges for 2001 and 2002 were below minimum hydrologic conditions for the gage as documented by the historical record. Drought emergencies were declared in the area in both 2001 and 2002. Unmonitored diversions for irrigation that routinely take place during the snowmelt season also affected conditions downstream. Annual cycles in Se concentration, load, and selenate (Se 6+ ) reached maxima in the spring during the period of maximum flow at the gaging station. Thirty-seven to 44% of annual flow occurred dur- ing the three-month high-flow season (April through June) in 2001 and 56% of annual flow occurred during that time period in 2002. Extrapolation from historical hydrographs for average and wet years and a limited data set of regional Se concentrations for 2001 and 2002 indicated potential for a 3.6- to 7.4-fold increase in Se loading because of increased seasonal flows in the Blackfoot River watershed. Supplementation data indicate that: (a) the difference between total Se and dissolved Se, as a measure of the contribution of particulate Se, was < 10% except at the peak of con- centration when total Se was 18% more than dissolved Se; (b) selenite (Se 4+ ) represented less than 10% of the dissolved species during all months of 2001; and (c) dissolved Se was approximately a 50:50 mixture of selenate and organic selenide (operationally defined Se 2- ) during summer 2001 (June through August). Ecological risk based on regional Se drainage occurred during both the high- and low-flow seasons. Seventy to 83% of the Se load occurred during the high-flow season. During early May of both years, dissolved-Se concentrations exceeded the criterion for the protection of aquatic life and the ecological threshold of 5 gL 1 Se at which sub- stantive risk occurs. During the majority of the three-month high-flow season, dissolved- Se concentrations exceeded the 2 gL 1 Se concern level for aquatic biota. The Se concentration in suspended material during high flow in 2002 was within the range of marginal risk to aquatic life (2-4 gg 1 Se, dry weight). Selenate was the major species during peak flows, with both selenate and organic selenide being major species during relatively low-flow periods in summer. A change in speciation to reduced Se may indicate elevated biotic productivity during summer months and could result in enhanced Se uptake in food webs. In addition to the magnitude of regional Se release in the Blackfoot River watershed, Se concentrations in individual source drains and waste-rock seeps, and those predicted by experimental column leaching of proposed mining overburden materials, also indicate that drainage options that currently meet existing demands for phosphate mining cause eco- logical risk thresholds to be exceeded. At times, the drinking-water Se standard (50 g L 1 Se) and the criterion for hazardous Se waste (1000 L -1 Se) ( US Department of the Interior, 1998 ; US Environmental Protection Agency, 1987 ) are also exceeded. For water-years 2001 and 2002, seasonal increased input of water in the mining area resulted in increased Se transport, suggesting a mechanism of contamination that involves a significant Se reservoir. Hence, recognition and monitoring of Se loading to the envi- ronment on a mass balance basis (i.e. inputs, fluxes and storage within environmental media, and outputs) are essential to evaluating how to control Se concentrations within environmentally protective ranges ( Presser and Piper, 1998 ). In areas where release of Se to aquatic systems is anticipated as a product of future expansion of phosphate mining, continuous monitoring of flow and development of seasonal Se loading patterns would help to model watersheds in terms of sources, flow periods, and environmental-Se con- centrations that most influence bioavailability. These data, in turn, could be linked to Se- bioaccumulation models specific to food webs and vulnerable species of the impacted areas to accurately project ecological effects. Gaging at this site on the Blackfoot River is planned to continue in order to establish a long-term (>10 year) record of hydrologic conditions.

Idaho↗

USGS2008

Hydrogeology and water quality of the Leetown area, West Virginia

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

West Virginia↗

USGS2011

Borehole geophysical investigation of a formerly used defense site, Machiasport, Maine, 2003-2006

The U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, collected borehole geophysical logs in 18 boreholes and interpreted the data along with logs from 19 additional boreholes as part of an ongoing, collaborative investigation at three environmental restoration sites in Machiasport, Maine. These sites, located on hilltops overlooking the seacoast, formerly were used for military defense. At each of the sites, chlorinated solvents, used as part of defense-site operations, have contaminated the fractured-rock aquifer. Borehole geophysical techniques and hydraulic methods were used to characterize bedrock lithology, fractures, and hydraulic properties. In addition, each geophysical method was evaluated for effectiveness for site characterization and for potential application for further aquifer characterization and (or) evaluation of remediation efforts. Results of borehole geophysical logging indicate the subsurface is highly fractured, metavolcanic, intrusive, metasedimentary bedrock. Selected geophysical logs were cross-plotted to assess correlations between rock properties. These plots included combinations of gamma, acoustic reflectivity, electromagnetic induction conductivity, normal resistivity, and single-point resistance. The combined use of acoustic televiewer (ATV) imaging and natural gamma logs proved to be effective for delineating rock types. Each of the rock units in the study area could be mapped in the boreholes, on the basis of the gamma and ATV reflectivity signatures. The gamma and mean ATV reflectivity data were used along with the other geophysical logs for an integrated interpretation, yielding a determination of quartz monzonite, rhyolite, metasedimentary units, or diabase/gabbro rock types. The interpretation of rock types on the basis of the geophysical logs compared well to drilling logs and geologic mapping. These results may be helpful for refining the geologic framework at depth. A stereoplot of all fractures intersecting the boreholes indicates numerous fractures, a high proportion of steeply dipping fractures, and considerable variation in fracture orientation. Low-dip-angle fractures associated with unloading and exfoliation are also present, especially at a depth of less than 100 feet below the top of casing. These sub-horizontal fractures help to connect the steeply dipping fractures, making this a highly connected fracture network. The high variability in the fracture orientations also increases the connectivity of the fracture network. A preliminary comparison of all fracture data from all the boreholes suggests fracturing decreases with depth. Because all the boreholes were not drilled to the same depth, however, there is a clear sampling bias. Hence, the deepest boreholes are analyzed separately for fracture density. For the deepest boreholes in the study, the intensity of fracturing does not decline significantly with depth. It is possible the fractures observed in these boreholes become progressively tighter or closed with depth, but this is difficult to verify with the borehole methods used in this investigation. The fact that there are more sealed fractures at depth (observed in optical televiewer logs in some of the boreholes) may indicate less opening of the sealed fractures, less water moving through the rock, and less weathering of the fracture infilling minerals. Although the fracture orientation remained fairly constant with depth, differences in the fracture patterns for the three restoration sites indicate the orientation of fractures varies across the study area. The fractures in boreholes on Miller Mountain predominantly strike northwest-southeast, and to a lesser degree they strike northeast. The fractures on or near the summit of Howard Mountain strike predominantly east-west and dip north and south, and the fractures near the Transmitter Site strike northeast-southwest and dip northwest and southeast. The fracture populations for the boreholes on or near the summit of Howard Mountain show more variation than at the other two sites. This variation may be related to the proximity of the fault, which is northeast of the summit of Howard Mountain. In a side-by-side comparison of stereoplots from selected boreholes, there was no clear correspondence between fracture orientation and proximity to the fault. There is, however, a difference in the total populations of fractures for the boreholes on or near the summit of Howard Mountain and the boreholes near the Transmitter Site. Further to the southwest and further away from the fault, the fractures at the Transmitter Site predominantly strike northeast-southwest and northwest-southeast.Heat-pulse flowmeter (HPFM) logging was used to identify transmissive fractures and to estimate the hydraulic properties along the boreholes. Ambient downflow was measured in 13 boreholes and ambient upflow was measured in 9 boreholes. In nine other bedrock boreholes, the HPFM did not detect measurable vertical flow. The observed direction of vertical flow in the boreholes generally was consistent with the conceptual flow model of downward movement in recharge locations and upward flow in discharge locations or at breaks in the slope of land surface. Under low-rate pumping or injection rates [0.25 to 1 gallon per minute (gal/min)], one to three inflow zones were identified in each borehole. Two limitations of HPFM methods are (1) the HPFM can only identify zones within 1.5 to 2 orders of magnitude of the most transmissive zone in each borehole, and (2) the HPFM cannot detect flow rates less than 0.010 + or - 0.005 gal/min, which corresponds to a transmissivity of about 1 foot squared per day (ft2/d). Consequently, the HPFM is considered an effective tool for identifying the most transmissive fractures in a borehole, down to its detection level. Transmissivities below that cut-off must be measured with another method, such as packer testing or fluid-replacement logging. Where sufficient water-level and flowmeter data were available, HPFM results were numerically modeled. For each borehole model, the fracture location and measured flow rates were specified, and the head and transmissivity of each fracture zone were adjusted until a model fit was achieved with the interpreted ambient and stressed flow profiles. The transmissivities calculated by this method are similar to the results of an open-hole slug test; with the added information from the flowmeter, however, the head and transmissivity of discrete zones also can be determined. The discrete-interval transmissivities ranged from 0.16 to 330 ft2/d. The flowmeter-derived open-hole transmissivity, which is the combined total of each of the transmissive zones, ranged from 1 to 511 ft2/d. The whole-well open-hole transmissivity values determined with HPFM methods were compared to the results of open-hole hydraulic tests. Despite the fact that the flowmeter-derived transmissivities consistently were lower than the estimates derived from open-hole hydraulic tests alone, the correlation was very strong (with a coefficient of determination, R2, of 0.9866), indicating the HPFM method provides a reasonable estimate of transmissivities for the most transmissive fractures in the borehole. Geologic framework, fracture characterization, and estimates of hydraulic properties were interpreted together to characterize the fracture network. The data and interpretation presented in this report should provide information useful for site investigators as the conceptual site groundwater flow model is refined. Collectively, the results and the conceptual site model are important for evaluating remediation options and planning or implementing the design of a well field and borehole completions that will be adequate for monitoring flow, remediation efforts, groundwater levels, and (or) water quality. Similar kinds of borehole geophysical logging (specifically the borehole imaging, gamma, fluid logs, and HPFM) should be conducted in any newly installed boreholes and integrated with interpretations of any nearby boreholes. If boreholes are installed close to existing or other new boreholes, cross-hole flowmeter surveys may be appropriate and may help characterize the aquifer properties and connections between the boreholes.

Maine↗

USGS2023

Selenium hazards in the Salton Sea environment—Summary of current knowledge to inform future wetland management

Quaternary marine and continental shales in the western United States are sources of selenium that can be loaded into the aquatic environment through mining, agricultural, and energy production processes. The mobilization of selenium from shales through agricultural irrigation has been recognized since the 1930s; however, discovery of deformities in birds and other wildlife using agricultural habitats during the 1980s spurred studies to determine the extent and effects of the contamination. Through these early studies, researchers determined that biota in the Salton Sea drainage basin was at risk from legacy selenium contamination in the Colorado River watershed. The Salton Sea and its surrounding managed and unmanaged wetlands provide vital inland habitat and trophic support for diverse assemblages of resident and migratory wildlife, and understanding regional selenium hazards for these trust species is a priority for many Federal and State agencies. The modern Salton Sea is a shallow, landlocked saline lake in Riverside and Imperial Counties (not shown) of California that is sustained by irrigation return and perennial river inflow. Changes in water transfer agreements under the 2003 Quantification Settlement Agreement (QSA) have resulted in reduced irrigation flow, declining lake levels, and the evolution of unmanaged wetlands in areas where drains and rivers no longer reach the Salton Sea. These wetlands provide additional habitat for some species of concern, but their potential to increase selenium hazards for trust species is largely unknown. From the 1980s to 2020, efforts to document selenium contamination and effects throughout the region have resulted in a considerable amount of selenium data from the Salton Sea and its surrounding drainage basin; however, no long-term (greater than 20 years), consistent sampling program has been established, and all data have been collected by different entities using a variety of protocols and analytical techniques. This lack of coordination has been previously documented in regional management plans and has led to difficulty in reliably assessing selenium hazards in the Salton Sea environment. This report provides a summary of the available disparate selenium information collected from water, sediment, and biota in the Salton Sea region since the 1980s and to identify data gaps that need to be filled to understand the potential effects of selenium on species of concern, including federally endangered desert pupfish ( Cyprinodon macularius ) and Yuma Ridgway’s Rail ( Rallus obsoletus yumanensis; formerly Yuma Clapper Rail, Rallus longirostris yumanensis ). Available data from the Salton Sea drainage basin show that water from the Colorado River has the lowest selenium concentration of all surface water sources. All other surface water flowing into the Salton Sea has elevated selenium concentrations due to evaporation and evapotranspiration that occurs in agricultural fields and associated water delivery infrastructure or leaching of selenium from irrigated farmland soils. The Salton Sea has lower selenium concentrations because of various biogeochemical processes that recycle selenium into the sediment or volatilize it to the atmosphere; however, these mechanisms are not well defined, and it is not clear if selenium cycling will change in response to possible changes in the oxidation state of the Salton Sea bottom waters as water levels decline. Agricultural drains have the highest average selenium concentrations, but few drains have been sampled since changes in irrigation practices have occurred (due to the 2003 QSA). Groundwater selenium concentrations are variable; some wells south of the Salton Sea have selenium concentrations as high as 300 micrograms per liter (µg/L), whereas selenium concentrations are below detection in other wells. Groundwater and surface-water geothermal discharge zones around the margins of the Salton Sea and in unmanaged wetlands have not been studied in detail, and published selenium measurements are not available for these surface features. Selenium concentrations in the sediment of the Salton Sea drainage basin are highest in wetland particulate organic matter and the Salton Sea lakebed, indicating that removal of selenium from the water to the sediment has been a primary mechanism for keeping selenium concentrations low in the water column. Sediment selenium concentrations in wetlands are lower than in the Salton Sea but higher than inflowing drains and rivers, indicating the lentic wetland sites also may be important sinks for selenium because of biogeochemical processes. Sediment selenium data have not been collected in agricultural drains since changes in irrigation practices occurred (due to the 2003 QSA), and it is unknown if selenium sequestration from the water column has changed in these systems. We divided biological data into broad taxonomic categories, including primary producers, invertebrates, herpetofauna, mammals, fishes, and birds to facilitate evaluation of selenium concentrations and spatiotemporal trends observed in the Salton Sea. Overall, selenium concentrations were substantially greater in algae samples compared to all vascular plant samples combined. Median selenium concentrations in several invertebrate taxa (Chironomidae, Formicidae, Corixidae, Corbiculidae and Nereididae, and Decapoda) exceeded the maximum suggested dietary threshold of 3.0–4.0 micrograms per gram (µg/g) dry weight (dw) for predators consuming invertebrates in aquatic food webs. The greatest number of samples were collected from fish, and selenium distributions among species and locations showed that the range for most samples was lower than the U.S. Environmental Protection Agency selenium criterion for aquatic life (8.5 µg/g dw whole body, 11.3 µg/g dw fillets). The median selenium concentrations for whole body fish were below the selenium criterion in most locations, except for bairdiella ( Bairdiella icistia ) from the Salton Sea and irrigation drains, a few individual tilapia spp. (family Cichlidae, including genera Tilapia, Oreochromis , and their hybrids) from the river and river outlets, and several western mosquitofish ( Gambusia affinis ) and sailfin molly ( Poecilia latipinna ) from irrigation drain outlets. For avian samples combined among years and locations, median selenium concentrations in livers from all families except waders and Ibis (family Threskiornithidae) were higher than levels expected to cause selenium toxicosis (10–20 µg/g dw), and all median egg concentrations were above or near 6.0 μg/g dw, which is a conservative threshold value for reproductive impairment. Most knowledge gaps we identified for water, sediment, and biota were interrelated, and the use of integrated approaches to address knowledge gaps can provide greater insight into the drivers behind selenium hazards. Integrated water, sediment, and biota studies could help identify cost-effective management solutions that serve multiple purposes. A comprehensive analysis of the hydrology, biogeochemistry, and food-web processes in wetlands and other habitats can inform predictive models to identify drivers of selenium bioavailability, uptake from the environment and subsequent trophic transfer, ultimately forming the basis for experimental habitat management manipulations to minimize selenium hazards to wildlife. Furthermore, a comprehensive, long-term sampling and analytical laboratory plan would enable comparison of data among different entities that are sampling at the Salton Sea. Such efforts are well suited to help fill knowledge gaps that preclude understanding of selenium hazards and future management options for biota using Salton Sea habitats, including newly formed wetlands throughout the region. All data compiled for this report are available in two U.S. Geological Survey data releases: Groover and others (2022) for water and sediment samples and De La Cruz and others (2022) for biological samples. The data releases include all publicly available data for selenium concentrations in water, sediment, and biological samples collected in and around the Salton Sea, including the Coachella and Imperial Valleys. The data releases also include previously unpublished data.

California↗

USGS2020

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

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

California, Nevada↗

USGS2013

Native Prairie Adaptive Management: a multi region adaptive approach to invasive plant management on Fish and Wildlife Service owned native prairies

Much of the native prairie managed by the U.S. Fish and Wildlife Service (FWS) in the Prairie Pothole Region (PPR) of the northern Great Plains is extensively invaded by the introduced cool-season grasses, smooth brome ( Bromus inermis ) and Kentucky bluegrass ( Poa pratensis ). Management to suppress these invasive plants has had poor to inconsistent success. The central challenge to managers is selecting appropriate management actions in the face of biological and environmental uncertainties. In partnership with the FWS, the U.S. Geological Survey (USGS) developed an adaptive decision support framework to assist managers in selecting management actions under uncertainty and maximizing learning from management outcomes. This joint partnership is known as the Native Prairie Adaptive Management (NPAM) initiative. The NPAM decision framework is built around practical constraints faced by FWS refuge managers and includes identification of the management objective and strategies, analysis of uncertainty and construction of competing decision models, monitoring, and mechanisms for model feedback and decision selection. Nineteen FWS field stations, spanning four states of the PPR, have participated in the initiative. These FWS cooperators share a common management objective, available management strategies, and biological uncertainties. Though the scope is broad, the initiative interfaces with individual land managers who provide site-specific information and receive updated decision guidance that incorporates understanding gained from the collective experience of all cooperators. We describe the technical components of this approach, how the components integrate and inform each other, how data feedback from individual cooperators serves to reduce uncertainty across the whole region, and how a successful adaptive management project is coordinated and maintained on a large scale. During an initial scoping workshop, FWS cooperators developed a consensus management objective: increase the composition of native grasses and forbs on native sod while minimizing cost. Cooperators agreed that decision guidance should be provided annually and should account for local, real-time vegetation conditions observed on the ground. Over the course of development, two prototypes of the decision framework were considered. The final framework recognized four alternative actions that managers could take in any given year: (1) Graze—targeted use of grazing ungulates to achieve defoliation, (2) Burn—application of prescribed fire as the single form of defoliation, (3) Burn/Graze—a combination treatment, and (4) Rest—no action. The study area included northern mixed-grass and tallgrass prairie. Native vegetation in mixed–grass prairie has a strong cool-season component and thus the dominant native species have a phenology similar to that of smooth brome and Kentucky bluegrass, making management of those species challenging. In contrast, tallgrass prairie has a strong warm-season native component, leading to an existence of cool-season windows, periods of time in the fall and spring when cool‐season invasive grass species are actively growing and vulnerable to damage via select management actions, but warm‐season grass species are not active and are thus less susceptible to damage via the same actions. This dichotomy between prairie types necessitated the development of separate but parallel decision support systems for mixed-grass and tallgrass biomes. Management units are parcels of native prairie that receive a single management treatment at any one time over their entire extent. At any particular time, the vegetation state of each management unit is characterized by the amount of cover of native grasses and forbs and the type of invasive grass that is dominant. In addition, each unit has a defoliation state which reflects the number of years since the last defoliation event and an index to how intensively the unit was managed during the previous 7 years. State-transition models are used to predict the state of a management unit in year t+1 from its state in year t and a prescribed management action that was applied between the two monitoring events. Alternative models are built around key uncertainties that make choice of a management action difficult. Three uncertainties revolve around whether the effect of management actions depends on (1) type of dominant invader, (2) past defoliation history, and (3) level of invasion. Two additional uncertainties are considered when choosing a management action for tallgrass units: (4) the effectiveness of grazing within the cool-season window as a surrogate for burning when smooth brome is the dominant invader, and (5) the differential effect of active management outside the window as compared to rest. Because data on the probability of transitioning from one state to another under the various models were lacking, expert opinion and elicitation were used to parameterize the models. In addition, cooperators participated in elicitation exercises to extract their beliefs regarding the value of having native prairie compared to the cost of achieving it. Quantifying the subjective expression of utility in this way allowed for mathematical representation of the management objective into an objective function. By maximizing the objective function, cumulative utility is maximized, leading to the identification of a sequence of decisions that will achieve the management objective. The NPAM system adopted a vegetation monitoring protocol that was rapid, inexpensive, and familiar to many of the cooperators. The monitoring protocol served three purposes: (1) determining current vegetation and defoliation states of each unit, (2) evaluating progress toward the management objective, and (3) assessing predictive performance of the alternative models. The management year runs from September 1 to August 31. Management can be applied anytime during that period and monitoring takes places from late June to mid-August. Cooperators enter vegetation data and management information into a centralized database by August 25 of each year. Given the current state of the system (vegetation and defoliation states) and the current understanding of the system (or the belief state), identifying the current best management decision is a matter of looking up the combination (that is, system state and belief state) in the appropriate (mixed-grass or tallgrass) optimal decision table. Given complete uncertainty at the outset of decision-making, initial assignment of equal belief weights to each model was believed reasonable. The decisions in the optimal decision table that correspond to the current belief state constitute the current optimal decision policy. By August 31 of each year, individual cooperators are provided with a recommended management action for each of their management units for the upcoming management year. Upon receiving the management recommendations for their units, managers consider the recommendation, along with other relevant information, and at some point during the year one of the management alternatives is carried out. This iterative cycle of making and implementing a management decision, predicting the response, monitoring the outcome, comparing predicted and observed outcomes, updating model weights, and recommending a management action for the next cycle is expected to result in an accumulation of weight on a representative model of system dynamics, thereby increasing understanding needed to effectively manage native prairies. The NPAM system is now entering its second full year of complete operation, and represents one of only a few fully implemented applications of adaptive management within the U.S. Fish and Wildlife Service. NPAM is truly unique in that it originated from the ground up as a result of the leadership and steadfastness of several refuge biologists and managers confronted with a common problem. These biologists recognized that working together across a large landscape presented perhaps the best opportunity for halting and reversing the invasion of native grasslands by non-native cool-season grasses. Importantly, the NPAM system encapsulates the collective thinking and experience of tens if not hundreds of individuals who have battled this vexing problem for much of their careers. The NPAM initiative is rooted in principles of adaptive management, thereby affording the opportunity for grassland managers to pursue management objectives while acquiring information to reduce uncertainty and improve future management. The project introduced a number of technical innovations that will serve as templates for conservation efforts throughout and beyond the U.S. Fish and Wildlife Service. First, NPAM is an on-the-ground implementation of active adaptive management—possibly the first of its kind in conservation management—in which recommended management actions result from a prospective analysis of future learning (Williams, 1996). Second, by the use of dynamic optimization, NPAM demonstrates how decisions can be made that take into account possible future transitions of the system. Third, NPAM demonstrates how models of partial controllability are an effective means of accommodating unpredictable circumstances that cause a manager to follow a different course than was intended. Finally, the database developed for NPAM is an unparalleled system that enables the rapid integration of data from the field for the generation of ‘just-in-time’ management recommendations. In all, NPAM provides an example of how a science-management partnership can be forged to achieve large-scale conservation objectives.

Minnesota;Montana;North Dakota;South Dakota↗

USGS2024

Simulation of groundwater flow in the Long Island, New York regional aquifer system for pumping and recharge conditions from 1900 to 2019

The U.S. Geological Survey has developed a transient, groundwater-flow model that simulates hydrologic conditions in the Long Island aquifer system as part of an ongoing (since 2016) multiyear, cooperative investigation with the New York State Department of Environmental Conservation. The goals of this investigation are to assist stakeholders and resource managers to evaluate the response of the hydrologic system to changes in future hydraulic stresses. Responses in the hydrologic system include changes in water levels in the hydrogeologic units; discharge to streams, coastal waters, and subsurface infrastructure; and the extent of saline groundwater in the aquifers. Hydraulic stresses include future water-supply management and changes in land use and infrastructure. The numerical model synthesizes a diverse set of physiographic, geologic, climatic, land-use, and historical population, water use, and infrastructure data to physically represent the Long Island aquifer system from land surface to bedrock and to simulate annual hydrologic conditions between 1900 and 2019. A three-dimensional hydrogeologic framework was developed from existing and recently collected borehole geologic and geophysical data collected as part of a companion drilling program. Water-transmitting properties of the principal aquifer sediments were defined in three dimensions from new and existing lithologic logs. The distribution of recharge from precipitation was estimated from landscape characteristics and climate data. Anthropogenic recharge from wastewater, leaky infrastructure, and storm runoff were estimated from population, infrastructure, and pumping data. Water-use data, including well locations, depths, and pumping rates, were obtained from historical sources and records and used to estimate pumping stresses continuously in time and space, at an annual average time scale. The data were incorporated into a three-dimensional numerical model using the U.S. Geological Survey finite difference modeling code MODFLOW 6; the model encompassed all of Long Island and surrounding surface waters and simulated historical hydrologic conditions from 1900 to 2019. The calibration process involved trial and error adjustments using prior knowledge to improve general fit to observations followed by an inverse calibration to update and optimize input parameters, using an iterative ensemble smoother algorithm implemented in PEST++ version 5.0. This resulted in a model that generally was in good agreement with observed, dynamically varying hydrologic conditions from 1900 to 2019. The calibrated model was used to develop two base-case models for scenario testing of future, hypothetical conditions where one represented average-annual conditions, and one represented average-seasonal conditions from 2010 to 2019. The model representing average-annual conditions was modified further to represent an alternate sea-level position of 6 feet above the North American Vertical Datum of 1988, and the model representing average-seasonal conditions was modified to represent the average seasonal effects of a 5-year drought imposed upon current hydrologic conditions. Recharge is the sole source of water to the aquifer system; groundwater discharges to coastal water and streams and is withdrawn by pumped wells. Model-estimated annual recharge ranged from about 11 inches in 1965 to 41 inches in 1983. On average, from 2010 to 2019, about 23 percent of water was pumped from wells, and about 47 and 27 percent discharged to coastal waters and streams, respectively; the remaining 4 percent was water that moved into storage in the aquifer matrix. Water levels on Long Island vary naturally during time in response to changes in recharge; the amount of variation is largest in the interior of the island, in areas with highest water table altitudes near groundwater divides and lowest near streams and the coastal waters. The total range of water table altitudes on Long Island between 1900 and 2019 ranged from near 0 to more than 70 feet in western parts of Long Island. The largest range in altitudes is in New York City and is associated with areas of large historical withdrawals between the 1920s and the late 1980s. Water table altitudes generally varied by less than 10 feet in eastern Suffolk County, where the aquifer is under more natural conditions. Saltwater intrusion is of great concern on Long Island, particularly in western Long Island where both the unconfined and confined parts of the aquifer system have been intruded in response to large-scale groundwater withdrawals; however, the volume of freshwater in the islandwide aquifer system only has changed by about 5 percent between 1900 and 2019. The decadal change in the freshwater volume was largest during the early and mid-20th century, corresponding to the largest historical pumping, but that volume change did not exceed 1 percent. The negligible change in freshwater volume suggests that saltwater intrusion as of 2019 was limited at an islandwide scale but continues to occur in local areas of Queens and Nassau Counties, adversely affecting current water supplies and limiting future water supplies for affected communities. The regional groundwater model developed for this investigation is a tool that can be used to help determine the viability of current and future water supplies at a regional scale and can be used to support development of additional models at finer scale to support more focused assessments of groundwater sustainability.

New York↗

USGS2021

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

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

Open-File Report↗

USGS2017

Groundwater quality in the Western San Joaquin Valley study unit, 2010: California GAMA Priority Basin Project

Water quality in groundwater resources used for public drinking-water supply in the Western San Joaquin Valley (WSJV) was investigated by the USGS in cooperation with the California State Water Resources Control Board (SWRCB) as part of its Groundwater Ambient Monitoring and Assessment (GAMA) Program Priority Basin Project. The WSJV includes two study areas: the Delta–Mendota and Westside subbasins of the San Joaquin Valley groundwater basin. Study objectives for the WSJV study unit included two assessment types: (1) a status assessment yielding quantitative estimates of the current (2010) status of groundwater quality in the groundwater resources used for public drinking water, and (2) an evaluation of natural and anthropogenic factors that could be affecting the groundwater quality. The assessments characterized the quality of untreated groundwater, not the quality of treated drinking water delivered to consumers by water distributors. The status assessment was based on data collected from 43 wells sampled by the U.S. Geological Survey for the GAMA Priority Basin Project (USGS-GAMA) in 2010 and data compiled in the SWRCB Division of Drinking Water (SWRCB-DDW) database for 74 additional public-supply wells sampled for regulatory compliance purposes between 2007 and 2010. To provide context, concentrations of constituents measured in groundwater were compared to U.S. Environmental Protection Agency (EPA) and SWRCB-DDW regulatory and non-regulatory benchmarks for drinking-water quality. The status assessment used a spatially weighted, grid-based method to estimate the proportion of the groundwater resources used for public drinking water that has concentrations for particular constituents or class of constituents approaching or above benchmark concentrations. This method provides statistically unbiased results at the study-area scale within the WSJV study unit, and permits comparison of the two study areas to other areas assessed by the GAMA Priority Basin Project statewide. Groundwater resources used for public drinking water in the WSJV study unit are among the most saline and most affected by high concentrations of inorganic constituents of all groundwater resources used for public drinking water that have been assessed by the GAMA Priority Basin Project statewide. Among the 82 GAMA Priority Basin Project study areas statewide, the Delta–Mendota study area ranked above the 90th percentile for aquifer-scale proportions of groundwater resources having concentrations of total dissolved solids (TDS), sulfate, chloride, manganese, boron, chromium(VI), selenium, and strontium above benchmarks, and the Westside study area ranked above the 90th percentile for TDS, sulfate, manganese, and boron. In the WSJV study unit as a whole, one or more inorganic constituents with regulatory or non-regulatory, health-based benchmarks were present at concentrations above benchmarks in about 53 percent of the groundwater resources used for public drinking water, and one or more organic constituents with regulatory health-based benchmarks were detected at concentrations above benchmarks in about 3 percent of the resource. Individual constituents present at concentrations greater than health-based benchmarks in greater than 2 percent of groundwater resources used for public drinking water included: boron (51 percent, SWRCB-DDW notification level), chromium(VI) (25 percent, SWRCB-DDW maximum contaminant level (MCL)), arsenic (10 percent, EPA MCL), strontium (5.1 percent, EPA Lifetime health advisory level (HAL)), nitrate (3.9 percent, EPA MCL), molybdenum (3.8 percent, EPA HAL), selenium (2.6 percent, EPA MCL), and benzene (2.6 percent, SWRCB-DDW MCL). In addition, 50 percent of the resource had TDS concentrations greater than non-regulatory, aesthetic-based SWRCB-DDW upper secondary maximum contaminant level (SMCL), and 44 percent had manganese concentrations greater than the SWRCB-DDW SMCL. Natural and anthropogenic factors that could affect the groundwater quality were evaluated by using results from statistical testing of associations between constituent concentrations and values of potential explanatory factors, inferences from geochemical and age-dating tracer results, and by considering the water-quality results in the context of the hydrogeologic setting of the WSJV study unit. Natural factors, particularly the lithologies of the source areas for groundwater recharge and of the aquifers, were the dominant factors affecting groundwater quality in most of the WSJV study unit. However, where groundwater resources used for public supply included groundwater recharged in the modern era, mobilization of constituents by recharge of water used for irrigation also affected groundwater quality. Public-supply wells in the Westside study area had a median depth of 305 m and primarily tapped groundwater recharged hundreds to thousands of years ago, whereas public-supply wells in the Delta–Mendota study area had a median depth of 85 m and primarily tapped either groundwater recharged within the last 60 years or groundwater consisting of mixtures of this modern recharge and older recharge. Public-supply wells in the WSJV study unit are screened in the Tulare Formation and zones above and below the Corcoran Clay Member are used. The Tulare Formation primarily consists of alluvial sediments derived from the Coast Ranges to the west, except along the valley trough at the eastern margin of the WSJV study unit where the Tulare Formation consists of fluvial sands derived from the Sierra Nevada to the east. Groundwater from wells screened in the Sierra Nevada sands had manganese-reducing or manganese- and iron-reducing oxidation-reduction (redox) conditions. These redox conditions commonly were associated with elevated arsenic or molybdenum concentrations, and the dominance of arsenic(III) in the dissolved arsenic supports reductive dissolution of iron and manganese oxyhydroxides as the mechanism. In addition, groundwater from many wells screened in Sierra Nevada sands contained low concentrations of nitrite or ammonium, indicating reduction of nitrate by denitrification or dissimilatory processes, respectively. Geology of the Coast Ranges westward of the study unit strongly affects groundwater quality in the WSJV. Elevated concentrations of TDS, sulfate, boron, selenium and strontium in groundwater were primarily associated with aquifer sediments and recharge derived from areas of the Coast Ranges dominated by Cretaceous-to-Miocene age, organic-rich, reduced marine shales, known as the source of selenium in WSJV soils, surface water, and groundwater. Low sulfur-isotopic values (δ34S) of dissolved sulfate indicate that the sulfate was largely derived from oxidation of biogenic pyrite from the shales, and correlations with trace element concentrations, geologic setting, and groundwater geochemical modeling indicated that distributions of sulfate, strontium, and selenium in groundwater were controlled by dissolution of secondary sulfate minerals in soils and sediments. Elevated concentrations of chromium(VI) were primarily associated with aquifer sediments and recharge derived from areas of the Coast Ranges dominated by the Franciscan Complex and ultramafic rocks. The Franciscan Complex also has boron-rich, sodium-chloride dominated hydrothermal fluids that contribute to elevated concentrations of boron and TDS. Groundwater from wells screened in Coast Ranges alluvium was primarily oxic and relatively alkaline (median pH value of 7.55) in the Delta–Mendota study area, and primarily nitrate-reducing or suboxic and alkaline (median pH value of 8.4) in the Westside study area. Many groundwater samples from those wells have elevated concentrations of arsenic(V), molybdenum, selenium, or chromium(VI), consistent with desorption of metal oxyanions from mineral surfaces under those geochemical conditions. High concentrations of benzene were associated with deep wells located in the vicinity of petroleum deposits at the southern end of the Westside study area. Groundwater from these wells had premodern age and anoxic geochemical conditions, and the ratios among concentrations of hydrocarbon constituents were different from ratios found in fuels and combustion products, which is consistent with a geogenic source for the benzene rather than contamination from anthropogenic sources. Water stable-isotope compositions, groundwater recharge temperatures, and groundwater ages were used to infer four types of groundwater: (1) groundwater derived from natural recharge of water from major rivers draining the Sierra Nevada; (2) groundwater primarily derived from natural recharge of water from Coast Ranges runoff; (3) groundwater derived from recharge of pumped groundwater applied to the land surface for irrigation; and (4) groundwater derived from recharge during a period of much cooler paleoclimate. Water previously used for irrigation was found both above and below the Corcoran Clay, supporting earlier inferences that this clay member is no longer a robust confining unit. Recharge of water used for irrigation has direct and indirect effects on groundwater quality. Elevated nitrate concentrations and detections of herbicides and fumigants in the Delta–Mendota study area generally were associated with greater agricultural land use near the well and with water recharged during the last 60 years. However, the extent of the groundwater resource affected by agricultural sources of nitrate was limited by groundwater redox conditions sufficient to reduce nitrate. The detection frequency of perchlorate in Delta–Mendota groundwater was greater than expected for natural conditions. Perchlorate, nitrate, selenium, and strontium concentrations were correlated with one another and were greater in groundwater inferred to be recharge of previously pumped groundwater used for irrigation. The source of the perchlorate, selenium, and strontium appears to be salts deposited in the soils and sediments of the arid WSJV that are dissolved and flushed into groundwater by the increased amount of recharge caused by irrigation. In the Delta–Mendota study area, the groundwater with elevated concentrations of selenium was found deeper in the aquifer system than it was reported by a previous study 25 years earlier, suggesting that this transient front of groundwater with elevated concentrations of constituents derived from dissolution of soil salts by irrigation recharge is moving down through the aquifer system and is now reaching the depth zone used for public drinking water supply.

California↗

USGS2010

Effects of selected low-impact-development (LID) techniques on water quality and quantity in the Ipswich River Basin, Massachusetts: Field and modeling studies

During the months of August and September, flows in the Ipswich River, Massachusetts, dramatically decrease largely due to groundwater withdrawals needed to meet increased residential and commercial water demands. In the summer, rates of groundwater recharge are lower than during the rest of the year, and water demands are higher. From 2005 to 2008, the U.S. Geological Survey, in a cooperative funding agreement with the Massachusetts Department of Conservation and Recreation, monitored small-scale installations of low-impact-development (LID) enhancements designed to diminish the effects of storm runoff on the quantity and quality of surface water and groundwater. Funding for the studies also was contributed by the U.S. Environmental Protection Agency’s Targeted Watersheds Grant Program through a financial assistance agreement with Massachusetts Department of Conservation and Recreation. The monitoring studies examined the effects of (1) replacing an impervious parking lot surface with a porous surface on groundwater quality, (2) installing rain gardens and porous pavement in a neighborhood of 3 acres on the quantity and quality of stormwater runoff, and (3) installing a 3,000-square foot (ft2) green roof on the quantity and quality of stormwater runoff. In addition, the effects of broad-scale implementation of LID techniques, reduced water withdrawals, and water-conservation measures on streamflow in large areas of the basin were simulated using the U.S. Geological Survey’s Ipswich River Basin model. From June 2005 to 2007, groundwater quality was monitored at the Silver Lake town beach parking lot in Wilmington, MA, prior to and following the replacement of the conventional, impervious-asphalt surface with a porous surface consisting primarily of porous asphalt and porous pavers. Changes in the concentrations of the water-quality constituents, phosphorus, nitrogen, cadmium, chromium, copper, lead, nickel, zinc, and total petroleum hydrocarbons, were monitored. Increased infiltration of precipitation did not result in discernible increases in concentrations of these potential groundwater contaminants. Concentrations of dissolved oxygen increased slightly in groundwater profiles following the removal of the impervious asphalt parking lot surface. In Wilmington, MA, in a 3-acre neighborhood, stormwater runoff volume and quality were monitored to determine the ability of selected LID enhancements (rain gardens and porous paving stones) to reduce flows and loads of the above constituents to Silver Lake. Flow-proportional water-quality samples were analyzed for nutrients, metals, total petroleum hydrocarbons, and total-coliform and Escherichia coli bacteria. In general, when all storms were considered, no substantial decreases were observed in runoff volume as a result of installing LID enhancements. However, the relation between rainfall and runoff did provide some insight into how the LID enhancements affected the effective impervious area for the neighborhood. A decrease in runoff was observed for storms of 0.2 inches (in.) or less of precipitation, which indicated a reduction in effective impervious area from approximately 10 percent to about 4.5 percent for the 3-acre area. Water-quality-monitoring results were inconclusive; there were no statistically significant differences in concentrations or loads when the pre- and post-installation-period samples were compared. Three factors were probably most important in minimizing differences: (1) the small decrease in effective impervious area, (2) the differences in the size of storms sampled for water-quality constituents before and after installation of the infiltration enhancing measures, and (3) small sample sizes. In a third field study, the characteristics of runoff from a vegetated “green” roof and a conventional, rubber-membrane roof were compared. The amount of precipitation and the length of the antecedent dry period were the two primary factors affecting the green roof’s water-storage capacity. The green roof retained more than 50 percent of the precipitation from storms with 0.04 to 1.0 in. of rain. Approximately 95 percent of the precipitation from one storm of nearly 2 in. was retained by the green roof. On the rubber-membrane roof, only a small, shallow puddle of insubstantial volume ever remained after a storm. Bulk precipitation from 10 storms was monitored for the same constituents (nutrients, metals, and total petroleum hydrocarbons) as the roof runoff, and the results were compared with those for roof-runoff samples. The use of fertilizers to help establish the vegetation during the study probably distorted any effect the plants and growing medium may have had on the retention of target analytes. As a result of the fertilizer and growing medium chemistry, median concentrations of total nitrogen, total phosphorus, cadmium, copper, and nickel in runoff from the green roof were greater than in the runoff from the conventional roof or in bulk precipitation. Concentrations of lead and zinc were greater in runoff from the conventional roof, probably a result of passage through the old, metal drainpipes. Simulations of the effects of LID on streamflow in the Ipswich River Basin were conducted with a previously calibrated Hydrological Simulation Program-FORTRAN (HSPF) precipitation-runoff model. Simulations were conducted at multiple spatial scales to evaluate the effects of (1) updated water withdrawals for the towns of Reading and Wilmington; (2) potential land-use changes at buildout (potential future development); (3) effective impervious area reductions upstream from the South Middleton streamgage to represent the effects of widespread implementation of LID retrofit techniques; (4) basin-scale water withdrawal reductions scaled up (expanded to the town level) from water-conservation pilot programs conducted by the Massachusetts Department of Conservation and Recreation; and (5) land-use change and LID techniques at a local scale, which is smaller than the HSPF subbasin. Effects on streamflow generally were evaluated by comparing results of two or more related simulations for selected reaches in the basin; thus, relative rather than absolute changes in simulated flow were the focus of the assessment. Simulations indicated that reduced withdrawals for the towns of Reading and Wilmington led to substantially higher medium and low flows in most of the reaches upstream from the South Middleton streamgage. Simulations of water-conservation measures resulted in negligible effects on streamflow. Overall, simulations indicated that spatial scale is an important factor in determining the effects of land-use change and LID practices on streamflow. Potential land-use changes at buildout had modest (percent differences of less than 20 percent) effects on streamflow in most subbasins because relatively little land in the basin was available for development (about 17 percent); moreover, most of the available open land is zoned for low-density residential development, and this land-use category was simulated to contain relatively little effective impervious area and to be similar hydrologically to the forested land in place prior to development. Results of the simulations conducted to evaluate widespread effective impervious area reductions upstream from the South Middleton streamgage indicated that the percentage of urban land use and associated effective impervious area was too small for a 50-percent reduction of effective impervious area to appreciably affect streamflow (percent differences of less than 20 percent) in most subbasins. In contrast, the results of the hypothetical local-scale simulations indicated that for smaller streams, where the percentage of urban land use and associated effective impervious area in the drainage area may be substantially higher, land-use change, development patterns, and LID practices potentially have much greater effects on streamflow. Modeling results also indicated that LID was potentially most beneficial for minimizing streamflow alteration when applied to dense urban development, largely because larger tracts of effective impervious area were available for reduction than were available for other land-use categories. For example, commercial-industrial-transportation land use is composed of 37 percent pervious area and 63 percent effective impervious area in the HSPF model, whereas low-density residential area is composed of 97.5 percent pervious area and only 2.5 percent effective impervious area. Field and modeling studies concurred in the assessment that LID enhancements would likely have the greatest effect on decreasing stormwater runoff when broadly applied to highly impervious urban areas. A measurable effect for small rainfall events (less than 0.25 inch) was determined in the small, highly pervious area that was monitored in this study, but the volume difference was not great.

Massachusetts↗

USGS2014

Dissolved-solids sources, loads, yields, and concentrations in streams of the conterminous United States

Recent studies have shown that excessive dissolved-solids concentrations in water can have adverse effects on the environment and on agricultural, domestic, municipal, and industrial water users. Such effects motivated the U.S. Geological Survey’s National Water Quality Assessment Program to develop a SPAtially-Referenced Regression on Watershed Attributes (SPARROW) model that has improved the understanding of sources, loads, yields, and concentrations of dissolved solids in streams of the conterminous United States. Using the SPARROW model, long-term mean annual dissolved-solids loads from 2,560 water-quality monitoring stations were statistically related to several spatial datasets that are surrogates for dissolved-solids sources and land-to-water delivery processes. Specifically, sources in the model included variables representing geologic materials, road deicers, urban lands, cultivated lands, and pasture lands. Transport of dissolved solids from these sources was modulated by land-to-water delivery variables that represent precipitation, streamflow, soil, vegetation, terrain, population, irrigation, and artificial drainage characteristics. Where appropriate, the load estimates, source variables, and transport variables were statistically adjusted to represent conditions for the base year 2000. The nonlinear least-squares estimated SPARROW model was used to predict long-term mean annual conditions for dissolved-solids sources, loads, yields, and concentrations in a digital hydrologic network representing nearly 66,000 stream reaches and their corresponding incremental catchments that drain the Nation. Nationwide, the predominant source of dissolved solids yielded from incremental catchments and delivered to local streams is geologic materials in 89 percent of the catchments, road deicers in 5 percent of the catchments, pasture lands in 3 percent of the catchments, urban lands in 2 percent of the catchments, and cultivated lands in 1 percent of the catchments. Whereas incremental catchments with dissolved solids that originated predominantly from geologic sources or from urban lands are found across much of the Nation, incremental catchments with dissolved solids yields that originated predominantly from road deicers are largely found in the Northeast, and incremental catchments with dissolved solids that originated predominantly from cultivated or pasture lands are largely found in the West. The total amount of dissolved solids delivered to the Nation’s streams is 271.9 million metric tons (Mt) annually, of which 194.2 million Mt (71.4%) come from geologic sources, 37.7 million Mt (13.9%) come from road deicers, 18.2 million Mt (6.7%) come from pasture lands, 13.9 million Mt (5.1%) come from urban lands, and 7.9 million Mt (2.9%) come from cultivated lands. Nationwide, the median incremental-catchment yield delivered to local streams is 26 metric tons per year per square kilometer [(Mt/yr)/km 2 ]. Ten percent of the incremental catchments yield less than 4 (Mt/yr)/km 2 , and 10 percent yield more than 90 (Mt/yr)/km 2 . Incremental-catchment yields greater than 50 (Mt/yr)/km 2 mostly occur along the northern part of the West Coast and in a crescent shaped band south of the Great Lakes. For example, the median incremental-catchment yield is 81 (Mt/yr)/km 2 for the Great Lakes, 78 (Mt/yr)/km 2 for the Ohio, and 74 (Mt/yr)/km 2 for the Upper Mississippi water-resources regions. Incremental-catchment yields less than 10 (Mt/yr)/km 2 mostly occur in a wide band across the arid lowland of the interior West that excludes areas along the coast and the extensive, higher mountain ranges. For example, the median incremental-catchment yield is 3 (Mt/yr)/km 2 for the Lower Colorado, 5 (Mt/yr)/km 2 for the Rio Grande, and 8 (Mt/yr)/km 2 for the Great Basin water-resources regions. Predicted incremental loads were cascaded down through the reach network, with loads accumulating from reach to reach. For most stream reaches, the entire incremental load of dissolved solids delivered to the reach was transported to either the ocean or to one of the large streams flowing along the U.S. international boundary without losses occurring along the way. The exceptions to this include streams in the southwestern part of the country, such as the Colorado River, Rio Grande, and streams of internally drained drainages in the Great Basin, where dissolved-solids loads decreased through streamflow diversion for off-stream use, or by infiltration through the streambed. Long-term mean annual flow-weighted concentrations were derived from the predicted accumulated-load and stream-discharge data. Widespread low concentrations, generally less than 100 milligrams per liter (mg/L), occur in many reaches of the New England, South Atlantic-Gulf, and Pacific Northwest water-resources regions as a result of moderate dissolved-solids yields and high runoff rates. Widespread moderate concentrations, generally between 100 and 500 mg/L, occur in many reaches of the Great Lakes, Ohio, and Upper Mississippi River water-resources regions. Whereas dissolved-solids yields are generally high in these regions, runoff rates are also high, which helps moderate concentrations in these regions. Widespread higher concentrations, generally greater than 500 mg/L, occur across a belt of reaches that extends almost continuously from Canada to Mexico in the Midwest, cutting through the Souris-Red-Rainy, Missouri, Arkansas-White-Red, Texas-Gulf, and Rio Grande water-resources regions. Although dissolved-solids yields are moderate to low in these areas, low runoff rates result in the high concentrations for these areas. In 12.6 percent of the Nation’s stream reaches, predicted concentrations of dissolved solids exceed 500 mg/L, the U.S. Environmental Protection Agency’s secondary, nonenforceable drinking water standard. While this standard provides a metric for evaluating predicted concentrations in the context of drinking-water supplies, it should be noted that it only applies to drinking water actually served to customers by water utilities, and it does not apply to all stream reaches in the Nation nor does it apply during times when water is not being withdrawn for use. Exceedance of 500 mg/L is more pronounced in certain water-resources regions than others. For example, about half of the reaches in the Souris-Red-Rainy region have concentrations predicted to exceed 500 mg/L, and between 25 and 37 percent of the reaches in the Missouri, Arkansas-White-Red, Texas-Gulf, Rio Grande, and Lower Colorado regions are predicted to exceed 500 mg/L. Development of stream-load data for use in the SPARROW model also provided long-term temporal trend information in dissolved-solids concentrations at the monitoring stations for their period of record, which was constrained between 1980 and 2009. For the 2,560 monitoring stations used in this study, long-term trends in flow-adjusted dissolved-solids concentrations increased over time at 23 percent of the stations, decreased at 18 percent of the stations, and did not change over time at 59 percent of the stations. Long-term trends show a strong regional spatial pattern where from the western parts of the Great Plains to the West Coast, concentrations mostly either did not change or decreased over time, and from the eastern parts of the Great Plains to the East Coast, concentrations mostly either did not change or increased over time. Results from the trend analysis and from the SPARROW model indicate that, compared to monitoring stations with no trends or decreasing trends, stations with increasing trends are associated with a smaller percentage of the predicted dissolved-solids load originating from geologic sources, and a larger percentage originating from urban lands and road deicers. Conversely, compared to stations with increasing trends or no trends, stations with decreasing trends have a larger percentage of the predicted dissolved-solids load originating from geologic sources and a smaller percentage originating from urban lands and road deicers. Stations with decreasing trends also have larger percentages of predicted dissolved-solids load originating from cultivated lands and pasture lands, compared to stations with increasing trends or no trends.

Scientific Investigations Report↗