Geology ReportsSearch

SEARCH · Geology Reports

Results for “Annual Review of Earth and Planetary Sciences”

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.

1,695 records · Page 53Linked to original sources

Use of a ground-penetrating radar system to detect pre- and post-flood scour at selected bridge sites in New Hampshire, 1996-98

Ground-penetrating radar was used to measure the depth and extent of existing and infilled scour holes and previous scour surfaces at seven bridges in New Hampshire from April 1996 to November 1998. Ground-penetrating-radar survey techniques initially were used by the U.S. Geological Survey to study streambed scour at 30 bridges. Sixteen of the 30 bridges were re-surveyed where floods exceeded a 2-year recurrence interval. A 300-megahertz signal was used in the ground-penetrating radar system that penetrated through depths as great as 20 feet of water and as great as 32 feet of streambed materials. Existing scour-hole dimensions, infilled thickness, previous scour surfaces, and streambed materials were detected using ground-penetrating radar. Depths to riprap materials and pier footings were identified and verified with bridge plans. Post data-collection-processing techniques were applied to assist in the interpretation of the data, and the processed data were displayed and printed as line plots. Processing included distance normalization, migration, and filtering but processing was kept to a minimum and some interference from multiple reflections was left in the record. Of the 16 post-flood bridges, 22 ground-penetrating-radar cross sections at 7 bridges were compared and presented in this report. Existing scour holes were detected during 1996 (pre-flood) data collection in nine cross sections where scour depths ranged from 1 to 3 feet. New scour holes were detected during 1998 (post-flood) data collection in four cross sections where scour depths were as great as 4 feet deep. Infilled scour holes were detected in seven cross sections, where depths of infilling ranged from less than 1 to 4 feet. Depth of infilling by means of steel rod and hammer was difficult to verify in the field because of cobble and boulder streambeds or deep water. Previous scour surfaces in streambed materials were identified in 15 cross sections and the depths to these surfaces ranged from 1 to 10 feet below the streambed. Riprap materials or pier footings were identified in all cross sections. Calculated record depths generally agree with bridge plans. Pier footings were exposed at two bridges and steel pile was exposed at one bridge. Exposures were verified by field observations.

New Hampshire

Land application of drill waste: A scope analysis

Drilling fluid waste land application, a process where drilling wastes are spread and tilled into the land surface, has become common in some petroleum-producing states, however, the potential benefits and risks of this practice are not well studied. Drilling fluids can be water- or oil-based and can have high concentrations of total soluble salts and total petroleum hydrocarbons. Comprehensive chemical characterization of these fluids is not well documented in the literature, and the extent of land application is largely unknown. We hypothesized that the land application of drill waste would fluctuate over time due to economic factors. To begin to understand the extent of historical and potential future land application, we analyzed data from over 5,800 drilling fluid land application permits collected by the Oklahoma Corporation Commission for years 2000, 2005, 2010, and 2015–2020. During the years studied, drilling fluid wastes were applied to more than 250,000 acres in Oklahoma, with over 54,000 thousand barrels (Mbbl) of liquids and nearly 21,000 Mbbl of solids applied. Land application is widespread (occurring in 59/77 counties), however recent drilling activity, land availability, and the economics of transportation have created conditions favorable for land application specifically in the Anadarko Basin. Land application can co-occur with sensitive areas, such as important groundwater and surface-water drinking sources and agricultural fields used for subsistence or feed crop production. Our approach for quantifying the extent of land application, along with further chemical characterization studies, can aid operators and land managers who are considering this practice in assessing the associated benefits and risks.

Oklahoma

Utilizing hydrophones to detect streambed mobilization in the Wild and Scenic reach of the Rio Chama

This paper describes a study that deployed underwater acoustic sensors (hydrophones) to detect streambed mobilization in the Wild and Scenic reach of the Rio Chama, New Mexico. Fine sediment accumulation in the gravel‑cobble bed affects brown trout spawning habitat, and understanding incipient motion thresholds is essential for planning effective environmental flushing flows released from El Vado Dam. Two hydrophone stations were operated during 2021–2022, capturing audio data during a high‑flow pulse in late 2021. Acoustic processing techniques, including spectral analysis and threshold‑based impact detection, were used to identify sediment‑generated noise associated with coarse sediment transport. Results indicate detectable sediment mobilization during the high‑flow pulse, demonstrating hydrophones’ potential as a passive monitoring tool for environmental flow planning and sediment‑management decision support.

Conference Paper

Hydraulic properties of the Madison aquifer system in the western Rapid City area, South Dakota

Available information on hydrogeology, data from borehole geophysical logs, and aquifer tests were used to determine the hydraulic properties of the Madison aquifer. From aquifer-test analysis, transmissivity and storage coefficient were determined for the Minnelusa and Madison aquifers, and vertical hydraulic conductivity (K v ') along with specific storage (Ss') for the Minnelusa confining bed. Borehole geophysical well logs were used to determine the thickness and location of the Minnelusa aquifer, the lower Minnelusa confining bed, and the Madison aquifer within the Madison Limestone. Porosity values determined from quantitative analysis of borehole geophysical well logs were used in analyzing the aquifer-test data. The average porosity at the two aquifer-test sites is about 10 percent in the Minnelusa aquifer, 5 percent in the lower Minnelusa confining bed, and 35 percent in the Madison aquifer. The first aquifer test, which was conducted at Rapid City production well #6, produced measured drawdown in the Minnelusa and Madison aquifers. Neuman and Witherspoon's method of determining the hydraulic properties of leaky two-aquifer systems was used to evaluate the aquifer-test data by assuming the fracture and solution-opening network is equivalent to a porous media. Analysis of the aquifer test for the Minnelusa aquifer yielded a transmissivity value of 12,000 feet squared per day and a storage coefficient of 3 x 10 -3 . The specific storage of the Minnelusa confining bed was 2 x 10 -7 per foot, and its vertical hydraulic conductivity was 0.3 foot per day. The transmissivity of the Madison aquifer at this site was 17,000 feet squared per day, and the storage coefficient was 2 x 10 -3 . The second aquifer test, which was conducted at Rapid City production well #5 (RC-5) produced measured drawdown only in the Madison aquifer. Hantush and Jacob's method of determining the hydraulic properties of leaky confined aquifers with no storage in the confining bed was used to evaluate the aquifer-test data by assuming the fracture and solution-opening network is equivalent to a porous media. The analysis of data from the RC-5 aquifer test showed that transmissivity was not equal in all directions. Hantush's method was used to determine the direction of radial anisotropy and magnitude of the major and minor axes of transmissivity. The major axis of transmissivity is at an angle of 42° east of north, and the transmissivity along this axis is about 56,000 feet squared per day. The minor axis of transmissivity is at an angle of 48° west of north, and the transmissivity along this axis is about 1,300 feet squared per day. The major axis of transmissivity intersects Cleghorn Springs, a large resurgent spring on the west edge of Rapid City. The shape of the potentiometric contours of the Madison aquifer near RC-5 agree with the orientation of the transmissivity ellipse. The average value of the storage coefficient from the isotropic analysis of the aquifer-test data was 3.5 x 10 -4 , and the average vertical hydraulic conductivity of the lower Minnelusa confining bed was 9.6 x 10 -3 foot per day.

South Dakota

Stream sediment sources in Medicine Creek, northern Missouri and southern Iowa

This report presents the results of a cooperative study by the U.S. Geological Survey and Missouri Department of Natural Resources to quantify sediment transport source contributions in the Medicine Creek drainage basin. Understanding relative source contributions provides valuable information for selecting the conservation practices that may be most effective in reducing sediment and sediment-associated nutrient transport in the Medicine Creek drainage basin and similar areas of the Lower Grand River drainage basin. Sediment samples were collected from potential contributing areas (source samples) and from fluvial-transported samples (target samples). Source sample types included streambanks, row crop fields, and a combined pastures and forests category. Samples were analyzed for particle size and quantity of carbon, nitrogen, stable isotopes of carbon and nitrogen, and 49 mineral elements as potential tracers. Results for the carbon stable isotope ratio of carbon-13/carbon-12 (δ 13 C) and concentrations of total carbon, total nitrogen, calcium, potassium, and copper were selected by discriminant function analysis as the best combination of multiple tracers to differentiate each source type. The discriminant function analysis poorly differentiated pastures and forests, so these source types were combined. The sources defined by the discriminant function analysis were then used in an unmixing model to apportion sources for each target sample. In the study area, transported sediment was predominantly bank sediment, with an overall average of 86.9 percent of suspended-sediment samples and depositional streambed samples attributed to bank material. Suspended-sediment samples from the mainstem of Medicine Creek were dominated by bank sediments (average of 95.8 percent), and depositional streambed samples from throughout the drainage basin had more variable source contributions with an average of 71.1 percent attributed to bank material. The relative importance of upland sources (row crop fields and the combined pastures and forests category) varied seasonally and with streamflow but was not related to land use or drainage basin size. Relative contributions from upland sources were greater in the summer through winter rather than spring and during lower streamflow, though this may be driven by the seasonality of streamflow. These results indicate management practices that reduce bank erosion could be effective strategies for managing the dominant source of sediment and sediment-associated phosphorus.

Iowa, Missouri

Transmissivity and water quality of water-producing zones in the intermediate aquifer system, Sarasota County, Florida

The intermediate aquifer system is an important water source in Sarasota County, Florida, because the quality of water in it is usually better than that in the underlying Upper Floridan aquifer. The intermediate aquifer system consists of a group of up to three water-producing zones separated by less-permeable units that restrict the vertical movement of ground water between zones. The diverse lithology, that makes up the intermediate aquifer system, reflects the variety of depositional environments that occurred during the late Oligocene and Miocene epochs. Slight changes in the depositional environment resulted in aquifer heterogeneity, creating both localized connection between water-producing zones and abrupt culmination of water-producing zones that are not well documented. Aquifer heterogeneity results in vertical and areal variability in hydraulic and water-quality properties. The uppermost water-producing zone is designated producing zone 1 but is not extensively used because of its limited production capability and limited areal extent. The second water-producing zone is designated producing zone 2, and most of the domestic- and irrigation-supply wells in the area are open to this zone. Additionally, producing zone 2 is utilized for public supply in southern coastal areas of Sarasota County. Producing zone 3 is the lowermost and most productive water-producing zone in the intermediate aquifer system. Public-supply well fields serving the cities of Sarasota and Venice, as well as the Plantation and Mabry Carlton Reserve well fields, utilize producing zone 3. Heads within the intermediate aquifer system generally increase with aquifer depth. However, localized head-gradient reversals occur in the study area, coinciding with sites of intense ground-water withdrawals. Heads in producing zones 1, 2, and 3 range from 1 to 23, 0.2 to 34, and 7 to 42 feet above sea level, respectively. Generally, an upward head gradient exists between producing zones 3 and 2. However, an upward head gradient between producing zones 2 and 1 does not consistently occur throughout Sarasota County, probably the result of greater ground-water withdrawals from producing zone 2 than from producing zone 1. The transmissivity of the intermediate aquifer system is spatially variable. Specific-capacity data from selected wells penetrating producing zones 2 and 3, were used to estimate transmissivity. Estimated transmissivity values for producing zones 2 and 3 range from about 100 to 26,000 feet squared per day and from about 1,300 to 6,200 feet squared per day, respectively. Because the capacity of specific water-producing zones is highly variable from site to site, estimating the performance of a specific water-producing zone as a water resource is difficult. Water samples collected during the study were analyzed for major-ion concentrations. Generally, bicarbonate type water from rock interaction occurs in northern Sarasota County; enriched calcium-magnesium-sulfate type water from deeper aquifers occurs in central Sarasota County; and sodium-chloride type water from saltwater mixing occurs in southern Sarasota County. In some areas of northern Sarasota County, the major-ion concentrations in water are lower in producing zone 2 than in producing zone 1. Major-ion concentrations in water are higher in producing zone 3 throughout the study area. A major objective of the study was to evaluate hydraulic and water-quality data to determine distinctions that could be used to characterize a particular producing zone. However, data indicate that both hydraulic and water-quality properties are highly variable within and between zones, and are more related to the degree of connection between and areal extent of water-producing zones than to aquifer depth and distance from the coast.

Florida

Geochemistry of quartz and fluid inclusions associated with gold, tungsten, and antimony in the Stibnite-Yellow Pine district, Idaho

Deposits in the historic Stibnite-Yellow Pine district in east-central Idaho were mined for Ag, Au, Hg, Sb, and W. Ore is hosted in breccia and fault zones along NE-striking faults that transect the Atlanta lobe of the Late Cretaceous Idaho batholith and metamorphosed Neoproterozoic to lower Paleozoic rocks. Geochronological studies have delineated five superimposed and zoned hydrothermal events in the mining district which range in age from Late Cretaceous to the Eocene. They formed the following consecutive assemblages: (1) vein and disseminated pre-ore milky quartz (± calcite), muscovite, and brassy pyrite with minor molybdenite and base metal sulfides; (2) vein and disseminated quartz, calcite, sericite, auriferous-pyrite, and arsenopyrite; (3) potassium feldspar, quartz, dolomite, and calcite with vein and disseminated scheelite ± pyrite and arsenopyrite; (4) Au-Ag veins that contain adularia; and (5) breccia, vein, and disseminated quartz, calcite, and stibnite with minor Ag sulfosalts and distal zones of cinnabar and realgar at Fern and Hermes mines. In this study, the pressure, temperature, and composition of ore-forming fluids were measured in Au-, W-, and Sb-bearing samples collected from drill core and surface exposures. Using optical petrography, scanning electron microscopy cathodoluminescence (SEM-CL) and electron microprobe cathodoluminescence (EMPA-CL) imaging, as well as electron microprobe (EMP) and laser ablation-inductively coupled plasma-mass spectrometry trace element analyses, four major overprinting generations of hydrothermal quartz are defined. Fluid inclusion assemblages (FIAs) in each generation of quartz were characterized by microthermometry and laser Raman spectroscopy. Quartz and fluid inclusion paragenesis are as follows: (Qi) altered and overprinted Qi with no coeval FIAs that formed before or during event 1; (Q1) early veins as well as disseminated quartz with dull gray SEM-CL and H 2 O-NaCl-CO 2 FIAs that formed at minimum temperatures and pressures of 281° to 358°C and 1346 to 1789 bar during event 2 (Au); (Q2) vein and disseminated quartz, often euhedral, with H 2 O-NaCl ± CO 2 ± CH 4 FIAs that formed at minimum temperatures and pressures of 183° to 213°C and 72 to 602 bar during event 3 (W); and (Q3) euhedral and disseminated quartz with H 2 O-NaCl FIAs that formed at minimum temperature and pressures of 156° to 161°C and 5 to 6 bar during event 5 (Sb). Stibnite hosts H 2 O-NaCl FIAs trapped at minimum temperatures and pressures of 150° to 189°C and 5 to 12 bar that are similar to those of Q3-hosted FIAs. The quartz paragenetic history in the district follows a progressive temporal path from Tienriched quartz deposited by a relatively deep and higher temperature aqueous-carbonic fluid to a Li, Al, and Sb enriched quartz precipitated from a lower temperature aqueous fluid. Gases in fluid inclusion extracts from ore and gangue minerals were analyzed by mass spectrometry. Noble gas isotope analysis of fluid inclusion extracts revealed three signatures. Pre-ore Qi has R/R A (the isotopic ratio of the sample 3 He/ 4 He to atmospheric 3 He/ 4 He (1.384 × 10 −6 )) values that reflect the composition of the crustal melt-derived Atlanta lobe of the Idaho batholith. Quartz (Q1) and auriferous arsenian pyrite of hydrothermal event 2 have 20 Ne/ 4 He and R/R A values that plot along the air-crust mixing line, and these gases were likely introduced by metamorphic fluids derived from metasedimentary rocks. Stibnite with associated quartz (Q3 ± Q2) have 20 Ne/ 4 He and R/R A values that extend from the air-crust mixing line to the air-mantle mixing line, showing a shift from a crustal signature to an overprinting mantle signature. Late calcite has 20 Ne/ 4 He and R/R A values that extend from air toward the mantle end member, which suggests that magmatic volatiles condensed into meteoric groundwater late in the paragenesis. A mantle volatile source is also supported by high 20 Ne/ 22 Ne ratios. The Stibnite-Yellow Pine district is similar to other Au-W-Sb districts around the world in ore fluid chemistries, with late-stage stibnite being a product of cooling. The fluid pressure-temperature-composition (P-T-X) estimates and the crustal source of He during hydrothermal event 2 gold mineralization are characteristic of orogenic Au deposits, and the observed inclusion reequilibration textures indicate that trapped fluid inclusions were deformed during transitions between the brittle-ductile environment. Scheelite was deposited during event 3 with Q2 quartz and is overgrown by Q3 quartz and stibnite. Our interpretation of the quartz texture, chemistry, and fluid inclusion P-T-X data suggest that Au, W, and Sb mineralization occurred as hydrothermal fluids evolved from H 2 O-NaCl-CO 2 at >300°C and depths of ~4 to 6 km near the ductile to brittle transition to H 2 O-NaCl at ~150°C and epithermal depths as the Idaho batholith was exhumed. Geochronology indicates that from the Late Cretaceous to early Eocene, overprinting hydrothermal events formed deep orogenic Au, intermediate W, and shallow epithermal Au-Ag and Sb. Without comprehensive geochronologic constraints, our fluid inclusion data and textural observations could be interpreted as a single evolving hydrothermal system, as has been done in many other mining districts with similar mineralization.

Idaho

Water-quality conditions and relation to drainage-basin characteristics in the Scituate Reservoir Basin, Rhode Island, 1982-95

The Scituate Reservoir Basin covers about 94 square miles in north central Rhode Island and supplies more than 60 percent of the State of Rhode Island's drinking water. The basin includes the Scituate Reservoir Basin and six smaller tributary reservoirs with a combined capacity of about 40 billion gallons. Most of the basin is forested and undeveloped. However, because of its proximity to the Providence, Rhode Island, metropolitan area, the basin is subject to increasing development pressure and there is concern that this may lead to the degradation of the water supply. Selected water-quality constituent concentrations, loads, and trends in the Scituate Reservoir Basin, Rhode Island, were investigated locate parts of the basin likely responsible for exporting disproportionately large amounts of water-quality constituents to streams, rivers, and tributary reservoirs, and to determine whether water quality in the basin has been changing with time. Water-quality data collected between 1982 and 1995 by the Providence Water Supply Board PWSB) in 34 subbasins of the Scituate Reservoir Basin were analyzed. Subbasin loads and yields of total coliform bacteria, chloride, nitrate, iron, and manganese, estimated from constituent concentrations and estimated mean daily discharge records for the 1995 water year, were used to determine which subbasins contributed disproportionately large amounts of these constituents. Measurements of pH, color, turbidity, and concentrations of total coliform bacteria, sodium, alkalinity, chloride, nitrate, orthophosphate, iron, and manganese made between 1982 and 1995 by the PWSB were evaluated for trends. To determine the potential effects of human-induced changes in drainage- basin characteristics on water quality in the basin, relations between drainage-basin characteristics and concentrations of selected water-quality constituents also were investigated. Median values for pH, turbidity, total coliform bacteria, sodium, alkalinity, chloride, nitrate, and iron were largest in subbasins with predominately residential land use. Median instantaneous loads reflected drainage-basin size. However, loads normalized by drainage area (median instantaneous yields) also were largest in residential areas where point and non-point sources are likely, and in areas of poorly drained soils. Significant trends in water-quality constituents from 1982 to 1995 in the Scituate Reservoir Basin indicate that the quality of the water resources in the basin may be slowly changing. Scituate Reservoir subbasins with large amounts of residential land use showed increasing trends in alkalinity and chloride. In contrast, subbasins distributed throughout the drainage basin showed increasing trends in pH, color, nitrate, and iron concentrations, indicating that these characteristics and constituents may be affected more by atmospheric deposition. Although changing, water-quality constituent concentrations in the Scituate Reservoir Basin only occasionally exceeded Rhode Island and USEPA water-quality guidelines and standards. Result of correlation analysis between pH, color, turbidity, and concentrations total coliform bacteria, sodium, alkalinity, chloride, nitrate, orthophosphate, iron, and manganese and land use, geology, wetlands, slope, soil drainability, and roads indicated that the percentage of wetlands, roads, and slope appear have the greatest effect on water-quality in the Scituate Reservoir Basin. The percentage of urban, residential, and commercial land use also are important, but to a lesser degree than wetlands, roads, and slope. Finally, geology appears to have the least effect on water quality compared to other drainage-basin characteristics investigated.

