USGS ScienceSearch

USGS · 70276589

Fifty years of riverine harmful algal bloom modeling: A global synthesis of approaches, challenges, and opportunities

Abstract

This systematic literature review critically examines 162 articles on harmful algal bloom (HAB) modeling in riverine systems to uncover persistent gaps, redefine critical challenges, and propose trackable opportunities to advance future modeling efforts. Articles largely focused on site-specific applications (93%) across more than 80 rivers worldwide. Most modeled systems were large, eutrophic rivers with flow modifications or obstructions. Geographic clustering of modeled systems was pronounced, with South Korea accounting for 26% of articles, followed by Europe (25%), United States (21%), and China (12%). Modeling approaches were led by process-based models (59%), though use of data-driven models (37%) increased over time, reflecting advances in computing and monitoring technology. Modeling endpoints varied widely across the articles with many focused on gross measures of algal abundance and fewer representing more refined endpoints like algal toxins or community composition. Furthermore, inconsistent units and taxonomic resolution hindered comparability between models. Datasets used for model development and calibration typically spanned 5 years, with weekly to monthly sampling at 1–10 sites, though durations and site counts were positively skewed. Quantitative metrics of model skill were often absent and included a diverse set of metrics when reported. Across all models, nutrients, light availability, streamflow, algal physiological processes, and water temperature emerged as key predictors, though algal processes were rarely incorporated in data-driven models. Scenario analyses primarily were conducted with process-based models and addressed flow management, whereas forecasting applications were less common and typically used data-driven models. After almost 50 years of riverine HAB modeling, persistent challenges include underrepresentation of benthic habitats, neglect of side-channel and backwater influences, insufficient documentation of river features, and weak linkages between modeled endpoints and potential harms. Addressing these gaps through reporting of contextual information, models from other aquatic settings, benchmark datasets, and community-driven tools could advance riverine HAB modeling towards increased transferability and ultimately operational forecasts.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Jennifer C. Murphy, Rebecca M. Gorney, Lisa V. Lucas, Jacob Aaron Zwart, Jennifer L. Graham. 2026. Fifty years of riverine harmful algal bloom modeling: A global synthesis of approaches, challenges, and opportunities. https://doi.org/10.1016/j.watres.2026.126240

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Aging of biosolids in fields prolongs the release of biosolid-borne PFAS and promotes the release of select PFCA in water

Land application of biosolids can introduce PFAS into croplands. Biosolids transported to surface waters through runoff may subsequently release PFAS into aquatic environments. The objective was to evaluate the effects of the aging of biosolids and the presence of sediment on the release of biosolid-borne PFAS to water at the water-sediment interface under anaerobic conditions. Through two independent rounds of microcosm experiments, we demonstrate that by day 136, most of the PFAS compounds (e.g., perfluorocarboxylic acids, or PFCAs) released <10 ng per g biosolids dry weight, except for 5:3 FTCA. Kinetically, the aging process substantially delayed the release of some long-chain PFAS (i.e., PFNA, PFDA, and PFOS) to water, due to the collapse of the porous structure of biosolids during aging. After day 136, the amount of long-chain PFAS released from aged biosolids reached levels similar to those released from the fresh biosolids due to the recovery of the pore structure during the release experiment. Mass balance analysis showed considerably higher, final masses of PFPeA, PFHxA, and PFOA than their initial masses, suggesting that precursors were biotransformed into these compounds during the release experiment, and biosolid aging led to elevated biotransformation and release of these three PFCAs. Sediment had little effect on the kinetics or extent of PFAS release from biosolids, except that the extent of FTCA release was reduced in the presence of sediment potentially due to adsorption. Overall, biosolid aging can delay the release of long-chain PFAS and promote the biotransformation to and release of PFCAs.

Water Research

Decoding the oxidative digestion mechanism for polystyrene nanoplastic detection in the Great Lakes using a customizable Raman spectral processing algorithm

Despite the concerns surging around nanoplastics (NPs) regarding their prevalence and bioavailability in freshwater systems, robust detection of NPs in complex environmental matrices is hindered by the lack of standardized sample pretreatment and a mechanistic understanding of oxidative digestion. Here, we systematically investigate the interaction between hydrogen peroxide (H 2 O 2 ) and polystyrene (PS) NPs during digestion in deionized (DI) water and four environmental matrices from in and around the Great Lakes Basin. To facilitate high-throughput analysis, we develop Pre_peak, a customizable Raman spectral processing algorithm that achieves >99% accuracy for both NP identification and interference rejection, allowing reliable NP quantification via pixel counting and systematic decoding of the oxidative digestion mechanisms. In DI water, varying H 2 O 2 doses from 0 to 30% has negligible effects on the recovery and Raman signal intensity of PS NPs over 24 hours of digestion. However, morphological changes and aggregation of PS NPs are observed when the H 2 O 2 dose exceeds 20%. Prolonged digestion further leads to progressive NP loss. In natural waters, the optimal dosage and digestion duration depend on matrix characteristics, including dissolved organic matter (DOM) and ion composition. This study provides mechanistic insights into NP–oxidant interactions and underscores the need for matrix-tailored digestion protocols to advance standardized NP detection in freshwater environments.

Great Lakes

Transit time modeling framework for predicting freshwater salinization in urban catchments

The salinity of inland freshwaters is rising globally, particularly in urban watersheds where winter road deicers are widely applied. Attributing stream salinity dynamics to specific sources and transport pathways remains challenging due to episodic salt inputs, engineered drainage, and strong coupling between hydrology and subsurface storage. We present a modeling framework that couples climate-driven deicer build-up and wash-off with transient transit time distribution theory to simulate salt transport through drainage, interflow, and groundwater pathways. Applied to an urban watershed in Northern Virginia (USA), the model reproduces ten years of high-frequency stream salinity measurements across daily-to-decadal timescales. The calibrated model implies an average deicer application of 206 tonnes Cl −1 yr −1 , or roughly one 20 kg bag of rock salt person −1 yr −1 when normalized by the 20,000 people living in the watershed. In winter months, higher infiltration routes a large fraction of snowmelt and deicers into shallow subsurface pathways, enhancing vadose-zone and interflow contributions to stream salinity. Limited subsurface storage capacity and seasonal hydrologic turnover flush excess chloride from the vadose zone and groundwater during subsequent summer storms. By linking climate-driven deicer inputs, hydrologic connectivity, and stream water age, the framework provides a transferable basis for diagnosing and managing freshwater salinization in urban watersheds.

Virginia