Geology ReportsSearch

SEARCH · Geology Reports

Results for “Journal of Water Process Engineering”

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,711 records · Page 53Linked to original sources

Deep groundwater total dissolved solids mapping in the Dakota Group, Williston Basin, USA

Growing concern about the quantity of available freshwater around the world has led to interest in surveying groundwater total dissolved solids (TDS) below water well depths. Deep TDS has not been systematically mapped, and there is much to learn about the distribution and controls on deeper groundwater. In sedimentary basins across the United States, groundwater resources often overlie hydrocarbon resources, providing an opportunity to use borehole geophysical data collected for hydrocarbons to characterize groundwater and pore space resources. This study adapts a recently developed subsurface geostatistical and geophysical modeling approach to continuously map groundwater TDS, porosity, and temperature in the Dakota Group of the Williston Basin—an undercharacterized regional aquifer system overlying deeper hydrocarbon reservoirs. Groundwater TDS in the Dakota Group ranges from approximately 4800 to 26,900 mg/L. TDS patterns are stratified with higher TDS in the lower and upper Dakota Group, and relatively lower TDS in the middle Dakota Group. The lower TDS in the middle zone may represent a preferential regional flow path for lower-TDS meteoric recharge from the west. The alternating pattern of TDS may also be evidence of higher-TDS inflows into the Dakota Group from underlying and potentially from overlying aquifers. Porosity is lower near the center of the Williston Basin and tends to be higher to the east, which may be related to grain size distributions. The new regional TDS and porosity modeling serves as a quantitative reference for water users and provides supporting evidence for hypotheses on Dakota Group recharge.

Montana, North Dakota, South Dakota

Multi-scale geophysical imaging of a hydrothermal system in Yellowstone National Park, USA

Little is known about the local plumbing systems that fuel Yellowstone’s famous hot springs, geysers and mud pots. A multi-method, multi-scale geophysical investigation was carried out in the Obsidian Pool Thermal Area (OPTA) to: (i) delineate the lateral extent of the hydrothermal area and associated surface features; (ii) estimate the dimensions of the upflow zone and identify its main controlling structures; (iii) assess fluids circulation pathways from depth to surface. Ground and airborne geophysical data were acquired to connect local and regional scales, from shallow to large depths. Maps of surface electrical resistivity show a strong correlation with hydrothermal features. At in-termediate depths, electrical resistivity permits delineating the upper limit of the upflow zone, while Poisson’s ratio highlights differences in subsurface fluid content. Combining these results with surface observations and topographic information, we speculate that differential mixing of hydrothermal and fresh water could explain the wide diversity of features observed at OPTA. Low electrical resistivity observed at large depths also suggest that a vast upflow zone, controlled by rhyolite flows and conjugate faults, underlies the OPTA. We speculate that hydrothermal fluids rise along fractures and reach the surface in topographic lows to form hydrothermal features. Our results show that synoptic, multi-scale geophysical measurements provide a roadmap for understanding where and how geologic heterogeneity, topography, fluid-gas separation, and the mixing of thermal and meteoric waters conspire to produce the wide variety of Yellowstone’s renowned hydrothermal features.

Wyoming

How do hydrological variability and human activities control the spatiotemporal changes of riverine nitrogen export in the Upper Mississippi River Basin?

Excessive nitrogen export from agricultural watersheds remains a critical water quality challenge, with the Upper Mississippi River Basin (UMRB) significantly contributing to downstream eutrophication and hypoxia in the Gulf. This study investigates the spatiotemporal dynamics of riverine nitrate plus nitrite (NO 3 – + NO 2 – -N) export across the UMRB at high spatial resolution (12-digit Hydrologic Unit Codes or HUC12 subwatershed scale) during 2001–2020 and quantifies the effects of anthropogenic activities and hydrological variability on riverine NO 3 – + NO 2 – -N export changes in the region between 2001–2005 and 2016–2020. Our results revealed hotspots of substantial increases in NO 3 – + NO 2 – -N yields across the UMRB, with distinct regional patterns in driving factors. Over the entire UMRB, NO 3 – + NO 2 – -N yields increased by 9.7 kg/ha/yr on average from 2001–2005 to 2016–2020, with anthropogenic activities contributing 4.8 kg/ha/yr and hydrological variability contributing 4.9 kg/ha/yr. The northern and western UMRB had combined influences from both anthropogenic activities and hydrological variability, while the east-central regions had predominantly hydrologically driven changes. Agricultural sources, including fertilizer, manure, and biological nitrogen fixation, collectively contributed over 80% of NO 3 – + NO 2 – -N loading throughout the basin. This framework for disentangling human and hydrological impacts provides critical insights for developing effective and targeted watershed management strategies to reduce nutrient losses and improve water quality.

