USGS Science⌕ Search

SEARCH · USGS Science

Results for “Final Management Report”

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

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

At least 451 records · Page 25Linked to original sources

Delineation of marsh types from Corpus Christi Bay, Texas, to Perdido Bay, Alabama, in 2010

Coastal zone managers and researchers often require detailed information regarding emergent marsh vegetation types (that is, fresh, intermediate, brackish, and saline) for modeling habitat capacities and needs of marsh dependent taxa (such as waterfowl and alligator). Detailed information on the extent and distribution of emergent marsh vegetation types throughout the northern Gulf of Mexico coast has been historically unavailable. In response, the U.S. Geological Survey, in collaboration with the Gulf Coast Joint Venture, the University of Louisiana at Lafayette, Ducks Unlimited, Inc., and the Texas A&M University-Kingsville, produced a classification of emergent marsh vegetation types from Corpus Christi Bay, Texas, to Perdido Bay, Alabama. This study incorporates about 9,800 ground reference locations collected via helicopter surveys in coastal wetland areas. Decision-tree analyses were used to classify emergent marsh vegetation types by using ground reference data from helicopter vegetation surveys and independent variables such as multitemporal satellite-based multispectral imagery from 2009 to 2011, bare-earth digital elevation models based on airborne light detection and ranging (lidar), alternative contemporary land cover classifications, and other spatially explicit variables. Image objects were created from 2010 National Agriculture Imagery Program color-infrared aerial photography. The final classification is a 10-meter raster dataset that was produced by using a majority filter to classify image objects according to the marsh vegetation type covering the majority of each image object. The classification is dated 2010 because the year is both the midpoint of the classified multitemporal satellite-based imagery (2009–11) and the date of the high-resolution airborne imagery that was used to develop image objects. The seamless classification produced through this work can be used to help develop and refine conservation efforts for priority natural resources.

Alabama, Louisiana, Mississippi, Texas↗

Preliminary hydrogeologic assessment near Tassi and Pakoon Springs, western part of Grand Canyon-Parashant National Monument, Arizona

Tassi and Pakoon Springs are both in the Grand Wash Trough in the western part of Grand Canyon-Parashant National Monument on the Arizona Strip. The monument is jointly managed by the National Park Service (NPS) and the Bureau of Land Management. This study was in response to NPS’s need to better understand the influence from regional increases in groundwater withdrawals near Grand Canyon-Parashant on the groundwater discharge from Tassi and Pakoon Springs. The climate of the Arizona Strip is generally semiarid to arid, and springs in the monument provide the water for the fragile ecosystems that are commonly separated by large areas of dry washes in canyons with pinyon and juniper. Available hydrogeologic data from previous investigations included water levels from the few existing wells, location information for springs, water chemistry from springs, and geologic maps. Available groundwater-elevation data from the wells and springs in the monument indicate that groundwater in the Grand Wash Trough is moving from north to south, discharging to springs and into the Colorado River. Groundwater may also be moving from east to west from Paleozoic rocks in the Grand Wash Cliffs into sedimentary deposits in the Grand Wash Trough. Finally, groundwater may be moving from the northwest in the Mesoproterozoic crystalline rocks of the Virgin Mountains into the northern part of the Grand Wash Trough. Water discharging from Tassi and Pakoon Springs has a major-ion chemistry similar to that of other springs in the western part of Grand Canyon-Parashant. Stable-isotopic signatures for oxygen-18 and hydrogen-2 are depleted in the water from both Tassi and Pakoon Springs in comparison to other springs on the Arizona Strip. Tassi Spring discharges from multiple seeps along the Wheeler Fault, and the depleted isotopic signatures suggest that water may be flowing from multiple places into Lake Mead and seems to have a higher elevation or an older climate source. Elevated water temperatures and a depleted stable-isotopic signature for Pakoon Springs suggest that the water may be traveling along a deep circulating flowpath, have multiple sources of water, been recharged at a high elevation, and (or) has an older climate source.

Arizona↗

Estimating freshwater productivity, overwinter survival, and migration patterns of Klamath River Coho Salmon

An area of great importance to resource management and conservation biology in the Klamath Basin is balancing water usage against the life history requirements of threatened Coho Salmon. One tool for addressing this topic is a freshwater dynamics model to forecast Coho Salmon productivity based on environmental inputs. Constructing such a forecasting tool requires local data to quantify the unique life history processes of Coho Salmon inhabiting this region. Here, we describe analytical methods for estimating a series of sub-models, each capturing a different life history process, which will eventually be synchronized as part of a freshwater dynamics model for Klamath River Coho Salmon. Specifically, we draw upon extensive population monitoring data collected in the basin to estimate models of freshwater productivity, overwinter survival, and migration patterns. Our models of freshwater productivity indicated that high summer temperatures and high winter flows can both adversely affect smolt production and that such relationships are more likely in tributaries with naturally regulated flows due to substantial intraannual environmental variation. Our models of overwinter survival demonstrated extensive variability in survival among years, but not among rearing locations, and demonstrated that a substantial proportion (~ 20%) of age-0+ fish emigrate from some rearing sites in the winter. Our models of migration patterns indicated that many age-0+ fish redistribute in the basin during the summer and winter. Further, we observed that these redistributions can entail long migrations in the mainstem where environmental stressors likely play a role in cueing refuge entry. Finally, our models of migration patterns indicated that changes in discharge are important in cueing the seaward migration of smolts, but that the nature of this behavioral response can differ dramatically between tributaries with naturally and artificially regulated flows. Collectively, these analyses demonstrate that environmental variation interacts with most phases of the freshwater life history of Klamath River Coho Salmon and that anthropogenic environmental variation can have a particularly large bearing on productivity.

Arcata Fisheries Technical Report↗

Local to landscape-level controls of water fluxes through Hawaiian forests: Effects of invasive animals and plants on soil infiltration capacity across substrate and moisture gradients

Given the potential effect of invasive plants and animals to water fluxes through forests, the invasive-driven degradation of native ecosystems is a topic of great concern for many downstream land and water managers. The infiltration rate determines the partitioning between runoff and infiltration into soil in Hawaiian forests and beyond. Thus, to explore the ecohydrological effects of plant and animal invasion in mesic and wet forests in Hawaii, we measured soil infiltration capacity in multiple fenced (i.e., ungulate-free)/unfenced and native/invaded forest sites along moisture and substrate age gradients across the islands of Hawai‘i and Kaua‘i. We also characterized forest composition and structure and soil characteristics at these sites to assess the direct and vegetation-mediated impacts of invasive species on infiltration capacity. Infiltration capacity is highly variable across forested sites and the wider landscape. Much of this variability is determined by a complex set of soil, vegetation, and disturbance factors that affect infiltration capacity at the immediate surrounding of measurement plots. Consequently, the effect of any given factor can be masked by variability in other factors. However, by controlling for variability in soil and vegetation conditions at a local plot level, we found that the presence of invasive species in forests has complex and sometimes non-intuitive effects on infiltration. Our final models showed that invasive ungulates negatively affect soil infiltration capacity consistently across the wide moisture and substrate age gradients considered. Additionally, because several soil characteristics known to be affected by ungulates were associated with local infiltration rates (e.g., soil organic matter, bare soil cover, soil depth), the long-term secondary effects of high ungulate densities in Hawaiian forests may be higher than effects observed in this study. These results provide clear evidence for land managers that ungulate control efforts likely improve ecohydrologic function to mesic and wet forest systems critical to protecting downstream and nearshore resources and maintaining groundwater recharge. Compared to ungulate effects, the effect of invasive plants on water infiltration capacity in Hawaiian forests appeared much more complex. In general, elements of forest structure including increased canopy, understory and floor cover, greater presence of large roots, and lower grass and bare soil covers were positively associated with water infiltration. Whether native or not, a plant species’ potential to alter infiltration rates in Hawaiian forests was likely to depend on its physiognomy and how it affects forest community structure. For instance, while the cover of native dominant tree ‘ōhi‘a, Metrosideros polymorpha, was found to be positively associated with infiltration capacity (perhaps as an indicator of overall forest integrity), invasive Himalayan ginger, Hedychium gardnerianum, was also positively correlated with infiltration capacity, possibly due to preferential flow channels created by the presence of large root mats. Few studies have conducted comprehensive integrated ecological and hydrological sampling in forests of high conservation value. While we show there are large benefits to understanding how conservation efforts may help shape water fluxes, we also found that the commonly used study design for infiltration studies used here and elsewhere (i.e., adjacent paired sites) could be modified to provide more accurate effects of invasion in future studies for ecosystems in Hawaii and beyond.

Hawaii↗

Species of conservation concern and environmental stressors: Local regional and global effects

Species conservation has traditionally been based on individual species within the context of their requisite habitat, which is generally defined as the communities and ecosystems deemed necessary for their presence. Conservation decisions are hampered by the fact that environmental stressors that poetically threaten the persistence of species can operate at organizational levels larger than the habitat or home range of a focal species. Resource managers must therefore simultaneously consider local, regional, and/or global scale stressors for effective conservation and management of species of concern. The wide ranging effects associated with global stressors such as climate change may exceed or exacerbate the effects of local or regional stressors, they still need to understand the direct and interactive effects of global stressors and ultimately how they affect the lands they manage. Conservation of species in southern Nevada is further complication by the fact that the region includes one of the largest and fastest growing urban centers in North America. To accomplish the goal of species conservation, resource managers must identify actionable management options that mitigate the effects of local and regional stressor in the context of the effects of global stressors that are beyond their control. Species conservation is typically focused on a subset often referred to as species of conservation concern that have either demonstrated considerable decline or are naturally rare or have limited distributions. Stressors can directly and indirectly impact species in a variety of ways and through a diversity of mechanisms. Some stressors have been more intense in the past (e.g., livestock grazing) whereas other are now only emerging as new stressors (e.g., solar energy development, climate change). The primary stressors affecting southern Nevada ecosystems are listed in table 2.1 and reviewed in detail in Chapter 2. This chapter addresses Dub-goal 1.4 in the SNAP Science Research Strategy which is to sustain and enhance southern Nevada's biotic communities to preserve biodiversity wand maintain viable populations (table 1.3; Turner and others 2009). We provide numerous examples of how stressors affect the range and/or habitat of select species of conservation concern. It is important to note that the species or groups discussed in this chapter by no means represent a comprehensive treatment of all species of conservation concern listed in Table 1.2 (Chapter 1). Rather, several species were chosen as examples for each southern Nevada ecosystem type to illustrate how stressors and linkages among them can affect species of conservation concern, keeping in mind that many of the species considered here are found in more than one ecosystem type. In addition, the stressors that may impact a species in one ecosystem may not be those that affect it in another ecosystem and different species in the same ecosystem may not be affected by the same suite of stressors. Finally, at the start of each ecosystem section we summarize key resource concerns, species used as examples, key stressors, and potential synergistic effects of those stressors relative to the species example.

Nevada↗

Forecasting stream habitat and Brook Trout responses to climate change in Catoctin Mountain Park

Anticipating and mitigating the effects of climate change is a fundamental challenge for natural resource conservation. In this report, we respond to research needs identified by Catoctin Mountain Park (CATO) for native Brook Trout (Salvelinus fontinalis) conservation and management as part of the US Geological Survey (USGS) Natural Resources Preservation Program in FY15-16. We addressed three overarching research questions: (1) How will anticipated changes in air temperature affect stream habitats? (2) How will changes to stream habitat affect the distribution of Brook Trout? (3) Which stream segments are most and least vulnerable to the effects of climate change? First, we surveyed Brook Trout abundance and fish community composition using electrofishing techniques within three watersheds: Owens Creek, upper Big Hunting Creek, and Blue Blazes Creek (a tributary to Big Hunting Creek). Second, we deployed a network of stream temperature gages to assess spatial variation in stream temperature and groundwater (GW) influence. Third, we used modeling techniques to forecast future stream temperatures that account for GW influences and air temperature scenarios. Fish sampling detected 13 species and 15,345 individual fish, the majority of which were Blacknose Dace (60%), Blue Ridge Sculpin (26%), and Brook Trout (6%). Brook Trout were not observed in Blue Blazes Creek and exhibited higher densities in Owens Creek than upper Big Hunting Creek (average densities = 19 fish/100 m and 4 fish/100 m, respectively). In contrast, Brown Trout were present in Blue Blazes Creek and exhibited greater density in Blue Blazes Creek than either Owens Creek or upper Big Hunting Creek (average densities = 3.0 fish/100 m, 0.3 fish/100 m, and 1.7 fish/100 m, respectively). Brown Trout occurred in sympatry with Brook Trout in Owens Creek and upper Big Hunting Creek, but appeared to have replaced Brook Trout in Blue Blazes Creek. Our fish surveys also revealed important locations for Brook Trout reproduction and young-of-year (YOY) dispersal within the Owens Creek watershed. Our study also revealed surprising differences in the distribution of Blue Ridge Sculpin among CATO streams. This species was abundant in Owens Creek (average density = 83 fish/100 m) but was less common in Blue Blazes Creek (average density = 12 fish/100 m) and was not detected in upper Big Hunting Creek. Histological examination of several specimens from Blue Blazes Creek by V. Blazer at the USGS Leetown Science Center revealed the presence of a novel parasite (Dermosystidium sp.) which has been linked to fish population declines elsewhere (Blazer et al. 2016). The parasite was not detected in Blue Ridge Sculpin samples from Owens Creek, and all trout appeared to be uninfected. Our survey results suggest that Blue Ridge Sculpin have been extirpated from upper Big Hunting Creek and have not recolonized from downstream source populations due to the fish passage barrier of Cunningham Falls. We recommend additional research to (1) evaluate the feasibility of reintroducing Blue Ridge Sculpin into upper Big Hunting Creek and (2) continue monitoring the distribution and potential spread of Dermocystidium in downstream waters. Stream temperatures ranged from 9.6 – 27.6 ºC during baseflow conditions in 2015 and 2016. Sites within upper Big Hunting Creek were consistently warmer than in Owens Creek or Blue Blazes Creek, suggesting an effect of headwater ponds outside CATO on upper Big Hunting Creek temperatures. For instance, in 2016 the maximum observed temperature in upper Big Hunting Creek was 27.6 ºC whereas Owens Creek reached a maximum of 23.7 ºC that year. Stream temperature data also revealed that 2016 was warmer than 2015 throughout the study area but did not exceed thermal tolerance limits for Brook Trout in either year. We estimated the influence of GW on stream temperatures using a statistical modeling approach based on the relationship between daily mean air temperature and stream temperature over time. Results indicated that effects of GW were generally stronger in the Owens Creek watershed than in Blue Blazes or upper Big Hunting Creek. However, we detected substantial spatial variation in GW influence among Owens Creek sites, with stream temperatures at some locations showing relatively little GW influence and others showing very strong influences (and correspondingly small influence of daily mean air temperatures). Although incoming lateral seeps were detected in upper Big Hunting Creek (D. Ferrier, Hood College, personal communication), the strongest effects of GW in the study area were due to GW upwelling within portions of the Owens Creek watershed (i.e., Tributary C in Figure 4) where we also observed high numbers of Brook Trout juveniles. Our results therefore identified potential high-priority areas for Brook Trout conservation in CATO. Finally, we modeled future stream temperatures based on scenarios characterizing GW sensitivity to air temperature and future air temperature increases. Stream temperature forecasts revealed important differences in habitat suitability for Brook Trout within and among watersheds. Big Hunting Creek sites were generally more sensitive to air temperature increases than sites in Owens Creek or Blue Blazes Creek. For instance, an increase in mean annual air temperature of 1.5 ºC (lowest level evaluated) exceeded thermal thresholds for Brook Trout in the majority of sites within that watershed, regardless of GW influence levels. In contrast, an air temperature increase of 1.5 ºC did not exceed thermal thresholds for Brook Trout in Owens Creek. However, modeled air temperature increases of 5 ºC resulted in a loss of Brook Trout thermal suitability throughout the study area. Model results revealed spatially patchy responses to air temperature increases that could provide an early-warning system for trout monitoring designs in CATO.

Maryland↗

A Natural History Summary and Survey Protocol for the Southwestern Willow Flycatcher

The Southwestern Willow Flycatcher (Empidonax traillii extimus) has been the subject of substantial research, monitoring, and management activity since it was listed as an endangered species in 1995. When proposed for listing in 1993, relatively little was known about the flycatcher's natural history, and there were only 30 known breeding sites supporting an estimated 111 territories rangewide (Sogge and others, 2003a). Since that time, thousands of presence/absences surveys have been conducted throughout the historical range of the flycatcher, and many studies of its natural history and ecology have been completed. As a result, the ecology of the flycatcher is much better understood than it was just over a decade ago. In addition, we have learned that the current status of the flycatcher is better than originally thought: as of 2007, the population was estimated at approximately 1,300 territories distributed among approximately 280 breeding sites (Durst and others, 2008a). Concern about the Southwestern Willow Flycatcher on a rangewide scale was brought to focus by Unitt (1987), who described declines in flycatcher abundance and distribution throughout the Southwest. E. t. extimus populations declined during the 20th century, primarily because of habitat loss and modification from activities, such as dam construction and operation, groundwater pumping, water diversions, and flood control. In 1991, the U.S. Fish and Wildlife Service (USFWS) designated the Southwestern Willow Flycatcher as a candidate category 1 species (U.S. Fish and Wildlife Service, 1991). In July 1993, the USFWS proposed to list E. t. extimus as an endangered species and to designate critical habitat under the Act (U.S. Fish and Wildlife Service, 1993). A final rule listing E. t. extimus as endangered was published in February 1995 (U.S. Fish and Wildlife Service, 1995); critical habitat was designated in 1997 (U.S. Fish and Wildlife Service, 1997). The USFWS Service released a Recovery Plan for the Southwestern Willow Flycatcher in 2002 (U.S. Fish and Wildlife Service, 2002), and re-designated critical habitat in 2005 (U.S. Fish and Wildlife Service, 2005).

