Spatio-Temporal Modeling for Forecasting High-Risk Freshwater Cyanobacterial Harmful Algal Blooms in Florida

Forecasting freshwater cyanobacterial harmful algal blooms for Sentinel-3 satellite resolved U.S. lakes and reservoirs

This forecasting approach may be useful for water managers and associated public health managers to predict near-term future high-risk cyanobacterial harmful algal blooms (cyanoHAB) occurrence. Freshwater cyanoHABs may grow to excessive concentrations and cause human, animal, and environmental health concerns in lakes and reservoirs. Knowledge of the timing and location of cyanoHAB events is important for water quality management of recreational and drinking water systems. No quantitative tool exists to forecast cyanoHABs across broad geographic scales and at regular intervals. Publicly available satellite monitoring has proven effective in detecting cyanobacteria biomass near-real time within the United States. Weekly cyanobacteria abundance was quantified from the Ocean and Land Colour Instrument (OLCI) onboard the Sentinel-3 satellite as the response variable. An Integrated Nested Laplace Approximation (INLA) hierarchical Bayesian spatiotemporal model was applied to forecast World Health Organization (WHO) recreation Alert Level 1 exceedance >12 μg L-1 chlorophyll-a with cyanobacteria dominance for 2192 satellite resolved lakes in the United States across nine climate zones. The INLA model was compared against support vector classifier and random forest machine learning models; and Dense Neural Network, Long Short-Term Memory (LSTM), Recurrent Neural Network (RNN), and Gneural Network (GNU) neural network models. Predictors were limited to data sources relevant to cyanobacterial growth, readily available on a weekly basis, and at the national scale for operational forecasting. Relevant predictors included water surface temperature, precipitation, and lake geomorphology. Overall, the INLA model outperformed the machine learning and neural network models with prediction accuracy of 90% with 88% sensitivity, 91% specificity, and 49% precision as demonstrated by training the model with data from 2017 through 2020 and independently assessing predictions with data from the 2021 calendar year. The probability of true positive responses was greater than false positive responses and the probability of true negative responses was less than false negative responses. This indicated the model correctly assigned lower probabilities of events when they didn't exceed the WHO Alert Level 1 threshold and assigned higher probabilities when events did exceed the threshold. The INLA model was robust to missing data and unbalanced sampling between waterbodies.

One-Week-Ahead Prediction of Cyanobacterial Harmful Algal Blooms in Iowa Lakes

Cyanobacterial harmful algal blooms (CyanoHABs) pose serious risks to inland water resources. Despite advancements in our understanding of associated environmental factors and modeling efforts, predicting CyanoHABs remains challenging. Leveraging an integrated water quality data collection effort in Iowa lakes, this study aimed to identify factors associated with hazardous microcystin levels and develop one-week-ahead predictive classification models. Using water samples from 38 Iowa lakes collected between 2018 and 2021, feature selection was conducted considering both linear and nonlinear properties. Subsequently, we developed three model types (Neural Network, XGBoost, and Logistic Regression) with different sampling strategies using the nine selected variables (mcyA_M, TKN, % hay/pasture, pH, mcyA_M:16S, % developed, DOC, dewpoint temperature, and ortho-P). Evaluation metrics demonstrated the strong performance of the Neural Network with oversampling (ROC-AUC 0.940, accuracy 0.861, sensitivity 0.857, specificity 0.857, LR+ 5.993, and 1/LR- 5.993), as well as the XGBoost with downsampling (ROC-AUC 0.944, accuracy 0.831, sensitivity 0.928, specificity 0.833, LR+ 5.557, and 1/LR- 11.569). This study exhibited the intricacies of modeling with limited data and class imbalances, underscoring the importance of continuous monitoring and data collection to improve predictive accuracy. Also, the methodologies employed can serve as meaningful references for researchers tackling similar challenges in diverse environments.

Recent changes in cyanobacteria algal bloom magnitude in large lakes across the contiguous United States

Cyanobacterial blooms in inland lakes produce large quantities of biomass that impact drinking water systems, recreation, and tourism and may produce toxins that can adversely affect public health. This study analyzed nine years of satellite-derived bloom records and compared how the bloom magnitude has changed from 2008-2011 to 2016-2020 in 1881 of the largest lakes across the contiguous United States (CONUS). We determined bloom magnitude each year as the spatio-temporal mean cyanobacteria biomass from May to October and in concentrations of chlorophyll-a. We found that bloom magnitude decreased in 465 (25 %) lakes in the 2016-2020 period. Conversely, there was an increase in bloom magnitude in only 81 lakes (4 %). Bloom magnitude either didn't change, or the observed change was in the uncertainty range in the majority of the lakes (n = 1335, 71 %). Above-normal wetness and normal or below-normal maximum temperature over the warm season may have caused the decrease in bloom magnitude in the eastern part of the CONUS in recent years. On the other hand, a hotter and dryer warm season in the western CONUS may have created an environment for increased algal biomass. While more lakes saw a decrease in bloom magnitude, the pattern was not monotonic over the CONUS. The variations in temporal changes in bloom magnitude within and across climatic regions depend on the interactions between land use land cover (LULC) and physical factors such as temperature and precipitation. Despite expectations suggested by recent global studies, bloom magnitude has not increased in larger US lakes over this time period.

Satellites quantify the spatial extent of cyanobacterial blooms across the United States at multiple scales

