Hydrogeochemical Assessment of Groundwater and Suitability Analysis for Domestic and Agricultural Utility in Southern Punjab, Pakistan

: Groundwater is a critical water supply for safe drinking water, agriculture, and industry worldwide. In the Khanewal district of Punjab, Pakistan, groundwater has severely deteriorated during the last few decades due to environmental changes and anthropogenic activities. Therefore, 68 groundwater samples were collected and analyzed for their main ions and trace elements to investigate the suitability of groundwater sources for drinking and agricultural purposes. Principal component analysis (PCA) and cluster analysis (CA) were employed to determine the major factors inﬂuencing groundwater quality. To assess the groundwater’s appropriateness for drinking and irrigation, drinking and agricultural indices were used. The pH of the groundwater samples ranged from 6.9 to 9.2, indicating that the aquifers were slightly acidic to alkaline. The major cations were distributed as follows: Na + > Ca 2+ > Mg 2+ > K + . Meanwhile, the anions are distributed as follows: HCO 3 − > SO 42 − > Cl − > F − . The main hydrochemical facies were identiﬁed as a mixed type; however, a mixed magnesium, calcium, and chloride pattern was observed. The reverse ion exchange process helps in exchanging Na + with Ca 2+ and Mg 2+ ions in the groundwater system. Rock weathering processes, such as the dissolution of calcite, dolomite, and gypsum minerals, dominated the groundwater hydrochemistry. According to the Weight Arithmetic Water Quality Index (WAWQI), 50% of the water samples were unsafe for drinking. The Wilcox diagram, USSL diagram, and some other agricultural indices resulted in around 32% of the groundwater samples being unsuitable for irrigation purposes. The Khanewal’s groundwater quality was vulnerable due to geology and the inﬂuence of anthropogenic activities. For groundwater sustainability in Khanewal, management strategies and policies are required.


Introduction
Groundwater is the principal water supply source for household, agricultural, and industrial purposes in many countries [1]. Approximately one third of the world's population relies on groundwater for drinking [2]. In some arid and semi-arid regions, groundwater is the only source of drinking water [3,4]. Groundwater resources have been experiencing increasing problems in recent years, including quality issues and the depletion of freshwater resources, attributed to both anthropogenic and geogenic factors [5,6]. N-fertilizer, manure applications, and pesticides have enhanced crop yields, but they have a negative impact on the groundwater quality [7].
Recharge sources, geological structure, hydrological setting, the mineralogy of the watersheds and aquifers, and water-rock interaction, including mineral dissolution, ion exchange, redox, and anthropogenic activities, are all potential factors that affect groundwater quality [6]. Water quality is just as important as its quantity when meeting basic human needs, since a shortage of water and sanitation creates a slew of public health issues [8,9]. Contamination and overexploitation are responsible for the poor quality of groundwater. Understanding the significant factors influencing the qualitative and hydrochemical properties of accessible groundwater is thus crucial for maintaining a sufficient supply for various purposes [10,11]. Analysis and plots can be alternative to the isotopic studies for identifying the major factors that characterize the hydrogeochemical mechanism, such as multivariate analysis, USSL, silicate weathering, Gibbs plots and saturation indices. This analysis was carried out in some previous and recent studies [4,8,12,13].
The groundwater quality and hydrochemical features within a region are the basis for analyzing the utilization and availability of regional water resources in recent years [14]. Many hydrochemical components play an essential role in human health [15]. As a result, a detailed understanding of the hydrogeochemical processes is required to evaluate the suitability of groundwater for drinking, agriculture, and other applications. Pakistan is a water-stressed country. In the last few decades, both surface and groundwater availability have declined substantially in the surrounding regions of Pakistan [16][17][18], while the demand for domestic groundwater grew from 5.20 million acres to 9.70 million acres. The significant rise in groundwater quality problems is due to over-exploitation in many areas of Pakistan. Groundwater with elevated arsenic and fluoride concentrations is a serious health issue, while the research on the health effects is still limited. Multiple causes have enhanced Pakistan's groundwater vulnerability in recent decades. TDS levels of 400-20,000 ppm have been reported in entrapped severely mineralized groundwater [19]. The groundwater in Sindh province (Tharparker, Nagarparker, and Umarkot) is polluted with arsenic, resulting in skin injuries on major body parts (arms, face, and feet), rough skin with black dots from dermal and oral contact, and dental and skeletal fluorosis from arsenic-and fluoride-contaminated water [20]. Arsenic and chromium concentrations were significantly higher in the groundwater in Lahore, Muzaffargarh, and Mailsi [21]. As a result of the high fluoride concentration in the groundwater, bone softening and deformation problems have been identified in some areas of Punjab province. Furthermore, high fluoride levels in the groundwater have been reported in Lahore city [22]. In different parts of the world, several hydrogeochemical methods have been used to assess the suitability of groundwater for drinking and agricultural purposes. These methods have recommended various measures to increase water quality [23].
The purpose of this study is to determine the hydrochemistry and suitability of groundwater in the Khanewal district for drinking and irrigation. This study will help us better understand the hydrogeochemical processes that cause ions to form in groundwater, as well as the factors that influence the groundwater quality in the area.

