This study reported the results of groundwater quality for drinking and irrigation purposes in Kurukshetra district, Haryana, India. Twelve parameters, namely pH, TDS, TH, sodium (Na+), magnesium (Mg2+), calcium (Ca2+), potassium (K+), chloride (Cl), sulfate (SO42−), bicarbonate (HCO3), nitrate (NO3), and fluoride (F) were evaluated and the Water Quality Index (WQI) was calculated. Based on the WQI, 31.57% and 68.43% of samples fall in excellent and good drinking water quality. The sodium adsorption ratio (SAR: 5.75–33.8), magnesium hazard ratio (MHR: 0.21–0.75), percent sodium (%Na: 31.34–72.84), residual sodium carbonate (RSC: 139–770), potential salinity (PS: 18.5–90), and Kelley Ratio (KR: 0.9–3.13) were calculated. The SAR, %Na, MHR, and KR indicated that most of the groundwater is suitable for irrigation. In addition, the inverse distance weighted (IDW) interpolation method was also employed to determine the spatial distribution of groundwater quality in the form of maps using ArcGIS software. The Pearson correlation matrix has been calculated to analyze the relationship between the water quality parameters. Gibbs plots have shown that evaporation and rock weathering are primary processes responsible for affecting the hydrochemistry of groundwater. Piper plots suggested that maximum groundwater samples are (HCO3–Na+), (Na+), and (HCO3) types.

  • The groundwater quality for drinking was compared with the BIS and WHO by evaluating the WQI in the Kurukshetra district of Haryana, India.

  • To analyze groundwater for irrigation, SAR, MHR, %Na, RSC, PS, and KR were calculated.

  • Piper trilinear, Gibbs diagrams, and spatial distribution maps were plotted to monitor groundwater's hydrochemistry.

  • The groundwater of the Kurukshetra can be used for drinking and irrigation after some primary treatments.

Graphical Abstract

Graphical Abstract
Graphical Abstract

Groundwater is being used as a primary source for industrial, irrigation, and drinking purposes in many parts of the world. Worldwide, 15% of groundwater is used for industrial work, 20% for irrigation purposes, and 65% for drinking purposes. About one-third population of the world depends upon groundwater for drinking purposes. Due to the excessive use of fertilizers and chemicals in agriculture, groundwater quality is being degraded continuously leading to an impact on human health (Rango et al. 2012; Zhang et al. 2020). In the past, several studies have been conducted around the world to assess the quality of groundwater and its suitability for drinking and irrigation purposes (Chung et al. 2015; Al-Barakah et al. 2017; Raza et al. 2017; Sheikh et al. 2017; Wu et al. 2018; Chaudhry & Sachdeva 2020; Adimalla 2021). Groundwater may contain harmful chemicals that affect human health seriously when consumed (Appelo & Postma 2004). A proper understanding of the science of groundwater geochemistry which uncovers the processes that influence the chemical configuration of groundwater is essential to monitor the water quality for drinking and irrigation purposes (Dassargues 2018). Therefore, it is strongly suggested to regularly evaluate groundwater quality to assess its fitness for irrigation and drinking purposes.

In the past few decades, groundwater contamination has been a severe problem worldwide. Rapid industrial growth, intensive urbanization, over usage of fertilizers in irrigation sectors, animal and human wastage, etc., are significant reasons for the deterioration of the groundwater quality. Therefore, several studies have been focused on groundwater quality studies for several usages (Amiri et al. 2015; Narsimha & Sudarshan 2018; Gharbi et al. 2019; Sood & Sharma 2020; Raheja et al. 2022; Wang et al. 2022). A study by Shan et al. (2021) evaluated the water quality using the Water Quality Index (WQI) in Bhindawas bird sanctuary, Haryana, India. Paul et al. (2019) studied groundwater quality in Tripura, Northeast India, using a Geographical Information System (GIS). Simsek & Gunduz (2007) studied the GIS-based irrigation water quality index (IWQI) using hydrochemical parameters like magnesium hazard ratio (MHR), infiltration hazard, electrical conductivity (EC), salinity hazard, sodium adsorption ratio (SAR), and others of the Simav Plain in Turkey. Gao et al. (2020) worked on the integrated weight WQI for evaluating groundwater quality in Xi'an, China. Karnena et al. (2022) investigated the suitability of groundwater using the WQI in the GIS environment covering a coastal region of Andhra Pradesh, India. They observed that groundwater quality varied from good to poor for drinking purposes and was free from fluoride content. Ram et al. (2021) studied groundwater quality by using the WQI under the GIS framework in Mahoba District, Uttar Pradesh, India, and found that the groundwater is safe for drinking purposes except at a few locations. The GIS-based inverse distance weighted (IDW) interpolation method is commonly used to obtain spatial distribution maps of groundwater quality (Şener et al. 2017; Balamurugan et al. 2020; Fang et al. 2020; Soujanya Kamble et al. 2020; Ram et al. 2021). It is a time-efficient and cost-effective technique to create various spatial distribution maps and projections. A GIS-based study is used for the spatial evaluation of various groundwater quality parameters. The literature review concluded that the assessment of groundwater quality is essential and offers basic information on groundwater for the sustainable development of any region. Most of the land in the Kurukshetra district of Haryana (state) is being used for agriculture. Further, this region is considered quite important for good-quality rice crops in the country. Therefore, it is very important to ascertain whether groundwater is suitable for irrigation or not. In addition to this, the groundwater of this study area is also frequently used for drinking purposes. As a result, the main aim of this study was to analyze the suitability of groundwater quality for drinking and irrigation purposes.

In the present study, multivariate statistical plots and spatial distribution maps developed in the GIS have been used to analyze physicochemical information of the groundwater of the Kurukshetra district of Haryana state, India. The groundwater quality was compared with guidelines laid down by the Bureau of Indian Standards (BIS 2015) and the World Health Organization (WHO 2011) by calculating the WQI as proposed in the literature. On the other hand, the suitability of groundwater for irrigation purposes was evaluated by computing indices such as SAR, MHR, percent sodium (%Na) or sodium hazard, residual sodium carbonate (RSC), potential salinity (PS), and Kelley Ratio (KR). Furthermore, Piper trilinear and Gibbs diagrams were also plotted for the classification of hydrochemistry of groundwater.