Upper Mississippi River Basin

Light absorbing particles deposited to snow cover across the Upper Colorado River Basin, Colorado Rocky Mountains, 2013-16: Interannual variations from multiple natural and anthropogenic sources

Atmospheric particulate matter (PM) as light-absorbing particles (LAPs) deposited to snow cover can result in early onset and rapid snow melting, challenging management of downstream water resources. We identified LAPs in 38 snow samples (water years 2013–2016) from the mountainous Upper Colorado River basin by comparing among laboratory-measured spectral reflectance, chemical, physical, and magnetic properties. Dust sample reflectance, averaged over the wavelength range of 0.35–2.50 μm, varied by a factor of 1.9 (range, 0.2300–0.4444) and was suppressed mainly by three components: (a) carbonaceous matter measured as total organic carbon (1.6–22.5 wt. %) including inferred black carbon, natural organic matter, and carbon-based synthetic, black road-tire-wear particles, (b) dark rock and mineral particles, indicated by amounts of magnetite (0.11–0.37 wt. %) as their proxy, and (c) ferric oxide minerals identified by reflectance spectroscopy and magnetic properties. Fundamental compositional differences were associated with different iron oxide groups defined by dominant hematite, goethite, or magnetite. These differences in iron oxide mineralogy are attributed to temporally varying source-area contributions implying strong interannual changes in regional source behavior, dust-storm frequency, and (or) transport tracks. Observations of dust-storm activity in the western U.S. and particle-size averages for all samples (median, 25 μm) indicated that regional dust from deserts dominated mineral-dust masses. Fugitive contaminants, nevertheless, contributed important amounts of LAPs from many types of anthropogenic sources.

Colorado

An early Holocene wet period in the southwestern United States

Multiple generations of spring-fed streams traversed ∼800 km 2 of the Las Vegas Valley in southern Nevada between ca. 10.9 ka and 8.5 ka, depositing an extensive tufa network. The scale of this network and diversity of tufa morphologies is novel in North America and offers an opportunity to obtain quantitative paleoclimate data for the region during the early Holocene. We determined isotopic compositions and estimated past temperatures using clumped isotope data from early Holocene tufa on the valley floor (698 m) as well as tufa forming today at higher elevation in the nearby Spring Mountains at Cold Creek Spring (1856 m). Modern and fossil tufa yielded comparably low δ 18 O values, implying that source waters for both were derived from high-elevation winter precipitation. Clumped isotope temperatures of modern tufa average 15.8 ± 2.5 °C, aligning with mean summer temperatures of the emergent spring water, and indicate equilibrium conditions of tufa formation. The early Holocene tufa yielded similar clumped isotope temperatures, averaging 15.2 ± 3.9 °C, meaning it precipitated at temperatures that occur at much higher elevations today. The Las Vegas tufa record, combined with nearby and temporally correlative paleospring and lacustrine records, suggest that cool/wet conditions prevailed throughout the Mojave Desert during the early Holocene. These records also demonstrate that spring ecosystems responded to millennial-scale hydroclimate variations that supersede climate change driven solely by insolation. The previously unrecognized pattern of ecosystem response to hydroclimate documented here may assist in understanding climate drivers for the early Holocene and provide critical information for the fate of groundwater-dependent ecosystems in the southwestern United States.

Nevada

Suspended-sediment characteristics in the Housatonic River Basin, western Massachusetts and parts of eastern New York and northwestern Connecticut, 1994-96

