Geology Reports⌕ Search

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 14Linked to original sources
USGS2025

Widespread occurrence of former anhydrite phenocrysts in Laramide-age magmas related to porphyry-skarn Cu mineralization at Santa Rita and Hanover-Fierro, New Mexico, USA

Reports of magmatic anhydrite are relatively rare, with only ~30 occurrences documented worldwide so far. However, magmatic anhydrite saturation is difficult to recognize because anhydrite decomposes rapidly in near-surface environments. In most cases, only anhydrite inclusions shielded within other phenocryst phases were able to survive. Alternatively, since anhydrite phenocrysts preserved in fresh volcanic rocks are characteristically intergrown with apatite phenocrysts, the former presence of anhydrite phenocrysts can be recognized based on the occurrence of lath-shaped cavities that show a strong spatial association with apatite phenocrysts. These cavities can be either empty or filled with low-temperature, secondary minerals such as zeolites, carbonates, or microcrystalline silica. A systematic search for the occurrence of such cavities, combined with optical and Raman-spectroscopic identification of anhydrite inclusions preserved within apatite, hornblende and quartz phenocrysts, demonstrates that most of the Laramide-age magmas associated with the Santa Rita and Hanover-Fierro porphyry-skarn Cu (Zn, Mo, Au, Pb) deposits were saturated in magmatic anhydrite. The anhydrite typically coexisted with monosulfide solid solution (MSS), suggesting oxygen fugacities of ~2.0±0.5 log units above the fayalite-magnetite-quartz buffer. The magmas range from andesitic to rhyodacitic in composition, and from shortly pre-mineralization (~61 Ma) to shortly post-mineralization (~57 Ma) in age. In three samples with particularly well-recognizable former anhydrite phenocrysts, their modal abundance could be quantified based on high-resolution scans of polished hand specimens. The observed modal anhydrite abundances of 0.63–1.8 vol% translate into minimum magma sulfur contents of 0.20–0.56 wt% S. The highest sulfur content of 0.56 wt% S is difficult to reconcile with available anhydrite solubility models, but it could be reproduced in an anhydrite solubility experiment performed at 950 °C and 1.15 GPa on a natural latite containing 13.1 wt% dissolved H2O. The sample with the second-highest sulfur content of 0.26 wt% S requires ~10 wt% H2O in the silicate melt, and, consequently, a minimum pressure of ~0.5 GPa. Taken together, the results suggest that the magmas of the Central Mining District were extremely hydrous and thus originated from great depth. Indeed, their major element compositions and reconstructed H2O and S contents agree well with experimentally observed and numerically predicted compositions of residual silicate melts after 50–70 wt% crystallization of ordinary arc basalts at high pressure and high oxygen fugacities.

New Mexico↗
USGS2026

Altitude of the potentiometric surface and depth to water in the Mississippi River Valley alluvial aquifer, spring 2024

Potentiometric-surface and depth-to-water maps for spring 2024 were created for the Mississippi River Valley alluvial aquifer using groundwater-altitude data from 1,151 wells completed in the Mississippi River Valley alluvial aquifer and from the altitude of the top of the water surface in area rivers from 160 U.S. Geological Survey and U.S. Army Corps of Engineers streamgages. The potentiometric-surface and depth-to-water maps for 2024 were created to support investigations to characterize the Mississippi River Valley alluvial aquifer as part of the U.S. Geological Survey Water Availability and Use Science Program. Sufficient data were available to map the potentiometric surface and depth to water of the Mississippi River Valley alluvial aquifer for spring 2024 for about 83 percent of the aquifer area. The potentiometric contours ranged from 0 to 330 feet (ft) above the North American Vertical Datum of 1988. The regional direction of the groundwater gradient was generally to the south-southwest, except in areas of groundwater-altitude depressions, where the groundwater gradient direction was into the depression, and near rivers, where the groundwater gradient direction was either from the aquifer to the river or from the river into the aquifer. Groundwater depressions in the potentiometric-surface map are in the lower one-half of the Cache region and in most of the Grand Prairie and Delta regions. Depth to water by well in the Mississippi River Valley alluvial aquifer in spring 2024 ranged from 0.27 ft above land surface to 145.65 ft below land surface.

