Introduction

Groundwater is the main source for irrigation, domestic and drinking purposes in many basins of the Indian peninsula (Siva Prasad and Venkateswara Rao 2018). Groundwater plays an important role for supplying the water to irrigation in Subledu Basin even though there are many tanks and canals existing in the basin. Due to overirrigation, more groundwater is withdrawn; as a result, the groundwater quality also gets degraded in addition to depletion of groundwater quantity. In order to know the status of groundwater quality, groundwater samples were analyzed in the basin as no such study was carried out earlier in the basin. However, some of the studies carried out in the similar hydrogeologic settings in the Indian peninsula indicate that there are a fluoride contamination due to geogenic sources and nitrate pollution due to excess use of fertilizers are very much in vogue causing skeletal fluorosis and methemoglobinemia respectively (Mehdi et al. 2018). In order to verify the suitability of groundwater quality for the drinking and irrigation purposes, physicochemical analysis of water samples in the study area was carried out.

According to Vodela et al. (1997) and Moharir et al. (2017), the groundwater pollution in the world is a serious issue. In an area, the groundwater quality can be described in the form of a function which includes physical, chemical and microbiological parameters of the groundwater (Subramani et al. 2005; Schiavo et al. 2006; Magesh and Chandrasekar 2011; Krishna Kumar et al. 2011). In India, the main reason for groundwater pollution is over irrigation (Pawar and Shaikh 1995; Sujatha and Reddy 2003). The groundwater quality data give the information of groundwater recharge, discharge and storage (Walton 1970). Hydrogeochemistry information is also required to estimate the groundwater quality of an area where the groundwater is to be used for drinking and irrigation purposes (Srinivas et al. 2013). In India, nearly 80% of the rural population and around 50% of the urban population are dependent on groundwater for domestic purposes (Biswas et al. 2014). At several places, fluoride and nitrate concentrations in the groundwater of rural agricultural areas of Telangana state, India, were reported to be more than the permissible limits of World Health Organization (WHO) and Bureau of Indian Standards (BIS) (Alekhya et al. 2018; Shaik Arshad et al. 2015; Vijaya lakshmi et al. 2016; Annadasankar Roy et al. 2018).

By using the remote sensing and Geographical Information System (GIS) software, Srinivasa Rao and Jugran (2003) have prepared the groundwater quality maps as well as groundwater potential zones in the hard rock terrains of Chittoor. Inverse distance weighted (IDW) method of interpolation technique was used to generate spatial distribution maps of various irrigation water parameters (Gagandeep Singh et al. 2019).

Groundwater suitability for irrigation was assessed by salinity, chlorinity, sodicity and other parameters (Mills 2003; Nishanthiny et al. 2010). According to Chaitanya et al. (2018), spatial interpolation and statistical methods were used for mapping the groundwater quality parameters in basaltic hard rocks of Akola and Buldhana districts of Maharashtra, India. The consumption of nitrate and fluoride due to contaminated water causes the health risk, which is in the decreasing order of children, women and men. It depends upon sensitiveness to pollutants and smaller body weights of children, followed by women and men (Li et al. 2016; Ahada and Suthar 2017; Yang et al. 2018; Karunanidhi et al. 2020a, b).

The present work is aimed to find out the quality of groundwater for irrigation and drinking purposes apart from verifying geogenic and anthropogenic pollution of fluoride and nitrate, respectively. The spatial distribution maps of important chemical constituents are prepared by using the contour maps created by utilizing the IDW interpolation tool in the Geographical Information System (GIS) for the drinking water quality and irrigation water quality in order to know about the polluted and nonpolluted areas as done by earlier researchers (Jasrotia et al. 2018, 2019; Adimalla and Taloor 2020).

Description of the study area and methods

Location and geology

The study area, namely Subledu Basin, is located in Khammam district, Telangana, India, as shown in Fig. 1. The areal extent of Subledu basin is 130.29 km2 and falls in the toposheet No. E44N15 and E44N16 on the 1:50,000 scale. The toposheets are collected from the Survey of India. Geographically, the area exists between the latitudes 17°13' N to 17°24' N and the longitudes 79°51' E to 79°59' E.

