Impact of climate change on hydrochemical properties and quality of groundwater for domestic and irrigation purposes in arid environment: a case study of Al-Baha region, Saudi Arabia

Groundwater over-exploitation, rising temperature and declined precipitation, have resulted in groundwater quality deterioration, subsequently affecting human health and agricultural productivity. Therefore, to evaluate the suitability of groundwater for domestic and irrigation purposes, groundwater samples were collected from 88 locations in Sarat Al-Baha region, Saudi Arabia. The hydrochemical characteristics were assessed to evaluate the suitability of groundwater for drinking and irrigation purposes. The results revealed that most groundwater samples were medium to highly salinity hazards; nevertheless, no sodicity hazards were expected. By using the water quality indices, 70.5% of studied waters were class II. Nevertheless, 91% of collected samples were in acceptable concentrations of heavy metals and nitrate. The saturation index showed that the groundwater was largely under-saturated with most minerals, indicating a possibility of salinity increase due to increasing the concentration of iron, calcium, magnesium, sodium, chloride and sulfate due to the dissolution of under-saturated minerals. Due to the increased salinity levels, more than 90% of the collected groundwater samples showed low quality for irrigation due to higher potential for salinity problems. Furthermore, the development of salinity of groundwater in Sarat Al-Baha could be the result of the increased temperatures and decreased rainfall. Therefore, it is recommended to implement appropriate management practices to optimize groundwater usage and decrease the potential of further deterioration in its quality. The discharge of groundwater need to be controlled, strictly preventing the drilling of new wells. It is also important to assess the suitability of planted crops, and implement proper water management practices to sustain groundwater resources in Sarat Al-Baha region. The approach used in this study can be used to alike environment worldwide.


