Released December 03, 2025 15:00 EST

2025, Scientific Investigations Report 2025-5087

Dominick J. Antolino, Gerard J. Gonthier, Georgina M. Sanchez

In 2019, the U.S. Geological Survey and Wake County Environmental Services began a collaborative study to evaluate groundwater resources and long-term groundwater availability in the county’s fractured-rock groundwater system. Wake County, in central North Carolina, is experiencing rapid population growth, associated land development, and changing water use. Hydrogeologic data including groundwater levels, aquifer testing, borehole fracture flow measurements, water-quality samples, and groundwater age-dating tracers were collected, along with findings from previous investigations, to help inform a conceptual model of the flow system used to develop a modular three-dimensional finite-difference groundwater-flow model (MODFLOW) for simulating historical and future groundwater conditions from 2000 to 2070.

Hydraulic conductivity and transmissivity ranges were estimated from 17 slug tests and 21 borehole-flow measurements. Groundwater-quality analytical results from 19 sampling sites indicate that oxidation-reduction (redox) conditions varied within the regolith and bedrock and that minimal evaporation occurred before recharge entered the groundwater system. Age dating revealed mixtures of older and younger water, ranging from the 1940s to the 1990s—indicating variable flow pathways of recharge within permeable bedrock fracture zones.

To simplify the complex fractured-rock groundwater system, two layers representing the regolith and the fractured bedrock were used in the MODFLOW model. Model calibration included parameter estimation and provided a reasonable fit to observed groundwater levels and estimated stream base flows. The model forecast scenarios incorporated future climate-model data for two emissions scenarios with land cover change projections to simulate potential impacts to future groundwater levels, recharge, and base flows. Recharge and base flow projections were largely within historical ranges, with no apparent long-term trends, but did indicate a slight downward shift in median values—likely, in part, because of differences in spatial resolution of input climate datasets. Seasonal patterns were consistent with historical data, with projections of possible increases in future winter recharge. Model limitations are discussed, and additional monitoring and model refinement needs are highlighted to support decision making for local groundwater management.

Released December 03, 2025 11:50 EST

2025, Fact Sheet 2025-3048

Jenny H. Lagesse, Christopher J. Schenk, Jane S. Hearon, Sarah E. Gelman, Thomas M. Finn, Benjamin G. Johnson, Tracey J. Mercier, Phuong A. Le, Heidi M. Leathers-Miller, Andrea D. Cicero, Ronald M. Drake II

Using a geology-based assessment methodology, the U.S. Geological Survey estimated undiscovered, technically recoverable mean conventional and continuous resources of 4.7 trillion cubic feet of gas in the Mesaverde Group and Lance Formation in the Southwestern Wyoming Province, Wyoming, Utah, and Colorado.

Released December 02, 2025 16:56 EST

2025, Open-File Report 2025-1042

Anya Metcalfe, Morgan Ford, Lawrence E. Stevens, Theodore Kennedy

Management practices that enhance habitat complexity in dam tailwaters can increase biodiversity and improve ecosystem health, particularly in the face of climate change and warming water temperatures. In some instances, however, management practices may be intended to simplify habitat features to impede establishment of invasive species. In Glen Canyon National Recreation Area, a backwater known colloquially as “the Slough”, located 5 km downstream of Glen Canyon Dam, is being considered for removal to reduce breeding habitat for warmwater non-native fish. In this report, we assess the occupancy of and habitat suitability for dragonflies and damselflies (order: Odonata) at the Slough. U.S. Geological Survey staff conducted site visits to the Colorado River in Glen Canyon, the Slough, and another backwater (‘Frogwater’) on 11-13, and 26 September 2024. We characterized the physical habitat of the sampling sites by recording water temperatures, specific conductance, dissolved oxygen, turbidity, flow, depth, and benthic substratum size distribution, and we sampled the assemblages of aquatic macroinvertebrates as well as riparian macroinvertebrates using benthic and aerial collection methods, respectively. We describe three distinct benthic aquatic invertebrate communities in and around the Slough, two of which contained Odonata. We found no Odonata larvae in the mainstem at Frogwater, or in the Lower Slough where it joins the Colorado River mainstem. Using historic specimen data from Museum of Northern Arizona, we report eight species of damselflies (Coenagrionidae) and eight species of dragonflies among three families (Aeshnidae, Gomphidae, Libellulidae) from Glen Canyon between 1985-2024. We discuss the habitat requirements of Odonata larvae known to occur in the Slough, as well as their cultural and recreational values. We conclude that re-engineering the Slough to cool water temperatures will likely reduce larval Odonata habitat locally, but is unlikely to affect their diversity and abundance on a regional scale.

Released December 02, 2025 10:03 EST

2025, Preprint

Salme Ellen Cook, Liv M. Herdman

Coastal flooding and the associated damages due to storms are increasing with sea level rise around the world, with regional variability in the severity of impacts., Researchers and resource managers need to better understand and predict the future shifts in coastal flooding due to these processes to plan for resilient and sustainable communities. Here we present an analysis of long-term historical records of water levels, tides, and modeled present-day wave climatologies, to characterize the present-day inundation extent in Long Island Sound and Long Island, NY. To understand the potential future changes in inundation extent, we provide a similar analysis of future climate projections of non-tidal residuals (storm surge) for the year 2050 and compare these projections with our present-day results. We examine both the magnitude of relatively frequent events with a 0.99 annual exceedance probability to more extreme events with a 0.01 annual exceedance probability (or the 1 in 100-year event). This range of events is relevant for local managers to understand the spatial variability in coastal inundation, in addition to planning for larger more catastrophic events.

Released December 01, 2025 14:45 EST

2025, Open-File Report 2025-1051

Kendra L. Russell, William J. Andrews, Amy R. McHugh

Executive Summary 

A Decree of the Supreme Court of the United States entered June 7, 1954 (New Jersey v. New York, 347 U.S. 995), established the position of Delaware River Master within the U.S. Geological Survey. In addition, the Decree authorizes the diversion of water from the Delaware River Basin and requires that compensating releases from certain reservoirs owned by New York City be made under the supervision and direction of the River Master. The Decree stipulates that the River Master provide reports to the Court, not less frequently than annually. This report is the 65th annual report of the River Master of the Delaware River. The report covers the 2018 River Master report year, from December 1, 2017, to November 30, 2018.

During the report year, precipitation in the upper Delaware River Basin was 60.39 inches or 136 percent of the long-term average. On December 1, 2017, combined useable storage in the New York City reservoirs in the upper Delaware River Basin was 193.230 billion gallons or 71.3 percent of the combined useable storage capacity of 270.837 billion gallons. The reservoirs had a usable capacity of 99.5 percent on May 31, 2018. Combined storage remained high (above 80 percent combined capacity) and did not decline below 80 percent of combined capacity through November 30, 2018. River Master operations during the year were conducted as stipulated by the Decree and the Flexible Flow Management Program.

Diversions from the Delaware River Basin by New York City and New Jersey fully complied with the Decree. Reservoir releases were made as directed by the River Master at rates designed to meet the flow objective for the Delaware River at Montague, New Jersey, on 42 days during the report year. Interim Excess Release Quantity banks and conservation releases, designed to relieve thermal stress and protect the fishery and aquatic habitat in the tailwaters of the reservoirs, were also made during the report year.

Released December 01, 2025 10:37 EST

2025, Nature Microbiology (10) 3059-3068

Julia Kelliher, Chloe Mirzayi, Sarah R. Bordenstein, Aaron Oliver, Christina A. Kellogg, Eneida L. Hatcher, Maureen Berg, Petr Baldrian, Mashael Aljumaah, Cassandra Maria Miller, Christopher Mungall, Vlastimil Novak, Alexis Palucki, Ethan Smith, Nazifa Tabassum, Gregory Bonito, J. Rodney Brister, Patrick S. Chain, Jose Pablo Dundore-Arias, Joanne B. Emerson, Vanessa Moreira Fernandes, Roberto Flores, Antonio Gonzalez, Zoe A. Hansen, Scott A. Jackson, Ahmed M. Moustafa, Trent R. Northen, Nonia Pariente, Jennifer Pett-Ridge, Sydne Record, Linta Reji, Anna-Louise Reysenbach, Virginia I. Rich, Lorna Richardson, Simon Roux, Lynn M. Schriml, Reed S. Shabman, Maria A. Sierra, Matthew B. Sullivan, Punithavathi Sundaramurthy, Katherine M. Thibault, Luke R. Thompson, Scott W. Tighe, Ethell Vereen, STREAMS Consortium, Emiley A. Eloe-Fadrosh

The interdisciplinary nature of microbiome research, coupled with the generation of complex multi-omics data, makes knowledge sharing challenging. The Strengthening the Organization and Reporting of Microbiome Studies (STORMS) guidelines provide a checklist for the reporting of study information, experimental design and analytical methods within a scientific manuscript on human microbiome research. Here, in this Consensus Statement, we present the standards for technical reporting in environmental and host-associated microbiome studies (STREAMS) guidelines. The guidelines expand on STORMS and include 67 items to support the reporting and review of environmental (for example, terrestrial, aquatic, atmospheric and engineered), synthetic and non-human host-associated microbiome studies in a standardized and machine-actionable manner. Based on input from 248 researchers spanning 28 countries, we provide detailed guidance, including comparisons with STORMS, and case studies that demonstrate the usage of the STREAMS guidelines. STREAMS, like STORMS, will be a living community resource updated by the Consortium with consensus-building input of the broader community.

Released December 01, 2025 10:13 EST

2025, Frontiers in Earth Science (13)

Michelle R. Plampin, Steven T. Anderson, Stefan Finsterle, Ashton M. Wiens

Since the vast majority of carbon dioxide (CO₂) storage resources in the United States are in deep saline aquifers, optimizing the use of these saline storage resources could be crucial for efficient development of geologic CO₂ storage (GCS) resources and basin- or larger-scale deployment of GCS in the country. Maximum CO₂ injection rates can be enhanced by extracting brine from the CO₂ storage unit. However, disposal of the extracted brine is both a technological and economic challenge. The lowest-cost option would likely be reinjection of the extracted brine into another formation above or below the CO₂ storage unit. Therefore, it is important to estimate brine injectivity as it will constrain the potential to increase CO₂ injectivity at an injection site that has access to multiple geologic storage units where either CO₂ or brine can be injected. Using a simulation-optimization framework, coupled with a non-isothermal, multiphase CO₂-water-salt equation-of-state module, we developed a computationally efficient method for evaluating optimization of simultaneous CO₂ injection, brine extraction, and brine (re)injection at hypothetical injection sites deployed across a geologic basin. The Illinois basin is ideal for testing our methodology because it contains multiple geologic storage units with seals in between them to isolate injection of CO₂ in one unit from interfering with the injection of either brine or CO₂ in another unit above or below it. In addition, we investigated the relative effects of variation in key geologic parameters as well as two reservoir structures (hydrogeologic heterogeneity/anisotropy and homogeneity/isotropy) on CO₂ injectivities and enhancement of CO₂ injectivity through extracting brine. Results suggest that permeability, depth, and especially thickness of the storage unit could be the most influential parameters determining CO₂ injectivity. They also suggest that only injecting CO₂ into the storage unit with the greatest injectivity, enhancing that unit’s injectivity by extracting brine, and disposing of the produced brine in other suitable units could maximize total CO₂ injectivity in limited regions of the basin. At the majority of simulated injection sites, however, we found that injecting CO₂ into all of the accessible and suitable storage units was more likely to maximize the CO₂ storage resource.

Released December 01, 2025 09:37 EST

2025, Scientific Investigations Report 2025-5097

Kirsten E. Faulkner, Bryant C. Jurgens

More than 2 million Californians rely on groundwater from domestic wells for drinking-water supply. This report summarizes a 2022 California Groundwater Ambient Monitoring and Assessment Priority Basin Project (GAMA-PBP) water-quality survey of 33 domestic and small-system drinking-water supply wells in the Gilroy-Hollister Valley groundwater basin and the surrounding areas, where more than 20,000 residents are estimated to utilize privately owned domestic wells. The study area includes the Llagas subbasin in the north, the North San Benito subbasin in the south, and the surrounding uplands. The study was focused on groundwater resources used for domestic drinking-water supply, which are mostly drawn from shallower parts of aquifer systems rather than those of groundwater resources used for public drinking-water supply in the same area. This assessment characterized the quality of ambient groundwater in the aquifer before filtration or treatment, rather than the quality of drinking water delivered to the tap.

To provide context, the measured concentrations of constituents in groundwater were compared to Federal and California State regulatory and non-regulatory benchmarks for drinking-water quality. A grid-based method was used to estimate the areal proportions of groundwater resources used for domestic drinking wells that have water-quality constituents present at high concentrations (above the benchmark), moderate concentrations (between one-half of the benchmark and the benchmark for inorganic constituents, or between one-tenth of the benchmark and the benchmark for organic constituents), and low concentrations (less than one-half or one-tenth the benchmark for inorganic and organic constituents, respectively). This method provides statistically representative results at the study-area scale and permits comparisons to other GAMA-PBP study areas. In the study area, inorganic constituents in groundwater were greater than regulatory benchmarks (U.S. Environmental Protection Agency [EPA] or State of California maximum contaminant levels [MCLs]) for public drinking-water quality in 24 percent of domestic groundwater resources. The inorganic constituents present at concentrations greater than MCLs for drinking water were nitrate (as nitrogen), barium, chromium, and selenium. Total dissolved solids (TDS) or manganese were present at concentrations greater than the secondary maximum contaminant levels (SMCLs) that the State of California uses as aesthetic-based benchmarks in 48 percent of domestic groundwater resources. No volatile organic compounds or pesticide constituents were present at concentrations greater than regulatory benchmarks. Total coliform bacteria and enterococci were detected in 4 percent of domestic groundwater resources. Per- and polyfluoroalkyl substances (PFAS) were detected in 19 percent of domestic groundwater resources, and 10 percent had concentrations greater than recently enacted (April 2024) EPA MCLs.

