**Hydrogeochemical Assessment of Groundwater Quality of Mokopane Area, Limpopo, South Africa Using Statistical Approach**

#### **Mmasabata Dolly Molekoa 1, Ram Avtar 1,2,\*, Pankaj Kumar 3, Huynh Vuong Thu Minh 1 and Tonni Agustiono Kurniawan 4**

Received: 14 July 2019; Accepted: 7 September 2019; Published: 11 September 2019

**Abstract:** Despite being a finite resource, both the quality and quantity of groundwater are under tremendous pressure due to rapid global changes, viz. population growth, land-use/land-cover changes (LULC), and climate change. The 6th Sustainable Development Goal (SDG) aims to "Ensure availability and sustainable management of water and sanitation for all". One of the most significant dimensions of the SDG agenda is the emphasis on data and governance. However, the lack of good governance coupled with good observed data cannot ensure the achievement of SDG6. Therefore, this study strives to evaluate water quality status and hydrochemical processes governing it in the data-scarce Mokopane area of South Africa. Groundwater is the main source of fresh water supply for domestic usage, intensive agriculture, and mining activities in Mokopane. In this study, hydrogeochemical analysis of groundwater samples was employed to calculate the water quality index (WQI) and evaluate factors governing water quality evolution in the study area. Statistical and spatial analysis techniques were carried out to divide sampling sites into clusters and delineate principal factors responsible for determining water quality of the sampled groundwater. Results suggest that most of the physico-chemical parameters are within permissible limits for drinking water set by the World Health Organization (WHO), except for high fluoride in some samples. Na-HCO3 is the most abundant water type followed by Mg-HCO3, which indicates dominance of Na+, Mg2+, and HCO3 ±. Rock-water interaction is the prime factor responsible for fluoride enrichment in water. The alkaline nature of groundwater favors the release of exchangeable F− from minerals like muscovite. The WQI suggests that 80% of water samples fall into the good and excellent categories. Poor management of untreated domestic sewage and agricultural runoff is a main factor for the bad/very bad categories of water samples. As the area lacks any credible scientific/government work to report water quality and its management aspects, the findings of this study will definitely help both scientific communities and policy makers to do what is needed for sustainable water resource management in a timely manner.

**Keywords:** groundwater; geospatial analysis; hydrogeochemical assessment; multivariate statistical analysis; water quality index

#### **1. Introduction**

Groundwater is one of the most important finite natural resources for fresh water on Earth [1–3]. Its importance has increased as people depend on it for agriculture and industry, with domestic

use exponentially increasing [4–8]. Groundwater exploration and climate change in recent years have led to accelerated depletion of groundwater levels as well as deterioration of groundwater quality, especially in shallow aquifers [9]. Groundwater quality is as important as its quantity, and therefore it is essential to consider in depth the holistic approach towards its management [10]. According to Subramani et al. [11], ground water quality is relatively more important than surface water as it is less susceptible to contamination/pollution and is suitable for different purposes. The quality or physico-chemical characteristics of groundwater is highly dependent on anthropogenic activities like agriculture, mining, urban settlements, etc.; as well as natural processes like rock–water interaction, geological formations, and aerobic/anaerobic conditions of the aquifers [6,12–14]. It is reported that spatio-temporal variation in groundwater quality is also predominantly governed by the hydrogeochemical process occurring at the recharge area and groundwater–surface water interaction [15]. Also, climate change and frequent extreme weather conditions exaggerate the severity of water scarcity [16–18]. Groundwater quality assessment can indicate its suitability for different purposes like drinking and irrigation. The concentration of different chemical parameters beyond their permissible limit leads to serious health complications such as fluorosis, blue baby syndrome, and arsenicosis from fluoride, nitrate, and arsenic, respectively [19,20]. Weathering and leaching are the biggest drivers for high fluoride concentration in groundwater [21,22], whereas nitrate enrichment is mainly attributed to agricultural runoff with high usage of fertilizers and leaching from untreated sewerage systems [7,23]. Therefore, diligent monitoring is a vital part of groundwater quality management, and thus, for maintenance of healthy aquatic ecosystems and human health [24–26].

