Assessment Framework for Natural Groundwater Contamination in Arid Regions: Development of Indices and Wells Ranking System Using Fuzzy VIKOR Method

: Limited groundwater resources in arid regions have been found polluted for drinking purpose due to the presence of natural minerals and radioactive substances, in the sub-soils, higher than the drinking water quality standards. Municipalities in these regions are spending extensive resources to transport (from well ﬁelds) and treat this raw water to provide safe water to the community. Regular monitoring of various physical, chemical, and radioactive water quality parameters (WQPs) in raw water generates large datasets, which makes it di ﬃ cult to come up with convenient ﬁndings for both the decision-makers and general public. A hierarchical water quality assessment framework develops three sub-indices, an overall water quality index, and a system for ranking of groundwater wells. Fuzzy analytical hierarchy process (fuzzy-AHP) establishes the importance weights of di ﬀ erent WQPs and the sub-indices based on their impacts on human health, treatment processes, distribution system infrastructure, and irrigation applications. Fuzzy VIekriterijumsko KOmpromisno Rangiranje (fuzzy-VIKOR) method aggregates the WQPs’ performance for each well and ranks all the wells in a well ﬁeld based on their overall pollution levels, i.e., remoteness from the applicable standards. For evaluating the pragmatism of the framework, data of 11 WQPs were obtained for 39 wells operating in three di ﬀ erent well ﬁelds located along the boundaries and the central part of Buraydah, Qassim, Saudi Arabia. Spatial water quality maps showing physical, chemical, radioactive, and overall water quality assessment results revealed that the oldest well ﬁeld located in the middle of the city outperforms the other two more recently developed well ﬁelds with lesser anthropogenic activities in their catchments. These ﬁndings testify that the primary source of contamination in deep aquifers is the natural sub-soil condition. The water quality indices will be useful to demonstrate the current situation of groundwater quality in Qassim Region and will facilitate the decision-makers for deﬁning the intended uses of raw water sources (i.e., drinking, unrestricted irrigation, and restricted irrigation) and rehabilitation and renewal planning of the groundwater wells. The framework is applicable in the Kingdom of Saudi Arabia (KSA), Gulf Region, and elsewhere for groundwater quality assessment with desired modiﬁcations.


