USGS Publications Warehouse

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

Thu, 28 May 2026 14:13:31

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

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

Tue, 26 May 2026 18:25:26

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

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

Laske, Sarah; Young, Daniel; Bartz, Krista; von Biela, Vanessa; Carey, Michael — Wed, 27 May 2026 14:38:46

Objective

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

Methods

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

Results

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

Conclusions

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

Ecology of reintroduced Rocky Mountain bighorn sheep in Dinosaur National Monument

Wed, 27 May 2026 14:33:30

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

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

Lopez, Brianna; Shelton, Jenna; Marketti, Michael; Ritzel, Kate; Graham, Brandon — Tue, 26 May 2026 18:37:48

Introduction

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

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

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

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

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

Sweetkind, Donald S. — Tue, 26 May 2026 18:39:21

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

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

Harkness, Jennifer; Faulkner, Kirsten; Jurgens, Bryant — Tue, 26 May 2026 18:29:13

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

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

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

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

Atkinson, Allison — Tue, 26 May 2026 18:19:14

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

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

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

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

Remote sensing enables basin-scale inventories of coal mine methane

Wed, 20 May 2026 14:51:11

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

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

Tue, 26 May 2026 21:07:13

Executive Summary

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

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

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

System characterization report on Tanager

Kim, Minsu; Park, Seonkyung; Clauson, Jeff; Vrabel, Jim; Sampath, Ajit — Tue, 26 May 2026 13:18:54

Executive Summary

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

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

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

Carson, Jennifer N.; Benacquisto, Matthew T. — Wed, 27 May 2026 13:00:48

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

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

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

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

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

Tue, 19 May 2026 14:46:14

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

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

Mon, 18 May 2026 15:28:47

Study region

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

Study focus

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

New hydrological insights for the region

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

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

Clements, Kaylin; English, James; Wilkins, Emily; Moore, Megan; Schuster, Rudy — Thu, 14 May 2026 14:21:59

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

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

Fri, 15 May 2026 17:52:16

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

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

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

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

Tue, 19 May 2026 15:19:53

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

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

Medler, Colton; Anderson, Todd — Fri, 15 May 2026 17:44:04

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

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

Mon, 11 May 2026 17:07:27

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

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

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

Mon, 18 May 2026 15:50:01

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

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

Blommel, Caroline; Lamont, Margaret; Kendall, William — Wed, 20 May 2026 14:57:35

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

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

Russell, Aaron Daniel; Bair, Lucas; Meldrum, James; Hawbaker, Todd — Fri, 15 May 2026 13:16:54

Background

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

Aims

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

Methods

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

Key results

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

Conclusions

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

Interpretations

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

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

Wed, 13 May 2026 14:23:56

Context

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

Objectives

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

Methods

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

Results

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

Conclusions

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

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

Legleiter, Carl; Kinzel, Paul; Laker, Mark; Conaway, Jeff — Tue, 12 May 2026 13:47:57

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

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

Kenworthy, Megan — Mon, 11 May 2026 17:06:06

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

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

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

Water use in Louisiana, 2020

Robinson, Angela — Mon, 11 May 2026 17:02:38

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

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

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

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

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

Hillenbrand, Ian; Thomson, Kelly — Tue, 26 May 2026 14:38:17

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

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

Wed, 13 May 2026 14:10:53

Objective

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

Methods

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

Results

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

Conclusions

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

Los Planes watershed vegetation monitoring: Standard operating procedures

Wilson, Natalie R. — Mon, 11 May 2026 13:48:50

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

Analysis of alternative weir designs for improved passage of select fish at the U.S. Geological Survey streamgaging weir at Blackwells Mills, New Jersey

Suro, Thomas; Niemoczynski, Michal; Mulligan, Kevin — Mon, 11 May 2026 16:59:37

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

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

Regional conservation planning tool: A spreadsheet model to support spatial prioritization and resource allocation decisions

Fri, 8 May 2026 14:24:23

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

USGS 2025 critical minerals review

Jones, James; Gallegos, Tanya; Kunledare, Mojisola; Riggs, Charlotte — Mon, 18 May 2026 15:39:31

The United States Geological Survey (USGS) provides scientific information for the Department of Interior and the nation, consistent with its original mission expressed in the Organic Act of 1879 (43 U.S.C. 31): “the classification of the public lands and examination of the geological structure, mineral resources, and products within and outside the national domain.” Legislation such as the Energy Act of 2020 and the 2022 Infrastructure Investment and Jobs Act (43 USC 31l) and recent executive actions (Executive Orders 14154 , 14153, 14241 Secretary’s Orders 3417, 3418, 3422, 3436) underscore the importance of mineral resources and focus USGS activities on mapping and assessing mineral resources, with a particular focus on those presently identified as critical, both in ground and above ground in mine wastes.

This article reviews selected activities and accomplishments by the USGS Mineral Resources Program related to critical minerals in 2025. Highlights include a new List of Critical Minerals, a first-ever national mine waste inventory, international minerals partnerships, mineral resource assessment advancements, and national data collection activities and outcomes of the Earth Mapping Resources Initiative (Earth MRI). The selected contributions are not comprehensive but are intended to demonstrate USGS leadership in critical mineral mapping and assessment, the importance of domestic and global partnerships, and the breadth of research activities that are responsive to national needs and priorities.

An overview and participatory framework for choosing spatial boundaries in social–ecological systems modeling

Thu, 7 May 2026 15:15:49

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

Revisiting the utility of regional-scale, high-quality geophysical data in mineral exploration - A case study featuring the Mammoth Magnetic Anomaly, Pinal County, Arizona

Walter, Callum — Mon, 18 May 2026 15:10:18

Regional aeromagnetic surveys passively measure the total magnetic intensity (TMI) and are a foundational tool used in mineral exploration (Airo, 2015). With the increased global demand and the number of critical mineral resources required for manufacturing high-tech devices, developing high-quality, regional-scale geophysical surveys could aid critical mineral exploration efforts and geologic mapping. In 2019, the U. S. Geological Survey launched the Earth Mapping Resources Initiative (Earth MRI) to modernize the geologic and geophysical mapping of regions that have the potential to contain critical mineral resources within the United States. In support of planning Earth MRI geophysical surveys, Drenth and Grauch (2019) defined five aeromagnetic data quality rankings (rank 1 through rank 5) applying them to the airborne geophysical survey inventory of the United States (Johnson et al., 2021). Rank 1 aeromagnetic surveys are of the highest quality, meeting modern standards and allowing best practices for qualitative and quantitative interpretation; whereas rank 5 aeromagnetic surveys are of the lowest quality, being useful only for qualitative interpretation of broad features. Through the Earth MRI effort, 48 high-quality, regional-scale rank 1 and 2 airborne magnetic and radiometric geophysical surveys have been planned, collected, or publicly release through May 2025 (U. S. Geological Survey, 2025). Here, a portion of a rank 1 Earth MRI aeromagnetic survey in southeast Arizona is presented and compared to a legacy rank 5 aeromagnetic survey over the Mammoth Magnetic Anomaly (MMA), demonstrating how modern, high-quality aeromagnetic data improves our view of crustal geology, aiding mineral exploration.

The United States Magnetotelluric Array and the National Impedance Map

Fri, 1 May 2026 14:26:55

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

Understanding the occurrence and distribution of per- and polyfluoroalkyl substances (PFAS) in surface waters of the nontidal Passaic River Basin

Fri, 1 May 2026 16:47:30

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

Assessment of long-term trends in streamflow statistics within and near the Mobile Bay and Perdido Bay watersheds, United States, 1950–2022

Asquith, William; Crowley-Ornelas, Elena; Whaling, Amanda — Fri, 1 May 2026 16:45:34

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

Evaluating approximations of river channel shape using a national cross section database

Legleiter, Carl; Kinzel, Paul — Fri, 1 May 2026 14:00:47

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

Changes in suspended sediment concentration along tidal rivers of the Chesapeake Bay: The tidal freshwater “sediment shadow”

Noe, Gregory; Murphy, Rebecca; Krauss, Ken — Fri, 15 May 2026 13:15:07

Transport of terrigenic sediment from nontidal watersheds into estuaries has important impacts on coastal habitat quality, pollutant transport, and resilience to sea-level rise. However, relatively little is known about changes in suspended sediment as nontidal rivers encounter tide, transition into tidal rivers through the tidal freshwater zone (TFZ), and enter saline portions of estuaries. The goal of this paper is to identify spatial and temporal patterns in suspended sediment concentration (SS) changes across tidal and salinity gradients over multiple tidal rivers, using a robust monitoring long-term dataset from the Chesapeake Bay. The multiple TFZs in the Chesapeake Bay consistently have a “sediment shadow” shown by a local spatial minimum in SS compared to upstream nontidal and downgradient oligohaline river reaches. Similarly, freshwater inputs from nontidal rivers have diminishing influence on tidal SS temporal dynamics with distance downstream from the head-of-tide. Therefore, little of the contemporary watershed sediment load is likely transported past the TFZ except during extreme floods when some sediment may be delivered to saline portions of the estuary. Tidal freshwater and brackish portions of the estuary have spatially variable trends in SS over time, both increases and decreases. However, the more saline downstream ends of tidal rivers and the mainstem of the Chesapeake Bay have had a consistent average 25% decline in SS over the past decades. In summary, the presence of “sediment shadows” suggests watershed loads of sediment are currently mostly not transported through the TFZ into the saline estuary, and likely generate sediment deficits for tidal freshwater wetlands.

Hydrologic investigation of water level fluctuations at Moreau Lake, Moreau Lake State Park, town of Moreau, New York

Heisig, Paul M. — Fri, 1 May 2026 16:44:17

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

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

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

Aligning legacy NLCD land cover maps based on Landsat Collection 1 to Collection 2

Li, Congcong; Jin, Suming — Thu, 21 May 2026 14:40:45

The transition from Landsat Collection 1 to Collection 2 introduced significant improvements in radiometric and geometric accuracy. However, the improvements cause location misalignment between the existing Landsat-derived land cover products and the new collection. The legacy National Land Cover Database (NLCD) has been used as a cornerstone land cover source for a variety of research. Therefore, a method aligning the legacy NLCD product to Collection 2 is required to ensure its continuity and consistency of service. We developed a strategy to not only align legacy NLCD to match new Collection 2 geometric locations but also improve land cover labeling in the region that was affected by the geometric shifts. The method identifies boundary pixels of homogeneous land cover patches as potential problem areas that are likely impacted by geometric shifts and generates candidate labels from 3 × 3 window with the target pixel at the center and segmentation-derived majority label. Standard phenology patterns of each candidate land cover type are established based on the random samples except boundary pixels within a 1000-pixels × 1000-pixels processing window region. The phenological distance to each standard land cover type pattern is calculated through a penalty dynamic time warping (DTW) method for each target pixel in the boundary region. Finally, the method determines the most suitable label based on the phenological distance from the candidate labels. Both visual and accuracy assessment results demonstrate that the alignment preserves the overall land cover patterns in the original legacy NLCD product while reducing the spatial discrepancies between the Landsat Collection 2 and land cover. In addition, it enhances the accuracy of land cover labeling of boundary pixels. The overall accuracy (OA) was increased by 7% in the land cover boundary regions after alignment. The quality and confusion matrix comparison between the alignment results and the original legacy NLCD confirm the reliability of the method. Our alignment method has the potential to serve as a framework for aligning other Landsat-derived land cover products to future collections.

Shallow hydrogeologic framework of the Tully Valley mudboil area, Onondaga County, New York

Fri, 1 May 2026 16:42:15

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

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

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

Springtime formation of laminated soil carbonate rinds and changes in fluvial terrace soils on orbital timescales at Rio Mesa, Utah, USA

Fri, 1 May 2026 13:21:06

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

Modeling chronic wasting disease transmission risk in mule deer related to habitat characteristics

Fri, 1 May 2026 13:56:46

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

Characterizing the long-term (1981–2023) temperature and precipitation dynamics in the Trans-Mountain regions of Kazakhstan, Central Asia

Wed, 29 Apr 2026 14:27:38

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

Peak-, mean-, and low-streamflow regional-regression equations for natural streamflow in central and western Colorado, 2019

Kohn, Michael; Mast, M.; Gross, Tara — Mon, 18 May 2026 16:07:47

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

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

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

Uranium—Deposits, production and resources, market dynamics, and supply chain risks

Mihalasky, Mark — Fri, 24 Apr 2026 18:39:05

Introduction

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

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

Cyanobacterial bloom occurrence and emergency department visits for asthma or wheeze, Wisconsin, 2017–2019

Thu, 7 May 2026 15:56:42

Background:

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

Methods:

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

Results:

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

Conclusion:

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

Development of a two-stage lifecycle model to inform the trap-and-haul program for Oncorhynchus kisutch (coho salmon) in the Lewis River, Washington

Plumb, John; Perry, Russell — Thu, 23 Apr 2026 13:56:24

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

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

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

Fish body midline segmentation using binary search

Thu, 23 Apr 2026 15:07:26

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

A site prioritization tool for invasive species management: Integrating diverse spatial data to improve decision making

Tue, 5 May 2026 14:55:27

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

Non-native invasive beetle alters structure of a riparian bird community in a biodiversity hotspot

Kus, Barbara E.; Yee, Julie; Mendia, Shannon M. — Thu, 30 Apr 2026 15:17:01

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

Standardized method for logging drill core at the Idaho National Laboratory, Idaho

Dietz, Haley — Tue, 21 Apr 2026 14:52:03

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

Modeling future groundwater depletion to evaluate sustainability goals set under the Sustainable Groundwater Management Act in the critically overdrafted basins of the Central Valley, California, USA (2020–2070)

Wed, 20 May 2026 14:13:52

