USGS ScienceSearch

SEARCH · USGS Science

Results for “Applied and Environmental Microbiology”

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

Using a temporary emigration model to estimate abundance of stream fishes from hybrid removal surveys with and without block nets

Monitoring programs are often faced with a decision to allocate resources into either robust spatiotemporal coverage to estimate a population index (e.g., not true abundance) or confirming closed sampling conditions (e.g., with block nets) for an unbiased population estimate at the cost of spatiotemporal coverage. However, making accurate and precise abundance estimates at robust spatiotemporal scales is possible when combining open and closed sampling designs with integrated modeling techniques. We used simulations and a case study of backpack electrofishing surveys in the Santa Ana River, California to test the efficacy of an integrated abundance model (temporary emigration model, TE) to estimate abundance of fishes using removal sampling methods with a hybrid sampling design (sampling with and without block nets during removal sampling). We found that the TE model performed well under most modeling scenarios (sample size, amount of closure violation, number of samples collected during closure), although at least a few samples with block nets were necessary for all parameters to be estimable. When applied to fish surveys in the Santa Ana River, we found that catch of the fishes fit to the TE model (Santa Ana Sucker, Arroyo Chub, Channel Catfish, Largemouth Bass, Yellow Bullhead) showed little evidence that the closure assumption was violated when block nets were not used. Additionally, we found that the abundance of non-native fishes negatively affected the abundance of the native Santa Ana Sucker, which was also found to adversely affect the native fish’s access to critical habitat consisting of gravel and cobble substrate. Our results indicate that the TE model presents a viable solution to common sampling problems that impact many monitoring programs, where precise and accurate population estimates can be made at large spatiotemporal scales even when most samples violate the closure assumption.

California

Deterministic physics-based earthquake sequence simulators match empirical ground-motion models and enable extrapolation to data poor regimes: Application to multifault multimechanism ruptures

We use the deterministic earthquake simulator RSQSim to generate complex sequences of ruptures on fault systems used for hazard assessment. We show that the source motions combined with a wave propagation code create surface ground motions that fall within the range of epistemic uncertainties for the Next Generation Attenuation‐West2 set of empirical models. We show the model is well calibrated where there are good data constraints, and has good correspondence in regions with fewer data constraints. We show magnitude, distance, and mechanism dependence all arising naturally from the same underlying friction. The deterministic physics‐based approach provides an opportunity for better understanding the physical origins of ground motions. For example, we find that reduced stress drops in shallow layers relative to constant stress drop with depth lead to peak ground velocities in the near field that better match empirical models. The simulators may also provide better extrapolations into regimes that are poorly empirically constrained by data because physics, rather than surface shaking data parameterizations, is underlying the extrapolations. Having shown the model is credible, we apply it to a problem where observations are lacking. We examine the case of crustal faults above a shallow subduction interface seen to break coseismically in simulations of the New Zealand fault system. These types of events were left out of consideration in the most recent New Zealand national seismic hazard model due to the modeling complexity and lack of observational data to constrain ground‐motion models (GMMs). Here, we show that in the model, by breaking up the coseismic crustal and interface rupturing fault motions into two separate subevents, and then recombining the resulting ground‐motion measures in a square‐root‐of‐sum‐of‐squares incoherent manner, we reproduce well the ground‐motion measures from the full event rupture. This provides a new method for extrapolating GMMs to more complex multifault ruptures.

Seismological Research Letters

A targeted approach for mapping groundwater discharge to surface water and fish thermal refuge in four Lake Ontario tributaries

The duration, magnitude, and frequency of heatwaves are predicted to increase in the coming decades, a combination that can reduce the survival of many fish species. Across the world, there is broad interest in identifying thermal refuge for heat-intolerant fish species and exploring opportunities to enhance or protect these areas. Because deeper groundwater maintains a relatively constant temperature, groundwater-influenced areas along streams can provide cool-water refuge for fish during periods of extreme heat. A targeted approach was developed for identifying existing cold-water zones and areas of substantial groundwater discharge in four high priority Lake Ontario tributaries. Our approach included: (1) predicting where groundwater discharge is most likely with a simple geospatial model and (2) using model predictions to select field sites for intensive high-resolution study, including ground-based mapping of groundwater features (springs, seeps, tributaries) as well as drone-based optical and thermal infrared surveys. Results from field sites were used to both verify model performance and map different types and aerial extents of thermal anomalies. Geospatial modelling successfully predicted regions of widespread groundwater upwelling, later verified and mapped by field and drone surveys. Comparison of model and field survey results further highlighted specific geospatial layers, such as soil/bedrock types and topographic wetness index, as being particularly useful for predicting groundwater influence on streams in the study area. In addition, a comparison of geospatial model results with a model of fish abundances along the studied streams showed significant positive correlations for many heat-intolerant fish species over a wide geographic area. The approach developed in this study can be applied to other watersheds to highlight areas of probable groundwater discharge and could be used by fishery and water resource managers to support cold-water fish habitat management decision-making and resource conservation.

