Hydrogeochemical Investigation of Elevated Arsenic Based on Entropy Modeling, in the Aquifers of District Sanghar, Sindh, Pakistan

Arsenic (As) contamination in drinking groundwater is a common environmental problem in Pakistan. Therefore, sixty-one groundwater samples were collected from various groundwater sources in District Sanghar, Sindh province, Pakistan, to understand the geochemical behavior of elevated As in groundwater. Statistical summary showed the cations and anions abundance in decreasing order of Na+ > Ca2+ > Mg2+ > K+, and HCO3− > Cl− > SO42− > NO3−. Arsenic was found with low to high concentration levels ranging from 5 µg to 25 µg/L with a mean value of 12.9 µg/L. A major water type of groundwater samples was mixed with NaCl and CaHCO3 type, interpreting the hydrochemical behavior of rock–water interaction. Principal component analysis (PCA) showed the mixed anthropogenic and natural sources of contamination in the study area. Moreover, rock weathering and exchange of ions controlled the hydrochemistry. Chloro-alkaline indices revealed the dominance of the reverse ion exchange mechanism in the region. The entropy water quality index (EWQI) exposed that 17 samples represent poor water, and 11 samples are not suitable for drinking.


Introduction
Arsenic (As) is a toxic metalloid existing in groundwater from commonly anthropogenic and geogenic sources [1,2]. All over the world, more than 150 million people are exposed to excessively high As contamination of drinking water sources [3,4]. The prevalent range of health hazards of As in groundwater is a serious threat to Pakistan [5]. In several regions of Pakistan, heavy metal concentrations are detected, including Peshawar, Karachi, Hyderabad, DG Khan, Lahore, Malakand, Lower Dir, Muzaffargarh, and the Chitral district [6][7][8][9][10][11][12][13][14]. Arsenic is commonly dispersed in Sindh province in different regions such as Manchar lake, Jamshoro, Hyderabad, Jaccobabad, and Nagar Parker [15][16][17][18]. The carcinogenic consequence of As may lead to kidney failure, lung infections, hair loss, skin The groundwater samples were collected from district Sanghar, Sindh, province of Pakistan. The study area lies between 25° and 30° North latitude and 70-130° East longitude ( Figure 1). Sanghar district covers a total area of 987,313 hectares. The growing crops of the study area are rice, sugarcane, cotton, wheat, groundnuts, barley, and vegetables. The climate of the research area is colder in winter, hot, and drier in the summer season. There is a 2.2 million population in the study area. The primary occupation of the people of the Sanghar district is agriculture practices. Agriculture action is the main occupation, and more than 30% of people are engaged in this occupation to fulfill their daily needs. The second dominant income source of the area is unskilled labor that employs 22% of the local people to earn their livelihood.

Geology and Hydrology
The study area has rough topography from north to south. The Nara Canal irrigates the western section of the site is the leading water source for irrigation and drinking purposes. The study zone is bounded by alluvial pledges, primarily composed of silts, sand, and clay of Tertiary rocks [34]. The nature of soil differs significantly from one place to another sand and clay extent all over the study area. The research area is protected mainly by floodplains which consist of bar deposits. According to lithological studies, there are 250 to 450 feet thick sandy layers beneath the shallow aquifer [35]. They are satisfactory to intermediate micaceous sands having well-sorted bands and lenses of silt and clay. The Sanghar district covers a total area of 987,313 hectares. The growing crops of the study area are rice, sugarcane, cotton, wheat, groundnuts, barley, and vegetables. The climate of the research area is colder in winter, hot, and drier in the summer season. There is a 2.2 million population in the study area. The primary occupation of the people of the Sanghar district is agriculture practices. Agriculture action is the main occupation, and more than 30% of people are engaged in this occupation to fulfill their daily needs. The second dominant income source of the area is unskilled labor that employs 22% of the local people to earn their livelihood.

Geology and Hydrology
The study area has rough topography from north to south. The Nara Canal irrigates the western section of the site is the leading water source for irrigation and drinking purposes. The study zone is bounded by alluvial pledges, primarily composed of silts, sand, and clay of Tertiary rocks [34]. The nature of soil differs significantly from one place to another sand and clay extent all over the study area. The research area is protected mainly by floodplains which consist of bar deposits. According to lithological studies, there are 250 to 450 feet thick sandy layers beneath the shallow aquifer [35]. They are satisfactory to intermediate micaceous sands having well-sorted bands and lenses of silt and clay. The pledges are highly diverse, and the lithology varied so widely that the situation is challenging to match the strata discovered in two adjacent wells [36,37]. The proportion of sand and clay bands, on the other hand, is highly consistent over huge areas. In this aquifer, the depth of water tables ranges from 1.39 m to 12.76 m, having an average depth of 3.93 m [38]. The leading causes of salinity in the region include the under-watering of crops, low harvesting strength, and outflow from streams and adjacent canals. Chemical investigation shows that the studied area is modest to severe saltish, with saline-alkali soil accounting for 5.5% of the overall landmass.

