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.

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.

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.

Released May 06, 2026 13:04 EST

2026, Open-File Report 2026-1009

Allison M. Nguyen, Jonathan P. Rose, Anna C. Jordan, Giancarlo R. Napolitano, Daniel Macias, Elliot J. Schoenig, Gabriel A. Reyes, Brian J. Halstead

The giant gartersnake (Thamnophis gigas) is a semi aquatic snake endemic to the Central Valley of California. After losing 95 percent of its historic wetland habitat (Frayer and others, 1989), giant gartersnakes became state and federally listed as a threatened species (California Fish and Game Commission, 1971; U.S. Fish and Wildlife Service 1993, 1999). Continued monitoring of current populations and implementation of suggested management actions is necessary to recover the species. The Natomas basin in Sacramento, California, supports a population of giant gartersnakes persisting in restored marshes and rice agriculture. This annual report summarizes the giant gartersnake monitoring project for 2024, focusing on the apparent survival, abundance, density, and distribution of the giant gartersnakes and the connectivity of habitat throughout the Natomas basin. In 2024, 131 giant gartersnakes were captured 216 times at 44 sites by hand or trap. The catch-per-unit effort decreased from 2023 to 2024 but was similar to other years of the study. Estimates of occupancy increased between 2023 and 2024, although the trend of occupancy from 2011 through 2024 is still decreasing overall at a mean annual rate of 3 percent per year. Apparent survival was much higher at Betts-Kismat-Silva from 2018 to 2019 and from 2021 to 2022 than in other years, but this may be partly attributed to different sampling efforts over the years. Trapping effort was more consistent in the Sills tract, and apparent survival was slightly higher in later years (2022–23 and 2023–24). Giant gartersnake populations appeared to remain stable in 2024, but abundance, density, survival, and distribution is highly variable across different sites and years of the study. Continued monitoring of the populations would allow for better trend estimates over time and assessment of the effects of management activities. Giant gartersnake populations throughout the basin and on reserve lands would likely benefit from the following: (1) creating more managed marsh; (2) increasing the amount of emergent tule vegetation in existing marshes (for example, Cummings, Natomas Farms, and Lucich South); (3) continuing to flood existing marshes in early spring; (4) maintaining rice agriculture; and (5) continuing research into conservation actions that target the giant gartersnake, such as habitat and water management and translocation.

Released May 06, 2026 11:50 EST

2026, Fact Sheet 2026-3004

Rand Gardner, Justin E. Birdwell, Jason A. Flaum, Scott A. Kinney, Janet K. Pitman, Stanley T. Paxton, Andrea D. Cicero, Jenny H. Lagesse, Jeffrey D. Pepin, John W. Counts, Benjamin G. Johnson, Celeste D. Lohr, Katherine J. Whidden, Katherine L. French, Tracey J. Mercier, Heidi M. Leathers-Miller

Using a geology-based assessment methodology, the U.S. Geological Survey estimated undiscovered, technically recoverable mean resources of 3 million barrels of oil and 343.5 trillion cubic feet of gas in reservoirs of the Bossier Formation within the onshore United States and State waters of the Gulf Coast region.

Released May 06, 2026 10:02 EST

2026, Fisheries

Matthew D. Robertson, Holly Susan Embke, Abigail Lynch, Stephen R. Midway, Craig Paukert

Recreational fisheries are important global contributors to food security, socio-cultural practices, and local and regional economies. However, inland recreational fisheries are often overlooked by policymakers due to a limited understanding of the magnitude of participation, harvest, and economic impact. Here, we used the U.S. Inland Creel and Angler Survey Catalog and catch and effort model (CreelCatch) and several assumptions to provide an initial estimate of the magnitude of total inland recreational fisheries harvest in the conterminous USA. The CreelCatch model projected fishing harvest across lakes, ponds, and reservoirs based on fishing effort, water body area, and regional effects. We estimated that recreational lake fisheries in the conterminous USA likely harvest 236,000–671,000 tonnes of fish per year, 17–48 times greater than total inland fisheries harvest reported to the United Nations. Inland recreational fisheries may warrant greater consideration for their contribution to national scale socioeconomics and impacts on fish stocks and ecosystems.

Released May 06, 2026 09:33 EST

2026, Scientific Reports (16)

Jessica Z. LeRoy, Grace Loppnow, P. Ryan Jackson, G. Everett Lasher

Invasive carp pose ecological and economic risks to North American freshwater systems. This study uses the Fluvial Egg Drift Simulator to model the drift of invasive silver carp (Hypophthalmichthys molitrix) eggs and larvae after hypothetical spawning in Pools 1–10 of the Upper Mississippi River. Although adult invasive carps have been detected in this region, no reproduction has been confirmed as of this publication. A total of 450 spawning scenarios were simulated, representing 5 water temperatures, 9 flows, and 10 spawning locations in the tailwaters of lock and dam structures. The study examined egg and larval positions at two key developmental stages: hatching and gas bladder inflation, when larvae seek nursery habitat. Under a wide variety of flow conditions and water temperatures, eggs spawned upstream from Lake Pepin (Pool 4) are likely to settle in the lake before hatching, possibly increasing mortality rates. Eggs that survive passage through Lake Pepin reach gas bladder inflation within the study area, except in scenarios with lower temperatures and higher flows. Conversely, larvae spawned downstream from Lake Pepin generally drift out of the study area before reaching gas bladder inflation, except in cases of higher temperatures and lower flows. These findings inform ichthyoplankton sampling strategies and management actions aimed at reducing invasive carp populations in areas likely to support recruitment.

Released May 04, 2026 11:50 EST

2026, Scientific Investigations Report 2026-5002

Thomas P. Suro, Michal J. Niemoczynski, Kevin B. Mulligan

As the population of New Jersey continues to remain dense, the need for water supply will likely continue to be high, which can lead to water managers needing to make difficult decisions about managing drinking-water supply. Streamgaging weirs like the ones used by the U.S. Geological Survey (USGS) play a critical role in providing accurate and stable streamflow data, but their presence can affect the passage of diadromous fish species such as river herring (Alosa pseudoharengus [alewife], Alosa aestivalis [blueback herring], and Alosa sapidissima [American shad]). In some situations, weirs existing in rivers and streams are no longer used because they were part of a farm irrigation system or some type of industrial operation. The weir at the USGS streamgage 01402000 Millstone River at Blackwells Mills, New Jersey, was purposefully built as a hydraulic-control structure that provides a precise and stable control for the measurement of stage and computation of continuous streamflow. To satisfy the dual need of maintaining accurate streamflow data and providing improved fish passage for select species of fish during migration season, the USGS proposed the development and evaluation of two alternative weir designs that would meet the criteria established for successful passage of American shad, alewife, and blueback herring during their yearly migration. The designs were also required to maintain adequate control of the upstream pool elevation necessary for the precise computation of streamflow used by State agencies for municipal water-supply purposes for surrounding communities.

Two alternative weir design modifications were incorporated at the center of the Blackwells Mills weir and modeled using two-dimensional hydraulic modeling software and three-dimensional computational fluid-dynamics software to simultaneously evaluate conditions for passage of the target fish species and effects to streamflow computations at the streamgage. The models were calibrated to existing conditions around the weir location using surveyed-elevation data and recorded stage, streamflow, and velocity in the Millstone River. The alternative weir designs lowered the weir crest by 1.02 feet (ft) and the resulting simulations showed an effective increase in depth of 0.98 ft at the median streamflow of 251 cubic feet per second (ft³/s) and 0.96 ft at the 95-percent exceedance streamflow of 98 ft³/s. The alternative weir designs were also found to increase streamflow depth across the shallowest portions of the weir structure at the downstream anti-scour skirt by lowering the skirt about 4 inches, allowing for two or more body depths of water for American shad, alewife, and blueback herring at the median migration streamflow of 251 ft³/s. The alternative weir designs also reduced the highest stream velocities across the downstream weir sill and anti-scour skirt from about 9 to 10 feet per second, and the depth-averaged velocity to about 7 to 8 feet per second. The sensitivity of the weir with respect to the computation of streamflow was increased from about 1.8 cubic feet per second per hundredth foot to 1.6 cubic feet per second per hundredth foot for streamflows of about 10–100 cubic feet per second.

Released May 04, 2026 11:00 EST

2026, Open-File Report 2026-1010

Arthur J. Merschat, Mark W. Carter, Ashley S. Lynn, Benjamin R. Weinmann, William E. Odom, Ryan J. McAleer, Shannon A. Mahan, Kevin G. Stewart, Christopher S. Holm-Denoma, E. Allen Crider, Jr.

Introduction 

New bedrock and surficial geologic mapping in the Sparta East, Sparta West, and parts of the Glade Valley and Whitehead 7.5-minute quadrangles, North Carolina and Virginia, investigates the geologic framework and causative mechanisms of the August 9, 2020, Mw 5.1 earthquake near Sparta, North Carolina. The mapping documents (1) the coseismic surface rupture from the 2020 earthquake and related brittle structures in the bedrock; (2) the fault contact between the western Blue Ridge and eastern Blue Ridge; (3) lithostratigraphy in the Lynchburg Group, Ashe Metamorphic Suite, and Alligator Back Metamorphic Suite; (4) the nature of the contact between the Lynchburg Group, Ashe Metamorphic Suite, and Alligator Back Metamorphic Suite; and (5) surficial deposits.

Released May 04, 2026 09:20 EST

2026, Wildlife Society Bulletin

Anastasia Couvillon, Gregory J. Soulliere, David H. Gordon, Diane Eggeman, Mohammed A Al-Saffar, Dale D. Humburg, James E. Lyons

Prioritization is a central component of natural resource management because conservation needs routinely exceed available resources. Waterfowl and wetland conservation programs in North America are at the forefront of landscape-scale prioritization and transboundary management decisions due to the migratory nature of ducks, geese, and swans. The growing availability of geographic information systems (GIS) and geospatial technologies has accelerated the development of multi-objective landscape prioritization models, including applications of structured decision making and multi-criteria decision analysis to spatial planning for waterfowl and wetlands at the continental scale. However, regional managers and conservationists could benefit from flexibility in downscaling continental tools, selecting objectives, and assigning weights for rapid production of spatial prioritization models at smaller spatial scales without extensive computer coding or GIS analysis. We developed a spatial value model that prioritizes landscapes at sub-continental scales (e.g., states and provinces, bird conservation regions, etc.) and provides flexibility for users to select waterfowl conservation objectives of interest and weights. Our model can be used for direct downscaling of an existing continental geospatial model or further customized with region-specific geospatial data. We illustrate how regional prioritization can vary with the spatial scale selected by the user. The spatial value modeling framework and the downscaling tool presented here could increase the use of multi-criteria decision analysis and linear value modeling in spatial landscape prioritization, while also providing flexibility for selecting scales, objectives, and weights. Our spreadsheet tool was developed specifically for use by regional biologists, conservationists, and managers and does not require knowledge of GIS software (although results can be exported from the spreadsheet for spatial analysis using GIS). Together, the model outputs and the accompanying spreadsheet tool provide a bridge between continental waterfowl conservation and regional implementation, enabling rapid, stakeholder-driven, value-explicit prioritization.

Released May 04, 2026 09:09 EST

2026, Fisheries Management and Ecology

Kurt C. Heim, Jonah L. Withers, William Arden, Laurie Earley, David Minkoff, Theodore Castro-Santos

In the Boquet River (NY, USA) a low-head dam set above a ~200-m bedrock cascade was removed in 2015. We used radio-telemetry to assess landlocked Atlantic salmon passage at the remaining cascade (2020, 2022). Across years, 52% of males (13/25) attempted cascade passage whereas females made no discernable attempts (0/11). Attempt probability increased with stream discharge and decreased with fish size, though overall passage success was low (1/36). Shallow depths—likely owing to an artificially widened channel—appear to be limiting passage. Additionally, we transported fish upstream but observed high fallback (72%) that was associated with fish size and energetic status. Following dam removal, this cascade continues to limit upstream passage resulting in increased vulnerability to angling during migratory delay. Overall, we highlight the importance of follow-up studies after dam removal, and that further modifications at this site may be required to improve passage.

Released May 01, 2026 12:11 EST

2026, Fact Sheet 2026-3065

Gardner C. Bent, Brendan A. McCarthy, Luke P. Sturtevant, Meghan A. McCallister, Amanda L. Tudor, Ian P. Armstrong, Mark W. Poe, Alexander P. Graziano, Carl S. Carlson

Introduction

Currently (2026), many of the about 25,000 roadway crossing structures over rivers and streams in Massachusetts are undersized. Undersized culverts and bridges can be detrimental to fish and wildlife movement, habitat continuity, and the health of aquatic organisms. Undersized culverts also can lack the resiliency needed to withstand large floods, which could be worsened by potential increases in flood magnitude and frequency due to climate change. Improving culvert and bridge designs for stream crossing projects may improve aquatic organism passage, stream continuity, and resiliency during future floods by decreasing upstream overbank flooding, road flooding and erosion, and degradation of aquatic habitat.

The U.S. Geological Survey (USGS), Massachusetts Department of Environmental Protection (MassDEP), and University of Massachusetts Amherst began a series of cooperative studies in July 2019 to develop an automated geographic information system (GIS) hydraulic modeling tool for preliminary culvert designs for stream crossings. The USGS plans to provide preliminary culvert designs in the web-based StreamStats application, which enables municipalities and engineers to view potential designs and related information for stream crossing replacement projects in Massachusetts. This application can (a) provide information on hydrology, hydraulics, and ecological conditions at stream crossing sites, (b) provide users with potential culvert designs to improve aquatic organism passage and flood resiliency, and (c) assist MassDEP in implementing the Massachusetts Wetlands Protection Act regulations for stream crossing projects.

Released May 01, 2026 10:11 EST

2026, ISPRS International Journal of Geo-Information (15)

Christina D. Perella, Jelena Vukomanovic, Caleb R. Hickman, Adam J. Terando, Mitchell J. Eaton, Marie Schaefer

A common challenge when modeling social–ecological systems (SESs) is defining the spatial extent of the system. Boundaries that do not adequately capture both social and ecological processes and their interactions can lead to mischaracterization of the system, while expanding boundaries too widely can impact model complexity and required resources. Socially, boundaries can invoke and influence identity, culture, power, and sense of place. Boundary decisions benefit from flexible, iterative approaches and the expertise of local communities. Here, we use a structured database search supplemented with citation searching to identify and review the literature that addresses choosing or defining spatial boundaries in SESs mapping or modeling and, when applicable, how participatory methods were used in the research process. In a review of the resulting 79 studies, we discovered that pre-existing social or ecological boundaries were used most frequently (36 and 18 publications, respectively). Twenty-one publications combined social and ecological boundaries or data to create custom boundaries, and four studies used an alternative approach to conventional boundaries. Informed by the literature review, we present a general framework for defining boundaries at the outset of SES research. We then connect the framework to a specific case study based on a collaborative project with Tribal, university, and federal scientists to develop a social–ecological climate adaptation plan. We present guiding questions alongside candidate boundaries for our study system and explore the tradeoffs of these boundary options, which can function as a useful template for other social–ecological research collaborations.

Released May 01, 2026 09:02 EST

2026, Reviews of Geophysics (64)

Anna Kelbert, Paul A. Bedrosian, Adam Schultz, Gary D. Egbert, Louise Pellerin, Jeffrey J. Love, Andy Frassetto, Benjamin S. Murphy

The United States Magnetotelluric Array (USMTArray) data set, collected in the years 2006–2024, consists of more than 1,700 long-period magnetotelluric stations covering the entirety of the contiguous United States on a quasi-regular 70 km grid. Funding across multiple federal agencies was critical to sustaining this effort to its completion. Important components of the project included active guidance and participation from the MT community, the open and timely availability of all data, and the application of consistent instrumentation and robust data processing. Together with parallel advancement in the development of publicly available three-dimensional (3D) inversion codes, the USMTArray has revitalized the US magnetotelluric community and increased the visibility of magnetotellurics within the Earth-science community. Taken as a whole, these data are visualized as the National Impedance Map, which, together with a 3D synthesis conductivity model of the nation, reveals the electrical architecture of the contiguous US. USMTArray data are used by researchers worldwide for fundamental and applied studies, including investigations of continental architecture and evolution, estimation of hazards to critical infrastructure due to geomagnetic storms, and assessment of the nation's undiscovered geothermal and mineral resources. We here review the history and development of the project, discuss the challenges and successes in its execution, present the National Impedance Map and synthesis conductivity model, and highlight the breadth of research stemming from this rich data set.

Released April 30, 2026 15:25 EST

2026, Scientific Investigations Report 2026-5018

Molly L. Schreiner, Kristin M. Romanok, Jacob T. Gray, Eileen J. Brown, Brianna M. Williams, Maureen Kneser, Albert J. Capuzzi, Jason Boerner, Luke Giunta, Paul Serillo, John J. Trainor, Kelly L. Smalling

