Skip to main content
U.S. flag

An official website of the United States government

Dot gov

Official websites use .gov
A .gov website belongs to an official government organization in the United States.

Https

Secure .gov websites use HTTPS
A lock ( ) or https:// means you’ve safely connected to the .gov website. Share sensitive information only on official, secure websites.

Publications recently added to the Pubs Warehouse

(500 records max)
High-resolution magnetic survey using an unoccupied aerial vehicle to constrain buried lava flow geometry, volume, and eruptive history of Little Cones, Crater Flat, Nevada

Released June 05, 2026 14:42 EST

2026, Professional Paper 1890-R

Robert Van Alphen, Mel Rodgers, Rocco Malservisi, Charles B. Connor, Rachel Bakowski, Troy Berkey

Magnetic surveys are an important tool used to augment geologic mapping in distributed volcanic fields. Using magnetic anomalies, it is possible to model the geometry of shallowly buried volcanic features, such as conduits, sills, and lava flows. This subsurface mapping is important for understanding eruption dynamics and emplacement of lava flows, and it sometimes reveals buried volcanoes no longer visible at the surface. These data are critical to better interpret the numbers, styles, and magnitudes of eruptions in distributed volcanic fields and their associated volcanic hazards. New advances in unoccupied aerial vehicles (UAVs) offer an attractive middle range of resolution and aerial coverage between ground-based magnetic surveys and aeromagnetic surveys.

Here, we present the results of a UAV fluxgate magnetic survey of the Little Cones, Nevada, scoria cones, which have been the target of previous ground and aeromagnetic surveys. The magnetic anomalies at Little Cones are of interest because the surrounding alluvium conceals lava flows that erupted from Little Cones, making it very difficult to understand the volume and morphology of lava flows from geologic mapping alone. Nonlinear inversion of UAV-collected magnetic data were used to model the thickness and morphology of buried Little Cones’ lava flows with higher precision than achieved previously. The sequence of events and calculated flow characteristics are then interpreted. The total volume of Little Cones, including concealed lava flows, is approximately 0.016 cubic kilometer, and the initial sheet flow erupted in less than 24 hours. The findings presented herein demonstrate that UAV-based magnetic surveys are a reliable method of data collection and an efficient alternative to other survey methods, facilitating development of a three-dimensional perspective of distributed volcanic fields.

Global pegmatite-hosted lithium, cesium, and rubidium resources: A dataset for grade and tonnage modeling

Released June 04, 2026 08:57 EST

2026, Ore Geology Reviews (194)

Joshua Mark Rosera, Dalton M. McCaffrey, Niki E. Wintzer

Quantitative mineral resource assessments of potential undiscovered deposits can inform future mineral supply scenarios, but their accuracy is conditional on building robust grade and tonnage models of known deposits. This study presents an up-to-date global compilation and analysis of recently discovered and original, in-situ pegmatite-hosted Li, Cs, and Rb resources prior to historic production. Our analysis yields a median tonnage of 21.2 million tons (Mt) and grade of 1.12% Li2O, respectively, for global Li pegmatite deposits (n = 73). The grades and tonnages of Li pegmatite resources vary depending on the age of the bedrock host domain, pegmatite crystallization age, and primary ore mineralogy. Lithium pegmatite resources hosted in Archean to transitional Archean-Paleoproterozoic domains have the largest median tonnage (29.8 Mt; n = 38), and those hosted in Paleoproterozoic to Mesoproterozoic domains have smaller median tonnages (6.5 Mt; n = 16). Cesium deposits where pollucite is the primary ore mineral have a bimodal grade distribution, with modes of 2.40 and 0.035 wt% Cs2O for high- and low-grade deposits, respectively, while Rb deposits are more unimodal with a median grade of 0.247 wt% Rb2O. Pegmatite-hosted Cs and Rb resources have median tonnages of 7.6 and 6.3 Mt, respectively. Covariation between ore mineralogy and the degree of crustal enrichment in pegmatite-hosted deposits is diagnostic of petrogenetic differences, including melt source characteristics, magma evolution, or variable degrees of volatile solubility. The Li pegmatite compilation is suitable for fitting robust numerical models to support quantitative assessments. More well-defined Rb and Cs pegmatite resources are required for quantitative assessments, but these data provide useful information about original in-place resources for framing supply discussions.

Arizona Water Science Center activities at Lees Ferry, Arizona

Released June 03, 2026 14:15 EST

2026, Fact Sheet 2026-3002

Kathryn Anne Cooney

Introduction 

In 1921, the U.S. Geological Survey (USGS) established a streamgage on the Colorado River at Lees Ferry, Arizona, to monitor the river’s flow and level as it enters Grand Canyon. The following year, the seven States encompassing the Colorado River Basin (Arizona, California, Colorado, Nevada, New Mexico, Utah, and Wyoming) negotiated the 1922 Colorado River Compact to regulate distribution of the river’s waters between them. The compact divided the basin into two regions—the Upper Basin and the Lower Basin—and established the dividing point between them about one mile downstream from Lees Ferry, just below the confluence of the Colorado and Paria Rivers.

The Colorado River at Lees Ferry streamgage (USGS station 09380000) is one of the most important streamgages in the United States because it is used to measure how much water passes from the Upper Basin to the Lower Basin through Glen Canyon Dam. The dam, constructed between 1956 and 1966, generates hydropower and stores water in Lake Powell reservoir, which is used to provide Upper and Lower Basin states with the water allotted to them by the compact. Lower Basin states depend on releases from the dam to receive their allotments. The Lees Ferry streamgage, located less than 16 miles downstream from Glen Canyon Dam, produces publicly available, real-time water data that allows the Colorado River’s streamflow below the dam to be monitored.

Most years, the Colorado River runs dry before reaching its historical terminus at the Gulf of California in Mexico, so measuring and monitoring the river at Lees Ferry is critical for the Lower Basin ecosystems, agricultural resources, and municipal industries that rely on the river’s every drop. Additionally, Grand Canyon river guides and recreationalists depend on water level data from the Lees Ferry streamgage to determine when to run rapids and camp on sandbars. Streamflow and water-quality data collected at Lees Ferry are also important for monitoring the health of the Colorado River’s aquatic life because some species, including fish and macroinvertebrates, require certain water conditions to survive, reproduce, and spawn.

The Arizona Water Science Center is responsible for maintaining and collecting water data from the Lees Ferry streamgage. The Arizona Water Science Center is a branch of the USGS dedicated to providing high quality, impartial water data to resource managers and the public for their use in understanding and managing critical water resources in Arizona and the Southwest.

Black abalone surveys at Naval Base Ventura County, San Nicolas Island, California—2023 annual report

Released June 03, 2026 08:16 EST

2026, Open-File Report 2026-1015

Michael C. Kenner, Julie L. Yee

The U.S. Geological Survey monitors a suite of intertidal black abalone (Haliotis cracherodii) sites at San Nicolas Island, California, in cooperation with the U.S. Navy, which owns the island. The nine rocky intertidal sites were established in 1980 to study the potential effect of translocated sea otters on the intertidal black abalone population at the island. The sites were monitored from 1981 to 1997, typically annually or biennially. Monitoring resumed in 2001 and has been completed annually thereafter. Since 2018, the work has been carried out by the U.S. Geological Survey Western Ecological Research Center. The study sites became particularly important, from a management perspective, after a virulent disease decimated black abalone populations throughout southern California beginning in the mid-1980s. The disease, withering syndrome (Candidatus Xenohaliotis californiensis), was first observed on San Nicolas Island in 1992 and over the next few years, withering syndrome reduced the black abalone population on San Nicolas Island by more than 99 percent. In 2009, the black abalone subsequently was listed as endangered under the Endangered Species Act.

The subject of this report is the 2023 survey of the sites and the status of the measured population in comparison to long-term patterns (based on data collected since 1981) at San Nicolas Island. Between the years 2000 and 2023, the total monitored black abalone population at the island has grown from roughly 200 to more than 2,500 abalone following disease-related decline. Since it was first consistently measured in 2005, the average distance between adjacent black abalone has decreased substantially from approximately 50 centimeters to less than 15 centimeters, indicating that abalone are sufficiently close together at several of the sites to reproduce successfully. The total abalone count in 2023 was 2,570, which was 19.2 percent higher than in 2022 and the highest count since 1993. All nine sites had higher counts in 2023 than in the previous year. Over 25 percent of the black abalone counted in 2023 were classified as recruits, defined as having a shell length of 3 centimeters or less.

Summaries of goals, actions, and information needs by management entity

Released June 02, 2026 16:30 EST

2026, Scientific Investigations Report 2025-5018-B

Christine D. Miller Hesed, Heather M. Yocum, editor(s)

The grasslands in the North Central region are managed by a diverse group of Federal, State, and Tribal agencies; nongovernmental organizations; partnerships; and private landowners. This chapter highlights these various grassland management entities, provides background information on their mission and organizational structure, and describes some of their key grassland management activities, including the way in which each entity engages private landowners in grassland management. Each section also describes emerging challenges and opportunities and high-level information needs. The review and synthesis of grassland management-related documents identified specific information needs, which are listed in an appendix to provide additional detail for anyone looking to collaborate with grassland management entities on shared interests in grassland management or research.

Hydrogeology of the Tully Valley and characterization of mudboil activity, Onondaga County, New York

Released June 02, 2026 12:30 EST

1996, Water-Resources Investigations Report 96-4043

William M. Kappel, Donald A. Sherwood, William H. Johnston

Mudboil activity in the Tully Valley, in central New York, is causing turbidity in nearby Onondaga Creek, where it has caused a bridge to collapse; it also has threatened or damaged other structures and has caused extensive land subsidence. Mudboil activity was intermittent from its first reported appearance in the 1890's until the 1970's, when the rates of mudboil discharge and land subsidence began to increase. Historically, the water discharged from mudboils was reported as fresh, but chemical analyses in the late 1970's indicated an increase in specific conductance and chloride concentration.

Mudboil discharge is driven by artesian pressure in unconsolidated sediments that are confined by a 60-foot layer of silt and red clay. This process, once begun, has been self-propagating. Artesian pressures are about 20 feet above land surface over most of the valley floor but exceed 30 feet above land surface along Onondaga Creek where Rattlesnake Gulf and Rainbow Creek enter the Tully Valley. The source of artesian pressure is recharge from the Tully (Valley Heads) Moraine at the south end of the valley, and the alluvial fans of Rattlesnake Gulf and Rainbow Creek. The mudboils are found within a 300-foot-wide by 1,500-foot-long corridor along Onondaga Creek just upstream from the two alluvial fans, and in a 5-acre subsided area just west of that corridor.

Remediation efforts have entailed (1) diversion of flow from the tributary that feeds the subsided area, (2) installation of depressurizing wells at several locations, and (3) construction of a dam and settling impoundment to detain mudboil sediment that would normally discharge to Onondaga Creek. These efforts have been partly successful, but further work is needed to slow the mudboil activity, which is expected to persist in both areas. Mudboil activity is normally greatest during the early spring and late fall, when artesian pressures increase in response to seasonal ground-water recharge.

Suspended-sediment concentrations at the out-flow of the subsidence area ranged from 31,210 mg/L (milligrams per liter) in October 1991 to 17 mg/L after remediation efforts in the summer of 1993. Yearly average suspended-sediment loads to Onondaga Creek from the subsidence area for water years 1992, 1993, 1994, and 1995 were 29.8, 9.75, 1.41, and 1.80 tons per day, respectively. Sediment discharged from the mudboils initially was 30 to 60 percent clay and 80 to 100 percent silt-sized or smaller sediment, and the sand fraction never exceeded 20 percent. After the remediation projects, 50 to 80 percent was clay, and nearly all sediment was silt size or smaller.

