Next Article in Journal
Analysis of Siphonic Roof Drainage Systems with EPANET
Previous Article in Journal
Artificial Reefs around the World: A Review of the State of the Art and a Meta-Analysis of Its Effectiveness for the Restoration of Marine Ecosystems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling Cyanotoxin Production, Fate, and Transport in Surface Water Bodies Using CE-QUAL-W2

Department of Civil and Environmental Engineering, Portland State University, Portland, OR 97207-0751, USA
*
Author to whom correspondence should be addressed.
Environments 2023, 10(7), 122; https://doi.org/10.3390/environments10070122
Submission received: 28 April 2023 / Revised: 8 July 2023 / Accepted: 10 July 2023 / Published: 14 July 2023

Abstract

:
Cyanobacteria are frequently associated with forming toxic blooms. The toxins produced by cyanobacteria, cyanotoxins, are harmful to both humans and animals. Rising temperatures due to global climate change are expected to increase the occurrence of cyanobacteria, and it is vital that we protect our drinking water supplies and natural water resources. Modeling the production, fate, and transport of these toxins is an important step in limiting exposure to them and evaluating management strategies to mitigate their impact. The research provided here offers an overview of some of the main cyanotoxins of concern and presents preliminary models for the transport and fate of these toxins. Cyanotoxins can be either intracellular or extracellular, and a model for each was developed. The models were incorporated into the two-dimensional (longitudinal and vertical) hydrodynamic and water quality model CE-QUAL-W2. The toxin models were tested using a model of Henry Hagg Lake (Oregon, USA). The models were able to produce similar trends as found in published data, but since the toxin data available at Henry Hagg Lake was minimal, no direct comparisons between model results and field data were made. Four scenarios were conducted to test the functionality of the toxin models in CE-QUAL-W2. The predicted results from each test scenario matched the expected outcomes based on the parameters used in each scenario. Further applications of the toxin models to other water bodies with more consistent toxin data will help verify the accuracy of the models. This research provides a first step at modeling cyanotoxins using CE-QUAL-W2 and provides a framework to further develop the models through continued research of the cyanotoxins.

1. Introduction

Cyanobacteria are found throughout the world in both fresh and marine water, and since warmer temperatures have been shown to increase cyanobacterial growth, the occurrence of cyanobacteria could likely increase with the rise in global temperatures brought about by climate change [1,2,3,4]. While all bacterial and algal blooms can have harmful effects on the environment by reducing the availability of oxygen in water bodies, cyanobacterial blooms can also be harmful due to their ability to produce cyanotoxins, exposure to which may cause illnesses in humans [2,5,6]. The U.S. Environmental Protection Agency (EPA) Contaminant Candidate List 5 (CCL 5) lists cyanotoxins as contaminants in drinking water that may require future regulation [7]. By modeling the production and movement of these toxins, we can evaluate management strategies to mitigate their impact and help protect our drinking water supplies and natural resources.
Four of the main cyanotoxins studied (summarized in Table 1) are microcystins (MC), cylindrospermopsins (CYN), anatoxin-a (ATX), and saxitoxins (STX). Many publications have been produced in the last decade that provide reviews of the various cyanobacteria and cyanotoxins (see [5,8,9,10,11,12,13,14]).
One of the difficulties with predicting cyanotoxin concentrations is that not all toxic cyanobacteria may produce toxins in all environments, and some species include both toxic and nontoxic cyanobacteria strains [9,12,13]. Additionally, some strains may be capable of producing multiple different toxins. A study by Ballot et al. [15] detected both microcystin and anatoxin-a in one strain of Arthrospira fusiformis. The toxins can be present as intracellular toxins (toxins bound within a cell) or extracellular toxins (dissolved in the water). Cell lysis is one of the pathways by which cyanotoxins become extracellular, and some studies have also suggested that release from live cells may also occur. The active release is thought to occur for some CYN-producing species due to the high extracellular concentrations that have been observed [16,17], whereas cell lysis is likely the main cause of extracellular microcystins [18]. However, the overall quantity of toxins present in a water body appears to be largely a factor in the overall abundance of the particular cyanobacteria strains present producing the toxins [9].
Various models have been proposed for the production and degradation of cyanotoxins. For example, Harris and Graham [19] analyzed twelve different “linear and nonlinear regression models” for their ability to predict cyanobacteria and toxin concentrations in a reservoir in Kansas; however, the models were generally unable to predict peak concentrations. Mechanistic models have also been developed that relate cyanobacteria cellular concentrations with a constant toxin production coefficient or cell quota (mass of toxin per cell) to estimate the toxin concentration [20,21,22,23]. Table 2 shows a summary of some of these mechanistic models of cyanotoxin production and degradation.
The Hazen-Adams CyanoTOX (Cyanotoxin Tool for Oxidation Kinetics) tool has also been developed that predicts extracellular concentrations of various cyanotoxins after exposure to different oxidants [24,25]. The Version 3.0 of the CyanoTOX tool has been updated to include the release of intracellular toxins through lysing or leakage using assumed release rates (none, slow, moderate, fast, or instantaneous).
While the factors that affect cyanobacteria and cyanotoxin production are complex, the objectives of this research were to develop a general framework to model cyanotoxin production, distribution, and degradation in surface water bodies that could be incorporated into CE-QUAL-W2 [26,27], two-dimensional (longitudinal and vertical) hydraulic and water quality model of rivers, lakes/reservoirs and estuaries. Even though we are using the CE-QUAL-W2 framework to evaluate cyanotoxin production and transport, the framework is applicable to any water quality model. CE-QUAL-W2 has been used extensively to model the water quality of many different water bodies throughout the world. CE-QUAL-W2 is suited for long and narrow water bodies under stratified conditions and has been used in many different regions, such as those with temperate or tropical climates. A water quality and hydrodynamic model of Henry Hagg Lake in Oregon, USA, was used to test the cyanotoxin model framework.

2. Methods