Collection and Analysis of Samples
Groundwater samples (n = 61) were collected from various sources, including dug wells, boreholes, hand pumps, and tube wells, to evaluate the groundwater quality in Sanghar district, Sindh province ( Figure 1). The wells were pumped for more than 5 min before sampling to escape the effect of stagnant water [39]. Groundwater samples were then kept in polyethylene 1.5 L contamination-free bottles and immediately shifted to the research laboratory. The samples were subsequently tested in the standard water quality research laboratory of the Pakistan Council for Water Resources Research (PCRWR). The basic water quality parameters like TDS, EC, and pH were measured using an electrochemical analyzer and pH meter (Hac 44600-00, Loveland, CO, USA) [40]. An ultraviolet-visible (UV-VIS) spectrophotometer (Germany) was used to examine the samples for significant anions such as No 3 − and SO 4 2− . The titration method was used to analyze bicarbonate (HCO 3 − ) and chloride (Cl − ). A flame photometer (PFP7, Cambridge, UK.) was used to measure the main cations such as Na + , K + . To achieve the value for Ca 2+ and Mg 2+ , a volumetric titration with ethylene diamine tetra acetic acid (EDTA, 0.05 N) with <2% analytical error was used [41]. An atomic absorption spectrophotometer (AAS Vario 6, Analytik Jena, Jena, Germany) was used to measure the quantity of arsenic in groundwater samples. The charge balance error (CBE) was applied to verify the structures of groundwater analysis. The CBE is negative in water samples with a high concentration of anions while showing positive CBE; when cations, concentration increases. It was calculated using the following (Equation (1)).
Ionic absorptions are expressed in milliequivalent per liter (meq/L). Only those samples with less than ±5% CBE were accepted for further analysis following a standard protocol [7,42].

Assessment of Water Quality Entropy Water Quality Index (EWQI)
The water quality index (WQI) is a simple operative measure for defining the suitable drinking water used worldwide [43]. The old-style approaches of WQI suffer from some attentive limitations and could not deliver important evidence on groundwater quality. This study uses an advanced tool to quantify groundwater quality using entropy weight (EWQI) [44]. EWQI was measured via the following equation: where n denotes the sum of factors used to calculate the EWQI, and (ωj) indicates the jth A = P 1 + r n nt .
The parameter of entropy weight and q j shows the jth parameter quality rating scale. Equation (3) was used to estimate the entropy weight (ωj) for individual factors: Water 2021, 13, 3477 5 of 20 The following connections express the information of entropy (e j ): where (m) specifies the whole sum of samples, and p ij is the indexed value ratio for the j index for sample I and is considered using Equation (5): y ij is the consistent value of the jth parameter for the ith sample. The following equation is used to determine the normalizing construction function for the efficiency type: where C ij signifies the perceived value of the ith sample of jth parameter, and C j min and C j max signify the lowest and high standards of the jth parameter individually. Calculating a significant rating scale q j for each parameter is the second stage in EWQI estimation. The following equation is used to compute the q i value: while C j denotes each measured parameter concentration in mg/L, and S j represents the desired limit (mg/L) for every factor conferring to WHO standards.

Statistical Data Analysis
For understanding the data set, statistical analysis plays a vital role in representing various operations. Pearson correlation analysis was applied in SPSS (Armonk, NY, USA) to explain the relationship between water quality parameters. To interpret hydrochemical facies, a Piper diagram was organized over Aqua-Chem (version 2010.1). Saturation indices were measured over the geochemical simulation program PHREEQC (version 3.1), which determines the affinity of groundwater to dissolve a specific mineral. The geological map of the study area, and the water quality map for drinking purposes, were made using ARC GIS (version 9.3).

Groundwater Composition
The statistical summary of the physiochemical parameters of groundwater is shown in Table 1. Note: All parameters are represented in mg/L except EC, and As their unit is µs/cm and µg/L. The electrical conductivity value ranged from 305-5570 µs/cm, with a mean value of 1200 µs/cm, in all groundwater sources of district Sanghar. The high (EC) water samples reflect leakage or suspension of the aquifer constituents or other bases such as saline water bodies [45]. The pH of samples was in the range of 6.8 to 8.7, having a mean value of 7.6, showing the alkaline nature of groundwater sources. As a key water quality parameter, pH determination is compulsory due to its special effects on water chemistry, alkalinity, speciation, and solubility [46]. Total dissolved solids (TDS) were recorded, varied from 195-3564 mg/L with an average value of 739 mg/L, and were found inside the acceptable limit recommended by WHO. High TDS in groundwater indicated ion dissolution, which might be credited to progressively depleting salts and minerals over time [47]. Groundwater total hardness (TH) ranged from 13 to 1350 mg/L, with a mean of 316 mg/L. The rock-water interaction was mainly responsible for the significant random changes of the parameters mentioned above under extreme anthropogenic activities, which cause the solubilization of the minerals and salts [48]. Among cations, Na + exhibits a high mean concentration 125 mg/L, followed by Ca 2+ 61 mg/L, Mg 2+ 41.6 mg/L, and K + 2.81 mg/L in all groundwater sources and was within the recommended acceptable range of WHO. Under the anthropogenic influence, comprehensive interface and interaction among groundwater and adjacent rocks may result in higher Na + and K + concentrations [49]. The anions showed high dominance in comparison with cations. HCO 3 − was detected with high concentration varied from 90 mg/L to 590 mg/L, with a mean value of 233 mg/L. The elevation in HCO 3 − concentration is caused by the dissolution of calcite, carbonate, marble, and dolomite-bearing minerals [7,50,51]. The value of Cl − and SO 4 2− ranges from 20-850 mg/L, 20-1024 mg/L, respectively, with a mean value of 158.34 mg/L and 137.81 mg/L. NO 3 − values range from 0.1-27 mg/L having a mean value 0.97 mg/L. The concentrations of cations and anions were found in decreasing order of Na + > Ca 2+ > Mg 2+ > K + , and HCO 3 − > Cl − > SO 4 2− > NO 3 − , respectively. Arsenic had low to high concentration levels varied from 5-25 µg/L, with mean values of 12.9 µg/L in the study area. For the value of As in all groundwater samples, 28 samples (45%) had beyond the permissible limit out of sixty-one groundwater samples. As a result of natural and anthropogenic sources, elevated As concentration in groundwater sources has become a main environmental problem [52]. Water-logging and excessive pesticide use are all anthropogenic causes of arsenic enrichment [29,53].

Hydrochemical Facies
The hydrochemical facies represent the whole situation of groundwater clarifications interacting within a lithological configuration [54]. The richness of each ion as a proportion of total cations and anions was used to define the geochemical evolution of groundwater samples [55]. The Piper diagram (1944) provides a detailed graphical design of sample hydrochemistry and hydrochemical regimes [56]. Groundwater samples were plotted in a Piper diagram as shown in (Figure 2).
Maximum samples fell in the NaCl and CaHCO 3 type, which indicates that hydrochemical activities, anthropogenic actions, interface with the unsaturated region, improved resident time, and the rock-water interaction influences ion exchange in the study area. Few samples lie in the mixed type with CaMgCl, indicating that most depend on silicate weathering, dissolution of the carbonates, and ion-exchange processes in the aquifer. Similarly, regarding cations, the samples were lie in with Mg 2+ type and non-dominant type, Mg 2+ concentration persisted above the contents of other cations, indicating the supremacy of ion exchange and weathering of silicate. In contrast, in anions, most of the samples fall in with HCO 3 − type and No-dominant type in comparison with anions, depicting the eminence of carbonate weathering [57]. Most of the samples lie in both the cations and anions B zones mixed type, portraying the importance of silicate weathering and ion exchange [58]. The prevalence of Mg 2+ and HCO 3 − in both zones suggests that they are derived mainly through carbonate and sulfate solubilization. The dissolution of calcite, dolomite, and gypsum as the significant bases of these ions is generally indicated by a close relationship with the equilibrium line [57]. Maximum samples fell in the NaCl and CaHCO3 type, which indicates that hydrochemical activities, anthropogenic actions, interface with the unsaturated region, improved resident time, and the rock-water interaction influences ion exchange in the study area. Few samples lie in the mixed type with CaMgCl, indicating that most depend on silicate weathering, dissolution of the carbonates, and ion-exchange processes in the aquifer. Similarly, regarding cations, the samples were lie in with Mg 2+ type and non-dominant type, Mg 2+ concentration persisted above the contents of other cations, indicating the supremacy of ion exchange and weathering of silicate. In contrast, in anions, most of the samples fall in with HCO3 − type and No-dominant type in comparison with anions, depicting the eminence of carbonate weathering [57]. Most of the samples lie in both the cations and anions B zones mixed type, portraying the importance of silicate weathering  Table 2 shows the findings of principal component analysis (PCA) for groundwater parameters. Thus, following varimax rotation, PCA findings were attained to clarify the obtained sources that influenced the groundwater [7,59]. Four aspects, including PC1, PC2, PC3, and PC4 for groundwater parameters, were obtained with eigenvalues of 8.88, 1.68, 1.12, and 0.98 total variance of 68.33%, 12.99%, 7.08%, and 4.5%, respectively. , respectively. PC-1 had a high contribution of loading factors for the above parameters in PCA results. PC-1 demonstrated the geogenic and anthropogenic sources in the study area by showing a high contribution of moderate and strong positive loadings for all groundwater samples. The PC1 indicates the ionic configuration of groundwater and resultant from minerals dissolution, ion exchange, and weathering of host granitic rocks. Due to the high relationship between Na + and Cl − , PC1 describes the great difference in the data set, may be hard and salinity. The occurrence of carbonated rocks in the study area and the higher concentration of Na + and Cl − contributed to the high TDS values [60].