Fig. 1
figure 1

Location of the study area map of Subledu Basin and the groundwater samples

The area of the basin is basically underlain by granites. Gray granites are exposed in the middle parts of the basin. Pink hornblende granites are prominently observed in the villages, namely Sitarampuram, Tallacheruvu, and also can be observed at eastern parts, apart from other isolated places. Pink hornblende biotite gneiss exists in the southeast of the basin. The geological map of the study area is presented in Fig. 2, and it is collected from the Geological Survey of India (GSI) Southern Region, Hyderabad. The average annual rainfall of the basin is 868 mm (TSDPS 2020).

Fig. 2
figure 2

Geology map of the Subledu Basin

Topography and drainage

Figure 3 depicts the digital elevation model (DEM) of Subledu Basin. Minimum and maximum elevation range is 55 m to 188 m above the mean sea level. The longitudinal elevation is falling from north-west to south of the basin causing the main stream flow in that direction. Similarly, the latitudinal elevations are also falling from west to center of the basin as well as east to center of the basin causing lower order streams joining the main stream from both the east and west directions.

Fig. 3
figure 3

Subledu basin digital elevation model map

The drainage map of the Subledu Basin (Fig. 4) indicates that the area is well drained with dendritic pattern of drainage and has considerable number of water bodies in the form of tanks and act as recharging structures.

Fig. 4
figure 4

Drainage map of the Subledu Basin

Materials and methods

To assess the groundwater quality for drinking and irrigation purposes, 22 groundwater samples were collected in the month of May 2019. The samples were collected in one liter capacity of well-cleaned plastic bottles and sealed tightly. The samples were analyzed at Telangana State Groundwater Department, Hyderabad, using APHA (American Public Health Association 1998) methods. The results were compared with BIS (Bureau of Indian Standards 2012) values. Table 1 shows the Bureau of Indian Standards of drinking water as well as WHO (World Health Organization 2011) standards.

Table 1 Groundwater quality values of Bureau of Indian Standards (BIS) (IS 10500:2012) and WHO (2011) standards (Source: State Groundwater Department, Telangana)

The major anions and cations, namely carbonate, bicarbonate, nitrate, sulfate, phosphate, sodium, potassium, fluoride, chloride, calcium and magnesium, were determined apart from other parameters such as pH, electrical conductivity (EC) and total dissolved solids (TDS). By using the ArcGIS 10.2 software, the groundwater quality maps were also prepared and discussed. For creation of contour maps in the study area, the contour tool option was used in the spatial analyst tools of ArcGIS 10.2.

The instrument, namely Water Quality Analyzer (Elico PE 138), was used for determination of pH and electrical conductivity. Carbonate (CO3) and bicarbonate (HCO3) were determined by using acid titrimetric method. By using the spectrophotometer, the sulfate and phosphate ion concentrations were determined. For evaluation of total dissolved solids (TDS), the method of gravimetric evaporation was used. For determination of fluoride concentration, the ion-selective electrode method was used and the instrument used for the analysis was ion meter. The argentometric (AgNO3) titration method is used for chloride analysis. For the nitrate calculation, UV spectrometric method is followed. The total hardness of groundwater is estimated with EDTA titration.

In order to know the suitability of groundwater for irrigation, various derived parameters, namely sodium hazard, sodium percentage, residual sodium carbonate, magnesium hazard and Kelley’s ratio, are calculated. Sodium hazard for irrigation water is determined by finding the sodium adsorption ratio (SAR). The formula for SAR is as follows.

$${\text{SAR}} = \frac{{{\text{Na}}^{ + } }}{{\sqrt {\frac{{{\text{Ca}}^{ + 2} + {\text{Mg}}^{ + 2} }}{2}} }}$$
(1)

The formula used to calculate the sodium percentage (SP) is