Suspended-sediment concentrations, discharges, loads, and yields were determined for eight subbasins in the Housatonic River Basin in western Massachusetts, eastern New York, and northwestern Connecticut from April 1994 through March 1996. Suspended-sediment samples were collected at three continuous-record sediment stations and at four partial-record sediment stations. Suspended-sediment concentrations in samples collected during the period of study ranged from less than 0.5 to 3,400 milligrams per liter, and concurrent streamflows ranged from 0.03 to 126 cubic feet per second per square mile at the seven stations. Median suspended-sediment concentrations in samples collected at each station ranged from 7 to 61 milligrams per liter. Median streamflows during suspended-sediment sampling ranged from 1.86 to 5.88 cubic feet per second per square mile. Instantaneous suspended-sediment yields ranged from less than 0.005 to 185 tons per day per square mile, and medians ranged from 0.03 to 1.12 tons per day per square mile at the seven stations. Total suspended-sediment loads (mass) from April 1994 through March 1996 at the continuous-record sediment stations were 11,603 tons at Housatonic River near Great Barrington, 7,929 tons at Green River, and 54,347 tons at Housatonic River near Ashley Falls. Suspended-sediment load during January 1996 at the Green River station accounted for about 54 percent of the total suspended-sediment load for the Green River during the 2-year study. Suspended-sediment load on January 19 and 20, 1996, at the Green River station accounted for about 50 percent of the January 1996 suspended-sediment load, or about 27 percent of the total suspended-sediment load during the 2-year study. This large suspended-sediment transport was the result of rainfall and snowmelt on January 19 and 20, 1996--the equivalent of a 5- to 6-inch rain storm in the Green River subbasin. Total suspended-sediment loads during the 2-year study at the partial-record sediment stations were 3,052 tons at Williams River, 1,758 tons at Ironworks Brook, and 17,927 tons at Konkapot River. Suspended-sediment yields from April 1994 through March 1996 at the continuous-record sediment stations were 21 (tons/yr)/mi 2 at Housatonic River near Great Barrington, 78 (tons/yr)/mi 2 at Green River, and 58 (tons/yr)/mi 2 at Housatonic River near Ashley Falls. Suspended-sediment yields during the 2- year study at the partial-record sediment stations were 35 (tons/yr)/mi 2 at Williams River, 78 (tons/yr)/mi 2 at Ironworks Brook, and 147 (tons/yr)/mi 2 at Konkapot River. Suspended-sediment yields were estimated for two subbasins in the Housatonic River Basin--Schenob Brook at Sheffield, a partial-record sediment station, and the area adjacent to the Housatonic River between Great Barrington and Ashley Falls. The estimate of suspended-sediment yield at Schenob Brook at Sheffield of 82 (tons/yr)/mi 2 is comparable to the yield determined for the Green River and Ironworks Brook. The estimate of suspended-sediment yield for the area adjacent to the Housatonic River between Great Barrington and Ashley Falls was 395 (tons/yr)/mi 2 . This estimated suspended-sediment yield was 2.7 to 18.8 times greater than that estimated for any of the other subbasins. Several basin and land-use characteristics thought to affect suspended-sediment transport in the subbasins were compared to the suspended-sediment yields. The characteristics that seemed to affect suspended-sediment discharge were dams, which contributes to decreased yields; Hadley, Limerick, Linlithgo, Saco, and Winooski silt loam soils (high erodibility soils), which contributes to increased yields; stratified-drift deposits, which contributes to increased yields; and agricultural and open land, which contributes to increased yields. The effect of stratified-drift deposits on suspended-sediment discharge is thought to be greater when those deposits are of glaciolacustrine (generally clay, silt, and fine sand), rather than glaciofluvial (clay, silt, sand, gravel, and cobbles) origin. The silt loam soils, glaciolacustrine deposits, and agricultural and open land are interrelated, inasmuch as the silt loam soils generally are associated with glaciolacustrine deposits and agricultural activities.

Connecticut, Massachusetts, New York

Land-based nutrient flux to a fringing reef: Insights from Ofu Island, American Samoa

