Search USGS⌕ Search

SEARCH · Search USGS

Results for “Risk Analysis”

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

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

1,041 records · Page 58Linked to original sources

Determination of steroid hormones and related compounds in filtered and unfiltered water by solid-phase extraction, derivatization, and gas chromatography with tandem mass spectrometry

A new analytical method has been developed and implemented at the U.S. Geological Survey National Water Quality Laboratory that determines a suite of 20 steroid hormones and related compounds in filtered water (using laboratory schedule 2434) and in unfiltered water (using laboratory schedule 4434). This report documents the procedures and initial performance data for the method and provides guidance on application of the method and considerations of data quality in relation to data interpretation. The analytical method determines 6 natural and 3 synthetic estrogen compounds, 6 natural androgens, 1 natural and 1 synthetic progestin compound, and 2 sterols: cholesterol and 3--coprostanol. These two sterols have limited biological activity but typically are abundant in wastewater effluents and serve as useful tracers. Bisphenol A, an industrial chemical used primarily to produce polycarbonate plastic and epoxy resins and that has been shown to have estrogenic activity, also is determined by the method. A technique referred to as isotope-dilution quantification is used to improve quantitative accuracy by accounting for sample-specific procedural losses in the determined analyte concentration. Briefly, deuterium- or carbon-13-labeled isotope-dilution standards (IDSs), all of which are direct or chemically similar isotopic analogs of the method analytes, are added to all environmental and quality-control and quality-assurance samples before extraction. Method analytes and IDS compounds are isolated from filtered or unfiltered water by solid-phase extraction onto an octadecylsilyl disk, overlain with a graded glass-fiber filter to facilitate extraction of unfiltered sample matrices. The disks are eluted with methanol, and the extract is evaporated to dryness, reconstituted in solvent, passed through a Florisil solid-phase extraction column to remove polar organic interferences, and again evaporated to dryness in a reaction vial. The method compounds are reacted with activated -methyl--trimethylsilyl trifluoroacetamide at 65 degrees Celsius for 1 hour to form trimethylsilyl or trimethylsilyl-enol ether derivatives that are more amenable to gas chromatographic separation than the underivatized compounds. Analysis is carried out by gas chromatography with tandem mass spectrometry using calibration standards that are derivatized concurrently with the sample extracts. Analyte concentrations are quantified relative to specific IDS compounds in the sample, which directly compensate for procedural losses (incomplete recovery) in the determined and reported analyte concentrations. Thus, reported analyte concentrations (or analyte recoveries for spiked samples) are corrected based on recovery of the corresponding IDS compound during the quantification process. Recovery for each IDS compound is reported for each sample and represents an absolute recovery in a manner comparable to surrogate recoveries for other organic methods used by the National Water Quality Laboratory. Thus, IDS recoveries provide a useful tool for evaluating sample-specific analytical performance from an absolute mass recovery standpoint. IDS absolute recovery will differ and typically be lower than the corresponding analyte’s method recovery in spiked samples. However, additional correction of reported analyte concentrations is unnecessary and inappropriate because the analyte concentration (or recovery) already is compensated for by the isotope-dilution quantification procedure. Method analytes were spiked at 10 and 100 nanograms per liter (ng/L) for most analytes (10 times greater spike levels were used for bisphenol A and 100 times greater spike levels were used for 3--coprostanol and cholesterol) into the following validation-sample matrices: reagent water, wastewater-affected surface water, a secondary-treated wastewater effluent, and a primary (no biological treatment) wastewater effluent. Overall method recovery for all analytes in these matrices averaged 100 percent, with overall relative standard deviation of 28 percent. Mean recoveries of the 20 individual analytes for spiked reagent-water samples prepared along with field samples and analyzed in 2009–2010 ranged from 84–104 percent, with relative standard deviations of 6–36 percent. Concentrations for two analytes, equilin and progesterone, are reported as estimated because these analytes had excessive bias or variability, or both. Additional database coding is applied to other reported analyte data as needed, based on sample-specific IDS recovery performance. Detection levels were derived statistically by fortifying reagent water at six different levels (0.1 to 4 ng/L) and range from about 0.4 to 4 ng/L for 16 analytes. Interim reporting levels applied to analytes in this report range from 0.8 to 8 ng/L. Bisphenol A and the sterols (cholesterol and 3-beta-coprostanol) were consistently detected in laboratory and field blanks. The minimum reporting levels were set at 100 ng/L for bisphenol A and at 200 ng/L for the two sterols to prevent any bias associated with the presence of these compounds in the blanks. A minimum reporting level of 2 ng/L was set for 11-ketotestosterone to minimize false positive risk from an interfering siloxane compound emanating as chromatographic-column bleed, from vial septum material, or from other sources at no more than 1 ng/L.

Techniques and Methods↗

Magnitude of flood flows for selected annual exceedance probabilities in Rhode Island through 2010

Heavy persistent rains from late February through March 2010 caused severe widespread flooding in Rhode Island that set or nearly set record flows and water levels at many long-term streamgages in the State. In response, the U.S. Geological Survey, in partnership with the Federal Emergency Management Agency, conducted a study to update estimates of flood magnitudes at streamgages and regional equations for estimating flood flows at ungaged locations. This report provides information needed for flood plain management, transportation infrastructure design, flood insurance studies, and other purposes that can help minimize future flood damages and risks. The magnitudes of floods were determined from the annual peak flows at 43 streamgages in Rhode Island (20 sites), Connecticut (14 sites), and Massachusetts (9 sites) using the standard Bulletin 17B log-Pearson type III method and a modification of this method called the expected moments algorithm (EMA) for 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent annual exceedance probability (AEP) floods. Annual-peak flows were analyzed for the period of record through the 2010 water year; however, records were extended at 23 streamgages using the maintenance of variance extension (MOVE) procedure to best represent the longest period possible for determining the generalized skew and flood magnitudes. Generalized least square regression equations were developed from the flood quantiles computed at 41 streamgages (2 streamgages in Rhode Island with reported flood quantiles were not used in the regional regression because of regulation or redundancy) and their respective basin characteristics to estimate magnitude of floods at ungaged sites. Of 55 basin characteristics evaluated as potential explanatory variables, 3 were statistically significant—drainage area, stream density, and basin storage. The pseudo-coefficient of determination (pseudo- R2 ) indicates these three explanatory variables explain 95 to 96 percent of the variance in the flood magnitudes from 20- to 0.2-percent AEPs. Estimates of uncertainty of the at-site and regression flood magnitudes are provided and were combined with their respective estimated flood quantiles to improve estimates of flood flows at streamgages. This region has a long history of urban development, which is considered to have an important effect on flood flows. This study includes basins that have an impervious area ranging from 0.5 to 37 percent. Although imperviousness provided some explanatory power in the regression, it was not statistically significant at the 95-percent confidence level for any of the AEPs examined. Influence of urbanization on flood flows indicates a complex interaction with other characteristics that confounds a statistical explanation of its effects. Standard methods for calculating magnitude of floods for given AEP are based on the assumption of stationarity, that is, the annual peak flows exhibit no significant trend over time. A subset of 16 streamgages with 70 or more years of unregulated systematic record indicates all but 4 streamgages have a statistically significant positive trend at the 95-percent confidence level; three of these are statistically significant at about the 90-percent confidence level or above. If the trend continues linearly in time, the estimated magnitude of floods for any AEP, on average, will increase by 6, 13, and 21 percent in 10, 20, and 30 years' time, respectively. In 2010, new peaks of record were set at 18 of the 21 active streamgages in Rhode Island. The updated flood frequency analysis indicates the peaks at these streamgages ranged from 2- to 0.2-percent AEP. Many streamgages in the State peaked at a 0.5- and 0.2-percent AEP, except for streamgages in the Blackstone River Basin, which peaked from a 4- to 2-percent AEP.

Rhode Island↗

Understanding the golden eagle and bald eagle sensory worlds to enhance detection and response to wind turbines