$${\text{Sodium Percentage}} = \left[ {\frac{{{\text{Na}}^{ + } + {\text{K}}^{ + } }}{{{\text{Ca}}^{ + 2} + {\text{Mg}}^{ + 2} + {\text{Na}}^{ + } + {\text{K}}^{ + } }}{ }} \right]{\text{x}}100$$
(2)

According to Eaton 1950, the residual sodium carbonate can be calculated by using the formula and it can be expressed in meq/l as.

$${\text{Residual Sodium Carbonate }}\left( {{\text{RSC}}} \right) \, = \, \left[ {{\text{HCO}}^{ - }_{3} + {\text{CO}}^{ - }_{3} } \right] \, - \, \left[ {{\text{Ca}}^{ + 2} + {\text{Mg}}^{ + 2} } \right]$$
(3)

Magnesium hazard (MH) as proposed by Szabolcs and Darab (1964) is expressed as.

$${\text{MH\% }} = \left[ {\frac{{{\text{Mg}}^{ + 2} }}{{{\text{Ca}}^{ + 2} + {\text{Mg}}^{ + 2} }}} \right]{\text{x}}100$$
(4)

Kelley’s ratio is calculated by using the formula.

$${\text{Kelley}}'{\text{s}}\;{\text{Ratio}} = \left[ {\frac{{{\text{Na}}^{ + } }}{{{\text{Ca}}^{{ + 2}} + {\text{Mg}}^{{ + 2}} }}} \right]$$
(5)

In all the above relations, the concentrations of anions and cations are expressed in milliequivalents per liter (meq/l).

In order to know the relative ionic dominance in the groundwater and its suitability, the hydrogeochemical data are analyzed with the piper trilinear plot and USSL diagram by using the AquaChem 4.0 software.

Results and discussion

The data so analyzed are presented in various tables. Tables 2 and 3 present the cation analysis data and the anion analysis data of groundwater samples respectively along with standard deviation, mean, variance, maximum and minimum values. From Tables 2 and 3, the cation dominance order is \({\text{K}}^{ + } > {\text{ Na}}^{ + } > {\text{Mg}}^{ + 2} > {\text{Ca}}^{ + 2}\). Similarly, for anions the order of dominance is\({\text{HCO}}^{ - } {\text{ > NO}}_{{3}}^{ - } {\text{ > Cl}}^{-} {\text{ > SO}}_{{4}}^{{ - 2}} {\text{ > CO}}^{{ - 2}} {\text{ > F}}^{ - } {\text{ > O - PO}}_{{4}}^{{ - 3}}\).

Table 2 Chemical analysis of groundwater for cations and other parameters.
Table 3 Chemical analysis of groundwater for anions and other parameters.

Suitability of groundwater for drinking purposes

Total dissolved solids (TDS)

In natural groundwater, a large number of salts are found dissolved. They are formed out of anion and cation combination. The major ions among them are carbonate, bicarbonate, chloride, sulfate, phosphate, nitrate, calcium, magnesium, sodium, potassium, iron, manganese, etc. TDS is a measure of the amount of material dissolved in the water. According to Mitchell and Stapp (1992), the water flow depends on TDS concentration in the organism’s cells. Sometimes, death can be caused due to high or low concentration of TDS values. In order to estimate the drinking water quality, TDS plays a very important role. The water having more TDS results in salty taste. Prolonged intake of high TDS groundwater causes the kidney stones and develop heart problems (Ali and Ali 2018).

The TDS results varied from 469 mg/l at Tallacheruvu village to 1702 mg/l at Kistapuram village. The high values of TDS are due to excess usage of fertilizers for irrigation purposes, and the average value of TDS is 961 mg/l. Out of 22 groundwater samples analyzed, only at two villages, namely Tallacheruvu and Madiripuram, are within the BIS limits, rest all the TDS values have exceeded the BIS limit. For TDS, the standard deviation, mean and variance are 401 mg/l, 961.95 mg/l and 160,806 (mg/l)2, respectively.