This study, completed by the U.S. Geological Survey in cooperation with the North Jersey District Water Supply Commission (NJDWSC), was designed to characterize the occurrence and distribution of per- and polyfluoroalkyl substances (PFAS) in surface waters of the nontidal Passaic River Basin in New Jersey that have the potential to affect public-drinking-water quality. In 2025, 37 sites in the Wanaque, Ramapo, Pompton, and Passaic River watersheds were sampled in January, March, July, and September under base-flow conditions and a subset of sites was sampled during two rain events. Samples were analyzed for 40 individual PFAS and total organic carbon and a subset of samples was analyzed for 1,4-dioxane and trace elements. Fifteen PFAS were detected at least once, with individual concentrations ranging from 0.42 to 28 nanograms per liter (ng/L; median, 2.8 ng/L). Perfluorooctanoate (PFOA) and perfluorooctane sulfonate (PFOS) were widespread and detected in 100 and 97 percent of the samples, respectively. Concentrations of PFOA and PFOS ranged from 1.2 to 28 ng/L (median, 7.7 ng/L) and from 0.52 to 12 ng/L (median, 3.8 ng/L), respectively. Generally, concentrations were lower in the Wanaque and Ramapo River watersheds compared to the Pompton and Passaic River watersheds. Concentrations of PFOA and PFOS were highest in July and September when flows were low. During rain events, median concentrations of PFOS were elevated compared to those observed under base-flow conditions, indicating potential inputs from non-point sources. To understand potential drivers of PFAS concentrations, land cover and potential PFAS sources were summarized for each sampling site, and an accumulated wastewater model was used to estimate the percentage of wastewater from upstream municipal and industrial sources in all flowlines of the Passaic River Basin. Developed land, the number of potential sources, and the mean-annual accumulated wastewater percentage were highly correlated with PFAS concentrations and Deciduous Forests were negatively related to concentrations. Data provided by this study can be used by water purveyors and resource managers to make treatment and mitigation decisions to minimize PFAS in local surface waters used as drinking-water resources.

Released April 30, 2026 10:23 EST

2026, Scientific Investigations Report 2026-5142

William H. Asquith, Elena R. Crowley-Ornelas, Amanda R. Whaling

The U.S. Geological Survey, in cooperation with the Gulf Coast Ecosystem Restoration Council, assessed monotonic trends for a variety of streamflow statistics for 69 long-term U.S. Geological Survey streamgages within either the Mobile Bay or Perdido Bay watersheds that were active through at least at the end of calendar year 2019. Long-term data were defined for this investigation as having at least 50 years of cumulative record within the period since January 1, 1950, with a requirement for a complete record of streamflow during the 2010s (2010–19). The 69 streamgages have at least 54 years and as many as 73 years of daily mean streamflow data; the median period of record is 72 years; and 15 of the streamgages are identified as “major nodes” on the basis of the criteria described. The occurrence of statistically monotonic significant trends for the 69 streamgages at the 0.05 significance level is spatially shown for six statistics. For the major node streamgages, the study depicts (1) time-series graphics of annual mean, annual harmonic mean, decadal 10th, 50th, and 90th-percentile streamflows, and (2) a variation on Quantile-Kendall plots of Kendall’s tau and streamflow nonexceedance probabilities for each of the 365 days of a year. Trend assessment synthesis shows that, except for a few streamgages with relatively greater counts of statistically significant trends than others, the majority (about 93 percent) of individual trend tests indicate no trend in the streamflow and ecological metrics considered.

Released April 30, 2026 09:48 EST

2026, Forests (17)

Jesse Gray, Mandy L. Slate, Alyson Ennis, Courtney Peterson, John B. Bradford, Adam Roy Noel, Michael C. Duniway, Tara B. Bishop, Ian P. Barrett, Chris Domschke, Joel T. Humphries, Nicole N. Barger

Pinyon–juniper (PJ) woodlands, one of the most extensive mature and old-growth woodland types in the Western United States, provide critical ecological, cultural, and economic benefits but face increasing threats from climate change, altered disturbance regimes, invasive species, and pests. We developed the PJ Woodland Climate Adaptation Management Menu, a decision support tool designed to guide adaptive, climate-informed management of PJ ecosystems, particularly within the Colorado Plateau ecoregion. The menu was created through an iterative, collaborative process involving literature review, integration of strategies from existing adaptation frameworks, and extensive input from scientists, land managers, and community partners during workshops and focus groups. The menu links specific, evidence-based approaches to each of six broad strategies, including soliciting community input, mitigating disturbance, enhancing and maintaining biodiversity, conserving ecotones, timing actions for optimal outcomes, and accepting climate-driven changes when appropriate. It is intended for use with the Adaptation Workbook to help managers connect local goals and climate vulnerabilities to tailored management tactics. Hypothetical scenarios demonstrate the menu’s application to contrasting PJ woodland conditions, from die-off events to old-growth maintenance. Lessons learned during development underscore the value of early stakeholder engagement, cross-sector collaboration, and balancing diverse ecological objectives. This menu offers a flexible, transferable framework to strengthen climate resilience in PJ woodlands and serves as a model that could improve adaptation planning in other dryland forest ecosystems.

Released April 30, 2026 09:44 EST

2026, Polar Biology (49)

Eden Smith, Rose J. Swift, Anna Courtemanche, Feipeng Huang, Mary Margaret Pelton, Lauren Puleo, Josiah Simmonds, Matthew Waller, Hannah Walton, Casey Weissburg, Luke R. Wilde, Nathan R. Senner

Climate change is transforming the Arctic and sub-Arctic at a pace that threatens many taxa with population declines and extinction. However, some habitats–such as muskeg bogs–can serve as climatic refugia and lessen the effects of a changing climate on the species that rely on them. Hudsonian Godwits (Limosa haemastica) are a species of migratory shorebird that utilizes the muskeg bogs of Alaska and Canada to breed. Our study focused on a muskeg bog in Beluga, Alaska, USA to see if it had changed from 2009 to 2023, if the availability of Godwit nests sites in the bogs changed in concert, and if Godwit nest survival was affected by any of these changes. We found that, overall, the bog dried and became more vegetated, with the proportional cover of graminoids, shrubs, and forbs all increasing during our study. Godwit nest sites also changed, with the proportion of shrubs and graminoids around nests increasing over time. Nonetheless, these changes did not negatively impact Godwit nest survival. Instead, nest survival increased ~ 22% during our study period, and we observed no decline in the number of potential nests sites available to Godwits. Taken together, these results suggest that while muskeg bogs are changing, they are also currently acting as climate refugia for Godwits. However, it is unclear for how long muskeg bogs can continue to buffer Godwits and other species from the effects of climate change.

Released April 30, 2026 08:57 EST

2026, Water Resources Research (62)

Carl J. Legleiter, Paul J. Kinzel

Many hydrologic applications require basic information on the size and shape of river channels, but measuring cross section (XS) geometry in the field or via remote sensing can be costly and often provides only partial coverage. Given these challenges, we capitalized upon an existing data set of 46,971 XS from gaging stations to evaluate various approximations of channel shape. After screening and pre-processing these data, we fit four model types to each XS, including a new approach that involves Stacking PDFs (probability density functions) to Approximate River Channel Shapes (SPARCS). This framework produced depth estimates that closely matched field measurements, with typical cross-sectional area errors <1% and a median R² of 0.77 for comparison of observed and predicted depths. SPARCS model parameters can be interpreted in terms of channel characteristics: mean depth, asymmetry, bar convexity, and flatness of the bed. The model performed well for the XS included in the database, which was biased toward straight, uniform channels conducive to operational streamflow measurement. Neither model parameters nor accuracy were dependent on discharge. We also assessed the potential of SPARCS to fill in measurement gaps and found that although the model can help, the accuracy of inferred depths decreased as the observable fraction of the channel decreased. An important limitation of SPARCS is that mid-channel bars or multi-threaded morphologies cannot be produced. Graphical tools can help visualize how model parameters affect simulated river forms. SPARCS could facilitate satellite-based discharge estimation by providing prior information on channel shape.

Released April 30, 2026 07:52 EST

2026, Nature Communications (17)

Carl August Norlen, Sadikshya Sharma, Francisco J. Escobedo

Human settlements are increasingly being impacted by urban fires initiated by wildfires. Metrics such as area burned and number of structures destroyed are important, but research often overlooks the socio-ecological complexity of urban fires. We study the impacts of the 2025 Los Angeles fires on two communities at the neighborhood and residential parcel scales. Geospatial analyses and econometric modeling explore the relationships between urban morphology, socio-demographic factors, and home destruction. Here we show that socio-ecological characteristics and scale are key in parsing the dynamics of urban fires. Also, new socio-demographic populations are being affected and urban morphology metrics are more important than vegetation cover. Despite parallels with 19^th and early 20^th century urban conflagrations, understanding these re-emerging urban fires requires transdisciplinary approaches and unique metrics. Investigating the socio-ecological scales and dynamics of urban fires provides a valuable next step towards understanding and adapting to the risk associated with these disasters.

Released April 29, 2026 11:00 EST

2026, Scientific Investigations Report 2026-5132

Paul M. Heisig

The causes of water level fluctuations at Moreau Lake, within Moreau Lake State Park in the town of Moreau, New York, were investigated from 2016 to 2021 after lake water levels dropped between 2015 and 2016, raising concerns about the loss of a shallow swimming area at the park beach. Annual variation in precipitation records from the area did not account for the lake water level decline. Two possible causes for the low lake water levels were investigated: the increase in groundwater withdrawals from new residential development since about 2000 and seasonal changes (nongrowing and growing seasons) in precipitation.

Investigation of the potential effects of nearby groundwater withdrawals required the compilation and collection of well-log data, seismic surveys, and measurements of lake and groundwater levels, field chemical parameters, and water isotopes to define the hydrogeologic system and to estimate water use. The net result of this work was the determination that Moreau Lake is a “flow though” lake with no surface water outlet; groundwater enters the lake on the upgradient side and exits through the downgradient side, however, groundwater does not flow southward from the lake toward nearby groundwater withdrawals from the semiconfined aquifer, and thus groundwater withdrawals were unlikely to have an effect on lake water levels.

Investigation of the historic precipitation records during nongrowing (November through April) and growing (May through October) indicated that (1) nongrowing season precipitation from 2011–12 to 2015–16 was more deficient than any similar period during the past 78 years and (2) since about 2000, nongrowing seasons have been drier overall and growing seasons have been considerably wetter. Initiation of lake water level monitoring in 2016 provided an opportunity to compare seasonal precipitation with seasonal lake water level changes. Nongrowing season lake water levels are very sensitive to precipitation, such that high precipitation (40 percent above the seasonal median) resulted in a 5-foot rise in lake water level. In contrast, the growing season lake water levels are sensitive to dry conditions; for example, deficient rainfall (about 6 percent below the seasonal median) resulted in a decline of lake water levels of about 3.5 feet. However, lake water levels are insensitive to high growing season rainfall inputs (about 10 to 47 percent above the seasonal median); lake water levels consistently declined about by 0.8 feet above this range of seasonal excessive precipitation.

Released April 29, 2026 09:51 EST

2026, Scientific Reports (16)

Clint Otto, Alma Christa Schrage, Audrey Claire Lothspeich, Larissa L. Bailey, Tamara Smith, Robert Planman, Judy Cardin, Kristen S. Ellis, Bethany Dennis, Ralph Grundel

Understanding the distributions of rare species is necessary to guide monitoring and inform species recovery efforts. The rusty patched bumble bee (RPBB; Bombus affinis, Cresson) is an endangered species with an extant, known distribution centered around urban areas of the Midwestern United States. We tested a novel approach for finding undocumented RPBBs outside of urban centers and estimated the species occurrence at two scales that are relevant to management. We confirmed presence of RPBBs at 54% of the sampled 100 km² grid cells where the species was previously undocumented, expanding the species’ known distribution by 5700 km². After accounting for imperfect detection, our occupancy model estimated the number of occupied grid cells was 67 of 105 sampled grids, suggesting our methods were effective for finding undiscovered RPBB sites. Occupancy within 100 km² grids was positively related to the number of occupied neighboring units but was not related to the area of developed land within 100km² grid cells or smaller subunits (i.e. 3.14 ha patches or roadside transects). We highlight the utility of our approach for guiding future survey efforts by identifying an additional 145 grid cells where the occupancy status of RPBB is unknown but we predict a relatively high likelihood of RPBB occurrence. Our approach can be extended to find undiscovered RPBB sites in other areas and applied to other bee species where occurrence information is lacking outside of their core distribution.

Released April 29, 2026 09:32 EST

2026, Scientific Investigations Report 2026-5129

John H. Williams, Neil C. Terry, William M. Kappel, Paul M. Heisig, Robin L. Glas, Joshua C. Woda

Mudboils have been documented in the Tully Valley in southern Onondaga County, New York, since the late 1890s. Sediment-laden water from the mudboils flows into Onondaga Creek, which empties into Onondaga Lake at Syracuse 15 miles to the north. Turbidity from the mudboils has degraded the water quality of Onondaga Creek despite a series of mitigation efforts that began in the early 1990s. Turbidity mitigation actions presently (2025) being considered include creek relocation and offline sediment settling. In support of these proposed actions during 2021–23, the U.S. Geological Survey, in cooperation with the New York State Department of Environmental Conservation, U.S. Environmental Protection Agency, Onondaga Nation, Onondaga Environmental Institute, and Central New York Regional Planning and Development Board, collected and analyzed geologic, hydrologic, geophysical, and geotechnical data to characterize the shallow hydrogeology along four proposed creek-relocation paths and in the proposed offline settling basin area.

The investigation indicated that the four proposed creek-relocation paths, two east of Onondaga Creek and two west of Onondaga Creek, are underlain by sediments including muck, alluvium, mudboil deposits, alluvial-fan sand and gravel, and lacustrine fines. The proposed excavations would penetrate partially to fully saturated conditions: generally, the water table is shallow near the creek and deep on the alluvial fans. The shallowest excavation, about 5 feet below land surface, would be near the creek and primarily in alluvium, and the deepest excavation, as much as 30 feet below land surface, would be in the alluvial-fan deposits. Brackish waters would be penetrated by proposed channel excavations on the eastern side of Onondaga Creek in an area downgradient from a potentially leaking historical salt-exploration borehole and near the main mudboil area. Excavation in these areas likely would provide a continuous source of brackish groundwater to the relocated creek. Proposed channel excavations of muck, soft to very soft lacustrine fines, and mudboil-type sediments in mudboil and suspected mudboil areas would pose an excavation and slope stability challenge and would have the greatest potential to create new mudboils. Proposed channel excavations below the water table on the Rattlesnake Gulf and Rainbow Creek alluvial fans would intercept groundwater and make the constructed streambank susceptible to seepage-induced slope instability. The substantial water-level fluctuation in the sediments of both alluvial fans would aggravate the stability condition. In addition, excavation on the Rattlesnake Gulf alluvial fan would have the potential to affect water-supply springs at the toe of the fan.

The proposed offline settling basin area is in the northern part of the Rattlesnake Gulf alluvial fan. Natural and man-made diversions of Rattlesnake Gulf have resulted in saturated conditions in the general area of the proposed basin. The proposed offline settling basin would be excavated in, and berms would be constructed on, alluvial-fan deposits and lacustrine fines. In the proposed basin area, the alluvial deposits overlying the lacustrine fines are less than 10 feet thick. Excavation, berm construction, and loading of the saturated, soft to very soft lacustrine fines may be problematic and require soil strengthening.

Released April 29, 2026 09:31 EST

2026, Geochemistry, Geophysics, Geosystems (27)

Tyler E. Huth, Thure E. Cerling, David W. Marchetti, Amy L. Ellwein, Shannon A. Mahan, David R. Bowling, Benjamin H. Passey, Victor J. Polyak, Yemane Asmerom

Laminated soil carbonate rinds are a Quaternary paleoclimate archive whose isotope composition is linked to soil formation conditions. At Rio Mesa, Utah (USA), we investigated the fidelity of rind records in a river terrace setting by determining the seasonal timing of rind formation and testing for inter-record replication. We infer soil carbonate formed in the spring season, contrasting with our prior inference of summer formation at Teasdale, Utah, ≈200 km distant. This apparent discrepancy occurs because of differences in the timing of the largest annual infiltration (spring vs. summer). At Rio Mesa, modern soil data show that soil carbonate δ¹³C would have high values (−2 to 2‰ VPDB) regardless of seasonal activity of C₃ versus C₄ plants because respiration rate is a strong control. We accordingly suggest reassessment of published records interpreting soil carbonate δ¹³C only via C₃ versus C₄ plant abundance. Three rind δ¹³C and δ¹⁸O records generally replicated. Intriguingly, rind δ¹³C may inversely correlate with summer insolation, evidence for global-scale influence on soils. Rind δ¹⁸O is not as clearly correlated with published western USA paleoclimate records, potentially due to regional differences in climate and because rinds record soil-specific processes. Our results support the fidelity of the soil carbonate rind paleoarchive and suggest that because rind formation seasonality is intimately tied to infiltration seasonality, spatial transects of rind records might be used to delineate boundaries between areas dominated by spring and summer infiltration, permitting reconstruction of the geographic extent of large-scale hydrologic phenomena such as the North American Monsoon.

Released April 29, 2026 08:52 EST

2026, PLoS ONE (21)

Erica Meta Christensen, Nathan J. Kleist, David R. Edmunds, Julie A. Heinrichs, D. Joanne Saher, Ashley L. Whipple, Melia DeVivo, Cameron L. Aldridge

Chronic wasting disease (CWD) is a prion disease of cervids that spreads to uninfected individuals through direct transmission (contact with infected individuals), vertical transmission (from mother to offspring), or indirect transmission (exposure to contaminated environments). The risk of indirect transmission is unevenly distributed on the landscape, and risk levels are expected to be controlled by patterns of habitat use by infected and uninfected individuals as well as environmental properties that alter the length of time prions remain infectious and available for uptake. Despite evidence from controlled or laboratory studies identifying environmental properties likely to affect patterns of CWD prion locations on the landscape, it remains difficult to connect mechanisms to realized increased or decreased risk of disease transmission, and few studies have attempted to detect patterns of different CWD risk in different environments. Using data from GPS-collared mule deer in Wyoming that were CWD-tested annually, we constructed models predicting annual probability of disease transmission contingent on environmental properties extracted from GPS use points. We compared models that emphasized different pathways of disease transmission by including or excluding sets of covariates that described deer density, habitat selection, and covariates expected to affect prion persistence in the environment. Results indicated that key habitat characteristics often selected by mule deer, such as proximity to secondary roads, were also associated with higher risk of testing positive for CWD, which supports the hypothesis that disease risk was correlated to patterns of habitat use by deer. We also found increased risk associated with spatial properties that were not selected-for by deer, such as areas where topography collects moisture, suggesting that prion retention mechanisms also play a role in risk. Incorporating these spatially-varying risk factors into our understanding of CWD transmission and outbreak progression can support managers in designing data collection and disease management strategies.

