Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Data Series”

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 325 records · Page 18Linked to original sources

Growth history of oil reserves in major California oil fields during the twentieth century

Oil reserves in 12 of California's 52 giant fields (fields with estimated recovery > 100 million barrels of oil) have continued to appreciate well past the age range at which most fields cease to show significant increases in ultimate recovery. Most of these fields were discovered between 1890 and 1920 and grew to volumes greater than 500 million barrels in their first two decades. Growth of reserves in these fields accelerated in th e1950s and 1960s and is mostly explained by application of secondary and tertiary recovery technicques, primarily waterflooding and thermal recovery. The remaining three-fourths of California's giant fields show a pattern of growth in which fields cease to grow significantly by 20-30 years following recovery. virtually all of these fields have estimated ultimate recoveries less than about 500 million barrels and most are in the 100-200 million barrel range. Three of six offshore giant fields, all discovered between 1966 and 1981, have shown decreases in their estimated ultimate sizes within about the first decade after production began, presumably because production volumes ailed to match initial projections. The data suggest that: 1. Only fields that attain an estimated ultimate size of several hundred million barrels shortly after discovery and have geologic characterisics that make them suceptible to advanced recovery techniques are likely to show substantial late growth. 2. Offshore fields are less likely to show significant growth, probably because projections based on modern seismic reflection and reservoir test data are unlikely to underestimate the volume of oil in the field. 3. Secondary and tertiary recovery programs rather than field extensions or new pool discoveries are responsible for most of the significant growth of reserves in California. 4. field size data collected ove rmany decades provide a more comprehensive context for inferring reasons for reserve appreciation than shorter data series such as the Oil and Gas Integrated Field file (OGIFF) from the U.S. Department of Energy's Energy Information Administration (EIA). 5. Efforts to project future growth in California fields, and perhaps fields in other regions, should focus on evaulating the potential for enhanced recovery in fields with current estimated ultimate recoveries of about 250-500 million barrels. 6. By analogy with oil, attempts to project growth in gas reservoirs, in California and perhaps elsewhere, should focus on larger fields with lower permeability reservoirs where advances in recovery technology, such as perhaps horizontal drilling, are more likely to add substantial reserves.

Bulletin↗

Historical Ice-Out Dates for 29 Lakes in New England

Historical ice-out dates for 29 lakes in New England were compiled and are presented in this report. The length of record for the lakes ranges from 64 to 163 years, with an average of 108 years. Many lakes in New England had their latest recorded ice-out date in 1888. Ice-out dates for lakes are an important hydrologic data series for climate researchers and other interested parties.

Open-File Report↗

Hydrologic and landscape database for the Cache and White River National Wildlife Refuges and contributing watersheds in Arkansas, Missouri, and Oklahoma

A hydrologic and landscape database was developed by the U.S. Geological Survey, in cooperation with the U.S. Fish and Wildlife Service, for the Cache River and White River National Wildlife Refuges and their contributing watersheds in Arkansas, Missouri, and Oklahoma. The database is composed of a set of ASCII files, Microsoft Access® files, Microsoft Excel® files, an Environmental Systems Research Institute (ESRI) ArcGIS® geodatabase, ESRI ArcGRID® raster datasets, and an ESRI ArcReader® published map. The database was developed as an assessment and evaluation tool to use in examining refuge-specific hydrologic patterns and trends as related to water availability for refuge ecosystems, habitats, and target species; and includes hydrologic time-series data, statistics, and hydroecological metrics that can be used to assess refuge hydrologic conditions and the availability of aquatic and riparian habitat. Landscape data that describe the refuge physiographic setting and the locations of hydrologic-data collection stations are also included in the database. Categories of landscape data include land cover, soil hydrologic characteristics, physiographic features, geographic and hydrographic boundaries, hydrographic features, regional runoff estimates, and gaging-station locations. The database geographic extent covers three hydrologic subregions—the Lower Mississippi–St Francis (0802), the Upper White (1101), and the Lower Arkansas (1111)—within which human activities, climatic variation, and hydrologic processes can potentially affect the hydrologic regime of the refuges and adjacent areas. Database construction has been automated to facilitate periodic updates with new data. The database report (1) serves as a user guide for the database, (2) describes the data-collection, data-reduction, and data-analysis methods used to construct the database, (3) provides a statistical and graphical description of the database, and (4) provides detailed information on the development of analytical techniques designed to assess water availability for ecological needs.

Arkansas;Missouri;Oklahoma↗

Stratigraphic cross section of measured sections and drill holes of the Neslen Formation and adjacent formations, Book Cliffs Area, Colorado and Utah

