Next Article in Journal
Influence of Different Types of Small Hydropower Stations on Macroinvertebrate Communities in the Changjiang River Basin, China
Next Article in Special Issue
Distribution Characteristics and Spatial Differences of Phosphorus in the Main Stream of the Urban River Stretches of the Middle and Lower Reaches of the Yangtze River
Previous Article in Journal
Correction: Gisele Cristina Rubert et al., Evapotranspiration of the Brazilian Pampa Biome: Seasonality and Influential Factors. Water 2018, 10, 1864
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Graduate School of Environmental Science, Hokkaido University, Sapporo 060-0810, Japan
2
Faculty of Environmental Earth Science, Hokkaido University, Sapporo 060-0810, Japan
3
Natural Resources and Ecosystem Services, Institute for Global Environmental Strategies, Hayama 240-0115, Japan
4
Key Laboratory of the Coastal and Wetland Ecosystems (Xiamen University), Ministry of Education College of the Environment and Ecology, Xiamen University, Xiamen 361102, Fujian, China
*
Author to whom correspondence should be addressed.
Water 2019, 11(9), 1891; https://doi.org/10.3390/w11091891
Submission received: 14 July 2019 / Revised: 6 September 2019 / Accepted: 7 September 2019 / Published: 11 September 2019
(This article belongs to the Special Issue Water Quality Assessments for Urban Water Environment)

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.

1. Introduction

Groundwater is one of the most important finite natural resources for fresh water on Earth [1,2,3]. Its importance has increased as people depend on it for agriculture and industry, with domestic use exponentially increasing [4,5,6,7,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,13,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,17,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,25,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,28,29,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,32,33,34,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.
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.

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.

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, SO42, and PO43 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,43,44,45,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).
SI   =   l o g 10 ( K I A P / K S P )
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:
  • Calculating relative weight: this was calculated using Equation (2).
    W i   =   w i i n w i
    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.
  • Calculating Q value: this was calculated using Equation (3).
    Q i =   C i × 100 S i
    where Qi is quality rating, Ci is concentration of each parameter (mg/L), and Si is derived from the WHO water quality standard.
  • Finally, WQI was calculated using Equation (4).
    W Q I   =   W i × Q i
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].

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 > SO42− > CO32−> NO3 > F> PO43−, 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 SO42− 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 SO42− can also be attributed to the leaching of organic matter and agricultural runoff carrying unutilized SO42−. Concentration of NO3 and PO43− 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.
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).
KAl 2 [ AlSi 3 O 10 ] F 2 + 2 OH KAl 2 [ AlSi 3 O 10 ] [ OH 2 ] + 2 F
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).
NaAlSi 3 O 8 +   H 2 CO 3 +   H 2 O   Na + + HCO 3 + 2 H 4 SiO 4 + 2 Al 2 Si 2 O 5 ( OH ) 4
CaF 2 + 2 HCO 3 CaCO 3 + 2 F + H 2 O + CO 2
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 SO42− 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 SO42− is found in groundwater samples with shallow well depths. As anthropogenic activities are also potential sources for SO42−, it can be expected that leaching of untreated and unmanaged domestic sewage will contaminate shallow aquifers first, thus firmly supporting the findings.

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.

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 SO42− are not coming from anhydrite. Also, SO42− is mainly coming from anthropogenic sources rather than geogenic origin.

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.

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.
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, CO32−, HCO3, K+, Si, and Ca2+ were reported. Cluster 1 groundwater samples normally have high pH, CO32−, SS, and K; whereas cluster 2 has higher values of SS, Si, and CaH. Finally, cluster 3 has high EC, TDS, CO32−, HCO3, and Si.

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 SO42− 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 SO42− indicates that leaching of untreated sewage is the prime source of high SO42− 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.
Table A1. Calculated water quality index (WQI) for each borehole water sample.
BoreholesWQIWater Quality
143.18Excellent water
264.47Good water
366.85Good water
444.78Excellent water
5239.66Very poor water
641.97Excellent water
777.78Good water
844.04Excellent water
963.76Good water
10176.77Poor water
1144.81Excellent water
1256.34Good water

References

  1. USGS. How Important Is Groundwater? Available online: https://www.usgs.gov/faqs/how-important-groundwater?qt-News_science_products=0#qt-news_science_products (accessed on 11 July 2019).
  2. Avtar, R.; Tripathi, S.; Kumar Aggarwal, A. Assessment of Energy–Population–Urbanization Nexus with Changing Energy Industry Scenario in India. Land 2019, 8, 124. [Google Scholar] [CrossRef]
  3. Avtar, R.; Tripathi, S.; Aggarwal, A.K.; Kumar, P. Population-Urbanization-Energy Nexus: A Review. Resources 2019, 8, 136. [Google Scholar] [CrossRef]
  4. Ahmad, M.; Bastiaanssen, W.G.; Feddes, R. Sustainable use of groundwater for irrigation: A numerical analysis of the subsoil water fluxes. Irrig. Drain. J. Int. Comm. Irrig. Drain. 2002, 51, 227–241. [Google Scholar] [CrossRef]
  5. Avtar, R.; Kharrazi, A. Exploring Renewable Energy Resources Using Remote Sensing and GIS—A Review. Resources 2019, 8, 149. [Google Scholar] [CrossRef]
  6. Minh, H.V.T.; Avtar, R.; Kumar, P.; Tran, D.Q.; Ty, T.V.; Behera, H.C.; Kurasaki, M. Groundwater Quality Assessment Using Fuzzy-AHP in an Giang Province of Vietnam. Geosciences 2019, 9, 330. [Google Scholar] [CrossRef]
  7. Minh, H.V.T.; Kurasaki, M.; Van Ty, T.; Tran, D.Q.; Le, K.N.; Avtar, R.; Rahman, M.; Osaki, M. Effects of Multi-Dike Protection Systems on Surface Water Quality in the Vietnamese Mekong Delta. Water 2019, 11, 1010. [Google Scholar] [CrossRef]
  8. Le, N.K.; Jha, K.M.; Jeong, J.; Gassman, W.P.; Reyes, R.M.; Doro, L.; Tran, Q.D.; Hok, L. Evaluation of Long-Term SOC and Crop Productivity within Conservation Systems Using GFDL CM2.1 and EPIC. Sustainability 2018, 10, 2665. [Google Scholar] [CrossRef]
  9. Zahedi, S. Modification of expected conflicts between Drinking Water Quality Index and Irrigation Water Quality Index in water quality ranking of shared extraction wells using Multi Criteria Decision Making techniques. Ecol. Indic. 2017, 83, 368–379. [Google Scholar] [CrossRef]
  10. Neisi, A.; Mirzabeygi Radfard, M.; Zeyduni, G.; Hamzezadeh, A.; Jalili, D.; Abbasnia, A.; Yousefi, M.; Khodadadi, R. Data on fluoride concentration levels in cold and warm season in City area of Sistan and Baluchistan Province, Iran. Data Brief. 2018, 18, 713–718. [Google Scholar] [CrossRef]
  11. Subramani, T.; Elango, L.; Damodarasamy, S.R. Groundwater quality and its suitability for drinking and agricultural use in Chithar River Basin, Tamil Nadu, India. Environ. Geol. 2005, 47, 1099–1110. [Google Scholar] [CrossRef]
  12. Avtar, R.; Kumar, P.; Singh, C.K.; Sahu, N.; Verma, R.L.; Thakur, J.K.; Mukherjee, S. Hydrogeochemical Assessment of Groundwater Quality of Bundelkhand, India Using Statistical Approach. Water Qual. Expo. Health 2013, 5, 105–115. [Google Scholar] [CrossRef]
  13. Kumar, P.; Ram, A. Chapter 4: Integrating major ion chemistry with statistical analysis for geochemical assessment of groundwater quality in coastal aquifer of Saijo plain, Ehime prefecture, Japan. In Water Quality: Indicators, Human Impact and Environmental Health; Nova Publication: Haryana, India, 2013; pp. 99–108. ISBN 978-1-62417-111-6. [Google Scholar]
  14. Avtar, R.; Kumar, P.; Surjan, A.; Gupta, L.; Roychowdhury, K. Geochemical processes regulating groundwater chemistry with special reference to nitrate and fluoride enrichment in Chhatarpur area, Madhya Pradesh, India. Environ. Earth Sci. 2013, 70, 1699–1708. [Google Scholar] [CrossRef]
  15. Matthess, G. The Properties of Groundwater; Department of Earth Science, J. Wiley and Sons Inc.: New York, NY, USA, 1982; p. 406. [Google Scholar]
  16. Kumar, P.; Kumar, M.; Ramanathan, A.L.; Tsujimura, M. Tracing the factors responsible for arsenic enrichment in groundwater of the middle Gangetic Plain, India: A source identification perspective. Environ. Geochem. Health 2010, 32, 129–146. [Google Scholar] [CrossRef] [PubMed]
  17. Dat, T.Q.; Kanchit, L.; Thares, S.; Trung, N.H. Modeling the Influence of River Discharge and Sea Level Rise on Salinity Intrusion in Mekong Delta. In Proceedings of the 1st Environment Asia International Conference, Bangkok, Thailand, 23–26 March 2011; Volume 35, pp. 685–701. [Google Scholar]
  18. Ty, T.V. Scenario-based Impact Assessment of Land Use/Cover and Climate Changes on Water Resources and Demand: A Case Study in the Srepok River Basin, Vietnam—Cambodia. Water Res. Manag. 2012, 26, 1387–1407. [Google Scholar] [CrossRef]
  19. World Health Organization (WHO). Guidelines for Drinking-Water Quality, 3rd ed.; WHO: Geneva, Switzerland, 2006; Volume 1, Available online: http://www.who.int/water_sanitation_health/dwq/gdwq0506.pdf (accessed on 8 July 2019).
  20. UNICEF. UNICEF Handbook on Water Quality; United Nations Childrens Fund: New York, NY, USA, 2008. [Google Scholar]
  21. Alagumuthu, G.; Rajan, M. Chemometric studies of water quality parameters of Sankarankovil block of Tirunelveli, Tamilnadu. J. Environ. Biol. 2010, 31, 581–586. [Google Scholar]
  22. Singh, C.K.; Rina, K.; Singh, R.; Shashtri, S.; Kamal, V.; Mukherjee, S. Geochemical modeling of high fluoride concentration in groundwater of Pokhran area of Rajasthan, India. Bull. Environ. Contam. Toxicol. 2011, 86, 152–158. [Google Scholar] [CrossRef] [PubMed]
  23. Kumar, P.; Kumar, A.; Singh, C.K.; Saraswat, C.; Avtar, R.; Ramanathan, A.; Herath, S. Hydrogeochemical evolution and appraisal of groundwater quality in Panna District, Central India. Expo. Health 2016, 8, 19–30. [Google Scholar] [CrossRef]
  24. Chapman, D.V. Water Quality Assessments: A Guide to the Use of Biota, Sediments and Water in Environmental Monitoring; CRC Press: Boca Raton, FL, USA, 1996; ISBN 0-419-21590-5. [Google Scholar]
  25. Olajire, A.A.; Imeokparia, F.E. Water Quality Assessment of Osun River: Studies on Inorganic Nutrients. Environ. Monit. Assess. 2001, 69, 17–28. [Google Scholar] [CrossRef] [PubMed]
  26. Vasanthavigar, M.; Srinivasamoorthy, K.; Prasanna, M.V. Evaluation of groundwater suitability for domestic, irrigational, and industrial purposes: A case study from Thirumanimuttar river basin, Tamilnadu, India. Environ. Monit. Assess. 2012, 184, 405–420. [Google Scholar] [CrossRef] [PubMed]
  27. Hua, A.K. Land use land cover changes in detection of water quality: A study based on remote sensing and multivariate statistics. J. Environ. Public Health 2017, 2017. [Google Scholar] [CrossRef]
  28. Khan, A.; Khan, H.H.; Umar, R. Impact of land-use on groundwater quality: GIS-based study from an alluvial aquifer in the western Ganges basin. Appl. Water Sci. 2017, 7, 4593–4603. [Google Scholar] [CrossRef] [Green Version]
  29. Narany, T.S.; Aris, A.Z.; Sefie, A.; Keesstra, S. Detecting and predicting the impact of land use changes on groundwater quality, a case study in Northern Kelantan, Malaysia. Sci. Total Environ. 2017, 599, 844–853. [Google Scholar] [CrossRef] [PubMed]
  30. Machiwal, D.; Cloutier, V.; Güler, C.; Kazakis, N. A review of GIS-integrated statistical techniques for groundwater quality evaluation and protection. Environ. Earth Sci. 2018, 77, 681. [Google Scholar] [CrossRef]
  31. Gustafsson, P. High resolution satellite data and GIS as a tool for assessment of groundwater potential of semi-arid area. In Proceedings of the IXth Thematic Conference on Geologic Remote Sensing, Pasadena, CA, USA, 8–11 February 1993. [Google Scholar]
  32. Dabral, S.; Sharma, N. An Integrated Geochemical and Geospatial Approach for Assessing the Potential Ground Water Recharge Zones in Mahi-Narmada Inter Stream Doab Area, Gujarat, India. J. Environ. Earth Sci. 2013, 3, 134–144. [Google Scholar]
  33. Demir, Y.; Erşahin, S.; Güler, M.; Cemek, B.; Günal, H.; Arslan, H. Spatial variability of depth and salinity of groundwater under irrigated ustifluvents in the Middle Black Sea Region of Turkey. Environ. Monit. Assess. 2009, 158, 279–294. [Google Scholar] [CrossRef] [PubMed]
  34. Baalousha, H. Assessment of a groundwater quality monitoring network using vulnerability mapping and geostatistics: A case study from Heretaunga Plains, New Zealand. Agric. Water Manag. 2010, 97, 240–246. [Google Scholar] [CrossRef]
  35. Dash, J.; Sarangi, A.; Singh, D. Spatial variability of groundwater depth and quality parameters in the national capital territory of Delhi. Environ. Manag. 2010, 45, 640–650. [Google Scholar] [CrossRef] [PubMed]
  36. Abbasnia, A.; Yousefi, N.; Mahvi, A.H.; Nabizadeh, R.; Radfard, M.; Yousefi, M.; Alimohammadi, M. Evaluation of groundwater quality using water quality index and its suitability for assessing water for drinking and irrigation purposes: Case study of Sistan and Baluchistan province (Iran). Hum. Ecol. Risk Assess. Int. J. 2019, 25, 988–1005. [Google Scholar] [CrossRef]
  37. Busari, O. Groundwater use in parts of the Limpopo Basin, South Africa. In Proceedings of the 3rd International Conference of Energy and Development-Environment-Biomedicine, Athens, Greece, 29–31 December 2009; pp. 13–18. [Google Scholar]
  38. Statistics. South Africa Statistics. South Africa (Web) 2016. Available online: http://www.statssa.gov.za/ (accessed on 8 July 2019).
  39. Climate Data. Climate Mokopane. Available online: https://en.climate-data.org/africa/south-africa/limpopo/mokopane-953/ (accessed on 1 July 2019).
  40. Mineral Council in South Africa. Chamber of Mines. Available online: https://www.mineralscouncil.org.za/ (accessed on 18 August 2019).
  41. USGS. Platinum-Group Elements in Southern Africa—Mineral. Inventory and an Assessment of Undiscovered Mineral Resources; US Department of the Interior: Reston, VA, USA, 2010; p. 126.
  42. Xin, X.; Li, K.; Finlayson, B.; Yin, W. Evaluation, prediction, and protection of water quality in Danjiangkou Reservoir, China. Water Sci. Eng. 2015, 8, 30–39. [Google Scholar] [CrossRef] [Green Version]
  43. Minh, H.V.T.; Ngoc, D.T.H.; Ngan, H.Y.; Men, H.V.; Van, T.N.; Ty, T.V. Assessment of Groundwater Level and Quality: A Case Study in O Mon and Binh Thuy Districts, Can Tho City, Vietnam. Naresuan Univ. Eng. J. 2016, 11, 25–33. [Google Scholar]
  44. Lever, J.; Krzywinski, M.; Altman, N. Principal component analysis. Nat. Methods 2017, 14, 641. [Google Scholar] [CrossRef]
  45. Shrestha, S.; Kazama, F. Assessment of surface water quality using multivariate statistical techniques: A case study of the Fuji river basin, Japan. Environ. Model. Softw. 2007, 22, 464–475. [Google Scholar] [CrossRef]
  46. Busico, G.; Cuoco, E.; Kazakis, N.; Colombani, N.; Mastrocicco, M.; Tedesco, D.; Voudouris, K. Multivariate statistical analysis to characterize/discriminate between anthropogenic and geogenic trace elements occurrence in the Campania Plain, Southern Italy. Environ. Pollut. 2018, 234, 260–269. [Google Scholar] [CrossRef] [PubMed]
  47. Fetter, C.W. Applied Hydrogeology, 3rd ed.; Macmillan College Publishing Company: New York, NY, USA, 1994. [Google Scholar]
  48. Garrels, R.M.; Christ, C.L. Solutions, Minerals, and Equilibria; Harper & Row: New York, NY, USA, 1965. [Google Scholar]
  49. Tiwari, T.; Mishra, M. A preliminary assignment of water quality index of major Indian rivers. Indian J. Environ. Prot. 1985, 5, 276–279. [Google Scholar]
  50. Singh, S.; Hussian, A. Water quality index development for groundwater quality assessment of Greater Noida sub-basin, Uttar Pradesh, India. Cogent Eng. 2016, 3, 1177155. [Google Scholar] [CrossRef]
  51. Fu, L.; Wang, Y.-G. Statistical Tools for Analyzing Water Quality Data. In Water Quality Monitoring and Assessment; InTech: London, UK, 2012; ISBN 978-953-51-0486-5. [Google Scholar] [Green Version]
  52. Singh, C.K.; Shashtri, S.; Mukherjee, S.; Kumari, R.; Avatar, R.; Singh, A.; Singh, R.P. Application of GWQI to Assess Effect of Land Use Change on Groundwater Quality in Lower Shiwaliks of Punjab: Remote Sensing and GIS Based Approach. Water Res. Manag. 2011, 25, 1881–1898. [Google Scholar] [CrossRef]
  53. Shepard, D. A two-dimensional interpolation function for irregularly-spaced data. In Proceedings of the 23rd ACM National Conference, Las Vegas, NV, USA, 27–29 August 1968; ACM: New York, NY, USA, 1968; pp. 517–524. [Google Scholar]
  54. Bartier, P.M.; Keller, C.P. Multivariate interpolation to incorporate thematic surface data using inverse distance weighting (IDW). Comput. Geosci. 1996, 22, 795–799. [Google Scholar] [CrossRef]
  55. Geographic Information Technology Tranning Alliance Distance-Based Interpolation. Available online: http://www.gitta.info/ContiSpatVar/en/html/Interpolatio_learningObject2.xhtml (accessed on 4 September 2019).
  56. Subramanian, V.; Saxena, K. Hydrogeochemistry of groundwater in the Delhi region of India. In Proceedings of the Hamburg Symposium, Hamburg, Germany, 18–19 August 1983. [Google Scholar]
  57. Parkhurst, D.L.; Appelo, C. User’s guide to PHREEQC (Version 2): A computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations. Water Resour. Investig. Rep. 1999, 99, 312. [Google Scholar]
