A hydrogeochemical approach to evaluate groundwater quality in the vicinity of three tributaries of the Beas River, North-West India

The present study focused on the seasonal investigation of hydro-geochemical characteristics of groundwater samples collected from the vicinity of three tributaries of the Beas River, Punjab, India. Total 45 samples were analyzed during the pre- and post-monsoon season for physico-chemical parameters and heavy metals along with health risk assessment. Results revealed that the majority of samples were below the permissible limits set by the BIS and WHO. The relative abundance of major cations was Ca2+ > Mg2+ > Na+ > K+ and Ca2+ > Na+ > Mg2+ > K+, while that of the major anions was HCO3− > SO42− > Cl− > CO32− in the pre- and post-monsoon season, respectively. Groundwater was alkaline and hard in nature at most of the sites. Bicarbonate content exceeded the desirable limit having an average concentration of 337.26 mg/L and 391.48 mg/L, respectively, during the pre- and post-monsoon season. Tukey’s multiple comparison test was applied for finding significant differences among samples at p < 0.05. The dominant hydrochemical face of water was Ca–Mg–HCO3 type. US salinity (USSL) diagram indicated that during the pre-monsoon, 48.9% samples were C2S1 type and remaining 51.1% were C3S1 type while during the post-monsoon all samples were C3S1 type. It indicates that groundwater of the study area is at risk of salinity hazards in future and is not to be ignored. Such monitoring studies are recommended to design future safety plans to combat soil and human health risks.


Introduction
Groundwater is the most important source for drinking and irrigation purpose but due to intensive farming, industrial discharge, pollution of rivers, its quality is deteriorating in North-West regions of Punjab, India for survival (Talabi and Kayode 2019;Adimalla et al. 2020). The Beas River, falling in the Northwest of Punjab, India and its tributaries adjoining agricultural fields loaded with contaminants like heavy metals, bicarbonates, chlorides, sulfates etc. due to inefficient agricultural practices Srivastava et al. 2019;Setia et al. 2020). Use of groundwater is very high in India, particularly for irrigation purposes . Indiscreet consumption and speedy rise of this contamination in groundwater is threat to supply of good quality water (Kumar and Shah 2006;Gupta et al. 2018;Jain and Mujawar 2019). Groundwater contamination is difficult to assess than that of surface water. It needs regular monitoring as contamination will remain for decades in the subsurface environment (Jha et al. 2020). India is the biggest consumer of groundwater in the world, with an estimated groundwater use of around 230 km 3 /year (CGWB 2017). The higher use has resulted in depletion of groundwater more prominently in the Punjab state (Krishan et al. 2014). Groundwater is intimately linked to surface water present in the rivers and water received as rainfall also contributes to groundwater reserves. So, pollution of surface water also affects the quality of groundwater (Srivastava et al. 2019).
About 80% of diseases in human beings are directly linked to consumption of polluted water (Ganiyu et al. 2018;Ali and Ahmad 2019). The heavy metals are also one of the major contaminants of water. Exposure to contaminants in groundwater via drinking water and cooked food had led to direct or indirect accumulation in human body resulting in health hazards like respiratory diseases (Jaishankar et al. 5 Page 2 of 18 2014), osteoporosis (Rodriguez and Mandalunis 2018), hypertension (Scheelbeek et al. 2017), dermal problems (Sharma and Bhattacharya 2017), neurological (Bondy and Campbell 2017) and reproductive disorders (Agrawal et al. 2010), etc. According to the data compiled by the Central Groundwater Board (CGWB 2017), the total 2420 number of villages had been found to be affected by heavy metal contamination in India, out of which 2139 had fallen in Punjab .The present study was therefore undertaken to investigate the quality of groundwater by determining the various hydro-geochemical characteristics of groundwater in the vicinity of the three tributaries of Beas River falling in Bist-Doab region, namely Pong left main canal, Odhra rivulet and Kali bein and determining the suitability of groundwater for drinking and irrigation purpose.