From Fig. 5, it can be observed that the southern, western and northern parts of the basin are having more TDS concentration compared to the eastern part of the Basin.

Fig. 5
figure 5

Total dissolved solids distribution in the Subledu Basin

pH

The negative logarithmic of hydrogen ion concentration is called pH. pH tells how acidic or alkaline the water is, and its values are evaluated by using the pH meter for all the 22 samples. If the pH is not in the desirable limit of 6.5–8.5 in the potable water, it harms the mucous membrane present in eyes, nose, mouth, abdomen, anus, etc. (Ibrahim et al. 2015). In the study area, pH values are in the range of 7.53–8.35. All the pH values are within the BIS range. The villages Turkagudem and Pocharam are having maximum pH of 8.35 value, whereas Bachodu is having the minimum pH of 7.53. All the values of pH in the basin are more than 7 which indicates that the groundwater is toward alkaline in nature. The standard deviation, mean and variance for pH are 0.19, 8.09 and 0.03, respectively.

Electrical conductivity (EC)

The water can act as a conductor of electricity due to the presence of dissolved salts where as pure water cannot pass the electricity through it. The conductivity is caused in the water when the higher concentration of electrolytes are present. Conductance is measured in mhos or Siemens. In the present study, EC is varying from 733 to 2660 (µS/cm). The minimum EC of 733 (µS/cm) is found at Tallacheruvu village and the maximum EC of 2660 (µS/cm) is found at Kistapuram village. The high conductivity values indicate increase in salinization of the groundwater. The standard deviation, mean and variance for EC are 626.68 µS/cm, 1503 µS/cm and 392,736 (µS/cm)2, respectively.

Fluoride

The presence of fluoride in the groundwater is due to the fluoride-bearing rocks like fluorspar, cryolite, fluorapatite, etc. Due to the improper presence of fluoride, dental and skeletal problems can occur. If the fluoride concentration is less than 0.4 mg/l in the drinking water, then dental fluorosis occurs but if fluoride concentration is greater than 1.5 mg/l, skeletal fluorosis will occur (Yadav and Lata 2004). Karunanidhi et al. (2019) stated that the contamination of groundwater is caused by F-containing minerals due to the presence of the country rocks of hornblende–biotite gneiss and granites as well as by leaching of waste disposals and uncontrolled usage of agrochemicals.

In the basin, the fluoride is ranging from 0.41 to 2.84 mg/l. Except at three villages, namely Turkagudem (1.87 mg/l), Kistapuram (2.76 mg/l) and Pocharam (2.84 mg/l), the fluoride concentrations are below the 1.5 mg/l which are in permissible (BIS 2012) limit for consumption of groundwater. The standard deviation, mean and variance are 0.66 mg/l, 1.08 mg/l and 0.43 (mg/l)2, respectively, for the fluoride concentration. From Fig. 6, it can be observed that the fluoride concentrations are greater than 1.5 mg/l in the southern part of the basin where it is not suitable for drinking groundwater without treatment for the fluoride.

Fig. 6
figure 6

Fluoride distribution in the Subledu Basin

Chloride

The presence of chloride is caused in the water due to animal organic wastes. There are many reasons to cause the presence of chloride in the surface water, like effluents from water treatment plants, wastes from oil wells and agricultural wastes joining the surface water. The effect of leaching which is contaminated with the chlorides is also one of the reasons for increase in chloride concentration in the groundwater (Central Ground Water Board 2005).

The salty taste to the ground water is due to the presence of higher Cl, and it shows a laxative effect. Most of the Cl released from a nonlithological sources into the groundwater is due to anthropogenic sources of pollution caused by household wastewaters, (MPCA 2019; Subba Rao et al. 2019a). Chloride causes the corrosion in the metals and also changes the taste of food products (Hilbert et al. 2006; Malcolm et al. 2017). The minimum chloride concentration of 40 mg/l is found in the village Madiripuram, and maximum chloride concentration of 340 mg/l is found in the village Kistapuram. The high concentrations of chloride are found at villages, namely Kistapuram—340 mg/l, Bachodu—330 mg/l, Mekalthanda—320 mg/l, Maripeda Bunglow, Narasimhulagudem and Turkagudem—310 mg/l. Standard deviation, mean and variance of chloride are 105.49 mg/l, 162.72 mg/l and 11,128 (mg/l)2, respectively. From Fig. 7, it can be observed that high chloride concentrations are present at southern, western and some parts of north and northeast parts of the basin.

Fig. 7
figure 7

Chloride distribution in the Subledu Basin

Nitrate

The nitrate contamination in the groundwater is due to the excess use of fertilizers (Sunitha 2013). The animal, human and agricultural activities create the nitrate pollution in the groundwater. According to He and Wu (2019) and He et al. (2019), household wastewaters and septic tank leakages as well as the soils containing the nitrogen and animal wastes are the sources of NO3 in the groundwater. In the drinking water, generally histomorphological changes in thyroid are observed due to the presence of 250–500 mg/l of nitrates in the groundwater for drinking (Eskiocak et al. 2005). The high concentration of nitrate in the groundwater causes the methemoglobinemia when we drink directly. Not more than 45 mg/l is allowed for drinking water. In the study area, the standard deviation, mean and variance values of nitrate are 128.32 mg/l, 104.31 mg/l and 16,467 (mg/l)2, respectively.

In the present study, the nitrate concentrations range from 2 to 460 mg/l. Ten groundwater samples are below 45 mg/l and another 12 groundwater samples are above 45 mg/l. From Fig. 8, it can be observed that the high contamination of nitrates is in the western, southern and northeast parts of the basin. The village Medidapalle is having the low nitrate concentration of 2 mg/l, while the village Bachodu is having the high nitrate concentration of 460 mg/l.

Fig. 8
figure 8

Nitrate distribution in the Subledu Basin

Calcium, magnesium and total hardness

The Ca+2 is an essential inorganic element for proper growth of bones, which is released into the groundwater due to weathering and dissolution of plagioclase feldspars occurring in the parent rocks. The occurrence of ferromagnesium minerals (hornblende, hypersthene, sillimanite, garnet and biotite) in the basement rocks is the sources of Mg+2 in the subsurface water, in addition to the sources such as domestic wastewaters, septic tank leakages, etc. (Subba Rao et al. 2017a).

The calcium concentrations range from 16 to 112 mg/l. The village Sitarampuram is having the minimum calcium concentration of 16 mg/l, while the village Bachodu is having the maximum calcium concentration of 112 mg/l. For calcium, the standard deviation, mean and variance are 25.77, 45.81 mg/l and 664 (mg/l)2, respectively. The magnesium concentrations range from 15 to 165 mg/l. The minimum magnesium of 15 mg/l is found at Madiripuram village, while the maximum magnesium of 165 mg/l is found at Mekalthanda village. For magnesium, the standard deviation, mean and variance are 42.57, 78.90 mg/l 1812 (mg/l)2, respectively.

The hardness of water is sum of the polyvalent cations dissolved in it. The dominant cations are calcium and magnesium. The calcium carbonate (CaCO3) hardness is termed as hardness of the water. According to US Department of Interior and the Water Quality Association, the total hardness can be classified as shown in Table 4.

Table 4 Total hardness classification of groundwater quality (Prakash and Somashekar 2006).

In the study area, total hardness values range from 200 to 820 mg/l. The minimum total hardness of 200 mg/l is found at Madiripuram village, and the maximum total hardness of 820 mg/l is found at Bachodu Village. The total hardness has exceeded the limits of Bureau of Indian Standards at four villages. In the present study, all the groundwater samples are in hard to very hard category. From Table 4, we can say that Subledu Basin is facing a problem with very hard water. Standard deviation, mean and variance of the total hardness are 164.53 mg/l, 439.09 mg/l and 27,071 (mg/l)2, respectively.

Sodium and potassium