Principal Component Analysis (PCA) and Correlation Analysis of Groundwater Variables
Moreover, EC, TDS, and level specified were expected to influence by erosion of schistose rocks having sulfide minerals. The TDS, EC, Na + , and K + could be derived through bedrock leaching, implying that PC-1 has geogenic sources [7,28]. In contrast, the Cl − , NO 3 − , and SO 4 2− impurity is caused by anthropogenic sources such as animal manure, agriculture fertilizer, and atmospheric and soil sources [12,61,62]. The origins of SO 4 2− and HCO 3 − resulted from the dissolution of gypsum and calcium-bearing minerals and rock-water interaction [7,51]. Thus, PC1 accounted for mixed sources of geogenic and anthropogenic sources in the study area.
PC2 shows variables of groundwater water accounted for 12.92% variability with an eigenvalue of 1.68 ( Table 2). The high loadings values of groundwater were pH (r = 0.80), NO 3 − (r = −0.55), and As (r = 0.78). pH and NO 3 − levels are probably influenced by agriculture activities, while As can be credited to smelting activities of antimony ores in the study area [63]. Hence, PC2 exhibits anthropogenic sources in the study area. The groundwater variables PC3 and PC4 accounted for 7.08% and 4.5% variability, with eigenvalues of 1.12 and 0.98. PCA results support the hydrogeochemical processes. The significant correlation of PC3 and PC4 was K + and NO 3 − with coefficient (r) values of (0.91, 0.60). In the study area, both natural and anthropogenic activities play a significant role in groundwater pollution.
Furthermore, the groundwater pollution in the study area is caused by agricultural operations and household garbage [60]. The nitrate contamination resulted from agricultural practices, fertilizers, sewage, and animal manure [64]. Thus, the factors PC3, PC4 showed the anthropogenic causes in the research region.
The Pearson correlation constant (r) mostly finds how properly the groundwater samples are assembled in a straight line [65]. Both variables are highly dependent on one another, as seen by the positive correlation values. On the other hand, significant inverse relationships show that agriculture soil factors do not affect each other [62,66,67]. As a result, the inverse relationship indicates that the data set has little variability. Significant positive correlation pairs have been observed between different variables in groundwater sources.

