Search USGSSearch

SEARCH · Search USGS

Results for “Final 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 1,009 records · Page 56Linked to original sources

Water-quality constituent concentrations and loads computed using real-time water-quality data for the Republican River, Clay Center, Kansas, August 2018 through July 2023

Milford Lake, the largest reservoir by surface area in Kansas, has had confirmed harmful algal blooms every summer since reporting began in 2011, except 2018–19. Milford Lake has been listed as impaired and designated hypereutrophic under section 303(d) of the 1972 Clean Water Act. In 2014, the Kansas Department of Health and Environment established a total maximum daily load for eutrophication and dissolved oxygen impairments. In 2018, the Natural Resources Conservation Service funded the Regional Conservation Partnership Program for the Milford Lake Watershed to focus on best management practices in the Lower Republican River Basin. The U.S. Geological Survey, in cooperation with the Kansas Water Office, completed this study to assess and quantify water-quality constituent concentrations and loads for total nitrogen (TN), total phosphorus (TP), and suspended sediment (SS) using previously published models for the Republican River near Clay Center, Kansas (U.S. Geological Survey station 06856600), about 15 miles upstream from Milford Lake, during August 1, 2018, through July 31, 2023. TN, TP, and SS concentrations and loads were monitored because of their relation to water supply and water-quality issues in Milford Lake, including nutrient and sediment transport, taste-and-odor events, potentially toxic cyano-harmful algal bloom events, and subsequent downstream transport of contaminants. Data from this report can be used to evaluate changing conditions, provide science-based information for decision making, and help meet regulatory requirements. The study mean annual loads for TN and TP were greater than the reported mean annual total maximum daily load and exceeded the watershed reduction goals as well as Kansas nonpoint source reduction goals defined by the Watershed Restoration and Protection Strategy for the Lower Republican watershed. TN and TP annual loads during 2019–20 were greater than the defined mean annual total maximum daily load. During 2022, TN and TP annual loads were less than the Kansas nonpoint source reduction goal and during 2023 were less than the watershed reduction goal. SS loads were less than the mean annual sedimentation rate computed from the total maximum daily load for the entirety of the study period, and the study mean annual load was 72 percent less than the designed annual reservoir sedimentation rate for Milford Lake. Data collected during the study period represented a wide range of streamflow and water-quality conditions at the Clay Center site, ranging from low-flow with less frequent runoff during 2023 to high-flow with frequent runoff during 2018. Nutrient reduction goals were only met in the final 2 years of the study period when annual mean flow conditions were lower than normal, indicating that goals may be unattainable during average or high-flow conditions. In all years except 2019, the annual mean SS load was less than the 20-year sediment load reduction target. Although annual SS loads at the Clay Center site generally decreased over time, corresponding reductions in annual streamflow indicated that these reductions may primarily be related to less frequent runoff from the upstream basin. Continued water-quality monitoring and tracking of best management practices are necessary to understand the success of Regional Conservation Partnership Program efforts to reduce nutrient transport in the Milford Lake Watershed.

Kansas

Groundwater-flow model of the Treasure Valley, southwestern Idaho, 1986–2015