Analyses of water from upstream and downstream of the subsidence area, as well as from mudboil vents within that area, indicate that the source of water for some mudboils is a confined freshwater aquifer, whereas for others it is an underlying, brackish-water aquifer. Water from the freshwater aquifer has specific conductance values ranging from about 400 (μS/cm (microsiemens per centimeter at 25° Celsius) to almost 900 (μS/cm, dissolved chloride concentrations range from 37 to 430 mg/L, and dissolved-solids concentrations range from 215 to 463 mg/L. Specific conductance of water from the brackish-water aquifer ranges from 17,000 to 28,000 (μS/cm, chloride concentrations range from 2,000 to 7,100 mg/L, and dissolved-solids concentrations range from 4,200 to 12,800 mg/L.

The largest landslide in New York State in the last 75 years occurred at the foot of Bare Mountain, 1 mile downstream from the mudboil area, in April 1993 and was the fourth in a series of slides that have occurred at the base of this hill. Slope instability was reported as early as May 1990. After the slide, intermittent mudboil-like activity was observed at several springs within the backscarp of the slide; water from these springs ranged from fresh to brackish. The chemical similarity between water from some springs in the backscarp area and water in the lower (brackish) aquifer beneath the mudboil area may indicate a hydraulic connection between this aquifer and the surficial deposits.

Hydrologic changes in the valley during the last 100 years have been attributed to salt-solution mining in the upstream (southern) end of the valley. The removal of nearly 150 feet of salt from four evaporite beds in the Syracuse Shale of the Salina Group has caused the collapse of bedrock and unconsolidated deposits in and near the brine field, 3 miles south of the mudboil area. These collapses have created a hydraulic connection among bedding plane aquifers in the bedrock and increased the hydraulic connection with unconsolidated aquifers. The ground-water flow system after brine field closure in 1988 may have reached a new semiequilibrium, but mudboil activity will likely continue because artesian pressures remain. Whether mudboils were present before salt solution-mining began is unknown.

Trends in subdaily to daily rainfall in Florida, 1990–2022

Released June 01, 2026 09:17 EST

2026, Journal of Hydrometeorology (27) 847-865

Saira Haider, Michelle M. Irizarry-Ortiz, Jayantha T. Obeysekera, Ana C. Maran, Tarana Solaiman, Brett D. Johnston

Changing rainfall patterns and intensifying rainfall extremes affect urban infrastructure and can increase flash-flood risk. Understanding how climate change has altered rainfall can support state and local agencies as they adapt and build resiliency. In this study, rainfall data from 23 weather stations in Florida were used to examine temporal and spatial trends over the period 1990–2022. Subdaily to daily rainfall events of durations 1, 2, 3, 6, 12, and 24 h were examined. A variety of statistical methods were applied to examine annual and seasonal trends, including quantile regression, extreme value analysis, run theory using the Mann–Kendall test, Sen–Theil slope, and Poisson and negative binomial tests, and threshold exceedance rates using generalized additive models. Using subdaily rainfall data posed challenges, including equipment failures, limited documentation of the quality assurance and control process, and potential measurement interferences. Results indicated that over 1990–2022, there was a decrease in hourly rainfall extremes but an increase at moderate quantiles. Overall, the number of rainfall events increased, particularly at shorter durations, but the mean total rainfall per event decreased. Additionally, the annual number of daily rainfall extremes showed more decreases than increases.

Geothermal district energy systems coupled with seasonal underground thermal energy storage: A U.S. techno-economic screening by climate and geology

Released May 30, 2026 08:44 EST

2026, Renewable Energy (271)

Scott Mello, Hyunjun Oh, Whitney J. Trainor-Guitton, Ryan Cain Cahalan, Jeffrey D. Pepin, Erick R. Burns

In the United States, cooling-dominated commercial building loads can cause geothermal heat pump-based district energy systems to accumulate a long-term subsurface thermal imbalance, motivating the incorporation of seasonal underground thermal energy storage. We developed a transferable workflow to evaluate geothermal district systems that pair ground heat exchangers with seasonal underground thermal energy storage. Using standardized hourly loads for seven commercial buildings and a uniform cost framework, we simulated ten U.S. cities with a physics-based ground heat exchanger model, subsurface storage simulations, and economic assessment to isolate the roles of climate and hydrogeology. In cooling-dominated cities, underground thermal energy storage supplied the majority of annual cooling, cutting electricity use and summer peaks substantially while achieving levelized costs comparable to or below conventional chiller-boiler plants. In cooler climates, the storage share shrunk, required borefield size and costs rose, and levelized cost of energy increased nearly linearly with declining underground thermal energy storage fraction, indicating storage fraction as the primary economic lever. Sensitivity analysis showed capital risk dominated by borefield drilling and surface heating, ventilation, and air-conditioning and piping, with underground thermal energy storage costs secondary. This workflow provides a transparent foundation for site-specific design and screening of next-generation geothermal district energy systems.

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

Released May 29, 2026 08:56 EST

2027, Fuel (427)

Aaron M. Jubb, Justin E. Birdwell, Leslie F. Ruppert, Martha Stokes, Ashton M. Wiens, Thomas Headen, Tristan G. A. Youngs

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

Co-occurrence of pesticides and pharmaceuticals and personal care products (PPCPs) across Zostera marina (common eelgrass) communities

Released May 29, 2026 08:51 EST

2026, Marine Pollution Bulletin (231)

Alexandra G. Tissot, Janet C. Niessner, Elise F. Granek, Kimberly Brown, Michelle L. Hladik

Anthropogenic pressures are driving changes in eelgrass communities, which are altering baseline conditions in estuarine environments. Field detections have validated the transport of land-sourced pollutants to aquatic systems; however, studies rarely sample concurrently for pesticides, and pharmaceuticals and personal care products (PPCPs) across environmental compartments. Moreover, studies on contaminant uptake by eelgrass and associated species are even more limited. In collaboration with the Confederated Tribes of the Coos, Lower Umpqua and Siuslaw Indians (CTCLUSI), this study collected samples of water, eelgrass, clams, and sediment at sites of Tribal significance in Southern Oregon to test for organic contaminants (i.e., herbicides and pharmaceuticals). Paired sampling was conducted for analysis by the CTCLUSI in tandem with the United States Geological Survey (USGS) in order for the Tribe to develop analytical standards for future sampling efforts. Ten pesticides and eight pharmaceuticals were detected across the four sites, with the highest number of overall detections (27) at the Florence Marina site. The insecticide bifenthrin was most frequently detected across all media (0.012–1.565 μg/g organic carbon in sediment, 2.7–30 ng/g in organismal tissue) and the anti-diabetic agent metformin was the most detected PPCP in clam tissues (1.33–3.78 ng/g). Pesticides and PPCPs were observed to co-occur in eelgrass habitats, with numerous pesticide detections across media types. These findings demonstrate numerous routes of exposure for estuarine organisms which could be addressed with pharmaceutical disposal strategies or pesticide use restrictions near these habitats.

Range-wide relative abundance of the Appalachian grizzled skipper (Pyrgus centaureae wyandot) in the Eastern United States

Released May 28, 2026 10:15 EST

2026, Open-File Report 2026-1017

Nimish B. Vyas, Jennifer Selfridge, David Cuthrell, Robert Somes, Erin White, Judith Ratcliffe, J. Merrill Lynch, Laurie Hamon, Eileen Wyza, Betsy Leppo, Pete Woods, Anthony Tur, Donovan Drummey, Kathryn Nolan, Ellison Orcutt, Andrew Rapp, Leah Card, Jakob Goldner, Susan Olcott

The U.S. Fish and Wildlife Service has designated the Pyrgus centaureae wyandot (Appalachian Grizzled Skipper [AGS]) to be at-risk, based on its declining populations and the lack of information on its status. The objective of this study was to complete range-wide surveys to locate extant AGS colonies and to quantify the number of AGS observed at each location. From 2021–24, 284 surveys were done in 25 unique (that is, distinct) counties in 8 States in the Eastern United States — Maryland, Michigan, New York, North Carolina, Ohio, Pennsylvania, Virginia, and West Virginia. We found AGS in only two counties: Alleghany County, Virginia, and Greenbrier County, West Virginia. AGS were observed 180 times in these two counties. Our results can inform U.S. Fish and Wildlife decisions about where and how future AGS conservation efforts can be implemented.

Effects of repeat prescribed burning in dry coniferous forests in national parks of California

Released May 28, 2026 08:57 EST

2026, Fire Ecology (22)

Phillip J. van Mantgem, Micah C. Wright, Calvin A. Farris, Eamon Engber, Emma J. McClure, Anthony C. Caprio, MaryBeth Keifer

Background

Prescribed fire is a common approach to reduce fuels and mitigate fire hazards. The accumulation of live and dead fuels following initial treatment means that repeated application of prescribed fire could be used to maintain this benefit. However, the effect of repeated prescribed fires is not well documented in many dry coniferous forests in the western United States. Here, we present observations of changes in live trees and surface fuels following two prescribed fires in dry coniferous forests in national parks of California.

Results

Changes in forest structure and accumulation of surface fuels were similar over time following initial-entry and second-entry fires. An exception was that repeated fires were associated with substantial reductions in stem density. There were smaller changes in live tree basal area and stem biomass.

Conclusions

Our results indicate that following initial-entry fires, subsequent burning maintained reductions in surface fuel loads without major inadvertent losses of live tree basal area and stem biomass, implying the survival of large trees.

Future water constraints on United States lithium mining under climate change

Released May 28, 2026 08:43 EST

2026, Communications Earth & Environment

Jenna Nicole Trost, Nedal T. Nassar, Jennifer B. Dunn

Lithium is necessary for low-carbon technologies that combat climate change, but lithium extraction is water-intensive. Changes in temperature and precipitation arising from climate change are altering water distribution, which could further strain supplies for new mines and industry, farms, and households. Here we explored how climate change, water use, and mining siting could impact lithium mining in the United States. We analyzed whether there would be sufficient water available to support the single existing and 22 proposed U.S. lithium mines at mid-century under four socioeconomic-climate scenarios and five climate models. Though dependent on socioeconomic-climate scenario, climate model, and lithium deposit type, available water supply in most subbasins would likely be unable to support new mines’ water demands, or even non-mining water demands from other sectors. Water scarcity could hinder the ability of the United States to produce enough lithium to meet domestic demand thereby necessitating higher imports.

Geochemical, mineralogical, and isotopic evidence for multi-stage genesis of the Hicks Dome REE + Y-HFSE-fluorite deposit, Illinois, USA

Released May 27, 2026 09:02 EST

2026, Ore Geology Reviews (194)

Julia A. McIntosh, Allen K. Andersen, Mitchell M. Bennett, Jay M. Thompson, Craig A. Johnson, Albert H. Hofstra, Laurence Nuelle

Hicks Dome hosts breccias enriched in rare earth elements (REE), Y, Th, F, Ba, Ti, Nb, and Be, alongside spatially associated lamprophyre dikes (ca. 271 Ma). Hicks Dome is located within the Illinois–Kentucky Fluorspar District, which hosts fluorite, Pb–Zn, and barite resources. This study investigates the genetic relationships between Hicks Dome mineralization in breccias, alkaline magmatism, and Illinois–Kentucky Fluorspar District mineralization. Lamprophyre dikes are light REE–enriched with chondrite-normalized abundances decreasing from La to Lu. The Host Breccia exhibits middle and heavy REE–enriched patterns that mirror those of the principal REE–Th host minerals, including fluorapatite, xenotime, and thorite. Textural evidence suggests recrystallization of phosphates, sulfates, and Ti–Nb oxides in the Host Breccia. U–Pb geochronology constrains multiple mineralizing events, with ages of 277 ± 18 Ma from low-Th apatite interpreted as main-stage mineralization, and 121.6 ± 9.7 Ma from high-Th apatite indicating later overprinting. O–H–C stable isotope data provide evidence for multiple stages of fluid-rock interaction and fluid mixing: (1) early magmatic fluids dissolved limestone country rock, (2) mixing between magmatic fluids and basinal brines led to main-stage mineralization in the Host Breccia, and (3) late-stage mineralization occurred following mixing of meteoric water and basinal brine. These results indicate that heavy REEs, high field strength elements, and fluorine precipitated proximal to its alkaline magmatic source because of fluid–rock interactions and fluid mixing. Subsequent fluid mixing drove late-stage recrystallization and additional fluorite formation, a process that may be similar to mineralization in the Illinois-Kentucky Fluorspar District.

Status and trends of pelagic and benthic prey fish populations in Lake Michigan, 2025

Released May 27, 2026 08:46 EST

2026, Report

Ralph W. Tingley III, Timothy P. O’Brien, Charles P. Madenjian, Peter C. Esselman, Patricia Dieter, Kristy Phillips, Ben Turschak, Dale Hanson, Steven A. Farha

Fall bottom trawl (fall BT) and lakewide acoustic (AC) surveys are conducted annually to generate indices of pelagic and benthic prey fish densities in Lake Michigan. The fall BT survey has been conducted each fall since 1973 using 12-m trawls at depths ranging from 9 to 110 m at fixed locations distributed across seven transects; this survey estimates densities of seven prey fish species [i.e., Alewife (Alosa pseudoharengus), Bloater (Coregonus hoyi), Rainbow Smelt (Osmerus mordax), Deepwater Sculpin (Myoxocephalus thompsonii), Slimy Sculpin (Cottus cognatus), Round Goby (Neogobius melanostomus), Ninespine Stickleback (Pungitius pungitius)]. The AC survey has been conducted each late summer/early fall since 2004 (except 2020). The 2025 AC survey consisted of 26 transects [470 km total (292 miles)] covering bottom depths ranging from 5 to 259 m and 44 midwater trawl tows at 1.4 to 82.4 m fishing depth; this survey estimates densities of three prey fish species (i.e., Alewife, Bloater, and Rainbow Smelt). The data generated from these surveys are used to estimate various population parameters that are, in turn, used by state and tribal agencies in managing Lake Michigan fish stocks.  

For the AC survey, total biomass density of prey fish equaled 9.3 kg/ha, continuing a recent trend of index values above the long-term average of 5.4 kg/ha. For the fall BT, total biomass density of prey fish equaled 3.4 kg/ha, close to values observed since 2014 and well below historic numbers and those observed earlier in the 2000s. Over the period both surveys have been conducted (2004-2025), the total biomass density index had trended downward in the fall BT through the mid-2010s and appears to have stabilized at low values, while the AC survey biomass density index has remained relatively stable over the time series.  

Mean biomass of yearling and older (YAO) Alewife was 4.3 kg/ha in the AC survey and 0.45 kg/ha in the fall BT. Since 2014, annual survey results suggest that the catchability of YAO Alewife for the fall BT is substantially lower than the AC survey. The 2025 AC survey YAO Alewife biomass density estimate was 57% higher than the average from 2004-2024. The Alewife population of Lake Michigan appears to be composed mostly of young fish and the proportion of age-4 and older Alewife was ~5% in both surveys. Age-0 Alewife numeric density from the AC survey was 259 fish/ha in 2025, lower than the long-term mean (487 fish/ha). Biomass density of large (≥120 mm) Bloater was 3.5 kg/ha in the AC survey and 1.9 kg/ha in the fall BT. The density of small (<120 mm) Bloater was 540 fish/ha in the AC survey, the second highest value in the time series.  Meanwhile, small Bloater density estimated in the fall BT was only 6.1 fish/ha. Biomass density of large Rainbow Smelt (≥90 mm) was 0.69 kg/ha in the AC survey and 0.04 kg/ha in the fall BT survey. Numeric density of small (<90 mm) Rainbow Smelt was 541 fish/ha in the AC survey, the highest value in the time series, and 41 fish/ha in the fall BT. All four prey fish species indexed only by the fall BT had below-average biomass densities. Deepwater Sculpin biomass density was 0.21 kg/ha, which makes 15 of the past 16 years with biomass <1 kg/ha. Slimy Sculpin was estimated to be 0.03 kg/ha, an order of magnitude lower than the long-term average from the fall BT. Round Goby biomass density was 0.44 kg/ha and Ninespine Stickleback density was 0.20 kg/ha, the highest value since 2007. 

Moment magnitude for small earthquakes in the Delaware basin of west Texas and southeast New Mexico, USA

Released May 27, 2026 07:42 EST

2026, Seismological Research Letters

Sydney Gable, Yihe Huang, David R. Shelly, Justin L. Rubinstein

The Delaware Basin region of west Texas and southeast New Mexico has become one of the most prolific regions of seismic activity in the continental United States due to widespread hydraulic fracturing and wastewater disposal injection. In response to the increased number of earthquakes in this region, rapid and accurate characterization of earthquake sources is necessary to understand the evolution of seismic activity and level of seismic hazard associated with these earthquakes. This study re-evaluates earthquake magnitudes, estimating moment magnitude (MW) for small earthquakes in the Delaware Basin using 1) moment-rate spectra derived from S-wave coda envelopes, and 2) a relative magnitude method that relies exclusively on the ratio of waveform amplitudes between highly correlated waveform pairs. The coda-envelope method produces accurate MW estimates for small earthquakes (M 1.5 – 3) that are consistent with independent, waveform modeled moment magnitudes for events with MW > 3. Using the relative amplitudes method to extend these MW magnitudes to many other events, we successfully provide relative moment magnitude (MW,rel) values for 81% of the Texas Seismological Network catalog in the Delaware Basin region, and 45% of the USGS Induced Seismicity Project’s catalog of events in southeast New Mexico. The adoption and integration of the calibrated MW,rel method with current magnitude estimation methods offers valuable insights into the relationships between local and moment magnitude and will contribute to improved characterization of widespread induced seismicity.

Continuous and high-resolution longitudinal profiles of the water surface and riverbed elevation for 282 miles of the Colorado River from Lees Ferry to Pearce Ferry, Arizona, 2021

Released May 26, 2026 10:00 EST

2026, Scientific Investigations Report 2026-5010

Shannon L. Sartain, Matthew A. Kaplinski, Keith Kohl, Katherine A. Chapman, Nathaniel D. Bransky, Joel B. Sankey, Paul E. Grams

Longitudinal profiles of water surface and riverbed elevations capture key geomorphic characteristics that can be affected by water infrastructure and natural processes. Continuous water surface profiles of the Colorado River in Grand Canyon, a river influenced by two of the largest dams in the United States, have been measured infrequently. The water surface profile was first measured in 1923, 13 years before the completion of Hoover Dam, which impounded water into western Grand Canyon, and 40 years before the completion of Glen Canyon Dam, which affected streamflow and sediment supply for all of Grand Canyon. The water surface profile was next measured in 2000, 37 years after the completion of Glen Canyon Dam, although this profile did not include the segment affected by Hoover Dam. A continuous profile of riverbed elevations has never been published. Here, we present the first complete, coupled water surface and riverbed elevation profiles, collected in 2021 during a period of steady releases from Glen Canyon Dam. The profiles were constructed from positions and elevations measured by boat-based global navigation satellite systems and from bathymetry collected by multibeam sonar. Data collected by boat were supplemented by data from a photogrammetry-derived digital surface model that was created from concurrently collected aerial images. Independent measurements made by conventional total stations referenced to a common geodetic control network were used to evaluate accuracy of all measurements. The final water surface and riverbed elevation profiles improved the accuracy and precision reported for previous profiles. In this study, the mean absolute vertical accuracy of water surface elevations was 0.07 meter for 85 percent of river miles and 0.19 meter for 11 percent of river miles. For the remaining 4 percent of river miles, water surface elevations were interpolated between measured values. The profiles reported herein can be used for current assessment of Colorado River geomorphic conditions, quantification of changes in the river over time, and predictive modeling of river resources for potential future management scenarios.quantification of changes in the river over time, and predictive modeling of river resources for potential future management scenarios.

Streamflow and surface-water presence data availability across the conterminous United States: A review for headwater systems

Released May 26, 2026 08:49 EST

2026, Hydrological Processes (40)

Roy Sando, Kristin Jaeger, Christa Kelleher, John C. Hammond, Jay R. Christensen, Catalina Segura, Heather E. Golden, Frederick Y. Cheng, Admin Husic, C. Nathan Jones, Charles R. Lane, Li Li, D. Tyler Mahoney, Hillary McMillan, Adam N. Price, Erin C. Seybold, Adam Ward, Margaret Zimmer, Steven James Pestana

Water is essential for life on Earth, supporting ecosystems, human health, and economic activities. Hydrology relies on observational data, and this paper discusses regional and national datasets for the conterminous United States (CONUS) publicly available as of 2023, focusing on headwaters, defined as first- and second-order streams at 1:24000 scale. It identifies 72 primary and secondary datasets and 11 repositories and argues how better integration and accessibility of hydrological data can improve research. The paper distinguishes between datasets where streamflow was the primary data collection objective and those where it was secondary. This distinction highlights opportunities to consider data from efforts peripheral to hydrology but is still useful for understanding hydrologic conditions. The analysis reveals that out of about 118 000 active and inactive stream observation sites, about 6.6% and 25% are located on first- and second-order streams, respectively. This indicates a substantial data gap for headwater systems, which account for over 77% of stream length in CONUS. Federal agencies manage 72% of hydrologic monitoring sites across all stream orders, but only 34% of these are in headwater systems. Academic institutions operate about 2% of sites, with almost half (48%) in headwater systems, focusing on ecosystem research. State agencies also operate about 2% of sites, primarily on larger systems, with 19% on headwaters. Additionally, 23% of sites are managed by multiple agencies. Spatial patterns further reveal pronounced disparities among physiographic regions. Eastern and coastal provinces show relatively dense monitoring, while central and western regions show sparse coverage. These gaps reflect historical priorities, logistical constraints, funding limitations, and the high cost of continuous instrumentation. To address biases in monitoring networks, data collection could be enhanced with low-cost monitoring, community science, and remote sensing technologies. This study also notes the benefits of long-term monitoring and prioritizing retention of streamgages with longer records.

Effects of warming on growth and leaf colonization by litter mat-forming fungi in a wet tropical forest in Puerto Rico

Released May 24, 2026 07:56 EST

2026, Biotropica (58)

Ari E. Puentes, D. Jean Lodge, Deyaneira A. Ortiz-Iglesias, Tatiana Barreto-Vélez, Laura C. Rubio-Lebrón, Hieu P. Chu, Christine S. O'Connell, Sasha C. Reed, Tana E. Wood

Wet tropical forests are experiencing rising temperatures and increased frequency and intensity of extreme climatic events, such as cyclones, which can increase rates of soil erosion and surface runoff. Fungal litter mats, formed by agaric decomposer fungi, play a crucial role in stabilizing slopes, preventing erosion, and aiding nutrient cycling; however, little is known about how warming affects litter mat growth and function. We investigated two litter mat-forming fungi, Gymnopus johnstonii and Marasmius aff. crinis-equi, in warmed (+4°C above ambient) and control plots in the Luquillo Experimental Forest, Puerto Rico. Growth and time-to-leaf colonization were monitored over 6 weeks in spring (both species) and summer (G. johnstonii only). We hypothesized that warming would inhibit fungal mat growth and slow leaf colonization, particularly for G. johnstonii since it is drought sensitive. As expected, warming significantly reduced relative growth rates (RGR) in spring, though M. aff. crinis-equi showed slightly higher RGR than G. johnstonii. Leaf colonization was also delayed by 22% in warmed plots, with M. aff. crinis-equi colonizing leaves 4.3 times faster than G. johnstonii. There were significant seasonal differences in response to warming for G. johnstonii, with warming increasing RGR during the consistently wetter summer sampling period. Overall, warming led to significant inhibition of leaf colonization when conditions were dry, whereas there was a trend toward increased colonization in warm and wet conditions. Our findings suggest that warming, combined with drier conditions, is likely to suppress drought-sensitive fungal mat growth, reducing their ability to prevent nutrient and soil loss via erosion.

Lowland tropical forests remain a methane sink under warming and long-term hurricane disturbance recovery

Released May 23, 2026 09:14 EST

2026, Agricultural and Forest Meteorology (386)

Gabriele Larocca Conte, Lucia Zuvela, Rachel Cruz-Perez, Tatiana Barreto-Vélez, Nibia Becerra-Santillan, Sophia F. Campbell, Hieu P. Chu, Trung Dam, Iana F. Grullón-Penkova, Miriam Kleit, Deyaneira A. Ortiz-Iglesias, Laura C. Rubio-Lebrón, Molly A. Cavaleri, Sasha C. Reed, Debjani Sihi, Tana E. Wood, Christine S. O'Connell

Methane (CH4) is a potent greenhouse gas, and tropical forests account for roughly one–third of global atmospheric CH4 uptake by soils. Projected warming and more frequent hurricanes in these ecosystems may alter soil CH4 sink strength, as warmer and wetter soils enhance methanogenesis activity. We measured soil CH4 and CO2 efflux during the calendar summer months of 2023 and 2024 alongside continuous records of soil moisture, soil and air temperature, and precipitation in an in–situ warming experiment (TRACE) located in a lowland tropical forest in Puerto Rico, six to seven years after Hurricanes Irma and Maria (2017). The realized warming (∼1.95°C) enhanced soil respiration only in summer 2023 (p < 0.05), but net soil CH4 uptake was invariant in both campaigns (p > 0.05). Instead, sampling day and between–plot variability explained soil CH4 dynamics much more than treatment contrasts. Importantly, CH4 uptake was consistently coupled to CO2 efflux, suggesting tight linkages between methanotrophic and heterotrophic activities. Between treatments, CH4 and CO2 responses to soil temperature variation were less sensitive in warmed plots, which may suggest weak metabolic upregulation under elevated temperatures. Together, these findings indicate that lowland tropical soils remain CH4 sink even under warming and years after hurricane disturbance, with CH4 dynamics driven more by spatial and temporal variability than experimental warming. Long–term, high–resolution monitoring integrating soil biogeochemistry and microbial processes will be critical to determine whether the observed net CH4 uptake signal represents a sustainable or transient response under continued warming and disturbance.

Biochar modulates the dynamics of legacy nutrients in enhancing soil health and crop productivity

Released May 22, 2026 09:17 EST

2026, Land (15)

Manish Kumar, Shiv Bolan, Rakesh Kumar, Juhi Gupta, Dingjiang Chen, Hao Wu, Sarah M. Stackpoole, Nitika Chandel, Santanu Mukherjee, Manoj Chandra Garg, Srinithi Mayilswami, Kadambot H. Siddique, Nanthi Bolan

Most major crops in agricultural soils exhibit relatively low nutrient use efficiency for nitrogen (N), phosphorus (P), and potassium (K), often necessitating supplemental nutrient inputs to achieve sustainable yields. Furthermore, the increasing use of biowastes such as compost, manure, and biosolids, which frequently have nutrient ratios that do not match crop requirements, has contributed to excessive nutrient inputs and subsequent accumulation in soils. This situation has been further exacerbated by intensive farming practices involving multiple cropping cycles per season. Overuse of nutrients causes them to accumulate in the soil, creating a legacy nutrient pool. The application of biochar as soil amendment is considered a potential strategy to control legacy nutrients dynamics. The current review inspects the possible value of biochar in modulating legacy nutrient reserves in the soil, thereby increasing the bioavailability of nutrients and improving crop yield. This review discusses the search scope and synthesis approaches for the bibliometric methodological component through rigorous screening process (Preferred Reporting Items for Systematic Reviews and Meta-Analysis (PRISMA)), focusing on journal articles published in last 20 years that specifically address legacy nutrient management. The significance of the economic and environmental effects of legacy nutrients and the insufficient knowledge of how biochar application influences nutrient dynamics in soil highlight the necessity for additional research to address current gaps.

Tracking toxins: A pilot investigation of cyanotoxins in north-central Tennessee’s surface waters and wells

Released May 22, 2026 08:56 EST

2026, Toxins (18)

Kristi Lynn Hill, Andrea Jaegge, Devin M. Moore, Thomas D. Byl

Cyanobacterial toxins (cyanotoxins) threaten aquatic ecosystems and human health, yet the factors influencing their production and distribution in freshwater remain unclear. In north-central Tennessee, nutrient-rich runoff from agricultural and urban areas, combined with a karst landscape that supports drinking and recreational water use, heightens the need to understand cyanotoxin behavior. To examine cyanotoxin patterns, the U.S. Geological Survey and the Tennessee Department of Environment and Conservation monitored 18 sites, including two wells under the influence of surface water, every two weeks from September 2022 to November 2024. At least one cyanotoxin was detected at all sites, with the highest concentrations in deep reservoirs and lower levels in shallow systems. Most detections occurred during summer and fall, aligning with high temperatures and rapid-onset drought. Statistical analysis indicated that increased specific conductivity and pH raised the likelihood of detecting total microcystin, likely resulting from drought conditions and nutrient-laden runoff. Additionally, dissolved microcystin showed an inverse relationship with Cumberland River water levels, and principal component analysis showed that Secchi depth, chlorophyll a, pH, temperature, and conductivity explained most water quality variability. These results help increase understanding of cyanotoxin distribution and associated water quality conditions during detections to guide future freshwater cyanotoxin monitoring studies.

Waves, watersheds, and sediment in a coral reef embayment: Towards parsimonious models of accumulation and composition

Released May 22, 2026 08:55 EST

2026, Estuarine, Coastal and Shelf Science

Trent Biggs, Alex Messina, Curt D. Storlazzi

High sedimentation rates can damage coral reef ecosystems. Sedimentation rates are controlled by both sediment loads from watersheds and resuspension by waves and associated circulation patterns, but the outcomes are system specific and difficult to predict. The percent terrigenous (non-organic and non-carbonaceous) material in sediment is also often used as an indicator of watershed influence, but its dynamics are poorly understood. Sediment accumulation rates, particle size, and percent terrigenous were monitored quasi-monthly for one year (March 2014-April 2015) at nine sites in a coral reef-fringed embayment in American Samoa, where an aggregate quarry had increased sediment loads to the coast but mitigation reduced loads during the monitored period. Gross and net sediment accumulation rates were measured using sediment traps and SedPods (pods), respectively. Gross accumulation rates exceeded thresholds for impacts on coral health during at least one collection period at most sites, with more exceedances on the northern reef where water residence times and sediment availability are higher and corals show signs of sediment stress. Percent terrigenous of coarse sediment was higher in the traps and pods compared with the surrounding benthic sediment, indicating that some of the terrigenous sediment was advected through the bay without accumulating on the reef. The 95th percentile of hourly wave energy density (E95) taken from a global wave model (WaveWatch 3) was the best predictor of gross accumulation rates of both total and carbonate sediment in a log-log regression at most (n = 6) sites (R2 range 0.72-0.92), indicating a strong role of resuspension of benthic sediment. Gross accumulation rates of terrigenous sediment were not correlated with E95 and only correlated with SSY at the site nearest the stream mouth, indicating that most terrigenous sediment was not from resuspended benthic material but rather from a consistent watershed source. Percent terrigenous decreased with increasing wave energy due to high accumulation rates of carbonates during periods of high wave energy. Detection of the impact of sediment mitigation at the quarry on sediment accumulation was complicated by low wave energy in the period following mitigation. The use of gross accumulation rates and percent terrigenous as indicators of the magnitude and sources of sediment accumulation over time needs to account for wave-induced resuspension, which can be modelled with a simple power function using inputs from a global wave model.

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

Released May 22, 2026 08:54 EST

2026, Report

Pamela L Nagler, Jennie N. Duberstein, James Broska, Kamel Didan, Myles B. Traphagen

We produced a 30-m resolution binational land cover map of Bird Conservation Region 33 (BCR 33) for the U.S. North American Bird Conservation Initiative. The region covers large portions of the Sonoran and Mojave Deserts. The map can support the U.S. Fish and Wildlife Service (FWS) Migratory Bird Program’s recovery planning efforts and constitutes the first known binational land cover dataset spanning sections of the United States–Mexico border and using a consistent classification system for both countries. The mapped region includes 152 distinct land cover classes, covering a total area of 38,421,453 ha (148,345 mi2), of which 13,148,345 ha (52,706 mi2) are located in Mexico and 24,770,640 ha (95,639 mi2) in the United States.

We primarily used Landsat 8 (OLI) imagery, supplemented by limited ground surveys from two field campaigns, drone-based aerial data, and existing vegetation classification frameworks from both countries. The classification applied a data-fusion approach integrating 30-m Landsat 8 imagery, decadal phenology metrics from vegetation indices, and a random forest model trained mainly with datasets from a comprehensive national mapping project from the U.S. Geological Survey (USGS) GAP Analysis Project (GAP) and federal wildland fire agencies’ Landscape Fire and Resource Management Planning Tools (LANDFIRE) (GAP/LANDFIRE) [United States side] and the National Institute of Statistics and Geography (INEGI) [Mexico side] as well as land cover maps and opportunistic open-access and field observations.  

Mapping of the full BCR 33 region was carried out in two phases: 1) Phase I, the prototype map, covered a smaller portion of the transboundary area and identified 31 land cover classes, and 2) Phase II, the full BCR 33 map (refer to Figure 1), which resulted in 152 land cover classes. Using a Random Forest classifier, we achieved an overall prediction accuracy of 92% for the Phase I map and 87% for the Phase II full region map. This slight decrease can be attributed to working on a larger, more complex area with a greater number of land cover classes. No formal validation was conducted, aside from using a subset of the collected field observations and training data to assess model performance during and after training. The training sites were further verified using Google Earth (Google, 2026) imagery. Two undergraduate students who worked for over a year visually inspected imagery and open access public images to confirm each training site during model training using in-house developed, online, visual tools. A portion of this field training data was reserved for model validation, and the corresponding results are to be presented in later sections. 