The objective for this study was to measure the auditory and visual physiology of Golden and Bald Eagles in order to use eagle sensory capabilities to inform the design of potential deterrent stimuli that could be used to reduce eagle/turbine collisions with wind turbines. The rationale for this approach is that sensory systems of any organism will limit the capability of that organism to perceive aspects of the world around it. Moreover, species can differ dramatically in their sensory physiology so it is important to examine these characteristics in the species of concern, rather than relying on data from similar birds. Our project consisted of two main phases. The first phase was the acquisition and analysis of visual and auditory information from Golden and Bald Eagles in rehabilitation centers. This was performed in order to identify light and sound stimuli tuned to sensitive areas in the eagle’s sensory systems. The second phase of the project was to present these stimuli to both species of eagles in a behavioral experiment to identify which stimuli would be the most effective in changing the behaviors of the eagles. Results of phase one indicated that the visual system of the Golden Eagle strongly absorbs ultraviolet light, making it unlikely the Golden Eagle (and most likely the Bald Eagle) will detect ultraviolet light signals. The Golden and Bald Eagles have differences in the sensitivities of their visual systems to light within the eye, but mathematical models indicate that both species are able to detect indigo/blue and orange/red light produced by LEDs (light emitting diodes) very well. We also found that both species of eagles have a blind spot above their head. This blind spot is particularly large in Golden Eagles due to a pronounced brow ridge above the eyes. This blind spot will result in the inability of a Golden Eagle to see something in front of it when its head is pointed down during flight – as might happen while hunting (i.e. searching the ground for prey). As such, the blind spot may increase the chance of collision with wind turbines if the eagle is actively hunting. This problem is less pronounced in Bald Eagles because their blind spot is smaller than in the Golden Eagles and their foraging strategy is different. Results of phase one also indicated that the auditory systems of the Golden and Bald Eagles respond differently to a variety of sounds (static tones, static chords (i.e. stacked tones), and sounds with dynamic changes through amplitude modulation or frequency modulation). Both species’ auditory systems responded strongly to tones across a wide range of frequencies (0.5 – 5kHz ), however the Bald Eagles’ auditory system was much better at processing complex sounds with dynamic rapid changes in amplitude or frequency modulation than the Golden Eagle. All of these sounds were then played with two types of noise in the background (white or pink). White noise more closely resembles the sound of wind and pink noise more closely resembles wind turbines or other sources of anthropogenic noise. Most sounds were more strongly masked by pink noise than by white noise, but several sounds (especially sounds with rapid modulation changes) showed little or no masking, indicating these were good candidate signals. However, even though rapidly changing sounds are less subject to noise masking, they are also less strongly processed by the Golden Eagle auditory system. This tradeoff does not exist in Bald Eagles because individuals of this species are very good at processing rapidly changing sounds. Given that Golden Eagle populations are at greater risk than Bald Eagle populations, we suggest that the most efficient alerting sound stimuli used in deterrent systems should be complex sounds that do not change very rapidly. We identified candidate light (indigo/blue and orange/red LED lights) and sound (sinusoidal frequency modulated sound, linear frequency sweeps, amplitude modulated sound, and a mistuned harmonic stack) stimuli that both eagle species sensory systems are highly sensitive to. Results of phase two, in which we presented these stimuli to eagles in a behavioral experiment, indicated that eagles behaviorally responded to all the stimuli presented, but at varying degrees. The Golden Eagles, especially, elicited higher rates of visual exploratory behavior with a flashing blue light stimulus and all sound stimuli. Our results suggest using these stimuli in field-testing of light/sound eagle deterrent systems on wind turbines. The eagles showed lower rates of behavior over the course of an experiment, suggesting either that they habituated to our stimuli or were initially stressed by the setup of the behavioral tests. These results underscore the need to test for habituation effects. Nonetheless, habitation to the stimuli in these field tests would likely be reduced by the use of random presentations of the four sounds and if possible random presentation of the candidate lights.

Report↗

Least Bell's Vireos and Southwestern Willow Flycatchers at the San Luis Rey flood risk management project area in San Diego County, California: Breeding activities and habitat use—2021 Annual report

Executive Summary Surveys and monitoring for the endangered Least Bell’s Vireo ( Vireo bellii pusillus ; vireo) were done at the San Luis Rey Flood Risk Management Project Area (Project Area) in the city of Oceanside, San Diego County, California, between April 4 and August 4, 2021. We completed four protocol surveys during the breeding season, supplemented by weekly territory monitoring visits. We identified a total of 122 territorial male vireos; 111 were confirmed as paired and 8 were confirmed as single males. For the remaining three territories, we were unable to confirm pair status. Five transient vireos were detected in 2021. The vireo population in the Project Area decreased by 24 percent from 2020 to 2021. Vireo populations decreased across San Diego County, with a 14-percent decrease documented at Marine Corps Base Camp Pendleton (MCBCP); a 5-percent decrease on the Otay River; a 6-percent decrease on the middle San Luis Rey River; and a 44-percent decrease at Marine Corps Air Station (although this decrease was likely exaggerated by large-scale vegetation clearing that occurred prior to the 2021 breeding season). We used an index of treatment (Treatment Index) to evaluate the impact of on-going vegetation clearing on the Project Area vireo population. The Treatment Index measures the cumulative effect of vegetation treatment within a territory (since 2005) by using the percent area treated weighted by the number of years since treatment. We found that the Treatment Index for unoccupied habitat was more than two times that of occupied habitat, indicating that vireos selected less treated habitat in which to settle. We monitored vireo nests at three general site types: (1) within the flood channel where exotic and native vegetation removal has occurred regularly (Channel), (2) three sites next to the flood channel where limited exotic and native vegetation removal has occurred (Off-channel), and (3) three sites that have been actively restored by planting native vegetation (Restoration). Nesting activity was monitored in 85 territories, 8 of which were occupied by single males. Of the completed nests, 39 percent were successful, and nest success did not differ among the three sites. Clutch size was greater in the Channel than the Off-channel sites, and the proportion of hatchlings that fledged was greater in Off-channel sites than Channel and Restoration sites. There were no other nest-level differences detected among site types, nor were there any differences in territory-level measures of productivity (young fledged per pair, double-brooding) among the sites. Overall, breeding success and productivity were slightly lower in 2021 than in 2020, with 66 percent of pairs fledgling at least one young and pairs fledging an average of 1.9±1.7 young. To investigate if the cumulative years of treatment had an impact on vireo reproductive effort, we looked at the effects of the Treatment Index on reproductive parameters. Results from generalized linear models indicated that treatment did not have an effect on vireo nesting effort or the number of vireo fledglings per pair produced in 2021. Similarly, our analysis of nest survival for 2021 revealed no effect of Treatment Index on daily survival rate. Analysis of vegetation data collected at vireo nests from 2006 to 2021 did not indicate an effect of vegetation at the nest on daily survival rate. We also found no differences in nest-placement characteristics among site types or successful/unsuccessful nests. Red/arroyo willow ( Salix laevigata or Salix lasiolepis ) was the species most commonly selected for nesting by vireos in all three site types. Black willow ( Salix gooddingii ) and mule fat ( Baccharis salicifolia ) also were commonly used. Vireos used a wider variety of species for nesting in Channel and Off-channel sites (seven and eight species, respectively) compared with Restoration sites (three species).

California↗

Data summary report: Unregulated contaminants monitoring project

The Drinking Water Protection Section of the Minnesota Department of Health conducted reconnaissance monitoring of selected public water systems in Minnesota. Funding was obtained primarily from the Environment and Natural Resources Trust Fund. Sampling was conducted in 2019 and 2021. Laboratory analysis of samples was conducted for a variety of different contaminants of emerging concern (CECs), including selected pharmaceuticals, pesticides, PFAS, wastewater indicators and other parameters chosen for the physical and land use setting surrounding the sampling points. Sampling site and parameter selection were designed with several goals, as follows: Characterize occurrence and distribution of selected CECs in settings where such chemicals are most likely to be present; Determine if any such occurrences represent a public health concern; Compare results from coupled source water and finished (i.e., treated) water samples at public water system sites where such sampling is feasible; Assess if results from geologically vulnerable (sensitive subject to rapid recharge) and geologically non-vulnerable settings differ significantly. 306 samples were collected as part of the study, from three networks of public water systems differentiated on the basis of source water type (i.e., surface water or groundwater) and land use environment (agricultural and wastewater influenced). This report provides a preliminary, qualitative evaluation of the results. Additionally, more rigorous research will be conducted on these water quality data to evaluate the below findings in more detail. High-level findings from this assessment include the following: Very few samples exceeded health-based guidance for CECs; o When this occurred, MDH staff conducted follow up sampling at the system and provided technical advice about managing the situation. Only a fraction of the CECs analyzed were detected; o Of the 522 different CECs analyzed in the water samples, 161 were detected in one or more samples; o Additionally, most detections were at low levels; Among the CEC classes included in the analytical work, pesticides and PFAS were generally detected at a greater frequency than other CECs; o See Executive Summary Figure 1. The ten most commonly detected individual compounds include: o Tribromomethane, or bromoform, (a disinfection by-product) (70% of sites where analyzed); o norgestrel (a pharmaceutical) (69% of sites where analyzed); o lithium (68% of sites where analyzed); o Metolachlor SA (52%), Deethylatrazine (49%), atrazine (45%), and deisopropylatrazine (31%) (pesticides); o PFBA (44%) and PFHxS (27%) (PFAS compounds); and o 5-methyl benzotriazole (29%) (a benzotriazole). Some CECs were detected more frequently in samples collected from surface waters than those collected from groundwater sources; CEC concentrations were generally higher in vulnerable settings compared to nonvulnerable settings; Whether CECs were detected more frequently in the source water or finished water varied by CEC class. For example, o Benzotriazoles and pharmaceuticals were more frequently detected in source water samples than finished water samples; and o Tribromomethane, or bromoform, a common disinfection by-product, was more frequently found in finished water samples than in source water samples. This work prompted a series of programmatic changes and innovations: A response framework was established for helping the program and public water systems manage detections of unregulated CECs in drinking water; Results were forwarded to the program within MDH responsible for developing healthbased guidance in order to nominate specific compounds found in drinking water but for which limited or no risk advice is available; MDH is seeking support from the Clean Water Council to support the establishment of permanent capacity within the Drinking Water Protection Section to continue sampling efforts of this type.