New York

Exploring the science and data foundation for Federal public lands decisions

Public lands provide diverse resources, values, and services worldwide. Laws and policies typically require consideration of science in public lands decisions, and resource managers are committed to science-informed decision-making. However, it can be challenging for managers to use, and document the use of, science and data in their decisions. To better understand science and data use in Federal public lands decisions in the United States, we assessed the number, type, and age of documents cited in 70 Environmental Assessments (EAs) completed by the Bureau of Land Management (BLM) in Colorado from 2015–2019. We focused on the BLM, as they manage the largest area of public lands in the United States. We selected Colorado as our study area, as actions proposed on BLM lands in Colorado are representative of those across the nation. Fifty percent of citations were categorized as science and 23% as data. EAs contained an average of 17 citations (range 0–111), with documents analyzing effects of oil and gas development and recreation actions including the highest and lowest mean number of citations (41 and 6, respectively). Of individual resource analysis sections within EAs, 24% contained ≥1 science citation and 21% contained ≥1 data citation. Journal articles were the most cited type of document (26% of citations) followed by non-BLM inventories (13%). Forty-seven percent of citations were relatively recent (2010 or later); the oldest citation was from 1927. Commonly analyzed resources with the highest mean number of citations were socioeconomics, mineral resources, and noise. Fourteen of 33 commonly analyzed resources included <1 citation on average. Actions and resources with no or few citations represent opportunities for strengthening the transparent use of science and data in public lands decision-making.

Colorado

Simulation of groundwater flow to evaluate hydrogeologic controls on a PFAS plume, Coakley Landfill Superfund site, Rockingham County, New Hampshire

Per- and polyfluoroalkyl substances (PFAS), including perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS), have been detected at combined concentrations above 2,000 nanograms per liter (ng/L) at groundwater seep locations near the Coakley Landfill Superfund site, in North Hampton, New Hampshire. The landfill was active from 1972 to 1985. An impermeable cap was placed on the landfill in 1998. The adjacent area to the Coakley Landfill has many water supply wells, and transport of PFAS compounds to the wells is a concern. Fracture anisotropy in the underlying bedrock aquifer complicates the understanding of PFAS transport because groundwater preferentially travels along fractures that may not align with the prevailing groundwater flow direction. In 2018, the U.S. Environmental Protection Agency and the U.S. Geological Survey began an investigation of the groundwater flow from the Coakley Landfill site. This report describes the modification of a numerical groundwater-flow model for the local area around the Coakley Landfill and summarizes findings of the investigation. In addition, this report includes a brief description of PFOA and PFOS occurrence, a discussion of model construction, evaluation of model performance through calibration, and discussion of simulation results for two periods (before and after capping). Limitations are also discussed. Results show that simulated groundwater flow moves from the Coakley Landfill to the west and north. Advective transport modeling using particle tracking shows that groundwater from the landfill discharges primarily to streams to the west and north, and a small amount is transported to distal wells. Dilution of contaminants through advection and dispersion likely plays a role in whether PFAS compounds from the landfill will be detected above laboratory reporting levels at distal wells.

New Hampshire

Evidence of mineral alteration in a salt marsh subterranean estuary: Implications for carbon and trace element cycling

Subterranean estuaries (STE) in salt marshes are biogeochemically active zones where interactions between terrestrial groundwater and seawater drive complex cycling of carbon and trace elements, influenced by mineral dissolution. These systems, characterized by fine-grained organic-rich peat overlying permeable coastal aquifers, play a crucial role as a blue carbon sink, yet their geochemical dynamics remain poorly understood. We investigated dissolved trace elements, carbon, silica, and radium isotopes in a salt marsh STE (Sage Lot Pond, Waquoit Bay, MA) over seasonal and annual cycles. Our results reveal that groundwater and estuarine water circulation through marsh peat and aquifer sediments leads to enrichments of dissolved organic and inorganic carbon (DOC and DIC), Si, Ba, and Mn, with variable source/sink behavior of Fe and net removal of U. Submarine groundwater discharge dominated Ba fluxes, whereas pore water drainage from marsh peat acted as the main sink for U and source of Si. Fe cycling was variable, with terrestrial Fe largely removed as groundwater passed through the STE, consistent with Fe-sulfide and amorphous phase formation. Radium isotope ratios identified two distinct subsurface flow pathways, influenced by metal-oxide cycling and organic matter breakdown. Si production was decoupled from DIC, suggesting Si originates from mineral alteration, whereas DIC results from both mineral weathering and microbial respiration. Silicate mineral alteration, coupled with marsh pore water drainage, accounts for up to 16% of annual DIC exports (66 g C m −2 y −1 ), highlighting the importance of STEs in coastal carbon and trace element cycling, especially as marshes face environmental change.