The project developed an end-to-end, medium- and fine-resolution remote sensing–based data fusion mapping approach. This effort produced a map (Nagler et al., 2025) and the online tools to support a dynamic, live, online map for visualizing the transboundary vegetation communities in BCR 33. The toolset is currently hosted by the University of Arizona (UofA) Vegetation Index and Phenology (VIP) Lab to support FWS partners (https://vip.arizona.edu/viplab_data_explorer?LCM_BCR33). The online map is designed to allow rapid updates using new training, validation, or correction data, making it dynamic and maintainable. 

The approach we took established a framework for rapid updating and correction of land cover maps, as the model can be quickly retrained with new field observations, updated training data, or other sources. This enables dynamic mapping and change detection of the region’s vegetation. This framework is an advance in data fusion and crowdsourced mapping of complex, vulnerable regions, providing support to regional stakeholders and the wider user community. 

This transboundary map can inform the protection, conservation, and restoration of vegetation, habitat, and ecosystems, particularly for threatened and endangered species across the two nations using consistent and harmonized binational mapping systems. Beyond supporting land management decisions and stakeholders in the transboundary desert ecoregions, this BCR 33 mapping effort establishes a foundation for future rapid, low-cost, cross-border land cover mapping that can benefit and advance ecosystem management. 

Rearing method has limited effect on post-release movement of reintroduced age-0 Lake Sturgeon

Released May 22, 2026 08:41 EST

2026, Canadian Journal of Fisheries and Aquatic Sciences (83)

Jorden R. McKenna, Justin A. Chiotti, Christopher S. Vandergoot, Richard Kraus, Matthew D. Faust, Zak J. Slagle, Eric J. Weimer, Matthew D. Cross, William D. Hintz

Overfishing, habitat loss, and pollution caused the extirpation of Lake Sturgeon (Acipenser fulvescens) throughout much of the Great Lakes. A Lake Sturgeon reintroduction program using two rearing strategies began in 2018 in the Maumee River, a tributary of Lake Erie. We assessed the movement of streamside or traditionally reared age-0 Lake Sturgeon using acoustic telemetry to determine if rearing strategy affected river residency, movement, and the habitat area used. Tagged sturgeon generally left the Maumee River for Lake Erie on average 3–47 days after stocking and spent most of their time in the western basin of Lake Erie. The majority of sturgeon moved through nearshore areas along the south shore of Lake Erie. While we found no differences in post-stocking movements or habitat area used between the two rearing strategies, understanding how older life stages respond to rearing strategy is needed. Adding upstream stocking sites, using source water to raise eggs or larvae if excessive straying becomes evident, and increased acoustic receiver coverage are options to facilitate and evaluate successful recovery of Lake Sturgeon.

Effects of fipronil bait pellets on two cricetid species: Potential implications for plague mitigation and wildlife conservation

Released May 22, 2026 07:47 EST

2026, International Journal for Parasitology: Parasites and Wildlife (30)

David A. Eads, Marc R. Matchett, Travis M. Livieri, Richard A. Bowen, Airn E. Hartwig, Stephanie Porter, Mary L. Wright, Jason Fly, Madisen Hartlaub, Phillip Dobesh, Paul Roghair, Eddie Childers, John P. Hughes, Michelle L. Hladik, Gregory P. Dooley, Brian J. Smith, Rachel A. LaCasse, Kristy Bly, Dean E. Biggins

We evaluated the effects of fipronil bait pellets on two cricetids that commonly occupy colonies of black-tailed prairie dogs (Cynomys ludovicianus; BTPDs): western deer mice (Peromyscus sonoriensis) and northern grasshopper mice (Onychomys leucogaster). In one experiment, bait pellets (0.96 mg fipronil/bait) were applied at 75 baits/ha to three 1.44-ha plots on a BTPD colony. Mouse abundance declined by 70% from before to 6-10 d after treatment. In a second experiment, bait pellets (0.46 or 1.52 mg fipronil/bait) were applied at 125 baits/ha to four plots (0.85-1.86 ha) on two BTPD colonies; two non-treated plots were baselines (1.09 and 2.06 ha). From before to 11-15 d after treatment, mouse abundance declined by 51%- 67% on the treated plots vs. a decline of 9% on the non-treated plots. Mouse survival from before to 11-15 d after treatment was 51% lower on the treated plots. In a third experiment, bait pellets (0.84 mg fipronil/bait) were applied at 125 baits/acre on two 1.44-ha plots on a BTPD colony; two 1.44-ha non-treated plots were baselines. Mouse survival from before to 30-44 d after treatment was 45% lower on the treated plots; the abundance of deer mice on the treated plots remained similar from before to 30-44 d after treatment, perhaps due to juvenile recruitment and/or immigration. In a laboratory experiment, 33 deer mice offered one bait pellet (0.86 mg fipronil/bait) consumed 27% of their bait, on average (range = 0-100%). Over 3 d, deer mouse mortality was estimated at 53%; mortality increased with fipronil dose, which averaged 11 mg fipronil/kg body mass (range = 3-46 mg/kg). Brain samples were available from 31 deer mice; all tested positive for fipronil sulfone, the primary mammalian metabolite of fipronil, at 19 to 61,205 ng fipronil sulfone/g. Additional experiments could determine if these findings scale up to larger landscapes.

Critical minerals memory match game

Released May 21, 2026 18:00 EST

2026, General Information Product 266

Danielle A. Olinger

An educational information packet about the 2025 List of Critical Minerals, which includes a memory match game about select critical minerals and how they are used.

Indicators of mercury concentration in Lake Trout: Can fish location and appearance provide information to anglers to reduce their exposure?

Released May 21, 2026 09:34 EST

2026, North American Journal of Fisheries Management (46) 663-679

Sarah M. Laske, Daniel Young, Krista K. Bartz, Vanessa R. von Biela, Michael P. Carey

Objective

People are exposed to mercury (Hg) through the consumption of fish. State and federal governments provide broad, often-generalized food safety guidance to reduce exposure; however, numerous rural fishing areas lack testing and location- or species-specific guidance. The aim of this study was to provide tangible, visible, or easily measured characteristics of Lake Trout Salvelinus namaycush that could convey information on Hg exposure to people harvesting and consuming fish where no location-specific guidance exists.

Methods

We investigated potential indicators of Lake Trout total Hg (THg) concentrations in muscle across 10 lakes in Alaska's national parks. Potential indicators, including lake, lake zone (i.e., littoral, pelagic, profundal), fish length, head size, body condition, and general appearance, were evaluated by competing linear mixed-effects models.

Results

Lake Trout THg concentrations ranged widely from 22 to 1,306 ng/g wet weight. Much of the variation (48%) in THg concentrations was attributed to differences among individual lakes, but the interaction of the fish's lake zone, body length, and head size accounted for an additional 21%. Predicted THg concentrations increased with Lake Trout length and head : body proportion, but the rate of THg concentration increase with length varied by head : body proportion and lake zone.

Conclusions

Given the overwhelming evidence of high lake-to-lake variability in Lake Trout THg concentrations, we find support for use of lake-specific guidance when data are available. When lake-specific THg concentrations are not available, the best potential way to reduce exposure is to harvest and consume Lake Trout with mean predicted THg concentrations that are within state and federal safe consumption guidelines. This included Lake Trout from surface waters (i.e., pelagic or littoral zone) that are ≤70 cm in length; if harvesting fish from deep waters (i.e., profundal zone), lower THg concentrations were found in Lake Trout with heads ≤25% of their body length. The indicators—lake zone, length, and head size—of Lake Trout THg concentrations can provide harvesters with additional information in the absence of data for specific lakes.

Ecology of reintroduced Rocky Mountain bighorn sheep in Dinosaur National Monument

Released May 21, 2026 09:27 EST

2026, Wildlife Monographs

Sarah Louise Carroll, Elizabeth P Flesch, Salix Scoresby, Emily Spencer, Rachel S. Crowhurst, Clinton W. Epps, Nathan L. Galloway, William Michael Janousek, Tabitha A. Graves

Translocations have been widely used to restore and conserve bighorn sheep (Ovis canadensis) populations in North America. Some translocations have been successful, but many populations remain small and genetically isolated. Population structure can influence the viability and long-term success of reintroductions. Social ungulates often function as interconnected subpopulations (metapopulations); however, few studies evaluate subpopulation sizes, connectivity, and genetic diversity within metapopulations. To address this gap, we conducted a comprehensive study of a reintroduced Rocky Mountain bighorn sheep (Ovis canadensis canadensis) population in Dinosaur National Monument in Colorado and Utah, USA, between 2006–2020. We analyzed global positioning system (GPS) radio-collar data, genetic samples, and results of health testing to evaluate abundance, distribution, genetic structure and diversity, habitat use, movement and connectivity, and presence of or exposure to respiratory pathogens. We integrated these analyses to evaluate the outcomes of a reintroduction effort that began in 1952, over 70 years ago, and to inform management decisions in Dinosaur National Monument. We also provide a framework for evaluating metapopulation processes, including a non-invasive approach that links genetic structure with Bayesian spatial capture-recapture analyses to estimate subpopulation sizes. Despite models indicating continuous suitable habitat, we found a spatially structured population with at least 4 subpopulations with constrained connectivity. Evidence from step selection and density analyses suggested that movement among subpopulations may be limited by semi-permeable barriers including rivers and human disturbance, which could contribute to maintenance of spatial structure over time. In 2006, antibody to Mycoplasma ovipneumoniae was detected in all geographically and genetically distinct subpopulations. Widespread clinical signs of disease and a confirmed exposure to M. ovipneumoniae in 2019 indicate a long-term disease challenge. Proximity to domestic sheep creates repeated opportunities for introduction of new M. ovipneumoniae strains. We estimated abundance in 2019 at 109 (95% CrI = 87–133), composed of subpopulations ranging from 18–39 animals (95% CrIs from 11–50). Genetic diversity was relatively high compared to other reintroduced and native Rocky Mountain bighorn sheep populations, which is likely a consequence of multiple translocations from different sources. Three of 4 subpopulation centers generally aligned with the locations of original translocation release sites. Persistence in the presence of pathogens may be facilitated by metapopulation structure and moderately high genetic diversity. Conversely, metapopulation structure can also facilitate pathogen persistence. Our approach offers a path to advance understanding of the population ecology of reintroduced bighorn sheep and can inform effective conservation and management of their populations.

Modeling the seasonality of wind-driven hydrocarbon waves in Titan’s polar lakes

Released May 21, 2026 09:25 EST

2026, JGR Planets (131)

Charlene E. Detelich, Una G. Schneck, Alexander G. Hayes, Milan Curcic, Rose Elizabeth Palermo, Andrew D. Ashton, J. Taylor Perron, Juan M. Lora, Jordan Steckloff

Titan, the only body in the solar system aside from Earth with standing liquids on its surface, has polar hydrocarbon lakes and seas. As Titan’s atmosphere generates light winds, there should be waves on the surface of these lakes and seas, yet, direct wave observations are scant. We introduce and use PlanetWaves, an open source 4D spectral wave model, to study Titan’s waves and create seasonal maps of wave shape and propagation on Ontario Lacus and Ligeia Mare. Titan’s modeled waves grow up to 30 times larger than terrestrial waves for the same wind speed, are seasonally present and are largest in the spring and summer when winds are strongest. Average daily winds almost never exceed the wave generation threshold of 0.5–0.7 m/s. Average storm winds (∼1.5 m/s) generate waves 15–48 cm in height with a period ranging 6–10.5 s while maximum storm winds (∼4 m/s) generate waves 2.7–3.2 m in height with a period up to 32 s. Titan’s waves become fetch-independent at ∼40 km for average storm winds occurring ∼1% of a Titan year and ∼100 kilometers for maximum storm winds occurring 2-3 times per Titan decade. On Ontario Lacus, storm winds blow nearly parallel to the eastern shore, potentially driving wave modification of the smooth eastern shoreline. On Ligeia Mare, waves rarely propagate toward a hypothesized wave modified shoreline suggesting that another process, such as tectonics, may contribute to a straight shoreline morphology.

The Great Lakes Geologic Mapping Coalition—Working collaboratively to understand the geology of the Great Lakes Region

Released May 20, 2026 13:25 EST

2026, Fact Sheet 2026-3010

Brianna Lopez, Jenna L. Shelton, Michael Marketti, Kate Ritzel, Brandon L. Graham

Introduction 

The Great Lakes Geologic Mapping Coalition (GLGMC), commonly referred to as the “Coalition,” is a partnership between the U.S. Geological Survey (USGS), the U.S. States of Illinois, Indiana, Michigan, Minnesota, New York, Ohio, Pennsylvania, and Wisconsin and the Canadian province of Ontario. The member States receive funding for geologic mapping work from the USGS National Cooperative Geologic Mapping Program (NCGMP), whereas Ontario participates as a nonfunded partner. The mission of the GLGMC is to produce three-dimensional (3D) geologic maps that depict unconsolidated sediments and near-surface bedrock in the Great Lakes region of North America. Geologic maps are the basis of most earth science investigations and help support resource exploration (energy, minerals, groundwater), natural hazard mitigation, infrastructure development, and land-use planning, all of which can be used to advance economic development and strengthen national security in the Great Lakes region. 

During the last few million years, the Great Lakes region has experienced repeated glacial advances and retreats, leaving behind extensive sediments, abundant natural resources, and widespread effects on the underlying bedrock geology (Swezey and others, 2022). Linked by shared histories of past glaciations, industrial agriculture, and legacy automotive, coal, steel, and manufacturing industries, the GLGMC member States collaborate to improve the understanding of the 3D distribution of the sediments overlying the region’s bedrock (fig. 1). Developing a comprehensive subsurface 3D framework of this glaciated terrain can provide earth science data to policymakers at all levels. These insights facilitate informed decisions on the exploration, use, and protection of vital resources, such as critical minerals, industrial materials, and aquifers, thereby supporting economic prosperity and the well-being of the citizens of this region.

Since its inception in 1998, the Coalition has completed more than 100 geologic mapping projects across the Great Lakes region. Each project aims to deliver geologic maps, 3D datasets, and other information that improves understanding of the geology of the Great Lakes region, with an emphasis on economic and water resources. Key deliverables include 3D geologic maps and models typically portraying sediment thickness, often derived from top-of-bedrock and borehole data. These products are developed through a combination of fieldwork, subsurface modeling, and the collection and analysis of rock and sediment cores.

To support Coalition goals, member States collaborate with scientists working on related STATEMAP, EDMAP, and FEDMAP projects. Coalition scientists also engage with Tribal Nations in the Great Lakes region to ensure that Tribal interests pertaining to Coalition work are addressed. Through this collaboration, the Coalition unites the efforts of State, Federal, and Tribal Nation stakeholders to advance geologic data production and enhance understanding of the geologic resources of the Great Lakes region.

Assessment of undiscovered conventional oil and gas resources in the Berkine, Illizi, Hamra, Murzuq, and Erdis Kufra Basins of northern Africa, 2026

Released May 20, 2026 11:55 EST

2026, Fact Sheet 2026-3009

Christopher J. Schenk, Tracey J. Mercier, Cheryl A. Woodall, Thomas M. Finn, Heidi M. Leathers-Miller, Michael E. Brownfield, Phuong A. Le, Stephanie B. Gaswirth, Marilyn E. Tennyson

Using a geology-based assessment methodology, the U.S. Geological Survey estimated undiscovered, technically recoverable mean conventional resources of 419 million barrels of oil and 5.5 trillion cubic feet of gas in the Berkine, Illizi, Hamra, Murzuq, and Erdis Kufra Basins of northern Africa.

Factors affecting benthic macroinvertebrate health in the City of Roanoke, Virginia, 2020–2023

Released May 20, 2026 09:14 EST

2026, Report

Samuel Adam Miller, Marcus F Aguilar, Logan Helsley, Sally Entrekin

Major waterways in the City of Roanoke (City) have failed to meet Virginia’s aquatic life designated use since 1996. Segments of the upper Roanoke River lack healthy benthic macroinvertebrate communities which prompted a total maximum daily load (TMDL) study by the Virginia Department of Environmental Quality (VDEQ) to identify the most probable stressor(s) causing the impairment. Excess fine sediment was identified as the most probable stressor impairing benthic macroinvertebrates on portions of the Roanoke River in 2006, and a watershed implementation plan published in 2016 required communities within the impaired watershed to implement projects that would reduce the load of fine sediment entering the Roanoke River. Additional benthic macroinvertebrate sampling and stream habitat assessments along the Roanoke River and Tinker Creek (a tributary to the Roanoke River that flows through the City) revealed continued impaired conditions, and subsequent stressor identification analysis was completed in 2023. Samples collected downstream of the City on the Roanoke River and Tinker Creek generally showed more impaired conditions relative to samples collected at locations upstream of the City. Based on this evaluation, sediment and sediment-bound polychlorinated biphenyls (PCBs) were identified as probable stressors while specific conductance, total nitrogen, and sediment metals were possible stressors in Tinker Creek; however, only a sediment TMDL target was identified to address impaired benthic macroinvertebrate communities. In the Roanoke River upstream of the Niagara Dam, sediment and total phosphorus were identified as probable stressors, sediment polycyclic aromatic hydrocarbons and sediment PCB were considered possible stressors; however, the TMDL target was only for total phosphorus.  

The City partnered with the U.S. Geological Survey (USGS) in 2016 to continuously monitor water quality and streamflow conditions on a major tributary of Tinker Creek, Lick Run, and by 2020, four similar monitoring stations were installed on the Roanoke River and Tinker Creek near the locations of benthic macroinvertebrate sampling. Monitored parameters included streamflow and/or gage height (water level), water temperature, pH, dissolved oxygen, specific conductance, and turbidity. Turbidity is a measure of the relative clarity of the water and was previously used to model suspended-sediment concentrations at the monitoring stations. The City also contracted Kirk Environmental, LLP (KE) to collect benthic macroinvertebrate samples and stream habitat assessments near the locations of the water-quality monitoring stations. Identified benthic macroinvertebrates were used to calculate the Virginia Stream Condition Index (SCI), a multi-metric index composed of eight biological attributes that represent elements of the structure and function of the benthic macroinvertebrate community that measure diversity, composition, and tolerance to pollution.  

Study objective: In this report, benthic macroinvertebrate samples and stream habitat assessment scores collected at four locations on the Roanoke River and Tinker Creek by KE and the VDEQ between 2020 and 2023 were compared to measured water-quality and streamflow conditions prior to sampling to evaluate patterns between benthic macroinvertebrate health, water quality, and hydrology.  

Methodology for construction of a three-layer geologic model of the conterminous United States using land surface, top of bedrock, and top of basement

Released May 19, 2026 11:25 EST

2026, Data Report 1220

Donald S. Sweetkind

This report describes the methodology used for the construction of a digital three-layer geologic model of the conterminous United States by mapping the altitude of three surfaces: land surface, the top of bedrock, and the top of basement. These surfaces are mapped through the compilation and synthesis of published stratigraphic horizons from numerous topical studies. The mapped surfaces create a three-layer geologic model with three geomaterial-based subdivisions: unconsolidated to weakly consolidated sediment; layered consolidated rock strata that constitute bedrock; and crystalline rocks that are described as “basement,” consisting of either igneous, metamorphic, or highly deformed rocks. The data compilation and synthesis are highly dependent on the definition of the informal terms “bedrock” and “basement,” which may describe different ages or types of rock in different parts of the conterminous United States. This report presents the conceptualization of the three mapped layers, describes the datasets used, and summarizes the decisions made while compiling the three-layer model from the various sources. This digital dataset was created as part of efforts by the U.S. Geological Survey to develop subsurface geologic data in geospatial form as part of a broad directive to develop two-dimensional and three-dimensional geologic information at detailed, national, and continental scales. This digital dataset partly fulfills the goal of the U.S. Geological Survey’s National Cooperative Geologic Mapping Program to construct a national-scale three-dimensional geologic model.

Status and understanding of groundwater quality in the San Joaquin Valley Kern County subbasin domestic-supply aquifer study unit, 2022—California GAMA Priority Basin Project

Released May 19, 2026 10:38 EST

2026, Scientific Investigations Report 2026-5012

Jennifer S. Harkness, Kirsten E. Faulkner, Bryant C. Jurgens

The quality of water accessed by domestic wells (here referred to as domestic groundwater resources) in the San Joaquin Valley Kern County subbasin (basin number 5-022.14) was assessed as part of the California Groundwater Ambient Monitoring and Assessment (GAMA) Program Priority Basin Project (GAMA-PBP), in cooperation with the California State Water Resources Control Board. Kern County is at the southern end of the San Joaquin Valley in California, and about 30,000 residents are estimated to use privately owned domestic wells for drinking water. Domestic wells typically draw from shallower parts of the aquifer system than public-supply wells and can be more vulnerable to effects from surface activities. Kern County is host to a highly productive agricultural industry, with Bakersfield as the main urban center. The Kern River runs through Bakersfield from the southern Sierra Nevada and intersects the Kern Water Bank, one of the largest groundwater banking operations in California, at the Kern River Intertie. The section of the Kern River running through the Kern Water Bank is dry most years. Kern County also encompasses some of the most productive oil and gas basins in California, with extensive underground and surface disposal of oil-field wastewater.

This study was based on data collected from 33 sites sampled by the U.S. Geological Survey for the GAMA-PBP in 2022. To provide context for the water quality assessment, measured concentrations were compared to regulatory and non-regulatory health-based and aesthetic benchmarks. A grid-based method was used to estimate the proportions of the groundwater resources used for domestic-supply wells that have water-quality constituents below (low relative concentration), approaching (moderate relative concentration), or above (high relative concentration) benchmark concentrations. At least one measured constituent with a regulatory benchmark was categorized as having a high relative concentration in 72 percent of the aquifer area used for domestic groundwater resources. Inorganic constituents were detected at high concentrations in 45 percent of the domestic groundwater resources, and the constituents detected above regulatory benchmarks were arsenic, nitrate, and uranium. At least one organic constituent was detected at high concentrations in 41 percent of the domestic groundwater resources, and the constituents exceeding regulatory benchmarks were the fumigants 1,2,3-trichloropropane (1,2,3-TCP), 1,2-dibromo-3-chloropropane (dibromochloropropane [DBCP]), 1,2-dibromoethane (EDB), and the per-and polyfluoroalkyl substance (PFAS) perfluorooctanesulfonate. The disinfection by-product chloroform, the fumigant 1,2-dichloropropane, the herbicides atrazine and hexazinone, and the herbicide degradates 2-chloro-6-ethylamino-4-amino-s-triazine, 2-chloro-4,6-diamino-s-triazine, 4-hydroxychlorothalonil, and metolachlor sulfonic acid were detected in more than 10 percent of domestic groundwater resources, but concentrations did not exceed regulatory benchmarks.

Land use, groundwater age (fraction of modern water and mean age), and geochemical environment (oxic or anoxic conditions, pH, alkalinity) were associated with the distribution of high relative concentrations of inorganic and organic constituents. Young, oxygenated water is recharged along the Kern River and adjacent recharge ponds, or as irrigation water in the agricultural areas. High concentrations of nitrate and volatile organic compounds occurred in the oxic water in urban and agricultural areas. The fumigants 1,2,3-TCP, DBCP, and EDB were reported throughout the agricultural areas, whereas chloroform, tetrachloroethene, and PFAS were associated with urban land use. High uranium concentrations were associated with young, modern groundwater in agricultural areas with low pH and high bicarbonate. Total dissolved solids increased with distance from the Kern River, as the contributions of fresh, oxic water decreased. High concentrations of arsenic were present in older anoxic or alkaline groundwater away from areas of recharge. Overall, groundwater age, redox conditions, and the source of recharge as a result of different land uses contribute to large aquifer-scale portions of domestic groundwater resources that exceed health-based benchmarks for nitrate, uranium, and fumigant concentrations.

Precipitation-based flood-inundation maps for the East Fork Little Blue River and tributaries at Lee’s Summit, Missouri, 2024

Released May 19, 2026 09:33 EST

2026, Scientific Investigations Report 2026-5017

Allison A. Atkinson

The U.S. Geological Survey, in cooperation with the City of Lee’s Summit, Missouri, assessed flooding of the East Fork Little Blue River and tributaries for varying precipitation magnitudes and durations, varying antecedent runoff conditions, and projected climate-change conditions. The precipitation scenarios were used to develop a library of flood-inundation maps for a 2.95-mile reach of the East Fork Little Blue River and tributaries within the city.

A two-dimensional U.S. Army Corps of Engineers Hydrologic Engineering Center–River Analysis System (HEC–RAS; ver. 6.5) rain-on-grid model was calibrated to selected runoff events representing a range of antecedent runoff conditions and hydrologic responses. Lowest adjacent grades for structures within the nearby study area were incorporated into the terrain, and depth grids and water-surface elevation grids were developed for the study area. Simulated velocities at selected bridge locations were also developed from the model. The model was calibrated using water-surface elevation data collected from water-level loggers (pressure transducers) and streamflow measurements and water-surface elevation measurements made at a reference point during runoff events. The calibrated HEC–RAS model was used to simulate streamflows from design rainfall events of 15-minute to 24-hour durations and ranging from a 100- to 0.1-percent annual exceedance probability (1-year to 1,000-year recurrence intervals). Flood-inundation maps were produced for depths at a reference location of 3 to 16 feet, or a depth exceeding the 0.1-percent annual exceedance probability interval precipitation. The results of each precipitation duration-frequency value were represented by a 1-foot-increment inundation map based on the generated peak streamflow from that rainfall event and the corresponding water-surface elevation at the East Fork Little Blue River reference location.

Within the HEC–RAS model, 240 scenarios were developed from the design rainfall events with each of 3 antecedent conditions. Additional scenarios were created to simulate the effects of projected precipitation scenarios on the 100-year recurrence interval, 24-hour storm and the 100-year recurrence interval, 6-hour storm. All simulation results were assigned to a flood-inundation map condition based on the generated peak flow and corresponding water-surface elevation at the East Fork Little Blue River reference location.

The flood-inundation maps are shown on a web mapping application made available to the public through the City of Lee’s Summit (hyperlink will be added when available). The flood-inundation maps are tied to real-time precipitation data obtained from the Automated Surface Observing System weather station at the Lee’s Summit Municipal Airport, accessible at https://mesonet.agron.iastate.edu/request/download.phtml?network=MO_ASOS. The availability of these maps, along with information regarding observed rainfall, could help provide emergency management personnel and residents with information that is critical for flood-response activities, such as evacuations and road closures, and for postflood recovery efforts.

Remote sensing enables basin-scale inventories of coal mine methane

Released May 19, 2026 09:28 EST

2026, Environmental Science and Technology (60) 14924-14933

Elise Penn, Daniel J. Jacob, Daniel M. Bon, Kate Howell, Kelly O’Neill, Tia Scarpelli, Zichong Chen, Robert A. Field, C. Ozgen Karacan, Elfie Roy, Daniel Cusworth

Underground coal mines are important global sources of methane, but emission estimates are uncertain. We show that emission estimates for individual mines from aircraft remote-sensing surveys in the United States agree within 40% with direct measurements used for national emission reporting (IPCC Tier 3 estimate). Such direct measurements are unavailable in most countries, which rely on estimated emission factors (EFs) applied to coal-production rates. We find that EFs from IPCC Tier 1 and the Model for Calculating Coal Mine Methane (MC2M) methods overestimate U.S. emissions 3-fold due to incorrect dependence on mine depth. An IPCC Tier 2 method using measured basin-specific mine gas content agrees with direct emission measurements but does not account for gob well emissions and requires gas content data that are generally unavailable. We show that aircraft remote sensing for a small sample of mines can successfully estimate basin-specific EFs for ventilation shafts and gob wells, enabling estimates of basin- and national-scale emissions. We discuss how the method can be applied with satellite remote sensing to quantify coal emissions worldwide.

Whatever it takes— Shaping the L&O Letters Early Career Publication Honor to deliver true benefit

Released May 19, 2026 09:14 EST

2026, Limnology and Oceanography Letters (11)

Rita M. Franco-Santos, Bridget Deemer, Laura J. Falkenberg, Mary R. Gradoville, Scott Hotaling, Erin K. Peck

There are clear advantages for those who openly share their research. Publishing Open Access (OA) articles can increase author visibility (McCabe and Snyder 2014), improve productivity metrics (i.e., more diverse and higher citation rates; Huang et al. 2024; Piwowar et al. 2018), widen collaborative networks (Tai and Robinson 2018), and help secure future funding and/or comply with funder mandates (Herrmannova et al. 2019; Larivière and Sugimoto 2018; McKiernan et al. 2016 and references therein). These benefits can be vital for students and early career researchers (ECRs) trying to advance and thrive in academia. However, publishing papers in Gold OA journals such as L&O Letters comes at a notable financial cost, as these journals require that the corresponding author (or their organization or funder) pay a fee to make their published article immediately freely available to the public. These article processing charges can be prohibitively expensive (Fontúrbel and Vizentin-Bugoni 2021; Mekonnen et al. 2022; Ross-Hellauer et al. 2022). While Read and Publish agreements and waiver programs may be available to help cover these costs, these programs often exclude independent authors as well as those affiliated with ineligible or non-participating institutions (e.g., publisher waivers using the Research4Life eligibility criteria for access currently only allow authors from one of 13 South American countries/territories to publish free of charge). Besides the financial barrier, authors from underrepresented groups can also face a myriad of other publishing roadblocks, such as linguistic challenges (for speakers of English as a foreign language; Amano et al. 2023; Franco-Santos 2024; Ramírez-Castañeda 2020) and geopolitical-scientific bias (e.g., science conducted in the Global South being seen as less impactful and innovative than that conducted in the Global North; Ghosh 2022; Smits et al. 2025). For context, Global South (GS) and Global North (GN) are not geographic determinations (i.e., South and North hemispheres), but geopolitical classifications regarding a nation's level of development (underdeveloped, developing, or developed). For example: Australia and Brazil are both located in the southern hemisphere, but the former is considered as a Global North (developed) country and the latter as a Global South (underdeveloped or developing) country.

When a subset of researchers is unable to openly publish their work, the diversity of voices represented in OA literature can decline (Williams et al. 2023). Loss of diversity is a loss to science, as diversity increases productivity, innovation, and scientific impact (refer to opening quote; Freeman and Huang 2014; AlShebli et al. 2018; Tomillo et al. 2022). To partially address the above-mentioned challenges and enable underfunded ECRs to publish their work in OA format, the biennial L&O Letters Early Career Publication Honor (ECPH) was established in 2020 by the ASLO Raelyn Cole Editorial (RCE) Fellows (Hotaling et al. 2022). Below we reflect on the benefits, outcomes, and scientific impact of the 2022 call and introduce the articles it helped publish in L&O Letters, which are bundled in this section of the ECPH Virtual Issue. Articles published in L&O Letters during other calls are available in their respective sections. We also refer the reader to six articles published by ECRs in other journals (not included in this Virtual Issue) whose content originally warranted their leading author an ECPH in 2022.

Tringa flavipes (Lesser Yellowlegs) from separate breeding sites subdivides the Prairie Pothole Region in space and time during southbound migration

Released May 19, 2026 07:54 EST

2026, Ornithological Applications

Rosalyn E. Bathrick, James A. Johnson, Daniel R. Ruthrauff, Katherine Christie, Anna Courtemanche, Callie Gesmundo, Laura Anne McDuffie, Nathan R. Senner

Some staging regions support multiple groups of the same migratory species, each of which may use the region differently. Characterizing the ways, in which separate groups use such regions can therefore help to identify vulnerabilities during this sensitive period of the annual cycle. The Prairie Pothole Region (PPR) is a massive wetland complex in the northern Great Plains of North America used by ∼11 million shorebirds during migration. The PPR has been heavily modified by agriculture and is experiencing varied effects of global climate change, threatening the health of the shorebirds that rely on it. Here, we used 6 seasons of southbound tracking data of Tringa flavipes (Lesser Yellowlegs)—a long-distance migratory shorebird species with an estimated population decline of 63% over the last 4 decades—from 9 sites across their breeding range to explore differences in migratory behavior within this important staging region. We found that 75% of tracked individuals used the region during southbound migration, and T. flavipes from different breeding sites detoured 110–875 km from their most direct migratory route to access the PPR. Individuals that arrived later stayed longer and made more stops within the region than those that arrived early. Individuals originating from different breeding sites also displayed spatial and temporal segregation within the region: T. flavipes from southwest and central Alaska relied heavily on the northwestern PPR, while those from Canada used the central and southeastern portions of the PPR. Finally, timing of use varied among groups, but the southeastern PPR became increasingly important over the course of the southbound migratory window, as other wetlands likely dried out. Our study highlights the portions of the PPR of critical importance to migrating T. flavipes and the diversity of ways, in which different groups from within the same species can use a single staging region.

ECCOE Landsat quarterly calibration and validation report—Quarter 4, 2025

Released May 18, 2026 11:29 EST

2026, Open-File Report 2026-1014

Md Obaidul Haque, Nahid Hasan, Ashish Shrestha, Rajagopalan Rengarajan, Mark Lubke, Daniel Steinwand, Paul Bresnahan, Jerad L. Shaw, Kathryn Ruslander, Esad Micijevic, Michael J. Choate, Cody Anderson, Jeff Clauson, Kurt Thome, Amit Angal, Raviv Levy, Jeff Miller, Cibele Teixeira Pinto

Executive Summary 

The U.S. Geological Survey Earth Resources Observation and Science Calibration and Validation (Cal/Val) Center of Excellence (ECCOE) focuses on improving the accuracy, precision, calibration, and product quality of remote- sensing data, leveraging years of multiscale optical system geometric and radiometric calibration and characterization experience. The ECCOE Landsat Cal/Val Team continually monitors the geometric and radiometric performance of active Landsat missions and makes calibration adjustments, as needed, to maintain data quality at the highest level.

This report provides observed geometric and radiometric analysis results for Landsats 8 and 9 for quarter 4 (October–December) of 2025. All data used to compile the Cal/Val analysis results presented in this report are freely available from the U.S. Geological Survey EarthExplorer website: https://earthexplorer.usgs.gov.

One specific activity that the ECCOE Landsat Cal/Val Team closely monitored was a Landsat 9 safehold anomaly. On October 17, 2025, Landsat 9 experienced a Solar Array Drive Assembly potentiometer fault. The onboard fault response put both the Operational Land Imager sensor and the Thermal Infrared Sensor into safe mode. Additionally, the Thermal Infrared Sensor focal plane assembly was turned off, but the cryocooler remained on. On October 20, 2025, the Solar Array Drive Assembly recovery commanding was successfully performed to put the spacecraft into nadir viewing mode. The following day, Operational Land Imager activation and recovery started, including focal plane assembly warmup. After reaching nominal operational temperatures and achieving thermal stability, science imaging resumed on October 23, 2025. Additional information about the Landsat 9 safehold anomaly is here: https://www.usgs.gov/landsat-missions/news/landsat-9-returns-normal-operations-following-brief-safehold.

System characterization report on Tanager

Released May 18, 2026 11:04 EST

2026, Open-File Report 2021-1030-W

Minsu Kim, Seonkyung Park, Jeff Clauson, Jim Vrabel, Ajit Sampath

Executive Summary 

This report addresses the system characterization of the Tanager satellite hyperspectral sensor created by Planet Labs PBC. and is part of a series of system characterization reports produced and delivered by the U.S. Geological Survey Earth Resources Observation and Science Cal/Val Center of Excellence. These reports present and detail the methodology and procedures for characterization; present technical and operational information about the Tanager hyperspectral sensor; and provide a summary of test measurements, data retention practices, data analysis results, and conclusions.

This report summarizes the sensor performance of the Tanager based on the U.S. Geological Survey Earth Resources Observation and Science Cal/Val Center of Excellence system characterization process. In summary, we determined that the Tanager exhibits a band-to-band geometric error ranging from -0.074 to 0.097 pixels. Compared to the Landsat Operational Land Imager, geometric offsets ranged from -5.980 meters (-0.20 pixels) to 11.348 meters (0.40 pixels). Radiometric comparisons showed offsets between -0.004 and 0.056 with slopes from 0.830 to 1.066. Spectral shifts are found between 0.65 and 0.75 nanometers. Finally, spatial performance evaluation yielded a PSF full width at half maximum of 1.27 to 1.75 pixels, a relative edge response of 0.802 to 0.651, and a modulation transfer function at Nyquist of 0.488 to 0.253.

Bird migration and energetics simulations incorporating oil spill effects

Released May 18, 2026 10:08 EST

2026, Ecological Modelling (519)

Benjamin M West, Mark L. Wildhaber, Wayne E. Thogmartin, Michael J. Hooper

Oil spills are well-known for causing acute mortality of birds, but sublethal and delayed impacts are less understood. Focusing on the mallard (Anas platyrhynchos), we used simulation modeling to explore how sublethal oiling may affect avian survival and breeding ground body condition. We used empirically informed migration and energetics simulations to model hypothetical spills occurring in northern Arkansas, USA occurring in either January to simulate thermoregulatory stress or March to simulate pre-migration effects. We modeled trace and lightly oiled female mallards (≤5% or 6 to 20% of feather area oiled, respectively), incorporating oiling-induced energetic effects on thermoregulation, flight, and energetic gain. We found that mortality was generally higher for simulated spills occurring in January versus March. In the simulations, mallards lost body mass due to oiling, but surviving individuals could partially recover body mass before arriving at the breeding grounds. Including oiling-induced energetic gain effects in simulations increased mortality as well as increased overall variability of simulation results. This modeling effort identified an important gap in knowledge regarding oiled bird energetics, specifically a need to better quantify oiling-induced energetic gain changes. Although the model is currently limited to a specific species and geographic area, it serves as a proof-of-concept for future research and modeling efforts aimed at understanding more broadly the impacts of oil spills on avian populations.

Simulating past and future refugia for temperate trees in northern Italy

Released May 17, 2026 09:07 EST

2026, Ecography

Azzurra Pistone, Paul Henne, Petra Boltshauser-Kaltenrieder, Willy Tinner, Christoph Schworer

During the Quaternary, trees responded to the climatic changes of glacial–interglacial cycles with large-scale range shifts. Over cold glacials, temperate tree species contracted their ranges and survived in areas known as refugia. Several studies point to the Euganean Hills (Colli Euganei), in Veneto, northern Italy, as one of the northernmost European refugia of temperate tree species during the Last Glacial Maximum (LGM, ca 23 000–19 000 calibrated years BP). Using LandClim, a spatially explicit, dynamic forest landscape model, we demonstrate that climate conditions during the LGM likely allowed temperate tree species to persist in the Euganean Hills. The identified refugial locations lie at intermediate to high elevations and in sheltered valleys within the hilly complex. Therefore, the combined palaeoecological and modelling evidence suggests that today's temperate forests of the Euganean Hills have a full glacial legacy.

Simulations under future climate conditions suggest a collapse of the sub-mediterranean and oro-mediterranean deciduous forests that are prevalent today and the expansion of thermo-mediterranean evergreen forests (with e.g. Quercus ilex, Q. suber, Olea europaea and Pinus sp.). Specifically, the extrazonal population of oro-mediterranean Fagus sylvatica, which is unique to the Po Plain and likely persisted locally through several glacial–interglacial cycles, is predicted to sharply decline and face local extinction, underscoring a conservation hazard.

Distribution, abundance, breeding activities, and restoration efforts for the Southwestern Willow Flycatcher at Marine Corps Base Camp Pendleton, California—2025 Annual Report

Released May 15, 2026 12:38 EST

2026, Open-File Report 2026-1016

Suellen Lynn, Scarlett L. Howell, Barbara E. Kus

Executive Summary 

The purpose of this report is to provide the Marine Corps with an annual summary of the distribution, abundance, and breeding activity of the endangered Southwestern Willow Flycatcher (Empidonax traillii extimus; flycatcher) and to present results of management actions implemented to attract flycatchers and enhance flycatcher habitat at Marine Corps Base Camp Pendleton (MCBCP, or Base). Surveys for the flycatcher were done on Base between May 6 and July 23, 2025. All MCBCP’s historically occupied riparian habitat (core survey area) was surveyed for flycatchers in 2025. None of the non-core survey areas were surveyed in 2025.

No resident flycatchers were detected on Base in 2025. The one resident (female) present in 2024 did not return to the territory she occupied in 2024, and she was not detected within the historically occupied habitat surveyed in 2025.

Eight transient Willow Flycatchers of unknown subspecies were observed on two of the five drainages surveyed in 2025: Las Flores Creek and the Santa Margarita River. No Willow Flycatchers were detected at Fallbrook, Pilgrim, or San Mateo Creeks. Transients in 2025 occurred in mixed willow and riparian scrub habitats, dominated by multiple willow species (Salix spp.). Exotic vegetation was recorded in most flycatcher locations and was dominant (cover of exotics greater than 50 percent) in more than half of all transient locations. The most common exotic plant in habitat used by flycatchers was poison hemlock (Conium maculatum). All six of the flycatchers that were observed closely enough to determine banding status were unbanded.

Two measures were initiated in recent years to attract and retain resident breeding flycatchers on MCBCP: conspecific attraction using flycatcher song broadcasts and installation of artificial seeps to enhance flycatcher habitat. We surveyed plots with and without speakers that broadcast flycatcher vocalizations throughout the breeding season and detected two transient Willow Flycatchers within 20 meters of one speaker in 2025. We set up permanent vegetation sampling points surrounding artificial seeps and nearby sites without artificial seeps (Reference sites) to determine the effects of surface-water enhancement by seep pumps. Vegetation cover was highest near the ground and decreased with increasing height. Woody vegetation made up most of the cover at all height categories. Soil saturation in 2025 was higher at the sites near seeps than at the Reference sites and was associated with higher native herbaceous cover and lower non-native cover.

Continuous stream discharge, salinity, and associated data collected in the lower St. Johns River and its tributaries, Florida, 2023

Released May 15, 2026 10:57 EST

2026, Open-File Report 2026-1012

Jennifer N. Carson, Matthew T. Benacquisto

The U.S. Army Corps of Engineers, Jacksonville District, deepened the St. Johns River channel in Jacksonville, Florida, to accommodate larger, fully loaded cargo vessels. The U.S. Geological Survey (USGS), in cooperation with the U.S. Army Corps of Engineers, monitored stage, discharge, and (or) water temperature and salinity at 26 continuous data collection sites in the St. Johns River and its tributaries. 

This report contains information collected during the 2023 water year, from October 2022 to September 2023. Data at each site were compared for the length of the project, 8 years so far, and on a yearly basis to show the annual variability of discharge and salinity.

The countywide annual rainfall for the 2023 water year was below the average yearly rainfall in four of the five counties. Annual mean discharge at 9 of the 10 tributary monitoring sites was lower for the 2023 water year than for the 2022 water year, and the annual mean flow at Broward River below Biscayne Boulevard near Jacksonville, Florida (USGS site number 02246751), was the lowest recorded at that site for the 8 years of data collection. The annual mean discharge for each of the main-stem sites was higher for the 2023 water year than for the 2022 water year and was above the average for the 8 years of data collected so far.

Among the tributary sites, annual mean salinity was highest at Clapboard Creek above Buckhorn Bluff near Jacksonville, Fla. (USGS site number 302657081312400), the site closest to the Atlantic Ocean, and was lowest at Durbin Creek near Fruit Cove, Fla. (USGS site number 022462002), the site farthest from the ocean, for all years. Annual mean salinity data from the main-stem sites indicate that salinity decreased with distance upstream from the ocean, which was expected. Annual mean salinity for the 2023 water year was higher than or equal to that of the 2022 water year for all main-stem and tributary sites, except at St. Johns River at Dancy Point near Spuds, Fla. (USGS site number 294213081345300), which was lower. Three main-stem monitoring stations (USGS site numbers 295856081372301, 02245340, and 301057081414800) and six tributary monitoring stations (USGS site numbers 300803081354500, 022462002, 301204081434900, 02246459, 02246518, and 02246804) either had the highest annual mean salinities or tied with the highest annual mean salinities at their respective sites since data collection began.

Multi-proxy thermal history of basin heating during Cordilleran orogenesis in the Magallanes-Austral retroarc foreland basin, Patagonian Andes

Released May 15, 2026 09:56 EST

2026, Basin Research (38)

Rebecca A. VanderLeest, Julie C. Fosdick, Theresa Maude Schwartz, E.G. Hyland, M. Mastalerz

Resolving thermal histories in sedimentary basins is crucial for interpreting orogenic growth, basin burial, and tectonic processes during Cordilleran orogenesis. In the Magallanes–Austral Basin, Patagonian Andes, we integrate new (U-Th)/He thermochronology, vitrinite reflectance (%Ro), calcite-cement clumped isotope data and thermal history modelling to resolve the origin of the regionally extensive Paleogene unconformity (51°S–50°S). Thermal history modelling results require post-depositional heating of Palaeocene (Danian–Selandian) strata below the unconformity and suggest maximum burial temperatures of 87°C–101°C (55–52 Ma) and 89°C–92°C (18–16 Ma). For lower Eocene strata above the unconformity, Miocene burial temperatures (89°C–92°C) are consistent with calcite cement formation temperatures (~62°C–92°C) from carbonate clumped isotopes. Our results indicate that basin burial and heating between ca. 60 and 52 Ma were likely driven by shallowing of the subducting Farallon plate and enhanced plate coupling preceding arrival of the Farallon–Phoenix mid-ocean ridge. Subsequent basin inversion and cooling from ca. 52 to 44 Ma correspond with subduction of this mid-ocean ridge. Refined thermal models, constrained by expanded thermochronometric and organic maturation datasets, indicate that up to ~1.7–2.0 km of proximal foreland basin strata were removed during uplift and erosion across the Paleogene basin margin. A return to basin subsidence beginning ca. 44 Ma may reflect dynamic subsidence after passage of the mid-ocean ridge and renewed coupling between the fold-thrust belt and foreland basin system. Neogene thermal histories document continued subsidence, localized hot orogenic fluid flow along stratigraphic boundaries, followed by a final phase of basin inversion and cooling at ca. 18–16 Ma, which we attribute to regional uplift associated with Chile ridge subduction. Altogether, this study demonstrates that multiple thermal indices when analysed and modelled can provide clarity for tectonic and stratigraphic events that affect foreland basins.

Spawning habitat suitability models for Lake Erie cisco (Coregonus artedi) during the historical period of pre- and post-population declines 1877–1957

Released May 15, 2026 09:41 EST

2026, Journal of Great Lakes Research

Katelyn King, Cory Brant, Arthur Cooper, Gust Annis, Matthew Herbert, Karen M Alofs

Coregonine fishes play a key role in the food webs and fisheries of the Laurentian Great Lakes and are a major focus of basin-wide conservation efforts. In Lake Erie, management goals prioritize rebuilding spawning populations of cisco (Coregonus artedi). However, the historical distribution of cisco spawning habitat and the environmental conditions that influence early life-stage success remain poorly defined. We used a novel database of historical coregonine spawning observations as well as novel habitat variables to describe historical conditions to model and determine where and what habitat was historically most suitable for spawning cisco in Lake Erie. The environmental predictors that produced the best model included reefs, distance to rivers, historical substrate, coefficient of variation of ice duration, fetch, and circulation. The highest suitability occurred in areas of high reef probability, near river mouths, in rocky and sandy substrate, and in areas of low variability in historical ice, fetch, and circulation. Suitable spawning habitat is predicted mostly around reefs in the western basin as well as along the coast and near rivers lake-wide. Our model identifies important habitat features and allows managers to envision relevant scales and locations at which to focus restoration efforts.

Melanoma and other melanistic lesions in brown bullhead Ameiurus nebulosus from waterbodies in the northeastern United States and Canada: Identification of risk factors

Released May 15, 2026 09:19 EST

2026, Journal of Fish Diseases

Vicki S. Blazer, P. Emerson, M. Bodnar, Thomas Jones, D. R. Russel, M. Pehrson, Cheyenne R. Smith, Danielle M. Cleveland, Mark J. Henderson, Patricia Mazik

Melanistic lesions, including non-raised black areas due to proliferations of melanocytes and melanomacrophages in the dermis and epidermis, as well as raised black areas consistent with melanoma, are described in brown bullhead (BBH) Ameiurus nebulosus from three water bodies in the northeastern United States and Quebec, Canada. First observed in the Vermont portion of Lake Memphremagog, Vermont, USA and Quebec, Canada, the prevalence of melanistic lesions during 2014–2020 was greater than 30% in BBH 200 mm and longer. In 2023, seven sites throughout the lake were assessed, and prevalence ranged from 18% to 42%. In Hermon Pond, Maine, the prevalence was 29% in 2024, and in Village Pond, New Hampshire, lesions occurred in 22% of BBH in 2025. Compared to skin from visibly normal BBH, skin with melanistic lesions had significantly higher concentrations of seven metals, including arsenic, a known carcinogen and zinc. Lesions associated with oxidative damage, such as the accumulation of ceroid/lipofuscin, were also observed in the gill, spleen and kidney tissue of both affected and visibly normal BBH. The progression of lesions, observed by histopathology, ranged from inflammation, signs of oxidative damage, proliferation and necrosis of club cells, and the presence of melanomacrophages and melanocytes in the epidermis to invasive melanoma and suggests chronic exposure of BBH to environmental initiators and promoters of carcinogenesis.

Feathers and flu: Identifying data gaps in avian influenza host dynamics to prioritize wildlife conservation

Released May 15, 2026 09:14 EST

2026, Wildlife Monographs

Johanna Harvey, Matthew Gonnerman, Shenglai Yin, Cody M. Kent, Joshua Cullen, Jeffery D. Sullivan, Jonathan Dain, Nichola J. Hill, Diann Prosser, Jennifer Mullinax

Highly pathogenic avian influenza viruses (HPAIV) have had disastrous, worldwide effects on wild birds and domestic poultry since the emergence of the A/goose/Guangdong/1/1996 (Gs/GD/96) lineage. The currently circulating H5N1 clade 2.3.4.4b has an expanded set of susceptible hosts, including many migratory wild birds, and is associated with higher transmission rates, increased susceptibility among wild bird hosts, and a greater number of wildlife reservoirs. Certain wild bird life-history strategies and behaviors have been suggested to explain avian hosts’ susceptibility and exposure to HPAIV. These biological traits include gregariousness, such as colonial nesting and mixed flock foraging, predation or scavenging on wild birds, and association with aquatic habitats. Variation in host infection responses (e.g., infectability, shedding rates and duration, mortality rate, antibody development) informs the overall infection risk across avian species, yet the specific role of biological traits is often inconsistent and unclear across taxa. Moreover, the interactions and potential compounding effects among these biological traits remain largely unknown. To develop a more holistic understanding of cumulative risk across bird species, we integrate existing information on infection risk factors (i.e., susceptibility, immunological response, and behavioral traits) into a qualitative multivariate analysis. This approach enabled us to examine how infection risk factors relate to biological traits (e.g., phylogeny, physiology, behavior, species range) and to begin disentangling their complex interactions. We quantified and summarized these risk factors across host species and qualitatively ranked species by their viral responses along a proposed HPAIV response continuum, guided by expectations of traits and metrics associated with competence or vulnerability to HPAIV. In doing so, we aimed to better understand how viral responses and biological traits synergistically interact to influence cumulative risk across wild bird species. This work broadly expands on the previous avian influenza literature, which has focused on Anseriformes and Charadriiformes as primary viral reservoirs. We tie our findings to effective disease management responses with links to risk components, including descriptions of potential surveillance strategies applied to research and One Health goals, as well as a fuller understanding of how resources may be better deployed for rapid response when spillovers do inevitably occur. Additionally, we identified numerous areas where vital epidemiological information is lacking to best characterize the spread of these viruses. Ultimately, this improved understanding will help identify and inform disease management needs and decision making.

Baseflow and snowmelt sustained streamflow in the Upper Colorado River Basin, 1986-2020

Released May 15, 2026 08:34 EST

2026, Environmental Research: Water (2)

Olivia L. Miller, Matthew P. Miller, Patrick C. Longley, Noah M. Schmadel, Daniel R. Wise, Morgan C. McDonnell, Jay R. Alder

The Upper Colorado River Basin (UCRB) faces substantial water availability limitations. Although most streamflow originates as snowmelt, the partitioning of snowmelt between surface runoff and groundwater recharge and subsequent groundwater discharge to streams is highly uncertain. On average, over half of the streamflow in the UCRB is estimated to originate from groundwater discharge to streams, highlighting the importance of baseflow in sustaining surface water. However, the historical patterns of baseflow and streamflow, along with their variability over space and time and their specific sources, remain unknown at the basin scale. This study addresses those gaps by characterizing the sources and transport pathways of both baseflow and streamflow in the UCRB at a seasonal timestep from 1986 to 2020, including the lagged delivery of subsurface water to streams beyond the current season, using coupled models of baseflow and streamflow. Between 1986 and 2020, on average 63% of UCRB streamflow originated from baseflow. About half of this baseflow took longer than one season to reach streams, and outside the snowmelt season, baseflow was the dominant source of streamflow. Snowmelt was a key source of both baseflow and streamflow. Current season snowmelt contributed 33% of streamflow via runoff, and 22% of the 29% of streamflow that originated as current season baseflow via subsurface flow to streams. Over the study period, baseflow index (BFI) declined in headwaters and increased at mid-elevations. Springtime increases in BFI demonstrate the increasingly important role baseflow plays in water supply. Identifying the sources, locations, and timing of water that contributed to the UCRB outlet can inform management of water resources in the basin.

Quantifying southern sea otter (Enhydra lutris nereis) reactions to a quadcopter drone in central California

Released May 15, 2026 08:25 EST

2026, Marine Mammal Science (42)

Colleen Young, Julie L. Yee, Gena Bentall, Michelle M. Staedler, Lilian P. Carswell, Margaret Daly

Drones are useful for wildlife research and management, but they can cause disturbance and harassment to wildlife. Sea otters (Enhydra lutris) are candidates for drone-based observation and monitoring but are vulnerable to disturbance. No studies have evaluated drone effects on sea otter behavior, but based on prior disturbance studies, we hypothesized: (1) sea otters would exhibit behaviors indicating higher reactivity in the presence of drones than when drones were absent and (2) drone disturbance to sea otters would be greater when drones were closer. At two sites in Monterey Bay, CA, we conducted 37 observational sessions, recording behavior codes for focal sea otters during a baseline (no drone) period and three consecutive drone flights. Data were analyzed using ANOVA and ordinal logistic regression models. At both locations, focal sea otters had higher behavior codes during drone trials compared to baseline, and behavior codes increased with descending drone altitude. Pup presence, group size, flight trial number, and gull presence were significant covariables. We calculated multipliers to predict drone-mediated behavioral responses at a range of drone altitudes. Our findings can inform best practices for a variety of uses of drones around sea otters, including population monitoring, oil spill response, and drone photography/videography.

VegET evapotranspiration for Africa: Continental-scale simulation, multi-product evaluation, and drought assessment

Released May 15, 2026 08:10 EST

2026, Journal of Hydrology: Regional Studies (66)

Komlavi Akpoti, Naga Manohar Velpuri, Mansoor Leh, Stefanie Kagone, Kirubel Mekonnen, Afua Owusu, Mulugeta Tadesse, Paranamana T. Prabhath, Lahiru Madushanka, Tharindu Perera, Gabriel Edwin Lee Parrish, Vinay Nangia, Souleymane Sy, Jan Bliefernicht, Samuel Guug, Abdulkarim Seid, Gabriel B. Senay

Gabriel B. Senay, editor(s)

Study region

Continental Africa, encompassing diverse climatic zones—tropical, arid, and temperate—and spanning major transboundary river basins such as the Nile, Niger, Congo, Volta, and Zambezi River Basins. The region exhibits pronounced hydroclimatic gradients and heterogeneous land use systems ranging from rainfed croplands and rangelands to dense tropical forests and irrigated schemes.

Study focus

Actual evapotranspiration (ETa) is a central component of the terrestrial water balance, governing the redistribution of water and energy between the land surface and the atmosphere. Accurate estimation of ETa at continental scale is critical for hydrological monitoring, water resource management, and climate adaptation, as well as for quantifying water, energy, and carbon fluxes that underpin sustainable development. In this study, we applied the agro-hydrologic VegET v2 model to simulate a new, high-resolution, continental-scale ETa dataset for Africa (2000–2021). The model results were benchmarked against four widely used remote sensing-based products—MODIS16 v6.1, SSEBop v6.1, WaPOR v3, and GLEAM v4.1a—across major climate zones, land use types, and River Basins, providing a comprehensive multi-product evaluation of evapotranspiration dynamics across the continent.

New hydrological insights for the region

Validation against eddy covariance flux tower observations at eight representative sites confirmed that VegET v2 accurately reproduces the seasonal dynamics of observed ETa, achieving a correlation (r) of 0.8 and an RMSE of 25 mm month⁻¹ —accuracy that is comparable to or higher than accuracies of satellite-based products MODIS16, SSEBop, and GLEAM. This study represents one of the first Africa-wide hydrological simulations of ETa, extending the VegET model beyond basin-scale applications. Intercomparisons reveal that VegET aligns closely with MODIS16, SSEBop, and GLEAM in humid and tropical regions (r = 0.80–0.90; RMSE < 20 mm month⁻¹), while greater discrepancies appear in arid and semi-arid zones, where WaPOR tends to overestimate ETa (RMSE ≥ 28 mm month⁻¹). Despite these differences, VegET effectively captures spatial and temporal ETa variability across rainfed croplands, forests, and savannas, supporting its utility in regional water balance assessments, water accounting, and drought monitoring. A key application of VegET v2 is the Evapotranspiration Deficit Index (ETDI), derived by integrating VegET-based ETa with potential evapotranspiration (PET) to quantify water stress. ETDI successfully captured major drought episodes across Africa, including persistent Sahelian and southern African dry spells, the 2020–2021 winter drought in the Maghreb, and the 2018–2019 austral summer drought in southern Africa, while identifying positive anomalies over central Africa indicative of recurrent wetness. These results underscore VegET’s capability as a hydrologically consistent, operational tool for continental ETa monitoring and drought assessment, offering support for basin-scale water balance studies, food security planning, and climate resilience across Africa’s diverse hydrological environments.

Evaluation of stream capture related to groundwater pumping, middle Humboldt River Basin, Nevada

Released May 14, 2026 10:40 EST

2026, Professional Paper 1906

Kyle W. Davis, William G. Eldridge, Kip K. Allander, David E. Prudic, Murphy A. Gardner, Michael T. Pavelko, Cara A. Nadler

Historical, future, and potential stream capture from groundwater pumping in the middle Humboldt River Basin (MHRB), Nevada, is estimated using a calibrated numerical groundwater flow model. The model was developed to estimate (1) stream capture, which is the change in flux between the groundwater system and the Humboldt River and tributaries, and (2) change in streamflow, which is the change in streamflow estimated for the Imlay gage on the Humboldt River (U.S. Geological Survey streamgage 10333000). Historical stream capture for water years (WYs) 1961–2015 is estimated using recorded and estimated groundwater pumping during that period. Future (predictive) stream capture was based on historical stresses (WYs 1961–2015) using a scenario that simulated non-mine pumping from WY 2015 at a uniform rate for 100 years into the future. Potential stream capture throughout the middle Humboldt River Basin from groundwater pumping during varying durations of time are presented in a series of capture maps. Maps also are presented that show the potential to capture from groundwater evapotranspiration, as well as the storage changes for pumping duration of 100 years.

Estimates of historical stream capture from the mainstem Humboldt River during the early 1960s are less than 400 acre-feet per year (acre-ft/yr) when groundwater withdrawals and pumping rates were relatively small compared to more recent times. In the late 1980s and early 1990s, groundwater withdrawals increased and estimated historical stream capture also increased from about 4,000 acre-ft/yr in the late 1980s and early 1990s to as much as 18,800 acre-feet (acre-ft) in WY 1998. In WY 2015, estimated historical stream capture declined to about 13,000 acre-ft because of decreasing groundwater withdrawals and lower streamflow during the drought of WYs 2012–15, resulting in less stream water available for capture. Stream capture was estimated for 100 years into the future based on WY 2015 non-mine pumping rates and mine-dewatering activity through WY 2015. Stream capture is forecast to increase to about 23,000 acre-ft/yr, and streamflow in the Humboldt River could decrease by as much as 19,000 acre-ft/yr.

Pumping for mine-dewatering and the associated discharge of that water affects streamflow in the Humboldt River at Imlay, Nevada (U.S. Geological Survey streamgage 10333000). Historically, from WYs 1991 to 2015, streamflow was greater at Imlay gage during active mine-dewatering from mine-water discharge operations and increased by as much as 105,000 acre-ft in WY 1998. The increase was attributed mostly to the discharge of groundwater from mine-related dewatering operations directly into the mainstem Humboldt River or its tributaries, with some of this increase associated with return flows from discharge to rapid infiltration basins. Results indicate that streamflow at Imlay gage is expected to decrease by as much as 1,600 acre-ft/yr 30 years after mine-related pumping and discharge are discontinued. The streamflow reductions at the Imlay gage are expected to then decrease to around 500 acre-ft/yr, 100 years after mine-related pumping and discharge are discontinued.

Potential capture maps were produced for pumping durations of 10, 25, 50, and 100 years. Capture map results indicate that areas of greater potential stream capture occur adjacent to the Humboldt River and for upstream tributaries areas north of the Humboldt River.

Predictable seismic cycles result from structural rupture barriers on oceanic transform faults

Released May 14, 2026 09:22 EST

2026, Science (392) 718-723

Jianhua Gong, Wenyuan Fan, Jeffrey J. McGuire, Mark D. Behn, Jessica M. Warren, Emily Roland, M. S. Boettcher, J. A. Collins, Y. Liu, C. R. German

Earthquakes of magnitude (M) >5.5 on oceanic transform faults (OTFs) repeatedly rupture the same locked patches, sometimes quasiperiodically. These patches are separated by “barriers” that halt earthquake propagation and slip mostly aseismically. However, the physical processes governing this systematic behavior remain unclear. We analyzed two barriers along the Gofar transform fault that have arrested ~15 M6 earthquakes over the past three decades. Ocean bottom seismometer data indicate that the barriers hosted intense microseismicity before the mainshocks and comprise multistrand faults and transtensional stepovers with 100- to 400-m lateral offset. These characteristics contradict earthquake rupture termination models invoking velocity-strengthening friction or large geometric steps and instead point to damage-enhanced porosity and dilatancy-strengthening mechanisms. By isolating rupture segments, the barriers regulate the quasiperiodic recurrence of OTF earthquakes.

Syn-magmatic subsidence during the early stages of continental rifting in the Mesoproterozoic—A reanalysis of legacy data for the Midcontinent Rift, western Lake Superior

Released May 14, 2026 08:37 EST

2026, Geosphere

V. J. Grauch, Laurel G. Woodruff, Samuel J. Heller, Esther K. Stewart

The Midcontinent Rift system (ca. 1.1 Ga) is a 2000-km-long series of elongated volcanic and sedimentary troughs and associated intrusive centers exposed chiefly in the Lake Superior region of North America. The rift system represents a long history of intense magmatism and subsequent sedimentation that was arrested by far-field tectonic events before sea-floor spreading was established. The premature cessation preserved a record of processes related to the beginning of continental rifting.

The rift system under Lake Superior has been long studied using seismic-reflection data collected as part of the Great Lakes International Multidisciplinary Program on Crustal Evolution (GLIMPCE). We reexamine GLIMPCE Line C by developing a detailed velocity model for time to depth conversion constrained by other legacy data. We corroborate the model and develop a geologic interpretation using gravity and magnetic modeling and ties to geology mapped onshore.

We recognize superposed subsiding sedimentary and volcanic basins for the southern half of the Line C depth section. This interpretation differs from previous paradigms that show major crustal faults that bound half-grabens or full grabens. We conclude that high-velocity (6.9 km/s) intrusive zones rather than major crustal faults border the sides of the basins. We speculate that the volcanic basin represents the initiation of seaward dipping reflectors.

The syn-magmatic subsidence can be explained by dike injection and volcanic loading. Discrete lava basins throughout the region likely subsided at different times in a disorganized manner along the rift trend, raising questions about the long-term role of lithospheric thinning and melt generation.

Practical guidance for engaging end-users and experts in developing scientific tools

Released May 13, 2026 11:55 EST

2026, Scientific Investigations Report 2026-5137

Kaylin R. Clements, James J. English, Emily J. Wilkins, Megan A. Moore, Rudy Schuster

This report provides actionable guidance for scientists developing scientific tools that inform on-the-ground decision making. Scientific tools, in the context of this report, are technology or protocols that help practitioners collect and analyze their own data, and information products and web tools that practitioners could use to inform decisions. Engaging end-users and fellow experts is fundamental to the creation of useful scientific tools. Scientists can use clear and specific direction on action steps and activities to effectively engage with end-users and fellow experts during development. Our study explores lessons learned from six U.S. Geological Survey projects that designed and implemented engagement activities with end-users and experts to coproduce scientific tools for natural resource managers. U.S. Geological Survey teams engaged end-users and experts across the United States from Federal, State, and local governments; universities; Tribes; territories; and nongovernmental organizations in designing and developing scientific tools intended to support end-users in their work. An online survey with 98 participants measured satisfaction across several indicators of successful engagement, including engagement activity frequency, sufficient opportunities to provide feedback, feedback implementation, inclusion of necessary perspectives, and functionality of the tool for end-users. Semistructured interviews were held with project leads, during which the project leads reviewed a summary of the survey results. The project leads reflected on the engagement efforts used in their project, then described lessons learned from the engagement experience and participant feedback. Common themes for ensuring effective engagement identified through thematic analysis included engaging end-users during product conceptualization; establishing clear roles and expectations; considering who end-users are and how end-users may use the tool; recruiting participants through your network, boundary spanners, and leadership; understanding individual use cases; communicating how feedback was integrated into the product; and strategically using virtual meeting tools. This guide shares practical steps and exercises for planning and facilitating effective engagement based on lessons learned from project leads and case study summaries of each project.

Timing, uncertainty, and opportunity cost: Lessons for ecosystem modification on the Colorado River

Released May 13, 2026 09:43 EST

2026, Journal of Environmental Economics and Management (139)

Pierce Donovan, Lucas Bair, Matthew N. Reimer, Michael R. Springborn, Charles B. Yackulic

While conservation goals have long been pursued through traditional species-augmenting actions, a broader set of episodic ecosystem modification (EEM) actions, such as hydropower dam releases, prescribed fire, and beach nourishment, is garnering attention. EEM actions face several implementation challenges stemming from high opportunity costs, delayed effect mechanisms, reliance on monitoring for deployment timing, and outcome uncertainty due to infrequent use. In this paper, we study the use of EEM actions in the form of designer flows—ecologically-motivated releases of water into regulated river segments—to maintain a viable population of a threatened native fish species in the Colorado River. We demonstrate how the cost-effectiveness of EEM actions can be hampered by the complex and delayed effects on species viability, but enhanced through targeted monitoring for timing deployment and experimentation for reducing uncertainty about effectiveness.

Integrating mark-recapture, catch, and expert habitat assessments to quantify recent increases in humpback chub abundance over a 200 km long river segment of the Colorado River in western Grand Canyon

Released May 13, 2026 09:20 EST

2026, Canadian Journal of Fisheries and Aquatic Sciences (83) 1-13

Maria C. Dzul, David R. Van Haverbeke, Kirk Young, Charles B. Yackulic, Pilar Rinker, Michael D. Yard

Humpback chub, Gila cypha, were historically distributed throughout large portions of the Colorado River basin and were federally listed in 1967. In the Grand Canyon segment of the Colorado River, located below Glen Canyon Dam, chub abundances continued to decline through the early 2000s. Recently, catch has increased substantially, especially in the western Grand Canyon. Here, we integrate mark-recapture and catch data of subadult and adult humpback chub, with expert assessments of habitat suitability and an underlying model of spatial autocorrelation, to estimate abundance in western Grand Canyon from 2010 to 2024, a time of rapid population increase and expansion. Our model suggests that adult abundance grew ∼160 fold during this 15-year period, with a median adult population abundance of 70 000 (40 000–200 000; 95% credible interval) in 2024. Our approach identifies years with high population growth and indicates that the spatial distribution has changed over time. We test the sensitivity of our results to movement into sampling reaches during sampling with baited hoop nets. Despite rapid population growth, the resilience of humpback chub in western Grand Canyon is unknown.

Storm surge barriers reduce seaward sediment supply to lagoonal estuaries

Released May 13, 2026 08:51 EST

2026, Earth's Future (14)

David K. Ralston, Philip M. Orton, John C. Warner, Shima Kasaei

Numerical simulations with realistic forcing of fixed infrastructure for a proposed storm surge barrier for a lagoonal estuary, Jamaica Bay (New York, USA), are analyzed during typical forcing conditions to assess alterations to flow and sediment transport with the barrier open. Lagoonal estuaries are shallow and have modest watershed freshwater and sediment inputs, so sediment delivery is primarily from offshore by tidal transport. The storm surge barrier infrastructure across the inlet channel reduces cross-sectional area and increases tidal velocities, increasing frictional and form drag. The overall reduction in tidal amplitude is about 1%, but the quarterdiurnal M4 component decreases by 11%. The salinity and stratification in the estuary are only slightly modified by mixing by stronger velocities near the barrier. Sediment transport in the inlet scales approximately with tidal velocity cubed and net landward transport is driven by flood-dominant tidal asymmetry. Additionally, tidal asymmetry in the jet flow through barrier openings causes a divergence in sediment transport within several kilometers. The alterations to the tidal currents reduce sediment import to the bay by 20% for fine sand; transport of sediment with slower settling velocities is less affected, with reductions of 3% for medium silt and <1% for fine silt. The study examined tidal exchange with an open barrier, but the overall impact also depends on barrier operations during major storm events. The impacts of barrier infrastructure on lagoonal estuaries are distinct from other estuary types due to their modest freshwater input, predominance of tidal transport, and offshore sediment supply.

Salinas Valley integrated hydrologic and reservoir operations models, Monterey and San Luis Obispo Counties, California

Released May 12, 2026 10:30 EST

2026, Scientific Investigations Report 2026-5005

Wesley R. Henson, Randy Hanson, Scott Boyce, Joseph Hevesi, Marisa M. Earll, Deidre M. Herbert, Elizabeth R. Jachens

The area surrounding the Salinas Valley groundwater basin in Monterey and San Luis Obispo Counties of California is a highly productive agricultural area, contributes substantially to the local economy, and provides a substantial portion of vegetables and other agricultural commodities to the Nation. This region of California provides about half of the Nation’s lettuce, celery, broccoli, and spinach each year. Thus, this agricultural area provides substantial volumes of agricultural products not just for California but for the United States.

Changes in population and increased agricultural development, which includes a shift toward more water-intensive crops, and climate variability, have put increasing demand on both surface-water and groundwater resources in the valley. This situation has resulted in water management challenges in the Salinas Valley that generally relate to the distribution of the water supply throughout the basin. Where and when the water is present in the surface and subsurface does not coincide with where and when the water is needed. Historically, to deal with the distribution issue, water has been used conjunctively in the valley. Conjunctive use is a water management strategy that coordinates surface-water and groundwater use to maximize water availability. Groundwater is used throughout the Salinas Valley to meet water demands when surface-water supplies are insufficient. The availability of surface water is constrained by climate. Precipitation and streamflow vary seasonally and year to year. Although there are two reservoirs in the Salinas Valley to capture and store water during wet periods, the only conveyance of reservoir water to coastal agricultural areas is the Salinas River. Increasing demand for groundwater and surface-water resources throughout the Salinas Valley has resulted in undesirable effects from unsustainable water use, such as surface-water depletion, groundwater-level declines, storage depletion in the principal aquifers, and seawater intrusion. To address these escalating issues, local communities, water management agencies, and groundwater sustainability agencies are evaluating how to sustainably manage both their surface-water and groundwater resources. To meet water demands and reduce the undesirable effects of unsustainable water use, continued conjunctive management of surface water and groundwater would ideally incorporate strategies to deal with increases in demand and climate variability.

To evaluate the challenging water management issues in the Salinas Valley, the U.S. Geological Survey, Monterey County Water Resources Agency, and the Salinas Valley Basin Groundwater Sustainability Agency developed a comprehensive suite of models that represent the Salinas Valley hydrogeologic system called the Salinas Valley System Model. The geologic framework is known as the Salinas Valley Geologic Framework and was developed to characterize the subsurface using various topographic and geologic data sources, including information on hydrogeologic units, their surfaces and extents, geologic structures, lithology, and elevations from borehole data and cross sections, as well as details on faults and existing models. The surface-water model is called the Salinas Valley Watershed Model and simulates the Salinas River watershed. Monthly surface-water inflows into the integrated hydrologic model domain were simulated using the Salinas Valley Watershed Model. The historical model uses historical climate data, water and land use data, and reservoir releases to simulate agricultural operations, including landscape water demands, diversions, and reclaimed wastewater. The operational model adds an embedded reservoir operations framework to the simulation of the historical model that allows specified operational rules to simulate reservoir releases and changes in reservoir storage. The operational model assumes current reservoir operations and constant land use, which differs from historical conditions. Thus, the operational model is a hypothetical baseline model that can be used by local water managers to evaluate and quantify potential benefits of water supply projects. Together, the geologic framework, watershed, historical, and operational models form a tool that can be used to simulate irrigated agriculture and associated reservoir operations of the integrated hydrologic system of the Salinas Valley.

Toxicity of synergized permethrin residues in cattle dung to two temperate dung beetle species after application of common livestock pour-on treatment

Released May 12, 2026 10:15 EST

2026, Agriculture, Ecosystems, and Environment (408)

Michael C. Cavallaro, Michelle L. Hladik, Rodrigo Soares, Mikaela Anderson, W. Wyatt Hoback

Essential to pasture health, dung beetles (Coleoptera: Scarabaeidae) provide key ecosystem services across natural and managed rangeland habitats. Insecticide residues in livestock dung can negatively impact dung beetle populations, and synergized pyrethroid products are commonly used to combat resistant pest fly populations. Here, permethrin residues were measured by GC-MS/MS in fresh cattle feces on Days −2 (pretreatment), 4, 8, 16, and 30 after the label rate application of a formulated pour-on treatment (a.i. 5% permethrin, 5% piperonyl butoxide [PBO]). Mean (± SE) measured permethrin concentrations were the highest on Day 4 at 1400 ± 360 ng of permethrin/g of dung (dry weight) with a maximum concentration of 2200 ng/g. Approximately, 99% of applied permethrin was excreted by Day 16, with no detection by Day 30. Field-collected dung was used in a 48-hour toxicity test and with three treatment groups (control [Day −2], low risk [Day 16], and high risk [Day 4]). Two temperate dung beetle species were tested: Onthophagus pennsylvanicus Harold and Canthon chalcites Haldeman. Mean (± SE) mortality of O. pennsylvanicus was 28 ± 5% and 58 ± 13% for low and high risk treatments, respectively. Mean (± SE) mortality of C. chalcites was lower than O. pennsylvanicus with 10 ± 4% and 40 ± 10% for low and high risk treatments, respectively. PBO was detected on Days 4 and 8, and the permethrin:PBO ratio was 10:1 on Day 4, i.e., high risk treatment. Data presented highlight episodic risks of pour-on products and support threshold-based, integrated pest management approaches.

Top Elevation of Glacial Till and Thickness of the Big Sioux Aquifer Delineated From Electrical Resistivity Tomography Surveys Near Sioux Falls, South Dakota, 2022 and 2025

Released May 12, 2026 09:48 EST

2026, Scientific Investigations Report 2026-5023

Colton J. Medler, Todd M. Anderson

The City of Sioux Falls, South Dakota, requested the U.S. Geological Survey perform electrical resistivity surveys on three parcels of land north of the city. Electrical resistivity data were collected along a total of 22 transects during March 14–18, 2022, and November 17–21, 2025. Results from electrical resistivity surveys were used to delineate the top of glacial till deposits for the purpose of characterizing the Big Sioux aquifer near the city. Delineating geologic contacts provides important information on groundwater storage, flow dynamics, well design and placement, contaminant transport, groundwater–surface-water interactions, and regional water modeling. The top elevation of glacial till and the thickness of the Big Sioux aquifer varied among the three survey areas. The interpreted top elevation of glacial till in the North survey area decreases from east to west toward a slough, with elevations ranging from 1,403 to 1,418 feet (ft). The estimated thickness of the Big Sioux aquifer in the North survey area increased from east to west, with thicknesses ranging from 23 to 38 ft. The top elevation of glacial till in the Well 72 survey area generally decreases from northwest to southeast. Top elevations of the glacial till in the Well 72 survey area ranged from 1,400 to 1,409 ft along the southern end of transect W72_2. The estimated thickness of the Big Sioux aquifer in the Well 72 survey area was greatest along a southeast to northwest trending channel, with thicknesses ranging from 28 to 40 ft. The top elevation of glacial till in the Nose survey area generally decreases west toward the Big Sioux River. Top elevations of the glacial till in the Nose survey area ranged from 1,362 to 1,395 ft. The estimated thickness of the Big Sioux aquifer in the Nose survey area ranged from 33 to 70 ft.

Calcareous nannofossil assemblage changes in the Surprise Hill core and their implications for floral response to the Paleocene-Eocene Thermal Maximum across the Salisbury Embayment of Virginia, USA

Released May 12, 2026 09:47 EST

2026, Marine Micropaleontology (204)

Masayuki Utsunomiya, Jean Self-Trail, D. Clay Kelly, Xiaodong Zhang, Kristina Frank Gardner, James C. Zachos

We present Paleocene-Eocene calcareous nannofossil biostratigraphy and paleoecology for the Surprise Hill core, U.S. Atlantic Coastal Plain, Virginia. Calcareous nannofossil datums ranging from Zone NP3 to NP14 were identified. The Danian-aged Brightseat Formation rests unconformably atop the Lower Cretaceous Potomac Group at 211.4 m and disconformably underlies the Aquia Formation at 208.8 m. The absence of Zone NP7 suggests a hiatus is present in the Aquia Formation (Zones NP5 – NP9a). The contact between the Marlboro Clay and the overlying Nanjemoy Formation (Zones NP10 – NP14) at 189.5 m is truncated. The Paleocene-Eocene transition is marked by a shift from glauconitic sands of the Aquia Formation to pelitic muds of the Marlboro Clay at 202.7 m. A 3–3.5‰ negative δ13C excursion of benthic foraminifer and a thin dissolution interval (201.6–202.5 m) are recorded in the basal Marlboro Clay. Nannofossil response to the Paleocene-Eocene Thermal Maximum (PETM) include (1) a bloom in taxa with affinities for changing salinity conditions just prior to the PETM basin wide (Hornibrookina australis arca), (2) a decline in taxa with ecological affinities for cool, eutrophic waters (Chiasmolithus bidens) during PETM, (3) fluctuations in mesotrophic to eutrophic, opportunistic taxa (e.g., Neochiastozygus junctus) during PETM, (4) successive turnovers in species of Toweius spp. during core-PETM and its recovery. Our findings suggest that overall nannofossil assemblages in the southernmost portion of the Salisbury Embayment responded similarly to assemblages from South Dover Bridge, but had differing response to local changes in nearshore paleoecology.

Watershed Continuum Monitoring Approach: Combining multiple water quality patterns along stream and river flowpaths to track sources, pathways, and processing of pollutants

Released May 12, 2026 09:25 EST

2026, Ecological Engineering (229)

Sujay Kaushal, Ashley Mon, Stanley Grant, Paul M. Mayer, Aaron J. Porter, Andrew J. Sekellick, Jason Hamilton Chase, Shantanu Bhide, John D. Jastram, Tammy Newcomer-Johnson, Sydney A. Shelton, Alexis M. Yaculak, Joseph T. Malin, Carly Marcella Maas, Nicholas Salanitri, Daniel J. Silberstein, Steven P. Hohman, Ashley B. Dann, Weston M Slaughter, Megan A. Rippy, Ahmed Monofy, Ruth R. Shatkay, Jenna E. Reimer, Madeleine Seppi, Randi Noel, Julianna Mussa, Bennett Kellmayer, Gwendolyn Sivirichi, Melissa Grese, Walter L.M. Boger, Jeffrey G. Chanat, Shuiwang Duan, Kenneth T. Belt

There is a growing need to improve and expand water quality monitoring approaches to more accurately track the sources, fate, and transport of multiple chemicals and pollutants holistically and quantify the effects of best management practices (BMPs) at the watershed scale. An overarching question raised by scientists, environmental managers, and the general public is: how far can water quality impacts from disturbances or benefits from watershed management and restoration propagate along stream and river flowpaths? Many studies using the classic watershed approach focus on analyzing changes in water quality over time at one or a few sampling stations, whereas theories such as the River Continuum Concept focus on predicting shifts in energy sources and biological communities along rivers but have not been directly applied to water quality. We propose to merge these concepts to create a Watershed Continuum Monitoring Approach (WCMA) that combines both spatial and temporal monitoring in order to better detect and quantify trends and transitions in multiple water quality indicators along flowpaths. Specifically, an array of multiple water quality indicators are analyzed at multiple downstream points along a watershed flowpath over time. These multiple water quality indicators are analyzed together for making comparisons to infer hydrological, biological, and geochemical processes controlling sources, transport, and attenuation of pollutants (e.g., analagous to stream tracer studies at the watershed scale). The WCMA leverages the natural expansion of watershed areas along a flowpath, which reflect transitions in land use, land cover, and environmental management across spatial and temporal dimensions for making direct comparisons across different stream reaches and spatial trend analysis. WCMA facilitates monitoring of multiple water quality indicators together, and identifcation of hot spots in sources and attenuation of pollutants or mixtures of pollutants. We illustrate practical applications of the WCMA to analyze water quality trends, transitions, and tradeoffs (i.e., a tradeoff occurs when one pollutant is reduced but another is directly or indirectly increased downstream). We explore case studies that quantify: (1) downstream reductions in concentrations of multiple pollutants along a stream flowing to a major drinking water source due to engineered and nature-based solutions, (2) downstream reductions in multiple pollutants and water quality tradeoffs along streams experiencing stormwater BMPs and stream restoration, (3) comparisons in downstream reductions of multiple pollutants and nutrient uptake along streams draining into major drinking water sources based on types of stream restoration, (4) comparisons of downstream pollutant reductions along streams experiencing riparian forest conservation vs. stream restoration, and (5) mapping and visualizing hot spots of increasing water quality problems such as hypoxia, contaminant mobilization, and freshwater salinization that extend downstream to tidal rivers of the Chesapeake Bay. We explore future applications of WCMA for tracking decreasing trends in salinity, E. coli, and other pollutants of emerging concern. WCMA can holistically inform progress towards achieving multiple water quality goals and also be used as a screening tool for selecting monitoring sites and targeting management in strategic locations. Overall, WCMA enables the simultaneous quantification and comparison of sources and transport and attenuation rates for different chemicals and pollutants across a broader range of watershed sizes and flowpath lengths, which is critical for understanding ecological, hydrological, geochemical, and biogeochemical processes along human-impacted streams and rivers.

Effects of wildfire on soil hydraulic properties in the western Oregon Cascades

Released May 12, 2026 08:22 EST

2026, JGR Biogeosciences (131)

Cedric Pimont, Evan A. Thaler, Brian A. Ebel, Kevin D. Bladon

Wildfires can substantially impact the hydrology of forested watersheds, increasing the risk of hydrologic hazards such as flash floods and debris flows. Soil hydraulic properties related to infiltration are a key control in determining the timing and magnitude of these hydrogeomorphic events. In our study, we collected 445 soil cores from burned (216 cores) and unburned (229 cores) reference catchments and analyzed them for soil hydraulic properties 10 months after the 2022 Cedar Creek Fire in Oregon, USA. We observed significantly greater field-saturated hydraulic conductivity (Kfs), sorptivity (S), and wetting front potential (Ψf) in burned soils relative to unburned soils, with median ratios of 5.7, 4.4, and 5.0, respectively. Among low-, moderate-, and high burn severity groups, soil hydraulic properties were not statistically different. Reductions in median soil bulk density with increasing burn severity suggested an expansion of pore sizes, which may have been partially responsible for increasing Kfs and S. Additionally, in some burned soil samples, the increase in soil hydraulic properties may have been partially related to a concurrent reduction in “natural background” water repellency that is characteristic of dry, unburned soils in the Western Cascades. We observed no evidence of spatial autocorrelation in Kfs using semivariogram analysis. Principal component analysis paired with a k-means cluster analysis suggested that soil physical properties explained variations in soil hydraulic properties better than landscape attributes. Although there is a lack of regional results for comparison, our results trend in the opposite direction from drier, lower net primary productivity regions that are typically studied for post-wildfire soil hydraulic properties.

Hydrogeologic framework and conceptual groundwater-flow model of the panhandle and northwest parts of the High Plains (Ogallala) aquifer in Oklahoma, 1998–2022

Released May 11, 2026 11:05 EST

2026, Scientific Investigations Report 2026-5009

Amy S. Morris, Colin A. Baciocco, Isaac A. Dale, Chloe Codner, Ethan A. Kirby, Grant M. Graves, Derrick L. Wagner, Eric G. Fiorentino, Alan LePera, Jon E. Sanford, Lara Joy

This study was conducted by the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, to update the hydrogeologic framework and conceptual flow model for the panhandle and northwest parts of the High Plains (Ogallala) aquifer in Oklahoma, which together compose the Ogallala aquifer focus area. The study included the construction of a potentiometric surface, and available geologic and hydrologic data were used to evaluate saturated thickness of the aquifer. The water budget for the updated conceptual groundwater-​flow model was based on estimated inflows and outflows for the 1998–2022 study period.

Saturated thickness of the Ogallala aquifer averaged 127 and 116 feet for the panhandle and northwest parts, respectively. Groundwater withdrawals from the Ogallala aquifer for 1998–2022 averaged 422,054 and 39,645 acre-feet per year (acre-ft/yr) for the panhandle and northwest parts, respectively. Recharge, the primary inflow, was estimated at 0.63 inch per year for the 1998–2022 study period, with the panhandle part of the Ogallala aquifer receiving 175,068 acre-ft/yr and the northwest part of the Ogallala aquifer receiving 49,376 acre-ft/yr. Additional inflows included irrigation return flows, estimated at 8,111 and 642 acre-ft/yr for the panhandle and northwest parts, respectively, of the Ogallala aquifer. Net lateral groundwater flows, considered to be aquifer outflows, were estimated to account for 31,908 acre-ft/yr for the Ogallala aquifer focus area. Streambed seepage, which was an outflow of 5,535 acre-ft/yr, was only present in the northwest part of the Ogallala aquifer. Vertical leakage and saturated-zone evapotranspiration were considered negligible outflows. These findings provide a revised conceptual groundwater-flow model water budget for the Ogallala aquifer focus area in Oklahoma.

Hazard potential of compound flooding from rainfall, storm surge, and groundwater in coastal New York and Connecticut

Released May 11, 2026 10:42 EST

2026, Natural Hazards and Earth System Sciences (26) 2169-2188

Robin L. Glas, Liv M. Herdman, Salme Ellen Cook, Archi Howlader, Kristina Kirkyla Masterson

Compound flood events, the co-occurrence of multiple flood drivers, can result in flood hazard potential exceeding that of any single driver alone. To evaluate compound flooding in a semi-urbanized coastal area, historical records dating back to 1970 are used to study the co-occurrences of high precipitation, storm surge, and shallow groundwater conditions along the coastlines of New York and Connecticut. Joint return periods for coincident precipitation-surge events were computed using statistical dependence models and compared to the assumption of independence as a ratio, referred to here as a return period adjustment. Results indicate distinct seasonality where compound events in the area disproportionately occur in the cold season between October and April. Return period adjustments range from a factor of 1 to almost 9, demonstrating the range in precipitation-storm surge dependence across the study area. Across all 24 station triad locations, groundwater levels were elevated during times of precipitation- surge co-occurrence, reflecting the tendency for coastal storms and shallow groundwater conditions to co-occur seasonally. The result is a pseudo-trivariate compound flood hazard score and corresponding hazard map that integrates dependence between daily precipitation-surge events and overall monthly groundwater levels (as a precondition) into a relative compound hazard score. The location with the highest compound flood hazard score is on the south shore of Long Island, as well as locations across coastal Connecticut where groundwater levels compound the co-occurrence of heavy precipitation and storm surge.

Quantitative mineral resource assessment of lithium pegmatite deposits in the southern Appalachian orogen

Released May 11, 2026 10:40 EST

2026, Natural Resources Research

Joshua Mark Rosera, Kelsey Elizabeth Crocker, Laura Pianowski, Jacob T. Murchek, Ashton M. Wiens, Margaret M. Sanders, Lucas Leonidus Evart, Jacob DeAngelo, Graham W. Lederer, Joshua A. Coyan

The first quantitative mineral resource assessment for undiscovered lithium pegmatite deposits in the southern Appalachian region of the United States was conducted. Permissive tracts for lithium pegmatite deposits were delineated by integrating lithological, tectonic, geochemical, geophysical and mineral occurrence data. Lithium pegmatite prospectivity of the tracts was ranked with simplified mappable criteria, including proximity to Paleozoic felsic intrusions and major lithotectonic structures, stream sediment geochemical anomalies, and pegmatite occurrence data. The geospatial data and permissive tracts were used to estimate the number of undiscovered lithium pegmatite deposits. These estimates were integrated into probabilistic simulations along with a new global lithium pegmatite grade and tonnage dataset to quantify potential contained undiscovered lithium resources. An economic filter was applied to convert the probabilistic estimates of contained lithium into recoverable material. The identified lithium pegmatite resources for the Carolina Lithium and Kings Mountain deposits, North Carolina, contain 1589 thousand tons (kt) of Li2O. The median contained undiscovered resource for the southern Appalachian orogen was estimated to be 2240 kt Li2O. At 90% confidence, the region contains at least 130 kt Li2O, and 10,700 kt at 10% confidence. After applying economic filters, the median recoverable contained resource was 1430 kt Li2O, corresponding to approximately 201 years of current lithium imports for consumption in the United States. North and South Carolina are likely to contain most of these resources. Coarse data resolution and intra-state variations in the geological data contribute to uncertainty of undiscovered lithium pegmatite resources. Continued efforts to harmonize disparate geospatial datasets with updated or new information can improve the accuracy and precision of estimated undiscovered lithium pegmatite resources in the study area and at broader scales.

Not so fatal attraction: Captive female Burmese Python lures do not improve wild python detection

Released May 11, 2026 10:02 EST

2026, Southeastern Naturalist (25) 191-200

Alex D. Potash, Maggie Jones, Michael Kirkland, Jenna Cole, Kristen Hart, Robert A. McCleery

The ongoing invasion of Python bivittatus (Burmese Python; henceforth, Python) across the Greater Everglades Ecosystem (GEE) has led to near total collapse of the affected mammal community over the past few decades. Management efforts to eliminate Pythons and control their spread have been hampered by the Python's low detectability, which may be improved by using a lure. In controlled settings, male Pythons show an attraction to the scent from reproductively active females. To test the effectiveness of using reproductively active female Pythons as a lure for attracting wild male Pythons in the field, we conducted a paired experiment with wild-caught female Pythons in pens and empty control pens. We monitored Python visitation at all sites using camera traps, which resulted in >3,000,000 photographs that we filtered to 4 independent detections of Pythons using AI software. Python detection was low at sites with female Pythons (3 observations) and control sites (1 observation) over 90 days at 12 sites. Stress associated with captivity may have halted reproductive females from producing pheromones, eliminating the chemosensory cue that lures males. Identifying and implementing husbandry techniques to reduce stress in wild-caught female Pythons could improve the effectiveness of this technique. Little is currently known about the chemical ecology of Pythons, and pheromonal communication in particular, and further research in this area could aid in the identification and production of effective, low-cost lures to increase detection and removal of this invasive species.

Accounting for emigration reveals high survival and bimodal size at departure from a loggerhead sea turtle (Caretta caretta) foraging area

Released May 11, 2026 09:52 EST

2026, Marine Biology (173)

Caroline M. Blommel, Margaret Lamont, William L. Kendall

The life history of hard-shelled sea turtles includes several ontogenetic shifts in habitat use and these complex permanent emigration patterns can impact estimates of stage-specific population rates, including survival. We developed several multistate mark recapture models to estimate survival of adult and juvenile loggerhead turtles from a coastal bay in the northern Gulf of America (also commonly referred to as the Gulf of Mexico) while, in some cases, accounting for permanent emigration and transient individuals. Our mark-recapture dataset consisted of 228 individual turtles with 37 total recaptures from 2011 to 2024. Of the models we fit, those that incorporated emigration produced higher estimates for annual survival than models that did not, and higher estimates than what is commonly seen in the literature for loggerheads. All models suggested a major permanent emigration pulse at the typical size of sexual maturity (70 cm straight carapace length) and another major pulse at > 90 cm. This bimodal pattern of departure may reflect differences in size at sexual maturity among loggerheads, possible genetic variability within the assemblage, or both. To assess the models’ ability to effectively recover true parameter values, we developed a simulation study of 50 randomly generated independent data sets under our specified models of similar sample size to our study dataset. Simulation results suggested that models that accounted for permanent emigration and transient individuals produced relatively unbiased estimates of survival, while models that did not often underestimated survival rates. Mark-recapture studies that may exhibit emigration and suffer from low recapture rates would benefit from auxiliary data collection such as acoustic telemetry detections to better estimate true rates of emigration and survival. Obtaining unbiased estimates of true survival by accounting for processes like emigration can support effective conservation of endangered long-lived species like loggerheads.

Tropicalization of the temperate zone: Spatiotemporal variability of winter warming and declining freeze days across the United States

Released May 11, 2026 09:30 EST

2026, International Journal of Climatology

Vincent M. Brown, Derek T. Thompson, Buren B. DeFee, Michael Osland, Barry D. Keim

We investigate changes in cool-season and winter daily minimum (Tmin) and maximum (Tmax) temperatures, and the occurrence of freeze days, from 1952 to 2024 across the conterminous United States (CONUS). Emphasis is placed on the tropical-temperate transition zone (TTTz) in the southeastern CONUS. During winter, ~70% of the land area exhibited Tmin warming rates exceeding those of Tmax. The countywide coldest Tmin became milder across 57% of the CONUS, while the coldest Tmax showed little change and even cooled east of the Rocky Mountains in the central CONUS. Across the TTTz, 75% of freeze days occur within a ~25–100-day window, often fewer than 75 days in the southernmost areas. Approximately 80% of counties exhibited significant contractions in freeze-day concentration, with the largest and most spatially consistent changes occurring in the Southeast, primarily driven by later start dates. Roughly 85% of the CONUS experienced a significant decline in freeze days, with the largest relative declines in regions where average winter Tmin is above freezing, while parts of the Pacific Northwest showed no significant change. An analysis of freeze day isopleths (30, 45, 60 and 75 days) across 20-year periods showed that the mean latitude of freeze days has migrated poleward substantially. Between 101° W and 79° W in the TTTz, the 30 freeze-day isopleth for the late period (2005–2024) was, on average, 122 km (~1.1° latitude) farther north than in the early period (1952–1971). Generally, the largest latitudinal shifts and percentage losses in freeze days occurred across low-elevation, low-relief regions at lower latitudes (e.g., the Mississippi River Valley), with abrupt shifts occurring near topographic gradients. Regions with sharp elevational gradients (e.g., Balcones Escarpment, Ouachita Mountains and Tennessee Valley) exhibited smaller temporal changes, likely reflecting the barrier-like influence of higher terrain on the poleward retreat of freeze days.

Variability and consistency in wildfire susceptibility: Insights from a national compilation

Released May 11, 2026 09:19 EST

2026, International Journal of Wildland Fire (35)

Aaron Daniel Russell, Lucas Bair, James R. Meldrum, Todd Hawbaker

Background

Wildfire risk in the United States is rising and remains a land management priority. The quantitative wildfire risk assessment (QWRA) framework integrates fuels, topography, weather and values at risk to estimate the potential change in value from wildfire. Within this, response functions (RFs) represent how values respond to fire intensity. These are often based on expert judgment, but variation across assessments is unclear.

Aims

This study uses data from the US Geological Survey (USGS) Wildfire Hazard and Risk Assessment Clearinghouse to characterize consistency and variation across categories and contexts.

Methods

We applied descriptive statistics to summarize RFs, using tables, box-and-whisker plots and heat maps stratified by highly valued resource or asset (HVRA) category and spatial scale.

Key results

RFs and value definitions vary, especially for ecosystem-related resources. Some functions, such as for buildings in the wildland–urban interface (WUI), translate well across contexts, while others require more input.

Conclusions

Some functions are broadly transferable, while others need customization. This analysis provides references and starting points for improvement to RFs in QWRAs.

Interpretations

Expanding the clearinghouse and dataset and building more transparency in expert elicitation can build trust among communities, agencies and end-users, and can support efficient use of limited resources to mitigate wildfire risk.

Patterns of floodplain forest mortality and recruitment along the Upper Mississippi and Illinois Rivers: Associations with forest fragmentation and flood inundation

Released May 11, 2026 09:18 EST

2026, Landscape Ecology (41)

Nathan R. De Jager, Jason J. Rohweder, Molly Van Appledorn, Shelby A. Weiss, Matthew Trumper, Lyle J. Guyon

Context

Different rates of floodplain forest recruitment and mortality can reveal important changes in ecosystem processes that drive forest dynamics, resulting in net changes in forest cover, thereby influencing a wide range of river habitat and morphological characteristics.

Objectives

We evaluated characteristics of forest change areas in the Upper Mississippi River System.

Methods

An overlay technique was used to map patches of forest loss, gain, and persistence between 2010 and 2020 in relation to a series of explanatory variables.

Results

We quantified a net decline in forest cover ranging from 3.2 to 16.8% in the uppermost five study reaches, and a net increase in forest cover ranging from 0.5 to 4.6% in the southernmost three reaches. Patches of forest loss and persistence were similarly tall (> 15 m), dense (> 90% cover), silver maple (Acer saccharinum) dominated forests, whereas forest gain patches were short (< 15 m), less dense (< 66% cover) and more likely to be dominated by willow (Salix) species. Both forest loss and gain patches were smaller than forest persistence patches and were typically found in areas with low neighborhood forest density (< 50% forested 10 ha neighborhood). Areas that experienced more than three flood events per growing season, more than 100 consecutive days of inundation during a single flood event, and more than 60 mean total days of inundation per growing season from 2011 to 2020 showed a net loss of forest cover in all study reaches. In contrast, net increases in forest cover were restricted to areas that experienced less than a single flood event per growing season, less than 40 consecutive days of inundation during a single flood event and less than 30 mean total days of inundation per growing season from 2011 to 2020.

Conclusions

Forest mortality along these river reaches is associated with forest fragmentation and an increasingly wetter hydrological regime.

Temporal and spatial changes in seismic attenuation associated with inferred fluid migration in the 2016 central Apennines earthquake sequence

Released May 11, 2026 09:08 EST

2026, Bulletin of the Seismological Society of America

Luca Malagnini, Francesco Pio Lucente, Irene Munafo, Douglas S. Dreger, Thomas E. Parsons, Roland Burgmann

Prior work suggests that high‐frequency seismic attenuation acts as a highly sensitive proxy for crustal permeability and fluid mobility in fractured media. We test the hypothesis that the fault system responsible for the 2016–2017 Amatrice–Visso–Norcia–Capitignano sequence acted as an impermeable seal, compartmentalizing pressurized fluids until dynamic rupture triggered widespread fluid diffusion. By tracking across the sequence the spatiotemporal evolution of the S‐wave anelastic attenuation parameter, we identify large, positive low‐frequency attenuation anomalies emerging within the hanging wall following the Amatrice mainshock and strictly preceding subsequent large ruptures. Conversely, we observe weaker, negative anomalies in the footwall, anticorrelated in time with those of the hanging wall, revealing a massive asymmetry in fluid redistribution and permeability evolution across the fault system. Furthermore, aftershock migration rates reveal distinct linear alignments in a distance‐reduced time space, allowing us to explicitly track and quantify episodes of lateral and upward fluid migration. These physically consistent patterns suggest that stress‐driven fluid diffusion directly weakens adjacent fault patches, dictating the spatiotemporal migration of seismicity. We conclude that near‐real‐time monitoring of seismic attenuation may help detect fluid redistribution in active fault systems and may provide useful information for time‐dependent seismic hazard assessment.

Publisher correction: Strong nickel enrichment co-located with redox-organic interactions in Neretva Vallis, Mars

Released May 11, 2026 09:05 EST

2026, Nature Communications (17)

H. T. Manelski, R. C. Wiens, A. Broz, J. A. Hurowitz, M. Tice, S. M. Clegg, E. Dehouck, N. Randazzo, S. A. Connell, O. Forni, S. J. VanBommel, S. Scrhoder, L. Mandon, Travis S.J. Gabriel, C. Bedford, R. K. Martinez, E. A. Cloutis, A. Cousin, M. L. Cable

Correction to: Nature Communications https://doi.org/10.1038/s41467-026-70081-3, published online 31 March 2026

Refinement of a framework for Moving Aircraft River Velocimetry (MARV) and application to particle tracking along Alaskan rivers

Released May 11, 2026 08:46 EST

2026, Water Resources Research (62)

Carl J. Legleiter, Paul J. Kinzel, Mark Laker, Jeff Conaway

Information on river velocities enhances understanding flood hazards, evaluating habitat conditions, and predicting the transport of floating materials. In this follow-up study, we used data from two new sites, one with a more complex morphology and the other with a lower suspended sediment concentration, to provide further evidence that Moving Aircraft River Velocimetry (MARV) can yield accurate velocity estimates ( R2 up to 0.87 when compared to field measurements) for long segments of large, turbid rivers. The MARV workflow is packaged in freely available software and is robust to implementation details; neither buffering to mitigate edge effects nor a new approach to aggregating velocity vectors improved performance. MARV was not sensitive to parameters used to establish overlapping image sequences, but combining a long window with a short jump between consecutive windows was the optimal configuration. Although accuracy varied from one cross section to the next, agreement between remotely sensed velocities and those measured in the field was independent of position within a frame range. As an initial step toward application of the approach to help address practical problems, we showed how MARV can drive particle tracking models. Our first-order simulations suggest that channel morphology and flow velocity are the primary controls on travel time and particle fate, with diffusive processes playing a lesser role. Although MARV can be used to characterize an instantaneous flow field, a more comprehensive framework that accounts for other physical processes would be required to model specific types of events like oil spills.

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

Released May 11, 2026 08:11 EST

2026, Scientific Investigations Report 2026-5008

Philip T. Harte, Andrew L. Collins

Per- and polyfluoroalkyl substances (PFAS), including perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS), have been detected at combined concentrations above 2,000 nanograms per liter (ng/L) at groundwater seep locations near the Coakley Landfill Superfund site, in North Hampton, New Hampshire. The landfill was active from 1972 to 1985. An impermeable cap was placed on the landfill in 1998. The adjacent area to the Coakley Landfill has many water supply wells, and transport of PFAS compounds to the wells is a concern. Fracture anisotropy in the underlying bedrock aquifer complicates the understanding of PFAS transport because groundwater preferentially travels along fractures that may not align with the prevailing groundwater flow direction.