Released April 29, 2026 08:09 EST

2026, Geophysical Research Letters (53)

Lisa Victoria Luna, Benjamin B. Mirus, Brian D. Collins, Jonathan P. Perkins

Mid-range landslide outlooks can facilitate weather-related landslide preparedness and disaster response planning, but seasonal landslide activity remains poorly quantified at continental scales. Leveraging >55,000 reported landslides from across the United States (U.S.), we used circular statistics to quantify landslide seasonality in 67 National Weather Service County Warning Areas (CWAs). We found regional differences in landslide season timing and duration, with transitions between domains variably corresponding to climate class or river basin. We assessed differences in seasonality by movement type for slides, flows, and falls, detecting apparent, but uncertain, differences between slide and fall seasonalities in 27 of 35 (77%) of CWAs with both types reported. In the Pacific Northwest, where long records exist, we found a credible shift toward a later mean landslide season in western Washington from 1990 to 2020, but no trend in western Oregon. Our results can provide emergency planners a resource to assess seasonal landslide probability nationwide.

Released April 28, 2026 09:52 EST

2026, Report

Michael J. Gould, Frank T. van Manen, Bryn Karabensh, editor(s)

No abstract available.

Released April 28, 2026 09:15 EST

2026, Water (18)

Baktybek Duisebek, Gabriel B. Senay, Talgat Usmanov, Kudaibergen Kyrgyzbay, Janay Sagin, Yerbolat Mukanov, Kanat Samarkhanov, Xuejia Wang, Sulitan Danierhan, Xiaohui Pan

Mountain regions are highly climate-sensitive, yet long-term observational evidence of elevation and seasonal climate dynamics in Central Asia remains limited. This study examines spatiotemporal trends in temperature (Tmean, Tmax, Tmin, and diurnal temperature range [DTR]) and precipitation across Kazakhstan’s transmountain regions using 74 meteorological stations (1981–2023). Data were analyzed using the Mann–Kendall test and Sen’s slope estimator, stratified across six elevation zones from lowlands (<400 m) to high mountains (>1500 m). Results reveal a robust, spatially coherent warming signal across all zones. Annual Tmean increased at a median rate of ~0.30 °C decade⁻¹, peaking at 0.36 °C decade⁻¹ above 1500 m, corresponding to an absolute increase exceeding 1.5 °C. Warming exhibited strong seasonal and diurnal asymmetries. Spring warmed most rapidly, with Tmean increasing >0.60 °C decade⁻¹ (approaching 3 °C total). Winter warming was driven by Tmin increases (up to 0.44 °C decade⁻¹), causing widespread DTR contraction, whereas summer warming was driven by Tmax increases, expanding DTR at higher elevations. Tmin showed the strongest elevation amplification overall. In stark contrast, precipitation trends were weak, spatially heterogeneous, and largely non-significant. Annual changes ranged from −6.63 to +14.35 mm decade⁻¹, with seasonal tendencies indicating modest, non-significant winter/spring wetting and summer drying. Ultimately, the results demonstrate a profound decoupling between strong, elevation-dependent warming and weak precipitation changes. The acute amplification of temperature, particularly during spring and summer at high elevations, has severe implications for snowmelt timing, glacier mass balance, evapotranspiration demand, and long-term water security in Kazakhstan.

Released April 28, 2026 08:33 EST

2026, Preprint

Kelly Filer Robinson, Sarah Nelson Sells, Conor P. McGowan, Elise R. Irwin

Communication regarding the mission of the U.S. Geological Survey (USGS) Cooperative Research Units Program (CRU) can take many forms, yet clear and concise messaging for various audiences is critical to highlight program accomplishments and increase visibility. Before the work described in this report, CRU did not have a communication strategy; therefore, CRU leadership supported a structured decision-making (SDM) workshop to develop a comprehensive strategy for multiple audiences. The workshop was conducted in November 2024, in Nebraska City, Nebraska. The working group for this SDM process included CRU Program leadership, the CRU Communications Team lead, Unit scientists, a Unit administrative assistant, a representative of the Wildlife Management Institute (WMI), a member of the USGS Ecosystems Mission Area (EMA), Office of Communications and Publishing (OCAP) team, and the team lead for the CRU Program strategic planning process, as well as three facilitators who were also unit scientists as well as experts in SDM. Over the course of a week, the SDM team followed the PrOACT framework which identified the problem, objectives, alternatives, consequences, and tradeoffs to guide us toward a strategy for implementation of a set of actions for CRU communications. Results of the SDM workshop included the development of a problem statement, an objectives hierarchy, a suite of alternatives that were evaluated using a consequences table and a clear process for assessing tradeoffs among alternative communication actions and strategies. Through the evaluation of consequences of each action or campaign, the team developed both the assessment tool (for the future) and an immediate plan for communication product development and distribution. The consequences table for this problem was meant to be flexible to accommodate changes in CRU thematic priorities and can be easily updated with new objectives, measures, and alternatives. In addition, the weight placed on objectives may change as the Team moves forward; the ranking and scoring system used in the workshop can be easily updated. Overall, the working group identified three different actions or campaigns—Fact Sheets, Who Are We Campaign, and Alumni Campaign—that scored high in the prototype decision framework. However, the tradeoffs analysis indicated that each action(s) performed better on some objectives than others. The working group identified a need to therefore develop an implementation plan that is composed of individual actions that each target different objectives to potentially create a holistic and feasible communications strategy that performs well for all objectives. In addition, the SDM prototype developed a scalable, objective-based framework for effectively communication of the value and accomplishments of the CRU program.

Released April 27, 2026 09:07 EST

2026, Molecular Ecology Resources (26)

Junman Huang, Richard Flamio Jr., Nathan R. Campbell, Aaron J. DeLonay, Amy C. Buhman, Edward J. Heist

Hatchery supplementation is vital for conserving dwindling fish populations. Effective augmentation requires distinguishing hatchery-origin from wild individuals and accurately identifying species, particularly in systems where closely related species coexist. Genetic monitoring is key to quantifying genetic differences, but conventional markers do not distinguish hybrids, especially backcrosses. Misidentifying hybrids in hatchery programs compromises wild gene pools because hatchery broodstock contributes to numerous offspring being released into the wild. Here, we present a workflow for developing and evaluating the Genotyping-in-Thousands by sequencing (GT-seq) single nucleotide polymorphism (SNP) panel for North American river sturgeons (Scaphirhynchus spp.). This panel is designed to detect complex hybrid classes and to determine parent-offspring relationships. Our species identification panel (S-loci) contains 155 SNPs selected for high genetic differentiation (F_ST) between Pallid Sturgeon (S. albus) and Shovelnose Sturgeon (S. platorynchus), and the parentage assignment panel (P-loci) includes 112 SNPs with high heterozygosity within Pallid Sturgeon. Simulation analyses demonstrated that our GT-seq S-loci panel reliably classifies pure species, F1, F2 and backcross hybrids, even with up to 70% missing data. The P-loci panel achieves high-confidence parentage assignment with ≥ 80% typed loci, with performance influenced by the proportion of sampled parents. Overall, the novel Scaphirhynchus GT-seq panel developed in this study represents a robust and efficient tool for detecting hybridisation, assigning parentage and providing critical information for management decisions in ongoing Pallid Sturgeon conservation.

Released April 24, 2026 13:10 EST

2026, Scientific Investigations Report 2025-5047

Michael S. Kohn, M. Alisa Mast, Tara A. Gross

The U.S. Geological Survey (USGS), in cooperation with the Colorado Department of Transportation, developed peak-, mean-, and low-streamflow regional-regression equations for estimating various statistics for natural streamflow in hydrologic regions of central and western Colorado. The peak-streamflow regression equations were developed using data from 418 streamgages, consisting of 15,202 years of record and a mean of approximately 36 years of record per streamgage. The mean- and low-streamflow regional-regression equations were developed using data from 323 streamgages where daily streamflow data were collected year-round. The annual exceedance-probability discharges for each streamgage were computed using the USGS software program PeakFQ. Mean monthly and 7-day minimum and maximum streamflows were computed using the USGS software program SWToolbox. Streamflow-duration values were computed using an R script. The regional-regression equations were determined using data for the period of record for a given streamgage through water year 2019. Geographic information systems datasets were used to develop 55 basin and 42 climatic characteristics, which were evaluated as candidate explanatory variables in the regression analysis.

For the peak-streamflow regional-regression equations, the study area was divided into four hydrologic regions based on mean basin elevation, including the Plateau (less than 8,014 feet), Mid-Elevation (8,015 feet to 9,492 feet), Sub-Alpine (9,493 feet to 10,490 feet), and Alpine (greater than 10,490 feet) regions. For the peak-streamflow equations, the selection of basin and climatic characteristics was based on the 1-percent annual exceedance-probability discharge for each hydrologic region.

For the mean streamflow, streamflow-duration values, and 7-day minimum and maximum streamflows, the study area was divided into four hydrologic regions based on river basin, including the (1) Colorado-East Slope Headwaters, (2) Green River, (3) Rio Grande, and (4) San Juan-Dolores. For mean streamflows, basin and climatic characteristics were evaluated separately for the annual period and each month for each hydrologic region. Regional regression equations published in this report are available for use in the USGS web-based program StreamStats.

Released April 24, 2026 11:06 EST

2026, Book chapter, The effects of noise on aquatic life IV

Daniel P. Costa, Rachel R. Holser, Garrett T. Shipway, Arina B. Favilla, Birgitte I. McDonald, Daphne M. Shen, Amber R. Diluzio, Sarah H. Peterson, Joshua T. Ackerman, Daniel E. Crocker

Arthur N. Popper, Joseph A. Sisneros, Paul A. Lepper, Kathleen J. Vigness-Raposa, editor(s)

Noise exposure is a potential stressor for free-ranging marine mammals and is often studied in the absence of other environmental factors. Here, a multi-investigator, interdisciplinary effort was undertaken to examine the response of elephant seals to multiple stressors. An integrated physiological and ecological approach was taken, including immunology, stress physiology, toxicology, animal behavior, population biology, and life history theory, to examine the cumulative effects of exposure to multiple stressors in elephant seals. While we measured the response of individual animals, a population response can be predicted by incorporating these results into the long-term data on elephant seal demographics.

Released April 23, 2026 11:45 EST

2026, Fact Sheet 2026-3061

Brian N. Shaffer, Elisa Alonso, Michelle N. Johnston, Scott A. Kinney

Plain Language Summary

Metallurgical coal (met coal; consumed to produce coke for steelmaking) must meet specific chemical and physical specifications. In 2023, the conterminous United States produced 66 million short tons (mst) of met coal, consumed 15.85 mst domestically, exported 51.1 mst, and imported 0.7 mst. Most met coal was produced in the Appalachian Basin, but there are also resources that meet the specifications for met coal in the Western United States.

Released April 23, 2026 11:45 EST

2026, Fact Sheet 2025-3057

Mark J. Mihalasky

Introduction

Interest in nuclear power for the generation of electricity has risen with the increase in the need for more diverse baseload power, enhanced energy security, and the development of new technologies, such as small modular reactors (SMRs), which could provide power for remote areas, industrial applications, and artificial intelligence (AI) data centers. In 2024, the U.S. Department of Energy received $2.7 billion in congressional funding to bolster the domestic uranium production and nuclear fuel supply chain and address reliance on imports from foreign suppliers. In 2025, the U.S. Government issued several Executive and Secretary’s orders aimed at revitalizing the U.S. nuclear sector. If SMRs are to be as widely deployed in the United States and worldwide as envisioned, demand for uranium (nuclear reactor fuel) will likely increase.

After the Fukushima nuclear accident in 2011, the market spot price of uranium began a decline, followed by a decrease in U.S. and global uranium exploration and mine development expenditures that led to a uranium supply deficit until 2020, when prices started to recover, prompting a resurgence in uranium exploration and development. In January of 2024, the uranium spot price rose to a 17-year high $106 (U.S. dollars) per pound of U3O8 (triuranium oxide, commonly known as “yellowcake”), which is expected to increase uranium exploration, mine development, and uranium production domestically and worldwide.

Released April 23, 2026 10:15 EST

2026, Environmental Epidemiology (10)

Amy M. Lavery, Jordan Murray, Audrey F. Pennington, Blake Schaeffer, Bridget Seegers, Elizabeth D. Hilborn, Keith Loftin, Stephen Scroggins, Lorraine Backer

Background: 

Cyanobacterial harmful algal blooms (cyanoHABs) pose risks to human and animal health.

Methods: 

We investigated the relationship between cyanoHABs and asthma or wheeze-related emergency department (ED) visits near three Wisconsin cities (Green Bay, Madison, and Oshkosh) during 2017–2019. CyanoHAB exposure was approximated using the Cyanobacterial Assessment Network remotely sensed satellite indicator of cyanobacterial biomass, a chlorophyl algorithm (Chl_BS) aggregated by water-adjacent ZIP Code Tabulation Areas (ZCTA), and distance weighted from the nearest waterbody. Weekly counts of ED visits for asthma or wheeze were aggregated by ZCTA. Poisson generalized linear models estimated the association between the weekly number of ED visits and weekly Chl_BS, adjusting for maximum temperature, dewpoint, fine particulate matter (PM_2.5), month, and correlation within ZCTA.

Results: 

During 2017–2019, 7,057 ED visits for asthma or wheeze occurred in the study area (42 ZCTAs). Peaks in Chl_BS occurred between June and October, with higher values in Lake Winnebago and Lake Mendota compared to Green Bay. Chl_BS was not associated with ED visits for asthma or wheeze (adjusted rate ratio = 1.00, 95% confidence interval = 0.99, 1.00), and the presence of onshore winds did not change this result. Monthly aggregations of ED visits and Chl_BS showed a monotonic trend between increasing Chl_BS and ED visits during July–September.

Conclusion: 

This study demonstrates the utility of remote sensing data in environmental health research. Future studies could explore individual-level exposure and outcomes to refine health risks associated with cyanoHABs.

Released April 22, 2026 14:45 EST

2026, Open-File Report 2026-1004

John M. Plumb, Russell W. Perry

Restoration of salmon populations in the upper Lewis River Basin, Washington, depends on a trap-and-haul program owing to the Lewis River Hydroelectric Project (hereinafter referred to as “Project”) operated by PacifiCorp and Cowlitz Public Utilities District (hereinafter referred to as “Utilities”), which has been a barrier to salmon passage since the 1930s. Thus, sustaining the Oncorhynchus kisutch (Walbaum, 1792; coho salmon) population upstream from the Project currently depends on two fundamental factors: (1) the collection of upstream migrating adult coho salmon at Merwin Dam, the lowermost dam within the Project, and transporting them by truck to spawn above Swift Dam, the uppermost dam within the Project; and (2) the collection of out-migrating juvenile coho salmon at the downstream collection facility at Swift Dam for transport and release below the Project. The reintroduction program began once the downstream collection facility at Swift Dam was commissioned in late 2012, with the first year of transport data being collected in 2013. Over the past decade, the Utilities have been collecting data on juvenile outmigrants and adult fish returns at the dams. The need to construct a lifecycle model for Lewis River anadromous fish was identified by the Lewis River Aquatic Technical Subgroup, with the understanding that many years (more than 15 years) of data collection are needed to adequately measure the lifecycle production of salmon. The U.S. Geological Survey was contracted to develop and apply the model to past data at the Lewis River dams to help inform future data collection and provide a framework that can be updated annually to measure trap-and-haul program performance within a lifecycle context.

Because coho salmon can live as long as 5 years, estimating demographic parameters for coho salmon populations over their lifecycle requires at least 10 or more years of data collection. Over the past decade, PacifiCorp has been collecting data on fish collection efficiency and the numbers of adult and juvenile salmon transported around the Lewis River dams, making this an ideal time to formulate a lifecycle model that can guide future data collection efforts and provide preliminary information to resource managers. The goal of the statistical lifecycle model is to estimate annual production and survival during two critical life-stage transitions: (1) the freshwater production from escapement of adults released upstream from Swift Dam, and the collection of downstream migrating juveniles at the downstream passage facility at Swift Dam; and (2) the smolt-to-adult survival from the time of collection at Swift Dam to their return as adults. We used the Beverton-Holt stock-recruitment model to estimate juvenile production from the number of spawners (Beverton and Holt, 1957). This approach allowed us to test for density dependence at current spawner abundances while estimating annual productivity, defined as the number of juveniles produced per spawner at low spawner abundance. Productivity was then expressed as a function of the number of juveniles collected and transported downstream from the Project. Because juvenile fish collection efficiency (FCE) directly affects the number of juveniles that survive to continue downstream migration, FCE is a primary determinant of fish production. Consequently, the modeling framework is well suited to evaluate the performance of trap-and-haul programs within a lifecycle context.

