Search USGS⌕ Search

SEARCH · Search USGS

Results for “Field Guides”

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.

626 records · Page 35Linked to original sources

The Valmy thrust sheet: A regional structure formed during the protracted assembly of the Roberts Mountains allochthon, Nevada, USA

A synthesis of field, biostratigraphic, detrital zircon geochronologic, and remote sensing data across north-central Nevada, United States, defines a thick, regionally extensive sheet of Middle–Upper Ordovician Valmy Formation quartzite that structurally overlies deformed early Paleozoic units of the Roberts Mountains allochthon. Late Paleozoic regional unconformities that record tectonic disruptions have been recognized in the foreland of central and eastern Nevada and locally within the Roberts Mountains allochthon; these identify multiple, regional tectonic events between the Devonian–Mississippian initiation of the Antler orogeny and the Permian–Triassic Sonoma orogeny. However, few studies have documented the regional kinematic history of the Robert Mountains allochthon sensu stricto. In the Independence Mountains of northern Nevada, emplacement of the Roberts Mountains allochthon is restricted to the Mississippian. In the Tuscarora Mountains, the range west and southwest of the Independence Mountains, several deformation events have been identified, and emplacement of the thrust sheet containing the Valmy Formation is restricted to the Late Pennsylvanian–Early Permian. These structural and temporal relations, reflected in the Antler foreland basin adjacent to the Roberts Mountains allochthon and overlap sequences, suggest that the Roberts Mountains allochthon is a composite stratigraphic terrane assembled along the Cordilleran margin during two or more late Paleozoic contractional events. Valmy Formation deposits likely represent the development of coalescing submarine fans below or within bypass channels in a deep slope or rise environment. Petrographic characteristics, biostratigraphy, and detrital zircon U-Pb age populations of the Valmy Formation link it to coeval slope and rise turbidites of the Vinini Formation and shelfal Eureka Quartzite; Valmy Formation detrital zircon age populations are dissimilar to the rift-to-drift facies of the Neoproterozoic–Cambrian Prospect Mountain Quartzite. Throughout north-central Nevada, the Valmy Formation is in fault contact with units of the Roberts Mountains allochthon, including the Devonian–Mississippian Slaven Chert, Silurian–Devonian Elder Sandstone, and Cambrian(?)–Ordovician Vinini Formation, which were deformed prior to, or during, emplacement of the thrust sheet containing Valmy Formation quartzite. Our mapping and data synthesis, guided by regional quartz maps based on remote sensing (Advanced Spaceborne Thermal Emission and Reflection Radiometer [ASTER]) data, delineate similar structural relationships discontinuously for >200 km along strike of the Roberts Mountains allochthon. Exploration for concealed gold deposits within reach of drilling requires knowledge of the relative thicknesses of the Roberts Mountains allochthon and the Valmy Formation. Overall thicknesses of deformed Roberts Mountains allochthon units between the Valmy Formation and underlying carbonate rocks, which host large, world-class Carlin-type gold deposits, vary by hundreds of meters, but are generally less than 700 m in three of the areas studied here. Recognition of windows through and klippen of the Roberts Mountains allochthon is essential for identification of areas where deposits may be at or near the surface. Correspondingly, most ongoing exploration for Carlin-type gold deposits subjacent to the Roberts Mountains allochthon targets concealed deposits. The model proposed in this study is applicable to determining depth to rocks prospective for undiscovered deposits.

Nevada↗

Simulated water budgets and ground-water/surface-water interactions in Bushkill and parts of Monocacy Creek watersheds, Northampton County, Pennsylvania: A preliminary study with identification of data needs

This report, prepared in cooperation with the Department of Environmental Protection, Office of Mineral Resources Management, provides a preliminary analysis of water budgets and generalized ground-water/surface-water interactions for Bushkill and parts of Monocacy Creek watersheds in Northampton County, Pa., by use of a ground-water flow model. Bushkill Creek watershed was selected for study because it has areas of rapid growth, ground-water withdrawals from a quarry, and proposed stream-channel modifications, all of which have the potential for altering ground-water budgets and the interaction between ground water and streams. Preliminary 2-dimensional, steady-state simulations of ground-water flow by the use of MODFLOW are presented to show the status of work through September 2005 and help guide ongoing data collection in Bushkill Creek watershed. Simulations were conducted for (1) predevelopment conditions, (2) a water table lowered for quarry operations, and (3) anthropogenic changes in hydraulic conductivity of the streambed and aquifer. Preliminary results indicated under predevelopment conditions, the divide between the Bushkill and Monocacy Creek ground-water basins may not have been coincident with the topographic divide and as much as 14 percent of the ground-water discharge to Bushkill Creek may have originated from recharge in the Monocacy Creek watershed. For simulated predevelopment conditions, Schoeneck Creek and parts of Monocacy Creek were dry, but Bushkill Creek was gaining throughout all reaches. Simulated lowering of the deepest quarry sump to an altitude of 147 feet for quarry operations caused ground-water recharge and streamflow leakage to be diverted to the quarry throughout about 14 square miles and caused reaches of Bushkill and Little Bushkill Creeks to change from gaining to losing streams. Lowering the deepest quarry sump to an altitude of 100 feet caused simulated ground-water discharge to the quarry to increase about 4 cubic feet per second. Raising the deepest sump to an altitude of 200 feet caused the simulated discharge to the quarry to decrease about 14 cubic feet per second.Decreasing the hydraulic conductivity of the streambed of Bushkill Creek in the reach of large losses of flow caused simulated ground-water levels to decline and ground-water discharge to a quarry to decrease from 74 to 45 cubic feet per second. Decreasing the hydraulic conductivity of a hypothesized highly transmissive zone with a plug of relatively impermeable material caused ground-water levels to increase east of the plug and decline west of the plug, and decreased the discharge to a quarry from 74 to 53 cubic feet per second. Preliminary results of the study have significant limitations, which need to be recognized by the user. The results demonstrated the usefulness of ground-water modeling with available data sets, but as more data become available through field studies, a more complete evaluation could be conducted of the preliminary assumptions in the conceptual model, model sensitivity, and effects of boundary conditions. Additional streamflow and ground-water-level measurements would be needed to better quantify recharge and aquifer properties, particularly the anisotropy of carbonate rocks. Measurements of streamflow losses at average, steady-state hydrologic conditions could provide a more accurate estimate of ground-water recharge from this source, which directly affects water budgets and contributing areas simulated by the model.

Pennsylvania↗

Vessel surveys for marine birds in the Gulf of America (Gulf of Mexico): Protocol manual (2017–2025)