In 2018, the U.S. Environmental Protection Agency and the U.S. Geological Survey began an investigation of the groundwater flow from the Coakley Landfill site. This report describes the modification of a numerical groundwater-flow model for the local area around the Coakley Landfill and summarizes findings of the investigation. In addition, this report includes a brief description of PFOA and PFOS occurrence, a discussion of model construction, evaluation of model performance through calibration, and discussion of simulation results for two periods (before and after capping). Limitations are also discussed.

Results show that simulated groundwater flow moves from the Coakley Landfill to the west and north. Advective transport modeling using particle tracking shows that groundwater from the landfill discharges primarily to streams to the west and north, and a small amount is transported to distal wells. Dilution of contaminants through advection and dispersion likely plays a role in whether PFAS compounds from the landfill will be detected above laboratory reporting levels at distal wells.

Landscape connectivity and wildlife access to water across an international border: Barriers and opportunities for facilitating transboundary movement

Released May 08, 2026 08:50 EST

2026, Global Change Biology (32)

Bogdan Chivoiu, Erin L. Koen, Michael Osland, Christopher A. Gabler, Jerald T. Garrett, Ernesto Reyes, Stephanie A. Bilodeau, Mitch A. Sternberg, Miguel L. Villarreal, Eric K. Waller, Samuel N. Chambers, Jude A. Benavides, Robert S. Lawson, James Martinez