In 2014, California's Sustainable Groundwater Management Act (SGMA) mandated local agencies to devise and implement groundwater sustainability plans to address critically overdrafted conditions throughout the state's aquifers. However, the feasibility of these agencies' sustainability goals has not previously been assessed through a regional-scale, integrative lens. Here, we develop and analyze a novel, basin-wide database of 936 sustainability indicator wells located within Central Valley subbasins designated as critically overdrafted, most of which lie in the San Joaquin Valley. Our database shows 2040 groundwater elevation goals vary widely from 60 m above to 80 m below 2020 levels, with variability within and between adjacent subbasins. To evaluate the feasibility of achieving these goals, we coupled the database with a regional hydrologic model (Central Valley Hydrologic Model version 2) and simulated multiple future pumping scenarios. Results show that under increased groundwater demand, 60%–70% of indicator wells may fail to meet their 2040 goals. Even a 50% reduction from 2020 demand levels leaves nearly 40% of wells failing to meet their sustainability thresholds by 2040. Baseline models show that by 2070, up to 70% of wells could fail to meet their goals due to large-scale, spatially connected regions of groundwater depletion. This integrated framework, linking the first region-wide compilation of SGMA indicator wells with a regional groundwater model, demonstrates that many local sustainability goals may be unattainable with substantial (up to 50%) reductions in pumping. Additional management interventions, such as expanded recharge or coordinated demand reductions, may help achieve sustainability goals.

Assessment of groundwater quantity and quality contributions to Lake Huron

Kaemming, Bridget; Ford, Chanse; Martin, Sherry — Fri, 24 Apr 2026 18:36:33

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

Thematic accuracy assessment of the National Land Cover Database (NLCD2021) for the conterminous United States

Tue, 19 May 2026 14:40:53

The MultiResolution Land Characteristics (MRLC) consortium’s National Land Cover Database (NLCD) was established as an operational, stakeholder-oriented land cover monitoring program with the release of the NLCD2001 database. Both land cover and land cover change data became available with the release of NLCD2006. Here, we report the land cover and land cover change accuracies for NLCD2021, documenting Level II and I accuracies for the 2016, 2019, and 2021 land cover datasets and 2016–2021 land cover change. An additional objective, not addressed in previous NLCD assessments, is the estimation of the total variance of area estimates, where total variance includes sampling variance and variability in reference label assignment. We estimate the total variance of the area estimates for the stable and change classes targeted by the strata used in the sampling design. For the 16 Level II classes, the overall accuracy (OA) was 73% ± 1% when agreement was defined using only the primary reference class label and 89% ± 0.7% when agreement was defined as a match between the map label and either the primary or alternate reference label (± standard error, SE). For the eight Level I classes, the corresponding OA estimates were 81% ± 0.9% and 89% ± 0.7%. Level II and Level I OA tended to be higher for 2016 land cover and lower and more similar for 2019 and 2021 land cover. Land cover change user and producer accuracies (UA and PA) tended to be <50%. When the alternate reference label was included in the definition of agreement, exceptions of high class-specific accuracies (≥75%) were more prevalent for PA than UA. The consistency of reference label interpretations was greater for no-change classes than for change classes. The estimates of the total variance incorporating interpreter variability were often smaller than the standard variance estimates, indicating the possibility that the total variance estimator is highly unstable. Further study is needed to improve the utility of this total variance estimator for practical applications.

Quantitative mineral resource assessment of lithium pegmatite deposits in the northern Appalachian orogen, USA

Mon, 20 Apr 2026 15:29:50

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

A Bayesian hierarchical modeling approach for species diversity in ecology

Qian, Song; Dufour, Mark; Jaffe, Sabrina; Hilling, Corbin; Hintz, William — Fri, 24 Apr 2026 14:13:47

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

A practical decision tool for marine bird mortality assessments

Thu, 14 May 2026 14:01:51

Given the rise in anthropogenic, environmental, and disease events contributing to marine bird mortality, there is a critical need to improve the rigor of mortality assessments. Deficits in data collection and mortality estimation can hinder a manager’s ability to document the scale of events and assess population level impacts. Therefore, to inform decisions required during activities, such as conservation status assessments or harvest management, organizations may choose to incorporate mortality assessments into response plans. Resources, capacity, and assets to assess mortality vary across jurisdictions (federal, state, Indigenous, local, etc.), and clear guidance to support mortality assessments is often unavailable or not clearly addressed. Here, we present a decision support tool to help managers identify and evaluate survey options to assess bird mortality in a diverse array of scenarios. The objective of the decision tool is to improve data collection and availability, which will increase the ability to estimate mortality robustly, given situation-specific attributes and constraints. This decision tool is designed to guide the response when a mortality event is initially encountered and offers suggestions for assessment and reporting procedures in the absence of other guidance or to complement existing protocols. The decision tool is also meant to inform decision making for response determination and resource allocation. The tool facilitates examination of options for further assessment and monitoring, which users determine by examining questions pertaining to species prioritization, determination of mortality minimum spatial extent, and the potential magnitude of impacts on affected species. Finally, identification of appropriate survey methods that address imperfect detection when a complete census is not possible are determined by exploring location, spatial and temporal extent, and the type of species affected. Ultimately, this decision tool aims to facilitate and improve the standardization of mortality assessments, equipping managers with a practical resource to navigate the decision-making process for marine bird mortality estimation.

Hydrogeology, groundwater salinity distributions, and assessment of the effect of oil-production activities on groundwater in the Midway Valley area, western Kern County, San Joaquin Valley, California

Fri, 24 Apr 2026 14:10:27

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

Historical ice jams and associated environmental conditions on Osoyoos Lake

Sutfin, Nicholas; Breen, Stephen — Mon, 20 Apr 2026 17:42:39

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

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

Comparative assessment of STIC sensors, streamflow and rain gauges for quantifying river connectivity in intermittent systems

Cooper, Cienna; Rogosch, Jane; Smith, Nathan; Robertson, Clinton; Wilson, Wade — Mon, 4 May 2026 16:56:07

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

Logical data model for hydrographic data based on HY_Features concepts

Mon, 18 May 2026 14:33:53

This report describes background and design of the “hydrofabric data model” which defines logic for implementation of data schemas and software that deals with hydrologic geospatial data. As a “logical” data model, the hydrofabric data model specifies details necessary to support compatibility of data and software that satisfy diverse needs without unnecessarily restricting implementation details. The logic presented in this report is based on concepts defined in WaterML2 Part 3 Surface Hydrology Features Concepts and is designed to serve the needs of a range of hydroscience use cases.

Development of international community standards applicable to hydrofabrics began, prompted by the World Meteorological Organization Commission for Hydrology, in 2012 [5] . More than 10 years later, this report documents one aspect of a long-term research and development activity that traces its roots back that far.

This report describes terminology, use cases, and background as context preceding presentation of the logical model and discussion of its design. Three appendices document related data models, an example encoding of the hydrofabric data model, and an artificial schematic and tabular data example. The sections of the report can be accessed in the Clause 5 section.

Characterizing changes in postfire debris-flow hazard as burned areas recover

Fri, 17 Apr 2026 15:26:10

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

Analyses of meteorological and hydrological records support Tribal members’ accounts of changing climate on the Fort Apache Reservation, east–central Arizona

Mason, Jon — Mon, 20 Apr 2026 17:40:45

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

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

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

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

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

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

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

The global proliferation of aquatic, benthic Microcoleus: Taxonomy, distribution, toxin production, ecology, and future directions

Wed, 11 Feb 2026 15:20:31

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

Assessment of natural gas pipeline construction on stream temperature and turbidity in southwestern Virginia, 2017–25

Foster, Brendan; Maas, Carly; Flota, Alejandra — Mon, 27 Apr 2026 18:42:46

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

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

Proceedings of the Floodplain Vegetation Monitoring Workshop for the Long Term Resource Monitoring Element of the Upper Mississippi River Restoration Program, January 7–8, 2025, Moline, Illinois

Wed, 15 Apr 2026 14:24:34

Preface

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

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

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

Albuquerque Seismological Laboratory strategic vision

Wed, 15 Apr 2026 14:59:00

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

Opportunities and challenges in using Solid Phase Adsorption Toxin Tracking (SPATT) samplers for monitoring cyanotoxins in freshwater and estuarine environments

Tue, 28 Apr 2026 17:26:16

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

The Sedimentary Geochemistry and Paleoenvironments Project Phase 2 data release: An open data resource for the study of Earth's environmental history

Wed, 22 Apr 2026 15:01:35

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

Two-stage approach to automatic detection with machine learning for improved surveillance of the invasive Cuban treefrog

Huber, Kaitlin; Waddle, J. Hardin; Glorioso, Brad; Donovan, Therese — Mon, 4 May 2026 15:33:56

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

Initial condition uncertainty exerts a large and persistent influence on model simulations of ecosystem carbon dynamics in California

Selmants, Paul; Sleeter, Benjamin M.; Daniel, Colin — Wed, 22 Apr 2026 14:41:53

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

The Climate Hazards Center Infrared Precipitation with Stations, version 3

Tue, 19 May 2026 15:42:18

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

Occurrence of cyanobacteria and associated cyanotoxins in the Raritan Basin Water Supply Complex, New Jersey, August 2020 to August 2021

Mon, 27 Apr 2026 15:32:44

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

Arsenic and isotope concentrations in the lower Platte River valley of eastern Nebraska, early 1970s to 2023

Moser, Matthew; Cherry, Mikaela; Hall, Brent — Thu, 16 Apr 2026 17:22:01

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

Escherichia coli monitoring and assessment in 2022 and 2023 after beach restoration at Lake St. Clair Metropark Beach, Macomb County, Michigan

Lockmiller, Hayden; Byers, Victoria; Fogarty, Lisa R. — Thu, 16 Apr 2026 17:23:22

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

Incorporating data sets with multiple sources of uncertainty in integrated species distribution models

Lunt, Fiona; Scher, C.; Mummah, Riley; Miller, David — Tue, 14 Apr 2026 16:29:45

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

Estimating GPS-based social aggregation metrics using collar data

Fri, 17 Apr 2026 15:45:55

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

Semantic segmentation of light-toned veins in multimodal ChemCam data

Tue, 14 Apr 2026 15:24:12

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

Selected water-quality data from the Cedar River and Cedar Rapids well fields, Cedar Rapids, Iowa, 2017–22

Meppelink, Shannon; Kalkhoff, Stephen — Fri, 10 Apr 2026 18:18:25

The Cedar River alluvial aquifer is the source of drinking water in Cedar Rapids, Iowa. Production wells are completed in the alluvial aquifer approximately 40 to 80 feet below land surface. The City of Cedar Rapids and the U.S. Geological Survey have studied the groundwater-flow system and water quality of the aquifer in the vicinity of Cedar Rapids since 1992. Results of these studies documented hydrologic conditions, water quality, and geochemistry of the alluvial aquifer and interactions with the Cedar River. Water-quality samples were collected for studies involving well field monitoring, trends, source-water protection, groundwater geochemistry, surface-water–groundwater interaction, and pesticides in groundwater and surface water. Water quality was analyzed for dissolved major ions (boron, bromide, calcium, chloride, fluoride, iron, magnesium, manganese, potassium, silica, sodium, sulfate, and total dissolved solids), dissolved nutrients (ammonia as nitrogen, ammonia plus organic nitrogen as nitrogen, nitrite plus nitrate as nitrogen, nitrite as nitrogen, orthophosphate as phosphorus, and phosphorus), dissolved organic carbon, and selected pesticides. Physical characteristics (alkalinity, dissolved oxygen, pH, specific conductance, and water temperature) were measured on site and recorded for each water sample collected. This report presents the results of routine water-quality data-collection activities from October 2017 through September 2022. Methods of data collection, quality assurance, water-quality analyses, and statistical procedures are presented. Data include the results of water-quality analyses from quarterly sampling from monitoring wells, production wells, two water treatment plants, and the Cedar River at Blairs Ferry Road at Palo, Iowa, streamgage (U.S. Geological Survey station number 05464420), as well as monthly nutrient sampling from the Cedar River and Morgan Creek near Covington, Iowa, streamgage (U.S. Geological Survey station number 05464475).

Abundance and movement patterns of fish accessing a reconnected Lake Erie coastal wetland: Insights from high-resolution sonar data

Fri, 17 Apr 2026 17:19:20

Coastal wetlands of the Laurentian Great Lakes are complex ecosystems that provide essential biological services, including providing habitat for a suite of fish species. As restoration efforts for these coastal wetlands increase, there is a need to characterize how well restored areas support the life histories of wetland, riverine, and lake fishes. Most traditional survey methods (e.g. 24hr trap netting, visual surveys) are implemented over periods of time that make it difficult to describe short-term fluxes in fish activity, especially in waters with difficult access or high turbidity. To address this gap, we used acoustic sonar technology to evaluate fish movement in a recently reconnected coastal wetland on the southern shore of western Lake Erie. Data collected over four years (2011-2014) indicated that the 10 ha wetland was continuously utilized by millions of fish, with discernible fluctuations in usage patterns observed on both daily and annual scales. These insights add to our understanding of how fish assemblages respond to restored coastal wetland habitats and can inform management decisions that may impact fish access (e.g. conduct management activities that may be disruptive to fish populations during periods of reduced utilization). Additional study of short-term fish movements using high-resolution sonar and other technologies will reveal patterns that may enhance the effectiveness of restoration and management efforts in Great Lakes coastal wetlands.

Mineral chemistry perspective on remobilization of stored magma at Kamakai'a Hills, Southwest Rift Zone of Kilauea, Island of Hawai'i, USA

Downs, Drew; Sas, May — Mon, 13 Apr 2026 15:00:24