Study Area
Khanewal is a district of Punjab province in Pakistan situated in the lower Bari of the river Doab between the Sutlej and Ravi rivers (from 29 • 52 03 N to 30 • 44 30 N and from 71 • 32 19 E to 72 • 28 31 E) with an altitude of 128 m, as shown in Figure 1. According to the Pakistani census of 1998, the district had 2,068,490 people, with 17.42 % of them living in urban areas [12]. Four seasons-summer, winter, autumn, and spring-make up the region's climate. The temperature ranges from 21 • C to 51 • C during the summer season and from 1 • C to 27 • C during the winter season [12]. Khanewal is 134 m above sea level and has a desert climate, with annual precipitation of about 166 mm. Furthermore, the Köppen-Geiger climate classification classifies the area's climate as BWh.

Geology and Hydrogeology of the Study Area
Khanewal is located in the upper Indus plains, and the river in the Khanewal district has significant impact on the basis of the physical qualities in the region. The area's geology is characterized by quaternary loess deposits, quaternary flood plain deposits, and quaternary flood plain deposits (Figure 1). Silt and sand make up the loess deposits. The silt concentration, on the other hand, is estimated to be 60-90%. Loess is wind-blown sediment that preserves a complete record of interglacial glacial cycles on land. Loess ranges in thickness from a few centimeters to more than 91 m. Khanewal loess deposits are heterogeneous sediments that contain buried soils. On the other hand, the quaternary loess deposits are made up of mica, feldspars, quartz, and calcium carbonate. The alluvium overlies semi-consolidated Tertiary rocks or Precambrian age metamorphic igneous rocks. The recent Pleistocene alluvial complex contains unconsolidated silt, sand with gravel, and minor clay. These sediments from the Ravi River and its tributaries make up the soils. The alluvial deposit thickness along major streams ranges up to 30.48 m and probably averages 12.192 m. Along minor streams, the thickness ranges up to 13.716 m and probably averages 6.096 m [12]. The individual stratum's vertical and horizontal continuity is limited in diverse alluvial deposits. The soil is primarily reddish brown to grayish brown in color, medium-textured, and contains a lot of fine to very fine sand and very small amounts of clay and gravel [24]. The majority of the alluvial complex comprises fine to medium sand, silt, and clay. In silty or clayey sand, siltstone and mudstone pebbles can be found. The alluvial deposits form a unified, highly transmissive non-artesian aquifer despite their heterogeneity. For the most part, the transmissivity ranges from 954 to 4263 m 3 /day/m [25]. Above 300 feet, the strata are compacted but exceedingly productive. The water table depth in the area varies between 15 m and 95 m, with an average of 48.2 m. The water level contour map was created by calculating the individual hydraulic heads (Figure 1).

Sampling and Analysis
In the Khanewal district, groundwater samples (n = 68) were collected from different groundwater sources such as wells, boreholes, and hand pumps to investigate the hydrochemical properties and suitability for drinking and agriculture. Groundwater samples were collected in 1.5-L polyethylene bottles and washed with a 10% HCl solution. To collect samples from the area, a duplicate sampling technique was used. The basic water quality parameters like pH, water temperature, total dissolved solids, and electric conductivity were measured in situ after sampling. To preserve the groundwater samples, 2-3 drops of HNO 3 were added and acidified with a sequence of duplicate samples. The pH, electrical conductivity (EC), and total dissolved solids (TDS) were calculated using a pH meter and an electrochemical analyzer (Hac 44600-00, Loveland, CO, USA). Significant anions such as nitrate (NO 3 -N) and sulfate (SO 4 2− ) were measured using an ultraviolet-visible (UV-VIS) spectrophotometer (Analytik Jena, Jena, Germany), while the chloride (Cl − ), total alkalinity, and bicarbonates (HCO 3 − ) were measured using titration methods. The iron (Fe 2+ ) and major cations viz. Ca 2+ , Mg 2+ , Na + , and K + were analyzed using an atomic absorption spectrophotometer (Varian Spectra AA-240) under standardized operating conditions. The charge balance error (CBE) was calculated to confirm the authenticity of the water sample analysis. A positive CBE is associated with higher cation concentrations in these water samples, while a negative CBE is associated with higher anion concentrations [26]. The ionic levels were measured in milliequivalents per liter (meq/L) in this study. The CBE% for all groundwater samples was within ±5%, indicating that the analytical data were reliable [27].