Arkansas, Illinois, Kentucky, Louisiana, Mississip↗
USGS2024

Representation of surface-water flows using Gradient-Related Discharge in an Everglades Network

The Everglades Depth Estimation Network interpolates water-level gage data to produce daily water-level elevations for the Everglades in south Florida. These elevations were used to estimate flow vectors (gradients and directions) and volumetric flow rates using the Gradient-Related Discharge in an Everglades Network (GARDEN) application developed by the U.S. Geological Survey in cooperation with the U.S. Army Corps of Engineers. Flow rates in both the east-west and north-south directions were computed on a 400-meter square grid using modified parameters in the Manning’s equation. The frictional resistance parameter in the Manning’s equation was calibrated to measured flow rates at coastal creeks fed by Everglades Depth Estimation Network boundary flows. Levees and other features that act as barriers to flow were defined as “no-flow” grid cells where vectors were set to zero. The flow volume magnitudes were calibrated with 2020 daily values of coastal river flows, and verification was performed using 2021 data. Within a given day, the measured coastal river flows fluctuate more than the GARDEN boundary flows because of tidal and wind forcings. Because the GARDEN boundary flows were the upstream water source for the coastal rivers, calibration focused on matching average daily flow volumes rather than daily fluctuations. The Pearson’s correlation coefficient is 0.766 for the 2020 calibration period and 0.566 for the 2021 verification period. Applying GARDEN to periods with hydraulic-control-structure releases allows the propagation of structure flows to be seen in the daily flow-vector maps along with the multiday response of flows farther downgradient. Flow vectors may be overestimated near control structures because of difficulties in resolving the water gradient downstream from the structure. Flow vectors farther from the structure are more accurate than those near the structure.

Florida↗
USGS2025

Fiber-optic distributed temperature sensing of hydrologic processes—Diverse deployments and new applications by the U.S. Geological Survey

Fiber-optic distributed temperature sensing instruments harness the temperature-dependent properties of glass to measure temperature continuously along optical fibers by using precise pulses of laser light. In the mid-2000s, this technology was refined for environmental monitoring purposes such as snowpack-air exchange, groundwater/surface-water exchange, and lake-water stratification. Fiber-optic distributed temperature sensing has revealed unprecedented details about preferential flow processes; however, the method is labor intensive and requires specific training, resulting in limited use by the broader water community. With the ongoing national implementation of the U.S. Geological Survey Next Generation Water Observing System, there has been renewed interest in harnessing the unique spatiotemporal monitoring capabilities of fiber-optic distributed temperature sensing. This fact sheet briefly describes this technology, highlights uses by the U.S. Geological Survey, and discusses current applications and future opportunities.

Fact Sheet↗
USGS2026

Effects of a low-lipid diet on the gut microbiome and head kidney transcriptome of juvenile Chinook Salmon

Objective Pacific salmon Oncorhynchus spp. reared in production hatcheries are typically fed high-lipid, energy-dense diets to achieve large size and high body condition prior to release. In contrast, juveniles in natural environments tend to consume low-lipid, high-protein diets, and fish reared for research or conservation purposes are sometimes fed diets that are formulated to mimic natural diets and promote wild-like phenotypes. Understanding how these alternative diets affect fish health beyond growth and body condition could ultimately contribute to improving hatchery fish fitness. Methods In this work, we evaluated changes in the fecal microbiome and gene expression of juvenile Chinook Salmon O. tshawytscha on a standard high-lipid hatchery diet versus a low-lipid diet formulated to mimic the nutrition profile of natural-origin fish. To evaluate the time scale at which diet alters the fecal microbiome, we collected longitudinal samples over a 12-week period and switched the diets of a subset of fish twice during the experiment. We used 16S ribosomal RNA gene amplicon sequencing to characterize fecal microbiome differences between fish on the two diets as well as hatchery-reared fish at a production hatchery, hatchery fish that had been captured after release into a stream, and natural-origin, stream-reared fish of similar ages. Additionally, we conducted RNA sequencing on head kidney samples from laboratory-reared fish to evaluate changes in gene expression in this important immune organ. Results We found that the low-lipid diet and the hatchery diet resulted in microbiomes that differed from the microbiome of natural-origin fish and from each other and that diet-driven changes to the microbiome could occur in under 14 d. The low-lipid diet did not result in a microbiome that resembled the microbiome of naturally produced fish. Instead, the low-lipid diet resulted in a microbiome community that was distinct from those of fish reared on the hatchery diet and fish sampled from the wild. The RNA sequencing results indicated differential enrichment of pathways related to immunity, metabolism, and hormone synthesis between fish that were fed the two experimental diets. Conclusions The results suggest that additional environmental factors influence the microbiome more strongly than diet formulation or that the low-lipid diet has a smaller effect on the microbiome than a natural, ­invertebrate-based diet. Given that the gut microbiome and systemic immune function contribute significantly to disease resistance, our findings highlight the importance of understanding how diets fed to fish in captivity may affect fish health beyond growth and body condition metrics.