Submarine groundwater discharge (SGD) is a critical driver of nutrient transport in coral reef ecosystems, shaping water quality, primary productivity, and overall reef health. This study quantifies SGD fluxes and associated nutrient dynamics in two reef flat pools within the Ofu Unit of the National Park of American Samoa: Papaloloa and Fatuana. A multi-method approach integrating unoccupied aerial system-based thermal infrared (UAS-TIR) surveys, radon-based SGD measurements, multichannel electrical resistivity tomography (ERT), and discrete water sampling was used to assess SGD rates and nutrient contributions. UAS-TIR imagery revealed cooler sea surface temperatures in both pools, indicative of SGD, with the higher fluxes observed in Papaloloa. Radon measurements revealed a strong inverse correlation between SGD rates and tidal stage, with a more immediate SGD response at Papaloloa due to its highly permeable calcareous sand and gravel substrate. In contrast, a 2–3-hour lag in SGD response at Fatuana suggests discharge from a more inland aquifer that has lower diffusivity. Nutrient concentrations correlated with temperature and salinity, confirming SGD as the dominant nutrient transport mechanism, whereas isotopic analyses indicated inputs from both groundwater and potential anthropogenic sources. Despite lower SGD flux at Fatuana, higher algal cover suggests additional factors influencing algal proliferation, including substrate availability and hydrodynamic conditions. Excess nutrient inputs from SGD may contribute to algal overgrowth, which threatens Ofu’s thermally tolerant corals by increasing competition for space and light. These findings underscore the complexity of SGD-mediated nutrient dynamics in reef environments and emphasize the need for integrated hydrological and ecological assessments to support effective reef conservation and management strategies.

American Samoa, Ofu Island, Olosega Island

U.S. Geological Survey Groundwater Climate Response Network, 2024

As of October 2024, the U.S. Geological Survey (USGS) operated 588 sites across the United States and its territories as part of the Groundwater Climate Response Network (CRN). The CRN is comprised of wells selected to monitor the effects of climate variability, such as droughts, on groundwater levels nationwide. The CRN includes nearly 500 locations with real-time data and more than 100 sites with non-real-time data available to the public on the CRN web mapper and the USGS National Water Dashboard.

General Information Product

Mapping the resistivity structure of Walker Ridge 313 in the Gulf of Mexico using the marine CSEM method

A marine controlled source electromagnetic (CSEM) campaign was carried out in the Gulf of Mexico to further develop marine electromagnetic techniques in order to aid the detection and mapping of gas hydrate deposits. Marine CSEM methods are used to obtain an electrical resistivity structure of the subsurface which can indicate the type of substance filling the pore space, such as gas hydrates which are more resistive. Results from the Walker Ridge 313 study (WR 313) are presented in this paper and compared with the Gulf of Mexico Gas Hydrate Joint Industry Project II (JIP2) logging while drilling (LWD) results and available seismic data. The hydrate, known to exist within sheeted sand deposits, is mapped as a resistive region in the two dimensional (2D) CSEM inversion models. This is consistent with the JIP2 LWD resistivity results. CSEM inversions that use seismic horizons provide more realistic results compared to the unconstrained inversions by providing sharp boundaries and architectural control on the location of the resistive and conductive regions in the CSEM model. The seismic horizons include: 1) the base of the gas hydrate stability zone (BGHSZ), 2) the top of salt, and 3) the top and bottom of a fine grained marine mud interval with near vertical hydrate filled fractures, to constrain the CSEM inversion model. The top of salt provides improved location for brines, water saturated salt, and resistive salt. Inversions of the CSEM data map the occurrence of a ‘halo’ of conductive brines above salt. The use of the BGHSZ as a constraint on the inversion helps distinguish between free gas and gas hydrate as well as gas hydrate and water saturated sediments.

Louisiana

Approximate inland extent of saltwater intrusion at the base of the Biscayne aquifer, Miami-Dade County, Florida, 2022

