USGS Science⌕ Search

SEARCH · USGS Science

Results for “Scientific Reports”

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,027 records · Page 57Linked to original sources
USGS2009

Occurrence and Trends of Selected Chemical Constituents in Bottom Sediment, Grand Lake O' the Cherokees, Northeast Oklahoma, 1940-2008

After over 100 years of continuous activity, lead and zinc mining in the Tri-State Mining District (hereafter referred to as the TSMD) in parts of southeast Kansas, southwest Missouri, and northeast Oklahoma ended in the 1970s. The mining activity resulted in substantial historical and ongoing input of cadmium, lead, and zinc to the environment including Grand Lake O' the Cherokees (hereafter referred to as Grand Lake), a large reservoir in northeast Oklahoma. To help determine the extent and magnitude of contamination in Grand Lake, a one-year study was conducted by the U.S. Geological Survey in cooperation with the U.S. Fish and Wildlife Service. Bottom-sediment coring at five sites was used to investigate the occurrence of cadmium, lead, zinc, and other selected constituents in the bottom sediment of Grand Lake. Cadmium concentrations in the bottom sediment of Grand Lake ranged from 2.3 to 3.6 mg/kg (milligrams per kilogram) with a median of 3.5 mg/kg (5 samples). Compared to an estimated local background concentration of 0.6 mg/kg, the historical mining activity increased cadmium concentrations by about 280 to 500 percent. Lead concentrations ranged from 35 to 102 mg/kg with a median of 59 mg/kg (50 samples). Compared to an estimated local background concentration of 20 mg/kg, the historical mining activity increased lead concentrations by about 75 to 410 percent. The range in zinc concentrations was 380 to 986 mg/kg with a median of 765 mg/kg (50 samples). Compared to an estimated local background concentration of 100 mg/kg, the historical mining activity increased zinc concentrations by about 280 to 890 percent. With the exception of the most upstream coring site, the lead and zinc depositional profiles generally were similar in terms of the range in concentrations measured and the temporal pattern observed. Depositional profiles for lead and zinc indicated mid-core peaks followed by concentrations that decreased since about the 1980s. The depositional profiles reflect the complex interaction of several factors including historical mining and related activities, mine drainage, remediation, landscape stabilization, precipitation and associated runoff, and the erosion and transport of contaminated and clean sediments within the basin. Compared to sediment-quality guidelines, the Grand Lake samples had cadmium concentrations that were substantially less than the general probable-effects concentration (PEC) (4.98 mg/kg) and a TSMD-specific PEC (11.1 mg/kg). The PECs represent the concentration above which toxic biological effects are likely to occur. Likewise, all sediment samples had lead concentrations that were substantially less than the general PEC (128 mg/kg) and a TSMD-specific PEC (150 mg/kg). Zinc concentrations typically exceeded the general PEC (459 mg/kg), but were substantially less than a TSMD-specific PEC (2,083 mg/kg). Throughout the history of Grand Lake, lead and zinc concentrations in the deposited sediment did not approach or exceed the TSMD-specific PECs. As of 2008, legacy effects of mining still included the delivery of contaminated sediment to Grand Lake by the Spring and Neosho Rivers. The Neosho River, with its larger flows and less-contaminated sediment, likely dilutes the load of contaminated sediment delivered to Grand Lake by the Spring River. The information contained in this report provides a baseline of Grand Lake conditions with which to compare future conditions that may represent a response to changes in mining-related activity in the Grand Lake Basin.

Scientific Investigations Report↗
USGS2017

Three approaches for estimating recovery factors in carbon dioxide enhanced oil recovery

Preface The Energy Independence and Security Act of 2007 authorized the U.S. Geological Survey (USGS) to conduct a national assessment of geologic storage resources for carbon dioxide (CO 2 ) and requested the USGS to estimate the “potential volumes of oil and gas recoverable by injection and sequestration of industrial carbon dioxide in potential sequestration formations” (42 U.S.C. 17271(b)(4)). Geologic CO 2 sequestration associated with enhanced oil recovery (EOR) using CO 2 in existing hydrocarbon reservoirs has the potential to increase the U.S. hydrocarbon recoverable resource. The objective of this report is to provide detailed information on three approaches that can be used to calculate the incremental recovery factors for CO 2 -EOR. Therefore, the contents of this report could form an integral part of an assessment methodology that can be used to assess the sedimentary basins of the United States for the hydrocarbon recovery potential using CO 2 -EOR methods in conventional oil reservoirs.

Scientific Investigations Report↗
USGS2018

Compilation and assessment of resource values and hazards to inform transportation planning and associated land-use planning

Land-use planning has an important role in local, regional, State, and Federal land management, and planning efforts can benefit from consistent, spatially explicit information that can help guide priorities and decisions. The credibility and relevance of information used to inform planning activities depends on the availability of consistent information about the resources and values of interest or concern within the planning area. To support long-range transportation planning and other regional land-use planning efforts, the U.S. Geological Survey gathered, processed, interpreted, and compiled spatial datasets representing a wide range of information on terrestrial and aquatic ecosystem condition and importance, cultural (historical) features and places, and natural hazards. This report describes the spatial data compiled to represent natural landscape conditions (including social, cultural, and natural attributes) to estimate the potential importance of lands for wildlife, wild habitats, recreation, and conservation based on abundance of species, habitats, land and water conditions, and conservation designations. Abundance of resources, including the potential number of species, presence of important habitats and protected areas, and proximity to particular features or habitats, indicates the potential sensitivity of the natural landscape to land use, especially transportation networks. The source data, derived indices, and the methods for processing these data are described in this final report. The datasets referenced in the report are available from the U.S. Geological Survey ( https://www.sciencebase.gov/catalog/ and https://doi.org/10.5066/F7MW2F8W ) or the Central Federal Lands Highway Division of the Office of Federal Lands Highway ( https://flh.fhwa.dot.gov ).

Scientific Investigations Report↗
USGS2018

Flood-inundation maps for Cayuga Inlet, Sixmile Creek, Cascadilla Creek, and Fall Creek at Ithaca, New York

Digital flood-inundation maps for a 2.9-square-mile area of Ithaca, New York, were created in 2015–18 by the U.S. Geological Survey in cooperation with the City of Ithaca, New York, and the New York State Department of State. The flood-inundation maps depict estimates of the maximum areal extent and depth of flooding corresponding to selected flood frequencies for Cayuga Inlet, Sixmile Creek, Cascadilla Creek, and Fall Creek and selected water-surface elevations of Cayuga Lake. Flood profiles for the stream reaches were computed by combining a one-dimensional step-backwater model for the stream channels and a two-dimensional model for the overbank areas. The resulting hydraulic model was calibrated by using water-surface profiles from five observed storm events. The model was then used to compute 15 water-surface profiles for 5 flood frequencies (50-, 10-, 2-, 1-, and 0.2-percent annual exceedance probabilities, or 2-, 10-, 50-, 100-, and 500-year recurrence intervals) and 3 lake levels (representing average conditions, a 2-year-high condition, and a 100-year-high condition). The simulated water-surface profiles were then combined with a digital elevation model (derived from light detection and ranging data having 0.31‑foot vertical accuracy and 3.3-foot horizontal resolution) to delineate the maximum area flooded at each water level. Flood-inundation maps and geographic information system flood-extent polygons and depth grids are available in the data release associated with this report. These maps can provide emergency management personnel and residents with information that is critical for flood-management planning, flood-response activities, and postflood recovery efforts.

New York↗
USGS2019

Sand Creek characterization study for Oncorhynchus clarkii virginalis (Rio Grande Cutthroat Trout), Great Sand Dunes National Park and Preserve, Colorado

The Oncorhynchus clarkii virginalis (Rio Grande cutthroat trout, RGCT) has undergone extensive declines in distribution and population. The RGCT is the southernmost distributed subspecies of cutthroat trout. Native to the Rio Grande Basin in Colorado and New Mexico, the subspecies is also found in the headwaters of the Pecos River and Canadian River basins in New Mexico. Currently, RGCT populations represent approximately 12 percent of the historic distribution. There are many factors that have contributed to the decline of the RGCT including small population sizes; hybridization with non-native salmonids; competition with non-native salmonids; angling; and loss of habitat resulting from wildfire, stream drying, disease, increased water temperatures; and poor land management. The eastern side of Colorado’s Rio Grande Basin is also home to Great Sand Dunes National Park and Preserve and the Sand Creek watershed. This study was designed to (1) characterize current physical and biological conditions of waterbodies within the Sand Creek watershed, from headwaters to lower terminus near the dune field; (2) characterize the spatial extent of existing fisheries within the Sand Creek watershed to inform the scope of potential future reclamation efforts; and (3) evaluate key limiting factors for a future native RGCT reintroduction. Bathymetric profiles were completed for two lakes within the upper Sand Creek drainage to characterize the physical geometry of each lake and to estimate the total lake volume required for future piscicide treatment and (or) fish removal efforts. Physical and biological conditions evaluated included stream water temperature and intermittency, discharge, and the genetics and existing fish community distribution and composition within the Sand Creek watershed were key components of this study. A baseline established the geographic extent and biological constraints factored into future piscicide treatment planning and native trout reintroduction efforts. As a result of this work, the Sand Creek watershed can be broken up into several distinct categories: Lakes that are good candidates for reclamation and reintroduction of RGCT, lakes that are poor candidates for reclamation, streams that currently have fish and are good candidates for reclamation and reintroduction, streams that currently lack fish and may be good candidates for introduction of RGCT, and streams that currently lack fish and are not good candidates for introduction of RGCT. This characterization study report is intended to inform State and Federal managers of the likelihood that the Sand Creek watershed can support a sustainable population of RGCT should they be reintroduced.

Colorado↗
USGS2020

Storage capacity and sedimentation characteristics of the San Antonio Reservoir, California, 2018

The San Antonio Reservoir is a large water storage facility in Alameda County, California, and is a major component of the Hetch Hetchy Regional Water System (RWS). The RWS is a water-supply system owned and operated by the San Francisco Public Utilities Commission (SFPUC) and provides water for about 2.7 million people in the San Francisco, Santa Clara, Alameda, and San Mateo Counties. The San Antonio Reservoir is one of two RWS reservoirs in Alameda County and the third largest of the RWS reservoirs in the San Francisco Bay Area. The reservoir was formed by the James H. Turner Dam, which was completed in 1965. At the time of construction, the reservoir was estimated to have 50,500 acre-feet (acre-ft) of storage capacity. That early estimate was based on a 1963 pre-construction topographic map, which was drawn from aerial photographs. The capacity of the reservoir was later surveyed in 1994 and 2000. These two later surveys did not include the upper 18 feet (ft) of the reservoir, which represents roughly 30 percent of the overall storage volume. To determine the storage capacity and provide updated stage-capacity curves up to the spillway, the U.S. Geological Survey, in cooperation with the SFPUC, surveyed the bathymetry and shoreline of the reservoir in April 2018. The bathymetric survey was performed by making depth soundings using a boat-mounted, multibeam echosounder. At the time of the survey, the water level was between 13 and 14 ft below the spillway elevation. To measure capacity between the water line up to the spillway elevation, topography along most of the shoreline was surveyed from the boat using a terrestrial Light Detection and Ranging (LiDAR) scanner and in other areas by using ground-survey techniques. Location during bathymetric and topographic data collection was determined using a Global Navigation Satellite System-Real Time Network system. Vertical profiles of sound speed were collected periodically. The sound-speed profiles were used to spatially and temporally adjust the sound-speed calculations used to determine depth from the soundings. Approximately 125 kilometers (78 miles) of transects with a total of about 560 million depth soundings and topographic LiDAR points were collected (about 160 per square meter). In addition, approximately 500 topographic survey points were collected in shallow, wadable areas and on land near the upper reservoir area using a Global Navigation Satellite System receiver attached to a fixed length survey rod. Depth soundings, terrestrial LiDAR points, topographic survey points, and a digitized shoreline were merged and interpolated to generate a digital elevation model (DEM) of the reservoir. Gridded elevation data extracted from the DEM were then tabulated to determine total reservoir capacity and create reservoir stage-surface area and stage-storage capacity tables. Results of the reservoir capacity analysis indicated that the reservoir has 53,266 (plus or minus 140) acre-ft of storage capacity, which is an increase of 2,766 acre-ft (or 5.5 percent) greater than the original 1965 estimate; the increase is likely due to improved survey methods. Also, at the time of this 2018 survey, Intake #1 (the lowest intake) was not in operation. Intake #1 is estimated to be buried approximately 10 ft below the bed, whereas Intake #2 is about 20 ft above the bed. There are five intakes at different elevation levels; however, when consecutive lower intakes become inoperable due to sedimentation, the live storage capacity (capacity available for use) is reduced. At the time of this survey, the remaining live storage (above Intake #2) was approximately 52,363 acre-ft. The 2018 stage-capacity curve was compared to the original 1965 stage-capacity curve. Although overall, the changes indicate an increase in storage capacity, the change in volume at 372.7 ft North American Vertical Datum of 1988 (370 ft National Geodetic Vertical Datum of 1929, NGVD 29) shows a decrease of 733 acre-ft (the elevation of 370 ft NGVD 29 was used because it is the lowest elevation available for the 1965 stage-capacity curves). This finding agrees with the observed accumulation of sediment over Intake #1. That volume was converted to an annual sediment yield of 0.35 acre-ft per square mile (or 165 cubic meters per square kilometer), which is of the same order of magnitude as that found in other watersheds for the Coast Ranges in California. A decrease of 733 acre-ft between 1965 and 2018 thus represents a loss of 1.5 percent of the overall storage capacity in the reservoir. The updated stage-surface area and stage-capacity tables provided in this report and online ( https://doi.org/10.5066/P9KC9DU8 ) can be used by the SFPUC to improve reservoir operations and serve as an accurate baseline to monitor bathymetric changes in the future.

California↗
USGS2024

Development of a hydrogeologic visualization model for western Sarpy County, Nebraska

Population in western Sarpy County, Nebraska, has steadily increased over the last several decades and has led to increased groundwater use for domestic purposes. To meet the increase in demand, the Papio-Missouri River Natural Resources District is seeking to use all available sources of groundwater in western Sarpy County. Additionally, elevated groundwater nitrate plus nitrite as nitrogen concentrations were detected, indicating the need to better understand the groundwater quality of the area. Although the general geology of the area is understood, the area does not have detailed information on the extent of the various aquifers, particularly the Dakota aquifer. To characterize these aquifers, the Papio-Missouri River Natural Resources District invested in airborne electromagnetic surveys of the area to better understand the subsurface geology. Although these surveys improved understanding of the groundwater systems in the area, the Papio-Missouri River Natural Resources District wanted to integrate the subsurface information with available water-quality and groundwater-level data. In response, the U.S. Geological Survey, in cooperation with the Papio-Missouri River Natural Resources District, the Nebraska Natural Resources Commission, and the Nebraska Department of Natural Resources, assembled geologic, hydrogeologic and nitrate plus nitrite as nitrogen information for the selected area into a three-dimensional visualization computer software package called GeoScene3D. The completed GeoScene3D project was assembled to provide a visualization of the groundwater systems and associated water-quality results in Sarpy County and to provide the Papio-Missouri River Natural Resources District managers with information that can be used to make more informed groundwater resource-planning decisions in the future. This report details the development of a three-dimensional model created within GeoScene3D to visualize the subsurface, particularly the Dakota Sandstone in western Sarpy County.

Nebraska↗
USGS2024

Seismic techniques and suggested instrumentation to monitor volcanoes

Introduction Changes in the pressure or location of magma can stress or break surrounding rocks and trigger flow of nearby waters and gases, causing seismic signals, such as discrete earthquakes and tremor. These phenomena are types of seismic unrest that commonly precede eruption and can be used to forecast volcanic activity. Mass movements at the surface, including avalanches, debris flows, and lahars, may also generate seismic signals that are specifically addressed in chapter H , this volume (Thelen and others, 2024). Our focus in this chapter is to determine the levels of instrumentation recommended to produce high-quality, well-constrained seismic observations important for early warning of impending eruptions, detecting changes in ongoing eruptions, and characterizing other hazardous volcanic events. There are emerging techniques and new types of instrumentation, such as distributed acoustic sensing or rotational seismometers, that we do not consider here. These types of instrumentation show promise for monitoring but still require maturation before being considered more generally in volcano monitoring. Most of the capabilities mentioned below are universal for all types of volcanic systems, although some are best applied to stratovolcanoes with an apical single vent. In some settings, such as calderas or shield volcanoes, we must broaden coverage to include multiple possible storage regions or vent locations. As an example, Thelen (2014) discretized the long rift zones of shield volcanoes in Hawaiʻi as a set of evenly spaced “vents.” In this construct, each vent comes with recommendations, and several thousand network configurations were simulated to assess the effect on network quality levels and to determine the most efficient network design. The same process could be applied in a caldera setting or a volcanic field, where an evenly spaced grid of potential vents is considered. Localized recommendations for each unique system are beyond the scope of this report and left up to local experts to assess based on the conditions, restrictions, and requirements of each volcano.

Scientific Investigations Report↗
USGS2004

Seismic monitoring at Cascade volcanic centers, 2004 — Status and recommendations

The purpose of this report is to assess the current (May, 2004) status of seismic monitoring networks at the 13 major Cascade volcanic centers. Included in this assessment are descriptions of each network, analyses of the ability of each network to detect and to locate seismic activity, identification of specific weaknesses in each network, and a prioritized list of those networks that are most in need of additional seismic stations. At the outset it should be recognized that no Cascade volcanic center currently has an adequate seismic network relative to modern-day networks at Usu Volcano (Japan) or Etna and Stromboli volcanoes (Italy). For a system the size of Three Sisters, for example, a modern-day, cutting-edge seismic network would ideally consist of a minimum of 10 to 12 short-period three-component seismometers (for determining particle motions, reliable S-wave picks, moment tensor inversions, fault-plane solutions, and other important seismic parameters) and 7 to 10 broadband sensors (which, amongst other considerations, enable detection and location of very long period (VLP) and other low-frequency events, moment tensor inversions, and, because of their wide dynamic range, on-scale recording of large-amplitude events). Such a dense, multi component seismic network would give the ability to, for example, detect in near-real-time earthquake migrations over a distance of ~0.5km or less, locate tremor sources, determine the nature of a seismic source (that is, pure shear, implosive, explosive), provide on-scale recordings of very small and very large-amplitude seismic signals, and detect localized changes in seismic stress tensor orientations caused by movement of magma bodies. However, given that programmatic resources are currently limited, installation of such networks at this time is unrealistic. Instead, this report focuses on identifying what additional stations are needed to guarantee that anomalous seismicity associated with volcanic unrest will be detected in a timely manner and, in the case of magnitude = 1 earthquakes, reliably located.

California, Nevada, Oregon, Washington↗
USGS2009

Regression equations for estimation of annual peak-streamflow frequency for undeveloped watersheds in Texas using an L-moment-based, PRESS-minimized, residual-adjusted approach

Annual peak-streamflow frequency estimates are needed for flood-plain management; for objective assessment of flood risk; for cost-effective design of dams, levees, and other flood-control structures; and for design of roads, bridges, and culverts. Annual peak-streamflow frequency represents the peak streamflow for nine recurrence intervals of 2, 5, 10, 25, 50, 100, 200, 250, and 500 years. Common methods for estimation of peak-streamflow frequency for ungaged or unmonitored watersheds are regression equations for each recurrence interval developed for one or more regions; such regional equations are the subject of this report. The method is based on analysis of annual peak-streamflow data from U.S. Geological Survey streamflow-gaging stations (stations). Beginning in 2007, the U.S. Geological Survey, in cooperation with the Texas Department of Transportation and in partnership with Texas Tech University, began a 3-year investigation concerning the development of regional equations to estimate annual peak-streamflow frequency for undeveloped watersheds in Texas. The investigation focuses primarily on 638 stations with 8 or more years of data from undeveloped watersheds and other criteria. The general approach is explicitly limited to the use of L-moment statistics, which are used in conjunction with a technique of multi-linear regression referred to as PRESS minimization. The approach used to develop the regional equations, which was refined during the investigation, is referred to as the 'L-moment-based, PRESS-minimized, residual-adjusted approach'. For the approach, seven unique distributions are fit to the sample L-moments of the data for each of 638 stations and trimmed means of the seven results of the distributions for each recurrence interval are used to define the station specific, peak-streamflow frequency. As a first iteration of regression, nine weighted-least-squares, PRESS-minimized, multi-linear regression equations are computed using the watershed characteristics of drainage area, dimensionless main-channel slope, and mean annual precipitation. The residuals of the nine equations are spatially mapped, and residuals for the 10-year recurrence interval are selected for generalization to 1-degree latitude and longitude quadrangles. The generalized residual is referred to as the OmegaEM parameter and represents a generalized terrain and climate index that expresses peak-streamflow potential not otherwise represented in the three watershed characteristics. The OmegaEM parameter was assigned to each station, and using OmegaEM, nine additional regression equations are computed. Because of favorable diagnostics, the OmegaEM equations are expected to be generally reliable estimators of peak-streamflow frequency for undeveloped and ungaged stream locations in Texas. The mean residual standard error, adjusted R-squared, and percentage reduction of PRESS by use of OmegaEM are 0.30log 10 , 0.86, and -21 percent, respectively. Inclusion of the OmegaEM parameter provides a substantial reduction in the PRESS statistic of the regression equations and removes considerable spatial dependency in regression residuals. Although the OmegaEM parameter requires interpretation on the part of analysts and the potential exists that different analysts could estimate different values for a given watershed, the authors suggest that typical uncertainty in the OmegaEM estimate might be about +or-0.10 10 . Finally, given the two ensembles of equations reported herein and those in previous reports, hydrologic design engineers and other analysts have several different methods, which represent different analytical tracks, to make comparisons of peak-streamflow frequency estimates for ungaged watersheds in the study area.

Texas↗
USGS2026

Slip tendency and dilation tendency of Quaternary faults in the Great Basin region, Arizona, California, Idaho, Nevada, Oregon, Utah, and Wyoming

This report presents the results of slip tendency and dilation tendency analysis of Quaternary faults in the Great Basin region, in the states of Arizona, California, Idaho, Nevada, Oregon, Utah, and Wyoming. This effort was undertaken as part of two studies: the U.S. Department of Energy Geothermal Technologies Office-funded Innovative Geothermal Exploration through Novel Investigations of Undiscovered Systems (INGENIOUS) project, which is a multidisciplinary, multi-institution effort to develop new methodologies and best practices to accelerate the discovery of new, commercially viable hydrothermal resources; and the U.S. Geological Survey Geothermal Resources Investigations Project. The goal of this work is to identify faults and fault segments that are appropriately oriented to be stress-loaded for slip or to dilate under the ambient stress conditions. Both conditions may make such faults more likely to host permeability and fluid flow relative to other faults in the region. Slip tendency and dilation tendency were calculated for all fault segments in the study area from the U.S. Geological Survey Quaternary Fault and Fold Database. Results demonstrate that faults with high slip tendency and dilation tendency are common in the Great Basin region. The areas containing many high slip tendency and dilation tendency fault segments are associated with high strain rates. These areas also host a relatively high proportion of the known hydrothermal systems. This indicates that faults with high slip tendency and dilation tendency in areas with high strain-rate are potentially most prospective for hosting hydrothermal systems. The database published in coordination with this report contains the slip tendency and dilation tendency results data, including the fault geometry and stress data used to calculate these quantities, for both the Great Basin study area and the smaller INGENIOUS study area.

Arizona, California, Idaho, Nevada, Oregon, Utah, ↗
USGS2008

A regression model for computing index flows describing the median flow for the summer month of lowest flow in Michigan

In 2006, Michigan enacted laws to prevent new large capacity withdrawals from decreasing flows to the extent that they would functionally impair a stream's ability to support characteristic fish populations. The median streamflow for the summer month of lowest flow was specified by state decision makers as the index flow on which likely impacts of withdrawals would be assessed. At sites near long-term streamflow-gaging stations, analysis of streamflow records during July, August, and September was used to determine the index flow. At ungaged sites, an alternate method for computing the index flow was needed. This report documents the development of a method for computing index flows at ungaged stream sites in Michigan. The method is based on a regression model that computes the index water yield, which is the index flow divided by the drainage area. To develop the regression model, index flows were determined on the basis of daily flows measured during July, August, and September at 147 streamflow-gaging stations having 10 or more years of record (considered long-term stations) in Michigan. The corresponding index water yields were statistically related to climatic and basin characteristics upstream from the stations in the regression model. Climatic and basin characteristics selected as explanatory variables in the regression model include two aquifer-transmissivity and hydrologic-soil groups, forest land cover, and normal annual precipitation. Regression model estimates of water yield explain about 70.8 percent of the variability in index water yields indicated by streamflow-gaging station records. Index flows computed on the basis of regression-model estimates of water yield and corresponding drainage areas explain about 94.0 percent of the variability in index flows indicated by streamflow-gaging station records. No regional bias was detected in the regression-based estimates of water yield within seven hydrologic subregions spanning Michigan. Thus, the single regression model developed in this report can be used to produce unbiased estimates of index water yield and flow statewide. In addition, a technique is presented for computing prediction intervals about the index flow estimates.

Michigan↗
USGS2013

Potential depletion of surface water in the Colorado River and agricultural drains by groundwater pumping in the Parker-Palo Verde-Cibola area, Arizona and California

Water use along the lower Colorado River is allocated as “consumptive use,” which is defined to be the amount of water diverted from the river minus the amount that returns to the river. Diversions of water from the river include surface water in canals and water removed from the river by pumping wells in the aquifer connected to the river. A complication in accounting for water pumped by wells occurs if the pumping depletes water in drains and reduces measured return flow in those drains. In that case, consumptive use of water pumped by the wells is accounted for in the reduction of measured return flow. A method is needed to understand where groundwater pumping will deplete water in the river and where it will deplete water in drains. To provide a basis for future accounting for pumped groundwater in the Parker-Palo Verde-Cibola area, a superposition model was constructed. The model consists of three layers of finite-difference cells that cover most of the aquifer in the study area. The model was run repeatedly with each run having a pumping well in a different model cell. The source of pumped water that is depletion of the river, expressed as a fraction of the pumping rate, was computed for all active cells in model layer 1, and maps were constructed to understand where groundwater pumping depletes the river and where it depletes drains. The model results indicate that if one or more drains exist between a pumping well location and the river, nearly all of the depletion will be from drains, and little or no depletion will come from the Colorado River. Results also show that if a well pumps on a side of the river with no drains in the immediate area, depletion will come from the Colorado River. Finally, if a well pumps between the river and drains that parallel the river, a fraction of the pumping will come from the river and the rest will come from the drains. Model results presented in this report may be considered in development or refinement of strategies for accounting for groundwater pumping in the river aquifer connected to the Colorado River in the study area.

Arizona;California↗
USGS2017

Groundwater flux and nutrient loading in the northeast section of Bear Lake, Muskegon County, Michigan, 2015

Bear Lake in North Muskegon, Michigan, is listed as part of the Muskegon Lake area of concern as designated by the U.S. Environmental Protection Agency. This area of concern was designated as a result of eutrophication and beneficial use impairments. On the northeast end of Bear Lake, two man-made retention ponds (Willbrandt Pond East and Willbrandt Pond West), formerly used for celery farming, may contribute nutrients to Bear Lake. Willbrandt Ponds (East and West) were previously muck fields that were actively used for celery farming from the early 1900s until 2002. The restoration and reconnection of the Willbrandt Ponds into Bear Lake prompted concerns of groundwater nutrient loading into Bear Lake. Studies done by the State of Michigan and Grand Valley State University revised initial internal phosphorus load estimates and indicated an imbalance in the phosphorus budget in Bear Lake. From June through November 2015, the U.S. Geological Survey (USGS) did an investigative study to quantify the load of nutrients from shallow groundwater around the Willbrandt Ponds in an effort to update the phosphorus budget to Bear Lake. Seven sampling locations were established, including five shallow groundwater wells and two surface-water sites, in the Willbrandt pond study area and Bear Lake. A total of 12 nutrient samples and discrete water-level measurements were collected from each site from June through November 2015. Continuous water-level data were recorded for both surface-water monitoring locations for the entire sampling period. Water-level data indicated that Willbrandt Pond West had the highest average water-level elevation of all sites monitored, which indicated the general direction of flux is from Willbrandt Pond West to Bear Lake. Nutrient and chloride loading from Willbrandt Pond West to Bear Lake was calculated using two distinct methods: Dupuit and direct seepage methods. Shallow groundwater loading calculations were determined by using groundwater levels to first determine a flux of shallow groundwater, then nutrient concentrations to determine a load. It was determined that Willbrandt Pond East and Willbrandt Pond West contributed between 2 to 4 percent of the total annual phosphorus load to Bear Lake by way of shallow groundwater flow. Annual loads calculated for other constituents include orthophosphate (40–100 pounds per year [lb P/yr]), total nitrogen (200–830 lb/yr), chloride (12,700–32,100 lb/yr), and ammonia (130–670 lb N/yr). Study results indicated that mean groundwater and surface-water nutrient concentrations calculated in this study were higher than reported Michigan statewide values. The data collected in this study allow understanding of groundwater nutrient loading into Bear Lake in an effort to help inform future restoration and management decisions.

Michigan↗
USGS2018

Method comparisons for determining concentrations of metals in water samples used in studies of fish migratory histories

Signatures developed from metal concentrations in water and fish bony structures can be used to demonstrate migration of individual fish between connected water bodies. The U.S. Geological Survey (USGS), in cooperation with the National Park Service and the Missouri Department of Conservation, compared two protocols for collecting and analyzing water samples for concentrations of several metals commonly used to develop metal signatures. In 2015, paired seasonal water samples were collected in two study areas incorporating three National Park Service units; paired water samples were collected using USGS protocols and simpler research protocols. Metal concentrations obtained using USGS and research protocols were compared using t-tests, percent differences, and simple linear regression analyses. Graphical plots of median values and measured ranges were used to compare ratios of strontium to calcium (Sr:Ca) and barium to calcium (Ba:Ca) obtained using the different protocols among individual stations within the two study areas. For stations on the Mississippi and St. Croix Rivers, ranges in concentrations of calcium, barium, and strontium (obtained using USGS protocols) were compared between samples collected from 1995 through 2012 and samples collected in this study. Comparisons were used to evaluate the long-term stability of metal concentrations in the environment. Collectively, results presented in this report demonstrated that research protocols provided metal concentration data that were similar to data obtained using USGS protocols for all compared metals except manganese. Holding times of 6–33 weeks prior to filtration and analyses for samples collected using research protocols may have caused greater changes in manganese concentrations compared to other metals. Strontium, barium, and calcium are the metals most commonly used in studies of fish migration, and concentrations of these metals were similar using different protocols. However, rivers within each study area were more easily distinguished from each other using metal concentration data obtained using USGS protocols compared to data obtained using research protocols. Information presented in this report can be used to develop studies that use identified metal signatures in connected water bodies and bony fish structures to demonstrate fish migration.

Mississippi National River and Recreation Area, Na↗
USGS2022

Hydrologic and geochemical characterization of the Petaluma River watershed, Sonoma County, California

Executive Summary The objectives of the study are to (1) develop an updated assessment of the hydrogeology and geochemistry of the Petaluma valley watershed (PVW) and (2) develop an integrated hydrologic model for the PVW. The purpose of this report is to describe the conceptual model of the hydrologic, hydrogeologic, and water-quality characteristics of the PVW and a numerical groundwater-flow model of PVW.

California↗
USGS2007

Natural gases in ground water near Tioga Junction, Tioga County, north-central Pennsylvania: Occurrence and use of isotopes to determine origins, 2005

In January 2001, State oil and gas inspectors noted bubbles of natural gas in well water during a complaint investigation near Tioga Junction, Tioga County, north-central Pa. By 2004, the gas occurrence in ground water and accumulation in homes was a safety concern; inspectors were taking action to plug abandoned gas wells and collect gas samples. The origins of the natural-gas problems in ground water were investigated by the U.S. Geological Survey, in cooperation with the Pennsylvania Department of Environmental Protection, in wells throughout an area of about 50 mi2, using compositional and isotopic characteristics of methane and ethane in gas and water wells. This report presents the results for gas-well and water-well samples collected from October 2004 to September 2005. Ground water for rural-domestic supply and other uses near Tioga Junction is from two aquifer systems in and adjacent to the Tioga River valley. An unconsolidated aquifer of outwash sand and gravel of Quaternary age underlies the main river valley and extends into the valleys of tributaries. Fine-grained lacustrine sediments separate shallow and deep water-bearing zones of the outwash. Outwash-aquifer wells are seldom deeper than 100 ft. The river-valley sediments and uplands adjacent to the valley are underlain by a fractured-bedrock aquifer in siliciclastic rocks of Paleozoic age. Most bedrock-aquifer wells produce water from the Lock Haven Formation at depths of 250 ft or less. A review of previous geologic investigations was used to establish the structural framework and identify four plausible origins for natural gas. The Sabinsville Anticline, trending southwest to northeast, is the major structural feature in the Devonian bedrock. The anticline, a structural trap for a reservoir of deep native gas in the Oriskany Sandstone (Devonian) (origin 1) at depths of about 3,900 ft, was explored and tapped by numerous wells from 1930-60. The gas reservoir in the vicinity of Tioga Junction, depleted of native gas, was converted to the Tioga gas-storage field for injection and withdrawal of non-native gases (origin 2). Devonian shale gas (shallow native gas) also has been reported in the area (origin 3). Gas might also originate from microbial degradation of buried organic material in the outwash deposits (origin 4). An inventory of combustible-gas concentrations in headspaces of water samples from 91 wells showed 49 wells had water containing combustible gases at volume fractions of 0.1 percent or more. Well depth was a factor in the observed occurrence of combustible gas for the 62 bedrock wells inventoried. As well-depth range increased from less than 50 ft to 51-150 ft to greater than 151 ft, the percentage of bedrock-aquifer wells with combustible gas increased. Wells with high concentrations of combustible gas occurred in clusters; the largest cluster was near the eastern boundary of the gas-storage field. A subsequent detailed gas-sampling effort focused on 39 water wells with the highest concentrations of combustible gas (12 representing the outwash aquifer and 27 from the bedrock aquifer) and 8 selected gas wells. Three wells producing native gas from the Oriskany Sandstone and five wells (two observation wells and three injection/withdrawal wells) with non-native gas from the gas-storage field were sampled twice. Chemical composition, stable carbon and hydrogen isotopes of methane (13CCH4 and DCH4), and stable carbon isotopes of ethane (13CC2H6) were analyzed. No samples could be collected to document the composition of microbial gas originating in the outwash deposits (outwash or 'drift' gas) or of native natural gas originating solely in Devonian shale at depths shallower than the Oriskany Sandstone, although two of the storage-field observation wells sampled reportedly yielded some Devonian shale gas. Literature values for outwash or 'drift' gas and Devonian shale gases were used to supplement the data collection. Non-native gases fr

Pennsylvania↗
USGS2010

Summary of Hydrologic Data for the Tuscarawas River Basin, Ohio, with an Annotated Bibliography

The Tuscarawas River Basin drains approximately 2,600 square miles in eastern Ohio and is home to 600,000 residents that rely on the water resources of the basin. This report summarizes the hydrologic conditions in the basin, describes over 400 publications related to the many factors that affect the groundwater and surface-water resources, and presents new water-quality information and a new water-level map designed to provide decisionmakers with information to assist in future data-collection efforts and land-use decisions. The Tuscarawas River is 130 miles long, and the drainage basin includes four major tributary basins and seven man-made reservoirs designed primarily for flood control. The basin lies within two physiographic provinces-the Glaciated Appalachian Plateaus to the north and the unglaciated Allegheny Plateaus to the south. Topography, soil types, surficial geology, and the overall hydrology of the basin were strongly affected by glaciation, which covered the northern one-third of the basin over 10,000 years ago. Within the glaciated region, unconsolidated glacial deposits, which are predominantly clay-rich till, overlie gently sloping Pennsylvanian-age sandstone, limestone, coal, and shale bedrock. Stream valleys throughout the basin are filled with sands and gravels derived from glacial outwash and alluvial processes. The southern two-thirds of the basin is characterized by similar bedrock units; however, till is absent and topographic relief is greater. The primary aquifers are sand- and gravel-filled valleys and sandstone bedrock. These sands and gravels are part of a complex system of aquifers that may exceed 400 feet in thickness and fill glacially incised valleys. Sand and gravel aquifers in this basin are capable of supporting sustained well yields exceeding 1,000 gallons per minute. Underlying sandstones within 300 feet of the surface also provide substantial quantities of water, with typical well yields of up to 100 gallons per minute. Although hydraulic connection between the sandstone bedrock and the sands and gravels in valleys is likely, it has not been assessed in the Tuscarawas River Basin. In 2001, the major land uses in the basin were approximately 40 percent forested, 39 percent agricultural, and 17 percent urban/residential. Between 1992 and 2001, forested land use decreased by 2 percent with correspondingly small increases in agricultural and urban land uses, but from 1980 to 2005, the 13-county area that encompasses the basin experienced a 7.1-percent increase in population. Higher population density and percentages of urban land use were typical of the northern, headwaters parts of the basin in and around the cities of Akron, Canton, and New Philadelphia; the southern area was rural. The basin receives approximately 38 inches of precipitation per year that exits the basin through evapotranspiration, streamflow, and groundwater withdrawals. Recharge to groundwater is estimated to range from 6 to 10 inches per year across the basin. In 2000, approximately 89 percent of the 116 million gallons per day of water used in the basin came from groundwater sources, whereas 11 percent came from surface-water sources. To examine directions of groundwater flow in the basin, a new dataset of water-level contours was developed by the Ohio Department of Natural Resources. The contours were compiled on a map that shows that groundwater flows from the uplands towards the valleys and that the water-level surface mimics surface topography; however, there are areas where data were too sparse to adequately map the water-level surface. Additionally, little is known about deep groundwater that may be flowing into the basin from outside the basin and groundwater interactions with surface-water bodies. Many previous reports as well as new data collected as part of this study show that water quality in the streams and aquifers in the Tuscarawas River Basin has been degraded by urban, suburban, and rural

Scientific Investigations Report↗