Rapid global acceleration in the construction of physical barriers along international borders has greatly influenced biodiversity and animal movement. Physical barriers can fragment landscapes, hinder access to essential resources, impact long-distance migrations, and inhibit dispersal and gene flow. The effects of physical barriers on animal movement and landscape connectivity can be exacerbated in dryland environments where access to water is a limiting factor. In recent decades, the construction of border barrier infrastructure has accelerated along the international boundary between the United States and Mexico. Here, we used a landscape connectivity model to investigate the effects of barriers on wildlife access to the river in the Lower Rio Grande Valley. We used a modified omnidirectional connectivity model to compare access to the river for three large, terrestrial mammal species across three border barrier scenarios: (1) a landscape without border barriers; (2) a landscape with the existing barrier system; and (3) a potential future landscape with a continuous barrier system. The existing barrier system includes many discrete sections of barrier within tracts of the Lower Rio Grande Valley National Wildlife Refuge or on lands associated with the region's flood control system. Our results indicate that the existing border barriers can impede connectivity and wildlife access to the river in some areas, while some existing gaps between border barrier sections can serve as conduits for wildlife movement and river access. Our future scenario results show how a potential continuous border barrier system could further impede wildlife access to the river. We discuss management and landscape conservation options for enhancing wildlife access to water and riverine habitats. Collectively, our results illustrate the potential effects of border barriers on wildlife movement and access to water, providing information that can be used to better anticipate and lessen the ecological impacts of transboundary barriers.