Sampling and sample pretreatment
An array of 15 villages (5 villages/tributary) distributed across the study area was selected which are close to confluence of the selected tributaries of the Beas river, Punjab, India (Fig. 1). Total 90 groundwater samples (45 samples through each season) were collected from tubewells located in these villages (3 sampling sites i.e., 45 × 3 = 90) in June 2019 (pre-monsoon) and October 2019 (post-monsoon) from same locations to evaluate the seasonal variations. Samples were collected in clean prewashed polyethylene bottles around 3-4 times with same groundwater to be collected. Samples were collected after pumping the tubewells for about 10 min so that water previously stored in pipes can be removed . For estimation of heavy metals such as manganese (Mn), lead (Pb), molybdenum (Mo), nickel (Ni), selenium (Se), aluminum (Al), arsenic (As), copper (Cu), cadmium (Cd), iron (Fe) , cobalt (Co) and zinc (Zn), a portion of water sample was acidified with conc. HNO 3 (pH<2) and rest of the portion were used to analyze different physico-chemical parameters like pH, electrical conductivity (EC), total dissolved solids (TDS) and contents of anions and cations. The groundwater samples were transported in cooled ice box and stored in refrigerator at 5 °C until analysis.

Physico-chemical analysis of groundwater samples
Physico-chemical parameters of collected groundwater samples were analyzed by following the standard methods (Trivedy et al. 1987;APHA (2017); The portable kit was used for measuring pH and TDS at sampling sites. Carbonate (CO 3 2− ), bicarbonate (HCO 3 − ) and chloride (Cl − ) ions were determined by the titrimetric method while calcium (Ca 2+ ) and magnesium (Mg 2+ ) were quantified using versenate method.
First letter in name of each site represents the name of tributary and number represents the code of village as P 1 : Sadarpur, P 2 : Lamipur, P 3 : Jamalpur, P 4 : Bhagran, P 5 :  (Sheen et al. 1935) using colorimeter (make: Systronics; model: 115) at a wavelength of 420 mμ with a blue filter. Na + and K + contents in samples were analyzed using flame photometer (make: Systronics; model: 130).

Heavy metal analysis of groundwater samples
Heavy metals were determined in water samples by following standard method used by Oyekanmi et al. 2018 using microwave plasma-atomic emission spectrometer (make: Agilent; model-4200). Before determination of concentration of heavy metals, water samples were filtered using Whatman filter paper no. 42 .

Reagents and quality control
The chemicals and reagents used for analyzing physicochemical parameters and heavy metal contents in present study were of analytical grade. Stock standard solutions (Agilent technologies, USA) were used to draw calibration curves for determination of heavy metal content in samples. Double distilled water was used at regular intervals to run instrument for washing to avoid the deposition of any analyst inside it.

Water quality assessment
Groundwater samples were categorized into various classes for irrigation and drinking purposes on the basis of physico-chemical parameters such as total dissolved solids (TDS) and total hardness (TH) etc. The hardness of water is due to dissolution of salts of Ca 2+ and Mg 2+ . So, to determine TH, we have estimated Ca 2+ and Mg 2+ contents in water samples. Hem (1989) had given following formula for calculating TH of water: Groundwater samples can also be classified on the basis of indices such as electrical conductivity (EC), sodium adsorption ratio (SAR), residual sodium carbonate (RSC), Kelly's ratio (KR), magnesium ratio (MR), permeability index (PI), corrosivity ratio (CR) and sodium percentage (Na%). Suitability of groundwater samples for drinking and irrigation purposes was determined by keeping in consideration the standard limits given by WHO (2011) andBIS (2012).
Sodium adsorption ratio (SAR) is used for estimating sodium hazard to crops and to classify groundwater to check its suitability for irrigating crop fields (Gholami and Srikantaswamy 2009). SAR values of all groundwater samples were estimated by using the equation given by Richards 1954as: Todd (1995 had given following equation to calculate sodium percentage (Na%) in water: Magnesium ratio (MR) was evaluated using following formula (Palliwal 1972): Kelly's ratio (KR) was calculated by using equation given by Kelley (1946) which is expressed as follows: Residual sodium carbonate (RSC) is the measure of excessive amount of sodium bicarbonate and carbonate over Ca 2+ and Mg 2+ which has affect on physical properties like pH, EC and SAR of the soil by dissolution of organic matter (Naseem et al. 2010). It has been evaluated using the following equation (Ragunath 1987): (1) TH as CaCO 3 mg∕L = Ca 2+ + Mg 2+ meq∕L × 50 (3) Na% = Na + + K + × 100 Corrosivity ratio (CR) is used to check the corrosiveness of water. It was evaluated using following equation (Raman 1985): Permeability index (PI) is a qualitative measure of rate of vertical movement of water from the ground surface through unsaturated zone to the subsurface layers of the soil (Kumar 2016;Rawat et al. 2018). It is given as:

Health risk assessment
Health risk assessment is essential to check the suitability of groundwater for drinking purposes. The US Environmental Protection Agency (USEPA 2004) had given following two equations to measure the exposure dose (E) through ingestion and potential non-carcinogenic risk of Hazard Quotient (HQ): where E: chronic daily intake (mg/kg/day), CPW: concentration of contaminant (mg/L), IR: human ingestion rate (L/ day: 0.78 L/day for children and 2.5 L/day for adults), ED: exposure duration (years: 12 and 64 for children and adults, respectively), EF: exposure frequency (days/years: 365 days for children and adults), ABW: average body weight (Kg: 18.7 and 57.5 for children and adults, respectively), AET: average time (days: 4380 and 23360 for children and adults, respectively), HQ: Hazard quotient, RfD: reference dose of pollutant (Adimalla and Qian 2019).

Statistical analysis
Analysis of variance (ANOVA) was applied on physicochemical parameters and heavy metal content to assess variation among location and season of sampling site at p < 0.05 using SPSS 16.0 software. Further, Tukey's multiple comparison test was used to confirm significant differences among samples at p < 0.05.

Physico-chemical analysis
All the tributaries were in the Bist-Doab region where groundwater level falls systematically in various regions over the last two decades. The pH of the studied samples was within the maximum permissible limits given by WHO and BIS during both the pre-and post-monsoon season (Kamble and Saxena 2016). All the groundwater samples of both seasons were alkaline (average of 7.80 and 7.96) in nature. Dissolution of primary CO 3 − and silicate minerals after monsoon had led to increase in concentration of Ca 2+ , Na + , Mg 2+ and HCO 3 − which had raised the value of pH (Rouabhia et al. 2011). During each season, values of EC remained higher than the permissible value of 250 μS/cm (WHO 2011). All samples have shown higher EC (μS/cm) values than permissible limits, i.e., 750.80 (pre-monsoon) and 891 (post-monsoon). TDS in mg/L of all the samples was observed to be less than the permissible limits set by WHO (2011) andBIS (2012) with average value of 316.09 (pre-monsoon) and 406.18 (post-monsoon). It was noticed a small increase in EC and TDS in the post-monsoon season than in premonsoon season.
It has been reported that leaching and dissolution or mixing of salts due to geochemical or anthropogenic activities was responsible for this variation in the study area Subba Rao et al. 2017ao et al. 2017. During summer, the water present in the interspaces of soil get evaporated and salts were accumulated in upper layer of the soil which were leached back again during monsoon period had resulted high EC and TDS values .
Total hardness (TH) of all the samples in the premonsoon was below the permissible limit but in the postmonsoon season, 31 samples of groundwater have high TH values than desirable value of 200 mg/L as per BIS (2012). It was probably due to increase in Ca 2+ and Mg 2+ ions in the studied samples. The pre-monsoon value (mg/L) of Ca 2+ were ranged from 24.67 to 31.33 having average value 27.56 while in post-monsoon season it was from 31.33 to 41.33 having average value of 36.58. The Mg 2+ concentration (mg/L) in the post-monsoon ranged from 19.60 to 28.40 having average value of 24.00 while in the post-monsoon it was from 22.80 to 29.20 with mean concentration of 26.27. The CO 3 2− content was negligible in both the seasons of the study area while HCO 3 − was found beyond desirable limit in all the 45 samples having average concentration (mg/L) of 337.26 and 391.48 in the pre-and post-monsoon season, respectively. The carbonic acid helps in initiating acid-base reactions and Page 5 of 18 5 fastens the weathering of rocks like limestone (CaCO 3 ) and dolomite [CaMg(CO 3 ) 2 ] present in the study areas and similar results had been reported by Singh et al. 2013;Humbarde et al. 2014;Dev and Bali 2018. Due to this weathering process, the concentration of Ca 2+ and Mg 2+ were also increased after the monsoon (Sharma and Chhipa 2016;Sharma et al. 2019). The abundance of major cations was in the order of Ca 2+ >Mg 2+ >Na + >K + and Ca 2+ > Na + >Mg 2+ >K + , respectively, in the pre-and post-monsoon season. On the other hand, abundance of major anions was in the order of HCO 3 − > SO 4 2− > Cl − > CO 3 2− in both the pre-and post-monsoon seasons as similarly reported by Humbarde et al. 2014;Thakur et al. 2016. The mean concentration of chloride (Cl − ) in the study area was found to be 42.08 and 54.21(mg/L) in the pre-and post-monsoon season, respectively. The high concentration of Cl − could be due to the addition of sewage effluents directly into the tributary water and utilizing that water for irrigation (Saraswat et al. 2019). Leaching of salts along with rainwater was also reported as one of the main reasons behind high chloride concentration (Purushothaman et al. 2012;Pathak et al. 2014). Both the K + and Na + ions were found in concentrations (mg/L) lower than permissible limits having mean value of 11.80 and 12.20 for K + and 24.22 and 26.88 for Na + in the pre-and post-monsoon season, respectively. High concentration of K + was noticed in the post-monsoon season may be due to weathering of saltpetre during the rainy season ). The average value (mg/L) of SO 4 2− was 68.90 in the pre-and 78.85 during the post-monsoon, was lower than the permissible limit. The presence of high SO 4 2− ion during the post-monsoon may be due to leaching, breakdown of organic matter in soil, excess use of fertilizers and anthropogenic activities. Similar findings were documented by Ganiyu et al. 2018 andHejaz et al. 2020 in their studies.

Heavy metal analysis
The heavy metal content in samples of study areas during both the pre-and post-monsoon season is represented in Table 1. The average concentration of Fe, Al, Mn, Cu, As and Se in samples to be 0.30 mg/L, 0.02 mg/L, 0.02 mg/L, 0.03 mg/L, 0.19 μg/L and 0.78 μg/L, respectively during pre-monsoon and 0.33 mg/L, 0.03 mg/L, 0.03 mg/L, 0.03 mg/L, 0.25 μg/L and 0.84 μg/L, respectively, during the post-monsoon. Main sources of heavy metals in the present study area were due to indiscriminate use of agrochemicals and leaching of heavy metals along with the rainwater as reported by Sharma et al. 2019 andVirk 2019. Selenium was found in more concentration mainly at sites near Kali Bein, for the unknown reason but had predicted to be due to Se rich soil in that area (Bajaj et al. 2011). It has also been reported in groundwater due to weathering, leaching and oxidation of soluble selenium compounds (Virk 2018). Though, selenium has an essential role in human health but at the same time toxicity had been reported at higher concentrations (Dhillon and Dhillon 2003) in Punjab. There was no reason behind occurrence of As contamination in Punjab (Virk 2019). The aluminum was reported naturally as oxides, hydroxides and silicates or as a complex with organic matter. The samples with negligible concentrations of Zn, Cd, Co, Pb, Mo and Ni in both the seasons were observed in the presentstudy.3.3 Statistical analysis One-way ANOVA is used for assessment of variation during the pre-and post-monsoon season in physico-chemical parameters as well as for the heavy metal content. It is necessary to evaluate the level of difference over change in sampling sites for both the seasons (pre-and post-monsoon) and their interaction. The heavy metals which were below detection limits were not used for statistical analysis. Results had shown statistically significant difference for most of the parameters at p < 0.05 (Tables 2, 3 and 4) ). The significant difference was observed due to different sampling locations, leaching of agrochemicals, influence of river and tributary water, inputs via industrial and domestic waste in groundwater.

