3.1 Statistical Analysis of Nutrient Concentration and Flow
Nitrogen and TP loads varied significantly among sampling locations over the 20-year study period (N: F3,53 = 179.85, P < 0.001; TP: F3,53 = 28.48, P < 0.001). Post-hoc analysis indicated significant pairwise differences in N loads among all the sites investigated (i.e., Sites A and C; A and D; B and C; B and D; C and D [P < 0.001]), but not between sites A and B (P = 0.46). In addition, TP was significantly lower at site D than at sites A, B and C (P < 0.001). No further significant differences among sites were apparent in TP loads (P > 0.05).
Actual N load varied significantly from the allowable load when the sampling sites were considered collectively (F1,2088 = 10.15, P = 0.001) (Fig. 2). Similarly, TP actual load (Fig. 3) was significantly higher than the allowable load in the study area (F1,1953 = 137.12, P < 0.001). N loads varied significantly among the WWTWs over the study period (F1,19 = 8.43, P < 0.05), whereas TP did not vary significantly among WWTWs (F1,19 = 3.63, P = 0.07). Moreover, the actual N loads were significantly higher than allowable loads at both the Baviaanspoort (F1,767 = 12.42, P < 0.001) and Zeekoegat (F1,423 = 566.25, P < 0.001) WWTWs. A similar trend was apparent for TP, with actual TP being significantly higher at both the Baviaanspoort (F1,669 = 1800.57, P < 0.001) and Zeekoegat (F1,385 = 401.24, P < 0.001) WWTWs.
Tab. 2–6 present correlations between the nutrient concentrations and flow using Spearman’s rank test. With linear regression/correlation data, the p- value is good to consider, but it is more important to assess the R value. An R below 0.4 is considered a weak association. Notable, there was a weak, but significant positive correlation between N concentration and flow when the study area was considered collectively (rs= 0.11, n = 2162, p < 0.001). Similarly, TP correlated weakly, but significantly with flow rate estimates (rs= 0.20, n = 982, p < 0.001).
3.2 Flows into the RD System
Fig. 4 depicts the FDC of the cumulative frequency of 2001–2021 flow data recorded at five gauging stations: A2H027, A2H029, A2H054, A2H055, and A2R009. At A2H027, 90% of flow was ≥ 0.44 m3.s-1 and 10% was ≥ 1.76 m3.s-1. At A2H029, 90% of flow was ≥ 0.01 m3.s-1 whereas 10% was ≥ 0.34 m3.s-1. Relatively similar conditions were observed at A2H054 and A2H055, where 90% and 10% flows were ≥ 0.11 m3.s-1 and ≥ 1.14 m3.s-1, respectively. These findings suggested the dominance of low flows in the tributaries, with the Pienaars River (A2H027) being the major flow contributor.
A large volume of water from the Pienaars River into RD is probably effluent from the Baviaanspoort WWTW (30 Mℓ/day full operating capacity) into the river in the immediate upstream reach of A2H027 (DWAF 2008). The second largest contributor of flow to the RD is the Hartbeespruit River (A2H054), which drains the suburban western part of the RD and an industrial area through its tributary, the Moreleta River (A2H055). The third largest contributor is the Edendalespruit River (A2H029), which mainly flows through the agricultural area on the east side of the RD (Bosman and Kempster 1985). In addition, Zeekoegat WWTW, which discharges into the central point of the RD (monitoring site B), is operating at 50Mℓ/day (58% of full capacity) (Burkard and Van der Merwe-Botha 2017).
At A2R009, there was no flow 50% of the time, whereas for 10% of the time the flow was ≥ 4.71 m3.s-1; the site is at the dam wall, and it is associated with periodic dam spills and releases. As a result, high-flow (flooding) events downstream occur when the RD is filled to capacity. This condition is only likely to occur following substantial rainfall in the upstream reaches of the RD catchment during the wet seasons, whereas, in the dry season, waterbodies are susceptible to water quality impacts. This usually negatively impacts ecosystem health and may have significant cost implications as water may require additional treatment for beneficial use.
3.3 Nutrient Loading Capacity
In the early 1980s, the South African Department of Water Affairs and Forestry (now DWS) promulgated special standards (for the purification of effluent) of 1 mg/ℓ for P and 1.5 mg/ℓ for N as a long-term eutrophication control program (DWAF 1998; DWS 2023). However, this approach has been criticized on the grounds that it provides a blanket approach and ignores that there are differences in the nutrient-receiving capacity of waterbodies and that in some catchments the contribution from non-point sources was high enough that the removal of point sources would have ignored the effects on the trophic status of waterbodies (Pretorius 1983; Griffin 2017).
In the 1990s, DWAF adopted the approach of water quality objectives (as a management tool) that focused on the cumulative impacts on the receiving water resources. Such an approach informed the development of RQOs. In cases where catchment nutrient numerical targets are not yet established, national water quality standards may be used. The RQOs are numerical and narrative descriptors of conditions that need to be met in order to achieve the required management scenario as provided during the resource classification and are required for both the quality and flow of the water resource (Odume et al. 2018). Therefore, the results from monitoring sites A–D were compared against the RQOs instead of the specialized standards. Loads above the LDC illustrated the exceedance of RQOs, and those below the curve indicated compliance.
Fig. 5 illustrates the loading capacities for NO3 + NO2 and TP using the RQO targets of 1.0 mg/ℓ and 0.13 mg/ℓ, respectively. Site A depicted an increase in concentration with a decrease in flows for both NO3 + NO2 and TP, which resulted in frequent exceedances of RQOs under dry conditions. This suggested a constant contribution from point source pollution. Frequent exceedances of allowable concentration were noted in Site C between the 10th and 90th percentiles, which could be an indication that the site receives both point and non-point loads. In site D, water quality exceeded allowable concentrations for both NO3 + NO2 and TP occurred mostly under wet conditions, below the 60th percentile.
Unfortunately, the DWS flow data does not include the effluent flows from the WWTW. Therefore, the effluent load from Zeekoegat WWTW was estimated by adding data from gauging stations A2H027, A2H029, A2H054, and A2H055 and subtracting it from A2R009 data. The estimated effluent flow was plotted against the corresponding water quality concentrations collected at the discharge point (Fig. 6A), and the same approach was applied to site B (Fig. 7). In contrast, the effluent from the Baviaanspoort WWTW (Fig. 6B) was calculated by multiplying the flow from the nearest gauging station (A2H027) by the concentration collected at the discharge point. Importantly, Fig. 6 illustrates a random relationship between nutrient concentrations in WWTW effluent loads into the RD system. Notably, loads increased with discharge before the dilution effect occurred.
Zeekoegat WWTW recorded the highest loads of 3341.8 kg/d at 5 m3.s-1 for N and 3864.9 kg/d at 15 m3.s-1 for TP, and Baviaanspoort WWTW recorded the highest loads of 9193.30 kg/d at 3 m3.s-1 for N and 6487.46 kg/d at 8 m3.s-1 for TP. High effluent loads could be associated with the malfunction of the two WWTWs.
It was observed that high effluent loads were experienced before the system could assimilate the pollution load. In this regard, self-purification, also known as assimilative capacity, is a natural biogeochemical process occurring in any ecosystem and leading to the elimination or assimilation of nutrients or other pollutants by the natural activity of its resident biological communities (Breil et al. 2022). Generally, river flow facilitates assimilative capacity and mitigates the impact of pollution (Farhadian et al. 2015).
The US-EPA (2007) indicated that streams receiving effluent from WWTWs exert a significant influence on water quality at low flows. This was the case at sites A (Fig. 6) and B (Fig. 7), which receive effluent from WWTWs. For instance, it was observed in Fig. 7 that load exceedance for both N and TP occurred during low flows, with N reaching a high level of 1399.52 kg/d at approximately 5 m3.s-1 compared to a level of 225.30 kg/d observed at 15 m3.s-1 and TP reaching a level of 530.57 kg/d at 6 m3.s-1. Whereas the highest TP load of 637.29 kg/l was observed during high flows (15 m3.s-1). These results could be attributed to the discharge from the WWTW; however, a further contribution could be from non-point sources in the catchment.
When the actual load exceeds the allowable loading under dry or low flow conditions, a contravention of the water quality standard occurs (Teague et al. 2011). In this context, water quality standard exceedance near high flows (flow exceedance range <40%) was associated with rainfall events and pollution sources were classified as non-point sources; conversely, water quality standard exceedance near low flows (flow exceedance range > 60%) was associated with dry weather conditions, and pollution sources were classified as point sources. Thus, this study concurred with other studies (Bowes et al. 2008; Slaughter 2011; Slaughter and Hughes 2013) that water quality for sites where there are known point sources is rarely associated with a clear asymptotic decrease in P concentrations with increasing flows as point source load decreases with increasing flow due to dilution of constant input, whereas non-point load usually increases with overland flows.
Dlamini et al. (2019) assessed the compliance of RQOs at a wastewater facility by investigating the effects of nutrient (NO3- and P) loading capacity at the catchment level. The authors revealed that the management of nutrient loads during low flows is crucial, and that compliance status should be based on the flow regime for each season. The findings of the current study concurred with Dlamini et al. (2019) and, in addition, recommended that pollution be controlled at the source to minimize impacts on the water resource.
N and P are major nutrient parameters in water resources, and TP, TN, and biodegradable organic matter are considered critical pollution indicators and the best estimates for pollution loading (Cheruiyot and Muhandiki 2014). Nutrients can be transported to waterbodies through landscapes (Mockler and Bruen 2018). For instance, nitrate (NO3-) is transported to streams via sub-surface pathways, whereas P from diffuse sources is driven by storm events and transported through overland flows (Jordan et al. 2005; Kröger et al. 2007; Tesoriero et al. 2009; Mockler and Bruen 2018). Mining and utilization of P minerals in fertilizers and detergents have caused major changes to the global biochemical P cycle in water resources (DWS 2023). Rapid aquatic plant growth is a common effect of nutrient enrichment in freshwater ecosystems (Yao et al. 2018) and can lead to toxic algal blooms, which constrict water security (Wei et al. 2022).
Fundamental to trophic status assessment is the determination of the relationship between the level of nutrient loading (P and N in particular) and the in-lake conditions (Harding 2008). Thus, the attainment of water quality management objectives is key to avoiding or reducing water quality problems such as eutrophication. This study further recommended that P and N-special standards be updated according to the characteristics of water resources and that water quality be evaluated in terms of desired beneficial water use and per system or catchment, considering the impactors.
3.4 Estimated Nutrient Load Reduction
Tab.7 presents the estimated annual mean loads measured across the RD from 2001–2021. The growing trend on N was noted between 2015 and 2021, with the highest load of 414.74 t/a recorded in 2021. On the other hand, the TP load was observed to be fluctuating through the period of 2001–2021, with the highest load of 173.18 t/a being recorded in 2020.
Notably, the annual TN load exceeded the allowable load 38% of the time, whereas the annual TP load exceeded the allowable limit by 95%. Therefore, the RD was characterized by long-term variation in nutrient load, which could be related to the fact that the catchment was impacted by both point and non-point sources. However, the major contributors were point sources, as depicted in Fig. 6 and 7.
The findings in Tab. 7 were used to estimate the nutrient load reduction. Tab. 8 below presents the proposed total, median annual nutrient load limits (TMANL) for which the RD could reduce nutrient loads in order to meet the desired trophic state. Thus, instream nutrient concentrations must be improved to maintain aquatic ecosystem health, ensure ecological sustainability, and meet all water quality requirements for water users. This can be done by reducing N by 15.3% in order to meet the numerical limit of 1 mg/ℓ and TP by 66.27% to fulfil the 0.130 mg/ℓ requirements.
According to DWAF (1998), approximately 35–50% of the P in domestic wastewater and most of the P in industrial wastes comes from synthetic detergents, and 32% ends up in reservoirs (Quayle et al. 2010). This makes P the most manageable nutrient. While N removal is more complex, for instance, the conventional wastewater N removal process is an energy-intensive process due to the aeration demand (Jia and Yuan 2016). On the other hand, (Oberholster et al. 2019a) suggested a specific algae treatment for domestic wastewater effluent from pond systems as an alternative practice to improve the quality of the effluent. A reduction of 43.1% and 35.1% for TN in the two respective ponds was observed, indicating other possible effective and sustainable methods to be used to manage N loading in African countries. Thus, this current study recommended that each WWTW in the vicinity of the RD be upgraded to remove at least 8% of the N load and 30% of the TP load per facility per year. This can be achieved by the application of nature-based solutions to control pollution and reduce stormwater and runoff input, as well as employing low-cost green treatment technology options to reduce nutrient loads from domestic wastewater effluent, such as using a consortium of microalgae in municipal wastewater treatment pond systems, which requires less energy. Oberholster et al. (2021) presented compelling evidence in terms of the feasibility of using this technology in developing countries to reduce nutrient loads from domestic wastewater effluent.
In addition, rehabilitation programs for dams in developing countries and dams impacted by different sources of pollution, such as RD, would need to focus on the management of both point and non-point sources, taking into consideration the following actions:
- Enforcement of stringent regulatory actions to ensure that the surrounding WWTWs are complying with discharge permits.
- Adoption of the best management practices for non-point pollution sources such as nutrient tracing programs, nutrient and path management programs for fertilizers, and construction erosion and sediment control ordinances.
- The introduction of zero-phosphate detergents in African countries linked to awareness and education campaigns through citizen science.