Regression models for estimating suspended sediment concentrations and loads and comparison with acoustic surrogate model on the Snake River, Weiser, Idaho, 1977–2022

Released May 07, 2026 15:45 EST

2026, Scientific Investigations Report 2026-5007

Megan K. Kenworthy

The U.S. Geological Survey, in cooperation with Idaho Power, developed streamflow- based regression models to estimate suspended sediment concentration (SSC) and loads on the Snake River at Weiser, Idaho site (U.S. Geological Survey streamgage 13269000; hereafter referred to as “Snake at Weiser site”). This site sits upstream from the dams and reservoirs of the Hells Canyon Complex and the Hells Canyon National Recreation Area, where large sandbars along the Snake River that provide recreation and riparian habitat and host archaeological resources have declined since 1973. Analyses of samples from historical (1977- 2003) and modern (2017- 22) periods show that SSC has decreased over time, with median concentrations declining from 50 milligrams per liter (mg/L) to 28 mg/L. Mann- Kendall trend tests confirm statistically significant declines in total SSC and the fine and sand fractions of suspended sediment through the full period of record.

Regression models specific to each period outperformed models using the full dataset, suggesting changes in the sediment supply to this reach of the Snake River and highlighting the need for period- based approaches. Regression models for total SSC and fine sediment were more accurate than those for sand, which exhibited greater error and bias, likely reflecting a sand supply limited by upstream dams. The regression model for modern period total SSC and a previously developed acoustic surrogate model showed similar performance, indicating both methods are viable for estimating SSC and loads.