Differentiated magmas stored in the rift zones of Kīlauea have received more attention in recent years following eruption of andesite during the early phase of 2018 lower East Rift Zone activity. Despite this growing interest, some of the most voluminous eruptions of differentiated rift zone magmas remain poorly studied. One such eruption, and the most voluminous exposed differentiated flow field at Kīlauea, is the Kamakaiʻa Hills. This eruption took place in the Southwest Rift Zone of Kīlauea, a region that is hypothesized to contain a long-lived rift zone reservoir. The Kamakaiʻa Hills flow field encompasses >250 × 10⁶ m³ of basaltic andesite and basalt compositions with a mineral assemblage of orthopyroxene + clinopyroxene + plagioclase during its early ʻaʻā phase and clinopyroxene + plagioclase + olivine during its late pāhoehoe phase. To better understand storage conditions and magma accumulation, this study focuses on major, minor, and trace elements from the mineral assemblage present within the early ʻaʻā and late pāhoehoe phases. The diversity of clinopyroxene and plagioclase compositions within the early ʻaʻā and late pāhoehoe phases, as well as diverse compositions of plagioclase and orthopyroxene within the early ʻaʻā phase, suggest multiple magma bodies and limited pre-eruption magma mixing within the broader Kamakaiʻa Hills reservoir. Oscillatory zoning patterns (particularly in clinopyroxene) imply processes such as recharge events, magma mixing or mingling, or convection within a differentially cooling, chemically stratified reservoir over protracted time intervals, whereas only limited resorbed mineral textures indicate incomplete mixing of heat and chemically distinct magmas during the dike intrusion that triggered the eruption. Mineral-mineral and mineral-melt thermobarometry indicate predominantly shallow (≤2.5 km depth) crustal storage conditions of the cooled, differentiated magma (∼1100 °C and cooler for the basaltic andesites) to hotter temperatures for the basalts (all >1100 °C). Despite the known large standard errors estimated for mineral-melt and mineral-mineral barometry (10s to >100 MPa), the calculated pressures and depths broadly correspond with earthquake swarm depths beneath the Kamakaiʻa Hills, and drill core and fluid inclusion barometry storage depths of differentiated magmas within the lower East Rift Zone. The Kamakaiʻa Hills differentiated magmas have H₂O contents (∼0.5 wt%, using plagioclase-melt hygrometry) equivalent to typical Kīlauea basalts. Our data and interpretations demonstrate a complex, long-lived rift zone storage system that consisted of multiple magma bodies and was mobilized into eruption through intrusion of a hotter and more primitive summit-derived (uprift) magma.

Towards affordable wetland evapotranspiration monitoring using the Variance-Bowen Ratio method: Insights from three contrasting wetlands

Thu, 9 Apr 2026 13:30:55

Accurate measurement of evapotranspiration (ET) is essential for sustainable water management. Standard methods such as eddy covariance (EC) are costly, while alternatives such as surface renewal are cheaper but require calibration and complex data processing. This study evaluates the utility of the Variance-Bowen Ratio (VBR) method for estimating ET across three California’s wetlands. Using data from 2023, half-hourly latent heat flux (λE) and daily/monthly ET from VBR were compared with EC at one non-tidal (site A) and two tidal (sites B and C) wetlands. λE_VBR consistently underestimated λE_EC, with root mean squared errors (RMSE) of 61.2 W m^-2 at sites A, 106.1 W m^-2 at site B, and 137.2 W m^-2 at site C, largely due to storage fluxes across sites. Temporal integration improved VBR’s performance at tidal sites, where compensating water heat storage errors yielded low daily and monthly biases (site B: RMSE = 0.78 mm/d and 12 mm/month; r = 0.93; site C: RMSE = 0.90 mm/d and 13 mm/month; r = 0.93), with reduced major axis (RMA) regression slopes of 0.98 and ~0.91. In contrast, biomass heat storage at site A caused persistent biases (RMSEs = 0.97 mm/d and 23 mm/month; daily and monthly RMA slopes ~0.75; r = 0.85). These results highlight VBR’s limitations in environments with substantial storage fluxes. Despite this, VBR is cost-effective for estimating daily and monthly ET, with sensor costs at least tenfold lower than EC and simpler setup, making it suitable for ET monitoring in resource-limited and hard-to-access regions.

Long-term monotonic trends in water budget components in the contiguous United States: Insights from two hydrologic models

Goodling, Phillip; Foks, Sydney; Ayers, Jessica — Wed, 8 Apr 2026 13:47:35

Characterizing changes to water availability for domestic, industrial, agricultural, and other uses is essential to support water management. To better quantify these changes, the U.S. Geological Survey and National Science Foundation National Center for Atmospheric Research produced two hydrologic models simulating water budget components from 1980 to 2021 over the contiguous United States (CONUS). Both hydrologic models were driven by a common atmospheric forcing dataset and aggregated to common spatial and temporal scales, which enables a novel evaluation of congruency between the models. We present annual and seasonal trends in six water budget components (precipitation, evapotranspiration, streamflow, groundwater recharge, soil saturation, and snow water equivalent) based on the Mann–Kendall test for monotonic trend and Theil-Sen slope estimate for the water year 1983–2021 period for ~86,000 catchments in CONUS. Additional components and metrics from our analysis pipeline are available in an associated published dataset, which contains more than 46 million trend results. The water budget trends showed broad agreement with prior observational and modeling studies that indicate increasing trends in the northeast and decreasing trends in southwestern CONUS. We found the seasonal variability in water budget trends was greatest in the southern, central, and northwest CONUS. These findings support integrated trend assessments when coupled with trends in water quality and use.

Distribution and abundance of Least Bell’s Vireo (Vireo bellii pusillus) and Southwestern Willow Flycatcher (Empidonax traillii extimus) at the Hansen Dam Basin, Los Angeles County, California—2025 data summary

Lynn, Suellen; Kus, Barbara E. — Mon, 6 Apr 2026 13:36:04

Executive Summary

We surveyed for Least Bell’s Vireos (Vireo bellii pusillus; vireo) and Southwestern Willow Flycatchers (Empidonax traillii extimus; flycatcher) along Big Tujunga Creek in the Hansen Dam Basin in Los Angeles County, California, in 2025. Four vireo surveys were completed between April 17 and July 2, 2025, and three flycatcher surveys were completed between May 20 and July 2, 2025. We detected 62 territorial male vireos, 51 of which were confirmed as paired, and 2 transient vireos. Additionally, we detected 32 juvenile vireos during surveys. Seventy-seven percent of vireos were detected in habitat characterized as mixed willow, and 95 percent of vireos were detected in habitat with greater than 50-percent native plant cover. Most vireo territories were dominated by Goodding’s black willow (Salix gooddingii).

On May 20, 2025, we detected 18 transient Willow Flycatchers of unknown subspecies, none of which were confirmed to be paired, and no juveniles were detected. Mixed willow habitat was used by 78 percent of Willow Flycatchers, and all Willow Flycatchers were detected in habitat with greater than 50-percent native plant cover. Most Willow Flycatchers were detected in locations dominated by Goodding’s black willow.

Deep groundwater total dissolved solids mapping in the Dakota Group, Williston Basin, USA

Stephens, Michael; Hoogenboom, Bennett; Ball, Lyndsay; Chang, Will — Tue, 19 May 2026 15:33:56

Growing concern about the quantity of available freshwater around the world has led to interest in surveying groundwater total dissolved solids (TDS) below water well depths. Deep TDS has not been systematically mapped, and there is much to learn about the distribution and controls on deeper groundwater. In sedimentary basins across the United States, groundwater resources often overlie hydrocarbon resources, providing an opportunity to use borehole geophysical data collected for hydrocarbons to characterize groundwater and pore space resources. This study adapts a recently developed subsurface geostatistical and geophysical modeling approach to continuously map groundwater TDS, porosity, and temperature in the Dakota Group of the Williston Basin—an undercharacterized regional aquifer system overlying deeper hydrocarbon reservoirs. Groundwater TDS in the Dakota Group ranges from approximately 4800 to 26,900 mg/L. TDS patterns are stratified with higher TDS in the lower and upper Dakota Group, and relatively lower TDS in the middle Dakota Group. The lower TDS in the middle zone may represent a preferential regional flow path for lower-TDS meteoric recharge from the west. The alternating pattern of TDS may also be evidence of higher-TDS inflows into the Dakota Group from underlying and potentially from overlying aquifers. Porosity is lower near the center of the Williston Basin and tends to be higher to the east, which may be related to grain size distributions. The new regional TDS and porosity modeling serves as a quantitative reference for water users and provides supporting evidence for hypotheses on Dakota Group recharge.

Sampling and analysis plan for the water-quality monitoring program in Lake Koocanusa and upper Kootenai River, Montana, water years 2022–23

Fri, 3 Apr 2026 18:10:53

The U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency, collected water-quality samples and environmental data in Lake Koocanusa (also known as “Koocanusa Reservoir”), the Kootenai River, and the Tobacco River during water years 2022–23. The transboundary Lake Koocanusa is in southeastern British Columbia, Canada, and northwestern Montana, United States. It was formed by constructing Libby Dam on the Kootenai River 26 kilometers upstream from Libby, Montana. One of the lake sites and the Kootenai River site, in the Libby Dam tailwater (the outflow of the lake flow into the Kootenai River), were equipped with automated, high-frequency ServoSipper water samplers. At the lake site, these samplers were mounted to pontoon platforms during the summer, and a submersible ServoSipper sipper was deployed with ice buoys during the winter. Samples were automatically collected from multiple depths. At the Kootenai River site, these samplers were housed in the gage house. In water year 2022, discrete water-quality samples were collected every 4–6 weeks, year round, at all four lake sites in the Kootenai River between April and November. In water year 2023, discrete water-quality samples were collected at three lake sites and the Kootenai and Tobacco River sites every 4–6 weeks. The goal of this project was to collect multidepth, high-frequency vertical and temporal water-quality samples and data to understand the limnological and biological processes that control variations and trends in selenium concentrations and loads throughout Lake Koocanusa and in the Libby Dam tailwater at the southern end of the lake. This sampling and analysis plan documents the organization, sampling and data-collection scheme and design, pre- and post-collection processes, and quality-assurance and quality-control procedures of the Koocanusa/Kootenai water-quality monitoring program during water years 2022–23.

Seabed maps showing topography, ruggedness, backscatter intensity, sediment mobility, and the distribution of geologic substrates in quadrangle 3 of the Stellwagen Bank National Marine Sanctuary region offshore of Boston, Massachusetts

Valentine, Page; Cross, VeeAnn — Fri, 3 Apr 2026 17:31:25

The U.S. Geological Survey, in cooperation with the National Marine Sanctuary Program of the National Oceanic and Atmospheric Administration, has conducted seabed mapping and related research in the Stellwagen Bank National Marine Sanctuary (SBNMS) region since 1993. The area being mapped using geophysical and geological data includes the SBNMS and the surrounding region, which totals approximately 3,700 square kilometers (km²) and is subdivided into 18 quadrangles. The seabed is a glaciated terrain that is topographically and texturally diverse. Quadrangle 3, the subject of this scientific investigations map, has a mapped area of 185 km² and has water depths that range from about 30 meters (m) on the Stellwagen Bank crest to about 135 m in a basin east of South Ninety Bank, which lies off the eastern margin of Stellwagen Bank. Seven map types, each at a scale of 1:25,000, depict seabed topography, ruggedness, backscatter intensity, distribution of geologic substrates, sediment mobility, distribution of fine- and coarse-grained sand, and substrate mud content. These maps show the distribution of geologic substrates on the southeastern part of Stellwagen Bank, on adjacent banks and basins in deeper water to the east, in the eastern part of Race Point Channel to the south of the bank, and on the northern slope of Cape Cod. Interpretations of multibeam sonar bathymetric and seabed backscatter imagery, photographs, video imagery, and grain-size analyses were used to create the geology-based maps. Data from 309 stations were analyzed, including 279 sediment samples. The geologic substrate maps of quadrangle 3 show the distribution of 21 geologic substrates that represent a wide range of textures, such as rippled sand, immobile sand, immobile muddy sand, sand that partially veneers gravel, and boulder ridges. Mapped substrates are characterized by sediment grain-size composition, surface morphology, substrate layering, the mobility or immobility of substrate surfaces, and water depth range. This scientific investigations map portrays the major geological elements (substrates, topographic features, and processes) of environments in quadrangle 3. It is intended to provide a foundation for research into present and past sediment transport processes in a complex terrain, provide insights into the ecological requirements of invertebrate and vertebrate species that use the various substrates, and support seabed management in the region.

Methods for estimating daily upstream location of the freshwater-saltwater interface along the Maurice and Cohansey Rivers, New Jersey

Closson, Jennifer; Suro, Thomas; Niemoczynski, Lukasz — Thu, 16 Apr 2026 17:25:15

The Delaware River basin (DRB) provides drinking water to 15 million people in the surrounding area. Water is frequently withdrawn from the freshwater reaches of streams, above head of tide, in the DRB for use as public drinking water. During extended periods of low flow, saltwater can move upstream, which can threaten drinking-water supplies in the basin. Due to spatial patterns in bathymetry, tidal influences within the DRB, and varying weather conditions, it can be hard to predict the movement and upstream extent of the freshwater-saltwater interface, often defined as the salt-front. Although there is a relationship that predicts this location in the main stem of the Delaware River, there lacks a relationship for its tributaries, such as the Maurice and Cohansey Rivers in southwestern New Jersey. In this study, a relationship was developed between daily specific conductance (SC) at gage locations along the tidal river reaches of the Maurice and Cohansey Rivers to the daily upstream location of the salt-front. The study augmented existing real-time tide gage data with the collection of water temperature and specific conductance data to develop the relationship. Additionally, longitudinal profiles upstream of the selected tide gages were conducted during a range of high tide conditions to define the location of the salt-front. Equations were then developed that related the daily SC measured at the tide gage to the upstream location of the salt-front. The equations were used to estimate the daily upstream location of the salt-front for the period of July 15, 2021, to July 15, 2024. This work can aid in understanding the propagation of saltwater upstream, which can affect local communities and crop farmers along these tidal reaches of the DRB.

Thickness and other characteristics of overbank sediment deposited during an extreme flood in May 1978 along the Powder River, Montana

Moody, John; Meade, Robert — Fri, 3 Apr 2026 17:29:50

An extreme flood on the Powder River in southeastern Montana in May 1978 inundated its valley and deposited sediment on the floodplains and terraces at multiple heights. The recurrence interval for this flood was less than 1 percent in the reach between Moorhead and Broadus, Montana. Peak discharges at the U.S. Geological Survey streamgages at Moorhead and Broadus were 779 and 711 cubic meters per second (m³/s), respectively, the difference reflecting the water and sediment stored on the valley surfaces. Bankfull discharge depended on the height of the bank at the start of the valley transect and varied from 243 to 713 m³/s. Sediment-thickness and particle-size data were collected and analyzed in the autumn of 1978 by U.S. Geological Survey scientists at about 900 sites along 20 valley transects between Moorhead and Broadus, Mont. These transects were approximately orthogonal to the floodflow across the floodplain from near the edge of the channel to the high-water mark. Estimated maximum flood depths along these transects ranged from 0.9 to 4.2 meters.

Contrary to theory and controlled laboratory experiments, the distribution of sediment thickness and particle sizes along valley transects did not decrease systematically with distance from the main channel but were affected by the distribution of vegetation. Additionally, some water and sediment—primarily muds and silts—were conveyed by subsidiary channels (often connected to the main channel downriver from the valley transect) during the early stages of the flood before water overtopped the banks at the start of the valley transect. The vegetation created natural sediment traps in the recirculation and wake zones in the lee of trees and shrubs. Sediment that accumulated in these traps formed dunes and thus an undulating surface with many local maximums and minimums in sediment thicknesses. Sediment in the traps are referred to as lee dunes, which recorded flow conditions and a predominance of coarsening-upward sequence of particle sizes (mud to silt to sands) starting at the preflood surface. These sequences were associated with the rising limb of the hydrograph, and later as the flood began to recede, the lee dunes recorded a fining-upward sequence associated with the falling limb of the hydrograph.

Distribution and threats to the Arizona toad in Clark County

Stemp, Kenzi; Hossack, Blake — Wed, 15 Apr 2026 16:00:50

The Arizona Toad (Anaxyrus microscaphus) has been petitioned for protection under the U.S. Endangered Species Act (ESA) and is under evaluation for inclusion as a Covered Species under the Clark County Multi Species Habitat Conservation Plan Amendment (Clark County Department of Comprehensive Planning and USFWS, 2001; USFWS, 2015b). Although the species is locally abundant in some parts of its range, Arizona Toads have not been confirmed in Clark County since at least the 1980s. Extensive amphibian surveys in the late 1990s and examination of museum specimens reported no incidence of nonhybridized Arizona Toads in Clark County (Bradford et al., 2005), indicating a need to assess the current distribution of the species and availability of potential habitat for recovery efforts at the western extent of its historical range.

This document reports on the results of 2 years of Clark County Desert Conservation Program (DCP) funding to assess the status of Arizona Toad populations in the study area. This project supports data collection for an ongoing Species Status Assessment for the Arizona Toad by the U.S. Fish and Wildlife Service (USFWS) and provides data to support management actions to create or preserve suitable habitat for the Arizona Toad on Clark County Riparian Reserve properties.

Low streamflows in Massachusetts: Variability over space and time and relations with climatic and basin variables

Chamberlin, Catherine; Hodgkins, Glenn — Mon, 6 Apr 2026 14:19:40

Streamflows in Massachusetts have set record lows in recent years despite generally wetter conditions than during the drought of the 1960s, and the reasons for this are not known. To analyse potential drivers of low streamflows in Massachusetts, six low-flow metrics were computed at 107 streamgages. These metrics represent low-flow magnitude, magnitude normalized to median flows, and duration. Multiple linear regressions were used to analyse the variability of low flows over space and time. Potential explanatory variables were computed using climatic, land use, water use, and basin data. For all low-flow metrics, the ratio of precipitation to potential evapotranspiration (P/PET) in July–August explained the most variability, with decreasing P/PET largely explained by lower precipitation. Water/wetland area was a significant explanatory variable in all the normalized-magnitude and duration models, with greater area associated with lower normalized magnitudes and with shorter durations of low flows. Human influence (characterized by development, population, water use, and artificial water storage) had mixed effects. Trends from 1983 to 2022 in summer P/PET and human influence have been strongest in the eastern part of the state where the strongest decreases in flows are observed. Low flows in Massachusetts seem to be driven by a combination of low summer precipitation and human effects, though the specific mechanisms of human influence on flow likely vary between basins.

Fisheries research and monitoring activities of the Lake Erie Biological Station, 2025

Fri, 24 Apr 2026 13:54:40

Lake Erie has the most populated watershed of all the Great Lakes and has undergone dramatic anthropogenic changes. Since the 1800s, overexploitation of fish populations, habitat destruction, non-native species proliferation, industrial contamination, and changes in nutrient loading have impacted the fish community including declines in or extirpation of many native species (Regier et al. 1969, Hartman 1973; Leach & Nepszy 1976; Ludsin et al. 2001). Implementation of the Clean Water Act and Great Lakes Water Quality Agreement in the 1970s improved habitat conditions (Reutter 2019), which contributed to several strong percid year-classes (Vandergoot et al. 2019). These strong year-classes also benefited from more restrictive management practices that reduced harvest, ultimately rehabilitating Lake Erie percid stocks (Kayle et al. 2015, STC 2020). Historically, Lake Erie supported a cool water fish community dominated by percids and salmonids. Recently updated FCOs set forth a vision that “Lake Erie will consist of diverse fish communities that support ongoing societal benefits, including thriving commercial and recreational fisheries, improved fish habitat and desirable ecosystem performance, and reduced adverse impacts from invasive fish” (Francis et al. 2020). Today, mixed fisheries resulting from seasonally changing cool and warm water habitats have developed in Lake Erie, and the new FCOs reflect a desire to manage both predator and prey fish communities within them.

Although Lake Erie management agencies have traditionally focused on numerical indices of a few economically important species, aquatic ecosystem models are typically evaluated in terms of entire fish community biomass. As a result, our understanding of fish community structure and ecosystem dynamics from biomass-based models has been limited to short-term investigations and proxy measurements (e.g., length-weight conversion; FTG 2020). Therefore, many Lake Erie fish community databases are now incorporating biomass-based measurements.

In response, USGS revised the Lake Erie trawl program to provide biomass-based measurements for all encountered species (Table 1). The survey design change occurred in 2012, coincident with commissioning of a new research vessel and a change in bottom trawl gear. These modifications already altered the existing time series; therefore, the survey design was also expanded to include greater spatial coverage and increased sample size generating a new time series. The purpose of this report is to develop a comprehensive understanding of the long-term changes and fish community dynamics including population dynamics of key fishes of interest to management agencies, such as native percids and their prey. Here, we summarize survey results for the most recent series of West Basin trawl data from 2013 through 2025.

Note that a detailed description of the sampling process along with traditional numericallybased catch data (e.g., fish/ha) for individual species can be downloaded online (DuFour et al. 2026) or obtained for earlier years (https://doi.org/10.5066/F75M63X0; U.S. Geological Survey, Great Lakes Science Center 2019).

Extrinsic factors similarly affect nest survival of a threatened shorebird in natural and human-created habitats

Wed, 8 Apr 2026 15:02:20

Knowledge of factors that influence nest survival can inform effective conservation management for imperiled avian species. Habitat availability and quality are common priorities of conservation efforts, and climate and interspecific associations can also affect survival rates. In the lower Platte River system of eastern Nebraska, USA, Piping Plovers (Charadrius melodus, hereafter plovers) nest on river sandbars and different types of human-created off-river sites (i.e., sand and gravel mines, housing developments, and transitional sites) that are unique within the Northern Great Plains breeding population. However, off-river habitat may not be suitable for plover nesting long-term because of reduced habitat availability on both river sandbars and off-river sites. We evaluated plover nest survival at off-river and sandbar sites using data from 2008 to 2023 (n = 285). In addition, we examined the effects of extreme temperatures, proximity to conspecific and Interior Least Tern (Sternula antillarum athalassos, hereafter tern) nests, and temporal factors on nest survival. Plover nest survival did not differ between off-river sites and sandbars or amongst different off-river site types. Daily nest survival was 0.9818 (95% CI = 0.9729–0.9878) from 2008 to 2013 (unexclosed; n = 87) and 0.9950 (0.9918–0.9970) from 2014 to 2023 (exclosed; n = 198). Nest survival increased with the proportion of above average temperature days, increased with proximity to neighboring tern nests, and decreased with later nest initiation dates. Therefore, prioritizing early season nests, recognizing the benefits terns provide to plover nest survival, and monitoring the effect of climatic trends may aid future conservation efforts. Finally, as nest survival at off-river sites is comparable to sandbars, the predicted decline of habitat provided at off-river sites may reduce the overall breeding productivity and abundance of plovers in the lower Platte River system with ramifications to broader population viability.

Multiscale resource selection for a reintroduced elk population

Tue, 5 May 2026 13:42:44

Patterns of resource selection are driven by the decision-making processes of animals occurring at multiple scales from where to establish a home range (i.e., second order selection) to which resource patches to use within the home range (i.e., third order selection). Elk (Cervus canadensis) were reintroduced to southwestern Virginia, USA, from 2012 to 2014 following successful translocations onto reclaimed surface coal mines in the region. We sought to understand how elk have acclimated following their translocation using location data from GPS-collared adult female elk (n = 33) collected from 2019 to 2022 along with remotely sensed terrain and land cover data. We utilized continuous-time movement models paired with generalized linear mixed-effects modeling to describe seasonal resource selection at second and third orders. At both scales of selection and throughout the year, female elk selected reclaimed surface mines, conifer forests, ridgetops, and areas with lower terrain roughness, while avoiding mixed hardwood and oak (Quercus spp.) forests. Unmined open land was only selected at the third order during periods of forage scarcity (i.e., winter) and increased metabolic requirements (i.e., late gestation). Although surface coal mining leaves legacy environmental impacts on the landscape, management of these sites provides benefits to elk and maintains open habitat that is otherwise limited.

Challenges and opportunities for national-scale projections of future coastal landscape change

Mon, 18 May 2026 14:55:52

Local to global scale projections of future coastal landscape change are essential to improve land and resource management decisions that aim to prepare for and reduce risk exposure to impending coastal hazards. However, the availability of actionable knowledge is often limited due to the complexity of drivers of change, their consequences, and uncertainties that span disciplines. Varying spatial and temporal scales of physical processes that occur on landscapes with differing levels of coastal development and modification complicate decision-making. This paper reviews advances in observational approaches, modeling, and quantification of uncertainty to support the next generation of coastal change projections. We also highlight transdisciplinary opportunities for making this information more accessible and useful to support decision-making in the coastal zone. This includes consideration of different levels of user engagement to address specific needs, as well as resources and mechanisms that can support knowledge transfer broadly. A central recommendation is the development of a coordinated framework to support research integration including three components: accounting for relevant drivers, processes, and feedbacks; working with users to identify their information needs; and generating meaningful projections and uncertainty. Actions that would support such frameworks include the evaluation of models through coordinated research approaches and multi-model comparisons; ongoing investments in observations and data-driven approaches necessary to understand this complex environment; expansion of modeling capabilities using novel techniques; prioritization of co-development efforts between researchers and users; and continued training and support for transdisciplinary research.

Impact of fog on California waterfowl flight activity: Historical and modern insights into effects post-Clean Air Act

Fri, 1 May 2026 14:32:54

Since establishment of the Clean Air Act in the early 1970s, occurrence of the dense ‘Tule Fog’, historically prevalent throughout winter across California’s Central Valley, has substantially reduced. At the same time, waterfowl body masses have generally increased. Flight is metabolically expensive, and fog visually and navigationally impairs birds in flight, likely causing them to remain aloft for longer than usual periods. If less fog results in less flight and reduced energy expenditure, then fewer winter Tule fog events could contribute to increased body masses of California waterfowl since the mid-1980’s. Therefore, we aimed to assess the relationship between waterfowl flight and fog occurrence/density with historic (1991–93) and modern (2015–23) waterfowl tracking data in the Central Valley of California (CCV). Historic tracking data showed that the probability of flight increased with increasing fog density. Birds were significantly more likely to fly in fog than when there was no fog and most likely to fly in heavy fog. Modern data showed similar responses to fog with flight occurring significantly more during dawn fog events. This relationship between improved waterfowl body mass and fewer fog events may provide an opportunity to redirect scarce funding to focus on other population requirements such as improving habitats for nesting, molting and brood rearing that are currently lacking. Unanticipated benefits of the Clean Air Act should stand as a strong recommendation to maintaining this Act into the future.

Using large databases of groundwater chemistry in the northern Midwest USA: The effects of geologic and anthropogenic factors

Fri, 10 Apr 2026 18:49:38

Regional geochemical databases for the northern Midwest USA are being compiled to examine the various geogenic and anthropogenic factors that control the chemistry of groundwater. At the regional scale, variations are seen that are attributable to agricultural and urban effects, or to geologic factors. Examples of the former include enrichments of nitrate in groundwater, while examples of the latter mainly highlight geochemical differences between carbonate rocks and all other rock types in the region. This paper examines a few of these regional effects and the spatial scales at which they can be observed.