Techniques and Methods↗

Estimates of peak flood discharge for 21 sites in the Front Range in Colorado in response to extreme rainfall in September 2013

Extreme rainfall in September 2013 caused destructive floods in part of the Front Range in Boulder County, Colorado. Erosion from these floods cut roads and isolated mountain communities for several weeks, and large volumes of eroded sediment were deposited downstream, which caused further damage of property and infrastructures. Estimates of peak discharge for these floods and the associated rainfall characteristics will aid land and emergency managers in the future. Several methods (an ensemble) were used to estimate peak discharge at 21 measurement sites, and the ensemble average and standard deviation provided a final estimate of peak discharge and its uncertainty. Because of the substantial erosion and deposition of sediment, an additional estimate of peak discharge was made based on the flow resistance caused by sediment transport effects. Although the synoptic-scale rainfall was extreme (annual exceedance probability greater than 1,000 years, about 450 millimeters in 7 days) for these mountains, the resulting peak discharges were not. Ensemble average peak discharges per unit drainage area (unit peak discharge, [ Q u ]) for the floods were 1–2 orders of magnitude less than those for the maximum worldwide floods with similar drainage areas and had a wide range of values (0.21–16.2 cubic meters per second per square kilometer [m 3 s -1 km -2 ]). One possible explanation for these differences was that the band of high-accumulation, high-intensity rainfall was narrow (about 50 kilometers wide), oriented nearly perpendicular to the predominant drainage pattern of the mountains, and therefore entire drainage areas were not subjected to the same range of extreme rainfall. A linear relation (coefficient of determination [ R 2 ]=0.69) between Q u and the rainfall intensity ( I Tc , computed for a time interval equal to the time-of-concentration for the drainage area upstream from each site), had the form: Q u =0.26( I Tc -8.6), where the coefficient 0.26 can be considered to be an area-averaged peak runoff coefficient for the September 2013 rain storms in Boulder County, and the 8.6 millimeters per hour to be the rainfall intensity corresponding to a soil moisture threshold that controls the soil infiltration rate. Peak discharge estimates based on the sediment transport effects were generally less than the ensemble average and indicated that sediment transport may be a mechanism that limits velocities in these types of mountain streams such that the Froude number fluctuates about 1 suggesting that this type of floodflow can be approximated as critical flow.

Colorado↗

Estimating flood magnitude and frequency on gaged and ungaged streams in Maine

Accurate estimates of flood frequency and magnitude on rivers and streams in Maine are a key component of effective flood risk management, flood mitigation, and flood recovery programs for the State. Flood-frequency estimates are published here for 148 streamgages in and adjacent to Maine. Equations are provided for users to compute flood-frequency estimates at any location on a stream that does not have a streamgage. Estimates and equations are presented for peak flows with annual exceedance probabilities (AEPs) of 50, 20, 10, 4, 2, 1, 0.5, and 0.2 percent. AEPs correspond to flood recurrence intervals of 2, 5, 10, 25, 50, 100, 200, and 500 years, respectively. New estimates use a regional skew coefficient of 0.02 with a standard error of prediction of 0.30 developed specifically for Maine as a part of this work. Equations are designed for use at ungaged sites without substantial flow regulation or urbanization in Maine, with drainage areas between 0.26 and 5,680 square miles. The equations were developed from streamflows and basin characteristics at 124 unregulated streamgages using generalized least-squares regression techniques. Explanatory variables used in the equations for computing peak flows are drainage area, percentage of area in the basin that contains wetlands, and basin mean 24-hour rainfall intensities. The average standard error of prediction (ASEP) for these equations ranges from −31.5 to 45.9 percent for the 50-percent AEP and from −34.2 to 52.0 percent for the 0.2-percent AEP. Equations that use only drainage area are provided for use in cases where lower accuracy is acceptable. The ASEP for estimating peak flows with these simpler equations ranges from −40 to 66 percent for the 50-percent AEP and from −44 to 79 percent for the 0.2-percent AEP. Final peak flows at unregulated streamgages are computed as weighted averages between the at-station peak flows and peak flows computed at those same sites using the regression equations. Peak flow estimates and equations presented here are accessible in the U.S. Geological Survey StreamStats application. StreamStats is a web application that computes selected basin characteristics and estimates of peak flows and other available streamflow statistics for user-selected streams in Maine.

Maine↗

User's Manual for the New England Water-Use Data System (NEWUDS)

Water is used in a variety of ways that need to be understood for effective management of water resources. Water-use activities need to be categorized and included in a database management system to understand current water uses and to provide information to water-resource management policy decisionmakers. The New England Water-Use Data System (NEWUDS) is a complex database developed to store water-use information that allows water to be tracked from a point of water-use activity (called a 'Site'), such as withdrawal from a resource (reservoir or aquifer), to a second Site, such as distribution to a user (business or irrigator). NEWUDS conceptual model consists of 10 core entities: system, owner, address, location, site, data source, resource, conveyance, transaction/rate, and alias, with tables available to store user-defined details. Three components--site (with both a From Site and a To Site), a conveyance that connects them, and a transaction/rate associated with the movement of water over a specific time interval form the core of the basic NEWUDS network model. The most important step in correctly translating real-world water-use activities into a storable format in NEWUDS depends on choosing the appropriate sites and linking them correctly in a network to model the flow of water from the initial From Site to the final To Site. Ten water-use networks representing real-world activities are described--three withdrawal networks, three return networks, two user networks, two complex community-system networks. Ten case studies of water use, one for each network, also are included in this manual to illustrate how to compile, store, and retrieve the appropriate data. The sequence of data entry into tables is critical because there are many foreign keys. The recommended core entity sequence is (1) system, (2) owner, (3) address, (4) location, (5) site, (6) data source, (7) resource, (8) conveyance, (9) transaction, and (10) rate; with (11) alias and (12) user-defined detail subject areas populated as needed. After each step in data entry, quality-assurance queries should be run to ensure the data are correctly entered so that it can be retrieved accurately. The point of data storage is retrieval. Several retrieval queries that focus on retrieving only relevant data to specific questions are presented in this manual as examples for the NEWUDS user.

Open-File Report↗

Effects of climate change on plague exposure pathways and resulting disease dynamics

Introduction and Objectives: Sylvatic plague, a zoonotic flea-borne disease, caused by the bacterium Yersinia pestis , is relevant to the Department of Defense (DOD), because prairie dogs and other susceptible rodents are present on military installations in several western states. Arthropod-borne diseases, like plague, are thought to be particularly sensitive to local climate conditions. Expected changes in temperature and humidity over the next several decades will likely increase the geographical expansion of plague outbreaks in wildlife. Through a combination of field and laboratory work, along with data-driven modeling, we evaluated the potential effects of climate change on plague exposure pathways in prairie dogs and associated rodents to provide guidance to DOD partners regarding the potential for future outbreaks. Briefly, our specific objectives were to determine the relation between local climate conditions and the prevalence of plague and other pathogens while assessing the ecological roles of specific rodent hosts and vector species in plague dynamics, evaluate flea intensity on rodent hosts and in burrows in relation to local climate conditions, and develop models to predict the effects of climate change on plague dynamics. Technical Approach: Using data and samples collected during a large field study on the effectiveness of vaccination to manage plague in prairie dogs, we assessed rodent/flea assemblages, pathogen prevalence in fleas, and determined how local climate conditions influence flea development rates and relative abundance. Live animals (prairie dogs and some small rodents) were trapped to collect fleas and other samples on 46 prairie dog plots in 6 western states, many sites near DOD lands. At seven additional locations on a latitudinal gradient, fleas were collected from burrows several times per year to assess seasonality and effects of local climate conditions on flea abundance. These data were then used to develop predictive models that could be used to test specific hypotheses. Results: We determined that flea developmental rates, on-host flea abundance, species composition of the flea community, and burrow temperatures varied across a latitudinal gradient. Rodent and flea community composition and abundance differed geographically and were highly specialized. Flea-switching between prairie dogs and short-lived rodents was rare. Flea development rates, on-host flea abundance, and burrow temperatures increased with increasing ambient temperature. Although relative humidity can affect flea development, burrow humidity was uniformly high (~85%) across sampling sites and seasons. A large increase in the number of fleas found on a prairie dog colony, coupled with a greater number of infested burrows, could have substantial effects on plague dynamics in the western United States as the climate warms. In addition to affecting flea load, climate change may also influence body condition of prairie dogs by reducing the amount of forage. This may result in animals being more tolerant of high flea loads (less engaged in grooming behavior) and more vulnerable to disease.

Arizona, Montana, South Dakota, Texas, Utah, Wyomi↗

Hydrologic Record Extension of Water-Level Data in the Everglades Depth Estimation Network (EDEN) Using Artificial Neural Network Models, 2000-2006

The Everglades Depth Estimation Network (EDEN) is an integrated network of real-time water-level gaging stations, ground-elevation models, and water-surface models designed to provide scientists, engineers, and water-resource managers with current (2000-present) water-depth information for the entire freshwater portion of the greater Everglades. The U.S. Geological Survey Greater Everglades Priority Ecosystem Science provides support for EDEN and the goal of providing quality assured monitoring data for the U.S. Army Corps of Engineers Comprehensive Everglades Restoration Plan. To increase the accuracy of the water-surface models, 25 real-time water-level gaging stations were added to the network of 253 established water-level gaging stations. To incorporate the data from the newly added stations to the 7-year EDEN database in the greater Everglades, the short-term water-level records (generally less than 1 year) needed to be simulated back in time (hindcasted) to be concurrent with data from the established gaging stations in the database. A three-step modeling approach using artificial neural network models was used to estimate the water levels at the new stations. The artificial neural network models used static variables that represent the gaging station location and percent vegetation in addition to dynamic variables that represent water-level data from the established EDEN gaging stations. The final step of the modeling approach was to simulate the computed error of the initial estimate to increase the accuracy of the final water-level estimate. The three-step modeling approach for estimating water levels at the new EDEN gaging stations produced satisfactory results. The coefficients of determination (R2) for 21 of the 25 estimates were greater than 0.95, and all of the estimates (25 of 25) were greater than 0.82. The model estimates showed good agreement with the measured data. For some new EDEN stations with limited measured data, the record extension (hindcasts) included periods beyond the range of the data used to train the artificial neural network models. The comparison of the hindcasts with long-term water-level data proximal to the new EDEN gaging stations indicated that the water-level estimates were reasonable. The percent model error (root mean square error divided by the range of the measured data) was less than 6 percent, and for the majority of stations (20 of 25), the percent model error was less than 1 percent.

Open-File Report↗

Northern Great Plains