Massachusetts

Earthquake magnitude and source parameter estimation with a distributed acoustic sensing dataset in the Gorda subduction zone

Distributed acoustic sensing (DAS) systems offer a cost‐effective way to create large‐scale strainmeter arrays for seismological applications using fiber‐optic cables. DAS‐based strain measurements are known to be influenced by various factors, bringing into question their general reliability for accurate earthquake characterization. A 15‐km‐long DAS deployment in northern California was operational within 3 days of the 2022 M w 6.4 Ferndale earthquake and ran continuously throughout the aftershock sequence. We utilize these aftershock data to validate DAS‐based strain measurements in two ways. We first test the accuracy of DAS‐based magnitude estimates from peak dynamic strains by comparing them with magnitude and attenuation scaling relations derived independently from traditional borehole strainmeter (BSM) data. We demonstrate that DAS‐based magnitudes are comparable to BSM‐based magnitudes when corrections for variations in site response along the fiber‐optic cable are properly made. Magnitude errors are spatially correlated, potentially because of factors such as finite‐fault effects (e.g., stress drop) or more complex, unmodeled path attenuation or because of wave propagation effects in heterogeneous media. We then apply more advanced source characterization methodology to the DAS data using a time‐domain empirical Green’s function (EGF) deconvolution approach to measure details of the moment rate history. The EGF approach using DAS data depends on careful treatment of distorting factors such as anthropogenic sources of noise and optical phase wrapping but successfully isolates source spectra for moderate‐magnitude earthquakes: source spectral ratios obtained from DAS data, broadband seismometer data, and BSM data in the same region show consistent results, revealing differences in directivity and spectral shape among earthquakes. Although further research is needed to refine source‐time‐function estimation techniques for DAS data, particularly for larger magnitude events, these case studies demonstrate the clear potential of DAS for earthquake source characterization.

California

The 3D Elevation Program—Supporting Rhode Island’s economy

