Evaluation of the Groundwater Quality Using the Water Quality Index and Geostatistical Analysis in the Dier al-Balah Governorate, Gaza Strip, Palestine

: Groundwater contamination is a major problem in the Gaza Strip. In this study we investigate the groundwater quality in the Dier al-Balah Governorate. Water samples were collected from 19 municipal wells in April 2009 and April 2014 and analyzed for physio-chemical parameters (pH, TDS, Ca 2 + , Mg 2 + , Na + , K + , Cl − , SO 42– , HCO 3 − and NO 3 − ). The aim of the research is to determine the groundwater quality and to produce groundwater quality maps using the water quality index (WQI) method and geostatistical analysis. The results show that all water samples are very saline due to the intrusion of Mediterranean seawater in the coastal aquifer. Di ﬀ erences in chemical composition between 2009 and 2014 indicate that about 1% more seawater was mixed with the groundwater in this period. The majority of the observed chemical parameters of all wells are well above the WHO water quality standards and all WQI values indicate that the water quality is problematic. The spatial variation of the WQI scores is modelled by a deterministic component expressing a linear dependence on the distance to the coastline and a stochastic residual described by an exponential variogram with a practical range of 3000 m. The mapping of the WQI scores and derived water quality classes is achieved through regression-kriging. The results indicate that the groundwater in a large area along the coastline is unsuitable for human consumption and comparison of the maps of 2009 and 2014 shows that this region further expanded by about 700 m inland in a period of 5 years. The results of this study are worrying, but they also contribute to a better understanding of the factors that determine the groundwater quality and can help authorities and stakeholders with sustainable development.


Introduction
Groundwater is an important source of fresh water for human consumption, irrigation and industrial use in many countries of the world. However, residential, industrial, commercial, agricultural and other anthropogenic activities together with natural conditions often lead to a deterioration in groundwater quality [1,2]. That is why an assessment of the quality of groundwater is of great importance for society. Water quality assessment includes an evaluation of the physical, biological and chemical properties of water in relation to the natural quality, intended use and human effects that can influence the health of aquatic systems [3].
More than 90% of the population in the Gaza Strip is dependent on the municipal drinking water network, while the remaining 10% of the population rural areas is dependent on private

Study Area
The Dier al-Balah Governorate is located in the center of the Gaza Strip (Figure 1a). The governorate is about 10 km long and 5.5 km wide, with a total land area of 55.2 km 2 and covers 15% of the Gaza Strip. Located in a coastal zone, transitional between the temperate Mediterranean climate and arid climate of the Negev and Sinai deserts, the Gaza Strip experiences a semi-arid climate with two well-defined seasons: a wet season from October to April and a dry season from May to September. The average daily temperature ranges from 27 • C in summer to 13 • C in winter and the average annual rainfall is about 335 mm/y [19]. In 2017 around 273,200 people lived in the Dier al-Balah Governorate resulting in a population density of 4949 persons per km 2 . In recent decades, the demand for water for domestic purposes has increased noticeably as a result of the high population growth.
the Dier al-Balah Governorate belongs to Pliocene-Pleistocene age and consists of Kurkar Group deposits of calcareous sandstone, silt, clay, unconsolidated sands and conglomerates [21]. The aquifer is the only natural water resource and is used for various purposes, such as drinking water, crop irrigation and industrial uses [10].

