Released January 27, 2026 15:30 EST

2026, Professional Paper 1890-N

Marissa E. Mnich, Christopher D. Condit

A detailed characterization of the >3,000 square kilometer (km²) Springerville volcanic field, located on the southern tip of the Colorado Plateau in Arizona, United States, with its more than 501 volcanic units and widely distributed >420 cinder cones and lava flows, provides constraints toward an integrated petrogenetic model for the field. Large-volume effusive tholeiitic eruptions at 2–1.5 mega-annum (Ma) transitioned to more numerous, smaller volume alkali olivine basalt (AOB) events at 1.5–1.0 Ma, with increasing abundances of evolved alkalic rocks (EARs), and a final 1.0–0.3 Ma period dominated by smaller volume, more explosive alkalic eruptions.

Early large-scale melting in a relatively enriched lithospheric mantle (EM) source generated large-volume effusive tholeiitic magmas. Depths of tholeiite magma generation average about 90 kilometers (km) across the field, but depths for individual units decreased southward, consistent with lithospheric thinning toward the Colorado Plateau margin. Early and middle-stage transitional basalts, alkali olivine basalts (AOBs), and basanites originate from a progressively deeper (>100 km) region in a prevalent mantle (PREMA)-like asthenospheric source produced by increasingly smaller degrees of melting, as low as about 2 percent. The chemical signature of the basanites is consistent with small degrees of melting in a carbonated, asthenospheric source to depths of about 140 km. As heat waned, the last phase of volcanism was dominated by more explosive EARs derived at shallower lithospheric pressures but that have isotopic and trace element similarities to the deeper asthenospheric magmas. This suggests mixing between deeper basanitic and shallower tholeiitic magmas. With waning heat, eruptions became more localized along alignments, likely related to boundaries between blocks of Proterozoic crust with differing properties that affected magma ascent.

The petrogenetic patterns are consistent with a variety of processes. Basin and Range Province extension, melting, and heat-induced weakening progressively eroded the Colorado Plateau’s thicker lithosphere, giving rise to relatively high degrees of partial melting from shallower (<90 km) sources that produced the early large volumes of effusive tholeiitic magma. Possible lithosphere delamination and removal, and a resulting steep boundary step with the asthenosphere, may have given rise to mantle edge convection, resulting in the mixing of basanitic and tholeiitic magmas. With plate motion, shear-driven upwelling likely gave an eastward component to convection, affecting relative amounts of melting at the field’s western boundary and corresponding with volcanism moving eastward during times of peak vent production by eruption of alkalic lavas, giving rise to more explosive and lower volume eruptions.

Released January 27, 2026 11:56 EST

2026, Open-File Report 2025-1024

Rebecca K. Mann, Michael C. Duniway, Jo Ellen Hinck

Mining uranium from breccia-pipe deposits in the greater Grand Canyon region has occurred since the mid-1900s. However, possible ecosystem contamination with harmful levels of radionuclides may have occurred due to mining activities in the 21st century. In response, a 20-year Federal moratorium on new mining claims in the Grand Canyon watershed was initiated in 2012, to allow time to evaluate the potential effects of uranium exploration and mining on human health, wildlife, and water resources. This moratorium, nor the 2023 designation of the “Baaj Nwaavjo I’tah Kukveni–Ancestral Footprints of the Grand Canyon National Monument,” precludes operation or development of mining claims predating 2012.

Vegetation is a core ecosystem component that may be affected by uranium mining (for instance, through uptake and storage of radionuclides from the air or soil) or may act as a vector of exposure to wildlife, livestock, and humans (for instance, via their consumption of contaminated plant tissues). To provide baseline information about the plant communities associated with uranium mines in the Grand Canyon region, the U.S. Geological Survey surveyed an approximately 200-meter-wide buffer surrounding four breccia-pipe deposits, each in a unique stage of mine development, and at one reference area (a livestock water tank) that underwent ground disturbance but contains no mineral deposits. We sectioned the buffer zones into 0.65–4.52 hectare plots, within which we (1) inventoried all plant species, (2) measured percent cover of plant species, plant functional groups, and ground surface types (dark cyanobacteria, lichen, moss, bedrock, rock, embedded litter, duff, plant bases, and bare soil) using line-point intercept, and (3) measured length and frequency of gaps between perennial plant canopies using canopy gap intercept. We found that plant composition at the mines and the reference area differed from one another but were all characteristic of expected regional vegetation patterns. We provide this data summary as potential baseline information for future research and management efforts.

Released January 27, 2026 08:18 EST

2026, Open-File Report 2026-1059

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

Executive Summary 

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

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

Released January 26, 2026 08:50 EST

2026, Fact Sheet 2025-3058

Judith Z. Drexler, Jake F. Weltzin

Introduction 

The San Francisco Bay and Sacramento–San Joaquin Delta Estuary (hereafter, Bay-Delta) is the largest estuary on the West Coast of the United States. The Bay- Delta covers more than 1,600 square miles and drains a watershed of more than 75,000 square miles, which is greater than 40 percent of California. The region surrounding the Bay- Delta is home to about 10 million people, and its habitats (fig. 1) support more than 800 plant and animal species. The waterways of the Bay- Delta are the central hub of California’s extensive freshwater delivery system, supplying water to more than 27 million Californians and 4 million acres of farmland in the Central Valley.

The U.S. Geological Survey (USGS) is the primary science agency of the U.S. Department of the Interior, providing a broad range of Earth, water, biological, and mapping data and expertise to inform natural resource management across the country. This fact sheet focuses on research conducted by the USGS in the Bay- Delta region, mostly within the past 5 years. The fact sheet is organized across five major themes: water flow and water quality, fish and wildlife, wetland restoration, invasive species, and hazards. In each of these areas, the USGS works closely with Federal, State, and local agencies; academia; and non- governmental organizations to provide objective, evidence-based science. The data and knowledge gained from USGS research and monitoring are publicly available and are routinely used by agencies, including the U.S. Bureau of Reclamation, U.S. Army Corps of Engineers, U.S. Fish and Wildlife Service, California Department of Water Resources, California Department of Fish and Wildlife, California State Water Resources Control Board, Delta Stewardship Council, and many additional organizations and stakeholders to inform management of the Bay- Delta region and its resources.

Released January 23, 2026 07:45 EST

2026, Journal of Soil and Water Conservation

Brittany Rose Hanrahan, Matthew W. Diebel, Rebecca B. Carvin, Edward G. Dobrowolski, Matthew John Hardebeck, Andrew Jason Kowalczk, Chad A. Toussant, Matthew J. Komiskey

The Laurentian Great Lakes are a vital freshwater resource in the United States, and nonpoint source (NPS) nutrient pollution, specifically phosphorus (P) and nitrogen (N), from agricultural land use continues to negatively impact water quality throughout the Great Lakes basin. One focus of the Great Lakes Restoration Initiative (GLRI), a mechanism to coordinate conservation efforts in the Great Lakes that began in 2010, is reducing NPS nutrient pollution through the implementation of conservation practices in priority watersheds (Genesee River, Fox River, Maumee River, and Saginaw River). As part of GLRI efforts, the objective of the study presented here was to evaluate the effects of conservation implementation, specifically increasing vegetative cover on fields and in primary flowpaths through perennial or cover crop planting and grassed waterways, on surface-runoff water quality at 12 agricultural fields (six paired and six unpaired) located in priority watersheds. We determined the percentage difference in mean event response variables between the periods before and after conservation implementation at individual sites, describing patterns across sites to synthesize lessons learned from these GLRI evaluations. Generally, we found that mean event flow-weighted concentration (FWC) and yield (kilograms per hectare) decreased for suspended sediment (SS) and nitrate (NO₃^–-N) across many sites. Mean event FWC and yield for total P (TP) showed mixed results across sites, while mean event FWC and yield for orthophosphate generally increased across sites. These results indicate that perennial or cover crop planting and grassed waterways effectively reduce SS and NO₃^–-N losses in surface runoff from agricultural fields, but mitigating TP and dissolved P losses remains a challenge.

Released January 22, 2026 09:04 EST

2026, Palaeogeography, Palaeoclimatology, Palaeoecology (686)

John A. Barron, Thomas J. DeVries, Jason J. Coenen

Diatom records from the East Pisco Basin (EPB) of southern Peru and offshore Ocean Drilling Program (ODP) Hole 682 A reveal stepwise increases in the primary productivity of the Humboldt Current during the middle and late Miocene. Although diatoms are present back through the late middle Eocene, successively enhanced diatom production occurs during the Miocene in four steps. The first step between 14.2 and 13.8 Ma marks the onset of diatom deposition in the Pisco-0 sequence. A second step at ∼12.9 Ma coincides with a major drop in global sea level during which diatom deposition ceased in the EPB but continued in ODP 682 A as well as in Chile, Ecuador, and Colombia. Beginning at 10.4 to 10.0 Ma, a major increase in diatom sedimentation rates in both the EPB and in ODP 682 A signals an abrupt intensification of the Humboldt Current productivity, which may be related to the closure of the Central American Seaway to deeper water circulation. A fourth step coincides with the transgressive deposition of the diatom-rich Pisco-2 sequence at 8.4 Ma and is enhanced by the onset of the global Late Miocene Biogenic Bloom (8.0–4.5 Ma). Between 7 and 6 Ma, common subtropical diatoms in ODP 682 A contrast with high diatom deposition rates in the EPB, suggesting variable El Niño-La Niña conditions. During this time, fossil vegetation evidence from southern Peru supports the presence of seasonal periods of enhanced rainfall typical of El Niño conditions. The Miocene trend of increasing sea mammal diversity in the EPB parallels these steps, further supporting stepwise enhancement of primary productivity in the Humboldt Current during the Miocene.

Released January 22, 2026 08:44 EST

2026, Remote Sensing (18)

Nicholas Brimhall, Kelvyn K. Bladen, Tom Kerby, Carl J. Legleiter, Cameron Swapp, Hannah Fluckiger, Julie E Bahr, Makenna Roberts, Kaden Hart, Christina L. Stegman, Brennan Bean, Kevin Moon

Remote sensing enables large-scale, image-based assessments of river dynamics, offering new opportunities for hydrological monitoring. We present a publicly available dataset consisting of 281,024 satellite and aerial images of U.S. rivers, constructed using an Application Programming Interface (API) and the U.S. Geological Survey’s National Hydrography Dataset. The dataset includes images, primary keys, and ancillary geospatial information. We use a manually labeled subset of the images to train models for detecting rapids, defined as areas where high velocity and turbulence lead to a wavy, rough, or even broken water surface visible in the imagery. To demonstrate the utility of this dataset, we develop an image segmentation model to identify rivers within images. This model achieved a mean test intersection-over-union (𝐼𝑜𝑈) of 0.57, with performance rising to an actual 𝐼𝑜𝑈 of 0.89 on the subset of predictions with high confidence (predicted 𝐼𝑜𝑈 > 0.9). Following this initial segmentation of river channels within the images, we trained several convolutional neural network (CNN) architectures to classify the presence or absence of rapids. Our selected model reached an accuracy and F1 score of 0.93, indicating strong performance for the classification of rapids that could support consistent, efficient inventory and monitoring of rapids. These data provide new resources for recreation planning, habitat assessment, and discharge estimation. Overall, the dataset and tools offer a foundation for scalable, automated identification of geomorphic features to support riverine science and resource management.

Released January 21, 2026 19:43 EST

2026, Scientific Investigations Map 3542

Gregory J. Walsh, Sean P. Regan, Phillip S. Geer, Arthur J. Merschat, Kaitlyn A. Suarez, Ryan J. McAleer, Matt S. Walton, Jr., E. Allen Crider, Jr.