Introduction High-resolution elevation data are critical to applications of landscape modeling and planning, both of which have a significant effect on Rhode Island’s economy. In these and other enterprises, program managers, while aiming to strike a balance between accuracy and cost, strive to obtain the best available elevation data to help them address a range of issues. Programs focused on climate change, environmental management, transportation design and asset management, aviation navigation and safety, riverine ecosystem management, wildlife habitat characterization and management, shellfish aquaculture, and the management and mapping of forests, parks and recreation areas, soils, wetlands, and impervious surfaces are also among the critical applications that meet the State’s management needs and depend on light detection and ranging (lidar) data that provide a highly detailed three-dimensional (3D) model of the Earth’s surface and aboveground features. The 3D Elevation Program (3DEP) is managed by the U.S. Geological Survey (USGS) in partnership with Federal, State, Tribal, U.S. territorial, and local agencies to acquire consistent lidar coverage at quality level 2 or better to meet the many needs of the Nation and Rhode Island. The status of available and in-progress 3DEP baseline lidar data in Rhode Island is shown in figure 1. 3DEP baseline lidar data include quality level 2 or better, 1-meter or better digital elevation models, and lidar point clouds, and must meet the Lidar Base Specification version 1.2 ( https://www.usgs.gov/3dep/lidarspec ) or newer requirements. The National Enhanced Elevation Assessment identified user requirements and conservatively estimated that availability of lidar data would result in at least $178,560 in new benefits annually to the State. The top 10 Rhode Island business uses for 3D elevation data, which are based on the estimated annual conservative benefits of 3DEP, are shown in table 2.

Rhode Island

Conceptualization and simulation of groundwater flow and groundwater availability in the Boone and Roubidoux aquifers in northeastern Oklahoma, 1980–2017

Oklahoma Groundwater Law (Oklahoma Statute § 82-1020.5) requires that the Oklahoma Water Resources Board conduct hydrologic investigations to determine the maximum annual yield for the State’s groundwater basins. The Boone and Roubidoux aquifers (also known as the Springfield Plateau aquifer and Ozark aquifer, respectively) are bedrock aquifers that extend from northeastern Oklahoma into Kansas, Arkansas, and Missouri. At present (2024), the Oklahoma Water Resources Board has yet to legally issue orders for the final determination of maximum annual yields for the Boone and Roubidoux aquifers. To support determination of a maxi­mum annual yield, the U.S. Geological Survey, in coopera­tion with the Oklahoma Water Resources Board, developed a hydrogeologic framework, a conceptual groundwater-flow model, and a calibrated numerical groundwater-flow model for the Boone and Roubidoux aquifers. Three types of groundwater-availability scenarios were simulated by using the calibrated numerical model. These scenarios were used to (1) estimate equal-proportionate-share groundwater withdrawal rates (groundwater withdrawal applied equally over the aquifer), (2) quantify the potential effects of projected groundwater withdrawals on groundwater storage over a 50-year period, and (3) simulate the poten­tial effects of a hypothetical 10-year drought. For the Boone aquifer, equal-proportionate-share groundwater withdrawal rates were 1.10, 0.98, and 0.96 acre-feet per acre per year for the 20-, 40-, and 50-year scenarios, respectively. For the Roubidoux aquifer, equal-proportionate-share groundwater withdrawal rates were 1.76, 1.34, and 1.25 acre-feet per acre per year for the 20-, 40-, and 50-year simulations, respectively. For the 50-year scenarios, stream seepage was minimally affected. Over the 10-year drought scenario, groundwater storage in the Boone and Roubidoux aquifers decreased by 660,451 acre-feet (6.7 percent) and 508,472 acre-feet (1.0 per­cent), respectively.

Oklahoma

The 2025 Puerto Rico and Virgin Islands U.S. National Seismic Hazard Model Update: Ground motion model selection and comparison

We evaluate, select, and describe the ground-motion models (GMMs) used in the 2025 update of the U.S. National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI). We identify the most appropriate models that align with GMM selection criteria for use in the PRVI region to improve the accuracy of seismic hazard assessments. The update incorporates globally applicable GMMs suited for the active crustal and subduction earthquakes in the region. We include region-specific adjustments to these GMMs derived from local site response analyses derived from ground motion records. The unadjusted and regionally-corrected GMMs are combined to create a robust model for predicting median ground motion. The model integrates epistemic uncertainty through a median ground motion logic tree that accounts for variations in magnitude and distance. This study compares the GMMs selected for the 2025 PRVI NSHM, including both as-provided and regionally adjusted NGA-West2 and NGA-Subduction models, with those used in the 2003 PRVI NSHM. We evaluate how changes in model selection, weighting, aleatory variability, and epistemic uncertainty influence seismic hazard estimates. Trends with distance, magnitude, and spectral period are analyzed to evaluate how the scaling behavior of the newer GMMs differs from that of earlier models. Relative to the GMMs used in the 2003 NSHM for this region, the 2025 models generally predict lower ground motions. Comparisons with additional GMMs indicate that the adjustments applied for PRVI are consistent with regional-specific modifications developed elsewhere globally. The increase in aleatory variability and epistemic uncertainty in the 2025 update results in a notable increase in hazard levels from these wider uncertainty bounds. These changes can result in as much as a 10%–20% variation in probabilistic ground motion at the 2% in 50 years exceedance level for hazard maps computed across the region for representative site classes and periods.

Puerto Rico, Virgin Islands

Don’t Let Negatives Hold You Back: Accounting for Underlying Physics and Natural Distributions of Hydrothermal Systems When Selecting Negative Training Sites Leads to Better Machine Learning Predictions

Selecting negative training sites is an important challenge to resolve when utilizing machine learning (ML) for predicting hydrothermal resource favorability because ideal models would discriminate between hydrothermal systems (positives) and all types of locations without hydrothermal systems (negatives). The Nevada Machine Learning project (NVML) fit an artificial neural network to identify areas favorable for hydrothermal systems by selecting 62 negative sites where the research team had confidence that no hydrothermal resource exists. Herein, we compare the implications of the expert selection of negatives (i.e., the NVML strategy) with a random sample strategy, where it is assumed that areas outside the favorable structural ellipses defined by NVML are negative. Because hydrothermal systems are sparse, it is highly probable that, in the absence of a favorable geological structure, hydrothermal favorability is low. We compare three training strategies: 1) the positive and negative labeled examples from NVML; 2) the positive examples from NVML with randomly selected negatives in equal frequency as NVML; and 3) the positive examples from NVML with randomly selected negatives reflecting the expected natural distribution of hydrothermal systems relative to the total area. We apply these training strategies to the NVML feature data (input data) using two ML algorithms (XGBoost and logistic regression) to create six favorability maps for hydrothermal resources. When accounting for the expected natural distribution of hydrothermal systems, we find that XGBoost performs better than the NVML neural network and its negatives. Model validation was less reliable using F1 scores, a common performance metric, than comparing probability estimates at known positives, likely because of the extreme natural class imbalance and the lack of negatively labeled sites. This work demonstrates that expert selection of negatives for training in NVML likely imparted modeling bias. Accounting for the sparsity of hydrothermal systems and all the types of locations without hydrothermal systems allows us to create better models for predicting hydrothermal resource favorability.

Geothermal Resources Council Transactions

Capturing the uncertainty of seismicity observations in earthquake rate estimates: Implications for probabilistic seismic hazard analysis and the USGS National Seismic Hazard Model

The rate of earthquakes in a region is a fundamental input to Probabilistic Seismic Hazard Analysis. We present a Monte Carlo method for computing that rate from seismicity catalogs while including a range of data and analysis uncertainties. This method is applied to regions for which the b value is assumed to be spatially invariant. Each region is broken down into epochs for which each epoch is estimated to have a uniform magnitude of completeness ( ⁠⁠ M c ). The distribution of earthquake rates for M ≥ M c is determined for each epoch by considering the Poisson likelihood of rates given the number of observed earthquakes with M ≥ M c ⁠ . We use a Monte Carlo process to include the uncertainty in b , ⁠ M c , and individual event magnitudes. The result for each epoch is the joint distribution of the Poisson rate of earthquakes with magnitudes larger than the minimum value used to calculate hazard ( ⁠⁠ M 1 ) and the Gutenberg–Richter b values, which control the extrapolation to other magnitudes. The rate for each region is either the duration‐weighted average over the epochs or, to better capture temporal variations, we also consider mixture models. The mixture models also provide an avenue to allow temporal variations in b values. To implement this joint distribution in a logic tree, we use the mean and 95% confidence branches, each of which is parameterized with an M ≥ M 1 rate and b value. We explore different ways of defining those branches, as well as non‐Gutenberg–Richter branches, and their impact on hazard estimates. The mean hazard, but not the fractiles, is robust with respect to these choices. To illustrate these new methods, we use synthetic data and catalogs from recent U.S. Geological Survey National Seismic Hazard Models for the Central and Eastern United States and for Puerto Rico and the U.S. Virgin Islands.

central and eastern United States, Puerto Rico, U.

Comparing subduction ground-motion models to observations for Cascadia

We evaluate Cascadia subduction ground-motion models (GMMs), considered for the 2023 US National Seismic Hazard Model (NSHM) update, by comparing observations to model predictions. The observations comprise regional recordings from intraslab earthquakes, including contributions from 2021 and 2022 events in southern Cascadia and global records from interface earthquakes. Since the 2018 NSHM update, new GMMs for Cascadia have been published by the Next Generation Attenuation (NGA)-Subduction Project that require independent evaluation. In the regional intraslab comparisons, we highlight a characteristic frequency dependence for Cascadia data, with short periods having lower ground motions and longer periods being comparable to other subduction zones. We evaluate differences in northern and southern Cascadia and find that the NGA-Subduction GMMs developed using southern Cascadia data perform better in this region than the model that did not consider these data. We compare ground-motion variability in Cascadia with the NGA-Subduction model predictions and find differences at short periods ( T = 0.1 s) due to the use of global versus regional data in the development of these models. Moreover, the within-event component of aleatory variability from the GMMs overpredicts the standard deviation of Cascadia recordings at very short periods ( T < 0.05 s). Using global interface earthquakes as a proxy to evaluate the Cascadia GMMs, we find long-period overprediction from a simulation-based GMM and some of the empirical GMMs. When comparing recent observations, we find a similar misfit to GMMs and the 2010 and 2022 Ferndale earthquakes. Finally, we observe different basin amplification factors arising in different subsets of the data, which indicate that differences in basin factors between empirical GMMs could arise from the data selection choices by the developers. As part of evaluating the regional basin terms, we apply basin amplification factors from the magnitude 9 Cascadia earthquake simulations to the empirical GMMs for interface earthquakes. The comparisons presented in this study indicate that the NGA-Subduction GMMs for Cascadia perform well relative to observations and older subduction GMMs.

British Columbia, California, Oregon, Washington

Pesticide occurrence in shallow groundwater in three regions of agricultural land use: Baldwin County, the Wiregrass region, and the Tennessee River valley region of Alabama, 2009–20

As part of a cooperative investigation between the U.S. Geological Survey and the Alabama Department of Agriculture and Industries, a network of 22 groundwater wells were sampled from 2014 through 2020 for about 230 pesticide and pesticide degradate compounds. Wells were located in three regions of intensive agricultural land use in Alabama: Baldwin County, the Wiregrass region, and the Tennessee River valley region. Metolachlor sulfonic acid, a degradate of the herbicide metolachlor, was the most frequently detected compound, occurring in about 70 percent of the samples. Three other compounds, metolachlor, atrazine, and 2-chloro-4-isopropylamino-6-amino-s-triazine, were also detected in more than half of the samples. Metolachlor and its degradates accounted for 33 of the 50 greatest compound concentrations study-wide, including the maximum pesticide concentration across all compounds (62,500 nanograms per liter). The frequency and magnitude of detections of many specific pesticide compounds varied among the three regions, but all detected pesticide concentrations were well below the U.S. Environmental Protection Agency maximum contaminant levels and applicable human health benchmarks. Sample results were combined with results of previous (2009–13) sampling to provide a continuous time-series of data for 2009–20. More than half of the 289 pesticide compounds analyzed during 2009–20 were not detected in any samples. Only four compounds were detected at great enough frequency throughout the 10 sampling years to evaluate patterns of change through time. Metolachlor and its degradate, metolachlor sulfonic acid, were frequently detected in all regions. Atrazine and its degradate, 2-chloro-4-isopropylamino-6-amino-s-triazine, were also detected in wells from all regions, but the variability and magnitude of concentrations were greatest in the Tennessee River valley region. No apparent temporal pattern in concentrations was found.

Alabama

Navigating uncertainty and competing objectives: Spring chinook salmon recovery in the upper Willamette river

Globally, anadromous fish populations are threatened with extinction due to multiple factors, including river impoundments that block migration, widespread alteration of physical habitat, temperature, and flow regimes, commercial and recreational harvest, changes in biological communities, and long-term climate trends. In this synthesis, we describe how physical and biological factors, policy and law, and public resource allocation affect management of spring-run Chinook Salmon ( Oncorhynchus tshawytscha ) in the Upper Willamette River (UWR), Oregon, USA. Efforts to restore salmon populations often involve trade-offs, for example, water management decisions must account for flood control, agricultural and municipal water use, and management of other species. UWR Chinook Salmon recovery efforts are also governed by multiple laws, policies, and action agencies that in some cases have incongruent objectives. Competing and conflicting objectives, combined with uncertainty regarding the outcomes of proposed management actions, create significant challenges for decision makers. Resource limitations, regulatory constraints, and incongruent objectives have created challenges for managers in the UWR, but opportunities exist for improved adaptive management and collaboration, which could build trust and engagement among stakeholders and ultimately support the successful implementation of decision tools. A suite of decision support tools that have been used in the Willamette River basin and elsewhere across the range of Pacific salmon have the potential to enhance decision processes. These tools can be effective when decision makers and stakeholders commit to long-term, collaborative utilization of the tools for decision making, and when available resources and regulatory environments are conducive to implementation and adaptive refinement of both the decision-support tools and management plans. The general biological and institutional principles discussed for UWR Chinook Salmon in this work apply more generally to populations of anadromous fishes across their range, especially to populations in highly regulated river systems.

Oregon

A methodology to estimate CO2 and energy gas storage resources in depleted conventional gas reservoirs

Depleted hydrocarbon reservoirs are subsurface geological structures capable of sequestering vast quantities of carbon dioxide (CO 2 ) as well as storing other energy gases for later usage, such as natural gas, and potentially hydrogen (H 2 ). Here we outline a methodology to quantify multi-gas storage resources in depleted conventional gas reservoirs for usage in assessments by the United States Geological Survey (USGS) at the scale of sedimentary basins. The methodology consists first of quantifying accessible pore volume in a depleted reservoir for natural gas storage using up to three equations. Input data are derived from commonly reported or estimated reservoir parameters and natural gas production volumes, and equations may be combined in linear models to improve pore volume estimates. Storage estimates from these equations are tested and validated for 31 reservoirs in the Michigan Basin Province, USA that were previously converted to underground gas storage facilities and have known (federally reported) natural gas storage capacities. Secondly, natural gas storage capacities can be transformed via fluid substitution calculations to estimate the storage resources for non-native fluids, applied here for, CO 2 , H 2 , and methane-H 2 blends, accounting for molecule-specific deviations from ideal gas behavior at reservoir pressures and temperatures as well as differing storage efficiencies. Importantly, the storage of non-native fluids may not be appropriate in all depleted gas reservoir settings due to potential risks like leakage, in particular in the case of H 2 storage, requiring additional knowledge of caprock sealing capacity. Given this caveat, we demonstrate the fluid substitution method for natural gas reservoirs of the Northern Niagaran Reef and Southern Niagaran Reef USGS plays in the Michigan Basin Province, as these trends of Silurian pinnacle reefs are capped with tight-sealing evaporite facies. The deterministic equations outlined from this methodology can be incorporated into future probabilistic USGS gas storage assessments for CO 2 , H 2 , and natural gas resources in the United States.

Michigan

Assessment of water levels, nitrate, and arsenic in the Carson Valley Alluvial Aquifer and the development of a data visualization tool for the Carson River Basin, Nevada

Residents of Carson Valley, Douglas County, Nevada, rely on the basin-fill alluvial aquifer underlying the valley for drinking water. Since the 1980s, groundwater levels and water-quality data have been collected to monitor the status of the aquifer system and to assist in planning efforts to address current (2024) and future demand. The U.S. Geological Survey (USGS), in cooperation with Douglas County, Nevada, evaluated trends in water levels, nitrate, and arsenic concentrations from a network of monitoring and domestic wells in Carson Valley. This work also assessed the monitoring well network to determine the suitability of wells for characterizing the occurrence of arsenic in the groundwater. Monitoring of constituents, such as nitrate and arsenic concentrations, is needed to assess changes in contaminant distribution and to evaluate the effect that changing land use and groundwater pumping has on their occurrence and transport. Results of the trend analysis indicate water levels are declining ( p <0.05) in 17 of 26 selected monitoring wells (65 percent). Areas with the largest change in water levels, with more than 20 feet of declines, were within the community areas of Johnson Lane, Ruhenstroth, South Agricultural, East Valley, and Fish Springs. Variations in water levels measured in wells from the Central Agricultural, Minden, Foothill, Alpine County (one well), and Gardnerville Ranchos areas show periods of increase and decrease over time, but they also maintain long-term declining trends ( p <0.05). Increases in nitrate concentrations in groundwater samples collected from 9 out of 14 selected monitoring wells (64 percent) are statistically significant ( p <0.05) within the Ruhenstroth, Gardnerville Ranchos, East Valley, Genoa, and Johnson Lane community areas. Samples collected from a well in Indian Hills/Jacks Valley indicated a decreasing trend in nitrate concentration over time. Nitrate concentrations in samples collected from wells in East Valley, Genoa, Johnson Lane, and Indian Hills/Jack Valley were consistently low (less than 3 milligrams per liter [mg/L]) and stable. Nitrate concentrations from selected wells in Johnson Lane and Garnerville Ranchos exceeded the U.S. Environmental Protection Agency (EPA) maximum contaminant level (MCL) of 10 mg/L (as nitrogen) and have trends that are increasing over time. In 2022, a sample collected from Johnson Lane had a concentration (7.3 mg/L) below the MCL with an increasing trend over time. Temporal trend analyses for groundwater arsenic concentrations in Carson Valley could not be done because of a lack of temporal data. However, using available historical data, arsenic concentrations seem to be greater in groundwater from wells located on the eastern and northern areas of the valley than in wells located on the western or southern areas. Groundwater arsenic concentrations exceed 5 micrograms per liter (μg/L) in most samples collected from wells in Johnson Lane, Airport, Central Agricultural, and East Valley areas and in many cases exceed the U.S. Environmental Protection Agency (EPA) MCL of 10 μg/L. Data indicate that groundwater from domestic wells screened at deeper intervals are likely more vulnerable to elevated arsenic concentrations than shallower wells. A groundwater network evaluation for Carson Valley identified potential modifications in the sampling locations and frequency to better understand the effect of groundwater pumping in communities where municipal and domestic demand are greatest, potentially enhancing understanding of contaminant transport in these areas. Potential modifications to the active well network include reducing the frequency of sample collection from existing network wells (6 out of 11) that have consistently shown low and stable nitrate concentrations, adding wells in areas where data are sparse, and increasing the number of wells in areas with elevated groundwater nitrate concentrations. Including the analysis of arsenic in samples from the active groundwater monitoring well network will provide more detail on the temporal and spatial variability of arsenic concentrations. A visualization tool for the Carson River Basin was developed to provide access to discrete and near real-time hydrologic and water-quality data. The Carson River Basin Hydro Mapper (CBH; U.S. Geological Survey, 2023b) shows active and historical discrete water levels measured by the USGS and the State of Nevada Division of Water Resources, discrete groundwater nitrate and arsenic concentration data collected by the USGS, near real-time streamflow, and surface water levels for select waterbodies. The hydrologic data in the CBH provides resource managers, the public, and the scientific community with an easily accessible tool to present and communicate the most up-to-date information available about local and basin-wide water resources.

Nevada

Hydrologic investigations and a preliminary conceptual model of the groundwater system at North Penn Area 1 Superfund Site, Souderton, Montgomery County, Pennsylvania

The U.S. Geological Survey (USGS) conducted hydrogeologic investigations, reviewed existing data, and developed a preliminary conceptual model of the groundwater system as part of technical support of the U.S. Environmental Protection Agency (EPA) at the North Penn Area 1 Superfund Site (hereafter, the NP1 Site) located within the Borough of Souderton in Montgomery County, Pennsylvania. Field work and monitoring took place during 2012–18. The area is underlain by sedimentary formations that form a fractured-rock aquifer used for drinking water and industrial supply. The EPA placed the Site on the National Priorities List in 1989, identifying tetrachloroethylene (PCE) and trichloroethylene (TCE) as contaminants of concern. During 2012–18, the USGS conducted field activities that included drilling an 82-foot (ft)-deep monitoring well (MG 2220) in 2016, reconstructing a 208-ft-deep former industrial production well (MG 668 [Granite Knitting Mill]), and collecting borehole geophysical and video logs and water levels from those and five additional wells, which ranged in depth from about 50 to 200 ft below land surface. Continuous water levels were collected during 2014–17, and a synoptic set of water levels were measured in April 2018 in the seven wells. The borehole geophysical logs (caliper, acoustic televiewer, natural gamma, single-point resistance, vertical flow, and fluid temperature and resistivity) and borehole video logs in the seven wells were evaluated to assess potential for lithologic correlation and to identify and describe water-bearing features, which included both low- and high-angle fractures and other openings oriented along dipping bedding planes, joints, or possible faults. Borehole geophysical logs collected by USGS in 1992 in a 300-ft-deep former production well near the Site were also evaluated. Few to no distinctive features were identified on geophysical logs (natural gamma and single-point resistance) that could be used for correlation, thus limiting this approach to determining local geologic structure. Extensive fracturing in the upper 62 ft of monitoring well MG 2220 indicates that the well was likely drilled through a zone of faulting, and other evidence of faulting is present in the area near the Site. Assessment of continuous water levels showed hydraulic connections among some wells as indicated by rising or falling water levels in response to changes in pumping rates at nearby wells. A map of water levels measured in April 2018 indicates potential for groundwater flow generally toward the stream to the south and southwest of the Site, but the limited water-level data are insufficient to describe vertical groundwater gradients or lateral gradients in any detail. Review of 1999–2022 volatile organic compound (VOC) monitoring data collected by the Pennsylvania Department of Environmental Protection for five monitoring wells indicates that the highest groundwater concentrations of PCE and TCE were found in samples from extraction well MG 2201 (S-1) downgradient from, and nearest to, the previously identified Site contaminant source area, and these concentrations fluctuated through time. PCE concentrations were higher than TCE concentrations in samples from all five monitoring wells and were much higher than TCE concentrations in samples from extraction well MG 2201 (S-1). Temporally variable recharge is a possible factor affecting observed fluctuations in PCE concentrations in groundwater samples from well extraction MG 2201 (S-1), as indicated by a general inverse relation between PCE concentrations and water levels in a nearby long-term observation well. The PCE concentration of 1,830 micrograms per liter (μg/L) in a May 2018 water sample from monitoring well MG 2220 was more than four times the PCE concentration of 444 μg/L in a December 2017 sample from the nearby extraction well MG 2201 (S-1), which is open to fewer fractures. Low concentrations of VOCs were measured in surface water at two stream sites downgradient from wells with the highest groundwater VOC concentrations at the Site, indicating that discharge of contaminated groundwater to the stream is likely. Development of a conceptual model of the groundwater system was constrained by limited data. In areas with no pumping, groundwater-flow directions generally are thought to be controlled by topography and geologic structure (bedding orientation) and likely to the south and southwest of the Site, with local flow directions affected by orientations of fractures, joints, and local faults. Additional investigations that could help improve the conceptual model of the groundwater system and help delineate the extent of groundwater contamination and its transport are discussed.

Pennsylvania