Correlation Analysis
Correlation analysis is a statistical tool for determining how closely two variables are related. Pearson correlation analysis was used to distinguish the link between several water chemical parameters in this study.

Hierarchical Cluster Analysis
Hierarchical clustering is a clustering approach that categorizes distinct datasets based on their similarity, obtained by computing the distance between each category's data points and all of the data points. The Euclidean distance was used in this study to compute the distance between various types of data points. Variance analysis was then used to examine the distances between distinct clusters using Ward's linkage approach [8,28]. The final result is displayed as a dendrogram. In addition, the number of clusters is governed by the phenom line.

Principal Component Analysis
Principal component analysis (PCA) is a popular method for reducing data dimensionality. PCA can keep the majority of the information while reducing the data's dimensionality [8,29]. In this study, principal component analysis (PCA) was used to determine whether the hierarchical cluster analysis was reasonable. SPSS Statistics (Version 20.0) has been used to find the PCA.

Geochemical Modeling
The tendency of groundwater to dissolve or precipitate a particular mineral can be determined using saturation indices [30]. In addition, the measurement of the mineral balance is helpful in predicting the dissolved mineral reactivity in water [31]. Saturation indices have been computed using a geochemical simulation program called PHREEQC Interactive (version 3.4).

Water Type
A Durov diagram is commonly used to represent the hydrogeochemical type and relative content of main anions and cations in distinct samples, and it can also show some potential geochemical processes that might assist in the understanding and assessment of groundwater quality [32]. The Durov diagram was prepared with the help of Grapher (version 14).

Groundwater Classification
The Gibbs diagram is a common tool to determine the relationship between lithology and hydrochemistry in aquifers [8]. According to this diagram, the dissolved hydrochemical compositions are divided into three distinct areas. Both the cation and anion ratios viz. TDS were plotted against Na + / (Na + + Ca 2+ ) and Cl − / (Cl − + HCO 3 − ). We used Microsoft Excel 2020 to draw the Gibbs diagram.

Ion Exchange Processes
The chloro-alkaline indices (CAI1 and CAI2) were used to analyze the ion exchange reaction [33]. These indices have a significant influence on the chemistry and evolution of water [34]. The CAI1 and CAI2 were calculated using the formulas below: All units are in milliequivalents per liter (meq/L). When both indices had positive values, direct ion exchange occurred (Equation (3)). When CAI1 and CAI2 had negative levels, reverse ion exchange occurred (Equation (4)):

Water Quality Analysis for Drinking
Water Quality Index (WQI) The weighted arithmetic water quality index (WAWQI) values were calculated to assess the groundwater's suitability for drinking purposes [35]. The WHO (2011) standard for drinking water was used to calculate the WAWQI. The WAWQI was computed by using Equation (5):

Groundwater Quality Evaluation for Irrigation
Understanding the salinity levels of water, which affect the soil structure and agricultural output, is the primary goal of determining the suitability of water for irrigation. To determine the qualities of irrigation water, the following criteria were applied.
The formula for its calculation is as follows:

Kelly's Ratio
Kelly's ratio is a crucial requirement for determining the quality of groundwater used in agriculture [40]:

Magnesium Hazard
Magnesium hazard analysis is another approach suggested by Szabolcs and Darab for determining the agricultural suitability of groundwater [41]. Soil modification increases as the Mg 2+ content relative to Ca 2+ grows, allowing clay particles to scatter and causing damage to the soil structure by reducing hydraulic conductivity:

Permeability Index
The Permeability Index (PI) is a primary factor for assessing water's appropriateness for irrigation, according to [42]. Class 1 (>75%, suitable), class 2 (25-75%, good), and class 3 (less than 25%, unsuitable) are the three classifications for PI: Additionally, Grapher (version 14) software was used to create a U.S. salinity laboratory diagram and a Wilcox diagram by linking the sodium absorption ratio and electrical conductivity and the percentage of sodium and electrical conductivity, respectively, to define and analyze the suitability of groundwater for irrigation. Table 1 shows the groundwater statistical results compared with the World Health Organization guidelines. The pH ranged from 6.9 to 9.2 with a mean value of 7.79, which suggests that the nature of the groundwater was slightly acidic to alkaline. However, 66 groundwater samples with 97% contributions were within the WHO's recommended guidelines [43]. The pH value of the groundwater has no direct impact on human health. However, groundwater in an acidic medium (pH = 6) may accelerate heavy metal solubility and leaching, whereas a higher pH is expected to boost calcite mineral precipitation [44]. The alkalinity of the research area was high due to weathering and anthropogenic activity. Bicarbonate, carbonate, and hydroxide are the most well-known alkalinity components. The breakdown of minerals produces these chemicals [45]. The turbidity value varied between 0 and 7 Nephelometric Turbidity Units (NTUs) with an average value of 2.52 NTUs, and only three samples showed high turbidity. The higher turbidity in the groundwater system of the area belonged to a shallow-depth well due to contamination exposure compared with the deep aquifer and mid-depth aquifer. Thus, surface recharge, water runoff, weathering processes, and industrial effluents can further increase the turbidity in the water system of the area. Moreover, high turbidity also results from shallow and poorly constructed wells that can cause high turbidity [29]. The most acceptable EC limit in drinking water is 400 µS/cm (WHO, 2006). The EC values ranged between 350 µS/cm and 1525 µS/cm with a mean value of 783.6 µS/cm. The groundwater ion concentrations were directly associated with the EC, resulting in higher salinity and dissolved concentrations. However, 92.6% of the samples exceeded the permissible limits. The TDS levels ranged from 137.0 to 1414 mg/L in the study area with a mean of 707.3 mg/L. Water with more than 1000 mg/L dissolved solids typically gives an unpleasant taste or renders the water unsuitable. The high TDS value affects the water's taste, hardness, and corrosive properties [46]. The maximum permissible guideline value of TDS recommended is 1000 mg/L [43].