A modeling framework was developed for cyanotoxin production and degradation so that a predictive model could answer questions about toxic cyanobacterial blooms and their impact on water quality. Once the algorithms were developed, they were added, tested, and evaluated in the model CE-QUAL-W2 [26,27]. The growth and death kinetics of cyanobacteria are treated the same as algae in CE-QUAL-W2 [27]. A simplified equation of the sources and sinks of cyanobacteria in a batch reactor (no sources or sinks associated with inflows and outflows, nor settling since the reactor is well-mixed) based on those used in CE-QUAL-W2 are shown in Equation (1) [27].
d a d t = k g k r k e k d a
where a is the concentration of cyanobacteria, k g is the growth rate, k r is respiration rate, k e is excretion rate, and k d is the death rate. Figure 1 shows some of the hypothesized sources and sinks of intracellular and extracellular cyanotoxins in a batch reactor based on the general growth and death kinetics of cyanobacteria as well as processes specific to cyanotoxins.
It is assumed that during cellular respiration, when the cells release carbon dioxide [28], the cells may use the toxin for cell processes, thereby decreasing the amount of intracellular toxin while not adding to the amount of extracellular toxin. It is also assumed that during excretion, when the cells release nutrients and organic carbon [28], the cells may also release toxins to the water column, thereby adding to the extracellular toxin concentration and decreasing the intracellular concentration. Internal decay, leakage, and active release are other hypothesized processes that might occur. Referring to Figure 1, it is assumed that for a specific toxin in a batch reactor, the sources and sinks of intracellular toxin concentration can be modeled by Equation (2), and the sources and sinks of extracellular toxin concentration can be modeled by Equation (3)
d C i n t r a d t = k g k r k e k d a β k l e a k C i n t r a k d e c a y _ i n t r a C i n t r a k a c t i v e C i n t r a
d C e x t r a d t = k e + k d a β + k l e a k C i n t r a k d e c a y _ e x t r a C e x t r a + k a c t i v e C i n t r a
where C i n t r a is intracellular toxin concentration, k g a β is the increase in intracellular toxin during growth, k r a β is the loss of intracellular toxin during respiration, ± k e a β is the loss of intracellular toxin (an increase in extracellular toxin) to cell excretion, ± k d a β is the loss of intracellular toxin (an increase in extracellular toxin) to death, β is the ratio of intracellular toxin mass to the mass of dry weight organic matter, ± k l e a k C i n t r a is the loss of intracellular toxin (an increase in extracellular toxin) due to leakage from the cell, k d e c a y _ i n t r a C i n t r a is the loss of intracellular toxin to internal decay, ± k a c t i v e C i n t r a is the loss of intracellular toxin (an increase in extracellular toxin) due to active toxin release from the cell, C e x t r a is the extracellular toxin concentration, and k d e c a y _ e x t r a C e x t r a is the loss of extracellular toxin to extracellular decay. It is assumed that the value of β represents the intracellular toxin production rate (100 percent of the production rate) and that the extracellular toxin is only released from the intracellular toxin produced.
Since the distinction between various processes nor the rates of these different processes are not well known or measured, many processes were combined and simplified when applying this model to CE-QUAL-W2 as shown in Equations (4) and (5) for intracellular and extracellular toxin concentrations, respectively.
C i n t r a = i = 1 n   a l g a l   g r o u p s C T P Φ a β
d C e x t r a d t = K a m Φ a β C T P + k r e l e a s e C i n t r a k d e c a y C e x t r a
where C i n t r a is the intracellular toxin concentration (mg toxin L−1), Φ a is the cyanobacteria concentration (mg dry weight organic matter L−1), β is the ratio of intracellular toxin mass to the mass of dry weight organic matter (mg toxin mg dry weight−1), C e x t r a is the extracellular toxin concentration (mg toxin L−1), K a m is the cyanobacteria mortality rate (day−1), k r e l e a s e is the combined toxin release rate of intracellular to extracellular representing any possible methods of excretion, leakage, and active release (day−1), k d e c a y is the extracellular toxin decay rate (day−1), and CTP (from CyanoToxin Produced) is described below. Separate toxin mass balances are used for each toxin of interest. In contrast to Equations (2) and (3), these simplified equations assume that the intracellular toxin concentration will not decrease as a result of toxin release (as determined by k r e l e a s e ), but rather the cell will manufacture any toxin lost to keep the ratio of intracellular toxin to dry-weight organic matter constant.
The CTP term mentioned above is the fraction of the algae group that produces the specific toxin. CTP represents which fraction of the modeled algae groups produce a toxin. For example, if a user selected one algae group to represent all the cyanobacteria in the reservoir, but only half of the group consisted of predicted microcystin producers, then CTP for Microcystin (CTP_MC) would be set to 0.5. If a group is not a toxin producer, then the CTP values for that group would be set to zero.
The intracellular toxin is calculated based on the concentration of cyanobacteria present in a model cell and the fraction of cyanobacteria producing that toxin. The intracellular concentration is summed for all the cyanobacteria groups that produce that toxin. There are no decay or other rates associated with the intracellular toxin concentration other than the scaling with the biomass algae concentration. The extracellular rate equation is a function of the mortality rate of cyanobacteria, the release rate of intracellular toxin, and the extracellular decay.
Table 3 shows the literature values for dry weight ratios (toxin mass per dry weight biomass), toxin cell quotas (mass of toxin per cell), ratios of toxin to chlorophyll a, and percent of toxin as extracellular for specific species producing each of the four toxins. See Buratti et al. [8], Chorus and Welker [9], Cirés and Ballot [10], and Testai et al. [14] for similar tables with additional values. Table 4 shows the degradation rates of cyanotoxins exposed to various conditions. The taxonomy of some of the cyanobacteria species in Table 3 has been revised [9].
Table 5 lists suggested ranges of values to use for each parameter from Equations (4) and (5) based on the literature values and model tests. CTP has been described previously, CTB (from CyanoToxin Biomass) represents the β term, CTR (from CyanoToxin Release) represents the k r e l e a s e term and CTD (from CyanoToxin Decay) represents the k d e c a y term.

3. Results and Discussion

