USGS Publications Warehouse

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

Fri, 5 Jun 2026 13:55:33

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

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

Miller, Samuel; Aguilar, Marcus; Helsley, Logan; Entrekin, Sally — Fri, 29 May 2026 14:23:57

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

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

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

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.

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.

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.

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

Wed, 13 May 2026 14:33:57

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

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

Pimont, Cedric; Thaler, Evan; Ebel, Brian A.; Bladon, Kevin — Thu, 14 May 2026 13:32:52

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

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

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.

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.

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.

Water scarcity and infrastructure risk of amplified seasonal sediment transport

Tue, 19 May 2026 14:28:32

Climate warming and deglaciation are reshaping hydrological seasonality in cold–dry regions, threatening the long-term sustainability of agriculture, ecosystems and local communities. However, existing evidence is limited to runoff seasonality. Changing sediment-transport seasonality, a more sensitive component, is emerging as a substantial yet under-recognized threat to water infrastructure. Leveraging monthly observations from the upper Tarim River from the 1960s to 2000s, we show that a warmer and wetter climate has intensified sediment-transport seasonality, with a 43% increase in summer sediment fluxes. Over half of this amplification stems from more frequent extreme sediment transport, particularly events triggered by high sediment supply rather than high discharge. Supported by a state-of-the-art river change dataset, we show that enhanced sediment seasonality and extreme sediment transport have largely contributed to increased river mobility since 2000. Sediment-driven changes are pushing riverine processes towards greater unpredictability and pose growing threats to water infrastructure and water security in vulnerable cold–dry regions.

An automated geographic information system-based hydraulic modeling tool for developing preliminary culvert designs for stream crossings in Massachusetts

Fri, 1 May 2026 19:08:48

Introduction

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

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

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.

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.

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.

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.

Spatially consistent but temporally divergent changes in nitrate and phosphorus loads and yields in Illinois watersheds, 1997–2022

Wed, 29 Apr 2026 15:01:54

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

Spatial heterogeneity of salt marsh vulnerability to sea-level rise: Dual controls of hydrological setting and salinity regime

Thu, 23 Apr 2026 14:59:21

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

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.

A novel drive-point multilevel system to investigate PFAS and other contaminants of global concern in the hyporheic zone of a wastewater effluent dominated stream

Meyer, J.; Mianecki, A.; Occhi, E.; Kolpin, Dana W.; LeFevre, G. — Wed, 29 Apr 2026 14:14:57

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

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.

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.

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

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

Kappel, William; Sherwood, Donald; Johnston, William — Wed, 3 Jun 2026 13:43:37

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

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

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

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

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

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

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