Introduction
Elevated temperature and declined precipitation rate ultimately reduce the recharge rates of groundwater, subsequently deteriorating the quality of groundwater particularly in arid regions. Consequently, the groundwater quality and vital ecosystem water requirements are threatened, particularly during water scarcity periods (Treidel et al. 2012;Turral et al. 2011;Wu et al. 2020;Obianyo 2019;Várallyay 1994). As the Kingdom of Saudi Arabia (KSA) is located in arid environment, it is suffering from accelerated groundwater deterioration due to limited precipitation and warming climate (Al-Omran et al. 2016;Aly et al. 2016a, b;Kløvea et al. 2014;Alghamdi et al. 2021).
Climate change affects (both directly and indirectly) groundwater quantity and quality in many complex ways (Alghamdi et al. 2021;Kløvea et al. 2014;NOAA 2017).
For example, Xiaolong and Boufadel (2017) used field measurements and numerical modeling to investigate the effect of evaporation and rainfall rates on subsurface salinity and groundwater dynamics. They found that increased evapotranspiration rates increased pore water salinity in the surface 10 cm layer. In addition, climate change can greatly affect recharge rates, depth of groundwater, and water consumption for human activity and irrigation (Ludwig and Moench 2009). Furthermore, low recharge rates in dry regions increase groundwater depletion, affecting both quantity and quality of groundwater (Treidel et al. 2012).
As an essential constituent of the hydrologic cycle, groundwater is highly vulnerable to climate changes i.e., increased temperature, especially with the rapid changes in terms of increased temperature observed nowadays in comparison to historical patterns. It is expected that these rapid changes will have far-reaching consequences for Saudi's water resources (Al-Barakah et al. 2020).
It is mostly significant to recognize that alarming problems with regard to increased groundwater depletion and groundwater quality deterioration are starting to already exist in Saudi Arabia in response to observed climate change (Aly 2015). Climate change may enlarge or accelerate long standing stresses determined by human population growing, temperature increase and diminishing groundwater reservoir. These stresses are more severe in arid and semi-arid areas where aquifer recharge rates are low, and the demand for groundwater for irrigation purposes is higher (Turral et al. 2011). Therefore, temporal monitoring of groundwater quality is of critical importance in changing climate under arid climatic conditions to avoid potential hazardous impacts on human health and agricultural production. Groundwater recharge depends on precipitation, evapotranspiration, and land cover. Warmer summer temperatures can reduce the amount of groundwater due to the increase of evaporation and more groundwater abstraction (Jyrkama and Sykes 2007). In addition, aquifer type and properties affect groundwater recharge. Shallow and unconfined aquifers react promptly to elevated temperature; however, confined aquifers respond slowly (Lee et al. 2006). Therefore, deep and confined aquifers will be less delicate to the direct impact of temperature (Wada et al. 2012). During groundwater recharge, an increase in the risk of contaminants leaching is anticipated (Okkonen et al. 2010). Although groundwater can be contaminated by natural sources; however, human activity is usually the main causes of groundwater contamination worldwide (Okkonen et al. 2010). The high population density, usually led to high groundwater vulnerability to contamination. This is due to the human activity accompanied by chemicals and wastes released to the environment; consequently, to the groundwater (Alghamdi et al. 2021;Okkonen et al. 2010). When groundwater contaminated, its restoration and clean up is not easy and not affordable (Alghamdi et al. 2021). Numerous studies presented the direct relationship between agriculture and nitrate concentration in the groundwater (Al-Omran et al. 2016;Aly et al. 2016a, b;Hansen et al. 2012). Aly et al. (2016a, b) and Bloomfield et al. (2006) reported that the excessive use of nitrate fertilizers and pesticide in agricultural areas may led to their leaching into groundwater. MacQuarrie et al. (2001) said that the septic tanks can also contribute to groundwater nitrate pollution. Al-Barakah et al. (2020) demonstrate that the wells near urban area in Al-Baha region, KSA were found contaminated by E. coli and nitrate due to leaching of watewater from septic tanks into groundwater.
The groundwater hydrochemical parameters comprising anions, cations, nitrates, heavy metals, pH, dissolve salts, … etc (El-Sayed et al. 2012). These parameters are requested to determine the water quality indices (WQI) and spatiotemporal monitoring of the groundwater (Al- Barakah et al. 2017). The WQI is an efficient tool to assess water quality and aid policy-makers since it integrates water parameters into one number representing the water suitability and quality (Al-Barakah et al. 2020). The WQI have been properly calculated and used by different investigator such as Backman et al. (1998) which used the WQI for assessing groundwater quality in Finland, Ketata-Rokbani et al. (2011) applied the WQI to assess the groundwater quality in Tunisia, Al-Omran et al. (2015) used the WQI to assess water quality in Riyadh, KSA, and finally, Aly et al. (2015) applied the WQI in the groundwater of the city of Hafar Albatin, KSA.
Water management decisions regarding groundwater use a dominant role in dictating the extent of climate impact on state of aquifers and streams. If a considerable increase in the groundwater net volume withdrawal is used to address global warming problems, there are likely to be severe and detrimental effects on groundwater storage and quality. In KSA, groundwater resources are considered the main source of irrigation requirements for plants (Al-Omran et al. 2016). However, increased groundwater depletion and deterioration in the quality of groundwater in KSA have been reported over the past two decades owing to extreme temperatures (Aly et al. 2016a, b;Aly 2015), which necessitates increased effort to monitor and assess the quality of groundwater in KSA . Al-Baha region, located in the southwest of KSA, is the smallest provinces of KSA and due to the arid conditions, the groundwater in Al-Baha is considered a very valuable natural resource. The Al-Baha region is covered by crystalline basement rocks of the Arabian shield. The crystalline basement has overlain by basalts, volcanic rocks, in some areas. The groundwater occurs within the weathered and fractured crystalline rocks, igneous and metamorphosed igneous rocks (Al-Barakah et al. 2020). The water flow in the direction of lower areas namely called wadis (where there are many wadis in the region). The wadis are mostly consists of the unconsolidated wadi deposits/ alluvial deposits (Mahmoud and Alazba 2015). In Al-Baha regions and alike environment, the contaminants move quickly through fractures in rocks into the groundwater. Rock fracture is considered problem with respect to contaminations due to the fractures are usually spaced randomly and do not follow hydraulic gradient or the land surface contours line (Mahmoud and Alazba 2015).
Until recently, Al-Baha region was considered as one of the least developed regions in Saudi Arabia. However, over the past few years, there has been an increased interest in the socioeconomic development of the region, particularly as Al-Baha enjoys a natural variation of terrestrial ecosystems, moderate climate and natural resources. However, AL-Baha region is facing a continuous increase in temperature, while significant decline in rate of precipitation. These climatic variations are likely to influence the quality and volume of groundwater to a substantial level. Al-Baha region was characterized by an increase in mean temperatures by 1.1-1.6 °C, and a decrease in annual rainfall by 24-41% during the years 1995-2019. However, very limited studies are available for the assessment of groundwater quality for irrigation and drinking purposes in this region under rapidly changing climatic conditions (Al-Barakah et al. 2020). Therefore, current study was conducted to (1) evaluate the suitability of groundwater in Al-Baha region; 2) evaluate the sutability of groundwater in Al-Baha region for domestic and irrigation purposes using water quality index (WQI) analyses; and (2) characterize the hydrochemical composition and mineral saturation indices of groundwater in Al-Baha region under changing climatic conditions.