The CE-QUAL-W2 model with the toxin code was applied to a previously developed and calibrated model of Henry Hagg Lake in Oregon, USA [53]. The model simulated eight years of continuous water quality from 1 January 2013 through 31 December 2020, and the model was calibrated to temperature, nutrient, chlorophyll a, and organic matter data from this time period. The first two years of the model simulation (2013 and 2014) were used in this paper to show annual trends. The climate of Oregon near Henry Hagg Lake shows distinct seasonal patterns with hot and dry (summer) weather typically from July through September and cold and wet (winter) weather typically from December through February.
There were minimal toxin data available at the lake, so the model was only used to test the proof of concept of the toxin equations within the code. Two toxins (MC and CYN) were included in the model simulations to test the model’s ability to predict multiple toxin concentrations and to test the toxin dynamics between MC and CYN. The model was set up using three algal groups where group 1 represented cooler temperature algal species (diatoms), group 2 represented warmer temperature algal species (greens), and group 3 represented cyanobacteria.
Four model scenarios were conducted to test the functionality of the code updates. The first model scenario tested the toxin decay in the reservoir without any toxin production by the algal or cyanobacterial groups within the lake. This was achieved by adding an initial toxin to the lake and setting all the CTP (fraction of group producing a toxin) values to zero for all the groups. The second model scenario tested toxin decay and production by cyanobacteria death only and for only one group (group 3). The third scenario tested toxin decay and production by cyanobacteria death and release for only group 3. The fourth scenario tested toxin decay and production with all three groups set as toxin producers to test the ability of the model to add together intracellular toxins for multiple groups. The third scenario was chosen as a representative model that would likely be used in further applications of the CE-QUAL-W2 model, and additional results are provided for this scenario. This scenario was chosen because it modeled all four of the toxin parameters for two toxins typically of concern in lakes and from only the cyanobacteria group. Table 6 summarizes the control file parameters chosen for each scenario based on Equations (4) and (5).
Total toxin concentrations, as predicted by CE-QUAL-W2, at the surface of the lake for each of the four scenarios are shown in Figure 2, Figure 3, Figure 4 and Figure 5.
The results from each scenario matched the predicted outcomes based on the toxin parameters chosen for each scenario. The results for the first scenario showed a faster decay for MC than CYN, which agrees with the higher decay rate for MC chosen than CYN. The results for the second scenario showed peaks of extracellular concentrations in the summer, with concentrations close to zero in the winter months. CYN was given a larger CTP value and a slower decay rate which resulted in higher CYN concentrations than MC concentrations. The results for the third scenario showed peaks at similar time periods to the second scenario, but the overall concentrations were higher. Since both toxins in the third scenario had release rates, the concentrations in this scenario were higher than in the second scenario. The fourth scenario showed an even greater increase in concentrations for both toxins compared to the other scenarios since all three groups were turned on as toxin producers. Figure 6 shows the concentration and percent of intracellular and extracellular MC, and Figure 7 shows the concentration and percent of intracellular and extracellular CYN for the third scenario.
The model results showed peak MC concentrations of about 0.4 ng mL−1 in the summer months and peak CYN concentrations of about 1 ng mL−1 in the summer months. The relative fractions of the intracellular and extracellular components of each toxin were similar to the literature values reported for microcystins and cylindrospermopsins. The intracellular fraction of MC is often found at a higher percentage than the extracellular fraction. Figure 6 shows that the model predicted approximately equal fractions of intracellular and extracellular MC during the summer months and predicted the intracellular toxin to be about 90 percent of the total toxin concentration during the winter months. In contrast to MC, CYN is often found at higher extracellular concentrations than intracellular concentrations. Figure 7 shows that the model predicted extracellular CYN to be about 70 percent of the total toxin concentration in the summer and fall. This agrees with CYN being found at high extracellular concentrations. The model predicted that most of the toxin would be present as the intracellular toxin in the winter for both MC and CYN. This is likely due to increased flows and more mixing in the winter that led to the dilution and removal of the extracellular toxin that exceeded the amount of toxin being produced during those months.
Toxin data were available at Henry Hagg Lake in the spring of 2019 for microcystin and cylindrospermopsin [54]. All the recorded values for CYN were classified as non-detects, and all but three measurements for MC were recorded as non-detects (0.15, 0.2, and 0.43 ng mL−1). The model results showed peak MC concentrations of about 0.4 ng mL−1 in the summer months and peak CYN concentrations of about 1 ng mL−1 in the summer months. The peak microcystin concentrations, as predicted by the model, were in the approximate range of those measured in the lake, but since there was minimal field data, it is difficult to adequately compare these values. Even though no cylindrospermopsin concentrations were detected in the lake, the toxin was still included in the model simulations to test the model’s ability to predict multiple toxin concentrations and to test how the dynamics vary between MC and CYN.

4. Conclusions

The goals of this research were to develop models for the transport and fate of cyanotoxins in surface water bodies and to incorporate the models into the two-dimensional (longitudinal and vertical) hydrodynamic and water quality model CE-QUAL-W2. The models were developed with various parameters that can be adjusted depending on the type of toxin being produced or the species producing the toxin, such as the leakage and loss rates, decay rates, and the ratio of intracellular toxin to dry-weight organic matter (β). These parameters allow for flexibility in applying the model to different cyanobacteria species and different toxins. While the model allows for flexibility in choosing values for each parameter, the values are fixed for each toxin. Future research could help to develop each of the model parameters as functions of changing environmental conditions such as light, temperature, and pH; however, as mentioned previously, the quantity of toxins present in a water body is due mainly to the abundance of toxin-producing strains.
The toxin models were incorporated into the CE-QUAL-W2 code and tested on a previously developed model of Henry Hagg Lake. Four scenarios were used to test the functionality of the toxin models in CE-QUAL-W2. The model results from each scenario matched the predicted outcomes based on the toxin parameters chosen for each scenario. There were limited field data available for the toxins, so the model was only used to test the proof of concept of the toxin equations within the code, and no direct comparisons between model results and field data were made. The model results of the relative fractions of intracellular and extracellular microcystin and cylindrospermopsin concentrations were similar to what would be expected based on the literature data for these two toxins. The model results showed that intracellular concentrations of both toxins were higher than the extracellular concentrations in the winter months. This is likely due to higher flows and more mixing in the winter leading to dilution of the extracellular concentration. Further applications of the model to water bodies with more consistent toxin data will help verify the accuracy of the CE-QUAL-W2 model in predicting toxin concentrations. The toxin models developed as part of this project are a first step in modeling toxin fate and transport in surface water bodies using CE-QUAL-W2.

Author Contributions

