Next Article in Journal
Biotransformation of 17β-Estradiol through a Denitrifying Sludge
Previous Article in Journal
Students’ Performance in Face-to-Face, Online, and Hybrid Methods of Teaching and Assessment in Anatomy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evolution Mechanism of Arsenic Enrichment in Groundwater and Associated Health Risks in Southern Punjab, Pakistan

1
Key Laboratory of Groundwater Resources and Environment, Ministry of Education, Jilin University, Changchun 130021, China
2
College of New Energy and Environment, Jilin University, Changchun 130021, China
3
Soil, Water and Environment Discipline, Khulna University, Khulna 9208, Bangladesh
4
School of Environmental Studies, China University of Geosciences, Wuhan 430074, China
5
Institute of Environmental Sciences and Engineering, School of Civil and Environmental Engineering, National University of Science and Technology, Islamabad 44000, Pakistan
6
School of Chemistry and Chemical Engineering, Shaanxi Normal University, Chang’an West Street 620, Xi’an 710119, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2022, 19(20), 13325; https://doi.org/10.3390/ijerph192013325
Submission received: 21 August 2022 / Revised: 29 September 2022 / Accepted: 12 October 2022 / Published: 15 October 2022
(This article belongs to the Section Water Science and Technology)

Abstract

:
Arsenic (As) contamination in groundwater is a worldwide concern for drinking water safety. Environmental changes and anthropogenic activities are making groundwater vulnerable in Pakistan, especially in Southern Punjab. This study explores the distribution, hydrogeochemical behavior, and pathways of As enrichment in groundwater and discusses the corresponding evolution mechanism, mobilization capability, and health risks. In total, 510 groundwater samples were collected from three tehsils in the Punjab province of Pakistan to analyze As and other physiochemical parameters. Arsenic concentration averaged 14.0 μg/L in Vehari, 11.0 μg/L in Burewala, and 13.0 μg/L in Mailsi. Piper-plots indicated the dominance of Na+, SO42−, Ca2+, and Mg2+ ions in the groundwater and the geochemical modeling showed negative saturation indices with calcium carbonate and salt minerals, including aragonite (CaCO3), calcite (CaCO3), dolomite (CaMg(CO3)2), and halite (NaCl). The dissolution process hinted at their potential roles in As mobilization in groundwater. These results were further validated with an inverse model of the dissolution of calcium-bearing mineral, and the exchange of cations between Ca2+ and Na+ in the studied area. Risk assessment suggested potential carcinogenic risks (CR > 10−4) for both children and adults, whereas children had a significant non-carcinogenic risk hazard quotient (HQ > 1). Accordingly, children had higher overall health risks than adults. Groundwater in Vehari and Mailsi was at higher risk than in Burewala. Our findings provide important and baseline information for groundwater As assessment at a provincial level, which is essential for initiating As health risk reduction. The current study also recommends efficient management strategies for As-contaminated groundwater.

Graphical Abstract

1. Introduction

Water is a critical resource for human life, but chemical and microbial contamination frequently challenges its safety [1,2,3,4]. Specifically, groundwater is one of the most commonly used water resources for drinking [5,6,7]. Therefore, there is always a need to provide and supply clean groundwater for drinking [8,9]. Arsenic (As) is a naturally occurring element with high toxicity and carcinogenic potential [10,11]. Groundwater As pollution has been a severe worldwide health and environmental problem for many years [12,13]. Drinking As-contaminated water has been linked to various health problems in Asian countries [14]. Chronic exposure to As through drinking groundwater may lead to health risks, such as kidney cancer, skin lesions, neurological diseases, and liver damage, and death may occur in the most severe cases [15]. As release in aquifers is caused by geochemical changes in underground sediments and is usually ascribed to the process of the reductive dissolution of arsenic-bearing minerals in aquifers [16,17,18]. A few examples include the water–rock interactions, sorbent/desorbent processes, residence time along the flow channel, oxidative/reductive dissolution processes, and the variation in the chemical composition of groundwater containing arsenopyrite (FeAsS), and calcium carbonate (CaCO3), all of which can release As in groundwater systems [19,20,21]. As also enters into aquifers from human activities like mining and smelting, industrial waste disposal, fertilizer use, agricultural operations, and wastewater irrigation [18]. As in the aquatic environment is mainly found in the forms of arsenite (As(III)) and arsenate (As(V)). As(III) is more widespread in reduced environments, whereas As(V) is abundant in oxidized environments [22,23].
As-contaminated groundwater threatens millions of people’s health and lives, especially in densely populated places. Over 200 million people in 105 countries are affected, and, in Asia, the health of about 150 million people is at risk through the consumption of As-contaminated groundwater [24,25,26]. Although some countries have lower arsenic concentrations in the groundwater (<10 μg/L), many others, including Bangladesh [27], China [28], Egypt [29], India [30], Indonesia [31], Nepal [32], and the Philippines [33] have over than 50 μg/L of As in groundwater aquifers. Pakistan is also one of the As-contaminated countries, facing water scarcity and groundwater vulnerabilities [34]. A recent study on groundwater pollution in Pakistan emphasized the problem of As contamination, especially in the Punjab and Sindh provinces connected to the Indus River and its tributaries [35,36,37]. In these contaminated regions, arsenic concentration was estimated 2580 μg/L in Lahore and Kasur (Punjab) [38,39], 905 μg/L in Muzaffargarh (Punjab) [40], 201 μg/L in Chichawatni (Punjab) and 29.30 μg/L in Dera Ghazi Khan (Punjab) [35,41]. In addition, arsenic concentration was as high as 1516 μg/L in Jamshoro (Sindh) [42], 158 μg/L in Manchar lake (Sindh) and 2580 μg/L in Tharparkar (Sindh) [43,44]. In Hyderabad, over 40 people died in 2004 because of polluted drinking water with high levels of As [45]. About 70% of Pakistan’s groundwater and surface water is contaminated with biological, organic and inorganic contaminants [46]. According to the current state of As contamination in Pakistan’s drinking water, about 47 million people live in places where more than 50% of well water contains >10 μg/L of As [46] and 17% surpasses 50 μg/L [47]. Only 26% of Pakistan’s population has access to safe drinking water [48]. The health risk map for As concentration in 1187 instances of groundwater suggests that the Indus Basin is severely contaminated by As [49]. The elevated As concentration may be attributable to oxidative desorption driven by many physicochemical parameters [47]. Numerous reports have explored As contamination in groundwater in various places, including rural and urban areas [24], peri-urban areas [50], rivers [47] and health facilities [51], while some earlier studies established As distribution on a small scale without comparisons between tehsils in Punjab’s groundwater and health risk assessments.
Considering the above-mentioned circumstances, this study explored As distribution, formation and evolution mechanisms in groundwater through geochemical modeling and multivariate analysis. We aimed to unravel the As contamination situation in Punjab’s groundwater and conduct an assessment on the health risks caused by the consumption of As-contaminated groundwater.

2. Materials and Methods

2.1. Study Area

Vehari, Burewala and Mailsi are three sub-districts (tehsils) of the Punjab Province in Pakistan. The area is bounded by the Sutlej and Ravi rivers and is located at 83–170 m asl (altitude) between 30°04′19″ N and 72°35′28″ E, as illustrated in Figure 1. The total population is 3 million and the climate is hot in the summer, from April to August, rising to a maximum of 50 °C, and cold in winter, from November to February, dropping to a minimum of about 5 °C. Dust storms are typical throughout the summer, and the annual average rainfall is around 125 mm, most of which falls during the monsoon season (July and August) [52]. Electric pumps and tube wells are used to exploit groundwater.
The Indus and Jhelum rivers pass around the study area. In contrast, the Ravi and Satluj rivers run through the study area, the rainwater and rivers serves as recharge sources for the aquifer. The south Indus River forms the plain alluvial deposition and its five major tributaries, which contain sediments from the Pleistocene and Holocene, are carried by the Ravi and Sutlej rivers. Unconsolidated alluvial deposits form the aquifer with varying amounts of sand with a large proportion of fine sand and silt and low organic matter. Since the late tertiary period, these deposits have been formed through the Indus River and its tributaries in the wide alluvial plain, which stretches from the Himalayan foothills to the Arabian Sea. The minerals, including anhydrite, aragonite, calcite, dolomite, goethite, gypsum, halite, and hematite, were detected in the mineralogical assessment [24,53]. The tributaries of the Indus River deposited a 400-m thick layer of sediment throughout the Pleistocene. Groundwater originated in interbedded silt and alluvial sand layers that are well-distributed across the Punjab area. This region’s residents depend on groundwater for drinking, livestock, and agriculture.

2.2. Groundwater Samples

In total, 510 groundwater samples were collected from shallow aquifers (<35 m) of three tehsils, namely Vehari, Burewala, and Mailsi from April–May 2020 and filtered to 0.45 μm for further analysis. A global positioning system (GPS) was used to record the groundwater sample’s location. (APHA et al. 2005) The American Public Health Association’s standard methods were followed [54]; 250 and 100 mL glass bottles were used to collect samples for the analysis of cations and anions, respectively. A drop of nitric acid (HNO3) was immediately added to the samples to adjust pH (<2.0) for cation analysis.

2.3. Hydrogeochemical Processing and Health Risk Assessment

All analyses followed the methods prescribed by the APHA 2005 [54]. The samples were analyzed at the Pakistan Council of Research in Water Resources (PCRWR) Islamabad, Pakistan. The calibrated portable meters (RS232C/Meter CON 110) were used to measure water pH, turbidity, electric conductivity and total dissolved solids on-site. Arsenic and major cations (Na+, Ca2+, Mg2+, K+, Fe2+) were analyzed by inductively coupled plasma atomic emission spectrometry (ICP-AES/27). An ion chromatograph (IC/P60) was used to measure SO42−, Cl, NO3−, and F, while HCO3 and total alkalinity were measured using the acid-base titration method. Finally, the samples were further examined for accuracy by computing percentage charge balance errors (CBE), as seen in Equation (1).
% CBE = cations anions cations + anions   × 100
Here, the unit of cations and anions is meq/L and the physicochemical analysis with a %CBE within ±5% are considered perfect for further research.
SPSS v16.0 (SPSS Inc., Chicago, IL, USA) was used for the statistical analysis of groundwater physicochemical parameters. XLSTAT v19 was used for Pearson’s correlation and principal component analysis (PCA) to determine the factors influencing groundwater quality. Piper, 1994 [55], was incorporated, using Golden Software Grapher 18.3 (Golden, CO, USA). Origin v19 (Sydney, Australia) plotted the ionic ratios to describe the dominant ionic water facies displaying the water chemistry in the study area.
According to USEPA (2005), the chronic daily intake (CDI), hazard quotients (HQ), and cancer risk (CR) were evaluated [56,57]. Two models were used for the health risk assessment of exposure to uncover As-contaminated groundwater. Carcinogenic and non-carcinogenic health hazards were determined by CDI, as described in Equation (2). Non-carcinogenic risk was estimated by HQ, as listed in Equation (3), and carcinogenic risk was calculated by CR, as described in Equation (4).
CDI = C   × IR BW
HQ = CDI RfD
CR = CDI × CSF
Here, C represents the As concentration in groundwater (mg/L), and IR represents daily water consumption, 0.5 L/day for children and 2 L/day for adults. BW represents body weight, 20 kg for children and 70 kg for adults. RfD refers to the oral reference dose of As, 0.3 mg/kg/day for children and adults, and CSF denotes the cancer slope factor for As, 1.5 mg/kg/day for children or adults [57]. Groundwater samples with (HQ < 1.0) were considered safe to drink, and groundwater samples with (CR > 10−4) exhibited considerable cancer risk [57].

2.4. Geochemical Modeling

PHREEQC version 3.0 software package was used to calculate the mineral saturation level [58]. The saturation index was used to illustrate the representative mineral phases, including anhydrite, aragonite, calcite, dolomite, goethite, gypsum, halite, and hematite, which are essential characteristics in the study area. The saturation indices (SI) that explain the thermodynamic tendency of minerals to precipitate or dissolve were calculated using Equation (5).
SI = Log IAP   Ksp   = LogIAP LogK Sp
Here, IAP represents the ion activity product for the solution of dissociated types, and Ksp represents the equilibrium solubility product at groundwater sample temperature for the chemicals involved. The hydrochemical equilibrium state was at SI = 0, whereas negative (SI < 0) and positive (SI > 0) values indicated the under-saturation and oversaturation of the minerals, respectively.
The inverse groundwater geochemical evolution modeling over two flow pathways along the (A to C) and (X to Z) were carried by PHREEQC (Figure 1). According to the groundwater flow pattern and hydrogeological conditions, the selected path I flows from A to B and B to C, and the selected path II flows from X to Y and Y to Z. Major ions (Na+, Mg2+, Ca2+, Fe2+, K+, HCO3, SO42–, Cl, and NO3) were used in the model simulation. The mineral phases include anhydrite, aragonite, calcite, dolomite, goethite, gypsum, halite, and hematite. Due to the wide-ranging distribution of the clay minerals in the study area, the exchange of cation is also considered. For significant tendencies, our model moreover considers the CaX2/NaX, CO2(g), and H2O(g) phases [59].

3. Results and Discussion

3.1. Characteristics of Groundwater Geochemistry

Groundwater hydrogeochemical properties in three tehsils are listed in Table 1 and compared to the WHO 2022 recommended drinking water quality standard [60]. EC ranged from 308 to 4550 μS/cm (averaging 1570 μS/cm), 85 to 4400 μS/cm (averaging 1231 μS/cm), and 226 to 2690 μS/cm (averaging 1185 μS/cm) in Vehari, Burewala, and Mailsi, respectively, exceeding recommended guidelines. Such huge changes in EC are usually caused by human activity and mineral dissolution in groundwater. Our findings imply that geochemical reactions, rock–water interactions, and anthropogenic sources influence the groundwater chemistry in the study area [61]. However, higher TDS in Vehari ranged from 234 to 3148 mg/L (averaging 1072 mg/L), Burewala 277 to 3173 mg/L (averaging 896.68 mg/L), and Mailsi 359 to 3298 mg/L, (averaging 896 mg/L). Water pH values in Vehari, Burewala, and Mailsi ranged from 6.85 to 7.61 (averaging 7.17), 6.78 to 7.11 (averaging 6.98), and 7.2 to 8.35 (averaging 7.66), respectively. The chemical composition of groundwater is pH-dependent, the groundwater chemistry fluctuates with changes in pH, and this variance is mostly driven by different chemicals, such as fertilizer for agriculture activities [5]. In the study area, the groundwater pH was alkaline, which might be caused by the presence of HCO3 and rock weathering. The turbidity in Vehari (26.37 NTU) and Burewala (15 NTU) were higher, whereas the turbidity in Mailsi (1.22 NTU) was lower than the WHO guidelines. Elevated EC and TDS levels in groundwater are caused by electrolytes and high salinity, which are often associated with semi-arid and dry climatic conditions and are responsible for the increased turbidity in Vehari and Burewala. Furthermore, the alkaline environment may increase conductivity over time by hastening the dissolving process. The obtained higher EC values here might be explained by the presence of more dissolved salts in groundwater. High TDS levels are most likely caused by wastewater discharged into pits from residential and dyeing units, ponds, and lagoons. The average value of total hardness in Vehari, Burewala, and Mailsi was 368.47, 348.85, and 414.58 mg/L, respectively. The higher concentration of bicarbonate in the groundwater may be attributed to rock weathering, atmospheric sources, and anthropogenic activities. The average HCO3 contents were: Vehari, 345.08 mg/L; Burewala, 299.20 mg/L; and Mailsi (324.23 mg/L). The majority of HCO3 comes from geological sources, such as carbonate dissolution and carbonate cement formations. Both anion and cations observed consistent sequences—cations: Na+ > Ca2+ > Mg2+ > K+ > Fe2+; and anions: HCO3 > SO42− > Cl > NO3 > F. As concentrations in Vehari, Burewala, and Mailsi ranged to 45 μg/L (averaging 14.0 μg/L), 52 μg/L (averaging 11.0 μg/L), and 89 μg/L (averaging 13.0 μg/L), respectively. They all exceeded the permitted limit of 10 μg/L. More precisely, As concentrations in Vehari and Mailsi was significantly higher than Burewala because of high bicarbonate (HCO3) content and pH [62]. In this study, the high concentrations of HCO3 and alkaline water quality indicated an oxidative condition in aquifers. Most areas in Southern Punjab belong to tributary areas, where groundwater has high minerals and As concentrations, which are not suitable for drinking. Thus, Vehari and Mailsi face the critical challenge of the most widespread As contamination. Although a prior investigation indicated a significant Fe2+ content in Jamshoro [42] and compared to other As-affected places the southeast Asia region [63], Fe2+ iron in our study was lower than that report, while the high HCO3 may dissolve the carbonate minerals, possibly explaining the higher As concentration in groundwater.

3.2. Hydrogeochemical Origin of Groundwater

3.2.1. Geochemical Classifications of Groundwater

Piper (1994) proposed a convenient method to classify and compare water chemistry types based on an ionic composition by plotting the hydrogeochemical data on a trilinear diagram [55]. The first type of the triangle in Piper’s plot represents anions and the second type represents cations (Figure 2). The groundwater type is dominated by Na+, Ca2+, and Mg2+ amongst cations, while HCO3, SO42−, and Cl amongst anions. Cation concentrations were in the order of Na+ > Ca2+ >Mg2+ > K+ > Fe2+, and the anions were ranked as HCO3 > SO42− > Cl > NO3− > F. Rasool et al. reported a similar pattern of anion and cation concentrations [64]. In addition, groundwater showed mixing behavior to Ca2+, SO42−, and Mg2+, which explained the geochemical formation with calcium, magnesium, and bicarbonate-containing mineral phases [65]. The Na+ and Cl type of groundwater indicated that rock–water interaction and agricultural activities played significant roles in the study area. In Piper’s plot, Na+−SO42− and Ca2+ and Mg2+ types of waters indicated As release by the sedimentary rocks into groundwater. Several important mechanisms result in As mobilization in groundwater, such as the dissolution of calcium, the dissolution of salt minerals, and desorption due to high pH [66,67].