Study area location, geology, and GIS analyses
Al-Baha region is located between latitudes 19° 27′ 17″ and 20° 49′ 75″ N and longitudes 40° 46′ 30″ and 42° 10′ 10″ E (Fig. 1). The total area of Al-Baha region is 11221 km 2 . Numerous geological features are recognized in the KSA (Fig. 1). Al-Baha region is located in the Arabian Shield in the western part of KSA. The Arabian Shield contains crystalline basement covered by the basalts (Al- Barakah et al. 2020). In this study, the major physical, chemical, and hydrological properties of the soil and groundwater samples were geo-referenced to its geographic locations in ArcGIS 10.4 (ESRI 2015) to draw thematic maps of the spatial distribution of soil and groundwater properties in Sarat Al-Baha region. The geostatistical analyst function within ArcGIS was used to interpolate values between measurement locations. A 30-m digital elevation model (DEM) was generated for the study area. The Spatial analyst function was used to calculate main flow directions based on the generated DEM. Most of the flow directions in study area were found in the north, west, and east directions, with approximately 46% of all flow in these three directions ( Fig. 1) (Mahmoud and Alazba 2015). The Al-Baha has six administrative governorates (Al-Aqiq, Belgrashi, Al-Mandq, Al-Qora, Al-Mkhwah, and Qolwah) in addition to Al-Bahaa city (Fig. 2). Al-Baha region is characterized by diverse soils, varied plant cover, and different patterns of topography and climate. Since Al-Baha region is located on the Arabian Shield, its soils are mostly consisting of crystalline basement composed of Precambrian continental crust (Mohammed et al. 2011). The weather is generally moderate to hot in summer and warm to cold in winter, with average temperatures ranging between 12 and 23 °C. The average annual precipitation rate is about 150-200 mm, with moisture between 50 and 70% (Al- Barakah et al. 2020).
In this study, we focused our investigation on the area representing the middle part of Al-Baha region, which is known as Sarat Al-Baha, and is consisting of Al-Qora, Al-Mandq, Al-Baha, and Belgrashi Governorates (Fig. 2). The area of Sarat Al-Baha represents only 26.7% of the total area of Al-Baha region. Nevertheless, Sarat Al-Baha has the largest population percentage in the entire Al-Baha region. In addition, Sarat Al-Baha area includes the largest cultivated area, forest and mountainous areas, and, therefore, is a major area of interest for future socioeconomic development in Al-Baha region. With groundwater being a major resource for drinking and irrigation in Al-Baha region, it is, therefore, of crucial importance to assess its quality as impacted by climatic changes in the region. We obtained historical daily climatological reports for Al-Baha region from the General Authority for Meteorology and Environmental Protection in KSA. The data included average daily temperature and precipitation for 25 years  in Al-Baha region.