Conceptualization, B.G. and S.W.; data curation, B.G.; investigation, B.G.; methodology, B.G. and S.W.; project administration, S.W.; software, S.W.; supervision, S.W.; validation, B.G.; visualization, B.G.; writing—original draft preparation, B.G.; writing—review and editing, S.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

All relevant data are included within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, L.; Giesy, J.P.; Adamovsky, O.; Svirčev, Z.; Meriluoto, J.; Codd, G.A.; Mijovic, B.; Shi, T.; Tuo, X.; Li, S.-C.; et al. Challenges of using blooms of Microcystis spp. in animal feeds: A comprehensive review of nutritional, toxicological and microbial health evaluation. Sci. Total Environ. 2021, 764, 142319. [Google Scholar] [CrossRef] [PubMed]
  2. Fristachi, A.; Sinclair, J.L.; Hall, S.; Hambrook Berkman, J.A.; Boyer, G.; Burkholder, J.; Burns, J.; Carmichael, W.; DuFour, A.; Frazier, W.; et al. Occurrence of Cyanobacterial Harmful Algal Blooms Workgroup Report. In Cyanobacterial Harmful Algal Blooms: State of the Science and Research Needs; Hudnell, H.K., Ed.; Springer: New York, NY, USA, 2008; pp. 45–103. ISBN 978-0-387-75865-7. [Google Scholar]
  3. O’Neil, J.M.; Davis, T.W.; Burford, M.A.; Gobler, C.J. The rise of harmful cyanobacteria blooms: The potential roles of eutrophication and climate change. Harmful Algae 2012, 14, 313–334. [Google Scholar] [CrossRef]
  4. Paerl, H.W.; Gardner, W.S.; Havens, K.E.; Joyner, A.R.; McCarthy, M.J.; Newell, S.E.; Qin, B.; Scott, J.T. Mitigating cyanobacterial harmful algal blooms in aquatic ecosystems impacted by climate change and anthropogenic nutrients. Harmful Algae 2016, 54, 213–222. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Huang, I.-S.; Zimba, P.V. Cyanobacterial bioactive metabolites—A review of their chemistry and biology. Harmful Algae 2019, 83, 42–94. [Google Scholar] [CrossRef]
  6. United States Environmental Protection Agency (US EPA) Cyanobacteria and Cyanotoxins: Information for Drinking Water Systems, 2014, EPA-810F11001. Available online: https://www.epa.gov/sites/default/files/2014-08/documents/cyanobacteria_factsheet.pdf (accessed on 24 February 2020).
  7. United States Environmental Protection Agency (US EPA) Contaminant Candidate List 5—CCL 5. Available online: https://www.epa.gov/ccl/contaminant-candidate-list-5-ccl-5 (accessed on 5 December 2022).
  8. Buratti, F.M.; Manganelli, M.; Vichi, S.; Stefanelli, M.; Scardala, S.; Testai, E.; Funari, E. Cyanotoxins: Producing organisms, occurrence, toxicity, mechanism of action and human health toxicological risk evaluation. Arch. Toxicol. 2017, 91, 1049–1130. [Google Scholar] [CrossRef]
  9. Chorus, I.; Welker, M. Toxic Cyanobacteria in Water: A Guide to Their Public Health Consequences, Monitoring and Management, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2021; ISBN 978-1-00-308144-9. [Google Scholar]
  10. Cirés, S.; Ballot, A. A review of the phylogeny, ecology and toxin production of bloom-forming Aphanizomenon spp. and related species within the Nostocales (cyanobacteria). Harmful Algae 2016, 54, 21–43. [Google Scholar] [CrossRef]
  11. Corbel, S.; Mougin, C.; Bouaïcha, N. Cyanobacterial toxins: Modes of actions, fate in aquatic and soil ecosystems, phytotoxicity and bioaccumulation in agricultural crops. Chemosphere 2014, 96, 1–15. [Google Scholar] [CrossRef]
  12. Merel, S.; Walker, D.; Chicana, R.; Snyder, S.; Baurès, E.; Thomas, O. State of knowledge and concerns on cyanobacterial blooms and cyanotoxins. Environ. Int. 2013, 59, 303–327. [Google Scholar] [CrossRef]
  13. Pearson, L.A.; Dittmann, E.; Mazmouz, R.; Ongley, S.E.; D’Agostino, P.M.; Neilan, B.A. The genetics, biosynthesis and regulation of toxic specialized metabolites of cyanobacteria. Harmful Algae 2016, 54, 98–111. [Google Scholar] [CrossRef]
  14. Testai, E.; Scardala, S.; Vichi, S.; Buratti, F.M.; Funari, E. Risk to human health associated with the environmental occurrence of cyanobacterial neurotoxic alkaloids anatoxins and saxitoxins. Crit. Rev. Toxicol. 2016, 46, 385–419. [Google Scholar] [CrossRef]
  15. Ballot, A.; Krienitz, L.; Kotut, K.; Wiegand, C.; Pflugmacher, S. Cyanobacteria and cyanobacterial toxins in the alkaline crater lakes Sonachi and Simbi, Kenya. Harmful Algae 2005, 4, 139–150. [Google Scholar] [CrossRef]
  16. Preußel, K.; Chorus, I.; Fastner, J. Nitrogen Limitation Promotes Accumulation and Suppresses Release of Cylindrospermopsins in Cells of Aphanizomenon Sp. Toxins 2014, 6, 2932–2947. [Google Scholar] [CrossRef] [Green Version]
  17. Preußel, K.; Wessel, G.; Fastner, J.; Chorus, I. Response of cylindrospermopsin production and release in Aphanizomenon flos-aquae (Cyanobacteria) to varying light and temperature conditions. Harmful Algae 2009, 8, 645–650. [Google Scholar] [CrossRef]
  18. Rohrlack, T.; Hyenstrand, P. Fate of intracellular microcystins in the cyanobacterium Microcystis aeruginosa (Chroococcales, Cyanophyceae). Phycologia 2007, 46, 277–283. [Google Scholar] [CrossRef]
  19. Harris, T.D.; Graham, J.L. Predicting cyanobacterial abundance, microcystin, and geosmin in a eutrophic drinking-water reservoir using a 14-year dataset. Lake Reserv. Manag. 2017, 33, 32–48. [Google Scholar] [CrossRef]
  20. Grover, J.P.; Roelke, D.L.; Brooks, B.W. Modeling of plankton community dynamics characterized by algal toxicity and allelopathy: A focus on historical Prymnesium parvum blooms in a Texas reservoir. Ecol. Modell. 2012, 227, 147–161. [Google Scholar] [CrossRef]
  21. Grover, J.P.; Baker, J.W.; Roelke, D.L.; Brooks, B.W. Current Status of Mathematical Models for Population Dynamics of Prymnesium parvum in a Texas Reservoir. J. Am. Water Resour. Assoc. 2010, 46, 92–107. [Google Scholar] [CrossRef]
  22. Huo, X.; Chang, D.-W.; Tseng, J.-H.; Burch, M.D.; Lin, T.-F. Exposure of Microcystis aeruginosa to Hydrogen Peroxide under Light: Kinetic Modeling of Cell Rupture and Simultaneous Microcystin Degradation. Environ. Sci. Technol. 2015, 49, 5502–5510. [Google Scholar] [CrossRef]
  23. Jähnichen, S.; Long, B.M.; Petzoldt, T. Microcystin production by Microcystis aeruginosa: Direct regulation by multiple environmental factors. Harmful Algae 2011, 12, 95–104. [Google Scholar] [CrossRef]
  24. American Water Works Association AWWA Resources on Cyanobacteria/Cyanotoxins. Available online: https://www.awwa.org/Resources-Tools/Resource-Topics/Source-Water-Protection/Cyanobacteria-Cyanotoxins (accessed on 28 November 2022).
  25. Stanford, B.; Crafton, L.; Owen, C. A Digital Twin Model for Cyanotoxin Removal at Drinking Water Plants. Available online: https://www.hazenandsawyer.com/articles/a-digital-twin-model-for-cyanotoxin-removal-at-drinking-water-plants (accessed on 5 December 2022).
  26. Wells, S.A. CE-QUAL-W2 Hydrodynamic and Water Quality Model. Available online: https://cee.pdx.edu/w2/ (accessed on 20 December 2021).
  27. Wells, S.A. CE-QUAL-W2: A Two-Dimensional, Laterally Averaged, Hydrodynamic and Water Quality Model, Version 4.5; User Manual; Department of Civil and Environmental Engineering, Portland State University: Portland, OR, USA, 2021. [Google Scholar]
  28. Chapra, S.C. Surface Water-Quality Modeling; Waveland Press, Inc.: Long Grove, IL, USA, 2008; ISBN 978-1-57766-605-9. [Google Scholar]
  29. Halinen, K.; Jokela, J.; Fewer, D.P.; Wahlsten, M.; Sivonen, K. Direct Evidence for Production of Microcystins by Anabaena Strains from the Baltic Sea. Appl. Environ. Microbiol. 2007, 73, 6543–6550. [Google Scholar] [CrossRef] [Green Version]
  30. Long, B.M.; Jones, G.J.; Orr, P.T. Cellular Microcystin Content in N-Limited Microcystis aeruginosa Can Be Predicted from Growth Rate. Appl. Environ. Microbiol. 2001, 67, 278–283. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Jähnichen, S.; Ihle, T.; Petzoldt, T.; Benndorf, J. Impact of Inorganic Carbon Availability on Microcystin Production by Microcystis aeruginosa PCC 7806. Appl. Environ. Microbiol. 2007, 73, 6994–7002. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Orr, P.T.; Jones, G.J. Relationship between microcystin production and cell division rates in nitrogen-limited Microcystis aeruginosa cultures. Limnol. Oceanogr. 1998, 43, 1604–1614. [Google Scholar] [CrossRef]
  33. Oh, H.-M.; Lee, S.J.; Jang, M.-H.; Yoon, B.-D. Microcystin Production by Microcystis aeruginosa in a Phosphorus-Limited Chemostat. Appl. Environ. Microbiol. 2000, 66, 176–179. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Kosol, S.; Schmidt, J.; Kurmayer, R. Variation in peptide net production and growth among strains of the toxic cyanobacterium Planktothrix spp. Eur. J. Phycol. 2009, 44, 49–62. [Google Scholar] [CrossRef] [Green Version]
  35. Preußel, K.; Stüken, A.; Wiedner, C.; Chorus, I.; Fastner, J. First report on cylindrospermopsin producing Aphanizomenon flos-aquae (Cyanobacteria) isolated from two German lakes. Toxicon 2006, 47, 156–162. [Google Scholar] [CrossRef]
  36. Cirés, S.; Wörmer, L.; Timón, J.; Wiedner, C.; Quesada, A. Cylindrospermopsin production and release by the potentially invasive cyanobacterium Aphanizomenon ovalisporum under temperature and light gradients. Harmful Algae 2011, 10, 668–675. [Google Scholar] [CrossRef]
  37. Davis, T.W.; Orr, P.T.; Boyer, G.L.; Burford, M.A. Investigating the production and release of cylindrospermopsin and deoxy-cylindrospermopsin by Cylindrospermopsis raciborskii over a natural growth cycle. Harmful Algae 2014, 31, 18–25. [Google Scholar] [CrossRef]
  38. Harada, K.; Nagai, H.; Kimura, Y.; Suzuki, M.; Park, H.-D.; Watanabe, M.F.; Luukkainen, R.; Sivonen, K.; Carmichael, W.W. Liquid Chromatography/Mass Spectrometric Detection of Anatoxin-a, a Neurotoxin from Cyanobacteria. Tetrahedron 1993, 49, 9251–9260. [Google Scholar] [CrossRef]
  39. Ballot, A.; Fastner, J.; Lentz, M.; Wiedner, C. First report of anatoxin-a-producing cyanobacterium Aphanizomenon issatschenkoi in northeastern Germany. Toxicon 2010, 56, 964–971. [Google Scholar] [CrossRef]
  40. Gagnon, A.; Pick, F.R. Effect of nitrogen on cellular production and release of the neurotoxin anatoxin-a in a nitrogen-fixing cyanobacterium. Front. Microbiol. 2012, 3, 1–7. [Google Scholar] [CrossRef] [Green Version]
  41. Aráoz, R.; Nghiêm, H.-O.; Rippka, R.; Palibroda, N.; de Marsac, N.T.; Herdman, M. Neurotoxins in axenic oscillatorian cyanobacteria: Coexistence of anatoxin-a and homoanatoxin-a determined by ligand-binding assay and GC/MS. Microbiology 2005, 151, 1263–1273. [Google Scholar] [CrossRef]
  42. Negri, A.P.; Jones, G.J. Bioaccumulation of paralytic shellfish poisoning (PSP) toxins from the cyanobacterium Anabaena circinalis by the freshwater mussel Alathyria condola. Toxicon 1995, 33, 667–678. [Google Scholar] [CrossRef] [PubMed]
  43. Velzeboer, R.M.A.; Baker, P.D.; Rositano, J.; Heresztyn, T.; Codd, G.A.; Raggett, S.L. Geographical patterns of occurrence and composition of saxitoxins in the cyanobacterial genus Anabaena (Nostocales, Cyanophyta) in Australia. Phycologia 2000, 39, 395–407. [Google Scholar] [CrossRef]
  44. Ongley, S.E.; Pengelly, J.J.L.; Neilan, B.A. Elevated Na+ and pH influence the production and transport of saxitoxin in the cyanobacteria Anabaena circinalis AWQC131C and Cylindrospermopsis raciborskii T3. Environ. Microbiol. 2016, 18, 427–438. [Google Scholar] [CrossRef]
  45. Tsuji, K.; Naito, S.; Kondo, F.; Ishikawa, N.; Watanabe, M.F.; Suzuki, M.; Harada, K. Stability of Microcystins from Cyanobacteria: Effect of Light on Decomposition and Isomerization. Environ. Sci. Technol. 1994, 28, 173–177. [Google Scholar] [CrossRef]
  46. León, C.; Boix, C.; Beltrán, E.; Peñuela, G.; López, F.; Sancho, J.V.; Hernández, F. Study of cyanotoxin degradation and evaluation of their transformation products in surface waters by LC-QTOF MS. Chemosphere 2019, 229, 538–548. [Google Scholar] [CrossRef]
  47. Edwards, C.; Graham, D.; Fowler, N.; Lawton, L.A. Biodegradation of microcystins and nodularin in freshwaters. Chemosphere 2008, 73, 1315–1321. [Google Scholar] [CrossRef]
  48. Chiswell, R.K.; Shaw, G.R.; Eaglesham, G.; Smith, M.J.; Norris, R.L.; Seawright, A.A.; Moore, M.R. Stability of Cylindrospermopsin, the Toxin from the Cyanobacterium, Cylindrospermopsis raciborskii: Effect of pH, Temperature, and Sunlight on Decomposition. Environ. Toxicol. 1999, 14, 155–161. [Google Scholar] [CrossRef]
  49. Smith, C.; Sutton, A. The Persistence of Anatoxin-a in Reservoir Water; Report No FR0427; Foundation for Water Research: Marlow, UK, 1993. [Google Scholar]
  50. Stevens, D.K.; Krieger, R.I. Stability studies on the cyanobacterial nicotinic alkaloid anatoxin-a. Toxicon 1991, 29, 167–179. [Google Scholar] [CrossRef]
  51. Pereira, P.; Dias, E.; Franca, S. Persistence of Paralytic Shellfish Toxins in Freshwater Environments. In Proceedings of the 10th International Conference on Harmful Algae, St. Pete Beach, FL, USA, 21–25 October 2002; Steidinger, K.A., Landsberg, J.H., Tomas, C.R., Vargo, G.A., Eds.; Florida Fish and Wildlife Conservation Commission, Florida Institute of Oceanography, and Intergovernmental Oceanographic Commission of UNESCO: St. Petersburg, FL, USA, 2004; pp. 166–168. [Google Scholar]
  52. Jones, G.J.; Negri, A.P. Persistence and degradation of cyanobacterial paralytic shellfish poisons (PSPs) in freshwaters. Water Res. 1997, 31, 525–533. [Google Scholar] [CrossRef]
  53. Wells, S.A.; Garstecki, B.; Berger, C.J.; Mansell, S. Tualatin River and Hagg Lake 2013–2020 Water Quality Model Development and Calibration. Department of Civil and Environmental Engineering. Portland State University: Portland, OR, USA, 2022; Unpublished work. [Google Scholar]
  54. Mansell, S.; Clean Water Services, Hillsboro, OR, USA. Personal communication, 2021.
