Assessment of groundwater quality and its controlling processes in Bemetara District of Chhattisgarh State, India

Groundwater withdrawal at very fast rate poses threat on existing groundwater resources in different parts of the world. This reduction in groundwater levels significantly disturbed the natural aquifer flow rate and thereby different hydrogeochemical processes, which may further impair the groundwater quality. The groundwater quality in rural area of Chhattisgarh State is degraded, and the problem of saline water poses health risk to people. In this research investigation, suitability of groundwater of Bemetara District, Chhattisgarh, India, has been evaluated for drinking purpose through water quality index (WQI) method and principal component analysis (PCA). Total 116 groundwater samples were collected during the pre-monsoon (June 2019) and post-monsoon season (December 2019) and analyzed for physicochemical parameters. Total dissolved solids ranged from 250 to 10,440 mg/L and 289 to 3583 mg/L during pre-monsoon and post-monsoon, respectively, and 55% of the total samples exceeded acceptable BIS limit in pre-monsoon, while about 66% samples exceeded in post-monsoon season. SO42− concentrations varied from 3 to 5734 mg/L during pre-monsoon and 4.5 to 2002 mg/L during post-monsoon, respectively. Total 28% samples in pre-monsoon and 18% samples in post-monsoon season exceeded the maximum permissible BIS limit (400 mg/L) of SO42− ion in the study area. On the basis of WQI, the quality of groundwater varies from “Excellent water” to “Good water” category. The groundwater of northeastern part of the district is not suitable for drinking, and therefore, it is recommended to treat this groundwater before human consumption with special reference to SO42− contamination. PCA inferred that four components are sufficient to explain the variance in chemistry of groundwater that is mainly governed by dissolution of gypsum mineral, other rock–water interaction and anthropogenic activities. Further, water quality was improved in the direction of groundwater flow in the study area, establishing a direct relationship between groundwater flow and water quality of the Bemetara District. This study provides very useful database to design sustainable groundwater management plan for the district.


Introduction
Untreated waste discharge leads to degradation of the surface water quality, and therefore, water supply for different sectors like agricultural, industrial and domestic needs has been fulfilled by groundwater resources. But the rate of usage of groundwater has resulted into declining groundwater levels, which is very critical to available resources, and this decline level reaches up to 80 m in zone of depression (Chen et al. 2005). This reduction in groundwater levels may significantly altered the groundwater flow conditions (Wang et al. 2008;Zhang et al. 1997Zhang et al. , 2000Fan 1998;Xia et al. 2004). Therefore, different studies have been taken to investigate the aquifer flow conditions in respect to the sustainable utilization of groundwater resources (USGS 1999; CGWB Report 2019). This decline in groundwater level leads to change in the aquifer hydrogeochemical processes (CGWB Report 2019).
Groundwater quality depends on important process like atmospheric precipitation, runoff and inland surface water. Moreover, the groundwater quality is degraded due to disposal of industrial waste and mining activities (Rodell et al. 2009;Chopra and Gopal 2014;Malyan et al. 2019). Further, groundwater quality in a region is influenced by physical and chemical parameters that are strongly affected by natural processes such as water chemistry in the recharged area, water intermixing, groundwater recharge, aquifer discharge and recharge, and water flow path (Singh et al. 2017;Azhdarpoor et al. 2019;Soleimani et al. 2018). The spatial distribution and zoning of NO 3 − and F − concentration and their health risk assessment in drinking groundwater of Shiraz metropolitan area in the southwest of Iran and Behbahan City were studied by application of Monte Carlo simulation, sensitivity analysis and geographic information system (Badeenezhad et al. , 2021. The most convenient method to describe the quality of drinking water resources is the Water Quality Index (WQI). The first attempt to develop a WQI was made in 1948, when scientific community found a correlation between pollution load and certain group of organisms (fish, plant and benthic community) (Alves et al. 2014). Later, Horton developed WQI technique in 1965 (Horton 1965). After that, National Sanitation Foundation (NSF) of United States developed WQI in 1970 which is widely used (Brown et al. 1970). Therefore, WQI is not new tool, but has been extensively used across the world to determine the water quality (Abbasi and Abbasi 2012).
Principal component analysis (PCA) is a useful tool to investigate the chemical relationship between different water quality parameter and thereby predicting the dominating parameters (Sharma and Jain 2006). PCA is a multivariate statistical technique that has been widely used to reduce the dimensionality of large data (Vega et al. 1998;Duan et al. 2016;Zhang et al. 2016). The goal of PCA is to describe the majority of the data sets in a few principal components (PCs), and these PCs are the linear combination of observed data with maximum variations with minimum loss of the actual information (Baghanam et al. 2020). PCA has been applied to explain water quality variable in several studies (Vermonden et al. 2009;Daou et al. 2016;Zhang et al. 2016;Abdelaziz et al. 2020;Chai et al. 2021). Huang et al. (2007) used PCA to explain the storm water quality data for pattern recognition and identification of pollution sources from different urban surface-type catchments. Sharma and Jain (2006) evaluated groundwater quality of Jodhpur District, Rajasthan (India), using multivariate technique and concluded that F − , NO 3 − , pH and K + have significant influence on the quality of the aquifer. Marin Celestino et al. (2019) used PCA to investigate hydrogeochemical processes in a wastewater-irrigated region (central Mexico) and reported that groundwater chemistry was dominated by three processes: salinization, mineralization and groundwater contamination. Chai et al. (2021) investigated the source assessment of pollution in the Fen River for different seasons. Giakwad et al. (2020) evaluated the groundwater quality of western coast of Maharashtra, India, with the use of PCA technique. Abdelaziz et al. (2020) studied groundwater quality index based on PCA in Wadi El Natrun, Egypt, and used PCA to reduce the complexity of the data and identified the group of parameters (Na + , SO 4 2− , Cl − , strontium, Ca 2+ and molybdenum) that control the groundwater quality.
The groundwater quality in rural area of the Chhattisgarh State is degraded, and the problem of saline water (EC ~ 2000-4500 µS/cm) poses a threat to the people. There is a serious problem of saline water in 113 villages of district Bemetara of Chhattisgarh State, and 1200 to 2600 mg/L of TDS was observed in the groundwater of problematic villages of district Bemetara, which may cause health problem, viz. digestion, high blood pressure, heart attack and kidney problems (CGWB Report 2015). In some parts, the value of SO 4 2− is observed up to 800 ppm which causes gastrointestinal disorders among the inhabitants of the area (Mukherjee and Gupta 2010). According to the groundwater quality assessment by Central Ground Water Board, a central agency of government of India, the SO 4 2− concentration found up to 763 ppm in Bemetara village of Bemetara District in (CGWB Report 2015. For the alternate sources, the residents of these villages are using contaminated water from ponds, rivers and drains in the area. In a study of health risk evaluation of uranium in groundwater of Bemetara District, uranium levels in water samples range from 1.15 to 83.5 µg/L and 0.68 to 96.08 µg/L during premonsoon and post-monsoon, respectively, with few samples exceeding the safe limit of 30 µg/L prescribed by WHO (2011), and there is no harmful effect by radiological risk, but chemical risk can affect human health (Sahu et al. 2020). Dahariya et al. (2020) studied the contamination, sources and environmental hazard of groundwater in Bemetara District of Chhattisgarh, India, and reported WQI (406 ± 82) values which clearly demarcate groundwater unsuitability for drinking purposes.
In view of the above scenario of degraded groundwater quality in the district Bemetara of Chhattisgarh State, the aim of the present investigation is (1) to monitor the quality of groundwater for drinking purpose using WQI; (2) to monitor spatial and seasonal variation of important water quality parameters using the application of GIS software; and (3) to identify several factors responsible for degradation of quality of groundwater using principal component analysis.

Study area
Bemetara District is newly formed district of Chhattisgarh State, India, and covering area of 2854.81 km 2 (Fig. 1). It lies in between 21° 22′ and 22° 03′ North latitude and 81° 07′ and 81° 55′ East longitude. Bemetara District has huge quantity of mineral deposits, namely sandstone, limestone (low grade), dolomite and quartzite. Dolomite and limestone mineral were found high in whole district. The study area has a dry and wet tropical climate. The temperature varies from 10 to 48 °C, where the maximum temperature is reached in the month of May and June and minimum temperature fall in January. Bemetara District has flat topography, and totally six rivers flow in the direction of slope of district (north to east), namely Shivnath, Kharun, Surahi, Haff, Sakari and Phonk rivers. Bemetara District geologically comes under Meso-to Neoproterozoic rock sequence of Chhattisgarh supergroup. This Chhattisgarh Fig. 1 Map showing the location of sampling sites in the study area supergroup is divided further into different groups; our study area falls under Raipur group which comprises four types of geological formation (Maniyari formation, Hirri formation, Tarenga formation and Chandi formation). The reddish brown and purple non-calcareous shale containing gypsum form the characteristics of Maniyari formation, under most of the district area occurred (District Survey Report Bemetara Chhattisgarh 2016).

Hydrogeology of the study area
The hydrogeological formation of study area mainly consists of arenaceous-argillaceous-calcareous rocks and is enriched by limestone/dolomite and calcareous shale. The groundwater in these formations occurs under water table, semi-confined and confined conditions. The weathered, cavernous and fractured part of the formation constitutes the aquifer in the area and has great potential in regards to groundwater yield and thereby development groundwater in the district. Geology and hydrogeology of the study area shown in Figs. 2 and 3, respectively. The gypsum karsts occurring in the Maniyari formation of this province are more productive. Though gypsum is more soluble than calcite, their alternative assemblage with thinly laminated shale provides special condition that favors dissolution of gypsum laminae causing roof collapses to create larger openings. However, all the formations in the district are productive (CGWB Report 2015).

Water sampling and analytical techniques
In total, 116 samples were collected from groundwater sources of 51 locations of the district Bemetara, viz. open wells, dug wells, borewells and handpumps, during premonsoon (June 2019) and post-monsoon (Dec. 2019), which are extensively being used for drinking water purpose and analyzed for physicochemical parameters using standard methods (APHA 2005). Groundwater samples were taken from one open well from each of 51 locations, but borewell samples were also collected from seven locations along with open well. Before collecting samples, handpumps/borewells were pumped for 5 min to get represented water sample and the sampling bottle was rinsed with the same water. In situ parameters were analyzed on site like pH and electrical conductivity using Hach, USA make HQ40d portable handheld multimeter. Other parameters like major cation and anion were analyzed using Metrohm ion chromatograph. The ionic balance error (IBE) test was performed (Eq. 1) (Freeze and Cherry 1979) and was below 5% for all samples, which support the data accuracy and reliability of the analysis. Total alkalinity and HCO 3 − were determined by autotitrator of SI analytical instrument, a Xylem brand.

Water quality index (WQI)
WQI is an important method that is used for the evaluation of quality of water for drinking purposes (Subba Rao 1997;Avvannavar and Shrihari 2008;Mishra and Patel 2001;Badeenezhad et al. 2020). BIS (2012) and WHO (2011) set standard limit to different water quality parameters for drinking purposes, and some of these are incorporated into calculation of WQI. Here, we took ten parameters (Table 1) for calculating WQI and each parameter is assigned a weight (w i ) depending upon the importance of overall quality of water.
The assigned weight ranges from 1 to 5, where 5 is the highest and 1 is the lowest weight (Srinivasamoorthy et al. 2008;Vasanthavigar et al. 2010). In the next step, the relative weight (W i ) is calculated by Eq. 2 as follows: (1) (2)  where W i is the relative weight, w i is the assigned weight of each parameter and n is the total number of parameters. Chemical parameters and their calculated relative weight (W i ) are given in Table 1.
The next quality rating value (q i ) is calculated by Eq. 3: where q i is the quality rating, C i is the concentration of parameter in water sample (ppm) and S i is the BIS value for each parameter in ppm. And finally WQI is calculated by Eqs. (4) to (5): where SI i = subindex of ith parameter; q i = rating based on the concentration of ith parameter; n = number of parameters.
Five WQI classes and water types can be categorized accordingly as given in Table 2.

Data processing
An univariate and multivariate normal distribution are required for the best outcomes of the statistical multivariate methods like PCA (Zhou et al. 2007;Oppong and Agebedra 2016;Marin et al. 2018). Shapiro-Wilk's (Shapiro and Wilk 1965) and Royston's tests (Royston 1983) were used to verify the univariate and multivariate normality conditions, respectively. The Spearman's rank correlation was used in this multivariate analysis because water quality data were non-normal distribution, and this correlation method is best suitable for reducing deviation of variables (Marin Celestino et al. 2019). Based on Royston's test, the dataset has a non-normal distribution. To achieve a normal-like distribution, the original set of variables was transformed using a logarithmic transformation (natural logarithm). To approach the best conditions of the multivariate analysis, feature scaling on the database was done using standardization (or Z-score normalization). Standardization minimizes the variance in variables and protects dissimilarity metrics such as the Euclidean distance from being severely influenced (Davis and Sampson 1986). Each variable was normalized to its Z score, which was determined using Eq. (6): where Z i is the standardized Z score, X i is each variable's value, and mean and S are each variable's mean value and standard deviation, respectively. The Kolmogorov-Smirnov (K-S) test was used to assess how well the transformed variables were adjusted to the normal log distribution (Rizvi et al 2015;Muangthong and Shrestha 2015;Marin Celestino et al. 2019;Castillo et al. 2021). The precision and acceptability of the data for PCA were assessed by using the Kaiser-Meyer-Olkin (KMO) and Bartlett's sphericity tests. KMO is a metric for determining the sampling's adequacy by identifying the proportion of shared variation that might be linked by unknown factors (Marin Celestino et al. 2019;Castillo et al. 2021). The KMO and Bartlett's sphericity test results are described in the results and discussion section.

Principal component analysis (PCA)
In PCA, the standardized and normalized data as discussed in earlier section are used. Factors are produced through an eigenvalue analysis of the correlation matrix. These factors are vector which shows interorthogonality within a multidimensional space defined by a number of variables in the analysis (Sharma and Jain 2006). Unlike the original variable, the factors are uncorrelated with each other. They are described by means of their correlation with (or "loading" on) the original variables and ranked in order of the amount of the total variance they explain. A loading close to one indicates a strong relationship between factor and the variable, whereas a zero loading indicates no relationship (Davis and Sampson 1986;Sharma and Jain 2006). Among the whole factors, first few factors explain the bulk of total variance and the remaining factors are not used in the analysis. These retained factors were then rotated using Varimax method. Varimax rotation tried to attain a simple structure, whereby factor loading is close to one or zero. This helps in the interpretation of the factors that either does or does not include a particular variable. Principal component and factor analysis are applied on the groundwater of the district Bemetara, and the whole analysis has been carried out using software RStudio Vs 1.4.1106.

Hydro-chemical characteristics of groundwater of the study area
In total, 116 samples were collected from the study area during the year 2019-2020 from the abstraction sources in collaboration of Water Resources Department (WRD), Govt. of Chhattisgarh, Raipur. The minimum, maximum and average values of all the chemical parameters for both the seasons are given in Table 3.
The pH of all collected samples is alkaline in nature and varies from 6.2 to 7.7 and 6.1 to 7.1 during pre-monsoon season (PRS) and post-monsoon season (PMS), respectively. Almost all samples are found within BIS (2012) and WHO (1996) standard limits of 6.5-8.5.
In Bemetara District, total dissolved solid (TDS) in groundwater ranges from 250 to 10,440 mg/L and 289 to 3583 mg/L during PRS and PMS, respectively. About 55% samples exceeded the acceptable limit (ALT) but lies below the prescribed maximum permissible limit (MPL) of 2000 mg/L in PRS and about 66% samples were found above ALT but within 2000 mg/L in PMS. The maximum value of TDS 10440 mg/L was observed in the groundwater of the village Kunra of the block Nawagarh. The spatial distribution map of TDS for PRS and PMS is shown in Fig. 4. CO 3 − , HCO 3 − and hydroxides are mainly responsible for alkalinity in water system. The alkalinity value varied from 83 to 280 mg/L and 52 to 415 mg/L during PRS and PMS, respectively. About 50% samples exceeded ALT but under the prescribed MPL of 600 mg/L in PRS, and about 59% samples were found above ALT but are under 600 mg/L in PMS. None of the sample crossed the prescribed limit of 600 mg/L during PRS and PMS, respectively (Fig. 5).
In groundwater of the study area, the values of Ca 2+ range from 26 to 569 mg/L and 26 to 648 mg/L during PRS and PMS, respectively. About 24% samples crossed the prescribed MPL of 200 mg/L during the PRS and 17% samples crossed the prescribed MPL during the PMS. The less number (7%) of samples crosses the MPL in the PMS because of the dilution effect as compared to PRS. The value of Mg 2+ varies from 11 to 488 mg/L during PRS and 12 to 259 mg/L during PMS. Only 12% samples exceed prescribed MPL 100 mg/L during PRS and 7% exceed above 100 mg/L during PMS.
The value of Cl − ranges from 10 to 1080 mg/L and 12 to 652 mg/L during PRS and PMS, respectively. More than 90% samples lie within the prescribed ALT of 250 mg/L during both the seasons. NO 3 − content in the study area ranges from 0 to 194 mg/L and 0 to 569 mg/L during PRS and PMS, respectively. About 93% of the samples of the study area fall within MPL of 45 mg/L and 7% of samples even crossed the MPL during PRS and about 67% of the samples of the study area fall within the prescribed MPL of 45 mg/L and 33% of samples even crossed the prescribed MPL during PMS. F − in the groundwater of Bemetara District ranges from 0.06 to 2.4 mg/L and 0 to 1.04 mg/L during PRS and PMS, respectively. Almost all samples lie under the prescribed ALT of 1.0 mg/L during both the seasons.
SO 4 2− is generally found as soluble salts of Ca 2+ , Mg 2+ and Na + in the groundwater. During the PRS, the concentration of SO 4 2− ranges from 3 to 5734 mg/L and 4.5 to 2002 mg/L during PMS. BIS sets 200 mg/L as the ALT and 400 mg/L as MPL for SO 4 2− in drinking water. In Bemetara District about, 55% of the samples are below the prescribed ALT of 200 mg/L and 28% samples exceeded the value of MPL of 400 mg/L during the PRS, while 67% samples lie below the ALT and 19% sample exceeded prescribed MPL of 400 mg/L during PMS. The spatial distribution of SO 4 2− in the study area is shown in Fig. 7 for PRS and PMS. The northeastern area of the study area is highly contaminated as shown in the spatial distribution maps for both the seasons, and the maximum value of SO 4 2− 5735 mg/L was observed in the groundwater of the village Kunra of the block Nawagarh.

Depthwise variation of groundwater quality of the study area
Depthwise distribution of groundwater samples has been arranged in three categories, i.e., shallow aquifer (0-20 m), medium aquifer (21-40 m) and deep aquifer (> 40 m). In PRS, about 31% samples from the shallow aquifer (SA), no sample from medium and 2% sample from deep aquifer lie under the ALT of TDS, i.e., 500 mg/L (BIS 2012), while about 49% samples from the SA, 2% samples each from medium and deep aquifer exceeded limit of TDS, i.e., 2000 ppm (BIS 2012), and only 14% samples from SA exceeded the MPL of TDS. In case of PMS, about 26% samples from the SA, no sample from the medium and deep aquifer fall within the ALT, while about 64% samples from SA, 2% samples from deep aquifer exceeded the ALT but under the MPL of TDS and 9% samples from SA exceeded the MPL of TDS (Table 4).
During PRS, about 48% samples from the SA, 2% samples from medium and no sample from deep aquifer fall within ALT of alkalinity, i.e., 200 mg/L (BIS 2012), while about 46% samples from the SA, 0% samples from medium aquifer and 3% samples from deep aquifer lie under the limit of MPL of alkalinity, i.e., 600 mg/L (BIS 2012) and no samples from SA, medium aquifer and deep aquifer exceeded the MPL of alkalinity. In case of PMS, about 41% samples from the SA, no sample each from medium and deep aquifer fall within the ALT, while only 57% samples from SA and 2% samples from deep aquifer lie above ALT but under the MPL of alkalinity and none of samples from SA, medium aquifer and deep aquifer exceeded the MPL of alkalinity (Table 5).
In PRS, about 12% samples from the SA, 0% samples each from medium and deep aquifer fall within the ALT of hardness, i.e., 200 mg/L (BIS 2012), while about 48% samples from the SA, 0% samples from medium aquifer and 2% samples from deep aquifer lie below the MPL of hardness, i.e., 600 mg/L (BIS 2012), and only 38% samples from SA exceeded the MPL of hardness. In case of PMS, about 5% samples from the SA, no sample from the medium and deep aquifer fall within the ALT, while about 66% samples from SA, 2% samples from deep aquifer crossed the ALT but are under MPL of hardness and 27% samples from SA exceeded the MPL of hardness (Table 6).
In PRS, about 53% samples from the SA, 0% samples from medium and 2% samples from deep aquifer fall within the ALT of SO 4 2− , i.e., 200 mg/L (BIS 2012), while about 17% samples from the SA, 0% samples from each medium aquifer and deep aquifer lie under MPL of SO 4 2− , i.e., 400 mg/L (BIS 2012), and 26% samples from SA, 2% samples from medium aquifer and 0% samples from deep aquifer exceeded the MPL of SO 4 2−. In PMS, about 66% samples from the SA, no sample from the medium and 2% samples from deep aquifer fall within the ALT, while only 14% samples from SA lie in between ALT and MPL of SO 4 2− and 19% samples from SA exceeded the MPL of SO 4 2− (Table 7). From the above discussion, it may be concluded that almost all collected groundwater samples belong to shallow aquifer and the significant amount of collected contaminated with higher SO 4 2− contamination. Further, water quality at different depths at the same site of few locations has been studied and it was observed that higher concentrations of different water quality parameters were generally observed at higher depths below the ground, which is due to more residence time of groundwater in the deeper aquifer (Fig. 8).

Correlation among the chemical parameters
Correlation among the physicochemical parameters for both the seasons has been evaluated using correlation matrix and given in i.e., dissolution of dolomite or gypsum mineral. Further, SO 4 2− plotted against the Ca 2+ , Mg 2+ , Na + and K + (Fig. 10) and the best relationship was observed between Ca 2+ and SO 4 2− (maximum r 2 ), further supporting the fact that the source of SO 4 2− in the groundwater of the study area may be CaSO 4 , i.e., gypsum, which is present in Maniyari shale formation of the region.

Principle component analysis
PCA is a multivariate method that reduces the dimensions of a dataset that contain different variables which are interrelated to each other. PCA was used to study the origin of major salts in the groundwater of Bemetara. Before applying PCA, two tests were performed to analyze the statistical interrelation among the parameters. The KMO test value was 0.5, and Bartlett's sphericity value was 0.00 that confirmed the data are appropriate and suitable for PCA (Zhang et al. 2020). Eigenvalues and cumulative contribution of all principal components for both the seasons are given in Table 10, and varimax rotated component loadings are given in Table 11. Different components incorporated in the explication of sources for both the seasons are presented in Table 12.
In pre-monsoon season, only four principal components (PCs) were taken into account with eigenvalue more than   (Table 10). The first principal component (PC1) was mainly characterized by SO 4 2− , EC, Mg 2+ , Ca 2+ , Cl − and F − in groundwater that may be attributed to dissolution of gypsum mineral and weathering of Cl − bearing rocks (Sharma and Kumar 2020;Mullaney et al. 2009). In addition, the Fmay be due to rock-water interactions in the aquifer, and therefore, this factor (F1) may be considered to represent the local geogenic process. Also, PC1 shows 50.49% of the total variation in the data. The second principal component (PC2) accounts17.80% of the total variation in the dataset and had loading on NO 3 − , HCO 3 − and alkalinity. NO 3 − in groundwater is mainly contributed from extensive usage of chemical fertilizers (Zhang et al. 2020), and therefore, this factor (F2) was considered to represent the anthropogenic source. The third principal component (PC3) characterized by K + , Mg 2+ , Cl − and NO 3 − accounts for 10.72% of the total variation in the data. K + , Mg 2+ and Cl − were resulted from the interaction of groundwater and rock materials. The NO 3 − in groundwater was associated with chemical fertilizers application in farming and could consider this factor (F3) mixed type source, i.e., geogenic and anthropogenic. The fourth principal component (PC4) had been loaded with pH, Na + and NO 3 − which accounts for 8.08% of the total variation in the dataset. Na + may be observed in livestock and domestic waste (Zhang et al. 2020) and could enter subsurface system where there are inappropriate management activities of waste. Therefore, this factor (F4) was called as nonpoint source of pollution mainly from agriculture.
In the post-monsoon season, four PCs can explain 87.61% of total variation with individual contributions of 54.30%, 16.43%, 11.13% and 5.76%, respectively (Table 10). EC, hardness, Na + , Ca 2+ , Mg 2+ and SO 4 2− had high loading in PC1. These ions are related to natural geogenic processes as same explained in factor (F1) in pre-monsoon season. In PC2 loading element were pH, HCO 3 − , SO 4 2− and NO 3 − , accounting for16.43% of the total variance. NO 3 − and SO 4 2− could be resulted from different agricultural activities  (Zhang et al. 2020) and consider similar as to source (F4) in the pre-monsoon season, i.e., nonpoint source of pollution. HCO 3 − , Cl − , NO 3 − and F − were loading elements that contribute 11.13% of the total variance. NO 3 − related to anthropogenic source, whereas F − and Cl − are mainly attributed to geogenic source. Therefore, this factor (F3) is considered as mixed type source as explained in the pre-monsoon season, i.e., geogenic and anthropogenic. The fourth principal component (PC4) was mainly characterized by Na + and Cl − , explaining 5.76% of the total variance. As explained earlier in the pre-monsoon data, this factor (F1) can be regarded as local geogenic processes and salinity.

Variation of water quality along groundwater flow
The investigation of aquifer water flow conditions is very crucial for sustainable water resource management of any region. The movement of groundwater flow interacts with aquifer rock material and carries chemical species in dissolved form (Rakhmatullaev et al. 2010(Rakhmatullaev et al. , 2012Huneau et al. 2011).
Groundwater level observations from the mean sea level (MSL) were taken from the all selected sampling points in the study area for both the seasons. For the identification of the direction of water flow, spatial contour maps of water levels have been prepared for both the seasons (Fig. 11). Groundwater moves from northeast direction to northwest direction and from northeast to southeast direction during both the seasons of the study area. Further, River Kharun and River Sheonath are flowing along the southeastern boundary, further supporting the direction of groundwater flow direction (Fig. 3). However, from the spatial distribution map of SO 4 2− (Fig. 7), it is evident that the concentration increases from northwest to northeast and southwest to northeast part of the study area. Therefore, it may be inferred from the above discussion that the SO 4 2− decreases with the direction of flow of groundwater (i.e., NE to NW). Also, among all four blocks of the district (Nawagarh, Bemetara, Berla and Saja), two blocks, i.e., Nawagarh and Bemetara, have high SO 4 2− contamination in the study area. High concentration in these blocks may be attributed to existing Maniyari shell formation comprising reddish brown and purple non-calcareous shale with gypsum interbands.
Water quality improves from northeast to northwest part of the district during both the seasons (Fig. 9). As we compared the trend of groundwater flow and water quality (WQI), it was observed that water quality was improved along the flow of groundwater. Therefore, a direct relationship was found     Geogenic factor (F4) Na + and Cl − Fig. 11 Groundwater level (msl) in the study area during pre-and post-monsoon (2019) between groundwater flow and water quality of the Bemetara District.

Conclusion
For any scheme of drinking water supply in an area, it is mandatory to have the status of water quality of the water resources being used for drinking water supply. Therefore, the evaluation of water quality becomes essential for development of water resource strategies for sustainable water use and to provide database for future planning. Spatial distribution maps were prepared to identify degraded water quality zones, possible sources of pollution and specific parameters not conforming to drinking water quality standards. BIS for drinking water have been violated for physicochemical parameters, viz. TDS, total hardness, Ca 2+ , Mg 2+ , SO 4 2− and NO 3 − . The groundwater quality changes from region to region along the depth of water table, and generally higher concentrations were observed in deeper aquifers. The source of SO 4 2− in the groundwater of the study area may be attributed to dissolution of gypsum as evident from relationship between Ca 2+ and SO 4 2− (r 2 > 0.8). WQI at different locations has been computed to check the quality of groundwater for drinking purpose and 60% of the samples fall from "Excellent water" to "Good water" category in PRS and about 62% of samples in PMS. About 9% samples were found to be unsuitable for drinking purpose. Spatially, it is revealed that the drinking water sources existing in the northeastern area of the district were found to be contaminated with high SO 4 2− concentration which is not fit for direct public consumption. Multivariate analysis gives insight to different possible potential sources contributing to groundwater pollution in the area and inferred that four components are sufficient to explain the variance in groundwater chemistry mainly controlled by dissolution of gypsum mineral, other rock-water interaction and anthropogenic activities. Water quality was improved in the direction of groundwater flow in the study area, establishing a direct relationship between groundwater flow and water quality of the Bemetara District. This study provides very useful database for policymaker and state government to design sustainable groundwater management plan for the district.