Groundwater sampling and analysis
Groundwater samples were collected from 88 locations in 1.0 L air tight plastic bottles during August 2020 (Fig. 3). Five replicte samples were collected from same well, which were placed in ice-boxes and moved to the laboratories in the Department of Soil Sciences at King Saud University, Riyadh for analyses. A portion of the groundwater sample was acidified by concentrated HNO 3 to a pH < 2 for accurate determination of heavy metal concentration (APHA 1998).
The samples were analyzed for pH using pH meter (pH meter-CG 817), and total dissolved salts (TDS) by measuring electrical conductivity (EC) in dS m −1 at 25 °C (Test Kit Model 1500_20 Cole and Parmer). Magnesium and calcium ions were determined by titration by versenate as described by Sparks (1996). The potassium and sodium ions were measured by flame photometer (Corning 400) (Matiti 2004). Bicarbonates and carbonates ions were determined by acid titration (Sparks 1996). Chloride ion was determined by titration with silver nitrate (Matiti 2004). Sulfate was measured by the turbidity method (Tabatabai 1996), while the phenoldisulfonic acid method was used for the determination of nitrate (NO 3 − ) (APHA 1998). Heavy metal concentrations Blanks and standard solutions examination were frequently run to check for probable inaccuracies in the analysis. The accuracy of the chemical analysis was estimated  (Appelo and Postma 1996). IBE values equal to ±3 % were considered satisfactory. Samples with an IBE out of this range were re-estimated. More than 97% of the analyzed water samples were in the acceptable range with regard to the IBE values, indicating the suitability of the obtained results to carry out geochemical modeling and assessment of groundwater quality in the study area.

Calculation of water quality index (WQI)
WQI was computed to assess the quality of groundwater for drinking and irrigation purposes in the study area. Computation of WQI was carried out according to Al-Barakah et al. (2017). Ten measured parameters (i.e.: pH, TDS, Fe, Zn, Cu, Cd, Cr, Pb, As, and NO 3 -N) were included in the calculation of WQI. Each parameter was allocated a weight (wi) conferring to its significance in the quality of groundwater (Table 1). Parameters with the largest significance (i.e.: As, Pb, Cd, and NO 3 -N) were given a weight of 5. Whereas, parameters with the lowest significance (i.e.: Fe, Zn, and Mn) were given a weight of 1 (Ramakrishnalah et al. 2009). The remaining parameters (i.e.: pH, TDS, Cu, and Cr) were given a weight value between 1 and 5 (Table 1). A Relative where RWi is the relative weight of the i th parameter, Wi is the weight assigned to i th parameter, and n is the total number of measured parameters. RWi values of measured parameters are presented in Table 1.
A Quality Rating Scale (QRSi) was calculated for each measured parameter (i) according to the equation: where Ci is the measured value of the ith parameter (mg L −1 ), and Si is the WHO acceptable standard limit of the ith parameter (Alexakis 2011;WHO 2011).
A parameter index value (PIi) was calculated by multiplying the relative weight (RWi) by the quality rating scale (QRSi) for each measured parameter. The final value of the WQI was calculated as the sum of all computed parameter index values (PIi) according to the equation: where PIi is the parameter index value of the ith parameter, and n is the total number of measured parameters. According to Ramakrishnalah et al. (2009), the quality of groundwater, for drinking purposes, based on the values of the WQI is classified into five classes: excellent, good, poor, very poor, and unsuitable for drinking, with value ranges of <50, 50.1-100, 100.1-200, 200.1-300, and >300, respectively.

Geochemical modeling
Interactions amongst the groundwater and the surrounding rocks are the main processes controlling the chemical features of groundwater. The water's deviation from equilibrium with respect to dissolved minerals is defined by the saturation index (SI). The SI of a mineral is calculated using the following equation: PIi where IAP is the ion activity product of the dissociated chemical species in solution and k t is the equilibrium solubility product of the chemical involved (Alexakis 2011). The PHREEQC model (Ramakrishnalah et al. 2009) was used to obtain the SI of groundwater samples in Sarat Al-Baha region with respect to the main mineral phases. Saturation indices (SI) are usually used to assess mineral dissolution in groundwater. If the SI of the mineral is equal to zero, the mineral is in equilibrium with groundwater. SI values smaller than zero indicate that the mineral is under-saturated, whereas SI values greater than zero indicate the mineral is over-saturated.

Historical climatic conditions in Sarat Al-Baha
Records of climatic data over the past 25 years  in Al-Baha region revealed mean temperature periodic cycles approximately every 2-3 years, whereas, annual rainfall showed slightly longer periodic cycles of 4-5 years (Fig. 4). During the period 1995-2019, average temperature in Al-Baha region was 23.2 °C. However, a clear evidence of the impact of climate change on mean annual temperatures in the region was the observed gradual increase in mean temperatures by 1.1-1.6 °C that was noticed between the years 1995 and 2019 (Fig. 4A). In contrast, there was a gradual decrease in total annual rainfall that ranged between 24 and 41% during the same period (excluding the extremely dry period that took place from 2010 to 2014) (Fig. 4B). Previously, Al-Baha region, which indicates that Al-Baha is a fragile ecosystem, highly vulnerable to climate change and subsequent environmentally, induced water resource degradation (Al-Barakah et al. 2020). Improper water management will further increase the adverse impacts of climate change, and will accelerate problems associated with the degradation of groundwater quality in Al-Baha region (Aly et al. 2016a, b).

Characterization of groundwater resources in Sarat Al-Baha
The digital elevation model (DEM) of Sarat Al-Baha region showed considerable variations in topography, with an average elevation of 1882 m above sea level (Fig. 5A). Areas in the northern parts of Al-Baha governorate and in the southern parts of Al-Qora and Al-Mandaq governorates showed the highest elevation that reached 2446 m above sea level. On the other hand, the northern parts of Belgrashi and Al-Qora governorates showed the minimum surface elevation of 1318 m above sea level. The large variability in elevation in Sarat Al-Baha region emphasizes the unique nature of the hydrogelology in the region. It is, therefore, anticipated that the quantity and quality of groundwater in Al-Baha depend mainly on rainfall rate. Consequently, we can say that the Al-Baha groundwater has a meteoric origin, which means that the water comes from precipitation in the form of rain with a short residence time of water flow through the vadose zone to the aquifers (Craig 1961;Good et al. 2014). The maximum depth to groundwater in Sarat Al-Baha was 190 m, which was observed in Al-Mandaq governorate. However, the average well depth was found to be 48.3 m from the soil surface (Fig. 5B). These results suggest that most groundwater in Sarat Al-Baha are present in shallow, and as such will likely to be more sensitive to be adversely impacted by climate change in the region (Lee et al. 2006;Wada et al. 2012). Most studied groundwaters in Al-Baha region have normal range pH (6.7-8.4); On the other hand, their salinity were moderatly high and ranged bewteen ( <1 50-3000 mg/L) (Fig. 5C, D) (Al-Barakah et al. 2020; Ayers and Westcot 1985). Groundwater in Al-Baha region is considered the main source for drinking, irrigation, and livestock purposes. Therefore, the quality of groundwater in terms of its salinity and mineral content is extremely important to maintain sustainable development in the region. In particular, salinity is highly significant when groundwater is used for irrigation to reduce the hazards of soil salinization (Aly 2020). Approximately, 40% of the tested wells showed total dissolved salts (TDS) values between 500 and 1000 ppm (Table 2), indicating acceptable levels for the cultivation of most crops and low potential for soil salinity problems. However, the majority of the collected groundwater samples showed higher salinity levels ( > 1000 ppm) (Figure 5D), making this groundwater unsuitable for drinking. Furthermore, using groundwater having this level of salinity for irrigation could have adverse impacts on soil quality. Secondary salinization is particularly anticipated under the arid and hot environment of Sarat Al-Baha region. The observed increase in mean temperatures in Al-Baha region, and the accompanied decrease in average annual precipitation over the past 25 years, could have contributed to the increasing rates of salination observed in groundwater. Other reasons could include excessive depletion of groundwater and agricultural expansion in the region (Wu et al. 2020;Aly 2015).

Evaluation of groundwater quality for drinking
To facilitate the assessment of the quality of groundwater in Sarat Al-Baha region, a map showing all sampling locations along with well numbers is presented in Fig. 6. The calculated values of the WQI were interpolated using the geostatistical analyst function inside ArcGIS 10.4 (Fig. 7). Calculations of WQI were used to evaluate the quality of groundwater for drinking. We found that 20.5% of collected groundwater samples were excellent for drinking (class I), 70.5% were good for drinking (class II), 8.0% were poor water for drinking (class III), and approximately 1.0% were very poor for drinking (class IV) (Fig. 7). The area showing very poor groundwater quality for drinking was around well number 75 in Al-Mandq governorate. This result is  1995-1999, 2000-2004, 2005-2009, 2010-2014, and 2015-2019 in agreement with the finding of Al-Barakah et al. (2020) whose concluded that the groundwater near urban area with the absent of sewere system are subjected to E. coli and nitrate contamination due to leaching of watewater from septic tanks (around 5 m in depth) into groundwater. It is worth mentioning that arsenic (As) concentration in groundwater collected from well 75 showed the maximum concentration (0.163 mg L −1 ) in the entire study region. The WHO (2011) recommended standard for As concentration in drinking water is 0.01 mg L −1 (Table 1), which can explain the very poor quality of groundwater obtained from well 75 based on WQI value. Al-Farraj et al. (2013) reported that the main source of As released to the groundwater systems in KSA is nature geothermal and is linked to processes similar to those involved in the boron release. Furthermore, Taleb (2009) said that high levels of As occur due to mixing of cold waters with geothermal fuids and are often associated with a high total dissolved solids (TDS). Our results are slightly different from Al- Barakah et al. (2020) who used WQI to assess groundwater quality for drinking in Al-Baha region in 2016, and concluded that approximately 97% of groundwater samples were class I and class II in quality. The slight deterioration in groundwater quality between the years 2016 and 2020 could be attributed to increased groundwater salinity (Sharaf 2013). This study suggest that the increasing of groundwater salinity in Al-Baha region is due to the increase in mean temperatures and gradual decrease in total annual rainfall. Since the groundwaters in Al-Baha are shallow and located in unconfined aquifers, consequently, react promptly to climate change (Lee et al. 2006). The climate change can affect recharge rates, depth of groundwater, and water consumption for human activity and irrigation (Ludwig and Moench 2009). This can led to natural groundwater salinity increase due to mineral dissociation increase (such as halite mineral) in groundwater aquifer, and groundwater deterioration due to anthropogenic activiies (such as an increase water consumption and leaching of saline drainge water to the groundwater system) (Al-Barakah et al. 2020;Alghamdi et al. 2021).