Minnesota↗

Ungulate migration in a changing climate—An initial assessment of climate impacts, management priorities, and science needs

Executive Summary Migratory behavior among ungulates in the Western United States occurs in response to changing forage quality and quantity, weather patterns, and predation risk. As snow melts and vegetation green-up begins in late spring and early summer, many migratory ungulates leave their winter range and move to higher elevation summer ranges to access high-quality forage and areas with vegetative cover for protection during fawning. Ungulates remain on these ranges until the fall when increasing snowfall and decreasing temperatures trigger them to migrate back to their lower elevation winter ranges. While researchers have begun to assess the effects of physical barriers such as roads and energy infrastructure on migration, less attention has been paid to understanding how changing climate conditions might affect ungulate movements and range habitats. Does earlier spring green-up make ungulates leave their winter ranges sooner? Do persistent drought conditions reduce the carrying capacity of seasonal range habitats or lead to shifts in migration pathways? These and other questions remain largely unanswered but could have cascading effects on ungulate population dynamics and migratory behavior. In February 2018, the Secretary of the Interior signed Department of the Interior Secretarial Order 3362 (SO3362), “Improving Habitat Quality in Western Big-Game Winter Range and Migration Corridors.” The order, which focuses on elk, mule deer, and pronghorn in 11 Western States, directs the Bureau of Land Management (BLM), the U.S. Fish and Wildlife Service (FWS), the National Park Service (NPS), and the U.S. Geological Survey (USGS) to partner with State wildlife agencies on their priorities and objectives for identifying and conserving ungulate migration corridors and winter-range habitat. The USGS Climate Adaptation Science Centers (CASCs) were established to help managers of the Nation’s fish, wildlife, waters, and lands understand the effects of climate change and adapt to changing conditions. To support the recent Department of the Interior (DOI) emphasis on ungulate migration corridors and winter-range habitat, this report assesses current information on how climate change could affect elk, mule deer, and pronghorn migration. The report synthesizes the drivers of migration, outlines what is known about how climate change might affect these drivers, and summarizes management priorities and science needs related to ungulate migration corridors and range habitat. A review of the literature on ungulate migration shows that the core drivers of spring migration are the timing of spring green-up and snowmelt, and the core driver of fall migration is winter severity. After exploring what is known about how these drivers affect or could be affected by climate change, several pathways through which ungulate migration could be altered were identified: (1) ungulates alter migration timing to better track plant phenology or in response to changes in winter conditions; (2) ungulates change their migration route or distance traveled during migration to accommodate changes in environmental conditions; and (3) ungulate populations that are currently migratory may begin to demonstrate interannual variability in whether they migrate, depending on environmental conditions and density-dependence, and may remain resident for sets of consecutive years. Through discussions with managers, physical barriers to movement such as roads and fences were identified as a core concern. In addition, the primary research needs of States are the acquisition and analysis of data on ungulate movements, to refine delineation of winter range, summer range, and corridors, and to support a better understanding of how ungulates use these habitats. When it comes to understanding climate effects, managers were more concerned with understanding the vulnerability of winter- and summer-range habitats than the vulnerability of migration corridors because of the influence of summer and winter forage on ungulate condition and reproductive success. Managers were also concerned about how forage quality and quantity might change because of stressors such as drought, wildfire, and invasive species and how they might need to alter habitat-treatment strategies as a result. More baseline data are needed before effective projections of ungulate migration, at a West-wide scale under climate change, can be made. These data needs include (1) more clearly defined corridors and seasonal range habitats; (2) a comprehensive understanding of the ecological drivers of migration across ungulate species and populations; and (3) the identification of environmental thresholds for key variables that influence migration, above which ungulates alter migratory behavior. The CASCs have several opportunities to play a role in addressing these needs. The CASCs could initiate projects to identify past and potential future changes and trends in key variables known to affect ungulate migration, such as plant phenology, forage quality, or winter severity. However, it would be difficult to use this information to determine what those trends mean for ungulate migration due to the lack of knowledge about environmental thresholds for ungulates. Additional projects would be required to compare multiple years of movement data with key variables to define thresholds. Once available, information on environmental thresholds could be integrated with projections of key variables to forecast the likelihood that the migration routes or the distance traveled could change—another area in which the CASCs could contribute. A more immediate role for the CASCs would be to carry out synthesis projects. One such project could summarize the “state of the science” on the drivers of ungulate migration. Although there are dozens of population- and location-specific studies on this topic, collating this information could help highlight trends in migration drivers that span species and geographies: a necessary first step toward determining the extent to which migration drivers could be affected by climate change. A second project could focus on what is known about how climate variability and change affect ungulate life-histories, population dynamics, and migration in the Western United States. The goal of this effort could be to identify knowledge clusters and information gaps that require further investigation. Together, these synthesized products could focus future scientific activities on the most pressing issues of ungulate migration and climate change in the Western United States.

Arizona, California, Colorado, Idaho, Montana, New↗

Statistics for stochastic modeling of volume reduction, hydrograph extension, and water-quality treatment by structural stormwater runoff best management practices (BMPs)

The U.S. Geological Survey (USGS) developed the Stochastic Empirical Loading and Dilution Model (SELDM) in cooperation with the Federal Highway Administration (FHWA) to indicate the risk for stormwater concentrations, flows, and loads to be above user-selected water-quality goals and the potential effectiveness of mitigation measures to reduce such risks. SELDM models the potential effect of mitigation measures by using Monte Carlo methods with statistics that approximate the net effects of structural and nonstructural best management practices (BMPs). In this report, structural BMPs are defined as the components of the drainage pathway between the source of runoff and a stormwater discharge location that affect the volume, timing, or quality of runoff. SELDM uses a simple stochastic statistical model of BMP performance to develop planning-level estimates of runoff-event characteristics. This statistical approach can be used to represent a single BMP or an assemblage of BMPs. The SELDM BMP-treatment module has provisions for stochastic modeling of three stormwater treatments: volume reduction, hydrograph extension, and water-quality treatment. In SELDM, these three treatment variables are modeled by using the trapezoidal distribution and the rank correlation with the associated highway-runoff variables. This report describes methods for calculating the trapezoidal-distribution statistics and rank correlation coefficients for stochastic modeling of volume reduction, hydrograph extension, and water-quality treatment by structural stormwater BMPs and provides the calculated values for these variables. This report also provides robust methods for estimating the minimum irreducible concentration (MIC), which is the lowest expected effluent concentration from a particular BMP site or a class of BMPs. These statistics are different from the statistics commonly used to characterize or compare BMPs. They are designed to provide a stochastic transfer function to approximate the quantity, duration, and quality of BMP effluent given the associated inflow values for a population of storm events. A database application and several spreadsheet tools are included in the digital media accompanying this report for further documentation of methods and for future use. In this study, analyses were done with data extracted from a modified copy of the January 2012 version of International Stormwater Best Management Practices Database, designated herein as the January 2012a version. Statistics for volume reduction, hydrograph extension, and water-quality treatment were developed with selected data. Sufficient data were available to estimate statistics for 5 to 10 BMP categories by using data from 40 to more than 165 monitoring sites. Water-quality treatment statistics were developed for 13 runoff-quality constituents commonly measured in highway and urban runoff studies including turbidity, sediment and solids; nutrients; total metals; organic carbon; and fecal coliforms. The medians of the best-fit statistics for each category were selected to construct generalized cumulative distribution functions for the three treatment variables. For volume reduction and hydrograph extension, interpretation of available data indicates that selection of a Spearman’s rho value that is the average of the median and maximum values for the BMP category may help generate realistic simulation results in SELDM. The median rho value may be selected to help generate realistic simulation results for water-quality treatment variables. MIC statistics were developed for 12 runoff-quality constituents commonly measured in highway and urban runoff studies by using data from 11 BMP categories and more than 167 monitoring sites. Four statistical techniques were applied for estimating MIC values with monitoring data from each site. These techniques produce a range of lower-bound estimates for each site. Four MIC estimators are proposed as alternatives for selecting a value from among the estimates from multiple sites. Correlation analysis indicates that the MIC estimates from multiple sites were weakly correlated with the geometric mean of inflow values, which indicates that there may be a qualitative or semiquantitative link between the inflow quality and the MIC. Correlations probably are weak because the MIC is influenced by the inflow water quality and the capability of each individual BMP site to reduce inflow concentrations.

Scientific Investigations Report↗

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↗

Geohydrology of the Cross-Florida Barge Canal area, with special reference to the Ocala vicinity