Data Collection and Procedures
In April 2009 and April 2014, water samples were collected from 19 municipal wells in the Dier al-Balah Governorate (Figure 1b) by the first author on behalf of the Coastal Municipalities Water Utility, responsible for the production of drinking water. The depth of the wells varies from 60 m to 90 m. Groundwater quality parameters were analyzed using standard procedures and suggested precautions were taken to prevent contamination [22]. During the groundwater sampling, temperature, pH and electrical conductivity (EC) were measured in situ according to WHO guidelines [3,23]. Total dissolved solids (TDS) was derived from EC using a calibrated linear relationship. Concentrations of major cations (Ca 2+ , Mg 2+ , Na + and K + ) and anions (HCO3 -, SO4 2-, Cland NO3 -) were determined in the laboratory of the Coastal Municipalities Water Utility according to standard methods prescribed by the American Public Health Association [22]. For quality assurance and quality control, duplicate tests were performed throughout the laboratory analyses.
An indication of the accuracy of the water analysis can be obtained using the charge-balance equation based on the electroneutrality condition. The deviation from equality is expressed as: = ∑cations ∑anions ∑cations ∑anions 100% (1) where is the charge-balance error, and ∑cations and ∑anions are the sums of cations and anions, respectively, expressed in equivalents per liter. The calculated charge-balance error must be within an acceptable limit of ±10% [24]. The major cation and anion compositions of the groundwater samples are presented in a piper diagram [25], with which hydrogeochemical facies can be categorized and trends in the data can be detected.

Water Quality Index
The derivation of WQI involves weight assignment for each parameter, normalization of the weights, standardization of the parameters, calculation of individual WQI scores and aggregation of the scores. Weights for the parameters (pH, TDS, Ca 2+ , Mg 2+ , Na + , K + , HCO3 -, SO4 2-, Cland NO3 -) given in Table 1 were adapted from previous studies [10,16]. The weights values range from one to five based on their importance for water quality assessment. Total dissolved solids and nitrate The Dier al-Balah Governorate is generally flat with an average topographic elevation of about 45 m above the mean sea level. The soil type is mainly sandy loess or loessal sand in addition to some sandy regosols along the coast [20]. Land use is dominated by agricultural land, urban and built-up areas, natural reserves and recreational areas in the coastal part [20]. The coastal aquifer in the Dier al-Balah Governorate belongs to Pliocene-Pleistocene age and consists of Kurkar Group deposits of calcareous sandstone, silt, clay, unconsolidated sands and conglomerates [21]. The aquifer is the only natural water resource and is used for various purposes, such as drinking water, crop irrigation and industrial uses [10].

Data Collection and Procedures
In April 2009 and April 2014, water samples were collected from 19 municipal wells in the Dier al-Balah Governorate (Figure 1b) by the first author on behalf of the Coastal Municipalities Water Utility, responsible for the production of drinking water. The depth of the wells varies from 60 m to 90 m. Groundwater quality parameters were analyzed using standard procedures and suggested precautions were taken to prevent contamination [22]. During the groundwater sampling, temperature, pH and electrical conductivity (EC) were measured in situ according to WHO guidelines [3,23]. Total dissolved solids (TDS) was derived from EC using a calibrated linear relationship. Concentrations of major cations (Ca 2+ , Mg 2+ , Na + and K + ) and anions (HCO 3 − , SO 4 2-, Cl − and NO 3 − ) were determined in the laboratory of the Coastal Municipalities Water Utility according to standard methods prescribed by the American Public Health Association [22]. For quality assurance and quality control, duplicate tests were performed throughout the laboratory analyses. An indication of the accuracy of the water analysis can be obtained using the charge-balance equation based on the electroneutrality condition. The deviation from equality is expressed as: where E is the charge-balance error, and cations and anions are the sums of cations and anions, respectively, expressed in equivalents per liter. The calculated charge-balance error must be within an acceptable limit of ±10% [24]. The major cation and anion compositions of the groundwater samples are presented in a piper diagram [25], with which hydrogeochemical facies can be categorized and trends in the data can be detected.

Water Quality Index
The derivation of WQI involves weight assignment for each parameter, normalization of the weights, standardization of the parameters, calculation of individual WQI scores and aggregation of the scores. Weights for the parameters (pH, TDS, Ca 2+ , Mg 2+ , Na + , K + , HCO 3 − , SO 4 2-, Cl − and NO 3 − ) given in Table 1 were adapted from previous studies [10,16]. The weights values range from one to five based on their importance for water quality assessment. Total dissolved solids and nitrate concentration are given the maximum weight of five because these are important indicators for assessment of the overall water quality, while potassium is given the minimum weight of one because it poses little risk to consumer health. Other chemical parameters are assigned a weight between one and five depending on their importance in determining the water quality. A factor of three is assigned to pH because of the indication of alkalinity and effect on hardness and to sodium, chloride and sulfate because of their indication of salinity and effect on taste. Calcium, magnesium and bicarbonate are given a weight of two because of their effect on hardness. The normalized weight of each parameter is obtained as: where w i is the normalized weight for parameter i, W i is the assigned weight for parameter i, and n is the total number of parameters. Water quality observations are then standardized by dividing the values by their corresponding water quality standard: where q i is the partial WQI score for parameter i, C i is the observed concentration for parameter i (pH-7 for acidity), and S i is the water quality standard for parameter i. For the water quality standards, we use the WHO 2011 [26] guidelines for drinking-water, as shown in Table 1. The overall WQI score is obtained by summation of the scores of each parameter multiplied by their normalized weight: A WQI value less than one indicates that the water can be used without any precaution, while higher values indicate poorer water quality. WQI scores make it possible to classify the suitability of water for human consumption into categories such as excellent, good, poor, very poor and unsuitable for use as shown in Table 2 [9,14,15]. Our method for deriving WQI scores does not differ from what is in the literature, the only minor difference is that we have not multiplied the values by one hundred as is traditionally done, but that actually serves no purpose.

Geostatistical Data Analysis
Mapping of the WQI estimates is achieved by geostatistical analysis in which an empirical variogram is derived, fitted by a variogram model, cross-validation of the model variogram and regression-kriging using the gstat library in R [27,28]. The spatial variation of WQI is modeled as: where x is the location vector, m(x) is the deterministic component of the water quality index and (x) is the stochastic residual. The empirical variogram of the stochastic residual is estimated based on the observations as where γ is the variogram or semi-variance, h is the lag distance and N(h) is the number of observed pairs within the lag distance interval. The empirical variogram is fitted by an exponential model: where n is the nugget, s is the sill and r is the practical range. The variogram model is verified by cross-validation whereby the data points are removed one by one and predicted by regression-kriging using the remaining data. Finally, maps of WQI are obtained by spatial interpolation using regression-kriging where β k are optimized deterministic model coefficients, q k (x) are external drift components, p is the number of external drift components, λ i are optimized kriging weights and N(x) is the number of data points selected within the neighborhood of x.

