Next Article in Journal
Editorial—Managing and Planning Water Resources for Irrigation: Smart-Irrigation Systems for Providing Sustainable Agriculture and Maintaining Ecosystem Services
Next Article in Special Issue
Optimizing Laboratory Investigations of Saline Intrusion by Incorporating Machine Learning Techniques
Previous Article in Journal
Estimating Soil Organic Carbon in Agricultural Gypsiferous Soils by Diffuse Reflectance Spectroscopy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Department of Hydrology and Hydraulic Engineering, Vrije Universiteit Brussel, Pleinlaan 2, 1050 Brussels, Belgium
2
Shellharbour City Council, Shellharbour, NSW 2529, Australia
*
Author to whom correspondence should be addressed.
Water 2020, 12(1), 262; https://doi.org/10.3390/w12010262
Submission received: 13 December 2019 / Revised: 7 January 2020 / Accepted: 14 January 2020 / Published: 16 January 2020
(This article belongs to the Special Issue Seawater Intrusion into Coastal Aquifers)

Abstract

:
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, Ca2+, Mg2+, Na+, K+, Cl, SO42–, HCO3 and NO3). 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. Differences 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.

Graphical Abstract

1. 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 sources [4]. The quality of groundwater in the Gaza Strip is influenced by many processes, including return flow from irrigation, wastewater leakage, upconing of underground brine water and sea water intrusion [5,6]. Shomar [4] found that only 10% of the municipal wells meet WHO drinking water standards; in particular the chloride, nitrate and fluoride concentrations are 2–9 times higher than WHO standards. Almasri and Ghabayen [7] reported that a large portion of nitrate concentrations observed in wells in the Gaza Strip did not meet the US Environmental Protection Agency drinking water standard for the years 1990 and 2000–2004. Chemical analysis of municipal wells conducted in 2009 by Aish [8] showed that total dissolve solids, chloride and sodium values exceeded WHO standards. Abbas et al. [9] sampled 58 municipal wells in 2010 and found that no groundwater in the Gaza Strip meets all WHO drinking water standards. Alastal et al. [10] derived groundwater quality maps for the Gaza Strip using six water quality parameters (chloride, nitrate, sulfate, calcium, magnesium, and alkalinity) and found that the fraction for which the groundwater quality is not good increased from 30% to 55% between 2000 and 2010. The Palestinian Water Authority [11] reported that only 12.4% of the wells comply with the WHO standard for nitrate and only 19.3% with the WHO standard for chloride. Abu-alnaeem et al. [12] reported that around 90% of all wells in the Gaza Strip are polluted by nitrate, mainly due to wastewater inputs.
To assess the suitability of groundwater for human consumption, it is essential to determine and evaluate its quality. Researchers have used different methods to express water resources quality. Traditionally, water engineering professionals compare individual chemical parameters with recommended allowable limits. In many regions with scarce water resource, however, the use of water at concentrations slightly above these limits is generally not harmful. Horton [13] proposed a water quality index (WQI) to describe the suitability water for human consumption in a single score that can be ranked into categories using terms such as excellent, good, poor, very poor and unsuitable for use, which are easy to understand for decision makers and consumers. Various methods have been proposed to derive WQI scores [14]. A weighted WQI score is usually used in which ratios of concentrations of water quality parameters and their recommended standard values are weighted and combined in a single number. Recent applications of the weighted WQI approach in groundwater quality studies have been presented [10,15,16,17,18]. All methods used to derive WQI scores are similar, the only difference being the number of parameters (observations) used and their corresponding weights.
The purpose of this study was: (1) to assess the groundwater quality in the Dier al-Balah Governorate of the Gaza Strip using groundwater samples taken from municipal wells in 2009 and 2014, (2) to evaluate the suitability of the groundwater for human consumption using the WQI method, (3) map the spatial distribution of the groundwater quality, and (4) evaluating the evolution of the groundwater quality for the 2009–2014 period.

2. Materials and Methods

2.1. 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 km2 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 km2. 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 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].

2.2. 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 (Ca2+, Mg2+, Na+ and K+) and anions (HCO3, SO42–, Cl and 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:
E =   cations   anions   cations +   anions   × 100 %
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.

2.3. 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, Ca2+, Mg2+, Na+, K+, HCO3, SO42–, Cl and 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 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:
w i = W i / i = 1 n W i
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:
q i = C i / S i
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:
WQI = i = 1 n w i q i
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.

2.4. 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:
WQI ( x ) = m ( x ) + ϵ ( x )
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
γ ( h ) = 1 2 N ( h ) ( i , j ) N ( h ) [ ϵ ( x i ) ϵ ( x j ) ] 2
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:
γ ( h ) = n + ( s n ) [ 1 exp ( 3 h / r ) ]
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
WQI ( x ) = k = 1 p β k q k ( x ) + i = 1 N ( x ) λ i ϵ ( x i )
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 .

3. Results and Discussion

3.1. 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 Table 3 and Table 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 and Table 4 indicate that the ionic dominance pattern is of the order Na+ > Ca2+ > Mg2+ > K+ for the cations and Cl > HCO3 > SO42– > NO3 for the anions, both in 2009 and 2014. 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% Ca2+ + Mg2+, 100% HCO3) and Eastern Mediterranean seawater (78% Na+, 100% Cl + SO42–) [29]. 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% of seawater mixing, while the increase in ion content from 2009 to 2014 is approximately 1% of seawater mixing.

3.2. 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 Table 3 and Table 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.

3.3. 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 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 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.
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.

4. 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–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.

Author Contributions

Conceptualization, M.E.B. and P.K.; methodology, M.E.B. and F.D.S.; software, M.E.B. and P.K.; validation, M.E.B., P.K. and F.D.S.; writing—original draft preparation, M.E.B. and P.K.; writing—review and editing, M.H. and F.D.S.; supervision, M.H. and F.D.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The first author acknowledges the support from Erasmus Mundus and the Vrije Universiteit Brussel for providing financial support for his PhD research. The authors are thankful to the Palestinian Water Authority, the Palestinian Meteorological Office, the Coastal Municipalities Water Utility, and the Ministry of Health for providing data and information.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Foster, S.; Hirata, R.; Gomes, D.; D’Elia, M.; Paris, M. Groundwater Quality Protection: A Guide for Water Utilities, Municipal Authorities, and Environment Agencies; World Bank: Washington, DC, USA, 2002; 114p. [Google Scholar]
  2. Nair, I.S.; Rajaveni, S.P.; Schneider, M.; Elango, L. Geochemical and isotopic signatures for the identification of seawater intrusion in an alluvial aquifer. J. Earth Syst. Sci. 2015, 124, 1281–1291. [Google Scholar] [CrossRef] [Green Version]
  3. Bartram, J.; Ballance, R. Water Quality Monitoring—A Practical Guide to the Design and Implementation of Freshwater Quality Studies and Monitoring Programs; World Health Organization, United Nations Environment Programme, E. & F.N. Spon: London, UK, 1996; 383p. [Google Scholar]
  4. Shomar, B. Groundwater of the Gaza Strip: Is it drinkable? Environ. Geol. 2006, 50, 743–751. [Google Scholar] [CrossRef]
  5. Metcalf and Eddy Inc. Coastal Aquifer Management Program (Camp); Integrated Aquifer Management Plan (Gaza Strip); The United States Agency for International Development and Palestinian Water Authority: Gaza, Palestine, 2000. [Google Scholar]
  6. Ghabayen, S.; McKee, M.; Kemblowski, M. Ionic and isotopic ratios for identification of salinity sources and missing data in the Gaza aquifer. J. Hydrol. 2006, 318, 360–373. [Google Scholar] [CrossRef]
  7. Almasri, M.; Ghabayen, S. Analysis of nitrate contamination of Gaza coastal aquifer, Palestine. J. Hydrol. Eng. 2008, 13, 132. [Google Scholar] [CrossRef]
  8. Aish, A.M. Drinking water quality assessment of the Middle Governorate in the Gaza Strip, Palestine. Water Resour. Ind. 2013, 4, 13–20. [Google Scholar] [CrossRef] [Green Version]
  9. Abbas, M.; Barbieri, M.; Battistel, M.; Brattini, G.; Garone, A.; Parisse, B. Water quality in the Gaza Strip: The present scenario. J. Water Resour. Prot. 2013, 5, 54–63. [Google Scholar] [CrossRef] [Green Version]
  10. Alastal, K.M.; Alagha, J.S.; Abuhabib, A.A.; Ababou, R. Groundwater quality assessment using water quality index (WQI) approach: Gaza Coastal aquifer case study. J. Eng. Res. Tech. 2015, 2, 80–86. [Google Scholar]
  11. PWA (Palestinian Water Authority). Water Resources Status Summary Report/Gaza Strip; Palestinian Water Authority: Ramallah, Palestine, 2016. [Google Scholar]
  12. Abu-alnaeem, M.F.; Yusoff, I.; Ng, T.F.; Alias, Y.; Raksmey, M. Assessment of groundwater salinity and quality in Gaza coastal aquifer, Gaza Strip, Palestine: An integrated statistical, geostatistical and hydrogeochemical approaches study. Sci. Total Environ. 2018, 615, 972–989. [Google Scholar] [CrossRef]
  13. Horton, R.K. An index number system for rating water quality. J. Water Poll. Cont. Fed. 1965, 37, 300–305. [Google Scholar]
  14. Tyagi, S.; Sharma, B.; Singh, P.; Dobhal, R. Water quality assessment in terms of water quality index. J. Am. Water Res. 2013, 1, 34–38. [Google Scholar] [CrossRef]
  15. Sahu, P.; Sikdar, P.K. Hydrochemical framework of the aquifer in and around East Kolkata wetlands, West Bengal, India. Environ. Geol. 2008, 55, 823–835. [Google Scholar] [CrossRef]
  16. Ketata, M.; Gueddari, M.; Bouhlila, R. Use of geographical information system and water quality index to assess groundwater quality in El Khairat deep aquifer (Enfidha, Central East Tunisia). Ara. J. Geos. 2012, 5, 1379–1390. [Google Scholar] [CrossRef]
  17. Rabeiy, R.E. Assessment and modeling of groundwater quality using WQI and GIS in Upper Egypt area. Environ. Sci. Pollut. Res. 2018, 25, 30808. [Google Scholar] [CrossRef]
  18. Kawo, N.S.; Karuppannan, S. Groundwater quality assessment using water quality index and GIS technique in Modjo River Basin, central Ethiopia. J. Afr. Earth Sci. 2018, 147, 300–311. [Google Scholar] [CrossRef]
  19. Ministry of Environmental Affairs (MEnA). Palestinian Environmental Action Plan; Ministry of Environmental Affairs: Gaza, Palestine, 2001; pp. 63–69. [Google Scholar]
  20. El Baba, M.; Kayastha, P.; De Smedt, F. Landfill site selection using multi-criteria evaluation in the GIS interface: A case study from the Gaza Strip, Palestine. Ara. J. Geos. 2015, 8, 7499–7513. [Google Scholar] [CrossRef]
  21. Ubeid, K.F. The nature of the Pleistocene-Holocene palaeosols in the Gaza Strip, Palestine. Geologos 2011, 17, 163–173. [Google Scholar] [CrossRef] [Green Version]
  22. APHA. Standard Methods for the Examination of Water and Wastewater, 21st ed.; American Public Health Association: Washington, DC, USA, 2005; 1368p. [Google Scholar]
  23. World Health Organization. Guidelines for Drinking-Water Quality. Vol. 2, Health Criteria and other Supporting Information, 2nd ed.; World Health Organization: Geneva, Switzerland, 1996; 973p. [Google Scholar]
  24. Domenico, P.A.; Schwartz, F.W. Physical and Chemical Hydrogeology; Wiley: New York, NY, USA, 1990; 824p. [Google Scholar]
  25. Piper, A.M. A graphic procedure in the geochemical interpretation of water-analyses. EoS Trans. Am. Geophys. Union 1944, 25, 914–928. [Google Scholar] [CrossRef]
  26. World Health Organization. Guidelines for Drinking-Water Quality, 4th ed.; World Health Organization: Geneva, Switzerland, 2011; 542p. [Google Scholar]
  27. Pebesma, E.J. Multivariable geostatistics in S: The gstat package. Com. Geos. 2004, 30, 683–691. [Google Scholar] [CrossRef]
  28. Gräler, B.; Pebesma, E.; Heuvelink, G. Spatio-Temporal interpolation using gstat. R J. 2016, 8, 204–218. [Google Scholar] [CrossRef]
  29. Nessim, R.B.; Tadros, H.R.Z.; Taleb, A.E.A.; Moawad, M.N. Chemistry of the Egyptian Mediterranean coastal waters. Egypt J. Aquat. Res. 2015, 41, 1–10. [Google Scholar] [CrossRef] [Green Version]
  30. Qahman, K.; Larabi, A. Evaluation and numerical modeling of seawater intrusion in the Gaza aquifer (Palestine). Hydrogeo. J. 2006, 14, 713–728. [Google Scholar] [CrossRef]
  31. Efron, S.; Fischbach, J.R.; Blum, I.; Karimov, R.I.; Moore, M. The Public Health Impacts of Gaza’s Water Crisis: Analysis and Policy Options; RAND Corporation: Santa Monica, CA, USA, 2018; 87p. [Google Scholar]
  32. World Bank. Securing Water for Development in West Bank and Gaza; International Bank for Reconstruction and Development: Washington, DC, USA, 2018; 30p. [Google Scholar]