The study area district Kurukshetra is situated in the northeast part of Haryana, India. It is located between 29°53′00″–30° 15′02″ N latitude and 76°26′27″–77°07′57″ E longitude (Figure 1) and has a geographical area of 1,530 km2 with an average temperature of 23.9 °C. Kurukshetra district has been divided into six blocks: Pehowa, Thanesar, Shahabad, Ladwa, Ismailabad, and Babain, with a total population of 9,64,231 as per the census 2011 and the forecasted population of Kurukshetra District in the year 2022 is 1,036,329 (estimates as per Ministry of Electronics & Information Technology 2019). The study area falls within two river basins, i.e., the Upper Yamuna Basin and Upper Ghaggar Basin. The average annual rainfall is about 582 mm, out of which 81% is received during the southwest monsoon and 19% is received during the non-monsoon period. The entire study area covers tropical arid brown soils.
Figure 1

A map showing groundwater sampling locations of northeast part of Haryana, Kurukshetra, India.

Figure 1

A map showing groundwater sampling locations of northeast part of Haryana, Kurukshetra, India.

Close modal

For all agricultural and drinking purposes, groundwater is being utilized as the main source of water in the entire district. The groundwater level lies between 20.18 and 32.64 m below ground level (bgl) in pre-monsoon and 21.80 and 34.41 m bgl in post-monsoon (CGWB 2013) and is being depleted at a rate of more than 1 m/year due to excessive water use for rice crop. As per CGWB (2013), the groundwater level is in the dark zone and groundwater depletion is a significant problem in the area. Moreover, the average groundwater level depth is more than 30 m. Therefore, it is essential to study the groundwater quality of aquifers in the Kurukshetra district.

Sampling and analysis

Nineteen groundwater samples were collected from different locations in the Kurukshetra district, Haryana, during February 2021 (Figure 1). The samples were collected in 1-L high-density polyethylene (HDPE) plastic bottles and were used further for physicochemical analysis in the laboratory. To avoid external impurities, groundwater was flushed for 5–7 min before filling the HDPE bottles. The bottled groundwater samples were stored in an icebox and transported to the laboratory. These samples were analyzed for pH, total dissolved solids (TDS), EC, total hardness (TH), sodium (Na+), magnesium (Mg2+), calcium (Ca2+), potassium (K+), chloride (Cl), sulfate (SO42−), carbonate (CO32−), bicarbonate (HCO3), nitrate (NO3), and fluoride (F) within 48 h of collection. The pH, TDS, and EC were measured in the field by using a pH/TDS/EC meter (Hanna HI 9811-5). While TH, Ca2+, and Mg2+ were determined by the titration method (EDTA). The Cl ion was estimated using the titration method with silver nitrate. The CO32− and HCO3 ions were estimated by titrating with sulfuric acid. The SO42−, NO3, and F ions were analyzed using a spectrophotometer (APHA 2005). The soluble K+ and Na+ ions were evaluated by using Flame Photometer (EI-380). To plot spatial distribution maps of various chemical constituents, ArcGIS software with IDW interpolation was applied. An IDW interpolation method helps to calculate a value for any unknown location and gives more weight to points closest to the prediction locations in the study area. The statistical summary and Pearson correlation matrix of the evaluated hydrochemical parameters have been reported in Tables 2 and 5. The charge balance error (CBE) was used to assuring the accuracy of the chemical analysis of the groundwater sample and can be calculated using the following equation:
(1)

The observed CBE is within the suggested limit of ±5% (Zhang et al. 2018), where all the anions (Cl, SO42−, HCO3, NO3, and F) and cations (Na+, Mg2+, Ca2+, and K+) are in milligrams per liter (mg/L).

Water quality index

The Bureau of Indian Standards (BIS 2015) suggested the use of the WQI for drinking water quality standards. To calculate the WQI, the following steps are followed:

In the first step, 12 chemical parameters were assigned weights () as in Table 1 (Ramakrishnaiah et al. 2009).

Table 1

Relative weight of hydrochemical parameters

Sr. No.Name of parametersWeight ()Relative weight =Permissible limits (BIS 2015)
pH 0.091 6.5–8.5 
TDS 0.091 500 
TH 0.091 200 
Ca2+ 0.068 200 
Mg2+ 0.068 30 
Na+ 0.091 200 
K+ 0.023 12 
HCO3 0.023 120 
Cl 0.11 250 
10 SO42− 0.11 200 
11 NO3 0.11 45 
12 F 0.11 1.5 
   = 44   
Sr. No.Name of parametersWeight ()Relative weight =Permissible limits (BIS 2015)
pH 0.091 6.5–8.5 
TDS 0.091 500 
TH 0.091 200 
Ca2+ 0.068 200 
Mg2+ 0.068 30 
Na+ 0.091 200 
K+ 0.023 12 
HCO3 0.023 120 
Cl 0.11 250 
10 SO42− 0.11 200 
11 NO3 0.11 45 
12 F 0.11 1.5 
   = 44   
Table 2

Statistical summary of measured/calculated hydrochemical parameters of groundwater samples of Kurukshetra district and comparison with drinking water standards (WHO 2011)

Sr. No.Name of parameterMinimum valueMaximum valueMean valueWHO limits
Number of Samples
Desirable limitMaximum allowable limitWithin desirable limitOut of maximum allowable limit
pH 7.85 8.96 8.48 6.5–8.5 9.5 12 – 
TDS (mg/L) 403 1,400 713.53 500 1,500 – 
TH (mg/L) 177 276 235.80 200 600 – 
Ca2+ (mg/L) 10 45 28.52 75 200 19 – 
Mg2+ (mg/L) 14 45 25.79 50 150 19 – 
Na+ (mg/L) 34 187 83.74 – 200 – – 
K+ (mg/L) 3.4 14.5 7.81 – 12 – 
CO32− (mg/L) 23 86 44.63 – – – – 
HCO3 (mg/L) 155 745 294.74 – – – – 
10 NO3 (mg/L) 1.5 22.6 6.08 10 45 15 – 
11 Cl (mg/L) 12 75 37 250 1,000 19 – 
12 SO42− (mg/L) 10 65 37.21 200 400 19 – 
13 F (mg/L) 0.35 0.98 0.61 1.0 1.5 19 – 
Sr. No.Name of parameterMinimum valueMaximum valueMean valueWHO limits
Number of Samples
Desirable limitMaximum allowable limitWithin desirable limitOut of maximum allowable limit
pH 7.85 8.96 8.48 6.5–8.5 9.5 12 – 
TDS (mg/L) 403 1,400 713.53 500 1,500 – 
TH (mg/L) 177 276 235.80 200 600 – 
Ca2+ (mg/L) 10 45 28.52 75 200 19 – 
Mg2+ (mg/L) 14 45 25.79 50 150 19 – 
Na+ (mg/L) 34 187 83.74 – 200 – – 
K+ (mg/L) 3.4 14.5 7.81 – 12 – 
CO32− (mg/L) 23 86 44.63 – – – – 
HCO3 (mg/L) 155 745 294.74 – – – – 
10 NO3 (mg/L) 1.5 22.6 6.08 10 45 15 – 
11 Cl (mg/L) 12 75 37 250 1,000 19 – 
12 SO42− (mg/L) 10 65 37.21 200 400 19 – 
13 F (mg/L) 0.35 0.98 0.61 1.0 1.5 19 – 

