Development of a Universal Water Quality Index (UWQI) for South African River Catchments

The assessment of water quality has turned to be an ultimate goal for most water resource and environmental stakeholders, with ever-increasing global consideration. Against this backdrop, various tools and water quality guidelines have been adopted worldwide to govern water quality deterioration and institute the sustainable use of water resources. Water quality impairment is mainly associated with a sudden increase in population and related proceedings, which include urbanization, industrialization and agricultural production, among others. Such socio-economic activities accelerate water contamination and cause pollution stress to the aquatic environment. Scientifically based water quality index (WQI) models are then essentially important to measure the degree of contamination and advise whether specific water resources require restoration and to what extent. Such comprehensive evaluations reflect the integrated impact of adverse parameter concentrations and assist in the prioritization of remedial actions. WQI is a simple, yet intelligible and systematically structured, indexing scale beneficial for communicating water quality data to non-technical individuals, policymakers and, more importantly, water scientists. The index number is normally presented as a relative scale ranging from zero (worst quality) to one hundred (best quality). WQIs simplify and streamline what would otherwise be impractical assignments, thus justifying the efforts of developing water quality indices (WQIs). Generally, WQIs are not designed for broad applications; they are customarily developed for specific watersheds and/or regions, unless different basins share similar attributes and test a comparable range of water quality parameters. Their design and formation are governed by their intended use together with the degree of accuracy required, and such technicalities ultimately define the application boundaries of WQIs. This is perhaps the most demanding scientific need—that is, to establish a universal water quality index (UWQI) that can function in most, if not all, the catchments in South Africa. In cognizance of such a need, this study attempts to provide an index that is not limited to certain application boundaries, with a contribution that is significant not only to the authors, but also to the nation at large. The proposed WQI is based on the weighted arithmetic sum method, with parameters, weight coefficients and sub-index rating curves established through expert opinion in the form of the participation-based Rand Corporation’s Delphi Technique and extracts from the literature. UWQI functions with thirteen explanatory variables, which are NH3, Ca, Cl, Chl-a, EC, F, CaCO3, Mg, Mn, NO3, pH, SO4 and turbidity (NTU). Based on the model validation analysis, UWQI is considered robust and technically stable, with negligible variation from the ideal values. Moreover, the prediction pattern corresponds to the ideal graph with comparable index scores and identical classification grades, which signifies the readiness of the model to appraise water quality status across South African watersheds. The research article intends to substantiate the methods used and document the results achieved.


