Water Quality Indices and GIS-based evaluation of a decadal groundwater quality

Abstract The water quality data of fixed 44 wells were monitored during post- and pre-monsoon season 2005–2013 of Kanchipuram district of state Tamil Nadu, India. Water Quality Indices (WQIs) methodology was applied to categorize groundwater on the basis of rating. In total, 13 groundwater quality parameters: calcium (Ca2+), magnesium (Mg2+), sodium (Na+), potassium (K+), bicarbonate (HCO3 −), nitrate (NO3 −), chloride (Cl−), fluoride (F−), sulphate (SO4 2−), total hardness (TH), total dissolved solids (TDS), hydrogen ion concentration (pH) and electrical conductivity (EC) were taken for statistical analysis. Three different WQI models were computed using a decadal physicochemical water quality data mainly for drinking (WQI1 and WQI2) and irrigation (WQI3) need. The WQIs variability thematic maps were generated in the Geographical Information System environment through ArcGIS software. These spatial distribution maps of WQI provide a detailed overview. The work outlined that the groundwater of the area is impaired by the anthropogenic activities and needs proper management plan to control further contamination and pollution of the groundwater of the area. Further, the groundwater of few monitoring wells needs a proper treatment before its consumptive uses. The alternative water supply should be provided to local inhabitants for daily uses.