Journal of Aquatic Animal Health↗
USGS2026

A review and synthesis of post-wildfire shifts in hydrologic processes and streamflow generation mechanisms

Critical water supply watersheds in the western United States (WUS) are impacted by wildfires, with potential negative effects on water quality and quantity. Scientific understanding is currently insufficient to deliver estimates of wildfire consequences for water quantity that are regionally accurate. Regional variability in the directionality and magnitude of post-wildfire shifts in streamflow generation fuels uncertainty in estimates of wildfire effects on water supply. In this work we provide a narrative review of wildfire effects on hydrologic processes and the resulting changes in streamflow generation mechanisms with a focus on the WUS, incorporating other global regions when pertinent. A conceptual model summary of wildfire effects on streamflow generation emphasizes: (1) precipitation seasonality, (2) synchrony of precipitation and potential evapotranspiration, (3) net shifts in interception, evaporation, and transpiration relative to total annual precipitation, (4) vegetation changes, including compensatory uptake and type conversion, (5) degree of overlap in rainfall rates and infiltration, (6) fire extent and severity, (7) burn scar positioning (e.g. in headwaters or proximal to watershed outlet), (8) scale-dependent groundwater leakage, (9) near-surface water storage reduction, and (10) soil to groundwater connectivity. Ongoing gaps and challenges include separating the influences of precipitation variability, water withdrawals, and post-fire land management; compound and overlapping disturbances; and lack of pre-fire data. Notable future opportunities include: harnessing ever-improving gridded and remotely sensed precipitation and fire-effects data; linking geophysical, isotopic tracer, and geochemical signatures to diagnose hydrologic changes; leveraging physically based and data-driven model advancements; and analyzing streamflow generation recovery trajectories across diverse watersheds.

western United States↗
USGS2026

Decoding the oxidative digestion mechanism for polystyrene nanoplastic detection in the Great Lakes using a customizable Raman spectral processing algorithm

Despite the concerns surging around nanoplastics (NPs) regarding their prevalence and bioavailability in freshwater systems, robust detection of NPs in complex environmental matrices is hindered by the lack of standardized sample pretreatment and a mechanistic understanding of oxidative digestion. Here, we systematically investigate the interaction between hydrogen peroxide (H 2 O 2 ) and polystyrene (PS) NPs during digestion in deionized (DI) water and four environmental matrices from in and around the Great Lakes Basin. To facilitate high-throughput analysis, we develop Pre_peak, a customizable Raman spectral processing algorithm that achieves >99% accuracy for both NP identification and interference rejection, allowing reliable NP quantification via pixel counting and systematic decoding of the oxidative digestion mechanisms. In DI water, varying H 2 O 2 doses from 0 to 30% has negligible effects on the recovery and Raman signal intensity of PS NPs over 24 hours of digestion. However, morphological changes and aggregation of PS NPs are observed when the H 2 O 2 dose exceeds 20%. Prolonged digestion further leads to progressive NP loss. In natural waters, the optimal dosage and digestion duration depend on matrix characteristics, including dissolved organic matter (DOM) and ion composition. This study provides mechanistic insights into NP–oxidant interactions and underscores the need for matrix-tailored digestion protocols to advance standardized NP detection in freshwater environments.