These findings help to better quantify suspended sediment concentrations and loads upstream of the Hells Canyon Complex and provide resource managers with tools to better quantify sediment loads affecting reservoir storage and the maintenance of sandbars in the Hells Canyon National Recreation Area.

Fossil footprints and Ice Age ecosystems of White Sands National Park

Released May 07, 2026 13:50 EST

2026, Fact Sheet 2025-3046

Kathleen B. Springer, Jeffrey S. Pigati, David Bustos, Thomas M. Urban, Matthew R. Bennett

Introduction

In September 2021, National Park Service staff, U.S. Geological Survey scientists, and an international team of researchers revealed evidence in the form of human footprints at White Sands National Park, New Mexico, that showed people were present in North America between 23,000 and 21,000 years ago. This time was during the Last Glacial Maximum, when large ice sheets covered much of the continent. The results stunned the scientific community and sparked a global debate. The story of how the discoveries were made, how they upended traditional thought, and how they “rewrote the book” on the earliest phases of North American prehistory is a classic example of the process of science.

Water use in Louisiana, 2020

Released May 07, 2026 09:31 EST

2026, Scientific Investigations Report 2026-5135

Angela L. Robinson

The U.S. Geological Survey (USGS), in cooperation with the Louisiana Department of Transportation and Development, collected water-withdrawal and water-use data from a 2020 inventory of water withdrawals in Louisiana. In 2020, approximately 8,700 million gallons per day (Mgal/d) of water was withdrawn from groundwater and surface-water sources in Louisiana, which represented a 0.22-percent decrease from 2015. Total groundwater withdrawals were about 1,900 Mgal/d, an increase of 7.1 percent from 2015, and total surface-water withdrawals were about 6,800 Mgal/d, a decrease of 2.1 percent from 2015 to 2020.