Groundwater Quality and Hydrochemical Characteristics
The total hardness varied from 35.0 to 530 mg/L with an average of 215.7 mg/L. The water could be classified as very hard if the TH was more than 300.0 mg/L [47]. Ten samples from the study area were under the very hard category. The concentrations of the cations could be written as Na + > Ca 2+ > Mg 2+ > K + > Fe 2+ . Similarly, the concentrations of the anions were in the order of HCO 3 -> Cl -> SO 4 2-> F -> NO 3 -. The groundwater sodium (Na+) concentrations ranged from 15.0 mg/L to 324.0 mg/L with a mean of 90.87 mg/L. The maximum permissible limit for Na + is 200 mg/L, according to [43]. The most common causes of elevated levels of Na + in groundwater are salt deposit erosion and Na + -bearing rock minerals. Natural brackish water from some aquifers is also absorbed by halite mineral water [4,48]. Agricultural byproducts may be another source of sodium in groundwater in the study region [8]. The concentrations of Ca 2+ and Mg 2+ were found to be 16-34 mg/L and 1-73 mg/L, respectively, with mean values of 46.78 mg/L and 20.78 mg/L, respectively. The permissible limit for Ca 2+ is 100 mg/L, and for Mg 2+ , it is 50 mg/L. According to the guidelines in [49], four samples had increased Ca 2+ levels, and three samples had increased Mg 2+ levels. Figure 2a shows carbonate weathering releasing Ca ++ and Mg 2+ , and Figure 2b reflects the Na + K against the total cations. Potassium (K + ) concentrations in the groundwater ranged from 0.2 to 34 mg/L, with an average of 6.42 mg/L. According to [43], a maximum K + concentration of 12 mg/L is considered acceptable in the groundwater. In the study area, six samples exceeded the permissible limits. K + is produced in groundwater by weathering potassium-bearing rocks such as sylvite (KCl) and silicates, especially clay minerals [4]. K + may also be added to groundwater by fertilizers and the breakdown of waste materials [50]. The maximum permissible concentration of iron allowed in the groundwater according to [43] is 0.3 mg/L. The higher Fe 2+ concentration in the groundwater system was due to the ferruginous minerals of the Earth's crust.
The most dominant anion in the study area was HCO 3 − , with concentrations that varied from 100 to 510 mg/L and a mean of 243.7 mg/L. In our study area, 38% of the samples exceeded the permissible limit. The presence of organic matter in the aquifer that is oxidized to produce carbon dioxide, promoting mineral dissolution, results in the release of bicarbonate in the groundwater of the study area [51]. SO 4 2− concentrations ranged from 22.0 to 498.0 mg/L with an average of 117.55 mg/L. The permissible limit for SO 4 2− is 250 mg/L, and 5.9% of the samples were found above the allowable limit of the WHO. In addition to water-rock interaction, a high SO 4 2− content may also be contributed by anthropogenic sources in the study region [8,51].
The Cl − ion is an essential element found in most groundwater in the world. Groundwater Cl − depends on many variables, such as weathering, sedimentary leaching of rocks, and anthropogenic sources [8,52]. The Cl − content in the groundwater of the study area varied from 9.6 to 167 mg/L with a mean of 46.7 mg/L. The nitrates in the groundwater were between 0 and 4.2 mg/L with a mean value of 0.30 mg/L. Anthropogenic sources cause the pollution of NO 3 -N [53].
The F − concentration in the groundwater of the study area ranged from 0.01 to 3.95 mg/L, with a mean of 0.54 mg/L. F − is primarily attributed to the oxidation of minerals and the weathering of granite and igneous rocks. Fluoride is a vital component of the human body, and the fluoride content of potable water should not reach 1.5 mg/L [43]. Long-term exposure to fluoride-contaminated water causes dental and skeletal fluorosis [54,55]. In the KL-3 and KL-7 groundwater samples, the fluoride concentration was beyond the permitted limits. Fluoride concentrations are induced by the dissolution of fluorides comprising minerals from granite and gneiss rock [56].
Cation exchange mechanisms with greater alkaline pH, HCO 3 − , and lower Ca 2+ levels in groundwater increase F − pollution [55]. The formation of higher F − in groundwater is influenced by water-rock interactions and dry weather conditions. In acidic circumstances, clay minerals adsorbed F − , which was then desorbed in alkaline pH.
The Durov diagram is a useful graphical tool for representing hydrogeochemical data [57] (see Figure 3). In the groundwater, the mixed type was the most prevalent. A few samples were plotted in the D (sodium form) field, indicating that Na + K predominated in the groundwater. On the other hand, there was a noticeable difference in the groundwater anions. The E field samples belonged to the bicarbonate type, suggesting the HCO 3 − dominant anion. Those samples were primarily characterized by freshwater (TDS < 1000 mg/L) located in the aquifer recharge region. The other belonged to brackish water's sulfate and mixed water types, suggesting the combined effects of evaporation, water-rock contact, and human activities.

Silicate Weathering
The majority of the samples were above the aquiline of 1:1 in the scatter diagram of (Ca + Mg) vs. (HCO 3 + SO 4 ), indicating that carbonate weathering was a primary source of Ca 2+ and Mg 2+ [8]. The samples fell along the aquiline 1:1, implying that these ions resulted from carbonate and silicate weathering [58]. Carbonate weathering is responsible for ionic concentrations above the aquiline, while silicate weathering is responsible for those below the aquiline [59]. HCO 3 − is mainly derived from CO 2 in the soil zone during the weathering of source rock minerals. The soil zone in the subsurface environment has a higher CO 2 pressure due to the decomposition of organic matter and root respiration, which reacts with rainfall (H 2 O) to produce HCO 3 − [60]. However, the sample data would also plot above the 1:1 line of a Na + + K + vs. total cation graph, as silicate weathering is involved in the geochemical process. The data plot above the line confirms that silicate mineral weathering contributed to the groundwater chemistry ( Figure 4a).
However, evaporite dissolution, silicate, carbonate rock weathering, and the effect of rock weathering on the groundwater chemistry were investigated using Na-normalized molar ratios (Figure 4b). As shown in Figure 4b, most of the samples were found in the bottom left and middle sections, indicating that evaporite dissolution and silicate weathering were dominant in the formation of the groundwater chemistry. The dissolved solutes were regulated mainly by carbonate weathering, as the distributed sample in the middle and upper right sections illustrates. In the study area, the sodium is released into the groundwater by halite and silicate minerals. If the mole ratio of Na/Cl is <1, then halite dissolution may have occurred, or a mole ratio of Na/Cl > 1 may indicate silicate dissolution [8]. In igneous rock, silicate minerals mostly occurred, as 95% of the Earth's crust is composed of silicate minerals and the mantle. Silicon and oxygen are the primary components of silicate. The known minerals of silicon include quartz, mica, feldspar, pyroxene, olivine, amphibole, and clay minerals. In comparison, the common minerals include quartz, tridymite, coesite, cristobalite, keatite, and lechatelierite. The composition of the groundwater in our analysis showed that 85.3% of the samples had a Na/Cl molar ratio of ≥1, and 14.7% showed Na/Cl < 1. Therefore, the molar ratio of Na/Cl was ≥1, denoting silicate weathering in the groundwater samples of our analysis, and the ratio less than 1 showed non-silicate weathering. In addition, the groundwater Na/Cl ratio for this study varied between 0.45-8.1 and 2.16, respectively. However, the silicate minerals' range and average results were observed (1.03-8.10 and 2.42), and those for the non-silicate minerals were identified (0.45-0.93 and 0.70). As a result, the findings of this analysis indicate that both halite and silicate minerals caused groundwater contamination in the study region. Therefore, these findings are consistent with the studies mentioned above that were carried out in Pakistan [4,8,48].

Hydrogeochemical Evolutional Processes
The Gibbs diagram depicts three variables that influence the groundwater chemistry, including the dominance of evaporation, precipitation dominance, and weathering dominance [8,55,61]. According to Figure 5, most samples were plotted in the rock dominance region, indicating that rock dominance influenced most of the groundwater samples. The main contributor to mineral groundwater enrichment was the weathering of rock. The process by which the groundwater mixed dissolvable salts and minerals promoted the weathering of the parent rock. In addition, a long rock-water contact residence time also allows mineral dissolution [8]. Therefore, aquifer lithology and groundwater bedrock mineralogy have been highlighted as geochemical processes of importance in the study area.  Table 2 shows the correlation matrix of the physiochemical parameters of the study area for the groundwater samples. According to the Pearson correlation matrix, it was found that TDS had strong correlation with Na + (0.83), HCO 3 − (0.74), Cl − (0.60), and SO 4 2− (0.86). All ionic concentrations increased with the increase in the TDS value. These ions were indicated to be primarily intended for sedimentary rock weathering. A poor correlation between Na + and Cl − may have been due to their genetic association with low evaporation. A low to moderate correlation between SO 4 2− and Ca 2+ suggests that gypsum dissolution did not occur, but the excess of Ca 2+ was due to carbonate and the cation exchange [62], while a positive association between Ca 2+ and K + could be attributed to the ion exchange and clay minerals.

Mineral Phases
Measurements of the saturation could be used to estimate the subsurface mineralogy. Some minerals (calcite and dolomite) are typically present in equilibrium in subsurface and surface waters. SI calculations have shown in the current analysis that the carbonate minerals had varying degrees of saturation. The SI of the dolomite ranged from −1.19 to +2.66, as shown in Figure 6, and that of calcite ranged from 0.57 to 1.54. By contrast, the SI for gypsum ranged from −2.69 to −0.92. These findings indicate that the dolomite and calcite showed variation in different samples, with both supersaturation and unsaturation levels. These minerals may precipitate (SI > 1), and further dissolution (if SI < 0) is also possible. On the other hand, the Gypsum Value was below 0, indicating an unsaturation condition, and the findings also mean that Na + and Cl − may not be the primary source of halite. This research shows that silicate minerals led to groundwater pollution in the study region [63].

Ion Exchange Processes
The interaction between the Na + K-Cl and Ca + Mg cations is widely used to identify cation exchange processes (HCO 3 -SO 4 ) [64]. Furthermore, the chlor-alkali index (CAI) can be used to determine whether the cation exchange is reverse or direct [8]. If the CAI value is greater than 0, it indicates that the groundwater Ca 2+ was exchanged for Na + in the aquifer, whereas a CAI value >0 indicates reverse cation exchange. Positive CAI-1 and CAI-2 results suggest that Na + and K + ions are exchanged for Mg 2+ and Ca 2+ ions in water. Furthermore, when the Mg 2+ and Ca 2+ ion indices are negative, they are swapped with Na + and K + from rocks. This is an indication of an alkaline-chloride imbalance. The majority of the samples are shown in the lower left corner of Figure 7, showing reverse ion exchange. However, a small number of samples fell in the upper right corner, indicating direct ion exchange. The dominance of HCO 3 − +SO 4 2− vs. Ca 2+ + Mg 2+ could be related to cation exchange and silicate weathering. Only a few samples were discovered above the 1:1 line, demonstrating that reverse cation exchange is important in groundwater. Similarly, if cation exchange is the primary activity in groundwater, the linear relationship of Na + + K + Cl vs. (Ca 2+ + Mg 2+ ) − (HCO 3 − + SO 4 2− ) has a slope of −1 [65]. The linear fitting line had a slope of −1.002 (Figure 8a), being close to −1, indicating cation exchange in the groundwater system. In addition, in the groundwater, reverse cation exchange reactions should occur.
The lithogenic Na + in the groundwater was presumed and available from the meteoric origin, which was balanced by equivalent Cl − , and the equivalent Cl − was subtracted from Na + [4]. Cation exchange mechanisms control the ionic compositions of groundwater. The R 2 value of Na + against (Na + Cl − ) meq/L was 0.99 (Figure 8b), favoring Ca ++ and Mg ++ precipitation during ion exchange processes. The surface of the clay mineral absorbed Na + and mobilized F − in the groundwater. The results of the ion exchange process are relatively supported by the previous studies [8,66].

Cluster Analysis
Cluster analysis (CA) is a critical multivariate statistical technique often used to organize complex water data into various clusters better. In Ward's method, the squared Euclidean distance between the water variable was calculated for the similarity index. CA combined water observations that fell within the comparable cluster with the most similar geochemical results [8,67]. Thus, three cluster groups were constructed after Ward's method of measuring the Euclidean distance. C1 represents less polluted, C2 represents moderately polluted, and C3 indicates severely polluted.

Principal Component Analysis
The principal component analysis (PCA) results for the groundwater parameters are shown in Table 3. The PCA results were obtained after varimax rotation to elucidate the acquired bases that influenced the groundwater [8,67,68]. For the groundwater parameters, F1, F2, F3, and F4 were derived with eigenvalues of 5.86, 2.20, 1.48, and 1.37, respectively, resulting in a total variance of 36.60 %, 13.75%, 9.24%, and 8.54%, respectively.
The first factor (F1) represented 36.60% of the total variance with an eigenvalue of 5.86 as shown in the Table 3 indicates the ionic configuration of groundwater and resultant from mineral dissolution, ion exchange, and weathering of the host granitic rocks. The dissolution of gypsum-and calcium-carrying rocks, as well as the interaction of water and rock, initiate the source of SO 4 2− and HCO 3 − . F1 is also linked to water salinity, as determined by EC and major ion chemistry. Natural processes, such as carbonate (CaCO 3 ) and mineral bearing, contributed significantly to the F1 component. High levels of major ions like Ca 2+ , Mg 2+ , Na + , and Cl − may reflect anthropogenic sources which affect mineral water reactions [4,8]. The second factor (F2) accounted for 13.75% of the overall variance, with an eigenvalue of 1.84. The strong correlation coefficient for Mg 2+ and Ca 2+ , with correlation coefficients of 0.53 and 0.66, respectively, were observed, with water forms of HCO 3 − Ca and Mg having a percentage contribution of component F2. The high concentrations of Mg 2+ and Ca 2+ might have been due to extreme evaporation, scarce precipitation, and a shallow water table influenced by irrigation water [13].
F3 and F4 accounted for 9.24% and 8.54% of the variability with eigenvalues of 1.48 and 1.37, respectively. The PCA results support the hydrogeochemical processes, and F3 and F4 had a significant correlation with NO 3 − , pH, and HCO 3 − , with coefficient (r) values of 0.64, 0.51, and 0.63, respectively. Both natural and anthropogenic factors caused groundwater pollution in the study area.

Groundwater Suitability Evaluation for Drinking Purposes
The Water Quality Index (WQI) is a vital groundwater quality parameter used to evaluate groundwater's aptness for drinking [69]. The WQI is a method that calculates the effect of the groundwater parameters on the groundwater quality, and it is measured by using the WHO's drinking standards [43]. There are five categories according to the WQI criteria: excellent (<50), good (>50), bad (>100), very poor (>200), and water not suitable for drinking (>300) [70]. Table 4 shows the WQI classification of the study area, listing the different types of groundwater based on the WQI references. It depicts more than 50% of the groundwater samples as not being safe for drinking purposes. A spatial distribution map of Figure 10 provided more information about the different types of water in the study region, where 2.94% of the samples were excellent for drinking, while 39.7% were poor, 11.76% were good, 27.94% were very poor, and 17.64% were not suitable for drinking. Furthermore, Figure 10 demonstrates that the majority of the samples in the groundwater of the study area were inappropriate for drinking purposes.