Evaluation of groundwater quality for irrigation
The method of the Food and Agriculture Organization (FAO) (Ayers and Westcot 1985) was used to evaluate the suitability of Al-Baha groundwater for irrigation. Results indicate that 29.5% of groundwater samples showed acceptable salinity levels (EC < 0.7 dS m −1 ), and good quality for irrigation without any restriction. However, the majority of groundwater samples (68.2%) in Sarat Al-Baha had larger EC values (0.7-3.0 dS m -1 ), and were classified as having slight to moderate degree of restriction for irrigation. Groundwater samples with salinity levels exceeding 3.0 dS m −1 were classified as having severe restriction for irrigation. These samples represented only 2.3% of groundwater samples collected from the study region, and were mainly located in the eastern parts of Al-Baha governorate (e.g. wells 30, 49, 50, 51, and 52), and the western part of Al-Mandq governorate (e.g. well 80). The concentration of heavy metals and nitrate were within the acceptable safe limit for irrigation for all studied wells, with the exception of one sample (well 75 in Al-Mandq governorate), which had As concentration much larger than recommended level (Al-Barakah et al. 2017;Aly et al. 2013). The combined effect of alkalinity and salinity on the quality and classification of groundwater used for irrigation was evaluated using the salinity diagram of the U.S. Salinity Laboratory (Richards 1954). Based on the characteristics of the groundwater samples from Sarat Al-Baha, the sodium hazard was low (class S1) for all groundwater samples, indicating low potential for alkalinity problems (Fig. 8). In contrast, only 4.5% of groundwater samples showed low salinity hazard, giving groundwater in these areas a final classification of C1-S1. These areas were mainly found around wells 74 and 76 in Al-Mandq governorate. Groundwater samples showing medium salinity and low sodium hazards (i.e. class C2-S1) also represented 4.5% of total samples. These results indicate that only 9% of groundwater samples collected from Sarat Al-Baha had low potential for salinity problems, and were considered suitable for the irrigation of most field crops, vegetables and tree plants without the need for special practices to control salinity levels in the soil. The majority of groundwater samples (80.7%) had high salinity and low sodium hazards (i.e. class C3-S1), while approximately 10.3% had very high salinity and low sodium hazards (i.e. class C4-S1), indicating the lowest quality of groundwater for irrigation.
In 2016, the study of Al-Barakah et al. (2020) used the salinity diagram of the U.S. Salinity Laboratory (Ramakrishnalah et al. 2009) to classify groundwater in Al-Baha region, and found that the class C3-S1 (i.e. groundwater having high salinity and low sodium hazards) represented 60% (compared to 80.7% in our study) of all collected groundwater samples. Taking into consideration the variations that exist between the study of Al-Barakah et al. (2020) and this research in terms of the number of groundwater samples and their locations, we still can see a clear trend of increase in groundwater salinity in Al-Baha during the period 2016-2021. The fact that 10.3% of collected groundwater samples had very high salinity (i.e. class C4-S1) can be considered as an alarming warning of quick trend of deterioration in the quality of groundwater in Sarat Al-Baha region (Aly and Benaabidate 2010). Our findings indicate that groundwater in most of the sampling locations (91%) had excessive amount of salts, and hence it cannot be used for irrigation of most crops without special measures for salinity control such as the planting of salt-tolerant plants, proper water management, and the application of additional leaching requirements (Lauchli and Epstein 1990).

Hydrochemical characterization
Results of the chemical analyses of groundwater samples were plotted on a Piper trilinear diagram (Piper 1944) (Fig. 9). Piper trilinear diagram offers a suitable method to classify groundwater types based on the ionic composition (Aly and Benaabidate 2010). The diagram illustrates that all groundwater in Sarat Al-Baha were calcium-magnesium, and chloride-sulfate type water. Therefore, the anticipated predominant groundwater type in Sarat Al-Baha should be Mg/CaSO4, which is commonly due to the geology of the area that contains gypsum, anhydride, and dolomite minerals (Al-Barakah et al. 2020).
The concentration of ions present in groundwater were plotted on Gibbs's diagrams (Fig. 10). The diagram helps to determine main factors affecting the quality of groundwater. Results indicated that the chemical weathering of rock-forming minerals is affecting groundwater quality in Sarat Al-Baha. Most groundwater samples were located in the rock dominance section of the diagram, suggesting that rock-water interaction is the major source of dissolved ions for most groundwater in Sarat Al-Baha region. The rock-water interaction process includes chemical weathering, dissolution-precipitation of secondary carbonates, and ion exchange between water and clay minerals (Subba 2006). The increase in air temperatures and evaporation rates could significantly increases the ion concentrations formed by chemical weathering, leading to a higher salinity (Kumar et al. 2006).

Geochemical modeling
Variations in pH greatly affects the calculated values of SI. Inaccuracy in the determination of pH, and/or variations in pH of groundwater as water flows can result in an inaccurate calculation of SI. Therefore, it is recommended that SI should be taken in a wider range (−1 < SI < +1) to minimize the impact of pH variations on the calculation of SI values (Nasher et al. 2013 (Fig. 11). In contrast, SI values showed that all groundwater samples were over-saturated with respect to Hematite (Fe 2 O 3 ) (as indicated by the positive SI values) (Fig. 12). Our findings indicate that there is a possibility of increasing the concentration of iron, calcium, Fig. 7 Interpolated areal distribution of the Water Quality Index (WQI) in Sarat Al-Baha region magnesium, sodium, chloride and sulfate in groundwater in Sarat Al-Baha due to possible further dissolution of under-saturated minerals (Appelo and Postma 1996;Jalali 2006). This explain why groundwater salinity increases in Al-Baha region. From sake of brevity, it should've said that most of the flow directions in study area were found in the north, west, and east directions toward wadis (Fig. 1). The Al-Baha region, like other regions, located in the Arabian Shield comprises Tertiary-Quaternary basaltic lava flows, Precambrian-Cambrian basement rocks, Quaternary-Recent alluvial deposits, and Cretaceous-Tertiary sedimentary succession (Sharaf 2011). Many studies showed the presence of numerous recent local buried faults (Al-Shaibani 2008; Metwaly et al. 2021a;Metwaly at al. 2021b;Sharaf 2011). The water in the Arabian Shield generally stored in the subsurface wadi deposits and shallow fracture basement (Al-Shaibani 2008). The groundwater flow in Al-Baha is toward the north, west, and east directions, following mostly the surface drainage system. Hydraulic gradient varies greatly from one point to another depending on the pumping rates and cross-sectional area of the aquifer in addition to give a magnitude of the transmissivity (Sharaf 2011). The most important recharge mechanisms for the main sub-basaltic aquifer are the seasonal rainfall and the floods with short period through fractures. However, the main discharge mechanisms (approximately 576 m 3 /day) are the pumping of groundwater for irrigation, and domestic purposes (GAS 2019; Metwaly et al. b;Alghamdi et al. 2022). The shallow aquifer is a saturated thickness of 3-20 m and water is established under semiconfined and unconfined aquifers (Metwaly et al. 2021a, b).

Conclusions
Al-Baha region, a fragile agro-ecosystem, is highly vulnerable to climate change and subsequent environmentally, induced groundwater resource degradation. A clear evidence of mean annual temperatures gradual increase; on the other hand, a gradual decrease in total annual rainfall between the years 1995 and 2019. The groundwater in Sarat Al-Baha is the main resource for domestic and irrigation purposes. To evaluate the quality and hydrochemical characterizations of groundwater, 88 groundwater samples were collected from Sarat Al-Baha. Calculation of water quality indices showed that most groundwater in Sarat Al-Baha had acceptable concentrations of heavy metals and nitrate, and were good for drinking. Only 8% of the collected groundwater samples were classified as having poor and very poor quality for drinking, mainly because of increased concentrations of As and NO 3 -N. Assessing the quality of groundwater for irrigation using the salinity diagram of the U.S. Salinity Laboratory showed that only 9% of the collected groundwater had low (C1-S1) and medium (C2-S1) salinity hazards. The majority of groundwater samples (80.7%) had high salinity and low sodium hazards (i.e. class C3-S1), while approximately 10.3% had very high salinity and low sodium hazards (i.e. class C4-S1), indicating the lowest quality of groundwater for irrigation. Saturation index showed that groundwater samples were over-saturated with Hematite, and under-saturated with respect to other minerals (Halite, Anhydrite, Gypsum, Manganite, Jarosite-K, Hausmannite, Pyrolusite, Pyrochroite, Melanterite, and Sylvite). Therefore, there is a possibility of increasing the concentration of iron, calcium, magnesium, sodium, chloride and sulfate due to the dissolution of under-saturated minerals. Comparing our results to the findings of previous studies in the region showed a trend of groundwater deterioration, particularly with increased salinity levels. More than 90% of the groundwater samples collected showed low quality for irrigation. Since the groundwaters in Al-Baha are shallow and located in unconfined aquifers, consequently, react promptly to climate change. Thus, it is necessary for decision makers to formulate an action plan to minimize the effect of prospective climate change, and decrease the potential of further deterioration in the quality of groundwater. Therefore, it is imperative to control groundwater discharge, strictly prevent the drilling of new wells, and closely monitor the quality of groundwater in the region. It is also important to implement management practices to minimize agricultural water demand by changing the crop patterns and enhancing water productivity. In particular, policy-makers should enforce regulations and guidelines to control excessive use of fertilizers and pesticides in agricultural production, and implement proper water management to maintain socioeconomic development and sustain groundwater resources in Sarat Al-Baha region.