Hydro-Chemical Characterization
The results of in the in-situ tests and laboratory analyzes of the water samples from the 19 municipal wells in the Dier al-Balah Governorate in 2009 and 2014 are shown in Table 3; Table 4, respectively. These tables also provide the average and standard deviation for each parameter. The second-to-last column of Tables 3 and 4 shows the charge-balance error obtained with Equation (1). All charge-balance errors are within an acceptable limit of ±10%, except for the sample from well W11 in 2014, which is just above the limit. This small might be due to some ions missing in the analysis, such as for example organic acids that may be present due to contamination by industrial or household waste. Table 3. Physio-chemical analysis of the groundwater samples taken in 2009; also given are the charge-balance error E, the WQI score, the water quality classification, and the mean and standard deviation (SD).  Table 4. Physio-chemical analysis of the groundwater samples taken in 2014; also given are the charge-balance error E, the WQI score, the water quality classification, and the mean and standard deviation (SD).  Figure 2 shows the major ion composition in a piper plot. All water quality samples are clustered in the sodium-chloride category (Na + + K + > 50%, Cl − > 50%), from which we can conclude that the groundwater is very brackish. Also shown in Figure 2 is the theoretical mixing line between fresh groundwater (100% Ca 2+ + Mg 2+ , 100% HCO 3 − ) and Eastern Mediterranean seawater (78% Na + , 100% Cl − + SO 4 2-) [29].
Water 2020, 12, 262 7 of 14 All groundwater samples are clustered around this line, which suggest that the groundwater quality is likely to be influenced by mixing with seawater. There appears to be little difference in the major ion composition between the groundwater samples collected in 2009 and 2014. However, the Piper diagram shows only the relative distribution, while absolute concentration values indicate a clear increase in salinity over time. This is shown in Figure 3 where the total cation concentration is plotted against the chloride concentration, both expressed in milliequivalent per liter. The total cation content is used here because the relative distribution of cations is influenced by adsorption and desorption of cations on soil particles, while chloride is a conservative tracer. In the graph, samples taken from a particular well in 2009 and 2014 are linked to show the chemical evolution over time, which clearly indicates that for all wells both the total cation content and chloride concentration increase considerably. The solid blue line in the graph shows the chemical evolution of mixing fresh water with East-Mediterranean seawater; the ratio between the total cation content and the chloride content for Eastern Mediterranean seawater is 1:0.90 [29]. All the groundwater samples plot above this line because groundwater has a larger cation versus chloride ratio than East-Mediterranean seawater due to the presence of more bicarbonate and sulfate. However, the increase in ion content of the well samples from 2009 to 2014 is generally parallel to the sea water mixing line, which suggests that the groundwater is constantly mixed with seawater. A probable explanation for the chemical evolution of the groundwater is given by the dotted line, which shows the mixing of seawater with groundwater with an initial ion content of 20 meq/L without chloride. This line runs more or less through most of groundwater samples and indicates that the chemical evolution of the groundwater is strongly determined by the intrusion of seawater. The marks on this mixing line indicate increases by 1% sea water during the mixing process. Most groundwater samples plot along a range of 3-6% There appears to be little difference in the major ion composition between the groundwater samples collected in 2009 and 2014. However, the Piper diagram shows only the relative distribution, while absolute concentration values indicate a clear increase in salinity over time. This is shown in Figure 3 where the total cation concentration is plotted against the chloride concentration, both expressed in milliequivalent per liter. The total cation content is used here because the relative distribution of cations is influenced by adsorption and desorption of cations on soil particles, while chloride is a conservative tracer. In the graph, samples taken from a particular well in 2009 and 2014 are linked to show the chemical evolution over time, which clearly indicates that for all wells both the total cation content and chloride concentration increase considerably. The solid blue line in the graph shows the chemical evolution of mixing fresh water with East-Mediterranean seawater; the ratio between the total cation content and the chloride content for Eastern Mediterranean seawater is 1:0.90 [29]. All the groundwater samples plot above this line because groundwater has a larger cation versus chloride ratio than East-Mediterranean seawater due to the presence of more bicarbonate and sulfate. However, the increase in ion content of the well samples from 2009 to 2014 is generally parallel to the sea water mixing line, which suggests that the groundwater is constantly mixed with seawater. A probable explanation for the chemical evolution of the groundwater is given by the dotted line, which shows the mixing of seawater with groundwater with an initial ion content of 20 meq/L without chloride. This line runs more or less through most of groundwater samples and indicates that the chemical evolution of the groundwater is strongly determined by the intrusion of seawater. The marks on this mixing line indicate increases by 1% sea water during the mixing process. Most groundwater samples plot along a range of 3-6% of seawater mixing, while the increase in ion content from 2009 to 2014 is approximately 1% of seawater mixing.

