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 similar to 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.