The objectives of this study were to (1) gather and collate available data on adult and juvenile coho salmon at Merwin and Swift Dams; (2) quantify adult escapement, juvenile abundance, and the age at outmigration and adult return; (3) describe, formulate and fit the integrated population model to the data; and (4) summarize our findings, identify data gaps, and identify opportunities for future studies that could improve model estimation and inference. Our key findings were: (1) over and above the number of spawning females, FCE was the primary factor affecting productivity of coho salmon above Swift Dam; (2) smolt-to-adult return (SAR) rates were relatively high considering that harvest was included in the estimate, averaging about 4.5 percent and ranging as high as 12.9 percent; and (3) juvenile capacity upstream from Swift Dam was difficult to estimate due to the limited range in spawning females over the time series of data, suggesting the model may be improved by collecting data at higher spawner abundances. In addition, by including FCE in the model, we estimated that the median pre-collection productivity, defined as the number of juveniles produced per spawner when FCE=1, was 64 juveniles per spawner. Because the two-stage lifecycle model partitions factors that affect fish production in rivers versus the ocean, the model estimates may help inform fishery managers about the overall role that fish collection at Swift Dam plays in the recovery and sustainability of Lewis River coho salmon. By providing the model with (1) more years of data, (2) higher numbers of spawning females, and (3) data on age at juvenile migration in relation to age at adult return, greater certainty in the estimates of capacity and SAR can be attained. Ultimately, information provided by the model may assist in the evaluation and continued improvement of the current trap-and-haul program to support anadromous fishes in the Lewis River Basin.

Released April 22, 2026 11:50 EST

2026, Fact Sheet 2026-3005

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

 Using a geology-based assessment methodology, the U.S. Geological Survey estimated undiscovered, technically recoverable mean resources of 80.1 trillion cubic feet of shale gas in the Grand Erg/Ahnet Basin Province of Algeria.

Released April 22, 2026 10:00 EST

2026, Computers and Electronics in Agriculture (248)

Robert M.H. Sterling, Elsa Marie-Catherine Goerig, M Buzdalov, Theodore Castro-Santos, O. Akanyeti

Body and caudal fin locomotion is ubiquitous in aquatic vertebrates, and kinematic models describing it are used in robotics, biomechanics and fisheries research. This paper presents a new algorithm to translate continuous body midlines of fish into a series of interconnected segments by identifying favorable joint positions along the body. The algorithm employs binary search to generate parsimonious kinematic models, aiming at minimizing the number of segments yet keeping approximation error below a user-defined threshold. To achieve this, the algorithm maximizes the length of each segment by determining the most distal joint position through repetitive shrinking of the search space. Theoretical and empirical analysis using two different datasets show that the binary search algorithm is substantially faster when compared to segment growing algorithm, which employs linear search to generate its models. There is four-fold improvement in computation time when generating models with less than 10 segments, which are typically sufficient to describe fish and fish-inspired robot movements. Furthermore, the multi-segment models generated by the binary search algorithm matched the ground truth models obtained through dynamic programming in over 97% of cases, and on average, contained one fewer segment than those produced by the Ramer–Douglas–Peucker algorithm, which is widely used in curvature simplification tasks. Our findings suggest that the binary search algorithm provides a computationally efficient approach for generating compact kinematic models and may facilitate the analysis of large datasets with high temporal and spatial resolution.

Released April 22, 2026 09:57 EST

2026, Journal of American Water Resources Association (62)

Brock J.W. Kamrath, Jennifer C. Murphy, Lindsey Ayn Schafer, Hannah Lee Podzorski, Gregory F. McIsaac

Illinois contributes substantial nutrient loads to the Gulf of America, warranting watershed-scale assessment. This study estimated nitrate-nitrogen (nitrate-N) and total phosphorus (TP) loads and yields for 49 Illinois 8-digit hydrologic unit code (HUC8) watersheds draining to the Mississippi River Basin from 1997–2022, comparing recent (2018–2022) to baseline (1997–2011) conditions. Estimates included point and nonpoint source contributions, dissolved phosphorus, and water yields. During the recent period, nonpoint sources dominated nutrient export (82% nitrate-N, 78% TP), though point sources drove high yields in the Chicago area. Spatially, nonpoint source nutrient hotspots persisted with nitrate-N yields highest in east-central and northern Illinois and TP yields higher in southern and western Illinois. Temporally, statewide nitrate-N loads decreased 9%, while TP loads increased 27%. Nitrate-N yields increased in 22 HUC8s and decreased in 20, while TP yields increased in 32 HUC8s and decreased in 9. For both nutrients, baseline yields were negatively correlated with yield changes, indicating high-yielding watersheds tended toward larger decreases or smaller increases. Water yields increased 19% on average but were weakly correlated with nutrient yield changes (r = 0.23 and 0.20 for nitrate-N and TP). These results reveal spatially persistent yet temporally divergent nutrient export across Illinois, with contrasting nitrate-N and TP trajectories for nonpoint sources.

Released April 22, 2026 09:50 EST

2026, Geophysical Research Letters (53)

Dongxiao Yin, Zafer Defne, Neil K. Ganju, John C. Warner, David K. Ralston, Courtney K. Harris, Bin Li

Salt marsh vulnerability to sea-level rise (SLR) is typically assessed using point measurements of vertical accretion, neglecting three-dimensionality of geomorphic evolution and spatial variability. Recent studies suggest links between vertical and horizontal vulnerability, with differences between oligohaline and polyhaline marshes, yet these relationships remain untested in estuary-marsh systems. Here we combine geospatial analysis with hydrodynamic modeling to evaluate how unvegetated/vegetated marsh ratio (UVVR), a metric of marsh degradation, relates to elevation across hydrological regions and salinity regimes in the Albemarle-Pamlico Estuarine System, the largest lagoonal estuary in U.S. We show that at given normalized elevation, UVVR decreases across hydrological regions and salinity regimes from offshore to inland. UVVR-elevation relationship varies systematically with both hydrological setting and salinity regime, with hydrology exerting stronger influence. These findings challenge the assumption of a universal marsh deterioration trajectory and underscore the need to account for spatial heterogeneity when predicting responses to SLR.

Released April 22, 2026 09:38 EST

2026, Earth and Planetary Science Letters (686)

Xudong Zhao, Yifei Li, Huiping Zhang, Richard O. Lease, Ying Wang, Yuqi Hao, Zifa Ma, Hao Xie, Huan Kang, Jianguo Xiong, Peizhen Zhang

The modern drainage network of eastern Tibet is widely believed to have developed through a series of river capture and flow reversal events; however, the timing and mechanisms driving this reorganization remain contentious. Among these events, the river capture that formed the First Bend of the Yangtze River (YFB) stands out as both iconic and particularly debated. Here we present sedimentary provenance data from the Late Miocene–Quaternary Dali Basin, located south of the YFB, which indicate that a southward-flowing Jinsha River (i.e., the present-day upper Yangtze River) sourced sediment to the Dali basin at ∼7.4–6.4 Ma in a drainage configuration different from that of today. Because this interval postdates the initial establishment of a near-modern Jinsha River system prior to the Miocene, our results imply at least two discrete fluvial reorganizations occurred at the YFB—one preceding ∼7.4 Ma and another following ∼6.4 Ma. By integrating these findings with landscape evolution modeling, we infer that the initiation of rapid uplift of the Yulong-Haba Mountains and the Diancang Shan may have been responsible for these drainage reorganizations. These results underscore that Cenozoic drainage systems on the eastern Tibetan Plateau have evolved dynamically on a short timescale of ∼10⁵–10⁶-year, rather than remaining in a long-term stationary configuration on ∼10⁷-year timescales.

Released April 22, 2026 07:55 EST

2026, Canadian Journal of Earth Sciences (63) 1-15

Benjamin J. Drenth, Amy L. Radakovich, Amanda Kate Souders, George J. Hudak, Stacy Saari

The Neoarchean Mentor anorthosite intrusive complex (MAIC) lies within the Wawa subprovince in northwestern Minnesota, in a region where the Wawa, Quetico, and Wabigoon subprovinces are juxtaposed in close proximity. Archean rocks are entirely concealed, and interpretations are developed from aeromagnetic, gravity, and borehole samples. The MAIC includes both anorthosite (dense, weakly magnetized) and oxide-rich gabbro (strongly magnetized) lithologies. Anorthosite is coarse-grained to megacrystic, intensely altered, and locally brecciated. Pervasive epidote alteration enhanced the density of the anorthosite via introduction of higher-density mineral assemblages, explaining why the MAIC produces a significant gravity high. Oxide-rich gabbro forms a border phase of the MAIC and has potential for vanadium, chromium, and titanium mineralization, and produces a strong aeromagnetic high. The MAIC is interpreted to extend over an area of 640 km², making it the largest known anorthosite complex of the Superior Province, as measured by preserved areal extent. Modeling indicates the MAIC extends more than 6 km into the subsurface. A new Pb–Pb zircon age of 2737.2 ± 4.5 Ma is interpretated as the crystallization age of anorthosite within the MAIC, showing that the MAIC formed well before the ca. 2690 Ma Shebandowanian orogeny, and raising new questions about correlations with other parts of the Wawa subprovince. A low-density batholith, here informally called the Fertile batholith, is interpreted to intrude the southern part of the MAIC. A new Pb–Pb zircon age of 2701.1 ± 6 Ma is interpreted as the magmatic age of the Fertile batholith.

Released April 22, 2026 07:47 EST

2026, Ecological Informatics (95)

Janet S. Prevéy, Cameron J. Reimer, Peder S. Engelstad, Pairsa N. Belamaric, Terri Hogan, Jillian M. LaRoe, Colter J. Mumford, Jennifer L. Sieracki, Catherine S. Jarnevich

Resource managers are tasked with protecting natural areas from invasive species with limited resources. Further, invasive management goals can vary greatly based on different management priorities specific to management agencies or taxa of interest. The site prioritization tool for invasive species management addresses these challenges by creating a platform to view and combine diverse spatial data layers to estimate cumulative invasion risk based on user-specific needs. For this tool, we developed a human transport risk layer, estimating invasion risk based on proximity to human population centers and transportation corridors, and created maps of non-native species richness across the conterminous United States. The tool also includes spatial layers showing projected changes in key climate variables through the end of the century to identify areas where invasion risk may shift. Users can explore these layers to prioritize sites based on the invasive taxa of interest, likely invasion pathways, and disturbances that may elevate invasion risk. This interactive tool will allow managers to make the spatial comparisons needed to focus efforts on areas that are highly susceptible to invasion and efficiently target monitoring and suppression efforts.

Released April 21, 2026 10:11 EST

2026, Frontiers in Ecology and Evolution (14)

Barbara E. Kus, Julie L. Yee, Shannon M. Mendia

A serious emerging threat to southern California riparian ecosystems is the invasive shot hole borer (Euwallacea spp.; SHB), a non-native beetle that cultivates a pathogenic fungus that kills trees of 66 reproductive host species. We examined the response of the bird community at the Tijuana River, California, to a massive SHB infestation in 2015 using data from a Monitoring Avian Productivity and Survivorship (MAPS) station operated during 7 pre-infestation (2009-15) and 7 post-infestation (2017-23) years. Species richness did not change between pre- and immediate (2017-18) post-SHB periods, but average annual adult captures declined by 27%. Among the species making up ≥ 5% of the total individuals caught in any one year (n=15), abundance declined by up to 76% in 10 species, including those most abundant at the station (Bushtit (Psaltriparus minimus), Song Sparrow (Melospiza melodia), Common Yellowthroat (Geothlypis trichas), Orange-crowned Warbler (Leiothlypis celata), and Wilson’s Warbler (Cardellina pusilla)). Mean annual abundance increased slightly for the endangered Least Bell’s Vireo (Vireo bellii pusillus) and Northern Yellow Warbler (Setophaga aestiva) and doubled for House Finch (Haemorhous mexicanus) and Western Warbling-Vireo (V. swainsoni). We compared species trends at the Tijuana River to those at a nearby uninfested MAPS station on the Santa Margarita River to isolate the effect of SHB from other factors influencing annual abundance. The contribution of SHB to changes in abundance post-SHB was high (63-80%) for 7 declining species, moderate (22-45%) for 4 species, and weakly to moderately positive (18-40%) for 3 species. By 2019, the SHB infestation at the Tijuana River had abated and canopy cover was recovering through resprouting of mature willows (Salix spp.) and seedling establishment. Bird abundance tracked this regrowth, with all of the species strongly affected by SHB increasing between 2019-23. The rapid recovery of the Tijuana River habitat and the associated response by the bird community are encouraging signs that the threat of the invasive shot hole borer to regional biodiversity may not be as great as originally anticipated.

Released April 21, 2026 09:27 EST

2026, Seismological Research Letters

Nadine G. Reitman, Yuanshi Wang, Yu-Ting Kuo, Catherine Elise Hanagan, Alexandra Elise Hatem, Christopher B. DuRoss, Chun-Chi Chen, Dara Elyse Goldberg, Harriet Zoe Yin, Richard W. Briggs, Jessica A. Thompson Jobe, Sylvia R. Nicovich, Emerson Madelyn Lynch, Joseph Hoss Powell, William D. Barnhart, Robert G. Schmitt

The 2025 M_w 7.7 Mandalay, Burma (Myanmar), earthquake ruptured 475 km of the central Sagaing fault and is the longest continental strike-slip rupture on record. The observed rupture length is 1.6–4.7 times the value expected (100–300 km) from existing length-magnitude scaling relations for strike-slip earthquakes. The earthquake resulted from shallow dextral faulting and ruptured bilaterally with supershear speeds south of the epicenter, rupturing close to three major cities in Myanmar and exposing over six million people to violent or extreme shaking. We report on the surface rupture character, length, and slip distribution based on sub-pixel correlation of Sentinel-2 (10 m) and Planet Dove (3 m) optical images and visual analysis of SkySat and WorldView (0.3–0.5 m) optical images. The earthquake had moderate surface slip (average = 3.3 m, maximum = 5.6 m, 25–75% range = 3.0–4.0 m), narrow deformation zone width (1–10 pixels in sub-pixel correlation and up to 190 meters for the detailed surface rupture mapping), and simple fault geometry (no stepovers or large changes in strike, 87% of the rupture that was mapped in detail is single-stranded). We attribute the extreme length of the Mandalay earthquake to supershear rupture speed, simple fault geometry, narrow down-dip width, and moderate surface slip. Based on a compilation of 25 supershear strike-slip earthquakes (M_w 6.5–8.6; 1979–2025), we find that the rupture length of supershear earthquakes does not fit empirical scaling relationships for strike-slip earthquakes that predict length from magnitude. A length-magnitude scaling relationship based on supershear earthquakes has a best fit of log₁₀(surface rupture length) = 0.89 M_w – 4.44, indicating that supershear earthquakes tend to be longer than their subshear counterparts for any given magnitude and thus may expose a greater population to shaking.

Released April 21, 2026 09:10 EST

2026, Techniques and Methods 5-E1

Haley M. Dietz

The U.S. Geological Survey’s (USGS) Lithologic Core Storage Library (CSL) at the Idaho National Laboratory stores more than 120,000 feet of drill core that is accessible to the public for research and sampling. To effectively convey the physical and descriptive properties of the drill core, USGS staff at the Idaho National Laboratory Project Office log the drill core and publish the lithologic logs as data releases. The logs provide essential data on the lithology, texture, mineralogy, alteration, and other physical properties of the core, which serve as valuable information for researchers to guide their research and sampling efforts. To ensure consistent, quality, and dependable lithologic logs, this document outlines the procedures and expectations for logging drill core at the CSL. This document describes the processes for storing, photographing, and logging core, and includes a variety of resources, reference materials, and appendixes designed to standardize and aid the logging process. Following the procedures outlined in this document will produce consistent, detailed logs that facilitate dependable observations and serve as an easy reference for researchers and other interested parties.

Released April 20, 2026 14:45 EST

2026, Scientific Investigations Report 2026-5136

Bridget B. Kaemming, Chanse M. Ford, Sherry L. Martin

Lake Huron, one of the five Great Lakes, borders the United States and Canada, with Michigan as the only U.S. State on its shoreline. Like other freshwater lakes, it faces water-quality challenges from nutrients and chemicals applied across its drainage basin. Although past studies focused on surface-water sources, groundwater contributions remain less understood. To address this gap, the U.S. Geological Survey, as part of the Cooperative Science and Monitoring Initiative, classified drainage basins to Lake Huron into eight hydrogeologic zones based on bedrock rock type and glacial sediment transmissivity. Utilizing existing data and empirical field data, we quantified groundwater discharge and identified areas of concern for loading of chloride and nitrate to Lake Huron. Groundwater contributions, including indirect and shoreline discharge, ranged from 5.8 to 11.5 inches annually, totaling 1.9 cubic miles and 0.09 cubic mile, respectively. Hydrogeologic zones with higher glacial sediment transmissivity yielded greater indirect groundwater discharge. Chloride levels above the U.S. Environmental Protection Agency’s 250-mg/L recommendation were mainly in the Saginaw lowlands, whereas nitrate above the 10-mg/L standard was rare—found in only 11 wells. Together, the analysis of where groundwater discharge is occurring in the Lake Huron Basin and the identification of areas with potential groundwater-quality concerns can help prioritize areas that are critical to protecting the long-term health of Lake Huron.

Released April 20, 2026 11:08 EST

2026, Open-File Report 2026-1006

Sarah Timbie, Shelby Jo Weidenkopf, Daniel A. Grear

The New World screwworm (Cochliomyia hominivorax; NWS) is a parasitic blowfly that lays its eggs in open wounds of live, warm-blooded animals including livestock, wildlife, and potentially humans. The larvae consume living animal tissue, and if untreated, the infestation can lead to death. Although NWS was eradicated in the United States in 1966, it has been moving northward from its endemic range in South America during the past decade and could seriously threaten the health of U.S. wildlife populations, making detection, treatment, and surveillance of the disease far more difficult across this multi-sector disease system.
As the likelihood of NWS reintroduction to the United States increases, veterinarians, wildlife managers, and conservation specialists need to be informed and prepared to respond. The existing knowledge about NWS interactions with wildlife hosts is lacking, especially regarding North American species where the NWS has been eradicated for more than 50 years. To address this knowledge gap, we compiled an annotated bibliography that consolidates key information from the existing literature on NWS infestation in wild animals.

Released April 20, 2026 09:06 EST

2026, Hydrological Processes (40)

J. R. Meyer, A. L. Mianecki, E. Occhi, Dana W. Kolpin, G. H. LeFevre