Assessment for Irrigation Purposes
The suitability of the groundwater in the study area for irrigation was assessed by using Wilcox and USSL diagrams, as well as other agricultural indices such as the sodium percentage (Na%), sodium adsorption ratio (SAR), Permeability Index (PI), magnesium hazard (MH), and Kelly ratio (KR). Table 5 summarizes the irrigation indices of the groundwater samples statistically. According to SAR, the sodium hazard in the groundwater could be rated as low (SAR10), medium (10)(11)(12)(13)(14)(15)(16)(17)(18), moderate (18)(19)(20)(21)(22)(23)(24)(25)(26), or very high (SAR > 26). The SAR ranged from 1.65 to 65.19 in the study area, with an average of 12.87, indicating a low, medium, and high risk of sodium. The USSL diagram (Figure 11a) shows 45 samples belonging to the C 2 S 1 and C 3 S 1 groups, indicating that they could be used for irrigation with minimal Na + exchange. Nine samples belonged to the C 3 S 2 group, four to the C 3 S 4 group, five to the C 2 S 2 group, and five to the C 1 S 1 group. These types of water are not suitable for agricultural use.  The Na + % varied from 10.64% to 89.78% with a mean value of 45.66%, and 77.94 % of the samples were suitable for agricultural purposes. A high percentage of Na + causes deflocculation, soil tilt, and permeability. Groundwater containing less than 60% Na+ is considered acceptable for irrigation usage. Wilcox's diagram further defines groundwater for irrigation, where EC is graphed against the percentage of Na + used. As in Figure 11b, approximately 36.7% of the samples were good to permissible, followed by 22% percent belonging to doubtful to permissible, and 22% of the samples fell into doubtful categories for irrigation use, while 11.7% were unfit for irrigation use.
The study area findings were in the range of 34.7-70.3, with a mean value of 70.3. Most samples fell into class 1 under the Deneen classification sampling programs, suggesting that the groundwater was suitable for irrigation. A magnesium hazard less than 50 is considered acceptable for irrigation purposes, while one greater than 50 is unfit for irrigation, and the crop yield declines as soils become more alkaline [4,8]. With a mean value of 38.42, the magnesium value for the study area varies from 9.4 to 68.95. In 13 samples, the magnesium hazard was more than 50, so these were unsuitable for irrigation, while others were good for irrigation.
The Kelly ratio ranged from 0.12 to 8.6 with a mean value of 1.29. If the Kelly ratio is less than one, it is considered ideal for irrigation purposes [71]. Kelly's ratio was lower than one in 39 samples, which means that the water sample was suitable for irrigation, whereas the 29 samples with a ratio greater than 1 were not suitable for irrigation in the study area.

Conclusions
For the current study, a combination of statistical analysis, graphical techniques, and hydrochemical modeling was utilized to gain a full understanding of the groundwater chemistry and its suitability for drinking and agricultural uses in the designated research region. Statistical analysis shows that cations were abundant in the following order: Na + > Ca 2+ > Mg 2+ > K + > Fe 2+ . Meanwhile, anions were abundant in the following order: HCO 3 -> Cl -> SO 4 2-> F -> NO 3 -. The high concentrations of sodium and bicarbonate ions in the research illustrate the process of silicate weathering, which regulates the chemistry of key ions in groundwater. Furthermore, ion exchange and rock weathering processes, such as the dissolution of albite, carbonate, and gypsum minerals, control the rock-water interactions that regulate the groundwater hydrochemistry. The Gibbs diagrams showed rock weathering as the primary source of the hydrochemical compositions of the groundwater.
Furthermore, the Durov diagram revealed no dominating ions, although some cation and anion samples revealed dominance in the Na + -K + type, and some samples revealed observable chemistry of the types HCO 3 and Cl -. The geochemical modeling showed undersaturation of the evaporite mineral gypsum and supersaturation of carbonate minerals calcite and dolomite. The ion exchange method showed that reverse ion exchange was dominant. Multivariate analysis statistically proved the relation between different variables and their hydrogeochemistry in the groundwater. The WQI showed that 50% of the groundwater samples were unfit for drinking purposes. However, according to the agricultural indices, approximately 32% of the groundwater samples were unsuitable for irrigation. Anthropogenic activities and sedimentary rock weathering play a significant role in the depletion of the groundwater system in the study area. It is highly recommended that the local government and higher authorities control groundwater pollution in the surrounding vicinity of the Khanewal district.