Water Quality
The pH-values show that all samples have alkaline properties as expected for groundwater. There seems to be a slight increase from 2009 to 2014, but all values are well within the WHO standard. Although the pH usually has no direct influence on human health [26], it can have an effect on the hardness of water. There are also no real health guidelines for major ions or TDS in drinking water, but excessive salinity can have an effect on odor and taste, while high levels of calcium and magnesium can lead to scaling. All TDS values are well above the WHO standard, some are 5 to 6 times larger, making the water less suitable for domestic use. TDS increased slightly from 2009 to 2014, except for wells W17 and W18 for which TDS increased by approximately 30%; these wells are close to the coastline. TDS values for groundwater that exceed the WHO standard in Gaza have been reported by other researchers; Aish [8] reported TDS values between 680 mg/l and 3107 mg/l for 9 municipal wells sampled in 2009 and Abbas et al. [9] reported that 88% of the TDS values sampled from 58 municipal wells in 2010 exceeded the WHO standard.
For the major ions, excessive amounts of sodium, chloride and bicarbonate are found, two to four times higher than the WHO standards and for calcium and magnesium one to two times highr than the WHO standards in all wells, except well W11. The levels of potassium and sulfate are moderate and are all below or around the WHO standards. The major ion levels in all wells increased from 2009 to 2014 except potassium. In particular, chloride and sodium levels increased significantly in all wells, some by more than 30%, suggesting salinization due to seawater intrusion in the aquifer. Excessive concentrations of sodium and chloride above WHO standards in the Gaza Strip have also been reported by other researchers; Aish [8] reported sodium levels in the range of 109-680 mg/1 and chloride levels in the range 179-1231 mg/L for nine municipal wells sampled in 2009, Abbas et al. [9] reported that 67% of the 58 wells sampled in 2010 had sodium levels above the WHO standard and 71% above the WHO standard for chloride, Abu-alnaeem et al. [12] found that 75% of the 219 wells sampled in 2013-2014 were above the WHO standard for sodium and 80% above the WHO standard for chloride.
High levels of nitrate in drinking water are undesirable because of a potential reduction to nitrite that can be hazardous to health, especially for infants and pregnant women. Moreover, nitrate in groundwater is mainly due to contamination, for example leaching of fertilizers or manure,

Water Quality
The pH-values show that all samples have alkaline properties as expected for groundwater. There seems to be a slight increase from 2009 to 2014, but all values are well within the WHO standard. Although the pH usually has no direct influence on human health [26], it can have an effect on the hardness of water. There are also no real health guidelines for major ions or TDS in drinking water, but excessive salinity can have an effect on odor and taste, while high levels of calcium and magnesium can lead to scaling. All TDS values are well above the WHO standard, some are 5 to 6 times larger, making the water less suitable for domestic use. TDS increased slightly from 2009 to 2014, except for wells W17 and W18 for which TDS increased by approximately 30%; these wells are close to the coastline. TDS values for groundwater that exceed the WHO standard in Gaza have been reported by other researchers; Aish [8] reported TDS values between 680 mg/L and 3107 mg/L for 9 municipal wells sampled in 2009 and Abbas et al. [9] reported that 88% of the TDS values sampled from 58 municipal wells in 2010 exceeded the WHO standard.
For the major ions, excessive amounts of sodium, chloride and bicarbonate are found, two to four times higher than the WHO standards and for calcium and magnesium one to two times highr than the WHO standards in all wells, except well W11. The levels of potassium and sulfate are moderate and are all below or around the WHO standards. The major ion levels in all wells increased from 2009 to 2014 except potassium. In particular, chloride and sodium levels increased significantly in all wells, some by more than 30%, suggesting salinization due to seawater intrusion in the aquifer. Excessive concentrations of sodium and chloride above WHO standards in the Gaza Strip have also been reported by other researchers; Aish [8] reported sodium levels in the range of 109-680 mg/L and chloride levels in the range 179-1231 mg/L for nine municipal wells sampled in 2009, Abbas et al. [9] reported that 67% of the 58 wells sampled in 2010 had sodium levels above the WHO standard and 71% above the WHO standard for chloride, Abu-alnaeem et al. [12] found that 75% of the 219 wells sampled in 2013-2014 were above the WHO standard for sodium and 80% above the WHO standard for chloride.
High levels of nitrate in drinking water are undesirable because of a potential reduction to nitrite that can be hazardous to health, especially for infants and pregnant women. Moreover, nitrate in groundwater is mainly due to contamination, for example leaching of fertilizers or manure, wastewater disposal, leakage from septic tanks, etc., and as such it indicates the possible presence of other hazardous contaminants, such as bacteria or pesticides. Almost all nitrate levels in the well samples are above the WHO standard, some even three to four times higher; only the samples from wells W9, W14 and W15 in 2009 had a nitrate concentration within the permissible range. Large-scale pollution of groundwater by nitrate has also been reported in other studies; Shomar [4] found that 89% of the 94 wells sampled in 2002-2004 showed nitrate levels above the WHO standard, Almasri and Ghabayen [7] reported that 75% of the 2413 nitrate concentrations observed in 568 wells in 1990 and 2000-2004 exceeded the WHO standard and Abu-alnaeem et al. [12] reported that 90% of 219 wells were nitrate polluted in 2013-2014.
The WQI scores and associated water quality classes are given in the last columns of Tables 3 and 4. All WQI scores are all larger than one, except well W11 in 2009, which indicates that the water quality is problematic. The average WQI score was 2.51 in 2009 and rose to 2.87 in 2014, demonstrating a continuous deterioration in water quality. All water samples collected from wells close to the coast have WQI scores greater than two, indicating very poor water quality and about half of them even have values greater than three, indicating that the groundwater is unsuitable for human consumption. The inland wells generally have WQI scores less than 2.5 and can be classified as having poor or very poor water quality, with the exception of well W11 in 2009, the only groundwater sample with good quality.

Water Quality Mapping
The water quality of the well samples is clearly influenced by the proximity to the coastline, which clearly indicates a strong influence of salinization due to sea water intrusion. Figure 4 shows the observed WQI scores plotted against the distance to the coast. Also shown in the figure are the fitted linear regression lines through the data points, which suggest a significant decrease in WQI score with distance to the sea. The influence of the proximity of the coastline was also noted by other researchers [6,9,12].
Water 2020, 12, x FOR PEER REVIEW 10 of 14 wastewater disposal, leakage from septic tanks, etc., and as such it indicates the possible presence of other hazardous contaminants, such as bacteria or pesticides. Almost all nitrate levels in the well samples are above the WHO standard, some even three to four times higher; only the samples from wells W9, W14 and W15 in 2009 had a nitrate concentration within the permissible range.
Large-scale pollution of groundwater by nitrate has also been reported in other studies; Shomar [4] found that 89% of the 94 wells sampled in 2002-2004 showed nitrate levels above the WHO standard, Almasri and Ghabayen [7] reported that 75% of the 2,413 nitrate concentrations observed in 568 wells in 1990 and 2000-2004 exceeded the WHO standard and Abu-alnaeem et al. [12] reported that 90% of 219 wells were nitrate polluted in 2013-2014. The WQI scores and associated water quality classes are given in the last columns of Tables 3 and 4. All WQI scores are all larger than one, except well W11 in 2009, which indicates that the water quality is problematic. The average WQI score was 2.51 in 2009 and rose to 2.87 in 2014, demonstrating a continuous deterioration in water quality. All water samples collected from wells close to the coast have WQI scores greater than two, indicating very poor water quality and about half of them even have values greater than three, indicating that the groundwater is unsuitable for human consumption. The inland wells generally have WQI scores less than 2.5 and can be classified as having poor or very poor water quality, with the exception of well W11 in 2009, the only groundwater sample with good quality.