Previous studies indicate that cyanobacterial harmful algal bloom (cyanoHAB) frequency, extent, and magnitude have increased globally over the past few decades. However, little quantitative capability is available to assess these metrics of cyanoHABs across broad geographic scales and at regular intervals. Here, the spatial extent was quantified from a cyanobacteria algorithm applied to two European Space Agency satellite platforms-the MEdium Resolution Imaging Spectrometer (MERIS) onboard Envisat and the Ocean and Land Colour Instrument (OLCI) onboard Sentinel-3. CyanoHAB spatial extent was defined for each geographic area as the percentage of valid satellite pixels that exhibited cyanobacteria above the detection limit of the satellite sensor. This study quantified cyanoHAB spatial extent for over 2,000 large lakes and reservoirs across the contiguous United States (CONUS) during two time periods: 2008-2011 via MERIS and 2017-2020 via OLCI when cloud-, ice-, and snow-free imagery was available. Approximately 56% of resolvable lakes were glaciated, 13% were headwater, isolated, or terminal lakes, and the rest were primarily drainage lakes. Results were summarized at national-, regional-, state-, and lake-scales, where regions were defined as nine climate regions which represent climatically consistent states. As measured by satellite, changes in national cyanoHAB extent did have a strong increase of 6.9% from 2017 to 2020 (|Kendall's tau (τ)| = 0.56; gamma (γ) = 2.87 years), but had negligible change (|τ| = 0.03) from 2008 to 2011. Two of the nine regions had moderate (0.3 ≤ |τ| < 0.5) increases in spatial extent from 2017 to 2020, and eight of nine regions had negligible (|τ| < 0.2) change from 2008 to 2011. Twelve states had a strong or moderate increase from 2017 to 2020 (|τ| ≥ 0.3), while only one state had a moderate increase and two states had a moderate decrease from 2008 to 2011. A decrease, or no change, in cyanoHAB spatial extent did not indicate a lack of issues related to cyanoHABs. Sensitivity results of randomly omitted daily CONUS scenes confirm that even with reduced data availability during a short four-year temporal assessment, the direction and strength of the changes in spatial extent remained consistent. We present the first set of national maps of lake cyanoHAB spatial extent across CONUS and demonstrate an approach for quantifying past and future changes at multiple spatial scales. Results presented here provide water quality managers information regarding current cyanoHAB spatial extent and quantify rates of change.

Development of a Risk Characterization Tool for Harmful Cyanobacteria Blooms on the Ohio River

A data-driven approach to characterizing the risk of cyanobacteria-based harmful algal blooms (cyanoHABs) was undertaken for the Ohio River. Twenty-five years of river discharge data were used to develop Bayesian regression models that are currently applicable to 20 sites spread-out along the entire 1579 km of the river’s length. Two site-level prediction models were developed based on the antecedent flow conditions of the two blooms that occurred on the river in 2015 and 2019: one predicts if the current year will have a bloom (the occurrence model), and another predicts bloom persistence (the persistence model). Predictors for both models were based on time-lagged average flow exceedances and a site’s characteristic residence time under low flow conditions. Model results are presented in terms of probabilities of occurrence or persistence with uncertainty. Although the occurrence of the 2019 bloom was well predicted with the modeling approach, the limited number of events constrained formal model validation. However, as a measure of performance, leave-one-out cross validation returned low misclassification rates, suggesting that future years with flow time series like the previous bloom years will be correctly predicted and characterized for persistence potential. The prediction probabilities are served in real time as a component of a risk characterization tool/web application. In addition to presenting the model’s results, the tool was designed with visualization options for studying water quality trends among eight river sites currently collecting data that could be associated with or indicative of bloom conditions. The tool is made accessible to river water quality professionals to support risk communication to stakeholders, as well as serving as a real-time water data monitoring utility.

Satellite remote sensing to assess cyanobacterial bloom frequency across the United States at multiple spatial scales

Cyanobacterial blooms can have negative effects on human health and local ecosystems. Field monitoring of cyanobacterial blooms can be costly, but satellite remote sensing has shown utility for more efficient spatial and temporal monitoring across the United States. Here, satellite imagery was used to assess the annual frequency of surface cyanobacterial blooms, defined for each satellite pixel as the percentage of images for that pixel throughout the year exhibiting detectable cyanobacteria. Cyanobacterial frequency was assessed across 2,196 large lakes in 46 states across the continental United States (CONUS) using imagery from the European Space Agency's Ocean and Land Colour Instrument for the years 2017 through 2019. In 2019, across all satellite pixels considered, annual bloom frequency had a median value of 4% and a maximum value of 100%, the latter indicating that for those satellite pixels, a cyanobacterial bloom was detected by the satellite sensor for every satellite image considered. In addition to annual pixel-scale cyanobacterial frequency, results were summarized at the lake- and state-scales by averaging annual pixel-scale results across each lake and state. For 2019, average annual lake-scale frequencies also had a maximum value of 100%, and Oregon and Ohio had the highest average annual state-scale frequencies at 65% and 52%. Pixel-scale frequency results can assist in identifying portions of a lake that are more prone to cyanobacterial blooms, while lake- and state-scale frequency results can assist in the prioritization of sampling resources and mitigation efforts. Satellite imagery is limited by the presence of snow and ice, as imagery collected in these conditions are quality flagged and discarded. Thus, annual bloom frequencies within nine climate regions were investigated to determine whether missing data biased results in climate regions more prone to snow and ice, given that their annual summaries would be weighted toward the summer months when cyanobacterial blooms tend to occur. Results were unbiased by the time period selected in most climate regions, but a large bias was observed for the Northwest Rockies and Plains climate region. Moderate biases were observed for the Ohio Valley and the Southeast climate regions. Finally, a clustering analysis was used to identify areas of high and low cyanobacterial frequency across CONUS based on average annual lake-scale cyanobacterial frequencies for 2019. Several clusters were identified that transcended state, watershed, and eco-regional boundaries. Combined with additional data, results from the clustering analysis may offer insight regarding large-scale drivers of cyanobacterial blooms.