Water Quality Indices (WQIs) are widely used tool in evaluation of water quality; it summarizes a number of water quality parameters into a meaningful single numerical value (Singh et al. 2013aStambuk-Giljanovik, 2003). It is simple and easy to understand for decision makers about quality and possible uses of any water body (Jacintha et al., 2016). WQIs serve the understanding of water quality issues by integrating large complex data and generating a score that provide to prioritize the management plan and fund. The water resources are impaired by human induced activities and geogenic activities (Amin, Fazal, Mujtaba, & Singh, 2013;Singh, Singh, Kumar, Gupta, & Mukherjee, 2009). WQI reflects the composite influence of different water quality parameters and is being used to provide comparisons of water quality status for different locations and at different time and acts as a water quality status prediction tool (Singh et al., 2013aTiwari et al., 2014). Still there is no globally accepted composite WQI in literature. Few commonly used WQIs namely National Sanitation Foundation Water Quality Index Additive Model (NSFAWQI)(assessment and communication), National Sanitation Foundation Water Quality Index Multiplicative Model (NSFMWQI) (insensitive to low sub index values), Unweighted Multiplicative Water Quality Index (UMWQI) (need to minimize overall index and sub index values) and Canadian Council of Ministers of the Environment Water Quality Index (CCME-WQI; to manage large complex data and to evaluate success and failure of management strategies) and many more are referred by many researchers. WQIs are relying on normalizing, or standardizing parameters values by parameter according to expected concentrations, and divulge as "good" vs. "bad" concentration. The UMWQI model parameters are weighted according to their perceived importance to overall water quality and the index is calculated as the weighted average of all observations of interest (Liou, Lo, & Wang, 2004;Pesce & Wunderlin, 2000;Sargaonkar & Deshpande, 2003;Stambuk-Giljanovik, 2003;Tsegaye et al., 2006). Scientists are trying to develop a unique WQI by minimizing the overall index and sub index values. WQIs provide ways to introduce and determine that particular standard value and helps in prediction of harmful effects.
WQIs have been mostly used in evaluation of potable water; very limited study has been carried out to evaluate the water for agriculture purpose . EC and Na + play a vital role in evaluation of water status for irrigation. Higher EC in water creates a saline soil and deteriorates the soil health and reduces the crop yield. The salts apart from affecting the growth of plants also affect the soil structure, permeability and aeration which indirectly affect plant growth. Higher salt content in irrigation water causes an increase in soil solution osmotic pressure (Shaki & Adeloye, 2006;Thorne & Peterson, 1954). The saline soils are reclaimed by the use of gypsum.
Geographical information system (GIS) used for variety of application from micro-level planning to macro-level planning, implementation and regular monitoring (Srivastava et al., 2013). The GIS is rapidly growing field of research with an integration of different subjects from atmosphere to ocean, land, sea, water resources, health, defence, agriculture, demography, forestry to economics and many more. GIS is defined as a technique to capture, store, retrieve, analyse and predict (Gajbhiye et al., 2016). GIS is also used as database system to store and prepare 2D, 3D or 4D water quality status maps according to concentration values of different chemical constituents. Different researchers used the GIS in the field of groundwater quality mapping (Singh et al., 2013a(Singh et al., , 2013b. The integration of WQIs with GIS provides the detailed, quick and reliable information for decision makes to adopt or implement strategies related to water pollution and scarcity (Singh et al., 2013b). The integration of WQIs data with other physical or biological driver's database in GIS helps in improving the already polluted water bodies and suggests the decision maker to take necessary action in time. The water suitability of sampled well is evaluated based on recommended permissible limits for parameters as recommended by WHO (2006) for drinking use and food and agriculture organization (FAO) for agriculture purpose.
The main objective of the work was to evaluate the decadal water quality data for potable and irrigation uses. In this work, WQI 1 and Nemerow's Pollution Index (NPI) (WQI 2 ) index is used for appraisal of groundwater for domestic purposes. Whereas ratio of dissolved sodium as %Na (WQI 3 ), is used to evaluate groundwater for irrigation purpose.

Study area
Study area Kanchipuram district of Tamil Nadu state (11° 00′-12° 00′ North latitudes and 77° 28′-78° 50′ East longitudes; Figure 1(a)) with average elevation as 83.2 m above mean sea level. The land around study area is flat and slopes towards the south and east. Study area is bound by Bay of Bengal in the East. The 47% population is engaged in agriculture activities and paddy is dominant crop. The soil in the region is mostly clay, with some loam, clay and sand. This district is mainly depending on the seasonal rains, the distress conditions prevail in the event of the failure of rains. Northeast and Southwest monsoon are the major donors with 54 and 36% contribution each to the total annual rainfall. During normal monsoon, the district receives average annual rainfall of 1213.3 mm (Figure 1(b)).

Materials and methods
The data related to groundwater quality were acquired from Chennai Water Metro Board/Central Ground Water Board (CGWB) during the pre-and the post-monsoon seasons using fixed sampling (because of fixed tube wells) technique. Forty-four groundwater samples were acquired of May 2005 (pre-monsoon) and sampling activity was repeated during January 2006 (post-monsoon) and it is a regular monitoring. Out of 44 samples points, most of wells are situated near agricultural location while some points were along village and town. Information related to date of collection/sampling, vegetation and cropping system, location of nearby ponds etc. was also recorded.

Computation of WQI
For computation of WQI, 13 quality parameters have been chosen, analysis involved in two steps; a single analysis for each parameter and the WQI ( WQI 1 , WQI 2 and WQI 3 ). Standards of water quality recommended by the WHO (2006) have been used to calculate WQI.
The weighted arithmetic index method (Brown, McClelland, Deininger, & Tozer, 1970;Singh et al. 2013bSingh et al. , 2015Singh et al. , 2016 was used for the calculation of WQI. Further, quality rating or sub index (Q n ) was calculated using the following expression Equationsn (1-3): Q n = Quality rating for the nth water quality parameter; V n = Actual value of nth parameter; v i = Ideal value of this parameter; V s = Standard permissible value of the (Table 1); [Consider V i = 0 for all except pH where V i = 7 for pH] Unit weight was calculated by a value inversely proportional to the recommended standard value V s of the corresponding parameter. (2)

Variability analysis of WQI n
For variability analysis of WQI n (n = 1, 2, 3) during study period is segregated into three parts (i) average WQI n of the pre-monsoon, (ii) average WQI n of the post-monsoon and (iii) average WQI n (pre + post-monsoon) during study period. For creating special destitute map of WQI n of study area, IDW (Inverse Distance Weightage) interpolation method was used in GIS environment.

GIS and inverse distance weighted (IDW)
The boundary of the study area was digitized from Survey of India toposheet using ArcGIS software. The precise locations of monitoring wells were recorded using GARMIN 12-Channel GPS and later are imported in GIS platform. The computed WQIs obtained from the water quality parameter were used in GIS and thematic maps were created with the help of spatial analyst module available in ArcGIS 10 software. Inverse distance weighted (IDW) interpolation technique was used for generating the WQIs thematic maps. IDW is versatile, easy to use program, and fairly accurate under a wide range of conditions and capture the extent of local spatial variation . By applying IDW method, the property at each unknown location for which a solution is seeking can be expressed as mathematically by Equation 5 . The expression for IDW can be depicted using the Equations 6, 7, and 8: where λ i -is the weight for neighbor i (the sum of weights must be unity to ensure an unbiased interpolator); d(S o , S i ) -is the distance from the new point to a known sampled point; β -is a coefficient that is used to adjust the weights and n -is the total number of points in the neighborhood analysis.

Hydrogeochemistry
The descriptive statistical results of the pre-monsoon of the year (2005-2013) are given in Table 3(a). EC varied from 180 to 12240 dS/m with a mean of 1322.90 dS/m. The range, mean and SD of pH are 6.80-8.6, 8.37 and 0.40. The mean cations concentration followed a ; > 1 where W n = unit weight for the nth parameters; K = constant for proportionality; V s = standard value for the nth parameters (Table 1).

Nemerows pollution index (NPI or WQI 2 )
It is a simplified pollution index introduced by Neme (Mohan, Singh, Panday, Kumar, & Jain, 2007) which is also known as Nemerow's Pollution Index. It is given as by Equation (4):

Soluble sodium percentage (SSP or WQI 3 )
Excess Na + , combining with carbonate, leads to formation of alkali soils, whereas with chloride, saline soils are formed. Neither soil will support plant growth (Rao, 2006;Singh et al., 2013). Generally, %Na + should not exceed 60% in irrigation waters. Sodium content is generally expressed in terms of per cent sodium or soluble sodium percentage (WQI 3 ). Wilcox (1955) planned a method for rating the irrigation waters based on the percentage of sodium. The percentage of sodium in water can be calculated using Equation (5) as follow.
Sodium percentage is an important factor for studying sodium hazard. It is also used for adjudging the quality of water for agricultural purposes. High percentage sodium water for irrigation purpose may stunt the plant growth and reduces soil permeability (Joshi, Kumar, & Agrawal, 2009;Singh et al., 2015). the pre-monsoon (July) 2012, 25 percent samples were under excellence limit and these percentage increased up to 84 percent due to only chemical dilution process by rainfall (during 2005 good amount of rainfall receive by study area, Figure 1(b)). Similarly, during year 2008 pre-monsoon (July), 22 percent well's come under < 10 limit but after rainfall (good amount of rainfall events so more groundwater recharge, Figure 1(b)) and dilution process this percentage increase up to 43%during postmonsoon year 2009. Due to rainfall factor, year 2009 the pre-monsoon 9% increase up to 68% (  (2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013). The mean concentration of Mg 2+ (50.26 mg/l) in the pre-monsoon is higher than recommended WHO (2006) as 50 mg/l.
Results of the post-monsoon year (2006-2013) are given Table 3(b). The results of the descriptive statistics revealed that there is no statistically significant differences between the mean concentration of the parameters values between June (pre-monsoon) and January (post-monsoon), only few degree of decreasing except mean value of K + (increasing 0.93 mg/l with respect to the pre-monsoon value of K + (18.46 mg/l). The mean cations concentration is followed a descending order as: Na + (125.56 mg/l) > Ca 2+ (62.74 mg/l) > Mg 2+ (47.93 mg/l) > K + (19.39 mg/l). The mean concentration in groundwater during the post-monsoon year (2006-2013) followed a descending order as: 244.94, 223.54, 76.04, 28.70, 7.15 and 0.39 mg/l, respectively. The mean concentration of pH ranged between 8.37 ± 0.40 at pre-monsoon and 8.20 ± 0.43 at the post-monsoon (Table 3(a) and (b)).

Post-monsoon water quality index (WQI)
Table 4(a) shows the WQI value of pre-monsoon. The post-monsoon (month of January) years wise all the samples with their category are given in Table 4(b) and Table 5 ((a)-(c)). During the post-monsoon, concentration of chemical effect slow down the process (dilution process by rainfall) therefore more number of samples falls under first category according to WQI 1 with respect to the pre-monsoon WQI 1 . During the post-monsoon year 2006 largest percent (84 %, total 37 well's water) of samples come under excellent category because of during  in study area, therefore, almost for all well WQI 1 values has been reduced. Figure 2(b) and Table 4(b), it is observed that WQI 2 also capable to show dilution process effect over wells during post-monsoon years. During the post-monsoon time series only seven time different wells as well no. 6 (2008), 24 (2006), 27 (2007, 2010 and 2013), 34 (2007) and 36 (2011 and 2012) are comes under third category (20-30) of WQI 2 (Figure 2(b) and Table 5(a)) otherwise all wells come under excellence and good class according to WQI 2 scaling method. According to average values throughout year WQI 2 values exist under first category (Table 5(a)). Figure 2(c) and Table 5(b), reveal that after the pre monsoon due to rainfall factor mostly no. of wells come under first and second class of WQI 3 . Impotency of rainfall factor can be understand as that during the pre-monsoon total 44 no. of wells were (Table 5(c)) come under unfit or unsuitable category for irrigation but due to rainfall factor these 44 well reduce up to 8 (6 in poor and 2 in unfit class, Table 5(c)). On an average during the post-monsoon wells were exist under second class of WQI 3 (Table 5(c)).

Pre-monsoon water quality index (WQI)
Year wise all samples with their corresponding category are given in Table 4(a) and Table 6 Table 5(a), clearly reveal effect of rainfall phenomena over WQI, after rainfall any well's water not come under unsuitable category.
Table 5(b), most time of study period much no. of wells are found under excellence (< 10) and good (10-20) category some time few no. of wells come under third category (20-30) during the post-monsoon years but it is observed with this index that is after rainfall no well not exist in fourth and fifth category of WQI 3 it is only due to dilution process by rainfall.
Table 5(c), it is clearly observed effect of rainfall over WQI value due to dilution process during post monsoon. According to Table 5(c) after pre-monsoon, any well not exist in category 30-40 (forth) and > 40 (fifth). After rainfall leaching of Ca 2+ , Mg 2+ , Na + and K + are more chances from rock but its net result give high yield in dominator of Equation 5. Its result will be come as less WQI 3 value, due to this reason after rainfall mostly wells water come under suitable classes for irrigation purposes.      (Table 4(a)) reveal that most of study period (2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013), area come under good (40-70) category therefore groundwater of study area is suitable for irrigation purposes (Shaki & Adeloye, 2006).  (Table 6(b)), shows large number of samples comes under excellent category (<10) because WQI 2 index is based on ratio and it is very rare change for sum of ratio, cross first limit (<10) barrier due to this fact WQI 2 cannot touch > 40   2− = 612/400 and TDS = 7341/1000) parameters were so high that why WQI 2 value (66.70) is higher than 2007 WQI value (55.37). Throughout the study period (2005-2013) average value of WQI 2 for each comes under excellence (< 10) category (Table 4(a)).
From Figure 3(c) and Table 4(a), most time of study period much no. of wells are found under excellence (<40) and good (40-70) category some time WQI 3 hit to fourth (80-90, very poor category) category. WQI 3 is basically depend on sum of concentration of Ca 2+ , Mg 2+ , Na + and K + parameter, and which year this sum will be lowest this year WQI 3 will be highest or touch to fifth category of WQI 3 . During the pre-monsoon year 2005, well no. 25 show value of WQI 3 above 80 because of sum of three parameters (Ca 2+ (=6)+ Mg 2+ (=12)+Na + (=104) and K + (=6)) is 128 which is vey less than numerator value (104 × 100) of Equation 6, therefore ratio of numerator and denominator comes under the forth category limit (80-90, very poor, Figure 3(c) and Table 4(a)). Similarly, well no. 4 come under forth category during 2006, 2008, 2009, 2011, 2012 and 2013, and in these years numerator (28500, 38000, 34000, 36600, 35200 and 35200) value was so higher then denominator value (348, 457, 402, 456, 431.59 and 434.67). After post-monsoon concentration of chemicals in groundwater will likely to increase due to leaching processes. Well no. 12 and 25 are poor for irrigation point of views during 2009 and 2013 respectively, according to WQI 3 . Overall average value of WQI 3 at each year (2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013) comes under good (40-70) category during the pre-monsoon (Table 4(a)). TDS (1075 mg/l), F − (2.17 mg/l) and Na + (340 mg/l). In same manner well no. 23 also come under unsuitable limit (>100, Figure 3(a)) during 2009 because this year four parameters (pH = 8.8, EC = 1430 dS/m, Mg 2+ = 49~50 mg/l and F = 1.6 mg/l) out of 13 parameters cross the WHO (2006) standard limit and second reason is two parameters as pH (W pH = 0.134) and F -(W F = 0.759) out of four parameters (pH, EC, Mg 2+ and F − ) have high weight with respect to other parameters therefore WQI 1 having high value for well no. 23. Figure 3(a) and Table  4(a) also reveal that well no. 34 come under the fifth class (>100) because six chemical parameters EC (1880 dS/m), Ca 2+ (80 mg/l), Mg 2+ (83 mg/l), Cl − (383 mg/l), TDS (1019 mg/l) and F − (1.78 mg/l) are above the WHO (2006) standard values during the pre-monsoon year 2009 and during 2007 this well was also polluted because this time seven parameters (EC = 12240 dS/m, Ca 2+ = 380 mg/l, Mg 2+ = 1021 mg/l, Cl − = 4219 mg/l, Na + = 823 mg/l, SO 4 2− = 612 mg/l and TDS = 7341 mg/l) were above the standard limit WHO(2006). Well no. 35 was also unsuitable (Figure 3(a)) for drink during the premonsoon year 2013, because present of heavy weightage parameter (pH and F). Figure 3(b) and Table 4(a) reveal that according WQI 2 water of well no. 34 is unfit for drinking purposes because WQI 2 index is sum of normalized values of individual groundwater parameters therefore most of time WQI values close to 1 or near to 1, when numerator will be near to WHO (2006) and when numerator is so high with respect to denominator than WQI 2 value will be so high like well no. 34 during the pre-monsoon 2007, six parameters (EC = 10400/1400, Mg 2+ = 742/50, Na + = 1380/200, Cl − = 3368/250, SO 4 2− = 960/400 and TDS = 6839/1000) out of 13 parameter having high value of numerator with respect to denominator S n , standard value of parameter according to WHO (2006). Average value of WQI 1 for each year comes     part of study area come under second category of WQI 2 . According to WQI 2 scaling ways study area existing good area of Blue patch (or excellent category). WQI 2 is not capable to pick or representing accurately third, fourth, and fifth category of water because this method is based on normalization process, after normalization it is very rear for any well to achieve third fourth and fifth category. The map of WQI 2 for the post-monsoon eight years average (Figure 4(b)-B) is indicating third category of WQI 2 is increased while first category of WQI 2 in decreasing mod with respect to Figure 4(b)-A. Second category of WQI 2 is almost remaining constant which is reveal that some part of Blue colour shifted toward second class. From observation of Figure 4(b)-(C), we can conclude that due to rainfall, leaching process is high in study area therefore nominator of Equation 4, is become so higher than dominator (S n ) which give higher value of WQI 2 during post monsoon and this reason also explain to redaction in Blue colour patch in study area during post-monsoon. Figure 4(b)-C it can be conclude that by averaging process of (pre-post-monsoon) full time series data sets (2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013) for calculating WQI 2 , mathematically chemical concentration are rearrange such a ways that they reduces effect of high concentration values therefore Figure 4(b)-(C) shows large area of second category (or light green patches). Figure 4(b)-(C) reveal average effect of the pre-and the post-monsoon reduces the permissible category in to good category. Over all WQI 2 having capability to show/pick leaching process during the post-monsoon because WQI 2 give more importance to individual contamination value of parameter.

Variability of groundwater based on WQI 3
The nine year the pre-monsoon average WQI 3 of 44 wells are interpolated for irrigation suitability point of views represent by Figure 4(c)-A. Figure 4(c)-A, reveal that almost 50 percent part of study area can easily found suitable groundwater (blue colour patch (or excellences category)+light green patch (or good category)) for irrigation during pre-monsoon. Figure 4(c)-A indicating study area also having a high Na% area (red patches) while this unsuitable category lies into very poor category (Fig. 4(c)-B), due to rainfall more leaching of cations (Ca 2+ , Mg 2+ , Na + and K + ) therefore dominator of Equation 5 become strong in compare of numerator its result low value of WQI 3 during post-monsoon. From Figure4(c)-C, red patches (or high Na%) area increases in compare of the post-monsoon but less than the pre-monsoon due to net effect (pre + post monsoon) most part of study area easily found suitable irrigation water.

Discussions
There is less significant differences in WQIs values were observed among the 44 sampling sites over the study year (Table 4(a), (b)) shows poor to medium and good

Variability of groundwater based on WQI 1
The nine year the pre-monsoon average WQI 1 of 44 wells are interpolated and resultant as a distributed WQI 1 map represent by Figure 4(a)-A. Figure 4(a) (A), reveal that most part of study area comes under the second category (good class) which is illustrated in map by light green colour. Figure 4(a)-A clearly indicated that after 25-50 range 50-70 range (forth class of WQI 1 or permissible class) is dominate in study area which is illustrated by yellow colour patches. Blue or excellent class area in study area is very small in compare of yellow that mean during study period (2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013) less amount of excellent category was available in study area. Interpolated map of WQI 1 (Figure 4(a)-A) during study period, reveal that few area of study area are comes under groundwater polluted area. Maximum and minimum average value of WQI 1 during years 2005-2013 were 69.9 (under permissible limit) and 21.7 (excellent), respectively. Figure 4(a) (B), represents special destitute WQI 1 map of eight year (2005-2013) the pre monsoon average WQI 1 for 44 wells. According to WQI 1 scaling and IDW method, study area showing clear effect of rainfall factor over groundwater because during the post-monsoon Blue colour area (excellent/first class) increased and light green area (good/second class) decrease or may be merge into Blue coloured area due to dilution process by rainfall. From Figure 4(a)-A and (a)-B, clearly indicated that rainfall effect is more over yellow colour occupied area (permissible/third class). During the eight year, the post-monsoon period WQI 1 average value increased for class fourth and fifth, because every time these area GW's chemical parameter existing in higher ranges from WHO (2006) standard values and during rainfall leaching process is much from leather factory (near of big black patch). Figure 4(a)-C represents net effect of average WQI 1 the pre and the post-monsoon of year 2005-2013. Figure  4(a)-C, reveal that due to net effect over excellence category of groundwater is almost vanish or may be due to net chemical concentration effect Blue colour patch converted/merge into third category (or yellow color patches). Net yellow colours patch (permissible/third class) area are more in compare to average pre and the post-monsoon WQI 1 , this is indication of great changes in first category due to more pollution/contamination of GW. Second category area almost reaming constant but very poor (fourth) and unfit (fifth) category become reduces in compare of Figure 4(a)-A and (a)-B. Over all according WQI 1 scaling, study area wells water can be use for drinking before some degree of treatment.

Variability of groundwater based on WQI 2
The map (Figure 4(b)-A) is average of nine years the pre-monsoon WQI 2 for study area. Figure 4(b)-A, clearly indicates that during last nine the pre-monsoon year most   and provides a composite numerical value at a defined scale. The WQIs has potential for decision making and management in flag contaminate of concern, predicts potentially harmful conditions, guides prioritization management efforts and funds, assess overall impacts of water quality interventions, communicate impacts of management and policy decisions, monitor water resource and health. The WQI output numeric values at defined scale suggest the water quality status with respect to excellent to severely polluted. In this study, physicochemical parameters have been used in process of spatial and temporal pattern of WQI n for the study area. During observation we found that WQI 1 is more reliable than WQI 2 because for WQI 1 weightage factor (W s ) will be constant always if each time chemical parameter of groundwater is constant in number while WQI 1 can show variability due to individual parameter concentration. Advantage of WQI 2 is scaling, it is capable to show the effect of leaching of particular (or individual) parameter. WQI 3 index can be used as irrigation parameter due to its easiness. Even though the WQIs tool is quick and reliable but they do not replace analysis of detailed water quality and other environmental data. These are dependent on these detailed data to drive the algorithms behind the composite numeric values. For improved prediction of water quality the WQIs can be used together with such indices of biological integrity, in order to give a more comprehensive view of water resources health. Spatial thematic map of WQIs shows the utility of GIS technique. and also excellent at few monitoring locations. Slight variations were observed due to seasonality.
According to the classifications for the different water uses, high values indicate that the water requires treatment before use for irrigation or other purposes. Longterm data-sets are required to confirm the environmental effect of water pollution. It should be noted that several authors have reported the polluting effect of rain and climatic conditions on the dynamics and water quality of rivers and reservoirs. The rainfall changes the overall WQI value due to dilution effect. The water which has higher value of WQI needs treatment before irrigation or other purposes.
The lowest WQI values were observed in the northern part of the study area, while the highest water quality values were found in the southern part. The water quality of any study area is notably affected by anthropogenic activities such as agricultural production, the industry, and urban growth (Singh et al., 2012. Many researchers suggested that anthropogenic activities in their study area have been responsible for changes in the groundwater quality (Singh, Tiwari, Panigary, and Mahato, 2013;Singh et al. 2015;Thakur et al., 2013;Thakur et al., 2015), sediment quality  and surface water quality .
Higher concentration of nitrate is observable in the rural part of the study area due to excess input of chemical or mineral fertilizer, improper manure management practices and improper operation and maintenance of septic systems (Singh et al., 2012). The higher value of WQI at study stations have been reported due to higher values of iron, nitrate, TDS, hardness, fluorides, bicarbonate and manganese in the groundwater (Ramakrishnaiah, Sadashivaiah, & Ranganna, 2009).
The calculated WQI showed fair water quality rating in post-monsoon season later it changed to medium in summer and winter seasons mainly due to human activities, soil process and climatic effects in the study area (Rajankar, Gulhane, Tambekar, Ramteke, & Wate, 2009). Singh et al. (2016) have studied the Bhopal lake water through WQI methodology and their finding suggests that generally, pre-monsoon samples have higher pollution level due to more evaporation of lake water, input of excess of domestic waste water and less dilution.

Conclusion
Groundwater is a natural resource, the evaluation of groundwater using the long-term water quality helps in developing new strategies. This study demonstrates that the use of WQIs and GIS methods together, they could provide useful information for water quality assessment. WQIs have potential in assessing and communicating overall impacts of existing, planned, or proposed water quality interventions and management decisions. The WQIs approach helps in categorization the water based on their physiochemical parameters