The bedrock geology of the 7.5-minute Eagle Lake quadrangle, Essex County, New York, consists of deformed and metamorphosed Mesoproterozoic gneisses of the Adirondack Highlands unconformably overlain by weakly deformed lower Paleozoic sedimentary rocks of the Champlain Valley. The Mesoproterozoic rocks occur on the eastern edge of the Adirondack Highlands and represent an extension of the Grenville Province of Laurentia. Granulite facies Mesoproterozoic paragneiss, marble, and amphibolite hosted the emplacement of an anorthosite-mangerite-charnockite-granite (AMCG) suite, now exposed mostly as orthogneiss, at approximately 1.18–1.15 giga-annum (Ga, billion years before present). The earliest of four phases of deformation (D1) predated AMCG magmatism and is characterized by gneissosity, rarely preserved F1 isoclinal folds, and migmatite in the paragneiss host rocks. A sample of hornblende quartz syenite from the AMCG suite, collected from an abandoned railroad cut on Old Furnace Road, yielded a U-Pb zircon age of 1,149±10 million years before present. D2 deformation produced a composite penetrative gneissosity, migmatite, and isoclinal F2 folds. Towards the end of D2, felsic magmatism (including the regionally extensive Lyon Mountain Granite Gneiss, abbreviated “LMG”) spread by penetrative migration as semiconcordant alkali feldspar granite sheets subparallel to S2 into the previously deformed lithologies. The LMG crystallized at approximately 1.15 to 1.14 Ga and displays synkinematic F2 folds thus constraining the time of D2 deformation. Exhumation of the Marcy anorthosite began during D3 along a mylonitic extensional detachment, as a type of core complex. Protracted D3 produced F3 folds exhibited in regional domes and basins, such as the Hammondville antiform, reactivation of the S2 foliation, partial melting, metamorphism, metasomatism, iron ore remobilization, and intrusion of magnetite-bearing pegmatite both as layer-parallel sills and crosscutting dikes. D4 created NE- and NW-trending boudinage, local high-grade ductile shear zones, and crosscutting granitic pegmatite dikes. Kilometer (km)-scale lineaments readily observed in lidar data are Ediacaran mafic dikes and Phanerozoic brittle faults. Lower Paleozoic rocks are part of the Early Cambrian to Late Ordovician great American carbonate bank on the ancient margin of Laurentia. The Potsdam Sandstone preserves the Cambrian stratigraphy in outliers above the Great Unconformity. The Paleozoic rocks are weakly folded and block faulted. Parts of the quadrangle are covered by undifferentiated glacial deposits, but much of the quadrangle contains only a variably thick, veneer of unmapped glacial till over significant areas of exposed bedrock. The map also shows waste rock piles and locations of historical mining operations. This study was undertaken to improve our understanding of the bedrock geology in the Adirondack Highlands, establish a modern framework for 1:24,000-scale bedrock geologic mapping in the Adirondack Mountains, and provide a modern context for historical mines. This Scientific Investigations Map of the Eagle Lake 7.5-minute quadrangle consists of a map sheet, an explanatory pamphlet, and a geographic information system database that includes bedrock geologic units, faults, outcrops, and structural geologic information. The map sheet includes a bedrock geologic map, a correlation of map units, a description of map units, an explanation of map symbols, and two cross sections. The explanatory pamphlet includes a discussion of the geology.

Released January 21, 2026 09:14 EST

2026, Ecology and Evolution (16)

C. M. Costello, L. L. Roberts, D. D. Bjornlie, M. D. Cameron, J. G. Clapp, Mark Haroldson, G. V. Hilderbrand, K. Joly, W. Kasworm, J. M. Nicholson, T. Radandt, M. S. Sorum, J. E. Teisberg, Frank T. van Manen, M. A. Vinks

Timing of grizzly bear (Ursus arctos) parturition during hibernation has been explained by ancestral traits (delayed implantation, altricial young, obligate maternal denning), but the ultimate driver underlying precise timing has not been fully explored. Capitalizing on an observed latitudinal increase in denning duration among four populations in interior North America, we tested two alternative hypotheses. First, that birth timing results from a physiological cue that synchronizes implantation with the onset of hibernation, allowing females to forgo reproduction should they lack adequate fat stores. Alternatively, that parturition is optimally timed relative to den exit to balance an energetic tradeoff between minimizing lactation time to protect the mother and maximizing developmental time to increase cub survival. Using parturition dates previously predicted from accelerometer data (27 Dec–28 Feb), we classified 115 females according to apparent litter survival when first visually observed after den exit: 57% successful (with cubs), 22% unsuccessful (alone), and 21% unknown (not observed). The number of days between birth and den exit showed no association with latitude (p = 0.29). It averaged 103 days among successful females but only 77 days among unsuccessful females (p < 0.001) owing to later births and earlier exit. With each increasing degree of latitude, birth date increased by 1.0 and number of days between den entry and birth increased by 2.5 (p < 0.001). Implantation dates were not centered on den entry dates (p < 0.001). These results supported the energetic tradeoff hypothesis and suggested natural selection has favored a consistent number of days between parturition and den exit under average body conditions and shifts toward later or earlier births for females with lower or higher levels of bodily stored energy, respectively. This flexible tradeoff may support resilience to climate change and present a possible mechanism explaining reduced natality and cub survival in high-density populations.

Released January 21, 2026 07:00 EST

2026, Open-File Report 2025-1057

Suellen Lynn, Alexandra Houston, Barbara E. Kus, Shannon M. Mendia

Executive Summary 

The purpose of this report is to provide the Marine Corps with a summary of abundance, breeding activity, demography, and habitat use of endangered Least Bell’s Vireos (Vireo bellii pusillus) at Marine Corps Base Camp Pendleton, California (MCBCP or Base). The report presents results of vireo surveys and monitoring in 2024 and summarizes a subset of data collected from 2020 through 2024. Surveys for the Least Bell's Vireo were completed at MCBCP between April 4 and July 9, 2024. Core survey areas and a subset of non-core areas in drainages containing riparian habitat suitable for vireos were surveyed two to four times. We detected 542 territorial male vireos and 17 transient vireos in core survey areas. An additional 102 territorial male vireos and 2 transients were detected in non-core survey areas. Transient vireos were detected on 5 of the 10 drainages/sites surveyed (core and non-core areas). In core survey areas, 87 percent of vireo territories were on the four most populated drainages, with the Santa Margarita River containing 67 percent of all territories in core areas surveyed on Base. In core areas, 77 percent of male vireos were confirmed as paired; 76 percent of male vireos in non-core areas were confirmed as paired.

The number of documented Least Bell’s Vireo territories in core survey areas on MCBCP decreased 3 percent from 2023. In five core survey area drainages, the number of territories increased by at least two, and in two core survey area drainages, the Santa Margarita River and Las Flores Creek, the number of vireo territories decreased by at least nine between 2023 and 2024. The number of vireo territories at Marine Corps Air Station, Camp Pendleton did not change from 2023 to 2024. The proportion of surveys during which Brown-headed Cowbirds (Molothrus ater) were detected decreased to 0.03 from a peak of 0.45 in 2022. Cowbirds were detected in April and June in 2024.

Most core-area vireos (58 percent, including transients) used mixed willow (Salix spp.) riparian habitat. An additional 9 percent of birds occupied willow habitat co-dominated by Western sycamores (Platanus racemosa). Riparian scrub dominated by mule fat (Baccharis salicifolia), sandbar willow (S. exigua), or blue elderberry (Sambucus mexicana) was used by 33 percent of vireos. Habitat dominated by non-native vegetation was used by 1 percent of vireos.

Since 2020, the number of vireos detected in each of the non-core survey groups was greater than expected, based on the change in vireo numbers in core survey areas. Although, the number of vireo territories on Base decreased from 2020–24, from approximately 1,224 to approximately 960, the trend in vireo territory numbers on Base since 2005 has been positive.

In 2019, MCBCP began operating an artificial seep along the Santa Margarita River; then, in 2021, two additional artificial seeps became operational. The artificial seeps pumped water to the surface during daylight hours starting in mid-April and ending in August each year and were designed to increase the amount of surface water to enhance Southwestern Willow Flycatcher (Empidonax traillii extimus) breeding habitat. Although this enhancement was designed to benefit flycatchers, few flycatchers have inhabited MCBCP, including the seep areas, within the past several years; therefore, vireos were selected as a surrogate species to determine effects of the habitat enhancement. This report presents the fifth year of annual monitoring and analyses summarizing all 5 years of vireo and vegetation response to the artificial seeps.

In 2020, we established four study sites along the Santa Margarita River, two surrounding and extending downstream from existing and proposed seep pumps at the Old Treatment Ponds and along Pump Road and two Reference sites in similar habitat downstream from the Seep sites. Seep pumps began operating at the Old Treatment Ponds in 2020 and along Pump Road in 2021. In 2023, seep pumps at the Pump Road Seep site did not function, and we recategorized that study site as Intermediate. We sampled vegetation at Seep, Intermediate, and Reference sites to determine the effects of surface-water enhancement by seep pumps. In 2024, vegetation cover was highest near the ground and decreased with increasing height. Woody vegetation made up most of the cover at all height categories. We determined that Seep and Intermediate sites differed from each other in addition to differing from Reference sites, which likely is, in part, because seep-pump operation at the Intermediate site was inconsistent compared to the Seep site. Soil saturation in 2024 was high at the Intermediate site and was associated with high native herbaceous cover and low non-native herbaceous cover. Sites differed, with the Intermediate site having more upper canopy cover in general, the Seep site having more low woody cover, and the Reference sites having more mid-canopy non-native vegetation cover.

Soil saturation significantly increased from 2020 through 2024 at the Seep site and was significantly higher at Seep and Intermediate sites than at their paired Reference sites in all years. Soil saturation likely was increased by the supplemental surface water at the Seep site. However, soil saturation at the Intermediate site was not clearly associated with seep pumps but likely affected by soil saturation at the site before seep-pump installation and flooding from high precipitation. Canopy height increased at the Intermediate site from 2020 through 2024 and increased with increasing soil saturation at the Intermediate and Reference sites. The canopy at the Seep site was shorter than at the Intermediate and Reference sites and decreased from 2020 through 2024 because tall trees were damaged and killed by shothole borer beetles (Euwallacea spp.).

We used Redundancy Analysis to discover associations among vegetation types, plant species, and other environmental variables (soil saturation, site, precipitation, and seep operation, defined as the site and year seep pumps were operating). These associations explained less than 15 percent of the variability in the vegetation, with the remaining 85 percent of variation unexplained. Generally, as soil saturation increased, understory vegetation increased and non-native cover decreased in the mid-and upper canopy. Non-native herbaceous plant species decreased in wetter soil.

The Seep site was characterized by more understory and less canopy, contrasting with the Intermediate site, which was characterized by less understory and more higher canopy cover. The addition of surface water via seep pumps or precipitation was associated with more vegetation near the ground. Higher early winter precipitation was associated with taller canopy and more woody vegetation in the upper canopy. We also created a Redundancy Analysis model isolating the components of Southwestern Willow Flycatcher habitat, as identified by Howell and others (2018). In this model, increased soil saturation resulted in increased cover of stinging nettle (Urtica dioica) and black willow (Salix gooddingii) below 3 meters (m), total cover 3–6 m, and black willow above 6 m. Cover of poison hemlock (Conium maculatum) and stinging nettle below 3 m was higher at the Seep site and lower at the Intermediate site.

Vireo territory density among the Seep, Intermediate, and Reference sites was similar before the seep pumps were installed. However, vireo territory density at Seep and Intermediate sites combined was significantly higher than at Reference sites after the seep pumps were installed.

We banded and resighted color banded vireos as part of a long-term evaluation of vireo survival, site fidelity, between-year movement, and the effect of surface-water enhancement on vireo return rate and between-year movement. We banded 164 Least Bell's Vireo nestlings during the 2024 season.

In 2024, we resighted 31 Least Bell's Vireos on Base that had been banded before the 2024 breeding season, and we were able to identify 25 of them. Of the 25 that we could identify, 24 were banded on Base and 1 was originally banded on the San Luis Rey River. Adult birds of known age ranged from 1 to 9 years old.

Base-wide survival of vireos was affected by sex, age, and year. Males had significantly higher annual survival than females (60 percent versus 47 percent, respectively). Adults had higher annual survival than first-year vireos (61 percent versus 11 percent, respectively). The return rate of adult vireos to Seep, Intermediate, or Reference sites was not affected by the original banding site (Seep versus Intermediate versus Reference).

Most returning adult vireos, predominantly males, showed strong between-year site fidelity. Of the adults present in 2023, 92 percent (all males) returned in 2024 to within 100 m of their previous territory. The average between-year movement for returning adult vireos was 0.4±0.03 kilometers (km). The average movement of first-year vireos detected in 2024 that fledged from a known nest on MCBCP in 2023 was 2.4±3.1 km.