Figure 1. Hypothesized sources and sinks of intracellular and extracellular cyanotoxins.
Figure 1. Hypothesized sources and sinks of intracellular and extracellular cyanotoxins.
Environments 10 00122 g001
Figure 2. CE-QUAL-W2 toxin and algae predictions for scenario 1 at the surface of Henry Hagg Lake. Microcystins: MC, Cylindrospermopsins: CYN, cooler temperature algal species: Group1, warmer temperature algal species: Group2, cyanobacteria: Group3.
Figure 2. CE-QUAL-W2 toxin and algae predictions for scenario 1 at the surface of Henry Hagg Lake. Microcystins: MC, Cylindrospermopsins: CYN, cooler temperature algal species: Group1, warmer temperature algal species: Group2, cyanobacteria: Group3.
Environments 10 00122 g002
Figure 3. CE-QUAL-W2 toxin and algae predictions for scenario 2 at the surface of Henry Hagg Lake. Microcystins: MC, Cylindrospermopsins: CYN, cooler temperature algal species: Group1, warmer temperature algal species: Group2, cyanobacteria: Group3.
Figure 3. CE-QUAL-W2 toxin and algae predictions for scenario 2 at the surface of Henry Hagg Lake. Microcystins: MC, Cylindrospermopsins: CYN, cooler temperature algal species: Group1, warmer temperature algal species: Group2, cyanobacteria: Group3.
Environments 10 00122 g003
Figure 4. CE-QUAL-W2 toxin and algae predictions for scenario 3 at the surface of Henry Hagg Lake. Microcystins: MC, Cylindrospermopsins: CYN, cooler temperature algal species: Group1, warmer temperature algal species: Group2, cyanobacteria: Group3.
Figure 4. CE-QUAL-W2 toxin and algae predictions for scenario 3 at the surface of Henry Hagg Lake. Microcystins: MC, Cylindrospermopsins: CYN, cooler temperature algal species: Group1, warmer temperature algal species: Group2, cyanobacteria: Group3.
Environments 10 00122 g004
Figure 5. CE-QUAL-W2 toxin and algae predictions for scenario 4 at the surface of Henry Hagg Lake. Microcystins: MC, Cylindrospermopsins: CYN, cooler temperature algal species: Group1, warmer temperature algal species: Group2, cyanobacteria: Group3.
Figure 5. CE-QUAL-W2 toxin and algae predictions for scenario 4 at the surface of Henry Hagg Lake. Microcystins: MC, Cylindrospermopsins: CYN, cooler temperature algal species: Group1, warmer temperature algal species: Group2, cyanobacteria: Group3.
Environments 10 00122 g005
Figure 6. Scenario 3 toxin predictions for microcystin: (top) intracellular and extracellular percent of the total toxin concentration, (bottom) intracellular and extracellular toxin concentrations.
Figure 6. Scenario 3 toxin predictions for microcystin: (top) intracellular and extracellular percent of the total toxin concentration, (bottom) intracellular and extracellular toxin concentrations.
Environments 10 00122 g006
Figure 7. Scenario 3 toxin results for cylindrospermopsin: (top) intracellular and extracellular percent of the total toxin concentration, (bottom) intracellular and extracellular toxin concentrations.
Figure 7. Scenario 3 toxin results for cylindrospermopsin: (top) intracellular and extracellular percent of the total toxin concentration, (bottom) intracellular and extracellular toxin concentrations.
Environments 10 00122 g007
Table 1. Summary of the cyanotoxins microcystins, cylindrospermopsins, anatoxin-a, and saxitoxins [9,11].
Table 1. Summary of the cyanotoxins microcystins, cylindrospermopsins, anatoxin-a, and saxitoxins [9,11].
ToxinMicrocystinsCylindrospermopsinsAnatoxin-aSaxitoxins
Toxin TypeHepatotoxinCytotoxinNeurotoxinNeurotoxin
OccurrenceMostly found intracellularly. Extracellular portion likely due mainly to cell lysis.High extracellular fractions observed. Extracellular portion likely due to active release from cells in addition to cell lysis.Extracellular portion likely due mainly to cell lysis.STX congeners may undergo reactions that can change their toxicity.
Table 2. Cyanotoxin models in the literature.
Table 2. Cyanotoxin models in the literature.
EquationParametersReferences
d C d t = μ c N c ε c k c C C = c y a n o t o x i n   c o n c e n t r a t i o n , μ g L 1
μ c = c y a n o b a c t e r i a   g r o w t h   r a t e ,   d a y 1
N c = c y a n o b a c t e r i a   c o n c e n t r a t i o n ,   c e l l   L 1
ε c = c y a n o t o x i n   p r o d u c t i o n   c o e f f i c i e n t ,   μ g   c e l l 1
k c = f i r s t   o r d e r   c y a n o t o x i n   d e c a y   r a t e ,   d a y 1
[20,21]
d M d t = p d X d t d M M M = M C   c o n c e n t r a t i o n , f g   m L 1
X = c y a n o b a c t e r i a   c o n c e n t r a t i o n ,   c e l l   m L 1
p = M C   p r o d u c t i o n   c o e f f i c i e n t , f g   c e l l 1
d M = f i r s t   o r d e r   i n t r a c e l l u l a r   M C   d e p l e t i o n   r a t e , d a y 1
[23]
d C e d t = q d N d t k O H , M C C e C O H
d C i d t = q d N d t
d C t d t = d C e d t + d C i d t
C e = e x t r a c e l l u l a r   M C   c o n c e n t r a t i o n ,   μ g   L 1
C i = c e l l b o u n d   M C   c o n c e n t r a t i o n , μ g   L 1
C t = t o t a l   M C   c o n c e n t r a t i o n ,   μ g   L 1
q = M C   c e l l   q u o t a ,   μ g   c e l l 1
k O H , M C = s e c o n d   o r d e r   r e a c t i o n   r a t e   c o n s t a n t   b e t w e e n   O H   r a d i c a l   a n d   M C s , M 1 s 1
C O H = O H   c o n c e n t r a t i o n , M
N = c y a n o b a c t e r i a   c o n c e n t r a t i o n   o f   i n t a c t   c e l l s ,   c e l l   L 1
[22]
Table 3. Literature values for species producing cyanotoxins.
Table 3. Literature values for species producing cyanotoxins.
SpeciesRatio of Toxin to Dry Weight
(μg Toxin g−1 DW)
Cell Quota
(fg Toxin Cell−1)
Ratio of Toxin to Chlorophyll a
(µg Toxin µg−1 chl a)
Percent Extracellular of Total Toxin (Intracellular + Extracellular)References
Microcystins
Anabaena1300–3900 [29]
Microcystis aeruginosa1180–6470 0.59 [30]
Microcystis aeruginosa 18–240.70–0.81 a21–47 [31]
Microcystis aeruginosa1500–9500 b56–165 [32]
Microcystis aeruginosa555–1113 [33]
Planktothrix agardhii1170–446044–343 0–62[34]
Planktothrix rubescens320–451027–857 0–23[34]
Cylindrospermopsins
Aphanizomenon flos-aquae2300–6600 [35]
Aphanizomenon flos-aquae 8–58[17]
Aphanizomenon ovalisporum90–63702–1910.01–0.5323–64[36]
Cylindrospermopsis raciborskii 10–25 b 14–50[37]
Anatoxin-a
Anabaena circinalis1396 [38]
Anabaena flos-aquae1107–13,013 [38]
Aphanizomenon sp.1562 [38]
Aphanizomenon issatschenkoi 100 21–39[39]
Aphanizomenon issatschenkoi6–1683 3–47[40]
Oscillatoria sp.4000 [41]
Oscillatoria sp.2713 [38]
Saxitoxins
Anabaena circinalis1580 c [42]
Anabaena circinalis7–2553 c<1–1105 c [43]
Anabaena circinalis 36–38 d[44]
Cylindrospermopsis raciborskii 53–63 d[44]
DW, dry weight. a Value obtained by dividing toxin cell quota by chlorophyll a cell quota. b Value estimated from figure in the cited reference. c Total saxitoxins (paralytic shellfish poisoning toxins). d Saxitoxin congener only, value estimated from ratio of intracellular to extracellular toxins.
Table 4. Influences of different parameters on cyanotoxin degradation.
Table 4. Influences of different parameters on cyanotoxin degradation.
ParameterToxin Degradation Rates aReferences
Microcystins
Sunlight and range of water-extractable pigment concentrations from 0–5 mg mL−10.006–0.159 day−1 (1–85% left after 29 days as estimated from figure in the cited reference)[45]
Artificial equivalence to approximately 20 days of natural sunlight in surface water0.071 day−1 (24% left after 144 h of light, equivalent to 20 days of natural sunlight)[46]
254 nm UV light at 250 W m−2 in surface water1326 day−1 (1% left after 5 min)[46]
Chlorination by 3.0 mg L−1 NaOCl in surface water798 day−1 (33% left after 2 min)[46]
Biodegradation in water samples0.039–0.173 day−1 (half-lives of 4–18 days)[47]
Cylindrospermopsins
Natural sunlight4.2–11.1 day−1 (half-lives of 1.5 and 4 h for algal extract solutions)
0.046–0.063 day−1 (half-lives of 11 and 15 days for natural water samples)
[48]
Range of artificial light from 9–42 μE m−2 s−1 with initial toxin concentrations of 1 and 4 mg L−10.005–0.025 day−1 (42–84% left after 5 weeks)[48]
Range of pH values from 4–10 with initial toxin concentrations of 1 and 4 mg L−10.004–0.005 day−1 (75–81% left after 8 weeks)[48]
Range of temperatures from 4–35 °C with initial toxin concentrations of 1 and 4 mg L−10.008–0.019 day−1 (77–89% left after 14 days)[48]
257 nm UV light at 300 mW m−210.3 day−1 (concentration decreased from 1.5 to 1.3 mg L−1 after 20 min)[48]
257 nm UV light at 400 mW m−20.924 day−1 (half-life of 18 h)[48]
Artificial equivalence to approximately 20 days of natural sunlight in surface water0.032 day−1 (53% left after 144 h of light, equivalent to 20 days of natural sunlight)[46]
254 nm UV light at 250 W m−2 in surface water27 day−1 (1% left after 250 min)[46]
Chlorination by 3.0 mg L−1 NaOCl in surface water2318 day−1 (4% left after 2 min)[46]
Anatoxin-a
Normal pH conditions and microbes0.139 day−1 (half-life of 5 days) [49]
Natural sunlight at pH range of 6–12 in aqueous toxin solution3.0–10.4 day−1 (half-lives ranging between 96 and 330 min)[50]
Natural sunlight at pH of 9 in algal lysate solution1.4 day−1 (half-life of 690 min)[50]
Headspace purged with nitrogen at pH of 90.069 day−1 (half-life of 10 days)[50]
Headspace purged with oxygen at pH of 90.139 day−1 (half-life of 5 days)[50]
Saxitoxins
Temperature of 20 °C and pH values of 7 or 9 in sterile water0.011 day−1 (total saxitoxins)[51]
Temperature of 30 °C and pH value of 7 in sterile water0.022 day−1 (total saxitoxins)[51]
Temperature of 20 °C in culture medium under different sterilization and deproteinization conditions0.089–0.240 day−1 (total saxitoxins)[51]
Sterile water at 25 °C0.015–0.033 day−1 (individual congeners)[52]
Irrigation drain water at 25 °C0.025–0.075 day−1 (individual congeners)[52]
a Unless provided directly, the decay rate was estimated using the first-order decay equation c = c0e−kt, where c is the toxin concentration after time (t), c0 is the initial toxin concentration, and k is the decay rate [28].
Table 5. Suggested model parameter ranges for CE-QUAL-W2 toxin models.
Table 5. Suggested model parameter ranges for CE-QUAL-W2 toxin models.
ParameterMCCYNATXSTX
CTP, fraction of cyanobacteria concentration producing toxinWater body dependent, determined by species present
CTB, ratio of intracellular toxin to dry weight biomass (mg-toxin mg-DW−1)0.0005–0.010.0001–0.010.0001–0.0150.0001–0.005
CTR, release rate day−1Approx. 0–1 times excretion rateApprox. 0–2 times excretion rate Approx. 0–1 times excretion rateApprox. 0–1 times excretion rate
CTD, extracellular decay day−10.01–0.20.001–0.10.05–1.00.01–0.1
DW, dry weight.
Table 6. Summary of toxin model test scenarios for Henry Hagg Lake.
Table 6. Summary of toxin model test scenarios for Henry Hagg Lake.
ParameterScenario 1Scenario 2Scenario 3Scenario 4
MCCYNMCCYNMCCYNMCCYN
CTP, fraction of cyanobacteria producing toxin, fraction000.10.30.10.30.10.3
CTB, ratio of intracellular toxin to dry weight biomass, fraction000.010.0050.010.0050.010.005
CTR, release rate, day−100000.010.030.010.03
CTD, extracellular decay rate, day−10.10.050.10.050.10.050.10.05
Initial concentration10 ng mL−10 ng mL−10 ng mL−10 ng mL−1
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Garstecki, B.; Wells, S. Modeling Cyanotoxin Production, Fate, and Transport in Surface Water Bodies Using CE-QUAL-W2. Environments 2023, 10, 122. https://doi.org/10.3390/environments10070122

AMA Style

Garstecki B, Wells S. Modeling Cyanotoxin Production, Fate, and Transport in Surface Water Bodies Using CE-QUAL-W2. Environments. 2023; 10(7):122. https://doi.org/10.3390/environments10070122

Chicago/Turabian Style

Garstecki, Bernadel, and Scott Wells. 2023. "Modeling Cyanotoxin Production, Fate, and Transport in Surface Water Bodies Using CE-QUAL-W2" Environments 10, no. 7: 122. https://doi.org/10.3390/environments10070122

APA Style

Garstecki, B., & Wells, S. (2023). Modeling Cyanotoxin Production, Fate, and Transport in Surface Water Bodies Using CE-QUAL-W2. Environments, 10(7), 122. https://doi.org/10.3390/environments10070122

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop