Next Article in Journal
Determination of Metal Content of Waste Mobile Phones and Estimation of Their Recovery Potential in Turkey
Next Article in Special Issue
Seasonal Variation of Water Quality in Unregulated Domestic Wells
Previous Article in Journal
Sleep Quality among Police Officers: Implications and Insights from a Systematic Review and Meta-Analysis of the Literature
Previous Article in Special Issue
A Hazard Assessment Method for Waterworks Systems Operating in Self-Government Units
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hydrogeochemical Characterization and Suitability Assessment of Groundwater: A Case Study in Central Sindh, Pakistan

1
School of Environmental Studies, China University of Geosciences, Wuhan 430074, China
2
College of Environmental Science & Engineering, Guilin University of Technology, Guilin 541000, China
3
Department of Dermatology, University of Health Sciences, Lahore 54000, Pakistan
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2019, 16(5), 886; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16050886
Submission received: 27 December 2018 / Revised: 1 March 2019 / Accepted: 5 March 2019 / Published: 11 March 2019
(This article belongs to the Special Issue Drinking Water and Health Risks)

Abstract

:
Groundwater is the most important water resource, on which depends human geo-economic development and survival. Recent environmental changes and anthropogenic activities render groundwater severely vulnerable. Groundwater in Central Sindh, Pakistan, is facing a similar situation. Hydrogeochemical characteristics of the groundwater in the said region were investigated by analyzing 59 groundwater samples via agricultural and drinking indices, using various statistical methods and graphical approaches to identify factors affecting groundwater. Major reactions occurring in the groundwater system were quantified by hydrogeochemical modeling. A statistical summary reveals the abundance of cations is Na+ > Ca2+ > Mg2+ > K+, while the abundance of anions is HCO3 > Cl > SO42. Groundwater chemistry is mainly of rock dominance. Correlation analysis and graphical relationships between ions reveal that ion exchange and rock weathering such as the dissolution of halite, albite, and dissolution of carbonate minerals are important rock–water interactions, governing the evolution of groundwater chemistry. Hydrochemical facies are predominantly of mixed CaMgCl and Na-Cl type, with few samples of Ca-HCO3 type, which constitutes fresh recharged water. Based on the Water Quality Index (WQI), 28.82% samples were found to be unsuitable for drinking. A United States Salinity Laboratory (USSL) diagram, Wilcox diagram, and other agricultural indices indicate that majority of the groundwater samples fall within the acceptable range for irrigation purposes.

1. Introduction

The importance of freshwater resources and its provision for every form of life is inevitable [1]. Being the elixir of life, water resources ensure the stability of the ecosystem and the status of human health [2]. Out of the total volume of water on Earth (i.e., 1600 × 106 km3), approximately 3.0 × 106 km3 of groundwater and 0.1 × 106 km3 of surface water is available as freshwater resources, which is 0.1938% of the total volume of water [3]. These approximations compel us to take strict initiatives for preservation of the available water resources to ensure its continuous and adequate availability for all living beings [4]. The key risks to freshwater resources are over-exploitation and over-consumption [5]. Since water scarcity is prevalent around the globe [6,7], groundwater exploitation emerges as an alternative to fulfill elevating water demands [8,9].
In the past few decades, global climatic change [10], urbanization, economic development [11], rapid population growth, and intensive industrialization renders the groundwater vulnerable and eventually deteriorating its quality, consequently risking environmental sustainability and preservation of life [12]. Groundwater pollution adversely affects the aqueous and terrestrial ecosystem [13], causing severe damage to human health [14,15]. Millions of people lose their lives being deprived of potable water [16] and using contaminated water [17]. Both anthropogenic and natural sources are posing significant threats to the quality of groundwater [18]. The use of pesticides, insecticides and fertilizers [19], waste from mining activities [20], industrial effluents [21], over-pumping [22], and interruptions in the river and canal networks are the prominent factors attributed to the anthropogenic activities. Moreover, natural processes include dissolution of rocks [23] and evaporation—particularly in shallow aquifers [24], which leads to water table rise and salt deposition, and seawater intrusion [25]. Considering all the adverse effects, it becomes vital to focus on groundwater management and protection from further deterioration.
Pakistan is a water-stressed country [26], not only facing a severe shortage of water [27] but is also threatened by its deteriorating quality. An abrupt decrease in both surface and groundwater resources has been noticed in the past few decades [28]. Out of the total annual water flows, Pakistan can only store 10%, which is only sufficient to meet the demand of the country for 30 days [29]. There has also been a drop in per capita availability of water from 5260 cubic meters in 1951 to 908 cubic meters in 2017, marking a significant drop of 82.73%. By 2025, with a population growth rate of 2.5% per annum, an additional population of 120 million will have to be fed. Hence, groundwater in Pakistan is mostly used for agricultural purposes besides being used for domestic and industrial purposes as well. Consequently, the demand of groundwater for domestic purposes would increase from 5.20 million acre feet to 9.70 million acre feet [30]. This drop is primarily credited to rapid population growth [31].
Sindh is the second most populated province in Pakistan. Physiographically, Sindh is bounded by the Kirthar and Laki ranges on the western side, Thar and cholistan deserts on eastern, and the Arabian sea to the south. The central part is extended from north-east to south-west, dominated by the lower Indus plain, Indus piedmont plain and Indus deltaic plain [32]. The available groundwater resource in Sindh is about 5 million acre foot (MAF) and has ample potential for irrigation however, the use of groundwater is comparatively lesser (4.3 billion cubic meter) than surface water because of two primary reasons: firstly, most of the area is lying on saline or brackish water; secondly, canal command areas are being provided with surface irrigation supplies [33]. The river Indus, having an influent behavior, loses water to the underlying aquifer, as it lies on a slight ridge, which slopes away in a lateral direction up to Larkana District. A part of the flow drains towards the desert in the east, whereas another flows towards the Khirthar Hills. During the harvest season of winter crops, the flow in the river below Sukkur Barrage becomes negligible, so the river starts receiving groundwater, especially from the left bank [33]. Studies have revealed contamination of groundwater by a variety of contaminants in Sindh province. Some of the water bodies are severely deteriorated by microbial contamination [34]; 53% of the area is affected by the calamitous outcomes of salinity and water-logging [35]; excessive fluoride exceeding the World Health Organization (WHO) and local thresholds has been reported in Nagar Parkar area [36]; and recent physiochemical analyses reveal elevated arsenic concentration in Matiari, Khairpur and Jamshoro districts [37]. Multiple approaches in the past few decades have been employed to address the hydrogeochemical characterization and groundwater quality assessment for drinking and irrigation purposes [38,39,40,41].
Categorically, the lower Indus plain comprises 14 cultivatable irrigation areas. The situation in the central Sindh, Rohri canal command area is vulnerable to surface and groundwater contamination which would ultimately risk the potability and agricultural utilization of groundwater. The literature regarding the evaluation of water quality in the study area is scarce. The present research will be explicitly helpful for the identification of hydrogeochemical characteristics and the processes governing the evolution of groundwater, and quality assessment for domestic and irrigation purposes. The findings may facilitate a clear understanding to address the adverse effects, and the solution is advocating the groundwater quality at canal command level in Sindh province.

2. Description of the Study Area

2.1. Location and Climate

The study area (Rohri canal command area) is located in the middle of the lower Indus plain and is one of the 14 annuls commonly called Canal Command Areas [42]. It stretches from longitudes 67°54′ to 69°20′ E and latitudes 24°44′ N to 27°16′ N covering an area of about 11,639 km2 (Figure 1). The climate of the area is mostly hot and dry, and can be categorized as an arid subtropical zone, i.e., hot in summer and cold in winter. The temperature rises to as high as 53 °C in summer, while in winter it drops to low as 2 °C [43]. The average annual precipitation is 260 mm [44]. The most cultivated crops in the area are wheat and cotton and are grown within the delta plain of Sindh and in the annually inundated lands [35].

2.2. Geology and Hydrogeology

The area has a rough topography, and the land slope is 0.000095 from north to south. Soils are generally coarse in texture along former meander floodplains, whereas the soil texture is smoother and finer across the flooded cover plain. The area is covered mostly by meander floodplains and consists of bar deposits, channel in-fills and channel scars (Figure 1). Fine-textured soils constitute 36% of the area and the medium soil constitute 60%. Soil surface salinity occurs in small scattered patches throughout the research area [45]. The primary reason of prevalent salinization in the area is under-watering the crops, low cropping intensity, and seepage from canals and lateral channels. Chemical analysis reveals that almost 15% of the land under study is moderate-to-highly saline with saline-alkali land constituting over 5.5% of the total land mass [46].
Quaternary sediments are a dominant feature of the geology of the area, which can often be hundreds of meters thick. Lithological studies indicate that there is a 200 to 400 feet thick sandy layer present below the surface. They are well sorted, fine to medium micaceous sands with bands and lenses of silt and silty clay [35]. The deposits are extremely variable, and the variation in lithology is so intense that it is difficult to correlate the strata found in two adjacent boreholes. However, the relative percentage of sand and clay bands is remarkably constant over large areas. The depth to water table variates from 1.56 m to 11.93 m in this aquifer with an average depth of 4.53 m. The direction of groundwater flow is defined by the individual measurement of hydraulic head which are combined to obtain the contour maps of the water level (Figure 1). The hydraulic heads in the current study area range from 10m to 28 m. In upper part, the flow direction is more south westward, while in middle and lower part, the flow direction is west ward, more toward and along river Indus. Sand is generally the predominant material making a fairly unified, transmissive and non-artesian aquifer. Groundwater yields from these sediments are typically around 50–300 m3/h. The storage factor is about 0.2 and transmissibility is 2.58 cm2/s. Vertical permeabilities are low in comparison to lateral ones but are just enough to ensure that there is a transmission of recharge from layer to layer [47].

3. Materials and Methods

3.1. Sample Collection and Analysis

In order to assess the groundwater quality of the Rohri canal command area, a total of 59 groundwater samples were collected from wells, boreholes, and hand-pumps. Before sampling, wells were pumped for several minutes to avoid the effects of stagnant water. Samples were then stored in rinsed and properly washed 2.5 L glass bottles. These samples were then analyzed in the standard water quality laboratory of the Pakistan Council of Research in Water Resources (PCRWR). A pH meter and electrochemical analyzer (Hac 44600-00, Loveland, CO, USA) were used to measure the pH, electrical conductivity (EC), and total dissolved solid (TDS). An ultraviolet-visible (UV-VIS) spectrophotometer (Analytik Jena, Jena, Germany) was used to analyze samples for major anions, i.e., nitrate (NO3-N) and sulfate (SO42), while titration method was employed to analyze chloride (Cl) and bicarbonates (HCO3) [48]. To measure major cations like Na+, K+, and Fe2+, a Flame photometer (PFP7, Cambridgeshire, UK) was used. Volumetric titration with ethylene diamine tetra acetic acid (EDTA, 0.05 N) with <2% analytical error was used to analyze Ca2+ and Mg2+. Total alkalinity in the samples was measured by acid titration using methyl-orange. Atomic absorption spectrophotometer (AAS Vario 6 Analytik Jena, Jena, Germany) was used to determine arsenic.
To validate the quality of water analyses, charge balance error (CBE) was computed. Water samples having a higher concentration of cations show positive CBE, while negative CBE is credited to higher concentrations of anions [49]. CBE was computed by Equation (1);
CBE = [ Σ cations Σ anions ] [ Σ cations + Σ anions ] × 100
where ionic concentrations are expressed in milliequivalent per liter (meq/L). According to the standard protocols, only those water samples were accepted that had less than ±5% CBE [50].
The Water Quality Index (WQI) values were computed to assess the suitability of groundwater for drinking purposes [51]. Three computing steps calculate WQI:
  • The first step involves the assignment of weight (wi) to each of the nine parameters (pH, TDS, Na, Mg, Ca, Cl, SO4, HCO3, and K) based on their relative importance to the overall quality of groundwater (Table 1). Total dissolved solids, chloride and sulfate are given a maximum weight of 5 due to their significant role in assessment, while bicarbonate is given a minimum weight of 1 because of its insignificant importance. The other parameters (pH, Na, Mg, Ca and K) are assigned weights between 1 and 5 based on their significant role while assessing the evaluation of groundwater quality.
  • The second step involves the computation of relative weight (Wi) of each parameter (Equation (2)).
    W i = w i i = 1 n w i
    where Wi is the relative weight, wi is the weight of each parameter and n is the number of parameters.
  • The third step is based on computation of the quality rating scale (qi) for each parameter (Equation (3)).
    q i = C i S i × 100
    where qi is the quality rating, Ci is the concentration of each parameter in mg/L, and Si is the WHO standard for each parameter in mg/L (Table 1).
The sub-index (SI) for each parameter is first calculated (Equation (4)) for the computation of WQI which requires the summation of sub-indices of all parameters in each sample (Equation (5)).
S I i = W i × q i
W Q I = i = 1 n S I i
where SIi and Wi is the sub-index and relative weight of the ith parameter, while qi is the rating based on the concentration of ith parameter.
For evaluating suitability for irrigation purposes, some commonly used indices like the sodium absorption ratio (SAR), residual sodium carbonate (RSC), sodium percentage (Na%), permeability index (PI), magnesium hazard (MH), Kelley’s ratio (KR) and potential salinity (PS) were computed [52]. Summary of the indices has been presented in Table 2.

3.2. Statistical Data Analysis

Statistical analyses play a pivotal role in interpreting the data set by representing various operations [53]. Regression along with Pearson correlation analysis was performed in SPSS (v23, SPSS Inc, Armonk, NY, USA) to delineate the relationships among water quality parameters. Piper diagram was prepared through Aqua-Chem (version 2010.1, Waterloo Hydrolgeologic, Kitchner, Ontario, Canada) to interpret hydrochemical facies. Saturation indices were calculated through geochemical modeling program PHREEQC (version 3.1) which determine the tendency of groundwater to dissolve or precipitate a particular mineral [54].

4. Results and Discussion

4.1. Chemical Characteristics of Groundwater

4.1.1. Hydrochemical Parameters Statistics

Hydrochemical parameters were computed in accordance with WHO protocols [55]. Table 3 presents the statistical summary of analyzed hydrochemical parameters.
The value of turbidity ranges from 0 to 165 NTU (Nephelometric Turbidity Units) with a mean value of 7.82, which is relatively elevated. The high turbidity is due shallow and poorly built wells [56]. The pH value varies from 6.5–8.2 having a mean value of 7.41, depicting neutral to slightly alkaline nature of groundwater. The mean values of EC, total dissolved solids (TDS), and total hardness (TH) are 1570.97 µS/cm, 993.92 mg/L, and 421.12 mg/L respectively, and is higher than the WHO standard apart from TDS. Dissolution of minerals and soluble salts, and evaporation of shallow groundwater are the contributing factors of elevated salinity across the study area [46]. The concentration of Na+ and K+ ranges from 17 mg/L to 638 mg/L and 0 to 25 mg/L, with average values of 165.20 mg/L and 3.53 mg/L, respectively. Mainly, the dissolution of minerals containing Na+ and K+ with surrounding rocks leads to higher concentrations [57]. Ca2+ and Mg2+ vary from 14 mg/L to 220 mg/L and 10 mg/L to 175 mg/L respectively, with mean values of 76.46 mg/L and 55.41 mg/L, which are mainly due to carbonate minerals. The concentration of HCO3 ranges from 59 mg/L to 950 mg/L, with an average concentration of 308.90 mg/L. SO42− and Cl varies from 20 mg/L to 600 mg/L and 16 mg/L to 779 mg/L respectively, with a mean values of 169.88 mg/L and 209.76 mg/L. Dissolution of gypsum, sulfate-bearing minerals and halite cause this rise in concentration [58]. According to the descriptive statistics of hydrochemical parameters, the abundance of cations is Na+ > Ca2+ > Mg2+, while the abundance of anions is HCO3 > Cl > SO42−.
Fluoride concentration varies from 0 mg/L to 2.0 mg/L with an average value of 0.42 mg/L, while most of the samples are in the permissible limits. Fluoride concentration in its extreme ranges portrays adverse effects on human health [59]. The value of NO3-N ranges from 0 mg/L to 9.90 mg/L having a mean value of 0.94 mg/L. The higher concentration of NO3-N is attributed to the excessive use of fertilizer and irrigation by wastewater but in the study area, all the samples have been found to be in the acceptable range. The concentration of Fe varies from 0 mg/L to 1 mg/L having a mean value of 0.09 mg/L.
Arsenic concentration in the research area ranges between 0 µg/L to 250 µg/L, while 15 (25.42%) samples were beyond the permissible range (10 µg/L). Elevated concentration of As in groundwater as a consequence of natural or anthropogenic sources has become a major environmental concern [60]. Climatic and geogenic sources have resulted in the varying concentration of As in groundwater, which is mainly because of moderately saline, alkalinity, and anoxic conditions [61]. It is agreed that groundwater contaminated with geogenic arsenic is more profound in alluvial aquifers [62,63]. The primary composition of these alluvial aquifers is sandstone, sand, silt, and gravel that remain in a flood plain or river channel over a longer period of time. The aquifer of the lower Indus is of similar composition [64]. As compared to other parts of the country, arsenic poisoning in the lower Indus is rather high [65], which is causing adverse effects on human health [66]. Natural mobilization of arsenic in groundwater is enhanced because of the pH-based reductive dissolution of iron hydroxide (FeOH) and competitive sorption of bicarbonate minerals in the presence of microorganisms along with evaporative enrichment. The anthropogenic causes of arsenic enrichment include water logging, coal-mining, and excessive use of pesticides [64].

4.1.2. Hydrochemical Facies

Hydrochemical facies depict the overall scenario of the interaction of groundwater solutions within a lithological structure. They are quite beneficial in interpreting the groundwater transition and pattern of its flow [67]. Piper diagram [68] presents a comprehensive graphical representation of the hydrochemistry of samples and their hydrochemical regimes.
In order to elucidate the chemical differences, all the samples were plotted in a piper diagram (Figure 2). Most of the samples fall in Zone 4 (mixed CaMgCl type) and Zone 2 (NaCl type), where the rock–water interaction significantly interprets hydrochemical behavior, anthropogenic activities, interaction with the unsaturated zone, increased resident time, ion exchange and reverse ion exchange. Few samples lie in Zone 1 (CaHCO3 type), constituting fresh recharged water. With respect to cations, groundwater samples fall in Zone D (Na + K type) and Zone B (mixed type), depicting the prominence of silicate weathering and ion exchange. With respect to anions, most of the samples fall in Zone E (HCO3 type) and Zone B (mixed type), while few samples fall in Zone G (Cl type), representing the dominance of carbonate weathering and dissolution of evaporite while the dissolution of gypsum is minimal.

4.2. Sources of Major Ions and Hydrogeochemical Evolution

4.2.1. Correlation Analysis

According to the Pearson correlation matrix (Table 4), a strong correlation between Na+ and Cl (0.92), SO42− (0.85), Mg2+ (0.60), HCO3 (0.75) exists, which depicts the derivation of Na+ from silicate weathering except the dissolution of halite and mirabilite, while the strong correlation of HCO3 with Ca2+ (0.80) and Mg2+ (0.77) illustrate the dissolution of carbonate minerals (calcite and dolomite). A strong correlation between Ca2+ and SO42− (0.73) suggests the dissolution of gypsum. The association between TDS and major constituent ions is beneficial in explaining the hydrochemical processes. TDS has been found to be strongly correlated with Na+ (0.92), Ca2+ (0.86), Cl (0.94), Mg2+ (0.84), SO42− (0.90) and HCO3 (0.86). All the ionic concentrations increase with the increase in TDS value (Figure 3).

4.2.2. Silicate Weathering

The contribution of silicate weathering plays a pivotal role in controlling the major ion chemistry of groundwater [69]. It is hard to quantify the silicate weathering products due to the incongruent nature of degradation of silicates, which yields a variety of solid phases (mostly clays) with other dissolved species [70]. Sodium which is the most abundant cation in the study area is derived from silicate weathering and by the dissolution of halite. If halite dissolution is the sole source of sodium (Equation 6), then the molar ratio of Na+/Cl would be equal to one, while the deviation from 1:1 line (Figure 4a) is attributed to the contribution of silicate weathering and ion exchange, e.g., the dissolution of albite (Equation (7)) might be responsible for the increment of Na+ in groundwater.
This fact is well supported by the saturation index of halite, which is under saturated with respect to water (Figure 5).
NaCl Na + + Cl .
2 NaAlSi 3 O 8 + 2 CO 2 + 11 H 2 O Al 2 Si 2 O 5 ( OH ) 4 + 4 H 4 SiO 4 + 2 Na + + 2 HCO 3
The influence of silicate weathering and carbonate dissolution is further illustrated by Na-normalized Ca2+ versus Mg2+ plot (Figure 6a) and Na-normalized Ca2+ versus HCO3 plot (Figure 6b). The Na-normalized Ca2+ versus Mg2+ plot shows that most of the Mg2+ is derived from carbonate dissolution except the silicate weathering, while Na-normalized Ca2+ versus HCO3 plot demonstrates that samples tend to fall close to silicate weathering rather than carbonate and evaporite dissolution [71].
If silicate weathering is a prominent source of sodium, then the most abundant anion should be HCO3 [72]. The abundance of bicarbonate among anions justifies this statement. The estimation of the ratio between Na+ and total cations (TZ+) would be helpful in understanding silicate weathering (Figure 4b). The majority of the samples fall near Na+ = 0.5 TZ+ line, which indicates the contribution of silicate weathering in the concentration of sodium in groundwater.
The contribution of dissolution of carbonates (calcite and dolomite) and sulfates (gypsum) is illustrated by the association between Ca2+ + Mg2+ and SO42− + HCO3 (Figure 4c). Samples falling on equiline (Ca2+ + Mg2+ and SO42− + HCO3) indicate dissolution of carbonates and sulfate as the sole source of Ca2+, Mg2+, and SO42− in groundwater (Equation 8, 9), while samples lying above the 1:1 line indicate less contribution of Ca2+ and Mg2+ by the dissolution of carbonates. The oversaturation of calcite and dolomite with respect to water indicated by the saturation indices of these minerals proves this statement (Figure 5).
CaCO 3 + CO 2 + H 2 O Ca 2 + + 2 HCO 3
CaMg ( CO 3 ) 2 + 2 CO 2 + 2 H 2 O Mg 2 + + Ca 2 + + 4 HCO 3
According to Figure 4d, samples falling along the 1:1 line indicate the dissolution of gypsum as a contributor of Ca2+ and SO42− ions in groundwater (Equation (10)). Moreover, under-saturation of gypsum with respect to water (Figure 5) allows it to dissolve in groundwater. However, samples falling below the equiline demonstrate the contribution of SO42− by the dissolution of Glauber’s salt (NaSO4·10 H2O) (Equation (11)), which also elevates the concentrations of Na+ ions in the groundwater.
CaSO 4 · 2 H 2 O Ca 2 + + SO 4 2 + 2 H 2 O
Na 2 SO 4 · 10 H 2 O 2 Na + + SO 4 2 + 10 H 2 O

4.2.3. Ion Exchange

Chloro-alkaline indices (CAI-1 and CAI-2) proposed by Schoeller [73] provide information about the ion exchange reactions, which have significant impacts on groundwater chemistry and its evolution [74,75]. The following formula can calculate the two indices:
CAI 1 = Cl ( Na + + K + ) Cl
CAI 2 = Cl ( Na + + K + ) HCO 3 + SO 4 2 + CO 3 2 + NO 3
where concentrations are in meq/L.
The direct ion exchange occurs when both indices give positive values expressed by Equation (14), while Equation (15) indicates the reverse ion exchange as a result of negative values of CAI-1 and CAI-2.
2 Na + + CaX 2 Ca 2 + + 2 NaX
Ca 2 + + 2 NaX 2 Na + + CaX 2
where soil exchange is expressed by X [76].
The possible role of ion exchange in hydrochemistry and evolution of groundwater can be studied by bivariate diagram between (Figure 7b) (Na+ + Cl + K+) and [(Ca2+ + Mg2+) – (SO42− + HCO3)] [77]. (Na+ + Cl + K+) indicates the increment of Na+ in the system by the processes other than halite dissolution, while [(Ca2+ + Mg2+) – (SO42− + HCO3)] represents the increment of cations Ca2+ and Mg2+ that are credited to processes excluding carbonate or silicate weathering (dissolution of gypsum, calcite, and dolomite) [78]. A linear relation between the parameters with slope equal to −1 indicates the significance of ion exchange as an important factor controlling the groundwater chemistry and its evolution [79]. According to Figure 7b, the relationship between parameters shows linearity with a slope value of −0.931 (Equation (16)), which is very close to the theoretical value of −1, indicating the ion exchange between Na+, Ca2+, and Mg2+.
y = −0.931x +0.043 (R2 = 0.895)
Most of the samples are plotted in the lower left part of Figure 7a, indicating reverse ion exchange expressed by Equation (11), which leads to the increment of Na+ and decrement of Ca2+ in groundwater. However, few samples fall in the upper right part, which signify direct ion exchange (Equation (10)).

4.2.4. Groundwater Chemistry Formation Mechanism

Groundwater resource development and quality protection require a better understanding of the key factors governing groundwater chemistry [75]. Gibbs (1970) proposed his well-known Gibbs diagram for studying groundwater chemistry formation mechanism. The two sub diagrams (Figure 8a,b) include relationship of TDS with weight ratio of (Na+ + K+) versus (Na+ + K+ + Ca2+) and TDS with weight ratio of Cl versus (Cl + HCO3). According to Gibbs’ diagram, three factors influence the groundwater chemistry, which include rock dominance, evaporation dominance and precipitation dominance [80,81]. According to Figure 8, most of the samples lie in the rock dominance zone, indicating that rock weathering is the primary source that controls the groundwater chemistry and its evolution. Parent rock weathering facilitates the process by which dissolvable salts and minerals become incorporated with groundwater. Moreover, long residence time of rock–water interaction also aids the mineral dissolution [57]. Figure 8 also shows that few samples fall in evaporation dominance zone. The influence of evaporation is attributed to the shallow groundwater depth (typically less than 3 m), which results in intense evaporation in the alluvial plain.

4.3. Groundwater Quality Assessment

The suitability of groundwater for drinking and irrigation is assessed on the basis of hydrochemical analysis data.

4.3.1. Assessment of Groundwater Quality for Drinking Purposes

WQI is one of the most comprehensive ways to address the overall quality of groundwater [82]. The composite influence of groundwater chemical parameters on overall groundwater quality is provided by the WQI, which is calculated by considering the standards recommended for drinking purposes by the WHO [55].
WQI values are classified into five categories: excellent (<50), good (>50), poor (>100), very poor (>200) and water unsuitable for drinking (>300). Table 5 illustrated the categories of groundwater quality according to WQI while the spatial distribution of water types is presented in Figure 9. The WQI values range from 35.15 to 231.46. According to Table 5, 9 samples (15.25%) are considered “excellent”, 33 samples (55.93%) as “good” while 17 samples (28.82%) were classified as “water unsuitable for drinking”.

4.3.2. Assessment of Groundwater Quality for Irrigation Purposes

The suitability of groundwater for irrigation is assessed through United States Salinity Laboratory (USSL) diagram (1954) and Wilcox diagram [83] along with some indices, such as SAR, RSC, %Na, PI, MH, KR and PS. A statistical summary of irrigation quality indices of groundwater samples is presented in Table 6.
SAR or Alkali hazard is an important tool to determine the groundwater suitability for irrigation purposes. Osmotic activity is reduced by higher salinity which hinders water to reach the branches and leaves of the plants, leading to their inferior production [84]. According to SAR, groundwater is classified as low (SAR < 10), medium (10 < SAR ≤ 18), high (18 < SAR ≤ 26) and very high (SAR > 26) sodium hazard. In the current study, the value of SAR ranges from 0.65 to 16.21, with an average of 3.41 depicting a low or medium sodium hazard. The USSL diagram (Figure 10) shows that 44 samples (73.78%) fall in C2S1 and C3S1 categories which can be used for irrigation purpose with little harm of Na+ exchange. Four samples belong to C3S2 category, six samples to C4S2 category, four samples to C4S3 category, while only one samples fall in C4S4 category. Water falling under these categories cannot be recommended for any agricultural practices.
RSC is an important index for the evaluation of groundwater suitability for irrigation. It can be estimated by subtracting alkaline earth metals (Ca2+ + Mg2+) from carbonates and bicarbonates (CO32− + HCO3). Adsorption of sodium in soil is attributed to the high value of RSC [85]. Groundwater is considered unsuitable if RSC value exceeds 2.5. According to RSC, 98.30% samples of the present study have been found to be suitable for irrigation purposes.
Irrigation water with high sodium percentage reduces the permeability of the soil, which consequently affects the plant growth, thus making it an important index for the evaluation of groundwater on agricultural scale [86]. Groundwater having less than 60% Na is considered suitable for irrigation use. %Na value ranges from 15.64 to 86.32 having a mean value of 41.70. 91.52% samples fall in excellent to permissible limits. According to Wilcox diagram (Figure 11), 42 samples (71.19%) lie in excellent to permissible zone, five samples in permissible to doubtful, while seven samples in doubtful to unsuitable and only five samples are in unsuitable zone.
The influence of ion contents in groundwater on the permeability of the soil is evaluated by Permeability Index (PI) proposed by Doneen [87]. Based on PI, values less than 25 are considered unsuitable for irrigation. In the current study, all samples have been found suitable for irrigation on PI basis.
Another tool for assessing the agricultural suitability of groundwater is MH, proposed by Szabolcs and Darab [88]. Sodication in the soils increases with the subsequent increment of Mg2+ content relative to Ca2+, which consequently results in the dispersion of clay particles hence damaging the soil structure by decreasing the hydraulic conductivity. MH value less than 50 is considered suitable for irrigation. Table 6 shows that 69.49% of samples have been found in the suitable range.
KR is another important indicator for the evaluation of groundwater for agricultural suitability proposed by Kelly [89]. A Kelly’s ratio value greater than 1 indicates excess level of sodium in groundwater, while less than 1 Kelly’s ratio is considered suitable for irrigation [90]. In the current study, as per KR, 42 samples (71.19%) fall in the suitable range for irrigation purposes.
The agricultural suitability of groundwater is also assessed by PS which is defined as the chloride concentration plus half of the sulfate concentration. Doneen [91] proposed that agricultural suitability is not dependent on the concentration of soluble salts. Every successive irrigation helps in the precipitation and accumulation of low soluble salts, whereas the salinity of the soil is enhanced by the concentration of highly soluble salts. Groundwater is considered suitable, when PS values are less than 10. For the present study, PS values show that 43 samples (72.89%) are suitable for irrigation purposes.

5. Conclusions

In order to carry out the current study, an amalgamation of statistical analyses, graphical techniques and hydrochemical modeling were used for comprehensive understanding of the groundwater chemistry, its evolution and suitability for drinking and agricultural purposes in the delimited area of research.
Statistical analyses demonstrate that the abundance of cations is in the order: Na+ > Ca2+ > Mg2+ > K+, while the abundance of anions is in the order: HCO3 > Cl > SO42−. The abundance of Na+ and HCO3 proves that silicate weathering is the eminent phenomenon controlling the major ion chemistry of the groundwater in the research area. This is also supported by the correlation analysis and graphical relationships between ions that reveal that ion exchange and rock weathering—like the dissolution of halite, albite, and dissolution of carbonate minerals—are important rock–water interactions, which govern the evolution of groundwater chemistry. Furthermore, rock dominance has been found as the key natural factor governing the groundwater evolution as depicted by the Gibbs diagrams, while few samples falling in the evaporation dominance propose the influence of evaporation in shallow groundwater depth zone. Hydrochemical facies are predominantly of mixed CaMgCl and Na-Cl type with few samples of Ca-HCO3 type, constituting fresh recharged water. Saturation indices computed by geochemical modelling indicate that the aqueous phase is under-saturated with respect to evaporites like halite and gypsum (negative values), while there is equilibrium to the over-saturated phase with respect to carbonates like calcite and dolomite (positive values).
Computed values of WQI shows that 71.18% samples are safe for drinking purpose while 28.82% fall in unsuitable range. The contributing factors towards the aforementioned range of groundwater samples include dissolution process and effective leaching of rock salt and gypsum-bearing rock formations. The majority of samples have been found suitable for irrigation purposes as per the USSL diagram, Wilcox diagram and other agricultural indices.

Author Contributions

Conceptualization, M.A.T. and Z.T.; Data curation, A.S.; Formal analysis, M.A.T., J.S. and M.F.; Investigation, M.A.T. and A.S.; Methodology, M.A.T., Z.T. and J.S.; Resources, A.S.; Software, M.A.T., A.S. and M.F.; Supervision, Z.T.; Validation, A.S. and J.S.; Visualization, M.A.T. and A.S.; Writing—original draft, M.A.T., J.S. and M.F.; Writing—review and editing, J.S. and M.F.

Funding

This research received no external funding.

Acknowledgments