Rhode Island

Seismic response of rock towers at the Trona Pinnacles (U.S.A.) to the 2019 Ridgecrest earthquake sequence: Theory, observations, and models

We analyze the seismic response of a class of fragile geologic features (FGFs), referred to as rock towers (RTs) at the Trona Pinnacles, a group of RTs in southern California that suffered strong shaking during the 2019 Ridgecrest earthquake sequence. FGFs, including RTs, may provide maximum constraints on past earthquake shaking intensity, and thereby support probabilistic seismic hazard assessments (PSHAs). In a rare case study, we explore the hypothesis that RT structural integrity is time dependent, as damage accumulates progressively. We develop finite‐element method (FEM) models of the RTs using photogrammetric shape models. We validate the models by comparing numerical simulations of their response to broadband ground shaking with low‐intensity seismic recordings obtained at the Pinnacles. Results of our simulations are in good agreement with the seismic recordings of actual earthquake aftershocks. We next use the results of the FEM models to analyze the response and evolution of RTs. Our analyses elucidate the influence of geometry over their seismic response, providing a rationale that may explain the rarity of slender RTs at Trona: high‐aspect‐ratio structures that respond in bending develop detrimental tensile stresses that crack the rock, whereas low‐aspect‐ratio ones’ response also includes shearing, which does not compromise material integrity as much as tension. Field measurements with a rebound hammer support this finding, suggesting that the material around the base of slender rocks has been weakened relative to other parts of the RT. We also study how to define simplified mechanical models (“archetypes”) to predict the natural frequencies of RTs. Results from our work illuminate the fundamental mechanisms of seismic response and progressive failure of RTs, and open new avenues of research to potentially incorporate these geologic features as long‐return period constraints on PSHA, in ways analogous to those of the widely used precariously balanced rocks.

California

Flood-inundation maps of the Blue River and major tributaries in and near Kansas City, Missouri, 2023–25

Digital flood-inundation maps for 35.5 miles of the Blue River, in and near Kansas City, Missouri, were created by the U.S. Geological Survey (USGS) in cooperation with the City of Kansas City. Maps were also developed for a combined 7.3 miles of Brush Creek (4.4 miles) and Indian Creek (2.9 miles), two primary tributaries to the Blue River in the study reach. The flood-inundation maps, available through the USGS Flood Inundation Mapping Program website at https://water.usgs.gov/osw/flood_inundation/ , depict estimates of the areal extent and depth of flooding, depth of water above first-floor elevations in selected buildings, water velocity, water-surface elevations, and time of travel corresponding to selected water levels or gage heights (hereafter referred to as “stages”) at nine reference USGS streamgages—five on the Blue River (USGS station numbers 06893590, 06893578, 06893500, 06893150, and 06893100), two on Brush Creek (USGS station numbers 06893562 and 06893557), one on Indian Creek (USGS station number 06893390), and one on the Missouri River (USGS station number 06893000). Near-real-time stages at these streamgages may be obtained from USGS Water Data for the Nation at https://doi.org/10.5066/F7P55KJN or the National Weather Service National Water Prediction Service at http://water.noaa.gov/ , which also forecasts flood hydrographs at seven of these sites (USGS station numbers 06893000, 06893590, 06893578, 06893500, 06893150, 06893557, and 06893390). Flood profiles were computed for eight map reaches by means of two-dimensional hydraulic models. The models were calibrated using the USGS-developed stage–streamflow relations at each reference streamgage. Two map reaches on the lower Blue River also include the effects of backwater from the Missouri River on flood inundation extent and hydraulic characteristics. The hydraulic models were used to compute water-surface profiles for flood stages at 1-foot intervals referenced to the streamgage datums. The profile stages ranged from the National Weather Service “Action Stage” or near bankfull to a stage exceeding the highest recorded water level at each streamgage. The simulated water-surface profiles were then combined with a digital elevation model (derived from light detection and ranging data having a nonvegetated vertical accuracy of a maximum 10-centimeter root-mean-square error) to delineate the area flooded at each water level and the associated hydraulic characteristics. The availability of these maps, along with information regarding current stage from the USGS streamgage and forecasted high-flow stages from the National Weather Service, will provide emergency management personnel, resource managers, and residents with information that could be critical for flood-response activities such as evacuations and road closures, as well as for postflood recovery efforts.

Kansas, Missouri

High resolution identification and quantification of diffuse deep groundwater discharge in mountain rivers using continuous boat-mounted helium measurements

Discharge of deeply sourced groundwater to streams is difficult to locate and quantify, particularly where both discrete and diffuse discharge points exist, but diffuse discharge is one of the primary controls on solute budgets in mountainous watersheds. The noble gas helium is a unique identifier of deep groundwater discharge because groundwater with long residence times is commonly enriched in helium. In this study, a portable mass spectrometer was used to measure longitudinal variation in dissolved helium concentrations in two mountainous rivers at high spatial resolution not feasible with traditional sampling techniques. Helium profiles were then simulated using a mass-balance model to quantify longitudinal variation in groundwater discharge to the receiving rivers. Results indicate helium concentrations were enriched by multiple orders of magnitude above atmospheric equilibrium in both rivers and that this persisted for up to 18 km below observed pulse inputs in the Colorado River. Helium mass-balance models match observed longitudinal patterns with the exception of sharp initial increases in helium observed in the rivers. Increased longitudinal groundwater discharge rates correspond to mapped geologic structures in both watersheds that likely transport deep geothermal water. Models show variable sensitivity to spatial assignment of input variables representing the groundwater source, illustrating the importance of collecting data from discrete groundwater discharges where possible. The methodology shows promise for field experiments designed to assess air–water exchange rates and to quantify total groundwater discharge from a combination of discrete and diffuse sources.

Colorado, Utah

A joint Gaussian process model of geochemistry, geophysics, and temperature for groundwater TDS in the San Ardo Oil Field, California, USA

Decline in availability of fresh groundwater has expanded interest in brackish groundwater resources; however, the distribution of brackish groundwater is poorly understood. Water resources in sedimentary basins across the United States often overlie oil and gas development. Mapping of groundwater total dissolved solids (TDS) using data from oil well geophysical logs has become an important technique for identifying fresh and brackish groundwater. Existing geophysical log analysis methods use porosity and temperature to relate formation resistivity to TDS. Typically, natural geothermal gradients are used to estimate temperature at the location of collected resistivity. However, in thermally enhanced oil fields, steam is injected into the subsurface to mobilize high viscosity oil, creating variable temperature distributions. Furthermore, TDS derived from resistivity also depends on the fractions of dominant ions. Typically, chloride and bicarbonate fractions must be determined. It is also necessary to model TDS across many geologic units with heterogenous porosity distributions. Collectively, each quantity used to estimate TDS (resistivity, porosity, temperature, bicarbonate fraction) varies in space and time, and available data points are rarely collocated. Here, we present a new method of mapping groundwater TDS that continuously models each quantity together with a joint Gaussian process. This method enables mapping fresh and brackish water with practically available data. We apply this method to the San Ardo Oil Field in Monterey County, California, where steam injection occurs. In some areas of the aquifer system overlying the oil zone, the temperature is ∼75 °C, roughly twice the natural background value. Groundwater TDS is typically <1,500 mg/L in the aquifer and increases with depth to ∼9,000 mg/L in the oil-producing zone. A low-permeability clay layer delineates the fresh and brackish water, likely by inhibiting surface recharge from penetrating the deeper zones, allowing higher-TDS connate water to remain in place. Weaker lateral TDS trends may be controlled by recharge patterns associated with the Salinas River. Our model reveals with high certainty that groundwater has freshened in one localized part of the oil-producing zone and suggests with less certainty that more widespread freshening has also occurred. The lowering of TDS was possibly from decades of low-TDS steam injection and the associated fluid production and disposal operations.

California

Assessing spatial variability of nutrients, phytoplankton, and related water-quality constituents in the California Sacramento–San Joaquin Delta at the landscape scale—2018 high resolution mapping surveys