Contaminants found in treated wastewater discharged to streams, including pharmaceuticals and per- and polyfluoroalkyl substances (PFAS), are of global concern due to their deleterious effects on aquatic ecosystems and potential impacts to human health. Hyporheic zones have strong potential for contaminant attenuation. Assessing this potential requires collection of physical and biogeochemical data within the hyporheic zone. This study tested the applicability of a novel drive-point multilevel system (DP-MLS) for quantifying head profiles and characterizing contaminant concentrations in the hyporheic zone of a temperate region effluent dominated stream (EDS). DP-MLS, each with 4 ports, were installed in the stream bed at two sites, DS-1 and DS-2, 0.2 and 4.7 km downstream of the effluent outfall, respectively. Head profiles were measured and groundwater collected for analysis of pharmaceuticals and PFAS temporally over two years. The DP-MLS withstood rapid changes in stage, ice formation, and floating debris. Vertical hydraulic gradients (VHG) were generally upward but varied in magnitude indicating heterogeneity in hydraulic conductivity and variability in flow conditions. Upward VHG were also about 2X larger at DS-1 than at DS-2. Contaminant concentration profiles consistently showed penetration of pharmaceuticals and PFAS to 1 m below the bed at DS-2 while there was less penetration, lower groundwater concentrations, and more temporal variability in concentrations at DS-1. Integration of the physical and chemical data suggests weaker upwelling conditions at DS-2 are more easily reversed during periods of high stream stage, which could facilitate migration of wastewater contaminants into the bed. However, further studies incorporating other transport processes and reach scale dynamics are required to fully characterize these exchanges. Overall, this study demonstrates the efficacy of these novel DP-MLSs for characterization of the hyporheic zone and provides new insights into the occurrence, composition, and persistence of wastewater derived contaminants in the hyporheic zone of a well-studied EDS.

Released April 18, 2026 10:21 EST

2026, Natural Resources Research

Niki E. Wintzer, Christopher S. Holm-Denoma, Jacob Evan Poletti, Dalton M. McCaffrey, Stanley Paul Mordensky, Erik Roger Tharalson, Collin Cronkite-Ratcliff

Lithium demand is projected to increase more than 48 times by 2040 due to electric vehicle production and other energy storage needs. Most lithium production is outside of the USA, thereby increasing supply chain vulnerability. The combined end use importance and heightened supply risk of lithium make this lightest metallic element a critical commodity to the USA. To mitigate this supply risk, the US Geological Survey is actively assessing lithium deposits in the USA. Herein, we detail an assessment for lithium-mineralized pegmatites in the US northern Appalachian Mountains. Permissive tracts were generated by cross-referencing tectonic and geologic maps and mineral occurrence data with mappable criteria derived from generalized and region-specific lithium pegmatite ore deposit models; tracts were then ranked as having high, medium, or low permissibility. Available geophysical and geochemical data were found to be of minimal utility for this deposit type at the scale of the assessment. The number of undiscovered deposits were estimated and integrated into probabilistic simulations, which included an expanded and updated global grade and tonnage model of pegmatite-hosted lithium ore. The estimated total amount of undiscovered resources for the northern Appalachian Orogen has a median value of 1,410,000 metric tons of Li₂O when considering moderate correlation across sub-regions. At a confidence level of 90%, a resource of at least 90,000 metric tons of Li₂O remains undiscovered, and at a 10% confidence level, a resource of as much as 7,380,000 metric tons Li₂O remains undiscovered. After applying an up-to-date economic filter to convert median contained lithium to recoverable material, a correlated total of 900,000 metric tons of Li₂O may be economically extractable, equating to enough Li₂O to provide the current annual US lithium supply deficit (presently obtained through net imports) for 127 years at 2025 rates of apparent consumption. This period of provision will inevitably shorten with projected increasing consumption rates, emphasizing that further research could be completed to better delineate regions of high lithium resource potential and support exploration and domestic production.

Released April 18, 2026 08:04 EST

2026, Geophysical Research Letters (53)

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

Freshwater salinization impacts the availability of water for human use and ecosystem needs worldwide. It has been estimated that total dissolved solids (TDS) in the Colorado River Basin cause $350 million/year in damages and substantial resources are devoted to reducing TDS loading to streams. This study describes the development and application of coupled watershed models that enable TDS source tracking through the subsurface and across the landscape at a seasonal timestep for 35 years in the Upper Colorado River Basin. Results indicate that, on average, 75% of TDS loading to streams originates as baseflow, and 50% of loading is lagged in delivery by longer than one season. Snowmelt was identified as a dominant process controlling the transport of lagged TDS to streams. This approach informs when and where TDS mitigation efforts may be effective in a watershed that serves as a critical water supply for the southwestern United States.

Released April 17, 2026 10:04 EST

2026, The Cryosphere (20) 2169-2179

Kajsa Holland-Goon, Randall Ray Bonnell, Daniel McGrath, W. Brad Baxter, Tate Meehan, Ryan Webb, Christopher F. Larsen, Hans-Peter Marshall, Megan A. Mason, Carrie Vuyovich

Snow is a vital component of high-latitude terrestrial systems, but environmental factors (e.g., permafrost) and complex vegetation challenge the accurate measurement of key snowpack properties. We evaluated local-scale ground-penetrating radar (GPR) and large-scale airborne lidar retrievals of snow depth collected during the NASA SnowEx 2023 campaign in tundra and boreal forest environments in Alaska along 44 short (3–12 m) transects. Compared to in situ observations, we identified modest biases for GPR snow depths (bias <0.03 m in tundra, +0.06 m in boreal forests) and larger biases for lidar snow depths in the boreal forests (–0.16 m). At the Upper Kuparuk-Toolik tundra site, lidar snow depths exhibited a small bias (–0.02 m), whereas the bias was much larger at the Arctic Coastal Plain tundra site (+0.19 m). For most sites, biases were primarily related to sub-snow vegetation, tussocks, and seasonally dynamic ground. However, we identified vertical alignment issues with the Arctic Coastal Plain lidar snow depth dataset that likely contributed to the higher bias. The complex ground surface and sub-snow vegetation in these environments present a challenge to established snow depth measurement methods, which needs to be considered when evaluating novel remote sensing approaches.

Released April 17, 2026 09:11 EST

2026, Ecological Informatics (95)

Song S. Qian, Mark Richard Dufour, Sabrina Jaffe, Corbin David Hilling, William D. Hintz

Species diversity is the foundation of many ecological disciplines. This metric is often approximated using species richness and evenness, even though actual richness likely exceeds observations due to imperfect sampling methods. Estimating the “true” species richness, which includes identifying the number of missing species, has intrigued ecologists for decades. We adopted a parametric model that appeared in Fisher et al. (1943), which models the numbers of individuals from different species as random samples from a negative binomial distribution, and developed a Bayesian computational approach to directly estimate the distribution model parameters. The model parameters represent species abundance and evenness, and can be used to derive species richness. We evaluated our parametric approach using (1) a simulation study and (2) three historical data sets. Furthermore, we illustrated the hierarchical modeling approach to combine data from multiple parallel studies using a biannual fishery survey data set. Our parametric model formulation is computationally efficient, and the hierarchical structure facilitates embedding diversity estimation into broader application, such as assessing spatial and temporal trends in species diversity associated with environmental stressors. Additionally, because the two parameters of the negative binomial distribution model represent species abundance and evenness of a community, this parametric approach facilitates a deeper understanding of the ecological systems under study. The negative binomial distribution model works with a wide range of species frequency distribution types. As a result, our emphasis on a parametric model can help us characterize the structure of an ecosystem and provide a greater depth of ecologically meaningful information.

Released April 17, 2026 08:57 EST

2026, PLOS Water (5)

Janice M. Gillespie, Riley Gannon, Lyndsay B. Ball, John G. Warden, Rhett R. Everett, Michael J. Stephens

This study seeks to determine the effects of oil field produced water disposal operations and well mechanical integrity issues on groundwater quality in oil fields in the southwest San Joaquin Valley, California. Whereas previous studies used groundwater wells to study shallow aquifers outside the oil fields, this study demonstrates that future approaches may use oil well geophysical logs to map groundwater head gradients, create salinity profiles and document changes in salinity over time in oil field areas with sparse groundwater well data and at depths greater than 330 m. We also incorporate an analysis of well histories to determine potential effects of compromised wellbore seals on changes in aquifer quality that cannot be explained by water disposal practices. Water quality in the aquifers is naturally brackish across most of the area, with better quality groundwater occurring in the eastern part. Geophysical logs are used to determine salinity variations within aquifers including the depth at which TDS exceeds 10,000 mg/L. This depth ranges from 366 m in the northwest to approximately 1,500 m in the southeast. Oil well porosity logs are used to determine water table elevations. These logs indicate the water table slopes south-southeast, showing the predominant groundwater flow direction is from oil field disposal areas toward better quality groundwater east of the oil fields. Geophysical logs show formation resistivity near some disposal facilities has decreased over time, indicating the salinity of the aquifer has increased due to disposal of saline produced water in injection wells and ponds. Oil well history analysis suggests that increased salinity over time in water-saturated sand intervals >1.5 km from disposal facilities may be caused by mechanical failures and/or incomplete borehole seals in poorly constructed or abandoned wellbores prevalent throughout the study area—particularly wells drilled prior to 1930.

Released April 17, 2026 07:49 EST

2026, Journal of Volcanology and Geothermal Research (474)

Shaul Hurwitz, R. Blaine McCleskey, Maxwell L. Rudolph, Sara Peek, David A. Roth, Melissa Schott-Atkins, Michael Manga, Kiernan F. Folz Donahue, Mara H. Reed, Jefferson D.G. Hungerford

The iconic Old Faithful Geyser in Yellowstone National Park, USA, has attracted a significant amount of research because of the relative regularity and impressive size of its eruptions. Numerous studies have included observations, measurements, and analyses that informed models of geyser eruptions. However, fundamental quantities, including the associated mass and heat discharged, remain poorly constrained. In April 2025 we measured the volume of water from 45 Old Faithful Geyser eruptions using a portable flume in an outflow channel and specific conductance measurements in the Firehole River. We used high-speed video to perform velocimetry, measured changes in water chemistry to calculate the volume of water evaporated along the outflow channels, and used thermodynamic calculations to estimate the volume of water erupted as steam and to quantify the geyser's heat output. The calculated average volume of water discharged by Old Faithful Geyser in each eruption is 27.9 ± 9.4 m³, with no relation between eruption volume and the length of the preceding eruption interval. Video analysis of the eruption's liquid-dominated phase yields similar volumes of 21–30 m³. The calculated heat flow from the geyser is 2.2–2.4 MW and the average annual discharge of chloride, fluoride, and arsenic are 63 tons, 3.9 tons, and 241 kg, respectively. Average annual silica deposition rate on the geyser cone and along the outflow channels is 7 tons. This study provides a methodology for future studies at geysers worldwide and a baseline for monitoring future activity changes at Old Faithful.

Released April 16, 2026 16:18 EST

2026, Fact Sheet 2026-3064

Sean Patrick Gaynor, Eric D. Anderson, Kyle A. Eastman, Karen Lund, Chris Gammons, Heather A. Lowers, Jay M. Thompson

Plain Language Summary 

U.S. Geological Survey research, in collaboration with Montana Technical University and Montana Bureau of Geology and Mines, is providing key critical mineral information that may have potential for critical mineral production of several mining districts in the Boulder Batholith region, to better understand the abundance and distribution of natural resources within this region. Continued research can be used to show the potential for previously undiscovered critical mineral resources in southwestern Montana and in other parts of the United States.

Released April 16, 2026 14:10 EST

2026, Open-File Report 2026-1002

Jason P. Pope, Alison D. Gordon, Ryan S. Frederiks

Computation of detailed groundwater flow budgets for subdivisions of the Virginia Coastal Plain aquifer system has enabled quantification and more thorough understanding of groundwater flow within this important water resource. A zone budget analysis based on previously published groundwater models of the Virginia Coastal Plain and Virginia Eastern Shore indicates that groundwater conditions vary substantially throughout the Coastal Plain aquifer system because of local variations in hydrogeology and historical and ongoing variations in groundwater use and management. Decades of substantial groundwater withdrawal from the Coastal Plain aquifer system have altered groundwater flow from predevelopment conditions. Rates of sustainable withdrawal are limited because the downward groundwater flow rate into confined aquifers is a relatively small part of the total groundwater budget for the aquifer system compared to the rate of recharge at the land surface.

Analyses of groundwater budgets from the Virginia Coastal Plain model indicate that groundwater flow is generally outward from the surficial aquifer to rivers and coastal waterbodies and downward through a series of underlying aquifers and confining units to the Potomac aquifer, which is the deepest aquifer and the source of most groundwater withdrawals. Downward flow into the Potomac aquifer is estimated to be only 7 percent of total net precipitation-derived net recharge at the land surface but makes up about 66 percent of inflow to the aquifer in Virginia, with much of the remaining inflow occurring laterally from outside of defined groundwater budget regions in Virginia. For several decades prior to 2010, high rates of withdrawal from the Potomac aquifer resulted in substantial decline in groundwater storage in the aquifer and in most overlying aquifers and confining units. From 2010 to 2023, rates of withdrawal substantially lower than the historical maximum resulted in small net increases in groundwater storage in the confined aquifer system for most regions of the Virginia Coastal Plain. Nevertheless, for the same period, groundwater storage for the entire model domain continues to incrementally decline, indicating that storage recovery in Virginia is offset by a continued decrease in storage in areas beneath the Chesapeake Bay or adjacent areas of Maryland and North Carolina. Withdrawals from the Potomac aquifer have induced substantial downward flow which is a large part of groundwater budgets for confined aquifers such as the Potomac. For the most recent simulated conditions (2023) downward groundwater flow continues, but because vertical flow rates are a function of the difference between water pressure in the upper surficial systems and lower confined units, rates of downward flow are lower than those in earlier decades as the confined water levels partially recover from larger groundwater withdrawals in the past. Geographically, groundwater flow is generally inward from perimeter regions of the Virginia Coastal Plain toward central regions with the largest withdrawal rates. Groundwater inflow from coastal regions could be contributing to saltwater intrusion, even though that was not measured in this study.

Analyses of groundwater budgets from the Virginia Eastern Shore peninsula, a geographic region of the Virginia Coastal Plain, indicate that groundwater flow for that isolated aquifer system is generally outward from the surficial aquifer to coastal water bodies and downward into the confined Yorktown-Eastover aquifer system, which is the source of most withdrawals. Downward groundwater flow into the confined Yorktown-Eastover aquifer system is estimated to be less than 2 percent of total recharge and less than 9 percent of net recharge at the water table but makes up more than 93 percent of all inflow to the confined aquifer system. Decades of substantial but relatively consistent groundwater withdrawals have induced greater downward flow rates into the confined aquifer system but also have resulted in loss of groundwater from storage. For the most recent simulated period (2023), estimated storage loss accounts for slightly under 7 percent of withdrawals from the confined aquifer system. The reported withdrawal rate for this period from the confined Yorktown-Eastover system is near the highest reported rate for the Virginia Eastern Shore, which means that the storage depletion is expected to continue, even though groundwater levels appear to be relatively stable. Estimated groundwater flow rates upward from the confining unit underlying the Yorktown-Eastover system and low rates of inflow from coastal water bodies underscore ongoing concerns about up-coning and lateral intrusion of salty groundwater.

Released April 16, 2026 12:50 EST

2026, Scientific Investigations Report 2026-5003

Nicholas A. Sutfin, Stephen J. Breen

Ice jams occur regularly at the southern outlet of Osoyoos Lake, which spans the border between the State of Washington and British Columbia, Canada. In recent winters, ice jams caused (1) decreases in downstream discharge that may adversely affect salmon spawning habitat and (2) short-duration lake-level rise that can interfere with lake level management agreements. In response, water managers sought to understand the environmental conditions associated with the historical ice-jam occurrences on Osoyoos Lake. Researchers compiled datasets of discharge, lake level, and air temperature from four meteorological and three hydrologic stations near Oroville, Washington, to determine “ice-jam” or “non-ice-jam” days from 1942 to 2024.

After confirming known ice jams since 1994 using Landsat 8–9 and Sentinel–2 satellite imagery along with discharge, lake level, and air temperature data, researchers designated ice-jam days. They conducted statistical analyses to examine environmental conditions associated with ice-jam occurrences on Osoyoos Lake. Statistical tests indicated significant differences in wind speed, wind direction, and air temperature between ice-jam and non-ice-jam days. A linear discriminant-analysis model correctly predicted 12 of 13 historical ice-jam days since 1994 and determined that ice jams are more likely under westerly and northwesterly winds near or above 10 kilometers per hour (km/h) and minimum temperatures near or below –9.4 degrees Celsius (°C). An analysis of historical discharge suggests that ice jams have occurred since at least the 1940s, but 13 ice jam days occurred in the past decade (2014–2024), exceeding any previous decade. The daily minimum air temperature in the Osoyoos Lake region has increased at a rate of 0.021 °C per year since the 1940s, but ice jams usually occur in winters with colder average temperatures.

Released April 16, 2026 09:54 EST

2026, Environmental Reviews (34) 1-7

Steve J. Cooke, Abigail Lynch, David Tickner, Robin Abell, Morgan L. Piczak, Angela H. Arthington, Michele Thieme, Denielle M. Perry, J. Robert Britton, Tatenda Dalu, Kim Birnie-Gauvin, Steve J. Ormerod, Fernanda Ayaviri Matuk, Rajeev Raghavan, John P. Smol

No abstract available.

Released April 16, 2026 09:50 EST

2026, Water Biology and Security

Cienna R. Cooper, Jane S. Rogosch, Nathan G. Smith, Clinton R. Robertson, Wade M. Wilson