This manual provides contextual background and detailed instructions for survey procedures that were used to collect, enter, archive, and analyze digitally-entered field data collected for marine birds during three large-vessel survey projects in the northern Gulf of America (also Gulf of Mexico, hereafter “Gulf”). Projects using this protocol included the Gulf of Mexico Marine Assessment Program for Protected Species (GoMMAPPS: 2017-2018), the Vessel Surveys for Abundance and Distribution of Marine Mammals and Seabirds (VSAD: 2023-2024), and certain repositioning legs of the Atlantic Marine Assessment Program for Protected Species (AMAPPS: 2021, 2023). Based substantially on Ballance and Force (2011), the protocol used in the Gulf is based on a standardized 300 m wide strip transect using software architecture available in program SeeBird (v.4.3.7), a data collection program developed by Robert Holland at the National Oceanic and Atmospheric Administration (NOAA) Southwest Fisheries Science Center, La Jolla, CA. Survey protocols relied on the digital entering of standardized metadata linked to sightings of seabirds (e.g., species, number, location, distance bin, association[s], behavior, flight direction, flight height, angle from ship, age, sex, plumage, and optional comments) as well as other protected species. SeeBird has additional modules that may be used to input data on mixed-species seabird flocks and flying fish (only the latter was used in the Gulf). SeeBird includes “user-driven” commands issued through the keyboard, touch screen, and/or mouse, many of which are aided by drop-down menus that prompt the user options that are available for input. In addition, the program features “event–driven” inputs updated at intervals selected by the user, and that also record exact time and location of such edits. We describe how SeeBird was adapted in the Gulf both to run directly via serial connections from the ship’s navigation system and to operate independently from the ship through a satellite-derived GPS feed. This user manual also gives information about: procedures used to document flying fish abundance on Gulf surveys; how to use the Beaufort scale for standardizing wind force and sea state environmental conditions; how to make and use a hand-held rangefinder for delimiting strip transect bin distances; additional functionalities in SeeBird for which user guides are not available; file-naming, file-archiving, and file location conventions used on Gulf marine bird surveys; protocols used in the Gulf for describing relative abundance of the macro-alga Sargassum as an environmental covariate; the standardized four-letter taxa codes (summarized in a SpeciesCode.txt file) used in the Gulf for denoting species, other taxonomic, and event categories as inputs to SeeBird; an Excel example spreadsheet that has successfully imported SeeBird’s ASCII text file; procedural steps to import, edit, and archive into Excel ASCII text files from SeeBird; initialization instructions for starting up program SeeBird on a Panasonic ToughBook™ that relies on integrated, satellite-derived GPS, and for writing “S” (seabirds + other taxa) and “E” (flying fish) ASCII files to a folder or subfolder other than the default root location on the laptop that contains program SeeBird; how to trouble-shoot miscellaneous start-up problems with program SeeBird; how to add new four-letter codes (into a SpeciesCode.txt file) when additional taxa or objects of interest are encountered. The primary objective of this manual is to archive a permanent, replicable history of survey methodology to facilitate continuity with comparable future projects that may be conducted in the Gulf. A secondary objective is to describe why and how program SeeBird was applied and adapted specifically for operating in Gulf survey conditions. As of late 2025, user guides for software program SeeBird were no longer readily accessible, so this manual serves also to address that gap.

Northern Gulf of America↗

Geotechnical characterization of TriNet sites: A status report