Total water withdrawals, in million gallons per day, in 2020 for the various categories of use were as follows: public supply, 720; industry, 2,100; power generation, 4,100; rural domestic, 39; livestock, 7.0; rice irrigation, 930; general irrigation, 250; and aquaculture, 590. From 2015 to 2020, Louisiana’s total withdrawals for public supply increased by 1.4 percent, industry decreased by 2.3 percent, power generation decreased by 4.9 percent, rural domestic decreased by 1.2 percent, livestock increased by 11 percent, rice irrigation increased by 13 percent, general irrigation increased by 12 percent, and aquaculture increased by 20 percent.

About 51 percent (approximately 960 Mgal/d) of all groundwater withdrawn was from the Chicot aquifer system and 24 percent (approximately 450 Mgal/d) was withdrawn from the Mississippi River alluvial aquifer. Since 2015, withdrawals from the Chicot aquifer system increased by 13 percent, and withdrawals from the Mississippi River alluvial aquifer increased by 18 percent. About 72 percent (4,900 Mgal/d) of all surface water withdrawn was from the Mississippi River main stem. This value represents a 1.1-percent decrease in withdrawals from 2015 to 2020.

All water-withdrawal and water-use data presented in this report should be considered estimates. Because of rounding, totals and percentages presented in the tables, figures, and text in the report may differ slightly from totals or percentages calculated individually.

Reconstructing ancient sedimentary source-to-sink systems – Examples from southern Laurentia’s Proterozoic accretionary orogens

Released May 07, 2026 09:29 EST

2026, GSA Bulletin

Ian William Hillenbrand, Kelly David Thomson

Provenance analysis is a powerful tool for investigating sediment delivery networks, constraining magmatic histories, and reconstructing the tectonic evolution of orogenic belts and basins. Basin analysis studies increasingly use detrital zircon (DZ) U-Pb forward mixture modeling to enhance provenance interpretations by quantifying the relative contributions of different sources. Forward mixture modeling requires significant a priori knowledge that limits deep-time applications. This challenge is overcome with an inverse mixture modeling approach non-negative matrix factorization to reconstruct the number and age distributions of paleo-source regions of Proterozoic metasedimentary rocks in the southwestern United States. This analysis indicates eight reconstructed end-member distributions representing unique sediment sources: two multi-modal end members characterized by ages older than ca. 1.8 Ga from cratonic Laurentia, five unimodal age distributions between ca. 1.80 Ga and 1.65 Ga consistent with Paleoproterozoic arc magmatic sources, and a ca. 1.6−1.5 Ga end member likely derived from exotic cratons in supercontinent Nuna (Columbia). Sediments deposited between ca. 1.80 Ga and 1.73 Ga yield heterogeneous age distributions suggesting multiple arc-backarc systems and several phases of slab roll back, contraction, and accretionary orogenesis, including input from pre−1.8 Ga Laurentian cratons. Homogenization of DZ signatures during the Yavapai orogeny (ca. 1.72−1.68 Ga) reflect crustal assembly as well as the uplift of Paleoproterozoic arcs in the orogenic hinterland. Detrital zircon age distributions from strata deposited during the Mazatzal orogeny (ca. 1.65−1.60 Ga) suggest the Mazatzal Province is a continental arc constructed on older crust. Mesoproterozoic samples are consistent with multiple basins derived from local recycling and long-distance sediment transport. Collectively, these data record the tectonic transition from the episodic accretion of disparate crustal domains to an increasingly integrated continental margin. These results provide new insights into the Proterozoic tectonic and paleogeographic evolution of the southwestern United States at basin to orogen scales and highlight the power of inverse DZ modeling to extract geologically meaningful quantitative mixture models from sedimentary records alone, offering a powerful tool for deep-time tectonic and basin analysis.

Riverine pesticide trends in the United States: Assessing a decade of national-scale monitoring

Released May 07, 2026 09:13 EST

2026, Environmental Science & Technology Water (ES&T Water)

Megan E. Shoda, Sara E. Breitmeyer, Elise Danica Hinman, Sarah M. Stackpoole

Pesticides in freshwater systems can compromise water availability by degrading water quality, with implications for human health and aquatic life. Despite recognition of the need for national-scale monitoring and analysis, few studies have documented long-term trends in surface water pesticide contamination across the US. This study addresses that need by analyzing temporal trends and acute and chronic benchmark exceedances for aquatic life and human health from 81 river sites sampled from 2013 to 2022 using an analytical method targeting 80 pesticides. The majority (79%) of single site and pesticide combinations had too few pesticide detections to estimate trends. When detections were more frequent, increasing trends in concentration were twice as common as decreasing trends. Increasing pesticide concentrations were common in primary drainages of the Mississippi River Basin. Aquatic life benchmarks were exceeded by 19 pesticides, and exceedances were geographically widespread, with both acute and chronic aquatic life benchmark exceedances at 62% of sites. The herbicides atrazine and metolachlor and the insecticide imidacloprid were identified as the greatest threats to surface water availability based on their trends and aquatic life benchmark exceedances. These findings demonstrate the need for continued monitoring and trend analysis, driver investigation, and management strategies to protect freshwater resources.

Acoustic measurements and modeling of sub-bottom layers on the summit of the Atlantis II Seamount

Released May 07, 2026 09:05 EST

2026, JASA Express Letters (6)

Tzu-Ting Chen, Matthew A. Milone, Jason Chaytor, James H. Miller, Gopu R. Potty, William S.. Hodgkiss, Ying-Tsong Lin

Sub-bottom profiler images taken from the summit plateau of the Atlantis II Seamount reveal distinct seabed stratigraphy, including marine sediment, limestone, and basalt layers. Acoustic reflection data also show arrivals reflecting from this sub-bottom structure. A wavenumber integration model with elastic geoacoustic properties is able to reproduce the arrival pattern of seabed reflections and particularly the phase inversion of the sub-bottom return. The reflection model suggests that the limestone layer is eroded with high porosity and possesses a lower compressional velocity than a well-cemented layer. The model results also highlight the necessity of incorporating elastic effects for realistic geoacoustic characterization.

Life history traits and population dynamics of Freshwater Drum across large river gradients

Released May 07, 2026 09:05 EST

2026, Transactions of the American Fisheries Society

Kristen L. Bouska, Levi E. Solomon, Andrew Bartels, Steven A. DeLain, Eric J. Gittinger, Travis Kueter, Kristopher A. Maxson, John L. West, James T. Lamer, Hae H. Kim, Quinton Phelps

Objective

Monitoring and assessment of nongame native fishes is limited, but conservation interest in these species is growing. Freshwater Drum Aplodinotus grunniens are a wide-ranging species that serve important functional roles and could serve as an indicator for similar but less common species. Our overall objectives were to quantify and compare population dynamic rates and life history of Freshwater Drum among study reaches in the upper Mississippi and Illinois rivers and relate these metrics to hypothesized environmental and anthropogenic factors.

Methods

We integrated recently collected age data with monitoring data to estimate age and size distributions, growth curves, maturation schedules, mortality rates, and young-to-adult ratios of Freshwater Drum in six study reaches spanning 1,500 km of river. Principal component analyses and linear regression were used to relate environmental and anthropogenic gradients (latitude, commercial harvest, hydrologic dynamics, primary productivity) to life history traits and population dynamic rates.

Results

We found latitudinal gradients in life history traits and population dynamic rates whereby Freshwater Drum in upstream, higher-latitude study reaches generally exhibited later maturity, slower growth, smaller maximum size, and lower mortality rates compared with those in lower-latitude study reaches. Further, young-to-adult ratios positively corresponded with chlorophyll-a concentration. No clear relationships were apparent between population dynamic rates and hydrologic variation or commercial harvest.

Conclusions

Latitude is an important structuring component of life history traits and population dynamics of Freshwater Drum in the upper Mississippi and Illinois rivers likely due to both temperature seasonality and disturbance regimes. The presence of demographic structure in a widespread, common species such as Freshwater Drum suggests similar patterns likely exist in other long-lived native fishes.

Los Planes watershed vegetation monitoring: Standard operating procedures

Released May 07, 2026 08:44 EST

2026, Report

Natalie R. Wilson

This is a description of survey procedures for short term vegetation monitoring at Natural Infrastructure in Dryland Stream (NIDS) structure sites and control sites a ranch in the Los Planes, La Paz, Baja California Sur. This study design was modified from USGS Short Term Vegetation Response Study (Wilson et al. 2021) with the goal to quantify changes in species abundance/cover, structure, and composition. The Society of Ecological Restoration identifies 3 major ecosystem attributes of importance when assessing restoration projects, such as the installation of NIDS (Society for Ecological Restoration International Science & Policy Working Group 2004; Ruiz-Jaen and Mitchell Aide 2005). These attributes are vegetation structure, diversity, and ecological processes. Our protocol can be used to directly quantify vegetation structure and diversity and by collecting data over several years we can indirectly assess the ecohydrological processes associated with NIDS (Norman, Lal, et al. 2022).