In intermittent stream systems, including those occurring in Texas, USA, the severity of low-flow conditions, duration of seasonal disconnection, and frequency of no-flow events have been amplified by drought. Documentation of these no-flow events is necessary to evaluate ecosystem health. However, many intermittent reaches remain un-gauged given that perennial river sec-tions are often prioritized for gauge placement. Our objectives were to 1) document stream flow using Stream Temperature, Intermittency, and Conductivity (STIC) loggers to determine the frequency and duration of no-flow events in intermittent tributaries of the Colorado River, Texas and 2) compare logger data to publicly available data from streamflow discharge and precipitation gauge networks to understand differences among these data types for drying event characterization. We use these comparisons to summarize benefits and limitations of the application of in-stream data loggers. STIC loggers were deployed at 19 sites, one in each pool and riffle habitat of a stream reach. STIC loggers recorded a measurement of relative conductance every six hours from June 2022 to March 2024, which was used to determine the presence or absence of flow connectivity in a reach. No-flow duration among intermittent reaches varied between 37 and 270 days across tributaries during an ongoing drought in the study area. Overall, logger data was more precise than discharge data for characterizing no-flow events or precipitation data when documenting presence of water in the stream channel due to runoff. Lack of discharge gauges in intermittent tributaries left large sections of stream reaches undocumented and resulted in mischaracterization of flow patterns. Drought severity across the tributaries did not follow longitudinal patterns that would be expected by the climatic precipitation gradient of the study area. More research is needed to determine if factors such as population size affect severity. Likewise, precipitation data did not correlate well with logger water presence data, lacking consideration for groundwater recharge, soil hydrophobicity, and surface compaction. This study shows that to monitor no-flow events, detailed spatial datasets are necessary and that STIC loggers are useful tools that provide data to fill spatial information gaps and facilitate more accurate flow characterization and water presence data in intermittent systems.

Released April 16, 2026 09:40 EST

2026, Science (395) 285-295

John J.Y. He, Ryan S. Crow, John R. Douglass, Christopher S. Holm-Denoma, Jorge A. Vazquez, Brian F. Gootee, Marsha I Lidzbarski, Laura Pianowski, Harrison J. Gray, Emma Heitmann, Phil Pearthree, Kyle House, Shannon Dulin

The timing and mechanism of the integration of the Colorado River and incision of the Grand Canyon remain among geology’s enduring controversies. A key question is the configuration of the upper Colorado River watershed between 11 and 6 million years ago. In this study, we present new evidence from zircon uranium-lead geochronology for the arrival of distinctive Colorado–Green River sediment in the Bidahochi basin by 6.6 million years ago derived from the Browns Park Formation. This is coeval with an order-of-magnitude increase in depositional rate, an increase in carbonate strontium isotope (⁸⁷Sr/⁸⁶Sr) ratios, the appearance of large fish species characteristic of fast-flowing waters, and other sedimentological changes. This evidence is consistent with the Colorado River supplying water and sediment to the Bidahochi basin before spillover integration of the river through the Grand Canyon.

Released April 16, 2026 08:18 EST

2026, Geosphere

Andrew Paul Graber, Matthew A. Thomas, Jason W. Kean, Jonathan Michael King, Jaime Kostelnik

Emergency assessments of postfire debris-flow hazards that are performed by the U.S. Geological Survey (USGS) provide estimates of debris-flow likelihood and rainfall triggering conditions that are used for evaluating and managing runoff-generated debris-flow hazards in recently burned areas throughout the western United States. Although the immediate postfire period, within roughly one year after fire, is typically the most susceptible to runoff-generated debris flows, the hazard evolves in time and space as the burned area recovers. The recovery trajectory a given burned area will take depends on local climate and weather and can be difficult to predict. Some burned areas recover quickly, whereas others experience debris flows for multiple years after fire. As a result, extending our ability to update debris-flow likelihood estimates and rainfall thresholds based on observed recovery of the burned area would be beneficial. We present a method for multi-year runoff-generated debris-flow hazard assessment that leverages the USGS “M1” debris-flow likelihood model and integrates updated, satellite-derived, normalized burn ratio data to estimate vegetation recovery. We predict recovery-aware rainfall thresholds and validate them against a multi-year debris-flow hazard prediction and could be adapted for use with other debris-flow models that incorporate burn severity data.

Released April 15, 2026 15:25 EST

2026, Scientific Investigations Report 2026-5140

Jon P. Mason

The Fort Apache Reservation in east–central Arizona, home to the White Mountain Apache Tribe of the Fort Apache Reservation, Arizona, contains several climate zones because of the large variation in surface elevation within the reservation. This study was carried out in cooperation with the White Mountain Apache Tribe of the Fort Apache Reservation, Arizona, to raise awareness of how the changing climate affects the Fort Apache Reservation. This report documents the evaluation of existing multidecadal meteorological and hydrological datasets for the Fort Apache Reservation, used to evaluate the effects of a changing climate on the reservation. In this evaluation, near-surface air temperature, snow depth, snow water equivalent, precipitation, and streamflow datasets were analyzed for monotonic trends indicative of changing climatic conditions during specified periods of time. The results of these trend analyses were then compared with the Tribal community's memories of the changing climate.

Trend analysis of near-surface air temperatures from a U.S. Historical Climatological Network station on the Fort Apache Reservation at Whiteriver, Arizona, indicated that mean annual air temperatures have increased by an average of 2.48 degrees Fahrenheit from 1980 to 2023. Records from the same station also indicated that average monthly maximum temperatures recorded for March increased by 5.39 degrees Fahrenheit for the same time period.

Annual precipitation at the five precipitation stations used in this study decreased greatly from the 1980s to 2023. The largest total decrease was 10.07 inches, or 34.7 percent. However, only one of the two precipitation stations with longer term data available prior to 1980 had a significant negative trend when data from the entire period of record, from 1901 to 2023, were analyzed.

Trend analyses show a decrease in the annual maximum snow water equivalent and an earlier disappearance of the snowpack at two Natural Resources Conservation Service snow telemetry stations in the mountainous region just east of the Fort Apache Reservation from 1981 to 2023. Based on the trend analyses, the average annual maximum snow water equivalent has decreased by more than 40 percent at both stations, and the average date when the snowpack was fully melted at the stations in the spring has moved earlier in time from late April to early April or late March. However, a statistically significant trend was not determined for the early April snow water equivalent measured at a nearby Natural Resources Conservation Service snow course across its period of record, indicating that the history of mountain snowpack in this area is not fully understood. Analysis of snowfall data from a National Oceanic and Atmospheric Administration Cooperative Observer Program network station on the Fort Apache Reservation at McNary 2N, AZ (station 025412) indicated that, on average, the measured total annual snowfall at the station decreased 42.4 percent from 1935 to 2023.

Streamflow data from six U.S. Geological Survey streamgages on the Fort Apache Reservation were analyzed for trends. For most streamflow gages, statistically significant trends were not determined for tested parameters when the entire streamflow period of record was used for stations with records going back to at least the 1960s. However, when the data from 1980 to 2023 was tested, most of the streamflow parameters had statistically significant negative trends. All six streamgages showed a decrease in average annual runoff of at least 50 percent from 1980 to 2023; one streamgage showed an 81.8 percent decrease.

A similar statistical finding was observed in the analysis of the annual spring snowmelt peak from one of the six streamgages used in the study and located in an area receiving measurable amounts of snowmelt runoff. When data from the entire period of record (1958–2023) was used, no trend in streamflow was determined; however, a significant negative trend was determined from 1980 to 2023, indicating a decrease in average annual springtime runoff of 62.6 percent. Statistical analysis on the timing of the annual spring snowmelt peak at the same streamgage indicated the snowmelt peak is happening on average about 12 days earlier now (2023) than it did in the past. The trend results for the timing of the annual spring snowmelt peak were the same and statistically significant for both periods tested (1958–2023 and 1980–2023). Two of the streamflow records from the Fort Apache Reservation were compared to the Palmer Hydrological Drought Index computed for Arizona Climate Division 4 (East Central) by the National Centers for Environmental Information. The comparison showed that the streamflow records generally tracked the Palmer Hydrological Drought Index.

In interviews, Tribal community members living on the Fort Apache Reservation described the changes in climate that they observed during their lifetimes. Common themes reported were that air temperatures have become warmer, and the weather is less predictable with changes in seasonal patterns. Drier conditions, lower snowfall, shorter winters, and lower river levels were also reported. These community member observations align with the results of this study.

Released April 15, 2026 10:20 EST

2026, Professional Paper 1890-J

Thomas W. Sisson, Andrew T. Calvert

Cenozoic volcanic rocks of the Arabia Plate cover about 140,000 square kilometers across a distance of about 3,000 kilometers from southern Yemen to southeastern Turkey. The majority of volcanic products are alkali basalts that erupted in restricted areas, commonly over periods of a million or more years, building mafic lava fields, each known in Arabic as a “harrat.” Harrat volcanism commenced following the Oligocene flood-lava effusions that blanketed the (now) Ethiopian highlands, southern Sudan, and western Yemen, and overlapped the latest Oligocene to early Miocene initial riftings of the Red Sea and Gulf of Aden, but the majority of harrat volcanism has been since approximately 13–10 million years ago. Persistent harrat magmatism in restricted locations led to the development of intermediate and evolved magmas (hawaiites, mugearites, benmoreites, trachytes, comendites, and phonolites) mainly through intracrustal crystallization-differentiation. Most of these intermediate and evolved magmas erupt at sites of the greatest aggregate volcanic relief, reflecting sites of the greatest overall magmatic fluxes. Production of fractionated magmas at these sites negates “monogenetic” as an appropriate descriptor of harrat volcanism. This chapter summarizes the geologic, eruptive, and tectonic history and aspects of the petrogenesis of the Cenozoic Arabian alkalic province. Particular emphasis is placed on results of a joint study of Ḩarrat Rahat adjacent to the city of Al Madīnah al Munawwarah, Kingdom of Saudi Arabia, published as U.S. Geological Survey Professional Paper 1862 and Saudi Geological Survey Special Report SGS–SP–2021–1. A goal of this chapter is to provide an introduction to those unfamiliar with this vast, enigmatic, and fascinating region of distributed continental volcanism.

Released April 15, 2026 10:15 EST

2026, Ornithology

Erin R. Stewart, Courtney J. Conway

Unequal representation of focal research areas can arise during the initial stages of project development when investigators make decisions about what, when, and where to study. Regarding when research is conducted, publications on vertebrates are strongly skewed toward breeding-stage studies, leaving sizeable gaps in our knowledge pertaining to behavior and demography in nonbreeding stages. Here, we quantified how the focus on different annual cycle stages has changed over the past 64 years in the field of ornithology, and tested 3 hypotheses to explain underlying correlates of the focus on different stages. We surveyed field-based studies published in Ornithology from 1960 to 2024 and documented the annual cycle stages examined, regions of study, migratory strategies of focal species, and author affiliation for each paper. The tendency to conduct breeding-stage research in ornithology became more, rather than less, pronounced over the 64-year interval, and breeding-stage research was more common when field work focused on migratory species. Therefore, investigators, authors, and editors could likely increase representation of other annual-cycle stages by supporting, conducting, and publishing more studies in the tropics and more studies using remote tracking technologies. More nonbreeding and year-round studies are necessary to fully understand the ecology and evolution of species and develop the most effective strategies for conservation.

Released April 15, 2026 08:14 EST

2026, Water Research (294)

Laura T. Kelly, Daniel G Beach, Joanna R. Blaszczak, Keith Bouma-Gregson, Sydney M. Brown, Haowu Cheng, Janette L. Davidson, Jutta Fastner, Marcus Francis, Andrea Garcia Jimenez, Laurel Genzoli, Ramesh Goel, Diego Gonzalez, Kim M. Handley, Sabine Hilt, Jean-Francois Humbert, Rob Jamieson, Lindsay Johnston, Pilar Junier, Janice Lawrence, Pearse McCarron, Sven Meissner, Jacob Mormando, Jonathan Puddick, Catherine Quiblier, Nagasaijanani Rajpirathap, Charlotte Schampera, Andy Selwood, Karen Shearer, Abeer Sohrab, Rosalina Stancheva, Cecilio Valadez-Cano, Jordan M. Zebrecky, Susanna A. Wood

There have been sporadic reports of aquatic, benthic Microcoleus proliferations in freshwater rivers, lakes, and reservoirs for four decades, with reports increasing in frequency over the last twenty years, suggesting a possible rise in their global distribution, frequency, and intensity. Microcoleus can produce anatoxins which are neurotoxic, and ingestion of toxic mats has caused hundreds of dog fatalities and raised serious human and ecological health concerns. This review synthesizes and evaluates current knowledge on Microcoleus distribution, taxonomy, toxin production, toxicity, ecology, environmental drivers, and biotic interactions. Toxin-producing Microcoleus have been reported in at least 18 countries, though many regions have not conducted toxin testing, suggesting a broader but under-reported distribution. Proliferations occur across diverse habitats, including cobble-bedded streams, large sandy rivers, reservoirs, and lakes. Microcoleus proliferations also occur on macrophytes, both in lakes and rivers. Genomic analyses currently classify anatoxin-producing Microcoleus into distinct species, with all known anatoxin-producers isolated from freshwater ecosystems. Anatoxin concentrations vary widely over space and time, within and among waterbodies. While studies on environmental drivers remain limited, research in cobble-bedded rivers suggests that moderate enrichment of dissolved inorganic nitrogen and low dissolved reactive phosphorus concentrations in the water column promote proliferation. Metagenomic approaches have revealed unique nutrient acquisition and storage strategies used by Microcoleus. Key knowledge gaps remain around the environmental and ecological triggers of proliferation, toxin production, genomic diversity and microbial interactions. Addressing these gaps through coordinated, global studies using robust datasets and consistent methods is critical to improve prediction, monitoring, and mitigation of this increasingly widespread public and ecological health threat.

Released April 14, 2026 15:46 EST

2026, General Information Product 265

Holly F. Weiss-Racine, Joseph A. Bard, Jessica L. Ball, Carolyn L. Mastin

Introduction 

Eruptions at Mount Rainier produce lava flows, plumes of airborne volcanic ash, and avalanches of hot rock, ash, and gas—pyroclastic flows—that rush down the steep, ice-covered slopes of the volcano. Hot rock and ash ejected during an eruption can melt large quantities of snow and ice, forming huge, fast moving mudflows called lahars that travel 30+ miles, all the way to Puget Sound. Very large lahars can also form when weak and water-saturated rock high on the volcano collapses with or without volcanic activity. Learn more inside!

Released April 14, 2026 14:40 EST

2026, Open-File Report 2026-1061

Meredith Leigh McPherson, Joshua B. Logan, Kristen Alkins, Sarah Groff, Gerry A. Hatcher, Ann E. Gibbs, Susan Cochran, Curt D. Storlazzi

Executive Summary 

Coral bleaching events have become increasingly common across the Hawaiian Archipelago since 1996 because of more frequent and intense marine heatwaves. The most significant bleaching event to date occurred from 2014 to 2015, which resulted in catastrophic state-wide coral loss. Bleaching events with less severe effects also occurred in 1996 and 2019. To understand the long-term effects of repeated bleaching events, along with other anthropogenic factors such as water quality, storms, sewage runoff, and coastal development, on coral reefs on the Kona Coast of the Island of Hawaiʻi, the U.S. Geological Survey, in collaboration with the National Park Service, collected underwater imagery in the early 2000s (baseline survey) and again in 2022 (resurvey). These images were captured within and adjacent to the National Historic Parks (NHP) and National Historic Sites (NHS) of Kaloko-Honokōhau NHP (KAHO), Puʻuhonua o Hōnaunau NHP (PUHO), and Puʻukohola Heiau NHS (PUHE). Imagery was classified for live coral cover and dominant type (four coral types, rubble, macroalgae, and two bottom substrate types). Change of percent live coral cover was determined for all sites. Change of coral and non-coral dominant types were calculated by aggregating classifications for each park into coral and non-coral. Net coral cover decreased between the baseline and resurvey period across all three parks, though PUHE exhibited the greatest loss of live coral cover. Across all three parks, the occurrence of lower coral cover classes (0–20 percent) increased and higher coral cover classes (greater than 50 percent) decreased. Furthermore, the total occurrence of non-coral dominant type classifications (rubble, macroalgae, sand, and volcanic pavement) increased by approximately 25 percent across all three parks, with PUHE experiencing a nearly 90-percent increase in the occurrence of non-coral types. There was little to no effect of water depth on change of live coral cover, indicating that marine heatwave driven bleaching events and additional anthropogenic influences affected the entire reef across all water depths from the lower fore reef to the reef flat.

Because coral loss was more severe at PUHE and PUHO than KAHO, creating a monitoring framework that utilizes periodic underwater camera surveys and fixed diver transects by the National Park Service would contextualize the periodic spatial surveys to the fixed transects that have greater temporal resolution. Similarly, increased frequency of spatial surveys would allow for the National Park Service to continue monitoring changes to critical nearshore habitats and marine resources relevant to National Park jurisdiction.

Released April 14, 2026 12:18 EST

2026, Scientific Investigations Report 2026-5011

Brendan M. Foster, Carly Marcella Maas, Alejandra Logan Flota

Despite the extensive natural gas pipeline network in the United States that intersects streams and other sensitive habitats, few case studies use a comparative upstream-to-downstream approach to evaluate potential short- and long-term effects of pipeline stream crossings from pre-construction through post-restoration. In 2017, the U.S. Geological Survey, in cooperation with the Virginia Department of Environmental Quality, deployed real-time continuous stream monitoring stations upstream and downstream from six proposed Mountain Valley Pipeline stream crossings in southwestern Virginia. Water temperature and turbidity data collected at the upstream and downstream sites were compared across three periods—before stream crossing construction, during stream crossing construction, and after stream crossing construction—to determine potential influences from the pipeline stream crossing. Additionally, the monitoring network was used to notify regulators of potentially anomalous conditions throughout the entire monitoring period.