Great Lakes↗
USGS2025

Detecting hidden sedimentary geothermal systems in the Upper Colorado River Basin

Geothermal resources exist in sedimentary rock where circulation of water for efficient extraction or storage of heat is possible. Except in rare instances where hot water is expressed at the land surface, sedimentary geothermal resources are hidden, so the identification of these systems is optimally accomplished using predictive subsurface modeling. An integrated approach using detailed paleogeographic interpretations, subsurface geologic mapping, and numerical modeling has produced regional geologic and temperature models for the Upper Colorado River Basin, a large watershed in central North America that contains many sedimentary basins. These models identify areas of hidden sedimentary geothermal resource potential in low temperature (<90°C), moderate temperature (90–150°C), and high temperature (>150°C) fairways across the study area. These models incorporate maps of key horizons in outcrop and the subsurface to create a robust structural framework that can be used to target favorable geology for natural or engineered permeability. This framework is populated with lithologies derived from detailed palaeogeographical maps and over 40,000 bottom hole temperature (BHT) values were used to create a calibrated three-dimensional (3D) temperature model across the region. The resulting maps serve as a regional sedimentary geothermal play fairway screening tool for evaluating different grades of sedimentary geothermal resources and for identifying areas of interest where more detailed, prospect-scale studies can be undertaken.

Arizona, Colorado, New Mexico, Utah, Wyoming↗
USGS2025

Natural source zone depletion of crude oil in the subsurface: Processes controlling mass losses of individual compounds

At many petroleum hydrocarbon spill sites, residual spilled product forms a long-term source of groundwater contamination. The phrase source zone natural depletion is used to refer to the mass loss rates. Overall mass lost under environmental conditions was analyzed using conservative biomarker concentrations for a 1979 oil spill in northern Minnesota, USA. After 40–41 years, an average of 50% of the mass was lost with values ranging from 22% to 57% depending on location. It is also important to understand the composition changes in the source. To understand controls on the losses of individual compounds, concentrations of volatile hydrocarbons in oil samples were compared with aqueous solubilities, and pore-space oil saturations. The results of the comparison show that losses of the oil compounds were controlled by pore-space oil saturations, solubility, and susceptibility to degradation under methanogenic conditions. Compounds that degrade under methanogenic conditions, including toluene, o -xylene, and n -alkanes are more depleted compared to benzene, ethylbenzene, and m - and p -xylene for which losses are dominated by dissolution. These rates and compound-specific behaviors form a foundation for improved modeling approaches and risk analyses.

Minnesota↗
USGS2024

Hyperspectral (VNIR-SWIR) analysis of roll front uranium host rocks and industrial minerals from Karnes and Live Oak Counties, Texas Coastal Plain