The TriNet project, launched in 1997, created an improved, real-time seismic monitoring network in Southern California. Planning of the network began in 1995 ( e.g. , Heaton et al. , 1996 ), building on the success of the earlier TERRAscope network, which included 24 digital broadband and strong-motion instruments throughout Southern California ( e.g. , Kanamori et al. , 1993 ). At the end of the five-year TriNet project the network comprised 150 real-time digital broadband stations and another 400 strong-motion sensors, 50 of which were also real-time. This network is now recording digital broadband data for Southern California earthquakes at an unprecedented rate, data that are already proving valuable for investigations of earthquake sources and regional wave propagation, as well as earthquake response. In this report we describe an ongoing effort aimed at a full geotechnical characterization of the newly installed TriNet sets. Shallow geologic structure is known to play a substantial role in controlling the ground motions recorded at any site. Documentation of amplified ground motion at soft-sediment sites can be found among even early macroseismic observations of strong ground motions ( e.g. , Drake, 1815 ). Seismic waves are also now known to be strongly affected by deep basin structure as well ( e.g. , Frankel et al. , 1991 ; Field, 2000 ; Joyner, 2000 ). It will be necessary to understand these effects to exploit fully the rich data set being recorded at TriNet sites. Clearly, the nature of site conditions and site response at the recording sites must be understood for studies focused on ground motions and hazard from future large earthquakes, but it is also necessary to understand these effects to conduct earthquake source studies of both large and small earthquakes. Our multifaceted site characterization project involves geological/geotechnical site investigations, database development, and investigation of empirical amplification factors determined from broadband and strong-motion data recorded to date. Our goal is to complete a first-order geologic site characterization by 2005 and then employ appropriate methods to obtain direct constraint on shallow shear-wave velocity structure at each site. In the latter effort we will endeavor to find the most accurate and cost-effective methods to quantify geotechnical parameters at sites that have ranges of geologic site conditions and cultural settings. Our long-term goal is to obtain direct estimates of the average shear-wave velocity in the upper 30 m, Vs30, at each site. If velocity information is available to greater depths, for example at sites characterized under the ROSRINE project ( http://geoinfo.usc.edu/rosrine/ ), this information will be included in the database as well. The original motivation for using Vs30 to characterize near-surface velocity was pragmatic, determined by such factors as the typical reach of a drill rig in a single day. The parameter has become the accepted standard with which many seismic recording sites are characterized, however. In a recent study using mainshock and aftershock recordings at sites in the Los Angeles region, Wald and Mori ( 2000 ) observed good correlation between Vs30 and amplification at 1-7 Hz, albeit with significant scatter. The purpose of this report is threefold. First, we describe our ongoing efforts and present site characterization results collected to date for 62 broadband and strong-motion stations in and around the greater Los Angeles metropolitan region. Second, we present preliminary results that illustrate how improved geologic site characterizations can improve the correlation between site conditions and site response. Finally, we invite feedback from the community to guide our future investigations, in particular with respect to our ongoing database development efforts.

California↗

Fluorine, fluorite, and fluorspar in central Colorado

Fluorine (F) is a widespread element that was deposited in a variety of rocks, minerals, and geologic environments in central Colorado. It occurs as a trace element, as a major component of the mineral fluorite (CaFs), and as a major economic source of fluorine in fluorspar deposits, which are massive concentrations of fluorite. This study has compiled available geochemical analyses of rocks, both unmineralized and mineralized, to determine the distribution of fluorine in specific age-lithologic categories, ranging from 1.8-giga-annum (Ga) metamorphic rocks to modern soils, throughout central Colorado. It also draws upon field studies of fluorine-rich mineral deposits, including fluorspar deposits, to decipher the nearly two-billion-year-long geologic history of fluorine in the study area, with implications for mineral-resource evaluations and exploration. The resulting compilation provides an important inventory of the naturally occurring levels and sources of fluorine that ultimately weather, erode, and become part of surface waters that are used for domestic water supplies in densely populated areas along the Colorado Front Range. Most commonly, fluorine is a trace element in virtually all rocks in the region. In the 3,798 unmineralized rocks that were analyzed for fluorine in the study area, the average fluorine content was 1,550 parts per million (ppm). The median was 640 ppm, nearly identical to the average crustal abundance of 650 ppm, and some high-fluorine rocks in the Pikes Peak area skewed the average to a value much greater than the median. Most unmineralized age-lithologic rock suites, including Proterozoic metamorphic rocks, 1.7- and 1.4-Ga granitic batholiths, Cambrian igneous rocks, Phanerozoic sedimentary rocks, and Laramide and Tertiary igneous rocks, had median fluorine values of 400 to 740 ppm fluorine. In all suites, however, a small number of analyzed samples contained more than 1 percent (10,000 ppm) fluorine. The 1.1-Ga plutonic rocks related to the Pikes Peak batholith had a mean fluorine content of 1,700 ppm, and primary magmatic fluorite and fluorite-bearing pegmatites are common throughout that igneous mass. Fluorine was deposited in many types of economic mineral deposits in central Colorado, and it currently is a significant trace element in some thermal springs. In the fluorspar deposits, fluorine contents were as high as 37 percent. Some fluorine-rich porphyry systems, such as Jamestown, had fluorine values that ranged from 200 ppm to nearly 37 percent fluorine, and veins in other deposits contained hydrothermal fluorite, although it was not ubiquitous. For the 495 samples from non-fluorspar mining districts (and excluding Jamestown), however, the median fluorine content was 990 ppm. This is above the crustal average but still relatively modest compared to the fluorspar deposits, and it indicates that the majority of the mineralizing systems in central Colorado did not deposit large amounts of fluorine. Nevertheless, the fluorine- and fluorite-rich mineral deposits could be used as guides for the evaluation and discovery of related but concealed porphyry and epithermal base- and precious-metal deposits. The Cenozoic geologic history of central Colorado included multiple periods during which fluorine-bearing rocks and mineral deposits were exposed, weathered, and eroded. This protracted history has released fluorine into soils and regoliths, and modern rainfall and snowmelt interact with these substrates to add fluorine to the hydrosphere. This study did not evaluate the fluorine contents of water or make any predictions about what areas might be major sources for dissolved fluorine. However, the abundant data that are available on fluorine in surface water and ground water can be coupled with the results of this study to provide additional insight into natural sources of fluorine in domestic drinking water.

Scientific Investigations Report↗

The Glacier National Park: A popular guide to its geology and scenery

The Glacier National Park includes that part of the Front Range of the Rocky Mountains lying just south of the Canadian line, in Teton and Flathead counties, Mont. It is bounded on the west by Flathead River (locally called North Fork), on the south by the Middle Fork of Flathead River and the Great Northern Railway, and on the east by the Blackfeet Indian Reservation. Although this part of the Rocky Mountains has been known since Lewis and Clark crossed the continent in 1805-6, the region later made a park appears not to have been visited by white men until 1853, when Cut Bank Pass was crossed by A. W. Tinkham, one of the Government engineers engaged in exploring a route for the Pacific railroad. Tinkham, who was encamped in the Bitterroot Valley, was ordered to examine Marias Pass, but in traversing Middle Fork of Flathead River along the line of the present railroad he was evidently misled by the large size of the valley of Nyack Creek and ascended that instead of keeping to the right up the main stream. He reported the pass impracticable for railroad construction, and so this region dropped out of public attention for a long time. The next explorers to enter the region were a group of surveyors who, under the direction of American and British commissioners, established the international boundary line along the forty-ninth parallel from the Pacific coast to the main summit of the Rocky Mountains. This party reached the area now included in the park in the summer of 1861, and the stone monument shown in Plate I, B, which they erected on the Continental Divide west of Waterton Lake, still marks a point on the boundary between the United States and Canada. The land on the west side of the range formed a part of the public domain which, until the erection of the park, was open to settlement, but the land on the east originally belonged to the Blackfeet Indians and the white men had no rights upon it. About 1890 copper ore was found near the heads of Quartz and Mineral creeks, and a great boom for this region followed. Many prospectors drifted in, expecting to reap rich rewards from the discovery of mineral deposits and the general development of the region. Several of the main trails were built about this time, and considerable money was spent in prospecting, in opening mines, and in providing machinery to handle the large output of copper ore that was expected. The copper-bearing veins were found to extend through the range to the east side, but prospecting in that part of the mountains was not possible, for the land was included in the Indian reservation. This situation produced a growing discontent among the prospectors, who began to have a strong feeling that the Government should come to their relief by acquiring the coveted land and placing it at their disposal. The urgent demand of the prospectors and promoters was felt in Congress, and a bill was passed providing for the purchase from the Indians of this supposed mineral land for $1,500,000. In accordance with this act, a treaty with the Blackfeet Indians was signed at Browning, Mont., September 26, 1895, and approved by the Senate on June 10, 1896, by which the west line of the reservation was removed from the Continental Divide and was fixed along the eastern points of the spurs of the mountain range, as shown on the accompanying map, and the land so acquired was thrown open to mineral entry only. Under the stimulus of the new territory acquired, active prospecting was carried on for a time, but copper ore was found only in small quantities, and gradually the prospectors and miners drifted away to newer or more promising fields, and the region reverted to its original condition. For a long time it was visited only by hunters in search of big game and by summer visitors who, in order to escape the heat of the plains, were willing to undergo the privations and discomforts of the rude hotels then to be found in the region. Although these mountains had ceased to interest prospectors in search of mineral wealth, they still possessed a fascination for the lovers of natural scenery and almost everyone who visited the region was impressed with its alpine beauty. Again it was brought to the attention of Congress, which was urged to secure it for the use of the people by setting it aside as a national park. On May 11, 1910, President Taft signed the bill creating the Glacier National Park, and this most interesting and beautiful region thus became a permanent playground for the American people.

Bulletin↗

Implementation plan for the southern Pacific Border and Sierra-Cascade Mountains provinces

Introduction The National Cooperative Geologic Mapping Program (NCGMP) is publishing a strategic plan titled Renewing the National Cooperative Geologic Mapping Program as the Nation’s Authoritative Source for Modern Geologic Knowledge (Brock and others, in press). The plan provides a vision, mission, and goals for the program during the years 2020–2030, which are: Vision .—Create an integrated, three-dimensional, digital geologic map of the United States. Mission .—Characterize, interpret, and disseminate a national geologic framework model of the Earth through geologic mapping. Goal .—Focus on geologic mapping as a core function of the U.S. Geological Survey (USGS) within the long-term vision of adequately mapping the Nation’s geologic framework in three dimensions. In order to achieve the goals outlined in the strategic plan, the NCGMP has developed an implementation plan. This plan will guide the annual review of projects carried out by USGS staff (FEDMAP) described in the plan and the development of the annual FEDMAP prospectus that will ensure the effective application of the NCGMP strategy. This publication describes the implementation plan of the NCGMP strategy for the southern Pacific Border and Sierra-Cascade Mountains provinces, as defined by Fenneman (1917, 1928, and 1946). This implementation plan focuses on the geology of California and a sliver of Nevada surrounding Lake Tahoe. The southern Pacific Border and Sierra-Cascade Mountains provinces encompass the varied landscapes of the high Sierra Nevada, the Central Valley, and Coast Ranges in northern and central California and the Peninsular Ranges, Continental Borderland, Los Angeles Basin-San Gabriel-San Bernardino valleys, western and central Transverse Ranges, and northernmost Salton Trough in southern California. Societal demands create a need for earth-science data in each of these landscapes. The broader San Francisco Bay area, Central Valley, Los Angeles-San Gabriel-San Bernardino lowlands, and the coastal lowlands that border the Peninsular Ranges are densely populated (about 30 million people) areas at high risk of natural hazards. The mountains of the Sierra Nevada, Peninsular Ranges, and Transverse Ranges, and the coast all provide numerous recreational opportunities that attract visitors from around the world, whereas previously these ranges attracted people to mine their resources. The agricultural capacity of the Central Valley is a critical resource for the Nation that is increasingly water limited. The southern. Pacific Border and Sierra-Cascade Mountains provinces, at the edge of the North American continent, were profoundly influenced by subduction zone tectonics during the Mesozoic and early Cenozoic (ongoing in northernmost California) and subsequently by the inception, development, and present activity of the San Andreas transform margin system. Although the geology of this region is the poster child of fundamental conceptual models of subduction zone complexes, forearc basins, ophiolite obductions, magmatic arcs, and suspect terranes, as well as hosting one of Earth’s most notorious continental transform faults—the San Andreas Fault—important questions that have important societal consequences remain to be answered. Most of California’s population reside in these provinces and live within 30 miles of an active fault (according to www.earthquakeauthority.com ) yet new faults continue to be discovered, highlighting the importance of deformation off the main San Andreas Fault. Bedrock, surficial, and three-dimensional (3D) geologic maps depicting stratigraphic structure and depth to crystalline basement rocks provide critical context and information for understanding fault rupture, distributed deformation, fault connectivity, and history in addition to providing crucial data that enable forecasting of shaking amplitude and length from hypothetical earthquake scenarios. The tectonic evolution of California produced not only stunning mountains, with associated hazards from landslides and active volcanoes, but also fertile valleys that make California the top agricultural producer in the country in terms of cash receipts (according to www.ers.usda.gov/faqs ). These valleys lie atop large basins that not only store groundwater but, in many cases, host oil and gas fields, contributing to the fourth highest hydrocarbon production by State in the country in 2016 (according to https://www.aei.org/carpe-diem/animated-chart-of-us-oil-production-by-state-1981-2017 ). Water is a key resource increasingly stressed by growing agricultural, industrial, and residential needs. Warmer and drier conditions have led to an increased reliance on extracting groundwater resources, whose availability and quality are dictated at the first order by the 3D spatial distribution of bedrock and Quaternary surficial deposits. Thus, assessment of this critical resource is inextricably tied to knowledge of the surficial and subsurface geologic structure and material types.

California, Nevada↗

A Framework for Long-term Ecological Monitoring in Olympic National Park: Prototype for the Coniferous Forest Biome

This report is the result of a five-year collaboration between scientists of the U.S. Geological Survey Forest and Rangeland Ecosystem Science Center, Olympic Field Station, and the natural resources staff of Olympic National Park to develop a comprehensive strategy for monitoring natural resources of Olympic National Park. Olympic National Park is the National Park Serviceʼs prototype monitoring park, representing parks in the coniferous forest biome. Under the umbrella of the National Park Serviceʼs prototype parks program, U.S. Geological Survey and Olympic National Park staffs are obligated to: develop strategies and designs for monitoring the long-term health and integrity of national park ecosystems with a significant coniferous forest component. design exportable monitoring protocols that can be used by other parks within the coniferous forest biome (i.e., parks having similar environments), and create a demonstration area and ʻcenter of excellenceʼ for assisting other parks in developing ecological monitoring programs. Olympic National Park is part of the North Coast and Cascades Network, a network of seven Pacific Northwestern park units created recently by the National Park Serviceʼs Inventory and Monitoring Program to extend the monitoring of ʻvital signsʼ of park health to all National Park Service units. It is our intent and hope that the monitoring strategies and conceptual models described here will meet the overall purpose of the prototype parks monitoring program in proving useful not only to Olympic National Park, but also to parks within the North Coast and Cascades Network and elsewhere. Part I contains the conceptual design and sampling framework for the prototype long-term monitoring program in Olympic National Park. In this section, we explore key elements of monitoring design that help to ensure the spatial, ecological, and temporal integration of monitoring program elements and discuss approaches used to design an ecosystem-based monitoring program. Basic monitoring components include ecosystem drivers, (e.g., climate, atmospheric inputs, human pressures), indicators of ecosystem integrity (e.g., biogeochemical indicators), known threats (e.g., impacts of introduced mountain goats), and focal or ʻkeyʼ species (e.g., rare or listed species, Roosevelt elk). Monitoring system drivers and key indicators of ecosystem integrity provide the long-term baseline needed to judge what constitutes ʻunnaturalʼ variation in park resources and provide the earliest possible warning of unacceptable change. Monitoring effects of known threats and the status of focal species will provide information useful to park managers for dealing with current park issues. In Part I we describe the process of identifying potential indicators of ecological condition and present conceptual models of park ecosystems. In addition we report results from several workshops held in conjunction with Olympic National Park aimed at identifying potential indicators of change in the parkʼs ecosystem. First, we describe the responses of Olympic National Park staff to the generic question, “What is the most important resource to monitor in Olympic National Park and why?” followed by the responses from resource and land managers from areas adjoining the park. We also catalogue the responses of various expert groups that we asked to help identify the most appropriate system drivers and indicators of change in the Olympic National Park ecosystems. Results of the workshops provided the justification for selecting basic indicators of ecosystem integrity, effects of current threats to park resources, and focal resources of parks to detect both the currently evident and unforeseeable changes in park resources. We conclude Part I by exploring several generic statistical issues relevant to monitoring natural resources in Olympic National Park. Specifically we discuss trade-offs associated with sampling extensively versus sampling intensively in smaller geographic regions and describe a conceptual framework to guide development of a generic sampling frame for monitoring. We recommend partitioning Olympic National Park into three zones of decreasing accessibility to maximize monitoring efficiency. We present examples of how the generic sampling frame could be used to help ensure spatial integration of individual monitoring projects. Part II of the report is a record of the potential monitoring questions and indicators identified to date in our workshops. The presentation is organized according to the major system drivers, components, and processes identified in the intermediate-level working model of the Olympic National Park ecosystem. For each component of the park system, we develop the need and justification for monitoring, articulate park management issues, and describe key resources and ecosystem functions. We also present a pictorial conceptual model of each ecological subsystem, identify monitoring questions, and list potential indicators for each monitoring question. We conclude each section by identifying linkages of indicators to other ecological subsystems in our general ecosystem model, spatial and temporal contexts for monitoring (where and how often to monitor), and research and development needs. Part II represents the most current detailed listing of potential indicators—the material for subsequent discussions of monitoring priorities and selection of indicators for protocol development. Collectively, the sections of this report contain a comprehensive list of the important monitoring questions and potential indicators as well as recommendations for designing an integrated monitoring program. In Part I, Chapter 6 we provide recommendations on how to proceed with the important next steps in the design process: establishing priorities among the many possible monitoring questions and indicators, and beginning to research and design effective long-term monitoring protocols.

Information and Technology Report↗

Preliminary map showing the thickness of glacial deposits in Ohio

This map was compiled as part of a U.S. Geological Survey project to portray the character and thickness of glacial deposits east of the Rocky Mountains. It is a preliminary map, compiled at a scale of 1:250,000 from sources of differing type and quality (fig. I, table 1). County drift-thickness maps produced by the Ohio Geological Survey are the most reliable source of data used. These county maps, either published or in manuscript form, were available for 36 of the 72 counties that contain glacial or glacially related deposits more than 50 ft thick. An earlier, reconnaissance drift-thickness map of northwestern Ohio (Ohio Department of Natural Resources, 1960) provided data for parts of 12 additional counties. Published maps of bedrock topography were utilized for 9 counties and for parts of 6 counties. The computation of drift thickness in these counties is less reliable then the drift-thickness maps published by the Ohio Geological Survey because the bedrock topography maps, most of which are at least 20 years old, were compiled from a more limited set of data. In the remaining counties, only well-log data were available. Drift-thickness data in these counties are the least reliable; the determination of thickness trends was commonly guided by the configuration of preglacial drainage channels (as drawn by Stout and others, 1943) and by topography. Throughout the glaciated areas of Ohio, and especially near the glacial margin, a network of drainage channels is buried by glacial drift. Some of these channels arc of a preglacial age while others were formed during interglacial periods; some were cut into the bedrock by rivers— either preglacial rivers or meltwater streams flowing along the glacial margin—while others, such as the lower Cuyahoga River valley, may have been overdeepened by the erosive power of glacial ice. The thickest deposits of drift in the state occur in the lower Cuyahoga River valley near Lake Erie, and above the buried valley in west-central Ohio that has been historically considered a part of the Teays Valley. The thickest drift reliably measured in Ohio to date is 602 ft in a drill hole near Cleveland. In contrast to the extreme variations in drift thickness encountered in the vicinity of buried channels, drift on the upland arcus is generally thinner and the variations in thickness are much less pronounced. Worthy of note, however, are three large areas where the drift sheet is relatively thick. In northwestern Ohio, a large volume of drift was deposited along the flanks of the Erie ice lobe (fig. 2) near the interlobate position with the Saginaw lobe to the northwest; drift thickness there exceeds 200 ft. Thick drift was also deposited in a roughly cast-west band across the Miami lobe. The mechanism that produced this band of thick drift is not obvious, but it may have been influenced in part by bedrock topography. Bedrock control of drift thickness is more clearly indicated to the cast of Columbus, along the eastern flank of the Scioto lobe, where ice slow was resisted by rocks of the Allegheny plateau. The edge of the plateau, or the Allegheny escarpment, is obscured by glacial deposits but its likely position (Fenneman, 1938; Stout and others, 1913; Dove, 1960; and Root and others, 1961) is shown on the map. Southward from the ice margin's reentrant position in southern Richland County, ice flowing eastward from the Scioto lobe encountered the topographically higher plateau, which constrained the ice and caused drift to accumulate in significant thicknesses just to the west of the escarpment.

Ohio↗

Coordinated bird monitoring: Technical recommendations for military lands

The Department of Defense (DoD) is subject to several rules and regulations establishing responsibilities for monitoring migratory birds. The Sikes Act requires all military installations with significant natural resources to prepare and implement Integrated Natural Resources Management Plans (INRMPs). These plans guide the conservation and long-term management of natural resources on military lands in a manner that is compatible with and sustains the military mission. An INRMP also supports compliance with all legal requirements and guides the military in fulfilling its obligation to be a good steward of public land.The management and conservation of migratory birds is addressed in installation INRMPs. The National Environmental Policy Act (NEPA) requires federal agencies to evaluate and disclose the potential environmental impacts of their proposed actions. More recently, DoD signed an MOU (http://www.dodpif.org/downloads/EO13186_MOU-DoD.pdf) for migratory birds, under Executive Order 13186, with the US Fish and Wildlife Service (USFWS) in July 2006 and a Migratory Bird Rule (http://www.dodpif.org/downloads/MigBirdFINALRule_FRFeb2007.pdf) was passed by Congress in February 2007. The Migratory Bird Rule addresses the potential impacts of military readiness activities on populations of migratory birds and establishes a process to implement conservation measures if and when a military readiness activity is expected to have a significant adverse impact on a population of migratory bird species (as determined through the NEPA process). The MOU states that for nonmilitary readiness activities, prior to initiating any activity likely to affect populations of migratory birds DoD shall (1) identify the migratory bird species likely to occur in the area of the proposed action and determine if any species of concern could be affected by the activity, and (2) assess and document, using NEPA when applicable, the effect of the proposed action on species of concern. By following these procedures, DoD will minimize the possibility for a proposed action to unintentionally take migratory birds at a level that would violate any of the migratory bird treaties and potentially impact mission activities. In addition, implementing conservation and monitoring programs for migratory birds supports the ecosystem integrity necessary to sustain DoD's natural resources for the military mission.Non-compliance with the procedural requirements of the MBTA could result in a private party lawsuit under the Administrative Procedures Act (APA). A lawsuit filed under APA involving a Navy bombing range is the basis for a court ruling that unintentional take of migratory birds applies to federal actions. Ensuring the necessary data is available to adequately assess impacts of a proposed action will help avoid lawsuits or help ensure such lawsuits have no grounds. The data gathered in a bird monitoring program will provide the best scientific data available to assess the expected impacts of a proposed action on migratory bird species through the NEPA process. This report presents recommendations developed by the U.S. Geological Survey (USGS) for the Department of Defense (DoD) on establishing a "Coordinated Bird Monitoring (CBM) Plan." The CBM Plan is intended to ensure that DoD meets its conservation and regulatory responsibilities for monitoring birds (Chapter 1). The report relies heavily on recommendations in the report, "Opportunities for improving avian monitoring" (http://www.nabci-us.org/aboutnabci/monitoringreportfinal0307.pdf), by the U.S. North American Bird Conservation Initiative (U.S. NABCI Monitoring Subcommittee, 2007) and on a review of 358 current DoD bird monitoring programs carried out as part of this project (Chapter 2). This report contains 12 recommendations which, if followed, would result in a comprehensive, efficient, and useful approach to bird monitoring. The recommendations are based on the entire report but are presented together at the end of Chapter 1. DoD has agreed to consider implementing these recommendations; however, final decisions will be based upon such factors as the availability of resources and military mission considerations. These recommendations from USGS can be summarized into 6 major themes: A major report on monitoring was released in 2007 by the U.S. North American Bird Conservation Initiative (http://www.nabci-us.org/main2.html). DoD can be consistent with this report by establishing policy that monitoring will be explicitly acknowledged as an integral element of bird management and conservation (Recommendation 1). The design of monitoring and assessment programs for birds should include the following steps: Preparation of a document describing the program's goals, objectives, and methods similar to a format we provide (Recommendation 2, Chapter 4). Selection of field methods using an "expert system" developed in this project (Recommendation 3, Chapter 5) or another well-documented system. Preparation and storage of metadata describing the monitoring program in the Natural Resources Monitoring Partnership (NRMP), and other appropriate databases Recommendation 4, Chapter 6). Entry of the survey data using eBird (http://ebird.org/content/dod) or the Coordinated Bird Monitoring Database (CBMD) and long-term storage of the data in the CBMD and the Avian Knowledge Network (AKN; Recommendation 5, Chapter 6; http://www.avianknowledge.net/). Submission of major results from the monitoring program for publication in a peer reviewed journal (Recommendation 6). The DoD Legacy Resource Management Program (Legacy; https://www.dodlegacy.org), Environmental Security Technology Certification Program (ESTCP; http://www.serdp.org/), and Strategic Environmental Research and Development Program (SERDP; http://www.serdp.org/) should be encouraged to continue their significant contributions to the foundations of bird monitoring (Recommendation 7, Chapters 1 and 3). Appropriate monitoring should be conducted to identify species of concern on installations. A year-round, one-time survey of birds on installations with habitat for migratory birds would provide the most information to assist compliance with the MOU, the Final Rule, and the NEPA analyses of proposed actions. However, less intensive survey efforts can still be conducted to yield useful information. We describe how various levels of survey effort might be organized and conducted. In addition, continuing surveys, as feasible, would further assist in documenting effects of military readiness and non-readiness activities on species of concern (SOC) (Recommendation 8, Chapter 7). Participation in well-designed, large-scale surveys [(e.g., North American Breeding Bird Survey (BBS; http://www.pwrc.usgs.gov/bbs/), Monitoring Avian Productivity and Survivorship (MAPS; http://www.birdpop.org/maps.htm)] on land that DoD manages or on lands where the results will be of high interest to DoD, will provide DoD and other NABCI members with information important to bird conservation (Recommendation 9, Chapter 8). Review and implementation of the CBM Plan should involve both higher level management and installation-level natural resources managers (Recommendation 11), be implemented through cooperative partnerships (Recommendation 12), and be followed on U.S territory lands and Army Corps of Engineers projects (Recommendation 10).Additional recommendations that pertain to implementing the DoD CBM Plan are discussed in Chapter 9.

Open-File Report↗

Uranium-bearing coal in the central part of the Great Divide basin, Sweetwater County, Wyoming

Field work leading to this report was done by the U.S. Geological Survey for the Division of Raw Materials of the U.S. Atomic Energy Commission. Nearly 24 townships were mapped in the central part of the Great Divide Basin, Sweetwater County, Wyoming. Fourteen of these townships contain outcrops of uranium-bearing coal. Thirty coal beds were mapped, but only seven of them have uranium-bearing coal reserves as defined in this report. Coal beds 2.5 or more feet thick are considered in calculating coal reserves, and of these, only beds containing 0.003 or more percent uranium are considered in calculating reserves of uranium in coal. Reserves of uranium in coal ash include those beds 2.5 or more feet thick that contain 0.015 or more percent uranium in coal ash. Measured and indicated coal reserves total about 700,000,000 short tons which contain about 2,600 short tons of uranium in the coal, or about 2,400 short tons of uranium in the coal ash. Strippable reserves, defined as reserves in beds beneath 60 or less feet of overburden, are about 250,000,000 short tons of coal containing about 1,100 short tons of uranium in coal, or about 600 tons of uranium in coal ash. The thickest coal beds underlie a relatively narrow belt that trends northwest and coincides approximately with the axis of the Red Desert syncline. The coal beds contain the most uranium on the east flank of the syncline near the southwesternmost edge of the Battle Spring formation (new). This formation is of early and middle Eocene age and consists predominantly of very coarse-grained arkosic sandstone which is highly permeable. It intertongues southwestward with the Tess permeable Green River and Wasatch formations. The Green River formation consists from youngest to oldest of the Morrow Creek and Laney shale members and the Tipton and Luman (new) tongues. The Wasatch formation interfingers with the Green River formation and consists from youngest to oldest of the Cathedral Bluffs, Niland, and Red Desert tongues. The latter two are here recognized for the first time and contain all the coal beds in the Wasatch formation. The Morrow Creek member of the Green River formation is of middle Eocene age. The Laney member and Cathedral Bluffs tongue of the Green River and Wasatch formation, respectively intertongue and are in part equivalent, and are of early or middle Eocene age. The other units are of early Eocene age. A broad gentle arch which is either the eastern extension of the Wamsutter arch, or a separate arch en echelon to it, separates the Washakie Basin the southeastern part of the map area from the Red Desert syncline in the northeastern part. The Red Desert syncline plunges gently northwest into the nearly circular structural Niland basin. The south flank of the Wamsutter(?) arch dips southeastward at an average rate of about 230 feet per mile. The northern flank of the arch dips northeastward at an average rate of about 140 feet per mile. North of the Niland basin, the structure is more complex and is dominated by a northwestwardtrending graben which parallels and includes the Cyclone Rim in the northern part of die map area. Schroeckingerite deposits in the northeastern part of the map area lie within this graben in a sequence of arkosic sandstone and clay shale stratigraphically equivalent to the lower part of the Cathedral Bluffs tongue of the Wasatch formation and to the uppermost part of the Tipton tongue of the Green River formation. Weakly radioactive tuffaceous sandstone beds of the Browns Park formation that probably once blanketed the entire Great Divide Basin are preserved a short distance north and northwest of the schroeckingerite deposits. The geologic settings of the schroeckingerite and uranium-bearing coal are similar and their source of uranium is probably the same. The uranium probably was leached from tuffaceous beds in the Browns Park formation and carried to its present site of deposition by groundwater whose circulation was guided by structure and changes in facies.

Wyoming↗

Natural offshore oil seepage and related tarball accumulation on the California coastline — Santa Barbara Channel and the southern Santa Maria Basin; source identification and inventory

Oil spillage from natural sources is very common in the waters of southern California. Active oil extraction and shipping is occurring concurrently within the region and it is of great interest to resource managers to be able to distinguish between natural seepage and anthropogenic oil spillage. The major goal of this study was to establish the geologic setting, sources, and ultimate dispersal of natural oil seeps in the offshore southern Santa Maria Basin and Santa Barbara Basins. Our surveys focused on likely areas of hydrocarbon seepage that are known to occur between Point Arguello and Ventura, California. Our approach was to 1) document the locations and geochemically fingerprint natural seep oils or tar; 2) geochemically fingerprint coastal tar residues and potential tar sources in this region, both onshore and offshore; 3) establish chemical correlations between offshore active seeps and coastal residues thus linking seep sources to oil residues; 4) measure the rate of natural seepage of individual seeps and attempt to assess regional natural oil and gas seepage rates; and 5) interpret the petroleum system history for the natural seeps. To document the location of sub-sea oil seeps, we first looked into previous studies within and near our survey area. We measured the concentration of methane gas in the water column in areas of reported seepage and found numerous gas plumes and measured high concentrations of methane in the water column. The result of this work showed that the seeps were widely distributed between Point Conception east to the vicinity of Coal Oil Point, and that they by in large occur within the 3-mile limit of California State waters. Subsequent cruises used sidescan and high resolution seismic to map the seafloor, from just south of Point Arguello, east to near Gaviota, California. The results of the methane survey guided the exploration of the area west of Point Conception east to Gaviota using a combination of seismic instruments. The seafloor was mapped by sidescan sonar, and numerous lines of high -resolution seismic surveys were conducted over areas of interest. Biomarker and stable carbon isotope ratios were used to infer the age, lithology, organic matter input, and depositional environment of the source rocks for 388 samples of produced crude oil, seep oil, and tarballs mainly from coastal California. These samples were used to construct a chemometric fingerprint (multivariate statistics) decision tree to classify 288 additional samples, including tarballs of unknown origin collected from Monterey and San Mateo County beaches after a storm in early 2007. A subset of 9 of 23 active offshore platform oils and one inactive platform oil representing a few oil reservoirs from the western Santa Barbara Channel were used in this analysis, and thus this model is not comprehensive and the findings are not conclusive. The platform oils included in this study are from west to east: Irene, Hildago, Harvest, Hermosa, Heritage, Harmony, Hondo, Holly, Platform A, and Hilda (now removed). The results identify three “tribes” of 13 C-rich oil samples inferred to originate from thermally mature equivalents of the clayey-siliceous, carbonaceous marl, and lower calcareous-siliceous members of the Monterey Formation. Tribe 1 contains four oil families having geochemical traits of clay-rich marine shale source rock deposited under suboxic conditions with substantial higher-plant input. Tribe 2 contains four oil families with intermediate traits, except for abundant 28,30-bisnorhopane, indicating suboxic to anoxic marine marl source rock with hemipelagic input. Tribe 3 contains five oil families with traits of distal marine carbonate source rock deposited under anoxic conditions with pelagic but little or no higher-plant input. Tribes 1 and 2 occur mainly south of Point Conception in paleogeographic settings where deep burial of the Monterey Formation source rock favored generation from all three members or their equivalents. In this area, oil from the clayey-siliceous and carbonaceous marl members (Tribes 1 and 2) may overwhelm that from the lower calcareous-siliceous member (Tribe 3) because the latter is thinner and less oil-prone than the overlying members. Tribe 3 occurs mainly north of Point Conception, where shallow burial caused preferential generation from the underlying lower calcareous-siliceous member or another unit with similar characteristics. It is very desirable to be able to clearly distinguish the naturally occurring seep oils from the anthropogenically derived platform oils. Within the “training set” of oils and tars (388 samples), the biomarker parameters are sometimes sufficient to allow unique discrimination of individual platform oils. More often however, platform samples and seep samples with sources geographically close to each other are too similar to each other, with respect to the biomarker parameters, to definitively differentiate them on that basis alone. In some cases other parameters can be helpful. These other parameters are related to the degree of biogeochemical degradation or weathering that the oils or tars have experienced. These components include the typical oil distribution of n-alkane hydrocarbons and isoprenoids pristane and phytane. All of the platform oils in our sample set contain these components. On the other hand, the seep oils or tars have been exposed to significant biodegradation while in the near subsurface. The majority, but not all of seep oils or tars have been biodegraded up to or beyond the loss of n-alkanes and isoprenoids. Seep oils found in the vicinity of Coal Oil Point or Arroyo Burro are apparently the least weathered and are particularly likely to retain significant n-alkanes and isoprenoids. Therefore the combination of chemometric fingerprinting and the presence or absence of n-alkanes and isoprenoids help to differentiate anthropogenic production oils versus natural seeps oils and tars. The differentiation is not always definitive because of the close chemical similarity of some samples and the variability in the biodegradation progression. This is the case near Coal Oil Point, and near Platform A (Dos Cuadros Field) where seep oils and Platform Holly and Platform A oils are genetically very similar and cannot be definitively distinguished after a period of a few days of weathering. In contrast, oils from the Point Conception platforms can be distinguished on the basis of chemometric fingerprinting alone. In the middle of this spectrum are oils from Platforms Harmony, Heritage, and Hondo, where it is expected that oil weathering would take on the order of two weeks to a month to produce tarballs similar to those seen near Point Conception. In this case there is a much greater degree of weathering needed to proceed from produced oil to the biodegraded tar characteristic of tarball stranded on the beach. Tar deposition on beaches was monitored as part of cooperative with the County of Santa Barbara Energy Division and the U.S. Geological Survey during 2001-2003. We found tar deposition varies on a seasonal basis. In general, tarballs accumulate at a faster rate or remain longer on all beaches during the summer and fall months. The reasons for this are unclear based on our limited observations, however we speculate that factors such as prevailing winds and currents combined with more quiescent wave conditions favors the accumulation and preservation of tarballs on the beach during the summer and fall months. In contrast, winter storms, with much greater wave action remove beach sand and other materials, and stormy seas tend to break up oil that might weather into tarballs. Natural seepage is affected by the spring/neap tidal cycle; however, the link to tar deposition is unclear. Longer periods of monitoring are needed to address the variability in the data and provide a more robust statistical analysis.