Physical and chemical factors from natural and anthropogenic sources that could affect the groundwater quality were evaluated using results from statistical testing of associations between constituent concentrations and potential explanatory variables. In this study, relevant physical factors include well construction characteristics, groundwater age, site proximity to groundwater recharge or discharge zones, and potential sources of contamination. Relevant chemical factors include the initial chemistry of the recharge water, the mineralogy of the aquifer sediments, and the subsequent shifts in chemistry as biologic and geologic reactions alter groundwater in the subsurface.

Nitrate concentrations were correlated to agricultural land use, distance from the boundary of the Gilroy-Hollister Valley groundwater basin, and the proportion of modern (post-1950s) water captured by the well. Denitrification under anoxic redox conditions can mitigate some nitrate derived from fertilizer application. Total dissolved solids primarily were derived from water-rock interactions with soils and aquifer materials in the study area, but there were high concentrations where agricultural practices contributed additional TDS. Mineralogy of aquifer sediments and rocks also affect barium, selenium, boron, and chromium concentrations in the Gilroy-Hollister Valley groundwater basin. PFAS were positively correlated with urban land use and the proportion of modern water captured by the well.

Released December 01, 2025 07:18 EST

2025, Data Report 1217

Brian G. Prochazka, Peter S. Coates, Cameron L. Aldridge, Michael S. O’Donnell, David R. Edmunds, Adrian P. Monroe, Steve E. Hanser, Lief A. Wiechman, Michael P. Chenaille

Greater sage-grouse (Centrocercus urophasianus; hereafter sage-grouse) are at the center of State and national land-use policies largely because of their unique life-history traits as an ecological indicator for the health of sagebrush ecosystems. This updated population trend analysis provides State and Federal land and wildlife managers with the best available science to help guide management and conservation plans aimed at benefiting sage-grouse populations and the ecosystems they inhabit. This analysis relied on previously published population trend modeling methodology from Coates and others (2021, 2022a) and incorporates population lek count data for 1960–2024. Included in this report are methodological updates to lek count data aggregation, state-space model forecasting, and targeted annual warning system signals, which are detailed under individual Modification sections. State-space models estimated a 2.9-percent average annual decline in sage-grouse populations between 1966 and 2021 (Period 1, six population oscillations) across their geographical range. The average annual decline among climate clusters for the same number of oscillations ranged between 2.2 and 3.4 percent. Cumulative declines were 41.2, 64.1, and 78.8 percent range-wide in Period 5 (19 years), Period 3 (35 years), and Period 1 (55 years), respectively.

Released November 27, 2025 09:18 EST

2025, Oecologia (208)

Anthony M. Pagano, Karyn D. Rode, Kerry L. Nicholson, William B. Leacock, Craig A. Stricker, Charles T. Robbins

Omnivores often face tradeoffs between selecting for spatially dispersed energy-dense vertebrate prey versus densely distributed herbivorous resources that have limited energetic value per unit intake. Arctic grizzly bears (Ursus arctos) are large omnivores within a resource-limited ecosystem that are known to exhibit smaller body masses and occur at lower densities than grizzly bears in other regions of North America. We evaluated the energy balance of Arctic grizzly bears during a portion of the fall hyperphagic period in two ecologically differing regions on Alaska’s northern Arctic coast by monitoring mass change, food intake, activity, and energy expenditure of 12 individuals over 17–22 days. Bears in coastal areas were more carnivorous than bears in the foothills that were predominantly herbivorous and frugivorous. Carnivory was associated with greater movement, body fat, and energy expenditure and two of four carnivorous bears lost mass. Overall, the mean body fat of the bears in this study was 34% lower than other grizzly bear populations in North America in the fall. Furthermore, the bears in this study exhibited relatively small changes in body mass (x̄= 3%, range =−2 to 11%) that were 60% lower than other grizzly bear populations which typically gain substantial mass in the fall in preparation for denning. Our results, while representing a snapshot from a small number of bears during the fall hyperphagic period, are consistent with previous studies and indicate limited availability of energy-dense food resources during this time for grizzly bears in this region of the Arctic.

Released November 26, 2025 12:05 EST

2025, Open-File Report 2025-1054

Andrea C. O'Neill, Sharon F. Batiste, Daniel D. Buscombe, Joseph Burgess, Kara S. Doran, Ann E. Gibbs, Rachel E. Henderson, Julia L. Heslin, Catherine N. Janda, Mark A. Lundine, Joseph F. Terrano, Jonathan A. Warrick, Kathryn M. Weber

Satellite-derived shorelines (SDS) have the potential to help researchers answer critical coastal science questions and support work to predict coastal change by filling in the spatial and temporal gaps present in current field-based and remote-sensing data collection methods. The U.S. Geological Survey conducted comparison analyses of traditionally sourced shorelines and SDS in diverse coastal landscapes to determine how SDS could be used in ongoing and future work across varied coastal environments and provided some initial findings that could be used for implementation. Using CoastSeg, a browser-based program for SDS detection and mapping, SDS for the period 1984–2023 for multiple locations across the United States were compared to shoreline positions from traditionally sourced shoreline data. In this report, the authors present these comparisons alongside lessons learned and challenges encountered when building SDS workflows in different coastal locations. Results show that individual SDS have larger uncertainty and yet produced similar linear trends to sparser, traditionally sourced shoreline data; because SDS methods provide orders of magnitude more data than traditional shoreline-detection methods, they can be used to evaluate shoreline behaviors. Refining average scalar slopes used in tidal corrections did not result in substantial decreases in uncertainty. Using lessons from this work to outline needs for regional implementation, initial setup time would be considerable, being on the order of weeks. However, once complete, shoreline detections and analyses are fast (on the order of minutes to hours) and achievable using a desktop computer.

Released November 26, 2025 11:01 EST

2025, Advances in Water Resources (206)

Gábor Fleit, Marian Muste, Sándor Baranya, Dongsu Kim, Amanda Whaling, Tate McAlpin, Hojun You

This paper introduces the use of acoustic Doppler current profiler (ADCP) measurements as input for the Acoustic Mapping Velocimetry (AMV) method, a technique for characterizing the dynamics of riverine bedforms. The performance of this new approach, ADCP-AMV, is compared with input from a multibeam echosounder through a field study conducted on the Mississippi River (USA). A virtual ADCP tool has been created to support the ADCP-AMV measurements with optimal data density predictions. To the authors’ knowledge, this is the first time ADCP measurements have been used in conjunction with the AMV dune-tracking method. Subsequently, the paper discusses the coupling of ADCP-AMV measurements with ancillary data extracted from the ADCP. These ancillary data are processed using previously developed protocols to characterize hydrodynamics and the suspended sediment distribution in the water column. This paper emphasizes the capability of ADCPs to characterize open-channel river hydromorphodynamic parameters with high spatiotemporal resolution. Recommendations to accurately and efficiently acquire these multi-variable measurements and derived datasets are discussed.

Released November 26, 2025 10:21 EST

2025, Freshwater Biology (70)

Freya Elizabeth Rowland, Michael J. Vanni, Nicole M. Hayes, Clifford E. Kraft

1. Although nitrogen and phosphorus deficiency of algal blooms have been the focus of substantial attention, organic nutrients can limit algal growth in aquatic systems. Growing evidence indicates thiamine (vitamin B₁) can influence the community of primary producers in marine systems, but comparatively little is known about the effect of thiamine on freshwater algal productivity.
2. We conducted 106 nutrient deficiency experiments with water from 39 Ohio lakes of varying trophic status during the growing seasons (April–October) of 2008–2009. Specifically, we tested the response of phytoplankton biomass (as chlorophyll a, chl-a) relative to controls to added nitrogen (N), phosphorus (P), thiamine (Th), or combinations of N + P and N + P + Th. Next, we compared the chl-a growth response of treatment/control to published thresholds based on frequentist approaches and compared the conclusions with Bayesian model results that focused on probability of a response.
3. Although N + P addition was consistently associated with the largest chl-a response, we found evidence of a thiamine influence on phytoplankton growth in some experiments. The Bayesian approach suggested thiamine may become more limiting as the growing season progresses. By late in the growing season, there was an 85% probability of a positive algal growth response to thiamine addition.
4. Understanding the role of thiamine or other overlooked nutrients is not likely to alter the prevailing understanding of nutrient deficiency in freshwater ecosystems. However, we present evidence that freshwater phytoplankton may experience thiamine deficiency and suggest limnologists consider thiamine when exploring resource deficiencies.

Released November 26, 2025 09:23 EST

2025, Economic Geology

Andreas Audétat, Jia Chang, Sean Patrick Gaynor

The depth level at which porphyry Cu–forming magmas fractionated and exsolved mineralizing fluids is actively debated. In the classic model, extensive magma fractionation occurs in large, upper crustal magma chambers, and concomitant fluid exsolution leads to forceful expulsion of residual magmas in the form of porphyry dikes, stocks, and breccia pipes, which subsequently serve as pathways for the mineralizing fluids. In contrast, some recent studies highlighting the role of deep crustal magma fractionation in the production of fertile magmas essentially deny the existence of upper crustal magma chambers at the time of mineralization. To address this, we conducted a detailed thermobarometric investigation of 13 intermediate to felsic, porphyritic intrusive rocks related to porphyry-skarn Cu mineralization at Santa Rita and Hanover-Fierro, New Mexico, United States, representing two premineralization magmas (61–60 Ma), seven synmineralization magmas (60–58 Ma), and four late- to postmineralization magmas (58–57 Ma).

For each sample, the pressure of last magma crystallization before final magma ascent to the current exposure level was reconstructed based on Al-in-hornblende barometry of small hornblende inclusions trapped within quartz phenocrysts and through titanium-in-quartz (TitaniQ) thermobarometry of the host quartz phenocrysts themselves. Since quartz is one of the last crystallizing magmatic minerals, and no significant phenocryst growth could have occurred in small dikes and stocks after final magma emplacement, quartz phenocrysts and their contained hornblende inclusions record the depth of last magma crystallization before final magma ascent. When present, hornblende phenocrysts and hornblende inclusions within other major phenocrysts were also analyzed. Both quartz and hornblende barometers return consistent average pressures of 3.2 ± 0.4 kbar for the entire suite of pre- to postmineralization magmas, corresponding to depths of 11 to 14 km. The synmineralization magmas return even more consistent average pressures of 3.1 ± 0.2 kbar, corresponding to a depth of 12 ± 1 km.

The volume of the mineralizing porphyry dikes and stocks at the emplacement level is far too small to have provided all the fluids and metals required to form the observed ore deposits. Therefore, the majority of the ore-forming fluids must have originated from the magmas that crystallized at 12 ± 1 km depth. The ore deposits, conversely, formed at ~5-km paleodepth. This implies that most of the mineralizing fluids traveled an average vertical distance of ~7 km from their magmatic source to the eventual site of ore precipitation. The relatively unaltered nature and low veining degree of deeper parts of mineralized porphyry dikes and stocks suggest that the fluid transport through these intrusive bodies occurred mostly at near-solidus conditions by means of fluid percolation along grain boundaries.

In summary, our results suggest that (1) a large, upper crustal pluton exists ~7 km beneath the Santa Rita and Hanover-Fierro deposits; (2) abundant phenocryst crystallization occurred at this depth level; and (3) this pluton was the main source for the exsolution of ore-forming fluids. However, the investigated rocks have elevated whole-rock Sr/Y ratios, indicating magma fractionation at deep crustal levels. As a result, our preferred model is a combination of the two end-member models introduced above, with most magma fractionation having occurred in the deep crust and with residual, intermediate to felsic melts having ascended and accumulated at 11 to 14 km paleodepth, where they continued to crystallize with comparatively little crystal-liquid separation, before some of these magmas ascended further to shallow levels and quenched to porphyries.

Released November 26, 2025 08:25 EST

2025, Scientific Investigations Report 2025-5099

Brock J.W. Kamrath, Courtney N. Lauderback, Jennifer C. Murphy

This report describes a cooperative study by the U.S. Geological Survey and Missouri Department of Natural Resources that evaluated temporal changes in total nitrogen (TN) and total phosphorus (TP) concentrations in the Lower Grand River hydrologic unit. The study focused on trends since 2010, when the basin was designated as a priority drainage basin of the Mississippi River Basin Healthy Watersheds Initiative (MRBI). At three local drainage basins within the Lower Grand hydrological unit (MRBI sites), stream nutrient trends were evaluated using flow-adjusted (FA) TN and TP concentrations for water years 2011 through 2023. FATN concentration trends were not statistically significant for any MRBI site. One site (site 2) showed a statistically significant increasing trend in FATP concentration, indicating a possible increase in phosphorus sources in parts of the basin. Overall, streamflow variability appeared to be the dominant factor affecting nutrient concentrations at MRBI sites. At five regional drainage basins, including the Grand River and nearby rivers with data from 1994 through 2023 (long-term sites), annual flow-normalized (FN) TN and TP concentrations were evaluated for trends before (water years 2000–10) and during (water years 2010–23) the MRBI. For water years 2010 through 2023, annual FNTN and FNTP concentrations decreased in the Grand River, as well as in the Nodaway and Chariton Rivers, which were not targeted by the MRBI. The Grand River (site 9) reversed from increasing to decreasing FNTP concentrations after 2010, with a 26-percent reduction. Annual FNTN and FNTP concentrations also decreased at the Missouri River sites. While nutrient reductions in the Grand River may reflect the effects of implemented conservation practices, similar trends in nearby, nontargeted rivers and the absence of strong decreasing trends at MRBI sites suggest that broader regional factors, instead of or in addition to MRBI efforts, may have contributed to nutrient reductions in the Grand River.

