Groundwater quality assessment for domestic and agricultural purposes using GIS, hydrochemical facies and water quality indices: case study of Rafsanjan plain, Kerman province, Iran

This paper investigates the suitability of groundwater for domestic and agricultural use in the Rafsanjan plain, southwest part of the Daranjir–Saghand basin, Iran. Fifty-five groundwater samples were collected and analyzed by six methods including the water quality index (WQI), Schoeller diagram, irrigation water quality (IWQ) parameters, Piper diagram, US salinity diagram and Wilcox diagram. The spatial distribution maps of chemical parameters and groundwater quality indices were plotted using the IDW method in GIS. The results showed a low concentration of major ions in the southeastern part and a high concentration from the central part towards northwestern part of the plain. The concentration of major ions in groundwater was strongly affected by groundwater flow, geological setting and the existence of the evaporative layers in the studied area. Moreover, results revealed that most of samples exceeded the acceptable limits recommended by the WHO and ISIRI1053 standards for domestic and agricultural purposes. In most of the wells, groundwater was classified into saline and very hard categories. The analyses based on WQI values indicated that above 87% of water samples were unsuitable for drinking purposes. IWQ parameters expressed that 85%, 67%, 32%, 51%, 43% and 50% of samples had EC > 3000, Na% > 60%, MAR > 50%, KR > 1, SAR > 9 and Cl− > 350, respectively, which were unsuitable for irrigation use. The dominant hydrochemical facies of water was Na–Cl–SO4 type, and 63% and 22.8% of samples were categorized as C4S4 and C4S3 class, with very high-salinity–high-sodium hazards and very high-salinity–high-sodium hazards, respectively. It indicated that most irrigated lands in this study area were affected by different levels of salinity and sodicity hazards that caused decreases in plant growth and crop productivity. The results can assist decision-makers and planners in prioritizing groundwater resources management in the region.