California↗

Simulation of groundwater flow and brine discharge to the Dolores River in the Paradox Valley, Montrose County, Colorado

Salinity, or total dissolved solids (TDS), of the Colorado River affects agricultural, municipal, and industrial water users and is an important concern in the Western United States. In the Paradox Valley of southwestern Colorado, natural discharge of sodium-chloride brine to the Dolores River from the underlying core of a salt-valley anticline accounts for about 6 percent of the salinity load to the Colorado River. Formation of the Paradox Valley began during the Miocene, and subsequent erosion exposed the Pennsylvania Paradox Formation in the core of the anticline where a cap rock, collapse features, breccia, and sodium-chloride saturated brine developed at the top of the exposed salt diapir. The discharge of brine to the Dolores River is affected by these dissolution features, along with seasonal hydrologic conditions and density-dependent flow between older dense brine and the younger fresh groundwater in the overlying alluvial aquifer. To reduce TDS concentrations in the Dolores River through the Paradox Valley, the Bureau of Reclamation has pumped brine from a series of shallow wells adjacent to the river since July 1996. The pumped brine is collected and piped to a deep disposal well where it is injected into the Mississippian Leadville Limestone at a depth of about 4,570-meters below land surface. The pumping and injection operation is collectively known as the Paradox Valley Unit (PVU), and by 2015, the PVU had substantially reduced TDS concentrations in the Dolores River by about 70 percent. Since 2019, injection-pressure limits and related seismic activity have constrained deep-well injection and thus brine pumping at the PVU. In cooperation with the Bureau of Reclamation, the U.S. Geological Survey developed a MODFLOW-6 three-dimensional, variable-density groundwater flow and TDS transport model of the Paradox Valley to evaluate the effects of PVU pumping operations on brine discharge to the Dolores River and to guide additional research. The finite-difference model grid consists of 76 rows and 48 columns oriented from northwest to southeast in alignment with valley topography and groundwater-flow directions in the near-surface freshwater alluvial aquifer. A 7-layer hydrogeologic framework was developed from existing datasets to represent the alluvial aquifer, cap rock, collapse breccia, and groundwater flow and TDS transport from the underlying Paradox Formation salt to the Dolores River. The model represents a 33-year transient calibration period from 1987 through 2020 that includes pre-PVU conditions from 1987 through June 1996 and post-PVU conditions from July 1996 through 2020. A 1,000-year simulation of groundwater flow and coupled TDS transport computed the initial conditions for the subsequent 33-year transient simulation. Observations of precipitation, streamflow, evaporation, agricultural land use, and PVU brine pumping rates were used to specify appropriate boundary conditions to the model representing time-varying recharge, tributary streamflow, groundwater underflow, evapotranspiration (ET), and PVU pumping. Values for average monthly streamflow and TDS concentration at the upstream streamgage, the Dolores River at Bedrock (USGS streamgage 09169500), were specified as model input where the Dolores River enters Paradox Valley. Observed pumping from the PVU, water levels and TDS concentrations in groundwater, and streamflow and estimated TDS concentrations at the downstream streamgage, the Dolores River near Bedrock (USGS streamgage 09171100), were calibration targets that constrained the manual calibration of model parameters representing aquifer hydraulic conductivity, storage, streambed conductance, recharge, and (ET). Two primary model-calibration targets were the match between observed and simulated TDS mass flux from PVU pumping wells and the match between estimated and simulated TDS mass flux to the Dolores River. The simulated TDS mass withdrawn by pumping wells is calculated by the model as the product of the assigned pumping rate and simulated groundwater TDS concentrations. Because actual pumping rates were assigned as simulated values, the total simulated PVU pumping for the 33-year calibration is within 0.5 percent of the observed values. However, simulated concentrations and thus mass flux of TDS withdrawn by the PVU pumping wells were consistently about 26 percent less than observed values for all the simulated time periods (33-year simulation, pre-PVU, and post-PVU). The representation of brine inflow was explored through additional modeling to evaluate the effect of the simulated brine source on groundwater TDS concentrations. Results indicated that a saturated-salt constant-flux brine source best replicated the magnitude and transient pattern observed for TDS mass flux from PVU pumping wells. The simulated TDS mass flux to the Dolores River is compared to estimates based on observed streamflow and specific conductance (SC) data for the downstream streamgage. The calibrated model provided a close fit of simulated to measured streamflow at the downstream streamgage, and the calibrated model fit to estimated TDS concentrations at the downstream streamgage was reasonable. The greatest differences between simulated and estimated values occurred during drought periods from June 2000 to March 2003, May 2012 to June 2013, and October 2013 to October 2014, when simulated TDS concentrations in the river were greater than estimated concentrations. In general, simulated TDS mass flux to the river for the pre-PVU period is in good agreement with estimated values (2-percent difference), but the model overestimated TDS mass flux to the river by about 41 percent during the post-PVU period. The model uncertainty with respect to TDS mass flux to the river indicates other processes or model parameters not well represented by the model are affecting the system, especially during drought. During model calibration, the most sensitive parameters were identified as vertical hydraulic conductivity of the alluvial aquifer, conductance of the Dolores River streambed, ET extinction depth and rate, and recharge rate. Five 5-year scenarios of conditions for 2021–25 were simulated to assist evaluation of alternative strategies to manage the discharge of brine into the Dolores River. The first scenario simulates no PVU pumping and serves as a base case for comparison to the other scenarios. Two scenarios simulate the effects of varying withdrawal timing at an annual rate about one-third less than during 2010 through 2018. During high-flow spring snowmelt runoff periods when brine discharge is naturally minimized, PVU pumping does not substantially affect salinity in the Dolores River, and comparison of these two scenarios indicates that scheduling brine withdrawals during times of low river stage is nearly as effective at reducing TDS mass flux to the river as pumping brine year-round. Cessation of pumping during periods of high river stage may be advantageous for system maintenance, brine injection, and seismic-risk reduction. The fourth scenario tested the effect of reducing irrigation-return flow on brine discharge and predicted a slight reduction of TDS mass flux to the Dolores River, but not as great a reduction as that of using the PVU to remove brine. The fifth scenario simulated 5 years of drought conditions without PVU pumping and indicates brine discharge during drought about 15 percent greater than during average hydrologic conditions. Results from scenario 5 are consistent with the calibrated model results and indicate that aquifer properties and ET processes and parameters may be affecting simulation results during drought. The Paradox Valley groundwater model provides a reasonable overall match to observed conditions in the Dolores River. The model is useful for evaluating relative differences between brine management scenarios to inform PVU operational decisions and to identify gaps in data and process understanding. Representation of the brine source, hydraulic-conductivity parameters, and recharge and ET processes were identified as potential areas for additional field and modeling research. Additional research in the Paradox Valley might include field-data collection that provides additional information on the hydrogeologic framework, groundwater levels, groundwater TDS concentrations, stream characteristics, and aquifer properties. Additional modeling efforts could benefit from applying advanced tools for model development, calibration, and visualization including parameter-estimation and sensitivity analysis. Statistical evaluation of known model uncertainties such as hydraulic conductivity, streambed conductance, representations of the brine source, recharge, and ET could improve the match between simulated and estimated TDS mass flux from PVU pumping wells and to the Dolores River further informing model predictions and system understanding for the Paradox Valley.

Colorado↗

Hydrogeology and water quality of the Leetown area, West Virginia

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

West Virginia↗