Figure 1. Map of the Gaza Strip with indication of the Dier al-Balah Governorate (a); map of the Dier al-Balah Governorate (study area) with the location of the municipal wells used in this study (b).
Figure 1. Map of the Gaza Strip with indication of the Dier al-Balah Governorate (a); map of the Dier al-Balah Governorate (study area) with the location of the municipal wells used in this study (b).
Water 12 00262 g001
Figure 2. Piper plot showing the major ion composition of the groundwater samples and the theoretical mixing line between fresh groundwater and seawater.
Figure 2. Piper plot showing the major ion composition of the groundwater samples and the theoretical mixing line between fresh groundwater and seawater.
Water 12 00262 g002
Figure 3. Plot of the total cation concentration versus the chloride concentration and the theoretical mixing line between fresh groundwater and seawater (solid line) or between groundwater with an initial ion content of 20 meq/L and no chloride and seawater (dotted line, the markings indicate intervals of 1% increase in sea water mixing).
Figure 3. Plot of the total cation concentration versus the chloride concentration and the theoretical mixing line between fresh groundwater and seawater (solid line) or between groundwater with an initial ion content of 20 meq/L and no chloride and seawater (dotted line, the markings indicate intervals of 1% increase in sea water mixing).
Water 12 00262 g003
Figure 4. Observed water quality index plotted against the distance to the coastline and the fitted linear regression lines.
Figure 4. Observed water quality index plotted against the distance to the coastline and the fitted linear regression lines.
Water 12 00262 g004
Figure 5. Empirical variogram of the stochastic residual of the water quality index and fitted exponential variogram model.
Figure 5. Empirical variogram of the stochastic residual of the water quality index and fitted exponential variogram model.
Water 12 00262 g005
Figure 6. Cross-validation results with estimated WQI scores versus observations and the 1:1 identity line with 95% confidence bands (conf.).
Figure 6. Cross-validation results with estimated WQI scores versus observations and the 1:1 identity line with 95% confidence bands (conf.).
Water 12 00262 g006
Figure 7. Maps with the spatial distribution of (a) WQI scores based on data from 2009, (b) WQI scores based on data from 2004, (c) groundwater quality classification for 2009 and (d) groundwater quality classification for 2014.
Figure 7. Maps with the spatial distribution of (a) WQI scores based on data from 2009, (b) WQI scores based on data from 2004, (c) groundwater quality classification for 2009 and (d) groundwater quality classification for 2014.
Water 12 00262 g007
Table 1. List of physico-chemical parameters with corresponding unit, WHO standard ( S i ), WQI weight ( W i ), and normalized weight ( w i ).
Table 1. List of physico-chemical parameters with corresponding unit, WHO standard ( S i ), WQI weight ( W i ), and normalized weight ( w i ).
Chemical ParameterUnitStandardWeightRelative Weight
SiWiwi
pH-6.5–8.530.103
TDSmg/L50050.172
Ca2+mg/L7520.069
Mg2+mg/L5020.069
Na+mg/L20030.103
K+mg/L1210.034
HCO3mg/L12020.069
SO42–mg/L25030.103
Clmg/L25030.103
NO3mg/L4550.172
Total 291.000
Table 2. Classification of the water quality according to the WQI score.
Table 2. Classification of the water quality according to the WQI score.
WQIWater Quality
RangeClass
<0.5Excellent
0.5–1Good
1–2Poor
2–3Very poor
>3Unsuitable
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 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).
WellpHTDSCa2+Mg2+Na+K+HCO3SO42–ClNO3EWQIWater Quality
IDmg/Lmg/Lmg/Lmg/Lmg/Lmg/Lmg/Lmg/Lmg/L%ScoreClass
W17.4032411901056705.727326511261393.763.15Unsuitable
W27.1027371611034807.9574267920200−7.973.18Unsuitable
W37.3029121821325248.1421252987191−0.143.22Unsuitable
W47.60349322612564014.532320412691831.493.54Unsuitable
W57.1031081601416007.52683801140147−0.993.18Unsuitable
W67.602877131867755.52312721028549.342.62Very poor
W77.0027791741286205.12862419731138.642.82Very poor
W87.1025971601396609.13705118211327.042.96Very poor
W97.90195386574804.0233204662384.541.82Poor
W107.802198111715204.2246209759625.112.11Very poor
W118.0071562311301.220160172505.430.92Good
W127.6019321125846011.1323186655992.332.12Very poor
W137.50211480354402.626618671761−6.021.94Poor
W147.30163184383102.83606555242−5.961.57Poor
W157.60200274534004.125422763243−3.531.81Poor
W167.20359816512378017.44275111160190−0.643.77Unsuitable
W177.502100138733606.72481286811022.052.12Very poor
W187.502233121604555.12731566811203.412.27Very poor
W197.502436150973707.7278407698120–4.372.49Very poor
Mean7.452456135875096.93082498231101.242.51
SD0.2870045371634.090126265555.080.75
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).
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).
WellpHTDSCa2+Mg2+Na+K+HCO3SO42–ClNO3EWQIWater Quality
IDmg/Lmg/Lmg/Lmg/Lmg/Lmg/Lmg/Lmg/Lmg/L%ScoreClass
W17.9534092071177606.032831212791442.803.43Unsuitable
W27.7629121761165606.35983311032211−6.913.44Unsuitable
W37.5034931981476606.85783531179236−2.893.91Unsuitable
W47.65351424316876015.239727313892333.734.01Unsuitable
W57.5433951791627407.728141312721712.643.57Unsuitable
W67.713031155988635.62713291259725.442.97Very poor
W77.5329681851376745.731327311241225.803.06Unsuitable
W88.1630031891467205.440953510221693.433.42Unsuitable
W97.77186296595094.2246241708473.431.90Poor
W107.512359132866345.7323277853796.422.45Very poor
W117.66774.264441641.4216711916210.161.04Poor
W127.612170136745036.1382262794118−1.562.45Very poor
W138.20246498535233.730421889383−5.802.35Very poor
W148.10172997583983.7394103623520.011.81Poor
W157.09236691675304.536027777755−1.842.22Very poor
W167.90390616913686020.04565541326204−1.844.11Unsuitable
W178.002954153924617.5293154998126−3.722.77Very poor
W188.102975133836306.13111829771333.062.87Very poor
W197.6826181691065036.8311442872132–2.262.81Very poor
Mean7.8027321511036036.83562959771291.062.87
SD0.3074947391704.2101129293624.510.81

Share and Cite

MDPI and ACS Style

El Baba, M.; Kayastha, P.; Huysmans, M.; De Smedt, F. Evaluation of the Groundwater Quality Using the Water Quality Index and Geostatistical Analysis in the Dier al-Balah Governorate, Gaza Strip, Palestine. Water 2020, 12, 262. https://doi.org/10.3390/w12010262

AMA Style

El Baba M, Kayastha P, Huysmans M, De Smedt F. Evaluation of the Groundwater Quality Using the Water Quality Index and Geostatistical Analysis in the Dier al-Balah Governorate, Gaza Strip, Palestine. Water. 2020; 12(1):262. https://doi.org/10.3390/w12010262

Chicago/Turabian Style

El Baba, Moustafa, Prabin Kayastha, Marijke Huysmans, and Florimond De Smedt. 2020. "Evaluation of the Groundwater Quality Using the Water Quality Index and Geostatistical Analysis in the Dier al-Balah Governorate, Gaza Strip, Palestine" Water 12, no. 1: 262. https://doi.org/10.3390/w12010262

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop