USGS Science⌕ Search

SEARCH · USGS Science

Results for “Cooperative Investigations Report”

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

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

At least 1,333 records · Page 74Linked to original sources

Hydrogeologic and geochemical characterization and evaluation of two arroyos for managed aquifer recharge by surface infiltration in the Pojoaque River Basin, Santa Fe County, New Mexico, 2014–15

In order to provide long-term storage of diverted surface water from the Rio Grande as part of the Aamodt water rights settlement, managed aquifer recharge by surface infiltration in Pojoaque River Basin arroyos was proposed as an option. The initial hydrogeologic and geochemical characterization of two arroyos located within the Pojoaque River Basin was performed in 2014 and 2015 in cooperation with the Bureau of Reclamation to evaluate the potential suitability of these two arroyos as sites for managed aquifer recharge through surface infiltration. The selected reaches were high-gradient (average 3.0–3.5 percent) braided channels filled with unconsolidated sand and gravel-sized deposits that were generally 30–50 feet thick. Saturation was not observed in the unconsolidated channel sands in four subsurface borings but was found at 7–60 feet below the contact between the unconsolidated channel sands and the bedrock. The poorly to well-cemented alluvial deposits that make up the bedrock underlying the unconsolidated channel material is the Tesuque Formation. The individual beds of the Tesuque Formation are reported to be highly heterogeneous and anisotropic, and the bedrock at the site was observed to have variable moisture and large changes in lithology. Surface electrical-resistivity geophysical survey methods showed a sharp contrast between the electrically resistive unconsolidated channel sands and the highly conductive bedrock; however, because of the high conductivity, the resistivity methods were not able to image the water table or preferential flow paths (if they existed) in the bedrock. Infiltration rates measured by double-ring and bulk infiltration tests on a variety of channel morphologies in the study reaches were extremely large (9.7–94.5 feet per day), indicating that the channels could potentially accommodate as much as 6.6 cubic feet per second of applied water without generating surface runoff out of the reach; however, the small volume available for storage in the unconsolidated channel sands (about 410 acre-feet in the east arroyo and about 190 acre-feet in the west arroyo) and the potential for the infiltrating water to preferentially flow over the bedrock contact and out of the reach present a challenge for storing water. Although a detailed assessment of the infiltration rate of the Tesuque Formation is beyond the scope of this investigation, one double-ring infiltrometer test was conducted on an outcrop, resulting in an estimated infiltration rate of about 4 feet per day. The shallow groundwater observed in this investigation was determined to be recharged locally on the basis of groundwater elevations and geochemical and isotopic signatures. The channel sands and shallow bedrock were observed to be weathered, indicating contact with oxic groundwater following deposition. This observation was supported by whole-rock elemental analysis and mineralogy of several core samples. The downward groundwater gradient between the shallow wells and those wells screened at greater depths suggests that the shallow groundwater is recharged by local precipitation and has the potential to migrate to the deeper aquifer units. The two age-dating tracers measured in this investigation, however, demonstrate that the shallow groundwater flow paths are very slow and that the deeper flow paths are likely part of a larger regional system. The composition of the shallow, native groundwater suggests that storing water diverted from the Rio Grande is not likely to leach constituents of concern that would cause the stored water to exceed health-based U.S. Environmental Protection Agency Maximum Contaminant Levels.

New Mexico↗

Updates to the Madison Lake (Minnesota) CE–QUAL–W2 water-quality model for assessing algal community dynamics

A previously developed CE–QUAL–W2 model for Madison Lake, Minnesota, simulated the algal community dynamics, water quality, and fish habitat suitability of Madison Lake under recent (2014) meteorological conditions. Additionally, this previously developed model simulated the complex interplay between external nutrient loading, internal nutrient loading from sediment release of phosphorus, and the organic matter decomposition of the algal biomass. However, the partitioning of Cyanophyta within the modeling framework was simplified to one group and did not account for how different Cyanophyta populations are affected by light conditions, use of nitrogen, temperature growth ranges, and differences in settling rates. Properly capturing Cyanophyta dynamics is important given the potential risks posed by potential large algal blooms. For example, when Cyanophyta form large blooms, recreational activities can become restricted in certain areas because of thick algal scums or algal mats, in addition to the possible production of a class of toxins, known as cyanotoxins, capable of threatening human health, domestic animals, and wildlife. Therefore, we updated the model to partition the Cyanophyta into a group that fixed nitrogen and a second, more buoyant Cyanophyta group that did not independently fix nitrogen. The U.S. Geological Survey, in cooperation with the St. Croix Watershed Research Station (Science Museum of Minnesota) with support from the Environmental and Natural Resources Trust Fund of Minnesota (Legislative-Citizen Commission on Minnesota Resources), updated the Madison Lake CE–QUAL–W2 model to address the shortcomings of simulating Cyanophyta in the previously developed model and better characterize Cyanophyta into two groups. In addition to updating the Cyanophyta group differentiation, the part of the model that handles the simulation of algal community dynamics was updated while preserving model predictive capabilities for nutrients, water temperature, and dissolved oxygen. The calibration and validation of the model was done under recent meteorological conditions with large and persistent Cyanophyta blooms (2014 and 2016). Overall, the model simulations predicted the persistently large total phosphorus concentrations in the hypolimnion of Madison Lake and key differences in nutrient concentrations between 2014 and 2016. The Cyanophyta bloom persistence throughout the summer was also simulated by the model in 2014 and 2016, a critical goal of the model update. Finally, monthly total phosphorus budgets were calculated for the updated Madison Lake model for 2014 and 2016.