Executive Summary This study examined the abundance and distribution of nutrients and phytoplankton in the tidal aquatic environments of the Sacramento–San Joaquin Delta (Delta) and Suisun Bay, comprising three spatial surveys conducted in May, July, and October of 2018 that used continuous underway high frequency sampling and measurements onboard a high-speed boat to characterize spatial variation across the extent of the Delta. The method used involves simultaneously collecting information about the concentration and spatial distribution of all major nutrient forms with analogous information about the major classes of phytoplankton and associated water-quality conditions. The results showed substantial variation across space and time, providing an unprecedented snapshot of the dynamic environmental processes that shape the ways nutrients interact with and affect aquatic habitats in the Delta. The purposes of this study were to improve our understanding of how hydrodynamics, landscape features, and aquatic primary productivity interact to drive nutrient cycling and transport in the Delta and to provide insights into the underlying processes most directly responsible for the conditions at the time of this study, and thus into the range of conditions that may be expected following the wide array of prospective future changes to the Delta. One major anticipated change at the time of this study was the planned upgrade to the Sacramento Regional Wastewater Treatment Plant, but the study also informs our understanding of potential effects from other changes to the Delta, such as those caused by other nutrient-management actions, flow actions, large-scale wetland restoration, drought, flood, levee failure, and changes to water management. Nutrient loading is the primary driver of nutrient concentrations in the Delta, but several other major drivers interact to shape their distribution and effects: geomorphology, hydrodynamics, landscape features, and aquatic productivity. Hydrodynamics affect timescales of transport and dilution of nutrient loads in the Delta. During transit through the system, channel geometry, tidal mixing, and water exports affect hydrodynamics in diverse ways that influence water-residence and transport times, thereby markedly affecting the range of times during which natural internal cycling can alter nutrient concentrations and forms. Channel geometry and location shape tidal energy and river currents into these observed dynamics. Interactions with Delta aquatic landscapes such as herbaceous tidal marsh, submerged aquatic vegetation, and large expanses of intertidal or subtidal sediments (all highly productive landscapes) exert demand on available nutrient supplies but can also simultaneously transform and generate nutrients. Finally, while phytoplankton require nutrients to sustain production and thus are a potential nutrient sink, the amount and form of nutrients also can influence the occurrence of harmful algal blooms (HABs) that adversely affect aquatic organisms as well as affect the occurrence of beneficial algal blooms that result in production of algae that are favorable for imperiled Delta pelagic aquatic food webs. The surveys revealed a complex mosaic of spatial variation, with nutrient concentrations varying from near zero to well above concentrations considered eutrophic; nutrient concentrations were more often related to the extent of hydrologic transport and mixing than to specific geographic locations or to specific landscape features. Similarly, the surveys identified phytoplankton abundance ranging from near detection to the level of large phytoplankton blooms, with large variation in phytoplankton community composition. Although the study occurred during a period of low bloom activity, phytoplankton productivity appeared to be the strongest potential sink for inorganic nutrients in the Delta, indicating that it is a larger control on nutrient concentrations and distribution than previously understood. Cycling and transformation within the water column only appeared to substantially lower total nutrient concentrations at the longest estimated transport timescales. Contrary to expectations, we did not observe substantial nutrient depletion near landscape-scale features such as open-water habitats, submerged aquatic vegetation beds, extensive wetlands, or exposed sediments, indicating that these habitat types did not act as major sinks for nutrients in the Delta during these surveys. These results indicated that nutrient reduction efforts may have the greatest effect on pelagic phytoplankton productivity in the more productive reaches of the Delta and estuary, but these effects are unlikely to be magnified by changes to nutrient loss within the Delta over conceivable changes in flow conditions, Delta water management actions, or large-scale wetland restoration activities. Nevertheless, local processes were shown to cause substantial loss, and thus integrating of nutrient effects with other indicators of aquatic habitat conditions will help inform planning future actions at specific sites. Finally, we note that the primary contribution of this study was intended to be the survey data themselves. Aside from the results highlighted in this report, the surveys are a benchmark against which future environmental change may be evaluated, including changes to nutrient management or water exports, drought, large-scale wetland restoration, and climate change. Further, although we highlight some of the main findings from the surveys in this report, the necessarily limited scope precludes examination of many topics for which these surveys may be highly informative. To facilitate the utility of these data to stakeholders, managers, and researchers, we have released the data online (Bergamaschi and others, 2020) and created an online data exploration portal ( https:​//ca.water​.usgs.gov/​bay-​delta/​2018-​delta-​wide-​mapping-​surveys.html ) where users may query the surveys in a variety of ways to test hypotheses, examine relationships, assess spatial trends, and download data. The data exploration portal is intended to be an immersive experience that allows users to gain greater understanding of the complex interactions that shape Delta aquatic environments. This report is intended as a companion to the portal, allowing the reader to challenge and further explore the highlighted findings. This study was a collaboration between the U.S. Geological Survey and the Delta Regional Monitoring Program, with additional funding provided from U.S. Geological Survey Cooperative Matching Funds Program.

California

Revisiting an enigma on California's north coast: The Mw6.5 Fickle Hill earthquake of 21 December 1954

Many earthquakes occur along the North Coast of California in the vicinity of the Mendocino Triple Junction (MTJ), where the Pacific, Gorda, and North American (NA) plates meet, and on the adjacent plate boundaries. The MTJ marks the nexus of the Mendocino and San Andreas faults with the Cascadia subduction zone (CSZ). Historically, most large earthquakes around the MTJ have been within the offshore Gorda plate and its subducted portion beneath the NA plate. North of the MTJ, active faults mapped in the NA plate are part of the CSZ fold‐and‐thrust belt. Although some events have been detected in the NA plate, no large historic events have been associated with mapped surface faults. The 21 December 1954 M w 6.5 earthquake in Humboldt County is one possible exception. Using published data from catalogs and articles, unpublished data from Berkeley’s archives, and S ‐ P times interpreted from two U.S. Coast and Geodetic Survey (USCGS) accelerometers, we determine a probability cloud for the earthquake’s hypocenter using NonLinLoc. The highest probability location lies beneath Fickle Hill just east of the city of Arcata, California, at 40.87° N, 124.03° W, and ∼11 km depth. Using P ‐wave polarities from Berkeley stations and the digitized waveforms from the accelerometers, we find that the focal mechanism most consistent with the data indicates thrust movement with strike, dip, and rake of 350°, 10°, and 90°, respectively, at a depth of 14 km. Given the depth uncertainties of both this event and the megathrust, this implies that the earthquake most likely took place on the subduction interface rather than on the mapped faults in the Mad River fault zone that trend 322° and dip to the northeast. The revisited intensity in the epicentral region also supports a location beneath Fickle Hill to the east of the city of Arcata, California.

California, Idaho, Nevada, Oregon, Washington

ShakeAlert® and schools: Incorporating earthquake early warning in school districts in Alaska, California, Oregon, and Washington

The U.S. Geological Survey-managed ShakeAlert® earthquake early warning system is the first public alerting system in the United States to provide rapid mass notification when an earthquake is detected. Although public alert delivery via mobile phones began in California in 2019 followed by Oregon and Washington in 2021, little is known about what might drive widespread implementation in at-risk institutional settings such as schools. For example, there has been limited research on how to best integrate earthquake early warning into existing emergency plans, alert systems, and drills to keep school children and personnel safe in an earthquake. To address this gap, in the spring of 2022, every school district superintendent in Alaska, California, Oregon, and Washington was sent a 15-min online survey. The survey assessed superintendent knowledge of ShakeAlert, preferences for alert messaging, and perceived opportunities and barriers to incorporating the system in schools. The results showed that superintendents had low awareness of ShakeAlert but held positive perceptions of the system's potential to enable life-saving protective actions. A major barrier to adoption included the perceived financial cost of implementing and maintaining the system. There were some statistically significant differences in state responses, and future research could investigate the specific needs of each state based on school district size and composition, hazard exposure, and earthquake experience. Together these findings can help inform targeted strategies to increase ShakeAlert adoption in schools and ultimately improve the safety of school children and staff during earthquakes.

Alaska, California, Oregon, Washington

Airborne radiometric data map alteration of porphyry copper systems in the Elkhorn district, MT

The Elkhorn district in southwestern Montana is known for past production of high-grade base and precious metals from skarn and carbonate replacement deposits. The district also hosts a porphyry Cu-Mo deposit. A radiometric survey flown over the district and surrounds provides maps of surficial concentrations of potassium, thorium, and uranium. Handheld gamma ray spectrometer measurements on altered and unaltered rocks in the district and surrounding region provide a means to ground truth the airborne data. We demonstrate that the airborne and ground-based measurements show similar map patterns. We calculate ratios of the radioelement concentrations and present a potassium enhancement map that combines potassium and ratios of potassium to thorium and uranium. The results highlight the rocks previously mapped as having K-feldspar-biotite and quartz-sericite alteration assemblages in the Elkhorn district and map additional hydrothermal systems in the region, including known alteration in the Radersburg district. The data enhancement techniques can be used as a screening tool for mapping additional porphyry copper systems.

Montana

Evolution of permeability and strength recovery of shear fracture under hydrothermal conditions

Geothermal energy is a clean and renewable resource that depends on the ability to move water through hot rock. In many locations, the ability to move water through rock requires the presence of extensive natural or human-made systems of fractures. However, these fracture systems are influenced by a variety of complex processes that occur at the temperature and pressure conditions found in geothermal reservoirs. To make geothermal systems more efficient, it is important to understand how these fractures evolve over time in response to these processes. This study explored how fractures in rock change under high-temperature and pressure conditions, like those found in geothermal reservoirs, through laboratory experiments and numerical simulations. In general fractures, tend to close due to pressure from the surrounding rock preventing the flow of water. One set of laboratory experiments revealed that fractures in granitic rock can weaken over time at elevated temperatures. These weak fractures can then slip slowly potentially keeping the fracture open for a longer period. This study also showed that at the highest temperature examined (250 °C) large differences between the chemical composition of the rock and the fluid moving through it could improve water flow temporarily likely due to the water dissolving parts of the rock. Flow rate predicted by models shows overall good agreement with the flow rate measured in our experiments but with some long time-scale variations that may be due to mechanical wear of the fracture surface. A second set of experiments found that, over the temperature range examined, the closure of fractures was caused by brittle failure driven by high stresses at points where fracture surfaces were in contact. This result suggests that fracture closure in geothermal reservoirs could be counteracted by maintaining high pore pressures. Additionally, we also tested a new method using electrical resistance to measure fracture closure in controlled systems with high precision.

Final Technical Report

Public water-supply systems and water use in Tennessee, 1988

This report summarizes the results of a study conducted by the U.S. Geological Survey, in cooperation with the Tennessee Department of Environment and Conservation (TDEC), Division of Water Supply in 1988. Data gathered during an inventory by the TDEC were collated to determine water use, supply sources, population served, and design and storage capacities of the systems. The inventory was limited to systems that were active on June 30, 1988. Results of a survey of the systems conducted by the Tennessee Department of Health and Environment during 1988 were a primary source of data for this report. Data from computer and manual files maintained by the Tennessee Department of Health and Environment and the U.S. Geological Survey also were used. The Division of Water Supply, TDEC, surveyed 541 public water-supply systems. These systems served 81 percent of the population of the State, or 3.95 million people. The gross per capita use statewide for public-supplied water was 179 gallons per day. Total water withdrawals for public supply increased about 39 percent from 510 million gallons per day (Mgal/d) in 1980, to 708 Mgalld in 1988. During the same period, the population increased about 7 percent. Surface-water withdrawals accounted for 63 percent (446 Mgal/d) of the total water withdrawn in the State. All of these withdrawals occurred in the Tennessee (56 percent or 249 Mgal/d) and the Ohio (44 percent or 197 Mgalld) hydrologic regions. Ground water supplied 262 Mgal/d or 37 percent of the total water withdrawn by public-supply systems statewide. Of that amount, 79 percent, or 208 Mgalld, was used in western Tennessee.

Tennessee