Introduction
The water quality index (WQI) is the most popular method of exhibiting the water quality of surface water bodies. WQI models are better known for delivering a comprehensive and explicit representation of water contamination for both surface water basins and groundwater reservoirs. The appraisal concept is concise and more straightforward, leading to wide acceptance across the water science community [1]. WQI provides a single numeric value that expresses the status of water quality through the integration of multiple microbiological and physico-chemical parameters [1,2]. Water quality index scores are classified using a diverse array of rating scales, but the frequently used grading system ranges from zero (bad quality) to one hundred (excellent quality) [2][3][4][5]. WQI scores are dimensionless [6], and can be further categorized using descriptive ranks associated with terms like "poor", "marginal", "fair", "good" and "excellent" [3][4][5]7]. Water quality indices (WQIs) are typically used by water authorities, water scientists, policymakers and the general public for decision-making, delineating spatial and temporal trends, tracing contamination sources, appraising regulatory guidelines and environmental policies and, most importantly, for suggesting future recommendations [6,8].
The main objective of WQIs is to convert multiple parameter data into information that is understandable by both technical and non-technical personnel. The ability of WQIs to synthesize complex scientific data into simple and easily understood formats makes them the most fundamental and indispensable elements of water quality monitoring agenda. Therefore, they are universally acknowledged as a "lifeline" for water quality studies, and their development continues as an ongoing affair [7]. Despite their range of applications and the variety of WQIs developed this far, there is still no definite and commonly acceptable methodology for developing water quality indices [6,9]. Instead, numerous techniques and approaches exist in WQI formation, but the conventionally employed method involves (a) the determination of relevant water quality variables, (b) the establishment of sub-indices, (c) the generation of significant weightage coefficients, (d) the aggregation of sub-indices and, lastly, (e) the attribution of a water classification schema [7,8,[10][11][12]. Each step has alternative methods to consider, which mean that it is extremely important to decide the most suitable method for each scenario. Notwithstanding technical knowledge of WQIs, the developer should apply due diligence and avoid subjective judgements and bias in the process of establishing WQIs, otherwise the WQI will inherit abnormalities and be considered dysfunctional [7]. For the current study, an index model for water pollution control and river basin planning functions has been established using expert opinion in the form of the participation-based Rand Corporation's Delphi Technique and extracts from the existing literature. The process yielded thirteen input variables, namely NH 3 , Ca, Cl, Chl-a, EC, F, CaCO 3 , Mg, Mn, NO 3 , pH, SO 4 and turbidity (NTU). In addition to the parameter selection, expert opinion was also applied to develop significant ratings and parameter weightage coefficients. The universal water quality index (UWQI) model is an increasing scale index founded on the weighted arithmetic sum method with resultant values ranging from zero (very bad quality) to one hundred (good quality). The overall classification is centered on five categories, with the Class 1 rank denoting "good water quality" and Class 5 rank representing very bad water quality. Following the review by Banda and Kumarasamy [7], it has been noted that most WQIs are designed for particular region and are source specific, thus creating a gap and ample scope to develop a universally acceptable WQI. However, it is a demanding task and is extremely difficult to develop a water quality model that is globally acceptable, and hence the current studies only focus on national boundaries-that is, they focus on a model only applicable to South African river catchments. Though this prospect is seemingly problematic to deal with, it is pertinent and recommended that water quality experts embark on developing a unified model that can be utilized across the globe. However, the immediate mission is to develop nationally acceptable water quality indices and break the barrier of region-specific models [7]. Moreover, this study attempts to break such barriers through the development of a universal index that is applicable to most river catchments in South Africa, thereby promoting a standardized way of monitoring and comparing the water quality of various watersheds at a national level, which might eventually assist in the prioritization of water resources across all nine provinces in South Africa. Umgeni Water Board (UWB) provided the water quality dataset used to test the UWQI. The data are from six sampling stations located in four different catchments under the jurisdiction of Pongola-Mtamvuna Water Management Area (WMA) in KwaZulu-Natal Province, South Africa. The four watershed regions are Umgeni, Umdloti, Nungwane and Umzinto/uMuziwezinto River catchments. The UWQI is earmarked for national application, but it is far-reaching and beyond the scope of the study to test the model against data from all 148 catchment regions in South Africa. Nevertheless, the four catchments used are adequate to ascertain the functionality of the model and the process is a step towards the ultimate goal of testing the model against most, if not all, the catchment areas in South Africa. The model responded steadily to the variation in parameter values and managed to indicate spatial and temporal changes in water quality for the four catchment areas considered for the study. It should be noted that the UWQI is formed autonomously without being linked to a particular dataset or specific region. The methods used are exclusively independent from such associations and the UWB dataset is entirely for testing purposes, upon which tasks can be performed using any other available data.

Research Data
The water quality monitoring process demands substantial amount of resources; therefore, the current study could not collect samples specific to the research work. Alternatively, water quality data from Umgeni Water Board (UWB) assisted in testing the functionality of the model. The dataset was comprised of 416 samples, tested monthly for a period extending to four years. The water quality records are from six sampling stations located in four different catchment areas, namely Umgeni, Umdloti, Nungwane and Umzinto/uMuziwezinto River catchments.
The UWB data were sampled in accordance with standard methods prescribed by the Department of Water and Sanitation (DWS), and analyzed according to international standards in an ISO 9001 accredited laboratory owned and operated by UWB [13]. The research dataset from UWB satisfactorily provided all the required thirteen water quality parameters, and these are ammonia, calcium, chloride, chlorophyll-a, electrical conductivity, fluoride, hardness, magnesium, manganese, nitrate, pondus Hydrogenium, sulphate and turbidity. Testing the model with data from these four river catchments supports the objective of establishing a UWQI applicable to a greater part of the country, if not the whole of South Africa. More than the availability of data from UWB, the economic significance of KwaZulu-Natal Province [14,15], the distinctiveness of its inter-basin arrangements, the scope of the transfer schemes involved and the extensive water demand [16,17] encouraged the choice of the study area, which falls under Pongola-Mtamvuna water management area (WMA) [18,19]. The project data were adequate to examine the model and complement the objective of developing a universally acceptable water quality model.

Universal Water Quality Index (UWQI)
Various methods are documented in the literature and, among them, there is no one distinctive method regarded as the supreme and favorable method for developing water quality indices (WQIs). Each method has its own considerable problems and the universal water quality index (UWQI) was formulated using the conventional method of establishing water quality indices. Moreover, the technique involves four common steps, which are (1) selecting water quality variables, (2) establishing relative weightage coefficients (3) forming sub-index rating curves and sub-index functions and (4) deriving the appropriate aggregation or indexing model [8,[10][11][12]. The methods employed for the development of the UWQI are selected based on a couple of reasons. Firstly, they eliminate individual bias through the incorporation of objective and subjective opinions from water quality scientists through appraisal questionnaires. Secondly, compared to other available techniques, the chosen methods are both practical, convenient and easy to implement in electing variables and generating weightage coefficients [20]. Lastly, the methods are proven and have been performed in various WQI studies [21][22][23][24][25][26][27][28][29][30][31][32][33][34][35].

Selection of Water Quality Variables
These steps and procedures were performed cautiously with cognizance of the fact that the model should widen its application boundaries and target to become a nationally accepted water quality monitoring tool. Based on this fact, a fixed set of parameters were established using expert opinions. The advantage of a fixed set of variables is that the model can be applied in various catchments without the possibility of altering the structure and functionality of the model [9], thereby permitting stakeholders to fairly compare the water quality of different sites and develop a more informed national prioritization schedule without prejudice. Further to this, expert opinion has the advantage of promoting the acceptability of the model, in the sense that most of the experts engaged are also the targeted end users of the model; therefore, the idea that they were involved in the process of developing the UWQI may eventually facilitate acceptance through a sense of ownership. Nevertheless, this alone does not warrant the usefulness of the model, the authors exercised enormous care and great attention to ensure that the most significant variables were incorporated in the UWQI. Importantly, the authors had to optimize the ideal number of parameters necessary to provide a meaningful water quality index value.
Following the Rand Corporation's Delphi Technique, a panel of thirty water specialists from government parastatals, the private sector and academia was established. Delphi Questionnaires were circulated to the participants and they were asked to consider twenty-one water quality parameters for their possible inclusion in the UWQI. The panelists were instructed to designate each variable as: "Include" and "Exclude" and further assign a relative significance rating against each variable elected as "Include." The rating scale ranged from one to five, whereby "scale 1" denoted the highest significance and "scale 5" represented a comparatively low significance. In addition to the prescribed twenty-one parameters, the experts were allowed to add, at most, five more variables if desired. A total of twenty-one questionnaires were returned out of the thirty questionnaires circulated. The Rand Corporation's Delphi Technique is described in detailed by Horton [21], Brown et al. [22] and Linstone and Turoff [36,37] and applied in several studies, which include those by Nagels et al. [33], Kumar and Alappat [34] and Almeida et al. [35].
Complementarily, the existing literature on WQIs was used to select the most significantly used water quality variables. Thirty-seven studies were considered and each variable was designated as "Include" if it corresponded to the twenty-one parameters considered for the Delphi Questionnaires; otherwise, it was designated as "Not Included." Furthermore, the formerly assigned significance rating was adopted as the relative significance rating for each parameter that was "Included" in the study in question. The rating was based on a scale ranging from one to five, with "Scale 1" representing the lowest significance and "Scale 5" for relatively high significance. If a different significance rating scale was used in the existing study, the original rating values were equivalently transformed to match the preferred rating scale using Equation (1), as follows: where y is the new rating; a, b, are minimum and maximum values of the targeted significance scale rating; x min , x max , are minimum and maximum possible ratings in the specified significance scale; x i is the ith rating value of the specified scale. Finally, a holistic ranking order was derived from the combined effect of the two aforementioned methods, upon which a rejection rationale was employed to eliminate redundant variables that are not commonly monitored across South African river catchments [20].

Establishing Weight Coefficients
Each parameter has different effects on water classification; hence, weighting factors are used to reflect their influence on the index model. These mathematical tools are assigned to each water quality variable based on the level of significance of the overall index value [9,38]. Parameter significance ratings assigned on Delphi questionnaires and those extracted from the literature on WQI publications were considered as preliminary significance ratings (see Section 2.2.1). Parameter significance ratings (b i ) were then established by aggregating the preliminary ratings from the two methods. Relative weight coefficients (w i ) are directly proportional to the significance ratings and they were established from dividing the parameter significance rating value (b i ) by the sum of all ratings ( b i ) using Equation (2) [3]. The weight coefficients are presented as decimal figures with a total sum of unity, the reason being that the combined effect of the water quality parameters should not exceed one hundred percent [3]. If this does occur, the aggregation of sub-indices will be compromised, and the water quality index will be deemed dysfunctional. Expert opinion techniques were employed primarily to produce comparative weights which minimize prejudice and uphold the integrity of the index model.
where b i is the assigned significance rating of the ith water parameter (one being the lowest rating and five the highest rating); w i is the final weight coefficient for the ith water parameter (decimal value); n is the total number of the rated water quality parameters.

Formation of Sub-Indices
Considering that water quality parameters are monitored in different scientific units, sub-indices are applied to convert the different units of measure into a single common non-dimensional scale [7]. This is a common practice and the conventional method involves sub-index rating curves which are later transformed into mathematical functions commonly known as sub-indices. For practical purposes, fixed key points of the rating curves were graphically established with reference to the permissible concentration limits. Straight-line graphs were used to converge the plotted points and produce a series of linear graphs, which were further converted into linear sub-index functions. Target Water Quality Ranges (TWQRs), as prescribed by DWAF [39][40][41], were consulted in the process.

Aggregation Formula
A scenario-based analysis was used to modified and align the model with local conditions to develop the final universal water quality index (UWQI), which is an improved version of the weighted sum method. The model equation integrates sub-index values of selected parameters in relation to the assigned weights and obtains the overall water quality status, which is presented as a unitless number ranging from 0 to 100. The rationale employed is based on solving multiple systems of equations [42], where key points of the rating curves were used to generate a series of m equations, with two unknown variables (x, z) and n water quality parameters in the form where WQI 1, . . . , m are the ideal water quality index values corresponding to the key points of the rating curves; x 1, . . . , m are the equation denominators (first unknown variable); SI m1, . . . , mn are the corresponding sub-indices; w 1, . . . , n are relative weight coefficients for the thirteen water quality parameters and z 1, . . . , m are the equation exponentials (second unknown variable). The first step was to find the optimum values of x and z; thereafter, the closest x-value was rounded off and substituted into the same set of equations to find the corresponding optimum z-value, which becomes the final exponential factor of the UWQI.
Using the thirteen selected water quality variables, weightage coefficients and sub-indices, the improved weighted sum method proved to be the most appropriate and relevant method to develop a UWQI for assessing water quality in South African river catchments. Hypothetically, this advocates the readiness of the UWQI model and deems the study a success. Such a milestone provides a tool that can be adopted at the national level to help solve the challenges experienced by water quality professionals. The structure of the universal water quality index model is represented in Figure 1.

Water Classification
In the interest of simplifying the interpretation of water quality index (WQI) values, mostly to accommodate non-technical individuals, an index categorization schema was established. The classification mechanism is based on an increasing scale index and the advantage of this system is that it is identical to a normal percentage hierarchy [3]; therefore, the public can easily relate to its function and interpretation. The UWQI model yields WQI values between zero and one hundred. Accordingly, the WQI scores are categorized using classes ranging from one to five, with "Class 1" representing water of the highest degree of purity with a possible maximum score of one hundred and, vice versa, "Class 5" denotes water quality of the poorest degree with index scores nearing or equal to zero. In order to close gaps identified in some of the existing classification scales [7], appropriate mathematical functions with logical linguistic descriptors, which include, but are not limited to, ʺgreater than", ʺless than", and ʺequal to" were used to appraise WQI scores and we respectively assigned them to the corresponding category.

Background and Specific Considerations
A substantial increase in population and improper disposal of wastewater have a significant influence on diminishing water quality in rivers and other surface water reservoirs. As a consequence, routine water quality assessment and pollution control measures are necessary to preserve and restore the healthiness of surface water bodies [43]. On the same basis, this study attempts to put forward a practical and standardized tool that can be used for monitoring surface water quality across all South African river catchments. Despite the fact that the current study targets all South African river catchments, a specific dataset from a distinct Water Service Authority (WSA) was considered to ascertain the appropriateness of the proposed model. It would be a far-reaching and considerable effort to test the model against water quality data from all the Water Boards (WBs) in South Africa. On these grounds, water quality data from Umgeni Water Board (UWB) was deemed appropriate to establish the effectiveness of the developed water quality model. The selection of UWB does not devalue the purpose

Water Classification
In the interest of simplifying the interpretation of water quality index (WQI) values, mostly to accommodate non-technical individuals, an index categorization schema was established. The classification mechanism is based on an increasing scale index and the advantage of this system is that it is identical to a normal percentage hierarchy [3]; therefore, the public can easily relate to its function and interpretation. The UWQI model yields WQI values between zero and one hundred. Accordingly, the WQI scores are categorized using classes ranging from one to five, with "Class 1" representing water of the highest degree of purity with a possible maximum score of one hundred and, vice versa, "Class 5" denotes water quality of the poorest degree with index scores nearing or equal to zero. In order to close gaps identified in some of the existing classification scales [7], appropriate mathematical functions with logical linguistic descriptors, which include, but are not limited to, "greater than", "less than", and "equal to" were used to appraise WQI scores and we respectively assigned them to the corresponding category.

Background and Specific Considerations
A substantial increase in population and improper disposal of wastewater have a significant influence on diminishing water quality in rivers and other surface water reservoirs. As a consequence, routine water quality assessment and pollution control measures are necessary to preserve and restore the healthiness of surface water bodies [43]. On the same basis, this study attempts to put forward a practical and standardized tool that can be used for monitoring surface water quality across all South African river catchments. Despite the fact that the current study targets all South African river catchments, a specific dataset from a distinct Water Service Authority (WSA) was considered to ascertain the appropriateness of the proposed model. It would be a far-reaching and considerable effort to test the model against water quality data from all the Water Boards (WBs) in South Africa. On these grounds, water quality data from Umgeni Water Board (UWB) was deemed appropriate to establish the effectiveness of the developed water quality model. The selection of UWB does not devalue the purpose of the study, rather it is the beginning of a long-term undertaking to demonstrate that the developed model is indeed universal and applicable to most, if not all, South African river catchments.
Umgeni Water Board is a Water Service Authority responsible for the water and sanitation affairs of KwaZulu-Natal Province in South Africa [44,45]. UWB falls under the jurisdiction of Pongola-Mtamvuna Water Management Area (WMA) which has four primary drainage regions labelled T, U, V and W. Among the four regions, primary drainage basin U was considered for the current study. Further to this, only four secondary drainage regions were selected and these are Umgeni, Umdloti, Nungwane and Minto River catchments, which are identified by the Department of Water and Sanitation (DWS) as U20, U30, U70 and U80, respectively. Umgeni River catchment is the largest of the four; consequently, it is regarded as the primary study area, and henceforth it is considered more significant than the other three catchments.

Umgeni River Catchment
Umgeni River catchment is a sub-humid drainage basin located along the Indian Ocean coastline in KwaZulu-Natal Province in the Republic of South Africa [14,46,47]. Having a diversified land usage and multiple water supply systems, Umgeni basin is regarded as one of the most complex drainage regions in the country. The basin is subdivided into twelve quaternary drainage regions, also known as quaternary catchments (QCs). Seven of the QCs are situated in the upper most part of Umgeni basin, three are in the middle and two are in the lower part of the secondary drainage region [13,46]. Umgeni River catchment plays a major role in the economic development of the country; it serves South Africa's biggest trading port and the second largest province in terms of population and economic sizes [14]. This is why Umgeni is considered one of the most significant river catchments in South Africa.
Significantly, Umgeni basin addresses the water needs of the Durban-Pietermaritzburg business corridor and act as the primary source of water supplied to the Port of Durban, which is the biggest trading port in Africa and contributes significantly to the Gross Domestic Product (GDP) of South Africa [14,15]. Considering the social and economic activities in KwaZulu-Natal, the province is regarded to be a highly ecologically disturbed region [13], and this describes the motivation for the adoption of Umgeni catchment as the main study area. The current activities and projected developments in Umgeni River catchment have extraordinary effects on the national water resources and require a comprehensive water management monitoring model that focuses on protecting the water reserves. It is therefore important to develop a water quality index model that can be adopted to better understand the dynamics of water quality changes in Umgeni River catchment and South Africa as a whole. The model will provide institutional support in delineating water quality concerns across various river catchments and provide substantial information to water technocrats and decision-makers.
Umgeni River catchment has a surface area nearly 4432 km 2 , with Umgeni River being the major water channel of the drainage basin [13,15,19,48]. The 232-kilometre-long river originates from the Drakensberg mountains and flows eastwards towards the Indian Ocean, with four main cardinal tributaries, namely Lions, Karkloof, Impolweni and Umsunduzi Rivers [19,47]. Lions River is the most contributing tributary in the upstream of Midmar Dam and it serves as the transfer channel, conveying water resources from the adjacent Mooi River Catchment [13]. The basin land cover is characterized as heterogeneous, mostly consisting of urban areas, natural forest, commercial sugarcane plantations, small-scale and commercial agricultural farms and the Port City of Durban [13][14][15]46]. Notably, Umgeni River supports the livelihood of informal settlers residing along the river course. They rely on the river for various household activities, irrigation and livestock production [49]. The rainfall pattern of Umgeni basin is seasonal, with rains concentrated in the summer months (October to March). The amount of precipitation is highly variable, increasing from the western side to the eastern part of the river catchment. Highest rainfall occurs in coastal areas with a range of 1000 mm/y to 1500 mm/y [15,47]. The inland parts of Umgeni basin generally receive rainfall ranging from 800 mm/y to 1000 mm/y [15,46,50]. The average annual temperature ranges from 12 • C to 22 • C; leading to evaporation rates between 1567 mm/y and 1737 mm/y [13]. Four major dams are used to regulate and reserve the water resources in Umgeni drainage region, and these are Midmar, Albert Falls, Nagle and Inanda [13,16]. Midmar Dam supplies Pietermaritzburg and some portions of Durban, whereas Albert Falls, Nagle and Inanda Dams cater for the greater part of the Durban metropolitan area [19,46,47]. In addition to the four major dams, there is also Henley Dam, situated south of Midmar Dam along the Msunduzi River, a tributary of the Umgeni River. Apart from that, there are about 300 farm dams utilized for irrigating nearly 185 km 2 of commercial farms in the Umgeni catchment area [46].

Umdloti River Catchment
Umdloti catchment is situated north-east of Umgeni basin, adjacent to Nagle and Inanda Dams. The catchment has an estimated area of 597 km 2 with Umdloti River as the main watercourse of the basin [51]. The river source is found in the Noodberg area and stretches for nearly 88 km, flowing eastwards toward the Indian Ocean. The river estuary is approximately 25 km northeast of central Durban [48,52]. A considerable portion of the catchment is utilized for commercial farming, dominated by sugarcane and banana plantations with minimal vegetable and citrus farming. Apart from these, other establishments include residential areas, Verulam Town, game reserves, Hazelmere Dam and Hazelmere wastewater treatment plant [52]. In a similar manner to Umgeni basin, the catchment experiences summer rainfall, with mean annual precipitation ranging between 800 mm and 1125 mm. Temperatures vary from 9 • C in winter to 38 • C in summer [52]. Hazelmere Dam is the major water impoundment in Umdloti catchment [51]. The dam was established to service the domestic, industrial and agricultural needs of the Durban area, including the new Durban International Airport [48,52].

Nungwane River Catchment
Located southwest of Umgeni drainage region, Nungwane River catchment has mean annual precipitation of 938 mm/y and annual evaporation close to 1200 mm/y. The significant impoundment in the quaternary catchment is the Nungwane Dam, situated along the Nungwane River, which is a tributary of Lovu River [53]. The impoundment was built in 1977, with a catchment area of 58 km 2 , and raw water from Nungwane Dam is treated at Amazimtoti water treatment plant and then supplies eThekwini Municipality [53].

Umzinto/uMuziwezinto River Catchment
Umzinto River catchment, also known as uMuziwezinto River catchment, lies further south of Nungwane Dam. According to Umgeni Water [53], the river basin receives rainfall averaging 985 mm/y, with an evaporation rate of 1200 mm/y. In 1983, Umzinto Dam was constructed along Umzinto/uMuziwezinto River, with a catchment area of about 52 km 2 . Together with EJ Smith Dam, raw water from Umzinto Dam is treated at Umzinto water treatment plant (WTP) and distributed to Ugu District Municipality [53,54]. Both dams, EJ Smith and Umzinto, supply raw water to the operation of Umzinto WTP [53,54].

Sampling Locations
Umgeni Water Board (UWB) established water sampling stations to enhance water quality monitoring and the stations are strategically positioned to provide a holistic understanding of water affairs within the service area of KwaZulu-Natal Province. Instead of establishing new research-based sampling stations, the current studies utilized water quality data collected by UWB. At least one or more stations were considered for each of the four drainage basins discussed in the preceding sections. Further details of the selected sampling stations are included in Table 1 and Figure 2.  The economic importance of the Umgeni Basin, the uniqueness of its inter-basin arrangements, the magnitude of the transfer schemes involved and the extensive water demand require comprehensive water resource management. All of these factors distinctively motivated the identification and selection of the Umgeni River catchment as the main study area. Beyond that, three additional catchments were incorporated into the study in order to further examine the model and complement the objective of developing a universally acceptable water quality model.

Research Dataset
Regular water quality sampling and analysis is a costly and demanding task, hence acquiring large volumes of water quality data is often a challenge and requires significant amount of financial resources [55,56]. Given that the authors could not gather their own samples, water quality data from Umgeni Water Board (UWB) were used to attest the functionality of the model. The dataset from UWB is for six sampling stations that fall under the jurisdiction of four different catchment areas. It contains 416 monthly samples for a period of four years, ranging from June 2014 to July 2018 and the dataset statistics The economic importance of the Umgeni Basin, the uniqueness of its inter-basin arrangements, the magnitude of the transfer schemes involved and the extensive water demand require comprehensive water resource management. All of these factors distinctively motivated the identification and selection of the Umgeni River catchment as the main study area. Beyond that, three additional catchments were incorporated into the study in order to further examine the model and complement the objective of developing a universally acceptable water quality model.

Research Dataset
Regular water quality sampling and analysis is a costly and demanding task, hence acquiring large volumes of water quality data is often a challenge and requires significant amount of financial resources [55,56]. Given that the authors could not gather their own samples, water quality data from Umgeni Water Board (UWB) were used to attest the functionality of the model. The dataset from UWB is for six sampling stations that fall under the jurisdiction of four different catchment areas. It contains 416 monthly samples for a period of four years, ranging from June 2014 to July 2018 and the dataset statistics are presented in Table 2. Hardness (CaCO 3 ) tests were recorded on a quarterly basis and, where possible, using the measured values of calcium (Ca) and magnesium (Mg); alternatively, the estimation of missing CaCO 3 values was achieved through Equation (4).
where all parameter concentration levels are in mg/L. The method is common practice and is prescribed in the literature [3,39,40,57,58]. Source: Umgeni Water Board. Notes: 1 sampling station identification number, 2 minimum measured water quality values, 3 mean/average of measured water quality values, 4 maximum measured water quality values, 5 standard deviation, 6 coefficient of variation as a percentage, and 7 water quality variables measured in mg/L, except for chlorophyll-a (µg/L), electrical conductivity (µS/m), pondus Hydrogenium (unitless), and turbidity (NTU).
The samples obtained from UWB were adequate and contributed significantly towards the success of the current study. Umgeni water quality data were considered based on availability; other than being a priority and limiting the number of WBs used for testing, the models do not devalue the significance of the study. The rationale used in developing the universal water quality index is completely independent of the dataset used for testing the functionality of the model; nevertheless, as an ongoing project and in support of the current studies, it is recommended that additional data from other WBs, if not all, be considered and documented separately.

Water Quality Variables and their Relative Weightage Coefficients
With the aid of the Delphi method and the existing literature, a fixed set of thirteen physico-chemical parameters were found to be adequate and appropriate to analyze and compare water quality status among different sites. A fixed system requires enormous care, attention, experience and proficiency to ensure that the most significant variables are incorporated in the WQI. Expertise is required to delineate what might be regarded as too few or too many variables, the ability to optimize the ideal number of parameters necessary or just enough to calculate a meaningful water quality index value [7]. Therefore, the study involved expert opinions through a group of selected professionals and extracts from similar studies. Accordingly, the most appropriate variables considered for inclusion in the UWQI are ammonia (NH 3 ), calcium (Ca), chloride (Cl), chlorophyll-a (Chl-a), electrical conductivity (EC), fluoride (F), hardness (CaCO 3 ), magnesium (Mg), manganese (Mn), nitrate (NO 3 ), pondus Hydrogenium (pH), sulphate (SO 4 ) and turbidity (Turb).
Considering that, in the current study, water quality parameters are viewed to have different influences on the overall classification of water, some variables are considered greater than others; therefore, weights were established to appropriately reflect the diversity of each parameter. The comparative scale used is biased towards the level of influence and significance towards the overall index value [9,38]. As represented in Table 3, two sets of parameter significance ratings obtained from the participation-based Delphi method and extraction from the existing literature are used to derive the weight ratings based on a common scale of influence, ranging from one (lowly rated) to five (highly rated). Given the Delphi significance rating (c i ) and the literature significance rating (d i ), the parameter weight rating (b i ) is given by (c i + d i )/2, whereas the final weight coefficient (w i ) is transformed by dividing the relevant parameter weight rating (b i ) by the sum of all weight ratings (w i = b i / b i ) and the index weight coefficients are represented as a decimal number with a sum equal to one (w 1 + w 2 + w 3 + . . . + w n = 1). In principle, this theory governs the model from computing index values in excess of one hundred percent, otherwise the aggregation process will be compromised and will jeopardize the scientific steadiness of the model [3]. Notes: parameters are listed alphabetically. Using final weight coefficients in Table 3, the following order of importance is achieved:

Formation of Parameter Sub-Index Rating Curves and Sub-Index Functions
Given the fact that the identified model input variables are assessed using different units of measurement, sub-indices were then developed to transform the measurement units into a common unitless scale. Moreover, the adopted indexing model can only aggregate parameters with a common scale, which became more necessary in order to harmonize the parameter values using a standardized non-dimensional scale. In relation to the permissible water quality parameter concentrations prescribed by DWAF [39][40][41], fixed key points of the rating curves were established (see Table 4) and converged with straight-line graphs. Thereafter, the linear equations associated to the straight-line graphs were collectively transformed into linear sub-index functions. The advantage with this technique is that sub-index functions are able to interpolate the limits in between water classification categories using the linear regression method. Examples of the final sub-index curves are included in Figure 3, whereas examples of sub-index functions are represented mathematically from Equation (5) to Equation (7).
where SI a,b, . . . , c are sub-index functions for (a) ammonia, (b) calcium, (c) chloride; x a,b, . . . , c are the observed water quality reading of the respective water quality parameter. Due to the nature of the article, only sub-index rating curves and mathematical functions for NH 3 , Ca and Cl are presented herein; the rest are documented elsewhere.
where SIa,b, … ,c are sub-index functions for (a) ammonia, (b) calcium, (c) chloride; xa,b, … ,c are the observed water quality reading of the respective water quality parameter. Due to the nature of the article, only sub-index rating curves and mathematical functions for NH3, Ca and Cl are presented herein; the rest are documented elsewhere.

Weighted Indexing Model (UWQI)
The mathematical structure and application of indexing models is normally governed by the degree of accuracy perceived and the type of weightage coefficients, which might be equally or unequally defined. Various aggregation methods exist, and each technique has its own formidable challenges; thus, the index developer has to decisively select the most appropriate and relevant indexing model, preferably one with fewer complications that might adversely influence the final index value. Otherwise, defining the best and absolute aggregation model is close to impossible. Since there is no supreme and favorable technique of formulating water quality indices (WQIs), various aggregation methods were tried and tested. The modified weighted sum (additive) method was found to be the most appropriate for the

Weighted Indexing Model (UWQI)
The mathematical structure and application of indexing models is normally governed by the degree of accuracy perceived and the type of weightage coefficients, which might be equally or unequally defined. Various aggregation methods exist, and each technique has its own formidable challenges; thus, the index developer has to decisively select the most appropriate and relevant indexing model, preferably one with fewer complications that might adversely influence the final index value. Otherwise, defining the best and absolute aggregation model is close to impossible. Since there is no supreme and favorable technique of formulating water quality indices (WQIs), various aggregation methods were tried and tested. The modified weighted sum (additive) method was found to be the most appropriate for the development of a universal water quality index for monitoring South African watersheds. The modified weighted sum (additive) method is represented in Equation (8) [26][27][28].
Upon conducting a scenario-based analysis, the modified weighted sum equation has been further improved to align it with local conditions and the developed final universal water quality (UWQI) model is presented in Equation (9): where UWQI is the aggregated index value ranging from zero to hundred, with zero representing water of poor quality and hundred denoting water of the highest quality; s i is the sub-index value of the ith water quality parameter obtained from the sub-index linear functions and the values range from zero to hundred, similar to WQI values; w i is the weight coefficient value for the ith parameter represented as decimal number and the sum of all coefficients is one, (w 1 + w 2 + w 3 + . . . + w n = 1); n is the total number of sub-indices-in this case, n = 13.

Scenario-Based Model Validation Analysis
Scenario-based analysis helps to identify potential data-processing gaps, which, in turn, enlighten us as to the necessary precautions that are imperative in order to minimize the impact, or perhaps eliminate the problem completely. To determine such precautions, ideal sets of predictive variables have been established under a variety of scenarios to calculate specific water quality variables. Considering increments of five scores, nine probable scenarios have been examined to demonstrate the model's ability to predict scores of all ranges, from class one (excellent) to class five (the worst). The nine forecasts use three-level grading, comprised of (i) worst-case scenario, 0 ≤ Index ≤ 10, (ii) base-case scenario, 45 ≤ Index ≤ 55 and, lastly, (iii) best-case scenario, 90 ≤ Index ≤ 100. Purposefully, the groupings provided a complete change of circumstances with each scenario, thereby widening the range of analysis and including a considerable array of possibilities. With reference to permissible concentration limits and developed linear sub-index functions, definite assumptions about all nine cases have been carefully considered. Accordingly, parameter values corresponding to each scenario have been established and applied to perform the analysis. The results of the scenario-based analysis are presented in Table 5 and Figure 4. Table 5. Comparison of modified weighted arithmetic water quality index and the developed universal water quality index using scenario-based analysis to establish the functionality and predictive capacity of the models.  While we do not wish to devalue the efforts by House [26][27][28], the modified weighted arithmetic model could not sufficiently satisfy the expected analytical results. Although the predictive pattern is recommended, there is a significant lag between the calculated results and the ideal case, especially with the base-case scenarios (45  Index  55). Henceforth, the model was further improved to suit our local conditions. In view of the analysis results, it is evident that the proposed UWQI is robust and technically  While we do not wish to devalue the efforts by House [26][27][28], the modified weighted arithmetic model could not sufficiently satisfy the expected analytical results. Although the predictive pattern is recommended, there is a significant lag between the calculated results and the ideal case, especially with the base-case scenarios (45 ≤ Index ≤ 55). Henceforth, the model was further improved to suit our local conditions. In view of the analysis results, it is evident that the proposed UWQI is robust and technically stable. The degree of variation from the ideal values is negligible; the prediction pattern followed the ideal graph with corresponding values in terms of both WQI scores and classification. This, therefore, pronounces the competence of the UWQI to be used as an evaluation tool for monitoring South African river catchments.

Sample
Good water quality index scores nearing or equal to one hundred are achieved when the surface water shows the virtual absence of threats or impairments, conditions very close to pristine (natural) levels. On the other hand, index values close or equal to zero are recorded when almost all the water quality parameters depart from desirable concentration levels. Ideal parameter concentrations for each possible level of contamination are documented in Tables 4 and 6.

Evaluation of Water Quality
Umgeni water quality data have been evaluated using the proposed universal water quality index (UWQI) model documented in Equation (9). Based on the UWQI, Tables 6 and 7 and Figure 5 indicate spatial and temporal water quality variations among the six sampling sites. In order to demonstrate further the ability of the suggested UWQI, Table 6 also includes WQI scores calculated using the ideal values derived from the key points of the rating curves.    High levels of NO3 are recorded during the summer periods and, considering the socio-economic developments surrounding the sampling stations (Figure 2), the source of contamination might be anthropogenic activities, especially wastewater discharge, among others. NO3 is a naturally occurring ion [59,60] that is widespread and is regarded as the most significant contaminant in water [61,62]. Nitrate itself is a low-toxicity compound, but when endogenously converted to nitrite (NO2), it becomes toxic to The results show that water quality in the region can be categorized as "acceptable water quality", with the lowest WQI score of 75.99 (class two) recorded at station 5 (Umzinto Dam). In this case, turbidity and Mn are the main contributors to the deterioration of the water quality, with concentration levels of 13.20 NTU and 1.05 mg/L, respectively. Sampling station 2 recorded the highest surface water quality with an index of 95.15 (class one) during the summer of 2017. NO 3 is the principal pollutant factor responsible for the minimum WQI scores for stations 2, 3, 4 and 6, with NO 3 concentrations of 0.99, 1.31, 4.50 and 1.77 mg/L, respectively (see Table 6).
High levels of NO 3 are recorded during the summer periods and, considering the socio-economic developments surrounding the sampling stations (Figure 2), the source of contamination might be anthropogenic activities, especially wastewater discharge, among others. NO 3 is a naturally occurring ion [59,60] that is widespread and is regarded as the most significant contaminant in water [61,62]. Nitrate itself is a low-toxicity compound, but when endogenously converted to nitrite (NO 2 ), it becomes toxic to human health and the aquatic environment [59,60], thus exemplifying the need for regular water quality monitoring to identify water quality trends over time and space [63].
High levels of turbidity are evident during the summer seasons at stations 1, 2, 5 and 6, with corresponding values of 97.20, 66.70, 13.20 and 13.30 NTU. Together with NO 3 , turbidity contributes significantly towards the deterioration of water quality among these sites. Sources of turbidity are diverse and include, but are not limited to, reservoir drawdown flushing, algal blooms (eutrophication), wastewater discharge, industrial effluent, exceptional rainfall events, soil erosion and decomposition of organic matter [64,65]. Chl-a concentrations at stations 2 and 3 exceed targeted water quality levels in summer, with values of 20.49 and 19.50 µg/L, respectively. Soluble nutrients, especially phosphorus and nitrogen, are the key determinants promoting algae blooms (eutrophication), which contribute significantly towards increased levels of chlorophyll-a [66,67]. Such enriching nutrients often originate from anthropogenic activities, which include wastewater discharge and fertilizer runoff [3,68,69].
Narrow variations in WQI are observed for stations 1 (77.98-88.08) and 4 (83. 30-94.92). The two stations are located upstream of the catchment and the rest of the sampling sites are situated downstream of the drainage region, towards Durban-Pietermaritzburg business corridor. WQI results indicate that surface water quality varies more with the increase in socio-economic activities along the river's watercourse, with station 2 having the largest variation (77.87-95. 15).
Testing the model with data from various river catchments promotes the objective of establishing a universal water quality index suitable for use across the catchment areas in South Africa. Noticeably, the UWQI model responded steadily to the highs and lows of each water quality parameter value, with the index output graphs showing the variations. This advocates the readiness of the UWQI to interpret water quality data and provide a simple non-dimensional score that is justifiable and repeatable. Such success fulfils the objective of developing a universal WQI and, more importantly, presents a "yardstick" that can be applied in most, if not all, of the distinct watersheds in South Africa. This accomplishment is a critical milestone not only to the authors, but also to most of the stakeholders directly or indirectly involved in water quality science.

Index Categorisation Schema
A five-class WQI categorization schema has been adopted for this study ( Table 8). The schema is an increasing scale that is identical to normal percentage hierarchy, offering a better understanding of water classification scales. The ranking mechanism is similar to WQI classifications suggested for the Boyacioglu Index (Turkey) and Vaal WQI (South Africa) [3,5]. Table 8. Index score classification for the universal water quality index for south african river catchments.

ID Water Quality Classification
Description of Rank and Classification Index Score In a similar manner to the methods used by and Sutadian et al. [20], Abrahão et al. [70], Rabee et al. [71] and Rubio-Arias et al. [72], appropriate mathematical functions with logical linguistic descriptors such as less than, equal to and greater than have been assigned to each categorization class. In this way, the categorization schema can accommodate all possible index scores regardless of the decimal value. More importantly, the established categorization schema aids in closing gaps identified in the literature and presents a progressive approach that will contribute significantly towards water quality index development. Such an academic contribution reflects on the efficiency of the model and can be attributed to the success of the current study.

Conclusion and Future Directions
Over four hundred water quality samples from six sampling stations located in four different river catchments are evaluated using UWQI, and Table 7, together with Figure 5, provides a summary of the trend analysis. The spatial and temporal changes in water quality for Umgeni Water Board are evident over a period of four years, with a varying sequence comprising of index scores as high as 95.15 (class one), with an average of 87.78 (class two) and the lowest score being 75.99 (class two) across the six sites. The best surface water quality was recorded at station 2 during the summer period of 2017, whereas the lowest water quality was recorded at station 5 during the month of August 2014. The main pollution contributors are NO 3 (station 2, 3, 4 and 6), turbidity (station 1, 2, 5 and 6), Chl-a (station 2 and 3) and, lastly, Mn at station 5. The sources of pollution may be associated with anthropogenic activities, considering the socio-economic developments surrounding the affected sampling stations. Otherwise, the rest of the water quality parameters are virtually within permissible levels. There is the need for regular water quality appraisal to monitor concentration levels against pollution control regulations and record the variation in trends, especially for sampling stations located within the Durban-Pietermaritzburg business corridor. The application of UWQI can perform sustainable water resource functions for river basin management.
The study opens a path for unified WQIs to be considered in South Africa. As this is the first attempt to demonstrate the use of nationally applicable indices, it is highly expected that the study will have an impact on methods of developing future water quality indices, contributing to our understanding of index models and supplementing our knowledge of water quality science. It is important to conduct research into unified WQIs formed based on multivariate statistical approaches. Further research is required to better understand the performance of objective methods on nationally applicable indices and address the effects of subjectivity on traditional methods of establishing WQIs. As an ongoing study, additional data from other river catchments should be considered and evaluated using the suggested UWQI. This will further demonstrate the universality of the model and perhaps provide guidance on necessary modification requirements. Nevertheless, an initial step towards the ultimate goal has been achieved, which is the development of a universal water quality index (UWQI).