VNIR-SWIR (400–2500 nm) reflectance measurements were made on the surfaces of various cores, cuttings and sample splits of sedimentary rocks from the Tertiary Jackson Group, and Catahoula, Oakville and Goliad Formations. These rocks vary in composition and texture from mudstone and claystone to sandstone and are known host rocks for roll front uranium occurrences in Karnes and Live Oak Counties, Texas. Spectral reflectance profiles, 569 in total, were reduced to 125 representative spectral signatures, which were analyzed using the U.S. Geological Survey's (USGS) Material Identification and Characterization Algorithm (MICA). MICA uses an automated continuum-removal procedure together with a least-squares linear regression to determine the fit of observed sample spectral absorption features to those of reference mineral standards in a spectral library. The reference minerals include various clay, mica, carbonate, ferric and ferrous iron minerals and their mixtures. In addition, absorption feature band-depth analysis was done to identify rock surfaces exhibiting absorption features related to uranium and zeolite minerals, which were not included in the command files used to execute MICA. Rocks from each of the four geologic units produced broadly similar spectral signatures as a result of comparable mineral compositions, but there were some notable differences. For example, Ca- and Na-montmorillonite was matched most frequently to the spectral absorption features in 2-μm (∼2000–2500 nm) wavelengths, while goethite occurred often at 1-μm (∼400–1000 nm) wavelengths. The latter is related to limonitic iron-staining in and around oxidized zones of the uranium roll front as described in previous papers. Rocks of the Jackson Group differed from those of the Catahoula, Oakville and Goliad units in that the former exhibited spectral features we interpret as being due to the presence of lignite-bearing mudstone layers. Goliad rocks exhibit spectral features related to dolomite, gypsum, anhydrite, and an unidentified green clay mineral that is possibly glauconite . Jackson Group rocks also exhibit weak but well-resolved absorption features at 964 and 1157 nm related to either or both zeolite minerals clinoptilolite and heulandite. These zeolite minerals and a few spectra exhibiting hydrous silica absorption features are indicative of alteration of volcanic glass in tuffaceous mudstone and claystone layers. A few sample spectra exhibited strong absorption features at around 1135 nm related to the uranium mineral coffinite. Both the 1135 nm coffinite and 1157 nm zeolite absorption features overlap somewhat, potentially making them difficult to distinguish without additional hyperspectral field, laboratory or remote sensing data. The results of this study were compared to mixtures of minerals described for ore, gangue and alteration minerals in deposit models for sandstone-hosted uranium, sedimentary bentonite and sedimentary zeolite. Use of these spectra can help facilitate mapping of both waste materials from the legacy mining of the above commodities, as well as future exploration and resource assessment activities.

Texas↗
USGS2024

Assessment of fresh groundwater discharge and saline surface-water intrusion at Operable Unit 2, North Chevalier Field Disposal Area (Site 11), Naval Air Station Pensacola, Florida, 2018–22

Site 11 is a former landfill at North Chevalier Field Disposal Area in Operable Unit 2 at Naval Air Station Pensacola, in northwest Florida. Site 11 is adjacent to Bayou Grande, a shallow, tidally influenced, saline estuary of the Pensacola Bay watershed. Federal and Florida regulators have expressed concern that contaminants detected in groundwater beneath the inland parts of Site 11 may discharge to Bayou Grande. In 2017, the Department of Defense, U.S. Navy, Naval Facilities Engineering Systems Command Southeast asked the U.S. Geological Survey to assess the occurrence of fresh groundwater discharge to Bayou Grande at Site 11 and to delineate to the extent practicable the location of groundwater discharge. Between 2018 and 2022, the U.S. Geological Survey used a multiple-lines-of-evidence approach that included a visual method and three physical methods based on the temperature difference between groundwater and surface water to assess groundwater discharge. One of the physically based methods also used the difference in specific conductance between fresh groundwater and brackish to saline surface water. Combined, the data indicate that fresh groundwater from across Site 11 discharges primarily along the shoreline of the northern and northeastern part of Site 11. The data also indicate that saline surface water from Bayou Grande intrudes tens of feet into the shallow aquifer beneath Site 11. The combined data indicate that the interface between fresh groundwater and saline surface water changes over space and time. Any new monitoring wells proposed for installation near the shoreline of Site 11 should include approaches to monitor the changes in the location of the freshwater/saltwater interface. Care would need to be taken to collect any groundwater samples at the correct season and tidal period to provide the highest probability of collecting a representative sample of Site 11 groundwater unaffected by saltwater.

Florida↗
USGS2025

Rapid emplacement of the Keaiwa Lava Flow of 1823 from the Great Crack in the Southwest Rift Zone of Kilauea volcano