The next step is to calculate relative weights as per the following equation.
(2)
where is the weight, is the relative weight, and n is the total number of parameters.
In the third step, the water quality rating has been estimated by using the following equation:
(3)
where is the quality rating, is the concentration of each chemical parameter, and is the limits of (BIS 2015) standards.
Finally, (water quality sub-index) and the WQI were calculated using the following equations:
(4)
(5)

Irrigation groundwater quality evaluation

RSC, SAR, MHR, percent sodium (%Na) or sodium hazard, PS, and KR were calculated for assessing the groundwater quality for irrigation purposes. To calculate these chemical indices, the following equations were used (Richards 1954; Wilcox 1955; Kelley 1963; Doneen 1964; Ragunath 1987; Gharbi et al. 2019; Gugulothu et al. 2022).

Sodium adsorption ratio

It is an essential indicator for evaluating groundwater quality for irrigation purposes. If groundwater has a low Ca2+ ion concentration and a high Na+ ion concentration, the soil structure may be damaged due to the dispersion of clay particles (Richards 1954; Subba Rao 2006). It can be calculated by using the following Equation (6).
(6)

Magnesium hazard ratio

The MHR is used to estimate groundwater quality for irrigation purposes. The ratio of magnesium cation to the sum of cations (Ca2+ + Mg2+) is known as the MHR and is defined by Equation (7).
(7)

Percent sodium (%Na)

Wilcox (1955) introduced the %Na ratio to measure irrigation groundwater quality. It is the ratio of Na+ ion to the sum of all cations (Na+ + K+ + Ca2+ + Mg2+) and it is calculated by Equation (8).
(8)

Residual sodium carbonate

It is estimated by using Equation (9) as under:
(9)

If the RSC value <1.25, then water is useful for irrigation; if the RSC lies between 1.25 and 2.25, then the condition is doubtful, if the value exceeds 2.25, then groundwater is unfit for irrigation.

Potential salinity

Doneen (1964) proposed the use of PS for groundwater quality for irrigation purposes and defined it by the following Equation (10).
(10)

If the PS value <3, then it can be used for irrigation; otherwise, it remains unsuitable.

Kelley ratio

The ratio of Na+ ion to the sum of cations (Ca2+ + Mg2+) is known as KR and can be calculated by using Equation (11).
(11)

If the KR value <1, then groundwater can be used for irrigation. If the KR value >1, groundwater is not used for irrigation.

Hydrogeochemistry

A summary of 19 groundwater samples in terms of various parameters is provided in Table 2. The pH or hydrogen ion plays an important role in assessing groundwater quality. The value of pH ion value ranges from 7.85 to 8.96 in the study areas with a mean of 8.48, signifying a mildly alkaline nature of groundwater (Table 2). The TDS lies between 403 and 1,400 mg/L with an average of 713.53 mg/L. If the value of TDS is below <500 mg/L, then it is considered potable and freshwater (BIS 2015). A higher value of TDS might contribute to bad taste in groundwater (Monteiro De Oliveira et al. 2021). This classification indicates that 86% of groundwater samples are suitable for drinking purposes (<1,000 mg/L). Davis & De (1966) proposed to divide TH into four different groups. If TH is below 75 mg/L, water will be considered slightly hard; moderately hard if it is between 75 and 150 mg/L; hard if the value is between 150 and 300 mg/L, and very hard if its value exceeds 300 mg/L. The TH concentration variation suggests that 100% of the study area falls under 150–300 mg/L (hard water) as per these classifications.

Cation chemistry

In the cation compositions of groundwater (Ca2+, Mg2+, Na+, and K+), calcium (Ca2+) ions range from 10 to 45 mg/L with a mean value equal to 28.52 mg/L. Due to agricultural drainage, a higher amount of Ca2+ was observed in Baronda village of Kurukshetra district. The desirable and allowable limit of Ca2+ for drinking purposes is 75 and 200 mg/L. In the case of magnesium (Mg2+) ion, the desirable and allowable limits as per drinking purposes are 50 and 150 mg/L. Analysis of the results (Table 2) suggested that 100% of groundwater samples in the study area fall under desirable limits. The concentration of Na+ ion lies between 34 and 187 mg/L, indicating that 100% of groundwater samples lie within maximum allowable limits (WHO 2011). Moreover, a higher amount of Na+ intake will cause hypertension and osteoporosis problems. The level of potassium (K+) ion is relatively high in the Kurukshetra district as its value lies between 3.4 and 14.5 mg/L with a mean value of 7.81 mg/L. A few groundwater samples, i.e., 16%, falls under the maximum allowable limits of K+ ion (Table 2).

Anion chemistry

Chloride (Cl) is considered the most essential anion present in groundwater for irrigation. As per WHO (2011), a maximum allowable limit of Cl is 1,000 mg/L as suggested for drinking purposes. A high intake of Cl anion in drinking water can cause heart and kidney problems (Patel et al. 2022). The Cl content of the Kurukshetra district ranged from 12 to 75 mg/L with an average of 37 mg/L. A higher level of Cl ion was found in the Samalkhi village of the study area. Sulfate (SO42−) is an important anion occurring in groundwater. It may enter groundwater through the seepage of industrial effluents (APHA 2005). The SO42− ion in the study area ranges from 10 to 65 mg/L with a mean of 37.21 mg/L (Table 2). If the concentration of SO42− lies between 300 and 400 mg/L, it causes a bitter taste in groundwater, and if it is more than 1,000 mg/L, it can cause stomach disorders (Patel et al. 2022). However, a high intake of SO42− ion in water may cause probable dehydration and laxative effects (Li et al. 2019). The carbonate (CO32−) ion range varies from 23 to 86 mg/L with an average of 44.63 mg/L. The value of bicarbonate (HCO3) ion lies in a range from 155 to 745 mg/L. Neither BIS nor WHO has recommended any limits (both desirable and allowable) for drinking purposes for CO32− and HCO3 ions. The F anion values ranged from 0.35 to 0.98 mg/L with a mean of 0.61 mg/L. The F anion is categorized into three different classifications. If the F value is below 1 mg/L, it is most suitable. If the range is between 1 and 1.5 mg/L, it is not fit for drinking (WHO 2011). Continuous consumption of higher F concentration may cause teeth problems and skeletal fluorosis. The nitrate (NO3) is formed by domestic and industrial discharge, biochemical, and fertilizer plants. The NO32− anion values ranged from 1.5 to 22.6 mg/L with a mean of 6.08 mg/L. As per BIS (2015), the maximum allowable limit of NO3 is 45 mg/L, and all the groundwater samples were are falling within the allowable limit.

Classification of groundwater quality for drinking purpose

Table 3 provides the WQI range, classifications, and percentage of all the groundwater samples tested. The calculated WQI values lie between 46.75 and 77.31, with an average of 55.33. The WQI over this study area is divided into five different categories: excellent water quality, if the WQI is <50; good water quality, if the WQI lies between 50 and 100; poor water quality, if the WQI is between 100 and 200; inferior water quality, if the WQI lies between 200 and 300, and unsuitable for the drinking purposes if the WQI values >300 as suggested by Ramakrishnaiah et al. (2009). It can be noted from Table 3 that 68.43% of groundwater samples are found to be good for drinking purposes, and 31.57% of groundwater samples fall in the excellent water quality category. The spatial distribution map of the WQI showed that good drinking water quality was distributed in most of the study area (Figure 2). It is noted from Figure 2 that the eastern part of the study area has excellent groundwater quality. The maximum value of the WQI was observed in Thana village and the minimum value in Berthala village of the study area (Table 3).
Table 3

Classification of the WQI, percentage of groundwater, and water type in the Kurukshetra district

Sample No.Sampling locationLatitudeLongitudeWQI valueWQI rangeType of groundwater% of samples
S1 Baronda 29.9556 77.0522 48.25 < 50 Excellent 31.57 
S2 Berthala 30.1292 76.9833 46.75 
S3 Ishaq 30.0167 76.4917 49.17 
S4 Ram Sharan Majra 30.0017 77.0003 49.92 
S5 Shahbad 30.1667 76.8667 48.41 
S6 Dabkhera 30.0006 76.7536 49.6 
S7 Jogna Khera 29.9736 76.8278 50.36 50–100 Good 68.43 
S8 Kaulapur 30.025 76.8958 60.02 
S9 Pehowa 29.9792 76.5833 53.53 
S10 Hatira 29.95 76.7625 52.3 
S11 Bateri 29.9208 76.6083 51.15 
S12 Samalkhi 30.235 76.8278 58.19 
S13 Sirsala 30.0003 76.5019 59.95 
S14 Jandaula 29.8889 76.5728 60.48 
S15 Santokhpur 30.0292 76.725 61.46 
S16 Bachki 30.0394 76.7556 60.67 
S17 Bodhni 30.0542 76.5583 60.57 
S18 Malikpur-Singhpura 30.00 76.7514 53.34 
S19 Thana 29.9083 76.5083 77.31 
 – – – – 100–200 Poor 
 – –  – 200–300 Very poor 
 – – – – >300 Not suitable 
Sample No.Sampling locationLatitudeLongitudeWQI valueWQI rangeType of groundwater% of samples
S1 Baronda 29.9556 77.0522 48.25 < 50 Excellent 31.57 
S2 Berthala 30.1292 76.9833 46.75 
S3 Ishaq 30.0167 76.4917 49.17 
S4 Ram Sharan Majra 30.0017 77.0003 49.92 
S5 Shahbad 30.1667 76.8667 48.41 
S6 Dabkhera 30.0006 76.7536 49.6 
S7 Jogna Khera 29.9736 76.8278 50.36 50–100 Good 68.43 
S8 Kaulapur 30.025 76.8958 60.02 
S9 Pehowa 29.9792 76.5833 53.53 
S10 Hatira 29.95 76.7625 52.3 
S11 Bateri 29.9208 76.6083 51.15 
S12 Samalkhi 30.235 76.8278 58.19 
S13 Sirsala 30.0003 76.5019 59.95 
S14 Jandaula 29.8889 76.5728 60.48 
S15 Santokhpur 30.0292 76.725 61.46 
S16 Bachki 30.0394 76.7556 60.67 
S17 Bodhni 30.0542 76.5583 60.57 
S18 Malikpur-Singhpura 30.00 76.7514 53.34 
S19 Thana 29.9083 76.5083 77.31 
 – – – – 100–200 Poor 
 – –  – 200–300 Very poor 
 – – – – >300 Not suitable 
Figure 2

A spatial distribution map of the WQI of groundwater of Kurukshetra district.

Figure 2

A spatial distribution map of the WQI of groundwater of Kurukshetra district.

Close modal

Classification of groundwater quality for irrigation purposes

Agriculture is the prime occupation in the study area, with extensive rice, wheat, and sugarcane farming. Keeping in view of extensive use of groundwater for water-intensive crops like rice and sugarcane, the estimation of groundwater quality is necessary. To assess the groundwater quality for irrigation uses, chemical indices like SAR, MHR, %Na, PS, and RSC were used in this study.

The spatial distribution map of SAR of Kurukshetra district is shown in Figure 3. The spatial distribution pattern of SAR suggested that the central portion along NW–SW in the district shows the presence of medium-to-high sodium in groundwater. The SAR value ranges from 5.75 to 33.86, with a mean value equal to 16.44 (Table 4). The highest value of SAR was recorded in Thana village, and the lowest value was recorded in Baronda village of Kurukshetra district. The EC value was recorded between 458 and 1,250 μS/cm with a mean of 666.26 μS/cm (Table 4). In present study, six samples (S1, S2, S5, S7, S10, S12) fall under excellent class (low sodium), i.e. SAR values <10. The seven samples (S4, S6, S8, S9, S11, S17, S18) fall in good class (medium sodium) and five water samples (S3, S13, S14, S15, and S16) have high sodium content. The samples that are graded as low sodium and medium sodium are used for irrigation. Based on the concentration of Na+ ion reacts with soil, which leads to an increase and decrease in its impermeability. Karanth (1987) stated that a higher %Na value could be a source of deflocculation, causing damage to the tilth and soil permeability. The %Na concentration in this study area lies between 31.33 and 72.83%, with a mean value of 55.36%. The lowest value of %Na was recorded in Baronda village and the highest was recorded in Bachki village in the Pehowa block of the study area (Figure 4). The %Na map clearly shows that it is maximum in the central part and SW portion of the study area. It can be seen from Table 5 that three samples (S3, S16, and S19) exceeded the limit of 60. A total of 15 groundwater samples have %Na values <60, which can be graded as good quality for irrigation use.
Table 4