3.2.2. Ionic Ratios and Mineral Phases of Groundwater

The main factors determining arsenic concentration in the groundwater of the study area were identified using ionic ratios. The plots [(Ca2 + + Mg2+) − (HCO3 + SO4)] and [(Na + K) − Cl] were employed to identify the exchange of cations in aquifers. [(Ca2+ +Mg2+) − (HCO3 + SO4)] illustrated that in the system, Ca2+ and Mg2+ amounts were removed/added separately from calcite, dolomite, and gypsum, whereas [(Na + K) − Cl] illustrated the sodium removed/added in the system separately to the chloride salts [68]. The ionic exchange of Na+, Ca2+, and Mg2+ within the system was indicated by the linear slope of −1. Therefore, our findings suggested the strong effect of cation exchange on significant correlations in Vehari (−0.95, slope = −0.93), Burewala (−0.96, slope = −0.99), and Maisli (−0.91, slope = −0.96) (Figure 3a). The importance of cation exchanges was further emphasized by the (Ca2+ + Mg2+ vs. SO4 + HCO3) plot (Figure 3b). Most groundwater samples were below the line of 1:1, showing the effects of silicate weathering and cation exchange in the study area [68]. The carbonate dissolution is the main supplier of sodium, magnesium, calcium, and bicarbonate ions into the groundwater. The plots revealed that silicate weathering for the contamination of the groundwater was influenced by dissolved carbonate minerals (HCO3/Na+) vs. (Ca2+/Na+) (Figure 3c) and (Mg2+/Na+) vs. (Ca2+/Na+) (Figure 3d).
The groundwater hydrogeochemistry was also supported by the saturation index (SI) and always affected by the minerals [69]. The SI estimation facilitates understanding the reaction pathways and the measurement of mineral dissolution and precipitation. In the geochemical simulation model (Table S1), aquifer conditions were undersaturated (SI < 0) with calcium carbonate and rock salt minerals, including aragonite (CaCO3), calcite (CaCO3), dolomite ((CaMg(CO3)2), and halite (NaCl). These mineral phases had negative SI values and were unlikely to precipitate, possibly playing a vital role in releasing As into aquifers due to their dissolution [70,71]. In contrast, the SI was positive for anhydrite (CaSO4), gypsum (CaSO42H2O), and iron oxide mineral phases, including goethite (FeOOH) and hematite (Fe2O3). These minerals tended to participate in groundwater (Figure 4a–c). Hasan et al. (2009) found that iron oxides in the sediments of the flood plain in Bangladesh inhibited As mobility in groundwater [72].

3.2.3. Correlation Analysis

In Pearson’s analysis, correlation was marked as weak (r < 0.3), moderate (r 0.5), and strong (r ≥ 0.7) [73]. From the correlation matrices for Vehari, Burewala, and Mailsi (Tables S2–S4), we were able to understand the geochemical process in the study area. In Vehari, As exhibited a negative correlation with EC, HCO3, SO42, Ca2+, and Fe2+ (p < 0.05), TDS, Mg2+, K+, and hardness (p < 0.01) (Table S2). Such correlations highlighted the influence of pH on As concentration in groundwater, matching well with a previous study [49]. In contrast, groundwater As in Burewala exhibited significant positive correlations with Ca2+ and hardness (p < 0.01) but negative correlations with NO3 and F (p < 0.05) (Table S3). In Mailsi, As was positively correlated with turbidity (p < 0.01) and Mg2+ (p < 0.05) but negatively correlated with Ca2+ (p < 0.05) and NO3 (p < 0.05) (Table S4). In all three tehsils, EC had a significant positive correlation with TDS, cations, and anions at the 1% level but not with Fe2+. The correlation of EC with other parameters indicated the higher possibilities of ion exchange in the aquifers [74].

3.2.4. Principle Component Analyses (PCA)

To identify the probable As sources in aquifers, PCA was used as a multivariate method lowering the number of independent variables to limited fundamental components [73]. As listed in Table 2, the three main components (PC1 to PC3) explained 59.15% of the variance in groundwater in Vehari. PC1 explained 38.04% of the variance and mainly consisted of EC (0.964), TDS (0.981), SO42− (0.936), Cl (0.867), Na+ (0.856), Mg2+ (0.816), HCO3 (0.707), and hardness (0.662), indicated the process of As geogenic in groundwater. PC2 contributed to 12.83% of total variance, mainly consisted of Ca2+ (0.831) and hardness (0.71). Thus, calcium was closely linked with hardness and mostly originated in sedimentary rocks. PC3 explained 8.2% of the total variance, predominated by As (0.514), pH (0.376), and turbidity (0.329). It hinted at geogenic origin and enrichment of As in Vehari.
In Burewala, the first three components explained 63.15% of the total variance. PC1 contributed 44.49% of the total variance, mainly comprising EC (0.957), TDS (0.962), alkalinity (0.837), HCO3 (0.821), Cl (0.844), SO42– (0.856), Ca2+ (0.592), Mg2+ (0.833), and Na+ (0.859). This indicated that mineral weathering and the rock–water interaction determine cation and anion distribution. PC2 was 10.88% of the total variance, including Ca2+ (0.706), hardness (0.527), and As (0.527), which explained As release by the sedimentary rocks in aquifers. PC3 explained 7.77% of the total variance and consisted of NO3 (0.648). This might be explained by the inappropriate usage of nitrate fertilizers, such as urea for irrigation purposes and anthropogenic activities, hinting at their significant impacts on groundwater quality [48].
Among the three main components explaining 62.86% of total variance in Mailsi, PC1 contributed 44.51% of the total variance and consisted of EC (0.869) and TDS (0.972). This results met well with the positive links between EC, TDS and hardness, evidencing the roles of water–rock interactions in substantial HCO3 due to the dissolution of carbonate minerals [39]. PC2 explained 11.33% of the total variance and included turbidity (0.745) and As (0.783), indicating that arsenic was mobilized geologically [42]. PC3 explained 7.01% of the total variance and consisted of pH (0.512) and Mg2+ (–0.502).
Taking all results in the three regions together, EC and TDS were positively linked with both cations and anions in the groundwater. Previous studies have used PCA and the hierarchical cluster analysis (HCA) to identify potential As sources in the groundwater. Bibi et al. (2015) suggested that As in groundwater was geogenic because of the weathering of sulphide minerals [75]. Brahman et al. (2013) also used PCA to investigate and report high concentrations of SO42– (398–409 mg/L), low Fe2+ concentration (0.097–0.132 mg/L), and high pH values (7.9–9.2) in the groundwater of Tharparkar (Sindh) [44]. These conditions indicate that As might be released in the oxidizing conditions of the quaternary sediment [76,77]. In this study, we observed similar levels of EC and TDS, anions, and cations in groundwater, hinting at the geogenic sources of As release in groundwater.

3.2.5. Groundwater and Arsenic Evolution along Flow Path

PHREEQC was used for geochemical inverse modeling [58], and two predefined flow lines were introduced in model to understand the effects of hydrogeochemical evolution on As mobilization along groundwater flow paths in this study. Six groundwater samples were used for geochemical inverse modeling (A-B-C and X-Y-Z) based on flow paths I and II, respectively. Briefly, four hydrogeochemical processes included: (1) the dissolution of calcium and salt bearing-minerals, (2) chemical weathering, (3) exchange of cations between Ca2+ and Na+, and (4) CO2 accounting for the interaction of atmosphere or microbial respiration resulting from organic decomposition. Modeling results (Figure 1 and Table S5) illustrated that anhydrite, gypsum, and iron oxide minerals precipitated, while other rock salt minerals, including goethite, hematite, halite, and NaX, remained in the recharge zone. Cation exchange was also involved in this process, Na+ (0.03124 mmol/L) was dissolved in groundwater being released from the soil and rock surface through the dominance of Ca2+, Na+, K+, and HCO3 flow of groundwater. When groundwater moved from B to C, where anhydrite, goethite, and hematite precipitated, the cation ion exchange Ca2+ and Na+ were released into groundwater. Generally, the increase of As from A to C was attributed to the dissolution process of calcium-bearing minerals.
In Path II X–Y, anhydrite precipitation and the cation exchange CaX2 (0.01642 mmol/L) was specified by the Na+ + K+-HCO3 facies. As the groundwater moved from Y to Z, the carbonate minerals were dissolved through the stronger exchange of cations, Ca2+ interchanged with Na+. The As concentration in groundwater was elevated along with the flow path X–Z, as evidenced by the increased HCO3, Cl, and TDS. This result suggested that the additional sources of As were introduced into groundwater [59]. The increasing HCO3, Cl, and TDS concentrations in groundwater might be related to the dissolution of minerals [70,71]. Notably, the dissolving of calcium carbonate and salt mineral (halite) is evident from the groundwater conditions. Our results are well supported by our saturation indices (SI < 0) results, which indicated the dissolution of calcium carbonate and rock salt minerals, possibly playing a vital role in releasing As into groundwater. This implies that the groundwater condition is affected by the higher dissolution of calcium carbonate and salt minerals regarding the As release into groundwater.

3.3. Mechanism of As Enrichment and Geospatial Distribution

Iron concentration is low in groundwater when groundwater oxygenation is available [78]. In this study, rivers pass through the study area, and the clay-covered aquifer top layer allows aerobic activities in the groundwater underneath. The SI of different minerals from the PHREEQC model suggested that the majority of groundwater was either supersaturated or near saturation in terms of anhydrite, gypsum, and iron oxide mineral phases, including goethite (FeOOH) and hematite (Fe2O3), which precipitated As in groundwater. In contrast, the negative SI of aragonite, calcite (CaCO3), dolomite ((CaMg(CO3)2)), and halite (NaCl) explained the evolution and mobilization of As in groundwater, and As release and mobility was associated with the dissolution of calcium carbonate, bicarbonate, and salt minerals [70,71]. Additionally, the fewer amounts of minerals such as FeCO3 react with natural calcium carbonate minerals, possibly causing their dissolution in the aquifer and releasing As into groundwater [70,79]. The presence of natural organic compounds and microbial activity may also be the reason to contaminate the groundwater [80,81,82,83]. Anthropogenic sources, the dissolution of CaCO3, and coal mining are probable sources of sulfur [70,84,85]. The positive correlations between As and pH in all three tehsils hinted at the critical roles of pH in As release in the study area [49], and the increase in the mineralization rates were also assumed to play a major role in the evolution mechanism of As [86].
The groundwater contamination might be interrelated with the geochemical processes, as shown in Figure 5a. The spatial distribution of As was mapped through the inverse distance weight (IDW) interpolation (Figure 5b). The fluctuating As concentrations in groundwater was attributed to different reasons, mainly due to the higher urbanization, industrialization, and coal mining in the Vehari and Mailsi tehsils compared to Burewala. It can also be clearly seen that the river crosses almost through middle of the Vehari and Mailsi; therefore, there are higher As concentrations in the Vehari and Mailsi study area than in Burwala. Baloch et al. (2021) reported an average As concentration of 1.53 μg/L in Sakrand Sindh, which is lower than in this study and others reported high concentrations in Punjab [73,87].

3.4. Health Risk Assessment

Health risks from groundwater As contamination were estimated by CDI, HQ, and CR for children and adults by considering drinking as the main exposure to As (Table S6). In this study, the CDI of groundwater in Vehari had an average value of 3.50 × 10−4 and 1.0 × 10−4 μg/(kg·day) for both children and adults, respectively. The average HQ was 1.17 and 0.30 for both children and adults, respectively, and the cancer risks were 5.27 × 10−4 and 1.51 × 10−4. For groundwater Burewala, the average CDI, HQ, and cancer risks were 2.71 × 10−4 (children) and 7.75 × 10−5 μg/(kg·day) (adults), 0.9 (children) and 0.25 (adults), and 4.07 × 10−4 (children) and 1.16 × 10−4 (adults), respectively. The groundwater in Mailsi exhibited slightly higher CDI (3.40 × 10−4 for children and 1.36 × 10−4 for adults), HQ (1.13 for children and 0.45 for adults), and cancer risk (5.10 × 10−4 for children and 2.04 × 10−4 for adults). Overall, the groundwater in Vehari and Mailsi exhibited more risk than that in Burewala. The carcinogenic risk was (CR > 10−4) in all three tehsils, indicating considerable cancer risks for adults and children (Figure 6a). There was significant non-carcinogenic risk for children (HQ > 1), which was higher than adults (Figure 6b).
The CDI in this study was similar to those reported in Punjab, averaging 3.4 × 10−4 μg/(kg·day) [88], slightly higher than previous investigations in Khyber Pakhtunkhwa [0 to 5.6 × 10−4 μg/(kg·day)] [89] and Lahore [0.11–3.7 × 10−4 μg/(kg·day)] [90]. However, they were much lower than those reported in Bangladesh [50–500 × 10−4 μg/(kg·day)], which were caused by the elevated arsenic levels in Bangladesh’s groundwater [91], and Turkey [0.023–5.2 × 10−4 μg/(kg·day)] [92]. In this study, carcinogenic risk indices showed non-neglectable health threats to both children and adults regarding As exposure, while non-carcinogenic risks were only obvious for children in Vehari and Mailsi. Shakoor et al. (2016) also reported a high HQ (0.12–18.5) in Punjab [88]. These findings suggested that people in study area who ingest As-contaminated groundwater face significant cancer risks, such as stomach cancer, spontaneous abortions, diabetes, and gastric and esophageal complications. The As level in groundwater needs constant monitoring and remediation to protect the health of local people.

4. Conclusions

In this study, we used statistical analyses and hydrogeochemical modeling to understand groundwater hydrogeochemistry, the evolution mechanism of As, and related health risks in Vehari, Burewala, and Mailsi (Punjab, Pakistan); levels were 14.0 μg/L, 11.0 μg/L, and 13.0 μg/L, respectively. Piper plots indicated that groundwater Na+-SO42−, Ca2+, and Mg2+ types were more dominant and ionic ratios showed that the As contamination in groundwater was influenced by the dissolution of calcium carbonate. Furthermore, geochemical modeling showed negative saturation indices with calcium carbonate and salt minerals, including aragonite (CaCO3), calcite (CaCO3), dolomite [CaMg(CO3)2], and halite (NaCl). Findings suggested their critical roles in As mobilization/release in aquifers, also evidenced by multivariate analysis. The current study proposed the clear concept of the As evolution mechanisms in the aquifer, as evidenced by the inverse model, namely calcium-bearing minerals, chemical weathering, the exchange of cations between Ca2+ and Na+, and anthropogenic activities. Geomorphology and mineralogical geochemistry are important in order to understand how the dissolution of mineral phases is determined the As-mobilization in the groundwater and how the agrochemical utilization in agricultural lands should be properly treated to protect groundwater resources in Pakistan. The spatial distribution of As followed the order of Vehari > Mailsi > Burewala. The possible cancer risk (CR > 10−4) was significant for both children and adults, while the non-carcinogenic risk hazard quotient was only significant for children. It is strongly recommended that decision-makers in Pakistan should adopt immediate management actions and develop long-term strategies for protecting groundwater resources from As pollution.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijerph192013325/s1.

Author Contributions

M.Y.J.B. and W.Z.: Writing—original draft, conceptualization, methodology, formal analysis, investigation, data curation. W.Z.: Supervision, resources, conceptualization, methodology, funding acquisition, writing—review and editing. D.Z.: Methodology, investigation, review and editing. B.A.A.S.: Formal analysis, review and editing. J.I., S.L., J.C., M.A.F. and A.P.: review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Key R&D Program of China (2019YFC1804804).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the first or corresponding authors.

Acknowledgments

We would like to thank the anonymous reviewers and editors for their valuable comments.

Conflicts of Interest

The authors declare that they have no competing financial interest or personal relationship that could have appeared to influence the work reported in this paper.

References

  1. Del Olmo, G.; Husband, S.; Sánchez Briones, C.; Soriano, A.; Calero Preciado, C.; Macian, J.; Douterelo, I. The microbial ecology of a Mediterranean chlorinated drinking water distribution systems in the city of Valencia (Spain). Sci. Total Environ. 2021, 754, 142016. [Google Scholar] [CrossRef] [PubMed]
  2. Zhang, W.; Zhu, Y.; Gu, R.; Liang, Z.; Xu, W.; Jat Baloch, M.Y. Health Risk Assessment during In Situ Remediation of Cr(VI)-Contaminated Groundwater by Permeable Reactive Barriers: A Field-Scale Study. Int. J. Environ. Res. Public Health 2022, 19, 13079. [Google Scholar] [CrossRef]
  3. Zhang, W.; Chai, J.; Li, S.; Wang, X.; Wu, S.; Liang, Z.; Baloch, M.Y.J.; Silva, L.F.; Zhang, D. Physiological characteristics, geochemical properties and hydrological variables influencing pathogen migration in subsurface system: What we know or not? Geosci. Front. 2022, 13, 101346. [Google Scholar] [CrossRef]
  4. Xia, X.; Teng, Y.; Zhai, Y.J.S. Biogeochemistry of Iron Enrichment in Groundwater: An Indicator of Environmental Pollution and Its Management. Sustainability 2022, 14, 7059. [Google Scholar] [CrossRef]
  5. Zhai, Y.; Han, Y.; Lu, H.; Du, Q.; Xia, X.; Teng, Y.; Zuo, R.; Wang, J. Interactions between anthropogenic pollutants (biodegradable organic nitrogen and ammonia) and the primary hydrogeochemical component Mn in groundwater: Evidence from three polluted sites. Sci. Total Environ. 2022, 808, 152162. [Google Scholar] [CrossRef]
  6. Zhai, Y.; Zheng, F.; Li, D.; Cao, X.; Teng, Y. Distribution, Genesis, and Human Health Risks of Groundwater Heavy Metals Impacted by the Typical Setting of Songnen Plain of NE China. Int. J. Environ. Res. Public Health 2022, 19, 3571. [Google Scholar] [CrossRef]
  7. Talpur, S.A.; Noonari, T.M.; Rashid, A.; Ahmed, A.; Jat Baloch, M.Y.; Talpur, H.A.; Soomro, M.H. Hydrogeochemical signatures and suitability assessment of groundwater with elevated fluoride in unconfined aquifers Badin district, Sindh, Pakistan. SN Appl. Sci. 2020, 2, 1038. [Google Scholar] [CrossRef]
  8. Xia, X.; Teng, Y.; Zhai, Y.; Zheng, F.; Cao, X. Influencing factors and mechanism by which DOM in groundwater releases Fe from sediment. Chemosphere 2022, 300, 134524. [Google Scholar] [CrossRef]
  9. Zhai, Y.; Cao, X.; Xia, X.; Wang, B.; Teng, Y.; Li, X. Elevated Fe and Mn Concentrations in Groundwater in the Songnen Plain, Northeast China, and the Factors and Mechanisms Involved. Agronomy 2021, 11, 2392. [Google Scholar] [CrossRef]
  10. Palma-Lara, I.; Martínez-Castillo, M.; Quintana-Pérez, J.; Arellano-Mendoza, M.; Tamay-Cach, F.; Valenzuela-Limón, O.; García-Montalvo, E.; Hernández-Zavala, A. Arsenic exposure: A public health problem leading to several cancers. Regul. Toxicol. Pharmacol. 2020, 110, 104539. [Google Scholar] [CrossRef]
  11. Sodhi, K.K.; Kumar, M.; Agrawal, P.K.; Singh, D.K. Perspectives on arsenic toxicity, carcinogenicity and its systemic remediation strategies. Environ. Technol. Innov. 2019, 16, 100462. [Google Scholar] [CrossRef]
  12. Li, P.; Karunanidhi, D.; Subramani, T.; Srinivasamoorthy, K. Sources and consequences of groundwater contamination. Arch. Environ. Contam. Toxicol. 2021, 80, 1–10. [Google Scholar] [CrossRef]
  13. Muhammad Yousuf Jat, B.; Chun, S.; Shakeel Ahmed, T.; Javed, I.; Kulvinder, B. Arsenic Removal from Groundwater using Iron Pyrite: Influences Factors and Removal Mechanism. J. Earth Sci. 2022, 32, 1–39. [Google Scholar]
  14. Rahaman, M.S.; Rahman, M.M.; Mise, N.; Sikder, M.T.; Ichihara, G.; Uddin, M.K.; Kurasaki, M.; Ichihara, S. Environmental arsenic exposure and its contribution to human diseases, toxicity mechanism and management. Environ. Pollut. 2021, 289, 117940. [Google Scholar] [CrossRef]
  15. Rehman, U.; Khan, S.; Muhammad, S. Ingestion of arsenic-contaminated drinking water leads to health risk and traces in human biomarkers (hair, nails, blood, and urine), Pakistan. Expo. Health 2020, 12, 243–254. [Google Scholar] [CrossRef]
  16. Smith, E.; Smith, J.; Smith, L.; Biswas, T.; Correll, R.; Naidu, R. Arsenic in Australian environment: An overview. J. Environ. Sci. Health Part A 2003, 38, 223–239. [Google Scholar] [CrossRef]
  17. Ali, W.; Aslam, M.W.; Feng, C.; Junaid, M.; Ali, K.; Li, S.; Chen, Z.; Yu, Z.; Rasool, A.; Zhang, H.; et al. Unraveling prevalence and public health risks of arsenic, uranium and co-occurring trace metals in groundwater along riverine ecosystem in Sindh and Punjab, Pakistan. Environ. Geochem. Health 2019, 41, 2223–2238. [Google Scholar] [CrossRef]
  18. Hernández-Flores, H.; Pariona, N.; Herrera-Trejo, M.; Hdz-García, H.M.; Mtz-Enriquez, A.I. Concrete/maghemite nanocomposites as novel adsorbents for arsenic removal. J. Mol. Struct. 2018, 1171, 9–16. [Google Scholar] [CrossRef]
  19. Aftabtalab, A.; Rinklebe, J.; Shaheen, S.M.; Niazi, N.K.; Moreno-Jiménez, E.; Schaller, J.; Knorr, K.-H. Review on the interactions of arsenic, iron (oxy)(hydr) oxides, and dissolved organic matter in soils, sediments, and groundwater in a ternary system. Chemosphere 2022, 286, 131790. [Google Scholar] [CrossRef] [PubMed]
  20. Savage, K.S.; Bird, D. Inverse geochemical modeling of pit lake evolution in a high-arsenic, alkaline pit lake. In Proceedings of the GSA Abstract Program, Denver, CO, USA, 27–30 October 2002. [Google Scholar]
  21. Dai, Z.; Samper, J.; Ritzi, R., Jr. Identifying geochemical processes by inverse modeling of multicomponent reactive transport in the Aquia aquifer. Geosphere 2006, 2, 210–219. [Google Scholar] [CrossRef]
  22. Shakoor, M.B.; Bibi, I.; Niazi, N.K.; Shahid, M.; Nawaz, M.F.; Farooqi, A.; Naidu, R.; Rahman, M.M.; Murtaza, G.; Lüttge, A. The evaluation of arsenic contamination potential, speciation and hydrogeochemical behaviour in aquifers of Punjab, Pakistan. Chemosphere 2018, 199, 737–746. [Google Scholar] [CrossRef]
  23. Ali, W.; Rasool, A.; Junaid, M.; Zhang, H. A comprehensive review on current status, mechanism, and possible sources of arsenic contamination in groundwater: A global perspective with prominence of Pakistan scenario. Environ. Geochem. Health 2019, 41, 737–760. [Google Scholar] [CrossRef]
  24. Shahid, M.; Khalid, M.; Dumat, C.; Khalid, S.; Niazi, N.K.; Imran, M.; Bibi, I.; Ahmad, I.; Hammad, H.M.; Tabassum, R.A.; et al. Arsenic level and risk assessment of groundwater in Vehari, Punjab Province, Pakistan. Expo. Health 2018, 10, 229–239. [Google Scholar] [CrossRef]
  25. Tabassum, R.A.; Shahid, M.; Dumat, C.; Niazi, N.K.; Khalid, S.; Shah, N.S.; Imran, M.; Khalid, S.J. Health risk assessment of drinking arsenic-containing groundwater in Hasilpur, Pakistan: Effect of sampling area, depth, and source. Environ. Sci. Pollut. Res. 2019, 26, 20018–20029. [Google Scholar] [CrossRef]
  26. Waqas, H.; Shan, A.; Khan, Y.G.; Nawaz, R.; Rizwan, M.; Rehman, M.S.-U.; Shakoor, M.B.; Ahmed, W.; Jabeen, M. Human health risk assessment of arsenic in groundwater aquifers of Lahore, Pakistan. Hum. Ecol. Risk Assess. Int. J. 2017, 23, 836–850. [Google Scholar] [CrossRef]
  27. Islam, M.; Reza, A.; Sattar, G.S.; Ahsan, M.; Akbor, M.; Siddique, M.; Bakar, A. Distribution of arsenic in core sediments and groundwater in the Chapai Nawabganj district, Bangladesh. Arab. J. Geosci. 2019, 12, 99. [Google Scholar] [CrossRef]
  28. Zhang, Y.; Xu, B.; Guo, Z.; Han, J.; Li, H.; Jin, L.; Chen, F.; Xiong, Y. Human health risk assessment of groundwater arsenic contamination in Jinghui irrigation district, China. J. Environ. Manag. 2019, 237, 163–169. [Google Scholar] [CrossRef]
  29. Embaby, A.; Redwan, M.J. Sources and behavior of trace elements in groundwater in the South Eastern Desert, Egypt. Environ. Monit. Assess. 2019, 191, 686. [Google Scholar] [CrossRef]
  30. Bhattacharya, P.; Welch, A.H.; Stollenwerk, K.G.; McLaughlin, M.J.; Bundschuh, J.; Panaullah, G.J. Arsenic in the Environment: Biology and Chemistry; Elsevier: Amsterdam, The Netherlands, 2007; pp. 109–120. [Google Scholar]
  31. Rochaddi, B.; Atmodjo, W.; Satriadi, A.; Suryono, C.A.; Irwani, I.; Widada, S. The Heavy Metal Contamination in Shallow Groundwater at Coastal Areas of Surabaya East Java Indonesia. J. Kelaut. Trop. 2019, 22, 69–72. [Google Scholar] [CrossRef]
  32. Mueller, B.; Hug, S.J. Climatic variations and de-coupling between arsenic and iron in arsenic contaminated ground water in the lowlands of Nepal. Chemosphere 2018, 210, 347–358. [Google Scholar] [CrossRef]
  33. Canencia, O.P.; Walag, A.M.P. Coal combustion from power plant industry in Misamis Oriental, Philippines: A potential groundwater contamination and heavy metal detection. Asian J. Microbiol. Biotechnol. Environ. Sci. 2016, 18, 55–59. [Google Scholar]
  34. Fida, M.; Li, P.; Wang, Y.; Alam, S.; Nsabimana, A.J.E. Water Contamination and Human Health Risks in Pakistan: A Review. Exposure Health 2022. [Google Scholar] [CrossRef]
  35. 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] [Green Version]
  36. Malik, A.H.; Khan, Z.M.; Mahmood, Q.; Nasreen, S.; Bhatti, Z. Perspectives of low cost arsenic remediation of drinking water in Pakistan and other countries. J. Hazard. Mater. 2009, 168, 1–12. [Google Scholar] [CrossRef]
  37. Rabbani, U.; Mahar, G.; Siddique, A.; Fatmi, Z. Risk assessment for arsenic-contaminated groundwater along River Indus in Pakistan. Environ. Geochem. Health 2017, 39, 179–190. [Google Scholar] [CrossRef]
  38. Abbas, Z.; Mapoma, H.W.T.; Su, C.; Aziz, S.Z.; Ma, Y.; Abbas, N. Spatial analysis of groundwater suitability for drinking and irrigation in Lahore, Pakistan. Environ. Monit. Assess. 2018, 190, 1–16. [Google Scholar] [CrossRef]
  39. Farooqi, A.; Masuda, H.; Firdous, N. Toxic fluoride and arsenic contaminated groundwater in the Lahore and Kasur districts, Punjab, Pakistan and possible contaminant sources. Environ. Pollut. 2007, 145, 839–849. [Google Scholar] [CrossRef]
  40. Nickson, R.; McArthur, J.; Shrestha, B.; Kyaw-Myint, T.; Lowry, D. Arsenic and other drinking water quality issues, Muzaffargarh District, Pakistan. Appl. Geochem. 2005, 20, 55–68. [Google Scholar] [CrossRef]
  41. Malana, M.A.; Khosa, M.A. Groundwater pollution with special focus on arsenic, Dera Ghazi Khan-Pakistan. J. Saudi Chem. Soc. 2011, 15, 39–47. [Google Scholar] [CrossRef] [Green Version]
  42. Baig, J.A.; Kazi, T.G.; Arain, M.B.; Afridi, H.I.; Kandhro, G.A.; Sarfraz, R.A.; Jamal, M.K.; Shah, A.Q. Evaluation of arsenic and other physico-chemical parameters of surface and ground water of Jamshoro, Pakistan. J. Hazard. Mater. 2009, 166, 662–669. [Google Scholar] [CrossRef]
  43. Arain, M.; Kazi, T.; Baig, J.; Jamali, M.; Afridi, H.; Shah, A.; Jalbani, N.; Sarfraz, R. Determination of arsenic levels in lake water, sediment, and foodstuff from selected area of Sindh, Pakistan: Estimation of daily dietary intake. Food Chem. Toxicol. 2009, 47, 242–248. [Google Scholar] [CrossRef] [PubMed]
  44. Brahman, K.D.; Kazi, T.G.; Afridi, H.I.; Naseem, S.; Arain, S.S.; Ullah, N. Evaluation of high levels of fluoride, arsenic species and other physicochemical parameters in underground water of two sub districts of Tharparkar, Pakistan: A multivariate study. Water Res. 2013, 47, 1005–1020. [Google Scholar] [CrossRef] [PubMed]
  45. Arain, M.; Kazi, T.; Jamali, M.; Jalbani, N.; Afridi, H.; Shah, A. Total dissolved and bioavailable elements in water and sediment samples and their accumulation in Oreochromis mossambicus of polluted Manchar Lake. Chemosphere 2008, 70, 1845–1856. [Google Scholar] [CrossRef] [PubMed]
  46. Shahid, M.; Niazi, N.K.; Dumat, C.; Naidu, R.; Khalid, S.; Rahman, M.M.; Bibi, I. A meta-analysis of the distribution, sources and health risks of arsenic-contaminated groundwater in Pakistan. Environ. Pollut. 2018, 242, 307–319. [Google Scholar] [CrossRef]
  47. Fatima, S.; Hussain, I.; Rasool, A.; Xiao, T.; Farooqi, A. Comparison of two alluvial aquifers shows the probable role of river sediments on the release of arsenic in the groundwater of district Vehari, Punjab, Pakistan. Environ. Earth Sci. 2018, 77, 382. [Google Scholar] [CrossRef]
  48. Memon, M.; Soomro, M.S.; Akhtar, M.S.; Memon, K.S. Drinking water quality assessment in Southern Sindh (Pakistan). Environ. Monit. Assess. 2011, 177, 39–50. [Google Scholar] [CrossRef]
  49. Podgorski, J.E.; Eqani, S.A.M.A.S.; Khanam, T.; Ullah, R.; Shen, H.; Berg, M. Extensive arsenic contamination in high-pH unconfined aquifers in the Indus Valley. Sci. Adv. 2017, 3, e1700935. [Google Scholar] [CrossRef] [Green Version]
  50. Shah, A.H.; Shahid, M.; Khalid, S.; Shabbir, Z.; Bakhat, H.F.; Murtaza, B.; Farooq, A.; Akram, M.; Shah, G.M.; Nasim, W. Assessment of arsenic exposure by drinking well water and associated carcinogenic risk in peri-urban areas of Vehari, Pakistan. Environ. Geochem. Health 2020, 42, 121–133. [Google Scholar] [CrossRef]
  51. Murtaza, B.; Amjad, M.; Shahid, M.; Imran, M.; Shah, N.S.; Abbas, G.; Naeem, M.A. Compositional and health risk assessment of drinking water from health facilities of District Vehari, Pakistan. Environ. Geochem. Health 2020, 42, 2425–2437. [Google Scholar] [CrossRef]
  52. Machado, I.; Falchi, L.; Bühl, V.; Mañay, N. Arsenic levels in groundwater and its correlation with relevant inorganic parameters in Uruguay: A medical geology perspective. Sci. Total Environ. 2020, 721, 137787. [Google Scholar] [CrossRef]
  53. Ahmad, N.; Ahmad, M.; Rafiq, M.; Iqbal, N.; Ali, M.; Sajjad, M.I. Hydrological modeling of the Lahore-Aquifer, using isotopic, chemical and numerical techniques. Back Issues J. Sci. Vis. 2002, 7, 16. [Google Scholar]
  54. Water Environmental Federation; APHA Association. Standard Methods for the Examination of Water and Wastewater; American Public Health Association: Washington, DC, USA, 2005; 21p. [Google Scholar]
  55. Piper, A.M.J.E. Transactions American Geophysical Union, A graphic procedure in the geochemical interpretation of water-analyses. Eos Trans. Am. Geophys. Union 1944, 25, 914–928. [Google Scholar] [CrossRef]
  56. United States Environmental Protection Agency. Arsenic, Inorganic; Integrated Risk Information System (IRIS), CASRN 7440-38-2). 1998. Available online: http://www.epa.gov/iris (accessed on 2 August 2022).
  57. USEPA. Guidelines for carcinogen risk assessment. In Risk Assessment Forum; EPA/630/P-03 F.; US Environmental Protection Agency: Washington, DC, USA, 2005. [Google Scholar]
  58. Parkhurst, D.L.; Appelo, C.A.J. User’s guide to PHREEQC (Version 2): A computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations. Water-Resour. Investig. Rep. 1999, 99, 312. [Google Scholar]
  59. Sharif, M.; Davis, R.; Steele, K.; Kim, B.; Kresse, T.; Fazio, J.A. Inverse geochemical modeling of groundwater evolution with emphasis on arsenic in the Mississippi River Valley alluvial aquifer, Arkansas (USA). J. Hydrol. 2008, 350, 41–55. [Google Scholar] [CrossRef]
  60. WHO. Guidelines for Drinking-Water Quality: Incorporating the First and Second Addenda; WHO: Geneva, Switzerland, 2022. [Google Scholar]
  61. Adimalla, N.; Qian, H.; Tiwari, D.M. Groundwater chemistry, distribution and potential health risk appraisal of nitrate enriched groundwater: A case study from the semi-urban region of South India. Ecotoxicol. Environ. Saf. 2021, 207, 111277. [Google Scholar] [CrossRef]
  62. Smedley, P.L.; Kinniburgh, D.G. A review of the source, behaviour and distribution of arsenic in natural waters. Appl. Geochem. 2002, 17, 517–568. [Google Scholar] [CrossRef] [Green Version]
  63. Rahman, M.M.; Naidu, R.; Bhattacharya, P. Arsenic contamination in groundwater in the Southeast Asia region. Environ. Geochem. Health 2009, 31, 9–21. [Google Scholar] [CrossRef]
  64. Rasool, A.; Xiao, T.; Farooqi, A.; Shafeeque, M.; Masood, S.; Ali, S.; Fahad, S.; Nasim, W. Arsenic and heavy metal contaminations in the tube well water of Punjab, Pakistan and risk assessment: A case study. Ecol. Eng. 2016, 95, 90–100. [Google Scholar] [CrossRef]
  65. Iqbal, J.; Su, C.; Rashid, A.; Yang, N.; Baloch, M.Y.J.; Talpur, S.A.; Ullah, Z.; Rahman, G.; Rahman, N.U.; Sajjad, M.M. Hydrogeochemical assessment of groundwater and suitability analysis for domestic and agricultural utility in Southern Punjab, Pakistan. Water 2021, 13, 3589. [Google Scholar] [CrossRef]
  66. Singh, R.; Singh, S.; Parihar, P.; Singh, V.P.; Prasad, S.M. Arsenic contamination, consequences and remediation techniques: A review. Ecotoxicol. Environ. Saf. 2015, 112, 247–270. [Google Scholar] [CrossRef]
  67. Baloch, M.Y.J.; Talpur, S.A.; Talpur, H.A.; Iqbal, J.; Mangi, S.H.; Memon, S.J. Effects of Arsenic Toxicity on the Environment and Its Remediation Techniques: A Review. J. Water Environ. Technol. 2020, 18, 275–289. [Google Scholar] [CrossRef]
  68. Xing, S.; Guo, H.; Zhang, L.; Wang, Z.; Sun, X.J. Silicate weathering contributed to arsenic enrichment in geotherm-affected groundwater in Pliocene aquifers of the Guide basin, China. J. Hydrol. 2022, 606, 127444. [Google Scholar] [CrossRef]
  69. Wen, X.; Wu, Y.; Wu, J.J.E.G. Hydrochemical characteristics of groundwater in the Zhangye Basin, Northwestern China. Environ. Earth Sci. 2008, 55, 1713–1724. [Google Scholar] [CrossRef]
  70. Anawar, H.M.; Akai, J.; Sakugawa, H.J.C. Mobilization of arsenic from subsurface sediments by effect of bicarbonate ions in groundwater. Chemosphere 2004, 54, 753–762. [Google Scholar] [CrossRef]
  71. Bhowmick, S.; Nath, B.; Halder, D.; Biswas, A.; Majumder, S.; Mondal, P.; Chakraborty, S.; Nriagu, J.; Bhattacharya, P.; Iglesias, M.J. Arsenic mobilization in the aquifers of three physiographic settings of West Bengal, India: Understanding geogenic and anthropogenic influences. J. Hazard. Mater. 2013, 262, 915–923. [Google Scholar] [CrossRef]
  72. Hasan, M.A.; Bhattacharya, P.; Sracek, O.; Ahmed, K.M.; Von Brömssen, M.; Jacks, G.J. Geological controls on groundwater chemistry and arsenic mobilization: Hydrogeochemical study along an E–W transect in the Meghna basin, Bangladesh. J. Hydrol. 2009, 378, 105–118. [Google Scholar] [CrossRef]
  73. Jat Baloch, M.Y.; Zhang, W.; Chai, J.; Li, S.; Alqurashi, M.; Rehman, G.; Tariq, A.; Talpur, S.A.; Iqbal, J.; Munir, M. Shallow Groundwater Quality Assessment and Its Suitability Analysis for Drinking and Irrigation Purposes. Water 2021, 13, 3361. [Google Scholar] [CrossRef]
  74. Sultana, J.; Farooqi, A.; Ali, U.J. Arsenic concentration variability, health risk assessment, and source identification using multivariate analysis in selected villages of public water system, Lahore, Pakistan. Environ. Monit. Assess. 2014, 186, 1241–1251. [Google Scholar] [CrossRef]
  75. Bibi, S.; Farooqi, A.; Ramzan, M.; Javed, A. Health risk of arsenic in the alluvial aquifers of Lahore and Raiwind, Punjab Province, Pakistan: An investigation for safer well water. Toxicol. Environ. Chem. 2015, 97, 888–907. [Google Scholar] [CrossRef]
  76. Rinklebe, J.; Shaheen, S.M.; Yu, K.J.G. Release of As, Ba, Cd, Cu, Pb, and Sr under pre-definite redox conditions in different rice paddy soils originating from the USA and Asia. Geoderma 2016, 270, 21–32. [Google Scholar] [CrossRef]
  77. Rasool, A.; Xiao, T.; Farooqi, A.; Shafeeque, M.; Liu, Y.; Kamran, M.A.; Katsoyiannis, I.A.; Eqani, S.A.M.A.S. Quality of tube well water intended for irrigation and human consumption with special emphasis on arsenic contamination at the area of Punjab, Pakistan. Environ. Geochem. Health 2017, 39, 847–863. [Google Scholar] [CrossRef] [PubMed]
  78. Bauer, M.; Blodau, C. Mobilization of arsenic by dissolved organic matter from iron oxides, soils and sediments. Sci. Total Environ. 2006, 354, 179–190. [Google Scholar] [CrossRef] [PubMed]
  79. Lengke, M.F.; Sanpawanitchakit, C.; Tempel, R.N.J.T.C.M. The oxidation and dissolution of arsenic-bearing sulfides. Can. Miner. 2009, 47, 593–613. [Google Scholar] [CrossRef]
  80. Anawar, H.M.; Akai, J.; Komaki, K.; Terao, H.; Yoshioka, T.; Ishizuka, T.; Safiullah, S.; Kato, K.J. Geochemical occurrence of arsenic in groundwater of Bangladesh: Sources and mobilization processes. J. Geochem. Explor. 2003, 77, 109–131. [Google Scholar] [CrossRef]
  81. Cui, J.; Jing, C. A review of arsenic interfacial geochemistry in groundwater and the role of organic matter. Ecotoxicol. Environ. Saf. 2019, 183, 109550. [Google Scholar] [CrossRef]
  82. Wang, Y.; Zhang, G.; Wang, H.; Cheng, Y.; Liu, H.; Jiang, Z.; Li, P.; Wang, Y.J. Effects of different dissolved organic matter on microbial communities and arsenic mobilization in aquifers. J. Hazard. Mater. 2021, 411, 125146. [Google Scholar] [CrossRef]
  83. Casentini, B.; Pettine, M.; Millero, F.J. Release of arsenic from volcanic rocks through interactions with inorganic anions and organic ligands. Aquat. Geochem. 2010, 16, 373–393. [Google Scholar] [CrossRef]
  84. Arain, G.M.; Aslam, M.; Majidano, S.A.; Khuhawar, M.Y. A preliminary study on the arsenic contamination of underground water of Matiari and Khairpur Districts, Sindh, Pakistan. J. -Chem. Soc. Pak. 2007, 29, 463. [Google Scholar]
  85. Bandara, U.; Diyabalanage, S.; Hanke, C.; Van Geldern, R.; Barth, J.A.; Chandrajith, R. Arsenic-rich shallow groundwater in sandy aquifer systems buffered by rising carbonate waters: A geochemical case study from Mannar Island, Sri Lanka. Sci. Total Environ. 2018, 633, 1352–1359. [Google Scholar] [CrossRef]
  86. Kumar, A.; Singh, C.K. Arsenic enrichment in groundwater and associated health risk in Bari doab region of Indus basin, Punjab, India. Environ. Pollut. 2020, 256, 113324. [Google Scholar] [CrossRef]
  87. Rasool, A.; Farooqi, A.; Xiao, T.; Masood, S.; Kamran, M.A. Elevated levels of arsenic and trace metals in drinking water of Tehsil Mailsi, Punjab, Pakistan. J. Geochem. Explor. 2016, 169, 89–99. [Google Scholar] [CrossRef]
  88. Shakoor, M.B.; Niazi, N.K.; Bibi, I.; Murtaza, G.; Kunhikrishnan, A.; Seshadri, B.; Shahid, M.; Ali, S.; Bolan, N.S.; Ok, Y.S. Remediation of arsenic-contaminated water using agricultural wastes as biosorbents. Crit. Rev. Environ. Sci. Technol. 2016, 46, 467–499. [Google Scholar] [CrossRef]
  89. Muhammad, S.; Shah, M.T.; Khan, S. Arsenic health risk assessment in drinking water and source apportionment using multivariate statistical techniques in Kohistan region, northern Pakistan. Food Chem. Toxicol. 2010, 48, 2855–2864. [Google Scholar] [CrossRef]
  90. Farooqi, A.; Masuda, H.; Kusakabe, M.; Naseem, M.; Firdous, N. Distribution of highly arsenic and fluoride contaminated groundwater from east Punjab, Pakistan, and the controlling role of anthropogenic pollutants in the natural hydrological cycle. Geochem. J. 2007, 41, 213–234. [Google Scholar] [CrossRef] [Green Version]
  91. Karim, M.M. Arsenic in groundwater and health problems in Bangladesh. Water Research 2000, 34, 304–310. [Google Scholar] [CrossRef]
  92. Caylak, E. Health risk assessment for arsenic in water sources of Cankiri Province of Turkey. Clean–Soil Air Water 2012, 40, 728–734. [Google Scholar] [CrossRef]