The Cross-Florida Barge Canal route commences at Palatka on the St. Johns River, about 75 miles upstream from the Atlantic Ocean, and extends 110 miles southwestward across Peninsular Florida into deep water in the Gulf of Mexico near Yankeetown. The canal will be equipped with five locks, each 600 feet long and 84 feet wide, and the channel will be a minimum of 12 feet deep and 150 feet wide. From near Ocala northeastward, the canal channel will replace much of the natural channel of the Oklawaha River, and will be excavated into beds of the so-called shallow sand aquifer of Miocene age and younger, which overlies limestone of the Floridan aquifer. Westward from Ocala, most of the canal will be excavated below the potentiometric surface into limestone and dolomite of the Floridan aquifer. Water levels of Rodman, Eureka, and Inglis Pools will be controlled by dams and spillways with the limited exchange of water between the pools and the aquifers. The water levels in the Summit Pools will fluctuate with the natural changes in the ground-water level of the Floridan aquifer, although the stage of the pool will be controlled partly by the stage held in the Eureka Pool. A dynamic inflow-outflow relationship will exist between the Summit Pool and the Floridan aquifer. The Floridan aquifer in the canal area is 1,000 to 1,200 feet thick and consists of limestone and dolomite of middle Eocene Miocene age, including from older to younger, the Lake City, Avon Park, and Ocala limestones plus permeable sandy, dolomitic limestone in the lower part of the Hawthorn Formation. It is possible that most of the flow to the two major springs in the area occurs in the upper 100 feet or so of the aquifer in the Ocala Limestone. The aquifer is underlain by the Oldsmar limestone of early Eocene age and is overlain by sand, clayey sand, clay and shell beds of Miocene through Holocene age, in thickness from a few feet to 300 feet. The permeable beds overlying the Floridan aquifer constitute the shallow aquifer, while the poorly permeable ones act as confining beds where the Floridan aquifer is under artesian conditions. A north-south line drawn separating the head of Silver Springs on the west from the Oklawaha River on the east marks the approximate western limit of a continuous blanket of materials of Miocene-Pliocene(?) age covering the rocks of the Floridan aquifer. East of the line, much of the aquifer is under artesian conditions, particularly in the Oklawaha River valley, although in some areas east of the valley, direct recharge through thick permeable Miocene-Pliocene(?) sands occurs. West of the line, only scattered remnants of a once continuous Miocene-Pliocene(?) cover remain. Lack of the cover is a result of erosion on the crest and flank of the Ocala Uplift, a broad northwest-southeast trending anticlinal upwarp, the axis of which is crossed by the canal route in the Dunnellon area. Over most of this area the Floridan aquifer is unconfined and receives direct recharge through a cover of a few tens of feet of sand and clayey sand of Quaternary age. Tensional stresses during the structural evolution of the Ocala Uplift produced an intersecting system of fractures and normal faults in rocks of the Floridan aquifer. The fractures and faults are important controls for orientation of solution channels and, therefore, for development of ground-water circulation patterns. When the system surface streams, which once drained the Barge Canal area, eroded the poorly permeable Miocene-Pliocene(?) cover from the flanks of the Ocala Uplift, surface runoff was reduced and precipitation began to directly infiltrate the underlying limestones. Now only principal streams remain, such as the Oklawaha and Withlacoochee Rivers and a few short tributaries, while one of the most highly developed subsurface drainage systems in the world has evolved in cavernous limestone of the Floridan aquifer. Two of the larger freshwater springs in the world now discharge from the Floridan aquifer in the canal area. Silver Springs near Ocala discharges an average 531 mgd (million gallons per day) down the 4-mile long Silver River, which flows on poorly permeable beds to the Oklawaha River. Rainbow Springs near Dunnellon discharges on average 468 mgd from numerous orifices in the bed of the 5-mile-long Rainbow River, which flows into the Withlacoochee River. The heads of the springs have migrated to their present positions partly because of a tendency of ground-water levels to decline as permeability in the aquifer is increased due to removal of limestone by solution, and because of mechanical erosion of the limestone in the vicinity of the spring heads. Also, points of principal spring discharges have shifted in the past due to changes in ground-water levels in response to changes in sea level. The subsurface drainage system is continuing to evolve today, as evidenced in part by frequent occurrence of new sinkholes and by the presence of significant amounts of calcium bicarbonate in the spring waters. Rodman Pool, at the east end of the canal, is separated from the Floridan aquifer by poorly permeable materials. The pool's operating water level will be only a few feet above the potentiometric surface at the downstream end, and at or slightly below the potentiometric surface at the upstream end. Little exchange of water between the Rodman Pool and the Floridan aquifer is expected. Eureka Pool, just upstream from Rodman Pool, will also be separated from the Floridan aquifer by poorly permeable beds. However, the stage of the pool will be about 15 feet higher than the natural potentiometric surface at the pool's downstream end, and some seepage into the Floridan aquifer is anticipated through faults and leaky parts of the poorly permeable beds, with a consequent rise in ground-water levels in areas adjacent to the lower end of the pool. Possibilities for particulate contamination of the aquifer will tend to be minimized because of the filtering capacity of the materials through which water must pass to reach the aquifer, although the natural filter will not preclude movement into the aquifer of contaminants which might become dissolved in the pool waters. No significant interchange of water between the pool and the aquifer is expected at the upstream end of the Eureka Pool. Present construction plans indicate an operating stage for Eureka Pool which will range between 38 and 40 feet above mean sea level, although it is possible to dredge the pool deep enough to permit a range in stage of 36 to 40 feet. A backwater effect extending up Silver River from the Eureka Pool is expected to regulate the stage at the head of Silver Springs between 39 and 44 feet above mean sea level if the pool ranges between 36 and 40 feet. If Eureka Pool ranges only between 38 and 40 feet, the range of stage at the head of the springs should be about 41 to 44 feet above mean sea level. From Inglis Lock west, the canal will have direct connection with the Gulf, and canal stage will fluctuate with the Gulf tide. Since the canal stage will be slightly lower than the adjacent ground-water levels along much of the reach, there will be some ground-water inflow to the canal. No significant changes in the existing ground-water regime are expected in the vicinity of Inglis Pool, the first step up in the canal east of the Gulf. Existing ground-water and surface-water levels in the area will not change appreciable, and the natural stage and flow of Rainbow Springs, which will flow by way of Rainbow River into Inglis Pool, should not be affected by canal operations. A possible adverse effect of the canal on the Inglis Pool area could result if sea water is locked up from the Gulf through Inglis Lock. However, the high step of 25 feet at the lock, flushing action of continuous flow from Inglis Pool to the lower reaches of the Withlacoochee River, and use of possible preventive locking procedures should minimize the problem. The potential for adverse effects on the ground-water regime is greatest in the area of the Summit Pool. Through most of the length of the pool, the canal channel will be excavated into limestone of the Floridan aquifer to depths of 12 to 27 feet below the potentiometric surface. Changes that will take place in the ground-water flow system in the Silver Springs drainage area, once the canal is completed, were estimated by flow-net analysis. Variation in aquifer transmissivity was determined by calculating transmissivity in 25 different flow cells surrounding Silver Springs. Transmissivity in the 25 cells averages about 15,600,000 gpd/ft [2,090,000 ft2/day (feet squared)], but transmissivity in the six cells through which the Summit Pool passes ranges from 9,000,000 to 44,000,000 gpd/ft (1,210,000 to 5,900,000 ft2/day). Transmissivity was used to compute static stage of the Summit Pool under given ground-water level conditions. Had the canal existed in May 1968 and had the stage of Eureka Pool been held at 36 feet at the time, the static stage in Summit Pool would have been about 42.1 feet above mean sea level. Thus, a conceptual model of the changes in the potentiometric surface wrought by the finished canal was drawn, and zones of ground-water inflow and outflow were delineated. Most outflow from the Summit Pool to the aquifer should be limited to one 4-mile-long zone along the north side of the pool, about 5 miles south of Silver Springs. It is estimated that a water volume equivalent to about 8 percent of the daily flow of Silver Springs will enter the Summit Pool each day from the southern one-third of the Silver Springs drainage area. A like amount will reenter the aquifer at the main zone of outflow and move toward Silver Springs at an estimated average velocity of about 200 feet per day, if something close to the natural static stage of the pool is maintained by return pumpage of the lockage losses. At a velocity of 200 feet per day, water from the Summit Pool would discharge at Silver Springs about 140 days later. However, any estimate of velocity in the highly cavernous limestone aquifer in the area should be used with caution, because difficult to measure changes in porosity and thickness of major zones of flow may cause large variations in velocity. If all lockage losses are returned to the Summit Pool by pumping from Eureka Pool, no net loss from the Silver Springs drainage area, except for some evaporation from the water surface in the canal and possible leakage around locks, will result from canal operations. The zone of outflow from the Summit Pool to the aquifer will be in a natural potentiometric trough, and the zone of inflow will be in a potentiometric ridge area. The equilibrium water level in the pool will tend to be about 1 foot higher than the altitude of the lowest level in the pre-canal potentiometric trough, and about 2 feet lower than the highest level on the pre-canal potentiometric ridge. West of the Silver Springs drainage area just east of Dunnellon Lock, in the area of a local potentiometric high, the water level in the Summit Pool is expected to be about 15 feet below the natural potentiometric surface. In most areas 2 to 3 miles away from the Summit Pool, effects of the canal on the natural potentiometric surface should be slight. The stage of the Summit Pool, judging from the 36-year record for ground-water level changes and the anticipated indirect effect of the controlled stage in Eureka Pool, should have a maximum of about 10.5 feet with a maximum water level of about 51.5 feet above mean sea level and a minimum of about 41.0 feet above mean sea level. Of particular importance in the Summit Pool is an implementation of well planned construction and operational procedures designed to minimize risks of ground-water contamination.