Released November 25, 2025 15:50 EST

2025, Scientific Investigations Report 2025-5062

Kristina Gutchess, Natasha Scavotto, Amanda Dondero, Joshua Woda, Neil Terry, Kirk Smith, John Williams

The U.S. Geological Survey, in cooperation with the New York State Department of Transportation, evaluated the effects of different deicing salt application rates on surface water, groundwater, and highway runoff quality near State highways in northern New York. Three reaches of State highways were tested with different deicing treatments between October 2019 and November 2022: a salt-sand mixture (Treatment A), a salt mixture applied at a lower rate (Treatment B), and a control mixture consistent with typical deicing salt amounts and application rates. Data on pavement conditions and the quality of surface water, highway runoff, and groundwater were collected. Surface electromagnetic data were also collected. Surface-water and groundwater quality downgradient from the State highways were compared with water quality at upgradient locations. The percentage of snow or ice coverage was used to evaluate the effectiveness of the salt applications.

This report provides an overview of the transport of deicing salt. The Treatment B watershed had deicing mixture applied more frequently than other highway reaches, which caused it to have the highest annual total chloride application. Despite differences in chloride application, flow-weighted mean chloride concentrations in highway runoff were comparable across treatments. Chloride concentrations were elevated in surface water and groundwater downgradient from highways relative to chloride concentrations upgradient from highways. A chloride mass balance, calculated for one treatment watershed, indicated that groundwater affected by legacy deicing practices may be contributing additional chloride to surface water. Spatial patterns from electromagnetic surveys show a shallow saline plume alongside the highway in that area.

Differences in winter severity and pavement-surface conditions drove deicing salt applications in the treatment areas. This study found that several factors affect chloride loads in the watersheds, including variable winter conditions, adaptive snow and ice management, legacy management practices, and area-specific aquifer and groundwater conditions.

Released November 25, 2025 09:50 EST

2025, Molecular Ecology

Jesse B. Parker, Sean Hoban, Laura Thompson, Scott E. Schlarbaum

The central–marginal hypothesis (CMH) predicts reduced genetic diversity and increased differentiation in range-edge populations due to ecological marginality and limited gene flow. Deviations from this pattern, however, can result from historical demographic processes, variation in reproductive strategies or interspecific hybridization. The genus Quercus, known for hybridization and long-distance pollination, offers an excellent model to examine the spatial patterns of genetic diversity, structure and introgression across species distributions. Here, we investigate these dynamics in Quercus bicolor Willd., a widespread eastern North American oak. Using RADseq, we genotyped 142 individuals from 12 sites at the fragmented trailing range edge and nine sites from the range core. To detect introgression, we incorporated reference data from six sympatric white oak species. We reveal extensive introgression, particularly from Q. lyrata Walt., in nearly all southern edge populations, but none in core populations despite sympatry with closely related congeners. Southern populations also showed increased genetic structure and differentiation, but not reduced diversity or increased inbreeding, even when only examining non-admixed individuals. Regression analyses reveal relationships between introgressed ancestry and heterozygosity, inbreeding and differentiation, indicating that introgression may buffer range-edge populations against genetic erosion by introducing novel alleles. Hindcast, current and forecast ecological niche models demonstrate temporally changing degrees of overlap between the geographic range of Q. lyrata and Q. bicolor and suggest higher hybridization potential in the future. These findings offer mixed support for the CMH while underscoring the evolutionary relevance of introgression in shaping genetic landscapes at range margins with significant implications for conservation.

Released November 24, 2025 13:01 EST

2025, Scientific Investigations Report 2025-5088

Scott A. Olson

This report provides estimates of flood discharge at selected annual exceedance probabilities (AEPs) for streamgages in and adjacent to Vermont and equations for estimating flood discharges at AEPs of 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent (recurrence intervals of 2-, 5-, 10-, 25-, 50-, 100-, and 500-years, respectively) for ungaged, unregulated, rural streams in Vermont with drainage areas between 0.47 and 851 square miles. The equations were developed using generalized least-squares regression and flood-frequency and drainage-basin characteristics from 156 streamgages. Flood-frequency analyses were completed using data through the 2023 water year. The drainage-basin characteristics used as explanatory variables in the regression equations are drainage area, percentage of wetland area, and basin-wide mean of the average annual precipitation. The average standard errors of prediction used to estimate flood discharges at the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent AEP with these equations are 34.9, 37.1, 38.2, 41.6, 43.8, 46.0, 49.1, and 53.2 percent, respectively.

Flood discharges at selected AEPs for streamgages were computed using the Expected Moments Algorithm. Techniques used to adjust an AEP discharge computed from a streamgage record with results from the regression equations and to estimate flood discharge at a selected AEP for an ungaged site upstream or downstream from a streamgage using a drainage-area adjustment are both described. The final regression equations and the flood-discharge frequency data used in this study will be available in StreamStats. StreamStats is an internet-based application that provides automated regression-equation solutions for user-selected sites on streams.

Released November 23, 2025 08:22 EST

2025, Ecology and Evolution (15)

Gavin G. Cotterill, Douglas A. Keinath, Tabitha A. Graves

Occupancy models have become increasingly popular for species monitoring and assessment, in part, because detection/non-detection data are readily obtained using a variety of methods. Multispecies occupancy models (MSOMs) can yield more accurate parameter estimates than single-species models (SSOMs) with less data through their hierarchical structure, making MSOMs an attractive option when species are hard to detect or when data collection is constrained, leading to sparse datasets. Such constraints may arise from limited sampling resources, but also occur in rare species monitoring or where preliminary results are desired to inform adaptive management. Further, experimental habitat treatments often impose spatial constraints on sampling based on the scale of their implementation. Whether a MSOM outperforms SSOMs depends on the volume of data, characteristics of the ecological community, research goals of a study and how these factors align with modeling assumptions. We performed a simulation study of hypothetical pollinator communities under varying sampling intensities for scenarios in which experimental habitat treatments produced different community-level effects. We fit occupancy models to simulated datasets and assessed model performance. At lower sampling intensities (< 20 spatial replicates and < 4 temporal replicates), MSOM community-level treatment effect estimates were biased. Even at twice this sampling intensity, SSOMs yielded more accurate species-specific effect estimates in treatment effect scenarios with high variance. In some cases, MSOMs can pull species in the tails of distributions too far toward the community mean effect, which risks incorrect conclusions concerning whether treatments help or harm individual species. When quantifying species-specific effects is the main objective, particularly for rarely observed species, SSOMs are more robust to outliers across a range of community response scenarios. Researchers can use this information to inform study design, guide simulation studies and decide whether the higher precision of MSOMs outweighs risks of improperly estimated effects for some species.

Released November 22, 2025 08:02 EST

2025, Journal of Great Lakes Research

Ralph W. Tingley III, Darryl W. Hondorp, Benjamin A. Turschak, Steven A. Pothoven, Amanda Susanne Ackiss, Jory Jonas, William W. Fetzer, Benjamin Scott Leonhardt, Andrew Edgar Honsey, Jeff Elliott, Lindsie Ann Egedy, Cory Brant, Lynn Benes, Kendra Kozlauskos, Renee Renauer-Bova, Ann J. Ropp

Characterizing fish movements is required for understanding habitat use, energy flow, and trophic structure and can inform fisheries management. Drowned river mouth (DRM) lakes are productive inland habitats in the Laurentian Great Lakes basin used by migratory fishes. Despite recognition of their ecological connections to the Great Lakes, the value of DRM lakes as seasonal habitats is not well understood for many fishes. One such species, cisco Coregonus artedi, has recently expanded in Lake Michigan from near extirpation to higher relative abundances in the northeastern portion of the lake. Cisco are recreationally harvested in some DRM lakes during winter, but little is known about cisco movement patterns and ecology. In winter 2022 and 2023, we collected cisco from three DRM lakes along the eastern shores of Lake Michigan (Lake Charlevoix, Portage Lake, Muskegon Lake) to characterize genetics, morphometrics, and diets. We also implanted telemetry tags in 20 cisco collected in Lake Charlevoix to examine movement patterns and determine DRM lake residency (i.e., seasonal vs. year-round). We found no consistent genetic or morphometric differentiation across DRM lakes, suggesting that recolonization began from a single stock. Fish were the only diet item found in cisco guts collected during winter months. Movement patterns from Lake Charlevoix indicated strong spawning site fidelity to Grand Traverse Bay as well as non-spawning site fidelity. However, given the presence of cisco in southern DRM lakes and some site-specific differences in morphometrics, managers could benefit from further research to determine whether spawning occurs in southern Lake Michigan.

Released November 21, 2025 11:55 EST

2025, Fact Sheet 2025-3052

Christopher J. Schenk, Marilyn E. Tennyson, Tracey J. Mercier, Phuong A. Le, Andrea D. Cicero, Ronald M. Drake II, Sarah E. Gelman, Jane S. Hearon, Benjamin G. Johnson, Jenny H. Lagesse, Heidi M. Leathers-Miller

Using a geology-based assessment methodology, the U.S. Geological Survey estimated undiscovered, technically recoverable mean conventional resources of 67 million barrels of oil and 56 billion cubic feet of gas in the Santa Maria Basin Province of California.  

Released November 21, 2025 09:23 EST

2025, Seismological Research Letters

Ruth A. Harris, Michael Barall, Grace Alexandra Parker, Evan Tyler Hirakawa

he Hayward, Calaveras, and Rodgers Creek faults in the San Francisco Bay region of California have a high probability of producing a large earthquake in the next decades. Although these faults creep, the creep is insufficient to keep up with their relatively rapid slip rates on their deepest sections, so they have been storing tectonic strain since their last large earthquakes, with the Hayward’s and Rodgers Creek’s more than 150 yr ago. We do not know what the next large Hayward–Calaveras–Rodgers Creek earthquakes will look like or how strongly they will shake the San Francisco Bay region. Harris et al. (2021) used the 3D dynamic (spontaneous) rupture method to simulate large earthquakes on these creeping faults. In this article, we examine the resulting simulated long‐period (T > 1 s) ground shaking from 0 to 50 km distance, for earthquakes nucleating on the Hayward fault and earthquakes nucleating on the Rodgers Creek fault. We compare these simulated long‐period ground motions with the Boore et al. (2014) well‐established empirically based ground‐motion model suitable for the slowest material velocity in our 3D velocity structure. We find that the simulated long‐period ground motions from the creeping‐fault earthquake scenarios produce a reasonable agreement with the empirical expectations if frictional cohesion is included only where it is appropriate.

Released November 20, 2025 12:37 EST

2025, JGR Earth Surface (130)

Rose Elizabeth Palermo, Jennifer L. Miselis, Daniel J. Ciarletta, Emily A Wei

Barrier island resilience to climate impacts depends on sediment redistribution between the subaqueous shoreface and subaerial barrier during sea-level rise and storms. However, autogenic interactions between the upper and lower shoreface and their influence on the subaerial barrier are poorly characterized. Here, we explore the influences of various shoreface components on barrier morphology using a model of barrier and shoreface evolution under sea-level rise, the Articulated Barrier Shoreface (ABSF) Model. This reduced-complexity model divides the shoreface into upper and lower shoreface panels that respond independently to sea-level rise and deviations from the equilibrium slope. We couple the ABSF with the Lorenzo-Trueba & Ashton, 2014, https://doi.org/10.1002/2013jf002941 model (LTA), a barrier island evolution model driven by overwash and sea-level rise. Through this coupled framework, we examine the influences of upper and lower shoreface slopes, their respective depths, and sensitivity to wave climate on long-term barrier evolution. Results show that the relative depths of the upper and lower shoreface toes influence barrier response to rising seas, alongside overwash flux and closure depth. Notably, the lower shoreface response to sea-level change lags that of the upper shoreface over decades, diminishing the resilience of the barrier over centennial timescales by slowing the overall barrier response. In fact, the ABSF model predicts barriers will drown faster and more than predicted with a linear shoreface. Results highlight the shoreface as an important sediment reservoir for barrier islands and that differences in upper and lower shoreface responses can reduce barrier resilience to sea-level rise due to limited lower shoreface sediment accessibility.

Released November 20, 2025 11:30 EST

2025, Professional Paper 1890-I

Carmen Jaimes-Viera, Amiel Nieto-Torres, Ana Lillian Martin Del Pozzo, Aurelie Germa, Chuck Connor, Michael H. Ort, Paul Layer, Jeff Benowitz

An analysis of 1,375 volcanoes in the Michoacán-Guanajuato (1,148 volcanoes in a 26,200 square-kilometer area) and Sierra Chichinautzin (227 volcanoes in a 3,500 square-kilometer area) volcanic fields in central Mexico identified patterns in the spatial and temporal distribution of past eruptions. A cluster agglomerative hierarchical method and kernel analysis confirmed that the Michoacán-Guanajuato volcanic field comprises four volcanic fields (Valle de Santiago, Uruapan, Apatzingán, and Pátzcuaro volcanic fields) controlled by different fault systems, indicating that it is not a single volcanic field but rather a group of volcanic fields (a “superfield”), each of which has distinct characteristics.