According to Subba Rao et al. (2017a, 2019a), the primary sources of Na+ are due to weathering of the plagioclase feldspars. Higher Na+ content reflects the secondary source, which is a result of influence of anthropogenic sources like household wastewaters, leakage of septic tanks and irrigation practice.

The sodium concentrations vary from 30  to 287 mg/l and the minimum value is found in the village Maripeda and the maximum value is observed at Maripeda Bunglow. In the present study, highest concentration of potassium is 366 mg/l and it is found at Hasnabad and the minimum value is 2 mg/l and is found at Maripeda, Pocharam and Eluvarigudem villages. Due to high consumption of potassium, organs such as kidneys and heart may not function properly. The standard deviation, mean and variance of the sodium are 79.66 mg/l, 131.81 mg/l and 6347 (mg/l)2, respectively, and for the potassium the standard, mean and variance values are 75.30 mg/l, 26.36 mg/l and 5671 (mg/l)2, respectively. Five samples in the case of sodium and one sample in the case of potassium have exceeded the permissible limit.

Sulfate and phosphate

Sulfate is generally found in natural water due to discharges from domestic sewage, effluents released from the industries and emissions from automobiles. The sulfate concentrations range from 9 mg/l at Hasnabad village to 51 mg/l at Narasimhulagudem village. In the present analysis of groundwater samples, all the sulfate values are below the permissible limit of BIS; hence, there is no problem with the sulfate. The sulfate’s standard deviation, mean and variance are 13.12 mg/l, 25.36 mg/l and 172.23 (mg/l)2, respectively.

As per WHO, the permissible value of phosphate concentration in the groundwater is 0.1 mg/l. Since no BIS standards exist for phosphate in India, according to Rajmohan and Elango (2005), the soils holding capacity of phosphate is low when it is less soluble. The maximum value of phosphate is found to be 0.42 mg/l at Hasnabad village, and the minimum value is found to be 0.04 mg/l at Tallacheruvu; only one sample is exceeding the permissible limit.

Integrated groundwater quality map for drinking water purposes

Figure 9 shows the integrated map of groundwater quality, and it is prepared using the IDW tool in ArcGIS 10.2 software. This map shows the parts of study area where the groundwater is suitable and not suitable for drinking purposes. It is prepared based on nitrate, fluoride, total hardness, total dissolved solids and chloride. Cool to hot colors indicate degree of suitability of groundwater for drinking purposes. Dark-blue color and light-blue color areas are suitable for drinking purposes, while all other areas are not suitable for drinking purposes. In the southern and western parts of the basin as well as few parts in the middle of the basin, the groundwater quality is poor compared to north, east and few parts in the west, south and eastern boundaries of the study area.

Fig. 9
figure 9

Integrated groundwater quality map of the Subledu Basin for drinking purposes

Suitability of groundwater for irrigation

According to Wagh et al. (2019) and Deepali et al. (2020a), the groundwater quality used for irrigation purposes plays an important role for growth and production of the crop. Physical properties such as soil permeability are very sensitive to type of exchangeable ions. By using salinity and sodium hazards, the suitability of groundwater for irrigation can be determined (Subba Rao 2008; Chen et al. 2019; Karunanidhi et al. 2020a, b). The crop yields depends on Salinity hazard, and it can reduce the crop yield due to excess of salt content in the groundwater. Similarly, the sodium hazard can also affect the soil permeability and reduce the infiltration rates of water and air into soil structure (Aravinthasamy et al. 2020). If the salinity and sodium hazards are greater, the higher is the risk of suitability of groundwater quality for irrigation purpose (Subba Rao 2017; Tolera et al. 2020).

Sodium adsorption ratio (SAR)

The permeability and infiltration rates decrease due to high levels of sodium concentrations. The sodium concentration can also effect the less sensitive crops (Zaman et al. 2018; Sarath Prasanth et al. 2012). With reference to SAR (Todd 1980), the irrigation water can be classified as presented in Table 5. According to Todd classification, almost all the samples in the study area are of excellent quality for irrigation.

Table 5 Irrigation water quality classification based on SAR (Todd 1980).

The SAR values are changing from 0.68 to 5.58 (Table 3). With the usage of water of high SAR values continuously for irrigation, there will be breakdown in the soil structure physically (Nagarajan et al. 2009). The standard deviation, mean and variance of the SAR are 1.35, 2.67 and 1.82, respectively. By using the SAR and specific conductivity, we can draw the US Salinity Laboratory (USSL) diagram (Fig. 10).

Fig. 10
figure 10

Groundwater quality classification in the Subledu Basin using US salinity diagram

Irrigation water quality based on salinity

Classification of irrigation water based on salinity is shown in Table 6 for the groundwater of the study area. Only three samples were not in permissible category. All other 19 samples were suitable for irrigation.

Table 6 Classification of irrigation water based on salinity (EC) (Richards 1954, Ntokozo Malaza 2017).

For irrigation water classification, US Salinity Laboratory (Richards 1954) diagram is also used in the Subledu Basin. In order to ensure the groundwater quality for irrigation purpose, United States Soil Salinity Laboratory staff (USSLS) diagram can be  used (Subba Rao 2017; Wagh et al. 2018; Deepali et al. 2020a). The graph (Fig. 10 ) is drawn by taking the SAR on Y- axis and the specific electrical conductance on X- axis. According to Hem (1985), this diagram is used for finding the salinity problems and exchange of ions. Based on the salinity hazard, it is divided into four groups, low (C1), medium (C2), high (C3), and very high (C4), and it is shown in Table 6. In the same way, sodium alkali hazard is also divided into low (S1), medium (S2), high (S3) and very high (S4). In the present study area, maximum groundwater samples fall under the C3S1 category indicating that the values are in the high salinity and low alkali hazard range, and it is shown in Fig. 10. It means that there is no sodium hazard but salinity hazard is present.

Sodium percentage (SP)

In order to estimate the sodium hazard in irrigation water, we use the sodium percentage. Based on sodium percentage, only one groundwater sample is doubtful. Ten groundwater samples are within permissible limits, nine samples are good and two samples are excellent. The standard deviation, mean and variance of the sodium percentage are 13.62, 38.87 and 185, respectively. Based on sodium percentage, the classification of groundwater samples is shown in Table. 7.

Table 7 Classification of irrigation water based on sodium percentage (Wilcox 1955).

Residual sodium carbonate (RSC) index

The high-concentrated bicarbonate and carbonate ions accompanied with alkaline earths are the possible cause of precipitation of Ca and Mg (Raghunath 2000). The residual sodium carbonate equals to the sum of carbonates and bicarbonate concentrations minus calcium and magnesium concentrations. The standard deviation, mean and variance of residual sodium carbonate are 3.15 mg/l, 0.7 mg/l and 9.96 (mg/l)2, respectively.

Residual sodium carbonate values range from 10.07 to 6.07 meq/l. Table 8 shows the irrigation water classification of groundwater in the study area based on RSC values. According to the this classification, we can say that 17 groundwater samples are having no hazard and three groundwater samples are in medium hazard and two are in high hazard for irrigation. The US Department of Agriculture stated that the water having RSC value more than 2.5 meq/l is not suitable for irrigation purposes (Ntokozo Malaza 2017).

Table 8 Classification of irrigation water based on residual sodium carbonate (Eaton 1950; Wilcox et al. 1954).

Magnesium hazard (MH)

The crop yield can be decreased due to the presence of magnesium in groundwater (Gowd 2005). If the magnesium hazard value is greater than 50%, then the water is not suitable for irrigation (Siva Prasad and Venkateswara Rao 2018). Magnesium hazard values in the study area range from 30.27% to 94.51%. Only five samples are less than 50%, and rest of the samples have exceeded the 50%. According to magnesium hazard, only 22.72% of groundwater in the study area is suitable for irrigation and the rest 77.28% is not suitable. The standard deviation, mean and variance of the magnesium hazard are 19.39, 70.05 and 376, respectively.