Florida↗

Simulation of groundwater flow and streamflow depletion in the Branch Brook, Merriland River, and parts of the Mousam River watersheds in southern Maine

Watersheds of three streams, the Mousam River, Branch Brook, and Merriland River in southeastern Maine were investigated from 2010 through 2013 under a cooperative project between the U.S. Geological Survey and the Maine Geological Survey. The Branch Brook watershed previously had been deemed “at risk” by the Maine Geological Survey because of the proportionally large water withdrawals compared to estimates of the in-stream flow requirements for habitat protection. The primary groundwater withdrawals in the study area include a water-supply well in the headwaters of the system and three water-supply wells in the coastal plain near the downstream end of the system. A steady-state groundwater flow model was used to understand the movement of water within the system, to evaluate the water budget and the effect of groundwater withdrawals on streamflows, and to understand streamflow depletion in relation to the State of Maine’s requirements to maintain in-stream flows for habitat protection. Delineation of the simulated groundwater divides compared to the surface-water divides suggests that the groundwater divides in the headwater areas do not exactly correspond to the surface-water divides. Under both pumping and non-pumping conditions, groundwater flows from the headwaters of the Branch Brook watershed into the Mousam River watershed. Pumping in the Mousam River watershed captures a small amount of groundwater from the Branch Brook basin. The cumulative effect of groundwater withdrawals on base flows in two rivers in the study area (Branch Brook and the Merriland River) was evaluated using the groundwater flow model. Streamflow depletion in the headwaters of Branch Brook was 0.12 cubic feet per second (ft 3 /s) for the steady-state simulation, or about 10 percent of the average base flow at that location. Downstream on Branch Brook, the total streamflow depletion from all the wells was 0.59 ft 3 /s, or 3 percent of the average base flow at that location. In the Merriland River downstream from the Merriland River well, the total amount of streamflow depletion was 0.6 ft 3 /s, or about 7 percent of the average base flow. The groundwater model was used to evaluate several different scenarios that could affect streamflow and groundwater discharging to the rivers and streams in the study area. The scenarios were (1) no pumping from the water-supply wells; (2) current pumping from the water-supply wells, but simulated drought conditions (25 percent reduction in recharge); (3) current recharge, but with increased pumping from the large water-supply wells; and (4) drought conditions and increased pumping combined. Simulations of increased pumping in the water-supply wells resulted in streamflow depletion in the headwaters of Branch Brook increasing to 16 percent of the headwater base flow. Simulated increases in the pumping in the coastal plain wells increased the amount of streamflow depletion to 6 percent of the flow in Branch Brook and to 8 percent of the flow in the Merriland River. The additional stress of a drought imposed on the model (25 percent less recharge) had a substantial impact on streamflows, as expected. If the simulated drought occurred simultaneously with an increase in pumping, the base flows would be reduced 48 percent in the headwaters of Branch Brook, compared to the no-pumping scenario. Downstream in Branch Brook, the total reduction in flow would be 29 percent of the simulated base flows in the no-pumping scenario, and in the Merriland River, the reduction would be 33 percent of the base flows in the no-pumping scenario. The study evaluated two different methods of calculating in-stream flow requirements for Branch Brook and the Merriland River—a set of statewide equations used to calculate monthly median flows and the MOVE.1 record-extension technique used on site-specific streamflow measurements. The August median in-stream flow requirement in the Merriland River was calculated as 7.18 ft 3 /s using the statewide equations but was 3.07 ft 3 /s using the MOVE.1 analysis. In Branch Brook, the August median in-stream flow requirements were calculated as 20.3 ft 3 /s using the statewide equations and 11.8 ft 3 /s using the MOVE.1 analysis. In each case, using site-specific data yields an estimate of in-stream flow that is much lower than an estimate the statewide equations provide.

Maine↗

Simulation of groundwater conditions and streamflow depletion to evaluate water availability in a Freeport, Maine, watershed

In order to evaluate water availability in the State of Maine, the U.S. Geological Survey (USGS) and the Maine Geological Survey began a cooperative investigation to provide the first rigorous evaluation of watersheds deemed "at risk" because of the combination of instream flow requirements and proportionally large water withdrawals. The study area for this investigation includes the Harvey and Merrill Brook watersheds and the Freeport aquifer in the towns of Freeport, Pownal, and Yarmouth, Maine. A numerical groundwater- flow model was used to evaluate groundwater withdrawals, groundwater-surface-water interactions, and the effect of water-management practices on streamflow. The water budget illustrates the effect that groundwater withdrawals have on streamflow and the movement of water within the system. Streamflow measurements were made following standard USGS techniques, from May through September 2009 at one site in the Merrill Brook watershed and four sites in the Harvey Brook watershed. A record-extension technique was applied to estimate long-term monthly streamflows at each of the five sites. The conceptual model of the groundwater system consists of a deep, confined aquifer (the Freeport aquifer) in a buried valley that trends through the middle of the study area, covered by a discontinuous confining unit, and topped by a thin upper saturated zone that is a mixture of sandy units, till, and weathered clay. Harvey and Merrill Brooks flow southward through the study area, and receive groundwater discharge from the upper saturated zone and from the deep aquifer through previously unknown discontinuities in the confining unit. The Freeport aquifer gets most of its recharge from local seepage around the edges of the confining unit, the remainder is received as inflow from the north within the buried valley. Groundwater withdrawals from the Freeport aquifer in the study area were obtained from the local water utility and estimated for other categories. Overall, the public-supply withdrawals (105.5 million gallons per year (Mgal/yr)) were much greater than those for any other category, being almost 7 times greater than all domestic well withdrawals (15.3 Mgal/yr). Industrial withdrawals in the study area (2.0 Mgal/yr) are mostly by a company that withdraws from an aquifer at the edge of the Merrill Brook watershed. Commercial withdrawals are very small (1.0 Mgal/yr), and no irrigation or other agricultural withdrawals were identified in this study area. A three-dimensional, steady-state groundwater-flow model was developed to evaluate stream-aquifer interactions and streamflow depletion from pumping, to help refine the conceptual model, and to predict changes in streamflow resulting from changes in pumping and recharge. Groundwater levels and flow in the Freeport aquifer study area were simulated with the three-dimensional, finite-difference groundwater-flow modeling code, MODFLOW-2005. Study area hydrology was simulated with a 3-layer model, under steady-state conditions. The groundwater model was used to evaluate changes that could occur in the water budgets of three parts of the local hydrologic system (the Harvey Brook watershed, the Merrill Brook watershed, and the buried aquifer from which pumping occurs) under several different climatic and pumping scenarios. The scenarios were (1) no pumping well withdrawals; (2) current (2009) pumping, but simulated drought conditions (20-percent reduction in recharge); (3) current (2009) recharge, but a 50-percent increase in pumping well withdrawals for public supply; and (4) drought conditions and increased pumping combined. In simulated drought situations, the overall recharge to the buried valley is about 15 percent less and the total amount of streamflow in the model area is reduced by about 19 percent. Without pumping, infiltration to the buried valley aquifer around the confining unit decreased by a small amount (0.05 million gallons per day (Mgal/d)), and discharge to the streams increased by about 8 percent (0.3 Mgal/d). A 50-percent increase in pumping resulted in a simulated decrease in streamflow discharge of about 4 percent (0.14 Mgal/d). Streamflow depletion in Harvey Brook was evaluated by use of the numerical groundwater-flow model and an analytical model. The analytical model estimated negligible depletion from Harvey Brook under current (2009) pumping conditions, whereas the numerical model estimated that flow to Harvey Brook decreased 0.38 cubic feet per second (ft 3 /s) because of the pumping well withdrawals. A sensitivity analysis of the analytical model method showed that conducting a cursory evaluation using an analytical model of streamflow depletion using available information may result in a very wide range in results, depending on how well the hydraulic conductivity variables and aquifer geometry of the system are known, and how well the aquifer fits the assumptions of the model. Using the analytical model to evaluate the streamflow depletion with an incomplete understanding of the hydrologic system gave results that seem unlikely to reflect actual streamflow depletion in the Freeport aquifer study area. In contrast, the groundwater-flow model was a more robust method of evaluating the amount of streamflow depletion that results from withdrawals in the Freeport aquifer, and could be used to evaluate streamflow depletion in both streams. Simulations of streamflow without pumping for each measurement site were compared to the calibratedmodel streamflow (with pumping), the difference in the total being streamflow depletion. Simulations without pumping resulted in a simulated increase in the steady-state flow rate of 0.38 ft 3 /s in Harvey Brook and 0.01 ft 3 /s in Merrill Brook. This translates into a streamflow-depletion amount equal to about 8.5 percent of the steady-state base flow in Harvey Brook, and an unmeasurable amount of depletion in Merrill Brook. If pumping was increased by 50 percent and recharge reduced by 20 percent, the amount of streamflow depletion in Harvey Brook could reach 1.41 ft 3 /s.