Water Quality Mapping
The water quality of the well samples is clearly influenced by the proximity to the coastline, which clearly indicates a strong influence of salinization due to sea water intrusion. Figure 4 shows the observed WQI scores plotted against the distance to the coast. Also shown in the figure are the fitted linear regression lines through the data points, which suggest a significant decrease in WQI score with distance to the sea. The influence of the proximity of the coastline was also noted by other researchers [6,9,12]. The distance to the coast is therefore an external drift component that explains the deterministic part of the spatial model for WQI, expressed by Equation (5). The stochastic residuals of the model can be estimated from Figure 4 as the differences between the WQI values and the corresponding values derived by linear regression. These residuals make it possible to determine the experimental variogram with Equation (6), for which we used a lag of 500 m. The results are shown in Figure 5, which indicates a spatial correlation of the residuals up to a distance of about 3000 m. The empirical variograms of the data for 2009 and 2014 do not differ much, so one suitable variogram model can be  The distance to the coast is therefore an external drift component that explains the deterministic part of the spatial model for WQI, expressed by Equation (5). The stochastic residuals of the model can be estimated from Figure 4 as the differences between the WQI values and the corresponding values derived by linear regression. These residuals make it possible to determine the experimental variogram with Equation (6), for which we used a lag of 500 m. The results are shown in Figure 5, which indicates a spatial correlation of the residuals up to a distance of about 3000 m. The empirical variograms of the data for 2009 and 2014 do not differ much, so one suitable variogram model can be fitted to all the data with the following parameters: exponential type, zero nugget, a sill of 0.23 and a practical range of 3000 m. This variogram model is verified by cross-validation using regression-kriging with the distance to the coastline as a drift component. The results are given in Figure 6, which shows the estimated WQI scores versus the observations. The 1:1 identity line and the 95% confidence bands are also shown. The latter are derived as ±1.96 times the standard deviation of the residuals, which can be estimated as the square root of the sill. The results show that all estimated values are within acceptable range, which proves that the variogram model is acceptable.
Water 2020, 12, x FOR PEER REVIEW 11 of 14 practical range of 3000 m. This variogram model is verified by cross-validation using regression-kriging with the distance to the coastline as a drift component. The results are given in Figure 6, which shows the estimated WQI scores versus the observations. The 1:1 identity line and the 95% confidence bands are also shown. The latter are derived as ±1.96 times the standard deviation of the residuals, which can be estimated as the square root of the sill. The results show that all estimated values are within acceptable range, which proves that the variogram model is acceptable.  Water quality index maps for 2009 and 2014 derived by regression-kriging are shown in Figure  7. The maps clearly indicate a regional trend in water quality as a result of sea water intrusion, in addition to some small-scale variations as a result of local contamination by domestic, agricultural or industrial waste. The regional trend shows a decrease in WQI by about 0.4 per km distance from the sea. Local increases in WQI scores are noted around well W16 and around wells W3 and W4, the first due to high sulfate levels and the later due to high levels of bicarbonate.  practical range of 3000 m. This variogram model is verified by cross-validation using regression-kriging with the distance to the coastline as a drift component. The results are given in Figure 6, which shows the estimated WQI scores versus the observations. The 1:1 identity line and the 95% confidence bands are also shown. The latter are derived as ±1.96 times the standard deviation of the residuals, which can be estimated as the square root of the sill. The results show that all estimated values are within acceptable range, which proves that the variogram model is acceptable.  Water quality index maps for 2009 and 2014 derived by regression-kriging are shown in Figure  7. The maps clearly indicate a regional trend in water quality as a result of sea water intrusion, in addition to some small-scale variations as a result of local contamination by domestic, agricultural or industrial waste. The regional trend shows a decrease in WQI by about 0.4 per km distance from the sea. Local increases in WQI scores are noted around well W16 and around wells W3 and W4, the first due to high sulfate levels and the later due to high levels of bicarbonate. Comparison of the maps of 2009 and 2014 shows a clear increase of the WQI scores over time by about 0.36 on average or about 0.07 per year, especially in the regional trend that expresses the intrusion of salt water. Qahman and  Water quality index maps for 2009 and 2014 derived by regression-kriging are shown in Figure 7. The maps clearly indicate a regional trend in water quality as a result of sea water intrusion, in addition to some small-scale variations as a result of local contamination by domestic, agricultural or industrial waste. The regional trend shows a decrease in WQI by about 0.4 per km distance from the sea. Local increases in WQI scores are noted around well W16 and around wells W3 and W4, the first due to high sulfate levels and the later due to high levels of bicarbonate. Comparison of the maps of 2009 and 2014 shows a clear increase of the WQI scores over time by about 0.36 on average or about 0.07 per year, especially in the regional trend that expresses the intrusion of salt water. Qahman and Larabi [30] also reported a continuous decline in groundwater quality as a result of over-exploitation of groundwater resulting in seawater intrusion.
Water 2020, 12, x FOR PEER REVIEW 12 of 14 Larabi [30] also reported a continuous decline in groundwater quality as a result of over-exploitation of groundwater resulting in seawater intrusion. The groundwater quality classification maps for 2009 and 2014 are also shown in Figure 7. For both years there are only three important water quality classes: poor, very poor and unsuitable for human consumption. In 2009, about 35% of the groundwater in the governorate was unsuitable for human consumption, which increased to 53% in 2014. All unsuitable groundwater is along the coast and the figures show that the boundary between unsuitable and very poor groundwater quality shifted by around 700 m inland in a span of 5 years. If this trend continues, there will soon be no safe groundwater left.