Minnesota↗

Paleomagnetic correlation of surface and subsurface basalt flows in the central and southwestern part of the Idaho National Laboratory, Idaho

The U.S. Geological Survey, in cooperation with the U.S. Department of Energy, used paleomagnetic data from 22 coreholes to construct 3 fence diagrams of subsurface basalt flows in the southern part of the Idaho National Laboratory. These diagrams provide comprehensive descriptions of the horizontal and vertical distribution of basalt flows and sediment layers beneath the surface, aiding geological studies and contributing valuable data to numerical models of groundwater flow and contaminant transport. The correlations established though these diagrams include spatial correlations between basalt flows found in multiple coreholes. Correlations were identified by matching average paleomagnetic inclinations and confirming or denying these correlations using petrology, geochemistry and radiometric ages. The fence diagrams aid in identifying potential locations of subsurface vents, volcanic vents that have been buried by more recent volcanic activity, associated to subsurface basalt flows. By tracing the subsurface flows and analyzing where the greatest thickness occurs, the locations of buried vents can be inferred. Some subsurface flows exhibit correlations across several coreholes and may indicate yet unidentified surface or buried vents, thereby enhancing our understanding of the volcanic history and subsurface geology of the region.

Idaho↗

Estimates of ground-water recharge, base flow, and stream reach gains and losses in the Willamette River basin, Oregon

Precipitation-runoff models, base-flow-separation techniques, and stream gain-loss measurements were used to study recharge and ground-water surface-water interaction as part of a study of the ground-water resources of the Willamette River Basin. The study was a cooperative effort between the U.S. Geological Survey and the State of Oregon Water Resources Department. Precipitation-runoff models were used to estimate the water budget of 216 subbasins in the Willamette River Basin. The models were also used to compute long-term average recharge and base flow. Recharge and base-flow estimates will be used as input to a regional ground-water flow model, within the same study. Recharge and base-flow estimates were made using daily streamflow records. Recharge estimates were made at 16 streamflow-gaging-station locations and were compared to recharge estimates from the precipitation-runoff models. Base-flow separation methods were used to identify the base-flow component of streamflow at 52 currently operated and discontinued streamflow-gaging-station locations. Stream gain-loss measurements were made on the Middle Fork Willamette, Willamette, South Yamhill, Pudding, and South Santiam Rivers, and were used to identify and quantify gaining and losing stream reaches both spatially and temporally. These measurements provide further understanding of ground-water/surface-water interactions.

Water-Resources Investigations Report↗

Evaluation of ground-water and boron sources by use of boron stable-isotope ratios, tritium, and selected water-chemistry constituents near Beverly Shores, northwestern Indiana, 2004