Groundwater quality for drinking purpose
The water samples have been evaluated for suitability for drinking purpose based on the guidelines given by WHO (2011). It was noticed that values of water quality parameter and number of samples exceeding desirable limits (DL) and maximum permissible limits (MPL) set by the WHO (2011) and their health implications are given in Table 5. The suitability of studied water samples for drinking was examined from its TDS, TH and ionic concentration which is explained under following headings:

Total dissolved solids (TDS)
Total dissolved solids are regarded as one of the major determinants for studying suitability of drinking water. As per WHO, TDS value higher than 1000 mg/L along with high pH value resulted in unpleasant taste. Groundwater samples were classified into six classes on the basis of TDS ) as: excellent, good, fair, poor and unacceptable. TDS value (mg/L) ranged from 272.33 to 378.33 with an average value of 316.09 and from 320.00 to 475.00 with an average value of 406.18 for the pre-and post-monsoon samples, respectively. Out of the total 45 samples collected during the each season, 15 samples were under excellent class while remaining 30 samples were under good class for the pre-monsoon while in the post-monsoon season, all samples were under good class (Table 6). Similar findings were observed by Krishan et al. 2014 andRao et al. 2017.

Health risk assessment
Health risk assessment characterizes the adverse health effects caused by exposure of pollutants on children and adults (Adimalla and Qian 2019). Drinking is the most prominent factor of ingestion of the heavy metals through contaminated groundwater (Chen et al. 2016). In the present study, Fe exceeded the desirable limits (DL) (WHO 2011). Fe concentration was more than DL in 23 water samples and 26 samples, collected, respectively, in the pre-monsoon and post-monsoon season. On the other hand, concentration of As and Se was found lower than DL during both seasons. The health risk assessment of these heavy metals was determined by calculating the hazard quotient (HQ). The calculations were done on the basis of standards given by USEPA (2004). The results of the health risk assessment are shown in Table 8. The HQ values for As during the pre-monsoon were ranged from 0 (B 3 ) to 6.12E-02 (O 5 ) for children and 0 (B 3 ) to 6.38E-02 (O 5 ) for adults while in the post-monsoon, it was ranged from 1.39E-03 (B 3 ) to 7.93E-2 (O 5 ) for children and from 1.45E-03 (B 3 ) to 8.26E-02 (O 5 ) for adults. On the other hand, results had shown that the HQ values for Fe during the post-monsoon ranged from 6.36E-03 (B 1 ) to 2.64E-02 (P 3 ) for children and from 6.63E-03 (B 1 ) to 2.75E-02 (P 3 ) for adults while it ranged from 9.53E-03 (B 2 ) to 3.10E-02 (P 3 ) for children and from 9.94E-03 (B 2 ) to 3.23E-02 (P 3 ) for adults during the pre-monsoon. While in case of Se, the HQ values ranged from 3.67E-03 (O 1 ) to 1.15E-02 (B 5 ) for children and from 3.83E-03 (O 1 ) to 1.20E-02 (B 5 ) for adults during the pre-monsoon while during the Table 2 One-way ANOVA and Tukey's multiple comparison tests for water quality parameters of samples collected from the proximity of tributaries of Beas River, India for pre-monsoon period (June, 2019) All values in mg/L except pH (unitless), EC(μS/cm), RSC (meq/L), Mg ratio (meq/L), SAR (meq/L), As (μg/L) and Se (μg/L) Similar alphabets in a row indicates no significant difference between the means at p < 0.05 NS, not significant at p < 0.05, *Indicates significant at p < 0.05  Victor et al. 2018). It was noticed that the HQ value for As, Fe and Se was found less than 1 in all the samples which showed that the contamination was below risk level in the study area for both the adults and children.