Classification of irrigation water based on SAR and EC

Name of parameterWater typeRangeNumber of samples% of samples
SAR Low sodium <10 31.58 
Medium sodium 10–18 36.84 
High sodium 18–26 26.32 
Very high sodium >26 5.26 
EC Low saline <250 
Medium saline 250–750 16 84.21 
High saline 750–2,250 15.79 
Very high saline >2,250 
Name of parameterWater typeRangeNumber of samples% of samples
SAR Low sodium <10 31.58 
Medium sodium 10–18 36.84 
High sodium 18–26 26.32 
Very high sodium >26 5.26 
EC Low saline <250 
Medium saline 250–750 16 84.21 
High saline 750–2,250 15.79 
Very high saline >2,250 
Figure 3

SAR distribution.

Figure 3

SAR distribution.

Close modal
Figure 4

%Na distribution.

Figure 4

%Na distribution.

Close modal
Ragunath (1987) proposed that the use of RSC may be used for groundwater quality for irrigation purposes. The RSC values greater than 2.25 suggested the unsuitability of groundwater for irrigation purposes as mentioned in Table 5. The spatial distribution map of the RSC of the Kurukshetra district is shown in Figure 5. The lowest value of RSC was recorded in Pehowa block, and the highest value was recorded in Thana village of Kurukshetra district. Tiwari & Manzoor (1988) have suggested that negative values of RSC show that Ca2+ and Mg2+ concentrations are higher in groundwater and can further deteriorate the quality of soil, and will decrease the yield of the crop. Moreover, Mg2+ concentration has more severe effects than Ca2+ concentration on the crop. Therefore, the MHR parameter is another significant ratio for estimating irrigation groundwater quality (Szabolcs & Darab 1964; Paliwal 1972; Ragunath 1987). In this study, MHR values are ranging from 29.17 to 75 mg/L with a mean value of 47.77 mg/L. it is evident from Table 5 that most of the groundwater samples (62.63%) are suitable for irrigation, and 37.37% (S7, S8, S14, S16, S17, and S19) of samples are not suitable for irrigation, which are having MHR >50. The maximum value of MHR was observed in Dab Khera village of Ladwa block, and the minimum value was in Jandaula village of Thanesar. The spatial distribution map showing the variation in MHR values of groundwater is shown in Figure 6. The PS is estimated using Equation (10) and values range from 18.5 to 90 mg/L with an average value equal to 55.61 mg/L. The spatial distribution map of the PS of groundwater was prepared as shown in Figure 7. The pattern of PS revealed that groundwater is not suitable for irrigation. Kelley ratio/index (KR) is also applied in the study area (Kelley 1963). According to the KR, 21.05% (S3, S15, S16, and S19) of samples are not suitable for irrigation purposes, whereas 78.95% of groundwater samples are having KR is <1 (Table 5). The minimum value of KR is observed from Baronda village and the maximum value is observed in Ishaq village in Pehowa block in Kurukshetra district (Figure 8).
Figure 5

RSC distribution.

Figure 5

RSC distribution.

Close modal
Figure 6

MHR distribution.

Figure 6

MHR distribution.

Close modal
Figure 7

PS distribution.

Figure 7

PS distribution.

Close modal
Figure 8

KR distribution.

Figure 8

KR distribution.

Close modal
The cation exchange process recognizes Na2+ and Cl concentrations (Sethy et al. 2016) in water samples. Figure 9(a) shows scatter plots of groundwater concentration of Cl versus Na2+ + Mg2+ in mg/L. It is clear from Figure 9(a) that all the groundwater samples of the Kurukshetra district are laying below the best-fit line. It means that ion exchange and reverse ion exchange processes (Sethy et al. 2016) are taking place. However, Figure 9(b) represents scatter plots of groundwater concentration of Ca2+ + Mg2+ versus HCO3 + CO32− + SO42− in mg/L. Mostly, groundwater lies above the best-fit line representing the dissolution of calcite, gypsum, and dolomite in the groundwater (Kammoun et al. 2022).
Figure 9

(a) Scatter plots of Cl versus Na2++ Mg2+ of groundwater concentration; and (b) scatter plots of HCO3 + CO32− + SO42− verses Ca2+ + Mg2+ of groundwater concentration.

Figure 9

(a) Scatter plots of Cl versus Na2++ Mg2+ of groundwater concentration; and (b) scatter plots of HCO3 + CO32− + SO42− verses Ca2+ + Mg2+ of groundwater concentration.

Close modal

Statistical analysis of hydrochemical parameters

Pearson, Kendall, and Spearman are the three essential methods for correlation. The Kendall method is used to measure either nonlinear or linear correlation, the Pearson method measures linear correlation, and Spearman measures nonlinear correlation. The Pearson correlation matrix has been used in this study to establish a relation between groundwater quality parameters. The values of the Pearson correlation coefficient range from −1 to +1. When the value is 0, there is no correlation between the two parameters. When the coefficient is +1, there is a perfect positive linear correlation. When the coefficient is −1, there is a perfect negative linear correlation. Overall, when the correlated value lies between 0.50 and 0.70, it specifies a moderate correlation. When the coefficient value >0.7, it indicates a strongly correlated (Giridharan et al. 2008). Table 6 presents the Pearson correlation matrix of groundwater concentration of 13 hydrochemical parameters. It can be observed from Table 6 that TDS has a strong positive correlation with CO32−, Na+, and HCO3 ions, whereas Na+ is well correlated with CO32− and HCO3 (Table 6).

Table 5

Irrigation chemical indices and percentage of groundwater samples

ParametersValuesWater typeNumber of samples% of samples
%Na (Wilcox 1955<20 Excellent 
20–40 Good 10.53 
40–60 Permissible 13 67.90 
60–80 Doubtful 21.57 
>80 Unsuitable 
PS (Doneen 1964<3 Suitable 
>3 Unsuitable 19 100 
MHR (Ragunath 1987<50 Suitable 12 62.63 
>50 Unsuitable 37.37 
RSC (Ragunath 1987<1.25 Good 
1.25–2.25 Doubtful 
>2.25 Unsuitable 19 100 
KR (Kelley 1963<1 Suitable 15 78.95 
>1 Unsuitable 21.05 
ParametersValuesWater typeNumber of samples% of samples
%Na (Wilcox 1955<20 Excellent 
20–40 Good 10.53 
40–60 Permissible 13 67.90 
60–80 Doubtful 21.57 
>80 Unsuitable 
PS (Doneen 1964<3 Suitable 
>3 Unsuitable 19 100 
MHR (Ragunath 1987<50 Suitable 12 62.63 
>50 Unsuitable 37.37 
RSC (Ragunath 1987<1.25 Good 
1.25–2.25 Doubtful 
>2.25 Unsuitable 19 100 
KR (Kelley 1963<1 Suitable 15 78.95 
>1 Unsuitable 21.05 
Table 6

Pearson correlation matrix of groundwater samples of major ions

 
 

Hydrochemical classification of groundwater

(a) Piper trilinear diagram

Piper diagram has been plotted to understand groundwater geochemistry by plotting the cations and anions (Piper 1944). This graph normally helps in revealing chemical relationships among groundwater in more concrete terms (Walton 1970). Water samples are divided into four basic types of diamonds. The anions (Ca2+, Mg2+, Na+, and K+) are expressed as a ternary plot in the lower rightward, while cations (HCO3, SO42−, Cl, and CO32−) are expressed as a ternary plot in the lower leftward. Both anions and cations are expressed by diamond plots in the middle. It can be inferred from Figure 10 that the groundwater described is rich in Ca2+ + Mg2+ and Cl−+SO42−, resulting in permanent hardness. It also shows that the maximum groundwater samples are HCO3 -Na+, Na+, and HCO3 type.
Figure 10

A Piper diagram presenting the geochemistry of groundwater samples in the Kurukshetra district.

Figure 10

A Piper diagram presenting the geochemistry of groundwater samples in the Kurukshetra district.

Close modal

(b) Gibbs diagram

In order to understand hydrochemical faces and aquifer lithology (Gibbs 1970; Pei-yue et al. 2010; Zhang et al. 2020), the Gibbs diagram is plotted (Figure 11). The ratios of cations [(Na+ + K+)/(Na+ + K+ + Ca2+)] or anions [Cl/(Cl + HCO3)] are taken on the horizontal plane and TDS are plotted on the vertical plane. There are three main natural mechanisms forming the water chemistry in these diagrams, i.e., precipitation dominance, rock dominance, and evaporation dominance. As shown in Figures 11(a) and 11(b), the groundwater samples were mainly in the zone of rock dominance and also toward pointing evaporation-dominant zone. It means that leaching and rock weathering are the major mechanisms controlling groundwater chemistry in the study area, which may be related to the local semi-arid climate with little rainfall and higher evaporation.
Figure 11

Gibbs diagrams: (a) and (b) present the features of rock-dominant.

Figure 11

Gibbs diagrams: (a) and (b) present the features of rock-dominant.

Close modal

This work is based on the hydrochemistry of significant anions and cations of groundwater of the Kurukshetra district of Haryana state in India. The main objective is to better understand groundwater chemical composition to estimate its suitability for drinking and irrigation uses. The conclusions drawn from this work are as follows:

  • The groundwater of Kurukshetra district, Haryana state, is found to be mildly alkaline and hard. The ions, namely Ca2+, Mg2+, Cl, SO42−, and F ions are within desirable limits (WHO 2011). However, HCO3 anion and Na+ cation are most dominant in the groundwater. It is observed that the K+ ion is relatively higher and is laying between 3.4 and 14.5 mg/L. A total of 16% of groundwater samples exceeded limits for K+ ions as per the recommendation of WHO (2011).

  • The calculated values of WQI lie between 46.75 and 77.31. The value of the WQI for drinking purposes revealed that 31.57% of groundwater samples belong to excellent water quality and 68.43% of groundwater samples are in a good category.

  • On the basis of values of groundwater salinity, SAR, RSC, MHR, sodium hazard, PS, and KR, it can be concluded that 15.79, 31.58, 100, 37.37, 21.57, 100, and 21.05% of the study area are not suitable for irrigation uses.

  • Piper plots exhibited that maximum groundwater samples are HCO3 -Na+, Na+, and HCO3 types. Gibbs plots revealed that evaporation and rock weathering are the major techniques influencing the hydrochemistry of groundwater.

In India, approximately 85% of groundwater is used for drinking and irrigation purposes and is the biggest user of groundwater around the Globe. This is a severe insinuation to the sustainability of living beings, irrigation, and commercial growth. The approach used in this work can be efficiently helpful to other areas having similar conditions and is beneficial for groundwater management authorities in ensuring safe water for the stakeholders.

H.R. is grateful to the Ministry of Education, Government of India, for providing the fellowship to do this research. The author conveys his gratitude to the National Institute of Technology, Kurukshetra, for giving research facilities in completing this work.

All relevant data are included in the paper or its Supplementary Information.

The authors declare there is no conflict.

Al-Barakah
F. N.
,
Al-jassas
A. M.
&
Aly
A. A.
2017
Water quality assessment and hydrochemical characterization of Zamzam groundwater, Saudi Arabia
.
Applied Water Science
7
(
7
),
3985
3996
.
https://doi.org/10.1007/s13201-017-0549-x
.
Amiri
V.
,
Sohrabi
N.
&
Dadgar
M. A.
2015
Evaluation of groundwater chemistry and its suitability for drinking and agricultural uses in the Lenjanat plain, central Iran
.
Environmental Earth Sciences
74
(
7
),
6163
6176
.
https://doi.org/10.1007/s12665-015-4638-6
.
APHA
2005
Standard Methods for the Examination of Water and Waste Water
, 21st Edn.
American Public Health Association
,
Washington, DC
.
Appelo
C. A. J.
,
Postma
D.
,
2004
Geochemistry, Groundwater and Pollution
(
Appelo
C. A. J.
&
Postma
D.
eds.).
CRC Press
.
https://doi.org/10.1201/9781439833544.
Balamurugan
P.
,
Kumar
P. S.
,
Shankar
K.
,
Nagavinothini
R.
&
Vijayasurya
K.
2020
Non-carcinogenic risk assessment of groundwater in southern part of Salem district in Tamilnadu, India
.
Journal of the Chilean Chemical Society
65
(
1
),
4697
4707
.
https://doi.org/10.4067/S0717-97072020000104697
.
BIS
2015
Indian Standard Drinking Water–Specification (Second Revision)
.
Bureau of Indian Standards (BIS)
,
IS 10500, New Delhi
, pp.
2
6
.
CGWB
2013
Ground Water Information Booklet of Kurukshetra District, Haryana
.
Central Ground Water Board, Ministry of Water Resources Government of India, North Western Region
,
Chandigarh
.
Chaudhry
A. K.
&
Sachdeva
P.
2020
Groundwater quality and non-carcinogenic health risk assessment of nitrate in the semi-arid region of Punjab, India
.
Journal of Water and Health
18
(
6
),
1073
1083
.
https://doi.org/10.2166/wh.2020.121
.
Chung
S. Y.
,
Venkatramanan
S.
,
Kim
T. H.
,
Kim
D. S.
&
Ramkumar
T.
2015
Influence of hydrogeochemical processes and assessment of suitability for groundwater uses in Busan City, Korea
.
Environment, Development and Sustainability
17
(
3
),
423
441
.
https://doi.org/10.1007/s10668-014-9552-7
.
Dassargues
A.
2018
Hydrogeology
.
CRC Press
.
https://doi.org/10.1201/9780429470660.
Davis
S. N.
&
De Wiest
R. J. M.
1966
Hydrogeology
.
Wiley
,
New York
.
Doneen
L. D.
1964
Notes on Water Quality in Agriculture
.
Water Science and Engineering Paper 4001, Department of Water Science and Engineering, University of California
,
Davis
.
Fang
Y.
,
Zheng
T.
,
Zheng
X.
,
Peng
H.
,
Wang
H.
,
Xin
J.
&
Zhang
B.
2020
Assessment of the hydrodynamics role for groundwater quality using an integration of GIS, water quality index and multivariate statistical techniques
.
Journal of Environmental Management
273
,
111185
.
https://doi.org/10.1016/j.jenvman.2020.111185
.
Gao
Y.
,
Qian
H.
,
Ren
W.
,
Wang
H.
,
Liu
F.
&
Yang
F.
2020
Hydrogeochemical characterization and quality assessment of groundwater based on integrated-weight water quality index in a concentrated urban area
.
Journal of Cleaner Production
260
,
121006
.
https://doi.org/10.1016/j.jclepro.2020.121006
.
Gibbs
R. J.
1970
Mechanisms controlling world water chemistry
.
Science, New Series
170
(
3962
),
1088
1090
.
Giridharan
L.
,
Venugopal
T.
&
Jayaprakash
M.
2008
Evaluation of the seasonal variation on the geochemical parameters and quality assessment of the groundwater in the proximity of River Cooum, Chennai, India
.
Environmental Monitoring and Assessment
143
(
1–3
),
161
178
.
https://doi.org/10.1007/s10661-007-9965-y
.
Gugulothu
S.
,
Subbarao
N.
,
Das
R.
&
Dhakate
R.
2022
Geochemical evaluation of groundwater and suitability of groundwater quality for irrigation purpose in an agricultural region of South India
.
Applied Water Science
12
(
6
),
1
13
.
https://doi.org/10.1007/s13201-022-01583-w
.
Kammoun
A.
,
Abidi
M.
&
Zairi
M.
2022
Hydrochemical characteristics and groundwater quality assessment for irrigation and drinking purposes: a case of Enfidha aquifer system, Tunisia
.
Environmental Earth Sciences
81
(
2
),
1
15
.
https://doi.org/10.1007/s12665-021-10163-1
.
Karanth
K. R.
1987
Role of Hydrogeology in Integrated Water Resource Development and Management
.
Environmental Centre, Groundwater News
, pp.
23
30
.
Karnena
M. K.
,
Konni
M.
,
Dwarapureddi
B. K.
,
Satyanarayana
Y.
&
Saritha
V.
2022
GIS-based approach qualitative features of sub-surface water from coastal district in Andhra Pradesh
.
Applied Water Science
12
(
3
),
1
15
.
https://doi.org/10.1007/s13201-021-01506-1
.
Kelley
W. P.
1963
Use of saline irrigation water
.
Soil Sciences
95
,
385
391
.
Li
P.
,
He
X.
&
Guo
W.
2019
Spatial groundwater quality and potential health risks due to nitrate ingestion through drinking water: a case study in Yan'an City on the Loess Plateau of northwest China
.
Human and Ecological Risk Assessment
25
(
1–2
),
11
31
.
https://doi.org/10.1080/10807039.2018.1553612
.
Ministry of Electronics and Information Technology
2019
Unique Identification Authority of India, Government of India. Available from: https://uidai.gov.in/
Monteiro De Oliveira
E. C.
,
Caixeta
E. S.
,
Santos
V. S. V.
&
Pereira
B. B.
2021
Arsenic exposure from groundwater: environmental contamination, human health effects, and sustainable solutions
.
Journal of Toxicology and Environmental Health, Part B
24
(
3
),
119
135
.
https://doi.org/10.1080/10937404.2021.1898504
.
Paliwal
K. V.
1972
Irrigation with Saline Water
.
Water Technology Centre, Indian Agriculture Research Institute
,
New Delhi
, p.
198
.
Patel
S. V.
,
Khatri
N.
,
Chavda
P.
&
Jha
A. K.
2022
Potential health concerns due to elevated nitrate concentrations in groundwater of villages of Vadodara and Chhota Udaipur districts of Gujarat, India
.
Journal of Water and Health
20
(
1
),
227
245
.
https://doi.org/10.2166/WH.2021.160
.
Paul
R.
,
Brindha
K.
,
Gowrisankar
G.
,
Tan
M. L.
&
Singh
M. K.
2019
Identification of hydrogeochemical processes controlling groundwater quality in Tripura, Northeast India using evaluation indices, GIS, and multivariate statistical methods
.
Environmental Earth Sciences
78
(
15
),
1
16
.
https://doi.org/10.1007/s12665-019-8479-6
.
Pei-yue
L. I.
,
Hui
Q.
&
Jian-hua
W. U.
2010
Groundwater quality assessment based on improved water quality index in Pengyang County, Ningxia, Northwest China
.
E-Journal of Chemistry
7
(
S1
),
209
217
.
Piper
A. M.
1944
A graphic procedure in the geochemical interpretation of water analyses
.
Transactions, American Geophysical Union
25
(
6
),
914
.
https://doi.org/10.1029/TR025i006p00914
.
Ragunath
H. M.
1987
Groundwater
.
Wiley Eastern Ltd
,
New Delhi
, p.
563
.
Raheja
H.
,
Goel
A.
&
Pal
M.
2022
Prediction of groundwater quality indices using machine learning algorithms
.
Water Practice and Technology
17
(
1
),
336
351
.
https://doi.org/10.2166/wpt.2021.120
.
Ram
A.
,
Tiwari
S. K.
,
Pandey
H. K.
,
Chaurasia
A. K.
,
Singh
S.
&
Singh
Y. V.
2021
Groundwater quality assessment using water quality index (WQI) under GIS framework
.
Applied Water Science
11
(
2
),
1
20
.
https://doi.org/10.1007/s13201-021-01376-7
.
Ramakrishnaiah
C. R.
,
Sadashivaiah
C.
&
Ranganna
G.
2009
Assessment of water quality index for the groundwater in Tumkur taluk, Karnataka state, India
.
E-Journal of Chemistry
6
(
2
),
523
530
.
https://doi.org/10.1155/2009/757424
.
Rango
T.
,
Kravchenko
J.
,
Atlaw
B.
,
McCornick
P. G.
,
Jeuland
M.
,
Merola
B.
&
Vengosh
A.
2012
Groundwater quality and its health impact: an assessment of dental fluorosis in rural inhabitants of the Main Ethiopian Rift
.
Environment International
43
,
37
47
.
https://doi.org/10.1016/j.envint.2012.03.002
.
Raza
M.
,
Hussain
F.
,
Lee
J. Y.
,
Shakoor
M. B.
&
Kwon
K. D.
2017
Groundwater status in Pakistan: a review of contamination, health risks, and potential needs
.
Critical Reviews in Environmental Science and Technology
47
(
18
),
1713
1762
.
https://doi.org/10.1080/10643389.2017.1400852
.
Richards
L. A.
1954
Diagnosis and improvement of saline and alkaline soils
.
Soil Science Society of America Journal
18
(
3
),
348
.
https://doi.org/10.2136/sssaj1954.03615995001800030032x
.
Şener
Ş.
,
Şener
E.
&
Davraz
A.
2017
Evaluation of water quality using water quality index (WQI) method and GIS in Aksu River (SW-Turkey)
.
Science of the Total Environment
584–585
,
131
144
.
https://doi.org/10.1016/j.scitotenv.2017.01.102
.
Sethy
S. N.
,
Syed
T. H.
,
Kumar
A.
&
Sinha
D.
2016
Hydrogeochemical characterization and quality assessment of groundwater in parts of Southern Gangetic Plain
.
Environmental Earth Sciences
75
(
3
),
1
15
.
https://doi.org/10.1007/s12665-015-5049-4
.
Shan
V.
,
Singh
S. K.
&
Haritash
A. K.
2021
Evaluation of water quality and potential metal contamination in ecologically important Bhindawas bird sanctuary, India
.
Applied Water Science
11
(
1
),
1
9
.
https://doi.org/10.1007/s13201-020-01334-9
.
Sheikh
M. A.
,
Azad
C.
,
Mukherjee
S.
&
Rina
K.
2017
An assessment of groundwater salinization in Haryana state in India using hydrochemical tools in association with GIS
.
Environmental Earth Sciences
76
(
13
),
1
13
.
https://doi.org/10.1007/s12665-017-6789-0
.
Simsek
C.
&
Gunduz
O.
2007
IWQ index: a GIS-integrated technique to assess irrigation water quality
.
Environmental Monitoring and Assessment
128
(
1–3
),
277
300
.
https://doi.org/10.1007/s10661-006-9312-8
.
Sood
S.
&
Sharma
P.
2020
Unsupervised learning techniques in groundwater quality assessment of Mewat region, Haryana, India
.
Sustainable Water Resources Management
6
(
6
),
1
11
.
https://doi.org/10.1007/s40899-020-00476-7
.
Soujanya Kamble
B.
,
Saxena
P. R.
,
Kurakalva
R. M.
&
Shankar
K.
2020
Evaluation of seasonal and temporal variations of groundwater quality around Jawaharnagar municipal solid waste dumpsite of Hyderabad city, India
.
SN Applied Sciences
2
(
3
),
1
22
.
https://doi.org/10.1007/s42452-020-2199-0
.
Subba Rao
N.
2006
Seasonal variation of groundwater quality in a part of Guntur District, Andhra Pradesh, India
.
Environmental Geology
49
(
3
),
413
429
.
https://doi.org/10.1007/s00254-005-0089-9
.
Szabolcs
I.
&
Darab
C.
1964
The Influence of Irrigation Water of High Sodium Carbonate Content of Soils
. In
Proceedings of 8th International Congress of ISSS, Trans II
, pp.
803
812
.
Tiwari
T. N.
&
Manzoor
A.
1988
River pollution in Kathmandu valley (Nepal) suitability of river water for irrigation
.
Indian J Env Protn
8
(
4
),
269
274
.
Walton
W. C.
1970
Groundwater Resource Evaluation
.
McGraw Hill Book Co.
,
New York
.
Wang
X.
,
Xu
X. Q.
,
Gao
C. H.
,
Li
L. H.
,
Liu
Y.
,
Zhang
N.
,
Xia
Y.
,
Fang
X.
&
Zhang
X. G.
2022
Assessing the drinking water quality in the Inner Mongolia autonomous region from 2014 to 2018
.
Journal of Water and Health
00
(
0
),
1
10
.
https://doi.org/10.2166/wh.2022.217
.
WHO
2011
Guidelines for Drinking-Water Quality
, Vol.
38
(
3
).
WHO Chronicle
, pp.
104
108
.
https://doi.org/10.1248/jhs1956.35.307
.
Wilcox
L.
1955
Classification and Use of Irrigation Waters
.
USDA
,
Washington, DC
,
Circular, 969
.
Wu
C.
,
Wu
X.
,
Qian
C.
&
Zhu
G.
2018
Hydrogeochemistry and groundwater quality assessment of high fluoride levels in the Yanchi endorheic region, northwest China
.
Applied Geochemistry
98
,
404
417
.
https://doi.org/10.1016/j.apgeochem.2018.10.016
.
Zhang
Y.
,
Wu
J.
&
Xu
B.
2018
Human health risk assessment of groundwater nitrogen pollution in Jinghui canal irrigation area of the loess region, northwest China
.
Environmental Earth Sciences
77
(
7
),
273
.
https://doi.org/10.1007/s12665-018-7456-9
.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY 4.0), which permits copying, adaptation and redistribution, provided the original work is properly cited (http://creativecommons.org/licenses/by/4.0/).