Concentrations of boron greater than the U.S. Environmental Protection Agency (USEPA) 900 μg/L removal action level (RAL) standard were detected in water sampled by the USEPA in 2004 from three domestic wells near Beverly Shores, Indiana. The RAL regulates only human-affected concentrations of a constituent. A lack of well logs and screened depth information precluded identification of whether water from sampled wells, and their boron sources, were from human-affected or natural sources in the surficial aquifer, or associated with a previously defined natural, confined aquifer source of boron from the subtill or basal sand aquifers. A geochemically-based classification of the source of boron in ground water could potentially determine the similarity of boron to known sources or mixtures between known sources, or classify whether the relative age of the ground water predated the potential sources of contamination. The U.S. Geological Survey (USGS), in cooperation with the USEPA, investigated the use of a geochemical method that applied boron stable isotopes, and concentrations of boron, tritium, and other constituents to distinguish between natural and human-affected sources of boron in ground water and thereby determine if the RAL was applicable to the situation. Boron stable-isotope ratios and concentrations of boron in 17 ground-water samples and tritium concentrations in 9 ground-water samples collected in 2004 were used to identify geochemical differences between potential sources of boron in ground water near Beverly Shores, Indiana. Boron and δ 11 B analyses for this investigation were made on unacidified samples to assure consistency of the result with unacidified analyses of δ 11 B values from other investigations. Potential sources of boron included surficial-aquifer water affected by coal-combustion products (CCP) or domestic-wastewater, upward discharge of ground water from confined aquifers, and unaffected water from the surficial aquifer that was distant from human-affected boron sources. Boron concentrations in potential ground-water sources of boron were largest (15,700 to 24,400 μg/L) in samples of CCP-affected surficial aquifer water from four wells at a CCP landfill and smallest (27 to 63 μg/L) in three wells in the surficial aquifer that were distant from human-affected boron sources. Boron concentrations in water from the basal sand aquifer ranged from 656 μg/L to 1,800 μg/L. Boron concentrations in water from three domestic-wastewater-affected surficial aquifer wells ranged from 84 to 387 μg/L. Among the representative ground-water samples, boron concentrations from all four samples of CCP-affected surficial aquifer water and four of five samples of water from the basal sand aquifer had concentrations greater than the RAL. A comparison of boron concentrations in acid-preserved and unacidified samples indicated that boron concentrations reported for this investigation may be from about 11 to 16 percent less than would be reported in a standard analysis of an acidified sample. The stable isotope boron-11 was most enriched in comparison to boron-10 in ground water from a confined aquifer, the basal sand aquifer (δ 11 B, 24.6 to 34.0 per mil, five samples); it was most depleted in CCP-affected water from the surficial aquifer (δ 11 B, 0.1 to 6.6 per mil, four samples). Domestic-wastewater-affected water from the surficial aquifer (δ 11 B, 8.7 to 11.7 per mil, four samples) was enriched in boron-11, in comparison to individual samples of a borax detergent additive and a detergent with perborate bleach; it was intermediate in composition between basal sand aquifer water and CCP-affected water from the surficial aquifer. The similarity between a ground-water sample from the surficial aquifer and a hypothetical mixture of unaffected surficial aquifer and basal sand aquifer waters indicates the potential for long-term upward discharge of ground water into the surficial aquifer from one or more confined aquifers. Estimated δ 11 B values for acidified samples were depleted by 1.9 to 2.8 per mil in comparison to unacidified samples from the four wells sampled; those differences were small in comparison to the differences between δ 11 B values of representative sources of boron in ground water. Tritium concentrations ranged from 7.0 to 10.3 tritium units in six samples from the surficial aquifer and were less than 0.8 tritium units in three samples from the basal sand aquifer. Water from wells in the surficial aquifer represents predominantly modern, post-1972 recharge and sources of boron and other constituents. Water from the basal sand aquifer is associated with pre-1952 recharge from sources not affected by local boron inputs. Ground water from six wells (five domestic wells and one public-supply well) where the ground-water source was unknown had boron concentrations, boron isotope ratios, and tritium concentrations similar to water from the basal sand aquifer. Boron concentrations greater than the RAL were found in water from four of these six wells. The boron isotope and tritium data from these four wells indicate a natural source of boron in ground water; therefore, the RAL does not apply to boron concentrations in water from these wells. Water samples from two domestic wells where the ground-water source was unknown had boron concentrations less than the RAL and boron isotope ratios and tritium concentrations that were similar to domestic-wastewater-affected water from the surficial aquifer. The boron isotope ratio for a sample from one domestic well was similar to that of CCP-affected water from the surficial aquifer and detergent compositions; the boron concentration of that sample was less than the RAL. The classifications of differences among representative sources of boron in ground water and water samples from wells where the ground-water source was unknown generally agreed with distinctions based on strontium-87/strontium-86 ratios and concentrations of strontium, chloride, nitrate, and ammonia. This application of boron concentrations, boron isotope ratios, and tritium concentrations to classify differences in relation to potential sources of boron in ground water was able to distinguish between boron from natural sources and from human-affected sources that are subject to regulation.

Indiana↗

Status of shallow-aquifer mapping in the Northern Front Range Area, Colorado

Mapping of shallow aquifers in the northern Front Range area of Colorado has been completed as part of the U.S. Geological Survey Front Range Infrastructure Resources Project. The aquifer mapping was undertaken as part of a comprehensive effort to better define the mineral, energy, cartographic, biological, and water resources that are critical to the support and development of the area's infrastructure, such as streets. highways, airports, and buildings. The aquifer mapping was undertaken in cooperation with the Colorado Division of Water Resources and the Colorado Water Conservation Board. The shallow aquifers have been mapped in a 2,450-square-mile area extending as an approximately 30-mile-wide band from north of Fort Collins to the Arapahoe-Douglas County line south of Denver (fig. I). The shallow aquifer mapping in the Denver metropolitan area was published in 1996 as Hydrologic Investigations Atlas HA-736 (Robson. 1996). Shallow aquifer mapping in the Greeley-Nunn area was published as HA-746A (Robson, Arnold, and Heiny, 2000a); mapping in the Fort Collins-Loveland area was published as HA-746B (Robson, Arnold, and Heiny, 2000b); mapping in the Fort Lupton-Gilcrest area was published as HA-746C (Robson, Heiny, and Arnold, 2000c); and mapping in the Boulder-Longmont area was published as HA-746D (Robson, Heiny, and Arnold, 2000d). Each of the five atlases contains five map sheets at 1:50,000 scale showing: 1. The thickness and extent of the unconsolidated sediments (loose gravel. sand. silt, and clay) that overlie the bedrock formations in the area (fig. 2). 2. The altitude and configuration of the bedrock surface. 3. The altitude of the water table and direction of ground-water movement. 4. The saturated thickness of the shallow aquifers. 5. The depth to the water table in the shallow aquifers.

Fact Sheet↗

Irrigated acreage and other land uses on the Snake River Plain, Idaho and eastern Oregon

Prompted by the need for a current, accurate, and repeatable delineation of irrigated acreage on the Snake River Plain, the U.S. Geological Survey entered into a cooperative agreement with the Idaho Department of Water Resources Image Analysis Facility and the U.S. Bureau of Reclamation to delineate 1980 land use form Landsat data. Irrigated acreage data were needed as input to groundwater flow models developed by the U.S. Geological Survey in a study of the regional aquifer system underlying the Snake River Plain. Single-date digital multispectral scanner data analyzed to delineate land-use classes. Source of irrigation water (surface water, ground water, and combined) was determined from county maps of 1975 water-related land use, data from previous investigations, and field checking. Surface-water diversions for irrigation on the Snake River Plain began in the 1840's. With the stimulus of Federal aid authorized by the Desert Land Act, Carey Act, and Reclamation Act, irrigated area increased rapidly in the early 1900's. By 1929, 2.2 million acres were irrigated. Ground water became and important source of irrigation water after World War II. In 1980, about 3.1 million acres of the Snake River Plain were irrigate: 2.0 million acres with surface water, 1.0 million with ground water, and 0.1 million with combined surface and ground water. About 5.2 million acres (half of the plain) are undeveloped rangeland, 1.0 million acres (one-tenth) are classified as barren. The remaining land is a mixture of dryland agriculture, water bodies, wetland, forests, and urban areas.

Idaho;Oregon↗

Under the Golden Gate Bridge — Views of the sea floor near the entrance to San Francisco Bay, California

San Francisco Bay in Northern California is one of the largest and most altered estuaries within the United States. The sea floor within the bay as well as at its entrance is constantly changing due to strong tidal currents, aggregate mining, dredge disposal, and the creation of new land using artificial fill. Understanding this dynamic sea floor is critical for addressing local environmental issues, which include defining pollution transport pathways, deciphering tectonics, and identifying benthic habitats. Mapping commercial interests such as safe ship navigation and dredge disposal is also significantly aided by such understanding. Over the past decade, the U.S. Geological Survey (USGS), the National Oceanic and Atmospheric Administration (NOAA), and California State University, Monterey Bay (CSUMB) in cooperation with the U.S. Army Corps of Engineers (USACOE) and the Center for Integrative Coastal Observation, Research and Education (CICORE) have partnered to map central San Francisco Bay and its entrance under the Golden Gate Bridge using multibeam echosounders. These sonar systems can continuously map to produce 100 percent coverage of the sea floor at meter-scale resolution and thus produce an unprecedented view of the floor of the bay. This poster shows views of the sea floor in west-central San Francisco Bay around Alcatraz and Angel Islands, underneath the Golden Gate Bridge, and through its entrance from the Pacific Ocean. The sea floor is portrayed as a shaded relief surface generated from the multibeam data color-coded for depth from light blues for the shallowest values to purples for the deepest. The land regions are portrayed by USGS digital orthophotographs (DOQs) overlaid on USGS digital elevation models (DEMs). The water depths have a 4x vertical exaggeration while the land areas have a 2x vertical exaggeration.

California↗

Hydrogeologic investigation, framework, and conceptual flow model of the Antlers aquifer, southeastern Oklahoma, 1980–2022

The 1973 Oklahoma Groundwater Law (Oklahoma Statute §82–1020.5) requires that the Oklahoma Water Resources Board conduct hydrologic investigations of the State’s groundwater basins to support a determination of the maximum annual yield for each groundwater basin. Every 20 years, the Oklahoma Water Resources Board is required to update the hydrologic investigation on which the maximum annual yield determinations were based. The maximum annual yield allocated per acre of land is used to set the equal-proportionate share pumping rate. The maximum annual yield of 5,913,600 acre-feet per year and equal-proportionate-share of 2.1 acre-feet per acre per year currently (2025) in place for the Antlers aquifer were issued by the Oklahoma Water Resources Board on February 14, 1995. Because more than 20 years have elapsed since the 1995 final order for the Antlers aquifer was issued, the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, completed an in-depth hydrologic study that included a hydrogeologic framework and conceptual groundwater-flow model for the 1980–2022 study period. The results of an analysis of land use, long-term climate patterns, streamflow and base-flow patterns, historical groundwater use, as well as groundwater-level fluctuations across the Antlers aquifer are described. In addition, groundwater quality was analyzed for total dissolved solids concentrations and major ions for the Antlers aquifer. An updated hydrogeologic framework was developed that included refining the aquifer boundary in Oklahoma, the creation of new potentiometric surface and saturated thickness of fresh groundwater maps, one multiple-well aquifer test, slug tests, and an analysis of lithologic logs across the aquifer. A conceptual groundwater flow model and water budget were developed by incorporating estimates of recharge from precipitation, saturated-zone evapotranspiration, streambed seepage, lateral groundwater flows, vertical leakage, and withdrawals from groundwater wells.

Oklahoma, Texas↗

Reconnaissance investigation of petroleum products in soil and ground water at Longmire, Mount Rainier National Park, Washington, 1990

The removal of an underground waste-oil storage tank in Mount Rainier National Park, at Longmire, Washington, led to the discovery that soil surrounding the tank was saturated with unidentified petroleum hydrocarbons. Subsequent investigations by the National Park Service indicated that a petroleum product smelling like diesel oil was present in the unsaturated zone as far as 120 feet from the tank site. A study was conducted by the U.S. Geological Survey in cooperation with the National Park Service to determine the extent to which the petroleum hydrocarbons have affected the unsaturated zone and ground water in the Longmire area. Measurements of water levels in wells and of water-surface elevations of the Nisqually River and a wetland west of Longmire indicate that ground water does not flow from the maintenance area to the river or to the wetland. Waste oil and diesel oil were detected in soil samples from the site closest to the waste-oil storage-tank site. Diesel oil was also detected in samples from a site about 200 feet northwest of the storage-tank site. Organic compounds of undetermined origin were detected in soil samples from all of the other sites. Waste oil was not conclusively detected in any of the ground-water samples. Diesel oil was detected in water samples from the well closest to the storage tank and from a well about 200 feet west of the storage-tank site. Ground-water samples from all of the other wells contained organic compounds of undetermined origin.

Water-Resources Investigations Report↗

Bedrock and surficial geologic map of the Satan Butte and Greasewood 7.5’ quadrangles, Navajo and Apache Counties, northern Arizona

The geologic map of the Satan Butte and Greasewood 7.5’ quadrangles is the result of a cooperative effort of the U.S. Geological Survey (USGS) and the Navajo Nation to provide regional geologic information for management and planning officials. This map provides geologic information useful for range management, plant and animal studies, flood control, water resource investigations, and natural hazards associated with sand-dune mobility. The map provides connectivity to the regional geologic framework of the Grand Canyon area of northern Arizona. The map area encompasses approximately 314 km 2 (123 mi 2 ) within Navajo and Apache Counties of northern Arizona and is bounded by lat 35°37'30" to 35°30' N., long 109°45' to 110° W. The quadrangles lie within the southern Colorado Plateau geologic province and within the northeastern portion of the Hopi Buttes (Tsézhin Bií). Large ephemeral drainages, Pueblo Colorado Wash and Steamboat Wash, originate north of the map area on the Defiance Plateau and Balakai Mesa respectively. Elevations range from 1,930 m (6,330 ft) at the top of Satan Butte to about 1,787 m (5,860 ft) at Pueblo Colorado Wash where it exits the southwest corner of the Greasewood quadrangle. The only settlement within the map area is Greasewood, Arizona, on the north side of Pueblo Colorado Wash. Navajo Highway 15 crosses both quadrangles and joins State Highway 264 northwest of Ganado. Unimproved dirt roads provide access to remote parts of the Navajo Reservation.

Arizona↗

Hydrology of Eagle Creek Basin and effects of groundwater pumping on streamflow, 1969-2009

Urban and resort development and drought conditions have placed increasing demands on the surface-water and groundwater resources of the Eagle Creek Basin, in southcentral New Mexico. The Village of Ruidoso, New Mexico, obtains 60-70 percent of its water from the Eagle Creek Basin. The village drilled four production wells on Forest Service land along North Fork Eagle Creek; three of the four wells were put into service in 1988 and remain in use. Local citizens have raised questions as to the effects of North Fork well pumping on flow in Eagle Creek. In response to these concerns, the U.S. Geological Survey, in cooperation with the Village of Ruidoso, conducted a hydrologic investigation from 2007 through 2009 of the potential effect of the North Fork well field on streamflow in North Fork Eagle Creek. Mean annual precipitation for the period of record (1942-2008) at the Ruidoso climate station is 22.21 inches per year with a range from 12.27 inches in 1970 to 34.81 inches in 1965. Base-flow analysis indicates that the 1970-80 mean annual discharge, direct runoff, and base flow were 2,260, 1,440, and 819 acre-ft/yr, respectively, and for 1989-2008 were 1,290, 871, and 417 acre-ft/yr, respectively. These results indicate that mean annual discharge, direct runoff, and base flow were less during the 1989-2008 period than during the 1970-80 period. Mean annual precipitation volume for the study area was estimated to be 12,200 acre-feet. Estimated annual evapotranspiration for the study area ranged from 8,730 to 8,890 acre-feet. Estimated annual basin yield for the study area was 3,390 acre-ft or about 28 percent of precipitation. On the basis of basin-yield computations, annual recharge was estimated to be 1,950 acre-ft, about 16 percent of precipitation. Using a chloride mass-balance method, groundwater recharge over the study area was estimated to average 490 acre-ft, about 4.0 percent of precipitation. Because the North Fork wells began pumping in 1988, 1969-80 represents the pre-groundwater-pumping period, and 1988-2009 represents the groundwater-pumping period. The 5-year moving average for precipitation at the Ruidoso climate station shows years of below-average precipitation during both time periods, but no days of zero flow were recorded for the 11-year period 1970-80 and no-flow days were recorded in 11 of 20 years for the 1988-2009 period. View report for unabridged abstract.

Scientific Investigations Report↗

A three-dimensional numerical model of predevelopment conditions in the Death Valley regional ground-water flow system, Nevada and California

In the early 1990's, two numerical models of the Death Valley regional ground-water flow system were developed by the U.S. Department of Energy. In general, the two models were based on the same basic hydrogeologic data set. In 1998, the U.S. Department of Energy requested that the U.S. Geological Survey develop and maintain a ground-water flow model of the Death Valley region in support of U.S. Department of Energy programs at the Nevada Test Site. The purpose of developing this 'second-generation' regional model was to enhance the knowledge an understanding of the ground-water flow system as new information and tools are developed. The U.S. Geological Survey also was encouraged by the U.S. Department of Energy to cooperate to the fullest extent with other Federal, State, and local entities in the region to take advantage of the benefits of their knowledge and expertise. The short-term objective of the Death Valley regional ground-water flow system project was to develop a steady-state representation of the predevelopment conditions of the ground-water flow system utilizing the two geologic interpretations used to develop the previous numerical models. The long-term objective of this project was to construct and calibrate a transient model that simulates the ground-water conditions of the study area over the historical record that utilizes a newly interpreted hydrogeologic conceptual model. This report describes the result of the predevelopment steady-state model construction and calibration. The Death Valley regional ground-water flow system is situated within the southern Great Basin, a subprovince of the Basin and Range physiographic province, bounded by latitudes 35 degrees north and 38 degrees 15 minutes north and by longitudes 115 and 118 degrees west. Hydrology in the region is a result of both the arid climatic conditions and the complex geology. Ground-water flow generally can be described as dominated by interbasinal flow and may be conceptualized as having two main components: a series of relatively shallow and localized flow paths that are superimposed on deeper regional flow paths. A significant component of the regional ground-water flow is through a thick Paleozoic carbonate rock sequence. Throughout the flow system, ground water flows through zones of high transmissivity that have resulted from regional faulting and fracturing. The conceptual model of the Death Valley regional ground-water flow system used for this study is adapted from the two previous ground-water modeling studies. The three-dimensional digital hydrogeologic framework model developed for the region also contains elements of both of the hydrogeologic framework models used in the previous investigations. As dictated by project scope, very little reinterpretation and refinement were made where these two framework models disagree; therefore, limitations in the hydrogeologic representation of the flow system exist. Despite limitations, the framework model provides the best representation to date of the hydrogeologic units and structures that control regional ground-water flow and serves as an important information source used to construct and calibrate the predevelopment, steady-state flow model. In addition to the hydrogeologic framework, a complex array of mechanisms accounts for flow into, through, and out of the regional ground-water flow system. Natural discharges from the regional ground-water flow system occur by evapotranspiration, springs, and subsurface outflow. In this study, evapotranspiration rates were adapted from a related investigation that developed maps of evapotranspiration areas and computed rates from micrometeorological data collected within the local area over a multiyear period. In some cases, historical spring flow records were used to derive ground-water discharge rates for isolated regional springs. For this investigation, a process-based, numerical model was developed to estimat

Water-Resources Investigations Report↗

Characterization of the hydrogeologic framework, groundwater-flow system, geochemistry, and aquifer hydraulic properties of the shallow groundwater system in the Wilcox and Lorraine process areas of the Wilcox Oil Company Superfund site near Bristow, Oklahoma, 2022

The Wilcox Oil Company Superfund site (hereinafter referred to as “the site”) was formerly an oil refinery northeast of Bristow in Creek County, Oklahoma. Historical refinery operations contaminated the soil, surface water, streambed sediments, alluvium, and groundwater with refined and stored products at the site. The Wilcox and Lorraine process areas are where the highest concentrations of volatile organic compounds, semivolatile organic compounds, polycyclic aromatic hydrocarbons, and trace elements (including metals) (collectively hereinafter referred to as “contaminants”) were measured in a local shallow perched groundwater system within the alluvium (hereinafter referred to as the “alluvial aquifer”) at the site during previous site assessments. In order to understand the potential migration of contaminants through the soil and groundwater in these areas, the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency, investigated aquifer characteristics of the alluvial aquifer in the Wilcox and Lorraine process areas of the site to (1) document hydraulic conductivity and other aquifer characteristics of the alluvial aquifer that govern contaminant fate and transport, (2) describe the geospatial extent and concentration of the contaminants in the alluvial aquifer in the Wilcox and Lorraine process areas, and (3) describe the geochemical controls pertaining to oxidation and reduction governing the fate and transport and the degradation potential of contaminants in the groundwater. Various data were compiled and collected to evaluate the aquifer characteristics at the site including the hydrogeologic framework, groundwater-flow system, geochemistry, and hydraulic properties of the aquifer. A total of 20 new (2022) groundwater monitoring wells were installed at the site to collect data used to supplement groundwater-level altitude and groundwater-quality data collected from older, existing groundwater monitoring wells and piezometers. Data compiled and collected for the study were used to evaluate the characteristics of the alluvial aquifer at the site. These aquifer characteristics are defined by the hydrogeologic framework, groundwater-flow system, geochemistry, and hydraulic properties of the aquifer.

Oklahoma↗

Characterization of major lithologic units underlying the lower American River using water-borne continuous resistivity profiling, Sacramento, California, June 2008

The levee system of the lower American River in Sacramento, California, is situated above a mixed lithology of alluvial deposits that range from clay to gravel. In addition, sand deposits related to hydraulic mining activities underlie the floodplain and are preferentially prone to scour during high-flow events. In contrast, sections of the American River channel have been observed to be scour resistant. In this study, the U.S. Geological Survey, in cooperation with the U.S. Army Corps of Engineers, explores the resistivity structure of the American River channel to characterize the extent and thickness of lithologic units that may impact the scour potential of the area. Likely lithologic structures are interpreted, but these interpretations are non-unique and cannot be directly related to scour potential. Additional geotechnical data would provide insightful data on the scour potential of certain lithologic units. Additional interpretation of the resistivity data with respect to these results may improve interpretations of lithology and scour potential throughout the American River channel and floodplain. Resistivity data were collected in three profiles along the American River using a water-borne continuous resistivity profiling technique. After processing and modeling these data, inverted resistivity profiles were used to make interpretations about the extent and thickness of possible lithologic units. In general, an intermittent high-resistivity layer likely indicative of sand or gravel deposits extends to a depth of around 30 feet (9 meters) and is underlain by a consistent low-resistivity layer that likely indicates a high-clay content unit that extends below the depth of investigation (60 feet or 18 meters). Immediately upstream of the Watt Avenue Bridge, the high-resistivity layer is absent, and the low-resistivity layer extends to the surface where a scour-resistant layer has been previously observed in the river bed.

California↗

Documentation of linear regression models for computing water-quality constituent concentrations using continuous real-time water-quality data for the North Fork Ninnescah River and Cheney Reservoir, Kansas, 2014–21

Cheney Reservoir, in south-central Kansas, was constructed to provide a reliable municipal water supply for the city of Wichita, Kansas, and to provide downstream flood control, wildlife habitat, and recreation. Cheney Reservoir will continue to be important for municipal water supply use as needs increase with ongoing population growth and urban development. Advanced notification of changing water-quality conditions near water-treatment facility intakes and in source waters allows water-treatment facilities and resource planning officials to proactively monitor changing conditions. The U.S. Geological Survey (USGS), in cooperation with the City of Wichita, collected water-quality data at the North Fork Ninnescah River above Cheney Reservoir (USGS station 07144780) and Cheney Reservoir near Cheney, Kans. (USGS station 07144790), monitoring sites to update and develop regression models relating continuous water-quality constituents, streamflow, reservoir storage, and seasonal components to discretely sampled water-quality constituent concentrations of interest. Linear regression analysis was used to update and develop models for alkalinity, major ions, nutrients (nitrogen and phosphorus species), total and dissolved organic carbon, total suspended solids, suspended sediment, fecal indicator bacteria, and atrazine at the North Fork Ninnescah River site and total dissolved solids, major ions, hardness as calcium carbonate, nutrients (nitrogen and phosphorus species), chlorophyll a , and suspended sediment at the Cheney Reservoir site. New and updated models for both sites are applicable to the period of YSI EXO water-quality monitor and sensor deployment (November 14, 2015, through September 30, 2021, at the North Fork Ninnescah River site; October 1, 2014, through September 30, 2021, at the Cheney Reservoir site). Models and resulting water-quality information included in this report can be used in real time, potentially as guidance for water-treatment processes, and can be used to characterize changes in water-quality conditions over time in Cheney Reservoir and its contributing drainage basin provided that the deployed equipment, sensors, and location do not change.

Kansas↗

Escherichia coli monitoring and assessment in 2022 and 2023 after beach restoration at Lake St. Clair Metropark Beach, Macomb County, Michigan

Lake St. Clair Metropark Beach in Michigan has a history of closures because of elevated Escherichia coli ( E. coli ) concentrations in its recreational waters. To reduce closures, restoration projects were implemented in 2021 to deter waterfowl from congregating on the beach. In this study, the U.S. Geological Survey, in cooperation with the Michigan Department of the Environment, Great Lakes, and Energy and in collaboration with Huron-Clinton Metroparks and the Macomb County Health Department, monitored E. coli from 2022–23 in surface water, shallow groundwater, and sediment at Lake St. Clair Metropark Beach. Results were compared to data from a prerestoration (2018–19) study. A significant decrease in daily geometric mean E. coli concentrations in surface water was observed postrestoration, but the number of high concentration events increased. This resulted in more frequent beach closures postrestoration. Surface-sediment E. coli concentrations significantly decreased after restoration, and waterfowl populations generally decreased from 2021 to 2023, suggesting that the deterrence measures could be influencing E. coli concentrations in surface sediments and surface water. Groundwater E. coli concentrations were orders of magnitude higher than those in surface water and revealed no change correlated with restoration. Seepage measurements indicated that groundwater occasionally discharges into surface water, potentially providing a transport mechanism for E. coli to reach the lake. Continued monitoring and consideration of environmental factors could help to better understand the beach system.

Michigan↗

Simulation of flow and sediment mobility using a multidimensional flow model for the White Sturgeon critical-habitat reach, Kootenai River near Bonners Ferry, Idaho

In 1994, the Kootenai River white sturgeon (Acipenser transmontanus) was listed as an Endangered Species as a direct result of two related observations. First, biologists observed that the white sturgeon population in the Kootenai River was declining. Second, they observed a decline in recruitment of juvenile sturgeon beginning in the 1950s with an almost total absence of recruitment since 1974, following the closure of Libby Dam in 1972. This second observation was attributed to changes in spawning and (or) rearing habitat resulting from alterations in the physical habitat, including flow regime, sediment-transport regime, and bed morphology of the river. The Kootenai River White Sturgeon Recovery Team was established to find and implement ways to improve spawning and rearing habitat used by white sturgeon. They identified the need to develop and apply a multidimensional flow model to certain reaches of the river to quantify physical habitat in a spatially distributed manner. The U.S. Geological Survey has addressed these needs by developing, calibrating, and validating a multidimensional flow model used to simulate streamflow and sediment mobility in the white sturgeon critical-habitat reach of the Kootenai River. This report describes the model and limitations, presents the results of a few simple simulations, and demonstrates how the model can be used to link physical characteristics of streamflow to biological or other habitat data. This study was conducted in cooperation with the Kootenai Tribe of Idaho along a 23-kilometer reach of the Kootenai River, including the white sturgeon spawning reach near Bonners Ferry, Idaho that is about 108 to 131 kilometers below Libby Dam. U.S. Geological Survey's MultiDimensional Surface-Water Modeling System was used to construct a flow model for the critical-habitat reach of the Kootenai River white sturgeon, between river kilometers 228.4 and 245.9. Given streamflow, bed roughness, and downstream water-surface elevation, the model computes the velocity field, water-surface elevations, and boundary shear stress throughout the modeled reach. The 17.5 kilometer model reach was subdivided into two segments on the basis of predominant grain size: a straight reach with a sand, gravel, and cobble substrate located between the upstream model boundary at river kilometer 245.9 and the upstream end of Ambush Rock at river kilometer 244.6, and a meandering reach with a predominately sand substrate located between upstream end of Ambush Rock and the downstream model boundary at river kilometer 228.4. Model cell size in the x and y (horizontal) dimensions is 5 meters by 5 meters along the computational grid centerline with 15 nodes in the z (vertical) dimension. The model was calibrated to historical streamflows evenly distributed between 141.6 and 2,548.9 cubic meters per second. The model was validated by comparing simulated velocities with velocities measured at 15 cross sections during steady streamflow. These 15 cross sections were each measured multiple (7-13) times to obtain velocities suitable for comparison to the model results. Comparison of modeled and measured velocities suggests that the model does a good job of reproducing flow patterns in the river, although some discrepancies were noted. The model was used to simulate water-surface elevation, depth, velocity, bed shear stress, and sediment mobility for Kootenai River streamflows of 170, 566, 1,130, 1,700, and 2,270 cubic meters per second (6,000, 20,000, 40,000, 60,000, and 80,000 cubic feet per second). The three lowest streamflow simulations represent a range of typical river conditions before and since the construction of Libby Dam, and the highest streamflow simulation (2,270 cubic meters per second) is approximately equal to the annual median peak streamflow prior to emplacement of Libby Dam in 1972. Streamflow greater than 566 cubic meters per second were incrementally increased by 570 cubic meters per second. For each

Scientific Investigations Report↗