Mechanisms Controlling Groundwater Chemistry
The formation of groundwater chemistry, evaporation, rock-weathering, crystallization, and precipitation are significant contributors [68]. Gibbs (1970) [69] proposed a quite beneficial model for understanding the mechanisms that control groundwater chemistry. Hence, to monitor the influences of hydrogeological interfaces of the groundwater and groundwater data, Gibbs plots were considered. The subplots indicate the relationship of TDS with the weight ratio of Cl − versus (Cl − + HCO 3 − ) and a weight ratio of (Na + + K + ) versus (Na + + K + + Ca 2+ ), respectively, as shown in (Figure 3).
A high percentage of samples are found in the rock dominance region, demonstrating that rock weathering is the primary factor influencing groundwater chemistry and evolution. The process of dissolvable salts and minerals becoming integrated with groundwater is assisted by parent rock weathering.
Furthermore, the mineral dissolution is assisted by the long resident time of rockwater interaction [70], further enhanced under the anthropogenic influence. Meanwhile, few samples lie in the evaporation dominance zone in the study area, identifying the concentrations of Na + and Cl − and the TDS levels in groundwater [7,71]. No sample occurs in the precipitation zone, which was sodomized to be negligible in the study area. The hydrogeological setting and weathering mechanisms play an essential role in the configuration of groundwater chemistry, as shown in the diagram. In addition to rock dominance, cation exchange is a credible mechanism for governing water chemistry composition in the study area [7]. The controlling mechanisms of rocks dominance contribution were consistent with the previous study [51,72]. TDS with the weight ratio of Cl − versus (Cl − + HCO3 − ) and a weight ratio of (Na + + K + ) versus (Na + + K + + Ca 2+ ), respectively, as shown in (Figure 3). A high percentage of samples are found in the rock dominance region, demonstrating that rock weathering is the primary factor influencing groundwater chemistry and evolution. The process of dissolvable salts and minerals becoming integrated with groundwater is assisted by parent rock weathering.
Furthermore, the mineral dissolution is assisted by the long resident time of rockwater interaction [70], further enhanced under the anthropogenic influence. Meanwhile, few samples lie in the evaporation dominance zone in the study area, identifying the concentrations of Na + and Cl − and the TDS levels in groundwater [7,71]. No sample occurs in the precipitation zone, which was sodomized to be negligible in the study area. The hydrogeological setting and weathering mechanisms play an essential role in the configuration of groundwater chemistry, as shown in the diagram. In addition to rock dominance, cation exchange is a credible mechanism for governing water chemistry composition in the study area [7]. The controlling mechanisms of rocks dominance contribution were consistent with the previous study [51,72].

Silicate Weathering
The significance of silicate weathering in influencing the main ion chemistry of groundwater is essential [73]. The highest rich cation in the study area is sodium, produced via silicate weathering and the breakdown of halite. The Na + /Ca 2+ values (>1) indicate the silicate weathering, and the distribution of samples away from the 1;1 in the Na + versus Ca 2+ plot suggested weathering of silicate minerals. The presence of additional Na + from silicate weathering was indicated by samples away [7]. The deviation from the 1:1 line (Figure 4a) is accredited to the silicate weathering and ion exchange, such as the dissolution of albite might be responsible for increasing Na + in groundwater if halite dissolution is the sole source of sodium. Hydro-geochemistry was governed by the association