Usage of different tools and techniques (e.g., numerical modelling, geographical information systems (GIS)) by both scientific communities and decision makers for groundwater quality assessment and management has rapidly increased over the past few decades [27–30]. A GIS is one of the proven vital tools in studies of groundwater for different purposes (spatial distribution of contaminants, estimating the effect of land-use/land-cover (LULC) changes on water quality, overlay analysis for estimating groundwater potential, etc.) throughout the world [31–35]. In this regard, the use of the water quality index (WQI) can provide simplistic information from complex water quality data to the public and managers. WQI can be used for assessment of surface as well as groundwater quality. Moreover, various WQIs that were developed for irrigation and drinking purposes can help to minimize the adverse effects of water use [6,36]. With the above background, this study focuses on groundwater quality assessment in Mokopane in South Africa using laboratory analysis of water quality parameters well-supported by interpretation through different statistical and GIS tools. This integrated assessment has helped us to understand well the processes governing hydro-geochemical changes in spatial patterns, and hence, the direction for effective management of groundwater resources. The rationale behind choosing this study area is that most of the households in this area rely on groundwater as the sole source of fresh water. However, there are various local reports and complaints relating to groundwater quality deterioration credited to poor sanitation, mining activities, and extensive agriculture in the area [37]. Also, the majority of Mokopane communities rely on pit toilets which have the potential to adversely affect the groundwater quality. Despite the above facts, no detailed scientific study has been conducted to analyze the current status of groundwater quality in the Mokopane. Therefore, this study is an integrated approach (chemical, statistical, spatial analysis) to analyze groundwater quality and processes which are the most likely source of pollutants.

#### **2. Study Area**

#### *2.1. Site Description*

The study area is located in the northern part of Limpopo province of South Africa. Limpopo borders Namibia, Botswana, Zimbabwe, and Mozambique. Mokopane is a town located in the Waterberg district of Limpopo province (Figure 1) with a population of 328,905 in 2016 [38]. It is located at 24◦11 2 S latitude and 29◦0 46 E longitude. Figure 2 shows the climatic condition of the study area which falls within steppe climate with mean annual maximum and minimum temperatures of 23.4 ◦C and 13 ◦C respectively [39]. The average annual rainfall is 490mm which normally occurs in November, December, and January. Mokopane has less rainfall during the winter season. July and August are normally dry due to the winter climate.

**Figure 1.** Location of the study area sampling points and mining sites.

**Figure 2.** Climate data of average rainfall and temperature in Mokopane.

Economically, Mokopane is known for its intensive agricultural activities until the commencement of Mogalakwena Platinum mine and Ivanplats mine. The mines provide employment to a large portion of the Mokopane population. The Mogalakwena mine is located about 20 km northwest of Mokopane and Ivanplats mine inside Mokopane area. Geologically the study area has deposit of the Bushveld complex. It contains the world's largest reserves of platinum-group minerals (PGMs). The Bushveld Complex comprises the Rustenburg Layered suite, the Lebowa Granites, and the Rooiberg Felsics, that are overlain by the Karoo sediments [40,41]. The development of the mining sector in the study area triggered rapid development in informal settlements and, hence, poor domestic sewerage treatment in the area. Poor settlement and sewerage treatment can cause high risks to human health in the near future if these issues are not taken into consideration.

#### *2.2. Land-Use*/*Land-Cover Classification*

The LULC map of the study area was prepared using the Landsat 8 satellite image acquired on 16th March 2018. The area was classified into five classes using maximum likelihood classification algorithm in the ENVI 5.2 software (Figure 3). Five classes of LULC map were classified, viz. agriculture land, settlement, vegetation, bare land, and water bodies. The selection of classification classes was based on the field work and other existing maps. Figure 3 shows the LULC maps and locations of the water samples and Ivaplats mining area. Results show that most of the study area was covered with agriculture land (53%) followed by settlement (32%), vegetation (13%), bare land (2%), and water bodies (1%). Recently, the settlement area has increased in the study area because of the increase in mining activities resulting in high employment opportunities.

**Figure 3.** Land-use/land-cover (LULC) map based on Landsat 2018 by using maximum likelihood classification algorithm.

#### **3. Methodology**

#### *3.1. Flow Chart of Research Methodology*

Figure 4 shows the flowchart of methodology adopted in this study. Analyzed water quality data was further processed using piper diagram, speciation modeling, groundwater quality index (GWQI) development, statistical analysis, and was well supported by spatial image development.

**Figure 4.** Flowchart of the methodology employed in this study for the assessment of groundwater quality.

#### *3.2. Sample Collections and Hydrochemical Analysis*

A total of 12 groundwater samples were collected on 12th June 2018 and analyzed for 21 physico-chemical parameters at Capricorn Veterinary Laboratories. Sampling coordinates were recorded using the Global Positioning System (GPS III, Garmin) as shown in Figure 1. The average well depth of groundwater sampling wells is 103.4 m below ground level (mbgl). Figure 1 shows the locations of the three major mining areas (Ivanplats mine, Mogakawena platinum mine, and Anglo American platinum mine). The inline flow cell that ensures the exclusion of minimized fluctuations and atmospheric contamination was used on-site in order to measure the electrical conductivity (EC) and pH. Precision of > 5% achieved with the Model Beverly, MA, 01.915, known as the transportable Orion Thermo Water analyzing kit, was used for all kinds of on-site measurements. Collected groundwater samples were stored in thoroughly rinsed polyethylene bottles. In order to stabilize trace metals (pH~2), the collected samples for major cation analysis were acidified by 1% HNO3; while samples collected for nitrate were acidified with H3BO3. The concentration of HCO3 − was analyzed by acid titration (using Metrohm Multi-Dosimat); while other anions Cl−, NO3 −, SO4 2−, and PO4 3− were analyzed by DIONEX ICS-90 ion chromatograph with an error percentage of less than 2% using duplicates. Inductively coupled optical emission spectrometry (ICP-OES) was used to evaluate major cations using duplicates, with a precision of less than 2%. One replicate was used to check the accuracy of the instruments after the analysis of every five samples for each instrument. The analytical precision was checked by the normalized inorganic charge balance (NICB) for major ions. The fractional difference between the total cations and anions is defined as [(Tz+ − Tz−)/(Tz+ + Tz−)]. Here, Tz+ and Tz− represent the total milli-equivalent of cations and anions, respectively. The observed charge balance supports the quality of the data points, which is better than ±8%, and this charge imbalance was generally in favour of positive charge.

#### *3.3. Statistical Analysis*

Different multivariate statistical tools, viz. cluster analysis (CA) using agglomerative hierarchical clustering (AHC) and principal component analysis (PCA) were applied on the analyzed dataset to find out the interrelationship between the groundwater quality variables and to determine the major factors influencing the groundwater quality [7,12,14,42–46]. In addition, the Piper diagram was also used to classify water samples into water types [46,47].

#### 3.3.1. Saturation Index for Estimating Groundwater Mineralization

The saturation index (SI) was calculated in order to understand the origin of groundwater mineralization in the study area. According to Garrels and Christ [48], mineral saturation states with respect to selected minerals were computed as shown by Equation (1).

$$\text{SI} = \log\_{10}(\mathcal{K}\_{\text{IAP}} / \mathcal{K}\_{\text{SP}}) \tag{1}$$

where SI is the saturation index, *KIAP* is the ion activity product of a particular solid phase, and *KSP* is the solubility product of the phase. Based on the SI values, the saturation states were recognized as saturated (equilibrium; SI = 0), unsaturated (dissolution; SI < 0), and oversaturated (precipitation; SI > 0).

#### 3.3.2. Water Quality Index

This study has used the water quality index (WQI), which has been considered as one of the most reliable tools to classify water contamination levels for both ground and surface water [49]. WQI is based on a scale of 0–300 whereby lower values indicate good water quality, whereas higher values are an indication of contaminated water. Based on the WQI values of < 50, 50–100, 100–200, 200–300, and > 300, water samples were categorized as excellent, good, poor, very poor, and unsuitable for drinking, respectively. The following steps were taken in order to calculate WQI:

1. Calculating relative weight: this was calculated using Equation (2).

$$\mathcal{W}\_{\dot{i}} = \begin{array}{c} w\_{\dot{i}}\\ \frac{\sum\_{i}^{n} w\_{\dot{i}}}{\sum\_{i}^{n}} \end{array} \tag{2}$$

where *Wi* represents the relative weight of each parameter sampled, *wi* represents the weight of each parameter, and *n* represents the total number of parameters. All calculated values of *Wi* are shown in Table 1.

2. Calculating Q value: this was calculated using Equation (3).

$$Q\_i = \frac{C\_i \times 100}{S\_i} \tag{3}$$

where *Qi* is quality rating, *Ci* is concentration of each parameter (mg/L), and *Si* is derived from the WHO water quality standard.

3. Finally, WQI was calculated using Equation (4).

$$WQI = \sum W\_i \times Q\_i \tag{4}$$

Water quality classification and ratings were designed by Singh et al. [50]. In their paper, they have shown that excellent ≤ 50, good water = 50–100, poor water = 100–200, very poor = 200–300, and water unsuitable for drinking purpose ≥ 300. Table 1 shows the relative weights of chemical parameters.

#### 3.3.3. Spatial Analyses

In general, asymmetrical distribution of water quality occurred due to water samples being collected from many sites and at different time [51]. To reduce this asymmetry, statistical analysis is one of the methods often used to evaluate water quality variation in many studies in Vietnam, India, and Japan [7,12,52]. The various multiple statistical tools, namely, cluster analysis (CA), PCA, factor analysis, and the Piper diagram were useful to differentiate changes in water quality in a spatial pattern [7,12,52].

**Table 1.** The relative weight of chemical parameters.

#### *3.4. Mapping Spatial Distribution of Groundwater Quality*

Inverse distance weighting (IDW) interpolation is a geospatial method to see the spatial distribution of groundwater quality parameters. This interpolation method was developed by Shepard, Bartier, and Keller [53,54]. Although the IDW method contains a "bull's eye" which has an effect on concentric regions of the same value near the known points [55], IDW is based on the distance-weighting exponent which can precisely control the influence of the distance of the points [54]. The IDW interpolation method was used to see the spatial distribution of groundwater quality in Mokopane using ArcGIS 10.3 software.

#### **4. Results and Discussion**

#### *4.1. General Groundwater Chemistry*

A statistical summary of the groundwater quality is presented in Table 2. The pH values of the groundwater samples varied from 7.31 to 9.24 with an average value of 8.30, suggesting the alkaline nature of groundwater. The alkaline nature of groundwater also suggests enhanced dissolution due to the interaction between soil and rainwater [56]. In addition, the alkaline aquifer condition depicts the influence of an arid/semiarid climate with relatively low rainfall and high evaporation. The electrical conductivity of groundwater varied from 728.0–1895.0 μS/cm with an average value of 1022.17 μS/cm indicating high ionic activity in the area. This also suggests soil mineralization or salt deposition due to climatic conditions as supported by high pH. Looking into the analyzed water quality data for major ions, anionic abundance was in the order of HCO3 − > Cl− > SO4 2− > CO3 2−> NO3 − > F−> PO4 3−, while cationic abundance was found in the order of Na+ > Mg2+ > Ca2+ > K+. The highest average concentration of HCO3 − among anions is due to the presence of carbonaceous sandstones in the aquifers and weathering of carbonate minerals in rain water followed by subsequent precipitation of HCO3 − along with other cations. Higher Cl− and SO4 2− concentration also categorized as secondary salt mainly came from the leaching of sewage effluent, especially due to poor sewage management in the informal settlements. In addition, higher concentrations of SO4 2− can also be attributed to the leaching of organic matter and agricultural runoff carrying unutilized SO4 2−. Concentration of NO3 − and PO4 3− is not of concern as it is well below the permissible WHO limits for all groundwater samples. Concentration of F− varied from 0.24 to 3.39 mg/L. Sporadic fluoride contamination, i.e., >1.5 mg/L permissible limit set by WHO, for some of the groundwater samples is a matter of concern and puts the local people at moderate risk of fluorosis (tooth decay particularly in children and bone deformation) if consumed for longer periods of time and not given timely attention for its remediation. For cations, Na+ > Mg2+ > Ca2+ > K+ and the average milli-equivalent ratio of Na+ + K+/Mg2+ + Ca2+ was found to be 1.08, indicating the prominence of both alumina-silicate and carbonaceous weathering in the study area.

**Table 2.** Statistical summary of groundwater quality in the study area.

TH—Total hardness, EC—Electric Conductivity, TDS—Total dissolved solid, SS—suspended solid, mbgl—meter below ground level.

The concentrations of most of the parameters were found to be well within the permissible limits set by the WHO (2011) except for F− and Na+ for few samples. High concentration of Na+ in groundwater samples causes only scale formation in the pipelines due to high hardness of the water but does not cause any health effect to the consumers. Potential sources of high Na+ are alumino-silicate weathering like schist, granitic rocks containing Na+, and cation-exchange Na+/K+ with Ca2+/Mg2+. On the other hand, as high concentration of F− in the groundwater (i.e., >1.5 mg/L) causes fluorosis (dental, skeletal), further investigation was made about it. It was found that 75% of the water samples contain F− concentration less than 1.5 mg/L, hence the water is safe to drink or use for cooking. The rest of the water samples were categorized under low risk (17% with F− concentration between 1.5 to 2.9 mg/L) and high risk (8% with F− concentration between 3.0–5.0 mg/L).

The fluoride ion (F−) normally remains non-reactive for groundwater with pH ≤ 7.0 because of its low solubility and remains absorbed in clay or exchangeable materials. However, once the pH of groundwater goes beyond 7.0 (i.e., alkaline medium), the hydroxyl group (OH−) group replaces the exchangeable fluoride from silicate minerals like biotite, muscovite etc. resulting in fluoride enrichment. The whole chemical process can be explained for muscovite using Equation (5).

$$\mathrm{KAl\_2[AlSi\_3O\_{10}]F\_2 + 2OH^- \to KAl\_2[AlSi\_3O\_{10}][OH\_2] + 2F^-} \tag{5}$$

Another possible chemical reaction, where hydrolysis of alumina-silicate minerals produces bicarbonate ions, which can enhance fluoride dissolution, is calculated using Equations (6) and (7).

$$\text{NaAlSi}\_3\text{O}\_8 + \text{H}\_2\text{CO}\_3 + \text{H}\_2\text{O} \rightarrow \text{Na}^+ + \text{HCO}\_3^- + 2\text{H}\_4\text{SiO}\_4 + 2\text{Al}\_2\text{Si}\_2\text{O}\_5(\text{OH})\_4 \tag{6}$$

$$\rm CaF\_2 + 2HCO\_3^- \rightarrow CaCO\_3 + 2F^- + H\_2O + CO\_2 \tag{7}$$

IDW interpolation was used to show the spatial distribution of SO42− and F− as shown in Figure 5. The high concentration of F− was observed in boreholes 2, 3, and 7, where it exceeded the permissible limit of 1.5 mg/L. On the other hand, high SO42− concentration was observed in boreholes 8, 9, and 11. To further support its origin, speciation modelling for minerals related to fluoride was conducted. This is discussed later in the manuscript.

To further strengthen our findings, the depth profile of F− and SO4 2− is plotted and shown in Figure 6. It is found that F− concentration has a significant association with groundwater well depth. Its concentration is generally higher for groundwater samples from higher well depths. This supports the ideas of geogenic origin as presented above. On the other hand, the higher concentration of SO4 2− is found in groundwater samples with shallow well depths. As anthropogenic activities are also potential sources for SO4 2−, it can be expected that leaching of untreated and unmanaged domestic sewage will contaminate shallow aquifers first, thus firmly supporting the findings.

**Figure 5.** *Cont*.

**Figure 5.** Spatial distribution of sulphate and fluoride using inverse distance weighting (IDW) interpolation in the study area.