The Keaīwa Lava Flow of 1823 in the Southwest Rift Zone of Kīlauea volcano is unusual for its expansive pāhoehoe sheet flow morphology and lack of constructive vent topography, despite having a similar tholeiitic basalt composition to other lavas erupted from Kīlauea. This lava flow issued from a ∼10-km-long continuous fissure now known as the Great Crack, and has an unusually thin sheet flow morphology with margin thicknesses of ∼15–110 cm (average of 42 cm). Based on field observations of the lava flow at its fissure vent (e.g., drain-back features), we propose that the Great Crack formed, or at least significantly widened, just prior to and syn-eruptively with this 1823 eruption. The absence of pyroclastic cones or spatter ramparts indicates that the eruption consisted of a rapid outpouring of relatively degassed lava as the fissure unzipped. The rapidly moving lava flow overtopped pre-existing tumuli and scoria cones (e.g., Lava Plastered Cones) up to ∼10 m tall. Glass and whole-rock chemistry yield homogeneous compositions for the lavas erupted from the Great Crack, with glass compositions of 6.40 ± 0.10 wt% MgO and whole-rock compositions of 7.39 ± 0.07 wt% MgO. Lava pads erupted from a short western fissure system are richer in mafic minerals (e.g., olivine and clinopyroxene), and show slightly more MgO-rich whole-rock compositions (7.79 ± 0.05 wt%). MgO-in-glass thermometry on juvenile spatter yield eruption temperatures of 1153 ± 13°C that are typical of Kīlauea lavas. Thus, the extensive sheet-like lava flow morphology is not a direct consequence of unusual magmatic or rheological conditions (i.e., low viscosity). Instead, the flow morphology is associated with high effusion rates caused by sudden drainage of uprift magma as it erupted from the Great Crack. Lava flow modeling on a 2-m-resolution digital elevation model indicates that a minimum bulk effusion rate of ∼5800 m 3 /s (∼3500 m 3 /s dense rock equivalent) and a minimum flow velocity of ∼11 m/s are required for the lava flow to overcome the topography of the Lava Plastered Cones. This effusion rate is among the highest inferred for eruptions in Hawaiʻi and around the world. This study highlights a less frequent eruption style at Hawaiian volcanoes characterized by a sudden outpouring of lava from an unusual fissure system. Local eyewitness accounts indicate that the 1823 eruption was preceded by seismicity. Given the complex magmatic-volcanic-tectonic relations across Kīlauea, we speculate that the south flank could have slipped over one or more events that ultimately triggered unzipping of the Great Crack and passive release of briefly stored uprift magma. An eruption similar to 1823 at Kīlauea or Mauna Loa, with an eruptive timeframe that could be as short as an hour, with high effusion rates and rapid flow front velocities, would not easily allow for a timely response.

Hawaii↗
USGS2026

Continental-scale prediction of hydrologic signatures and processes

Understanding how dominant hydrologic processes and their drivers vary across diverse continental-scale landscapes is critical for hydrologic modeling and water management applications. Our research addresses this question by synthesizing large-sample watershed datasets, Caravan and GAGES-II, and developing random forest models to identify patterns in hydrologic function. We assessed dominant processes by examining hydrologic signatures – summary indicators of watershed function derived from hydroclimatic time series and random forest models across 14 146 gauged United States watersheds. The results reveal clear continental-scale gradients in hydrologic processes, including baseflow, overland flow, storage, and water balance losses. Our map of dominant processes highlights, for example, the transition from baseflow to fast responses and back to baseflow along the elevation gradient from the Appalachian spine, through the Piedmont, to the Eastern Coastal Plain; a distinct outer ring around the Great Lakes region; and sharp contrasts between coastal and inland processes in the West. Variable importance analysis from random forest models show that processes in the western U.S. are primarily controlled by climate, whereas in the eastern U.S., soil, geology, and topography play larger roles, with distinct human influences apparent in urban areas. Our approach of estimating dominant processes and their drivers facilitates extending process knowledge from research watersheds to the continental scale, assessing current hydrological understanding, and evaluating hydrological model structures.

conterminous United States↗
USGS2026

Recent scientific contributions by the U.S. Geological Survey in the San Francisco Bay and Sacramento–San Joaquin Delta Estuary

Introduction The San Francisco Bay and Sacramento–San Joaquin Delta Estuary (hereafter, Bay-Delta) is the largest estuary on the West Coast of the United States. The Bay-Delta covers more than 1,600 square miles and drains a watershed of more than 75,000 square miles, which is greater than 40 percent of California. The region surrounding the Bay-Delta is home to about 10 million people, and its habitats (fig. 1) support more than 800 plant and animal species. The waterways of the Bay-Delta are the central hub of California’s extensive freshwater delivery system, supplying water to more than 27 million Californians and 4 million acres of farmland in the Central Valley. The U.S. Geological Survey (USGS) is the primary science agency of the U.S. Department of the Interior, providing a broad range of Earth, water, biological, and mapping data and expertise to inform natural resource management across the country. This fact sheet focuses on research conducted by the USGS in the Bay-Delta region, mostly within the past 5 years. The fact sheet is organized across five major themes: water flow and water quality, fish and wildlife, wetland restoration, invasive species, and hazards. In each of these areas, the USGS works closely with Federal, State, and local agencies; academia; and non-governmental organizations to provide objective, evidence-based science. The data and knowledge gained from USGS research and monitoring are publicly available and are routinely used by agencies, including the U.S. Bureau of Reclamation, U.S. Army Corps of Engineers, U.S. Fish and Wildlife Service, California Department of Water Resources, California Department of Fish and Wildlife, California State Water Resources Control Board, Delta Stewardship Council, and many additional organizations and stakeholders to inform management of the Bay-Delta region and its resources.

California↗
USGS2025

Toward a new framework to evaluate process-based model configurations and quantify data worth prior to calibration

Model criticism, discrimination, and selection methods often rely on calibrated model outputs. Because calibration can be computationally expensive, model criticism can first be undertaken by assessing model outputs obtained from limited prior parameter ensembles. However, such prior-based methods are often heuristic and do not formalize the notion of balancing model consistency with data and model complexity (i.e., model adequacy). We present a new framework to discriminate among candidate models prior to calibration that formalizes prior-to-calibration model adequacy into a metric to implicitly balance prior model output data coverage with model complexity represented by prior output (co)variance. The prior model adequacy metric “Mahalanobis distance deviation” quantifies the deviation of (a) the set of squared Mahalanobis distances of data from a prior model output distribution from (b) the set of squared Mahalanobis distances of data from their own distribution. A new data worth metric “discernment value” is also presented which quantifies the value of data for screening less-adequate models prior to calibration. Discernment value is calculated from the change in variance of a weighted average of prior model outputs from all candidate models due to less-adequate model outputs receiving lower weight. The framework is demonstrated using a one-dimensional groundwater flow model with eight possible configurations. A synthetic data network is used to test the framework. Results show the framework identifies the candidate models most similar to the true model used to create the synthetic data. Discernment values show variation in the value of different data types and locations for screening less-adequate models.

Water Resources Research↗
USGS2025

Emerging investigator series: Post-wildfire sediment geochemical characterization reveals manganese reactivity and a potential link to water quality impairment in the Gallinas Creek watershed, New Mexico

Water quality post-wildfire is often impaired by increased turbidity and elevated concentrations of elements such as manganese (Mn) and iron (Fe). Precipitation events exacerbate these issues, due in part to increased erosion and transport of sediment from hillslopes to surface water. Both Mn and Fe are major redox-active elements in sediments that drive a variety of biogeochemical cycles, precipitate adsorptive phases, and can themselves be drinking water contaminants. By investigating Mn and Fe sediment geochemistry in post-wildfire sediment deposits, related water quality hazards can be assessed. To establish and strengthen this connection, we analyzed the geochemistry of sediment deposits and surface water in the Gallinas Creek watershed, New Mexico over 1.5 years post-wildfire. Analyses included particle size analysis, water extractions, sequential extractions and aqua regia extractions to determine metal partitioning in sediment deposits. Data demonstrate Mn concentrations were distributed across labile and reactive fractions, such as the exchangeable and oxyhydroxide fractions, while Fe concentrations were mainly associated with the residual fraction. Manganese concentrations in aqua regia extractions and several fractions of sequential extractions were also strongly and significantly correlated with fine-grained sediment while the same pools of Fe concentrations were not. Dissolved Mn concentrations in surface water were elevated (>50 μg L −1 ) multiple times over the 1.5 years post-wildfire, highlighting a relationship between sediment geochemistry and water quality. This work shows Mn in sediments mobilized post-wildfire has an influence on water quality and highlights how further investigation into Mn sediment redox processes and mineralogy post-wildfire can inform risk assessments and resource management.

New Mexico↗
USGS2026

Uncertainty quantification of geologic energy storage in depleted gas reservoirs using material balance equations embedded in a hierarchical errors-in-variables model

