Environmental drivers behind the exceptional increase in cyanobacterial blooms in Okavango Delta, Botswana

The Okavango Delta region in Botswana experienced exceptionally intense landscape-wide cyanobacterial harmful algal blooms (CyanoHABs) in 2020. In this study, the drivers behind CyanoHABs were determined from thirteen independent environmental variables, including vegetation indices, climate and meteorological parameters, and landscape variables. Annual Land Use Land Cover (LULC) maps were created from 2017 to 2020, with ∼89% accuracy to compute landscape variables such as LULC change. Generalized Additive Models (GAM) and Structural Equation Models (SEM) were used to determine the most important drivers behind the CyanoHABs. Normalized Difference Chlorophyll Index (NDCI) and Green Line Height (GLH) algorithms served as proxies for chlorophyll-a (green algae) and phycocyanin (cyanobacteria) concentrations. GAM models showed that seven out of the thirteen variables explained 89.9% of the variance for GLH. The models showcased that climate variables, including monthly precipitation (8.8%) and Palmer Severity Drought Index- PDSI (3.2%), along with landscape variables such as changes in Wetlands area (7.5%), and Normalized Difference Vegetation Index (NDVI) (5.4%) were the determining drivers behind the increased cyanobacterial activity within the Delta. Both PDSI and NDVI showed negative correlations with GLH, indicating that increased drought conditions could have led to large increases in toxic CyanoHAB activity within the region. This study provides new information about environmental drivers which can help monitor and predict regions at risk of future severe CyanoHABs outbreaks in the Okavango Delta, Botswana, and other similar data-scarce and ecologically sensitive areas in Africa. Plain Language Summary: The waters of the Okavango Delta in Northern Botswana experienced an exceptional increase in toxic cyanobacterial activity in recent years. Cyanobacterial blooms have been shown to affect local communities and wildlife in the past. To determine the drivers behind this increased bloom activity, we analyzed the effects of thirteen independent environmental variables using two different statistical models. Within this research, we focused on vegetation indices, meteorological, and landscape variables, as previous studies have shown their effect on cyanobacterial activity in other parts of the world. While driver determination for cyanobacteria has been done before, the environmental conditions most important for cyanobacterial growth can be specific to the geographic setting of a study site. The statistical analysis indicated that the increases in cyanobacterial bloom activity within the region were mainly driven by persistent drier conditions. To our knowledge, this is the first study to determine the driving factors behind cyanobacterial activity in this region of the world. Our findings will help to predict and monitor areas at risk of future severe cyanobacterial blooms in the Okavango Delta and other similar African ecosystems.