Miami-Dade County is part of a densely populated urban corridor in southeastern Florida. The Biscayne aquifer serves as Miami-Dade County’s primary drinking water source and is characterized by highly permeable karstic limestone and carbonate sand. The aquifer’s coastal location and permeable nature make it susceptible to saltwater intrusion. Monitoring the current inland extent and the rate of movement of the saltwater front in the aquifer can inform management strategies for conserving the long-term sustainability of the county’s water supply. In the 1950s, the U.S. Geological Survey published a map of the inland extent of saltwater intrusion in the Biscayne aquifer and has continued to update this map to monitor changes over time, with the most recent update published in 2018. An updated map has been created showing the approximate inland extent of saltwater intrusion in the Biscayne aquifer in eastern Miami-Dade County in 2022, with the 2018 extent shown for comparison. The inland extent of saltwater intrusion was mapped through the interpretation of borehole electromagnetic induction logs and measurements of chloride and specific conductance in groundwater samples. The location of the saltwater interface at the base of the Biscayne aquifer was represented by the 1,000-milligram-per-liter isochlor. This report describes changes in the location of the saltwater interface from 2018 to 2022. By 2022, the saltwater interface had moved farther inland in both the northern and southern parts of the county, advancing by as much as 0.3 kilometer in the north and up to 0.8 kilometer in the Model Land Area to the south. However, it remained relatively unchanged from its 2018 position in the east-central part of the county.

Florida

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

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

Nebraska

Geologic framework and hydrostratigraphy of the Edwards and Trinity aquifers within parts of Bandera and Kendall Counties, Texas

The karstic Edwards and Trinity aquifers are classified as major sources of water in south-central Texas by the Texas Water Development Board. During 2019–23 the U.S. Geological Survey, in cooperation with the Edwards Aquifer Authority, mapped and described the geology and hydrostratigraphy of the rocks composing the Edwards and Trinity aquifers within parts of Bandera and Kendall Counties from field observations of the surficial expressions of the rocks. The thicknesses of the mapped lithostratigraphic and hydrostratigraphic units were also estimated from field observations in the study area. The Cretaceous rocks in the study area are part of the Trinity Group and Edwards Group. The groups, formations, and members are composed primarily of layers of marls, shales, and limestones. The limestones are composed of mudstone through grainstone, framestone and boundstone, dolomite, and argillaceous and evaporitic rocks. The principal structural feature in the study area is the Balcones fault zone. The Balcones fault zone is the result of late Oligocene and early Miocene extensional faulting and fracturing that was a result of the eastern Edwards Plateau uplift. In the Balcones fault zone, most of the faults in the study area are high-angle to vertical, en echelon, normal faults that are predominantly downthrown to the southeast. Hydrostratigraphically, the rocks exposed in the study area are those that contain the Edwards aquifer, the upper zone of the Trinity aquifer, and the middle zone of the Trinity aquifer. Descriptions of the hydrostratigraphic units, thicknesses, hydrologic function, porosity types, and field identification and observations are provided, including those for the informal Bandera and Love Creek hydrostratigraphic units of the Edwards aquifer, which were identified through the mapping for this study.

Texas

U.S. Geological Survey Groundwater Climate Response Network—2023

As of October 2023, the U.S. Geological Survey (USGS) operated more than 660 sites across the United States and its territories as part of the Groundwater Climate Response Network (CRN). The CRN is comprised of wells and springs selected to monitor the effects of climate variability, such as droughts, on groundwater levels and spring discharge nationwide. The CRN includes more than 550 locations with realtime data and more than 100 sites with non-real-time data available to the public on the CRN web mapper and the USGS National Water Dashboard.

General Information Product

Utilization of multiple geochronology techniques to constrain the age of laterization and mineralization of the world-class Mount Weld rare earth element deposit, Western Australia

Pervasive chemical weathering on stable cratons may form thick regoliths and elemental enrichment, but constraining the age of regolith formation is challenging. In this study we utilize multiple geochronological techniques on different minerals from the world-class Mount Weld rare earth element (REE) deposit, formed by lateritic weathering of a carbonatite, to constrain the age of formation and provide insight into landscape evolution. The oldest dates, ca. 100 to 38 Ma, are from Lu-Hf dating of churchite [HREE(PO 4 )·2(H 2 O)], a heavy REE phosphate mineral. Growth bands on individual minerals show a younging outwards. 40 Ar/ 39 Ar geochronology of cryptomelane [K(Mn 4+ ,Mn 2+ )₈O₁₆] yielded dates from ca. 40 to 27 Ma. Similarly, (U-Th)/He geochronology of goethite [FeO(OH)] yielded dates ranging from ca. 45 to 19 Ma. Integrating results into regional constraints, suggests 1) churchite formed by mineral saturation in a karst-like setting below the water table from ca. 100 to 40 Ma, 2) with minor uplift and erosion, cryptomelane and goethite formed at or near the water table between ca. 45 and 19 Ma, 3) after ca. 15 to 10 Ma chemical weathering within the profile had ended. Other studies document that the region experienced minimal uplift and a wet, warm climate from ca. 100 Ma to 15 Ma. These conditions and the high carbonate content of the carbonatite promote extensive chemical weathering, a deep weathering profile, and the preservation of the weathered section. This study highlights the use of multiple geochronological techniques utilizing different minerals to provide insight into how laterites form and to constrain the timing and history of the formation of this important mineral deposit.

