Identifying Groundwater and River Water Interconnections Using Hydrochemistry, Stable Isotopes, and Statistical Methods in Hanumante River, Kathmandu Valley, Central Nepal

: Interconnection between river water and groundwater plays an important role in maintaining water quantity and quality in hydrological systems. Furthermore, the exact interconnection is often difficult to observe and measure. This study attempts to explain river and shallow groundwater interconnection in urbanized areas of the Kathmandu Valley, Nepal. Isotopic ( δ D and δ 18 O) and chemical analyses were performed on river and groundwater samples, and the results were analyzed using statistical methods to identify areas of interconnection between river water and groundwater. Higher concentrations and positive strong correlations of Na + with K + , NH 4+ -N, Cl − , HCO 3 − , and PO 4 − -P, and a change of water type from Ca-HCO 3 during the wet season to Na-K-HCO 3 during the dry season indicate higher contamination in river water during the dry season. Hierarchical cluster analysis was used in grouping water samples into clusters on the basis of isotopic and chemical (Na+ and Cl − ) composition. Grouping of river and groundwater samples in one–one clusters from wet and dry seasons shows the presence of interconnection, indicating the contribution of river water in recharging shallow groundwater. These results imply that shallow groundwater found near rivers is chemically contaminated by polluted river water through bank infiltration, in both wet and dry seasons