In the Northern Great Plains, the timing and quantity of both precipitation and runoff have important consequences for water supplies, agricultural activities, and energy production. Overall, climate projections suggest that the number of heavy precipitation events (events with greater than 1 inch per day of rainfall) is projected to increase. Moving forward, the magnitude of year-to-year variability overshadows the small projected average decrease in streamflow. Changes in extreme events are likely to overwhelm average changes in both the eastern and western regions of the Northern Great Plains. Major flooding across the basin in 2011 was followed by severe drought in 2012, representing new and unprecedented variability that is likely to become more common in a warmer world. The Northern Great Plains region plays a critical role in national food security. Among other anticipated changes, projected warmer and generally wetter conditions with elevated atmospheric carbon dioxide concentrations are expected to increase the abundance and competitive ability of weeds and invasive species, increase livestock production and efficiency of production, and result in longer growing seasons at mid- and high latitudes. Net primary productivity, including crop yields and forage production, is also likely to increase, although an increasing number of extreme temperature events during critical pollination and grain fill periods is likely to reduce crop yields. Ecosystems across the Northern Great Plains provide recreational opportunities and other valuable goods and services that are ingrained in the region’s cultures. Higher temperatures, reduced snow cover, and more variable precipitation will make it increasingly challenging to manage the region’s valuable wetlands, rivers, and snow-dependent ecosystems. In the mountains of western Wyoming and western Montana, the fraction of total water in precipitation that falls as snow is expected to decline by 25% to 40% by 2100 under a higher scenario (RCP8.5), which would negatively affect the region’s winter recreation industry. At lower-elevation areas of the Northern Great Plains, climate-induced land-use changes in agriculture can have cascading effects on closely entwined natural ecosystems, such as wetlands, and the diverse species and recreational opportunities they support. Energy resources in the Northern Great Plains include abundant crude oil, natural gas, coal, wind, and stored water, and to a lesser extent, corn-based ethanol, solar energy, and uranium. The infrastructure associated with the extraction, distribution, and energy produced from these resources is vulnerable to the impacts of climate change. Railroads and pipelines are vulnerable to damage or disruption from increasing heavy precipitation events and associated flooding and erosion. Declining water availability in the summer would likely increase costs for oil production operations, which require freshwater resources. These cost increases will either lead to lower production or be passed on to consumers. Finally, higher maximum temperatures, longer and more severe heat waves, and higher overnight lows are expected to increase electricity demand for cooling in the summer, further stressing the power grid. Indigenous peoples in the region are observing changes to climate, many of which are impacting livelihoods as well as traditional subsistence and wild foods, wildlife, plants and water for ceremonies, medicines, and health and well-being. Because some tribes and Indigenous peoples are among those in the region with the highest rates of poverty and unemployment, and because many are still directly reliant on natural resources, they are among the most at risk to climate change (e.g., Gamble et al. 2016, Cozzetto et al. 2013, Espey et al. 2014, Wong et al. 2014, Kornfeld 2016, Paul and Caplins 2016, Maynard 2014, USGCRP 2017)

Report↗

Integrated geoscience studies in the Greater Yellowstone Area - Volcanic, tectonic, and hydrothermal processes in the Yellowstone geoecosystem

Yellowstone National Park, rimmed by a crescent of older mountainous terrain, has at its core the Quaternary Yellowstone Plateau, an undulating landscape shaped by forces of volcanism, tectonism, and later glaciation. Its spectacular hydrothermal systems cap this landscape. From 1997 through 2003, the United States Geological Survey Mineral Resources Program conducted a multidisciplinary project of Yellowstone National Park entitled Integrated Geoscience Studies of the Greater Yellowstone Area, building on a 130-year foundation of extensive field studies (including the Hayden survey of 1871, the Hague surveys of the 1880s through 1896, the studies of Iddings, Allen, and Day during the 1920s, and NASA-supported studies starting in the 1970s—now summarized in USGS Professional Paper 729 A through G) in this geologically dynamic terrain. The project applied a broad range of scientific disciplines and state-of-the-art technologies targeted to improve stewardship of the unique natural resources of Yellowstone and enable the National Park Service to effectively manage resources, protect park visitors from geologic hazards, and better educate the public on geologic processes and resources. This project combined a variety of data sets in characterizing the surficial and subsurface chemistry, mineralogy, geology, geophysics, and hydrothermal systems in various parts of the park. The sixteen chapters presented herein in USGS Professional Paper 1717, Integrated Geoscience Studies in the Greater Yellowstone Area—Volcanic, Tectonic, and Hydrothermal Processes in the Yellowstone Geoecosystem , can be divided into four major topical areas: (1) geologic studies, (2) Yellowstone Lake studies, (3) geochemical studies, and (4) geophysical studies. The geologic studies include a paper by Ken Pierce and others on the influence of the Yellowstone hotspot on landscape formation, the ecological effects of the hotspot, and the human experience and human geography of the greater Yellowstone ecosystem as influenced by the Yellowstone hotspot. Another paper by Paul Carrara describes the recent movement of a large landslide block dated by tree-ring analyses in the Tower Falls area. The section under Yellowstone Lake studies begins with a classic paper by J. David Love and others on ancestral Lake Yellowstone. Other papers in this section include results and interpretation of the high-resolution bathymetric, seismic reflection, and submersible studies by Lisa Morgan and others. Ken Pierce and others describe results from their studies of shorelines along Yellowstone Lake and their interpretation of inflation-deflation cycles, tilting, and faulting in the Yellowstone caldera. The influence of sublacustrine hydrothermal vent fluids on the geochemistry of Yellowstone Lake is described by Laurie Balistrieri and others. In Pat Shanks and others’ chapter, hydrothermal reactions, stable-isotope systematics, sinter deposition, and spire formation are related to the geochemistry of sublacustrine hydrothermal deposits in Yellowstone Lake. The geochemical studies section considers park-wide geochemical systems in Yellowstone National Park. In Bob Rye and Alfred Truesdell’s paper, the question of recharge to the deep thermal reservoir underlying the geysers and hot springs of Yellowstone National Park is discussed. Irving Friedman and Dan Norton report on the chloride flux emissions from Yellowstone in their paper questioning whether Yellowstone is losing its steam. Wildlife issues as addressed by examining trace-element and stable-isotope geochemistry are discussed in a chapter by Maurice Chaffee and others. In another chapter by Chaffee and others, natural and anthropogenic anomalies and their potential impact on the environment using geochemistry is reported. Pam Gemery-Hill and others present geochemical data for selected rivers, lake waters, hydrothermal vents, and subaerial geysers for the time interval of 1996–2004. The life cycle of gold deposits near the northeast corner of the park is discussed by Brad Van Gosen. Under the geophysical studies segment, Ray Kokaly and others use AVIRIS (Airborne Visible and Infrared Spectroscopy) data to map vegetation cover and microbial communities in Yellowstone National Park. Eric Livo and others report their results using AVIRIS data on hydrothermally altered rock and hot-spring deposits. In his final paper following a half century of scientific research, Irving Friedman presents data on monitoring changes in geothermal activity at Norris Geyser Basin using satellite telemetry. These papers summarize a near-decade-long effort by the USGS from the late 1990s to mid-2000s. In 2001, the USGS in cooperation with the National Park Service (Yellowstone National Park) and the University of Utah established the Yellowstone Volcano Observatory, the 5th volcano observatory in the United States.

Professional Paper↗

Application of FTLOADDS to Simulate Flow, Salinity, and Surface-Water Stage in the Southern Everglades, Florida

The Comprehensive Everglades Restoration Plan requires numerical modeling to achieve a sufficient understanding of coastal freshwater flows, nutrient sources, and the evaluation of management alternatives to restore the ecosystem of southern Florida. Numerical models include a regional water-management model to represent restoration changes to the hydrology of southern Florida and a hydrodynamic model to represent the southern and western offshore waters. The coastal interface between these two systems, however, has complex surface-water/ground-water and freshwater/saltwater interactions and requires a specialized modeling effort. The Flow and Transport in a Linked Overland/Aquifer Density Dependent System (FTLOADDS) code was developed to represent connected surface- and ground-water systems with variable-density flow. The first use of FTLOADDS is the Southern Inland and Coastal Systems (SICS) application to the southeastern part of the Everglades/Florida Bay coastal region. The need to (1) expand the domain of the numerical modeling into most of Everglades National Park and the western coastal area, and (2) better represent the effect of water-delivery control structures, led to the application of the FTLOADDS code to the Tides and Inflows in the Mangroves of the Everglades (TIME) domain. This application allows the model to address a broader range of hydrologic issues and incorporate new code modifications. The surface-water hydrology is of primary interest to water managers, and is the main focus of this study. The coupling to ground water, however, was necessary to accurately represent leakage exchange between the surface water and ground water, which transfers substantial volumes of water and salt. Initial calibration and analysis of the TIME application produced simulated results that compare well statistically with field-measured values. A comparison of TIME simulation results to previous SICS results shows improved capabilities, particularly in the representation of coastal flows. This improvement most likely is due to a more stable numerical representation of the coastal creek outlets. Sensitivity analyses were performed by varying frictional resistance, leakage, barriers to flow, and topography. Changing frictional resistance values in inland areas was shown to improve water-level representation locally, but to have a negligible effect on area-wide values. These changes have only local effects and are not physically based (as are the unchanged values), and thus have limited validity. Sensitivity tests indicate that the overall accuracy of the simulation is diminished if leakage between surface water and ground water is not simulated. The inclusion of a major road as a complete barrier to surface-water flow influenced the local distribution and timing of flow; however, the changes in total flow and individual creekflows were negligible. The model land-surface altitude was lowered by 0.1 meter to determine the sensitivity to topographic variation. This topographic sensitivity test produced mixed results in matching field data. Overall, the representation of stage did not improve definitively. A final calibration utilized the results of the sensitivity analysis to refine the TIME application. To accomplish this calibration, the friction coefficient was reduced at the northern boundary inflow and increased in the southwestern corner of the model, the evapotranspiration function was varied, additional data were used for the ground-water head boundary along the southeast, and the frictional resistance of the primary coastal creek outlet was increased. The calibration improved the match between measured and simulated total flows to Florida Bay and coastal salinities. Agreement also was improved at most of the water-level sites throughout the model domain.

Scientific Investigations Report↗

Status and trends of prey fish populations in Lake Michigan, 2008

The Great Lakes Science Center (GLSC) has conducted lake-wide surveys of the fish community in Lake Michigan each fall since 1973 using standard 12-m bottom trawls towed along contour at depths of 9 to 110 m at each of seven index transects. The resulting data on relative abundance, size structure, and condition of individual fishes are used to estimate various population parameters that are in turn used by state and tribal agencies in managing Lake Michigan fish stocks. All seven established index transects of the survey were completed in 2008. The survey provides relative abundance and biomass estimates between the 5-m and 114-m depth contours of the lake (herein, lake-wide) for prey fish populations, as well as burbot, yellow perch, and the introduced dreissenid mussels. Lake-wide biomass of alewives in 2008 was estimated at 8.27 kilotonnes (kt) (1 kt = 1000 metric tons), which was the smallest biomass estimate in the entire time series and 29% lower than the 2007 estimate. Lake-wide biomass of bloater in 2008 was estimated at 3.33 kt, which was the lowest estimate since 1977 and 38% lower than the 2007 estimate. Rainbow smelt lake-wide biomass equaled 0.89 kt, which was only 0.01 kt higher than 2007, which is the lowest estimate in the time series. Deepwater sculpin lake-wide biomass equaled 5.23 kt, which is the fourth straight year of declining biomass. The 2008 estimate is the second smallest in the time series, and 39% lower than the 2007 estimate. Slimy sculpin lake-wide biomass remained relatively high in 2008 (2.75 kt), increasing 25% over 2007. Ninespine stickleback lake-wide biomass equaled only 0.50 kt in 2008, which was 79% lower than the 2007 estimate. The final prey fish, exotic round goby, increased two orders of magnitude between 2007 and 2008, from 0.02 to 4.65 kt. Round gobies now represent 18% of the prey fish biomass. Burbot lake-wide biomass (0.91 kt in 2008) has remained fairly constant since 2002. Numeric density of age-0 yellow perch (i.e., < 100 mm) equaled 0.7 fish per ha, which is indicative of a relatively poor year-class. Lake-wide biomass of dreissenid mussels dropped precipitously in 2008, down to 9.47 kt, and a 96% decline from the 2007 biomass estimate. Overall, the total lake-wide prey fish biomass estimate (sum of alewife, bloater, rainbow smelt, deepwater sculpin, slimy sculpin, round goby, and ninespine stickleback) in 2008 was 25.62 kt, which was the lowest observed since the survey began in 1973.

Lake Michigan↗