Western Australia

Introduction and methods of analysis for peak streamflow trends and their relation to changes in climate in Illinois, Iowa, Michigan, Minnesota, Missouri, Montana, North Dakota, South Dakota, and Wisconsin

Flood-frequency analysis, also called peak-flow frequency or flood-flow frequency analysis, is essential to water resources management applications including critical structure design and floodplain mapping. Federal guidelines for doing flood-frequency analyses are presented in a U.S. Geological Survey Techniques and Methods Report known as Bulletin 17C. A basic assumption within Bulletin 17C is that for drainage basins without major hydrologic alterations, statistical properties of the distribution of annual peak streamflows (peak flows) are stationary; that is, the mean, variance, and skew are constant. The stationarity assumption has been widely accepted within the flood-frequency community; however, a better understanding of long-term climatic persistence and concerns about potential climate change and land-use change has caused a reexamination of the stationarity assumption. Flood-frequency analyses that do not incorporate observed trends and abrupt changes may result in a poor representation of the true flood risk. Bulletin 17C does not offer guidance for incorporating nonstationarities when estimating floods, and it describes a need for studies that incorporate changing climate or basin characteristics. In response to this need and a history of concern regarding nonstationarity peak flows in the region, this study was done by the U.S. Geological Survey, in cooperation with the Departments of Transportation of Illinois, Iowa, Michigan, Minnesota, Missouri, South Dakota, and Wisconsin; the Montana Department of Natural Resources and Conservation; and the North Dakota Department of Water Resources, to assess potential nonstationarity in peak flows in the north-central United States. This chapter summarizes the methods used to detect hydroclimatic changes in peak-flow data in the study region. A wide range of analyses and statistical approaches are applied to document the primary mechanisms controlling floods and characterize temporal changes in hydroclimatic variables and peak flow. Four periods were selected for analysis of peak flow, daily streamflow, and climate data. The periods are (1) a 100-year period, 1921–2020; (2) a 75-year period, 1946–2020; (3) a 50-year period, 1971–2020; and (4) a 30-year period, 1991–2020. The climate data consist of monthly time series estimates of temperature, precipitation, potential evapotranspiration, actual evapotranspiration, snowfall, soil moisture storage, snow water equivalent, and runoff on a 3.1-mile by 3.1-mile grid for the conterminous United States. Statistical and graphical analyses were used to investigate potential changes in hydrology and climate. The starting point for these analyses was the initial data analysis of peak flow described in Bulletin 17C, which includes plotting the peak flow and checking for autocorrelation, monotonic trends, and changes points. Analyses were added to examine additional features in the data. To examine potential causal drivers of changes, the climate data were analyzed graphically and statistically. Results are provided in a U.S. Geological Survey data release. The study limitations are documented for users of the results.

Iowa, Illinois, Michigan, Missouri, Minnesota, Mon

Corwin Springs Known Geothermal Resources Area, Park County, Montana