Silicate Weathering
The significance of silicate weathering in influencing the main ion chemistry of groundwater is essential [73]. The highest rich cation in the study area is sodium, produced via silicate weathering and the breakdown of halite. The Na + /Ca 2+ values (>1) indicate the silicate weathering, and the distribution of samples away from the 1;1 in the Na + versus Ca 2+ plot suggested weathering of silicate minerals. The presence of additional Na + from silicate weathering was indicated by samples away [7]. The deviation from the 1:1 line (Figure 4a) is accredited to the silicate weathering and ion exchange, such as the dissolution of albite might be responsible for increasing Na + in groundwater if halite dissolution is the sole source of sodium. Hydro-geochemistry was governed by the association between aquifer lithology and penetrating water in the sub-surface [74]. The samples distribution along the equilibrium line of Ca 2+ and Mg 2+ against HCO 3 − + SO 4 2− (Figure 4b) plot revealed the influence of carbonate dissolution for Ca 2+ and Mg 2+ distributions of the samples along the equilibrium line of Ca 2+ + Mg 2+ versus HCO 3 − + SO 4 2− plot recommended that the effect of carbonate dissolution for Ca 2 + Mg 2+ . In contrast, the HCO 3 − + SO 4 2− values above five meq/L expected silicate weathering in groundwater sources (Figure 4c). Distributions of the samples along the 1:1 equiline of Ca 2+ + Mg 2+ versus HCO 3 − plot proposed the effect of calcite/dolomite dissolution (Figure 4d). tribution along the equilibrium line of Ca 2+ and Mg 2+ against HCO3 − + SO4 2− (Figure 4b) plot revealed the influence of carbonate dissolution for Ca 2+ and Mg 2+ distributions of the samples along the equilibrium line of Ca 2+ + Mg 2+ versus HCO3 − + SO4 2− plot recommended that the effect of carbonate dissolution for Ca 2 + Mg 2+ . In contrast, the HCO3 − + SO4 2− values above five meq/L expected silicate weathering in groundwater sources (Figure 4c). Distributions of the samples along the 1:1 equiline of Ca 2+ + Mg 2+ versus HCO3 − plot proposed the effect of calcite/dolomite dissolution (Figure 4d). The hydrochemistry of gypsum, anhydrite, and halite was controlled by the dissolution of sulfate and chloride minerals suggesting the under-saturated conditions. However, the HCO3 − formation was caused by the collaboration of penetrating groundwater with organic carbon in the topsoil and modest limestone dissolution in the aquifer [75]. The dominance of Ca 2+ and Mg 2+ in groundwater was caused by Na + substituting Ca 2+ in the aquifer matrix at promising exchange sites (i.e., Ca 2+ /Mg 2+ -bearing clay) reported by [76] in the reverse ion exchange field.
Silicate weathering is a comparatively slow process and contributes a minimal role in natural settings, and excessive silicate weathering indicates a significant anthropogenic The hydrochemistry of gypsum, anhydrite, and halite was controlled by the dissolution of sulfate and chloride minerals suggesting the under-saturated conditions. However, the HCO 3 − formation was caused by the collaboration of penetrating groundwater with organic carbon in the topsoil and modest limestone dissolution in the aquifer [75]. The dominance of Ca 2+ and Mg 2+ in groundwater was caused by Na + substituting Ca 2+ in the aquifer matrix at promising exchange sites (i.e., Ca 2+ /Mg 2+ -bearing clay) reported by [76] in the reverse ion exchange field.
Silicate weathering is a comparatively slow process and contributes a minimal role in natural settings, and excessive silicate weathering indicates a significant anthropogenic effect in the study area. A rock-water interface controlled general hydrochemistry via gypsum, anhydrite, halite dissolutions, and silicate weathering; hence, ion exchange mechanisms played a minor role.
The Na + normalized Ca 2+ against a HCO 3 − plot shows (Figure 5a) and an Nanormalized Ca 2+ against Mg 2+ plot (Figure 5b) was used to investigate the impact of silicate weathering, carbonate solubilization, and evaporite dissolution. effect in the study area. A rock-water interface controlled general hydrochemistry via gypsum, anhydrite, halite dissolutions, and silicate weathering; hence, ion exchange mechanisms played a minor role.
The Na + normalized Ca 2+ against a HCO3 − plot shows (Figure 5a) and an Na-normalized Ca 2+ against Mg 2+ plot (Figure 5b) was used to investigate the impact of silicate weathering, carbonate solubilization, and evaporite dissolution. The Na + normalized Ca 2+ versus Mg 2+ plot shows that most of the Mg 2+ are fall closed and derived from silicate weathering rather than carbonate and evaporite dissolution. The Na + -normalized Ca 2+ versus HCO3 − plot demonstrates that samples tend to fall into the carbonate and evaporite dissolution, except for silicate weathering. Overall, the results showed that most groundwater samples fall within the carbonate and evaporite dissolution, followed by the silicate carbonate weathering zone. The Na + normalized Ca 2+ versus Mg 2+ plot shows that most of the Mg 2+ are fall closed and derived from silicate weathering rather than carbonate and evaporite dissolution. The Na + -normalized Ca 2+ versus HCO 3 − plot demonstrates that samples tend to fall into the carbonate and evaporite dissolution, except for silicate weathering. Overall, the results showed that most groundwater samples fall within the carbonate and evaporite dissolution, followed by the silicate carbonate weathering zone.

Saturation Indices for Minerals Phases
The speciation of geochemical modeling of saturation index (SI) values for groundwater samples, illustrated in Figure 6.

Saturation Indices for Minerals Phases
The speciation of geochemical modeling of saturation index (SI) values for groundwater samples, illustrated in Figure 6. According to the geochemical modeling results, the aquifer conditions for SO4 2− and CO3 containing minerals, including gypsum having a value (−2.80) and halite (−8.82), were generally unsaturated. The negative SI values specify that there could be no prospective role in the sorption/desorption of As by these two mineral phases in groundwater [77]. On the other hand, the positive SI values for mineral phases of dolomite (2.6) and calcite (4.4) were observed. It demonstrates that these minerals could contribute to groundwater aquifers as the leading basis of As in the study area [78], releasing As due to the rockwater edge reactions. In the current study, the possible formation of SO4 2− and CO3 containing minerals and their role of As discharge under auspicious (alkaline) conditions could be found in dealing with groundwater's pH data, which were observed to be in alkaline limit in the majority of water samples. Similarly, the aquifers of Tharparkar, Pakistan [79] observed a high saturation index value of calcite (SI < 0), implying its precipitation and dissolution in the As release process.