Estimating Locations of Perennial Streams in Idaho Using a Generalized Least-Squares Regression Model of 7-Day, 2-Year Low Flows

Many State and Federal agencies use information regarding the locations of streams having intermittent or perennial flow when making management and regulatory decisions. For example, the application of some Idaho water quality standards depends on whether streams are intermittent. Idaho Administrative Code defines an intermittent stream as one having a 7-day, 2-year low flow (7Q2) less than 0.1 ft3/s. However, there is a general recognition that the cartographic representation of perennial/intermittent status of streams on U.S. Geological Survey (USGS) topographic maps is not as accurate or consistent as desirable from one map to another, which makes broad management and regulatory assessments difficult and inconsistent. To help resolve this problem, the USGS has developed a methodology for predicting the locations of perennial streams based on regional generalized least-squares (GLS) regression equations for Idaho streams for the 7Q2 low-flow statistic. Using these regression equations, the 7Q2 streamflow may be estimated for naturally flowing streams in most areas in Idaho. The use of these equations in conjunction with a geographic information system (GIS) technique known as weighted flow accumulation allows for an automated and continuous estimation of 7Q2 streamflow at all points along stream reaches. The USGS has developed a GIS-based map of the locations of streams in Idaho with perennial flow based on a 7Q2 of 0.1 ft3/s and a transition zone of plus or minus 1 standard error. Idaho State cooperators plan to use this information to make regulatory and water-quality management decisions. Originally, 7Q2 equations were developed for eight regions of similar hydrologic characteristics in the study area, using long-term data from 234 streamflow-gaging stations. Equations in five of the regions were revised based on spatial patterns observed in the initial perennial streams map and unrealistic behavior of the equations in extrapolation. The standard errors of prediction for the final equations ranged from a minimum of +75.0 to -42.9 percent in the central part of the study area to a maximum of +277 to -73.5 percent in the southern part of the study area. The equations are applicable only to unregulated, naturally-flowing streams and may produce unreliable results outside the range of explanatory variables used for equation development. Extrapolation outside the range of available data was necessary, however, to predict perennial flow initiation points and transition zones along stream reaches. The map of perennial streams was evaluated by comparing predicted stream classifications with four independent datasets, including field observations by other government agencies. Overall, 81 percent of the comparison data points agreed with the USGS perennial streams model. Regions with the highest number of disagreements had a high percentage of mountainous and forested area with potential mountain front recharge zones, and regions with the highest agreements had a high percentage of low gradient, low elevation area. As a whole, the USGS model predicted a higher number of perennial streams than predictions made with the independent datasets. Some disagreements were due to poor site location coordinates, timing of the comparison site visits during unusually wet or dry years, discrepancies in classification criteria, and variable ground water contributions to flow in some areas. The Idaho Department of Environmental Quality Beneficial Use Reconnaissance Program (BURP) dataset is considered the most representative dataset for comparison because it covered a range of climate conditions and the number of sites visited were consistent from year to year during the study period. Eighty-five percent of BURP comparison data points agreed with the USGS perennial streams model. Although site-specific flow data may be needed to correctly classify streams in some areas, this information rarely is available and is not always practical to o

Scientific Investigations Report↗

Development of an aerial population survey method for elk (Cervus elaphus) in Rocky Mountain National Park, Colorado

Executive Summary Since the early 1990s, substantial effort and funding have been expended to conduct research to guide development of a 20-year Elk and Vegetation Management Plan for Rocky Mountain National Park (RMNP) in Colorado. One goal of the plan is to maintain the elk ( Cervus elaphus ) population size at the lower end of the natural range of variation. To implement management actions called for in the plan, accurate and reliable population estimates are needed, as well as a better understanding of the spatial and temporal distribution of elk. The previous aerial survey protocol and population estimation model used by the park had not been updated since the model’s initial calibration more than 15 years ago and the model was developed with an insufficient number (n=44) of observations. Thus we initiated research to reevaluate, update, and improve elk population estimation protocols for RMNP. We considered several alternative survey and analysis methods, and concluded that a hybrid population-estimation model using a simultaneous double-observer technique with sighting covariates was the most appropriate and effective aerial survey methodology for this population and the environmental conditions in RMNP. Instructional protocols for conducting these surveys in the future, along with datasheets, are provided in this report’s appendixes. To develop an improved method for aerial elk surveys, we used elk radio-collar location data from our study and other studies, and applied geographic information system analyses to define the survey area and develop effective, repeatable survey transect lines. We used telemetry data from radio-collared elk to inform our understanding of the temporal and spatial scale of elk movements across the park boundary, elk use of tree cover during potential survey hours, and elk use of different elevations within their range. Determining where elk were during surveys helped to fine-tune a survey design that improved spatial cover-age and decreased costs where possible, while standardizing the method to make it repeatable from year to year. We conducted aerial helicopter surveys to test our methodology in an adaptive, iterative process during three winters: 2007–2008, 2008–2009, and 2009–2010. We gained new information on each survey and used results to refine subsequent surveys. Our results confirm that elk movements were highly dynamic with respect to park boundary crossings; an average of six round trips from the park to Estes Park, Colorado, and back per month were taken by global positioning system-collared bull elk. We observed a strong diurnal temporal pattern of bull elk use of trees; elk were found in dense tree cover from 15:00−22:00 but not as often during morning and early afternoon hours when surveys were conducted. We used information on elk use of different altitudes to refine and establish a more efficient survey area. During the time of our study, a concurrent study in the park deployed 120 very high frequency radio collars on elk cows. We used those collar locations during our flights to increase sample size and to evaluate the level of precision we would gain by using “known fates analysis” (in which elk were known to be in the survey area, out of the survey area, or deceased based on radio-collar locations collected simultaneously during aerial surveys). Using radio-collar locations reduced bias by 1.1–8.8 percent and increased precision (that is, reduced the width of confidence intervals). This report provides final population estimates analyzed with and without radio-collar data to demonstrate what is gained by using marked individuals.

Colorado↗