Maine↗

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↗

Detailed study of selenium and selected constituents in water, bottom sediment, soil, and biota associated with irrigation drainage in the San Juan River area, New Mexico, 1991-95

In response to increasing concern about the quality of irrigation drainage and its potential effects on fish, wildlife, and human health, the U.S. Department of the Interior began the National Irrigation Water Quality Program (NIWQP) to investigate these concerns at irrigation projects sponsored by the Department. The San Juan River area in northwestern New Mexico was one of the areas designated for study. Study teams composed of scientists from the U.S. Geological Survey, the U.S. Fish and Wildlife Service, the Bureau of Reclamation, and the Bureau of Indian Affairs collected water, bottom-sediment, soil, and biological samples at 61 sites in the San Juan River area during 1993-94. Supplemental data collection conducted during 1991-95 by the Bureau of Indian Affairs and its contractor extended the time period and sampling sites available for analysis. Analytical chemistry performed on samples indicated that most potentially toxic elements other than selenium generally were not high enough to be of concern to fish, wildlife, and human health. Element concentrations in some water, bottom-sediment, soil, and biological samples exceeded applicable standards and criteria suggested by researchers in current literature. Selenium concentrations in water samples from 28 sites in the study area exceeded the 2-microgramper-liter (lg/L) wildlife-habitat standard. Vanadium concentrations in water exceeded the 100-Kg/L standard for livestock-drinking water at one site. In biota, selenium and aluminum concentrations regularly equaled or exceeded avian dietary threshold concentrations. In bottom sediment and soil, element concentrations above the upper limit of the baseline range for western soils were: selenium, 24 exceedances; lead, 2 exceedances; molybdenum, 2 exceedances;strontium, 4 exceedances; and zinc, 4 exceedances. Concentrations of total selenium in bottom-sediment and soil samples were significantly greater for Cretaceous than for non-Cretaceous soil types in the study area and were generally similar for habitats within and outside irrigation-affected areas. Mean and median total-selenium concentrations in samples from areas with Cretaceous soil types were 4.6 and 2.2 micrograms per gram (ps/g), respectively. Mean and median total-selenium concentrations in samples from areas with non-Cretaceous soil types were 0.6 and 0.15 pg/g, respectively. Samples from the study area had low concentrations of organic constituents. Organochlorine pesticides and polychlorinated biphenyls were detected in a few biological samples at low concentrations. Polycyclic aromatic hydrocarbon (PAH) compounds were not detected in whole-water samples collected using conventional water-sampling techniques. In tests involving the use of semipermeable-membrane devices to supplement conventional water assays for PAH's, low concentrations of PAH's were found at several locations in the Hammond Irrigation Supply Canal, but were not detected in the Hammond ponds at the downstream reach of the Hammond irrigation service area. PAH compounds do not appear to reach the San Juan River through the Hammond Canal. Data indicate that water samples from irrigation-drainage-affected habitats had increased mean selenium concentrations compared with samples from irrigation-delivery habitat. The mean selenium concentration in water was greatest at seeps and tributaries draining irrigated land (17 μg/L); less in irrigation drains and in ponds on irrigated land (61.tg/L); and least in backwater, the San Juan River, and irrigation-supply water (0.5 - 0.6 μg/L). Statistical tests imply that irrigation significantly increases selenium concentrations in water samples when a Department of the Interior irrigation project is developed on selenium-rich sediments. Water samples from sites with Cretaceous soils had significantly greater selenium concentrations than water samples from sites with non-Cretaceous soils. Water samples from Department of the Interior project irrigation-drainage sites developed on Cretaceous soils contained a mean selenium concentration about 10 times greater than those in samples from Department of the Interior project sites developed on non-Cretaceous soils. Selenium was much less concentrated in water than in bottom sediment, soil, or biota in the study area. The range in concentrations of dissolved selenium in water was less than 1 ptg/L to 37 1.1g/L (less than 1 to 37 parts per billion). The range in concentrations of total selenium in bottom sediment and soil was less than 0.1 to 23lig/g (less than 100 to 23,000 parts per billion). The range in concentration of selenium in biota was less than 0.1 to 24.0 fig/g (less than 100 to 24,000 parts per billion). Data indicated that bioaccumulation and leaching from soil were the important processes at the study area that lead to elevated levels of selenium. Other processes examined included: (1) evapoconcentration of selenium; (2) atmospheric deposition of aerosols containing selenium; and (3) contamination of surface water by point-source or non-point-source discharges. Selenium concentrations in biological samples were evaluated by a number of variables including: (1) media sampled (emergent and submergent plants, nektonic and benthic invertebrates, omnivore/herbivore and carnivore fish, and terrestrial and aquatic amphibians); (2) habitat (San Juan River main-stem reaches, backwaters, tributary reaches, irrigation delivery or drainage canals, and ponds); (3) irrigation project area and reference sites; and (4) soil type (non-Cretaceous or Cretaceous soils). Graphical techniques and nonparametric statistical tests were applied to determine the influence of selected physiographic variables on selenium concentrations in biological samples collected in the San Juan River area. Species of sucker and of smaller fish contained significantly higher selenium concentrations in the upstream portion of the river where a productive community of plants and animals is found that is associated with warming, nutrient-rich waters discharged from an upstream reservoir. Selenium concentrations in algae, odonates, and mosquitofish collected from both irrigation-drain and pond habitats underlain by Cretaceous soils were significantly greater than in those collected from similar habitats underlain by non-Cretaceous soils. Investigators conclude that the major factor affecting the variability of selenium accumulation in biota at aquatic habitats was the presence of underlying Cretaceous soils. Median selenium concentrations were less than 2 lAg/g for plant samples, less than 7 μg/g for invertebrate samples, and less than 6 lAg/g for whole-fish samples collected from aquatic habitats underlain by non-Cretaceous soils. Similar samples collected from aquatic habitats underlain by Cretaceous soils contained median selenium concentrations two to five times greater. Leaching of selenium from Cretaceous soils in the San Juan River area increases the accumulation of selenium concentrations in the biota and thereby increases the exposure and potential health risks associated with selenium to migratory birds, fish, and other wildlife that use these aquatic habitats extensively. Aquatic habitats presenting the greatest average exposure to excess selenium concentrations in the diets of resident wildlife are from consumption of plants, invertebrates, and fish at irrigation-drain habitats underlain by Cretaceous soils. Of the irrigation projects evaluated in the San Juan River area, the highest median selenium concentrations in algae, cattail leaves, odonate nymphs, mosquitofish, and leopard frog samples from the study area were collected from the east hogback irrigation drain.

New Mexico↗

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↗

Reservoir evolution, downstream sediment transport, downstream channel change, and synthesis of geomorphic responses of Fall Creek and Middle Fork Willamette River to water years 2012–18 streambed drawdowns at Fall Creek Lake, Oregon

Executive Summary Chapter A. Introduction Fall Creek Dam impounds Fall Creek Lake, a 10-kilometer-long reservoir in western Oregon and is operated by the U.S. Army Corps of Engineers (USACE) primarily for flood-risk management (or flood control) in late autumn through early spring months, as well as for water quality, irrigation, recreation, and habitat in late spring through early autumn. Since 2011 (water year [WY] 2012), Fall Creek Lake has been temporarily drawn down each year to facilitate downstream passage of juvenile spring Chinook salmon ( Oncorhynchus tshawytscha ) through the 55-meter (m) high dam. This annual dam operation is temporary, typically lasting about 1–2 weeks from WY 2012 through 2020 (drawdown operations in WY 2022–24 have increased to more than 6 weeks). Drawdown of the reservoir results in lake levels being lowered to the elevation near the historical, pre-dam streambed. The annual streambed drawdowns of WY 2012–18 have improved fish passage and led the USACE to formally adopt streambed drawdowns as part of annual operations at Fall Creek Dam. However, temporarily lowering the lake to streambed creates free-flowing conditions in the reservoir that result in the erosion and episodic export of predominantly sand and finer-grained sediments (less than 2 millimeters [mm]) to the lower gravel-bed reaches of Fall Creek and the Middle Fork Willamette River. The introduction of large volumes of sand and finer-grain sediment into the dam-regulated reaches downstream from Fall Creek Dam prompted questions about the geomorphic responses to annual streambed drawdowns within Fall Creek Lake and downstream reaches along Fall Creek and the Middle Fork Willamette River. The U.S. Geological Survey (USGS) in partnership with USACE initiated a comprehensive geomorphic and sediment transport investigation to assess the coupled processes of reservoir erosion, sediment evacuation from Fall Creek Lake, and patterns of sediment transport and deposition in reaches downstream from the Fall Creek Dam that have resulted from annual streambed drawdowns. The purpose of this report is to systematically describe the processes of sediment erosion, transport, and deposition at Fall Creek Lake and geomorphic interactions between reaches upstream and downstream from Fall Creek Dam that relate to dam operations. Specifically, this report focuses on evaluating geomorphic responses to streambed drawdowns from WY 2012 through 2018 and placing drawdown-induced geomorphic responses within the broader context of physiographic and historical conditions and dam operations of Fall Creek and Middle Fork Willamette Rivers. Key objectives for this study were to characterize changes in reservoir morphology and substrate at Fall Creek Lake, describe the character and temporal pattern of sediment transport downstream from Fall Creek Dam, characterize geomorphic changes in channel reaches downstream from the Fall Creek Dam, and relate these data to the annual streambed drawdowns of WY 2012–18. This study uses multiple independent monitoring and measurement approaches to assess site, reach, and river-scale geomorphic responses to drawdowns to inform dam and reservoir management. Patterns and processes of reservoir evolution were assessed with geomorphic mapping and volumetric analyses of topography through comparison of multiple digital surface models (DSMs). Just downstream from Fall Creek Dam, analyses of sediment export from the reservoir focused on suspended sediment but also incorporated bedload analyses to assess sediment sizes. Geomorphic assessments downstream from the dam used reach-scale and site-scale approaches to document changes in channel morphology and substrate, including site measurements of sand and finer-grained sediment deposition and in-channel bed-material, volumetric change analyses from comparison of digital elevation models (DEMs), and repeat geomorphic mapping. Findings from this study inform river management and dam operations by providing an understanding of (1) coupled upstream-downstream geomorphic responses to the Fall Creek Lake streambed drawdowns, (2) geomorphic responses of Fall Creek Lake streambed drawdowns in comparison to drawdowns at other large dams, (3) controls on reservoir erosion and downstream geomorphic responses, and (4) implications for future hydrogeomorphic changes that may result from continued drawdowns and monitoring activities to assess those changes. Chapter B. Reservoir Morphology and Evolution Related to Dam Operations at Fall Creek Lake To understand the volume and distribution of sediment accumulation in Fall Creek Lake since dam closure in 1965, decadal-scale sedimentation patterns (spanning approximately 1965–2016) are evaluated using a combination of storage curve analyses and geomorphic mapping. Short-term (drawdown event-scale) patterns of erosion, sedimentation, and sediment export downstream are evaluated using a combination of geomorphic mapping and change detection analyses that quantify the distribution and total volume of sediment erosion and deposition within Fall Creek Lake. Geomorphic mapping of reservoir topography and analyses of historical datasets reveals four categories of landforms and sediment processes within Fall Creek Lake related to lake level operations: lacustrine sedimentation expressed in the reservoir floor, fluvial erosion and deposition within historical stream channels during streambed drawdowns, channel-like features created by erosion within the reservoir floor during streambed drawdowns, and erosion on reservoir hillslopes. Where the reservoir floor is mapped for this study as pelagic (deep water), deposition up to 3 meters (m) thick by lacustrine processes and burial of pre-dam topography with deposits thinning toward the edges of the valley floor and upstream areas of reservoir are observed. Despite over 50 years of sediment accumulation since dam construction, the main stream channels of Fall and Winberry Creeks (or reservoir thalwegs) through the reservoir are well defined, though their distinct morphology is likely influenced by a long history of recurring historical drawdowns to or near streambed since dam construction. Unregulated streamflow and sediment transport through the reservoir primarily are confined to these channels during the streambed drawdown periods. Erosional channel-like features created by drawdowns are carved through underlying, unconsolidated reservoir floor sediments and are most prominent in the lower reservoir below minimum conservation pool (the low pool elevation during winter flood season); sediment generated from the formation of these drawdown channels is more likely to be transported through and out of the reservoir than sediment deposits along the reservoir hillslopes at the valley margins that are separated from main channels by areas of low-gradient reservoir floor. Morphologic changes in the lower reservoir topography between January 2012 and November 2016 indicate overall net erosion of about 129,500 cubic meters (m 3 ). The most prominent geomorphic changes occurred along the main channels of Fall and Winberry Creeks near the Fall Creek Dam where incision, lateral migration, and slumping banks resulted in vertical and lateral adjustments to channel position, whereas most changes fell below the detectable limit on higher-elevation reservoir floor surfaces except where erosion occurred along features mapped as drawdown channels. Chapter C. Sediment Delivery from Fall Creek Lake and Transport through Downstream Reaches USGS implemented a sediment monitoring program in WY 2013–18 to evaluate the quantity and character of reservoir sediment exported from Fall Creek Lake during streambed drawdowns. Turbidity and suspended sediments were monitored annually autumn through spring to span the WY 2013–18 streambed drawdowns; however, unequal monitoring timeframes each year reduced the ability to compare results and factors affecting sediment export from the reservoir difficult between years. These data were originally measured to develop regressions and compute suspended-sediment loads (SSL). Bedload sediment monitoring from a cableway at the Fall Creek streamgage was completed in the autumn-winter of WY 2013 and 2017. The limited number of samples and presumed variability in sediment supply from the reservoir precluded construction of streamflow and bedload discharge relations to compute more than instantaneous bedload. Sand and finer-grained silts and clays were transported from the reservoir in suspension, though some coarser grains (up to 32 mm) were also mobilized and transported downstream from the dam as bedload. Observations of increased sediment transport downstream from Fall Creek Dam coincided with lake levels approaching about 3 m (10 feet [ft] or elevation 690 ft) above the streambed regulating outlets. Suspended-sediment loads computed for the full monitoring periods WY 2013–18 at the Fall Creek streamgage, located 1.4 kilometers (km) downstream from Fall Creek Dam, range from 54,700 metric tons (t) in WY 2013 to 13,900 t in WY 2018. Although the total annual SSL varied from year to year, the overall seasonal patterns of suspended sediment transport throughout each year were similar during monitoring in WY 2013-18. Suspended-sediment loads were low prior to the drawdown, then increased rapidly as lake levels lowered and approached the streambed. In the weeks following the drawdown period, as pool levels were increased, SSL remained slightly elevated above pre-drawdown levels but generally declined through the following winter and spring except during streamflow-driven pulses of suspended-sediment transport. WY 2013 had the greatest total computed SSL for each streambed drawdown and partial-year monitoring period. SSL computed for the partial-year period have generally decreased since WY 2013 and have varied by about 6,800 t with the exception of WY 2014. WY 2014 SSL reflects anomalously low sediment export due to low streamflows and freezing conditions that stabilized reservoir floor deposits. Bedload measurements in the short 1.4-km reach between Fall Creek Dam and the Fall Creek streamgage showed an inverse correlation between bedload transport rates and discharge, which probably reflects diminishing supply of coarse-sized sediment. Sand was more abundant (60–100 percent) than gravel in bedload samples confirming sand and finer-grained sediment dominated sediment evacuated from the reservoir during streambed drawdowns at Fall Creek Lake. Chapter D. Geomorphic Responses to Fall Creek Lake Streambed Drawdowns Downstream from Fall Creek Dam In the days, weeks, and months following streambed drawdown operations at Fall Creek Dam through WY 2018, sites downstream from the dam displayed a variety of geomorphic responses to reservoir sediment delivery within the main channel and overbank areas. Evaluation of streambed elevations at two streamgages located 1.4 km downstream from the dam on Fall Creek and 16.3 km downstream from the dam on the Middle Fork Willamette River indicated the effects of drawdown sediment on bed elevations were modest and transient. Repeat particle size measurements (October 2015 and September 2016) at five sites along Fall Creek and the Middle Fork Willamette River showed similar grain-sized distributions that do not reveal substantial deposition of fine-grained sediment related to the WY 2016 streambed drawdown. Altogether, these findings indicate that transport capacity in the main, low-flow channels of Fall Creek and Middle Fork Willamette River during WY 2012–18 was sufficient to mobilize and evacuate reservoir sediments from streambed drawdowns or other bank material and tributary sources. However, other monitoring for this study indicate low-velocity zones in off-channel areas are prime locations for sand and finer-grain sediment deposition. Patterns of overbank sediment accumulation indicate that the magnitude and timing of overbank deposition on bars and low-elevation floodplain varies with proximity to the dam, geomorphic setting, streamflows, and other factors. Sand and finer-grained reservoir sediments carried as suspended-sediment load in the reaches downstream from Fall Creek Dam were deposited in overbank areas as observed with clay-horizon markers during WY 2016–17. Overbank deposition quantified with Geomorphic Change Detection (GCD) software evaluated landform-scale patterns of erosion and deposition using repeat light detection and ranging (lidar) surveys at two sites in the Upper Fall Creek reach and one site in the Jasper reach for 3 years (2012–15) and one site in the Clearwater reach for 6 years (2009–15). Deposition thickness and spatial patterns from the GCD analysis were variable; some sites had dispersed but measurable deposition while at others, deposition was highly localized and exceeded 1 m in depth. Patterns of overbank deposition illustrate interactions among bar morphology, local hydraulics, and suspended-sediment transport dynamics that can create patches of highly localized deposition. The measured deposition at the two Fall Creek GCD sites likely resulted from reservoir sediments released from Fall Creek Lake during streambed drawdowns in WY 2016 and 2017 because the limited sediment inputs from bank material (geomorphically laterally stable reach) or tributaries (no significant tributaries) provided few other sediment sources. On the Middle Fork Willamette River, observed patterns of overbank deposition could reflect sediment sourced from upstream tributaries, bank erosion, or Fall Creek Lake streambed drawdown operations. Despite the introduction of several thousand tons of reservoir sediment delivered from the Fall Creek Lake streambed drawdowns to below-dam river corridors, reach-scale mapping of channel features downstream from Fall Creek Dam shows minimal evidence of changes in channel planform or landforms that can be attributed to a drawdowns in WY 2012–16. On Upper Fall Creek reach, widespread increases in gravel bars or other in-channel sediment did not result from the five streambed drawdowns. The main changes attributable to sediment releases from Fall Creek Lake were localized increases in vegetated bar area, particularly on channel margin areas where sand and finer-grain sediment was deposited and rapidly colonized by vegetation. The area of mapped secondary water features decreased between 2005 and 2016, but that may be due to lower discharges depicted in the 2016 aerial photographs and less mapped area of inundation. Primary changes along the Lower Fall Creek reach include a 6.4 percent decrease in area of secondary water features between 2011 and 2016, and a nearly twofold increase in the area of unvegetated bars. Immediately downstream from the Fall Creek confluence, there were negligible changes in the location and areas of vegetated bars and the main wetted channel between 2005 and 2016, and local increases in bar area cannot be attributed solely to deposition of reservoir sediments from Fall Creek Lake because (1) areas along the Middle Fork Willamette River just upstream from the Fall Creek confluence display similar type and magnitude of changes and (2) some of the increases at the confluence area pre-date the drawdowns. The cumulative effect of sediment releases from Fall Creek Lake streambed drawdowns from WY 2012 to 2016 on downstream channel planform and landforms are modest compared to the river-scale transformations and planform changes that occurred in the decades following dam construction. Chapter E. Discussion of Geomorphic Responses of Fall Creek and Middle Fork Willamette River to Streambed Drawdowns at Fall Creek Lake Multiple aspects of Fall Creek Dam infrastructure and operations exert first-order controls on the magnitudes of reservoir erosion that occur during the streambed drawdowns and ultimately determine the sediment delivery to downstream reaches. Key aspects of the dam and its operations that are most relevant to assessing geomorphic responses to streambed drawdowns include the (1) dam infrastructure, including configuration and size of regulating outlets and their proximity to the streambed which dictates the capacity and competence of the river to deliver sediment to downstream reaches and mode of sediment transport as suspended-sediment load or bedload; (2) frequency of historical drawdowns and long-term, year-round dam operations and lake level management, which partly dictate reservoir morphology and locations and magnitudes of readily erodible materials; (3) dam operations and hydroclimatic conditions during the streambed drawdown (including length of the drawdown and streamflows entering the reservoir), which directly control the timing, duration and magnitude of reservoir erosion and sediment evacuation; and (4) dam operations following the streambed drawdown operation that regulate streamflows (and thereby sediment transport conditions) downstream of Fall Creek Dam which primarily reflect interactions between hydroclimatic conditions and flood control operations. Patterns of sediment erosion and evacuation observed in this study at Fall Creek Lake from WY 2012–18 suggest that reservoir erosion during annual streambed drawdowns may remain similar or decrease in future years assuming (1) annual streambed drawdown operations are implemented in similar manner as the WY 2012–18 drawdowns (in terms of duration, late autumn or early winter implementation, rate of pool-level lowering to reach streambed, and other factors), (2) streambed drawdowns coincide with similar conditions as were observed WY 2012–18 (similar sediment yield into reservoir, low reservoir inflows, limited precipitation, moderate air temperature), and (3) no major geomorphic changes in the main reservoir channels of Fall and Winberry Creeks occur (for example, channel avulsion). Under such conditions, it is hypothesized that the stream channel within the reservoir would achieve a quasi-equilibrium state with respect to annual influx and export of sediment and aided by the substantial amount of in-channel bedrock, will remain laterally stable without erosion across reservoir deposits. Patterns of sediment transport measured at the Fall Creek streamgage downstream from Fall Creek Dam provide insight into the potential effects of future streambed drawdowns at Fall Creek Lake. Analyses of suspended sediment measured in WY 2013–18 show a major reduction in suspended-sediment loads between WY 2013 and later years, indicating streamflows transporting sediment through the reservoir to downstream reaches during streambed drawdowns have become supply limited. The 6-year suspended-sediment monitoring and sampling program is insufficient to make predictions about future sediment transport conditions because of uneven monitoring periods and varying controls on reservoir sediment erosion. It is likely that future suspended-sediment loads will be variable but similar to those observed in WY 2015–18 if operational, climatic, and geomorphological factors remain similar to those monitored WY 2015–18. Suspended-sediment loads downstream from Fall Creek Lake will likely remain highest when regulating outlets are fully open and Fall Creek is free flowing with the reservoir fully drained with little to no residual pool. Over time, it is possible that the suspended-sediment loads would reflect mobilization of reservoir sediment deposited in the previous year rather than erosion of sediment deposited years or decades earlier. Bedload is likely to remain a small fraction of the total sediment load evacuated from the reservoir and is relatively modest compared with pre-dam bedload transport rates because most coarse sediment remains trapped by the dam. If sediment releases from Fall Creek Lake and ensuing streamflow conditions follow a similar pattern in the future as was assessed in this study spanning WY 2012–18, near-term geomorphic adjustments downstream of the dam are expected to be modest. Barring major operational, climatic, and geomorphological changes, local site-scale deposition on bars, overbank areas, or off-channel features that persists several months after the streambed drawdown will likely continue to be highly variable, ranging from negligible to several centimeters of deposition. At the landform-scale, low velocity areas nearest to Fall Creek Dam will likely continue to undergo rapid deposition immediately during and after a streambed drawdown event, similar to patterns observed for WY 2012–18. Some of the sediment entering these off-channel features and margin areas may be temporarily stored, then later remobilized and dispersed farther downstream. But if newly deposited sediment persists through the following spring, there is a greater likelihood that local vegetation will establish, reinforce deposited material, and trap sediment during later drawdowns. The reach-scale geomorphic changes may become more apparent if (1) streambed drawdowns continued for several decades, and geomorphic changes were measured at decadal scales or (2) the amount of sediment introduced to downstream reaches substantially increased and (or) sediment transport capacity decreased. The continued streamflow regulation of Fall Creek Dam after sediment releases provides an opportunity to strategically manage streamflows during and after the streambed drawdowns to minimize downstream sediment impacts and ensure other operational thresholds are satisfied. This study provides a comprehensive foundation of datasets and geomorphic analyses to inform dam operations at Fall Creek Lake, monitor sediment transport downstream, and consider operational schemes for future drawdowns. The datasets from this study also provide baselines of sediment transport and geomorphic conditions to assess future changes in reservoir and downstream environments. Future monitoring could be tailored to address specific questions regarding the long-term geomorphic effects of streambed drawdowns on fluvial habitats, flood hazards, cultural resources, or downstream water quality. Future monitoring activities could focus on the relevant geomorphic processes and spatial domains within the three categories used for this study: (1) reservoir erosion and net sediment evacuation, (2) sediment delivery to downstream reaches, including magnitude and temporal pattern of sediment transport, and (3) geomorphic responses of downstream reaches to sediment delivery. Specifically, high priority future monitoring activities could include: Repeat topographic or photographic surveys in the reservoir to characterize changes occurring within individual drawdowns, to quantify sediment export, to determine temporal changes in reservoir storage, and to identify locations of erosion and deposition. Continuous, year-round turbidity monitoring supplemented with suspended-sediment measurements at a streamflow-gaging station immediately downstream from the dam to quantify sediment export. Repeat geomorphic monitoring, mapping, or modeling in downstream reaches to track changes in channel and over bank features using a combination of site- and reach-scale monitoring approaches. This could support assessments of sediment deposition and ensuing vegetation encroachment on flood hazards and habitats and examine how sediment transport and depositional processes may be affected by different sediment supply, streamflow, or dam management scenarios.

Oregon↗