The storage potential for gas in geologic settings, such as depleted hydrocarbon reservoirs and solution-mined salt caverns, is becoming salient to future energy infrastructure planning. Technologies such as carbon capture, utilization, and storage, carbon dioxide-enhanced oil recovery, and natural gas and hydrogen storage help meet growing energy demands, reduce carbon emissions to meet climate goals, and provide energy security amid geopolitical uncertainties. Therefore, estimates of underground gas storage capacity could be useful for efficiently navigating the energy transitions. Material balance is a fundamental method in reservoir engineering for estimating original gas in place and potential storage capacity at the scale necessary for national assessments of subsurface pore space resources. However, the deterministic method cannot accommodate multiple data sources or quantify uncertainty in predictions. In this study, a method that embeds material balance equations within a hierarchical errors-in-variables model is proposed which allows the estimation of the distributions of reservoir properties needed for assessments. Uncertainties associated with these reservoir properties have traditionally been expert-elicited, whereas the uncertainty estimates from the proposed models are data-driven. Capacity and uncertainty estimates can be used in a probabilistic resource assessment, supplementing information traditionally used by assessors or even replacing this expert elicitation step when data are unavailable. Various regression models are compared in a case study of the Michigan Basin, a large contributor to the United States’ current natural gas storage capacity. In particular, errors-in-variables models help ameliorate regression dilution and can quantify uncertainty in predictions of pressure in addition to storage capacity. Overfitting is addressed by quantifying generalization error and model averaging in simple and stratified cross-validation against reported working gas capacity, representing the varying quality and quantity of available data. Incorporating a statistical framework into existing numerical methods in reservoir engineering can improve the quality of estimation, and in particular, this method brings rigor to uncertainty quantification as part of a larger effort by the U.S. Geological Survey to assess domestic energy gas storage resources in depleted hydrocarbon reservoirs.

Mathematical Geosciences↗
USGS2024

Plan to coordinate post-earthquake investigations supported by the National Earthquake Hazards Reduction Program (NEHRP)

Introduction This report presents a plan supported by the National Earthquake Hazards Reduction Program (NEHRP) to coordinate domestic and international post-earthquake investigations (herein called “the Plan”). Post-earthquake scientific and engineering investigations are undertaken to capture critical information to understand the causes and impacts of the event, lessons from which can substantially improve the Nation’s resilience after future earthquakes. NEHRP is the Federal Government’s coordinated nationwide program to reduce risks to life and property from earthquakes. The Plan describes the activation and coordination of the four designated NEHRP Agencies in the Federal Government: Federal Emergency Management Agency, National Institute of Standards and Technology, U.S. National Science Foundation, and U.S. Geological Survey (USGS). The Plan also describes coordination between NEHRP Agencies and other organizations that may participate in pre-event and post-earthquake investigations, including non-NEHRP Federal agencies; State, regional, local, Tribal, and territorial agencies; domestic nongovernmental organizations; academic institutions and affiliated organizations; private companies; foreign governmental agencies and nongovernmental organizations; and international organizations. The Plan delineates the coordination of NEHRP post-earthquake scientific and engineering investigations to document the direct, indirect, and cascading physical and societal impacts from fault rupture and ground shaking hazards and from secondary hazards such as landslides, liquefaction, and tsunamis. In addition, the Plan identifies pre-event activities necessary to ensure that post-earthquake investigations are executed effectively. The USGS is the lead NEHRP Agency for activating and coordinating NEHRP post-earthquake investigations and for implementing this Plan. The USGS also leads coordination of the NEHRP Agencies in completing the pre-event activities identified in the Plan. The Plan has new information and supersedes USGS Circular 1242, “The Plan to Coordinate NEHRP Post-Earthquake Investigations,” which was published in 2003. This second-generation Plan was developed with the assistance of the Applied Technology Council (ATC) of Redwood City, California, under USGS contract 140G0121P0309, ATC-155. A 17-member Project Review Panel provided guidance on plan development, and input was solicited from subject matter experts representing key stakeholder groups and from participants at a public workshop.

Circular↗