The Corwin Springs Known Geothermal Resources Area (KGRA) is contiguous to Yellowstone National Park along a part of the northern boundary of the park near Gardiner, Park County, Mont. The area contains two known sites of hot-spring activity--LaDuke Spring, 2.8 km southeast of the small resort community of Corwin Springs, and Bear Creek Spring, 2.6 km east of Gardiner. LaDuke Spring issues from brecciated quartzite and has a flow rate of 380 1/min, a surface water temperature of 65°C, a silica-geothermometer temperature of 68.7°C, and a Na-K-Ca geothermometer temperature of 76.8°C. Bear Creek Spring issues from limestone and has a flow rate of 4 1/min, a surface water temperature of 32°C, a silica temperature of 46.7°C, and a Na-K-Ca temperature of 87.2°C. The springs, which are actively depositing travertine, are on or near the trace of the Gardiner fault-a high-angle reverse fault which forms the southwestern boundary of the Beartooth uplift. Features that make the Corwin Springs KGRA potentially significant include: (1) Proximity of the Corwin Springs area to significant Pliocene and Pleistocene volcanism in Yellowstone National Park and related geothermal activity, (2) existence of a large potential geothermal reservoir in the northern extension of the Sepulcher Mountain graben, (3) localization of thermal activity along the Gardiner fault in the Corwin Springs area, (4) occurrence of negative gravity anomaly centered over the northern end of the Sepulcher Mountain graben, (5) existence of negative magnetic anomaly in the Corwin Springs area, (6) location of the area within the Intermountain seismic belt, (7) evidence of recent tectonism in the area, as demonstrated by Pleistocene and Holocene faulting, and (8) observed surface temperatures and estimated geochemical temperatures of the two known hot springs (features which do not indicate high subsurface temperatures).

Montana

Ground deformation and gravity for volcano monitoring

Introduction When magma accumulates or migrates, it can cause pressurization and related ground deformation. Characterization of surface deformation provides important constraints on the potential for future volcanic activity, especially in combination with seismic activity, gas emissions, and other indicators. A wide variety of techniques and instrument types have been applied to the study of ground deformation at volcanoes (sidebar, p. 2; Dzurisin, 2000, 2003, 2007). Geodetic instruments include continuously recording Global Navigation Satellite System (GNSS; of which the United States’ Global Positioning System is one example) stations (fig. D1), borehole tiltmeters, and interferometric synthetic aperture radar (InSAR) measurements (from satellites, occupied and unoccupied aircraft systems, and ground-based sensors). Additional geodetic measurements like continuous- and survey-mode gravity (fig. D2) can contribute substantially to interpreting these data. Borehole strainmeters (see chapter K , this volume, by Hurwitz and Lowenstern, 2024) also have outstanding utility for monitoring deformation, although because of cost and permitting challenges, we do not include them as part of standard volcano monitoring networks for U.S. volcanoes. Still other techniques like light detection and ranging (lidar), structure from motion, and optical satellite data can be used to derive gross topographic changes, which can be used to map volcanic deposits, infer eruption rates, and gain insights into the source processes associated with eruptive activity (see chapter G , this volume, on tracking surface changes caused by volcanic activity; Orr and others, 2024). Experience has shown that no single geodetic monitoring technique is adequate to detect and track the entire range of ground-motion patterns that occur at volcanoes, primarily because of the temporal and spatial diversity of volcano deformation (fig. D3). Similarly, the magnitude of surface deformation varies widely. Geodetic monitoring strategies should therefore include multiple techniques and instrument types to cover a wide range of spatial and temporal scales. In identifying recommendations for geodetic instrumentation for volcano monitoring networks, we attempted to maximize the diversity of instrument types to measure the full range of deformation signals and minimize their expense and number; thus, we do not include several well-known deformation-monitoring techniques in our recommendations. Extensometers, for example, measure strains over distances of a few meters and have an excellent record of success in detecting changes in preeruptive localized ground motion across existing cracks, including at Mount St. Helens, Washington (Iwatsubo and others, 1992), and Piton de la Fournaise, Réunion Island (Peltier and others, 2006). Despite being relatively inexpensive, extensometers are best used primarily when localized ground displacements (for example, ground cracks) need to be tracked, and are not necessary at all volcanoes. In considering volcano deformation monitoring strategies, two complicating factors are deserving of special attention. First, not all deformation is driven by subsurface magmatic activity—for example, at many large stratovolcanoes (for example, Mount Rainier), flank collapses and landslides are significant geologic hazards (Reid and others, 2001) that may occur even in the absence of magmatic activity. Monitoring the stability of volcanoes is thus another critical application of geodetic monitoring networks to inform hazard assessment. One of the most famous examples of edifice instability is the large flank collapse that initiated the May 18, 1980, eruption of Mount St. Helens. Deformation monitoring had detected a bulge on the north flank of the mountain in April 1980 that was expanding by several meters per day (Lipman and others, 1981). Given that flank collapses can happen at any time during a period of volcanic unrest (or even outside a period of unrest), the capability to assess edifice stability is critical. Second, although volcanoes are commonly treated as idealized structures that erupt from single points, like centralvent stratovolcanoes, many are characterized by long rift zones from which eruptions may originate, and distributed volcanic fields are characterized by broadly spaced vents. For example, linear dikes are common at Kīlauea, Mauna Loa, and between Mount Shasta and Medicine Lake in California. At Kīlauea, one of these linear dikes emerged more than 40 kilometers (km) away from the summit of the volcano during the lower East Rift Zone eruption in 2018. Other volcanic fields, like Lassen volcanic center, California, or the San Francisco Volcanic Field, Arizona, have many small vents spread over a wide area. Although the instrumentation guidelines presented in this chapter remain phrased for central-vent volcanoes, they should be modified as needed in the context of the eruptive characteristics of each individual volcanic system. Spatial analysis of geodetic network coverage could help to ensure adequate instrumentation in areas where volcanism can occur over a broad area as opposed to a central vent. As an example, consider the adjacent volcanoes Mount Shasta and Medicine Lake. If station locations are chosen based only on the distance from the centers of the volcanoes, then any geodetic anomalies between the two volcanoes—an area of potential volcanism as indicated by the presence of volcanic features—may remain undetected by ground-based instrumentation. The spatial analysis is accomplished via a grid of pressure point sources (Mogi, 1958) evenly distributed across the map area, at a depth of 5 km in this example (fig. D4). Each source is inflated until predicted deformations exceed the GNSS white noise uncertainty estimates at one site (Langbein, 2017; Murray and Svarc, 2017). This volume of detectable magma provides a measure of the quality of the coverage (fig. D4). The results indicate that, as of 2022, there is a large area between Mount Shasta and Medicine Lake volcano with existing mapped dikes in which a substantial amount of magma could intrude without being detected geodetically. Applying this style of analysis to individual volcanic systems can provide a guide for designing network geometry given the expected locations of future eruptions.