The authors would like to acknowledge the Pakistan Council for Research in Water Resources (PCRWR) and Federal Flood Commission, Ministry of Water and Power, Pakistan for the material support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Holland, R.A.; Scott, K.A.; Flörke, M.; Brown, G.; Ewers, R.M.; Farmer, E.; Kapos, V.; Muggeridge, A.; Scharlemann, J.P.; Taylor, G. Global impacts of energy demand on the freshwater resources of nations. Proc. Natl. Acad. Sci. USA 2015, 112, E6707–E6716. [Google Scholar] [CrossRef] [Green Version]
  2. Shen, J.; Lu, H.; Zhang, Y.; Song, X.; He, L. Vulnerability assessment of urban ecosystems driven by water resources, human health and atmospheric environment. J. Hydrol. 2016, 536, 457–470. [Google Scholar] [CrossRef]
  3. Döll, P.; Douville, H.; Güntner, A.; Schmied, H.M.; Wada, Y. Modelling freshwater resources at the global scale: Challenges and prospects. Surv. Geophys. 2016, 37, 195–221. [Google Scholar] [CrossRef]
  4. Andreoli, R.; Van Haaren, B. Ecosystemic Services Assessment: Application to Forests for the Preservation of Water Resources in Tropical Islands. QGIS Appl. Territ. Plan. 2018, 3, 169–237. [Google Scholar]
  5. Pellegrini, G.; Ingrao, C.; Camposeo, S.; Tricase, C.; Conto, F.; Huisingh, D. Application of water footprint to olive growing systems in the Apulia region: A comparative assessment. J. Clean. Prod. 2016, 112, 2407–2418. [Google Scholar] [CrossRef]
  6. Veldkamp, T.I.; Wada, Y.; de Moel, H.; Kummu, M.; Eisner, S.; Aerts, J.C.; Ward, P.J. Changing mechanism of global water scarcity events: Impacts of socioeconomic changes and inter-annual hydro-climatic variability. Glob. Environ. Chang. 2015, 32, 18–29. [Google Scholar] [CrossRef] [Green Version]
  7. Schewe, J.; Heinke, J.; Gerten, D.; Haddeland, I.; Arnell, N.W.; Clark, D.B.; Dankers, R.; Eisner, S.; Fekete, B.M.; Colón-González, F.J. Multimodel assessment of water scarcity under climate change. Proc. Natl. Acad. Sci. USA 2014, 111, 3245–3250. [Google Scholar] [CrossRef] [PubMed]
  8. Amin, A.; Iqbal, J.; Asghar, A.; Ribbe, L. Analysis of Current and Future Water Demands in the Upper Indus Basin under IPCC Climate and Socio-Economic Scenarios Using a Hydro-Economic WEAP Model. Water 2018, 10, 537. [Google Scholar] [CrossRef]
  9. Shang, H.; Wang, W.; Dai, Z.; Duan, L.; Zhao, Y.; Zhang, J. An ecology-oriented exploitation mode of groundwater resources in the northern Tianshan Mountains, China. J. Hydrol. 2016, 543, 386–394. [Google Scholar] [CrossRef]
  10. Cheeseman, J. Food security in the face of salinity, drought, climate change, and population growth. In Halophytes for Food Security in Dry Lands; Elsevier: Amsterdam, The Netherlands, 2016; pp. 111–123. [Google Scholar] [CrossRef]
  11. Hua, S.; Liang, J.; Zeng, G.; Xu, M.; Zhang, C.; Yuan, Y.; Li, X.; Li, P.; Liu, J.; Huang, L. How to manage future groundwater resource of China under climate change and urbanization: An optimal stage investment design from modern portfolio theory. Water Res. 2015, 85, 31–37. [Google Scholar] [CrossRef] [PubMed]
  12. Wichelns, D.; Qadir, M. Achieving sustainable irrigation requires effective management of salts, soil salinity, and shallow groundwater. Agric. Water Manag. 2015, 157, 31–38. [Google Scholar] [CrossRef]
  13. Kløve, B.; Ala-Aho, P.; Bertrand, G.; Gurdak, J.J.; Kupfersberger, H.; Kværner, J.; Muotka, T.; Mykrä, H.; Preda, E.; Rossi, P. Climate change impacts on groundwater and dependent ecosystems. J. Hydrol. 2014, 518, 250–266. [Google Scholar] [CrossRef] [Green Version]
  14. Wu, J.; Sun, Z. Evaluation of shallow groundwater contamination and associated human health risk in an alluvial plain impacted by agricultural and industrial activities, mid-west China. Expo. Health 2016, 8, 311–329. [Google Scholar] [CrossRef]
  15. Wongsasuluk, P.; Chotpantarat, S.; Siriwong, W.; Robson, M. Heavy metal contamination and human health risk assessment in drinking water from shallow groundwater wells in an agricultural area in Ubon Ratchathani province, Thailand. Environ. Geochem. Health 2014, 36, 169–182. [Google Scholar] [CrossRef] [PubMed]
  16. Bhowmick, S.; Pramanik, S.; Singh, P.; Mondal, P.; Chatterjee, D.; Nriagu, J. Arsenic in groundwater of West Bengal, India: A review of human health risks and assessment of possible intervention options. Sci. Total Environ. 2018, 612, 148–169. [Google Scholar] [CrossRef] [PubMed]
  17. Zhai, Y.; Zhao, X.; Teng, Y.; Li, X.; Zhang, J.; Wu, J.; Zuo, R. Groundwater nitrate pollution and human health risk assessment by using HHRA model in an agricultural area, NE China. Ecotoxicol. Environ. Saf. 2017, 137, 130–142. [Google Scholar] [CrossRef] [PubMed]
  18. Khatri, N.; Tyagi, S. Influences of natural and anthropogenic factors on surface and groundwater quality in rural and urban areas. Front. Life Sci. 2015, 8, 23–39. [Google Scholar] [CrossRef]
  19. Douglas, S.H.; Dixon, B.; Griffin, D. Assessing intrinsic and specific vulnerability models ability to indicate groundwater vulnerability to groups of similar pesticides: A comparative study. Phys. Geogr. 2018, 1–19. [Google Scholar] [CrossRef]
  20. Nagaraju, A.; Kumar, K.S.; Thejaswi, A. Assessment of groundwater quality for irrigation: A case study from Bandalamottu lead mining area, Guntur District, Andhra Pradesh, South India. Appl. Water Sci. 2014, 4, 385–396. [Google Scholar] [CrossRef]
  21. Kanagaraj, G.; Elango, L. Hydrogeochemical processes and impact of tanning industries on groundwater quality in Ambur, Vellore district, Tamil Nadu, India. Environ. Sci. Pollut. Res. 2016, 23, 24364–24383. [Google Scholar] [CrossRef] [PubMed]
  22. Hoogesteger, J.; Wester, P. Intensive groundwater use and (in) equity: Processes and governance challenges. Environ. Sci. Policy 2015, 51, 117–124. [Google Scholar] [CrossRef]
  23. Varol, S.; Davraz, A. Evaluation of the groundwater quality with WQI (Water Quality Index) and multivariate analysis: A case study of the Tefenni plain (Burdur/Turkey). Environ. Earth Sci. 2015, 73, 1725–1744. [Google Scholar] [CrossRef]
  24. Li, P.; Wu, J.; Qian, H. Hydrochemical appraisal of groundwater quality for drinking and irrigation purposes and the major influencing factors: A case study in and around Hua County, China. Arab. J. Geosci. 2016, 9, 15. [Google Scholar] [CrossRef]
  25. Tomaszkiewicz, M.; Najm, M.A.; El-Fadel, M. Development of a groundwater quality index for seawater intrusion in coastal aquifers. Environ. Modell. Softw. 2014, 57, 13–26. [Google Scholar] [CrossRef]
  26. Alamgir, A.; Ali Khan, M.; Manino, I.; Shaukat, S.S.; Shahab, S. Vulnerability to climate change of surface water resources of coastal areas of Sindh, Pakistan. Desalin. Water Treat. 2016, 57, 18668–18678. [Google Scholar] [CrossRef]
  27. Dalin, C.; Wada, Y.; Kastner, T.; Puma, M.J. Groundwater depletion embedded in international food trade. Nature 2017, 543, 700. [Google Scholar] [CrossRef] [PubMed]
  28. MacDonald, A.; Bonsor, H.; Ahmed, K.; Burgess, W.; Basharat, M.; Calow, R.; Dixit, A.; Foster, S.; Gopal, K.; Lapworth, D. Groundwater quality and depletion in the Indo-Gangetic Basin mapped from in situ observations. Nat. Geosci. 2016, 9, 762–766. [Google Scholar] [CrossRef] [Green Version]
  29. Farooqui, S.Z. Prospects of renewables penetration in the energy mix of Pakistan. Renew. Sustain. Energy Rev. 2014, 29, 693–700. [Google Scholar] [CrossRef]
  30. Raza, M.; Hussain, F.; Lee, J.-Y.; Shakoor, M.B.; Kwon, K.D. Groundwater status in Pakistan: A review of contamination, health risks, and potential needs. Crit. Rev. Environ. Sci. Technol. 2017, 47, 1713–1762. [Google Scholar] [CrossRef]
  31. Kirby, M.; Mainuddin, M.; Khaliq, T.; Cheema, M. Agricultural production, water use and food availability in Pakistan: Historical trends, and projections to 2050. Agric. Water Manag. 2017, 179, 34–46. [Google Scholar] [CrossRef]
  32. Alamgir, A.; Khan, M.A.; Schilling, J.; Shaukat, S.S.; Shahab, S. Assessment of groundwater quality in the coastal area of Sindh province, Pakistan. Environ. Monit. Assess. 2016, 188, 78. [Google Scholar] [CrossRef] [PubMed]
  33. van Steenbergen, F.; Basharat, M.; Lashari, B.K. Key challenges and opportunities for conjunctive management of surface and groundwater in mega-irrigation systems: Lower Indus, Pakistan. Resources 2015, 4, 831–856. [Google Scholar] [CrossRef]
  34. Waseem, A.; Arshad, J.; Iqbal, F.; Sajjad, A.; Mehmood, Z.; Murtaza, G. Pollution status of Pakistan: A retrospective review on heavy metal contamination of water, soil, and vegetables. BioMed Res. Int. 2014, 2014. [Google Scholar] [CrossRef] [PubMed]
  35. Qureshi, A.S.; McCornick, P.G.; Qadir, M.; Aslam, Z. Managing salinity and waterlogging in the Indus Basin of Pakistan. Agric. Water Manag. 2008, 95, 1–10. [Google Scholar] [CrossRef]
  36. Rafique, T.; Naseem, S.; Usmani, T.H.; Bashir, E.; Khan, F.A.; Bhanger, M.I. Geochemical factors controlling the occurrence of high fluoride groundwater in the Nagar Parkar area, Sindh, Pakistan. J. Hazard. Mater. 2009, 171, 424–430. [Google Scholar] [CrossRef] [PubMed]
  37. Memon, A.H.; Ghanghro, A.B.; Jahangir, T.M.; Lund, G.M. Arsenic contamination in drinking water of district Jamshoro, Sindh, Pakistan. Biomed. Lett. 2016, 2, 31–37. [Google Scholar]
  38. Dey, S.; Raju, N.J.; Ram, P.; Singh, J. Hydrogeochemical characterisation and evaluation of seasonal variation in groundwater chemistry in upper panda river basin, India. In Management of Water, Energy and Bio-resources in the Era of Climate Change: Emerging Issues and Challenges; Springer: Berlin, Germany, 2015; pp. 21–36. [Google Scholar] [CrossRef]
  39. Sethy, S.N.; Syed, T.H.; Kumar, A.; Sinha, D. Hydrogeochemical characterization and quality assessment of groundwater in parts of Southern Gangetic Plain. Environ. Earth Sci. 2016, 75, 232. [Google Scholar] [CrossRef]
  40. Gu, X.; Xiao, Y.; Yin, S.; Hao, Q.; Liu, H.; Hao, Z.; Meng, G.; Pei, Q.; Yan, H. Hydrogeochemical Characterization and Quality Assessment of Groundwater in a Long-Term Reclaimed Water Irrigation Area, North China Plain. Water 2018, 10, 1209. [Google Scholar] [CrossRef]
  41. Sebei, A.; Slama, T.; Helali, M.A. Hydrochemical characterization and geospatial analysis of groundwater quality in Cap Bon region, northeastern Tunisia. Environ. Earth Sci. 2018, 77, 557. [Google Scholar] [CrossRef]
  42. Aslam, K.; Rashid, S.; Saleem, R.; Aslam, R.M.S. Use of geospatial technology for assessment of waterlogging & salinity conditions in the Nara Canal Command area in Sindh, Pakistan. J. Geogr. Inf. Syst. 2015, 7, 438. [Google Scholar]
  43. Hussain, I.; Marikar, F.; Jehangir, W.A. Productivity and Performance of Irrigated Wheat Farms Across Canal Commands in the Lower Indus Basin; IWMI: Colombo, Srilanka, 2000; Volume 44. [Google Scholar]
  44. Salma, S.; Shah, M.; Rehman, S. Rainfall trends in different climate zones of Pakistan. Pak. J. Meteorol. 2012, 9, 37–47. [Google Scholar]
  45. Shahab, A.; Shihua, Q.; Rashid, A.; Hasan, F.U.; Sohail, M.T. Evaluation of Water Quality for Drinking and Agricultural Suitability in the Lower Indus Plain in Sindh Province, Pakistan. Pol. J. Environ. Stud. 2016, 25. [Google Scholar] [CrossRef]
  46. Qureshi, A.; Lashari, B.; Kori, S.; Lashari, G. Hydro-salinity behavior of shallow groundwater aquifer underlain by salty groundwater in Sindh Pakistan. In Proceedings of the Fifteenth International Water Technology Conference, Alexandria, Egypt, 28–30 May 2011. [Google Scholar]
  47. Jain, S.; Sharma, B.; Zahid, A.; Jin, M.; Shreshtha, J.; Kumar, V.; Rai, S.; Hu, J.; Luo, Y.; Sharma, D. A comparative analysis of the hydrogeology of the Indus-Gangetic and Yellow River basins. In Groundwater Governance in the Indo-Gangetic and Yellow River Basins; CRC Press: Boca Raton, FL, USA, 2009; pp. 63–84. [Google Scholar] [CrossRef]
  48. Federation, W.E.; Association, A.P.H. Standard Methods for the Examination of Water and Wastewater; American Public Health Association (APHA): Washington, DC, USA, 2005. [Google Scholar]
  49. Bozdağ, A. Combining AHP with GIS for assessment of irrigation water quality in Çumra irrigation district (Konya), Central Anatolia, Turkey. Environ. Earth Sci. 2015, 73, 8217–8236. [Google Scholar] [CrossRef]
  50. Singh, S.; Raju, N.J.; Ramakrishna, C. Evaluation of groundwater quality and its suitability for domestic and irrigation use in parts of the Chandauli-Varanasi region, Uttar Pradesh, India. J. Water Resour. Prot. 2015, 7, 572. [Google Scholar] [CrossRef]
  51. WHO. Guidelines for drinking-water quality: First addendum to the fourth edition. 2017. [Google Scholar]
  52. Vincy, M.; Brilliant, R.; Pradeepkumar, A. Hydrochemical characterization and quality assessment of groundwater for drinking and irrigation purposes: A case study of Meenachil River Basin, Western Ghats, Kerala, India. Environ. Monit. Assess. 2015, 187, 4217. [Google Scholar] [CrossRef] [PubMed]
  53. Sohail, M.T.; Delin, H.; Talib, M.A.; Xiaoqing, X.; Akhtar, M.M. An Analysis of Environmental Law in Pakistan-policy and Conditions of Implementation. Res. J. Appl. Sci. Eng. Technol. 2014, 8, 644–653. [Google Scholar] [CrossRef]
  54. Raj, D.; Shaji, E. Fluoride contamination in groundwater resources of Alleppey, southern India. Geosci. Front. 2017, 8, 117–124. [Google Scholar] [CrossRef] [Green Version]
  55. Edition, F. Guidelines for drinking-water quality. WHO Chron. 2011, 38, 104–108. [Google Scholar]
  56. Azis, A.; Yusuf, H.; Faisal, Z.; Suradi, M. Water turbidity impact on discharge decrease of groundwater recharge in recharge reservoir. Procedia Eng. 2015, 125, 199–206. [Google Scholar] [CrossRef]
  57. Selvakumar, S.; Ramkumar, K.; Chandrasekar, N.; Magesh, N.; Kaliraj, S. Groundwater quality and its suitability for drinking and irrigational use in the Southern Tiruchirappalli district, Tamil Nadu, India. Appl. Water Sci. 2017, 7, 411–420. [Google Scholar] [CrossRef]
  58. Wu, J.; Wang, L.; Wang, S.; Tian, R.; Xue, C.; Feng, W.; Li, Y. Spatiotemporal variation of groundwater quality in an arid area experiencing long-term paper wastewater irrigation, northwest China. Environ. Earth Sci. 2017, 76, 460. [Google Scholar] [CrossRef]
  59. Rafique, T.; Naseem, S.; Ozsvath, D.; Hussain, R.; Bhanger, M.I.; Usmani, T.H. Geochemical controls of high fluoride groundwater in Umarkot sub-district, Thar Desert, Pakistan. Sci. Total Environ. 2015, 530, 271–278. [Google Scholar] [CrossRef] [PubMed]
  60. Shakoor, M.B.; Niazi, N.K.; Bibi, I.; Rahman, M.M.; Naidu, R.; Dong, Z.; Shahid, M.; Arshad, M. Unraveling health risk and speciation of arsenic from groundwater in rural areas of Punjab, Pakistan. Int. J. Environ. Res. Public Health 2015, 12, 12371–12390. [Google Scholar] [CrossRef] [PubMed]
  61. Baig, J.A.; Kazi, T.G.; Mustafa, M.A.; Solangi, I.B.; Mughal, M.J.; Afridi, H.I. Arsenic exposure in children through drinking water in different districts of Sindh, Pakistan. Biol. Trace Elem. Res. 2016, 173, 35–46. [Google Scholar] [CrossRef] [PubMed]
  62. Su, C.; Zhu, Y.; Abbas, Z.; Huq, M.E. Sources and controls for elevated arsenic concentrations in groundwater of Datong Basin, Northern China. Environ. Earth Sci. 2016, 75, 570. [Google Scholar] [CrossRef]
  63. Huq, M.E.; Su, C.; Fahad, S.; Li, J.; Sarven, M.S.; Liu, R. Distribution and hydrogeochemical behavior of arsenic enriched groundwater in the sedimentary aquifer comparison between Datong Basin (China) and Kushtia District (Bangladesh). Environ. Sci. Pollut. Res. 2018, 1–14. [Google Scholar] [CrossRef] [PubMed]
  64. Shahab, A.; Qi, S.; Zaheer, M. Arsenic contamination, subsequent water toxicity, and associated public health risks in the lower Indus plain, Sindh province, Pakistan. Environ. Sci. Pollut. Res. 2018. [Google Scholar] [CrossRef] [PubMed]
  65. Sanjrani, M.; Mek, T.; Sanjrani, N.; Leghari, S.; Moryani, H.; Shabnam, A. Current situation of aqueous arsenic contamination in Pakistan, focused on Sindh and Punjab Province, Pakistan: A review. J. Pollut. Eff. Control 2017, 5, 207. [Google Scholar]
  66. Yunus, F.; Khan, S.; Chowdhury, P.; Milton, A.; Hussain, S.; Rahman, M. A review of groundwater arsenic contamination in Bangladesh: The millennium development goal era and beyond. Int. J. Environ. Res. Public Health 2016, 13, 215. [Google Scholar] [CrossRef] [PubMed]
  67. Jain, C.; Sharma, S.; Singh, S. Physico-chemical characteristics and hydrogeological mechanisms in groundwater with special reference to arsenic contamination in Barpeta District, Assam (India). Environ. Monit. Assess. 2018, 190, 417. [Google Scholar] [CrossRef] [PubMed]
  68. Piper, A.M. A graphic procedure in the geochemical interpretation of water-analyses. Eos Trans. Am. Geophys. Union 1944, 25, 914–928. [Google Scholar] [CrossRef]
  69. Subramani, T.; Rajmohan, N.; Elango, L. Groundwater geochemistry and identification of hydrogeochemical processes in a hard rock region, Southern India. Environ. Monit. Assess. 2010, 162, 123–137. [Google Scholar] [CrossRef] [PubMed]
  70. Saravanan, K.; Srinivasamoorthy, K.; Gopinath, S.; Prakash, R.; Suma, C. Investigation of hydrogeochemical processes and groundwater quality in Upper Vellar sub-basin Tamilnadu, India. Arab. J. Geosci. 2016, 9, 372. [Google Scholar] [CrossRef]
  71. Mukherjee, A.; Fryar, A.E. Deeper groundwater chemistry and geochemical modeling of the arsenic affected western Bengal basin, West Bengal, India. Appl. Geochem. 2008, 23, 863–894. [Google Scholar] [CrossRef]
  72. Rajmohan, N.; Elango, L. Identification and evolution of hydrogeochemical processes in the groundwater environment in an area of the Palar and Cheyyar River Basins, Southern India. Environ. Geol. 2004, 46, 47–61. [Google Scholar] [CrossRef]
  73. Schoeller, H. Qualitative Evaluation of Groundwater Resources. Methods and Techniques of Groundwater Investigations and Development; UNESCO: Paris, France, 1965; Volume 5483. [Google Scholar]
  74. Zaidi, F.K.; Nazzal, Y.; Jafri, M.K.; Naeem, M.; Ahmed, I. Reverse ion exchange as a major process controlling the groundwater chemistry in an arid environment: A case study from northwestern Saudi Arabia. Environ. Monit. Assess. 2015, 187, 607. [Google Scholar] [CrossRef] [PubMed]
  75. Aghazadeh, N.; Mogaddam, A.A. Assessment of groundwater quality and its suitability for drinking and agricultural uses in the Oshnavieh area, Northwest of Iran. J. Environ. Prot. 2010, 1, 30. [Google Scholar] [CrossRef]
  76. Anders, R.; Mendez, G.O.; Futa, K.; Danskin, W.R. A geochemical approach to determine sources and movement of saline groundwater in a coastal aquifer. Groundwater 2014, 52, 756–768. [Google Scholar] [CrossRef] [PubMed]
  77. Ahmed, M.; Samie, S.A.; Badawy, H. Factors controlling mechanisms of groundwater salinization and hydrogeochemical processes in the Quaternary aquifer of the Eastern Nile Delta, Egypt. Environ. Earth Sci. 2013, 68, 369–394. [Google Scholar] [CrossRef]
  78. Yidana, S.M.; Yidana, A. Assessing water quality using water quality index and multivariate analysis. Environ. Earth Sci. 2010, 59, 1461–1473. [Google Scholar] [CrossRef]
  79. Rajabpour, H.; Vaezihir, A.; Sedghi, M.H. The North Tabriz Fault, a barrier to groundwater flow in an alluvial aquifer northwest of Tabriz, Iran. Environ. Earth Sci. 2016, 75, 849. [Google Scholar] [CrossRef]
  80. Gibbs, R.J. Mechanisms controlling world water chemistry. Science 1970, 170, 1088–1090. [Google Scholar] [CrossRef] [PubMed]
  81. Salem, Z.E.; Atwia, M.G.; El-Horiny, M.M. Hydrogeochemical analysis and evaluation of groundwater in the reclaimed small basin of Abu Mina, Egypt. Hydrogeol. J. 2015, 23, 1781–1797. [Google Scholar] [CrossRef]
  82. Sadat-Noori, S.; Ebrahimi, K.; Liaghat, A. Groundwater quality assessment using the Water Quality Index and GIS in Saveh-Nobaran aquifer, Iran. Environ. Earth Sci. 2014, 71, 3827–3843. [Google Scholar] [CrossRef]
  83. Wilcox, L.V. The Quality of Water for Irrigation Use; United States Department of Agriculture, Economic Research Service: Washington, DC, USA, 1948.
  84. Marghade, D.; Malpe, D.; Zade, A. Geochemical characterization of groundwater from northeastern part of Nagpur urban, Central India. Environ. Earth Sci. 2011, 62, 1419–1430. [Google Scholar] [CrossRef]
  85. Eaton, F.M. Significance of carbonates in irrigation waters. Soil Sci. 1950, 69, 123–134. [Google Scholar] [CrossRef]
  86. Amiri, V.; Sohrabi, N.; Dadgar, M.A. Evaluation of groundwater chemistry and its suitability for drinking and agricultural uses in the Lenjanat plain, central Iran. Environ. Earth Sci. 2015, 74, 6163–6176. [Google Scholar] [CrossRef]
  87. Doneen, L.D. Notes on Water Quality in Agriculture; Department of Water Science and Engineering, University of California: Davis, CA, USA, 1964. [Google Scholar]
  88. Szabolcs, I.; Darab, C. The influence of irrigation water of high sodium carbonate content of soils. In Proceedings of the 8th International Congress of ISSS, Trans, II, Bucharest, Romania, 31 August–9 September 1964; Research Institute for Soil Science and Agricultural Chemistry of the Hungarian Academy of Sciences: Budapest, Hungary, 1964; pp. 803–812. [Google Scholar]
  89. Kelly, W. Use of saline irrigation water. Soil Sci. 1963, 95, 355–391. [Google Scholar] [CrossRef]
  90. Bouderbala, A. Assessment of groundwater quality and its suitability for agricultural uses in the Nador Plain, north of Algeria. Water Qual. Expo. Health 2015, 7, 445–457. [Google Scholar] [CrossRef]
  91. Doneen, L.D. Salination of soil by salts in the irrigation water. Eos Trans. Am. Geophys. Union 1954, 35, 943–950. [Google Scholar] [CrossRef]