In the Sierra Chichinautzin volcanic field, well-constrained isotopic ages were used to build a model of how the spatial distribution of the eruptions has changed over time. Two new 40Ar/39Ar ages from a locally recognized volcanic feature near the town of El Cantil, herein called El Cantil volcano (1,537±17 kilo-annum [ka]) and the volcanic feature at Cerro el Elefante (herein called El Elefante dome) (1,485±92 ka) belong to the oldest volcanic group identified in the Sierra Chichinautzin volcanic field, confirming the timing of the beginning of monogenetic volcanism in the region. Based on the volcanic groups identified in the Sierra Chichinautzin volcanic field, the youngest volcanism (less than 35 ka) is found only in the central-western sector of the field. Principal component analysis determined the directional trends of feeder dikes only for vents <10 ka in the Sierra Chichinautzin volcanic field. Possible magma migration paths through the crust were identified using seismic data from both volcanic fields using an earthquake catalog from 1973 to 2023, which includes 9,016 earthquakes in the Michoacán-Guanajuato volcanic field and 841 in the Sierra Chichinautzin volcanic field. The spatial distribution of the hypocenters does not highlight any trend that could be associated with superficial movement of magma in the Sierra Chichinautzin volcanic field. In the Michoacán-Guanajuato volcanic field, however, eight seismic swarms since 1997 have been detected. These swarms are interpreted to result from ascending magma. Strengthening monitoring systems and reinforcing mitigation measures to address volcanic hazards and risk are important means of preparing for future eruptions in both regions. Analysis such as those herein provide insights as to where an eruption might occur and may help mitigate volcanic hazards.

Released November 20, 2025 10:30 EST

2025, Fact Sheet 2025-3049

Laszlo P. Keszthelyi, Joshua A. Coyan, Lori M. Pigue, Kristen A. Bennett, Travis S.J. Gabriel

Rare earth elements (REEs) are a scarce but vital resource for our modern economies and lifestyles. Since the late 1990s, China has supplied the vast majority of the world’s refined REEs. Increasing global demand has broadened the search for REE deposits to unconventional places, including the Moon. Although most lunar rocks have very low REE concentrations, Apollo samples showed that one type of lunar rock containing potassium (K), REEs, and phosphorus (P)—known by the acronym KREEP—has high concentrations of REEs. Data from orbiting satellites have identified locations where substantial deposits of KREEP are likely. The viability of mining these deposits depends on the evolution of REE economics, the development of the Earth-Moon infrastructure, and the findings from future lunar mineral exploration missions.

Released November 20, 2025 10:03 EST

2025, Oikos

A. Donatelli, Mark Haroldson, Justin G. Clapp, P. Ciucci, Frank T. van Manen

Plasticity of diel activity rhythms may be a key element for adaptations of wildlife populations to changing environmental conditions. In the last decades, grizzly bears Ursus arctos in the Greater Yellowstone Ecosystem (GYE) have experienced notable environmental fluctuations, including changes in availability of food sources and severe droughts. Although substantial research has been conducted on grizzly bear diets, space use, and demographic parameters, studies on factors that may influence their diel activity patterns are lacking. We investigated diel activity of grizzly bears in the GYE as a function of anthropogenic landscape modification, maximum daily ambient temperature, drought severity, and bear density. Specifically, we used accelerometry readings of 169 bears (39 females, 130 males) from 2009 to 2022 to compute three complementary activity measures, hourly intensity of activity, daily active minutes, and active bout length, each used as a response variable within a Bayesian modeling framework. Grizzly bears generally exhibited bimodal diel activity, with crepuscular peaks and slight variations across seasons. Females with young (i.e. cubs or yearlings) were an exception, with more pronounced diurnal activity patterns, possibly as a strategy to avoid infanticide by dominant males. Landscape modification and maximum ambient temperature were the factors most strongly associated with activity patterns of grizzly bears, with greater nocturnality observed in lone females and males as these factors increased. Females with young were comparatively less affected. The GYE is changing because of increasing land development, human recreation pressures, and effects of climate change. Given their greater diurnal activity compared with other cohorts, female grizzly bears with dependent offspring may be more constrained in their ability to modify activity patterns. Our findings add to a growing body of research emphasizing the importance of the temporal dimension of wildlife behavior as a critical factor in assessing species adaptability and vulnerability in a changing world.

Released November 20, 2025 09:11 EST

2025, Journal of Hydrology: Regional Studies (62)

Eleanor E. Fahrney, David C. Mays, Connor P. Newman

Released November 20, 2025 08:55 EST

2025, Open-File Report 2025-1049

Faith A. Fitzpatrick, Shelby P. Sterner, James D. Blount, Jana S. Stewart

Urban stream rehabilitation plans can benefit from knowledge of the landscape setting and vegetative communities that were adjacent to streams prior to urbanization. Downstream to upstream connections of these characteristics can be relevant for native migratory fish species that have a range of preferred spawning habitats. Based on a need for more quantitative data on these potential connections, the U.S. Geological Survey assembled geomorphic characteristics, surficial geology, and pre-Euro-American settlement vegetation for 333 kilometers of stream segments in the Kinnickinnic River and Menomonee River subbasins of the Milwaukee River, Wisconsin. Channel slopes ranged from less than 0.3 percent to greater than 2 percent, covering at least two channel morphology and bedform types spanning low-energy irregular and pool-riffle complexes. Postglacial surficial geology ranged from coarse-grained outwash sand and gravel to lacustrine silt and clay, allowing for a range of stream substrate sizes. Presettlement riparian vegetation was mainly forest, including forested uplands, forested lowlands, and to a lesser extent, conifer-dominated wetlands in headwaters. This resulting framework of geomorphic habitat response units can be used for habitat rehabilitation projects for migratory native fish in other urban Great Lakes tributaries.

Released November 20, 2025 08:10 EST

2025, Electronics (14)

Nusrat Zahan, Sidike Paheding, Ashraf Saleem, Timothy C. Havens, Peter C. Esselman

In recent years, Single-Image Super-Resolution (SISR) has gained significant attention in the geoscience and remote sensing community for its potential to improve the resolution of low-quality underwater imagery. This paper introduces MIMAR-Net (Multiscale Inception-based Manhattan Attention Residual Network), a new deep learning architecture designed to increase the spatial resolution of input color images. MIMAR-Net integrates a multiscale inception module, cascaded residue learning, and advanced attention mechanisms, such as the MaSA layer, to capture both local and global contextual information effectively. By utilizing multiscale processing and advanced attention strategies, MIMAR-Net allows us to handle the complexities of underwater environments with precision and robustness. We evaluate the model on three popular underwater image datasets, namely UFO-120, USR-248, and EUVP, and perform extensive comparisons against state-of-the-art methods. Experimental results demonstrate that MIMAR-Net consistently outperforms existing approaches, achieving superior qualitative and quantitative improvements in image quality, making it a reliable solution for underwater image enhancement in various challenging scenarios.

Released November 19, 2025 11:56 EST

2025, Scientific Investigations Report 2025-5017

Libby M. Wildermuth, Jeffrey R. Kennedy, Jacob L. Conrad

Managed aquifer recharge is a widespread practice for storing water in the subsurface as groundwater. At a managed aquifer recharge facility in southern Arizona, groundwater-level and repeat microgravity data were collected to monitor aquifer response. These data were used to inform parameter identification for an unsaturated-zone flow model used to simulate the recharge process. The facility, the Southeast Houghton Artificial Recharge Project (SHARP), consists of 3 surface basins (about 27,600 square meters [6.8 acres] total surface area) where recycled water is distributed in recharge cycles lasting several months, with dry periods in between. During the study period, December 2020–December 2022, Tucson Water (the City of Tucson’s water utility) reported 6.56×10⁶ cubic meters of water (5,320 acre-feet) recharged.

Monitoring included groundwater-level observations at 3 monitoring wells and repeat microgravity measurements at as many as 22 locations (some stations were destroyed between surveys). Six gravity surveys were carried out using absolute- and relative-gravity meters. Large gravity increases, more than 250 microgals, were observed during the first repeat survey, 3.5 months after the start of recharge, but only in the immediate vicinity of the recharge basins. Data show that water moved downward to the water table, and storage changes in the unsaturated zone away from the facility were likely minimal. Gravity decreased at stations more than 1 kilometer from the facility, consistent with regional groundwater-level changes. Groundwater-level increases in wells adjacent to the recharge basins began 2 months after the second repeat gravity survey, and 5.5 months after recharge began.

Unsaturated-zone flow modeling was carried out using software that simulates water movement and parameter estimation. Model calibration was carried out by minimizing an objective function calculated from the differences between simulated and observed groundwater levels, and between simulated and observed repeat microgravity data. Including repeat microgravity data in the objective function reduced the uncertainty in estimated parameter values for saturated hydraulic conductivity and saturated water content. Modeling indicated that the unsaturated zone between the recharge basins and the water table does not become saturated even after 685 days of simulated infiltration. This gradual wetting may account for increasing infiltration rates over time, as hydraulic conductivity increases with increasing water content. Unsaturated-zone water content decreased rapidly between recharge cycles. Model-simulated groundwater mounding extended about 1 kilometer from the center of SHARP after the 685-day period following the onset of recharge.

Released November 19, 2025 08:15 EST

2025, Hydrology and Earth System Sciences (29) 6445-6460

Phillip J. Goodling, Jennifer H. Fair, Amrita Gupta, Jeffrey D. Walker, Todd Dubreuil, Michael J. Hayden, Benjamin H. Letcher

Despite their ubiquity and importance as freshwater habitat, small headwater streams are under-monitored by existing stream gage networks. To address this gap, we describe a low-cost, non-contact, and low-effort method that enables organizations to monitor relative streamflow dynamics in small headwater streams. The method uses a camera to capture repeat images of the stream from a fixed position. A person then annotates pairs of images, in each case indicating which image has more apparent streamflow or indicating equal flow if no difference is discernible. A deep learning modeling framework called streamflow rank estimation (SRE) is then trained on the annotated image pairs and applied to rank all images from highest to lowest apparent streamflow. From this result a relative hydrograph can be derived. We found that our modeled relative hydrograph dynamics matched the observed hydrograph dynamics well for 11 cameras at 8 streamflow sites in western Massachusetts. Higher performance was observed during the annotation period (median Kendall's Tau rank correlation of 0.75, with a range of 0.6–0.83) than after it (median Kendall's Tau of 0.59, with range 0.34–0.74). We found that annotation performance was generally consistent across the 11 camera sites and 2 individual annotators and was positively correlated with streamflow variability at a site. A scaling simulation determined that model performance improvements were limited after 1000 annotation pairs. Our model's estimates of relative flow, while not equivalent to absolute flow, may still be useful for many applications, such as ecological modeling and calculating event-based hydrological statistics (e.g., the number of out-of-bank floods). We anticipate that this method will be a valuable tool to extend existing stream monitoring networks and provide new insights on dynamic headwater systems.

Released November 18, 2025 12:22 EST

2025, Ecosphere (16)

Bradley J. Butterfield, Daniel Rodolphe Schlaepfer, Robert Al-Chokhachy, Jason B. Dunham, Jeremiah D. Groom, Clint C. Muhlfeld, Christian E. Torgersen, John B. Bradford

Landscapes encompass both aquatic and terrestrial ecosystems that experience the same climate but may respond to climate in divergent ways. For example, the time lag between seasonal dry-down of terrestrial soil moisture and decline in streamflow has important implications for species and ecosystem processes across the aquatic–terrestrial interface. How these lags between aquatic and terrestrial hydrology vary with climate and spatial location within watersheds remains largely unexplored. Here, we examine seasonal patterns of aquatic–terrestrial dry-down across seven watersheds in the northwestern USA, spanning a wide range of aridity. We compared daily streamflow data from USGS gages at watershed outlets with simulated daily soil moisture (1979–2020) from multiple locations within each watershed. In all watersheds, annual dry cycles progressed sequentially through the following features: evapotranspiration, precipitation, shallow soil moisture, deep soil moisture, and finally streamflow. Seasonal streamflow minima lagged behind soil moisture minima for shorter durations in more arid watersheds and drier years. Within watersheds, lag times varied spatially due to interactions between elevation and aridity, with short lags in low-elevation soils near streams in arid watersheds and longer lags in less arid watersheds. Collectively, these results indicate shorter lags between seasonal aquatic and terrestrial dry periods in drier watersheds and years, and show that these tighter linkages are spatially aggregated in drier watersheds. The co-occurrence of seasonally dry conditions in both aquatic and terrestrial systems under increasing aridification is likely to intensify stressors on ecosystems and services. Recognizing these patterns may be critical for predicting ecosystem vulnerabilities and informing adaptation strategies to mitigate the impacts of seasonally dry conditions.

Released November 18, 2025 12:06 EST

2025, Scientific Investigations Report 2025-5055

Moussa Guira, Jonathan P. Traylor, Andrew T. Leaf, Alec R. Weisser

The U.S. Geological Survey has developed a high-resolution inset groundwater-flow model in the Mississippi Delta as part of an interdisciplinary collaboration coordinated by the Mississippi Alluvial Plain project to provide a tool that stakeholders can use to support water-resource management decisions. Groundwater withdrawals from the Mississippi River Valley alluvial (MRVA) aquifer have been vital to support agricultural production in the region, but substantial groundwater-level declines near Shellmound, Mississippi, have caused concerns for long-term sustainability of the aquifer. To better understand the subsurface and try to mitigate the long-term groundwater-level declines, stakeholders have undertaken actions including a Groundwater Transfer and Injection Pilot (GTIP) project using a riverbank filtration-based managed aquifer recharge approach. The pilot project consisted of extracting groundwater near the Tallahatchie River and reinjecting it into the aquifer 3 kilometers west where water levels have substantially declined. A high-resolution airborne electromagnetic (AEM) survey was also completed to collect electrical resistivity data to support the GTIP project and the development of the groundwater model.

The inset groundwater-flow model was developed to (1) integrate the AEM data into the optimal layering configuration of the MRVA aquifer that the available observation data can support through calibration, and (2) assess the potential effect of the GTIP project on the groundwater levels. The AEM data were processed into three different layering configurations leading to the development of model A (18 layers), model B (16 layers), and model C (8 layers), all at a 100- x 100-meter cell spatial resolution using the U.S. Geological Survey modular finite-difference flow model 6 code with Newton-Raphson formulation. The model development process integrated recent advances in modeling, such as the incorporation of AEM data, the use of outputs from the soil-water-balance (SWB) model, and the Aquaculture and Irrigation Water-Use Model, and was facilitated by robust automation using the open-source python packages Modflow-setup and SFRmaker. Using Parameter Estimation ++ Iterative Ensemble Smoother, the three numerical groundwater-flow models (models A, B, and C) were calibrated against a set of observations, which included aquifer groundwater levels, streamflows, stream stage, and aquifer transmissivity. Results indicate that the detailed representation of MRVA aquifer layers in model A produced the best calibrated model by history matching, and the integration of data representing surficial connectivity played a key role in improving groundwater recharge and enhancing the ability of the model to match groundwater levels in the cone of depression. A forecast model simulated the managed aquifer recharge approach, and the results indicated that, given average irrigation and recharge conditions (2010–15), the GTIP project has the potential to induce groundwater-level increases of as much as 3 meters around the injection site, but a sustained increase would require repetition in subsequent years of water transfer at 2022 rates or above.

Released November 18, 2025 09:42 EST

2025, Fact Sheet 2025-3013

Seth D. Burgess

Introduction 

Western California is home to a variety of volcanic rocks. The locations, ages, and chemical compositions of these volcanic rocks help tell part of the fascinating story of California’s plate tectonic evolution over the past 40 million years. These volcanic rocks are a product of multiple tectonic processes, including subduction of divergent and transform plate boundaries beneath continental North America, opening of a slab window, creation and migration of a tectonic triple junction, and the birth and growth of the San Andreas Fault. This fact sheet explains these tectonic processes and discusses their role in shaping the volcanic history of western California over the past 40 million years. By studying the volcanic rock record in western California, geologists are able to piece together how regional volcanism and plate tectonics are linked in space and time. Recognizing this linkage helps scientists to understand possible future volcanism in the region, potential hazards associated with this volcanism, and the impacts these hazards may have on population and infrastructure. The U.S. Geological Survey California Volcano Observatory (CalVO) closely monitors the parts of western California with the greatest potential for volcanism.

Released November 18, 2025 08:55 EST

2025, Journal of Wildlife Diseases (61) 1010-1027

Bryan Daniels, Erik E. Osnas, Megan Boldenow, Robert Gerlach, Christina Ahlstrom, Sarah Coburn, Michael J. Brook, Michael Brubaker, Julian Fischer, David N. Koons, Angela Matz, Marin Murphy, Daniel Rizzolo, Laura Celeste Scott, David R. Sinnett, Jordan M. Thompson, Juliana Lenoch, Mia Kim Torchetti, David E. Stallknecht, Rebecca L. Poulson, Andrew M. Ramey

In 2021–22, clade 2.3.4.4b highly pathogenic avian influenza (HPAI) viruses were introduced by wild birds into North America, leading to geographically widespread disease. In response to HPAI outbreaks throughout late 2021 and early 2022, we recorded observations of sick and dead birds, estimated abundance of carcasses, collected swab and sera samples to detect viruses, and monitored bird nesting on the Yukon-Kuskokwim Delta region of Alaska to document potential effects of disease. Thirty-six reports of sick and dead birds were registered across the region. Nineteen carcasses were opportunistically collected for diagnostic testing, of which 12 were confirmed to be infected with clade 2.3.4.4b HPAI viruses. Carcass abundance estimates from line-distance sampling provided evidence that the most common species of dead birds from the western Yukon-Kuskokwim Delta region were Cackling Goose (Branta hutchinsii minima), Glaucous Gull (Larus hyperboreus), and Black Brant (Branta bernicla nigricans). Only one paired cloacal and oropharyngeal swab sample from a Northern Pintail (Anas acuta) tested positive for clade 2.3.4.4b HPAI virus, out of 464 live-captured duck and goose samples. Of 195 sera samples from waterfowl screened for antibodies reactive to influenza A viruses, antibodies were found in 41–98% of samples collected from Emperor Goose (Anser canagicus), Cackling Goose, Black Brant, and Spectacled Eider (Somateria fischeri). In addition, 15–98% of the same sera samples were reactive to a clade 2.3.4.4b H5 antigen. Fewer Black Brant and Emperor Goose nests were found on long-term study plots during 2022 than in previous years. Collectively, we found that HPAI viruses affected at least seven species of wild birds inhabiting the region during 2022. The full scope of impacts of HPAI at this location during 2022 is unknown, but our data indicate that acute effects to avian population health on the Yukon-Kuskokwim Delta region were likely modest.

Released November 18, 2025 08:33 EST

2025, Climate of the Past (21) 2263-2281

Jacob T. Anderson, Nicolas E. Young, Allie Balter-Kennedy, Karlee Prince, Caleb K. Walcott-George, Brandon L. Graham, Joanna Charton, Jason P. Briner, Joerg M. Shaefer

The lack of geological constraints on past ice-sheet change in marine-based sectors of the Greenland Ice Sheet (GrIS) following the Last Glacial Maximum limits our ability to assess (1) the drivers of ice-sheet change, and (2) the performance of ice-sheet models that are benchmarked against the paleo-record of GrIS change. Here, we provide new in situ ¹⁰Be surface exposure chronologies of ice-sheet margin retreat from the outer Scoresby Sund and Storstrømmen Glacier regions in eastern and northeastern Greenland, respectively. Ice retreated from Rathbone Island, east of Scoresby Sund, by ∼ 14.1 ka, recording some of the earliest documentations of terrestrial deglaciation in Greenland. The mouth of Scoresby Sund deglaciated by ∼ 13.2 ka, and retreated at an average rate of ∼ 43 m yr⁻¹ between 13.2 and 9.7 ka. Storstrømmen Glacier retreated from the outer coast to within ∼ 3 km of the modern ice margin between ∼ 12.7 and 8.6 ka at an average rate of ∼ 28 m yr⁻¹. Retreat then slowed or reached a stillstand as ice retreated ∼ 3 km between ∼ 8.6 ka to the modern ice margin at ∼ 8.0 ka. These retreat rates are consistent with late glacial and Holocene estimates for marine-terminating outlet glaciers across East Greenland, and comparable to modern retreat rates observed at the largest ice streams in northeastern, and northwestern Greenland.

Released November 18, 2025 08:28 EST

2025, Scientific Reports (15)

Shannon L. White, Amanda Higgs, Dewayne Fox

Urban landscapes are often overlooked in conservation planning, allowing human activities to take precedence in ecosystem management. However, even heavily modified environments can support diverse species profiles, but continued expansion of the human footprint could transform these biodiversity hotspots into ecological traps that serve as hidden catalysts for demographic declines. In the backdrop of one of the world’s most urbanized landscapes-New York City, USA—is a federally endangered population of shortnose sturgeon (Acipenser brevirostrum) that has been quietly recovering for several decades despite many demographic threats. Here, we identify a unique behavioral phenotype of shortnose sturgeon that occupies habitats in New York Harbor in late spring and fall, likely using the area to optimize bioenergetic processes. As this study highlights, urbanized environments can be a nexus for cryptic phenotypic diversity which, if overlooked, can disrupt eco-evolutionary processes and contribute to population and species loss.

Released November 17, 2025 12:20 EST

2025, Fact Sheet 2025-3042

Alma Lizette Anides Morales, Laura M. Norman, Thomas J. Mack

The Los Planes watershed of Baja California Sur, Mexico, and its underlying aquifer are experiencing groundwater decline owing to low average annual rainfall (28.1 centimeters per year) and rising water demand from population growth and agricultural activities. This decline in water availability can lead to desertification—a process that changes arable land to desert by degrading soil and vegetation—and can pose serious challenges to livelihoods that depend on the land.

To address these issues, a ranch in the Los Planes watershed has installed many natural infrastructures in dryland streams (NIDS) in channels for soil and water conservation. In 2022, the U.S. Geological Survey (USGS) began working with regional researchers and land managers to investigate the effects of NIDS on natural biological, geochemical, and physical processes and determine the efficacy of NIDS for water augmentation in the Los Planes watershed. The USGS also worked with local academic institutions and nonprofit organizations to create public educational opportunities focused on the area’s hydrogeology. These and other collaborative efforts with the U.S. Water Partnership and Innovaciones Alumbra aim at enhancing water resources in the Baja California Sur region and promoting water security and safeguarding community well-being.

Resume

La cuenca de Los Planes, ubicada en Baja California Sur, México, y su acuífero subyacente, están sufriendo una disminución de las aguas subterráneas debido a la baja precipitación media anual (28.1 centímetros por año) y la alta demanda de agua por parte de una población creciente y la actividad agrícola. Esta disminución de la disponibilidad de agua puede conducir a la desertificación—un proceso que por medio de la degradación del suelo y la vegetación convierte a la tierra cultivable en desierto—representando un serio desafío para los medios de vida de las personas.

Para abordar estos problemas, un rancho en la cuenca de Los Planes ha instalado numerosas obras de Infraestructura Natural en Arroyos de Tierras Áridas (INATS) para conservación del suelo y del agua. En 2022, el Servicio Geológico de los Estados Unidos (USGS, por sus siglas en inglés) comenzó a trabajar con investigadores regionales y gestores de tierras para estudiar los efectos de INATS en los procesos biológicos, geoquímicos y físicos, y determinar su eficacia en el aumento de los recursos hídricos en la cuenca de Los Planes. El USGS se ha asociado con instituciones académicas y organizaciones locales sin fines de lucro para crear oportunidades educativas públicas centradas en la hidrogeología de la zona. Estos y otros esfuerzos colaborativos con la Asociación del Agua de Estados Unidos (U.S. Water Partnership) e Innovaciones Alumbra, tienen como objetivo mejorar el uso de los recursos hídricos en la región de Baja California Sur, promover la seguridad hídrica y proteger el bienestar de la comunidad.

Released November 17, 2025 12:00 EST

2025, Preprint

Brady W. Allred, Sarah E. McCord, Timothy J. Assal, Brandon T. Bestelmeyer, Chad S. Boyd, Alexander C. Brooks, Samantha M. Cady, Michael C. Duniway, Samuel D. Fuhlendorf, Shane A. Green, Georgia R. Harrison, Eric R. Jensen, Emily J. Kachergis, Anna C. Knight, Chloe M. Mattilio, Brian A. Mealor, David E. Naugle, Dylan O’Leary, Peter J. Olsoy, Erika S. Peirce, Jason R. Reinhardt, Robert K. Shriver, Joseph T. Smith, Jason D. Tack, Ashley M. Tanner, Evan P. Tanner, Dirac Twidwell, Nicholas P. Webb, Scott L. Morford

Rangelands are extensive ecosystems, providing important ecosystem services while undergoing continuous change. As a result, improved monitoring technologies can help better characterize vegetation change. Satellite remote sensing has proven effective in this regard, tracking vegetation dynamics at broad and fine scales. We leveraged the spatial, spectral, and temporal resolution of Sentinel-2 satellites to estimate fractional cover and canopy gap across rangelands of the western United States. We produced annual, 10 m spatial resolution estimates of fractional cover and canopy gap size class for years 2018 to 2024. Fractional cover estimates include that of common plant functional types (annual forb and grass, bareground, littler, perennial forb and grass, shrub, tree) and select genera (including invasive annual grass species, pinyon-juniper species, and sagebrush species); canopy gap size classes include gap sizes 25 to 50, 51 to 100, 101 to 200, and greater than 200 cm. We make these data available as Cloud Optimized GeoTIFFs, organized as 75×75 km tiles covering the 17 western states of the United States.

Released November 17, 2025 09:31 EST

2025, North American Journal of Fisheries Management

Brock Huntsman, Matthew J. Young, Kai Palenscar, Kurt E. Anderson, William Ota, Jordan Mae-Jean Buxton, Justin K. Clause, Danielle L. Palm, Jeff Lee Gronemyer, Brett Mills, Kerwin Russell, Rebecca Christensen

Released November 17, 2025 09:05 EST

2025, Remote Sensing (17)

Mahesh Shrestha, Victoria Mary Scholl, Aparajithan Sampath, Jeffrey Irwin, Travis Kropuenske, Josip Adams, Matthew Alexander Burgess, Lance R Brady