Introduction
In arid and semiarid areas, groundwater is a vital resource for domestic, industrial and irrigation uses (Delgado et al. 2010, Yousefi et al. 2018. In these areas, over-exploitation of groundwater resources and successive long-term droughts in recent years caused decreases in groundwater levels and deterioration of groundwater quality by the accumulation of dissolved ions in water (Abbasnia et al. 2019, Chidambaram et al. 2022. Therefore, assessments of groundwater quality received more attention worldwide to provide safe water for various applications (Srinivas et al. 2013) and the hydrogeochemical behavior of groundwater was evaluated using various techniques such as binary plots, ionic ratios, Wilcox plots, Piper diagrams, Schoeller diagram, USSL diagrams and GIS techniques (Chidambaram et al. 2022). Previous studies assessed water quality for irrigation and drinking purposes by considering major parameters such as electrical conductivity (EC), total dissolved solids (TDS), total hardness (TH), pH, SAR and major chemical elements, e.g., Ca 2+ , Mg 2+ , Cl − , HCO 3 − , Na + , K + and SO 4 2− (Jamshidzadeh and Mirbagheri 2011, Pour et al. 2014, Alavi et al. 2016, Piroozfar et al. 2018, Abbasnia et al. 2019, Palmajumder et al. 2021, Rajmohan et al. 2021, Arota et al. 2022. Moreover, several studies applied GIS techniques, such as IDW and Kriging interpolation methods, in order to show the spatiotemporal variations of water quality that provide a better understanding of future water resources quality trends (Mosaferi et al. 2014, Alavi et al. 2016, Makki et al. 2021, Silva et al. 2021. Applying GIS in mapping of major parameters of groundwater quality has a significant impact on regional water management and future decisionmaking (Mohammadi et al. 2017, Chidambaram et al. 2022. This study was conducted in Iran, Kavir Daranjir-Saghand basin, which is classified as one of the countries facing limited water resources due to its location in the arid and semiarid belt and lack of rainfall. Thus, the longterm rainfall statistics in the last 5 decades have shown a downward trend in annual precipitation to 243 mm (Mohammadjani and Yazdanian 2014). Kavir Daranjir-Saghand is a second-order watershed with an area of 50,737 km 2 in the central part of Iran. This basin is divided into 12 sub-basins including Rafsanjan, Kavir Daranjir, Kerman Baghin, Bafgh, Zarand, Arnan Dahaj, Bardsir, Bahadoran, Kohbanan, Gharyatolarab, Siriz Toghroljerd and Saghand. Rafsanjan plain is the biggest sub-basin in the Kavir Daranjir-Saghand basin and is located in the southwest part of this region. It is classified as a semiarid region, with 100 mm annual precipitation and mean annual potential evapotranspiration of 3000 mm (Hosseinifard and Mirzaei Aminiyan 2015). With population growth and the development of pistachio orchards in recent years, the high contribution of water needed by domestic, agricultural and industrial use highly depends on groundwater resources in the study area. Therefore, monitoring the quality of groundwater resources by water quality index (WQI) and standards (WHO and ISIRI1053), Schoeller diagram, irrigation water quality (IWQ) parameters, Piper, USSL and Wilcox diagram using GIS interpolation methods is the aim of conducting this study. The results can support groundwater resources management and its suitability for different consumptions in the coming years.

Study area
Rafsanjan plain is the largest plain with an area of 12,513 Km 2 in the southwest part of the Kavir Daranjir-Saghand basin (covers about 25% of the basin), between longitude 290,000 to 470000E and latitude 331,000 to 3480000N in UTM projection system, and the area of the Rafsanjan aquifer is 4236 Sq. Km (Fig. 1a). The Rafsanjan plain consists of three main plains, namely Rafsanjan, Nough and Anar. The average height of the area is 1469 m above sea level. The high mountainous area is located in the west (Sarcheshmeh Mountains (Dehj-Sardoueye belt)) and east (Davaran Mountains) of this plain sloping down toward the central part of the plain (Fig. 1b).
The climate of the subjected area is semiarid, due to low annual rainfall (average 100 mm) and high annual potential evapotranspiration (average 3000 mm). Groundwater has been used for different purposes, such as domestic, agricultural and industrial usage in this area. In addition, irrigation of pistachio orchards is the main agricultural water use in this plain. There are 92,108 hectares of pistachio orchard and farmlands in the Rafsanjani plain ( Fig. 1c) that uses 492 million cubic meters of groundwater per year.

Geological and hydrogeological profile of the aquifer
Surficial geological setting of the area composed of three main units: (1) alluvial and river deposits near the mountains in the southwest, west and northeast, (2) clay, salty and mud flat in the central part toward north and northwest, and (3) sand dunes covered mostly east and northeast part of the plain (Fig. 2 a). However, in the Rafsanjan plain, we are dealing with a complex setting of sedimentary and volcanic rocks as shown in Fig. 2b; the volcanic rocks are located in the south and south west of the study area (Sarcheshmeh Mountains (Dehj-Sardoueye belt)).
The aquifer geology and geohydrology consist of quaternary sediments that build an unconfined alluvial aquifer. Moreover, any compressed clayey layer that could be a potential for a confined aquifer has not been determined based on the analysis of geoelectric sounding profiles. The aquifer mainly consists of sand and clay and its thickness is around 100 m in the southwest and central parts of the plain. The aquifer thickness decreases toward the northwest of the plain due to the reduction in the depth of bedrock that causes shallow water table depth and creates salt marches and swamps in these part of the plain (Asar 2019). Furthermore, the evaluation of 128 geoelectric sounding profiles in this study area confirms the existence of salt water along fine-grained sediments in the northwest of this area (Anar region) and uplifted bedrock (Miocene gypsum and salt sediment) that caused shallow water table depth.
In the southeast part of the plain toward Sarcheshmeh Mountain, deep coarse-grained sediments that are suitable for composing the aquifer are identified. In the southeast toward the Kabootarkhan region, shallow sediments forming The entrance of groundwater to the plain is located in the southeast from Kerman and Bardsir aquifer, and the groundwater flows from the southeast toward the northwest (Fig. 3a). The analyses of piezometric data (68 wells) indicate that the main recharge areas of the aquifer are located in the northeast and the southeast of the plain. The depth of groundwater ranges between 4.6 m and 177 m. This reaches 140 m in the margin of the plain and 60 to 40 m in the central parts, and then decreases to less than 20 m in the central parts toward the northwest of the plain (Fig. 3b).
The analysis of geological profiles of 7 exploration wells and 52 piezometric wells, which have been drilled in different years, shows the aquifer mainly consisting of gravel, silt, sand and clay in different depths and the amount of clay increases toward northeast part of the aquifer (Fig. 4 a to e).
Based on the analysis of sedimentological properties of logs of wells, a permeability map was produced. This map indicates very high to high permeability rates in the south, southwest and southeast of the plain due to the existence of coarse-grained sediment, and moderate to very low permeability rates in the north, northwest and northeast of the plain due to the existence of fine-grained sediment (Fig. 5).

Data gathering and analyses
In this study, the groundwater samples were collected from 55 different wells of Rafsanjan plain in the summer of 2020. Figure 2 shows the water sampling locations within the block. One-liter polyethylene bottles were used for sampling. Physicochemical parameters like pH and electrical conductivity (EC) of the water samples were taken by using pH meter and an EC meter, respectively. Total dissolved solids (TDS) were computed by multiplying the EC (µ Sm −1 ) by a factor of 0.64. Table 1 reveals the standard laboratory procedures that were followed to determine the concentration of major ions (Pansu and Gautheyrou 2007).
The suitability of water for drinking purposes was determined by applying standard parameters such as World Health Organization (WHO 2011) and Institute of Standards and Industrial Research of Iran (ISIRI1053) (Water treatment physicochemical Specification of Iran (Tehran) 1997); the water quality index (WQI) and plotting the standard reference diagram such as Schoeller (Schoeller 1977).
In order to produce the spatial distribution of each element in the study area, inverse distance weighting (IDW) method was applied in GIS. Table 3 shows the flocculation of physicochemical parameters including pH, EC, TDS and TH, in Rafsanjan plain, where groundwater samples were collected. Figure 6a shows the contour plots for pH variation in summer 2020 in Rafsanjani plain. The pH values ranged from 5.8 to 8.1 with a mean of 7. The maximum value of pH was observed at the southwestern corner of the situ in Jafarabad Deighe (L 25), and the minimum pH value was situated at the northwestern region in Hojjatabad Anar (L 1). The decrease in pH results of the groundwater samples in the northwestern region of the basin (Anar plain) was caused by the sulfide volcanic rocks of the Sarcheshmeh Mountains (Dehj-Sardoueye belt), which is the main source of groundwater recharge in the Anar plain. Then, the water with low pH values was transferred to the end of the basin based on the flow gradient direction (Khajehpour and Abbasnejad, 2007;Malakootian & Khashi 2014). On the other hand, the presence of copper mines in the way of the groundwater slope toward the Anar plain caused a decrease in the pH values in the northwestern areas (Dehghan and Abbasnejad, 2011). The lowering pH of water and the high concentration of heavy elements of arsenic, lead and cadmium, which was affected by the sulfide veins of copper mines in this region, had also been reported in previous studies (Dehghan and Abbasnejad, 2011). Metal mines (such as gold, silver and copper) are often rich in sulfide minerals. The precipitation of sulfides exposed to weather leads to the formation of sulfuric acid. Sulfuric acid dissolves heavy metals and harmful metalloids (such as arsenic) from minerals and forms an acidic pH solution with high concentrations of arsenic, cadmium, lead, copper, etc. (Iatan, 2021). The research conducted by Honarmand et al., (2017) showed that the volcanic rocks of the Sarcheshmeh Mountains and the mineral processing of Sarcheshmeh copper factory led to a decrease in pH values an high amount of heavy elements in the water samples of the studied area (Honarmand et al. 2017). The effect of copper mine-induced pH decrease was also reported in Falun, central Sweden by Ek et al., 2001 (Ek Sonon et al. (2007) showed that the presence of copper in groundwater samples of the Georgia region can lead to an increase in water acidity (Sonon et al. al., 2007).

Variations of pH
The oxidation of dark Jurassic shales rich in pyrite in the recharge points of the Nough plain affected the groundwater pH and led to the acidity water. In the outlet of this region, the presence of sodium alkaline cations affected by the evaporative layers caused the pH levels to 7.6 (Hassanzadeh and Abbasnejad, 2019).

Variation of EC and TDS
The variation of EC and TDS is presented in the contour plots of Figs. 6b and 6c. The mean value of EC and TDS is 8041.4 µS/m and 5436.6 mg/l, respectively. Figure 6b and 6c clearly shows that the low values of EC and TDS were for the water of southeastern corner of the region and increased toward the northwestern zone. There was a close similarity of the variation of EC and TDS with the subsurface water flow from southeastern to northwestern zone in the situ. The maximum of EC and TDS was referred to Ishaqabad (L 7) in the northwestern zone, with the values of 29,300 µS/m and 19,045 mg/l, respectively.
According to Fig. 4 and Fig. 5, the soil texture changed from coarse-grained sediment in the south and southwest, to the fine-grained sediment in north and northwest; this caused different permeability rates and water movement in this plain. The decrease in the speed of water flow in the salt and swampy corridors toward the northwest area increased the concentration of soluble salts (Asar 2019). In the research conducted by Dehghan and Abbasnejad (2011), it was also reported that the main factor in high concentration of soluble salts of groundwater in the    (2018), it was reported that the geological features in the northwestern parts of Rafsanjan plain had caused the water of many wells unusable for agricultural purposes due to excessive salts (Mirnezami et al. 2018).
On the other hand, the increase in the number of authorized and unauthorized wells and over-exploitation in the northwestern part of the plain had led to an increase in electrical conductivity (EC) (Salehi et al. 2022). The overexploitation of groundwater in areas with limestone and salt deposit layers (such as pistachio farming areas in Anar plain) had a significant effect on decreasing water quality (Rahnama andZamzam 2013, Mirnezami et al., 2018).
In the northwest of the aquifer toward the Nough plain, the presence of shallow sediments of fine-grained and clay sedimentary rocks had also reduced the water flow. The low water flow through the limestone formation, gypsum and salt deposit layers caused the dissolution of more salts and the high TDS values (Hassanzadeh and Abbasnejad, 2019). Therefore, the different geological and hydrogeological attributes led to a significant difference in the EC and TDS values in the various parts of Rafsanjan plain. Variation of TH Figure 6d shows the contour plot of total hardness (TH) of the study area. The TH varied from 70 mg/l in Jafarabad Deighe (L 25) to 10,600 mg/l in Ishaqabad (L 7), with an average of 1892.82 mg/l. From the maps, it is clear that for the entire studied area the low values of TH were for the water samples of the southeastern corner of the region and exceeded rapidly toward the northwestern zone. Like EC and TDS, the TH variations were on the same pattern of groundwater flow and were affected by the evaporative deposition subsurface layers. Table 4 presents the concentration and variation of major cations and anions like calcium (Ca 2+ ), magnesium (Mg 2+ ), sodium (Na + ) and chlorine (Cl − ), sulfate (SO 4 −2 ) and bicarbonate (HCO 3 − ), respectively. The maximum concentration of major cations was plotted for the northwestern areas of the plain (Figs. 7a, 4b, 4c), with 2200 (mg/l) for calcium in Aliabad Golshan Anar (L 8) 1296 (mg/l) for magnesium in Ishaqabad (L 7) and 5520 for sodium in Ishaqabad (L 7). The mean concentrations of major cation were 380.3, 226.1 and 1434.9 mg/l for calcium, magnesium and sodium, respectively. The highest value of anions, i.e., chlorine, was observed at the northwestern of the area (Fig. 7e). However, the maximum concentrations of two other anions were reported 1159 mg/l at Zarrin Dasht (L 44) and 3456 mg/l at Abbasabad Javadiyeh Fala (L 14) for bicarbonate and sulfate, respectively.

Variation of major ions
The previous research showed that concentration of cations and anions had a high correlation with the EC and TDS values (Hosseinifard et al., 2015;Aminiyan et al. 2016). Therefore, based on the results in "Variation of EC and TDS" part, the high concentration of ions in the northwestern areas (Anar plain) was strongly affected by the geological and hydrogeological features, i.e., fine soil texture, low permeability and water flow through gypsum layers and salt deposits. The reason for high concentration of ions in the Nough plain had explained by the previous research. Hassanzadeh and Abbasnejad (2019) reported that the low-speed water flow through the evaporative sediments caused an increase in the dissolution of soluble salts, especially chlorine and sodium concentration, in the northwestern areas of Nough plain (Hassanzadeh and Abbasnejad, 2019).

Drinking water standards
The physicochemical parameters of groundwater samples had been compared with drinking water standard that recommended by the WHO (WHO 2011) and ISIRI1053 (Water treatment physicochemical Specification of Iran (Tehran) 1997) for safe drinking usages. The ranges of pH values in groundwater samples, with mean value of 7.02, indicated that all samples were within acceptable limits ( Table 5). The EC values of all the groundwater samples were considerably higher than the acceptable limits in comparison with both drinking water standards (Fig. 8). Another significant factor for assessing the suitability of drinking water is total dissolved solids (TDS). According to WHO and ISIRI1053, TDS values were higher than 1500 mg/L resulted in an unpleasant taste. As, 50 samples were above acceptable limit (about 88%- Fig. 8) while only 5 samples were categorized in suitable class for drinking purposes (Table 5). Based on Sawyer and McCartly (1967), the total hardness (TH) of groundwater samples is classified into four classes: soft (75 mg/L), moderately hard (75-150), hard (150-300) and extremely hard (> 300). Accordingly, about 2%, 3%, 7% and 88% of the groundwater samples in the present study area were classified as soft, moderately hard, hard and very hard, respectively.
The suitability of water is generally evaluated by the concentration of major ions. Table 5 presents the comparison analysis for the concentration of major cations and anions in water samples and the drinking water standards acceptable limits. However, among all reported ions, the sodium and chlorine concentrations in most of water samples are at the critical limit and of serious concern for drinking purpose. They may be due to the source of the evaporative deposition layer in the present study area (Fig. 8).

Water quality index
The water quality index (WQI) is an effective tool for assessing the suitability of groundwater for drinking purposes. In addition, it is a worldwide method to evaluate the impact of human activities on natural groundwater quality (Vasanthavigar et al. 2010, Logeshkumaran et al. 2015. The calculation of the WQI was done by considering weights that range from 1 to 5 for each parameters, based on its relative importance (Logeshkumaran et al. 2015). Table 6 presents the WQI weights and relative weight of physiochemical parameters. According to Arota et al. (2022), the WQI is calculated in four steps (Eqs. 1-4): where Wi is the relative weight, wi is the weight of each parameter and n is the number of parameters. In the second step, the quality rating scale for each parameter is calculated by dividing its concentration in each water sample by where qi is the quality rating, Ci is the concentration of each parameter in each groundwater sample in (mg/L) and Si is the (WHO 2011) standard for each parameter in mg/L. In the third step, the subindex (SIi) is determined for each parameter (Eq. 3). Therefore, in the final step the sum of SI values gives the water quality index for each groundwater sample (Eq. 4): ( where SIi is the subindex of ith parameter, Wi is the relative weight, Qi is the rating based on the concentration of ith parameter and WQI is the water quality index of each sample. Drinking water quality is classified into five classes as excellent, good, poor, very poor water quality and unsuitable for drinking, based on WQI values (Palmajumder et al. (   (Table 7). 3.6% of the samples were classified as poor water quality based on two standards, and 1.8% and 3.6% of the samples were categorized in very poor water quality class according to WHO and ISIRI1053 standards, respectively. The rest of the samples were classified as unsuitable for drinking (WQI > 100) category, with 92.7% and 87.3% for WHO and ISIRI1053 standards, respectively (Table 7). Figure 9 shows the spatial distribution of the WQI for both standards. The samples with poor and very poor water quality were situated in the southeast and some parts of the west and northwest of the region (Fig. 9).

Schoeller diagram
Schoeller diagram (1977) is the most popular method to classify the suitability of drinking water quality (Alavi et al. 2016, Almodaresi et al. 2019, Panahi et al. 2021). The diagram is based on the concentration of major cations (Ca 2+ , Mg 2+ , Na + ), anions (Cl − , HCO 3− , SO 4 2+ ), total dissolved solids (TDS) and total hardness (TH). According to this diagram, the quality of drinking water is classified into six zones including good, acceptable, inappropriate, bad, temporarily drinkable in emergency and undrinkable (Choramin, Safaei et al. 2015). Figure 10 shows that most of the water samples fell in the limits of 'bad' and 'temporarily drinkable in emergency' zones. Water samples with good quality were located in the southeast corner of the region and  (Fig. 10c).

Assessment of water quality for agriculture
The suitability of water for agriculture purposes has been investigated by different studies (Choramin et al. 2015, Palmajumder et al. 2021, Arota et al. 2022, Singh et al. 2022. EC, SAR, Na%, MAR, KR, RSC, PI and Cl − were the most popular parameters to evaluate the suitability of groundwater quality for irrigation usage (Tripathi et al. 2012). Table 8 reveals the classification of groundwater samples according to these parameters. In addition, a graphical representation of the suitability of water samples for irrigation purposes was plotted based on the US salinity, Wilcox and Piper diagrams.

Salinity hazard (EC)
Irrigation with saline water was identified as the primary factor reducing plant growth and crop production (Baghalian, Haghiry et al. 2008). Electrical conductivity (EC) is an important factor for classifying the salinity hazard and the suitability of irrigation water. According to Wilcox (1955), the salinity of irrigation water was classified into five classes based on EC values (Table 8). Analyses of water samples showed that one groundwater sample was under good class, and 4 and 3 samples were in the permissible and doubtful class, respectively. 85% (n = 47) of groundwater samples fell under the unsuitable class, which covered most parts of studied area (Table 8). These results showing samples with low salinity levels were located in the east corner of the plain, and it increases from the center toward the northwest of the Rafsanjan plain.

Sodium absorption ratio (SAR)
In Table 8, the alkalinity hazard of irrigation water was defined by SAR parameter. Following these classifications, about 25% and 49% of groundwater samples were under excellent and good classes that were situated in the east and southeast parts. The alkalinity hazard increased, i.e., permissible (15% of the samples) and doubtful (11% of the samples) categories, toward the northwest. Since SAR reflects the effect of sodium relative to calcium and magnesium, high concentrations of sodium in the central and northeast areas related to the existence of evaporative deposition layers in this area.

Sodium percentage (Na%)
The high sodium content destroys the soil structure, decreases soil permeability and damages crop production (Prihar et al. 1985). Thus, it is introduced as essential parameter to assess irrigation quality (Singh et al. 2022). The sodium percentage (Na%) of collected samples was classified into five categories recommended by Wilcox (1955) guideline. Table 8 shows most of groundwater samples were under permissible (16%) and doubtful (34%) classes.

Magnesium adsorption ratio (MAR)
According to Paliwal (1972), MAR is classified into two classes, i.e., MAR > 50% are suitable and MAR < 50% are unsuitable for agricultural purposes. Based on this category, it can be concluded that about 56% (n = 32) were not suitable and safe for irrigation, as shown in Table 8.

Kelly's ratio (KR)
In 1957, Kelly represented a parameter that KR value < 1 meant suitable irrigation water, while those values greater than this ratio were considered unsuitable (Kelly 1957). Thus, out of 55 water samples in the present studied area, about 93% of samples were unsuitable for irrigation use (Table 8).

Residual sodium carbonate (RSC)
The residual sodium carbonate is a parameter to evaluate the ability of irrigation water to reduce free calcium and magnesium in the soil (Sagar 2015). The RSC value < 0 indicates a little risk of sodium accumulation due to offsetting levels of calcium and magnesium, while positive values indicate sodium accumulation in the soil due to trapping free calcium and magnesium by the bicarbonate and carbonate (Raju 2007). Based on concluded samples in the studied area, except in five locations, about 91% (n = 50) of the groundwater samples had negative RSC value with a low risk of sodium accumulation ( Table 8). The RSC geospatial map shows that the application of the groundwater for irrigation use has little impact on the deposition of free calcium and magnesium in almost all locations (Table 8).

Permeability Index (PI)
The long-term effect of irrigation water usage with high levels of Na + and HCO 3 − on soil permeability has been evaluated by Doneen (1964). In this region, PI varies from 30 to 98 with a mean value of 67 ± 12. Table 8 shows that 31% and 69% of samples come under good to suitable categories for irrigation, respectively. The good quality water samples were situated in the southwest and western areas, while the other parts of the studied area fell under the suitable class (Table 8).

Specific ion toxicity
The cause of toxicity problems in plants is high concentrations of certain ions, such as sodium and chloride, in water or soil. According to plant type and sensitivity, the significant uptake and accumulation of these ions in the plants lead to crop damage or yield reduction (Simsek and Gunduz 2007). However, the sensitivity of the permanent crops to this type of toxicity is more than the annual crops (Ayers and Westcot 1985). Therefore, since Rafsanjan plain is Iranian's foremost producer of pistachio, the evaluation of ion toxicity in water for irrigation purposes is essential.

Sodium
As sodium affects soil physical properties and plant survival, the assessment of sodium hazard is expressed separately. Table 8 presents the classification of irrigation water quality based on its SAR value, which generally defines as sodium hazard. Except for some regions in the east and southeast (n = 12), about 78% of water samples were under severe sodium conditions and were unsuitable for irrigation use and causing more problems such as the appearance of toxicity symptoms on the plant, damage to crop quality and productivity, and reduced permeability and water infiltration (Ayers andWestcot 1985, Flowers, Munns et al. 2015).

Chloride
Chloride is another ion that, at high levels, can cause toxicity problems for plants. The chloride toxicity symptoms are leaf burns, leaf tissue deaths and early leaf drop (Ayers and Westcot 1985). The classification of irrigation water based on its chloride value is shown in Table 8. Although waters with chloride values below 140 mg/l are suitable for irrigation use, about 50% of water samples in the present study area had chloride values > 350 mg/l and severe hazard conditions for irrigation purposes.

Piper diagram
Piper in 1944 proposed a trilinear diagram to determine the geochemical evolution of groundwater. This diagram represents the chemistry of water samples. According to this procedure, the cations, i.e., calcium, magnesium and sodium plus potassium, and anions, i.e., sulfate, chloride and carbonate plus hydrogen carbonate, are presented in two distinct triangle plots. The result is then projected onto a diamond, which is a matrix transformation of a graph of the anions (sulfate + chloride/total anions) and cations (sodium + potassium/total cations) (Piper 1944). Figure 11a shows the Piper diagram of the water samples. It indicates that about 93% of the samples fell in sodium-chloride-sulfate facies, while 5% and 2% of the samples fell within calcium-magnesium-sulfate-chloride and sodium-bicarbonate facies, respectively (Fig. 11a). In addition, in more samples sodium and chloride were as dominant cations and anions with a frequency of 98% and 93%, respectively (Fig. 11b). Thus, most of the locations in the studied area fell under sodium-chloride-sulfate facies due to the presence of subsurface evaporative layers (Fig. 11c). It denotes that almost whole water samples in this region were alkaline water type with prevailing chloride (Ravikumar, Somashekar et al. 2015).

US salinity diagram (USSL 1954)
The US salinity diagram (USSL) is another hydrochemistry plot to evaluate the water suitability samples for irrigation purposes using EC and SAR values (US Salinity Laboratory (USSL) 1954). According to this diagram, groundwater quality was classified into sixteen classes (Table 9).
Based on the plot of groundwater data (Fig. 12a), most of the water samples in the study area have high EC and SAR values. Table 10 presents that 65.45% of samples are C4S4 type and 23.64% are C4S3 type with very highsalinity-high-sodium and very high-salinity-high-sodium hazards, respectively. The rest of samples are classified as C4S2 (1.82%), C3S3 (1.82%), C3S2 (5.45) and C2S1 (1.82%). However, using very passion groundwater type for agricultural use causes different levels of soil salinity and sodicity hazards that applying salt leaching and reclamation techniques can be effective in improving plant growth and crop productivity (Bastani and Hosseininia 2018). Wilcox (1955) represented a diagram (using EC and Na%) with five classes including: excellent to good, good to permissible, permissible to doubtful, doubtful to unsuitable and unsuitable, in order to assess the water suitability for irrigation uses (Fig. 12b). This diagram reveals four, four and one groundwater samples in permissible to doubtful, doubtful to unsuitable and unsuitable classes, respectively. The remaining water samples had EC and TDS of more than 3500 µS/m and 35 mg/l, respectively, that they were unsuitable for irrigation purposes. The spatial distribution map shows that except for a few locations at the southeast corner of the studied area, all other parts were categorized as unsuitable (out of the Wilcox diagram) with high salinity values and sodium percentages ( Fig. 12b and Fig. 13).

Conclusions
The assessment of the hydrogeological status of the study area in Rafsanjan plain, Iran was performed based on different approaches, namely drinking water standards, water quality index (WQI) and Schoeller diagram for domestic usage and irrigation water quality (IWQ) parameters, Piper diagram, US salinity diagram and Wilcox diagram for irrigation purpose. Accordingly, the main findings were categorized below: • The average pH and EC values of groundwater were 7.02 and 8041.4 µS/m, respectively, which indicated neutral and saline water status.   • The GIS-based maps showed that variation of EC, TDS, TH and major ions in groundwater was strongly affected by groundwater flow and the geological structure of the studied area. A considerable increase in the mentioned parameters value was noted in the locations in the north and northeast of the studied area. Passing water flow through the evaporative deposition layers is the main reason for the rapid increase in the concentration of ions. • According to WHO and ISIRI1053 standards, the concentrations of cations and ions and the values of EC, TDS and TH were significantly higher than acceptable limits and were critical for drinking purposes. • WQI values classified the groundwater into good (1.8%), poor (3.6%), very poor (1.8%) and not suitable (92.7%) for drinking uses, based on WHO standards. According to ISIRI1053, the classification of groundwater for drinking indicated that 5.5%, 3.6%, 3.6% and 87.3%, of the samples in good, poor, very poor and not suitable classes, respectively. • The Schoeller diagram illustrated that most of the water samples fell within the limits of 'bad' and 'in emergency' zones. The spatial distribution map demonstrated that the groundwater was rated as acceptable in a few locations in the southeast corner of the region; at the same time, the water quality drastically decreased toward the north and northwest area of the plain. • IWQ parameters expressed that 85%, 67%, 32%, 51%, 43% and 50% of samples had EC > 3000, Na% > 60%, MAR > 50%, KR > 1, SAR > 9 and Cl − > 350, respectively, which were unsuitable for irrigation uses. According to RSC and PI, about 91% of the groundwater samples had negative RSC value with little risk of sodium accumulation, and 69% of the samples belonged to suitable permeability conditions for long-term irrigation usage. • The groundwater types were NaClSO4 (93%), CaMg-SO4Cl (5%) and NaHCO3(2%), which denoted alkaline water type with prevailing chloride in most sites. • Plots on US salinity diagram revealed that 63% of the water samples fell in very high-salinity-high-sodium hazard zones, while there were limited samples in the low-salinity or sodicity class. Thus, conducting some practices in order to reduce soil salinity and sodicity hazards can be effective in improving the conditions of irrigated lands and crop productivity. • According to the Wilcox diagram, only 7% of groundwater samples were permissible to doubtful for irrigation, and the rest were doubtful (7%) to unsuitable (85%) classes, which existed in the central to northeast regions. • Based on the results, it was concluded that almost all the groundwater samples were unsuitable for drinking and irrigation usage. Thus, some practices are recommended to reduce groundwater salinity in this area, such as: limit-ing over-exploitation of groundwater resources, controlling agricultural development and practices, using soil and water amendments, mixing high-salinity water with stored rainwater and applying modern irrigation methods and organic fertilizer (Etikala, Adimalla et al. 2021).