Introduction
Interconnection of river water and groundwater is a process of exchange between waters located on the river channel with those in the rocks/sediments under the surface.The exchange rate of water is controlled by hydraulic conductivities of the river channel and aquifer sediments; the relative stage of the river channel and nearby groundwater level; and geometry of the river channel within the alluvial plain [1,2].The presence of a clogging layer on an aquifer and a riverbed or bank can decrease or stop the water flow exchange [3].Interconnection of river water and groundwater also depends on the distance from a river channel, the geological conditions, and climatic factors [4].Understanding of groundwater and surface water interconnection is very important to develop effective water resource management and policy as it can change the water quality and quantity of both water systems [5,6].
Studies related to the interconnection between river water and groundwater has increased in most developing countries during the last few decades.Utmost studies have been carried out to the assess areas of interconnection and the presence of exchange flow between river and groundwater on a regional as well as local scale [7][8][9].Essentially, two types of exchange flow conditions are involved in the river and groundwater interaction: 1) the influent condition and 2) the effluent condition.Based on the condition of exchange flow, several studies specify that base flow in the river during the dry season is the result of effluent flow from shallow groundwater [5,[10][11][12][13].Furthermore, the condition of the exchange process can be affected by anthropogenic activities that alter the exchange processes, reduce connectivity, and lead to chemical or biological contaminations [2,14].Increased sewage load into the rivers running through the urbanized cities can transfer toxic contamination to surrounding shallow groundwater in influent reaches.The decline of water table in a nearby shallow aquifer due to over-extraction can increase groundwater recharge from polluted river water [15,16].The anthropogenic activities driven by increased urbanization and population growth affect river water quality, which adversely reflects on nearby groundwater quality.
Previous research conducted in the Kathmandu Valley reported the presence of an interconnection between river water and groundwater, showing a recharge of the groundwater by river water [17][18][19].Additionally, a number of previous studies on river water quality have also reported heavy contamination of the downstream section of major rivers [20][21][22][23], inducing the occurrence of various water-borne diseases such as diarrhea, cholera, and dysentery among the people of riverside areas [24].
Methods such as heat tracer, solute tracer, direct measurement of water flux, and environmental tracer methods including isotope and geochemistry have been used to determine the interconnections of groundwater and surface water [25][26][27][28][29][30].Similarly, numerical modeling, geophysical methods, and statistical methods have also been used to describe interaction processes [31][32][33][34].Basically, water flux measurement mostly used a direct method to get exchange flow conditions involved in interacting processes.The recharging source for groundwater or river water is dependent on the direction of exchange flow [35].Further, stable isotopes of hydrogen and oxygen along with Na and Cl ions have been widely used to determine river and groundwater interconnections.Similarity in isotopic values and chemical ions between nearby groundwater and river water indicate the presence of river and groundwater interconnection [13,36].Additionally, hierarchical cluster analysis (HCA), one of the multivariate methods, has been widely used to determine the interaction of river water and groundwater at the absence of a water flux dataset [37].HCA can be used to analyze regional-scale as well as local small datasets.Several previous studies used HCA as a major statistical method and suggested similarities in chemical and isotopic compositions between river water and their nearby groundwater representing the presence of river and groundwater interaction [9,34,37].
The present study will focus on the application of stable isotope values, chemical compositions of the river and groundwater samples, followed by statistical analysis for one of the contaminated rivers in urban areas of the Kathmandu valley.Thus, identifying the occurrence of the spatial and seasonal interconnectivity and possible contamination load from the river to river periphery groundwater are major goals of this study.

Study Area
The Hanumante River is a centripetal river and is one of the most polluted tributaries of the Bagmati River in the Kathmandu Valley, Central Nepal [20,38].It is the only river that drains from the eastern part of the Kathmandu Valley and confluences with the Manahara River (a major tributary of the Bagmati River) at Jadibuti (Figure 1).This sixth-order river extends up to 18.29 km, covering nearly 97 km 2 of watershed areas [39].The Godawari Khola, Tabyakhusi Khola, and Chakkhu Khola are major tributaries of this river.
The quality of the Hanumante River water is deteriorating as urbanization increases downstream, having lower dissolved oxygen (0-7 mg/L) and higher biological oxygen demand (3.5-79.9mg/L), chemical oxygen demand (128 mg/L), ammonia (0.4-25 mg/L), and phosphorous (0.09-1.71 mg/L), such that the river water is harmful for domestic purposes [38].Direct disposal of sewage and solid waste effluent from industries converts the Hanumante River into an open sewer during the dry season [40].The sole municipal drinking water supply organization, namely, Kathmandu Upatyaka Khanepani Limited (KUKL), cannot fulfill the total water demands of the Kathmandu Valley [41], compelling it to fulfill the water deficit by extracting groundwater from shallow as well as deep aquifers.Meanwhile, Gautam et al. [17] discussed the high possibility of shallow aquifer contamination by polluted river water in the peripheral part of the rivers in the Kathmandu Valley.Geologically, Plio-Pleistocene fluvial, fluvio-lacustrine, and fluvio-deltaic sediments comprise the Kathmandu Valley [42,43].The study area presented in Figure 1 is characterized by four formations, namely, the Gokarna Formation, Patan Formation, Thimi Formation, and Lower Terrace Deposit [44].The upstream section of the study area is composed of the Gokarna Formation, containing dark brown colored, laminated arkosic sand, silty clay, and peat.The middle section of the study area is covered by the Thimi Formation, which consists of sand, silt, clay, peat, and gravel composed of granite and gneiss derived from the Shivapuri Range.Similarly, the lower section of the study area is dominated by the Patan Formation, which contains deposits of fluvial-lacustrine composed of sand, silt, clay, and peat.The Lower Terrace Deposit along the river corridor consists of micaceous sand, pebbles, and granules [44].
To use chemical and isotopic analysis to investigate any interconnectivity between river and groundwater, samples were collected from 20 locations-10 from rivers and 10 from dug wells (Figure 1) during the wet (August 2017) and dry seasons (February 2018).

FieldMeasurement and Water Sampling
Groundwater was collected from dug wells which are located within 10 to 100 m from the river channel, with the depth ranging between 1.5 and 15.7 m.Sampling was carried out in two seasons-August 2017 (wet season) and February 2018 (dry season)-where the samples were collected from 20 locations, 10 from rivers and 10 from dug wells (Figure 1) in the consecutive seasons, respectively.Water samples were collected in 100 mL polyethylene bottles.Each bottle was rinsed three times with the same water before sample collection.Groundwater samples were collected after removing a quantity of water using an installed hand pump or with the help of rope and a plastic bucket.The collected water samples were stored at −4°C at a laboratory until the chemical and isotope analyses were performed Additionally, during the sample collection, well depth, water level depth, electrical conductivity (EC), dissolved oxygen(DO), pH, and water temperature were measured at each sampling location.
Well depth was measured using a measuring tape and was verified with the dug well owner.A water depth logger was deployed for the water table measurement.In situ parameters were measured by using portable devices, namely, a DO meter (Mettler Toledo SG3-ELK, Greifensee, Zurich, Switzerland) and a pH/EC meter (Mettler Toledo Duo, Greifensee, Zurich, Switzerland).The location of the water samples is shown in Figure 1, and Table 1presents the data measured during the field survey.

Chemical and Isotopic Analysis
The Interdisciplinary Center for River Basin Environment, University of Yamanashi (ICRE-UY), Japan provided laboratory facilities to carry out chemical and isotopic analyses.As per the laboratory procedure, collected water samples were first filtered through 0.2 µm filter paper to prepare final samples for further laboratory analyses.The dominant chemical ions, including cations (Na + , K + , NH4-N + , Ca 2+ , and Mg 2+ ) and anions (Cl − , NO3 − -N, PO4 − -P, and SO4 2− ), were determined by using ion chromatography (ICS-1100, Dionex, Waltham, MA, USA with an analytical error of 5%.The bicarbonate ion (HCO3 − ) concentrations were measured by using a titration method with 0.01N sulfuric acid.
The stable isotopes of hydrogen (δD) and oxygen (δ 18 O) were analyzed using cavity ring-down spectroscopy (L1102-i, Picarro, Santa Clara, CA, USA).VSMOW (Vienna Standard Mean Ocean Water) is the standard water used to calculate isotopic ratios (δ) of D and 18 O of water samples.The results were reported in parts per thousand (per mill deviation) with respect to these standards with precision 0.5‰ for δD and 0.1‰ for δ 18 O.The isotopic ratios of hydrogen and oxygen were calculated by using the formula given by Craig [45]: 16 O in sampled water (Rsample) and standard mean ocean water (Rstandard).

Statistical Analysis
Temporal variations of chemical variables were evaluated using a paired t-test for significant difference in parameters [46,47] within a 95% confidence level.Spearman's rho correlation analysis [48] was adopted to establish any relationship among different variables.Hierarchical cluster analysis (HCA) was used to examine any similarity in chemical as well as isotopic composition between river water and groundwater.Cluster analysis is useful in distinguishing water showing similar chemical or isotopic composition from dissimilar ones [37,49,50].HCA was performed based on Ward's linkage method [51] with squared Euclidean distances as a measure of similarity between samples [11,40].Statistical Package for Social studies version 25 (SPSS Inc., Chicago, IL, USA) was used for statistical analyses.

InSitu Parameters
Table 1 presents locations of sample points, the depth of water level at different wells, and the measured in situ parameters in wet and dry seasons.All sample wells had shallow water depth in the wet season, which may imply high recharge and lower extraction rates during wet seasons, with a maximum fluctuation of 3.60 m and a minimum fluctuation of 0.65 m at HW9 and HW6, respectively (Table 1).
In situ parameters, namely, temperature, pH, EC, and DO, in river water and groundwater were measured in wet and dry seasons.The temperature of river water ranged from 22.3 °C to 24.3 °C in the wet season, while in the dry season the range was 14.0-17.9°C.The temperature range of groundwater was 20.7-24.1°C in the wet season and 13.4-20.5°C in the dry season (Table 1).The pH value slightly decreased during the dry season in both river water and groundwater (Table 1).EC measured in groundwater ranged from 290 to 934 µS/cm in the wet season, and from 576 to 1323µS/cm in the dry season.Groundwater exhibited higher EC in the dry season relative to the wet season, except at HW1, HW2, and HW4 (Figure 1, Table 1).However, in the case of river water, the value of EC was low (164.2 to 247 µS/cm) in the wet season and abruptly increased by up to eight times (604 to 2060 µS/cm) in the dry season.DO was high in river water during the wet season and abruptly decreased below the value measured in groundwater during the dry season (Table 1).
Ca 2+ and HCO3 − were the most dominant ions of river water in the wet season, with values ranging from 7.3 to 17.0 mg/L and 24.4 to 73.2 mg/L, respectively (Figure 2).The major cations had an order of Ca 2+ > Na + > K + > Mg 2+ in the wet season, which changed to the order Na + > NH4 + > Ca 2+ > K + > Mg 2+ in the dry season.Similarly, for anions, HCO3 − was dominant, followed by SO4 2-and Cl -in the wet season, and Cl − and SO4 2− in the dry season.Except for NO3 − -N, all other parameters showed strong significant seasonal variation (p < 0.01).Concentrations of all these parameters increased in the dry season, but the rate of increment varied for different parameters.NH4 + -N and PO4 − -P concentration was insignificant (<1 mg/L) in the wet season and significantly increased in the dry season, ranging from 10.3 to 102.9 mg/L (for NH4 + -N) and from 2.4 to 31.7 mg/L (for PO4 − -P).Similarly, concentrations of Na + , K + , Cl − , and HCO3 -increased by more than ten times than in the wet season (Figure 2).Based on piper plot from the chemical analyses, the Hanumante River can be categorized as Ca-HCO3 type in the wet season and Na-K-HCO3 type in the dry season (Figure 3).Likewise, Ca 2+ and HCO3 − were the most dominant ions of groundwater in both wet and dry seasons.Concentration of Ca 2+ ranged from 11.7 to 72.5 mg/L in the wet season and from 28.4 to 93.7 mg/L in the dry season.In the case of HCO3 − , 67.1-305.0mg/L was the range in the wet season, which increased to 122.0-579.5 mg/L in the dry season.Statistical analysis (paired t-test within a 95% confidence level) showed no significant temporal variation in K + , NH4 + -N, Ca 2+ , HCO3 − , NO3 − -N, and SO4 2− in groundwater.However, for Na + , Mg 2+ , and Cl − , temporal variation was significant, with a p-value of 0.03.In general, groundwater showed lesser increments in concentration compared with river water during the dry season (Figure 2) and is classified as Ca-HCO3 type in both dry and wet seasons, except at HW5, HW9, and HW10.Water samples collected from HW5 and HW10 changed slightly from Ca-HCO3 (in the wet season) to Ca-SO4, and Na-Cl-SO4, respectively, during the dry season (Figure 3).Groundwater collected from HW9 in both seasons is of Ca-SO4 type.
Determination of water types using a piper diagram suggests the origin of the water [52].Ca-HCO3 type represents recent infiltration of freshwater, whereas Ca-SO4 and Na-K-HCO3 types indicate water exhibiting simple dissolution or mixing and ion exchange, respectively [53].Groundwater and river water from the wet season of Ca-HCO3 type thus represents recent rainfall infiltration or runoff as a major contributing source for groundwater recharge and river discharge.Changes in river water type during the dry season reflect the presence of different water sources for river discharge, including direct discharge of untreated sewage from municipal and industrial sources [16].
The Spearman's rho correlation matrix between different chemical parameters in river water and groundwater is presented in Table 2.In river water, DO has a strong negative correlation (r = −0.60 to −0.86) with all parameters except pH and NO3 − -N.The negative correlation of DO indicates that the presence of a higher concentration of chemical parameters decreases the amount of dissolved oxygen in river water.In contrast, EC has a strong positive correlation with parameters which have a negative correlation to DO, suggesting full dependence of conductivity of river water on dissolved ion concentrations [50].The higher value of EC and positive correlation with most ions also indicate higher anthropogenic contamination during the dry season [54] (Table 1 and Figure 2).There is also a strong positive correlation of Na + , K + , NH4 + -N, Ca 2+ , Mg 2+ , Cl − , HCO3 − , PO4 -P, and SO4 2− between each other (r = 0.67 to 0.98), indicating that the river water is highly influenced by anthropogenic pollution [11] such as direct discharge of municipal and industrial sewage and leachate of solid waste disposal near the river channel during the dry season [40].Strong positive correlation of PO4 − -P with SO4 2− , NH4 + -N, and NO3 − -N also suggests the influence of fertilizer and pesticides used in the cultivated land of river peripheral areas [54][55][56].
Correlations of chemical parameters of groundwater are similar to those of river water.As in the case of EC, it has strong positive correlation with Na + , K + , NH4 + -N, Ca 2+ , Mg 2+ , Cl − , and PO4 − -P.Strong positive correlation between Na + , Mg 2+ and Cl − (0.94) in groundwater indicates the influence of anthropogenic activities [57] because there was no evidence of halite deposits in the study area [58].The negative correlation between NH4 + -N and NO3 − -N represents nitrification of NH4 + -N into NO3 − -N [59].Positive correlation of K + , NH4 + -N, and PO4 − -P represents agricultural impact from the surrounding cultivated land.

Isotopic Composition
Fewer studies regarding isotopic analysis of meteoric water have been carried out on a local scale [60][61][62].The local meteoric water line (LMWL), established by Gajurel et al. [60] and Giri [61], has a similar slope and intercept value as the global meteoric water line (GMWL) reported by Craig [45].In this study, GMWL is used as a reference of meteoric water.The δ 18 O verses δD plot (Figure 4) presents variations in stable isotopic compositions of groundwater and river water during both wet and dry seasons.During the wet season, almost all river samples showed similar isotopic compositions of δ 18 O (−8.0 to −8.1‰) and δD (−54.3 to −55.7‰) except HR8, which had a heavier isotopic composition (Figure 1, Figure 4).Composition of isotopes slightly changed to heavier with wider value range during the dry season, showing a range from −6.91 to −7.73‰ for δ 18 O and −48.75 to −54.41‰ for δD.Research work by Yang et al. [63] in the Jiulong River also presents a narrow and wider value range for the wet and dry season, respectively.Wet season river samples are plotted near the GMWL and have a similar slope and intercept (R 2 = 0.86) as the GMWL (Figure4a), indicating recent meteoric water as a major source for river discharge [4,8,34,58].Ca-HCO3 water type, defined from piper diagram and lighter isotopic composition of rainfall during wet season [62,64], also suggests a similar water source.However, there was no evidence of evaporation in a previous study in the Kathmandu Valley [58,60]; dry season river samples plotted below the GMWL with lower slope (slope = 5.46 with R 2 = 0.79; Figure 4b) as compared with the GMWL may indicate a possibility of evaporation [26,65].Large spatial variation of δ 18 O and δD in water from dug wells was observed in both wet and dry seasons (Figure 4b).The δ 18 O of groundwater ranged from −6.1‰ to −8.2‰ in the wet season and from−6.68‰ to −8.85‰ in the dry season.However, the overall range of δ 18 O was similar in both seasons, with samples from HW3, HW5, HW6, HW7, HW8, and HW10 showing lighter composition in the wet season.In both seasons, HW1 had the heaviest isotopic value, whereas HW2 had the lightest isotopic composition.Groundwater samples plot below the GMWL, with a similar slope (R 2 > 0.9) in both wet and dry seasons (Figure 4b), indicating recent meteoric water as a major recharge source for these dug wells [12].Although meteoric water is a major source for dug wells, spatial variation is noticeable in groundwater since isotopic composition of precipitation is dependent on rainfall amount, elevation, and source of water vapor of rainfall [62,64,66].

Clustering of River Water and Groundwater
Isotopic composition is a reliable source to identify recharge sources of groundwater, while the concentration of Na + and Cl − can be used as an indicator of the presence of contamination through increased urbanization [36].Thus, for HCA, δD, δ 18 O, Na + , and Cl − are used as major parameters to identify similarity or interconnection between river water and groundwater [4,30].All 20 water samples from river water and groundwater of both seasons are managed separately for HCA.In the wet season, three clusters, namely A, B, and C, are observed in a dendrogram (Figure 5a).Cluster A consists of the water samples with lower concentration of Na + , and Cl -.It includes one dug well site (HW5) with all river sites and indicates minor influence of sewage discharge in wet season river water.The clusters B and C consist entirely of groundwater locations (Figure 5a) with higher concentration of Na + and Cl -and different isotopic compositions compared with that of river water locations.
Two major clusters, D and E, are categorized from a dendrogram of the dry season (Figure 5b).Cluster D consists only of river sites (eight sites), whereas Cluster E contains a combination of river and groundwater sites, indicating similarity in selected parameters [56,67].Cluster E has two subgroups, E1 and E2.Six dug well sites and one river site (HR1) are contained in Cluster E1.Similarly, four dug well sites and one river site (HR2) are grouped in Cluster E2.These two river samples, which are clustered with groundwater in Cluster E, are located at the uppermost section of the river (Figure 1), having lower concentrations of Na + and Cl − and classified as Ca-HCO3 type.However, in the case of the river sites of Cluster D, they have the highest concentration of Na + and Cl -, and it is categorized as aNa-K-HCO3 water type, which indicates pronounced effects of sewage discharge.Continued discharge of sewage in river water can develop an organic clogging layer on the riverbed [2], which will drastically reduce water exchange from the river channel to groundwater [3].The absence of groundwater samples grouping with river samples in Cluster D may indicate the possibility of clogging layer formation within the river bed.However, the clusters combining river water and groundwater, Cluster A from the wet season (Figure 5a) and Cluster E from the dry season (Figure 5b), show similarity in selected parameters and indicate interconnection between river and groundwater [9,11,37].
However, the current research has a limitation in the methodological approach where the river channel stage could not be developed comparing relative river water and groundwater level.The HCA along with hydraulic head between river water and groundwater can give a clearer view of interconnection between river and groundwater.

Identifying Areas of River Water and Groundwater Interconnection
Results from HCA imply interconnection between river water and groundwater in both wet and dry seasons.Cluster A from the wet season combines all river water with one groundwater (HW5), showing the similar water type Ca-HCO3.HW5 has shallow well depth (5 m) with water level depth at 1 m.It is located in permeable lower terrace deposits (Figure 1), generally used for cultivation.Many flooding events were recorded for the Hanumante River [68], with a major one recently in July 2018 [69] which inundated the entire area from Jagati to Madhyapur Thimi, along with HW2, HW5, and HW6 (Figure 1).However, HW5 is the only groundwater sample that was grouped with all river sites in Cluster A (Figure 5); HR5 is closest site to HW5 (80 m away; Figure 1).Identical isotopic compositions of δ 18 O (−8.0‰) and δD (−54.7‰)observed in HW5 and HR5 suggest possible recharge of groundwater in HW5.Similar concentration of Na + and Cl − with HR5, along with drastic dilution of these concentrations (Figure 2) at HW5, supports bank infiltration [4,18,19,70] as the mechanism recharging HW5 during the wet season (August of 2017).Similar types of results are presented by previous research conducted in developing countries for local scale [36] as well as for regional scale [37], indicating river water can recharge nearby groundwater.
In the dry season, dug wells HW1, HW2, HW4, HW6, HW7, and HW9 are grouped with the river site HR1 in Cluster E1, whereas dug wells HW3, HW5, HW8, and HW10 grouped with HR2 are clustered in E2.However, Cluster E from the dry season (Figure 5b) is suggestive of possible interaction between river water (HR1 and HR2) and groundwater.As HR1 and HR2 are located at the uppermost section of the river channel compared with grouped groundwater sites (Figure 1), there is a higher possibility of downstream groundwater recharge by HR1 and HR2 in the dry season.But differences in isotopic compositions of river water and nearby groundwater samples (Figure 4b) imply the possibility of other recharge sources besides river water.Further, groundwater samples exhibited similar trends of isotopic composition as that of the GMWL, suggesting dry season rainfall as one of the alternative sources of recharge for these dug wells.
Dug wells HW2, HW3, HW9, and HW10 are located within 20 m of HR2, HR3, HR9, and HR10, respectively (Figure 1).However, no single cluster formed includes this entire sampling site.These results indicate that interaction between river water and groundwater not only depends on distance from the river channel to groundwater sites, it also depends on well depth, water level depth, River water Groundwater topography of well location, water level in river channel, and sedimentological formation of the areas.

River Water Contribution to Groundwater
Dug wells recharged by river water have a mixture of water, including river water (RW) and original groundwater (GW).During the study period, only HW5 showed interconnection with river water during the wet season.The proportion of river water to groundwater can be estimated using a mass balance approach equation [18,36]: f = [(CS − CGW)/(CRW − CGW)] ×100% where CS is the Cl − , δD, or δ 18 O (mg/L or ‰) of mixed water (HW5 in the wet season); CRW is the Cl − , δD, or δ 18 O of river water (HR5 in the wet season), and CGW is the Cl − , δD,or δ 18 O of original groundwater (HW5 in the dry season) which has not been influenced by river water recharge.
The calculated fractional contribution is nearly 100% for Cl − , δD, or δ 18 O, indicating the contribution of river water (HR5) to recharge HW5 during the wet season (August 2017).This is the only well which had an identical isotopic value as in HR5 with much diluted water during the wet season.
The existence of interconnection between the Hanumante River and shallow groundwater shows that the surrounding shallow groundwater is contaminated by polluted river water.As river water deterioration continues with increasing urbanization, groundwater contamination is expected.It is very important to maintain river water quality for the improvement of peripheral shallow groundwater quality.

Conclusions
This study analyzes chemical and isotopic compositions of river water and groundwater to investigate interconnectivity between river water and shallow aquifers.Hydro-chemical parameters of river water exhibit significant temporal variations compared with groundwater.Groundwater and river water reveal a Ca-HCO3 type with similar trends of isotopic composition as that of the GMWL, confirming freshwater as a major source for groundwater recharge and river discharge during the wet season.Water types shifted from Ca-HCO3 to Na-K-HCO3, and strong positive correlations between Na + , K + , NH4 + -N, Ca 2+ , Mg 2+ , Cl − , HCO3 − , PO4-P, and SO4 2− indicate anthropogenic activities as a major source of contamination in dry season river water.Slight deviations of isotopic compositions from GMWL with a lower slope may suggest the possibility of evaporation in river water during the dry season.
Clusters formed by combinations of river water and groundwater, both in wet (Cluster A) and dry (Cluster E) seasons, indicate the presence of interconnectivity between river water and shallow groundwater.The wet season included one groundwater sample with all the river samples, including one well which was near the riverside (within 80 m).Identical isotopic composition and similar concentrations of Na + and Cl − in groundwater and river water suggest almost 100% of water recharge through bank infiltration comes from river to shallow groundwater during the wet season.Grouping of two upstream river samples with all downstream shallow groundwater implies that upstream river water is one of the sources of recharge for downstream shallow groundwater during the dry season.
This research concluded that polluted river water is one of the contamination sources for shallow groundwater around riverside areas.Such types of studies can be applied in other highly contaminated rivers to give a clear view about the contamination source for river peripheral groundwater.Higher extents of similar studies in rivers of the most urbanized areas are very important for groundwater management of the river peripheral area.

Figure 2 .
Figure 2. Bar diagram showing temporal and spatial variation of chemical parameters in river water (a and b) and groundwater (c and d).HR = Hanumante river water and HW = Hanumante groundwater.

Figure 3 .
Figure 3. Piper diagram of major ions of groundwater and river water in (a) wet season (b) dry season.HR = Hanumante river water, HW = Hanumante groundwater.

Figure 5 .
Figure 5. Dendrogram based on hierarchical clustering (Ward's method) in (a) wet season and (b) dry season.(c) Relation of Hanumante river water (HR) and Hanumante groundwater (HW) with cluster number in dendrogram.

Table 1 .
In situ measured parameters of groundwater and river water in the wet (7 August 2017) and dry (18 February 2018) seasons.Hanumante well water, HR = Hanumante river water, WLD = water level depth.EC = electrical conductivity, DO = dissolved oxygen.

Table 2 .
Correlation matrix of different chemical parameters.