Introduction
Climate change, agricultural withdrawals, and increasing municipal demands have led to a situation of water scarcity in the Kingdom of Saudi Arabia (KSA) and other Gulf countries.According to the map developed by the World Resource Institute, KSA lies in the list of extremely high water stressed countries based on ratio of withdrawals to the renewable supply of water (i.e., 80%) [1].In addition to quantity, limited groundwater resources in the country have been found polluted (for drinking purposes) due to the presence of high levels of natural compounds in the confined (deep) aquifers [2], while unconfined (shallow) aquifers have been affected by the anthropogenic activates, such as land disposal of wastewater and agricultural runoffs [3].Past studies reported higher concentrations of various water quality parameters (WQPs) than the drinking water quality standards (DWQS), such as total dissolved solids (TDS), iron (Fe), and radionuclides, in different parts of KSA, including Qassim region [1][2][3][4][5][6].
The top priority objective of the Water Directorate in KSA is to provide safe water to the community.Consequently, regular monitoring of raw water quality is mandatory to protect the public health and operate the treatment facilities at their highest efficiency.Raw water transmitted from the well fields (WF) is treated through several treatment processes, where each treatment process is designed to treat a maximum level of a given pollutant.In addition to the health impacts, very high concentrations of different WQPs may affect the efficiency of the treatment processes [7].However, there are numerous circumstances when raw groundwater is being used for agricultural and domestic uses in Qassim Region, such as private wells, fuel stations along highways, parks, and some mosques outside the boundaries of cities.Hence, knowing the overall state of raw water quality is essential for both the public and decision-makers working in concerned organizations.
For large-sized urban water supply systems, the groundwater source is essentially a well field (WF) drawing water through multiple wells.The collected water is then conveyed through large transmission main to the water treatment facility.Deep confined aquifers are generally free from microbial contamination; however, physical, chemical, and radioactive WQPs need to be regularly monitored.As a result, large datasets originate, which are difficult to evaluate and to come up with useful insights.A robust water quality index (WQI) can develop the spatial maps and rank the wells in a well field to facilitate the decision-makers (municipality managers, field engineers) and can inform the general public on the state of natural water resources in their regions.The same indices with appropriate performance levels can simultaneously help in defining the water quality for more than one use, e.g., drinking, unrestricted irrigation (URI), and restricted irrigation (RI).Nevertheless, the indices should be able to incorporate uncertainties due to monitoring inaccuracies, variations in DWQS and health-based guidelines, and subjectivity in expert opinion to linguistically define the resultant water quality.
In an old study conducted (in 1996) on assessment of water quality in Saq aquifer, Sharaf and Hussein [8] found the range of TDS from less than 500 to higher than 3500 mg/L in Qassim region, based on samples of well water.According to their study, the main water type in this region is Na-Ca-Cl-SO 4 .They reported high salt concentration and designated the water quality from "fair" to rarely "poor" in Qassim region in comparison with Tabuk and Hail regions.Soil characteristics and agricultural activities were mentioned as the main sources of pollution in the region.
In some studies, efforts were made to develop water quality indices by aggregating various physico-chemical water quality parameters in KSA.Mohammad et al. [9] developed a WQI, without considering the inherent uncertainties in the data, for 12 groundwater and desalination points in Riyadh, KSA.They developed the WQI, for drinking water quality assessment, based on temperature and dissolved oxygen variations at the sampling points instead of the values of 15 other physico-chemical WQPs measured.Aly et al. [10] and Al-Omran et al. [11] also employed simple weighted sum method (WSM) to aggregate various parameters for developing the WQI to check the suitability of ground water for domestic use in Hafar Albatin and Riyadh, KSA.The above-stated uncertainties were not considered in these studies.Other studies accommodated the uncertainties by using fuzzy logic in the development of groundwater quality index.Although a detailed review cannot be included in this text, a brief summary of some recent and relevant works has been outlined in Table 1.Most of the studies on development of fuzzy-based WQIs were carried out in India and Iran using fuzzy interference (or rule-based) systems.One of the possible reasons could be their primary reliance on groundwater in comparison with developed countries where treated surface water in mostly used for municipal water supplies.For example, surface water accounts for 74% of all water withdrawals in USA [12].
In the development of a water quality index, there are certain challenges which the fuzzy inference system (FIS) approach does not handle.Although FIS is simple to use after the inclusion of fuzzy logic toolbox in MATLAB, it does not precisely compare the monitored concentration of a WQP with its applicable standard value, as one of the membership functions is actually containing the standard somewhere in its range.Consider a situation when two or more samples have much lower concentrations of one or more WQPs than the standards while the rest are higher (with varying degrees) than the standards, here the minimum pollution level (e.g., "very low") will be defined for values equal to or less than DWQS.Now, all the samples with parameters values less than DWQS will achieve the same result, i.e., "very low" pollution.The same is true for other performance levels, e.g., "low", "medium", and "very high".This problem is more highlighted while ranking the alternatives with the same linguistic outcomes, e.g., two, three, or more samples with "medium" water quality.Allocating equal weights to all the parameters is another limitation in FIS; particularly for defining higher pollution levels when the higher concentration of a less important parameter (for instance, chlorides) also results in the same index which is otherwise expected to be due to the violation of more important WQPs, such as Iron.It can also be noticed in the Table 1 that most of the studies evaluated the water quality using a combined WQI; sub-indices can provide more insights on water quality and facilitate in decision-making for well ranking and infrastructure assets management (see results and discussions section for details).
Recently, Haider et al. [1] developed a water quality index using fuzzy analytical hierarchy process (fuzzy-AHP) for weight estimation and fuzzy technique for order of preference by similarity to ideal solution (TOPSIS) for parameter aggregation for the city of Buraydah.The main limitation of their study is that only three WQPs (i.e., radium, iron, and total dissolved solids) were used due to data limitations at the time of study.A total of 24 wells were ranked, based on WQI ranging from "very low" to "very high".Finally, based on the average values, an overall water quality index was developed.However, this WQI does not represent the overall situation of physical (e.g., turbidity and electrical conductivity), chemical (e.g., pH and hardness), and radioactive water quality in the groundwater.In their study, the water quality of a well was compared with the best and the worst water quality in the available data.Whereas the best alternative could itself be not the best performer, similarly the worst might also have few parameters which meet the DWQS.Such limitation can be resolved by considering many WQPs, WQIs, and using an approach in which the wells water quality is compared with the alternative having absolute best and worst water quality.Moreover, they [1] used simple fuzzy-AHP method for weight estimation, which does not account for the uncertainties that lie in the fuzzy range selected by the decision-makers.The α-cut-based fuzzy-AHP accounts for these challenges and gives more precise weights for the criteria.Application of fuzzy-AHP and fuzzy VIekriterijumsko KOmpromisno Rangiranje (fuzzy-VIKOR) methods can resolve the above-stated limitations in the existing WQI with the help of fuzzy triangular functions corresponding to different concentrations of WQPs in each well.The main objective of the proposed research is to develop a framework for assessing the impact of natural contamination of groundwater in arid regions.The framework is applied on the study area of Buraydah City, Qassim, KSA.A methodology based on hierarchical approach is developed to aggregate the performance levels of different WQPs for developing the sub-indices, i.e., physical, chemical, and radioactive WQIs.Finally, sub-indices are aggregated to develop an overall water quality index (WQIo) at the top of the hierarchy.Uncertainties in data limitations, subjectivity in decision-makers opinion, and monitoring errors are addressed by integrating fuzzy logic with the conventional multicriteria analysis.The proposed approach is implemented on 39 wells in three WFs located in the center and periphery of the city's boundaries.The indices will be useful for water directorates in KSA to assess the overall groundwater quality and take future planning decisions on rehabilitation and renewal of aged wells and installation of new wells.Moreover, the WQI will provide convenient and simplified information on the underground water quality to the public in urban areas, agricultural consumers, and other concerned agencies.

Study Area
The northeastern region of KSA primarily relies on the Saq aquifer for its municipal and agricultural water supplies.The groundwater in this region is confined in Paleozoic and Mesozoic sedimentary rocks formations with high concentrations of TDS and implementing membrane treatment to meet DWQS [22].Buraydah, located at 26 • 19 16 N and 43 • 57 32 E, is the capital of Qassim region of KSA (see Figure 1).The community water supply mainly relies on the Saq Aquifer.With the highest population share (42%) of the province, domestic and industrial water demands in the city are exponentially increasing.The oldest well field (i.e., WF-1) constructed 40 years ago, around the old city center, is now located in the center of existing periphery of the city.Around 38% of the wells analyzed in this study are located in well field 1.The WF-2, consisting of 42% of the wells, was developed around 15 years ago, shown at the bottom of Figure 1.Operations of the most recent well field-3 were started almost 10 years ago; in this field, 40% of the wells are just 5 years old.This well field consists of 20% of the wells evaluated in present study.The average depth of the wells in well fields 1 and 2 is 650 m, while the depth of wells in WF-3 varies between 650 m and 760 m.Municipalities perform detailed water quality tests for each well annually, however the parameters which can be checked using meters are performed on a monthly basis, e.g., TDS, electrical conductivity (EC), pH, and Turbidity.The wells are drawing water from a deep confined aquifer, no microbiological contamination was observed in any of the wells, and hence the bacteriological parameters were not included in the analysis.
The raw water transmitted from WF-1 and WF-2 is being treated at a common water treatment facility, while raw water from WF-3 has a new treatment plant.Turbidity, TDS, Iron, individual salts (e.g., Cl, NO 3 , and SO 4 ), and radionuclides are the common naturally occurring pollutants in the groundwater of the study area.The schematic of the two water treatment facilities in the study area (one for WF-1 and WF-2 and one for WF-3), showing different treatment processes, is illustrated in Figure 2 [7].Plain sedimentation occurs at raw water storage.In the old treatment plant, installed for WF-1 and WF-2, some percentage of raw water goes to the conventional process for removal of iron using oxidation followed by sand filtration, while the remaining goes to the more recently introduced iron removal process consisting of ion detention and ultrafiltration.The later system is being used in the new plant for WF-3.High removals of iron, manganese, and turbidity are achieved through sand or ultrafiltration.Reverse osmosis (RO) is the subsequent process for removal of salts.Finally, the RO product is blended with the filtered water to achieve a desired concentration of less than 500 mg/L.Ultrafiltration and RO are also being used for treating naturally occurring radionuclides in groundwater.
Several small and dispersed agricultural areas (farms) can be observed in the close vicinity around the study area boundary in Figure 1.Farmers in these areas are directly pumping raw groundwater for both restricted and unrestricted irrigation.Suitability of raw water quality for this purpose also needs to be assessed.
Water 2020, 12, 423 6 of 25 groundwater for both restricted and unrestricted irrigation.Suitability of raw water quality for this purpose also needs to be assessed.

The Assessment Framework
The methodology developed in present research to develop the ground water quality index is presented in Figure 3.After defining the boundaries of study area, the raw water quality data for each well were obtained from the water directorate of Buraydah city.Locations, flows, depths, and hours of operations for each well were also gathered.Important WQPs were selected based on their concentrations (higher than the applicable DWQS), health significance, unrestricted and restricted irrigation uses, and effects on treatment (i.e., filters and membranes) and distribution (i.e., chlorination practice, and pipe material) systems infrastructure.DWQS in Saudi Arabia and World Health Organization (WHO) health-based guidelines were used to define the minimum desirable limits for the selected WQPs.However, the framework can include additional WQPs if found problematic in future or for other areas to be investigated.groundwater for both restricted and unrestricted irrigation.Suitability of raw water quality for this purpose also needs to be assessed.

The Assessment Framework
The methodology developed in present research to develop the ground water quality index is presented in Figure 3.After defining the boundaries of study area, the raw water quality data for each well were obtained from the water directorate of Buraydah city.Locations, flows, depths, and hours of operations for each well were also gathered.Important WQPs were selected based on their concentrations (higher than the applicable DWQS), health significance, unrestricted and restricted irrigation uses, and effects on treatment (i.e., filters and membranes) and distribution (i.e., chlorination practice, and pipe material) systems infrastructure.DWQS in Saudi Arabia and World Health Organization (WHO) health-based guidelines were used to define the minimum desirable limits for the selected WQPs.However, the framework can include additional WQPs if found problematic in future or for other areas to be investigated.

The Assessment Framework
The methodology developed in present research to develop the ground water quality index is presented in Figure 3.After defining the boundaries of study area, the raw water quality data for each well were obtained from the water directorate of Buraydah city.Locations, flows, depths, and hours of operations for each well were also gathered.Important WQPs were selected based on their concentrations (higher than the applicable DWQS), health significance, unrestricted and restricted irrigation uses, and effects on treatment (i.e., filters and membranes) and distribution (i.e., chlorination practice, and pipe material) systems infrastructure.DWQS in Saudi Arabia and World Health Organization (WHO) health-based guidelines were used to define the minimum desirable limits for the selected WQPs.However, the framework can include additional WQPs if found problematic in future or for other areas to be investigated.
The hierarchical-based framework proposed in the middle of Figure 3 develops different WQIs.The overall WQI at the top of hierarchy will be more desired by the senior management, policymakers, and general public.While the sub-indices will be more useful for the operational manger to rank the Water 2020, 12, 423 7 of 24 existing wells and plan new wells.The indices can also be applied on treated water quality to evaluate the efficiency of treatment plants [7].The framework is robust enough to include additional WQPs, in future, depending on data availability, e.g., As and SiO 2 .It is worth mentioning that both the total coliforms and fecal coliforms were found absent in the groundwater of the study area and not included in present research.

Chemical Parameters
For chemical analysis of groundwater in the study area, the municipality regularly monitors pH, total hardness (TH), calcium (Ca), magnesium (Mg), iron (Fe), fluoride (F), chlorides (Cl), ammonia (NH3), nitrites (NO2), nitrates (NO3), and sulphates (SO4).Amongst these parameters, Ca, Mg, F, NH3, and NO2 were found well below the WHO and SAS in all the wells.As a result, pH, Fe, TH, Cl, NO3, Identification of study area, selection of water quality parameters of concern and their applicable standards followed by collection selected parameters for designated uses, e.g., drinking, unrestricted irrigation, and restricted irrigation.Fuzzy logic is used to deal with uncertainties, including (i) possible errors in measurement of WQPs, (ii) low frequency of water quality monitoring, (iii) missing data, (iv) vagueness in expert opinion in establishing importance weights, and (v) subjective nature of WQI.Development of a WQI is a two-step process.The first step is the estimation of importance weights of the criteria, i.e., WQPs in the proposed research.Analytical hierarchy process (AHP) is a well-known method [3,13] and found appropriate to determine the weights of the parameters arranged in the hierarchical form, as shown in Figure 3.In the second step, the transformed values of different parameters are aggregated.To deal with the above-stated uncertainties, Fuzzy-AHP [23] and Fuzzy-VIKOR [24] methods are used for developing the groundwater quality index for the study area.Finally, geospatial maps are developed to illustrate the state of groundwater quality in the study area.

Physical Parameters
TDS and electrical conductivity (EC) were selected to assess the concentrations of salts, while turbidity was used as a measure of suspended solids in raw well water.Although there are no health-based guidelines established by the World Health Organization (WHO), as per aesthetic guidelines, water having TDS level less than 600 mg/L tastes "Good", water with TDS levels ranged between 600 and 900 mg/L tastes "Fair", and "Poor" describes water with TDS higher than 900 mg/L.Concentrations higher than 1200 mg/L or lower than 100 mg/L make the water either "unacceptable for drinking" or "insipid" [7].As per the drinking water quality standards in KSA, the highest acceptable limit for TDS is 600 mg/L [25].Applicable irrigation water quality standards in KSA recommend TDS levels of 2000 mg/L for unrestricted irrigation and up to 2500 mg/L for restricted irrigation [26].According to the Food and Agriculture Organization (FAO), there is no restriction of use for TDS levels less than 450 mg/L, slight to moderate restriction for 450-2000 mg/L, and severe restriction for levels higher than 2000 mg/L [27].In a review conducted by Hussain et al. [28], TDS levels less than 1000 mg/L were recommended safe, 1000-2000 mg/L as marginal, and the levels higher than 2000 mg/L as hazardous for irrigation use.
EC is associated to the ionic content of water.There are essentially no health-based guidelines for this parameter.Various agencies have recommended maximum limits for EC in drinking water based on their experiences and end use of supplied water.Environmental protection agency of Ireland reported an indicator parametric value of 2500 µS/cm at 20 • C in drinking water, while as per Saudi Arabia Standards (SAS), the acceptable range is 160-1600 µS/cm at 20 • C.
Desired limit of turbidity as per most of the regulations is less than 1 NTU in drinking water; the same has been selected as the target concentration by the treatment facilities operating in Qassim Region.However, maximum guideline value of 5 NTU has also been recommended by WHO and KSA DWQS [25,29].The same value has been recommended for both restricted and unrestricted irrigation in the country [26].All the samples were collected, preserved, and transported following the Standard Methods for Examination of Water and Wastewater.TDS, EC, and pH were monitored using HACH multi-parameter meter.Turbidity was measured with the help of HACH turbidity meter.

Chemical Parameters
For chemical analysis of groundwater in the study area, the municipality regularly monitors pH, total hardness (TH), calcium (Ca), magnesium (Mg), iron (Fe), fluoride (F), chlorides (Cl), ammonia (NH 3 ), nitrites (NO 2 ), nitrates (NO 3 ), and sulphates (SO 4 ).Amongst these parameters, Ca, Mg, F, NH 3 , and NO 2 were found well below the WHO and SAS in all the wells.As a result, pH, Fe, TH, Cl, NO 3 , and SO 4 were selected to assess the chemical water quality of groundwater sources in the study area.Iron, chlorides (Cl), nitrates (NO 3 ), and sulphates (SO 4 ) were tested using HACH UV-VIS Spectrophotometer, Canada.The selection was done based on the monitoring frequency and significance of the parameters for the designated uses.
Although pH has no direct impact on drinking water quality, it is the most important WQP affecting the operations of treatment (e.g., ion detention, coagulation, and reverse osmosis) and distribution systems (e.g., formation of chloramines) [7,29].As per SAS, pH of drinking water should range between 6.5 and 8.5.The recommended standards in Saudi Arabia for irrigation are pH range between 6 and 8.4 [26].
In KSA, the presence of Fe in the groundwater has been frequently reported in past studies [1,30,31].In the most recent edition of WHO drinking water quality guidelines, no health-based limits were suggested for Fe in drinking water systems; nevertheless, accumulation of iron can occur in water supplied with values > 0.05 mg/L, which can lead to the growth of "Crenothrix", i.e., an iron-based bacterium, and >0.3 mg/L can result in staining of fixtures and clothes at the consumer end [32].As per both the SAS and WHO, the target value of Fe in drinking water supplies is 0.3 mg/L [25,29].For irrigation, the standards for iron are 5 mg/L for both restricted and unrestricted irrigation [26].
There are no health-based guidelines established for ammonia; however concentrations up to 1.5 mg/L can lead to odor and 35 mg/L can cause taste problems [29].The concentration of ammonia in Buraydah groundwater was found less than 1 mg/L, hence not included in water quality assessment.Although there are no health-related problems with hard water, due to utilization of more soap, taste complaints, and impacts on water infrastructure, hardness was included in the chemical parameters.The maximum tolerable concentration, as per both the SAS and WHO, is 500 mg/L for potable supplies; however, the values higher than 200 mg/L can lead to scaling problem in treatment units, water mains, and storage tanks.Conversely, very soft water with hardness less than 100 mg/L can be corrosive for water mains [26].
Similarly, the concentrations of natural chloride in deep aquifers higher than 250 mg/L may lead to objectionable taste in drinking water.To avoid soil salinity in KSA, the recommended chlorine concentration for irrigation use is 100 mg/L [26]; however, FAO reported chlorides higher than 350 mg/L (i.e., 10 me/L) should be avoided.Based on the epidemiological studies, NO 3 -N (nitrate as nitrogen) levels higher than 11 mg/L may lead to a risk of methaemoglobinaemia in infants.Both the SAS and WHO recommended 10 mg/L as the highest permissible concentration of NO 3 -N in drinking water; the same standards have been recommend for both the unrestricted and restricted irrigation.To avoid gastrointestinal effects from potable water, water sources with SO 4 higher than 500 mg/L need to be carefully monitored [25].As per both the SAS and WHO, the maximum allowable concentration of SO 4 is 400 mg/L, while the permissible level for unrestricted irrigation is 600 mg/L [26].

