Cyanotoxin mixture models: Relating environmental variables and toxin co-occurrence to human exposure risk

Toxic cyanobacterial blooms, often containing multiple toxins, are a serious public health issue. However, there are no known models that predict a cyanotoxin mixture (anatoxin-a, microcystin, saxitoxin). This paper presents two cyanotoxin mixture models (MIX) and compares them to two microcystin (MC) models from data collected in 2016–2017 from three recurring cyanobacterial bloom locations in Kabetogama Lake, Voyageurs National Park (Minnesota, USA). Models include those using near-real-time environmental variables (readily available) and those using additional comprehensive variables (based on laboratory analyses). Comprehensive models (R² = 0.87 MC; R² = 0.86 MIX) explained more variability than the environmental models (R² = 0.58 MC; R² = 0.57 MIX). Although neither MIX model was a better fit than the MC models, the MIX models produced no false negatives in the calibration dataset, indicating that all observations above regulatory guidelines were simulated by the MIX models. This is the first known use of Virtual Beach software for a cyanotoxin mixture model, and the methods used in this paper may be applicable to other lakes or beaches.