**Figure 6.** Scatter plot showing depth profile for F− and SO4 2−.

#### *4.2. Graphical Presentations of Water Quality Data*

In this study, physico-chemical data were further analyzed using the Piper diagram to get a deep insight into hydro-chemical processes controlling water quality of groundwater samples. The water type was calculated by plotting major cations and anions in the Piper diagram as shown in Figure 7. All water samples fell under three water facies, i.e., Na-HCO3 (50%), Mg-HCO3 (42%), and Na-Cl (8%). These water facies show chemical reactions during rock-water interaction within lithological context. Results indicate the dominance of weak acids (i.e., HCO3) over hard acids (i.e., SO4 and Cl). Also, ion exchange, dissolution of magnesite, semi-arid/arid climate, alkaline conditions as well as alumina-silicate/carbonate weathering are responsible for high concentrations of Na+, Mg2+ and HCO3 −. However, with poor management of pollution sources like untreated sewerage and agricultural runoff, there is a migration tendency from bicarbonate pole to the chloride pole. In addition, the effect of pit toilets still in practice for several local communities in the study area cannot be neglected. Henceforth, a further detailed investigation of microbial colonies (in particular, *E. coli* or fecal coliform) is a matter of research priority.

**Figure 7.** Piper diagram of 12 groundwater samples in the study area.

#### *4.3. Speciation Modelling for Selected Minerals*

Speciation modelling using the PHREEQC code [57] was performed to calculate the saturation index (SI) of selected minerals in the groundwater samples. The result is shown in Figure 8. A positive value of SI means water samples are saturated with respect to that mineral, while negative value indicates under-saturation of water samples with respect to those minerals. Here, six different minerals namely hydroxyapetite, fluorite, calcite, anhydrite, dolomite, and aragonite were selected which could possibly control the water chemistry in the study area. The results show that most of the groundwater samples are over-saturated with respect to the calcite, dolomite and aragonite minerals clearly depicting the source of HCO3 − through carbonate weathering. Looking in fluoride and hydroxyapatite, it was found that water samples were normally under-saturated except for a few samples. Looking for the saturation index for anhydrite, it was found that all water samples were unsaturated with respect to this mineral, indicating that both Ca2+ and SO4 2− are not coming from anhydrite. Also, SO4 2− is mainly coming from anthropogenic sources rather than geogenic origin.

**Figure 8.** Saturation index for selected minerals related to fluoride in groundwater.

#### *4.4. Groundwater Quality Index*

WQI was calculated for groundwater samples. The results are shown in Figure 9 and Table A1. Computed WQI values ranged from 41.97 to 239.66 with an average value of 80.37, which can be placed in four categories, namely very poor water (1), poor water (1), good water (5), and excellent water (5). Still, more than 80% of the water samples fell into the good and excellent water categories. The very poor water type of only one groundwater sample (i.e., GW5) is mainly due to the combined effect of relatively higher concentrations of TDS, EC, and NO3 − as well as high salt values. Also, looking at GW7 and GW10, WQI values were 77.78 and 176.77, respectively. These are marginal values and are susceptible to worsening into the poor and very poor categories, respectively. Therefore, timely actions for their management is very much advisable.

**Figure 9.** Spatial distribution of groundwater quality index (GWQI) using IDW interpolation in the study area.

#### *4.5. Spatial Analysis of Groundwater Quality*

Cluster analysis was used in order to classify water quality and the pollution sources for different samples. Here, with regard to the dendrogram cross section, all groundwater samples were divided into three groups (Figure 10a,b). Cluster 1 includes four groundwater samples i.e., GW2, GW3, GW7, and GW12. Cluster 2 includes five samples i.e., GW1, GW4, GW6, GW8, and GW11. Finally, cluster 3 includes the rest of the three samples including GW5, GW9, and GW10. The first group, which has four samples, is characterized in terms of where water quality is mainly governed by vegetation, agricultural activities, and close vicinity to settlements. We can compare the locations of water samples and the LULC map shown in Figure 3. It can be seen that group 1 samples are near vegetation, agricultural activities, and in close vicinity to settlements. Similarly, looking into group 2, water quality is mainly governed by hydro-chemical reactions going on both in agricultural land and mining areas, whereas water quality for group 3 is mainly being governed by hydrochemical processes going on in settlement areas.