From Fig. 11, it can be observed that the magnesium hazard is more at southern and eastern parts of the basin compared to northern and western parts.

Fig. 11
figure 11

Distribution of magnesium hazard in the Subledu Basin

Piper’s trilinear classification

In the Piper’s trilinear diagram (Piper 1944), it is having two triangles as shown in Fig. 12, one is used for plotting cations and another is used for plotting anions (Fig. 12). The resultant of the combination of the two triangles can be expressed in the shape of a diamond field.

Fig. 12
figure 12

Piper’s trilinear classification of groundwater in the Subledu Basin

Figure 12 shows the Piper’s trilinear classification of the 22 groundwater samples of the study area. From Fig. 12, we can say that most of the samples, i.e., 16 samples, are in the area No.5. From the legend, we can infer that this area shows the carbonate hardness (secondary alkalinity) exceeds 50 percent. The remaining groundwater samples, i.e., six samples, are in the area No.9, indicating that it is the neutral water. Therefore, the secondary alkalinity is the dominant hydrochemical facies of groundwater in the study area.

Kelley’s ratio

Kelley et al. (1940) stated that the high presence of sodium in the irrigation water causes the retardation of plant’s growth. It is expressed as Kelley’s ratio and if the ratio is less than 1, the groundwater quality is safe for irrigation. If it is greater than 1, then it is not safe for irrigation. In the present study area, the Kelley’s ratio is varying from 0.16 to 1.37. The standard deviation, mean and variance of the Kelley’s ratio are 0.31, 0.63 and 0.09, respectively.

In the study area, out of 22 groundwater samples 20 samples are in safe range for irrigation and two groundwater samples have exceeded Kelley’s ratio; they are Maripeda Bunglow and Hasnabad. Based on Kelley’s ratio, the groundwater samples are classified in Table 9.

Table 9 Classification of groundwater quality based on Kelley’s ratio (Kelley 1963).

Integrated groundwater quality map for irrigation

Figure 13 shows the suitability of groundwater for irrigation. In this map, residual sodium carbonate, EC and Kelley’s ratio are integrated and the map is prepared. For irrigation purpose, the groundwater is unsuitable at few places in the western middle parts of the basin as well as north and southern tips of the basin as shown in Fig. 13. The large rest of the basin is suitable for irrigation.

Fig. 13
figure 13

Groundwater quality map of the Subledu Basin for irrigation purposes

Conclusions

Twenty-two groundwater samples were analyzed in the Subledu Basin for chemical composition and their suitability for irrigation and drinking purposes. The fluoride concentrations have exceeded the BIS standards at three places and nitrate concentration at 12 villages; therefore, totally 13 groundwater samples are not suitable for drinking purpose. The integrated groundwater quality map for drinking purposes shows that in the large parts of southern and western parts of the basin as well as in the middle of the basin, the groundwater is not suitable for drinking. For irrigation purpose, Kelley’s ratio is more than 1 at two villages only. According to USSL salinity diagram, we can conclude that maximum groundwater samples fall under the C3S1 category that means the samples are having high salinity and low alkali hazard range. The cation dominance order is, \(K^{ + } > \, Na^{ + } > Mg^{ + 2} > Ca^{ + 2}\). Similarly, for anions the order of dominance is \({\text{HCO}}^{ - } {\text{ > NO}}_{{3}}^{ - } {\text{ > Cl}}^{-} {\text{ > SO}}_{{4}}^{{ - 2}} {\text{ > CO}}^{{ - 2}} {\text{ > F}}^{ - } {\text{ > O - PO}}_{{4}}^{{ - 3}}\). In the case of TDS except two, all the groundwater samples are having high TDS concentration and exceeded the BIS desirable limit. Based on total hardness results, all the groundwater samples in the basin are very hard water. From the Piper’s trilinear diagram, the secondary alkalinity is the dominant hydrochemical facies of groundwater in the study area. The integrated groundwater quality map for irrigation shows that the groundwater is suitable for irrigation except in few packets of western middle parts of the basin as well as north and southern tips of the basin.