Conclusions
Water samples collected from 19 municipal wells in April 2009 and April 2014 in the Dier al-Balah Governorate of the Gaza Strip indicate a sodium-chloride dominated water quality. Plotted in a piper diagram, the major ion composition appears to be clearly influenced by mixing with Mediterranean seawater. When the chloride concentration is plotted against the total cation concentration, it can also be concluded that the observed groundwater samples are mainly the result of mixing with 3-6% seawater, while in the period 2009-2014 there was an increase by about 1% more seawater mixing with the groundwater.
Almost all observed chemical parameters are above their WHO water quality standard, in particular for total dissolved solids, sodium, chloride and nitrate and for the wells close to the coastline. The differences between 2009 and 2004 show an increase in all major ion concentrations in In 2009, about 35% of the groundwater in the governorate was unsuitable for human consumption, which increased to 53% in 2014. All unsuitable groundwater is along the coast and the figures show that the boundary between unsuitable and very poor groundwater quality shifted by around 700 m inland in a span of 5 years. If this trend continues, there will soon be no safe groundwater left.

Conclusions
Water samples collected from 19 municipal wells in April 2009 and April 2014 in the Dier al-Balah Governorate of the Gaza Strip indicate a sodium-chloride dominated water quality. Plotted in a piper diagram, the major ion composition appears to be clearly influenced by mixing with Mediterranean seawater. When the chloride concentration is plotted against the total cation concentration, it can also be concluded that the observed groundwater samples are mainly the result of mixing with 3-6% seawater, while in the period 2009-2014 there was an increase by about 1% more seawater mixing with the groundwater.
Almost all observed chemical parameters are above their WHO water quality standard, in particular for total dissolved solids, sodium, chloride and nitrate and for the wells close to the coastline. The differences between 2009 and 2004 show an increase in all major ion concentrations in particular for sodium and chloride. The resulting scores for the water quality index indicate that the water quality is very problematic, since almost all samples classify from poor water quality to unsuitable for human consumption.
The observed water quality depends on the distance to the coastline. Thus, the spatial variation of the WQI scores can be described by a model with a deterministic drift depending on the distance to the coastline and a stochastic residual described by an exponential variogram and a practical range of 3000 m. This model is verified by cross-validation and mapping of the groundwater quality is achieved by regression-kriging. The results indicate that the groundwater in a large area along the coastline is not suitable for human consumption and that this zone considerably expanded inland over a period of five years (2009)(2010)(2011)(2012)(2013)(2014). It is therefore appropriate that the authorities concerned in the Dier al-Balah Governorate closely monitor the continued deterioration of groundwater, take the necessary precautions to prevent use of contaminated water for human consumption and develop alternatives or corrective measures to meet the water needs of the region.
Seawater intrusion in the Gaza Strip is a national problem, but mitigation measures are lacking. Possible solutions and mitigation measures as rain water harvesting, seawater desalination and wastewater treatment and reuse are well known and have been recommended in the literature (e.g., [31,32]). It is recommended that seawater intrusion be regularly monitored to improve groundwater management in the Gaza strip.