We monitored 47 Least Bell's Vireo pairs to evaluate the effects of surface-water enhancement on nest success and breeding productivity. Breeding productivity in 2024 was similar among Seep, Intermediate, and Reference sites (2.8, 3.0, and 3.0 young fledged per pair, respectively), and the percentage of pairs that fledged at least one young was not significantly different among sites (83, 91, and 96 percent, respectively). According to the best model, daily nest survival from 2020–24 was not related to site. Other measures of breeding productivity were also similar among Seep, Intermediate, and Reference site pairs.

Between 2020 and 2024, the number of vireo fledglings produced per pair increased with increasing native herbaceous cover under 3 m and decreasing cover of all herbaceous vegetation under 5 m and was not affected by precipitation, site, or seep operation. The number of vireo fledglings produced per egg was lower at the Seep and Intermediate sites than at the Reference sites and increased with decreasing late winter precipitation, cover of poison hemlock, black mustard, non-native vegetation above 2 m, and all vegetation over 2 m. Vireo pairs at Seep and Intermediate sites were less likely to fledge young than vireo pairs at Reference sites. All vireo pairs were more likely to fledge young with less cover of poison hemlock and more cover of poison oak.

From 2020 through 2024, vireos placed their nests in 24 plant species. The most used plants in all years were willows, mostly red (S. laevigata), or arroyo (S. lasiolepis). The fate of a vireo nest (whether it successfully fledged young or not) was not affected by placement in native or non-native vegetation, by site, or by year, but nests were more likely to be successful if they were placed in woody plants than in herbaceous plants. Successful nests were placed higher in the host plant and farther from the outer edge of the nest clump than unsuccessful nests.

Released January 20, 2026 15:17 EST

2026, Scientific Investigations Report 2025-5102

Kevin K. Johnson

The Cook County Precipitation Network is a set of 25 precipitation gages established within Cook County, Illinois, on approximately a 5- to 7-mile square grid and used by the U.S. Army Corps of Engineers to help account for diversions of water from Lake Michigan to the State of Illinois. The transition from the precipitation gage network operated by the Illinois State Water Survey to the precipitation gage network operated by the U.S. Geological Survey (USGS) was compared for periods of overlapping data. This transition took place from May through September during the 2019 water year. The USGS was able to establish replacement precipitation gages at 17 of the 25 sites by the conclusion of the overlapping operational period.

The double-mass curve method was used to compare the two networks by creating a graph of the cumulated data collected by the Illinois State Water Survey and the comparable data collected by the USGS. Breaks in the double-mass curve method are caused by a change in the relation between variables. The eight sites that were installed following the overlapping period have a gap in the recorded data; however, the slope of the line for each of the eight sites is nearly equivalent to the previous data. In general, the cumulated precipitation data from the two networks were similar. Three sites had greater than 8-percent difference in their cumulative data ratios, located at Cicero, Ping Tom Park at Chicago, and South Shore, Ill.

Released January 20, 2026 09:22 EST

2026, Ecotoxciology and Environmental Safety (310)

Lillian M. Pavord, Melissa K. Driessnack, Aaron B. Shiels, Steven Volker, Barnett A. Rattner, Jenifer McIntyre

To restore native biodiversity on island ecosystems containing invasive rodents, partial- and whole-island eradications generally rely on broadcast baiting with anticoagulant rodenticides (ARs). This approach can result in bait pellets entering aquatic environments, raising concerns about effects to non-target fish. Salmonids are a dominant group of fishes on many temperate islands targeted for rodent eradication, and AR toxicity data for salmonids are limited. Our goal was to determine if coho salmon (Oncorhynchus kisutch) are susceptible to coagulopathy and death via exposure to commonly used ARs. We assessed risk of ARs to coho using dose-response curves generated through intraperitoneal injections after determining that coho would not directly ingest the AR baits. Median lethal doses (96-h LD₅₀) estimated using 100 % corn oil carrier were 85.7 µg/g for brodifacoum and 54.0 µg/g for diphacinone. Acetone (30–41 %), used to dissolve ARs in corn oil, reduced the toxicity of diphacinone (LD₅₀ = 102.3 µg/g, p < 0.001) but not brodifacoum (LD₅₀ = 73.3 µg/g, p = 0.126) indicating that solvent choice can influence toxicity outcomes. Behavioral changes and onset of mortality differed between the two ARs, with diphacinone acting more rapidly. Tissue analysis supported a difference in toxicokinetics between the two ARs, with significant decreases in liver and muscle residues for diphacinone but not brodifacoum. Sublethal brodifacoum exposure (53.9 µg/g; LD₁₃) impaired blood clotting at 72- and 96- h but returned to baseline by 120 h. No clotting impairment was observed up to 144 h after diphacinone exposure (45.5 µg/g; LD₄), suggesting a non-coagulopathy mode of action. These findings will inform risk assessments when considering use of these ARs for rodent management near streams and shorelines and clearly demonstrate that brodifacoum causes coagulopathy in coho.

Released January 19, 2026 08:59 EST

2026, Quaternary Research

Karissa Cordero, Nathan Brown, Lauren N. Harrison, Shaul Hurwitz

Hydrothermal explosions are a significant geological hazard in some active volcanic systems; however, the timing and triggering mechanisms of these explosions are poorly constrained. This study applies luminescence dating techniques to hydrothermal explosion deposits in the Yellowstone Plateau volcanic field to constrain explosion chronologies and evaluate potential triggering mechanisms. We tested four luminescence dating techniques: K-feldspar post-infrared infrared stimulated luminescence (pIRIR225), quartz blue light optically stimulated luminescence (BLOSL), quartz blue thermoluminescence (BTL), and quartz red thermoluminescence (RTL). The pIRIR225 and RTL protocols produce consistent age estimates that agree with independent radiocarbon ages and with the timing of the Pinedale deglaciation. This study focuses on two craters, Mary Bay, along the northern shore of Yellowstone Lake, and Pocket Basin in Lower Geyser Basin. The mean pIRIR225 ages from Mary Bay deposits (11.99 ± 0.68 ka) agree with previous radiocarbon constraints. The mean pIRIR225 results from Pocket Basin deposits (13.44 ± 1.06 ka) suggest a history of explosion following Pinedale deglaciation, followed by recent hydrothermal alteration. Luminescence dating techniques are a promising tool for reconstructing the timing of hydrothermal explosions in the Late Pleistocene and Holocene, helping to constrain recurrence intervals of the largest hydrothermal systems, informing risk, and improving hazard assessments.

Released January 16, 2026 09:04 EST

2026, Water Resources Research (62)

Charles R. Kimsal, Enrique R. Vivoni, Osvaldo E. Sala, H. Curtis Monger, Owen P. McKenna

Ephemerally flooded playas are common in the southwestern United States and globally in drylands. Often formed in closed basins, playas are depressions which inundate infrequently from local precipitation and streamflow produced near the playa or from upland areas. Few studies have quantified the hydrologic connectivity between upland catchments and playas using observations. Here, we used rain gauge-corrected precipitation from weather radar and water level measurements in 18 playas of the Chihuahuan Desert to identify precipitation thresholds leading to playa inundation over a 6.4-year period. Geospatial data sets on topography, soil properties, and vegetation cover were employed to determine the controls on inundation. Only 9.4% of all precipitation events above 1 mm led to inundation, with 69.8% of all inundations occurring during the North American monsoon (NAM, July-September). Mean and standard deviations (Std) of runoff ratios at all playas were 2.74 ± 4.08% and 3.29 ± 5.19% for annual and NAM periods. At the annual scale, playa inundation occurred when mean precipitation thresholds of 18.3 ± 7.5 mm (event total) and 12.0 ± 4.5 mm/hr (60-min intensity) were exceeded. Across all playas, inundation occurrence and volume were related most strongly to precipitation metrics and catchment area, with secondary controls of soil and terrain properties. The explanatory power of the derived regressions describing the inundation response across the playas were significantly improved when considering their geological origin. As a result, the inundation response classification system could be applied to ephemeral playas in other arid and semiarid landscapes.

Released January 16, 2026 09:00 EST

2026, Preprint

Emily J. Wilkins, Dieta Hanson, Whitney Boone, Spencer A. Wood, Christian S.L. Crowley, Rudy Schuster

Federally managed public lands and waters receive about 1 billion recreational visits each year. Data on these visitors can aid in guiding policy decisions, managing resources effectively, and communicating the economic contributions of lands and waters. This report explores the methods used by agencies to collect data on aspects of recreational visitor use to Federal lands and waters (apart from visitation numbers, which are the focus of a companion publication). Aspects of recreational visitor use include visitor demographics, recreational activity participation, visitor satisfaction, visitor attitudes and experiences, trip characteristics, and economic contributions. We review practices used to understand aspects of visitor use across seven Federal agencies, revealing similarities such as the use of visitor intercept surveys and coverage of similar topic area, and differences in how survey programs are operationalized and how specific questions on visitor surveys are worded. We also evaluate emerging technologies, such as geolocated social media and mobile device location data, for their potential to aid in understanding aspects of visitor use. This report concludes with potential opportunities to enhance data collection and coordination, ensuring cost-effective data collection and informed decision-making.

Released January 16, 2026 08:04 EST

2026, Journal of Sedimentary Research (96) 1-23

Jason A. Flaum, Katherine L. French, Justin E. Birdwell, Kira K. Timm

Standardizing facies descriptions has proven key to integrating interpretations of depositional processes and environments from sedimentologic observations with geochemistry data for mudstone lithologies. Because of their fine-grained nature, high degree of compaction, and heterogeneous composition, standardizing methods for mudstone descriptions has proven difficult, but it is critical to formulating meaningful interpretations of the processes that govern the accumulation of organic-rich lithologies and their role in both petroleum systems and the global carbon cycle. In this study, we have developed a modified facies classification scheme for mudstone lithologies that incorporates sedimentologic and compositional observation at the hand-sample and thin-section scales with geochemical measurements, including bulk organic and inorganic geochemistry, to characterize these rocks and their variability more completely for improved interpretations of depositional environments during a low-order sea-level transgression. The facies described in this study are of the Cenomanian–Turonian Greenhorn Formation in the USGS #1 Portland Core drilled in Fremont County, Colorado. Strata of the Greenhorn Formation span Oceanic Anoxic Event 2 (OAE-2) and the preceding interval. Lithologies range from organic-rich argillaceous mudstones with varied sedimentary structures to organic-lean, highly bioturbated limestones. Six facies were identified, each differentiated by varied sedimentary structures and geochemical composition. These facies occur in a predictable stratigraphic stacking pattern that represents a low-order sea-level transgression with interpreted depositional environments ranging from terrigenous-dominated pro-delta and muddy continental shelf at the base of the interval to pelagic offshore marine at the top of the Greenhorn Formation. Though the facies are consistent with previous interpretations of depositional environments at this locale in the Cretaceous Western Interior Seaway during the Greenhorn cyclothem, the sedimentary processes governing the accumulation of organic-rich strata that have defined this interval are significantly revised. Variability in the proximity and intensity of bottom currents driven by storms and geostrophic flows were key to the accumulation of each facies, with significant sediment transport occurring even through deposition in the most oxygen-depleted bottom waters. The methodology and interpretations provided here are now being employed to basin-scale predictions of organic enrichment utilizing calibrated petrophysical methods. The approach and results from this study improve understanding of how organic and inorganic carbon was sequestered during perturbations to the global carbon cycle associated with events such as OAE-2.

Released January 15, 2026 10:31 EST

2026, Scientific Investigations Report 2025-5084

Benjamin V. Miller

Oregon Caves National Monument and Preserve in southwestern Oregon is a 4,554-​acre area managed by the National Park Service that is home to several cave systems, including Oregon Caves, which is the longest cave in Oregon, with 3.03 miles of mapped passages. Because of the interconnected nature of karst hydrologic systems, it is critical to understand the areas that can influence water quality and quantity in karst environments. Toward this goal, dye tracing was conducted by the U.S. Geological Survey from 2021 to 2024 to better understand the pathways that karst groundwater follows at Oregon Caves National Monument and Preserve and to delineate recharge areas for two caves, Oregon Caves and Cave Next Door. During the project, eight dye injections were conducted, delineating a 0.51-​square-​mile recharge area for Oregon Caves and a 0.69-​square-​mile recharge area for Cave Next Door. Additionally, the study helped to identify three resurgences associated with Oregon Caves that were previously unknown and showed that the recharge areas for the two caves were distinct from one another. The dye traces also illuminated some unique recharge characteristics of the karst at Oregon Caves, including a high variance in karst groundwater velocities, retention within the karst aquifers, and a significant diffuse-​flow component.