Most of the population of the Treasure Valley and the surrounding area of southwestern Idaho and easternmost Oregon depends on groundwater for domestic supply, either from domestic or municipal-supply wells. Current and projected rapid population growth in the area has caused concern about the long-term sustainability of the groundwater resource. In 2016, the U.S. Geological Survey, in cooperation with the Idaho Water Resource Board and the Idaho Department of Water Resources, began a project to construct a numerical groundwater-flow model of the westernmost portion of the western Snake River Plain aquifer system, called the Treasure Valley. The development of the model was guided by several objectives, including: to improve the understanding of groundwater and surface water interactions; to facilitate conjunctive water management; to provide a tool for water resources planning; and to provide a tool for water allocation. The model was constructed with a spatial scale and level of detail that aimed to meet these objectives while balancing the sometimes-competing goals of fast runtimes, numerical stability, usability, and parsimony. The Treasure Valley Groundwater Flow Model (TVGWFM) is a three-dimensional finite-difference numerical model constructed using MODFLOW 6 (Langevin and others, 2017, Documentation for the MODFLOW 6 Groundwater Flow Model: U.S. Geological Survey Techniques and Methods, book 6, chap. A55, 197 p., https://doi.org/10.3133/tm6A55 ). The model covers the westernmost portion of the western Snake River Plain and is discretized into a regular grid of 64 by 65 cells with a side length of 1 mile and 6 layers of varying depth and active area. A historical model period was developed consisting of 360 month-long stress periods for 1986–2015. The model builds upon previous modeling efforts by adding a transient period, incorporating new head and discharge observations to constrain parameters, incorporating information from the hydrogeologic framework model (HFM) of Bartolino (2019, Hydrogeologic framework of the Treasure Valley and surrounding area, Idaho and Oregon: U.S. Geological Survey Scientific Investigations Report 2019–5138, https://doi.org/10.3133/sir20195138 ) and incorporating refined estimates of evapotranspiration and irrigation classification of lands in the study area. The TVGWFM includes all significant components of recharge to and discharge from the aquifer. Inflows include canal seepage, irrigation and precipitation recharge, mountain-front recharge, rivers and stream seepage, and seepage from Lake Lowell. Outflows include discharge to agricultural drainage ditches, discharge to rivers and streams, pumping, and discharge to Lake Lowell. Each recharge or discharge component is represented separately using individual MODFLOW 6 packages. Parameter values were derived with a combination of trial-and-error steps and automated parameter estimation using PEST software (Doherty, J.E., 2005, PEST, model-independent parameter estimation–User manual: Watermark Numerical Computing, https://pesthomepage.org/documentation ). Parameter estimates were constrained with several types of observation data, including water levels, water level changes, vertical water level differences, drain discharges, change in drain discharges, river seepage, seepage from Lake Lowell, and change in seepage from Lake Lowell. Material properties from the hydrogeologic framework were also used to assign the minimum and maximum values of some parameters. A final parameter realization was reached that minimized residuals between the observed and modelled values for the various observation groups. Mean residuals for the observation groups were 15.4 feet (ft) for water levels, 0.2 ft for water level changes, 19.4 ft for vertical water level differences, −3.9 cubic feet per second (ft 3 /s) for drain discharges, 0.0 ft 3 /s for changes in drain discharge, 45.0 ft 3 /s for river seepage, −40.1 ft 3 /s for Lake Lowell seepage, and 126.3 ft 3 /s for changes in Lake Lowell seepage. The quality of the model’s fit to observations varied spatially, with notable areas of under- or over-simulation of water levels present to the northwest and southwest of Lake Lowell, in the foothills along the eastern model boundary, and near the City of Eagle. Trends were observed in the residuals of many of the observation groups, indicating that the model is missing or not fully reproducing some phenomena that are observed in the system. The TVGWFM can be used as a tool for water resource planning, for understanding the interactions of groundwater and surface water at a basin scale, and for facilitating conjunctive management, but may lack the precision needed for water rights administration at a local scale. Additional sources of uncertainty or limitations of the model are noted. The quantity and spatial distribution of canal seepage and infiltration of irrigation water recharge, the largest sources of recharge to the system, are unknown and approximated indirectly. There is poor understanding of how canal seepage and incidental recharge change as land is converted from agricultural (irrigated) to suburban (semi-irrigated). These uncertainties will affect any scenarios that investigate changes to land use or irrigation practices. Finally, the model has relatively high water-level residuals around and to the southwest of Lake Lowell and should not be used to estimate water level effects in that region. The model was built with multiple, broadly expressed objectives and did not optimize performance for specific uses. However, the model and the tools included in an associated data release provide ample flexibility to improve the model for future uses. Adjustments and improvements could be made by refining the model in an area of interest, collecting additional calibration data, applying more rigorous boundary conditions, or re-estimating model parameters to optimize model performance for a specific model forecast.

Idaho

Hydrogeologic characteristics of the St. Croix River basin, Minnesota and Wisconsin: Implications for the susceptibility of ground water to potential contamination

Population growth in the St. Croix River Basin in Minnesota and Wisconsin has intensified concerns of county resource managers and the National Park Service, which is charged with protecting the St. Croix National Scenic Riverway, about the potential for ground-water contamination in the basin. This report describes a previously developed method that was adapted to illustrate potential ground-water-contamination susceptibility in the St. Croix River Basin. The report also gives an estimate of ground-water-residence time and surface-water/ground-water interaction as related to natural attenuation and movement of contaminants in five tributary basins. A ground-water-contamination-susceptibility map was adapted from a state-wide map of Wisconsin to the St. Croix River Basin by use of well-driller construction records and regional maps of aquifer properties in Minnesota and Wisconsin. Measures of various subsurface properties were combined to generate a spatial index of susceptibility. The subjective index method developed for the State of Wisconsin by Schmidt (1987) 1 was not derived from analyses of water-quality data or physical processes. Nonetheless, it was adapted for this report to furnish a seamless map across state boundaries that would be familiar to many resource managers. Following this method, areas most susceptible to contamination appear to have coarse-grained sediments (sands or gravels) and shallow water tables or are underlain by carbonate-bedrock aquifers. The least susceptible areas appear to have fine-grained sediments and deep water tables. If an aquifer becomes contaminated, the ground-water-residence time can affect potential natural attenuation along the ground-water-flow path. Mean basin ground-water-residence times were computed for the Apple, Kettle, Kinnickinnic, Snake and Sunrise River Basins, which are tributary basins to the St. Croix Basin, by use of average aquifer properties of saturated thickness, porosity, and recharge rates. The Apple River Basin had the shortest mean ground-water-residence times (20–120 years), owing largely to the moderate saturated thickness and high recharge rate in the basin. The Kinnickinnic and Sunrise River Basins had the longest mean residence times (60–350 and 70–390 years, respectively) chiefly because of the relatively large saturated thickness of the basins. Owing to limitations of the residence-time calculations, actual ground-water-residence times will vary around the mean values within each basin and may range from days or weeks in karst carbonate aquifers to millennia in deep confined sandstone aquifers. Areas of relatively short residence time (less than the median residence time in each basin) were identified by use of ground-water-flow models for each of the five tributary basins. Results of simulations show that these areas, in which contaminants may have relatively less time for natural attenuation along the short flow paths, generally occur near streams and rivers where ground water discharges to the surface. Finally, the ground-water-flow models were used to simulate ground-water/surface-water interaction in the five tributary basins. Results of simulations show that some lakes and reservoirs leak surface water into the ground-water-flow system on their downgradient side, where the surface-water outflow has been restricted by a dam or a naturally constricted outlet. These locations are noteworthy because contaminated surface waters could potentially enter the ground-water-flow system at these locations.

Minnesota, Wisconsin

Hydrogeology, karst, and groundwater availability of Monroe County, West Virginia

Monroe County is in southeastern West Virginia, encompassing an area of 474 square miles. The area consists of karst and siliciclastic aquifers of Ordovician, Silurian, Devonian, and Mississippian age and is in parts of two physiographic provinces: the Valley and Ridge Province to the east of Peters Mountain, and the Appalachian Plateau Province to the west of Peters Mountain. This study was developed in response to inquiries from the Monroe County Commission requesting assessment of the water resources of the county to better understand the quantity of the county’s groundwater resources, for both current [2023] and future demand, and to provide information to support protection and management of the county’s valuable groundwater resources. Various products were developed for this study that provide knowledge with respect to water availability and contamination susceptibility of the karst aquifers within the county. U.S. Geological Survey (USGS) geologists conducted extensive geologic mapping in support of the project, producing (1) a countywide bedrock geologic map, (2) a countywide hydrogeologic map, and (3) a light detection and ranging (lidar)-derived countywide digital elevation model and associated sinkhole map. A significant part of this work was to map in detail the Greenbrier Group at the formation level, which prior to this study had only partially been completed. The report also includes (4) a description of the lithologic units identified as part of the geologic mapping process. U.S. Geological Survey hydrologists completed several additional products for the hydrology part of the effort, including development of (1) a countywide potentiometric surface (water-table) map, (2) a countywide base-flow stream assessment, (3) countywide water-budget estimates, (4) well log surveys for 15 wells to better understand subsurface controls on groundwater flow within the study area, (5) two groundwater tracer tests to better refine the groundwater divide from the northern and southern parts of the karst aquifer in Monroe County; and finally, based on all available data collected for the study including the potentiometric surface map, geologic map, current [2023] and legacy fluorometric groundwater tracer tests, and base-flow stream assessments, (6) groundwater-basin delineations were reassessed for principal groundwater basins within the Greenbrier aquifer. In Monroe County, four principal hydrogeologic settings produce large yields of water for residential, agricultural, and other uses. The most relied upon water-bearing zone with respect to current [2023] public water supply is from springs along Peters Mountain. These springs are derived from intervals of fractured sandstone and resultant alluvial deposits. Groundwater flows downslope through these permeable alluvial deposits and discharges at the contact with less permeable strata, such as the Reedsville Shale. The second most relied upon water-bearing zone in Monroe County is within the karstic Greenbrier Group aquifer, in which the basal Hillsdale Limestone overlies the less permeable Maccrady Shale. This geologic contact between the Hillsdale Limestone and Maccrady Shale is not only targeted as a source of water for agricultural supply but also is targeted as a source of water for residential supply. The third most relied upon water-bearing zone is composed of shallow perched aquifers within the Greenbrier Group. The discontinuous nature of these perched aquifers makes mapping their extent impossible, but they are related to permeable geologic strata, such as karstified limestones with solutionally enhanced permeability that overlies less permeable shale or chert bedrock. During geologic mapping of the county, several of these perched aquifers were documented in the Pickaway, Union, and Alderson Limestones. A fourth zone consists of springs from Ordovician carbonates at the base of Peters Mountain, which are influenced by sinking streams as well as upwelling along faults. In terms of water quantity, the most sustainable springs are those having deeper-sourced flows. Public supplies are a principal source of water used for residential and commercial supply in the region, accounting for 0.49 million gallons per day (Mgal/d) of fresh-water withdrawals (0.14 Mgal/d of groundwater and 0.35 Mgal/d of surface water) for residential and commercial use and serving 6,645 individuals (49.2 percent of the population). An estimated 6,861 people, (50.8 percent of the population) primarily rely on private wells or other unregulated sources, such as springs, and withdraw 0.55 Mgal/d of groundwater for their residential use. Public water supply in the region is primarily (71.4 percent) derived from springs and augmented by stream withdrawals (backup sources mainly during low-flow periods), with the remaining portion (28.6 percent) derived from groundwater withdrawals from wells. For rural residents, however, 100 percent of their withdrawals are derived from groundwater (wells or springs).

West Virginia

Novel observations of larval fire survival, feeding behavior, and host plant use in the regal fritillary, Speyeria idalia (Drury) (Nymphalidae)

Speyeria idalia is a prairie specialist that has experienced dramatic population declines throughout its range. Speyeria idalia is nearly extirpated from the eastern portion of its former range; however, populations within Kansas are relatively stable. We made several previously undescribed field observations of late-instar larvae and post-diapause female S. idalia in northeastern Kansas during 2014–2016. We report finding late-instar larvae at locations that were burned within weeks of detection. The observations of larvae shortly following a burn suggests that S . idalia larvae are capable of surviving fire and contradicts our current knowledge of this species. Additionally, we describe a feeding behavior characteristic of late-instar larvae. Larvae observed in the field and lab stripped leaves of host plants leaving only stems. This strip-style feeding behavior provided unique feeding evidence that was valuable to detecting the presence of larvae in the field. Finally, we documented larvae and post-diapause, egg depositing females using Viola sororia. The use of this relatively widespread and common plant by S. idalia populations in the Central Great Plains has only been implicitly documented but may have important conservation implications. These novel observations further our knowledge of the ecology of this imperiled species and provide timely information that may improve research and conservation management efforts directed toward S . idalia populations.

Journal of the Lepidopterists' Society

Significance of headwater streams and perennial springs in ecological monitoring in Shenandoah National Park

Shenandoah National Park has been monitoring water chemistry and benthic macroinvertebrates in stream ecosystems since 1979. These monitoring efforts were designed to assess the status and trends in stream condition associated with atmospheric deposition (acid rain) and changes in forest health due to gypsy moth infestations. The primary objective of the present research was to determine whether the current long-term macroinvertebrate and water-quality monitoring program in Shenandoah National Park was failing to capture important information on the status and trends in stream condition by not sufficiently representing smaller, headwater streams. The current benthic-macroinvertebrate and water-chemistry sampling designs do not include routine collection of data from streams with contributing watershed areas smaller than 100 hectares, even though these small streams represent the overwhelming proportion of total stream length in the park. In this study, we sampled headwater sites, including headwater stream reaches (contributing watershed area approximately 100 hectares (ha) and perennial springs, in the park for aquatic macroinvertebrates and water chemistry and compared the results with current and historical data collected at long-term ecological monitoring (LTEM) sites on larger streams routinely sampled as part of ongoing monitoring efforts. The larger purpose of the study was to inform ongoing efforts by park managers to evaluate the effectiveness and efficiency of the current aquatic monitoring program in light of other potential stressors (for example, climate change) and limited resources. Our results revealed several important findings that could influence management decisions regarding long-term monitoring of park streams. First, we found that biological indicators of stream condition at headwater sites and perennial springs generally were more indicative of lower habitat quality and were more spatially variable than those observed at sites on routinely monitored larger streams. We hypothesized that poorer stream condition observed in smaller streams was due to stream drying that occurs more frequently in headwater areas. We also found that biological and water-chemistry measures responded differently to landscape drivers. Variation in most biological endpoints was driven primarily by stream size and was only secondarily associated with bedrock geology. In contrast, water chemistry showed essentially the opposite pattern, with underlying geology explaining much of the variation and stream size being of secondary importance. Therefore, expanding the LTEM program to include headwater areas would yield substantially different biological information, whereas broad inferences regarding spatial patterns in water chemistry would probably not change. Although significant differences in community composition were observed among streams of different sizes, no taxa were unique to headwater sites. All taxa collected at the 45 headwater sites also had been collected at one or more LTEM sites during one or more years. This observation indicates that headwater sites in the park may be structured by biotic nestedness; consequently, focusing management efforts on preserving the species pool at the larger LTEM sites would likely result in the protection of most taxa parkwide. Finally, linkages (correlations) between water chemistry and biological measures of stream condition were signficantly stronger when assessed at the LTEM sites than when assessed at the springs or headwater sites, indicating that conditions at downstream sites may be better indicators of water-quality trends.

Virginia

Geology, geochemistry, and genesis of the Greens Creek massive sulfide deposit, Admiralty Island, southeastern Alaska

In 1996, a memorandum of understanding was signed by representatives of the U.S. Geological Survey and Kennecott Greens Creek Mining Company to initiate a cooperative applied research project focused on the Greens Creek massive sulfide deposit in southeastern Alaska. The goals of the project were consistent with the mandate of the U.S. Geological Survey Mineral Resources Program to maintain a leading role in national mineral deposits research and with the need of Kennecott Greens Creek Mining Company to further development of the Greens Creek deposit and similar deposits in Alaska and elsewhere. The memorandum enumerated four main research priorities: (1) characterization of protoliths for the wall rocks, and elucidation of their alteration histories, (2) determination of the ore mineralogy and paragenesis, including metal residences and metal zonation within the deposit, (3) determination of the ages of events important to ore formation using both geochronology and paleontology, and (4) development of computer models that would allow the deposit and its host rocks to be examined in detail in three dimensions. The work was carried out by numerous scientists of diverse expertise over a period of several years. The written results, which are contained in this Professional Paper, are presented by 21 authors: 13 from the U.S. Geological Survey, 4 from Kennecott Greens Creek Mining Company, 2 from academia, and 2 from consultants. The Greens Creek deposit (global resource of 24.2 million tons at an average grade of 13.9 percent zinc, 5.1 percent lead, 0.15 troy ounce per ton gold, and 19.2 troy ounces per ton silver at zero cutoff) formed in latest Triassic time during a brief period of rifting of the Alexander terrane. The deposit exhibits a range of syngenetic, diagenetic, and epigenetic features that are typical of volcanogenic (VMS), sedimentary exhalative (SEDEX), and Mississippi Valley-type (MVT) genetic models. In the earliest stages of rifting, formation of precious-metal-rich silica-barite-carbonate white ores began at low temperature in a shallow, subaqueous setting, probably a thin carbonate shelf on the flanks of the Alexander landmass. Epigenetic carbonate replacement textures in the footwall dolostones are overlain by stratiform silica-carbonate-barite-rich ores and indicate that early mineralization formed at and just beneath the paleo sea floor by mixing of a reduced, precious-metal-rich, base-metal-poor hydrothermal fluid with oxygenated seawater. As rifting intensified, the shelf was downfaulted and isolated as a graben. Isolation of the basin and onset of starved-basin shale sedimentation was concurrent with emplacement of mafic-ultramafic intrusives at shallow levels in the rift, resulting in an increasingly higher temperature and progressively more anoxic ore-forming environment. The formation of the main stage of massive sulfide ores began as the supply of bacterially reduced sulfur increased in the accumulating shales. As the main-stage mineralization intensified, shale sedimentation inundated the hydrothermal system, eventually forming a cap. Biogenic sulfate reduction supplied reduced sulfur to the base of the shales where mixing occurred with hot, base-metal-rich hydrothermal fluids. Ore deposition continued by destruction and epigenetic replacement of the early white ores in proximal areas and by inflation and diagenetic replacement of unlithified shale at the interface between the white ores and the base of the shale cap. Ore deposition waned as the shales became lithified and as the supply of bacterially reduced sulfur to the site of ore deposition ceased. The final stages of rifting resulted in the emplacement of mafic-ultramafic intrusive rocks into the Greens Creek system and extrusion of voluminous basaltic flows at the top of the Triassic section. Greenschist facies metamorphism during the Jurassic-Cretaceous accretion of the Alexander terrane to the continental margin resulted in recrystalli

Professional Paper

Resurvey of the Marble Canyon and Bridge Canyon dam sites in Grand Canyon National Park—Changes in sediment storage and evidence supporting the occurrence of bedrock incision through the mid-20th century

The Bureau of Reclamation developed an extensive plan for a network of dams, water tunnels, and hydropower plants in and around Grand Canyon, Arizona, in the 1940s through 1960s. The two largest of these planned dams were the Marble Canyon and Bridge Canyon Dams on the Colorado River. Though these dams were ultimately never built, Reclamation conducted extensive topographic, bathymetric, and subsurface exploration work at the sites proposed for these dams in the 1940s and 1950s. Resurveys of these dam sites were conducted between 1998 and 2021 to determine the changes in sediment storage at these dam sites caused by the upstream construction and operation of Glen Canyon Dam and by the recession of Lake Mead, the reservoir impounded by Hoover Dam. The resurveys of the Marble Canyon dam sites indicate that the post-1950s changes in sediment storage at these dam sites are broadly consistent with flux-based estimates of voluminous sand erosion from Marble Canyon since the 1963 closure of Glen Canyon Dam. These resurveys also suggest that the pre-dam longitudinal variation in sediment thickness over bedrock played a key role in determining the locations of the sand erosion induced by Glen Canyon Dam; more sand eroded from locations where more sand was present in the 1950s. The resurvey of the Bridge Canyon dam sites indicates that the Colorado River’s incision of the Lake Mead delta is regulated both by bed-sediment grain size and downstream hydraulic controls. Finally, analyses of bed-sediment thickness and sedimentological data at the dam sites, and observations of bed scour and gravel transport, suggest that sufficient bedrock was exposed to allow bedrock incision during commonly recurring pre-dam snowmelt floods that entrained small boulders into transport.

Arizona

2005 annual progress report: Elk and bison grazing ecology in the Great Sand Dunes complex of lands

In 2000 the U.S. Congress authorized the expansion of the former Great Sand Dunes National Monument by establishing a new Great Sand Dunes National Park and Preserve in its place, and establishing the Baca National Wildlife Refuge. The establishment of Great Sand Dunes National Park and Preserve and the new Baca National Wildlife Refuge in the San Luis Valley (SLV), Colorado was one of the most significant land conservation actions in the western U.S. in recent years. The action was a result of cooperation between the National Park Service (NPS), U.S. Fish and Wildlife Service (USFWS), Bureau of Land Management (BLM), U.S. Forest Service (USDA-FS), and The Nature Conservancy (TNC). The new national park, when fully implemented, will consist of 107,265 acres, the new national preserve 41,872 acres, and the new national wildlife refuge (USFWS lands) 92,180 acres (fig. 1). The area encompassed by this designation protects a number of natural wonders and features including a unique ecosystem of natural sand dunes, the entire watershed of surface and groundwaters that are necessary to preserve and recharge the dunes and adjacent wetlands, a unique stunted forest, and other valuable riparian vegetation communities that support a host of associated wildlife and bird species. When the National Park was initially established, there were concerns about overconcentrations and impacts on native plant communities of the unhunted segments of a large and possibly growing elk (Cervus elaphus) population. This led to the designation of the Preserve as a compromise solution, where the elk could be harvested. The Preserve Unit, however, will not address all the ungulate management challenges. In order to reduce the current elk population, harvests of elk may need to be aggressive. But aggressive special hunts of elk to achieve population reductions can result in elk avoidance of certain areas or elk seeking refuge in areas where they cannot be hunted, while removals of whole herd segments and abandonment or alterations of migration routes can occur (Smith and Robbins, 1994; Boyce and others, 1991). Elk may seek refuge from hunting in the newly expanded Park Unit and TNC lands where they might overconcentrate and impact unique vegetation communities. In these sites of refugia, or preferred loafing sites, elk and bison could accelerate a decline in woody riparian shrubs and trees. This decline may also be due to changes in hydrology, climatic, or dunal processes, but ungulate herbivory might exacerbate the effects of those processes. To address the questions and needs of local resource managers, a multi-agency research project was initiated in 2005 to study the ecology, forage relations, and habitat relations of elk and bison in the Great Sand Dunes–Sangre de Cristo–Baca complex of lands. Meetings and discussions of what this research should include were started in 2001 with representatives from NPS, USFWS, TNC, the Colorado Division of Wildlife (CDOW), and USDA-FS/BLM. The final study plan was successfully funded in 2004 with research scheduled to start in 2005. The research was designed to encompass three major study elements: (1) animal movements and population dynamics, (2) vegetation and nutrient effects from ungulate herbivory, and (3) development of ecological models, using empirical data collected from the first two components, that will include estimates of elk carrying capacity and management scenarios for resource managers.

Colorado

GIS-based identification of areas that have resource potential for critical minerals in six selected groups of deposit types in Alaska

Alaska has considerable potential for undiscovered mineral resources. This report evaluates potential for undiscovered critical minerals in Alaska. Critical minerals are those for which the United States imports more than half of its total supply and which are largely derived from nations that cannot be considered reliable trading partners. In this report, estimated resource potential and certainty for the state of Alaska are analyzed and mapped for the following six selected mineral deposit groups that may contain one or more critical minerals: (1) rare earth elements-thorium-yttrium-niobium(-uranium-zirconium) [REE-Th-Y-Nb(-U-Zr)] deposits associated with peralkaline to carbonatitic igneous intrusive rocks; (2) placer and paleoplacer gold (Au) deposits that in some places might also produce platinum group elements (PGE), chromium (Cr), tin (Sn), tungsten (W), silver (Ag), or titanium (Ti); (3) platinum group elements(-cobalt-chromium-nickel-titanium-vanadium) [PGE(-Co-Cr-Ni-Ti-V)] deposits associated with mafic to ultramafic intrusive rocks; (4) carbonate-hosted copper(-cobalt-silver-germanium-gallium) [Cu(-Co-Ag-Ge-Ga)] deposits; (5) sandstone-hosted uranium(-vanadium-copper) [U(-V-Cu)] deposits; and (6) tin-tungsten-molybdenum(-tantalum-indium-fluorspar) [Sn-W-Mo(-Ta-In-fluorspar)] deposits associated with specialized granites. This study used a data-driven, geographic information system (GIS)-implemented method to identify areas that have mineral resource potential in Alaska. This method systematically and simultaneously analyzes geoscience data from multiple geospatially referenced datasets and uses individual subwatersheds (12-digit hydrologic units) as the spatial unit of classification. The final map output uses a red, yellow, green, and gray color scheme to portray estimated relative potential (High, Medium, Low, Unknown) for each of the six groups of mineral deposit types, and it indicates the relative certainty (High, Medium, Low) of that estimate for each 12-digit hydrologic unit through color shading. Accompanying tables describe the data layers employed to score favorability for the presence of each mineral deposit group, the values assigned for specific analysis parameters, and the relative weighting of each data layer that contributes to estimated measures of potential and certainty . Core datasets used include the Alaska Geochemical Database, Version 2.0 (AGDB2); the Alaska Division of Geological & Geophysical Surveys (ADGGS) web-based geochemical database; the digital “Geologic Map of Alaska;” the Alaska Resource Data File (ARDF); and aerial gamma-ray surveys flown as part of the National Uranium Resource Evaluation (NURE) Program by the U.S. Department of Energy. Maps accompanying this report illustrate the scores for estimated mineral resource potential for the six deposit groups for the state of Alaska. Areas that have known potential, as well as new areas that were not previously known to have potential, for the targeted minerals and deposit groups are identified and described. Numerous areas in Alaska, some of them large, have high potential for one or more of the selected groups of deposit types within Alaska. Contributors Matthew Granitto, Timothy S. Hayes, James V. Jones, III, Susan M. Karl, Keith A. Labay, Jeffrey L. Mauk, Jeanine M. Schmidt, Nora B. Shew, Erin Todd, Bronwen Wang, Melanie B. Werdon, and Douglas B. Yager

Alaska

The United States Geological Survey: 1879-1989

The United States Geological Survey was established on March 3, 1879, just a few hours before the mandatory close of the final session of the 45th Congress, when President Rutherford B. Hayes signed the bill appropriating money for sundry civil expenses of the Federal Government for the fiscal year beginning July 1, 1879. The sundry civil expenses bill included a brief section establishing a new agency, the United States Geological Survey, placing it in the Department of the Interior, and charging it with a unique combination of responsibilities: 'classification of the public lands, and examination of the geological structure, mineral resources, and products of the national domain.' The legislation stemmed from a report of the National Academy of Sciences, which in June 1878 had been asked by Congress to provide a plan for surveying the Territories of the United States that would secure the best possible results at the least possible cost. Its roots, however, went far back into the Nation's history. The first duty enjoined upon the Geological Survey by the Congress, the classification of the public lands, originated in the Land Ordinance of 1785. The original public lands were the lands west of the Allegheny Mountains claimed by some of the colonies, which became a source of contention in writing the Articles of Confederation until 1781 when the States agreed to cede their western lands to Congress. The extent of the public lands was enormously increased by the Louisiana Purchase in 1803 and later territorial acquisitions. At the beginning of Confederation, the decision was made not to hold the public lands as a capital asset, but to dispose of them for revenue and to encourage settlement. The Land Ordinance of 1785 provided the method of surveying and a plan for disposal of the lands, but also reserved 'one-third part of all gold, silver, lead, and copper mines to be sold or otherwise disposed of, as Congress shall thereafter direct,' thus implicitly requiring classification of the lands into mineral and nonmineral. Mapping of the public lands was begun under the direction of the Surveyor-General, but no special provision was made for classification of the public lands, and it thus became the responsibility of the surveyor. There was,of course, no thought in 1785 or for many years thereafter of employing geologists to make the classification of the mineral lands, for geology was then only in its infancy.

Circular

Crotalus oreganus lutosus (Great Basin Rattlesnake). Defensive behavior/head hiding and tongue display

A variety of snake species employ head hiding as defensive behavior (Greene 1973. J. Herpetol. 7:143–161), but such behavior seems to be rarely observed in rattlesnakes. A recent report (Medica 2009. Herpetol. Rev. 40:95–97) presented observations on Crotalus scutulatus scutulatus and cited cases of similar behavior involving C. atrox , C. ruber , and C. viridis viridis . Here we report an encounter on 28 May 2016 with a C. oreganus lutosus (total length ca. 900 mm) observed 16 km S of Fallon, Churchill County, Nevada, USA. At 0925 h this snake was spotted by AKJ near the edge of the pavement on U.S. Hwy. 95. When approached, the snake presented typical active defensive behavior, with the rattle elevated (although not rattling), striking occasionally, and frequently flicking its conspicuous lavender and black tongue. As AKJ continued to take pictures the snake coiled more tightly (Fig. 1A), maintaining its rattle concealed and finally hiding its head beneath its coils (Fig. 1B). After photographing this pose AKJ made his first physical contact with the snake by scooping it up with a small shovel and placing it off the road. At this point the snake recovered its feisty attitude and would not retreat despite attempts to move it further away from the road. The snake did not rattle during the entire encounter (which lasted less than 5 min). Head hiding as a “last resort” in response to a persistent unavoidable threat may be more widespread within the family Crotalidae than previously documented.

Nevada

Local-area-enhanced, 2.5-meter resolution natural-color and color-infrared satellite-image mosaics of the Takhar mineral district in Afghanistan: Chapter Q in Local-area-enhanced, high-resolution natural-color and color-infrared satellite-image mosaics of mineral districts in Afghanistan

The U.S. Geological Survey (USGS), in cooperation with the U.S. Department of Defense Task Force for Business and Stability Operations, prepared databases for mineral-resource target areas in Afghanistan. The purpose of the databases is to (1) provide useful data to ground-survey crews for use in performing detailed assessments of the areas and (2) provide useful information to private investors who are considering investment in a particular area for development of its natural resources. The set of satellite-image mosaics provided in this Data Series (DS) is one such database. Although airborne digital color-infrared imagery was acquired for parts of Afghanistan in 2006, the image data have radiometric variations that preclude their use in creating a consistent image mosaic for geologic analysis. Consequently, image mosaics were created using ALOS (Advanced Land Observation Satellite; renamed Daichi) satellite images, whose radiometry has been well determined (Saunier, 2007a,b). This part of the DS consists of the locally enhanced ALOS image mosaics for the Takhar mineral district, which has industrial evaporite deposits. ALOS was launched on January 24, 2006, and provides multispectral images from the AVNIR (Advanced Visible and Near-Infrared Radiometer) sensor in blue (420–500 nanometer, nm), green (520–600 nm), red (610–690 nm), and near-infrared (760–890 nm) wavelength bands with an 8-bit dynamic range and a 10-meter (m) ground resolution. The satellite also provides a panchromatic band image from the PRISM (Panchromatic Remote-sensing Instrument for Stereo Mapping) sensor (520–770 nm) with the same dynamic range but a 2.5-m ground resolution. The image products in this DS incorporate copyrighted data provided by the Japan Aerospace Exploration Agency (©JAXA, 2008), but the image processing has altered the original pixel structure and all image values of the JAXA ALOS data, such that original image values cannot be recreated from this DS. As such, the DS products match JAXA criteria for value added products, which are not copyrighted, according to the ALOS end-user license agreement. The selection criteria for the satellite imagery used in our mosaics were images having (1) the highest solar-elevation angles (near summer solstice) and (2) the least cloud, cloud-shadow, and snow cover. The multispectral and panchromatic data were orthorectified with ALOS satellite ephemeris data, a process which is not as accurate as orthorectification using digital elevation models (DEMs); however, the ALOS processing center did not have a precise DEM. As a result, the multispectral and panchromatic image pairs were generally not well registered to the surface and not coregistered well enough to perform resolution enhancement on the multispectral data. Therefore, it was necessary to (1) register the 10-m AVNIR multispectral imagery to a well-controlled Landsat image base, (2) mosaic the individual multispectral images into a single image of the entire area of interest, (3) register each panchromatic image to the registered multispectral image base, and (4) mosaic the individual panchromatic images into a single image of the entire area of interest. The two image-registration steps were facilitated using an automated control-point algorithm developed by the USGS that allows image coregistration to within one picture element. Before rectification, the multispectral and panchromatic images were converted to radiance values and then to relative-reflectance values using the methods described in Davis (2006). Mosaicking the multispectral or panchromatic images started with the image with the highest sun-elevation angle and the least atmospheric scattering, which was treated as the standard image. The band-reflectance values of all other multispectral or panchromatic images within the area were sequentially adjusted to that of the standard image by determining band-reflectance correspondence between overlapping images using linear least-squares analysis. The resolution of the multispectral image mosaic was then increased to that of the panchromatic image mosaic using the SPARKLE logic, which is described in Davis (2006). Each of the four-band images within the resolution-enhanced image mosaic was individually subjected to a local-area histogram stretch algorithm (described in Davis, 2007), which stretches each band’s picture element based on the digital values of all picture elements within a 315-m radius. The final databases, which are provided in this DS, are three-band, color-composite images of the local-area-enhanced, natural-color data (the blue, green, and red wavelength bands) and color-infrared data (the green, red, and near-infrared wavelength bands). All image data were initially projected and maintained in Universal Transverse Mercator (UTM) map projection using the target area’s local zone (42 for Takhar) and the WGS84 datum. The final image mosaics for the Takhar area are provided as embedded geotiff images, which can be read and used by most geographic information system (GIS) and image-processing software. The tiff world files (tfw) are provided, even though they are generally not needed for most software to read an embedded geotiff image.

Takhar Mineral District

Plagioclase populations and zoning in dacite of the 2004-2005 Mount St. Helens eruption: Constraints for magma origin and dynamics

We investigated plagioclase phenocrysts in dacite of the 2004-5 eruption of Mount St. Helens to gain insights into the magmatic processes of the current eruption, which is characterized by prolonged, nearly solid-state extrusion, low gas emission, and shallow seismicity. In addition, we investigated plagioclase of 1980-86 dacite. Light and Nomarski microscopy were used to texturally characterize plagioclase crystals. Electron microprobe analyses measured their compositions. We systematically mapped and categorized all plagioclase phenocrysts in a preselected area according to the following criteria: (1) occurrence of zones of acicular orthopyroxene inclusions, (2) presence of dissolution surface(s), and (3) spatial association of 1 and 2. Phenocrysts fall into three main categories; one category contains four subcategories. The range of anorthite (An) content in 2004-5 plagioclase is about An 57-35 during the last 30-40 percent crystallization of plagioclase phenocrysts. Select microphenocrysts (10-50 μm) range from An 30 to An 42 . Anorthite content is lowest near outermost rims of phenocrysts, but zonation patterns between interior and rim indicate variable trends that correlate with textural features. Crystals without dissolution surfaces (about 14 percent of total) show steadily decreasing An content outward to the crystal rim (outer ~80 μm). All other crystals are banded as a consequence of dissolution; dissolution surfaces are band boundaries. Such crystals display normal outward An zoning within a single band that, following dissolution, is then overgrown abruptly by high-An material of the next band. Swarms of acicular orthopyroxene inclusions in plagioclase are characteristic of 2004-5 dacite. They occur mostly inward of dissolution surfaces, where band composition reaches lowest An content. The relative proportions of the three crystal types are distinctly different between 2004-5 dacite and 1980s dome dacite. We propose that crystals with no dissolution surfaces are those that were supplied last to the shallow reservoir, whereas plagioclase with increasingly more complex zoning patterns (that is, the number of zoned bands bounded by dissolution surfaces) result from prolonged residency and evolution in the reservoir. We propose that banding and An zoning across multiple bands are primarily a response to thermally induced fluctuations in crystallinity of the magma in combination with recharge; a lesser role is ascribed to cycling crystals through pressure gradients. Crystals without dissolution surfaces, in contrast, could have grown only in response to steady(?) decompression. Some heating-cooling cycles probably postdate the final eruption in 1986. They resulted from small recharge events that supplied new crystals that then experienced resorption-growth cycles. We suggest that magmatic events shortly prior to the current eruption, recorded in the outermost zones of plagioclase phenocrysts, began with the incorporation of acicular orthopyroxene, followed by last resorption, and concluded with crystallization of euhedral rims. Finally, we propose that 2004-5 dacite is composed mostly of dacite magma that remained after 1986 and underwent subsequent magmatic evolution but, more importantly, contains a component of new dacite from deeper in the magmatic system, which may have triggered the new eruption.

Washington

Approach to data exchange: the spatial data transfer standard

Significant developments have taken place in the disciplines of cartography and geography in recent years with the advent of computer hardware and software that manipulate and process digital cartographic and geographic data more efficiently. The availability of inexpensive and powerful hardware and software systems offers the capability of displaying and analyzing spatial data to a growing number of users. As a result, developing and using existing digital cartographic databases are becoming very popular. However, the absence of uniform standards for the transfer of digital spatial data is hindering the exchange of data and increasing costs. Several agencies of the U.S. government and the academic community have been working hard over the last few years to develop a spatial data transfer standard that includes definitions of standard terminology, a spatial data transfer specification, recommendations on reporting digital cartographic data quality, and standard topographic and hydrographic entity terms and definitions. This proposed standard was published in the January 1988 issue of The American Cartographer. Efforts to test and promote this standard were coordinated by the U.S. Geological Survey. A Technical Review Board was appointed with representatives from the U.S. government, the private sector, and the academic community to complete the standard for submittal to the National Institute of Standards and Technology for approval as a Federal Information Processing Standard. The proposed standard was submitted in February 1992 for final approval.

Conference Paper

Movements of translocated wolves in Minnesota

From Feb. 1975-May 1978, 104 wolves ( Canis lupus ) captured at or near Minnesota farms where depredations on livestock had been reported were translocated northward and eastward for 50-317 km into extensive forests; 3 others were released westward. Nine wolves were translocated twice, and 1 three times. Information on movements was obtained by radiotracking 17 wolves and by recovery of 16 others. All radio-tagged wolves left the release areas; adults left quickly, but pups generally remained longer, behaved less decisively, and settled nearby temporarily. Wolves released together did not remain together. Initial travel of most radio-tagged adults was between south and west, the general direction to their original location. Final directions were primarily to the west and northwest, due in part to physiographic barriers. Eight adults homed, 1 twice, to capture areas that were < 64 km from release sites. Nonhoming wolves were radiotracked and/or recovered 32-351 km from their capture sites and 23-302 km from their release sites. Translocation was largely unsuccessful at keeping problem wolves out of lifestock production areas. The problem of initial travel away from release sites in wolf reintroductions probably could be minimized by transport and release of 6-9 mo. old wolves.

Minnesota

Density and velocity relationships for digital sonic and density logs from coastal Washington and laboratory measurements of Olympic Peninsula mafic rocks and greywackes

Three-dimensional velocity models for the basins along the coast of Washington and in Puget Lowland provide a means for better understanding the lateral variations in strong ground motions recorded there. We have compiled 16 sonic and 18 density logs from 22 oil test wells to help us determine the geometry and physical properties of the Cenozoic basins along coastal Washington. The depth ranges sampled by the test-well logs fall between 0.3 and 2.1 km. These well logs sample Quaternary to middle Eocene sedimentary rocks of the Quinault Formation, Montesano Formation, and Hoh rock assemblage. Most (18 or 82%) of the wells are from Grays Harbor County, and many of these are from the Ocean City area. These Grays Harbor County wells sample the Quinault Formation, Montesano Formation, and frequently bottom in the Hoh rock assemblage. These wells show that the sonic velocity and density normally increase significantly across the contacts between the Quinault or the Montesano Formations and the Hoh rock assemblage. Reflection coefficients calculated for vertically traveling compressional waves from the average velocities and densities for these units suggest that the top of the Hoh rock assemblage is a strong reflector of downward-propagating seismic waves: these reflection coefficients lie between 11 and 20%. Thus, this boundary may reflect seismic energy upward and trap a substantial portion of the seismic energy generated by future earthquakes within the Miocene and younger sedimentary basins found along the Washington coast. Three wells from Jefferson County provide data for the Hoh rock assemblage for the entire length of the logs. One well (Eastern Petroleum Sniffer Forks #1), from the Forks area in Clallam County, also exclusively samples the Hoh rock assemblage. This report presents the locations, elevations, depths, stratigraphic, and other information for all the oil test wells, and provides plots showing the density and sonic velocities as a function of depth for each well log. We also present two-way traveltimes for 15 of the wells calculated from the sonic velocities. Average velocities and densities for the wells having both logs can be reasonably well related using a modified Gardner’s rule, with p=1825v(1/4), where p is the density (in kg/m3) and v is the sonic velocity (in km/s). In contrast, a similar analysis of published well logs from Puget Lowland is best matched by a Gardner’s rule of p=1730v(1/4), close to the p=1740v(1/4) proposed by Gardner et al. (1974). Finally, we present laboratory measurements of compressional-wave velocity, shear-wave velocity, and density for 11 greywackes and 29 mafic rocks from the Olympic Peninsula and Puget Lowland. These units have significance for earthquake-hazard investigations in Puget Lowland as they dip eastward beneath the Lowland, forming the “bedrock” beneath much of the lowland. Average Vp/Vs ratios for the mafic rocks, mainly Crescent Formation volcanics, lie between 1.81 and 1.86. Average Vp/Vs ratios for the greywackes from the accretionary core complex in the Olympic Peninsula show greater scatter but lie between 1.77 and 1.88. Both the Olympic Peninsula mafic rocks and greywackes have lower shear-wave velocities than would be expected for a Poisson solid (Vp/Vs=1.732). Although the P-wave velocities and densities in the greywackes can be related by a Gardner’s rule of p=1720v(1/4), close to the p=1740v(1/4) proposed by Gardner et al. (1974), the velocities and densities of the mafic rocks are best related by a Gardner’s rule of p=1840v(1/4). Thus, the density/velocity relations are similar for the Puget Lowland well logs and greywackes from the Olympic Peninsula. Density/velocity relations are similar for the Washington coastal well logs and mafic rocks from the Olympic Peninsula, but differ from those of the Puget Lowland well logs and greywackes from the Olympic Peninsula.

Washington

Methodology for estimating times of remediation associated with monitored natural attenuation

Natural attenuation processes combine to disperse, immobilize, and biologically transform anthropogenic contaminants, such as petroleum hydrocarbons and chlorinated ethenes, in ground-water systems. The time required for these processes to lower contaminant concentrations to levels protective of human health and the environment, however, varies widely between different hydrologic systems, different chemical contaminants, and varying amounts of contaminants. This report outlines a method for estimating timeframes required for natural attenuation processes, such as dispersion, sorption, and biodegradation, to lower contaminant concentrations and mass to predetermined regulatory goals in groundwater systems. The time-of-remediation (TOR) problem described in this report is formulated as three interactive components: (1) estimating the length of a contaminant plume once it has achieved a steady-state configuration from a source area of constant contaminant concentration, (2) estimating the time required for a plume to shrink to a smaller, regulatoryacceptable configuration when source-area contaminant concentrations are lowered by engineered methods, and (3) estimating the time needed for nonaqueous phase liquid (NAPL) contaminants to dissolve, disperse, and biodegrade below predetermined levels in contaminant source areas. This conceptualization was used to develop Natural Attenuation Software (NAS), an interactive computer aquifers. NAS was designed as a screening tool and requires the input of detailed site information about hydrogeology, redox conditions, and the distribution of contaminants. Because NAS is based on numerous simplifications of hydrologic, microbial, and geochemical processes, the program may introduce unacceptable errors for highly heterogeneous hydrologic systems. In such cases, application of the TOR framework outlined in this report may require more detailed, site-specific digital modeling. The NAS software may be downloaded from the Web site http://www.cee.vt.edu/NAS/ Application of NAS illustrates several general characteristics shared by all TOR problems. First, the distance of stabilization of a contaminant plume is strongly dependent on the natural attenuation capacity of particular ground-water systems. The time that it takes a plume to reach a steady-state configuration, however, is independent of natural attenuation capacity. Rather, the time of stabilization is most strongly affected by the sorptive capacity of the aquifer, which is dependent on the organic matter content of the aquifer sediments, as well as the sorptive properties of individual contaminants. As a general rule, a high sorptive capacity retards a plume.s growth or shrinkage, and increases the time of stabilization. Finally, the time of NAPL dissolution depends largely on NAPL mass, composition, geometry, and hydrologic factors, such as ground-water flow rates. An example TOR analysis for petroleum hydrocarbon NAPL was performed for the Laurel Bay site in South Carolina. About 500 to 1,000 pounds of gasoline leaked into the aquifer at this site in 1991, and the NAS simulations suggested that TOR would be on the order of 10 years for soluble and poorly sorbed compounds, such as benzene and methyl tertiary-butyl ether (MTBE). Conversely, TOR would be on the order of 40 years for less soluble, more strongly sorbed compounds, such as toluene, ethylbenzene, and xylenes (TEX). These TOR estimates are roughly consistent with contaminant concentrations observed over 10 years of monitoring at this site where benzene and MTBE concentrations were observed to decrease rapidly and are approaching regulatory maximum concentration limits, whereas toluene, ethylbenzene, and xylene concentrations decreased at a slower rate and have remained relatively high. An example TOR analysis for petroleum hydrocarbon NAPL was performed for the Laurel Bay site in South Carolina.

South Carolina