The use of Uncrewed Aerial Systems (UASs) for remote sensing applications has increased significantly in recent years due to their low cost, operational flexibility, and rapid advancements in sensor technologies. In many cases, UAS platforms are considered viable alternatives to conventional satellite and crewed airborne platforms, offering very high spatial, spectral, and temporal resolution data. However, the radiometric quality of UAS-acquired data has not received equivalent attention, particularly with respect to absolute calibration. In this study, we (1) evaluate the absolute radiometric performance of two commonly used UAS sensors: the Headwall Nano-Hyperspec hyperspectral sensor and the MicaSense RedEdge-MX Dual Camera multispectral system; (2) assess the effectiveness of the Empirical Line Method (ELM) in improving the radiometric accuracy of reflectance products generated by these sensors; and (3) investigate the influence of calibration target characteristics—including size, material type, reflectance intensity, and quantity—on the performance of ELM for UAS data. A field campaign was conducted jointly by the U.S. Geological Survey (USGS) Earth Resources Observation and Science (EROS) Center and the USGS National Uncrewed Systems Office (NUSO) from 15 to 18 July 2023, at the USGS EROS Ground Validation Radiometer (GVR) site in Sioux Falls, South Dakota, USA, over a 160 m × 160 m vegetated area. Absolute calibration accuracy was evaluated by comparing UAS sensor-derived reflectance to in situ measurements of the site. Results indicate that the Headwall Nano-Hyperspec and MicaSense sensors underestimated reflectance by approximately 0.05 and 0.015 reflectance units, respectively. While the MicaSense sensor demonstrated better inherent radiometric accuracy, it exhibited saturation over bright targets due to limitations in its automatic gain and exposure settings. Application of the ELM using just two calibration targets reduced discrepancies to within 0.005 reflectance units. Reflectance products generated using various target materials—such as felt, melamine, or commercially available validation targets—showed comparable agreement with in situ measurements when used with the Nano-Hyperspec sensor. Furthermore, increasing the number of calibration targets beyond two did not yield measurable improvements in calibration accuracy. At a flight altitude of 200 ft above ground level (AGL), a target size of 0.6 m × 0.6 m or larger was sufficient to provide pure pixels for ELM implementation, whereas smaller targets (e.g., 0.3 m × 0.3 m) posed challenges in isolating pure pixels. Overall, the standard manufacturer-recommended calibration procedures were insufficient for achieving high radiometric accuracy with the tested sensors, which may restrict their applicability in scenarios requiring greater accuracy and precision. The use of the ELM significantly improved data quality, enhancing the reliability and applicability of UAS-based remote sensing in contexts requiring high precision and accuracy.

Released November 17, 2025 08:22 EST

2025, North American Journal of Fisheries Management

Blake Logan, Mark Pegg, Kirk D. Steffensen, Jonathan J. Spurgeon

Objective

Knowledge of invasive Silver Carp Hypophthalmichthys molitrix population demographics and distributions may inform estimates of efforts necessary to achieve reductions in abundance and identify locations to conduct removal. Although extensively studied in other parts of their invasive range (e.g., Mississippi and Illinois rivers), less is known regarding Silver Carp population demographics in Great Plains rivers and streams. As such, this study characterized Silver Carp population demographics along an invasion gradient in a Great Plains river network containing multiple hydrologically unique river reaches and tributaries.

Methods

Boat and tote-barge electrofishing surveys were conducted within the lower Platte River basin in the spring, summer, and fall of 2022 and 2023. Lapilli otoliths were collected for assessment of age, growth, and annual mortality. Information on sex also was collected. A spatial assessment of differences in population demographics was performed between core and periphery regions of the population.

Results

There were 1,528 Silver Carp collected. A spatial difference in relative abundance was observed and was associated with changes in body condition and total length along the longitudinal gradient of the lower Platte River basin. Silver Carp sex ratios were male-skewed across the lower Platte River basin (1.6:1), particularly in reaches near the invasion front (10.1:1).

Conclusions

Silver Carp population demographics within the lower Platte River basin were aligned with an establishing population characterized by rapid individual growth and skewed sex ratios. Broadscale variation in population characteristics, including growth and size structure, suggested density-dependent processes. Silver Carp occurred throughout the study area, indicating that braided Great Plains streams are susceptible to invasion. This study provided insight into Silver Carp population demographics in the lower Platte River basin and may provide useful information for the development of Silver Carp management plans in similar Great Plains streams.

Released November 17, 2025 07:48 EST

2025, Open-File Report 2025-1046

David L. George, Charles M. Cannon

This report describes a modeling investigation by the U.S. Geological Survey (USGS) of hazards in the Toutle and Cowlitz River valleys posed by hypothetical outburst floods from Spirit Lake, Washington. A massive debris avalanche resulting from the collapse of Mount St. Helens’ north flank during the May 18, 1980, eruption blocked Spirit Lake’s natural outlet into the North Fork Toutle River. Lacking a natural outlet, subsequent runoff in the Spirit Lake watershed contributed to a rising lake level, elevating the potential for debris-dam breaching or catastrophic failure. The influence of highly erodible bed sediment in the upper North Fork Toutle River on downstream flood and debris-flow dynamics and extent is assessed in this study. Simulations of clear-water (non-erosive) outburst floods were used as a baseline and compared to erosive flows that entrain large volumes of material and transition into debris flows along their flow path, revealing the influence of entrainment on hazard extent. Clear-water floods were modeled with the shallow water equations. Erosive flows were modeled with a two-phase granular fluid model that accommodates mobilization and incorporation of sediment from the bed into the overlying flow and resultant changes in flow rheology across a wide range of solid concentrations, from dilute suspensions to dense-granular debris flows. Entrainment of bed material was found to substantially increase the total flow volume (total volume of transported water and sediment is approximately 150 percent of the water volume for non-erosive flows). Erosive flows are shown to exhibit higher flow-front speeds and faster downstream arrival times than non-erosive flows, consistent with volume amplification effects near the actively mobilizing flow front. However, the larger total volume of transported material does not necessarily lead to an enhancement of total volume throughput (cumulative discharge) or inundation extent (total affected area) for all locations along the entire flow path; while entrainment leads to the displacement of a larger volume of material overall, much of this dislocated material (water and sediment) deposits upstream from the distal extent of the flows. These results are consistent with energetic considerations of initial potential energy and granular shear resistance.

Released November 15, 2025 09:01 EST

2025, Carbon Capture Science and Technology (17)

C. Ozgen Karacan

Residual oil zones (ROZs) can offer significant oil resources via enhanced oil recovery (EOR) as well as subsurface carbon dioxide (CO₂) retention during injection. If injected CO₂ is anthropogenic, the ROZs can offer a substantial geologic storage potential. The ROZs below the oil/water contact (OWC) of main pay zones (MPZ) in conventional reservoirs or brownfields, are more commonly developed for CO₂ injection and oil production and reported in the literature. However, CO₂-EOR in greenfield ROZs, reservoirs without a MPZ present, have rarely been developed for CO₂-EOR operation. The Tall Cotton Field of West Texas, Permian Basin, which started production in 2015 (Phase 1) and expanded in 2017 (Phase 2) from the San Andres Limestone, is one of the first examples of greenfield ROZs developed for EOR by injecting CO₂.

This paper analyses EOR and CO₂ retention performance of Tall Cotton Field using allocated injection and production data from inverted 5-spot well patterns of Phase-1 and -2 developments. Production and injection data allocated to each of the 28 identified patterns (nine 20-acre patterns for Phase-1, three 20-acre and sixteen 10-acre patterns for Phase-2) were analyzed for historical and forecasted oil recovery using ratio-trend decline analysis, and for CO₂ retention performance of the patterns. The allocated data were further used to calculate injected reservoir pore volume and void replacement ratios (VRR) for the analysis period. Quantitative results indicated that oil recovery factors of the 5-spot patterns varied between 4–10 %, and 5–30 % between the end of injection and the forecast periods, respectively. Storage of CO₂, on the other hand, increased to a mean value of ∼7130 MMscf per pattern in Phase-1 and to a mean storage of 3700 MMscf per pattern in Phase-2 until the end of injection, followed by a decline after the end of injection and into the forecast period. Resulting CO₂ utilization factors ∼6–50 Mscf/bbl were estimated at the end of injection. Overall, presented results suggested that developing greenfield ROZs for CO₂-EOR can be as promising as brownfield ROZs and mature MPZs for EOR and underground storage of injected CO₂. For Tall Cotton Field, results suggest that Phase-2 patterns generally outperformed Phase-1 for oil recovery factors, while Phase-1 performed better in CO₂ retention performance metrics. This is the first study in the literature that reports a detailed CO₂-EOR performance analysis of a greenfield ROZ in the Permian Basin, which can potentially allow for comparison with MPZs and brownfield ROZs.

Released November 14, 2025 14:55 EST

2025, Open-File Report 2025-1037

Mark D. Kozar, Samuel H. Austin

Seven potential sources of water, consisting of free-flowing discharge from abandoned coal mines at six locations and one abandoned flooded underground coal mine air shaft, were sampled for chemical analysis to assess the quality of the groundwater emanating from the seven mine sources. The six free-flowing mine discharge sources were also assessed for discharge by current-meter measurements on two separate occasions. The U.S. Geological Survey assessed these seven sources to provide information to the City of Gary, West Virginia (W. Va.), and the City of Gary’s consulting engineer with groundwater-quality and flow data to allow them to assess the seven sites as potential alternate sources of water for the City of Gary to augment its existing supply.

For the six sites where discharge could be measured, discharge ranged from a minimum of 0.082 cubic feet per second (ft³/s) to a maximum of 3.685 ft³/s. Of the six sites measured, only two, Harmon Branch at Thorpe, W. Va. (USGS site 372201081303501) and the abandoned public-supply water wells near Havaco, W. Va. (USGS site 372358081344601), had discharge in excess of 1.00 ft³/s. Discharge from the abandoned public supply wells was 3.685 ft³/s on September 20, 2023, and 2.888 ft³/s on October 16, 2023, and discharge from Harmon Branch at Thorpe, W. Va., was 1.049 ft³/s on September 22, 2023, and 1.038 ft³/s on October 17, 2023. Discharge in the abandoned underground mine air shaft (USGS site 372224081340901) could not be assessed, but the air shaft drains an abandoned mine that likely contains water stored in approximately 1.7 square miles (mi²) of abandoned underground coal mines in the Pocahontas No. 3 coal seam, and possibly an additional 0.9 mi² of leakage from the overlying Pocahontas No. 4 coal seam. Discharge for the six sites measured for the study was measured during a period between September 20 and October 18, 2023, and corresponded to the 12th to the 15th percentile of flow-duration statistics for the Tug Fork downstream of Elkhorn Creek at Welch, W. Va. streamgage (USGS site 03212750).

Water-quality data for the seven sites sampled overall were acceptable with respect to drinking water standards. Of the 203 constituents analyzed, only a few failed to meet applicable U.S. Environmental Protection Agency (EPA) drinking water standards. Iron exceeded the 300 micrograms per liter (μg/L) secondary maximum contaminant level (SMCL) at only 1 of the 7 sites (14.3 percent) sampled. Iron concentrations ranged from a minimum of less than (<) 5.00 μg/L to a maximum of 724 μg/L with a median concentration of 7.62 μg/L. Manganese exceeded the 50.0 μg/L SMCL at 2 of the 7 sites (28.6 percent) sampled. Manganese concentrations ranged from a minimum of 1.93 μg/L to a maximum of 271 μg/L with a median concentration of 4.03 μg/L. No sites sampled exceeded the arsenic maximum contaminant level (MCL) of 10 μg/L. Arsenic concentrations ranged from a minimum of <0.100 μg/L to a maximum of 2.35 μg/L with a median arsenic concentration of 0.200 μg/L. None of the seven sites sampled for selenium for this study exceeded the EPA MCL of 50.0 μg/L. Selenium concentrations ranged from a minimum of <0.050 μg/L to a maximum of 5.26 μg/L with a median concentration of 3.21 μg/L.

All seven sites were sampled for volatile organic compounds (VOCs), semivolatile organic compounds (SVOCs), and polychlorinated biphenyls (PCBs), but most had concentrations below the detection limit. Of the 10 PCB compounds analyzed for the seven sites sampled, none contained detectable concentrations of PCBs or Aroclor compounds. Of the 44 SVOCs analyzed at each of the seven sites sampled, only 1 SVOC, acenaphthene, was detected, at a concentration of 0.02 μg/L. Of the 96 VOCs analyzed, from each of the seven sites sampled, only two were found at detectable concentrations. Trichloromethane was detected only at 1 of the 7 (14.3 percent) sites sampled at a concentration of 0.027 μg/L, and benzene was detected at the same site and 3 additional sites (4 of the 7 sites or 57.1 percent of the sites sampled) at concentrations of 0.028, 0.029, 0.021, and 0.035 μg/L, but none exceeded the EPA MCL for benzene of 5.00 μg/L.

Total coliform bacteria are ubiquitous in the environment, and their presence only suggests the potential for contamination by near-surface processes. Escherichia coli (E. coli) bacteria are derived from either human or animal fecal material and can be an indicator of potential contamination by pathogenic bacteria or viruses. Total coliform bacteria were detected at all 7 sites sampled at concentrations ranging from 17.5 to greater than (>) 2,420 most probable number per 100 mL (MPN/100 mL) of sample, with a median total coliform concentration of 1,553 MPN/100 mL. Escherichia coli bacteria were detected at 4 of the 7 sites sampled at concentrations ranging from <1 to 11.9 MPN/100 mL, with a median E. coli concentration of 5.1 MPN/100 mL.

Released November 14, 2025 10:10 EST

2025, Arctic Science (11) 1-3

Anthony M. Pagano, Stephen N. Atkinson, Louise C. Archer