Figure 1. (a) Groundwater sampling points of the study area. (b) Flow and hydrogeology map.
Figure 1. (a) Groundwater sampling points of the study area. (b) Flow and hydrogeology map.
Ijerph 19 13325 g001aIjerph 19 13325 g001b
Figure 2. Piper diagram showing groundwater chemistry in study area.
Figure 2. Piper diagram showing groundwater chemistry in study area.
Ijerph 19 13325 g002
Figure 3. Scattered ionic plots. (a) [(Ca2+ + Mg2+) − (HCO3 + SO42)] vs. [(Na+ + K)−Cl]. (b) [(Ca2+ + Mg2+ vs. HCO3 + SO42)]. (c) [(HCO3/Na+) vs. (Ca2+/Na+)]. (d) [(Mg2+/Na+) vs. (Ca2+/Na+)].
Figure 3. Scattered ionic plots. (a) [(Ca2+ + Mg2+) − (HCO3 + SO42)] vs. [(Na+ + K)−Cl]. (b) [(Ca2+ + Mg2+ vs. HCO3 + SO42)]. (c) [(HCO3/Na+) vs. (Ca2+/Na+)]. (d) [(Mg2+/Na+) vs. (Ca2+/Na+)].
Ijerph 19 13325 g003
Figure 4. Relationships between As and saturation indices in groundwater. (a) Vehari, (b) Burewala, and (c) Mailsi.
Figure 4. Relationships between As and saturation indices in groundwater. (a) Vehari, (b) Burewala, and (c) Mailsi.
Ijerph 19 13325 g004
Figure 5. (a) Sources of As mobilization. (b) Spatial distribution of groundwater As in study area.
Figure 5. (a) Sources of As mobilization. (b) Spatial distribution of groundwater As in study area.
Ijerph 19 13325 g005aIjerph 19 13325 g005b
Figure 6. (a) Carcinogenic risk (CR) of adults and children. (b) Non-carcinogenic risk hazard quotient (HQ) of adults and children.
Figure 6. (a) Carcinogenic risk (CR) of adults and children. (b) Non-carcinogenic risk hazard quotient (HQ) of adults and children.
Ijerph 19 13325 g006
Table 1. Statistical physicochemical parameters of groundwater in three tehsils. All parameters are in (mg/L), except for pH, As (µg/L), EC (µS/cm), and turbidity (NTU). Below detection limit (BDL).
Table 1. Statistical physicochemical parameters of groundwater in three tehsils. All parameters are in (mg/L), except for pH, As (µg/L), EC (µS/cm), and turbidity (NTU). Below detection limit (BDL).
ParametersVehari (n = 170)Burewala (n = 170)Mailsi (n = 170)WHO (2022) Standard
Min–MaxMean ±
Standard
Min–MaxMean ±
Standard
Min–MaxMean ±
Standard
EC308–45501569.94 ± 848.6385–44001230.78 ± 646.37226–26901185.059 ± 468.991000
TDS234–31481072.12 ± 514.8277–3173896.68 ± 429.54359–3298896.64 ± 385.681000
pH6.85–7.617.17 ± 0.176.78–7.156.98 ± 0.077.2–8.357.66 ± 0.236.5–8.5
Turbidity0.3–99126.37 ± 86.340.3–20215.98 ± 38.570.2–4.21.22 ± 1.124.0
Alkalinity2.2–61010.82 ± 46.562.8–13.25.98 ± 1.811.4–12.86.50 ± 1.68-
Hardness100–820368.47 ± 132.94105–820348.85 ± 123.71190–1540414.58 ± 152.34-
HCO3110–900345.08 ± 116.3580–660299.20 ± 92.7870–640324.23 ± 84.59-
Cl10–518107.94 ± 88.8310–50274.08 ± 57.828–336103.51 ± 53.84200–300
SO42–41–1300316.87 ± 219.8118–1432244.24 ± 185.3618–840167.49 ± 126.37250
Ca2+24–19286.31 ± 35.4026–20886.44 ± 32.598–32086.91 ± 40.77200
Mg2+6–9836.65 ± 16.5510–9532.25 ± 14.522–18048 ± 21.55150
Na+14–850195.29 ± 157.1713–620131.78 ± 109.0821–360107.97 ± 68.98200
K+2.6–42.69.53 ± 6.533.7–698.085 ± 6.241–749.49 ± 8.2912
NO3BDL-17.661.51 ± 3.070.01–15.821.8 ± 3.31BDL-141.85 ± 2.9450
FBDL-3.150.59 ± 0.370.18–1.350.46 ± 0.200.22–1.240.51 ± 0.171.5
Fe2+0.01–2.940.18 ± 0.390.01–3.920.19 ± 0.430.01–2.940.23 ± 0.400.3
AsBDL-4514.0 ± 11.85BDL-5211.0 ± 11.62BDL-8913.0 ± 15.310
Table 2. Factor loading for groundwater physicochemical parameters in the study area.
Table 2. Factor loading for groundwater physicochemical parameters in the study area.
TehsilVehariBurewalaMailsi
ComponentPC1PC2PC3PC1PC2PC3PC1PC2PC3
EC0.964−0.1990.0920.957−0.094−0.1850.8690.0630.195
TDS0.9810.0510.1310.9620.141−0.160.972−0.0740.027
pH0.069−0.0620.376−0.0370.125−0.14−0.340.2080.512
Turbidity0.0170.2450.329−0.0520.3520.0160.2170.7450.086
Alkalinity0.072−0.2020.0140.837−0.2150.1760.8140.191−0.222
HCO30.707−0.067−0.1080.821−0.2540.1820.8230.178−0.209
Cl0.867−0.2070.0970.844−0.057−0.1140.868−0.0480.132
SO42–0.936−0.1110.1930.8560.038−0.3560.862−0.0810.28
Ca2+0.3920.8310.1570.5920.7060.1730.779−0.334−0.025
Mg2+0.8160.301−0.1550.8330.1280.1930.6170.206−0.502
Na+0.856−0.4330.1210.859−0.307−0.3320.7920.130.383
K+0.3310.108−0.5560.301−0.1240.4810.573−0.0970.291
Hardness0.6620.710.0260.7930.5270.2080.88−0.104−0.307
NO30.125−0.064−0.6370.328−0.0060.6480.322−0.5020.062
Fe0.1540.093−0.0720.156−0.120.039−0.040.307−0.002
F0.466−0.643−0.1530.517−0.459−0.1480.2250.3030.313
As−0.288−0.2570.514−0.0020.663−0.3780.0250.783−0.12
Eigen Value6.462.181.407.561.851.327.561.921.19
% of Variance38.0412.838.2744.4910.887.7744.5111.337.01
Cumulative %38.0450.8859.1544.4955.3763.1544.5155.8562.86
Extraction method: principal component analysis. Rotation method: Kaiser normalization of Varimax. Bold values show higher loading value.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jat Baloch, M.Y.; Zhang, W.; Zhang, D.; Al Shoumik, B.A.; Iqbal, J.; Li, S.; Chai, J.; Farooq, M.A.; Parkash, A. Evolution Mechanism of Arsenic Enrichment in Groundwater and Associated Health Risks in Southern Punjab, Pakistan. Int. J. Environ. Res. Public Health 2022, 19, 13325. https://doi.org/10.3390/ijerph192013325

AMA Style

Jat Baloch MY, Zhang W, Zhang D, Al Shoumik BA, Iqbal J, Li S, Chai J, Farooq MA, Parkash A. Evolution Mechanism of Arsenic Enrichment in Groundwater and Associated Health Risks in Southern Punjab, Pakistan. International Journal of Environmental Research and Public Health. 2022; 19(20):13325. https://doi.org/10.3390/ijerph192013325

Chicago/Turabian Style

Jat Baloch, Muhammad Yousuf, Wenjing Zhang, Dayi Zhang, Baig Abdullah Al Shoumik, Javed Iqbal, Shuxin Li, Juanfen Chai, Muhammad Ansar Farooq, and Anand Parkash. 2022. "Evolution Mechanism of Arsenic Enrichment in Groundwater and Associated Health Risks in Southern Punjab, Pakistan" International Journal of Environmental Research and Public Health 19, no. 20: 13325. https://doi.org/10.3390/ijerph192013325

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