The results of this study indicate that pipeline stream crossing did not affect long-term or short-term upstream-to-downstream water temperature conditions or long-term upstream-to-downstream turbidity conditions in any of the six monitored streams. Some short-term anomalously elevated turbidity conditions were observed and attributable to pipeline stream crossing; however, the magnitudes and durations were not sufficient to alter the long-term turbidity regimes of the streams in which they were observed. The application of the monitoring network as a real-time alert system successfully alerted regulators to potentially anomalous conditions.

Released April 14, 2026 08:35 EST

2026, Planetary Science Journal (7)

Jean-Pierre Williams, Sarah Valencia, Kristen A. Bennett, Margaret E. Landis, Kerri L. Donaldson Hanna, Addison T. Dove, Patrick O'Brien, Brett W. Denevi, Justin Hagerty, Craig Hardgrove, Paul O. Hayne, Adam LaMee, Thomas H. Prettyman, Katherine A. Shirley, Matthew A. Siegler, Jessica M. Sunshine

Characterizing terrain surface properties is an essential step in assessing the feasibility of landing successfully at a location on a planetary surface. Slopes and terrain ruggedness index (TRI) values derived from high-resolution (2 m pixel⁻¹) digital terrain models provided important constraints in selecting the landing site for the upcoming Payloads and Research Investigations on the Surface of the Moon program as part of the Commercial Lunar Payload Services task order CP-21 mission. The selected landing site needed to balance safety requirements with the ability to achieve the science and exploration goals of the Lunar Vulkan Imaging and Spectroscopy Explorer payload. In this study, we compare several morphometric parameters in the context of the CP-21 landing site on Mons Gruithuisen Gamma, or the Gamma dome, and quantify the information they convey about lunar surface properties to assess their utility for future landing site evaluation. TRI was found to be a useful metric for assessing landing site safety. Metrics that better decouple slope and surface roughness, the vector ruggedness measure and the standard deviation of slope, provided additional information about surface characteristics and textures such as the degree to which roughness is isotropic.

Released April 13, 2026 11:56 EST

2026, Open-File Report 2026-1001

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

Preface 

In anticipation for increased funding made possible by the Water Resources Development Act of 2020, the Upper Mississippi River Restoration (UMRR) Program identified a need to conduct river-wide assessments of floodplain vegetation. In January 2025, we assembled a group of subject matter experts to perform the following tasks:

1. Review Upper Mississippi River Restoration’s current floodplain vegetation research portfolio,
2. Identify important features and goals for long-term floodplain vegetation monitoring,
3. Evaluate the suitability of existing datasets for system-wide vegetation assessments, and
4. Discuss emerging opportunities to learn about floodplain vegetation dynamics from local-scale restoration and management projects.

This document is a summarization of what occurred at the meeting and provides suggested next steps toward developing the capacity to conduct routine long-term monitoring and assessment of floodplain vegetation as part of the Upper Mississippi River Restoration Program.

Released April 13, 2026 11:50 EST

2026, Circular 1563

Adam T. Ringler, David C. Wilson, Robert Anthony, Corey I. Beutel, Andrew Holcomb, Charles R. Hutt, Tom Telesha

This circular presents a strategic outlook for the U.S. Geological Survey (USGS) Albuquerque Seismological Laboratory (ASL) for the next 10 years (2026–36). The ASL is a USGS field office in the Geological Hazards Science Center that operates portions of the Advanced National Seismic System and the Global Seismographic Network and focuses on fundamental research for instrumentation testing and data quality. The strategic outlook is categorized into two types of tasks: “Foundational Tasks” and “Aspirational Tasks.” Foundational Tasks are those that maintain the laboratory’s basic operations and services, including regional and global seismic monitoring, improving data quality, and providing instrument testing and support. A suite of Aspirational Tasks is also articulated; these can be considered priority targets of ASL that could improve ASL’s seismic monitoring capabilities and mission. Such tasks include improvements to remote stations, testing capabilities of nonseismic geophysical instruments, detection threshold monitoring, rapid aftershock deployments, and expanding seismic monitoring networks. This report was written with input from the USGS Geological Hazards Science Center, the USGS Earthquake Hazards Program (EHP), and colleagues with an interest in the work done by the ASL. Although the details of these tasks may change, this document can provide guidance on the overarching tasks at the ASL from 2026 to 2036 and an overview of the various components of the ASL and how they fit into the EHP and the Global Seismographic Network Program.

Released April 13, 2026 11:44 EST

2026, Scientific Investigations Report 2025-5093

Andrea Cecile Jaegge, Keith Bouma-Gregson, Thomas D. Byl, Kurt D. Carpenter, Victoria Christensen, Rebecca M. Gorney, Jennifer L. Graham, Heather A. Heckathorn, Hayley T. Olds, Pamela A. Reilly, Joshua J. Rosen, Michael D.W. Stouder

Cyanobacterial toxins (cyanotoxins) represent a substantial threat to drinking water supplies and safe recreational uses of freshwater resources in watersheds worldwide. Monitoring cyanotoxins can be difficult because toxin events are variable in both space and time, are not always persistent, can be moved easily by wind and currents, and may be degraded biotically or abiotically. Thus, monitoring programs that collect discrete samples on a monthly or even bimonthly interval can miss key events and underestimate cyanotoxin risk or if they capture a high-concentration event, can give a false impression that cyanotoxins are a widespread health hazard. The use of Solid Phase Adsorption Toxin Tracking (SPATT) samplers helps address this issue by providing a time-weighted average estimate of dissolved cyanotoxin occurrence and relative concentrations. SPATT samplers have been used as a complement to traditional monitoring programs and can help elucidate cyanotoxin dynamics. SPATT samplers have been used by six U.S. Geological Survey (USGS) Water Science Centers (New York, California, Oregon, Upper Midwest, New Jersey, and Lower Mississippi-Gulf) to monitor various cyanotoxins in waterbodies such as streams, rivers, lakes, waterfalls, estuaries, and drinking-water intakes. Despite their use across the USGS, there is little guidance available to ensure consistent approaches and data quality across the Bureau. This report summarizes best practices for SPATT deployment and analysis, synthesizes data and describes lessons learned from USGS studies, identifies priority knowledge gaps, and offers considerations for future targeted experiments to help improve data collection and interpretation.

Released April 13, 2026 10:21 EST

2026, Report

Sarah E. Whipple, Sarah K. Carter, Zachary Ormsby, Benjamin Hoke, Michael Powers, Ryan E. O'Dell, Rebecca Spitzer

This science strategy describes the scientific mission of the unit, science previously conducted on or near the unit, and the partners who have given considerable effort to produce this science. Further, it identifies priority science needs, outlines the unit’s plan to meet those science needs in coordination with partners, shares scientific protocols for conducting new research, and identifies systems of communication to help ensure science information generated from this research is shared throughout BLM, with partners, and with the public. These components of the science strategy are intended to help achieve the fundamental goals of the effort: to support and expand partnerships while identifying priority science information needs of unit staff and partners and clearly outlining pathways for conducting, communicating, and applying that science.

Released April 13, 2026 09:35 EST

2026, Chemical Geology (712)

Una Farrell, Hunter Olson, Maya Thompson, Michelle L. Abshire, Oyeleye Adeboye, Anne-Sofie Ahm, Lewis Alcott, Thomas Algeo, Ross Anderson, Arif Ansari, Lucas Bastos, Kohen Bauer, Brian Beaty, Justin E. Birdwell, Fred Bowyer, Jochen J. Brocks, Tessa Brunoir, James F. Busch, Donald Canfield, Fabrício Caxito, Chao-Chin Chang, Meng Cheng, Jean Clemente, David Cordie, Peter W. Crockford, Huan Cui, Celeste Cunningham, Tais Dahl, Janaina Rodrigues de Paula, Carol Dehler, Lucas Del Mouro, Keith Dewing, Dermeval Aparecido do Carmo, Stephen Dornbos, Nadja Drabon, Julie A. Dumoulin, Omabehere Ejeh, Emily Ellefson, Maya Elrick, Joseph Emmings, Bokanda Eric, Hao Fang, Gabriella Fazio, Henrique Fernandes, Katherine L. French, Robert R. Gaines, Richard Gaschnig, Timothy M. Gibson, Geoffrey J. Gilleaudeau, Karin Goldberg, Zheng Gong, Amy Hagen, Galen P. Halverson, Kalev G. Hantsoo, Emma R. Haxen, Miles A. Henderson, João-Pedro Hippertt, Malcolm S. W. Hodgskiss, Paul Hoffman, Edward Huang, Benjamin Johnson, Pavel B. Kabanov, Junyao Kang, C. Brenhin Keller, Brian Kendall, Julien Kimmig, Sara Kimmig, Michael A. Kipp, Andrew H. Knoll, Timmu Kreitsmann, Anurag Kulkarni, Alexandra Kunert, Marcus Kunzmann, Jiankang Lai, Richard O. Lease, Chao Li, Sen Li, Alex Lipp, Yang Liu, David Loydell, Xinze Lu, Katie Maloney, Kaarel Mänd, Alexie E. G. Millikin, N. Tanner Mills, Kento Motomura, Chiza N. Mwinde, Lyle Nelson, Nora M. Nieminski, Brennan O'Connell, Edel Mary O'Sullivan, Juliana Okubo, Jaden Olah, Frantz Ossa Ossa, Chadlin Ostrander, Kärt Paiste, Camille A. Partin, Egberto Pereira, Shanan Peters, Tiffany L. Playter, Susannah M. Porter, Simon W. Poulton, Sara B. Pruss, Zhen Qiu, Daven Quinn, Mariano Remirez, Sebastian Richiano, Sylvain Richoz, Kathryn Rico, Samantha R. Ritzer, Zachary Roney, Alan D. Rooney, William C. Rose, Elias Rugen, Swapan K. Sahoo, Shane D. Schoepfer, Judith A. Sclafani, Nathan D. Sheldon, Yanan Shen, Graham Shields, Pulkit Singh, Arvind Singh, Sarah P. Slotznick, Emily Smith, Haijun Song, Sam Spinks, Richard G. Stockey, Justin Strauss, Eva Stüeken, Zongyuan Sun, Dongjie Tang, Lidya Tarhan, Danielle Thomson, Nicholas Tosca, Rosalie Tostevin, Chenyi Tu, Maoli Vizcaíno, Yuxuan Wang, Changle Wang, Xiaomei Wang, Lucas Warren, Lucy Webb, Philip R. Wilby, Christina R. Woltz, Rachel Wood, Yuyang Wu, Xiuqing Yang, Inessa A. Yurchenko, Junpeng Zhang, Jessica Whiteside, Benjamin C. Gill, Akshay Mehra, Kimberly V. Lau, Noah J. Planavsky, David T. Johnston, Erik A Sperling

Geochemical data from sedimentary rocks are the primary source of information regarding Earth's surface evolution through time, including its air and water envelopes and interactions with life and deep Earth processes. The Sedimentary Geochemistry and Paleoenvironments Project (SGP) is a scientific consortium centered around open data and community-driven development of cyberinfrastructure tools and resources for sedimentary geochemistry and Earth history. Here we describe the SGP Phase 2 data release, which focused on incorporating Paleoproterozoic and Mesoproterozoic (2500–1000 million years ago) data and better accommodating carbonate data. This data release was built through the involvement of >200 researchers worldwide in academia, government, and industry, and provides the largest available public data resource for our user community in the academic fields of geochemistry, sedimentology, tectonics, paleontology, Earth history, and paleoclimate, as well as the petroleum and minerals industries. The dataset now encompasses 126,006 samples and 4,132,371 geochemical analyses. In addition to direct entry by SGP Team Members, we have ingested and incorporated datasets from the Geoscience Australia OZCHEM database, the Alberta Geological Survey, and the Deep-Time Marine Sedimentary Element Database (DM-SED) compilation. This paper details sampling in the Phase 2 dataset with respect to age, geography, lithology, and other geological characteristics, documents access via our search website and API, discusses possible issues and/or biases in the dataset that could impact analyses, describes plans for governance and stewardship of data from Indigenous lands, and serves as the citable reference paper for the data release.

Released April 13, 2026 08:27 EST

2026, Ecological Informatics (95)

Kaitlin Huber, J. Hardin Waddle, Brad M. Glorioso, Therese M. Donovan

The Cuban treefrog (Osteopilus septentrionalis), as an invasive species in the southern United States, presents a need for effective surveillance. Automated detection expedites processing of audio data for large-scale surveillance and monitoring programs. However, current available methods commonly used for anuran species have not been sufficient to detect Cuban treefrogs. Here, we present results from a two-stage method for automated detection that employs both cross-correlation template matching and secondary supervised learning classifiers. In the first stage, audio data are screened for initial detections using template matching, in which the detections contain both true and false positives. In the second stage, the false positives are screened out using classifier algorithms. We used this method to process 139,985 audio recordings, consisting of 596,046 total minutes, collected at 13 locations in Louisiana and Florida from 2014 to 2022. From the stage 1 template matching, we detected 83,191 Cuban treefrog signals across recordings. The stage 2 machine learning model was able to identify stage 1 false positive detections with a testing accuracy of 98.46% and a testing false positive rate of 1.116%. After pruning false positive detections, a total of 20,271 individual Cuban treefrog detections remained, distributed mainly across 3 sites in an area with known presence. Locations with presumed absence had an easily verifiable number of false positive detections (n = 109 across all other sites). The two-stage methodology utilizing both template matching and machine learning algorithms can be integrated into wildlife surveillance or monitoring programs for species with distinctive, conserved calls as an effective way to achieve sensitive species detection with a low incidence of false positives.

Released April 13, 2026 07:47 EST

2026, Genome (69)

Michaela Ray Grossklaus, David S. Pilliod, Stephen Frank Spear, Matthew B B. Laramie, Akhil Kholwadwala, Amanda Jean Boone, Yer Lor, Marissa Kaminski, Jeffrey G. Everett

Environmental DNA (eDNA) sampling is a noninvasive alternative to conventional methods of surveying insects that may be particularly useful for detecting pollinators. We developed a quantitative polymerase chain reaction (qPCR) assay to detect the DNA of Franklin’s bumble bee (Bombus franklini) from flower samples and conducted an initial test of the assay using samples collected within and around the historical range of the species. We further analyzed all samples using metabarcoding. Our qPCR assay successfully amplified B. franklini DNA and exhibited no cross-reactivity with nontarget bumble bee DNA during in silico and in vitro testing. We did not detect B. franklini DNA from field-collected flower samples using either qPCR or metabarcoding. However, metabarcoding analysis revealed DNA of at least 16 other bumble bee species. This finding underscores the potential utility of eDNA sampling for surveying bumble bees. Nondetection of B. franklini from field-collected flower samples may be due to the extreme rarity of the species; B. franklini is endangered and has not been observed in the wild since 2006. Our B. franklini assay is among the first bee-specific qPCR assays ever developed and provides proof of concept for additional assays that may improve detection rates of rare and endangered bees.

Released April 13, 2026 07:37 EST

2026, Environmental Research: Ecology (5)

Paul C. Selmants, Benjamin M. Sleeter, Colin J. Daniel

Uncertainties in terrestrial ecosystem models limit their predictive power. Efforts to reduce projection error have rarely focused on constraining uncertainty in the initial state of the ecosystem, however, despite evidence that matching model initial conditions to real-world observations reduces overall model bias. Here we use an integrated land change and carbon gain-loss model to evaluate the influence of initial condition uncertainty on simulations of California wildland ecosystems during the years 1985–2020. We generated 36 initial conditions scenarios by varying the source data used to initialize state variables and then ran simulations based on each of these scenarios under a constant set of historical conditions. We found that discrepancies in initial forest extent and initial forest age among scenarios generated wide uncertainty ranges in model estimates of terrestrial ecosystem carbon stocks and flux rates at the outset of the simulation period, but differences in initial forest composition had no impact. Over time, forest age became more homogeneous across model scenarios leading to exponential rates of decline in the uncertainty ranges of live biomass and dead wood carbon but little to no impact on uncertainties in litter and soil organic carbon. Uncertainties in individual carbon flux rates were consistent with uncertainties in their source pools. In contrast, model estimates of ecosystem carbon balance demonstrated a shift in system behavior not apparent in trends for individual carbon stocks and fluxes. Specifically, estimates of ecosystem carbon balance converged across scenarios for the first 20 years of the simulation period but then began to diverge at an accelerating rate, possibly due to weakened resilience to the increased frequency and severity of climate-driven disturbances. Our results demonstrate that uncertainty in the initial state of the system can have large and persistent impacts on the predictability of ecosystem carbon dynamics, and that ongoing shifts in external forcing by climate and climate-driven disturbances can exacerbate these impacts.

Released April 11, 2026 08:47 EST

2026, Scientific Data

Chris Funk, Pete Peterson, Laura Harrison, Robert Saldivar, Martin Landsfeld, Diego Pedreros, Shraddhanand Shukla, Andreas H. Fink, Frank Davenport, Seth H. Peterson, William Turner, Austin Sonnier, Michael Budde, Karyn Tabor, James Verdin, Disha Hauzaree, Mohamed Naim, Daniella Alaso, Gregory Husak