Figure 1. Location of the study area and sampling stations (Rohri Canal Command).
Figure 1. Location of the study area and sampling stations (Rohri Canal Command).
Ijerph 16 00886 g001
Figure 2. Piper plot for groundwater facies classification.
Figure 2. Piper plot for groundwater facies classification.
Ijerph 16 00886 g002
Figure 3. Scatter plot for (a) the main cations; (b) the main anions versus total dissolved solid (TDS).
Figure 3. Scatter plot for (a) the main cations; (b) the main anions versus total dissolved solid (TDS).
Ijerph 16 00886 g003
Figure 4. Relationship between ions (a) Na+ versus Cl (b) Na+ versus TZ+ (c) Ca2+ + Mg2+ versus SO42− + HCO3 (d) Ca2+ versus SO42.
Figure 4. Relationship between ions (a) Na+ versus Cl (b) Na+ versus TZ+ (c) Ca2+ + Mg2+ versus SO42− + HCO3 (d) Ca2+ versus SO42.
Ijerph 16 00886 g004
Figure 5. Variation of saturation indices of selected minerals.
Figure 5. Variation of saturation indices of selected minerals.
Ijerph 16 00886 g005
Figure 6. Plots of (a) Na-normalized Ca2+ versus Mg2+ (b) Na-normalized Ca2+versus HCO3.
Figure 6. Plots of (a) Na-normalized Ca2+ versus Mg2+ (b) Na-normalized Ca2+versus HCO3.
Ijerph 16 00886 g006
Figure 7. Plots showing (a) CAI-1 versus CAI-2; (b) (Ca2+ + Mg2+) − (HCO3 + SO42−) versus (Na+ + K+) − Cl−.
Figure 7. Plots showing (a) CAI-1 versus CAI-2; (b) (Ca2+ + Mg2+) − (HCO3 + SO42−) versus (Na+ + K+) − Cl−.
Ijerph 16 00886 g007
Figure 8. Gibbs plot showing major processes controlling groundwater chemistry.
Figure 8. Gibbs plot showing major processes controlling groundwater chemistry.
Ijerph 16 00886 g008
Figure 9. Water Quality Index map.
Figure 9. Water Quality Index map.
Ijerph 16 00886 g009
Figure 10. United States Salinity Laboratory (USSL) diagram for irrigation water classification (USSL 1954).
Figure 10. United States Salinity Laboratory (USSL) diagram for irrigation water classification (USSL 1954).
Ijerph 16 00886 g010
Figure 11. Plot of percentage sodium (%Na) versus EC (Wilcox 1948).
Figure 11. Plot of percentage sodium (%Na) versus EC (Wilcox 1948).
Ijerph 16 00886 g011
Table 1. Weight and relative weight of each chemical parameter.
Table 1. Weight and relative weight of each chemical parameter.
Chemical ParametersWHO Standard (mg/L) (Si)Weight (wi)Relative Weight (Wi)
pH6.5–8.530.097
TDS100050.161
Na20040.129
Mg15030.097
Ca20030.097
Cl25050.161
SO425050.161
HCO325010.032
K1220.065
Σ w i = 31 Σ W i = 1
Table 2. Summary of water quality indices for irrigation.
Table 2. Summary of water quality indices for irrigation.
IndicesAcronymFormula
Sodium absorption ratioSAR Na + / ( Ca 2 + + Mg 2 + ) / 2
Residual sodium carbonateRSC ( CO 3 2 + HCO 3 _ ) ( Ca 2 + + Mg 2 + )
Sodium percentage%Na [ ( Na + + K + ) / ( Na + + K + +   Ca 2 + + Mg 2 + ) ] × 100
Permeability indexPI [ ( Na + + HCO 3 _ ) / ( Ca 2 + + Mg 2 + + Na + ) ] × 100
Magnesium hazardMH [ ( Mg 2 + ) / ( Ca + 2 + Mg 2 + ) ] × 100
Kelly’s ratioKR [ ( Na + )   / ( Ca 2 + + Mg 2 + ) ]
Potential salinityPS Cl + SO 4 2
Table 3. Statistical analyses of chemical parameters (units of all parameters are mg/L, except pH, electrical conductivity (EC) µS/cm, Turbidity NTU and arsenic µg/L).
Table 3. Statistical analyses of chemical parameters (units of all parameters are mg/L, except pH, electrical conductivity (EC) µS/cm, Turbidity NTU and arsenic µg/L).
ParametersMinimumMaximumMeanStandard DeviationWHO StandardsNSBL *NSBL %
Turbidity01657.8227.2951016.95
pH6.708.207.410.336.5–8.500
EC24949501570.971061.9010003966.10
TDS1593168993.92677.5110002135.59
TH751270421.12252.983004067.80
Alkalinity0195.693.69- 0
Na+17638165.20151.792001830.51
K+0253.535.151246.78
Mg2+1017555.4134.3115011.69
Ca2+1422076.4649.6520035.08
Fe010.090.190.323.39
F020.420.501.546.78
Cl16779209.76186.322501932.20
SO42−20600169.88140.072501118.64
HCO359950308.90170.292503355.93
NO3-N09.900.941.601000
As025021.9548.31101525.42
* Number of samples beyond (WHO) limits.
Table 4. Correlation coefficient matrix of major physiochemical parameters.
Table 4. Correlation coefficient matrix of major physiochemical parameters.
ParameterpHECTDSTHNa+K+Mg2+Ca2+ClSO42−HCO3
pH1
EC−0.1761
TDS−0.1670.995 *1
TH−0.2270.891 *0.888 *1
Na+−0.0820.925 *0.923 *0.659 *1
K+0.0860.1400.1280.1520.1001
Mg2+−0.3010.847 *0.838 *0.971 *0.602 *0.1381
Ca2+−0.1660.863 *0.859 *0.977 *0.626 *0.1610.924 *1
Cl−0.0800.955 *0.942 *0.812 *0.917 *0.1300.771 *0.811 *1
SO42−−0.1900.899 *0.904 *0.790 *0.849 *0.0910.753 *0.730 *0.833 *1
HCO3−0.2150.858 *0.865 *0.830 *0.746 *0.1450.774 *0.800 *0.716 *0.659 *1
* Correlation is significant at the 0.01 level (2-tailed). Bold = strong correlation (>0.90).
Table 5. Classification of groundwater quality according to the Water Quality Index (WQI).
Table 5. Classification of groundwater quality according to the Water Quality Index (WQI).
WQIWater TypeNo. of SamplesPercentage of Samples
<50Excellent water915.25
50–100Good water3355.93
>100Unsuitable for drinking1728.82
Table 6. Statistical summary of irrigation quality indices of groundwater samples.
Table 6. Statistical summary of irrigation quality indices of groundwater samples.
IndicesMinimumMaximumMeanSDPermissible LimitUnsuitable Samples Suitable Samples %
SAR0.6516.213.410.37≤18-100
RSC−16.352.52−3.380.45≤2.5198.30
%Na15.6486.3241.701.98≤60591.52
PI31.1996.3259.671.76>25-100
MH33.9465.5346.830.75≤501869.49
KR0.186.310.890.11≤11771.19
PS1.3725.687.740.77≤101672.89

Share and Cite

MDPI and ACS Style

Talib, M.A.; Tang, Z.; Shahab, A.; Siddique, J.; Faheem, M.; Fatima, M. Hydrogeochemical Characterization and Suitability Assessment of Groundwater: A Case Study in Central Sindh, Pakistan. Int. J. Environ. Res. Public Health 2019, 16, 886. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16050886

AMA Style

Talib MA, Tang Z, Shahab A, Siddique J, Faheem M, Fatima M. Hydrogeochemical Characterization and Suitability Assessment of Groundwater: A Case Study in Central Sindh, Pakistan. International Journal of Environmental Research and Public Health. 2019; 16(5):886. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16050886

Chicago/Turabian Style

Talib, Muhammad Afnan, Zhonghua Tang, Asfandyar Shahab, Jamil Siddique, Muhammad Faheem, and Mehak Fatima. 2019. "Hydrogeochemical Characterization and Suitability Assessment of Groundwater: A Case Study in Central Sindh, Pakistan" International Journal of Environmental Research and Public Health 16, no. 5: 886. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph16050886

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