USGS Science⌕ Search

SEARCH · USGS Science

Results for “Ore Geology Reviews”

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.

830 records · Page 38Linked to original sources
USGS2026

High-resolution transboundary vegetation community maps of the Sonoran and Mojave Desert ecoregion to support critical landscape conservation planning and habitat management needs

We produced a 30-m resolution binational land cover map of Bird Conservation Region 33 (BCR 33) for the U.S. North American Bird Conservation Initiative. The region covers large portions of the Sonoran and Mojave Deserts. The map can support the U.S. Fish and Wildlife Service (FWS) Migratory Bird Program’s recovery planning efforts and constitutes the first known binational land cover dataset spanning sections of the United States–Mexico border and using a consistent classification system for both countries. The mapped region includes 152 distinct land cover classes, covering a total area of 38,421,453 ha (148,345 mi 2 ), of which 13,148,345 ha (52,706 mi 2 ) are located in Mexico and 24,770,640 ha (95,639 mi 2 ) in the United States. We primarily used Landsat 8 (OLI) imagery, supplemented by limited ground surveys from two field campaigns, drone-based aerial data, and existing vegetation classification frameworks from both countries. The classification applied a data-fusion approach integrating 30-m Landsat 8 imagery, decadal phenology metrics from vegetation indices, and a random forest model trained mainly with datasets from a comprehensive national mapping project from the U.S. Geological Survey (USGS) GAP Analysis Project (GAP) and federal wildland fire agencies’ Landscape Fire and Resource Management Planning Tools (LANDFIRE) (GAP/LANDFIRE) [United States side] and the National Institute of Statistics and Geography (INEGI) [Mexico side] as well as land cover maps and opportunistic open-access and field observations. Mapping of the full BCR 33 region was carried out in two phases: 1) Phase I, the prototype map, covered a smaller portion of the transboundary area and identified 31 land cover classes, and 2) Phase II, the full BCR 33 map (refer to Figure 1), which resulted in 152 land cover classes. Using a Random Forest classifier, we achieved an overall prediction accuracy of 92% for the Phase I map and 87% for the Phase II full region map. This slight decrease can be attributed to working on a larger, more complex area with a greater number of land cover classes. No formal validation was conducted, aside from using a subset of the collected field observations and training data to assess model performance during and after training. The training sites were further verified using Google Earth (Google, 2026) imagery. Two undergraduate students who worked for over a year visually inspected imagery and open access public images to confirm each training site during model training using in-house developed, online, visual tools. A portion of this field training data was reserved for model validation, and the corresponding results are to be presented in later sections. The project developed an end-to-end, medium- and fine-resolution remote sensing–based data fusion mapping approach. This effort produced a map (Nagler et al., 2025) and the online tools to support a dynamic, live, online map for visualizing the transboundary vegetation communities in BCR 33. The toolset is currently hosted by the University of Arizona (UofA) Vegetation Index and Phenology (VIP) Lab to support FWS partners (https://vip.arizona.edu/viplab_data_explorer?LCM_BCR33). The online map is designed to allow rapid updates using new training, validation, or correction data, making it dynamic and maintainable. The approach we took established a framework for rapid updating and correction of land cover maps, as the model can be quickly retrained with new field observations, updated training data, or other sources. This enables dynamic mapping and change detection of the region’s vegetation. This framework is an advance in data fusion and crowdsourced mapping of complex, vulnerable regions, providing support to regional stakeholders and the wider user community. This transboundary map can inform the protection, conservation, and restoration of vegetation, habitat, and ecosystems, particularly for threatened and endangered species across the two nations using consistent and harmonized binational mapping systems. Beyond supporting land management decisions and stakeholders in the transboundary desert ecoregions, this BCR 33 mapping effort establishes a foundation for future rapid, low-cost, cross-border land cover mapping that can benefit and advance ecosystem management.

Arizona, Baja California, California, Nevada, Sina↗
USGS2026

Competition between invasive round goby and native slimy sculpin in a coldwater laboratory setting

Slimy sculpin ( Cottus cognatus ) abundance has substantially declined in many regions of the Laurentian Great Lakes coincident with the arrival and range expansions of invasive round goby ( Neogobius melanostomus ). Previous studies have shown that round goby can outcompete native mottled sculpin ( C. bairdii ) for shelter and negatively affect native spoonhead sculpin ( C. ricei ) growth, but direct behavioral interactions of slimy sculpin and round goby have not been observed. We conducted a laboratory experiment in which a pair of slimy sculpin were placed in an aquarium that contained an artificial spawning shelter and one round goby (treatment, n = 10) or no other fish (control, n = 10) for at least 30 days at temperatures reflective of their potential offshore overlap. We measured slimy sculpin egg deposition, aggressive interactions, artificial shelter occupancy, growth, and survival. Slimy sculpin were less likely to occupy shelters and experienced increased intraspecific competition when a round goby was present. In contrast, round goby rarely chased slimy sculpin, slimy sculpin were more likely to chase round goby, and slimy sculpin growth, egg deposition, and survival were not significantly different between groups. Differences in density between treatments and the relative size of fish could confound our observations. However, the dominant slimy sculpin occupied the shelter at similar rates between experimental groups if the round goby was smaller. Our results indicate that larger round goby overlapping with slimy sculpin could affect slimy sculpin shelter and space use, but the Great Lakes environment may create more stress for slimy sculpin than we observed.

Lake Michigan, Lake Superior↗
USGS2025

Meet the people where they are: Assessing user needs for aftershock forecast products in El Salvador, Mexico and the United States

Aftershock forecasts can help communities reduce their seismic risk by conveying how many aftershocks can be expected following a large earthquake, and how the expected number of aftershocks and their corresponding ground shaking evolves over time and space. Prior work finds that graphical forecast products may communicate such information better than only text or numbers. To identify which visual products can serve multiple user groups, we held workshops with members of several professions, including emergency managers, engineers, critical infrastructure operators, public health specialists, science communicators, and more. We conducted these workshops in El Salvador, Mexico and the United States to understand which forecast products may be effective across different countries. In these workshops, users performed small-group activities to elicit the types of aftershock forecast information that would support decisions in their respective roles and how this information would optimally be displayed. Maps of shaking hazards were frequently requested across all professions and countries, even for dissimilar forecast uses. The design of these maps, and other product needs, were differentiated by profession, country and other factors. Other forecast products, including those showing forecasts about the magnitudes and time periods of future aftershocks, also served a variety of users, but for different types of decisions. We found a greater variation in user needs by profession than by country, and that user needs also vary with time, communication channel and other contextual factors. We discuss practical implications for user-centered visual communication of operational aftershock forecasts.

International Journal for Disaster Risk Reduction↗
USGS2026

Wavelet Inversion for SliP (WISP): Open-source earthquake slip modeling software

Models of the spatiotemporal evolution of earthquake slip, termed finite-fault models, are a critical component of rapid earthquake and tsunami response, earthquake forecasting, seismic ground-motion estimates, and studies of earthquake kinematics. Here, we detail a newly released finite-fault modeling software, Wavelet Inversion for SliP (WISP), in use at the U.S. Geological Survey’s National Earthquake Information Center (NEIC) and available to the public. WISP version 1.1.0 allows inversion of teleseismic body and surface waves, as well as local strong-motion, static and dynamic Global Navigation Satellite System, and satellite imagery (e.g., Interferometric Synthetic Aperture Radar) observations on single or multiple planar fault segments. The software is used in NEIC rapid response of earthquakes M w ≥ 7, generally resulting in a published model within the first few hours after the event origin time. The rupture location and dimensions are then used as inputs to downstream products to estimate earthquake shaking, predict loss, and model the likelihood of secondary hazards, namely landslides and liquefaction. WISP is also used in research studies to evaluate the characteristics of complex ruptures including multifault ruptures and earthquake doublets, among others. The WISP version 1.1.0 software release is composed of Python-wrapped FORTRAN code to accomplish the inversion procedure. A simple command line interface facilitates ease of use even for those with only a cursory knowledge of Python scripting. WISP version 1.1.0 includes a Jupyter Notebook tutorial demonstrating use of the software for modeling the 2015 M w 8.3 Illapel, Chile, earthquake. In parallel with the tutorial, we demonstrate the typical usage of the WISP software using the M w 8.3 Illapel earthquake example here.

Seismological Research Letters↗
USGS1998

Design, revisions, and considerations for continued use of a ground-water-flow model of the coastal plain aquifer system in Virginia

A digital numerical model of the ground-water-flow system in the Virginia Coastal Plain was constructed as part of the Regional Aquifer System Analysis program of the U.S. Geological Survey and subsequently revised based on reinterpretations of the hydrogeologic framework. The revised model has been incorporated by the Virginia Department of Environmental Quality as a means to evaluate the effects of existing and proposed ground-water withdrawals. The revised model, which is documented here, includes many design features of the original model. Principal modifications are (1) reassignment of some model layers to represent different aquifers, (2) replacement of applied recharge and the uppermost surface-water layer by constant-head cells to represent the unconfined aquifer, and (3) revision of the lateral extents and transmissivity and vertical leakance values assigned to model layers. Ground-water levels simulated by the revised model under steady-state conditions, using withdrawal rates representing the period 1978-80, generally are similar to the water levels simulated by the original model under transient conditions using the same withdrawal rates. Differences probably result largely from changes in transmissivity and vertical leakance values. T he revised model enables the continued evaluation of development of the large and complex aquifer system. Improved information on parts of some aquifers has been incorporated. Only large, regional-scale trends in ground-water water levels, flow directions, and flow rates are simulated, however, and not local trends and (or) short-term changes in water level and flow. Changing hydraulic stresses, knowledge of hydrogeologic conditions, and resource-management needs require periodic revision of the model to maintain its usefulness for the management of ground-water resources in the Virginia Coastal Plain.

Virginia↗
USGS2026

Distinguishing natural from mining-related metal sources by including streambank groundwater data in a stream mass loading study

Distinguishing stream metal loading caused by mine features from that caused by natural background sources remains challenging, yet this distinction is essential for making effective remedial decisions at many legacy mine sites. We combine a stream tracer injection and synoptic sampling study with data from shallow near-stream groundwater wells to estimate left-bank versus right-bank metal loading contributions at the 100-m spatial scale. The study was performed in the East Mancos River, a mountain headwater stream in Colorado, USA. The dominant source of elevated stream metal concentrations could be either groundwater infiltration through right-bank Doyle Mine waste piles or natural acid rock drainage from hydrothermally altered bedrock located mainly on the left bank. For the five metals of concern (Cu, Al, Zn, Cd, and Mn), we find that 15 % of the load contributed by diffuse groundwater inputs in the section potentially influenced by Doyle mine waste originates from the right bank. This right-bank potential mine contribution equates to only 3 % of the total watershed instream load for these metals. Furthermore, apparent 3 H/ 3 He groundwater ages in segments contributing most of the right-bank metal loading are sufficiently old (9–12 yr) to suggest that infiltration through the waste piles, located only 140–180 m from the stream, is unlikely. Estimated potential Doyle mine loading contributions can therefore be considered maximum values. Study results thus indicate that Doyle mine waste piles are a minor source of metal loading under low-flow conditions, and streambank groundwater data can provide valuable additional information in stream mass loading studies.

Colorado↗
USGS2026

Biome-scale spatial patterns of avian abundance reveal proactive conservation opportunities in North American grasslands

North American grassland birds have experienced steeper population declines than any other avian guild, yet conservation efforts remain largely reactive and fragmented. We used nearly four decades of North American Breeding Bird Survey data to identify biome-scale spatial patterns (clustering) of grassland bird abundance for the Great Plains. Our results reveal an ecological core in the north-central Plains where community-level abundance is either increasing by >100% or remains high and stable, providing a strategic roadmap for a “Defend the Core” conservation approach. This approach flips the script from reactive triage centered on isolated population fragments to a proactive strategy of maintaining large-scale ecosystem integrity. Conversely, we found that population losses are more spatially clustered than wins, reflecting the relentless, one-way movement of woody encroachment and agricultural conversion. This asymmetry supports prioritizing intact landscapes, as current restoration rates are often outpaced by the scale of habitat loss. Notably, we found that community-level spatial clustering is a more robust indicator of biome condition than trends of individual flagship species, suggesting that managing for ecosystem integrity provides a more effective multi-species umbrella. Given our results, there is an opportunity for operationalizing a Great Plains Conservation Design that is ecosystem-centric and rooted in the sustainability of the private-land cattle production that maintains these open spaces. By leveraging avian abundance as a biological sensor, managers and producers can deploy a shared vision that matches the spatial scale of the threats, moving from reactive triage to proactive defense of core working grasslands in North America.

Great Plains biome↗
USGS2026

Re-examining the growing degree day minimum for grass carp spawning initiation in Lake Erie: Have we missed the mark?

Understanding the timing of life history processes of invasive fish species, such as adult aggregations during spawning, allows natural resource managers to optimize targeted population control efforts. Temperature-based estimates of the timing of spawning initiation for grass carp ( Ctenopharyngodon idella ) exist from its native range, but have not previously been developed for grass carp reproduction in its invaded range in the Laurentian Great Lakes of North America. We accounted for differences among methods across studies and calculated cumulative growing degree days (GDD) for the dates of grass carp egg capture in Lake Erie tributaries (Sandusky, Maumee, and Huron rivers) during 2015–2024 and compared these values to the existing GDD base 15°C requirement estimate from the native range (633 GDD–15°C). The earliest egg detection in Lake Erie (Sandusky River) was 19 GDD–15°C or 361 GDD–5°C (latest egg detection was 2285 GDD–5°C from the Sandusky River), indicating that spawning was initiated up to 25 days earlier than would be predicted by the native range GDD. Moreover, spawning in Lake Erie tributaries occurred earlier than predicted for 58% of spawning events examined based on the native range GDD. The GDD values for spawning initiation in Lake Erie tributaries can be used to refine the timing of control efforts targeting pre-spawn and spawning adult fish, inform the timing of egg sampling, and provide linkages to existing early life stage models to infer the length of the growing season available for age-0 fish to reassess tributary suitability.

Lake Erie, Maumee, River, Sandusky river↗
USGS2025

The 3D National Topography Model Call for Action—Part 2: The Next Generation 3D Elevation Program

The three-dimensional (3D) National Topography Model initiative to integrate elevation and hydrography data includes the next generation of hydrography data from the 3D Hydrography Program and the next generation of elevation data from the 3D Elevation Program (3DEP). The first-ever collection of light detection and ranging (lidar) data for the nation (IfSAR for Alaska) provides a critical baseline reference, and the addition of multiple repeat elevation mapping projects as part of the next generation of 3DEP would substantially expand analysis capabilities. As the U.S. Geological Survey (USGS) is closing in on our goal of complete coverage with 3DEP data available or in progress for 98.3 percent of the Nation at the end of fiscal year 2024, the USGS is already transitioning to the next generation of 3DEP. Based on the 3D Nation Study results and input from a broad range of stakeholders, the USGS National Geospatial Program has finalized a new design for 3DEP that provides increased lidar quality levels and refresh rates. The new program is designed with more flexibility to meet changing user needs and take advantage of improvements in mapping technologies. The program will aim to expand the level of interagency coordination for topobathymetric lidar acquisition for inland rivers. The next generation of 3DEP will also aim to emphasize research, including advancing program design, products, and services and engaging and leveraging the evolving 3D industry. Research goals also include becoming more flexible in meeting user needs and taking advantage of evolving remote-sensing technologies. The program also plans to move from focusing on producing standard products to producing a concept of a 3D Nation Ecosystem with a variety of inputs, products, and services.

Circular↗
USGS1989

Techniques for simulating flood hydrographs and estimating flood volumes for ungaged basins in east and west Tennessee

A dimensionless hydrograph developed for a variety of basin conditions in Georgia was tested for its applicability to streams in East and West Tennessee by comparing it to a similar dimensionless hydrograph developed for streams in East and West Tennessee. Hydrographs of observed discharge at 83 streams in East Tennessee and 38 in West Tennessee were used in the study. Statistical analyses were performed by comparing simulated (or computed) hydrographs, derived by application of the Georgia dimensionless hydrograph, and dimensionless hydrographs developed from Tennessee data, with the observed hydrographs at 50 and 75% of their peak-flow widths. Results of the tests indicate that the Georgia dimensionless hydrography is virtually the same as the one developed for streams in East Tennessee, but that it is different from the dimensionless hydrograph developed for streams in West Tennessee. Because of the extensive testing of the Georgia dimensionless hydrograph, it was determined to be applicable for East Tennessee, whereas the dimensionless hydrograph developed from data on streams in West Tennessee was determined to be applicable in West Tennessee. As part of the dimensionless hydrograph development, an average lagtime in hours for each study basin, and the volume in inches of flood runoff for each flood event were computed. By use of multiple-regression analysis, equations were developed that relate basin lagtime to drainage area size, basin length, and percent impervious area. Similarly, flood volumes were related to drainage area size, peak discharge, and basin lagtime. These equations, along with the appropriate dimensionless hydrograph, can be used to estimate a typical (average) flood hydrograph and volume for recurrence-intervals up to 100 years at any ungaged site draining less than 50 sq mi in East and West Tennessee.

Tennessee↗
USGS2025

Projecting stream water quality using Weighted Regression on Time, Discharge, and Season (WRTDS): An example with drought conditions in the Delaware River Basin

Future water availability depends on understanding the responses of constituent concentrations to hydrologic change. Projecting future water quality remains a methodological challenge, particularly when using discrete observations with limited temporal resolution. This study introduces Weighted Regression on Time, Discharge, and Season for Projection (WRTDS-P), a novel, computationally efficient method that enables the projection of daily stream water quality under varying hydrologic conditions using commonly available discrete monitoring data. WRTDS-P model performance was validated using 39 sites in the Delaware River Basin (DRB) and four key constituents: specific conductance (SC), nitrate (NO 3 − ), magnesium (Mg 2+ ) and calcium (Ca 2+ ). Projections were tested against holdout data from the final 1 to 5 years of each time series, demonstrating robust predictive capability, with median Nash-Sutcliffe efficiencies of 0.67 for SC, 0.56 for NO 3 − , 0.65 for Ca 2+ , and 0.79 for Mg 2+ . Model uncertainty was correlated with indicators of hydrologic or geochemical mass-sinks, such as groundwater storage and adsorption in wetland soils. Drought scenario analyses for SC used ranges of reduced discharge including flows from the 1965 drought of record. Scenarios predicted widespread increases of SC, especially in southern DRB streams where baseline SC levels are already elevated. Fractional increases of SC were more uniformly distributed, indicating potential risk to sensitive ecosystems. Notably, drought-induced SC increases were positively correlated with interannual SC trends, indicating that hydrologic extremes could exacerbate ongoing salinization. This work provides a transferable and interpretable framework for projecting future water quality and assessing hydrologic risk to water resources and aquatic ecosystems.

Delaware, New Jersey, New York, Pennsylvania↗
USGS2026

On-demand global Landsat evapotranspiration product: Development, evaluation, and dissemination

Global actual evapotranspiration (ET) is one of the essential climate variables needed to understand and manage the relationships among food, energy, and water resources. The U.S. Geological Survey Earth Resources Observation and Science (EROS) Center launched a provisional ET product in 2020, offering on-demand, field-scale global coverage derived from Landsat data through the EROS Science Processing Architecture (ESPA) platform. The ESPA interface provides ET data for cloud-free Landsat overpasses starting in 1982 with Landsat 4 through the current Landsat 9. The ET data are delivered as a Provisional Level-3 Science product created using the Operational Simplified Surface Energy Balance (SSEBop) model. Landsat surface temperature and reference ET are the main model drivers along with vegetation index and net radiation for model parameterization. A large volume of Landsat-based ET orders (e.g., over 1,200,000 images from June 2020 through December 2025) around the world indicate increasing awareness and application of the ET data. The ESPA platform enables land and water resource managers and researchers to access a first-order ET product without requiring advanced knowledge of remote sensing technology or evapotranspiration modeling. We present the methodology and workflow of the on-demand Landsat ET product and its performance evaluations over diverse hydro-climatic settings. The product can help estimate field-scale consumptive water use and thus quickly and consistently assess historical water use, allocation, and budget to inform water management under changing environments. Future ET data aggregated to monthly and seasonal time scales are expected to enhance integration with decision-making tools and procedures.

Remote Sensing of Environment↗
USGS2025

International data gaps at the Center for Engineering Strong Motion Data

The Center for Engineering Strong Motion Data (CESMD) is utilized by seismologists, engineers, and disaster management professionals in the US and has historically achieved and distributed waveforms from across the globe for significant earthquakes. The increased access to the waveforms via Web API (Application Programming Interface) offers a unique opportunity to provide the community complete datasets, sampling a variety of tectonic environments and geologic conditions, increasing the number of available ground motion records for use in ground motion models (GMMs) and improving the accuracy of earthquake engineering evaluations. The objective of this study is to programmatically identify gaps in global event data from the past decade and backfill missing data gaps at CESMD. We first compare the CESMD catalog with the Advanced National Seismic System (ANSS) Comprehensive Earthquake Catalog identifying regions and time periods where strong-motion data is limited or inadequate. To backfill datasets at CESMD for significant events, we pinpoint regions and time intervals that lack information, creating a list of events for which we’d like to obtain data. An important facet of this work is identifying the source of data and metadata across earthquake repositories around the world and integrating these data repositories into our current strong-motion data processing workflow. In parallel with these newly processed datasets, we are developing a script to produce data origination citations to include provenance and attribution information to associate with respective datasets at CESMD. We showcase our methodology for identifying and filling data gaps at CESMD using three case studies (the 2018 Anchorage Alaska earthquake sequence, seismicity associated with the 2018 Hawaiian Kilauea volcano eruption, and several earthquakes in Turkey) and then outline our strategy to apply our data gap backfilling methods on an international scale.

Conference Paper↗
USGS2026

New measurements indicate that natural geologic methane emissions from microseepage in the Michigan Basin are likely negligible

The magnitude of natural geologic methane (CH 4 ) emissions to the atmosphere (including emissions of fossil CH 4 from offshore and onshore gas and oil seeps, diffuse microseepage, mud volcanoes, volcanic vents, and geothermal areas) is highly uncertain. The largest component of geologic emissions is thought to be microseepage, which is the diffuse flux of CH 4 from soils across large areas of productive hydrocarbon basins. The accuracy of existing bottom-up estimates of microseepage is limited by low spatial coverage of published microseepage measurements. We present the first soil—atmosphere CH 4 flux measurements from Michigan Basin, USA. Results from 335 measurements taken during summer and winter seasons across a large portion of the basin suggest microseepage is nonexistent in the sampled region. Even areas with predictive features for microseepage (e.g., underlying mature, organically rich source rocks, proven gas accumulations, faults, and lineaments) yield null or negative fluxes, suggesting that CH 4 emissions from microseepage are negligible throughout our study region. A Monte Carlo method was used to place an upper limit on the regional-mean microseepage flux, in which synthetic patchy microseepage distributions were generated and tested against our observations to assess the impact of possible emission hot spots that were missed by our sampling strategy. Our analysis finds it is very unlikely that regional-mean emissions are as high as assumed in a previous global microseepage study. The observed lack of seepage may be explained by groundwater flow, active methanotrophy, glacial sediments, and bedded salt deposits, which could inhibit vertical gas migration and release to atmosphere.

Michigan↗
USGS2026

Postglacial ecosystem development of a hydrothermal landscape in Yellowstone National Park

The Yellowstone geo-ecosystem has been the subject of much research, but the ecological history of the Yellowstone Plateau volcanic field and its iconic geyser basins is less known. In this investigation, paleoenvironmental analyses of sediment cores from lakes in Lower Geyser Basin were compared with regional records and paleoclimate model simulations to reconstruct the vegetation, wildfire, limnology, hydrothermal dynamics, and climate drivers since deglaciation, 15,000 to 14,000 years ago. Pollen data from Lower Geyser Basin lakes reveal the strong influence of infertile rhyolitic soils on vegetation history: an initial late-glacial steppe was replaced by lodgepole pine forest from 12.8 to 11.0 ka, with little change in forest composition or cover thereafter despite changing climate. This stability contrasts with the more dynamic vegetation response on nonrhyolite substrates in the Yellowstone region where nutrient and moisture availability is greater. Highest wildfire activity and low lake nutrient levels in Lower Geyser Basin occurred from 12 to 4 ka, when summers were substantially warmer and drier and fire-inducing vapor pressure deficits were 29 to 56% higher. The hydrothermal history, inferred from sedimentary arsenic and cesium abundances, was spatially and temporally variable but lake-forming hydrothermal events align with periods of abundant moisture. Thus, long-term changes in wildfire, limnology, and, to some extent, hydrothermal activity were governed by insolation-driven climate variations, whereas the vegetation response was muted and constrained by geologic processes. These findings suggest that warmer, drier conditions in the future could result in less hydrothermal activity yet little change in forest cover across the Yellowstone Plateau volcanic field despite more wildfires.

Idaho, Montana, Wyoming↗
USGS2026

Neutron scattering reveals fractionation of natural gas mixtures in unconventional petroleum reservoir pores: Perspectives on energy resource recovery and storage

In unconventional petroleum reservoirs hydrocarbon fluids are hosted by both mineral and organic matter pores. These pores can have diameters that range from microns to less than a single nanometer and, for unconventional reservoirs, there is evidence that small pores ( <20 nm diameter) may constitute a large proportion of the available space. Understanding subsurface volumes and how fluids behave in them can be helpful for predicting hydrocarbon production and storage in the subsurface. One area with knowledge gaps regarding hydrocarbon behavior in small pores is the possibility for mixtures to fractionate (i.e., unmix) based on pore size or pore type. Mixture fractionation as a function of pore size could impact recovery of hydrocarbons, drive compositional shifts during production, and limit fluid storage within candidate reservoirs. To investigate natural gas fractionation in small geologic pores, we applied total neutron scattering to probe methane-ethane mixtures at reservoir pressures (up to ≈30 MPa) and temperature (60°C) within a sample from the Upper Cretaceous Niobrara Formation. Neutron scattering data reveal only minor fractionation occurs between methane and ethane in 20-nm diameter sample mesopores. Increased fractionation is observed for sample micropores, with up to 72% (±1% at 1-sigma) methane found in 2 nm diameter pores following injection of a 50%-50% methane-ethane mixture. These data provide rarely available direct experimental observations of hydrocarbon mixture behavior under nanoconfinement in a sample from an important unconventional petroleum reservoir. Our results are discussed in the context of evaluating hydrocarbon resources in unconventional reservoir meso- and micropores, reconciling observed gas composition changes during production, and more broadly, understanding subsurface pore volumes within an energy storage framework.

Fuel↗
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

Factors contributing to pesticide contamination in riverine systems: The role of wastewater and landscape sources

Wastewater treatment plant (WWTP) discharges can be a source of organic contaminants, including pesticides, to rivers. An integrated model was developed for the Potomac River watershed (PRW) to determine the amount of accumulated wastewater percentage of streamflow (ACCWW) and calculate predicted environmental concentrations (PECs) for 14 pesticides in non-tidal National Hydrography Dataset Plus Version 2.1 stream segments. Predicted environmental concentrations were compared to measured environmental concentrations (MECs) from 32 stream sites that represented a range of ACCWW and land use to evaluate model performance and to assess possible non-WWTP loading sources. Statistical agreement between PECs and MECs was strongest for insecticides, followed by fungicides and herbicides. Principal component analysis utilizing optical fluorescence and ancillary water quality data identified wastewater and urban runoff sources. Pesticides that indicated relatively larger sources from WWTPs included dinotefuran, fipronil, carbendazim, thiabendazole, and prometon whereas imidacloprid, azoxystrobin, propiconazole, tebuconazole, and diuron were more related to urban runoff. In addition, PECs generally comprised a low proportion of MECs, which indicates possible dominant loading sources beyond WWTP discharges. Cumulative potential toxicity was higher for sites with greater ACCWW and/or located in developed areas. Imidacloprid, fipronil, and carbendazim accounted for the largest portion of predicted potential toxicity across sites. The chronic aquatic life toxicity benchmarks for freshwater invertebrates were exceeded for 82 % of the imidacloprid detections ( n = 28) and 47 % of the fipronil detections ( n = 19). These results highlight the ecological implications of pesticide contamination from WWTP discharges and also the potential legacy effects from accumulated soil and groundwater sources. Pesticide management strategies that mitigate both current and historical impacts may improve the health of aquatic ecosystems.

Potomac River watershed↗