Groundwater Quality and Suitability Assessment for Irrigation Using Hydrogeochemical Characteristics and Pollution Indices: A Case Study of North Al-Quwayiyah Governorate, Central Saudi Arabia

: The present work focuses on a comprehensive hydrochemical assessment of groundwater within a shallow aquifer located in the central region of Saudi Arabia. This aquifer serves as the principal source of groundwater supply for agricultural irrigation purposes. A total of twelve ground-water samples were systematically collected and subjected to thorough analysis to determine various physicochemical parameters. These parameters encompassed electrical conductivity (EC), pH, temperature, total dissolved solids (TDS), as well as concentrations of major ions. Analysis of major ion analysis was employed to elucidate the chemical attributes of groundwater within the research area. This approach facilitated an understanding of the geochemical evolution and the potential suitability of the groundwater for agriculture irrigation. The recorded range of electrical conductivity (EC) for the groundwater in this region falls between 3283 and 11,000 µ S/cm, with an average value of 6709.4 µ S/cm. The pH levels of the groundwater exhibited a slightly alkaline nature, ranging from 7.8 to 8.6. All sampled wells exhibited brackish water conditions (TDS > 1000 mg/L) based on TDS concentrations. The ﬁndings of this investigation demonstrated that the average ion concentration followed the order: Na + > Ca + > Mg 2+ > K + and Cl − > SO 42 − > HCO 3 − > NO 32 − . Applying Piper’s classiﬁcation system, the groundwater samples were classiﬁed into two predominant categories: Ca-Cl and mix-Ca-Mg-Cl. The predominance of Ca 2+ and Mg 2+ over alkalis indicated mixed ground-water facies inﬂuenced by processes of reverse ion exchange and extensive interactions between water and rock formations. The distinct chemical characteristics of the groundwater were attributed to a combination of factors, including the percolation of irrigation-returned water, chemical transformations within the vadose zone, and extensive interactions with subsurface lithology. Despite major ion concentrations surpassing the permissible limits outlined by the World Health Organization (WHO) for drinking water, the assessment of quality indices was conducted to ascertain the suitability of the groundwater for irrigation purposes. These quality indices encompassed the permeability index (PI), magnesium hazard (MH), Kelly’s ratio (KR), sodium-adsorption ratio (SAR), residual sodium carbonate (RSC), sodium percentage (Na%), and magnesium ratio (MR). A meticulous evaluation of groundwater quality for agricultural utilization indicated that all sampled groundwater sources were deemed suitable for irrigation purposes


Introduction
The quality of water has deteriorated due to the escalating utilization of groundwater reserves globally, leading to an increase in pollutant levels [1].Geochemical processes arising from interactions between aquifers and water play a pivotal role in governing groundwater quality, especially in regions with limited recharge from rainfall, where it serves as a vital resource for domestic and agricultural needs [2][3][4][5].
The Kingdom of Saudi Arabia (KSA) faces a scarcity of renewable water resources due to its arid climate [6].Groundwater is the primary source of agricultural and potable water within the KSA [7,8].Over the past three decades, groundwater extraction in the KSA has surged, reaching a volume of 17 billion m 3 /year [6], meeting 80% of the nation's water demand [9].
Recent years have witnessed a surge in agricultural projects reliant on groundwater resources, contributing to groundwater stress and degradation [6].Coupled with anthropogenic activities and rapidly changing climatic conditions, water shortages are becoming a looming concern [20].Groundwater recharge lags behind pumping rates [21], and excessive exploitation has led to a decline in groundwater levels and compromised quality, with substantial socio-economic consequences [22,23].
The principal source of groundwater in the region is the fractured Precambrian bedrock, overlain by deep aquifers in the sedimentary formations of the Arabian Shield, encompassing a significant sequence of Paleozoic to modern sedimentary deposits [24].Additionally, the shallow aquifers found mainly in valleys play a vital role in the Arabian Shield and coastal zones [15].In the central region of Saudi Arabia, thick sedimentary layers from the Mesozoic and Cenozoic periods serve as an abundant groundwater reservoirs [4].The Arabian Shield, western and southwestern coastal plains, and shallow alluvial aquifers (wadies) constitute significant water sources.High evaporation rates, reduced soil moisture due to the high temperatures, and unpredictable rainfall collectively undermine Saudi Arabia's water accessibility [24].
Assessment of water quality is imperative for its comprehensive management and effective utilization, especially amid growing climate impacts and community health concerns.Groundwater quality hinges on complex interactions involving water-soil-sediment dynamics, the flow system, rock types, and prevailing geochemical conditions encompassing dissolution, redox state, precipitation, leaching, ion exchange, etc. [25].Hydrochemistry in arid aquifers is influenced by several variables such as topography, soil chemistry, watermineral interactions, and subsurface mixing of chemically diverse groundwater [26,27].Primary natural processes impacting groundwater chemistry encompass rock-water interactions [25], evaporation/crystallization, and precipitation-induced dilution.Scarce rainfall and high evaporation rates in arid settings influence groundwater chemistry [28].
Water quality investigations have evolved from traditional approaches such as scatter plots to modern methodologies such as geographic information systems (GIS), fuzzy modeling, and machine learning approaches [29][30][31][32].The classical graphical approach involves ionic ratios and scatter plots to comprehend the geochemical processes [33].Contemporary methods encompass geostatistical simulation, multivariate statistical methods, geographic information systems, analytical hierarchical process (AHP), and algorithmic methods for machine learning (e.g., random forest modeling, artificial neural networks, etc.) enabling effective groundwater quality characterization and visualization.
This study focuses on the North of Al-Quwayiyah region in central Saudi Arabia.Despite limited prior investigations, this study aims to comprehensively assess groundwater chemical composition for potable and irrigation suitability using irrigation quality indicators.The irrigation water quality index identifies suitable groundwater extraction zones for irrigation and potential resource management.Spatial distribution maps further facilitate the assessment of groundwater suitability for varied applications.

The Study Area
The research area is situated within the province of Al-Riyadh in central Saudi Arabia, approximately 150 Km from the country's capital Riyadh.The study area is specifically located in the northern part of Al-Quwayiyah Governate (Figure 1).It has an approximate area of 695 km 2 and the topography exhibits varying altitudes, ranging from 740 m above sea level downstream of the area (the eastern part) to approximately 1270 m at the highest point of the catchment (the western part) (Figure 1).
This study focuses on the North of Al-Quwayiyah region in central Saudi Arabia.Despite limited prior investigations, this study aims to comprehensively assess groundwater chemical composition for potable and irrigation suitability using irrigation quality indicators.The irrigation water quality index identifies suitable groundwater extraction zones for irrigation and potential resource management.Spatial distribution maps further facilitate the assessment of groundwater suitability for varied applications.

The Study Area
The research area is situated within the province of Al-Riyadh in central Saudi Arabia, approximately 150 Km from the country's capital Riyadh.The study area is specifically located in the northern part of Al-Quwayiyah Governate (Figure 1).It has an approximate area of 695 km 2 and the topography exhibits varying altitudes, ranging from 740 m above sea level downstream of the area (the eastern part) to approximately 1270 m at the highest point of the catchment (the western part) (Figure 1).The predominant landscape of Saudi Arabia consists of arid desert conditions characterized by high temperatures.The study area demonstrates an exceedingly hot climate, particularly from May to September of each year when the average monthly temperature might exceed 30 • C (Figure 2a).
KSA is characterized by annual average rainfall ranging from 50 to 100 mm [42].Greater rainfall can be found north of Riyadh and along the Red Sea escarpment, where annual average maximum monthly rainfall rates in January, March, and April (Figure 2b).The average annual evaporation at Riyadh is 2900 mm, with daily peaks that can approach 18.5 mm/day in the summer [42].The predominant landscape of Saudi Arabia consists of arid desert conditions characterized by high temperatures.The study area demonstrates an exceedingly hot climate, particularly from May to September of each year when the average monthly temperature might exceed 30 °C (Figure 2a).
KSA is characterized by annual average rainfall ranging from 50 to 100 mm [42].Greater rainfall can be found north of Riyadh and along the Red Sea escarpment, where annual average maximum monthly rainfall rates in January, March, and April (Figure 2b).The average annual evaporation at Riyadh is 2900 mm, with daily peaks that can approach 18.5 mm/day in the summer [42].

Geological and Hydrogeological Se ing
KSA is tectonically situated on the Arabian Plate.It is made up of volcanic and volcano-sedimentary rocks that have undergone Late Precambrian metamorphism [42].The Arabian Shield and Arabian Platform are the two geological regimes that make up the plate's regional division (Figure 3a).These strata have an eastern dip towards the Arabian Gulf as a consequence of the Arabian Shield basements tilted [42].
The subsurface geology of the research area is composed of the following [42]: 1. Surficial, aeolian, alluvial, and colluvial rocks.Quaternary sediments constitute an unconfined aquifer in the western region and a confined aquifer in the eastern region.

Geological and Hydrogeological Setting
KSA is tectonically situated on the Arabian Plate.It is made up of volcanic and volcano-sedimentary rocks that have undergone Late Precambrian metamorphism [42].The Arabian Shield and Arabian Platform are the two geological regimes that make up the plate's regional division (Figure 3a).These strata have an eastern dip towards the Arabian Gulf as a consequence of the Arabian Shield basements tilted [42].
regional aquifer system (Figure 3b).The Saq Sandstone often dips to the east and is above.4. The Arabian shield complex from the Precambrian basement.
Except in areas where a reversal of the groundwater hydraulic gradient has occurred because of significant water-level declines, natural groundwater flow is towards the northeast within the study area [42].

Groundwater Sampling and Analysis
Twelve wells in North Al-Quwayiyah Governorate (Figure 4) were used to collect groundwater samples for irrigation.Groundwater samples were named starting with GW01 (located east of the study region) and ending with GW12 (located west of the study area).The distribution of these wells is a result of the numerous farming activities near the sample locations.During the collection of groundwater samples in the eastern portion of the research region, the average distance between the wells was found to be approximately 2 km, while the remaining samples had an average distance between the wells of 5-6 km.The subsurface geology of the research area is composed of the following [42]: 1.
Surficial, aeolian, alluvial, and colluvial rocks.Quaternary sediments constitute an unconfined aquifer in the western region and a confined aquifer in the eastern region.

2.
The Permian Khuff Formation's Ash Shiqqah is composed of huge conglomerates and shale.Marl, dolomite, sandstone, and gypsum with fossils are all found in the top section.The Khuff Formation is a discontinuous aquifer of lesser significance.

3.
Saq Sandstone, which dates from the Ordovician to the Cambrian.It represents a regional aquifer system (Figure 3b).The Saq Sandstone often dips to the east and is above.4.
The Arabian shield complex from the Precambrian basement.
Except in areas where a reversal of the groundwater hydraulic gradient has occurred because of significant water-level declines, natural groundwater flow is towards the northeast within the study area [42].

Groundwater Sampling and Analysis
Twelve wells in North Al-Quwayiyah Governorate (Figure 4) were used to collect groundwater samples for irrigation.Groundwater samples were named starting with GW01 (located east of the study region) and ending with GW12 (located west of the study area).The distribution of these wells is a result of the numerous farming activities near the sample locations.During the collection of groundwater samples in the eastern portion of the research region, the average distance between the wells was found to be approximately 2 km, while the remaining samples had an average distance between the wells of 5-6 km.In March 2023, groundwater samples were collected using acid-washed plastic bottles.Each sample, which had a volume of 1000 mL, was taken immediately from the well after it had been pumped.The bo les were acidified, closed, and kept at room temperature for analysis the next day.The chemical analysis was conducted using the recommended standard procedures [46].Based on the APHA standard method, volumetric titration procedures were employed to measure Ca 2+ , Mg 2+ , Cl − , and HCO3 − [46].Ca 2+ and Mg 2+ were quantified using the EDTA titrimetric method, Cl using the argentometric titrimetric method, and HCO3 using the phenolphthalein (P).A flame photometer (AFP 100, Hamburg, Germany) was used to measure K + and Na + , while a spectrophotometer (UV 1600 PC) was used to quantify SO4 2− and NO3 − .These analyses were conducted at the Al-Jazzar Engineering Consulting Company in Riyadh, Saudi Arabia.Groundwater samples were examined in situ using a HANNA portable multi-parameter device to determine pH, electric conductivity (EC), and total dissolved solids (TDS).
The Piper triangular diagram and the Gibbs diagram were used to represent the chemical analysis findings and describe the chemical composition of groundwater in the study area.Identifying the hydrogeochemical facies of groundwater was carried out using Piper's trilinear diagram [47].The two base triangles of the Piper diagram display milliequivalent percentages of cation and anion to express the main ions.The total of the cations and anions is determined to be 100%.The diagonal field, which indicates the whole ionic association, is projected with the suitable positions of the cations and anions for examination [48].The interaction between the chemical composition of water and aquifer geology was investigated using the Gibbs diagram.To represent hydrogeochemical processes in the geochemistry of groundwater, such as water-rock interaction, precipitation, and evaporation, Gibbs (1970) [49] developed two plots.For this reason, the weight proportions of Na + /(Na + + Ca 2+ ) for the cations and Cl − /(Cl − + HCO3 − ) for the anions were used for illustrating the diagrams based on the concentrations of TDS.In March 2023, groundwater samples were collected using acid-washed plastic bottles.Each sample, which had a volume of 1000 mL, was taken immediately from the well after it had been pumped.The bottles were acidified, closed, and kept at room temperature for analysis the next day.The chemical analysis was conducted using the recommended standard procedures [46].Based on the APHA standard method, volumetric titration procedures were employed to measure Ca 2+ , Mg 2+ , Cl − , and HCO 3 − [46].Ca 2+ and Mg 2+ were quantified using the EDTA titrimetric method, Cl using the argentometric titrimetric method, and HCO 3 using the phenolphthalein (P).A flame photometer (AFP 100, Hamburg, Germany) was used to measure K + and Na + , while a spectrophotometer (UV 1600 PC) was used to quantify SO 4 2− and NO 3 − .These analyses were conducted at the Al-Jazzar Engineering Consulting Company in Riyadh, Saudi Arabia.Groundwater samples were examined in situ using a HANNA portable multi-parameter device to determine pH, electric conductivity (EC), and total dissolved solids (TDS).

Groundwater Quality Indices
The Piper triangular diagram and the Gibbs diagram were used to represent the chemical analysis findings and describe the chemical composition of groundwater in the study area.Identifying the hydrogeochemical facies of groundwater was carried out using Piper's trilinear diagram [47].The two base triangles of the Piper diagram display milliequivalent percentages of cation and anion to express the main ions.The total of the cations and anions is determined to be 100%.The diagonal field, which indicates the whole ionic association, is projected with the suitable positions of the cations and anions for examination [48].The interaction between the chemical composition of water and aquifer geology was investigated using the Gibbs diagram.To represent hydrogeochemical processes in the geochemistry of groundwater, such as water-rock interaction, precipitation, and evaporation, Gibbs (1970) [49] developed two plots.For this reason, the weight proportions of Na + /(Na + + Ca 2+ ) for the cations and Cl − /(Cl − + HCO 3 − ) for the anions were used for illustrating the diagrams based on the concentrations of TDS.

Groundwater Quality Indices
The following indices are all necessary in this study to evaluate the quality of groundwater for irrigation [50].These indices include sodium adsorption ratio (SAR), sodium percentage (Na%), residual sodium carbonate (RSC), Kelly's ratio (KR), magnesium hazard (MH), and the permeability index (PI).

Sodium Adsorption Ratio (SAR)
The SAR value is calculated using the concentration of primary alkaline and earth alkaline cations, which is determined in units of meq/L in the water.Increased sodium levels in the soil due to the high levels of salt in irrigation water may change the soil's permeability and cause infiltration problems.It may be possible for the soil to spread out, making it challenging to cultivate.The following equation was used to calculate SAR Keesari et al., 2016 [51]: The estimation of Na% is crucial for managing the supply of water for agricultural applications since high sodium levels in soil and water impede plant growth by decreasing soil permeability.The following equation [52] is used to compute it: All cations and anions are measured in meq/L.The RSC ratio is yet another sign of the water's efficacy for irrigation.When Ca 2+ and Mg 2+ ions dissolve out of groundwater with a higher concentration of HCO 3 − and CO 3 2− , the excess NaHCO 3 and CaCO 3 that results could be harmful to the soil mechanism.The following equation [53] can be used to determine residual sodium carbonate: The Kelly's ratio (KR), which can be determined using the following formula [54], can be used to determine whether water has excessive levels of sodium.
Generally, Ca 2+ and Mg 2+ exist in equilibrium.High Mg 2+ concentrations occasionally affect the equilibrium, and excessive Mg may prevent plants from growing by altering the water to more alkaline.The following formula is used to calculate the MR content [55]: 3.2.6.Permeability Index (PI) In another way, crop yield is impacted when soil permeability is decreased by using mineral-rich water over a longer time.According to the accepted standards, the PI was calculated using the following formula [56]:

Spatial Distribution Maps
Regular monitoring and assessment of a wide range of physical and chemical properties are necessary for groundwater investigations.The collection, preservation, analysis, and visualization of both spatial and non-spatial data can be improved using geographical information systems (GIS).It is frequently used in the evaluation of hazards, spatial mapping, and water quality assessments [57][58][59].In this work, ArcGIS was used to construct spatial distribution maps of the groundwater quality variables using the inverse distance weighted (IDW) approach.

Statistical Approaches
To determine the interrelationships between the hydrochemical factors, two correlation techniques were used: Spearman's rank-order index and Pearson's index, which assess monotonic correlations and measure linear linkages between variables, respectively.In addition, certain ionic ratios between elements illustrate the hydrogeochemical process that may occur.

Major Ion Distribution
The pH ranges from 7.8 to 8.6 and the groundwater is slightly alkaline.The average electrical conductivity (EC) is 6709 µS/cm, with a range of 3283 to 11,000 µS/cm.Each sample exceeds the EC's maximum allowable limit of 1500 µS/cm.The substantial increase in EC values is related to anthropogenic effects and geochemical development of groundwater through rock-water interaction over an extended time.The total dissolved solids (TDS), which have an average value of 4025 mg/L, varied from 1970 to 6600 mg/L.According to the TDS characterization, all of the wells have brackish water (TDS > 1000) [60] (Table 1).There is not much variety in the research area.According to spatial distribution maps, pH, EC, and TDS values all show a similar pattern of trending toward groundwater flow (Figure 5).At samples GW08, GW09, and GW12, relatively high values of EC and TDS were found in the western regions of the basin while low values were found in the center part.The major ions in the groundwater are distributed as follows: Na + > Ca + > Mg 2+ > K + and Cl − > SO4 2− > HCO3 − > NO3 2− .Na + concentrations range widely from 297 to 1478 mg/L, with an average of 804.6 mg/L.Na levels in all samples were higher than the 200 mg/L acceptable limit for drinking water.Ca 2+ levels in the study area range from 320 to 1240 mg/L, with an average of 628.3 mg/L.All samples exceeded the permi ed limit of 200 mg/L according to WHO guidelines.The ion exchange of minerals from these region's rocks is the main cause of Ca 2+ in groundwater.Additionally, this might be brought on by the minerals CaCO3 and CaSO4 found in the soil, as well as gypsum anhydrides, dolomite, etc.The average Mg 2+ concentration was 200 mg/L, but the range was 36 to 300 mg/L.All samples had Mg 2+ values that are higher than the desired 100 mg/L threshold except for two samples.K + concentrations range from 1.8 to 24.6 mg/L, with an average of 12.7 mg/L.In the study location, the spatial distribution pa erns of Ca 2+ , Mg 2+ , Na + , and K + are shown in Figure 6.Except for sodium and magnesium, major cation concentrations all decrease towards the east.The highest sodium concentrations were found in the east (GW01) and west (GW12).Regarding magnesium, the eastern direction had the maximum concentration (GW01).
The content of Cl varies between 665 and 2924 mg/L, with an average of 1579.7 mg/L.Cl concentrations in all samples were higher than the permi ed limit of 250 mg/L.The average SO4 content was 630.83 mg/L, with a range of 430 to 1100 mg/L.SO4 concentration in all samples was higher than the allowed limit of 250 mg/L.The average HCO3 content was 222.5 mg/L, while the range was 101 to 401 mg/L.Figure 6 demonstrates a decrease in element concentrations (HCO3, SO4, and Cl) as groundwater flows eastward.While the western or middle regions of the basin have the highest concentrations.
The average NO3 2− concentration was 21.6 mg/L, with a range of 2.9 to 63 mg/L.NO3 2− is often low in deep and confined aquifers because of the generally decreasing situations The major ions in the groundwater are distributed as follows: Na + > Ca + > Mg 2+ > K + and Cl − > SO 4 2− > HCO 3 − > NO 3 2− .Na + concentrations range widely from 297 to 1478 mg/L, with an average of 804.6 mg/L.Na levels in all samples were higher than the 200 mg/L acceptable limit for drinking water.Ca 2+ levels in the study area range from 320 to 1240 mg/L, with an average of 628.3 mg/L.All samples exceeded the permitted limit of 200 mg/L according to WHO guidelines.The ion exchange of minerals from these region's rocks is the main cause of Ca 2+ in groundwater.Additionally, this might be brought on by the minerals CaCO 3 and CaSO 4 found in the soil, as well as gypsum anhydrides, dolomite, etc.The average Mg 2+ concentration was 200 mg/L, but the range was 36 to 300 mg/L.All samples had Mg 2+ values that are higher than the desired 100 mg/L threshold except for two samples.K + concentrations range from 1.8 to 24.6 mg/L, with an average of 12.7 mg/L.In the study location, the spatial distribution patterns of Ca 2+ , Mg 2+ , Na + , and K + are shown in Figure 6.Except for sodium and magnesium, major cation concentrations all decrease towards the east.The highest sodium concentrations were found in the east (GW01) and west (GW12).Regarding magnesium, the eastern direction had the maximum concentration (GW01).
The content of Cl varies between 665 and 2924 mg/L, with an average of 1579.7 mg/L.Cl concentrations in all samples were higher than the permitted limit of 250 mg/L.The average SO 4 content was 630.83 mg/L, with a range of 430 to 1100 mg/L.SO 4 concentration in all samples was higher than the allowed limit of 250 mg/L.The average HCO 3 content was 222.5 mg/L, while the range was 101 to 401 mg/L.Figure 6 demonstrates a decrease in element concentrations (HCO 3 , SO 4 , and Cl) as groundwater flows eastward.While the western or middle regions of the basin have the highest concentrations.
groundwater ecosystem.This concentration might point to sewage pollution and th age of fertilizers (gypsum, potassium fertilizers, etc.).It can be challenging to disting between municipal and/or agricultural contamination sources [63].In many instance land use pa ern can be used to identify the real source of NO3 pollution.Even if they a common ancestor, there is not always a strong association between them becaus percentage might be diluted by a variety of reactions in the soil horizon and vadose In contrast to the other ions, the NO3 2− has a different geographical distribution tern where high levels were observed toward groundwater flow in the eastern portio the basin (Figure 6) and comparatively low concentrations were found in wells locat the western and central regions.Figure 7 demonstrates the evidence of human act such as urban areas close to certain wells that recorded the highest values.These acti might be responsible for the increases in nitrate levels in the research region.The average NO 3 2− concentration was 21.6 mg/L, with a range of 2.9 to 63 mg/L.NO 3 2− is often low in deep and confined aquifers because of the generally decreasing situations at deeper depths.According to Trojan et al. [61] and White et al. [62], the concentration of NO 3 2− can be utilized as an effective indicator of how land use is interacting with the groundwater ecosystem.This concentration might point to sewage pollution and the usage of fertilizers (gypsum, potassium fertilizers, etc.).It can be challenging to distinguish between municipal and/or agricultural contamination sources [63].In many instances, the land use pattern can be used to identify the real source of NO 3 pollution.Even if they have a common ancestor, there is not always a strong association between them because the percentage might be diluted by a variety of reactions in the soil horizon and vadose zone.
In contrast to the other ions, the NO 3 2− has a different geographical distribution pattern where high levels were observed toward groundwater flow in the eastern portions of the basin (Figure 6) and comparatively low concentrations were found in wells located in the western and central regions.Figure 7 demonstrates the evidence of human activity, such as urban areas close to certain wells that recorded the highest values.These activities might be responsible for the increases in nitrate levels in the research region.Generally, the water quality may be suitable for other purposes, such as irrigation and agriculture, even while the element concentrations are higher than the permi ed level for drinking.However, recent investigations, whether conducted within the Kingdom of Saudi Arabia [64][65][66] or other places such as Egypt [67,68], India [69], Morocco [70], and Tunisia [71] demonstrated agreement with this issue.Therefore, groundwater indices are calculated to assess the quality of groundwater for irrigation.

Sca er Matrix Plot and Correlation Analysis
Using a sca er matrix (Figure 8a), it was possible to examine the relationships and distribution of different hydrochemical parameters in the study area.The sca er matrix revealed that some factors, such as TDS and some other parameters, especially EC, have linear correlations, whereas the remaining elements either show irregular or no relationships.Most of the parameters have irregular distributions, according to the matrix's presented histograms.The relationships between several factors that influence groundwater quality were studied and presented using two distinct correlation types: the Pearson correlation, which relies on normally distributed data for analysis, and the Spearman correlation, which does not.Generally, the water quality may be suitable for other purposes, such as irrigation and agriculture, even while the element concentrations are higher than the permitted level for drinking.However, recent investigations, whether conducted within the Kingdom of Saudi Arabia [64][65][66] or other places such as Egypt [67,68], India [69], Morocco [70], and Tunisia [71] demonstrated agreement with this issue.Therefore, groundwater indices are calculated to assess the quality of groundwater for irrigation.

Scatter Matrix Plot and Correlation Analysis
Using a scatter matrix (Figure 8a), it was possible to examine the relationships and distribution of different hydrochemical parameters in the study area.The scatter matrix revealed that some factors, such as TDS and some other parameters, especially EC, have linear correlations, whereas the remaining elements either show irregular or no relationships.Most of the parameters have irregular distributions, according to the matrix's presented histograms.The relationships between several factors that influence groundwater quality were studied and presented using two distinct correlation types: the Pearson correlation, which relies on normally distributed data for analysis, and the Spearman correlation, which does not.
Figure 8b,c shows the two correlation coefficients (r) that the two types of correlations yielded as results.Significant correlations were between 0.7 and 1 for both kinds of correlations, moderate relationships were between 0.5 and 0.7, weak correlations were from 0.5 to 0.3, while values less than 0.3 represent no correlations.Negative numbers signify negative correlations, whereas positive ones signify positive correlations.An important positive correlation suggested that there was a shared cause, whether it was artificial or natural.Figure 8b,c shows the two correlation coefficients (r) that the two types of correlations yielded as results.Significant correlations were between 0.7 and 1 for both kinds of correlations, moderate relationships were between 0.5 and 0.7, weak correlations were from 0.5 to 0.3, while values less than 0.3 represent no correlations.Negative numbers signify negative correlations, whereas positive ones signify positive correlations.An important positive correlation suggested that there was a shared cause, whether it was artificial or natural.
According to the findings of the two-correlation analysis, TDS exhibits substantial positive correlations with EC, Na, Cl, Ca, Mg, Cl, and SO4, demonstrating that these ions According to the findings of the two-correlation analysis, TDS exhibits substantial positive correlations with EC, Na, Cl, Ca, Mg, Cl, and SO 4 , demonstrating that these ions regulate the aquifer's groundwater salinity and electric conductivity, and that EC and TDS are strongly reliant.Na and Cl, Na and Ca, Na and SO 4 , Ca and SO 4 , and Ca and Cl.All show a highly positive association, with these relationships primarily resulting from natural events such as rock-water interaction and dissolution.The dissolution of rocks including gypsum, halite, and silicates.

Water Type
Using Piper's schematic groundwater facies were categorized.The ternary diagrams (Figure 9) display the enrichment trends of the Cl anion, Na + , Ca 2+ , and mix cation types.The plots of acceptable data on the diamond-shaped field divide the different types of groundwater into two categories: Ca-Mg-Cl and Ca-Cl water types.
ural events such as rock-water interaction and dissolution.The dissolution of rocks including gypsum, halite, and silicates.

Water Type
Using Piper's schematic [47], groundwater facies were categorized.The ternary diagrams (Figure 9) display the enrichment trends of the Cl anion, Na + , Ca 2+ , and mix cation types.The plots of acceptable data on the diamond-shaped field divide the different types of groundwater into two categories: Ca-Mg-Cl and Ca-Cl water types.
The majority of the sample belongs to the Ca-Cl water type, which demonstrates groundwater chemistry is accomplished by reactions in unsaturated zones, rock-water interaction, ion exchange, reverse ion exchange, and anthropogenic impacts.

Hydrochemical Process
To define groundwater properties and understand chemical modification, the findings of groundwater chemistry analyses were examined.To reveal numerous relationships that are useful in pointing to potential sources of changes in groundwater chemistry, it was necessary to convert ions to equivalent weights before plo ing.The main mechanisms governing the chemistry of groundwater were identified using the Gibbs plot.The findings on the Gibbs diagram (Figure 10) reveal that the evaporation process, which is the main cause of increased salinity, results in poor groundwater quality [72,73], and controls groundwater chemistry primarily.This is generally expected because raised salinity results from evaporation considerably increasing the ion concentrations created by chemical weathering of the rock.The majority of the sample belongs to the Ca-Cl water type, which demonstrates groundwater chemistry is accomplished by reactions in unsaturated zones, rock-water interaction, ion exchange, reverse ion exchange, and anthropogenic impacts.

Hydrochemical Process
To define groundwater properties and understand chemical modification, the findings of groundwater chemistry analyses were examined.To reveal numerous relationships that are useful in pointing to potential sources of changes in groundwater chemistry, it was necessary to convert ions to equivalent weights before plotting.The main mechanisms governing the chemistry of groundwater were identified using the Gibbs plot.The findings on the Gibbs diagram (Figure 10) reveal that the evaporation process, which is the main cause of increased salinity, results in poor groundwater quality [72,73], and controls groundwater chemistry primarily.This is generally expected because raised salinity results from evaporation considerably increasing the ion concentrations created by chemical weathering of the rock.The proportions and ionic relations of the main elements found in groundw determined from their ionic concentrations.High correlation (r = 0.82) in the al Ca + Mg ionic abundance plot suggests mix-type pa erns in concentrations.Alk abundant in most of the samples, and Ca + Mg concentrations were higher in th ing samples than they were in the alkalis (Figure 11a).According to Subba Rao ion-exchange processes are shown by an abundance of Na with Ca concentration 11b) and the dominance of Ca and Mg over Na and K.
Alkalis have a strong preference for bonding with Cl and SO4 (Figure 11c); th analysis, which has a positive correlation with Cl + SO4 (r = 0.97) and reveals th all alkalis are consumed when bonding with Cl and SO4, excludes the origin of anthropogenic sources.In addition, a strong association between alkalis with C suggests that rocks rather than anthropogenic activity are the source of this dev [75,76].Additionally, a plot of Na and Cl demonstrates a strong correlation (r = 0 gesting that the dissolution of halite may be the cause of groundwater origin (Fi The proportions and ionic relations of the main elements found in groundwater were determined from their ionic concentrations.High correlation (r = 0.82) in the alkalis with Ca + Mg ionic abundance plot suggests mix-type patterns in concentrations.Alkalis were abundant in most of the samples, and Ca + Mg concentrations were higher in the remaining samples than they were in the alkalis (Figure 11a).According to Subba Rao et al. [74], ionexchange processes are shown by an abundance of Na with Ca concentrations (Figure 11b) and the dominance of Ca and Mg over Na and K.
Alkalis have a strong preference for bonding with Cl and SO 4 (Figure 11c); the current analysis, which has a positive correlation with Cl + SO 4 (r = 0.97) and reveals that nearly all alkalis are consumed when bonding with Cl and SO 4 , excludes the origin of Na + from anthropogenic sources.In addition, a strong association between alkalis with Cl and SO 4 suggests that rocks rather than anthropogenic activity are the source of this development [75,76].Additionally, a plot of Na and Cl a strong correlation (r = 0.98), suggesting that the dissolution of halite may be the cause of groundwater origin (Figure 11f).
Water 2023, 15, x FOR PEER REVIEW 15 of 23 from industry, specialized fertilizers, etc., may result in SO4 enrichment into groundwater systems [77,78].Regardless of lithology or depth, the spatial distribution of NO3 supports the prevalent anthropogenic impact on groundwater.

Compatibility for Irrigation
The impacts of the water's mineral components on both plants and the soil determine whether groundwater is suitable for agriculture [79,80].Estimates of sodium adsorption ratio (SAR), relative quantities of Na + indicated by Na%, and residual sodium carbonate (RSC), Kelly's ratio (KR), magnesium hazard (MH), and the permeability index (PI) are among the general parameters for evaluating the irrigation water quality (Table 2).Figure 11d shows the relative abundance of the anionic facies, where Cl + SO 4 is significantly more prevalent than HCO 3 .First, the main affinities between ionic pairs were evaluated to determine ionic affinity correlations.
The plots of Ca + Mg with Cl + SO 4 have a strong association (r = 0.79), which may indicate that they came from source rocks including dolomites and gypsum (Figure 11e).In addition to geogenic processes, human activity such as sewage contamination, waste from industry, specialized fertilizers, etc., may result in SO 4 enrichment into groundwater systems [77,78].Regardless of lithology or depth, the spatial distribution of NO 3 supports prevalent anthropogenic impact on groundwater.

Compatibility for Irrigation
The impacts of the water's mineral components on both plants and the soil determine whether groundwater is suitable for agriculture [79,80].Estimates of sodium adsorption ratio (SAR), relative quantities of Na + indicated by Na%, and residual sodium carbonate (RSC), Kelly's ratio (KR), magnesium hazard (MH), and the permeability index (PI) are among the general parameters for evaluating the irrigation water quality (Table 2).Groundwater for irrigation is classified depending on its SAR value.SAR values less than 10 (SAR < 10) are classified as "Excellent", while those between 10 and 18 are considered "Good".When the SAR value is between 18 and 26, it is identified "Doubtful" When the SAR value of groundwater exceeds 26, it is called "unsuitable" (SAR > 26) [81,82].
When the SAR value is greater than 10, it is regarded as unsuitable for irrigation purposes according to the World Health Organization (WHO 2011) [83] and the Food and Agriculture Organization (FAO 2017) [84].Higher SAR values increase the danger of sodium salinity, which subsequently in turn impacts crop growth by decreasing the amount of soil water that is accessible and the proportion of the two main minerals, calcium, and magnesium.Our results showed that all groundwater samples, except for GW12 (SAR = 10.04), were appropriate for irrigation (Table 3).The SAR ranged from 3.95 to 10, averaging 6.98 (Table 2).Therefore, 91.6 of groundwater samples in the research area have SAR values in the "Excellent" range, making them acceptable for irrigation.In fact, there are certain crops that have already been cultivated.The farms there primarily produce dates and leaves, with the existence of other fruit trees such as lemon and orange.The spatial distribution map (Figure 12) displays areas that are by high values in the east (GW01) and in the center of the research region (GW09), which recorded values close to 10, as well as GW12 located in the western region.The Na% in the groundwater samples varied from 33.6 to 49.3%, with an average of 41.84% (Table 2).Na% values within ranges 20-40 and 40-60 were recognized as being good and permissible for irrigation [66].Four groundwater samples were found to be good for irrigation, while eight groundwater samples were permissible for irrigation uses (Table 4).Excessive salt concentrations in soil harm aeration, infiltration, and soil composition [85].As seen Figure 12, there is a similarity between Na% and SAR, where the direction of increase is towards the east (GW01), as well as GW09 in the middle of the study region.Agriculture suffers when water exceeds acceptable limits for carbonate and bicarbonate over alkaline earth, especially calcium and magnesium [79,86].Organic matter is essential for agricultural productivity; however, high HCO 3 concentrations cause an increase in pH levels, which dissolve organic matter.To discriminate between the various water classes for irrigation, residual sodium carbonate is used.RSC is considered to be a helpful method for figuring out whether groundwater is acceptable for irrigation through the determination of the carbonate and bicarbonate ratio [87].By precipitating calcium as CO 2 , the extra sodium ion opposes the excesses of calcium and magnesium and produces a negative RSC value in the process.Additionally, the positive value of RSC indicates the increased concentration of calcium and magnesium, which is caused by the interaction with HCO 3 to produce calcium bicarbonate and magnesium bicarbonate [88,89].RSC levels were categorized as excellent for irrigation purposes if they were less than 1.25, dubious if they were between 1.25 and 2.5, and unsuitable if they were more than 2.5 [66].RSC levels in the groundwater samples ranged from −84.21 to −17.85, and an average value of −44.44.The findings revealed that the groundwater samples were appropriate for irrigation (Table 5).In contrast to the spatial distribution maps of the other ratios, the highest readings for this index were found at GW05, GW06, and GW07 in the center of the research region (Figure 12).According to Kelley et al. [90], the values of Kelley's ratio might be used to effectively solve the sodium issue with irrigational water.When Kelley's index is greater than one, there is excessive sodium in the water.As a result, water with Kelley's index below one is appropriate for irrigation, whereas water with a ratio above one is not.The results showed that groundwater samples within the study area varied from 0.5 to 0.97 with an average of 0.72.This means 100% of KR values are less than one (Table 6), indicating that the groundwater was acceptable for irrigation.According to Figure 12, GW01 in the east appears to have recorded the highest value, followed by GW09 in the center of the research region.2, with an average of 34.85.Nine groundwater samples were thought to be below 50%, making them acceptable for irrigation while three groundwater samples were above 50% and therefore unsuitable for irrigation (Table 7).The greatest values were found in GW01 in the east, followed by GW05 and GW04 in the research area's center (Figure 12).The most probable causes of the excessive levels of sodium in the water are the utilization of agrochemicals in agricultural practices and lithological origins resulting in the dissolution of these minerals in the water.Findings varied from 37.17 to 50.97, with 44.37 on average (Table 2).PI can be divided into three categories, according to Doneen (1964) [92]: class I (>75%, appropriate), class II (25-75%, good), and class III (25%, unsuitable).It is advised to use water under classes I and II for irrigation.These results showed that the collected groundwater samples had PI levels between 25 and 75, which be regarded as good for irrigation (Table 8) and are consistent with previous research [69,70].Figure 12 illustrates the spatial distribution, particularly GW07 and GW09 in the middle and GW01 in the east had the greatest values.Although integrated methodologies were used to assess groundwater quality and its suitability for irrigation in the study region, there are still certain limitations that need to be considered in future research.For instance, the distribution of groundwater samples in the research area is not perfect, despite the use of geospatial modeling and other statistical methodologies to solve this issue.Temporal changes will also be considered in future studies, which may give more accurate results.The use of the extra methodologies explains the variances in the study area's hydrogeochemistry.To identify the degree to which the hydrogeochemical processes and the origin of groundwater have affected the research region, additional hydrochemical analyses such as metals and isotopic analysis should be performed.

Conclusions
In arid regions, groundwater stands as one of the most reliable water sources.Saudi Arabia exemplifies such an area, reliant on this enduring and consistent supply, particularly within ongoing agricultural initiatives.This study aimed to ascertain the chemical composition of this groundwater and its viability for irrigation purposes.A total of twelve groundwater samples were procured for subsequent chemical analysis.
The outcomes unveiled discernible ion patterns, with Na + > Ca + > Mg 2+ > K + and Cl − > SO 4 2− > HCO 3 − > NO 3 2− prevailing among the major ions.The groundwater exhibited characteristics of Ca-Cl and Ca-Mg-Cl types, demonstrating an average salinity of 4025 mg/L, potentially attributed to the evaporation process.The question of whether this groundwater is suitable for irrigation necessitates further investigation.To ascertain irrigation feasibility, various contamination indices were employed, including the permeability index (PI), magnesium hazard (MH), Kelly's ratio (KR), sodium adsorption ratio (SAR), residual sodium carbonate (RSC), and sodium percentage (Na%).
Based on SAR analysis, 91.6% of the samples fall within the low-sodium category, while the remaining samples are categorized as medium sodium, indicating excellent and good suitability for irrigation water.Employing the Na% criterion, only 33.3% of the evaluated water samples exhibit a favorable irrigation classification, with 66.7% being deemed permissible.The entirety of the investigated groundwater samples demonstrates suitability for irrigation, as per the residual sodium carbonate (RSC) parameter.In terms of Kelly's ratio (KR), all samples exhibit values below unity, signifying their suitability for irrigation.In consideration of the magnesium hazard (MH), 75% of the analyzed water samples are deemed suitable, whereas 25% are determined unsuitable for irrigation.All assessed samples belong to the favorable permeability index (PI) classes, thus indicating their viability for use in irrigation practices.
Overall, the results of this study indicate that a variety of factors, natural and human activity, affect groundwater quality in the studied area.Additionally, groundwater samples are suitable for irrigation purposes.But to maintain the long-term quantitative and qualitative purity of this resource, it is essential to use groundwater carefully, supported by a thorough, long-term investigation.

Figure 1 .
Figure 1.Location and digital elevation model of the study area.

Figure 1 .
Figure 1.Location and digital elevation model of the study area.

Figure 2 .
Figure 2. (a) Monthly high, average, and low temperatures within the study area [43].(b) Average monthly rainfall amount (mm) in the study area [43].

Figure 2 .
Figure 2. (a) Monthly high, average, and low temperatures within the study area [43].(b) Average monthly rainfall amount (mm) in the study area [43].

Figure 4 .
Figure 4.The Satellite Landsat-9 map shows the locations of twelve groundwater samples within the study area.

Figure 4 .
Figure 4.The Satellite Landsat-9 map shows the locations of twelve groundwater samples within the study area.

Water 2023 , 23 Figure 7 .
Figure 7.The recent map displays the urban region near the wells that recorded the high nitrate concentrations.

Figure 7 .
Figure 7.The recent map displays the urban region near the wells that recorded the high nitrate concentrations.

Figure 9 .
Figure 9. Piper diagram for determining the types of water.

Figure 9 .
Figure 9. Piper diagram for determining the types of water.

Water 2023 ,Figure 10 .
Figure 10.Gibbs diagram illustrates the main affecting the chemistry of grou

Figure 10 .
Figure 10.Gibbs diagram illustrates the main mechanism affecting the chemistry of groundwater.

Table 2 .
Summary of statical indices of groundwater quality.

Table 3 .
Groundwater quality for irrigation according to SAR.

Table 3 .
Groundwater quality for irrigation according to SAR.

Table 4 .
Groundwater quality for irrigation according to Na%.

Table 5 .
Groundwater quality for irrigation according to RSC.

Table 6 .
Groundwater quality for irrigation according to KR. Ca 2+ and Mg 2+ usually keep an equilibrium state; however, in soil systems, Ca 2+ and Mg 2+ react differently, and Mg 2+ deteriorates soil structure, especially in waters that are salt predominant and extremely saline[91].High levels of Mg 2+ are typically caused by exchangeable Na + in irrigation soils.In a state of equilibrium increased Mg 2+ concentrations in water will have a negative impact on soil quality, converting it to alkaline and negatively affecting agricultural yields.Magnesium has the potential to be more harmful to soil than calcium.Magnesium concentrations range from 9.68 to 53.19, as shown in Table

Table 7 .
Groundwater quality for irrigation according to MH.

Table 8 .
Groundwater quality for irrigation according to PI.