Figure 1. Location of the study area sampling points and mining sites.
Figure 1. Location of the study area sampling points and mining sites.
Water 11 01891 g001
Figure 2. Climate data of average rainfall and temperature in Mokopane.
Figure 2. Climate data of average rainfall and temperature in Mokopane.
Water 11 01891 g002
Figure 3. Land-use/land-cover (LULC) map based on Landsat 2018 by using maximum likelihood classification algorithm.
Figure 3. Land-use/land-cover (LULC) map based on Landsat 2018 by using maximum likelihood classification algorithm.
Water 11 01891 g003
Figure 4. Flowchart of the methodology employed in this study for the assessment of groundwater quality.
Figure 4. Flowchart of the methodology employed in this study for the assessment of groundwater quality.
Water 11 01891 g004
Figure 5. Spatial distribution of sulphate and fluoride using inverse distance weighting (IDW) interpolation in the study area.
Figure 5. Spatial distribution of sulphate and fluoride using inverse distance weighting (IDW) interpolation in the study area.
Water 11 01891 g005
Figure 6. Scatter plot showing depth profile for F and SO42−.
Figure 6. Scatter plot showing depth profile for F and SO42−.
Water 11 01891 g006
Figure 7. Piper diagram of 12 groundwater samples in the study area.
Figure 7. Piper diagram of 12 groundwater samples in the study area.
Water 11 01891 g007
Figure 8. Saturation index for selected minerals related to fluoride in groundwater.
Figure 8. Saturation index for selected minerals related to fluoride in groundwater.
Water 11 01891 g008
Figure 9. Spatial distribution of groundwater quality index (GWQI) using IDW interpolation in the study area.
Figure 9. Spatial distribution of groundwater quality index (GWQI) using IDW interpolation in the study area.
Water 11 01891 g009
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.
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.
Water 11 01891 g010
Figure 11. Result for biplot diagram for three clusters in the study area.
Figure 11. Result for biplot diagram for three clusters in the study area.
Water 11 01891 g011
Table 1. The relative weight of chemical parameters.
Table 1. The relative weight of chemical parameters.
ParameterParameter StandardQiWeighted Factor (wi)Wi
pH9.589.4740.07
EC15049.6740.07
TDS100048.440.07
SS254440.07
Ca15031.620.04
K1564.820.04
Mg20034.3420.04
Na20027.3830.06
Cl2005.1440.07
F1.51650.09
NH40.213050.09
NO313450.09
SO44000.5820.04
Si9.222.3930.06
PO4100.5020.04
COD752430.06
Sum 541
Table 2. Statistical summary of groundwater quality in the study area.
Table 2. Statistical summary of groundwater quality in the study area.
ParameterRangeAverageSt DevWHO Permissible Limit
(WHO, 2011)
pH 7.31–9.248.300.79-
Depth (mbgl)60.0–169.0103.4233.3-
Temp. (°C)25.40–26.325.80.25-
TH (mg/L)81.62–440.46269.52105.050
EC (µs/cm)728.0–18951022.17333.54-
TDS (mg/L)473.00–1232664.25216.84-
SS (mg/L)10.00–4221.5012.11-
HCO3 (mg/L)162.50–544.6346.67106.81-
CO32− (mg/L)0.00–119.241.3243.96-
Ca2+ (mg/L)2.50–53.0219.6519.00200
K+ (mg/L)1.46–9.723.452.35200
Mg2+ (mg/L)7.64–100.5453.7524.86150
Na+ (mg/L)47.32–329.76133.5875.15200
Mn (mg/L)0.01–0.480.110.151
Cl (mg/L)7.21–188.1983.2450.76600
F (mg/L)0.24–3.391.350.981.5
NH4+ (mg/L)0.19–0.890.480.3225
NO3 (mg/L)0.06–19.964.998.0950
PO43− (mg/L)0.01–0.010.011.81E-1850
SO42− (mg/L)2.01–203.1844.7958.44400
Si (mg/L)0.43–36.4610.1212.75-
TH—Total hardness, EC—Electric Conductivity, TDS—Total dissolved solid, SS—suspended solid, mbgl—meter below ground level.

Share and Cite

MDPI and ACS Style

Molekoa, M.D.; Avtar, R.; Kumar, P.; Minh, H.V.T.; Kurniawan, T.A. Hydrogeochemical Assessment of Groundwater Quality of Mokopane Area, Limpopo, South Africa Using Statistical Approach. Water 2019, 11, 1891. https://doi.org/10.3390/w11091891

AMA Style

Molekoa MD, Avtar R, Kumar P, Minh HVT, Kurniawan TA. Hydrogeochemical Assessment of Groundwater Quality of Mokopane Area, Limpopo, South Africa Using Statistical Approach. Water. 2019; 11(9):1891. https://doi.org/10.3390/w11091891

Chicago/Turabian Style

Molekoa, Mmasabata Dolly, Ram Avtar, Pankaj Kumar, Huynh Vuong Thu Minh, and Tonni Agustiono Kurniawan. 2019. "Hydrogeochemical Assessment of Groundwater Quality of Mokopane Area, Limpopo, South Africa Using Statistical Approach" Water 11, no. 9: 1891. https://doi.org/10.3390/w11091891

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop