1. Introduction
Globally, freshwater resources have been put under increasing pressure due to the rapid increase of the world population and their improvement of living standards [
1]. It was indicated that an estimation of more than 83% of land surfaces that are surrounding freshwater systems have been significantly influenced by the footprint of humanbeing as a response to anthropogenic activities [
1,
2]. Contaminants from ineffective waste management, pesticides and fertilizers from agricultural areas, pollution from urban, industrial and domestic wastewater can often released to ground and surface water [
3], and even the landfills or slag heap disposals may release pollutants seeping into nearby water resources [
4,
5]. The inland freshwater bodies are more vulnerable to such problems of pollution and contamination [
6] resulting in the decline of quality and availability [
7]. Among the loaded pollutants, some serve as nutrients for the enhancing growth of algae. Such nutrients can be loaded from point or non-point sources, point sources, for example, wastewater treatment plants can easily be recognized and given more attention [
8,
9]. Non-point sources including urban areas, cultivated lands, natural forests and pastures are more difficult to be recognized and attention should be given to them according to the amounts of nutrients they can potentially release. For example, it was found that the estimated loads of nitrogen, phosphorus, and suspended sediments from urban areas and cultivated land are 10 to 100 times greater than idle and forested lands in the great lakes catchments of USA and Canada [
10].
Eutrophication is a serious freshwater problem caused by the excessive growth of harmful algal blooms (HABs). In recent years, the factors enhancing such blooms have become interesting research subjects. Many studies have been conducted in this field, and most of them were related to the HABs to the loading of nutrients such as nitrogen and phosphorus into the water bodies [
11,
12].
In Africa, societies are highly dependent on their natural inland freshwater resources creating more pressure on them resulting in significant changes in their water quality [
13]. South Africa has over 4000 freshwater storage dams, among which around 700 are public dams controlled and managed by the governments used for domestic, irrigation and industrial water supply [
14], they are underpinning the economic and social development of the country. Despite of such a huge number of dams, South Africa is water stressed and scarce country facing critical challenges due to poor water management practices, inadequate infrastructure, and a relentless surge in water demands [
15]. Owing to the fact that most South African dams are located downstream of metropolitan and urban areas, they have become more enriched with nutrients [
4,
16]. The climatic conditions of South Africa associated with nutrient loads resulted in extensive and widespread eutrophication and cyanobacterial blooms in the inland freshwater bodies [
17,
18]. This will continue to increase the cost of using these valuable resources. Nitrogen and phosphorus are considered as the leading factors in accelerating the lakes’ and reservoirs' eutrophication [
19]. The Water Act of 1956, section 21(1) (a) and the Department of Water Affairs in 1980 set exceptional standards of the phosphate concentration for effluent discharged to some mentioned rivers or their tributaries, it promulgated that the phosphate concentration should not be higher than 1 mg/l [
15,
20]. An investigation was conducted by CSIR and the Department of Water Affairs from 1985 to 1988 to predict the impact of the phosphorus standard (1mg/l). Based on their findings, the eutrophication control aim was set to maintain the mean concentration of chlorophyll-a within receiving water bodies at the level that no conditions of severe nuisance would occur in more than 20% of the time. To achieve this aim it requires to keep the phosphorus concentration level in the reservoirs water less than 0.14 mg/l [
21].
Vaal Dam is one of multiple water bodies affected by eutrophication and cyanobacterial blooms in South Africa. It was constructed by the Department of Irrigation of the national government and functionally completed in October 1936 at 20 km downstream of the confluence of Vaal and Wilge rivers [
22]. It is the second biggest dam by surface area in South Africa (about 320 Km
2). The dam has undergone several rises to increase its storage capacity which ended up with the total capacity of approximately 2.603 × 10
6 m
3 [
22]. The dam catchment area is approximately 38000 km
2 holding various human activities including major agricultural activities (crop cultivation and cattle grazing), mining and some industrial activities [
23,
24] as well as many formal and informal settlements. In many areas within the Vaal Dam catchment, mine dewatering and urban/industrial treated effluent discharges find their way to the streams causing serious water quality issues [
23]. Such activities have direct or indirect effects on the dam water quality in terms of nutrient loading which enhances the growth of HABs. A study in 2013 stated that the Vaal Dam Reservoir was classified as a mesotrophic water body according to the South African Department of Water Affairs Classification System, the mean concentration of chlorophyll-a was 14.8 μg/l, mean total phosphorus concentration was 0.077 mg/l and the time percentage where chlorophyll-a exceeds 30 μg/l was found 17% [
25].
Harmful Algal Blooms (HABs) have major effects on water quality and their aquatic system function, therefore, monitoring of their distribution in space and time is very important for water resources managers to address the issues related to it. Moreover, it is very important to address the question of which factors enhance and control the HABs in the Vaal Dam Reservoir. Historically, many studies were conducted to assess different methods of detecting HABs and cyanobacteria in Vaal Dam Reservoir, for example, an attempt was made to help the managers of the drinking water treatment facility with advanced prediction of
Microcystis sp. concentration in the Vaal Dam water, by building a model using physical, chemical, and biological water quality records between 2000 and 2012. The model showed a promising result of estimating
Microcystis sp. in 7-days in advance [
25]. Other studies were conducted using remote sensing to estimate the spatial distributions of HABs in the dam reservoir using satellite imageries, different algorithms and band ratios have been tested for Lndsat-8 and Sentinel-2 data. These have successfully revealed HABs predicting indices with strong correlations to the in-situ data [26-28]. The latest remote sensing based published article in the Vaal Dam (Obaid et al., 2021) applied to high resolution sensors that successfully used both the blue-green and red-infrared OC algorithms in estimating Chl-a concentrations [
16,
29]. As known, the blue-green algorithms are usually used to retrieve Chl-a concentration in Case I water [
30,
31]. In Case II waters, various constituents absorb light in the blue region which might create uncertainty when using it to detect chlorophyll-a concentration in such productive waters [32-34]. The successful use of blue-green OC algorithm may be because of the strong signal of HABs from the high biomass concentration in the reservoir water column. However, all above mentioned studies were focused on developing ways to detect HABs and cyanobacterial blooms in the reservoir, none of them tried to reveal the factors enhancing such blooms in terms of nutrient loading and environmental factors.
This study aims to investigate the relationship between HABs, and the potential nutrients and environmental processes during the times of peak blooms within the last few decades using historical water quality records from the dam reservoir, this can be obtained by comparing the records of nutrients with the Chl-a concentrations measured in the dam to identify which nutrients enhancing algal blooms. Moreover, the study tries to conduct spatial distribution of HABs using Landsat-8 and Sentinel-2 satellite images captured in late summer (April) between 2013 to 2023, which will help to identify the frequently and most affected parts of the dam reservoir. The Remote sensing approach is useful in this aspect by giving a synoptic view of HABs spatial and temporal distribution and improving our understanding of its dynamic variations. It has the potential to allow effective monitoring of HABs with high resolution data for more precise studies.
2. The Study Area
The Vaal Dam holds water used to supply potable water to the Gauteng metropolitans and its surrounding areas. The Vaal Dam Catchment extends within Free State, Mpumalanga, and Gauteng provinces and drains by Vaal and Wilge river systems, (
Figure 1). The main rivers consist of many tributaries draining different areas in terms of land use land cover types and human activities. For example,
The Waterval River contributes approximately 111 × 10
6 m
3 of water annually to the Vaal River. It drains very active areas holding intense human activities such as agriculture, industry, mining, urban and rural settlements, mainly in the upper reach of the Waterval River, these activities have been shown to be responsible for the deterioration of the water quality and ecological integrity of the river system [
35]. Another active area with mining and industrial activities within the Vaal River site is Grootdraai Dam catchment, it’s found that the greatest area of concern was the region’s closeness to the downstream of the urban, industrial, mining and cultivated lands cover [
36]. These two active areas in terms of land use activities put the Vaal River under more concern about its water quality issues.
The Wilge River system drains areas dominated by agriculture and grasslands, it contributes a great deal of water to the dam from its catchment and from the Lesotho Highland Water Project [
37]. The Tugela-Vaal Water Project also contribute a good portion of water to Vaal Dam via the Wilge and Nuwejaarspruit Rivers [
22].
5. Discussion
The data were explored to uncover the important water quality properties and the perusal of the graphs reveals some extremely high values as well as some trends and seasonality in temperature and dissolved oxygen values. The distribution frequency plot of chlorophyll-a (
Figure 2) after excluding the outliers (very high values), shows that most of the data centered near the low range values, below 10μg/l. This reveals that most of the time, the chlorophyll-a concentration in this measurement point remains in its lower concentrations.This might be because of the location of the measuring point which lies in the Vaal River near the dam wall where the water runs towards the discharge gates, in such situation the algae might be mixed and washed out comparing to other parts of the reservoir where the conditions are ideal for algal growth. The distribution of TP, NH
4_N and KJEL_N are also showing the measures are centered relatively near the lower values after removing some stochastically high values. These indicate that their concentration in the reservoir water is, most times, within their low limits with some periods of extremely high to very high values which probably indicate that some types of pollution ended up reaching the reservoir by increasing the nutrient concentrations which apparently increases the capacity of the reservoir water to support high rates of biomass productivity during such high nutrient load times.
The productivity in the open water systems is a function of multi biophysical factors, such as light, temperature, nutrients, etc. In many aquatic systems, some limiting nutrients mainly TP and N are considered the main drivers of HABs [
39,
40]. However, the time series plot (
Figure 3) showed that the TP levels between 1990 and 2022 are relatively low except on some specific dates, a corresponding increase in Chl-a and KJEL_N concentrations at such specific dates were detected, which jumped to extremely high values. These ar clearer when we look closely at the decadal time series plots in
Figure 4b,c. Thus, this explains that the dam water HABs productivity is primarily driven by TP and KJEL_N when their concentrations suddenly jumped to high values. The water temperature time series in
Figure 3 and
Figure 4a and 4c show no significant change and its decadal trends show a slightly decreasing trend for the first decade (1990 to 2000), see
Figure 5a while showing no increasing or decreasing trend during the last decade (2010 – 2020) in
Figure 5c, the average temperature for the first and last decade was 17.94 and 18.04 °C, respectively which suggest a slight increase in the average water temperature during the last decade.
Chl-a decadal trends (
Figure 5a–c) showed a slight decrease in the first decade (1990 to 2000) and remained constant within the second decade (2000 to 2010) before it increased noticeably during the last decade (2010 to 2020). The average decadal concentrations were 4.75, 10.51 and 16.7 μg/l respectively, with a significant increase throughout the last three decades. The low concentration of Chl-a between 1990 and 2000 may be because of the implementation of the bioremediation project between (the late 1980’s to the early 1990’s) which reduced the Chl-a levels effectively [
41].
In general, except for the first decade, the mean Chl-a concentration is above the threshold for eutrophic systems (7 μg/l) [
16]. This situation of increasing the decadal average of Chl-a during the study period while the temperature remains with no significant change suggests that the productivity in the Vaal Dam is not limited to the changes in the temperature during the study period.
TP, KJEL_N and NH
4_N decadal trends followed the same trend behavior of the Chl-a in the first and last decades, they were associated with general upward and downward trends of Chl-a. The average decadal concentrations of TP were 0.1043, 0.1096 and 0.1119 for the first, second and third decade respectively. These concentrations show that the TP concentrations were low except on the above mentioned individual dates where it jumps to high levels greater than the hypertrophic TP threshold level which is 0.25 mg/l [
16]. For KJEL_N, the decadal average concentrations were 0.8033 and 1.1417 mg/l for the first and last decade, respectively, and NH
4_N were 0.0397 and 0.0431 mg/l, respectively. However, the corresponding behavior of the Chl-a, TP and KJEL_N decadal trends alongside the erratic high Chl-a values recorded at the same individual dates of very high TP and KJEL_N records, explain that they are the driving factors of the algal blooms within the Vaal Dam.
On the other hand, the NO
3NO
2_N concentrations showed a negative trend compared to the Chl-a decadal trends during the first and last decade and DO also show a negative trend during the last decade. The NO
3NO
2 average decadal concentrations were 0.2455 and 0.2248 mg/l respectively with the average concentration of DO 7.928 mg/l. The trend of the DO is well understandable because dissolved oxygen is usually consumed during excessive algal blooms which have been known through last few decades [
42], but this situation might not be applicable to the behavior of NO
3NO
2_N trends. The analysis of the results suggests that it has no direct relation with the Chl-a trend, but a paper [
43] focused on sources and forms of most important nitrogen substrate for blooms in eutrophic Lake Erie suggested that, the NO
3ˉ was the most important source of N except in late blooming stages where phytoplankton relays on recycled N derived from dissolved organic nitrogen. Their study further showed that the NO
3ˉ depletion was related to the consumption by phytoplankton during its blooms showing a negative relationship. In this study NO
3NO
2_N also showed a negative correlation with Chl-a, but the assumption of consumption by the HABs needs more detailed investigations. The regression between the targeted WQ parameters showed a strong correlation between Chl-a and TP, Chl-a and KJEL_N in great agreement with the trend analysis. It also showed a positive correlation between Chl-a and temperature while there was no correlation between Chl-a and NH
4_N which explains that NH
4_N is not a driving factor for HABs blooms in the Vaal Dam. Like the decadal trend plots, the regression showed a negative correlation between the Chl-a and DO & NO
3NO
2_N.
The distribution of HABs in the reservoir area was analyzed by the interpretation of Chl-a concentrations derived from satellite data, some extreme concentrations have clearly been seen in 2015 and following images where the concentration of Chl-a was much more than those in summer 2013 and 2014. This change is strongly correlated with some high in-situ measured values of Chl-a, TP and KJEL_N during April and May 2015, and January 2016. High productivity was also noticed in images of 2019, 2021, 2022 and 2023. These high productivities may relate to the periods where TP and KJEL_N levels increase due to the increase of human activities such as discharge of partially treated or untreated wastewater, or through runoff from urban areas around the dam, and agricultural areas. The satellite data show high productivity in the Vaal Dam for the period from 2019 up to 2023 while there are no historical records available on the website during this time which put some warnings of increasing nutrient loads on the dam reservoir recently. The geometry of the Dam and the fluxes of the Vaal and Wilge rivers draining into it are directly linked to HABs spatial dynamics.
Within the past decade (2010 – 2020), the media reported some waste found its way to the reservoir. For example, on July 23
rd 2015 a report mentioned that an uncontrolled sewage discharge was overwhelmed the Deneysville town which is located on the Vaal Dam bank, just next to the dams wall (
https://mg.co.za/article/2015-07-23-sewage-in-gautengs-drinking-water/), accessed March 7, 2023.
Author Contributions
Conceptualization, Altayeb Obaid; Data curation, Altayeb Obaid; Formal analysis, Altayeb Obaid; Investigation, Altayeb Obaid; Methodology, Altayeb Obaid; Resources, Altayeb Obaid; Software, Altayeb Obaid; Supervision, Elhadi Adam, Khalid Adem Ali and Tamiru Abiye; Validation, Altayeb Obaid, Elhadi Adam, Khalid Adem Ali and Tamiru Abiye; Visualization, Altayeb Obaid, Elhadi Adam, Khalid Adem Ali and Tamiru Abiye; Writing – original draft, Altayeb Obaid; Writing – review & editing, Altayeb Obaid, Elhadi Adam, Khalid Adem Ali and Tamiru Abiye.