The spring is a critical period when polar bears (Ursus maritimus Phipps, 1774) are thought to have peak access to seals and acquire the majority of their annual energy requirements during a period of hyperphagia. Pagano et al. (Pagano A.M., Atkinson S.N., and Archer L.C. 2025. Arctic Science.11:1-14. doi:10.1139/as-2024-0051) examined the intra-seasonal changes in body mass of 31 polar bears on the spring sea ice and found polar bears exhibited a feast or famine lifestyle. A lack of a relationship between changes in body mass and recapture date suggested that many bears had not entered their primary period of hyperphagia. Pilfold extended our discussion to conclude that our data show polar bear hyperphagia begins after the period of ringed seal (Pusa hispida Schreber, 1775) whelping, and discusses this in relation to previous work on the timing of polar bear seal kills. Here, we reassess whether our data provide information on the timing of polar bear hyperphagia. We find no relationships in our data to conclude when polar bear hyperphagia begins. Instead, our data highlight the conflicting pressures individuals face between the spring breeding season, when time spent foraging is often reduced to engage in mating behavior, and the spring hyperphagia period, when the bulk of annual energy requirements are met.

Released November 14, 2025 09:52 EST

2025, Fact Sheet 2025-3045

James S. Webber, Kaylyn S. Gootman, Kenneth Hyer, Peter J. Tango, Douglas L. Moyer

The lands and waters of the Chesapeake Bay watershed provide more than $100 billion in economic benefits- an amount that is expected to increase by achieving the region’s clean-water goals. Achieving those goals requires accurate and timely information about the health of the watershed’s rivers and streams. The Chesapeake Bay nontidal monitoring network (NTN), a partnership of local, state, and federal agencies, as well as other partners, was established in 2004 to provide this information. The U.S. Geological Survey analyzes data collected from NTN stations to provide monitoring-based information about the amount of nitrogen, phosphorus, and sediment entering the Chesapeake Bay through its nontidal rivers. Thus, data collected from the NTN inform watershed management by providing decision makers with information on which to base their restoration and conservation actions.

Released November 14, 2025 09:00 EST

2025, Fact Sheet 2025-3044

Dan Walters

Introduction 

Massachusetts extends from the mountains of the Appalachian system in the west of the State to the sandy beaches and rocky shorelines of the Atlantic coast in the east. Inland topographic data support a wide range of important activities, including geologic mapping, transportation planning, forest and wildlife management, quantifying ecological services, water supply protection, commonwealth-wide infrastructure planning, local site planning, and flood-plain management. Nearshore bathymetry can be used to support coastal portions of the Commonwealth by addressing the combined threats of ocean warming, strong storm surge, and rising sea levels. The maintenance and (or) expansion of Massachusetts ports (for instance, Boston, New Bedford) and Cape Cod sediment management depends upon the accurate mapping of bathymetry and the frequent influx of sediment and redeposition. Critical applications that address the broad range of requirements 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 Massachusetts. The status of available and in-progress 3DEP baseline lidar data in Massachusetts 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 $1.23 million in new benefits annually to Massachusetts. The top 10 Massachusetts business uses for 3D elevation data, which are based on the estimated annual conservative benefits of 3DEP, are shown in table 2.

Released November 13, 2025 09:07 EST

2025, Canadian Journal of Microbiology

Jeffery D. Sullivan, Rebecca L. Poulson, Glenn H. Olsen, Alicia Berlin, Zijing Cao, Deborah Carter, Josh Homyack, Jennifer Kilburn, Scott R. McWilliams, Joshua Osborn, Tori Mezebish Quinn, Hannah Schley, Matthew M. Weegman, Christopher A. Williams, David E. Stallknecht, Diann Prosser

Clade 2.3.4.4b Eurasian-origin H5N1 entered North America in late 2021 and spread across the continent. While studies have characterized the antibody response mounted by dabbling ducks following exposure, little data are available for diving ducks. This study sought to identify influenza A virus (IAV) infection and antibodies in Lesser and Greater Scaup captured in Maryland, Illinois, and Rhode Island. In Maryland, IAV seroprevalence increased from the 2021/2022 to 2022/2023 sampling season, with IAV antibody prevalence increasing for juvenile (38% to 80%) and adult (82% to 90%) Lesser Scaup. While adult Lesser Scaup sampled in Illinois in 2021/2022 had IAV antibody prevalence comparable to those sampled in Maryland (76% and 82%, respectively), they had higher antibody prevalence to both H5 (48% and 18%) and N1 (68% and 35%), potentially due to being sampled in March versus December and January. Our data suggest that Lesser Scaup had limited antibodies to highly pathogenic H5 IAV prior to the introduction of clade 2.3.4.4b H5N1 to North America, but relevant antibodies were widely observed in the months and year following. Our more limited data suggest similar trends may have occurred in Greater Scaup as well.

Released November 13, 2025 07:54 EST

2025, Molecular Ecology Resources (26)

Rodney T. Richardson, Grace Avalos, Cameron J. Garland, Regina Trott, Olivia Hager, Mark J. Hepner, Clayton D. Raines, Karen Goodell

Terrestrial environmental DNA (eDNA) techniques have been proposed as a means of sensitive, non-lethal pollinator monitoring. To date, however, no studies have provided evidence that eDNA methods can achieve detection sensitivity on par with traditional pollinator surveys. Using a large-scale dataset of eDNA and corresponding net surveys, we show that eDNA methods enable sensitive, species-level characterisation of whole bumble bee communities, including rare and critically endangered species such as the rusty patched bumble bee (RPBB; Bombus affinis). All species present in netting surveys were detected within eDNA surveys, apart from two rare species in the socially parasitic subgenus Psithyrus (cuckoo bumble bees). Further, for rare non-parasitic species, eDNA methods exhibited similar sensitivity relative to traditional netting. Compared with flower eDNA samples, sequenced leaf surface eDNA samples resulted in significantly lower rates of Bombus detection, and these detections were likely attributable to high rates of background eDNA on environmental surfaces, perhaps due to airborne eDNA or eDNA movement during rainfall events. Lastly, we found that eDNA-based frequency of detection across replicate surveys was strongly associated with net-based measures of abundance across site visits. We conclude that the COI-based metabarcoding method we present is cost-effective and highly scalable for quantitative characterisation of at-risk bumble bee communities, providing a new approach for improving our understanding of species habitat associations.

Released November 12, 2025 09:17 EST

2025, Science Advances (11)

Jinxun Liu, Benjamin M. Sleeter, Zhiliang Zhu, Mark A. Cochrane, Qiang Zhou, Bin Wang, Grant Domke, Paul Selmants, Lisamarie Windham-Myers, Qiuan Zhu, Tamara Wilson, Kristin Byrd, Eric Ward, Terry Sohl, Todd Hawbaker, Zhen Zhang, Christopher Soulard, Kimberly Wickland, Robert G. Striegl

The land carbon sink of the conterminous United States was evaluated using a bottom-up modeling framework and 30-meter land change data from 1985 to 2020. This cross-scale, cross-landscape, and cross-system approach tracked fractional land cover changes and applied regional model calibration. Results show average terrestrial and aquatic carbon sinks of +110 ± 37 and +19 ± 0.5 teragrams of carbon per year, respectively. The terrestrial carbon sink, showing no clear trend, peaked in the 1990s, with more years as a carbon source since 2000, contradicting recent national and global studies. Land change had the largest impact (−70 ± 5.5 teragrams of carbon per year), exceeding impacts of climate (−33 ± 48 teragrams of carbon per year), wildfire (−7.7 ± 2.4 teragrams of carbon per year), and erosion transport (−1.9 ± 0.13 teragrams of carbon per year). The positive CO2 fertilization effect (+69 ± 12 teragrams of carbon per year) was insufficient to maintain the carbon sink strength. Our framework reveals key paths of carbon loss, with implications for carbon budget and energy policies in the United States and beyond.

Released November 12, 2025 08:45 EST

2025, Geochemistry: Exploration, Environment, Analysis (25)

K. E. Miltenberger, Zachary Shephard, Rachel Lynn Mixon, Johanna Blake, Shaleene Chavarria, Douglas Yager

Basin morphometry, climate and geology control how a hydrological network evolves over time, controlling the efficiency of weathering of elements from geological materials, and ultimately the input of sediment and dissolved constituents to river systems. Exceedances to the Navajo Nation surface water quality standards for trace metals have been reported in the San Juan River watershed. Because metals are transported adsorbed to fine-grain sediment, the identification of areas with elevated sources of trace metals and/or areas with increased erosion and sediment transport potential is an important first step in protecting water quality. Physical factors such as elevation, slope, relief and stream order were used to quantify morphometric parameters that effect the contribution of trace metals into the stream network. By correlating these parameters with water quality data that were collected from tributaries along the San Juan River, we identified statistically significant regressions between morphometric parameters and total Al, Pb, U, Fe and Mn in surface water. Positive correlations with trace metals include tributary drainage basin perimeter, pour point elevation and total number of streams, while negative correlations include stream length ratio, ruggedness number and longest basin axis. Stream reach measurements within geological units that contain known trace metal constituents reveal that Gallegos Canyon and Desert Creek are the most susceptible to sediment mobilization and transport, while other tributary drainage basins, such as Desert, Recapture and Salt creeks, are associated with naturally elevated concentrations of Al, As, Pb and U.

Released November 11, 2025 08:32 EST

2025, Journal of Caribbean Ornithology (38) 59-66

Yvan G. Satgé, J. Brian Patteson, Bradford S. Keitt, Chris P. Gaskin, Patrick G.R. Jodice

Pterodroma hasitata, the Black-capped Petrel (locally known as Diablotin), is the only extant Pterodroma petrel nesting in the Caribbean. The species is listed as globally Endangered by the IUCN and was recently listed as endangered under the U.S. Endangered Species Act. Pterodroma hasitata show a phenotypic gradient, ranging from a darker, smaller form to a paler, heavier form, that is reflected in a strong genetic structure. This phylogenetic divergence suggests the existence of at least two distinct breeding populations. We report on pre-breeding movements of two male Pterodroma hasitata, one of each form, tracked by satellite from non-breeding areas in Gulf Stream waters of the western North Atlantic Ocean to breeding locations in Hispaniola in late 2019. Based on a combination of tracking locations, location error classes, battery voltage, and satellite communication schedules, we infer that the light-form petrel visited a nest in central Dominican Republic during 2 to 8 October and 9 to 15 October, and the dark form visited a nest in southeastern Haiti during 9 to 22 November and 29 November to 3 December. This information supports earlier suggestions that Pterodroma hasitata forms breed in allochrony and in allopatry, both of which may be a driver of speciation.

Released November 10, 2025 10:52 EST

2026, Journal of Hydrology (664)

Nidia Tobon-Velazquez, Gerd Masselink, T.J. O’Hare, Robert Bates, Ferdinand Oberle, Curt D. Storlazzi, D. C. Conley

This study examines the spatial variability of shallow groundwater on Dhigelaabadhoo Island using electromagnetic induction surveys, groundwater monitoring, and sediment analyses. The research reveals how variations in island morphology—such as differences in elevation, reef flat width, and sediment composition—affect the spatial distribution of groundwater lenses and the overall aquifer dynamics. Saltwater intrusion is especially pronounced in low elevated areas, with narrow reef flat plate and areas where higher hydraulic conductivity—driven by the presence of coarser sediments—is observed, whereas regions characterized by finer sediments, higher elevation, and wider reef flat plates tend to support more symmetric and less saline groundwater lenses. The geophysical investigations reveal that tidal oscillations alter groundwater movement by markedly changing water levels and conductivity, thereby underscoring the critical need to account for temporal dynamics in atoll coastal aquifer systems and the importance of integrating tidal dynamics into the aquifer zone. The findings highlight the significant role of intrinsic morphological and external hydrodynamic factors in shaping groundwater distribution on atoll islands, offering critical insights for sustainable freshwater resource management.

Released November 07, 2025 09:03 EST

2025, Seismological Research Letters

Jeffrey J. McGuire, Andrew J. Barbour, Zack J. Spica, Verónica Rodríguez Tribaldos, Zhongwen Zhan, Bradley P. Lipovsky, Robert J. Mellors, Ettore Biondi, Clara Yoon, Martin Karrenbach, Adam T. Ringler, James William Atterholt, Avinash Nayak, Theresa Marie Sawi, Loic Viens, Eileen R. Martin, Allen L. Husker, Paul Bodin, Morgan P. Moschetti, Qibin Shi, Nathaniel C. Miller, Prisha Puri

The use of fiber‐optic sensing systems in seismology has exploded in the past decade. Despite an ever‐growing library of ground‐breaking studies, questions remain about the potential of fiber‐optic sensing technologies as tools for advancing if not revolutionizing earthquake‐hazards‐related research, monitoring, and early warning systems. A working group convened to explore these topics; we comprehensively examined the application of fiber optics in various aspects of earthquake hazards, encompassing earthquake source processes, crustal imaging, data archiving, and technological challenges. There is great potential for fiber‐optic systems to advance earthquake monitoring and understanding, but to fully unlock their capabilities requires continued progress in key areas of research and development, including instrument testing and validation, increased dynamic range for applications focused on larger earthquakes, and continued improvement in subsurface and source imaging methods. A key current stumbling block results from the lack of clear data archiving requirements, and we propose an initial strategy that balances data volume requirements with preserving key data for a broad range of future studies. In addition, we demonstrate the potential for fiber‐optic sensing to impact monitoring efforts by documenting the data completeness in a number of long‐term experiments. Finally, we outline the features of a instrument testing facility that would enable progress toward reliable and standardized distributed acoustic sensing data. Overcoming these current obstacles would facilitate progress in fiber‐optic sensing and unlock its potential application to a broad range of earthquake hazard problems.

Released November 07, 2025 08:58 EST

2025, Nature Geoscience (18) 1077-1080

Lara Mani, Jenni Barclay, Carina Fearnley, Richard E.A. Robertson, Blaise Mafuko Nyandwi, Sara Barsotti, Amy Donovan, Wendy K. Stovall

Since the 1985 Nevado del Ruiz eruption that killed over 23,000 people in Armero, Colombia, risk communication has become central to volcanic crisis management. Despite the development of effective tools and programmes for volcanic risk communication, considerable challenges remain.

Released November 06, 2025 10:00 EST

2025, Science (390) 628-632

Frank J. Pavia, Jesse Farmer, Laura Gemery, Thomas M. Cronin, Jonathan Treffkorn, Kenneth A. Farley

Arctic sea-ice loss affects biological productivity, sustenance in coastal communities, and geopolitics. Forecasting these impacts requires mechanistic understanding of how Arctic sea ice responds to climate change, but this is limited by scarce long-term records. We present continuous 30,000-year reconstructions of sea-ice coverage from the Arctic Ocean based on measurements of two isotopes, thorium-230 and extraterrestrial helium-3, whose burial ratio changes with sea-ice coverage. We found that the central Arctic was perennially covered by sea ice during the last glaciation. Sea-ice cover retreated during the deglaciation approximately 15,000 years ago, culminating in seasonal sea-ice coverage in the warm early Holocene, before ice coverage increased into the late Holocene. Sea-ice changes closely correlate with biological nutrient consumption, supporting projections of a nutrient-starved central Arctic Ocean with continued sea-ice loss.

Released November 06, 2025 09:41 EST

2025, Herpetologica (81) 336-345

Brian Halstead, Daniel Antonio Macias, Casey D. Moss, Patrick M. Kleeman, Jonathan P. Rose

Estimating the distributions of cryptic species is essential for conservation, yet our understanding is hampered by animal behavior and imperfect detection. We developed and implemented a multiscale occupancy survey protocol to estimate the probability of occurrence, probability of being active on the surface, and detection probability of two range-restricted terrestrial salamanders, Scott Bar Salamanders (Plethodon asupak) and Siskiyou Mountains Salamanders (P. stormi), in interior northern California, USA. We established survey sites near locations of historical occurrence of these salamanders and surveyed each site on one to six visits in late fall 2023 and spring 2024. We compared models with different environmental variables for predicting salamander occurrence and surface activity. Much model selection uncertainty in the effects of covariates existed, but the model with most support indicated that Plethodon salamanders in interior northern California were more likely to occur at sites near recent historical occurrences that had longer growing seasons after controlling for elevation and aspect. Plethodon stormi or P. asupak surface activity was higher at night than during the day and was highest when substrates were cool and moist. Our survey protocol was successful for quantifying the effects of site and visit characteristics on P. stormi or P. asupak occurrence and availability; minor modifications will likely improve its utility for providing unbiased inference about salamander occurrence and surface activity.

Released November 06, 2025 09:37 EST

2025, PLoS ONE (20)

Nora Boross, Ardo Laszlo, Duane C. Chapman, Gergely Boros, Zoltán Vitál, Viktor Tóth, Nathan Thompson, Katy E. Klymus, Catherine A. Richter

Silver carp (Hypophthalmichthys nobilis), bighead carp (H. molitrix) and their hybrids, collectively known as bigheaded carps, have been introduced to Lake Balaton, Hungary. The current stock sizes are difficult to assess. We investigated environmental DNA (eDNA) techniques targeted for bigheaded carps, assessed the spatial distribution of eDNA in Lake Balaton, compared eDNA concentrations to environmental variables to assess potential habitat selection based on those variables, and provided an estimate of biomass of bigheaded carps relative to eDNA shedding rates per unit biomass observed in controlled experiments. Water samples were collected from 70 sites in an array across the lake. Biomass estimation was calculated using mean eDNA concentration obtained by quantitative PCR of the samples and previously determined eDNA shedding rates of bigheaded carps under controlled conditions in a laboratory. Concentration of eDNA was highly variable between sites, resulting in wide confidence intervals. Basins did not significantly differ in eDNA concentration, and there were no strong relationships between environmental variables and eDNA concentration, indications that bigheaded carps use the entire lake. The model provided an estimate of 4,830 metric tonnes (2,750–8,030 tonnes) of bigheaded carps in Lake Balaton, or 81.0 kg/ha. The eDNA method produced a value close to previous estimates by traditional means of total biomass of bigheaded carps in the lake, and like traditional methods, there was a broad confidence interval on the estimate of the mean. The results of the present study support the utility of aquatic eDNA analysis, and the need for further comparisons with fisheries methods and supporting data from laboratory studies.

Released November 06, 2025 07:48 EST

2025, Science of the Total Environment (1005)

Marco Milardi, Louisa E. Wood, Elizabeth A. Nyboer, Holly Susan Embke, Sui C. Phang, Abigail J. Lynch

Inland recreational fisheries are globally significant leisure pursuits, with well-documented benefits to human health and well-being, but also one of the principal drivers of non-native fish introductions to enhance fishing opportunities, whether for sport or sustenance. In this study, we assess the relative reliance of global inland recreational fisheries on non-native versus native species for harvest. We further examine how this reliance varies by economic and nutritional value as well as the climate vulnerability of the species involved. We demonstrate that, of the 1,325,851 t of inland recreational fishes recreationally harvested for consumption worldwide in 2021, non-native fish were a small proportion (4 %; 53,651 t). On a global scale, non-native fish contributed a net positive 38.2 % economic value to inland recreational harvest. However, they also contributed a net negative −21.9 % nutritional value to inland recreational harvest. Non-native fishes were also more climate vulnerable (i.e., higher average climate vulnerability index values) and thus proportionally increased overall estimates of climate vulnerability with a net positive of 70.9 %. Our results quantitatively demonstrate that non-native species play a more important role in inland consumptive recreational fisheries than their mere harvest volume would suggest. However, many nuances were seen on the continent and country scale, which reflect the complexity of fisher behavior, fish distribution and socio-economic factors. Our findings help unravel the complex effects of non-native species on human activities and underscore the need to evaluate their global impacts holistically.

Released November 05, 2025 10:43 EST

2025, Seismological Research Letters

Adam T. Ringler, Andrew Holcomb, E. Joshua Rigler, Spencer Wilbur, C. Balch, Corey Beutel, Brendan Ryan Geels, J. Guerra, A. Horton, Edward Kromer, Kristen A. Lewis, Jeffrey J. Love, Yolando Root, Claudia Kristina Rossavik, N. Shavers, John Spritzer, Tyler Storm, Alexandra Nicole Wernle, David C. Wilson

The U.S. Geological Survey’s Geomagnetism Program is collaborating with the Earthquake Hazards Program and Global Seismographic Network Program to densify magnetic field observations. This collaboration focuses on the installation of magnetometers, or magnetic variometers, at existing seismic stations. Along with improving the density of space weather observations for hazard monitoring, these data can be used to correct colocated magnetic field induced noise in seismic data. Such corrections are especially useful during time periods of large magnetic storms where the magnetic field‐induced instrument noise can be of similar amplitude to earthquake ground‐motion records.

Released November 05, 2025 10:14 EST

2025, Ecological Monographs (95)

Michael Dumelle, Rob Trangucci, Amanda M. Nahlik, Anthony R Olsen, Kathryn Irvine, Karen A. Blocksom, Jay Ver Hoef, Claudio Fuentes

In ecology and related sciences, missing data are common and occur in a variety of different contexts. When missing data are not handled properly, subsequent statistical estimates tend to be biased, inefficient, and lack proper confidence interval coverage. Missing data are often grouped into three categories: missing completely at random (MCAR), missing at random (MAR), and missing not at random (MNAR). We review each category and compare their benefits and drawbacks. We review several approaches to handling missing data including complete case analysis, imputation, inverse probability weighting, and data augmentation. We clarify what types of variables should accompany imputation methods and how those variables are influenced by the analysis methods. Additionally, we discuss missing data that lack a formal basis for measurement and hence are fundamentally different from MCAR, MAR, and MNAR missing data. Throughout, we introduce concepts and numeric examples using both simulated data and data from the United States Environmental Protection Agency's 2016 National Wetland Condition Assessment. We conclude by providing five considerations for ecologists and other scientists handling missing data.

Released November 05, 2025 09:10 EST

2025, Journal of Fish Biology

Ella K. Humphrey, Jonathan J. Spurgeon, Lizabeth Bowen, Robert E. Wilson, Shannon C. Waters-Dynes, Braxton M. Newkirk, Sarah A. Sonsthagen

Environmental and associated ecosystem change may affect the persistence of fish species based on their ability to adapt to changing conditions, including decreasing flows and rising water temperatures. Exceeding the thermal tolerances of stream fish will likely result in a loss of ability to maintain metabolic processes. We evaluated the critical thermal maximum (CTmax) of bigmouth shiner (Ericymba dorsalis) and analysed the expression of heat shock protein 70 messenger RNA (mRNA) (HSP70) to quantify a thermal stress response over a gradient of temperatures (25°C–31°C). E. dorsalis HSP70 mRNA expression was upregulated in response to temperatures >25°C, indicating a stress response. This study supports the existence of a thermal stress threshold for E. dorsalis. The frequency at which this threshold is exceeded may increase under forecasted future climate scenarios for Nebraska.

Released November 05, 2025 08:00 EST

2025, GSA Bulletin

Abdullah Aohali, Shanaka L. de Silva, Alejandro Cisneros de Leon, Charles Lewis, Axel K. Schmitt, Martin Danišík, Mark E. Stelten, Sujoy Mukhopadhyay, Robert Duncan, Frank C. Ramos

Determining accurate and precise ages for Quaternary volcanic centers is essential for reconstructing volcanic field histories, understanding magmatic processes, and assessing potential hazards or risk. Harrat Khaybar, western Saudi Arabia, is one of the youngest and potentially most active volcanic fields on the Arabian plate, has been active since ca. 1.7 Ma, and is characterized by a spectrum of rock compositions ranging from predominantly alkalic basalt to trachyte and comendite. Previous work in Harrat Khaybar utilizing ⁴⁰Ar/³⁹Ar incremental heating geochronology to constrain morphological preservation and superpositional relationships bracketed the volcanic activity into broad age groups in intervals of ∼150 k.y., and the youngest and most compositionally evolved volcanoes, including Jabal Abyad, Jabal Bayda, and Jabal Qidr, were assigned to the age groups between ca. 300 ka and present. Herein, we establish a detailed chronology of prehistoric silicic and historical basaltic eruptions at central Harrat Khaybar using four independent eruption age determination techniques: zircon double-dating (ZDD), which combines ²³⁸U-²³⁰Th disequilibrium or U-Pb with (U-Th)/He dating; zircon U-Pb dating; cosmogenic ³He dating; and cosmogenic ³⁶Cl geochronology. These were employed to accurately date six volcanic centers, including the comenditic Jabal Abyad, Jabal Bayda, Jabal Ibayl, and Jabal Alhayyirah, the trachytic Jabal Aluthmor, and the basaltic Jabal Qidr. Additionally, our previously published ⁴⁰Ar/³⁹Ar ages have been recalculated using isochron intercept (nonatmospheric) ⁴⁰Ar/³⁶Ar for the trapped Ar component. Our new results reveal that zircon rims from Jabal Abyad and Jabal Bayda define isochron ²³⁸U-²³⁰Th crystallization ages of 125 ± 4 ka and 144 ± 6 ka, concordant with ZDD eruption ages of 132 ± 4 ka and 149 ± 5 ka, respectively. Zircon U-Pb crystallization and (U-Th)/He eruption ages from Jabal Alhayyirah are concordant at 471 ± 14 ka and 458 ± 18 ka, respectively. Finally, zircons from the nearby Jabal Ibayl yield a U-Pb weighted mean crystallization age of 566 ± 16 ka concordant with the corresponding (U-Th)/He eruption age of 554 ± 12 ka, both of which are notably older than the previously proposed eruption ages of 300−150 ka. Recalculations of published ⁴⁰Ar/³⁹Ar ages for the youngest volcanoes at central Harrat Khaybar are now in excellent agreement with new geochronological data.

Our new age data reveal several new insights into the development of Harrat Khaybar. It is now clear that the comenditic eruptions do not belong to the same eruptive phase and indicate an extended history during which comendites have episodically punctuated the basaltic volcanism since at least 600 ka. The data indicate that Jabal Ibayl and Jabal Alhayyirah represent separate older volcanic events, whereas the younger Jabal Abyad and Jabal Bayda volcanoes appear to be coeval, and their spatial proximity implies that they share a magmatic lineage and maybe a common plumbing system. Zircon age spectra of the comendites reveal obvious xenocrysts and antecrysts indicating assimilation of basement and plutonic progenitors, but otherwise they define broad unimodal populations of crystallization ages that overlap within error with the respective eruption ages. We interpret this to indicate that zircon crystallization continued up to the time of eruption. The cosmogenic ages of Jabal Aluthmor and Jabal Qidr reveal that these centers erupted as recently as ca. 2000 and 760 years ago, respectively. Broadly coeval young silicic and basaltic eruptions at northern Harrat Rahat, the harrat immediately south of Harrat Khaybar, may imply a shared geodynamic forcing between the two adjacent volcanic fields.