Released January 15, 2026 08:35 EST

2026, Preprint

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

Computation of detailed groundwater flow budgets for subdivisions of Virginia’s Coastal Plain aquifer system has enabled quantification and more thorough understanding of groundwater flow within this important water resource. A zone budget analysis conducted on previously published groundwater models of the Virginia Coastal Plain and Virginia Eastern Shore shows that groundwater conditions vary substantially throughout the Coastal Plain aquifer system due to 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 fundamentally altered groundwater flow from pre-development conditions. Rates of sustainable withdrawal are limited because the downward groundwater flow rate into confined aquifers supplying groundwater is a relatively small portion of the total groundwater water budget for the aquifer system.

Analyses of groundwater budgets from the Virginia Coastal Plain model show that groundwater flow is generally outward from the surficial aquifer to rivers and coastal water bodies 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 currently 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 areas 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 2025, rates of withdrawal substantially lower than the historical maximum have 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 in 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. Downward groundwater flow continues under current conditions, but because vertical flow rates are a function of the difference between water pressure in the upper surficial systems and lower confined units, those rates 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. Estimated groundwater inflow from coastal regions could be contributing to saltwater intrusion, though that was not measured directly in this study.

Analyses of groundwater budgets from the Virginia Eastern Shore peninsula, a geographic region of the Virginia Coastal Plain, show 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 over 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. Currently, estimated storage loss accounts for slightly under 7 percent of withdrawals from the confined aquifer system. The current withdrawal rate from the confined Yorktown-Eastover system is near the highest reported rate for the 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 small rates of inflow from coastal water bodies underscore ongoing concerns about up-coning and lateral intrusion of salty groundwater.

Released January 15, 2026 08:27 EST

2026, Science (391) 294-299

David R. Shelly, Amanda M. Thomas, Kathryn Z. Materna, Robert J. Skoumal

Accurate tectonic models are essential for assessing seismic hazard and fault interactions. However, the plate configuration at the complex Mendocino triple junction, where the San Andreas Fault and the Cascadia subduction zone meet, remains uncertain. We analyzed fault slip associated with a recently identified zone of tectonic tremor and low-frequency earthquakes (LFEs) near the southern edge of the subducting Gorda slab. Based on tidal sensitivity and P-wave first motions, we show that the LFEs are generated by dipping, strike-slip motion. This suggests that a former Farallon slab fragment, now captured by the Pacific plate, is translating northward beneath westernmost North America. This geometry effectively extends the slab interface fault, challenging prevailing interpretations of slab window formation and creating a potential unaccounted earthquake hazard in this region.

Released January 14, 2026 11:50 EST

2026, Fact Sheet 2026-3059

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

Using a geology-based assessment methodology, the U.S. Geological Survey estimated undiscovered, technically recoverable mean continuous and conventional resources of 1.6 billion barrels of oil and 28.3 trillion cubic feet of gas in the Woodford and Barnett Shales of the Permian Basin Province.

Released January 14, 2026 11:01 EST

2026, Water Resources Research (62)

Molly Van Appledorn, Nathan R. De Jager, Jason J. Rohweder, Marcella Windmuller-Campione, Daniel Griffin

The hydrologic regime of the upper Mississippi River (UMR) has become wetter, with greater discharges, longer-lasting high-flow conditions, and seasonal shifts in these patterns over the past several decades. How these changes are expressed spatially as floodplain inundation area, frequency, depth, duration, and timing is not well understood. It is also unclear to what degree spatial patterns of submergence are represented by examining discharge data alone. We assessed changes in floodplain inundation characteristics from 1940 to 2022 in navigation pools 3–10 of the UMR using a geospatial model to simulate daily inundation depths. Inundation characteristics shifted significantly across pools, but the direction and magnitude of change varied by pool and metric. Characteristics summarized at the pool scale correlated with streamgage-derived proxies but the strength of the relationship varied. Within pools, variability in inundation trends highlighted the importance of spatially explicit modeling. Our study demonstrates that changes in discharge over 83 years have manifested across the UMR floodplain in ways that may have consequences for ecological patterns and processes. By mapping hydrologically sensitive areas, we can anticipate which areas may be susceptible to additional shifts in river discharge in a climatically uncertain future.

Released January 14, 2026 09:24 EST

2026, Annals of the American Association of Geographers

Roger Sayre, Madeline T. Martin, Nadia Naji, Kelly B. Sides, Jill Janene Cress, Kevin Butler, Keith Van Graafeiland, Deniz Karagulle, Charlie Frye, Sean Breyer, Dawn Wright, John Klavitter, Dena Spatz, David Will, Gregg Howald, Alex Wegmann, Charlotte K. Stanley, Nick Holmes

To support conservation-focused research and management we developed a new 30-m resolution polygon data layer of the nonlacustrine and nonriverine islands of the United States, with associated attributes describing key physical and conservation geography characteristics. Islands were grouped into a three-tiered hierarchy of island regions (twelve), island provinces (twenty-eight), and individual islands (19,023). Islands were classified as either continental or oceanic based on their physiographic position relative to the North America continental shelf, and estuarine versus nonestuarine depending on their location within or external to estuaries. For each island we assessed the diversity of terrestrial and coastal ecosystems, the number of threatened and endangered (T&E) species listed under the Endangered Species Act, the number of T&E species critical habitats, the number of migratory bird species listed under the Migratory Bird Treaty Act, the number of Key Biodiversity Areas, and the number of and management responsibility for protected areas. We conclude that the conservation importance of islands is disproportionate to their total area as, for example, islands contain 52 percent of the T&E species yet their total area is only 2 percent of the area of the continental mainland. Similarly, of the global total of 431 World Terrestrial Ecosystems, 201 (47 percent) occur on U.S. islands compared with 286 (66 percent) that occur on the U.S. continental mainland.

Released January 14, 2026 08:50 EST

2026, JGR Biogeosciences (131)

Jessica R. Lacy, Samantha C. McGill, Janet Thompson, Rachel Allen, Francis Parchaso, David Hart, Lukas T. WinklerPrins, Joseph K. Fackrell, Andrew W. Stevens

The erodibility of cohesive sediment is known to vary both spatially and temporally but the factors governing its variation are not well understood. We conducted a field investigation of the influence of hydrodynamic forcing, sediment properties, and benthic infauna on erodibility in the muddy shallows of San Pablo and Grizzly Bays in northern San Francisco Bay in summer 2019 and winter 2020. An erosion rate parameter M_c was determined from regressions between near-bed vertical turbulent sediment flux, as a proxy for erosion, and bed shear stress due to currents. During each 2-month study period, we measured benthic infauna abundance and dry bulk density, particle size distribution, percent organic carbon, chlorophyll a, pheophytin a, and carbohydrates carbon concentrations of surficial bed sediments five or six times. M_c increased with bed shear stress due to waves in both embayments. In San Pablo Bay, erodibility was approximately 50% lower during the winter than the summer deployment, whereas in Grizzly Bay, there was no significant difference. The factor most strongly related to the decrease in M_c in San Pablo Bay was increased abundance of the amphipod Ampelisca abdita. The observed reduction in erodibility may occur in many muddy estuaries because A. abdita is broadly distributed in the coastal waters of North America. Erodibility was also directly related to biomass of the invasive clam Potamocorbula amurensis. Erodibility did not depend on dry bulk density: bulk density did not vary seasonally in San Pablo Bay and was lower in winter than summer in Grizzly Bay.

Released January 14, 2026 08:39 EST

2026, Earth Surface Processes and Landforms (51)

Sanne M. Vaassen, Karin R. Bryan, Andrew Swales, Joel A. Carr, Conrad A. Pilditch

Many estuaries worldwide face increasing sediment loading caused by catchment land use change and intensification, creating subsequent adverse effects on estuarine ecosystems. Extreme weather events can disproportionately alter sediment pathways and loading. Although storm-driven sediment exchange has been widely examined at open coasts and inlets, key transport mechanisms within constricted, sheltered estuaries remain understudied.

This study presents an observational dataset capturing the impact of a 99th percentile water-level event (based on 20 years of records) on sediment transport pathways in a sheltered, barrier-enclosed estuary. This event, driven by a 3-day storm surge (>0.5 m) combined with a spring tide, was recorded during a 3-week field campaign.

Sediment transport pathways and riverine contributions were analysed, and observations revealed substantial changes in suspended sediment concentrations increasing from 18 mg/l to 70 mg/l during the event. The elevated water levels and resulting pressure gradient at the constricted study site entrance caused by the storm surge increased local flood dominance. Combined with higher flow velocities and resuspension, the storm led to a sixfold increase in sediment import at the estuary entrance and a 600-fold increase in sediment flux to the upper estuary.

The decoupling of peak suspended sediment concentrations from streamflow indicates that the resuspension of estuarine legacy sediment, rather than catchment inputs, dominated the system's response.

These findings challenge assumptions about estuarine sediment budgets and emphasise that incorporating high water-level surge events into models can enhance the prediction of long-term estuarine evolution. Given projected increases in storm frequency under climate change, understanding these episodic but highly consequential sediment pulses can support the assessment of wetland resilience and inform estuarine management strategies.

Released January 14, 2026 07:45 EST

2026, Frontiers in Conservation Science (6)

Samuel J. Brenkman, Jeffrey J. Duda, Rebecca M. McCaffery, Katie E. Kierczynski, Marshal S. Hoy, Trevor J. Kumec,  William Baccus, Caren Suzanne Goldberg, Carl O. Ostberg, Steven C. Fradkin

Across the western United States, introductions of non-native fish into historically fishless mountain lakes have impacted native biota. Understanding the impacts of fish introductions is essential for conservation in Olympic National Park, a Biosphere Reserve. We reconstructed fish plantings using records dating back to 1930, followed by environmental DNA (eDNA) surveys to estimate the current distribution of fish and amphibians in 117 remote mountain lakes. We used Bayesian multiscale occupancy models to determine how lake attributes and planting history related to fish and amphibian occupancy. The most frequently detected species were Brook Trout, Rainbow Trout, Cascades Frog, and Northwestern Salamander. eDNA sampling revealed 52 lakes with amphibians only, 45 with fish and amphibians, 14 with fish only, and 6 unoccupied. Of the 53 lakes with planting records, 38 had fish eDNA detected. Fish eDNA was also detected in 21 lakes lacking planting records, which could reflect incomplete records, unauthorized plantings, and false positive detections. Of the three species planted, Cutthroat Trout had the highest failure rate and did not become established in 23 of 28 historically planted lakes. In a subset of 9 lakes sampled for up to 7 years, those with known fish and amphibian presence showed consistent eDNA detections over time. The number of times a lake was stocked was the best predictor of occupancy for Brook and Rainbow trout, while higher occupancy for Brook Trout was also associated with lower elevations, lower solar radiation, and larger lake area. We did not observe widespread negative associations between amphibian occupancy and fish presence, although there was a negative relationship between fish presence and Rough-skinned Newt and Long-toed Salamander occupancy. Cascades Frog occupancy showed no relationship to fish presence or lake traits. Our results suggest mechanisms of fish persistence over time and highlight areas where native amphibians are impacted by introduced fish. These results can guide management options like targeted fish removals that benefit native fauna while still supporting recreational fishing. More broadly, our work demonstrates the value of combining historical records with contemporary surveys and the utility of eDNA for broad-scale surveys of species distribution in remote wilderness areas.

Released January 14, 2026 07:40 EST

2026, JGR Solid Earth (131)

Helene Le Mevel, Nathan Lee Andersen, Annika E. Dechert, Josef Dufek

From 2019 to 2024, gravity surveys were conducted at the Three Sisters volcanic cluster (TSVC), measuring 246 gravity sites using a spring relative gravimeter. We calculated the residual Bouguer anomaly and identified three main zones with negative anomalies, ranging from −4 to −8 mGal, located southwest and west of South Sister, within an area that has been uplifting for the past two decades. After inversion, we obtain a 3D density model of the subsurface and identify low-density bodies extending from the surface down to 3 km. We estimate a total of 15 km³ of crustal bodies with density close to 2 g/cm³ that could store up to ~5 km³ of water, forming an extensive hydrothermal system beneath the TSVC. We explore the possible combinations of melt compositions and temperatures that could create a bulk density close to our reference crustal density (2.5 g/cm³) using MELTS thermodynamic simulations. Our results indicate that a magmatic mush with as little as 15% partial melt of bulk rhyolitic composition or as much as 52%–57% partial melt of a bulk dacitic composition could be stored in a magmatic system under TSVC without generating a detectable gravity anomaly. Episodic magma injections at the base of the magmatic system, such as the 1998–2000 intrusion at ~6 km depth, would bring heat and gas to the hydrothermal system while maintaining a low melt fraction in the magmatic mush, as imaged at other Cascade volcanoes.

Released January 13, 2026 10:58 EST

2026, Animal Biotelemetry (14)

Richard Kraus, James Roberts, Mark Richard Dufour, Branden E. Kohler

Increasing application of acoustic telemetry for determining survival, migration and habitat use of fishes highlights the need to improve interpretation of tracks that end abruptly: when is fishing mortality, predation, or some other cause to be inferred? Significant technological advances have led to the development of tags that “sense” predation and can be used to infer information about the type of predator that consumed the tagged fish. However, growing evidence suggests that bird predation is not effectively quantified by the technology. We hypothesized that reduction in sound transmission from acoustic tags in the gut of a bird combined with short bird diving intervals would eliminate detections of acoustic telemetry tags from the surface and severely reduce detection efficiency at depth. We test this hypothesis indirectly with two experiments using cormorant carcasses containing tagged fish in which carcasses were either tethered to a mooring for several hours or lowered through the water to simulate diving behavior. Detection of tagged prey fish in the gut of bird carcasses was severely reduced or negated completely, supporting our hypothesis. By comparison, as expected, tagged fish that were not in the gut of bird carcasses were detected at a higher frequency. Depth and distance to passive moored receivers also affected detection probability of tagged fish with more detections at depth and when closer to the receiver. Our results emphasized the importance of accounting for avian predation of tagged fish in studies of prey species in surface waters. Further, while recent development of predation sensing tags has illustrated a few examples of bird predation, our results demonstrate that determining that a tagged fish has been consumed by a diving bird will be difficult and will likely require alternative methods or technologies.

Released January 13, 2026 09:43 EST

2026, Journal of Fish Biology

Michelle A. Briggs, Molly A. Webb, Christopher S. Guy, Todd M. Koel

Yellowstone cutthroat trout Oncorhynchus virginalis bouvieri in Yellowstone Lake are the focus of intensive conservation efforts due to the threat of predation by invasive lake trout Salvelinus namaycush. Suppression gillnetting has reduced the abundance of predatory lake trout, and the Yellowstone cutthroat trout population is recovering. Long-term monitoring indicates the size structure of the population shifted following lake trout invasion, suggesting that reproductive demographic rates of Yellowstone cutthroat trout may have changed. Length at 50% probability of maturity, as assessed using histological analysis of gonadal tissue, was 479 mm (95% confidence interval [CI] 467–490 mm) for females and 406 mm (95% CI 386–430 mm) for males, compared to 330 mm for males and females historically. Currently, age at 50% probability of maturity is 6.6 for females and 5.4 for males. The rate of skipped spawning was 3% for females and 38% for males. Mean absolute fecundity was 2897 ovarian follicles/individual at present compared to 1141 ovarian follicles/individual before lake trout invasion. Mean relative fecundity was 2157 ovarian follicles/kg. This research illustrates the plasticity in the reproductive strategies of fishes as a result of an invasive species. Understanding the reproductive biology of fish populations is vital for effective fisheries management, and these results are integral to a population model that can be used to develop new conservation benchmarks for Yellowstone cutthroat trout.

Released January 13, 2026 09:08 EST

2026, Frontiers in Ecology and the Environment

Brendan K. Hobart, Wynne Emily Moss, Max C Cook, R. Chelsea Nagy, Valerie J McKenzie

Novel stressors like climate change and biological invasions alter ecological communities, resulting in changes to ecosystem services and biodiversity (that is, ecological transformation). Most ecological transformation research focuses on plants, but animals are likely affected by and plausibly mediate the extent, impact, and pace of transformations. In western North America, where invasive annual grasses are transforming sagebrush-steppe ecosystems, we quantified how transformation in vegetation drives songbird community change. We hypothesized that transformation of vegetation communities via invasion alters wildlife communities by favoring generalists over specialists, albeit with extinction debts that may temporarily obscure transformations’ consequences. Although local-scale songbird diversity increased with grass invasion, we found that this shift was accompanied by ongoing reorganization of songbird communities at larger scales (across the sagebrush biome), driven by heterogeneous impacts of invasion among species and guilds. Through our biome-wide analysis, we were also able to identify high-priority regions for conservation of sensitive songbird species. Our research provides evidence that wildlife communities are transforming alongside vegetation communities and offers insight into the approaches required to quantify nascent community turnover.

Released January 13, 2026 08:44 EST

2026, Ecosystems (29)

Charles R. Lane, Ellen D’Amico, Jay R. Christensen, Heather E. Golden, Frederick Y. Cheng, John C. Hammond, Admin Husic, Kristin L. Jaeger, C. Nathan Jones, Christa A. Kelleher, Li Li, D. Tyler Mahoney, Hilary K. McMillan, Adam N. Price, Roy Sando, Catalina Segura, Erin C. Seybold, Adam S. Ward, Margaret Zimmer

Headwater streams play critical roles in hydrologic and biogeochemical processes and functions, yet their spatial distribution and land cover context remain poorly understood at continental scales, and no dedicated geospatial dataset exists. Building from a high-resolution conterminous United States (CONUS) hydrography network dataset, we quantified the spatial extent, density, and upstream catchment characteristics of headwater stream segments across the CONUS. We identified approximately 8.4 million kilometers of headwater streams, finding that 77% of the total stream network consists of headwaters, nearly double the total length represented in prior estimates. Stream density varied fivefold across regions, from < 1 km·km⁻² in arid basins to > 5 km·km⁻² in humid, forested areas. Over 73% of the CONUS landmass drains from headwater streams. The majority of headwater stream length occurred in forested and cultivated catchments across the CONUS, while substantial regional differences were evident for headwater stream distribution in other land cover classes (for example, wetlands, urban areas, shrublands, and herbaceous-dominated catchments). The dedicated and novel geospatial dataset, HELiOS (HEadwater streams and Low-Order Systems) is introduced for management and research use. The HELiOS dataset provides the first continental-scale, high-resolution characterization of headwater streams, offering new insights and opportunities for hydrologic modeling, ecological assessments, and environmental policy.

Released January 12, 2026 09:30 EST

2026, PeerJ (14)

Robert S. Cornman, Mark J. Hepner, Clint Otto

Background

Detailed studies of foraging behavior are needed for scientific management of the endangered rusty-patched bumblebee (Bombus affinis) in the disjunct and ecologically differentiated habitats it presently occupies. Current knowledge gaps hinder recovery planning but are challenging to redress through direct observation of rare interactions in the field.

Methods

We used genetic metabarcoding to characterize the taxonomic composition of pollen collected by B. affinis workers in the Appalachian mountains of Virginia and West Virginia from 2021–2023. We developed a custom sequence database of the regional flora and compared results for two independent genetic loci, internal transcribed spacer 1 and internal transcribed spacer 2 (ITS1 and ITS2).

Results

While ITS2 consistently detected more plant diversity, results from the two loci were broadly concordant with a few notable exceptions. The plant genera Hydrangea, Actaea, Rhododendron, Tilia, and (unexpectedly) Laportea were prominent in midsummer samples, with Rubus a consistent contributor in late spring and early summer. Pea flowers (family Fabaceae) were relatively infrequent but the genera Securigera and Trifolium were detected before the Hydrangea bloom and again in late summer afterwards. The diversity of forage plants was highest in late summer, driven primarily by various genera of Asteraceae. Comparing the current data with previous work indicates regional differentiation in forage plants between Appalachia and the upper Midwest, but also allows ‘consensus’ forage sources that are supported by multiple lines of evidence and shared between regions to be tabulated. These results should help managers focus survey efforts for this endangered species and plan habitat enhancements.

Released January 12, 2026 08:48 EST

2026, Preprint

Kathryn A. Thomas

The long-term viability of the iconic Joshua tree of the Mojave Desert is being evaluated. In 2022, we measured the abundance and heights of Joshua tree stems on 62 1000 m2 plots in the eastern Mojave Desert of California. The 2022 plots were represented by 33 plots in a population of western Joshua trees (Yucca brevifolia) and 29 plots in a population of eastern Joshua trees (Y. jaegeriana). Five plots had no Joshua trees; two of which in the western Joshua tree population were destroyed by fire. The 57 plots with Joshua trees supported 627 stems ranging from stems less than 25 cm to mature trees. The plots were examined by abundance and four size classes of the stems, as indicated by their height and indicative of their reproductive status. The western Joshua tree population had more stems (407) than the eastern population (220 stems) although the median difference was not significant. The western population had significantly more stems in pre-reproductive size classes, juvenile (> 25 cm to one-meter, p = 0.001) and sub_adult (> one- to two-meters, p = 0.001), than the eastern population but significantly fewer adult (> two-meters, p < 0.001) stems. The eastern population had significantly greater mean height of adult trees (302 cm +/- 75.3 cm) than the western population (263 cm +/- 52.9 cm) and a larger proportion of stems taller than two-meters (56%) than the western population (39%). In contrast, the western population has 77% of its population in younger size classes (> 25 cm to < 2 meters). These abundance and size class measures alone do not predict whether either population has sufficient natural stand regeneration for long term persistence, but the younger size class structure of the western population suggests greater long-term resilience than for the eastern population.

Released January 12, 2026 07:52 EST

2026, IEEE Transactions on Antennas and Propagation

Sepehr Eskandari, Asem Melebari, Paul J. Kinzel, Robert Russell Lotspeich, Jack R. Eggleston, Mahta Moghaddam

We provide an overview of an uncrewed aircraft system (UAS)-based software-defined radar (SDRadar) system for high-resolution geophysical observations. The radar transceiver is implemented on a Radio Frequency System-on-Chip (RFSoC) platform, along with an ultra-wideband Vivaldi antenna that has a starting operating frequency of 150 MHz, enabling the system to be used across different applications, including measurements of freshwater bathymetry. In addition to system design and subsystem performance assessments, this paper presents the results of field testing conducted along the Sacramento River near Glenn, California, USA. The radar-derived river depth measurements are compared with ground truth data collected from a crewed boat using an acoustic Doppler current profiler (ADCP). The results show good agreement, with a root mean square error (RMSE) of 0.079 m.

Released January 12, 2026 07:44 EST

2026, Hydrological Sciences Journal

Melanie K. Vanderhoof, William Keenan, Wayana Dolan, Heather E. Golden, Charles R. Lane, Jay R. Christensen, Kylen Solvik, Adnan Rajib

Globally, many wetlands and lakes are at risk for further loss, which can amplify downstream consequences of flood and drought events. We derived remotely sensed based time series of surface water storage (SW_storage) to determine when and where accounting for SW_storage dynamics improves predictions of river discharge. We trained four long short-term memory (LSTM) models, that differed in their inclusion of storage data and catchment characteristics, to simulate daily river discharge (2016–2023) for select watersheds across the conterminous United States. Adding SW_storage to a meteorology-only or meteorology-and-catchment characteristics model improved upon model Nash-Sutcliffe efficiency (NSE) in 80.6% of the watersheds. Residuals during low-flow (Q70) events decreased by 47.6% when adding storage to meteorological data. Improvements were most consistent in ecoregions with a greater abundance of non-floodplain lakes and wetlands. This effort represents the first exploration to train a multi-watershed LSTM on landscape-scale remotely sensed time series of SW_storage.

Released January 09, 2026 11:50 EST

2026, Fact Sheet 2025-3053

Rand Gardner, Jason A. Flaum, Seth S. Haines, Justin E. Birdwell, Scott A. Kinney, Brian A. Varela, Katherine L. French, Janet K. Pitman, Stanley T. Paxton, Tracey J. Mercier, Christopher J. Schenk, Heidi M. Leathers-Miller, Hannah D. Shook