Scientific Investigations Report

Iodine-129 in the eastern Snake River Plain aquifer at and near the Idaho National Laboratory, Idaho, 2021–22

Between the 1950s and 1980s, wastewater generated at the Idaho National Laboratory contained Iodine-129 ( 129 I); this wastewater was discharged directly into the eastern Snake River Plain (ESRP) aquifer through a deep disposal well, unlined infiltration ponds, or leaked from distribution systems below industrial facilities. During 2021–22, the U.S. Geological Survey, in cooperation with the U.S. Department of Energy and the Idaho Department of Environmental Quality Idaho National Laboratory Oversight Program, collected groundwater samples from 64 monitoring wells in the ESRP aquifer, 6 of which are part of a multilevel monitoring system, to determine the concentration of 129 I in the groundwater. These samples were analyzed by accelerator mass spectrometry as part of a long-term ongoing study to track trends and occurrences of this carcinogenic, long-lived radionuclide in the environment. Concentrations ranged from slightly above the locally determined background concentration of 5.4×10 −6 picocuries per liter, to just below the U.S. Environmental Protection Agency’s maximum contaminant level of 1 picocurie per liter. Discharge of wastewater containing 129 I has been discontinued to the aquifer, and long-term trends from a subset (n=15) of sampled wells show decreasing 129 I concentrations over the last three decades. Concentrations of 129 I in groundwater from monitoring wells near facilities at the Idaho National Laboratory are affected by episodic recharge from an ephemeral surface-water source and by the fracture-flow dominated hydrologic regime in the ESRP aquifer. The spatially focused sampling effort has also identified a low-level 129 I plume that affects long-term water quality near and downgradient from the Advanced Test Reactor Complex in the southwestern part of the facility that had not been clearly defined in previous sampling efforts, although the definition of the plume is somewhat limited by available data.

Idaho