Groundwater quality for irrigation purpose
The quality of the groundwater for the irrigation purposes has been given by the different experts. However, the most important parameters to determine suitability of the water samples were EC, SAR, Na%, RSC, MR, CR, PI and KR (Tripathi et al. 2012). The classification of groundwater samples for their suitability for irrigation purpose is given in Table 9. In addition, the suitability of the water samples for irrigation purpose was determined by plotting data on various scientific diagrams which included USSL diagram, Wilcox diagram and Piper Diagram which are explained in the following discussion. Wilcox (1955) had given classification of the groundwater on the basis of salinity hazard (EC) to check its suitability for irrigation purposes. He had classified the groundwater samples into five classes. After analyzing, it was observed that 22 groundwater samples were under good class and 23 samples were in permissible class during the pre-monsoon while during the post-monsoon all the samples collected from the study area were found in permissible class. Thus, results have shown that salinity get increased after the monsoon season which might had been due to the leaching of salts from upper layer of the soil into the groundwater table (Vinod et al. 2009;Ahamed et al. 2015). Value of the EC had been reported to rise with increase in TDS of groundwater . Table 3 One way-ANOVA and Tukey's multiple comparison tests for water quality parameters of samples collected from the proximity of tributaries of Beas River, India for the post-monsoon period (October, 2019) All values in mg/L except pH (unitless), EC (μS/cm), RSC (meq/L), Mg ratio (meq/L), SAR (meq/L), As (μg/L) and Se (μg/L) Similar alphabets in a row indicates no significant difference between the means at p < 0.05 NS, not significant at p < 0.05 *indicates significant at p < 0.05

Alkalinity hazard (SAR)
High concentration of Na + relative to Ca 2+ and Mg 2+ had been reported to reduce water supply to crop by decreasing soil permeability and the soil texture. SAR values in study the area were observed to be less than 10 in both seasons. So, the water samples collected in the both the seasons were under excellent category. Table 9 represents that the groundwater in study area was suitable for the irrigation purpose in the both the pre-monsoon and the post-monsoon season.

Sodium percentage (Na%)
It was reported that sodium reduces the permeability of soil by reacting with the ions (Ca 2+ and Mg 2+ ) present in the soil (Selvakumar et al. 2014). Na + ion was absorbed by the clay particles and destroyed the structure of the soil (Singh et al. 2015). So, classification of the water on the basis of the sodium percentage is an important parameter to check its suitability for irrigation purposes. Wilcox (1955) had classified the water samples into five classes on the basis of Na% to evaluate their suitability for the irrigation purpose (Table 9). After analyzing samples, it was observed that all the water samples collected during both the seasons were under good class and suitable for irrigation use.

Magnesium ratio (MR)
Groundwater was classified into suitable (>50%) and unsuitable class (<50%) based on MR value (Khodapanah et al. 2009). Mostly Mg 2+ and Ca 2+ are present in equilibrium concentration with respect to each other. But if Mg 2+ content is higher in water then it was observed that soil becomes alkaline and adversely affected the crop yield (Nagaraju et al. 2014;Kumar et al. 2017). After analysis, it was noticed that all the 45 samples collected from study area during the premonsoon were found suitable but during the post-monsoon season, 42 samples fall in the suitable class and remaining 3 samples fall in the unsuitable class (Table 9). The high amounts of Mg 2+ ion at few places were due to the high dissolution of dolomite mineral after monsoon as reported in earlier study ).

Kelly's ratio (KR)
The Kelly's ratio more than 1 is an indication of alkali hazard (Karanth 1987;Bhat et al. 2018). It was observed that all the samples collected during both the seasons were in the suitable category (Table 9). So, according to classification based on KR, the water of study area was found to be suitable for irrigation. Further it was observed that presence of low sodium content was good for permeability and structure of the soil (Ramesh and Elango 2012;Patel and Vadodaria 2013).

Residual sodium carbonate (RSC)
The presence of RSC in water had resulted darkening of upper soil surface, burning of foliage which had resulted in reduce crop yield (Toumi et al. 2015). It was observed that out of total 45 samples, 29 samples were fit and 16 samples were marginal during the pre-monsoon while during the post-monsoon season, 23 samples were fit and 22 samples fell in marginal class (Table 9). Results have indicated that HCO 3 − content gets increased with respect to Ca 2+ and Mg 2+ after the rainy season which could be due to excessive leaching of HCO 3 − ions from upper layer to lower layers of soil. Nishanthiny et al. (2010) and Bhat et al. (2018) also had evinced similar results in their studies.

Corrosivity ratio (CR)
Water is considered safe if its CR value is less than 1 while the value more than 1 denotes corrosive nature of  Table 5 Guideline values (WHO 2011) for drinking water along with number of samples exceeding maximum permissible limit (MPL) and desirable limit (DL) and health implications All standard limits in mg/L except pH (unitless), As (μg/L) and Se (μg/L) water as per Tripathi et al. 2012. The CR value indicated non-corrosive nature of all the 45 water samples collected during both the seasons as represented in Table 9. As a result of CR value, groundwater of this area is advised to be safely uplifted and transported by pipes.

US salinity diagram (USSL 1954)
US salinity diagram (USSL) was used to access the suitability of water samples for irrigation use (USSL 1954) in this studied area. It is based on EC and SAR values of water samples. On the basis of this diagram, groundwater was classified into: C1, C2, C3 and C4 types on the basis of salinity hazard (EC) and S1, S2, S3 and S4 types on the basis of sodium hazard (SAR) (Lokhande and Mujawar 2016). The groundwater data of study area were plotted using Grapher software (version: 16.2.354) which illustrates in Fig. 2 that out of total 45 samples collected in the pre-monsoon, 48.9% samples are C2S1 type and remaining 51.1% samples are C3S1 type having medium salinity-low sodium hazards and high salinity-low sodium hazards, respectively, while during the post-monsoon, all the samples were C3S1 type having high salinity-low sodium hazards. So, leaching of salts is required to eliminate salinity hazard so that all crops can be produced efficiently (Bhandari and Joshi 2013;Adimalla et al. , 2020.

Wilcox diagram (Wilcox 1955)
The Wilcox (1955) had given a diagram to evaluate the suitability of water for irrigation purpose by classifying into various classes on the basis of values of EC and Na% . The samples had been classified into five classes: excellent to good, good to permissible, permissible to doubtful, doubtful to unsuitable and unsuitable. Diagrams software (version-6.61) was used to plot the Wilcox diagram which illustrates that out of the total 45 samples collected during the pre-monsoon season, 71.1% samples are in excellent to good class and remaining 28.9% samples are in good to permissible class whereas, during the post-monsoon season, only 17.8% samples fall in excellent to good class and the remaining 82.2% samples are in good to permissible class as reported in earlier study . The results are represented in Fig. 3 illustrating that groundwater of the study area was good to be used for irrigation purposes during both the seasons.

Geochemical classification of groundwater (Piper diagram)
Piper (1944) had given trilinear diagram for geochemical analysis of groundwater. It has two equilateral triangle fields (i.e., cations, anions) and one central diamond-shaped field (Oki and Akana 2016). The major ion data collected was plotted in a Piper diagram using Grapher software (version-16.4.432) which is represented in Fig. 4a, b. During both the seasons, cation triangle indicates, all the samples fall in no dominant type region but the anion triangle shows that HCO 3 − was dominant among anions and all the samples are in HCO 3 − type zone. The resultant diamond-shaped field has shown that all the samples were Ca-Mg-HCO 3 type tributaries. Thus, there was no variation in geochemical faces in water samples during both the seasons and Kaur et al. (2017) had also observed similar results. The Ca 2+ ion was present in the form of Kankar in this region while Mg 2+ and HCO 3 − were supposed to be present due to the influence of Beas River (Purushothaman et al. 2013). High content of Ca 2+ had been reported due to dissolution of agrochemicals during the monsoon . The main source of Mg 2+ ion had been documented due to presence of dolomite. Further, HCO 3 − ion had been found in excess due to its leaching, bacterial oxidation of organic matter through aerobic and anaerobic complex chemical reactions and dissolution of CO 3 2− -rich minerals ).

Conclusions
The interpretation of hydro-geochemical analysis of the study area in the proximity of tributaries of Beas River, Punjab has revealed that the quality of groundwater was   good. The bulk of samples had lower values than permissible limits prescribed by WHO and BIS. Considering TDS as a parameter, the samples were under suitable class of drinking and irrigation purposes in all the studied locations. The abundance order of major cations was Ca 2+ >Mg 2+ >Na + >K + and Ca 2+ > Na + >Mg 2+ >K + , respectively, in the pre-monsoon and the post-monsoon season while abundance of major anions was in the order of HCO 3 − > SO 4 2− > Cl − > CO 3 2− in both the seasons. Among heavy metals, Fe, As, Al and Se were found abundant in the samples having average concentration 0.30 mg/L, 0.19 μg/L, 0.02 mg/L and 0.78 μg/L, respectively, during the pre-monsoon and 0.33 mg/L, 0.25 μg/L, 0.03 mg/L and 0.84 μg/L, respectively, during the post-monsoon. Health risk assessment of Fe, As and Se show that their concentrations were below risk level in the study area. From this study, it was concluded that the groundwater of the study area was suitable for drinking and irrigation purposes. The continuous monitoring of water in this area is recommended to cope with the rising contamination of groundwater and adverse effects on human health.
Acknowledgements Authors would like to thank Principal Khalsa College, Amritsar for providing necessary resources. Authors would also like to thank Dr. Kultar Singh from post-graduation department of chemistry for analytical determinations.
Author contributions JS performed the experiments, analyzed the data, prepared figures and/or tables and approved the final draft. SKS assisted in design the experiment, statistically analyzed the data, prepared figures and/or tables, prepared manuscript and approved the final draft. KS critically reviewed the drafts of the paper and approved the final draft. DS conceived and designed the experiments, authored or reviewed drafts of the paper.

Funding
The author(s) received no specific funding for this work.
Availability of data and material The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflict of interest Author's have no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.