Ion-Exchange
The chloro-alkaline indices (CAI-1 and CAI-2) (Schoeller 1965) were used to explore the ion exchange reaction. These indices have a substantial impact on water chemistry and evolution [80]. The following formulas were used to determine the CAI1 and CAI2: According to the geochemical modeling results, the aquifer conditions for SO 4 2− and CO 3 containing minerals, including gypsum having a value (−2.80) and halite (−8.82), were generally unsaturated. The negative SI values specify that there could be no prospective role in the sorption/desorption of As by these two mineral phases in groundwater [77]. On the other hand, the positive SI values for mineral phases of dolomite (2.6) and calcite (4.4) were observed. It demonstrates that these minerals could contribute to groundwater aquifers as the leading basis of As in the study area [78], releasing As due to the rock-water edge reactions. In the current study, the possible formation of SO 4 2− and CO 3 containing minerals and their role of As discharge under auspicious (alkaline) conditions could be found in dealing with groundwater's pH data, which were observed to be in alkaline limit in the majority of water samples. Similarly, the aquifers of Tharparkar, Pakistan [79] observed a high saturation index value of calcite (SI < 0), implying its precipitation and dissolution in the As release process.

Ion-Exchange
The chloro-alkaline indices (CAI-1 and CAI-2) (Schoeller 1965) were used to explore the ion exchange reaction. These indices have a substantial impact on water chemistry and evolution [80]. The following formulas were used to determine the CAI1 and CAI2: Water 2021, 13, 3477 14 of 20 All units are milliequivalent per liter (meq/L). The direct-ion exchange happens when both indices have positive values (Equation (10)). In contrast, the reverse ion exchange occurs when CAI-1 and CAI-2 have negative values (Equation (11)): Ca 2+ + 2NaX → 2Na + + CaX 2 (11) Cation exchange mechanisms are frequently identified using the interaction between (Na + + K + − Cl − ) and (Ca 2+ + Mg 2+ ) − (HCO 3 − − SO 4 2− ) [81]. Moreover, the Chlor-alkali index (CAI) can be used to identify whether it is reverse or direct cation exchange [82]. If the CAI value is more significant than <0, it donates cation exchange in which Ca 2+ in groundwater was exchanged for Na + in the aquafer, whereas CAI value > 0 shows reverse cation exchange. Positive CAI-1 and CAI-2 readings show that Na + and K + ions in water are exchanged with Mg 2+ and Ca 2+ ions. Moreover, Mg 2+ and Ca 2+ ions are exchanged with Na + and K + from rocks when their indices are negative. That is a sign of a chloro-alkaline imbalance, as shown in Figure 7a.  (10)). In contrast, the reverse ion exchange occurs when CAI-1 and CAI-2 have negative values (Equation (11)): 2 2 2Na CaX Ca 2NaX Ca 2NaX 2Na CaX Cation exchange mechanisms are frequently identified using the interaction between (Na + + K + − Cl − ) and (Ca 2+ + Mg 2+ ) − (HCO3 − − SO4 2− ) [81]. Moreover, the Chlor-alkali index (CAI) can be used to identify whether it is reverse or direct cation exchange [82]. If the CAI value is more significant than <0, it donates cation exchange in which Ca 2+ in groundwater was exchanged for Na + in the aquafer, whereas CAI value > 0 shows reverse cation exchange. Positive CAI-1 and CAI-2 readings show that Na + and K + ions in water are exchanged with Mg 2+ and Ca 2+ ions. Moreover, Mg 2+ and Ca 2+ ions are exchanged with Na + and K + from rocks when their indices are negative. That is a sign of a chloro-alkaline imbalance, as shown in Figure 7a. The distribution of multiple water sample sites from the study area was near the y = x line, showing cation exchange impacts of groundwater chemical composition. Meanwhile, as shown in Figure 7b, the bulk of water sample locations in the study area lie inside the CAI < 0 zones, showing the dominance of reverse ionization.
The slope value of −0.812 (Equation (12)), which is very close to the theoretical value of −1, indicates linearity in the relationship between parameters, showing ion exchange between Na + , Ca 2+ , and Mg 2+ (Figure 7b): The majority of the samples are plotted in the lower-left corner of Figure 7a, showing reverse ion exchange (Equation (11)), an increase of Na + , and a decrease in Ca 2+ in groundwater. On the other hand, few samples lie in the upper right corner, indicating direct-ion exchange (Equation (10)).