Radioactive Parameters
A soluble form of natural Radium (Ra) generally exists in earth metal.Different combinations of three primary isotopes of Ra (i.e., 224 Ra, 226 Ra, and 228 Ra) exist in groundwater [33].Low levels of isotopes of natural radium also exist in the groundwater of some regions of Saudi Arabia [1,[34][35][36].Radium isotopes pose lesser risk to human health in comparison with chemical and biological contaminants [25].The United Nations Scientific Committee on the Effects of Atomic Radiation (UNSCEAR) reported that around 90% of the radon present in water supply entered into human body through inhalation [37].WHO recommended 30 pci/L for combined Radium ( 226+228 Ra), which is significantly higher than the maximum contaminant level (MCL) of 5 pci/L (0.0075 mg/L) and MCL goal (MCLG) of 0 pci/L guided by the United States Environmental Protection Agency (USEPA) [38].The debate on setting one permissible limit or even the limits in a more rational range is still going on.If found higher than 5 mg/L at the source, a complete removal of radium is being achieved at the treatment facility of the study area.
There are no standards available for radionuclides in irrigation water in KSA.Uptake of radionuclides by the plants depends on the transfer factor (ranging from 0.01 to 0.006 for raw vegetables), which is essentially the rate between the concentrations (as Bq kg −1 ) of radionuclides in the plant tissue and the dry soil [39].This means that if raw water contains 100 pci/L, and even if the same concentration in assumed in soil, 1 pci/L will reach to the plant tissue.which is well below the acceptable limits for human health.Although no clear irrigation standards are available for radium, a lower transfer rate can be expected for crops.

Fuzzy Analytic Hierarchy Process
Weights of WQPs and sub-indices were estimated using Fuzzy-AHP which is based on the pairwise comparison using linguistic terms (e.g., 'extremely important' and 'moderate unimportance') to define priorities and/ or posteriority amongst the criteria.In present research, triangular fuzzy numbers (TFN) were used for the linguistic terms for effectively approximating the subjective judgments of six decision-makers from practice and academia.Fuzzy-AHP based on the α-cut approach is used in present research.The approach further addresses the uncertainties in the fuzzy ranges selected by the decision-makers (see Figure 4).
to define priorities and/ or posteriority amongst the criteria.In present research, triangular fuzzy numbers (TFN) were used for the linguistic terms for effectively approximating the subjective judgments of six decision-makers from practice and academia.Fuzzy-AHP based on the -cut approach is used in present research.The approach further addresses the uncertainties in the fuzzy ranges selected by the decision-makers (see Figure 4).
The step-by-step procedure of -cut-based Fuzzy-AHP method is given in the following [40].
Step 1: Develop the pairwise comparison matrix.K number of decision-makers were asked to complete the pairwise comparison matrix using the nine point ranting scale given in Table 2.The fuzzy reciprocal judgment matrix  ̃ was developed for each decision maker as: where i is the criteria and j is the number of the criteria in the matrix, j = 1, 2, …, n.
Complete fuzzy reciprocal matrix  ̃ is defined as: where ̃  represents the relative importance difference between the criteria i and j and is a triangular  Step 2: Perform consistency check. ̃ = [̃  ] is the fuzzy positive reciprocal matrix where ̃  = (  ,   ,   ).In -cut approach the consistency of this matrix is checked for each decision maker with the help of the following equation: where   presents the dimension of the matrix and is the maximum eigenvalue.The step-by-step procedure of -cut-based Fuzzy-AHP method is given in the following [40].
Step 1: Develop the pairwise comparison matrix.K number of decision-makers were asked to complete the pairwise comparison matrix using the nine point ranting scale given in Table 2.The fuzzy reciprocal judgment matrix A k was developed for each decision maker as: where i is the criteria and j is the number of the criteria in the matrix, j = 1, 2, . . ., n.
Complete fuzzy reciprocal matrix R k is defined as: where r ij represents the relative importance difference between the criteria i and j and is a triangular fuzzy numbers (TFN) as Step 2: Perform consistency check.R k = r ij is the fuzzy positive reciprocal matrix where r ij = α ij , β ij , γ ij .In α-cut approach the consistency of this matrix is checked for each decision maker with the help of the following equation: where λ max presents the dimension of the matrix and is the maximum eigenvalue.Subsequently, the consistency ratio (CR) was calculated using Equation ( 4): where RI is the random index and its values were selected from Table 3 depending on the number of WQPs (i.e., n) under each sub-index.Calculated CR values less than '1' are acceptable and shows the consistency in decision-makers' judgment.
Step 3: Estimate the fuzzy weights.
Apply Equation ( 5) on the positive matrix 'k' for each decision maker: can be calculated by setting ∝ = 1 while the lower bound and the upper bound can be estimated by setting ∝ = 0. Now, estimate the criteria weights for all the decision-makers using Equations ( 1) and ( 6): where w i is the criteria weight and the weight vector W = (w i ), i = 1,2, . . .,n.
By applying Equation ( 6) to l, m, and u bounds, corresponding weight vertices were calculated as Calculate the smallest possible constant S k l and the largest possible constant S k lu for minimizing the fuzziness of the weights with the help of Equation (7a,b): Subsequently, the lower and upper bounds of the weight vector were estimated using Equation (8a,b): Finally, develop the fuzzy weigh matrix for each decision-maker as: Step 4: Combine the judgment of all the decision-makers.In this step, the fuzzy weights matrices obtained from Equation ( 9) for each decision maker are integrated with the help of Equation ( 10): where W i denotes the combined fuzzified weight of the criteria i estimated by accumulating the judgments of all the decision-makers K. Instead of finding the crisp weights for each criteria (i.e., sub-indices), the combined fuzzy weights from Equation (10) are used in the aggregation of WQPs' performance for each well with the help of Fuzzy-VIKOR method as described in the following.The VIKOR multicriteria decision-making method, first developed by Serafim Opricovic can deal with conflicting criteria [41].This method ranks the alternatives (groundwater wells) by comparing closeness of each alternative (water quality) to the ideal alternative (desired water quality standards).Extension of the method to Fuzzy VIKOR can help in the situations with uncertain conditions [42].In the following the step-by-step procedure to Fuzzy-VIKOR method is outlined [24]: Step 1: Estimate the fuzzified weights using Fuzzy-AHP method, i.e., W i = (w li , w lm , w lu ), where w li , w lm , and w lu are the lower, medium, and upper limits of the criteria (WQPs) weights represented as TFNs.
Step 2: The positive triangular ideal solution ( f * i ) correspond to a hypothetical well with all the parameters meeting DWQS and the negative triangular ideal solution ( f o i ) correspond to a hypothetical well having all the parameter higher than the highest ranges of WQPs in the study area will be determined as: where I c represents the set of WQPs as the cost criteria.
Based on the discussion made in Section 2.2, f * i and f o i values for all the WQPs are presented in Table 4.The parameters which are less than the DWQS (e.g., nitrates), lowest observed values are considered as the positive ideal solution.
Step 3: Calculate the normalized fuzzy difference ( d ij , j = 1, . . ., J , i = 1, . . ., n) using the following equation: Step 4: Compute fuzzy weighted sum S and fuzzy operator MAX R by the following relationship: Water 2020, 12, 423 13 of 24 Step 5: For ranking of wells, compute the value of Q j using the following equation: where S * is the MIN S j , S or is MAX S r j , R * is the MIN R j , SR or is MAXj R r j , and v is the weight of the strategy. 1As per the WHO drinking water quality guidelines [17] and Saudi Arabian Standards (SAS) [25]; 2 unrestricted irrigation [43,44]; 3 restricted irrigation [43,44].
Step 6: Defuzzified all the fuzzy results to find out crisp S j and Q j values for all the alternatives, i.e., groundwater wells, using the following equation: Step 7: The groundwater pollution index (GWPI) for each well is essentially the crisp Sj scores.Subsequently, estimate the groundwater quality index (GWQI) using the following equation: Step 8: Prioritize all the alternatives (i.e., wells) in descending order starting from the highest crisp Qj value, which correspond to the worst performance in terms of source water quality.The well with highest priority needs to be given more importance in terms of water quality monitoring, rehabilitation or renewal planning.
Step 9: Find out groundwater quality indices by averaging the indices for all wells in each well field as: Note that (GWQI) j is the generic form of the index defined in the proposed methodology.In application, this will replace (WQI)p, WQIc, WQIr, and WQIo.
Step 10: Develop geospatial water quality maps for the study area using the water quality assessment scheme presented in Table 5.The table shows that WQIp (including TDS, Turbidity, and EC) is relatively more important to define the suitability of all the uses.As per the rationale for defining the performance levels, the water quality cannot be considered "High" if the TDS levels are not meting drinking water quality standards.Consequently, filtration followed by reverse osmosis (particle for groundwater sources) becomes inevitable to remove the salts and make the water potable.Another example is raw water with only high content of Fe with TDS levels less than 500 mg/L (i.e., drinking water quality standards), which will result in "high" WQIp but "Medium" or "low" WQIc and hence not suitable for drinking.These results will also come up with "Medium" or "Low" WQIo.However, the treatment system would be required to remove iron only through iron oxidation followed by filtration, but does not need reverse osmosis for TDS removal.
Moreover, significance of WQIc should not be overlooked to compare individual WQPs with the corresponding standards (e.g., nitrates, sulphates, hardness) as described in Table 5. WQIr has a primary significance for drinking water, because the levels in the study area are not too high to effect URI or RI (see discussion on transfer rate for raw vegetables in Section 2.3.2 and more details in [39]).Table 5 also shows that for RI, TDS is the primary parameter of concern.Low (L) 0.0 to <0.63 • For WQIp, it might be a borderline situation for RI application.

•
For WQIc, careful monitoring is required for evaluating individual chemical water quality parameters for RI application.

•
High level of treatment is required to meet drinking water quality standards (DWQS).
Medium (M) 0.6 to <0.96 • Good water quality with most of the parameters meeting standards.But, some parameters are slightly higher than the drinking water quality standards.In general, water quality is suitable for irrigation applications.

•
Careful monitoring is required as the water quality is deteriorating with time, due to aging of wells and lowering of groundwater level, which may affect the effectiveness of existing treatment facilities.

•
High level of treatment is required to meet drinking water quality standards (DWQS).
High (H) ≥0.96 to 1.0 • For WQIp, WQIc, and WQIr, excellent water quality with all the parameters equal or lower than the target values stated in DWQS.

• AND •
For overall water quality index (WQIo) of the well field, all parameters meet DWQS after mixing of drawn water from different wells, i.e., impact of some small increase, of few WQPs, from DWQS in some wells vanishes due to dilution.

Groundwater Quality Monitoring Results
Comprehensive water data for the selected physical, chemical, and radioactive WQPs were obtained for the year 2016 from Water Directorate of the City of Buraydah.For each well field, the calculated average (MEAN), minimum (MIN), maximum (MAX), standard deviation (SD), and coefficient of variation (CV) values for all the WQPs are presented in Table 6.It can be seen in Table 6 that even the minimum TDS levels in WF-1 are higher than drinking water quality standards (600 mg/L) given in Table 4. Around 35% CV shows significant variations among 40 year olds wells in this well field.Average EC values are also very high and close to the upper limit of SAS.Turbidity ranging between 7.2 and 27.2 with a SD of 7.5 is also several times higher than the SAS.These results show that the physical water quality parameters of raw water (WF-1) do not meet SAS and need to be reduced through appropriate treatment (i.e., filtration followed by reverse osmosis) before domestic use.However, pH values remained within the desired range.An average concentration of iron higher than 0.3 mg/L also rationalizes the need of existing treatment facilities in the study area, i.e., iron oxidation and sand filtration or ion detention, coagulation, and ultrafiltration [7].An average value of 339 mg/l shows that the hardness in WF-1 is not problematic, particularly after mixing of water (drawn from different wells) in the transmission main.Further, both the Cl and NO 3 higher than the SAS further justify the need membrane filtration.Although SO 4 was found less than 400 mg/L, a maximum value of 324 justifies the inclusion of this parameter in the chemical WQI.Finally, average radium higher than 30 pci/L shows the presence of natural radioactivity in sub-soil strata, which needs to be treated using ultrafiltration and reverse osmosis [45].
Average TDS levels (1423 mg/L) in WF-2 are higher than those in WF-1, with lesser coefficient of variation of 18.3% in comparison with 35.4% in the case of WF-1.Mean turbidity is also higher than in WF-1.Mean and minimum values for pH lie in the desired range of SAS with an occasional value as high as 10.5.The reason for such high value could be instrumental or sampling error.Such results support the use of fuzzy logic in the aggregation process of different WQPs.Average Fe in WF-2 is also higher than that in both the WF-1 and SAS.A very high CV of 98.5% shows large variations between different wells; these findings show dissimilarities of the iron content in the geological and soil formation of this well field.Average concentrations of TH, Cl, and NO 3 are also higher than both the required DWQS and those in the WF-1; however SO 4 levels ranging between 71 mg/L and 233 mg/L are less than those in WF-1 and meet SAS as well.Radium levels are slightly higher than in WF-1 and thus need to be removed through similar types of treatment.
It can be seen in Table 6 that although the mean concentration of TDS in WF-3 meet the SAS for drinking water, such borderline values with maximum value higher than the standards and 153 mg/L of standard deviation identify the presence of natural salinity in groundwater.Average TDS levels were found to be almost one-third of those in the WF-1 and half of the TDS levels in WF-2.As we know that WF-3 is the newest among all the well fields in the study area, these results are in agreement with the findings of Kent and Landon [46], where they reported a notable increase in TDS levels with time in four groundwater sub basins in San Bernardino County, California.Overall, the turbidity values are higher than those in WF-1 and WF-2.Similar to other well fields, pH in WF-3 also meets SAS.Mean Fe levels are higher than those in WF-1 but lower than those in the WF-2.The levels of TH and Cl are much lower than those in the other two well fields and their entire ranges are below the desired SAS.Although average nitrates are slightly higher (i.e., 10.4 mg/l with maximum value of 15.5 mg/L) than the SAS, the values are less than those in the other two well fields.Similar to TDS, higher concentrations of NO 3 can be expected with the aging of wells in future.Sulphates level are almost consistent with the rest of the study area.Radium levels are higher than in WF-1 but less than in WF-2.

Groundwater Quality Index using Fuzzy-VIKOR
Significant variations in the groundwater samples due to possible changing patterns of natural soil and geology of the study area validate the need of a fuzzy-based index.To estimate the importance weights of water quality parameters, pairwise matrices for physical, chemical, and the overall sub-indices were completed by six decision-makers.As only two parameters (i.e., 226 Ra and 228 Ra) define the WQIr, equal weights were allocated to both of them.Consistency ratios for 18 pairwise matrices were checked, using Equation (4), and found less than 1.Table 7 presents the average estimated weights, from Equation ( 10) with K = 6, using the methodology of fuzzy-AHP described in Section 2.3.1.As the fuzzified weights have to be subsequently used in fuzzy-VIKOR method, crisp weight are not presented in the table.
After the weight estimation process, the aggregation process using fuzzy-VIKOR was also extensively applied for development of water quality indices and ranking of wells.All the detailed results cannot be provided due to space limitations.As an example, the estimation of chemical water pollution index for WF-3 is presented in the following.The performance matrix given in Table 8 presents the values of fuzzified values of six chemical WQPs for the eight wells located in WF-3.The last two columns contain scenarios of two hypothetical wells.The first one is the "best case scenario" with all the parameters equal to or lower than the DWQS, while the second corresponds to the "worst case scenario" with all the parameters having highest possible values held by any other well in the field.Table 9 presents the fuzzy results for S j , R j , and Q j , using the values for positive triangular ideal solution ( f * i ) and the negative triangular ideal solution ( f o i ) given in Table 4 and Equations ( 11) to (13).Defuzzified scores using Equation ( 16) are also presented in Table 9.The defuzzified S j scores are essentially the groundwater pollution index (GWPI) defined in Equation (13).S j scores of 0.0 for "best case scenario" and 1.0 for "worst case scenario" validate the proposed methodology for the estimation of GWPI.Higher values of S j correspond to higher pollution with larger number of water quality parameters violating the DWQS.After calculating the water pollution indices, Equation (17) was used to estimate the water quality indices for each well.Finally, WQIo for each well field was calculated by averaging the WQIo of all the wells using Equation (18).

Discussion
Very interesting results can be noticed in Tables 8 and 9.Only one well in WF-3 (i.e., W3-77) obtained the "High" chemical water quality index (see second last row of Table 9), which corresponds to the case when "all the WQPs meet DWQS", as per the subjective rating described in Table 5.These results can be confirmed from the performance data for these wells in Table 8.All the remaining wells obtained "Medium" WQIc with higher levels, primarily due to higher concentration of NO 3 in addition to the highest concentration of Fe in the entire well field.In general, Fe concentration is less than the irrigation standard of 5 mg/L; it is important to compare NO 3 with the standards of URI and RI for subsequent applications.Check wells obtained "High" and "Low" ratings for the best and worst case scenarios, i.e., WQIc values of 1.0 and 0.0.These results also support the subjective rating defined in Table 5 and the robustness of the proposed framework.It can be seen in Figure 1 that all the three fields are located quite uniformly in and around the study area, i.e., top, middle, and bottom.Therefore, tentative boundaries of the areas corresponding to the three well fields were marked by assuming the midpoints between the two respective well fields.Subsequently, the indices were generated by averaging different water quality indices for all the wells in each well field using Equation (18).The results in terms of spatial groundwater quality variations are presented in Figure 5. Figure 5a-c presents the over physical, chemical, and radioactive water quality indices for each spatial boundary.The central area corresponds to WF-1, the bottom to WF-2, and the area at the top is related to WF-1.Finally, the overall WQI (i.e., WQIo) was calculated by aggregating all the three sub-indices, and the results are presented in Figure 5d.Importance weights of the three indices are given in the last three rows of Table 7.
In Figure 5, Area-2 and Area-3 located on the periphery of the main city's boundary are essentially important for the semi-urban and rural settings in the proximity of the study area.In these areas, most of the population use raw groundwater for general domestic use (excluding drinking) and agriculture.The findings presented in Figure 5b are helpful (if properly disseminated with appropriate details and limitations) for the consumers in these areas.For instance, "Medium" physical and chemical water quality confirms that water is generally suitable for agricultural applications in study area.
Another useful application of the water quality maps presented in Figure 5 is observing the changing in raw water quality.For example, "Medium" WQIp in the case of WF-3 (referring the top area of Figure 5a) is due to lesser TDS in comparison with the other two well fields.Increasing trends in TDS levels with time have been reported in literature [46].It means that the field engineers need to keep an eye on such changes to note when the water quality becomes unsuitable for URI, i.e., changing from "Medium" to "Low".It can be seen in Table 6 that the highest observed value in WF-2 is already higher than 2000 mg/L.It can also be seen that on the basis of WQIo, with highest value of 0.87 of the oldest WF-1 outperforms the other two well fields.
Type and cost of RO membranes vary with the concentration of TDS in raw water, so the municipality can take a more rational decision while selecting the type of membranes for this well field, instead of having the same for the entire study area.In addition, the backwashing period will be longer and the amount of RO reject water might be less and more suitable for reuse after appropriate treatment.The WQIc primarily represents the spatial water quality variations.Once the type of membrane is selected based on TDS levels in raw water, the variations in WQIc have practically no impact on treatment plants' operations.
Finally, all the wells in each well field are ranked based on their overall water quality using Equation (15), i.e., crisp Q j values.The well are ranked in descending order as lower values of Q j correspond to less pollution.The wells with higher levels of pollution are ranked on the top, giving highest priority for their careful monitoring and preference for renewal or replacement activities.The results for all the well fields are listed in Table 10.For instance, in WF-1, W1-13 is ranked on the top because its WQIp is "Medium" and WQIr is "Low".While W1-4 has been ranked at 14 with "Medium" performance level for all the indices.Some interesting wells also exist in the middle ranks, such as W1-5 with "High" WQIc, "Medium" WQIp, but "Low" WQIr.It is worth mentioning that the ranking has been done individually for each well field.If the same municipality is responsible for all the three well fields, managers need to establish the overall ranks of all the wells, which might be different than the ranks listed in Table 10.In the study area, WTP-1 receives raw water from WF-1 and WF-2 and the raw water from WF-3 is transmitted to WTP-2.Therefore, the decision-makers at WTP-1 need to identify the top-ranked wells collectively in both the well fields.
The application of the proposed framework for developing various WQIs and well ranking system on the case of Buraydah city shows its pragmatism to manage many operational difficulties by extracting useful information from the large water quality monitoring datasets.Firstly, the three sub-indices will be useful to comprehend the performance of each well in terms of its physical, chemical, and radioactive water quality, instead of having several parameters.Secondly, the spatial maps developed in Figure 5 are useful for all the field engineers, treatment plant managers, well field

Conclusions and Recommendations
Groundwater in arid environments is diminishing and its water quality has also been deteriorating due to anthropogenic activities and the presence of natural substances in sub-soil strata.Groundwater quality assessment has gained importance in these regions for conservation of this limiting natural resource and planning and management of well field infrastructure.A hierarchical-based water quality assessment framework has been proposed to develop various WQIs and a system for ranking of groundwater wells.Three sub-indices developed in the middle of the hierarchy, including WQIc, WQIr, deduce the physical, chemical, and radioactive water quality from 11 water quality parameters (3 physical, 6 chemical, and 2 radioactive).Subsequently, these sub-indices are aggregated using fuzzy-VIKOR method for assessing the overall water quality (i.e., WQIo) on the top of hierarchy.All the three sub-indices and WQIo elucidate the state of water quality in the well field and nearby areas with the help of spatial maps.
Application of the proposed framework on a case of 39 wells in Buraydah City (Qassim, KSA) found at least one or more WQPs higher than the drinking water quality standards in almost all the wells.A subjective rating scheme ranging from "High" to "Low" is proposed for linguistically defining different types of water quality indices.WQIp and WQIc were found "medium" for all the three well fields, while the WQIr varied from "low" to "medium".However, the overall index (WQIo) was found "medium" for all the well fields.These results show that raw water is not potable without appropriate treatment.These assessment results justify the need of existing treatment processes in the study area.However, raw water can be used for both the unrestricted and restricted irrigation purposes with careful monitoring of individual water quality parameters.WQIo states the overall quality through a single measure for top-level management and public, while the sub-indices would be more useful for technical-level decision-making.
Application of fuzzy logic accommodates various possible uncertainties such as sampling and measurement errors, missing data, and difference of decision-makers opinions in the subjective weight estimation based on pairwise comparison.As small differences in concentrations cannot affect the type or cost of treatment, the results of WQIc cannot affect the existing treatment operations to meet the applicable drinking water quality standards.
Water quality assessment results revealed that the oldest well field located in the middle of the study area outperforms the remaining more recently developed well fields.Relatively new well fields located around the periphery of the city are supposed to be less affected by anthropogenic activities.Such observations affirm that the primary source of contamination in deep aquifers is the natural sub-soil condition.
The ranking of wells prioritizes the rehabilitation and renewal planning needs for underperforming wells in each well field.The spatial groundwater quality maps developed in present study are useful for the top-level management of the municipality, policymakers, operational personnel (water quality monitoring and engineering), as well as general public.The proposed framework provides a more structured approach for groundwater quality assessment in arid environmental regions.Future research may include additional WQPs (such as heavy metals) and data collected for a longer period in the assessment process.Moreover, prioritization of wells is based on their water quality performance, future research can include the cost implication in decision-making processes as well.

Figure 1 .
Figure 1.Study area showing the three well fields constructed in different times.The boundary of the city shown in figure is tentative and was not defined by the municipality of Buraydah.

Figure 2 .
Figure 2. Flow diagram of south water treatment plant in study area (modified from [7]).

Figure 1 .
Figure 1.Study area showing the three well fields constructed in different times.The boundary of the city shown in figure is tentative and was not defined by the municipality of Buraydah.

Figure 1 .
Figure 1.Study area showing the three well fields constructed in different times.The boundary of the city shown in figure is tentative and was not defined by the municipality of Buraydah.

Figure 2 .
Figure 2. Flow diagram of south water treatment plant in study area (modified from [7]).

Figure 2 .
Figure 2. Flow diagram of south water treatment plant in study area (modified from [7]).
Estimation of importance weights of different water quality parameters for each level of hierarchy using Fuzzy-AHP based on -cut approach  Municipality managers  Experts from academia and practice  Fuzzy weights for water quality parameters and water quality indices Development of water quality sub-indices and overall water quality index using Fuzzy-VIKOR method  Municipality managers  Experts from academia and practice  Overall WQI (WQIo)  Physical WQI (WQIp)  Chemical WQI (WQIc)  Radioactive WQI (WQIr) Development of maps showing the spatial water quality variations in the study area and suitability for various designated uses.

Table 1 .
Summary of some relevant past studies on groundwater quality assessment.

Table 2 .
Fuzzy scales and triangular fuzzy numbers (TFN) used for linguistic variables.

Table 3 .
Randomly generated values of consistency index (RI).

Table 4 .
Positive triangular ideal solution ( f i *

Table 5 .
Linguistic scheme for defining the groundwater quality index (GWQI).

Table 6 .
Summary of groundwater quality monitoring results for year 2016 in the study area.

Table 7 .
Estimated weights of WQPs and sub-indices using fuzzy-AHP.

Table 8 .
Performance matrix for eight wells located in Well Field-3.