The Climate Hazards Center Infrared Precipitation with Stations (CHIRPS) data stream combines: (1) a high-resolution climatology, (2) thermal infrared (TIR) geostationary satellite observations, and (3) station observations. In the past, CHIRPS version 2 (CHIRPS2) has proven to be valuable for drought monitoring, hydrologic modeling, scientific studies and agricultural decision making. Version 3 (CHIRPS3) improves each of these components. The new version, CHIRPS3 extends to 60°S/N, adopts an improved variance-preserving TIR-to-precipitation estimation method, uses many more stations and station sources than the original CHIRPS2 product, and implements gauge-undercatch correction. In this paper, we evaluate the performance of satellite-only CHIRP3, CHIRP2, IMERG, PERSIANN- CCS, and GPI using high quality interpolated data in twelve regions with dense station coverage. CHIRP3 represents both the observed mean and variance more accurately than CHIRP2. A usage section in Morocco shows that CHIRPS3 better captures the observed rainfall variability when compared to CHIRPS2. This section also demonstrates how station data should be gauge-undercatch-corrected when validating CHIRPS3.

Released April 10, 2026 13:33 EST

2017, Article

Richard B. Wanty, Andrew H. Manning, Michaela R. Johnson, Stephen J. Kalkhoff, Jessica D. Garrett, Jean M. Morrison, Stefania Da Pelo, Jeffrey L. Mauk

Released April 10, 2026 11:27 EST

2026, Scientific Investigations Report 2026-5141

Viviana Mazzei, Keith A. Loftin, Emily Karwacki, Jose V. Lopez, Lauren E. Krausfeldt, Barry H. Rosen, Hidetoshi Urakawa

The Caloosahatchee River, located in southwest Florida, is a eutrophic and colored river that flows from Lake Okeechobee westward into its estuary and the Gulf of America. Cyanobacterial harmful algal blooms (HABs) are a documented problem along this freshwater-to-marine waterway where nutrient enrichment has been identified as a key factor in bloom occurrence but has not been experimentally tested in the river. This study is the first to test the effects of inorganic nutrient loading on phytoplankton assemblages in the Caloosahatchee River and the effects of different nutrient sources on phytoplankton dynamics at different times of the year. Three independent, in situ experiments were conducted to test the effects of daily, incrementally increased ammonium, nitrate, and phosphate loading on phytoplankton at different times of the year (summer, fall, winter). Over the 72-hour enclosure period, phytoplankton abundance metrics (cell concentration, chlorophyll-a, and phycocyanin), dissolved oxygen, and pH increased, and fluorescent dissolved organic matter and turbidity decreased in all treatments and controls. Increased phytoplankton abundance metrics relative to controls were observed after 72 hours of exposure to elevated ammonium and nitrate in summer and only ammonium in winter, suggesting periodic nitrogen limitation; however, no treatment effects on phytoplankton assemblage structure in terms of resemblance and diversity metrics were found. Increases in total cell concentrations were driven by elevated growth rates of already dominant taxa but not sufficiently to form a visible bloom. Cyanobacteria consistently dominated the phytoplankton, particularly Aphanocapsa and Merismopedia, whereas the common HAB-forming Microcystis maintained consistently low abundance. This study provides new information on the ecology of phytoplankton assemblages in the Caloosahatchee River and could be used by water resources managers to evaluate strategies for controlling cyanobacterial HABs in the river.

Released April 10, 2026 09:37 EST

2026, Scientific Investigations Report 2026-5128

Rebecca M. Gorney, Heather A. Heckathorn, Kyle R. Clonan, Pamela A. Reilly, Kathryn Cahalane, Bradley W. Bjorklund

Harmful algal blooms, particularly cyanobacteria harmful algal blooms (cyanoHABs), have emerged as a substantial global concern because of their detrimental effects on water quality and aquatic ecosystem health. CyanoHABs can produce cyanotoxins, which pose serious health risks to humans and wildlife, such as liver failure and respiratory distress. This is particularly concerning for water bodies that serve as drinking-water sources. Recent trends indicate an increase in the frequency and intensity of cyanoHABs globally. This study focuses on the Raritan Basin Water Supply Complex in New Jersey, where extensive monitoring was conducted from August 2020 to August 2021 to assess the presence of cyanobacteria and associated cyanotoxins. The research utilized a combination of discrete water-quality sampling, continuous monitoring, and solid phase adsorption toxin tracking (SPATT) to capture the dynamics of cyanotoxin occurrence and potential transport. Findings revealed a widespread presence of cyanobacteria and potential for cyanotoxin production, although actual cyanotoxin concentrations remained below drinking water and recreational thresholds. The study, conducted by the U.S. Geological Survey (USGS) in collaboration with the New Jersey Water Supply Authority (NJWSA) and the New Jersey Department of Environmental Protection (NJDEP), highlighted the limitations of traditional sampling methods and emphasized that continuous monitoring can support better understanding of how cyanoHAB conditions change over time and in different places. Genetic testing included quantitative polymerase chain reaction (qPCR) analyses, which demonstrated higher sensitivity, or increased findings of cyanobacteria compared to microscopy, indicating the potential for use in early warning systems. This research underscores that integrating various detection methods and hydrological data can enhance understanding of cyanotoxin dynamics in river systems.

Released April 10, 2026 07:58 EST

2026, PLOS Sustainability and Transformation (5)

David T.S. Hayman, Steve Unwin, Kelly Bateman, Casey Barton Behravesh, Charlotte Berg, Jemma Bergfeld, Cristina Casalone, Claire Cayol, Erin Davis, Sunday Ekesi, Johan Esterhuizen, Merid Getahun, Federica Giorda, Keith Hamilton, Damien O. Joly, Christa Kuhn, Jean-Claude Manuguerra, Daniel Masig, Anita Michel, Paolo Mulatti, Misheck Mulumba, Annah Njui, Richard Paley, Antonio Fernandez, Sascha Knauf, David Poumo Tchouassi, Youming Wang, Nathalie Vachiery, Jandouwe Villinger, Frank Y.K. Wong, Gongxun Zhong, B. Dharmaveer Shetty

Wildlife health is integral to functioning, complex ecosystems [1], directly and indirectly influencing the health of people, animals, plants, and the environment [2–4]. Healthy wildlife populations are essential for ecosystem services and are at the heart of the One Health approach [3,4], which aims to sustainably balance and optimize the health of people, animals, and ecosystems through multisectoral and transdisciplinary collaboration [5].

Despite its importance, wildlife health initiatives often operate in silos, limiting capacity to address transboundary threats such as emerging diseases, pollution, and environmental changes. Anthropogenic changes, including habitat loss, degradation, fragmentation, and unsustainable harvesting, exacerbate wildlife health challenges [6–9]. These pressures disrupt species biology and alter host-pathogen dynamics [10–12], underscoring the importance of coordinated collective action in addressing harmful effects on the health of wild animals. While local conservation efforts are vital, long-term success in safeguarding biodiversity requires a unified, global network. For instance, without harmonized surveillance and response systems, individual institutions cannot effectively track pathogens across borders or share diagnostic capabilities.

The World Organisation for Animal Health (WOAH) Collaborating Centre Network for Wildlife Health—WOAH-WildNet—was established to bridge these gaps. By fostering global collaboration, sharing resources, and enabling data exchange, WOAH-WildNet provides a transformative, systems-based approach to wildlife health, managing risks, and enhancing ecosystem resilience. Central to this mission is breaking down silos to promote intersectoral coordinated responses to complex wildlife health challenges.

Released April 09, 2026 15:06 EST

2026, Scientific Investigations Report 2026-5138

Matthew T. Moser, Mikaela L. Cherry, Brent M. Hall

The City of Lincoln, Nebraska, has been monitoring concentrations of arsenic in their source water and evaluating their options for treatment and removal since at least 2002. In 2022, the City of Lincoln, Nebr., with funding assistance from the Nebraska Water Sustainability Fund, began cooperating with the U.S. Geological Survey to examine arsenic concentrations in surface water and groundwater in the lower Platte River valley and the area around City of Lincoln Water System (LWS) well field. Arsenic data collected from the Platte River since 1974 were examined using the “weighted regression on time, discharge, and season” model, which compared the streamflow (also referred to as “discharge”), time of year, and season to estimate concentrations of arsenic. Annual mean arsenic concentrations modeled for more than 49 years at the Platte River at Louisville, Nebr., U.S. Geological Survey streamgage (station 06805500), indicated a significant increasing trend. Arsenic concentrations in the Platte River were seasonal, with the highest concentrations being observed during mid- to late summer. When seasonal patterns and streamflow were combined with arsenic concentrations in the Platte River during low streamflow conditions, groundwater contributions, which can have higher arsenic concentrations, make up a larger portion of the streamflow. Arsenic samples were collected from upstream rivers in 2022 and 2023 and were paired to analyze the arsenic contributions at the U.S. Geological Survey streamgage on the Platte River near Ashland, Nebr. (station 06801000), near the City of Lincoln well field. The arsenic concentrations from the streamgage on the Platte River near Ashland, Nebr., location, were higher than the U.S. Geological Survey streamgage on the Elkhorn River at Waterloo, Nebr. (station 06800500), and significantly lower than at the U.S. Geological Survey streamgage on the Platte River near Leshara, Nebr.(station 06796500), indicating that the Platte River usually contributes a higher concentration of arsenic than does the Elkhorn River as they join near Ashland, Nebr. During 1991–2023, six groundwater monitoring wells were analyzed to identify trends in arsenic concentrations. Two of the six wells had a positive trend during the 33-year period. One monitoring well did not reveal a long-term trend during this period but showed a trend during 2019–23, correlating to a period when the island in the middle of the Platte River was connected to the east bank of the river when manganese reducing conditions were present and groundwater levels were declining in the well. Across all wells the oxidation and reduction (redox) condition during the time of sampling was assessed. Mixed anoxic and (or) oxic redox condition was the most common redox process and the highest sampled arsenic concentrations in monitoring wells were observed in anoxic conditions driven by manganese reduction. Groundwater arsenic concentrations had seasonal variation around the City of Lincoln well field, with higher arsenic concentrations tending to be further south in comparison to samples collected further north. Isotope samples were collected and analyzed in surface water and groundwater around the LWS well field. The samples indicate that the proportion of surface water present in the LWS production wells can be higher in the spring and lower in the summer. With higher arsenic concentrations observed in the stream water during the summer period, the LWS source water can be affected by these elevated arsenic concentrations even though the proportion of surface water is lower.

Released April 09, 2026 12:25 EST

2026, Circular 1564

Sara Ernst

The 2024 annual report of the U.S. Geological Survey Woods Hole Coastal and Marine Science Center highlights accomplishments of 2024, includes a list of 2024 publications, and summarizes the work of the center, as well as the work of each of its science groups. This product allows readers to gain a general understanding of the focus areas of the center’s scientific research and learn more about specific projects and progress made throughout 2024, all while enjoying photographs taken in various environments and laboratories, and applicable maps and figures.

Released April 09, 2026 11:16 EST

2026, Scientific Investigations Report 2026-5134

Hayden A. Lockmiller, Victoria (Tori) Byers, Lisa R. Fogarty

Lake St. Clair Metropark Beach in Michigan has a history of closures because of elevated Escherichia coli (E. coli) concentrations in its recreational waters. To reduce closures, restoration projects were implemented in 2021 to deter waterfowl from congregating on the beach. In this study, the U.S. Geological Survey, in cooperation with the Michigan Department of the Environment, Great Lakes, and Energy and in collaboration with Huron-Clinton Metroparks and the Macomb County Health Department, monitored E. coli from 2022–23 in surface water, shallow groundwater, and sediment at Lake St. Clair Metropark Beach. Results were compared to data from a prerestoration (2018–19) study. A significant decrease in daily geometric mean E. coli concentrations in surface water was observed postrestoration, but the number of high concentration events increased. This resulted in more frequent beach closures postrestoration. Surface-sediment E. coli concentrations significantly decreased after restoration, and waterfowl populations generally decreased from 2021 to 2023, suggesting that the deterrence measures could be influencing E. coli concentrations in surface sediments and surface water. Groundwater E. coli concentrations were orders of magnitude higher than those in surface water and revealed no change correlated with restoration. Seepage measurements indicated that groundwater occasionally discharges into surface water, potentially providing a transport mechanism for E. coli to reach the lake. Continued monitoring and consideration of environmental factors could help to better understand the beach system.

Released April 09, 2026 09:40 EST

2026, Preprint

Phil 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.

Released April 09, 2026 09:24 EST

2026, Ecology and Evolution (16)

Fiona Lunt, C. Lane Scher, Riley Olivia Mummah, David A.W. Miller

Data integration methods aim to improve species distribution estimates by incorporating multiple sources of uncertainty across datasets. Two major sources of uncertainty are: (1) variation in sampling effort across space and within datasets, and (2) variation in reliability associated with data collection protocols or timing among datasets. Our goal was to evaluate how different approaches to address these uncertainties influence predictive performance of integrated models. We modeled distributions of four bird species using three datasets that differed in sampling design. We examined three strategies to reduce uncertainty: (1) filtering data, (2) incorporating functions that account for uncertainty in observation models, and (3) varying how datasets are integrated into a single estimate. We first examine methods to account for variable effort in observations, focusing on both spatial differences in sampling intensity and effort given to a single observation record. We then examine approaches to account for data sets with differing reliability. Sampling effort was best addressed through conservative filtering, including spatial thinning and excluding observations with highly variable effort. Next, we considered how to account for potential false positive detections—due to either misidentification or changes in distributions. We found that treating less reliable data as a covariate, an approach previously suggested for data integration that can greatly speed up model fitting, performed well. Other effective approaches included directly modeling false positive rates and complete exclusion of less reliable data sets. Our results provide insights into best practices in integrated modeling for handling uncertainty in integrated models. We demonstrate the flexible options available when using integrated models to address uncertainty.

Released April 09, 2026 08:40 EST

2026, PLoS ONE (21)

William Michael Janousek, Gavin G. Cotterill, Olivia J. Lobo, Eric K. Cole, Sarah R. Dewey, Tabitha A. Graves

Understanding social aggregation patterns in ungulate herds is essential for gaining behavioral insights, optimizing resource use, reducing human-wildlife conflict, and managing disease risk. As chronic wasting disease is the preeminent disease-related threat to cervid populations in North America, knowledge of contact between individuals and spatiotemporal patterns of aggregation provides opportunity to understand and potentially reduce disease risk while supporting sustainable population sizes. Herd density metrics, derived from global positioning system (GPS) data, can be used to inform management decisions. To effectively compare aggregation behavior within and between herds, aggregation metrics must be accurate. However, the consistency of metrics across different GPS collar sample sizes remains unclear and robust studies of big game require understanding how these factors may vary in different contexts. We examined the minimum sample size necessary for reliable calculations of three aggregation metrics: pairwise inter-animal distances, daily proximity rates, and kernel density estimate (KDE) areas. We used GPS collar data from the Jackson and West Green River elk herds (Cervus canadensis) in western Wyoming, USA, that differ in herd size and group structure (single versus multiple sub-groups), representing common practical contexts. Elk locations were acquired for the Jackson herd between 2016 and 2019 and from 2005 to 2010 for the West Green River herd. Herd-specific characteristics substantially influence the sample size necessary for accurate density measurements. As predicted, larger herds with many groups require more GPS collars than small herds with fewer groups. The sample size needed to accurately estimate aggregation varies by metric, with KDE areas, useful for indexing environmentally transmitted disease risk, generally requiring fewer samples, especially in high-density contexts. The required sample size also varies with seasonal changes in density. During periods of highest density, similar sample sizes are required to estimate inter-animal distances and proximity rates regardless of herd characteristics. Our results have implications for costs associated with studying big game herds, indicating fewer collars may be sufficient in some cases. These insights can aid researchers and managers in determining the appropriate number of GPS collars required for effective herd monitoring and informing relevant aggregation metrics for their management goals.

Released April 09, 2026 08:12 EST

2026, Scientific Reports (16)

Ana Lomashvili, Kristin Rammelkamp, Protim Bhattacharjee, Olivier Gasnault, Elise Clavé, Christoph H. Egerland, Susanne Schröder, Travis S.J. Gabriel, Ari Essunfeld, Stéphane Le Mouélic, Begüm Demir

Since the Mars Science Laboratory landed in 2012, the ChemCam instrument aboard the rover has collected in-situ laser-induced breakdown spectroscopy (LIBS) data and context images along more than 35 km of the Gale Crater traverse, providing valuable observations including diagenetic features such as light-toned veins. These veins are of particular scientific interest because they are interpreted as indicators of past fluid circulation on Mars and provide insights into the evolution of habitability on Mars. Their identification, however, currently relies on manual visual inspection of Remote Micro Imager (RMI) images, a process that is time-consuming and sensitive to differences in human interpretation. To address this issue, in this paper we introduce a novel pixel-level labeled, multimodal dataset of ChemCam observations specifically tailored for vein detection, along with customized U-Net models to integrate both textural (RMI) and chemical (LIBS) modalities. To further ensure trustworthy scientific use, we incorporate the Learn-Then-Test (LTT) framework to provide statistical control of the false discovery rate without requiring model retraining. The experimental results demonstrate that the proposed customized U-Net models trained on the developed dataset, combined with risk-controlled prediction, increases the efficiency of pixel-level vein identification through automation and produces statistically reliable predictions for multimodal ChemCam data.

Released April 09, 2026 07:47 EST

2026, Circular 1566

Yellowstone Volcano Observatory

The Yellowstone Volcano Observatory (YVO) monitors volcanic and hydrothermal activity associated with the Yellowstone magmatic system, carries out research into magmatic processes occurring beneath Yellowstone Caldera, and issues timely warnings and guidance related to potential future geologic hazards. YVO is a collaborative consortium that includes the U.S. Geological Survey (USGS), Yellowstone National Park, University of Utah, University of Wyoming, Montana State University, EarthScope Consortium, Wyoming State Geological Survey, Montana Bureau of Mines and Geology, and Idaho Geological Survey. The USGS component of YVO also has the operational responsibility for monitoring volcanic activity in the Intermountain West of the United States, including Arizona, New Mexico, Utah, and Colorado. This report summarizes the activities and findings of YVO during the year 2024, focusing on the Yellowstone volcanic system.