This study updates a stratigraphic cross section published as plate 2 in Kirschbaum and Hettinger (2004) Digital Data Series 69-G ( http://pubs.usgs.gov/dds/dds-069/dds-069-g/ ). The datum is a marine/tidal ravinement surface within the Cozzette Sandstone Member of the Iles Formation and the Thompson Canyon Sandstone and Sulphur Canyon Sandstone Beds of the Neslen Formation. One of the cores shown was included on the original cross section, and new core descriptions have been added to the upper part of the cored interval. A new core description (S178) is included in this report. Cores are stored in the U.S. Geological Survey Core Research Facility at the Denver Federal Center, Colorado. The following information has also been added to help define the stratigraphic framework: 1) At least five claystones interpreted as altered volcanic ashes have been identified and may give future workers a correlation tool within the largely continental section. 2) Thickness and general geometry of the Sego Sandstone, Buck Tongue of the Mancos Shale, and Castlegate Sandstone have been added to provide additional stratigraphic context. 3) The geometry in the Sego Sandstone, Buck Tongue of the Mancos Shale, and Castlegate Sandstone has been added to provide additional stratigraphic context. 4) Ammonite collections are from Gill and Hail. The zone of Didymoceras nebrascense projected into the East Salt Wash area is based on correlation of the flooding surface at the base of the Cozzette Member to this point as shown in Kirschbaum and Hettinger. 5) A leaf locality of the Denver Museum of Nature and Science is shown in its approximate stratigraphic position near Thompson Canyon. 6) A dinosaur locality of the Natural History Museum of Utah is shown in the Horse Canyon area measured section at the stratigraphic position where it was extracted.

Utah, Colorado↗

Weekly resolution particulate flux from a sediment trap in the northern Gulf of Mexico, 2008-2012

The U.S. Geological Survey anchored a sediment trap in the northern Gulf of Mexico to collect time-series data on sediment flux from 2008 to 2012. There are continuous measurements of total mass flux and organic carbon flux (ogC) at 7–14 day resolution from 2008 to 2012. The flux of calcium carbonate (CaCO 3 ), particulate nitrogen (nitro), and biogenic silica (Opal) were also measured from January-December, 2008. The mass flux ranged from 0.01 g m -2 day -1 (grams per square meter per day) to 2.50 g m -2 day -1 , with a mean mass flux of 0.20 g m -2 day -1 over the 5-year study period.

Gulf Of Mexico↗

Future wave and wind projections for United States and United-States-affiliated Pacific Islands

Changes in future wave climates in the tropical Pacific Ocean from global climate change are not well understood. Spatially and temporally varying waves dominate coastal morphology and ecosystem structure of the islands throughout the tropical Pacific. Waves also impact coastal infrastructure, natural and cultural resources, and coastal-related economic activities of the islands. Wave heights, periods, and directions were forecast through the year 2100 using wind parameter outputs from four atmosphere-ocean global climate models from the Coupled Model Inter-Comparison Project, Phase 5, for Representative Concentration Pathways (RCP) scenarios 4.5 and 8.5 that correspond to moderately mitigated and unmitigated greenhouse gas emissions, respectively. Wind fields from the global climate models were used to drive a global WAVEWATCH-III wave model and generate hourly time-series of bulk wave parameters for 25 islands in the mid to western tropical Pacific for the years 1976–2005 (historical), 2026–2045 (mid-century projection), and 2085–2100 (end-of-century projection). Although the results show some spatial heterogeneity, overall the December-February extreme significant wave heights, defined as the mean of the top 5 percent of significant wave height time-series data modeled within a specific period, increase from present to mid-century and then decrease toward the end of the century; June-August extreme wave heights increase throughout the century within the Central region of the study area; and September-November wave heights decrease strongly throughout the 21st century, displaying the largest and most widespread decreases of any season. Peak wave periods increase east of the International Date Line during the December-February and June-August seasons under RCP4.5. Under the RCP8.5 scenario, wave periods decrease west of the International Date Line during December-February but increase in the eastern half of the study area. Otherwise, wave periods decrease throughout the study area during other seasons. Extreme wave directions in equatorial Micronesia during June-August undergo an approximate 30° clockwise rotation from primarily west to northwest. September-November RCP4.5 extreme mean wave directions rotate counterclockwise by approximately 30 to 45° in equatorial Micronesia; September-November RCP8.5 extreme mean wave directions within equatorial Micronesia rotate clockwise by approximately 20 to 30°. Extreme wind speeds decreased within both scenarios, with the largest decreases occurring in the September-November season. Extreme wind directions under RCP4.5 rotated clockwise by more than 60° in equatorial Micronesia during the September-November season and by approximately 30° during June-August. RCP8.5 extreme wind directions rotated counterclockwise during September-November within the same region by 30 to 50° and clockwise by 30 to 40° at one island. The spatial patterns and trends are similar between the two different greenhouse gas emission scenarios, with the magnitude and extent of the trends generally greater for the higher (RCP8.5) scenario.

Micronesia, Pacific Ocean↗

Summary of oceanographic and water-quality measurements in Chincoteague Bay, Maryland and Virginia, 2014–15

U.S. Geological Survey scientists and technical support staff measured oceanographic, waterquality, seabed-elevation-change, and meteorological parameters in Chincoteague Bay, Maryland and Virginia, during the period of August 13, 2014, to July 14, 2015, as part of the Estuarine Physical Response to Storms project (GS2–2D) supported by the Department of the Interior Hurricane Sandy recovery program. These measurements provide time series data that quantify the response and can be used to better understand the resilience of this back-barrier estuarine system to storms. The Assateague Island National Seashore (National Park Service) and the Chincoteague National Wildlife Refuge (U.S. Fish and Wildlife Service) are on the east side of Chincoteague Bay.

Maryland, Virginia↗

World petroleum assessment 2000; compiled PowerPoint slides

The slides in this compilation have been produced for a number of presentations on the World Petroleum Assessment 20000. Many of the figures are taken directly form the publication "U.S. Geological Survey World Petroleum Assessment 2000" - Description and Results: USGS Digital Data Series DDS-60, 2000. Some of the slides are modifications of figures from DDS-60, some are new descriptive slides, and a few are new slides. Several of the slides appear to be duplicates, but in fact are slight modifications for format or content from the same image. Forty-one people participated in this effort as part of the World Energy Assessment Team. The full list of contributors is given ion DDS-60.

Open-File Report↗

Simulation of Streamflow and Water Quality to Determine Fecal Coliform and Nitrate Concentrations and Loads in the Mad River Basin, Ohio

The Hydrological Simulation Program Fortran (HSPF) was used to simulate the concentrations and loads of fecal coliform and nitrate for streams in the Mad River Basin in west-central Ohio during the period 1999 through 2003. The Mad River Basin was divided into subbasins that were defined either by the 14-digit Hydrologic Unit (HU) boundaries or by streamflow-gaging-station locations used in the model. Model calibration and simulation processes required the formation of nine meteorologic zones to input meteorologic time-series data and water-quality data. Sources of fecal coliform and nitrate from wastewater-treatment discharges and combined sewer overflow discharges (CSOs) within the City of Springfield were point sources simulated in the model. Failing septic systems and cattle with direct access to streams were nonpoint sources included in the study but treated in the model as point sources. Other nonpoint sources were addressed by adjusting interflow and ground-water concentrations in the subsurface and maximum storage capacities and accumulation rates of the simulated constituents on the land surface for each meteorologic zone. Simulation results from the calibrated model show that several HUs exceeded the water-quality standard of 1,000 colony-forming units per 100 mL for fecal coliform based on the maximum 30-day geometric mean. Most HUs with high fecal coliform counts were within or downstream from the City of Springfield. No water-quality standard has been set for instream nitrate concentrations; however, the Ohio Environmental Protection Agency (Ohio EPA) considered a concentration of 5 mg/L or greater to be of concern. Simulation results indicate that several HUs in the agricultural areas of the basin exceeded this level. The calibrated model was modified to create scenarios that simulated loads of fecal coliform and nitrate that were either reduced or eliminated from selected sources. The revised models included the elimination of failing septic systems, elimination of direct access of cattle to streams, decrease in fecal coliform loads from the CSOs and selected wastewater-treatment facilities, and decrease in nitrate loads from land surfaces. The fecal coliform source-reduction model decreased the fecal coliform concentrations below a target concentration of 1,000 colonies per 100 milliliters for all HU outlets and decreased the load at the mouth of the Mad River by 73 percent. The nitrate source-reduction model decreased some HU mean concentrations to 5 milligrams per liter or less and decreased the load at the mouth of the Mad River by 52 percent. Other reduction scenarios may be run by Ohio EPA with the intent of identifying a management strategy that will attain a target concentration for the Mad River Basin.

Scientific Investigations Report↗

Flood plain delineation for the Fremont River and Bull Creek, Hanksville, Utah

Predicted inundation maps for the Fremont River and Bull Creek in Hanksville, Utah, were developed using one-dimensional hydraulic models. Estimates of the 1-percent chance (100-year) peak streamflows were determined for the Fremont River and Bull Creek study reaches by using annual peak series data from streamflow-gaging stations and regional peak-flow regression equations. Surveyed topographic data for the study reaches were processed for use in the one-dimensional hydraulic models. The 1-percent chance (100-year) peak streamflows were simulated with hydraulic models to obtain predicted water-surface elevations. Water-surface elevations were then used to map the predicted inundation on a recent aerial photograph. The 1-percent chance (100-year) flood plain for the Fremont River in Hanksville, Utah, included some agricultural lands and did not encroach upon the town. The 1-percent chance (100-year) flood plain on the west side of Bull Creek was found to include a large portion of the town of Hanksville, Utah.

Utah↗

Effects of urbanization, construction activity, management practices, and impoundments on suspended-sediment transport in Johnson County, northeast Kansas, February 2006 through November 2008

The U.S. Geological Survey, in cooperation with the Johnson County, Kansas, Stormwater Management Program, investigated the effects of urbanization, construction activity, management practices, and impoundments on suspended-sediment transport in Johnson County from February 2006 through November 2008. Streamgages and continuous turbidity sensors were operated at 15 sites within the urbanizing 57-square-mile Mill Creek Basin, and 4 sites downstream from the other largest basins (49 to 66 square miles) in Johnson County. The largest sediment yields in Johnson County were observed downstream from basins with increased construction activity. Sediment yields attributed to the largest (68 acre) active construction site in the study area were 9,300 tons per square mile in 2007 and 12,200 tons per square mile in 2008; 5 to 55 times larger than yields observed at other sampling sites. However, given erodible soils and steep slopes at this site, sediment yields were relatively small compared to the range in historic values from construction sites without erosion and sediment controls in the United States (2,300 to 140,000 tons per square mile). Downstream from this construction site, a sediment forebay and wetland were constructed in series upstream from Shawnee Mission Lake, a 120-acre reservoir within Shawnee Mission Park. Although the original intent of the sediment forebay and constructed wetland were unrelated to upstream construction, they were nonetheless evaluated in 2008 to characterize sediment removal before stream entry into the lake. The sediment forebay was estimated to reduce 33 percent of sediment transported to the lake, whereas the wetland did not appear to decrease downstream sediment transport. Comparisons of time-series data and relations between turbidity and sediment concentration indicate that larger silt-sized particles were deposited within the sediment forebay, whereas smaller silt and clay-sized sediments were transported through the wetland and into the lake. Data collected at sites up and downstream from the constructed wetland indicated that hydraulic retention alone did not substantially reduce sediment loading to Shawnee Mission Lake. Mean-daily turbidity values at sampling sites downstream from basins with increased construction activity were compared to U.S. Environmental Protection Agency turbidity criteria designed to reduce discharge of pollutants from construction sites. The U.S. Environmental Protection Agency numeric turbidity criteria specifies that effluent from construction sites greater than 20 acres not exceed a mean-daily turbidity value of 280 nephelometric turbidity units beginning in 2011; this criteria will apply to sites greater than 10 acres beginning in 2014. Although numeric criteria would not have been applicable to data from sampling sites in Johnson County because they were not directly downstream from construction sites and because individual states still have to determine additional details as to how this criteria will be enforced, comparisons were made to characterize the potential of construction site effluent in Johnson County to exceed U.S. Environmental Protection Agency Criteria, even under extensive erosion and sediment controls. Numeric criteria were exceeded at sampling sites downstream from basins with increased construction activity for multiple days during the study period, potentially indicating the need for additional erosion and sediment controls and (or) treatment to bring discharges from construction sites into compliance with future numeric turbidity criteria. Among sampling sites in the Mill Creek Basin, sediment yields from the urbanizing Clear Creek Basin were approximately 2 to 3 times those from older, more stable urban or rural basins. Sediments eroded from construction sites adjacent to or surrounding streams appear to be more readily transported downstream, whereas sediments eroded from construction sites in headwater areas are more likely to

Scientific Investigations Report↗

Development of a precipitation-runoff model to simulate unregulated streamflow in the South Fork Flathead River Basin, Montana

This report documents the development of a precipitation-runoff model for the South Fork Flathead River Basin, Mont. The Precipitation-Runoff Modeling System model, developed in cooperation with the Bureau of Reclamation, can be used to simulate daily mean unregulated streamflow upstream and downstream from Hungry Horse Reservoir for water-resources planning. Two input files are required to run the model. The time-series data file contains daily precipitation data and daily minimum and maximum air-temperature data from climate stations in and near the South Fork Flathead River Basin. The parameter file contains values of parameters that describe the basin topography, the flow network, the distribution of the precipitation and temperature data, and the hydrologic characteristics of the basin soils and vegetation. A primary-parameter file was created for simulating streamflow during the study period (water years 1967-2005). The model was calibrated for water years 1991-2005 using the primary-parameter file. This calibration was further refined using snow-covered area data for water years 2001-05. The model then was tested for water years 1967-90. Calibration targets included mean monthly and daily mean unregulated streamflow upstream from Hungry Horse Reservoir, mean monthly unregulated streamflow downstream from Hungry Horse Reservoir, basin mean monthly solar radiation and potential evapotranspiration, and daily snapshots of basin snow-covered area. Simulated streamflow generally was in better agreement with observed streamflow at the upstream gage than at the downstream gage. Upstream from the reservoir, simulated mean annual streamflow was within 0.0 percent of observed mean annual streamflow for the calibration period and was about 2 percent higher than observed mean annual streamflow for the test period. Simulated mean April-July streamflow upstream from the reservoir was about 1 percent lower than observed streamflow for the calibration period and about 4 percent higher than observed for the test period. Downstream from the reservoir, simulated mean annual streamflow was 17 percent lower than observed streamflow for the calibration period and 12 percent lower than observed streamflow for the test period. Simulated mean April-July streamflow downstream from the reservoir was 13 percent lower than observed streamflow for the calibration period and 6 percent lower than observed streamflow for the test period. Calibrating to solar radiation, potential evapotranspiration, and snow-covered area improved the model representation of evapotranspiration, snow accumulation, and snowmelt processes. Simulated basin mean monthly solar radiation values for both the calibration and test periods were within 9 percent of observed values except during the month of December (28 percent different). Simulated basin potential evapotranspiration values for both the calibration and test periods were within 10 percent of observed values except during the months of January (100 percent different) and February (13 percent different). The larger percent errors in simulated potential evaporation occurred in the winter months when observed potential evapotranspiration values were very small; in January the observed value was 0.000 inches and in February the observed value was 0.009 inches. Simulated start of melting of the snowpack occurred at about the same time as observed start of melting. The simulated snowpack accumulated to 90-100 percent snow-covered area 1 to 3 months earlier than observed snowpack. This overestimated snowpack during the winter corresponded to underestimated streamflow during the same period. In addition to the primary-parameter file, four other parameter files were created: for a "recent" period (1991-2005), a historical period (1967-90), a "wet" period (1989-97), and a "dry" period (1998-2005). For each data file of projected precipitation and air temperature, a single parameter file can be used to simulate a s

Montana;Idaho↗

A precipitation-runoff model for simulating natural streamflow conditions in the Smith River watershed, Montana, water years 1996-2008

This report documents the construction of a precipitation-runoff model for simulating natural streamflow in the Smith River watershed, Montana. This Precipitation-Runoff Modeling System model, constructed in cooperation with the Meagher County Conservation District, can be used to examine the general hydrologic framework of the Smith River watershed, including quantification of precipitation, evapotranspiration, and streamflow; partitioning of streamflow between surface runoff and subsurface flow; and quantifying contributions to streamflow from several parts of the watershed. The model was constructed by using spatial datasets describing watershed topography, the streams, and the hydrologic characteristics of the basin soils and vegetation. Time-series data (daily total precipitation, and daily minimum and maximum temperature) were input to the model to simulate daily streamflow. The model was calibrated for water years 2002–2007 and evaluated for water years 1996–2001. Though water year 2008 was included in the study period to evaluate water-budget components, calibration and evaluation data were unavailable for that year. During the calibration and evaluation periods, simulated-natural flow values were compared to reconstructed-natural streamflow data. These reconstructed-natural streamflow data were calculated by adding Bureau of Reclamation’s depletions data to the observed streamflows. Reconstructed-natural streamflows represent estimates of streamflows for water years 1996–2007 assuming there was no agricultural water-resources development in the watershed. Additional calibration targets were basin mean monthly solar radiation and potential evapotranspiration. The model estimated the hydrologic processes in the Smith River watershed during the calibration and evaluation periods. Simulated-natural mean annual and mean monthly flows generally were the same or higher than the reconstructed-natural streamflow values during the calibration period, whereas they were lower during the evaluation period. The shape of the annual hydrographs for the simulated-natural daily streamflow values matched the shape of the hydrographs for the reconstructed-natural values for most of the calibration period, but daily streamflow values were underestimated during the evaluation period for water years 1996–1998. The model enabled a detailed evaluation of the components of the water budget within the Smith River watershed during the water year 1996–2008 study period. During this study period, simulated mean annual precipitation across the Smith River watershed was 16 inches, out of which 14 inches evaporated or transpired and 2 inches left the basin as streamflow. Per the precipitation-runoff model simulations, during most of the year, surface runoff rarely (less than 2 percent of the time during water years 2002–2008) makes up more than 10 percent of the total streamflow. Subsurface flow (the combination of interflow and groundwater flow) makes up most of the total streamflow (99 or more percent of total streamflow for 71 percent of the time during water years 2002–2008).

Montana↗

Trends in hydraulic fracturing distributions and treatment fluids, additives, proppants, and water volumes applied to wells drilled in the United States from 1947 through 2010: data analysis and comparison to the literature

Hydraulic fracturing is presently the primary stimulation technique for oil and gas production in low-permeability, unconventional reservoirs. Comprehensive, published, and publicly available information regarding the extent, location, and character of hydraulic fracturing in the United States is scarce. This national spatial and temporal analysis of data on nearly 1 million hydraulically fractured wells and 1.8 million fracturing treatment records from 1947 through 2010 (aggregated in Data Series 868) is used to identify hydraulic fracturing trends in drilling methods and use of proppants, treatment fluids, additives, and water in the United States. These trends are compared to the literature in an effort to establish a common understanding of the differences in drilling methods, treatment fluids, and chemical additives and of how the newer technology has affected the water use volumes and areal distribution of hydraulic fracturing. Historically, Texas has had the highest number of records of hydraulic fracturing treatments and associated wells in the United States documented in the datasets described herein. Water-intensive horizontal/directional drilling has also increased from 6 percent of new hydraulically fractured wells drilled in the United States in 2000 to 42 percent of new wells drilled in 2010. Increases in horizontal drilling also coincided with the emergence of water-based “slick water” fracturing fluids. As such, the most current hydraulic fracturing materials and methods are notably different from those used in previous decades and have contributed to the development of previously inaccessible unconventional oil and gas production target areas, namely in shale and tight-sand reservoirs. Publicly available derivative datasets and locations developed from these analyses are described.

Scientific Investigations Report↗

Selenium in the upper Blackfoot River watershed, southeastern Idaho, 2001-12

The upper Blackfoot River in southeastern Idaho receives runoff from 12 large phosphate mines. Waste shales that are removed to access the phosphate ore are highly enriched with selenium, resulting in elevated selenium in runoff from the mine waste dumps. In 2001, in cooperation with the Bureau of Land Management, the U.S. Geological Survey (USGS) began monitoring streamflow, selenium, and other water-quality parameters at a single location near the outlet of the upper Blackfoot River to the Blackfoot Reservoir. Water samples primarily were collected by a flow triggered, automated pump sampler, supplemented by manual point and equal-width integrated manual samples. The approach to monitoring concentrations and streamflow over time at a fixed location is ideal for evaluating temporal trends, but provides no information about the relative source contributions from the mine waste dumps draining into various tributaries. In 2001, the Idaho Department of Environmental Quality (IDEQ) began an annual, mid-May, synoptic survey of selenium concentrations and streamflow at 21 locations along the main stem Blackfoot River and its tributaries. Individually, neither the intensive USGS sampling at the outlet nor the IDEQ annual synoptic sampling provides a comprehensive view of selenium runoff in the Blackfoot River watershed. Together, the efforts are complementary; therefore, in this report, results are presented from both sampling efforts. The USGS collected time-series data from 2001 to 2012 at a fixed location, the Blackfoot River near the outlet of the reservoir, near Henry, Idaho (USGS streamgage 13063000). Dissolved selenium concentrations from 450 filtered samples collected at this site ranged from 0.5 to 11.4 micrograms per liter (μg/L). The State of Idaho chronic aquatic life criterion concentration of 5 μg/L was exceeded in 31 percent of the samples, with most exceedances occurring during May of each year. No exceedances of the selenium criterion were recorded in months other than April, May, or June. Concentrations of selenium in unfiltered and filtered samples were similar, and concentrations from samples collected by depth and width integrated methods were similar to those collected by grab (point) samples, indicating that the grab samples adequately represent selenium concentrations across the entire river cross section. In speciation analyses made during 2003 and 2004, the median percentage of total selenium as selenate was 81 percent, ranging from 17 to 98 percent, and the median percentage of total selenium as selenite was 19 percent, ranging from 2 to 83 percent of the total selenium. During the period of study, selenium concentrations had an upward trend during the lowflow season of August–October. Time trends were not obvious during other seasons. Selenium daily loads varied by more than a factor of 900 during the study period and ranged from 0.03 kilograms per day (kg/d) to more than 24 kg/d. Annual maximum daily loads of selenium varied over nearly a factor of 12, ranging from about 2 to 24 kg/d. For the annual spring synoptic samples collected by the IDEQ along the main stem Blackfoot River and major tributaries, selenium concentrations ranged from less than 2 to 870 μg/L in 176 samples. In most years, the synoptic sampling showed that the majority of the selenium loads passing the USGS streamgage at the outlet of the watershed could be attributed to a single tributary, East Mill Creek, which enters the Blackfoot River through Spring Creek. Selenium loads decreased by about half from East Mill Creek before reaching the Blackfoot River, suggesting that much selenium is at least temporarily removed from the water column by uptake by aquatic vegetation or by losses to sediment. Similar decreases in selenium loads occurred through the main stem Blackfoot River before reaching the outlet in low flow years, but not in high flow years.

Idaho↗

Simulation of climate change effects on streamflow, groundwater, and stream temperature using GSFLOW and SNTEMP in the Black Earth Creek Watershed, Wisconsin

A groundwater/surface-water model was constructed and calibrated for the Black Earth Creek watershed in south-central Wisconsin. The model was then run to simulate scenarios representing common societal concerns in the basin, focusing on maintaining a cold-water resource in an urbanizing fringe near its upper stream reaches and minimizing downstream flooding. Although groundwater and surface water are considered a single resource, many hydrologic models simplistically simulate feedback loops between the groundwater system and other hydrologic processes. These feedbacks include timing and rates of evapotranspiration, surface runoff, soil-zone flow, and interactions with the groundwater system; however, computer models can now routinely and iteratively couple the surface-water and groundwater systems—albeit with longer model run times. In this study, preliminary calibrations of uncoupled transient surface-water and steady-state groundwater models were used to form the starting point for final calibration of one transient computer simulation that iteratively couples groundwater and surface water. The computer code GSFLOW (Groundwater/Surface-water FLOW) was used to simulate the coupled hydrologic system; a surface-water model represented hydrologic processes in the atmosphere, at land surface, and within the soil zone, and a groundwater-flow model represented the unsaturated zone, saturated zone, and streams. The coupled GSFLOW model was run on a daily time step during water years 1985–2007. Early simulation times (1985–2000) were used for spin-up to make the simulation results less sensitive to initial conditions specified; the spin-up period was not included in the model calibration. Model calibration used observed heads, streamflows, solar radiation, and snowpack measurements from 2000 to 2007 for history matching. Calibration was performed by using the PEST parameter estimation software suite. Simulated streamflows from the calibrated GSFLOW model and other basin characteristics were used as input to the one-dimensional SNTEMP (Stream-Network TEMPerature) model. SNTEMP was used to simulate daily stream temperature in selected stream reaches in the watershed. The temperature model was calibrated to high-resolution stream temperature time-series data measured in 2005. The calibrated GSFLOW and SNTEMP models were then used to simulate effects of potential climate change for the years 2010 through 2100. An ensemble of climate models and emission scenarios was evaluated. Downscaled climate drivers for the simulation period showed increases in maximum and minimum air temperature. Scenarios of future precipitation, however, did not show a monotonic trend like temperature. Uncertainty in the climate drivers increased with time for both temperature and precipitation. Forecasts of potential climate change scenarios showed growing season length increasing by weeks, and both potential and actual evapotranspiration rates increasing appreciably, in response to increasing air temperature. Simulated actual evapotranspiration rates increased less than simulated potential evapotranspiration rates as a result of water limitation in the root zone during the summer high-evapotranspiration period. The hydrologic-system response to climate change was characterized by a reduction in the importance of the snowmelt pulse and an increase in the importance of fall and winter groundwater recharge. The less dynamic hydrologic regime is likely to result in drier soil conditions, with relatively less drying expected in groundwater-fed systems. Groundwater discharge in the current upper cold-water reaches of Black Earth Creek is expected to decrease; flooding in downstream reaches may appreciably increase. The magnitude of changes in forecasted flow and associated groundwater/surface-water interaction is dependent on the General Circulation Model and emission scenario chosen. Potential future changes in air temperature drivers were consistently upward regardless of General Circulation Model and emission scenario selected; thus, simulated stream temperatures are forecast to increase appreciably with future climate. However, the amount of temperature increase was variable. Such uncertainty is reflected in temperature model results, along with uncertainty in the groundwater/surface-water interaction itself. The estimated increase in annual average temperature ranged from approximately 3 to 6 degrees Celsius by 2100 in the upper reaches of Black Earth Creek and 2 to 4 degrees Celsius in reaches farther downstream. As with all forecasts that rely on projections of an unknowable future, the results are best considered to approximate potential outcomes of climate change given the underlying uncertainty.

Wisconsin↗

Simulation of hydrodynamics and water temperature in a 21-mile reach of the upper Illinois River, Illinois, 2020–22

This report describes the development of a CE-QUAL-W2 river hydrodynamics and temperature model of a 21-mile reach of the Illinois River including a 3-mile reach of a major tributary, the Fox River. Model outputs consist of streamflow, water velocity, water-surface elevation, and water-temperature time series that can be used to simulate summer conditions in years with and without extensive development of harmful algal blooms (HABs). These analyses may provide a better understanding of some complex factors contributing to HAB development along the Illinois River. Such an understanding may provide more accurate HAB timing and location predictions and may help determine potential mitigating activities to prevent or limit the size and duration of HABs. Using the observed and simulated hydrodynamic conditions in the Illinois River study reach, it was possible to compare and contrast streamflow, velocity, and temperature conditions in years with varying HAB distributions. Occurrences of extensive HABs were documented in the study reach in June 2020 and June 2021, but only a small HAB restricted to the Marseilles Lock and Dam pool occurred in the summer of 2022. The objective then was to find similarities in site conditions between 2020 and 2021 that may contrast with the conditions in 2022. Among the 3 years included in the study, the variability in simulated water temperature exceeded variability in observed streamflow and simulated velocities. The longest period of water temperatures greater than 27 degrees Celsius (°C) in the selected locations in June of the three analysis years was in the second half of June 2022, yet no study-area wide HAB was documented in 2022. Simulations indicated that after warm water temperatures were established in the reach in June 2022, a cooling period broke up the warming period. This period of cooling was greater in magnitude and duration downstream from the location of a localized HAB perhaps limiting the spread of the bloom. Residence times differed substantially in segments representing different channel features; values ranged from 0.28 to 17.3 (days per 500 meters of channel) between the main stem and backwater areas, respectively. Variation in average June residence times was also greater among different channel features than among different years in the study period. The HABs in 2020 and 2021 at Starved Rock Dam were documented when water temperatures were about 26 °C. River backwater areas at some locations did attain these temperatures 2 to 3 days before the conditions in the main stem. Residence times in the backwater areas, however, generally exceeded 9 days, thus limiting the exchange of water carrying algal biomass into the main channel. Hydrodynamic model calibration involved adjusting model parameters until observed and simulated daily water-surface elevations, daily streamflows, discrete velocities, and channel areas were similar. Temperature calibration was done with near-surface continuous time-series data and discrete vertical profile temperatures. Observed and simulated water temperatures generally were within 1 °C at all monitoring locations.

Illinois↗

Analysis of water use associated with hydraulic fracturing and determination of baseline water quality in watersheds within the shale play of eastern Ohio, 2021–23

The U.S. Geological Survey, in cooperation with the Ohio Department of Natural Resources, performed a two-part study to (1) assess water use and temporal trends and changes in streamflow, and to (2) characterize 2021–23 baseline water quality, as they relate to oil and gas extraction activities in selected eastern Ohio watersheds. Between calendar years 2010 and 2019, hydraulic fracturing water withdrawals totaling about 27,168 million gallons were reported at 643 locations in Ohio. In 2021, wells developed with hydraulic fracturing were the source of most of the oil and gas produced in Ohio. Daily streamflow time-series data from seven study gages and two reference gages were used to assess temporal trends and changes in streamflow. The study gages were in basins with reported water withdrawals for hydraulic fracturing. The reference gages, which have long periods of record and were subject to minimal streamflow regulation, were in nearby basins with no hydraulic fracturing water withdrawals. Trend slopes for the period of record annual minimum and median daily streamflows and for annual daily streamflow nonexceedance probabilities less than 0.9 were all uniformly positive at the study and reference gages. This trend indicates a consistently increasing pattern over the periods of record, except for high flows. In addition, analyses of annual streamflow statistics showed no general indication that low flows or extreme low flows at the reference or study gages have lowered, become more frequent, or lengthened in duration since 2010, when records for hydraulic fracturing water withdrawals began in Ohio. In fact, in almost all cases, the opposite was indicated. Nonexceedance percentiles of daily streamflows were compared between the full and pre-2012 periods of record for reference and study gages. The streamflows associated with nonexceedance percentiles in the lower quartile of daily streamflows determined for the full period of record were larger than or equal to those determined for the pre-2012 period of record for all study and reference gages. This indicates that low flows did not decrease during the post-2011 period of record when water was withdrawn for hydraulic fracturing. Water-quality data were collected eight times at each of eight sampling sites (six of which were colocated with the study gages). Sampling was done during a variety of flow conditions to assess baseline water quality. In 2021, the 8 sampling sites had drainage basins that were wholly or partially within 7 of the 10 most active counties in Ohio for oil and gas development. As part of the record of baseline conditions, water-quality data were used to assess (1) water types based on major-ion chemistry; (2) sources of salinity to streams; (3) exceedances of aquatic life use criteria; and (4) the correlations between water chemistry and drainage-basin characteristics, such as density of oil and gas wells, density of wastewater treatment plants, or the percentage of different types of land cover (agriculture, developed, forest). Seven of the water-quality sampling sites were designated as coal-mine impacted based on criteria developed for assessing mine-drainage impacts in Ohio. Mine drainage from historical coal mining in the region likely affected the quality of these streams and complicated the use of some constituents typically used as indicators of oil and gas influence. Based on major-ion chemistry, three main types of water were in the study area―sulfate (three sites), calcium-bicarbonate (one site), and mixed bicarbonate-chloride (four sites) type waters. One site had samples with a higher proportion of sodium and chloride ions than other stream samples, indicating potential contamination with oil-field brine or road salt. Binary mixing curves revealed that 11 samples from 4 of the sampling sites likely contained a component of brine. The results of the baseline assessment of surface-water quality in the study area showed no exceedances of Ohio Environmental Protection Agency aquatic life use criteria. Spearman’s rank correlation coefficients indicated no significant positive correlations with the density of vertical or horizontal oil and gas wells.

Ohio↗