Groundwater Quality Assessment
The entropy water quality index (EWQI) technique was used to evaluate the excellence of groundwater in the research area. EWQI has been widely used to assess groundwater quality [83].
However, EWQI values are classified into five categories: excellent (<25), good (25-50), medium (50-100), poor (100-150), and extremely poor (>150). According to EWQI, the suitability of groundwater samples is presented in Table 4. The majority of water samples (n = 32) are considered "good", with a contribution rate of 54.24%. The samples (n = 20) fell in the "medium" category with a percentage contribution of 33.90%, whereas the samples (n = 17, and n = 11) were classified with 27.86% and 18.03% contribution for "poor" and "extremely poor" categories as shown in Table 4. None of the samples were considered in excellent drinking water class. Overall, most samples showed medium to good water quality levels, indicating that groundwater sources of the studied regions are suitable for drinking purposes except (n = 28) samples in the poor and extremely poor category in district Sanghar. The suitability map of water quality is shown in Figure 8.

Conclusions
This study has been conducted for the first time in this area based on several combinations of statistical analysis, graphical techniques, and hydrochemical modeling applied to understand groundwater chemistry, evolution, and suitability for drinking use in the studied region. Statistical analysis shows that the abundance of cations and anions is in a decreasing route, of Na + > Ca 2+ > Mg 2+ > K + , and HCO3 − > Cl − > SO4 2− > NO3 − . The richness of Na + and HCO3 − demonstrates that silicate weathering is the dominant factor affecting the main ion chemistry of groundwater in the study area. However, correlation and graphical relation between ions also reveal that ion exchange, rock weathering, and dissolution of carbonate minerals play an essential role in governing groundwater chemistry. Moreover, the Gibbs diagram showed that rock dominance had been identified as a fundamental natural component that controls groundwater evolution. In contrast, a few samples fall into the evaporation dominance zone, suggesting the importance of evaporation in the shallow groundwater depth zone. According to the hydrochemical facies data, the groundwater samples were mixed with an NaCl and CaHCO3 type. A few samples had mixed with CaMgCl type, interpreting the hydrochemical behavior of rock-water interaction.
Moreover, geochemical modeling-derived saturation indices show that the waterlogged phase is under-saturated compared to evaporite, such as gypsum and halite having negative values. In contrast, the excessive stage with carbonates such as calcite and dolomite has positive values. In the current study, 45% of groundwater samples showed As contamination above WHO recommended guideline values. The EWQI results revealed that 17 samples were found in poor and 11 samples in extremely poor quality, with

Conclusions
This study has been conducted for the first time in this area based on several combinations of statistical analysis, graphical techniques, and hydrochemical modeling applied to understand groundwater chemistry, evolution, and suitability for drinking use in the studied region. Statistical analysis shows that the abundance of cations and anions is in a decreasing route, of Na + > Ca 2+ > Mg 2+ > K + , and HCO 3 − > Cl − > SO 4 2− > NO 3 − . The richness of Na + and HCO 3 − demonstrates that silicate weathering is the dominant factor affecting the main ion chemistry of groundwater in the study area. However, correlation and graphical relation between ions also reveal that ion exchange, rock weathering, and dissolution of carbonate minerals play an essential role in governing groundwater chemistry. Moreover, the Gibbs diagram showed that rock dominance had been identified as a fundamental natural component that controls groundwater evolution. In contrast, a few samples fall into the evaporation dominance zone, suggesting the importance of evaporation in the shallow groundwater depth zone. According to the hydrochemical facies data, the groundwater samples were mixed with an NaCl and CaHCO 3 type. A few samples had mixed with CaMgCl type, interpreting the hydrochemical behavior of rock-water interaction.
Moreover, geochemical modeling-derived saturation indices show that the waterlogged phase is under-saturated compared to evaporite, such as gypsum and halite having negative values. In contrast, the excessive stage with carbonates such as calcite and dolomite has positive values. In the current study, 45% of groundwater samples showed As contamination above WHO recommended guideline values. The EWQI results revealed that 17 samples were found in poor and 11 samples in extremely poor quality, with a percentage value of 27.86% and 18.03%. However, no groundwater sample was found in excellent water status for drinking and household motive. The range mentioned above of groundwater samples is influenced by the dissolving process and leaching of rock-salt and gypsum-bearing rock formations.