**Figure 10.** (**a**) Dendrogram showing spatial similarities of groundwater sampling locations using cluster analysis, (**b**) results from discriminant analysis crusted in 3 groups based on the distance metric. Ward agglomeration method and dissimilarity of euclidean distance was applied.

In addition to the dendrogram, similarities and dissimilarities among three groups were evaluated with the help of factor analysis. The result is shown in Figure 5. Here, the similarity distance among clusters indicates that the groundwater samples in different clusters with lesser gaps between them had similar characteristics with respect to groundwater quality due to the same background source and type of pollution sources. Results suggest that there is a significant difference between cluster 3 and the other clusters.

The result of PCA is shown in Figure 11. Here, both *x* and *y* axes represented by F1 and F2, explained 31.06% and 22.15% of the total variance, respectively. High variation for parameters like pH, EC, TH, CO3 2−, HCO3 −, K+, Si, and Ca2+ were reported. Cluster 1 groundwater samples normally have high pH, CO3 2−, SS, and K; whereas cluster 2 has higher values of SS, Si, and CaH. Finally, cluster 3 has high EC, TDS, CO3 2−, HCO3 −, and Si.

**Figure 11.** Result for biplot diagram for three clusters in the study area.

#### **5. Conclusions**

The integrated approach using statistical analysis, speciation modeling, and graphical analysis was applied to examine the hydrochemical characteristics of groundwater quality in the Mokopane area, South Africa. The study suggests that groundwater quality is in good condition, except for a few boreholes due to sporadic pollution from F− and SO4 2− in the study area. Mainly F− pollution is of major concern due to very high concentration at sporadic locations. Results of speciation modelling for six selected minerals also suggest that most of the groundwater samples are over-saturated with respect to the calcite, dolomite, and aragonite minerals clearly depicting the source of HCO3 − through carbonate weathering. Looking into fluoride and hydroxyapatite, it was found that water samples are normally under-saturated except for few samples. The mineral composition as suggested by speciation modelling and the alkaline nature of groundwater supports the chemical reaction for fluoride enrichment through weathering. The depth profile of SO4 2− indicates that leaching of untreated sewage is the prime source of high SO4 2− concentration especially in shallow aquifers. Overall, both natural processes (i.e., rock–water interaction, silicate weathering) and poor management of anthropogenic activities (i.e., mining, leaching from untreated sewage and agricultural runoff) have a combined effect on the evolution of water quality. Computed WQI values suggest that more than 80% of the water samples fell into the good and excellent water categories. Although few samples stood in the poor and very poor categories and the water quality may further worsen if timely actions for their management are not taken. In summary, this study gives a snapshot of the status quo of water quality which is of high scientific merit despite being based on a small number of samples in the data-scarce region of Mokopane area. This study will be useful as a baseline for the planning of sustainable water resources management in the mining area in the near future to help in implementation of SDGs. The limitation of this study is the limited number of borehole samples for investigation. Hence, a detailed examination of water quality along with soil samples on a larger spatio-temporal scale can be a subject of future study.

**Author Contributions:** Conceptualization, M.D.M., R.A., P.K., H.V.T.M.; methodology, M.D.M., R.A., P.K., H.V.T.M.; investigation, D.M., R.A.; writing, M.D.M., R.A., P.K., H.V.T.M.; writing—review and editing, M.D.M., R.A., P.K., H.V.T.M., T.A.K.

**Funding:** This research received no external funding.

**Acknowledgments:** The author is grateful to the Environmental Department of Ivanplats mine for providing data for this research paper, the Faculty of Environmental Earth Science (Hokkaido University) for facilities, and anonymous reviewers for their comments. The author further extends gratitude to the JICA ABE Initiative scholarship for the opportunity to study in Hokkaido University.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **Appendix A**

**Table A1.** Calculated water quality index (WQI) for each borehole water sample.

#### **References**

© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).