Building on a geology-based assessment of undiscovered, technically recoverable hydrocarbon resources within the Haynesville Formation, the U.S. Geological Survey estimated the water and proppant necessary for development of the remaining resources associated with the Haynesville Sabine Uplift Continuous Gas Assessment Unit. Additionally, projections have been made on the volume of wastewater expected as a byproduct of possible future development. This fact sheet presents an overview of the methodology, along with the inputs and results of the Haynesville water and proppant assessment.   

Released January 09, 2026 11:50 EST

2026, Scientific Investigations Report 2025-5105

Jill N. Densmore, John A. Izbicki, Meghan C. Dick, Sandra Bond

The U.S. Army Fort Irwin National Training Center (NTC), 120 miles northeast of Los Angeles in the Mojave Desert of California, obtains its potable water supply from the Bicycle Valley and Langford Valley groundwater basins; Langford Valley groundwater basin is further subdivided into the Langford Well Lake and Irwin groundwater subbasins. The Irwin groundwater subbasin contains younger, unconsolidated deposits that have a saturated thickness of as much as 200 feet (ft) and a lower aquifer within older unconsolidated deposits as thick as 650 ft. Groundwater recharge under predevelopment conditions (before 1941) occurred primarily from infiltration of intermittent streamflow in small washes that cross the Irwin groundwater subbasin. Since that time, groundwater recharge has increased because of growth of the NTC in recent years and as a result of other processes, including (1) infiltration of treated wastewater into the aquifer through ponds near the NTC wastewater treatment facility (WWTF) and (2) infiltration of imported water and treated wastewater used for landscape irrigation at base housing and athletic fields.

Water samples were collected from 17 wells and analyzed for field parameters, chemical constituents, and isotope composition in the Irwin groundwater subbasin between 2014 and 2019. These data were supplemented with water- chemistry data collected during 1993–95 and at other times if available. Between 1993–95 and 2015–19, median dissolved solids and nitrate concentrations in water from wells in the Irwin groundwater subbasin increased from 620 to 1,030 milligrams per liter (mg/L) and from 2.8 to 4.5 mg/L as nitrogen, respectively. After 2014, dissolved solids and nitrate concentrations in water from wells near the NTC WWTF decreased as a result of changes in source water quality attributable to reverse osmosis of treated drinking water delivered within the Irwin groundwater subbasin and to increased levels of treatment at the NTC WWTF. Based on delta oxygen-18 and delta deuterium isotope data, increases in dissolved solids concentrations in water from most wells were consistent with evaporation prior to recharge and mobilization of soluble salts from the unsaturated zone. Arsenic and fluoride concentrations in water from wells decreased between 1993–95 and 2015–19 as the basin filled with treated wastewater, but 2015–19 concentrations generally exceeded the California State Water Resources Control Board maximum contaminant levels of 10 micrograms per liter for arsenic and 2 mg/L for fluoride. Most groundwater in the Irwin groundwater subbasin has unadjusted carbon- 14 ages ranging from 18,400 to 12,350 years before present. However, water from well 10E3, located along the wash near the subbasin outflow in the southeastern part of the Irwin groundwater subbasin, contained measurable tritium, which is consistent with infiltration of intermittent streamflow and groundwater recharge from the wash after 1952. Chemical and isotopic data indicate that treated wastewater is present in water from most wells in the upper aquifer that underlies the Irwin groundwater subbasin. Wells were not sampled to adequately determine the extent of treated wastewater and changes in water quality within the lower aquifer that underlies the Irwin groundwater subbasin.

Released January 09, 2026 08:48 EST

2025, Arizona Geological Survey Contributed Report CR-25-A

Gordon B. Haxel, Carl E. Jacobson, Gabe S. Epstein

No abstract available.

Released January 09, 2026 08:45 EST

2026, Fisheries

Paul L. Angermeier, James R. Karr, Chris O. Yoder, Robert M. Hughes

Karr (1981), which introduced the index of biotic (or biological) integrity (IBI) has been cited more often (>4,500 times) than any other paper in Fisheries. In this essay, we reflect on the historical context of this seminal publication and its broad, continuing impact on the management of natural resources, especially freshwater ecosystems.

Released January 08, 2026 10:47 EST

2026, Environmental Pollution (392)

Adam Biales, M. S. Hu, D. C. Bencic, M. J. See, Susan T. Glassmeyer, E.T. Furlong, Julia M. Stelman, W. Huang, Dana W. Kolpin, Marc A. Mills, L. D. Brunelle, Angela L. Batt, S. Thomas Purucker

Wastewater treatment effluents (WWTE) present complex risks to aquatic ecosystems that are difficult to characterize using traditional methods. This study systematically evaluated the consistency and performance of transcriptomic-based approaches over time with repeated sampling and with differing experimental approaches (selection of reference condition, grab vs. composite sampling, deployed vs. laboratory exposed). RNA-seq was performed on larval fathead minnow (FHM) exposed in the laboratory to moderately hard reconstituted water (MRHW) or individual grab samples collected from an upstream site and a WWTE in the morning and afternoon over two successive days, as well as FHM deployed concurrently with grab sampling at the same sites. Composite transcriptional profiles were generated by pooling count data from grab sample exposures. The choice of comparator significantly affected results. The use of the upstream site as the reference consistently yielded fewer differentially expressed genes (DEGs) and minimal overlap compared to DEG sets from the other comparisons. Using MRHW as a comparator, DEG sets showed high consistency across grab samples, with morning samples demonstrating larger, highly consistent gene expression sets (96 % overlap) compared to afternoon samples, revealing clear and consistent within-day expression patterns. With the MHRW comparator, DEG sets from grab sample composites and deployments also overlapped substantially, indicating that transcriptional profiles accurately reflect WWTE composition regardless of exposure method. Comparisons with non-targeted (NTA) and targeted analytical datasets confirmed that gene expression interpretations aligned with effluent composition while highlighting limitations of relying solely on targeted analyte sets for connecting expression to specific chemicals. Though highly dependent on experimental design, these results demonstrate that transcriptomic-based methods provide significant utility for characterizing the bioactivity of complex environmental mixtures, and when paired with NTA datasets, have the potential to deliver a comprehensive assessment of their environmental risk.

Released January 08, 2026 07:38 EST

2026, Oecologia (208)

Kristoffer T. Everatt, Leah Andresen, Jennifer F. Moore, James E. Hines, Graham I.H. Kerley

Apex carnivores are known to regulate ecosystem structure and function, including via interactions with syntopic, competitively inferior carnivores. These effects may be dependent on relative carnivore density and resource availability or productivity. We investigated the functional effect of African lions as an apex carnivore on the presence of co-occurring large carnivore species across two adjoining National Parks that contrast in relative densities of carnivores and prey. We employed two-species occupancy models from track data to test statistical interactions between lions and the other syntopic large carnivore species, while accounting for each species’ habitat selection. We further investigated the influence of anthropogenic and natural variables on these co-occurrence dynamics. Our models revealed that the occurrence of each carnivore species was best predicted by access to their own key resources. We also found significant statistical interactions between lions and cheetahs, lions and leopards, and lions and spotted hyenas in resource-rich landscapes. Finally, we found limited support for the competition exclusion hypothesis between most species, with the exception of lion-African wild dog co-occurrence patterns. Species’ co-occurrence dynamics were all influenced by resource availability, with lion-leopard and lion-cheetah co-occurrence decreasing strongly with increasing resource availability. Most species co-occurrence declined with increasing occurrence of lions. The patterns revealed by this study improves predictions of how changes in resource availability and carnivore occurrence could impact carnivore community dynamics and the functional role of apex carnivores.

Released January 07, 2026 09:15 EST

2026, Global Change Biology (32)

Lucia I.A. Enes Gramoso, Dustin Carrol, Kyle C. Cavanaugh, Remi Bardou, Michael Osland, Tom Van der Stocken

Warming winter temperatures are driving range expansion of tropical, cold-sensitive mangroves into temperate ecosystems. Along the Atlantic coast of North America, the mangrove range limit is particularly sensitive to climate variability and historical data demonstrate that the mangrove-salt marsh ecotone on this coast has shifted recurrently during recent centuries. However, a comprehensive understanding of how this mangrove-salt marsh ecotone may shift in the future remains lacking. Here, we combine ensemble forecasting of mangrove distribution for the next century with high-resolution oceanographic dispersal simulations, phenological observations, and historical hurricane data to project future mangrove-salt marsh dynamics at the rapidly changing range limit in northeastern Florida (USA). We show that warming winter temperatures will drive continued poleward expansion of mangroves along North America's Atlantic coast, potentially reaching South Carolina by 2100. With ongoing climate change, suitable mangrove habitat is projected to expand beyond the current range limit, and dispersal simulations suggest successful colonization of these sites from established mangrove populations. Additionally, patterns in hurricane directionality and intensity and field reports of propagule presence reveal that these high-energy events may significantly contribute to future mangrove range expansion by facilitating long-distance, storm-driven propagule dispersal. The encroachment of mangroves in salt marsh-dominated latitudes is expected to substantially modify wetland ecosystem function and structure, emphasizing how the identification of newly colonizable habitat can inform conservation strategies and site-specific decisions on mangrove management.

Released January 07, 2026 07:52 EST

2026, Global Ecology and Conservation (65)

Michael J. Gould, Ted Swem, Guthrie S. Zimmerman, Brian A. Millsap, Jay V. Gedir, Fitsum Abadi

Species status assessments require an understanding of underlying population dynamics and important drivers of species demography. Large-scale assessments can be difficult due to challenges collating data obtained through different methods and different sources at multiple scales. Integrated population models (IPMs) provide a unified framework to combine multiple data sources and jointly estimate population parameters over a large spatiotemporal scale. We developed separate IPMs to estimate abundance and demographic rates for a northern (NMP) and southern (SMP) management population of peregrine falcons (Falco peregrinus) in North America from 2008 through 2019 (SMP) and 2020 (NMP). An outbreak of highly pathogenic avian influenza (HPAI) starting in 2021 led us to extend our modeling effort to assess its impact on these management populations by updating both IPMs using index data of population size through 2024 in a predictive framework. Survival probabilities differed drastically between first-year and after-first-year individuals in both management populations. After-first-year survival was nearly identical between the NMP and SMP, but first-year survival was lower in the SMP. Mean productivity was significantly lower in the NMP compared to the SMP, whereas the probability of breeding was similar in both management populations. Estimated total abundance for the NMP was substantially larger than the SMP, representing most of the North American peregrine population. Population growth was positive for both management populations, albeit at a slower rate for the NMP. The NMP declined from 2017 to 2018 coinciding with a drop in 2018 estimated productivity. When we extended the IPMs with updated count data through 2024, the NMP slightly declined but estimated abundance remained above levels at the start of the time series analyzed. The SMP grew at a similar rate to that predicted during the period informed by demographic data. We did not detect a continental-scale change in population size or trajectory in either management population associated with the arrival of HPAI in 2021. Further monitoring can support determination of whether the declines in the NMP were temporary, can enhance understanding of the underlying mechanisms, and can be used to guide the conservation and management of the peregrine falcon population in North America.

Released January 07, 2026 07:29 EST

2026, Fact Sheet 2025-3055

Faith Fitzpatrick, Collin Roland, Angus Vaughan, Zhenduo Zhu, David Soong, Rachel Sortor

The FluOil tool was developed to help with planning and early response for oil spills in rivers where subsurface oil-sediment interactions result in the formation of oil-particle aggregates (OPA). The turbulence and variable velocity associated with water flowing within a natural stream channel creates the conditions needed for an oil slick to break up into small droplets and mix in the water column, collide with sediment or organic detritus, and form OPA. This process is similar to what occurs due to wave action along a shoreline. The FluOil tool estimates how fast OPA travel downstream in rivers as well as when and where they may deposit. The FluOil tool relies on pre-existing channel hydraulic data along with user-specified OPA characteristics of size, settling velocity and critical shear stress to compute OPA transport. It is important to know where OPA are transported and accumulated for understanding potential impacts on drinking water intakes, burial of sensitive habitat beds, potential toxicity to benthic organisms, and prolonged sheening from resuspension and breakup. OPA tend to accumulate with fine-grained (silt and clay) sediment deposits (“mud” or “muck”) in backwater areas, oxbows, side channels, pools, and other slow-moving reaches of rivers during low flows. Deposited OPA can be resuspended during high flows, driving continued environmental impact concerns that may extend beyond typical oil spill response timelines.

Released January 06, 2026 11:26 EST

2026, Journal of Hydrology: Regional Studies (63)

Saalem Tilahun Adera, Ayman H. Alzraiee, Richard G. Niswonger, Enrique Triana, Derek W. Ryter, John A. Engott

Released January 06, 2026 09:53 EST

2026, Seismological Research Letters

Andrea L. Llenos, David R. Shelly, Allison Shumway

Earthquake catalogs are essential data inputs for seismic hazard modeling. Because earthquake magnitudes are reported in a variety of types (e.g., local magnitudes and moment magnitudes), magnitude conversion relationships must be used to convert the different magnitude types present in a catalog to a uniform magnitude type to avoid biases in the hazard computation. However, these conversion relationships are often uncertain and have been shown to sometimes perform poorly. Here, we investigate the sensitivity of the gridded seismicity component of the National Seismic Hazard Model (NSHM) to the catalog conversion equations in the Eastern United States. In the 2023 NSHM, magnitudes of various types were converted to moment magnitudes using equations developed by the Central and Eastern United States Seismic Source Characterization for Nuclear Facilities (CEUS‐SSCn), based on least‐squares (LS) regressions made using data from a catalog containing events up through 2008. We recompute these equations using events in the Advanced National Seismic System Comprehensive Earthquake Catalog with multiple magnitudes from 2000 to 2023. Although we prefer the use of orthogonal regressions for our datasets, LS regressions produce broadly similar results, with both approaches exhibiting large deviations from the CEUS‐SSCn conversions, especially at smaller magnitudes. We compare the spatial distribution of annual rates using three different models: (1) the 2023 NSHM conversions, (2) our updated conversions, and (3) no conversions. We find that the choice of conversions leads to substantial differences in the rate forecasts, which can greatly impact the seismic hazard model, particularly in regions with low‐seismicity rates such as the Eastern United States, where the hazard is dominated by gridded seismicity rather than a fault model.

Released January 06, 2026 08:09 EST

2026, Frontiers in Water (7)

John C. Hammond, Phillip J. Goodling, Jeremy Alejandro Diaz, Hayley R. Corson-Dosch, Aaron Joseph Heldmyer, Scott Douglas Hamshaw, Ryan R. McShane, Jesse Cleveland Ross, Roy Sando, Caelan Simeone, Erik A. Smith, Leah Ellen Staub, David Watkins, Michael Wieczorek, Kendall C. Wnuk, Jacob Aaron Zwart

Forecasts of streamflow drought, when streamflow declines below typical levels, are notably less available than for floods or meteorological drought, despite widespread impacts. We apply machine learning (ML) models to forecast streamflow drought 1–13 weeks ahead at 3,219 streamgages across the conterminous United States. We applied two ML methods (Long short-term memory neural networks; Light Gradient-Boosting Machine) and two benchmark models (persistence; Autoregressive Integrated Moving Average) to predict weekly streamflow percentiles with independent models for each forecast horizon. ML models outperformed benchmarks in predicting continuous streamflow percentiles below 30%. ML models generally performed worse than persistence models for discrete classification (moderate, severe, extreme) but exceeded the benchmark models for drought onset/termination. Performance was better for less intense droughts and shorter horizons, with predictive power for 1–4 weeks for severe droughts (10% threshold). This work highlights challenges and opportunities to advance hydrological drought forecasting and supports a new experimental forecasting tool.

Released January 06, 2026 08:06 EST

2026, Geosphere

Erin Todd, Jonathan Saul Caine, Michael Bizimis, Andrew R.C. Kylander-Clark, Robert Reece Hammond, Alicja Wypych

New geochemical, U-Pb geochronology, and Sr-Nd-Hf isotope data provide evidence for the tectonic evolution of the Seventymile terrane in interior Alaska, USA. Ultramafic and mafic rocks of the Seventymile terrane are thought to represent components of a dismembered ophiolite and provide unique constraints on regional terrane evolution and accretion. The Seventymile ophiolite represents fragments of the Devonian to Permian Slide Mountain Ocean (SMO) that separated allochthonous and parautochthonous continental fragments of western North America. It now occurs as multiple thrust sheets containing Permian mafic and ultramafic rocks overlying and/or possibly imbricated with amphibolite-facies supracrustal rocks of allochthonous Yukon-Tanana terrane and parautochthonous North America. Seventymile klippen contain variably serpentinized peridotite, primarily harzburgite, low-grade meta-mafic rocks, and minor oceanic sedimentary rocks (argillite, chert, limestone, and metasandstone). Mafic rocks include gabbro to diabase, typically as dikes, veinlets, or rare massive stocks intruding peridotite. Mafic rocks also include greenstones of the Seventymile assemblage in klippen structurally underlying, and in shear zone contact with, ultramafic klippen.

New trace element and radiogenic isotope data from mafic magmatic rocks associated with the Seventymile ultramafic bodies show evidence for a weakly subduction-modified mantle source, like the mantle source of normal mid-ocean-ridge basalt (N-MORB) or back-arc basin basalt (BABB). Seventymile assemblage greenstones are more heterogeneous. They range from N-MORB to enriched mid-ocean-ridge basalt (E-MORB) and ocean-island basalt (OIB), with a subset of samples indicative of continental arc affinity. These geochemistry results indicate that distinct tectonic environments are represented by at least two, and possibly three, lithological and structural units comprising the Seventymile terrane. Hf-Nd isotope systematics are consistent with a depleted MORB mantle (DMM)−like component that overlaps with Pacific MORB. Primary zircon is rare, but new in situ U-Pb data for gabbro and greenstone indicate ca. 274−272 Ma peak zircon and titanite crystallization. Scattered younger zircons define a ca. 255 Ma zircon peak and correspond to secondary crystallization associated with baddeleyite reaction of high-Si fluids during low-grade metamorphism. If Seventymile suites are contemporaneous, obduction associated with the closure of the SMO resulted in the stacking of ophiolitic packages representing distinct tectonomagmatic settings across the transition from pericontinental, to epicontinental, to distal ocean back-arc. Intrusions hosted in klippe of ultramafic rocks, plus the least subduction-modified greenstones underlying them, geologically and compositionally resemble Slide Mountain rocks of the Campbell Range formation in eastern Yukon and may provide a new piercing point across the Tintina fault.

Released January 05, 2026 09:12 EST

2026, Journal of Contaminant Hydrology (277)

Andrew H. Manning, Robert L. Runkel, Jean M. Morrison, Sara Warix, Richard B. Wanty, Katherine Walton-Day, Michael Snook

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

Released January 05, 2026 08:38 EST

2026, Newsletter

Katrina Rossos

No abstract available.

Released January 05, 2026 08:18 EST

2026, Remote Sensing (18)

Monica Palaseanu-Lovejoy, Jeffrey J. Danielson, Minsu Kim, Bryan Eder, Gretchen Imahori, Curt D. Storlazzi

This research focusses on three satellite-derived bathymetry methods and optical satellite instruments: (1) a stereo photogrammetry bathymetry module (SaTSeaD) developed for the NASA Ames stereo pipeline open-source software (version 3.6.0) using stereo WorldView data; (2) physics-based radiative transfer equations (PBSDB) using Landsat data; and (3) a modified composite band-ratio method for Sentinel-2 (SatBathy) with an initial simplified calibration, followed by a more rigorous linear regression against in situ bathymetry data. All methods were tested in three different areas with different geological and environmental conditions, Cabo Rojo, Puerto Rico; Key West, Florida; and Cocos Lagoon and Achang Flat Reef Preserve, Guam. It is demonstrated that all satellite derived bathymetry (SDB) methods have increased accuracy when the results are aligned with higher-accuracy ICESat-2 ATL24 track bathymetry data using the iterative closest point (ICP). SDB vertical accuracy depends more on location characteristics than the method or optical satellite instrument used. All error metrics considered (mean absolute error, median absolute deviation, and root mean square error) can be less than 5% of the maximum bathymetry depth penetration for at least one method, although not necessarily for the same method for all sites. The SDB error distribution tends to be bimodal irrespective of method, satellite instrument, alignment, site, or maximum bathymetry depth, leading to the potential ineffectiveness of traditional error metrics, such as the root mean square error. However, our analysis demonstrates that performing detrending where possible can achieve an error distribution as close to normality as possible for which error metrics are more diagnostic.

Released January 03, 2026 08:02 EST

2026, Ecotoxicology (35)

Vicki S. Blazer, Cheyenne Rose Smith, Heather L. Walsh, P.M. Mazik, M.R. Magee

The U.S. Great Lakes Areas of Concern (AOCs) have been designated due to significant environmental degradation, with multiple Beneficial Use Impairments (BUIs) requiring remediation before delisting. The “fish tumors or other deformities” BUI remains at the Sheboygan River and the Lower Green Bay and Fox River AOCs, necessitating evaluation in white sucker (Catostomus commersonii), a sentinel species. This study documented the prevalence of skin and liver neoplastic and preneoplastic lesions at these AOCs, assessed temporal trends at the Sheboygan River AOC from 2012 to 2021, and investigated potential risk factors influencing liver neoplasm development using histological and molecular approaches. Histopathological analyses revealed that skin neoplasms were present in 15.5% of fish from the Green Bay AOC and 21.5% from the Sheboygan River AOC, while liver neoplasms were detected in 7.0% and 6.0% of fish, respectively. At the Sheboygan River AOC, the prevalence of skin tumors declined from 32.6% in 2012 to 21.5% in 2021, and liver tumor prevalence decreased from 8.3% in 2012 to 6.0% in 2021. Hepatic transcript analyses identified significant expression of genes associated with contaminant exposure, oxidative stress, and immune response, suggesting ongoing exposure to environmental stressors. The presence of white sucker hepatitis B-like virus (WSHBV) in hepatocytes was confirmed, though no direct correlation between viral transcript abundance and neoplastic lesions was found. Findings indicate tumor prevalence has declined over time, however rates remain elevated compared to previous reference site assessments. The study underscores the potential multifactorial nature of carcinogenesis in wild fish populations, involving both chemical and biological stressors.

Released January 02, 2026 08:54 EST

2026, The Seismic Record (6) 1-12

Harriet Zoe Yin, Kate E. Allstadt, William D Barnhart, Samantha Ann Clapp, Paul S. Earle, Dara Elyse Goldberg, Alex R. Grant, Matt Herman, Kishor S. Jaiswal, Sara K. McBride, Adam T. Ringler, Max Schneider, Eric M. Thompson, Nicholas van der Elst, David Wald, Dun Wang, Charles Worden, William L. Yeck

The 29 July 2025 M_w 8.8 Kamchatka, Russia, earthquake was the sixth largest instrumentally recorded earthquake. This event was seismically well observed at regional and teleseismic distances, but publicly available near‐source data were sparse at the time of the event, presenting unique challenges for rapid source and impact characterization. The U.S. Geological Survey (USGS) National Earthquake Information Center provides global real‐time monitoring for earthquakes, including rapid response information products that estimate source characteristics, shaking, and the resulting impacts. We describe the USGS rapid response earthquake information products following the Kamchatka event and discuss their implications for ongoing hazards in the region. We describe potential improvements to our response workflows motivated by this event, including more rapid constraints on source geometries and the automated selection of fault geometries for finite‐fault inversions. The rapid response products together support the interpretation of a unilateral southwestward rupture with significant slip on the southwestern end of the rupture extent. The M_w 8.8–9.0 event in 1952, which ruptured a comparable extent of the Kuril–Kamchatka subduction interface, has many similarities to the 2025 rupture. This illustrates that slip deficits may remain following great earthquakes and highlights the usefulness of comparative studies between historic and modern events.

Released January 01, 2026 09:49 EST

2026, GSA Today (36) 4-11

Lauren N. Schaefer, Francis K. Rengers, Benjamin Mirus, Liam Toney, Kate E. Allstadt, Richard Wooten, Patrick Moore, Paula Madeline Burgi, Anne Witt, Eric Leland Bilderback, Jennifer Bauer, David Korte, Matthew Crawford

Between 23 and 27 September 2024, antecedent rain followed by Hurricane Helene produced one of the most damaging weather events in southern Appalachia history. The back-to-back storm events resulted in a maximum cumulative rainfall of 848 mm and hurricane-force wind gusts over 170 km/h in western North Carolina, eastern Tennessee, and southwestern Virginia. The resulting regional flooding, landslides, and tree blowdown caused over 100 fatalities, damaged or destroyed critical infrastructure and thousands of structures, and severed connectivity across the region. Over the next several weeks, a multi-agency landslide response produced a rapid hazard assessment and mapped 2217 landslides, 55% of which damaged infrastructure or property. Orographic uplift enhanced rainfall, resulting in concentrated landsliding along the ~250 km swath of the Blue Ridge escarpment in western North Carolina. Landslides initiated predominantly on windward-facing (southeast-facing) slopes, and localized clustering of initiation points indicated a strong influence of hillslope-scale meteorological and geomorphic factors. Many shallow landslides mobilized into larger, highly mobile, and damaging debris flows that graded into floods. Here, we put our preliminary observations in the context of historical storm-driven landslide events and open new avenues for investigating the nature and extent of landslides and their effects in southern Appalachia and similar environments.

Released January 01, 2026 09:36 EST

2026, Limnology & Oceanography Bulletin

Michael Frederick Meyer, Jorrit Mesman, Carolina C. Barbosa, Jonathan J Borelli, Johannes Feldbauer, Merritt Elizabeth Harlan, Robert T. Hensley, Burak Kuyumcu, Robert Ladwig, Isabella Oleksy, Rachel M. Pilla, Jacob Aaron Zwart, Matthew Biddle, Paul J. Kinzel, Carl J. Legleiter, Tadhg Moore, Tylar Murray, Lipa Nkwalale, Brandon Overstreet, Mehraz Rumman, Whitney M. Woelmer

No abstract available.

Released December 31, 2025 10:45 EST

2025, Bulletin of the Seismological Society of America

Andrew J. Michael, Andrea L. Llenos

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

Released December 31, 2025 10:31 EST

2025, The Seismic Record (5) 352-362

Edward W Woolery, William J. Stephenson, Kevin Woller, Alena L. Leeds, Noah Silas Lindberg, Jackson K. Odum, Cooper Cearley, Ron Counts

The Meadow Bank scarp (MBS) in southeastern Illinois is a linear geomorphic expression, ∼10 km long and ∼8 m high above a relatively flat landscape. It parallels an underlying northeast‐oriented Late‐Precambrian–Early‐Cambrian structural fabric, called the Wabash Valley fault zone, and is within an area of modern, historic, and paleo seismicity, called the Wabash Valley seismic zone. In addition, the MBS acts as a boundary of the Wabash River floodplain, as well as Pleistocene glacial outwash channels, which show evidence of frequent outburst flood events. To better understand the MBS’s equivocal origin in this complex geologic environment, we acquired a 917‐m‐long seismic‐reflection survey across its axis to assess the subsurface geologic configuration. The resultant image indicates a complex set of faults that antiformally fold and displace the top of Paleozoic bedrock by ∼12 m across the survey. Moreover, fault and/or fold deformation extends into the shallowest imaged Quaternary strata at ∼6 m below the ground surface. This suggests the MBS origin is related to underlying Quaternary reactivated, positively inverted faults rather than exclusively to glacial outburst flood erosion. These results provide rare paleoearthquake spatial constraints for central U.S. regional seismic hazard consideration.

Released December 31, 2025 07:02 EST

2025, Scientific Investigations Report 2025-5101

Molly Van Appledorn, Lucie Sawyer, John Delaney, Chanel Mueller, Leigh Youngblood, Jane Harrell, Brian Breaker, Chris Frans

The hydrology of the Upper Mississippi and Illinois Rivers is a fundamental driver of ecosystem patterns and processes across a large portion of the United States. Quantitative hydrologic data for the main stems of these rivers underlie numerous scientific investigations, statistical models, and decision-making processes for local, State, and Federal agencies involved in the Upper Mississippi River Restoration program. Although historical hydrologic data exist, data representing potential future conditions of the Upper Mississippi and Illinois Rivers lack the resolution necessary to anticipate biotic and abiotic responses to altered hydrology and to determine resilient management actions. A source of future hydrologic scenarios is the readily available LOCA–VIC–mizuRoute hydrologic data products (named for the chain of models the data are produced from—localized constructed analogs, Variable Infiltration Capacity macroscale hydrological model, and the mizuRoute hydrologic routing model—that we shorten further to LVM in this report) that include simulated discharges for historic and future timeframes. The objective of this study is to assess the reliability of the hydrologic data products for their use in Upper Mississippi River Restoration program applications. Key study questions are (1) do the hydrologic data products reproduce characteristics of hydrology necessary to support ecological modeling and restoration decision-making applications within the Upper Mississippi River Restoration program? and (2) are there geographic differences in the reliability of the hydrologic data products?

Seven characteristics of river hydrology were selected related to flow magnitude, seasonality, and regime for evaluation. The seven characteristics were calculated using observed and historical simulated hydrologic data at 19 U.S. Geological Survey streamgages throughout the basins of the Upper Mississippi and Illinois Rivers; two streamgages are located on the main stem of the Mississippi River and two streamgages are located on the main stem of the Illinois River. Statistical comparisons between observed and historical simulated characteristics indicated that the hydrologic data products did not reliably represent historical hydrologic conditions in the basin or main stem. The hydrologic data products we evaluated could not reliably capture the overall hydrologic regime or flow magnitudes; the latter is evidenced by substantial underestimates of discharge at most streamgages. Seasonal hydrologic characteristics were captured more reliably than flow magnitude, but overall correspondence was low for most streamgages. A weak latitudinal pattern in seasonal characteristics indicated the hydrologic data products poorly represent streamflow timing in snow-affected regions of the basin. Discrepancies in magnitude, seasonality, and regime indicate the potential for multiple sources of error. Because poor correspondence was present across all 19 streamgages, it was not possible to identify specific drivers of poor performance (that is, drainage area or geography). The modeling chain should be evaluated for biases associated with meteorologic forcing data, as well as hydrologic model formulation and calibration.

We conclude that the hydrologic data products we evaluated appear unsuitable for applications tied to habitat and ecosystem restoration and management in the Upper Mississippi and Illinois Rivers. Plans to develop a future hydrology dataset for the Upper Mississippi River Restoration program would benefit from ongoing work to improve global climate model output downscaling methods, to improve hydrologic models, to make use of innovations in machine-learning approaches for projecting hydrology, and other efforts. The framework developed herein to evaluate hydrometeorological outputs generated using global climate models for a specific water resources application is a transferrable approach that could be applied to other data products and river systems.

Released December 30, 2025 09:46 EST

2026, Water Resources Research (62)

Barbara Bekins, William Herkelrath

At many petroleum hydrocarbon spill sites, residual spilled product forms a long-term source of groundwater contamination. The phrase source zone natural depletion is used to refer to the mass loss rates. Overall mass lost under environmental conditions was analyzed using conservative biomarker concentrations for a 1979 oil spill in northern Minnesota, USA. After 40–41 years, an average of 50% of the mass was lost with values ranging from 22% to 57% depending on location. It is also important to understand the composition changes in the source. To understand controls on the losses of individual compounds, concentrations of volatile hydrocarbons in oil samples were compared with aqueous solubilities, and pore-space oil saturations. The results of the comparison show that losses of the oil compounds were controlled by pore-space oil saturations, solubility, and susceptibility to degradation under methanogenic conditions. Compounds that degrade under methanogenic conditions, including toluene, o-xylene, and n-alkanes are more depleted compared to benzene, ethylbenzene, and m- and p-xylene for which losses are dominated by dissolution. These rates and compound-specific behaviors form a foundation for improved modeling approaches and risk analyses.

Released December 30, 2025 09:41 EST

2026, North American Journal of Fisheries Management

Jordan Mae-Jean Buxton, Marissa L. Wulff, Brock Huntsman, Kai Palenscar, Brett Mills, Kerwin Russell, Alicia Ruan, Tevin Bui

Released December 29, 2025 10:13 EST

2025, Western North American Naturalist (85)

Jeffrey E. Lovich, George T. Jefferson

he North American porcupine (Erethizon dorsatum) is widely distributed throughout many ecosystems on the continent from northern Mexico to the Arctic Ocean, except for the southeastern portion of the United States. Habi­tats include the arid Desert Southwest region where modern records are generally sparse. The paleogeographic range in the Desert Southwest is not unlike the modern distribution of Erethizon, with some exceptions. Although there are early Pleistocene records of Erethizon for southern California, no late Pleistocene records (Rancholabrean North American Land Mammal Age [NALMA]) are known for Imperial, Inyo, Riverside, or San Bernardino Counties in southern California or in Mexico. This is surprising given late-Pleistocene records in Arizona, Nevada, and New Mexico and modern records from southern California and Mexico. Pleistocene habitats for porcupines were similar to those occupied today, just at different elevations due to differing climates. Modern preferred habitats in the Desert Southwest include conifer-clad high-elevation areas and associated pinyon pine / juniper belts, and occasionally riparian corridors in desert environments. The overall rarity of porcupines today in arid parts of the Desert Southwest is likely due to the combination of increasing aridity in the region during the Holocene, past persecution by humans, and increasing predator populations.

Released December 29, 2025 09:56 EST

2025, Ecology and Evolution (15)

Jennifer A. Grauer, Jacqueline L. Frair, Krysten L. Schuler, Manigandan Lejeune, David W. Kramer, Angela K. Fuller

Parasite-induced morbidity and mortality can alter the trajectories of incidental host populations. Yet, parasites rarely act in isolation and may be one of a multitude of biotic and abiotic stressors that collectively shape mortality risk in vertebrate populations. We quantified sources of mortality in a low-density population of moose (Alces alces) in New York State and investigated factors including parasite infection, nutritional limitation, and thermal stress influencing mortality risk in calf moose. We observed high rates of annual survival (0.81–0.92) in adult (n = 25) and calf (n = 27) moose monitored 2015–2018 and 2022–2024, respectively. The majority of cause-specific mortality was attributed to disease induced by giant liver fluke (Fascioloides magna; 75% in adults, 67% in calves). Calf mortality risk increased by 72% for every unit increase in giant liver fluke infection risk, measured as cumulative monthly proportion of wetlands used by moose, and decreased by 16% with each additional unit of nutritional energy available. The combination of flukes, coinfecting parasites, and available nutritional energy is important to calf survival in this population, highlighting the importance of managing multiple stressors for species conservation, although the effects are hard to disentangle given the high rates of survival observed. Identifying causes of mortality and mechanisms underlying increased mortality risk contributes to the continued conservation of moose in fluctuating populations and highlights the importance of managing parasite-induced disease.

Released December 29, 2025 09:51 EST

2025, California Fish and Wildlife Journal (111)

Jordan Mae-Jean Buxton, Marissa L. Wulff, Veronica L. Violette, Danielle L. Palm, Matthew J. Young, Frederick Feyrer

No abstract available.

Released December 29, 2025 08:52 EST

2025, Conservation Science and Practice

Joshua B. Mouser, Ashley A. Dayer, Serena Ciparis, Sara Bottenfield, Paul L. Angermeier

Agricultural conservation practices are often used to protect stream health while continuing food production. However, recovery of stream health is often not as rapid or extensive as planned. The efficacy of practices may be improved by promoting their continued use by landowners (i.e., persistence) after cost-share contracts with government agencies end. Persistence rates and their drivers are typically unknown. Therefore, we mail-surveyed 889 landowners to understand persistence in grazing practices in southwest Virginia. Survey responses were analyzed quantitatively using logistic regression and qualitatively via coding. We found that persistence rates for vegetative practices, cattle-exclusion fencing, off-stream watering structures, and pasture management were 74%, 84%, 94%, and 94%, respectively. Both the quantitative and qualitative analyses indicated that landowner cognitions (i.e., attitudes toward practices, environment, and agencies), environmental motivations, and practice durability influence persistence. Our results highlight that persistence could be encouraged by: (1) providing targeted messaging that demonstrates persistence benefits and aligns with landowners' motivations, (2) ensuring that responsibility to maintain practices is transferred during land-tenure changes, and (3) allocating more agency funding to practice maintenance.