Evaluation and numerical modeling of seawater intrusion in the Gaza aquifer (Palestine)

A numerical assessment of seawater intrusion in Gaza, Palestine, has been achieved applying a 3-D variable density groundwater flow model. A two-stage finite difference simulation algorithm was used in steady state and transient models. SEAWAT computer code was used for simulating the spatial and temporal evolution of hydraulic heads and solute concentrations of groundwater. A regular finite difference grid with a 400 m² cell in the horizontal plane, in addition to a 12-layer model were chosen. The model has been calibrated under steady state and transient conditions. Simulation results indicate that the proposed schemes successfully simulate the intrusion mechanism. Two pumpage schemes were designed to use the calibrated model for prediction of future changes in water levels and solute concentrations in the groundwater for a planning period of 17 years. The results show that seawater intrusion would worsen in the aquifer if the current rates of groundwater pumpage continue. The alternative, to eliminate pumpage in the intruded area, to moderate pumpage rates from water supply wells far from the seashore and to increase the aquifer replenishment by encouraging the implementation of suitable solutions like artificial recharge, may limit significantly seawater intrusion and reduce the current rate of decline of the water levels.     

Simulation of seawater intrusion into the Khan Yunis area of the Gaza Strip coastal aquifer

 The Gaza Strip coastal aquifer is under severe hydrological stress due to over-exploitation. Excessive pumping during the past decades in the Gaza region has caused a significant lowering of groundwater levels, altering in some regions the normal transport of salts into the sea and reversing the gradient of groundwater flow. The sharp increase in chloride concentrations in groundwater indicates intrusion of seawater and/or brines from the western part of the aquifer near the sea. Simulations of salt-water intrusion were carried out using a two-dimensional density-dependent flow and transport model SUTRA (Voss 1984). This model was applied to the Khan Yunis section of the Gaza Strip aquifer. Simulations were done under an assumption that pumping rates increase according to the rate of population growth, or about 3.8% a year. Model parameters were estimated using available field observations. Numerical simulations show that the rate of seawater intrusion during 1997–2006 is expected to be 20–45 m/yr. The results lead to a better understanding of aquifer salinization due to seawater intrusion and give some estimate of the rate of deterioration of groundwater. Received, September 1997 Revised, January 1998, July 1998 Accepted, August 1998     

Environmental isotope study of seawater intrusion in the coastal aquifer (Syria)

Seawater intrusion into the shallow aquifer in the Syrian coast, north of Latakia (Damsarkho, Ras Ibn Hani) and south of Tartous (Al Hamidieh, Ein Zarka) was studied using hydrochemical and isotopic techniques. The electrical conductivity (EC) distribution map north of Latakia revealed that mixing in this area is the consequence of a frontal intrusion of seawater within the fresh groundwater aquifer which, in turn, results from intensive pumping since the 1960s which has lowered the water table inland below sea level. In Ein Zarka, south of Tartous, in contrast, the EC distribution revealed that seawater intrusion is due to local up-coning as a result of intensive pumping. The deuterium and oxygen-18 relationship is that of a mixing line with a slope of 5.55, indicating an intrusion between freshwater and seawater. In addition, the relationship between oxygen-18 and chloride reveals that the mixing has a dominant role compared to evaporation process. The mixing ratios are estimated to be between 6 and 10% north of Latakia, while they do not exceed 3% south of Tartous. A tritium model was applied to compute the “mean transit time”, which is estimated to be around 10 years, on average, to reach the equilibrium that existed originally between the fresh groundwater and seawater, provided that severe pumping is completely halted and the aquifer is naturally recharged by rainfall and deep percolation of irrigation water, thereby allowing the restoration of the hydraulic gradient. This paper is dedicated to the memory of Dr. Y. Yurtsever.     

Geological structure of the coastal aquifer in the southern part of the Gaza Strip,Palestine

The current study introduces the geological subsurface cross-sections in the southern part of the Gaza Strip to show the structure of the aquifer in the area. The cross-sections give evidence of four subaquifers of the coastal aquifer in the southern part of the Gaza Strip. These cross-sections give the natural reasons for the deterioration of the groundwater in the study area. The results show presence of clay lenses that prevent the replenishment processes of the aquifer of fresh water from the rainfall and returns flow from agricultural activities. Lithological formation was evident as one of the natural causes which accelerate destroying process of the coastal aquifer. The results also show that the structure of the aquifer causes the increase of the groundwater salinity in the Gaza Strip. The cross-sections had shown the shortage of storage capacity of high quantities of fresh water in the coastal aquifer in these areas. The role of lithological formation was evident as one of the natural causes to accelerate the process of destroying the coastal aquifer.     

Experimental and simulation study on the impact of storage and recovery of coastal aquifer to seawater intrusion

Natural Hazards - The horizontal seismic responses in cold regions have not received adequate attentions. In seismic design for frozen ground, computations for unfrozen ground are generally adopted...     

GIS-based GALDIT method for vulnerability assessment to seawater intrusion of the Quaternary coastal Collo aquifer (NE-Algeria)

The overexploitation of groundwater in coastal aquifers is often accompanied by seawater intrusion, intensified by climate change and sea level rise. Heading long-term water quality safety and thus the determination of vulnerable zones to seawater intrusion becomes a significant hydrogeological task for many coastal areas. Due to this background, the present study focussed the established methodology of the GIS-based GALDIT model to assess the aquifer vulnerability to seawater intrusion for the Algerian example of the Quaternary coastal Collo aquifer. According to the result analysis overall, more than half of the total surface of the northern study area can be classified as highly vulnerable. Besides the coastline, the areas nearby the local wadis of Guebli and Cherka occur to be the most vulnerable in the region. In view of further map removal performance as well as single-parameter sensitivity analyses from a coupled perspective respectively the GALDIT parameters, distance from the shore (D) and aquifer hydraulic conductivity (A) have been found to be of key significance regarding the model results (mean effective weightings ~?18–19%). Overall, the study results provide a good approximation basis for future management decisions of the Collo aquifer region, including various perspectives such as identification of suitable settings for prospective groundwater pumping wells.     

A geophysical and geochemical approach for seawater intrusion assessment in the Acquedolci coastal aquifer (Northern Sicily)

Vertical electrical sounding (VES’) surveys and chemical analyses of groundwater have been executed in the coastal plain of Acquedolci (Northern Sicily), with the aim to circumscribe seawater intrusion phenomena. This urbanized area is representative of a more general problem interesting most of Mediterranean littoral areas, where intensive overdraft favors a heavy seawater intrusion through the coastline. Aquifer resistivity seems decisively to be conditioned by the chlorine and magnesium content in the main aquifer of the region. Schlumberger VES’, together with piezometric and chemical–physical information of groundwater, allowed us to perform hydrogeological and geophysical elaborations, to propose the occurrence of a relatively narrow belt marked by fresh–salt water mixing. In the considered plain, pollution risk studies have been already realized by authors with the proposal to identify—by parametric and synthetic zoning of significant hydrogeological elements—the most vulnerable sectors. In detail, an intrinsic vulnerability mapping has been already performed, applying the well-known SINTACS system. This paper intends to give—in this sector of Sicily—an example of integration of different methodologies, showing the role of geophysics to describe the degradation of aquifers on the whole as well as to assess their pollution risk better.     

Integrating an artificial intelligence approach with k-means clustering to model groundwater salinity: the case of Gaza coastal aquifer (Palestine)

Artificial intelligence (AI) techniques have increasingly become efficient alternative modeling tools in the water resources field, particularly when the modeled process is influenced by complex and interrelated variables. In this study, two AI techniques—artificial neural networks (ANNs) and support vector machine (SVM)—were employed to achieve deeper understanding of the salinization process (represented by chloride concentration) in complex coastal aquifers influenced by various salinity sources. Both models were trained using 11 years of groundwater quality data from 22 municipal wells in Khan Younis Governorate, Gaza, Palestine. Both techniques showed satisfactory prediction performance, where the mean absolute percentage error (MAPE) and correlation coefficient (R) for the test data set were, respectively, about 4.5 and 99.8% for the ANNs model, and 4.6 and 99.7% for SVM model. The performances of the developed models were further noticeably improved through preprocessing the wells data set using a k-means clustering method, then conducting AI techniques separately for each cluster. The developed models with clustered data were associated with higher performance, easiness and simplicity. They can be employed as an analytical tool to investigate the influence of input variables on coastal aquifer salinity, which is of great importance for understanding salinization processes, leading to more effective water-resources-related planning and decision making.     

Alternatives to tap water: a case study of the Gaza Strip, Palestine

Gaza Strip is a highly populated, small area in which the groundwater is the main water source. During the last few decades, groundwater quality has deteriorated to a limit that the municipal tap water became brackish and unsuitable for human consumption in most parts of the strip. To overcome this serious situation, several attempts were made to replace the tap water or to improve its quality by water desalination applying the reverse osmosis (RO) technology, bottled water, importing water, and storm water harvesting. Water desalination, which is widely applied in the strip, has environmental and health risks. Brine water resulting from the RO systems has adverse environmental effects, whereas the produced fresh water with very low chemical concentrations may have harmful effects on human health. Therefore, introducing international water quality standards to delineate the harmless minimum limits is necessary. Storm water harvesting can provide a partial solution for the water problems in the strip.     

Stochastic forecasts of seawater intrusion towards sustainable groundwater management: application to the Korba aquifer (Tunisia)

A stochastic study of long-term forecasts of seawater intrusion with an application to the Korba aquifer (Tunisia) is presented. Firstly, a geostatistical model of the exploitation rates was constructed, based on a multi-linear regression model combining incomplete direct data and exhaustive secondary information. Then, a new method was designed and used to construct a geostatistical model of the hydraulic conductivity field by combining lithological information and data from hydraulic tests. Secondly, the effects of the uncertainties associated with the pumping rates and the hydraulic conductivity field on the 3D density-dependent transient model were analysed separately and then jointly. The forecasts of the impacts of two different management scenarios on seawater intrusion in the year 2048 were performed by means of Monte Carlo simulations, accounting for uncertainties in the input parameters as well as possible changes of the boundary conditions. Combining primary and secondary data allowed maps of pumping rates and the hydraulic conductivity field to be constructed, despite a lack of direct data. The results of the stochastic long-term forecasts showed that, most probably, the Korba aquifer will be subject to important losses in terms of regional groundwater resources.     

A simulation‐optimization model to control seawater intrusion in coastal aquifers using abstraction/recharge wells

Seawater intrusion is one of the most serious environmental problems in many coastal regions all over the world. Mixing a small quantity of seawater with groundwater makes it unsuitable for use and can result in abandonment of aquifers. Therefore, seawater intrusion should be prevented or at least controlled to protect groundwater resources. This paper presents development and application of a simulation‐optimization model to control seawater intrusion in coastal aquifers using different management scenarios; abstraction of brackish water, recharge of freshwater, and combination of abstraction and recharge. The model is based on the integration of a genetic algorithm optimisation technique and a coupled transient density‐dependent finite element model. The objectives of the management scenarios include determination of the optimal depth, location and abstraction/recharge rates for the wells to minimize the total costs for construction and operation as well as salt concentrations in the aquifer. The developed model is applied to analyze the control of seawater intrusion in a hypothetical confined coastal aquifer. The efficiencies of the three management scenarios are examined and compared. The results show that combination of abstraction and recharge wells is significantly better than using abstraction wells or recharge wells alone as it gives the least cost and least salt concentration in the aquifer. The results from this study would be useful in designing the system of abstraction/recharge wells to control seawater intrusion in coastal aquifers and can be applied in areas where there is a risk of seawater intrusion. Copyright © 2011 John Wiley & Sons, Ltd.     

Use of time-domain electromagnetic (TDEM) method to investigate seawater intrusion in the Lebna coastal aquifer of eastern Cap Bon,Tunisia

A multi-method geophysical survey, using time-domain electromagnetic (TDEM) fields, is conducted on the north-eastern coastal aquifer of Cap Bon, Tunisia. The aim of the study is to investigate the inland extent saltwater intrusion, to estimate the fresh-saltwater transition zone, and to characterize the lithologic spatial variability of the underground layers. Moreover, the TDEM is used to determine both fluid resistivity and bulk resistivity which are required to calculate the porosity by applying the Archie’s equation. The survey of TDEM profiles is carried out along the left and right overbanks of the Lebna river from the coast until the Lebna dam. The interpretation and the analysis of these profiles allow us to better identify the flux exchanges between the river and the underlying aquifer. On the basis of the 1D and 2D longitudinal and transversal resistivity profiles as well as the existing boreholes data and groundwater quality, it seems that the freshwater-saltwater interface is moving forward inland. Furthermore, electromagnetic results combined to the Archie’s law show that the porosity values are at the ranges of 12–36 and 5–31% for the Quaternary outcropping and Pliocene formation, respectively.     

Development of an efficient surrogate model based on aquifer dimensions to prevent seawater intrusion in anisotropic coastal aquifers,case study: the Qom aquifer in Iran

Coastal aquifers are usually exposed to saltwater intrusion. Therefore, groundwater extracted from these aquifers should be regulated considering their dimensions and effective parameters. In this paper, optimum discharge from a large number of exploitation wells is evaluated according to variations of width, length, and anisotropy coefficient in the Qom aquifer near the salt lake in central Iran. First, the wells are divided into clusters to decrease the number of decision variables. Then, the location and discharge from each cluster is obtained using SEAWAT and charged system search (CSS) simulation–optimization model with the assumption of three-dimensional variable density flow. The maximum discharge considering various anisotropy rates is computed based on different values of lengths and widths of the aquifer. Finally, an M5-tree model is trained using the obtained samples to derive a linear relationship between input and output data. Based on the results, for various ranges of width and length of an aquifer with impermeable boundaries, different linear equations for optimum discharge are obtained. Also, it was found that for an aquifer with a small width, the critical discharge is a function of the length while the effect of the boundaries is negligible. Sensitivity analysis of the anisotropy coefficient reveals that with increasing the anisotropy rate, thickness and slope of the transition zone decrease and as the maximum discharge increases consequently. However, the sensitivity of the discharge to anisotropy rate is not remarkable. A comparison between the results of this study with those of the analytical method based on sharp interface assumption is carried out. For the critical condition, the best agreement between analytical equation (\(\overline {L} =0.87\overline {W} +0.62\)) and proposed method (\(\overline {L} =0.83\overline {W} - 1.41\)) is achieved for the anisotropic aquifer when the 50% isochlor is assumed as the measure of salt water intrusion.     

Mapping groundwater contamination risk using GIS and groundwater modelling. A case study from the Gaza Strip, Palestine

Increasing pressure on water resources worldwide has resulted in groundwater contamination, and thus the deterioration of the groundwater resources and a threat to the public health. Risk mapping of groundwater contamination is an important tool for groundwater protection, land use management, and public health. This study presents a new approach for groundwater contamination risk mapping, based on hydrogeological setting, land use, contamination load, and groundwater modelling. The risk map is a product of probability of contamination and impact. This approach was applied on the Gaza Strip area in Palestine as a case study. A spatial analyst tool within Geographical Information System (GIS) was used to interpolate and manipulate data to develop GIS maps of vulnerability, land use, and contamination impact. A groundwater flow model for the area of study was also used to track the flow and to delineate the capture zones of public wells. The results show that areas of highest contamination risk occur in the southern cities of Khan Yunis and Rafah. The majority of public wells are located in an intermediate risk zone and four wells are in a high risk zone.     

Groundwater beneath the urban area of Khan Younis City,southern Gaza Strip (Palestine): assessment for multi-domestic purposes

The present study research investigation is aimed to assess the groundwater quality for the urban area in Khan Younis City, southern Gaza Strip, for multi-domestic purposes. The physicochemical analysis of the groundwater wells shows the major ions in the order of Na⁺ > Mg²⁺ > Ca²⁺ > B³⁺ > K⁺ and Cl^? > HCO₃ ^? > SO₄ ^2? > NO₃ ^2? > F^? > PO₄ ^3?. Groundwater quality is classified as very hard-brackish water type. Ninety-five percent of the wells are classified as saline water type with high NO₃ ^2? concentrations. Based on water quality index (WQI), the groundwater falls into one of three categories: fair water (10%), poor water (15%), very poor (45%), and worst (30%). The high WQI values are because of high Na⁺, Cl^?, SO₄ ^2?, and NO₃ ^2? concentrations, while synthetic pollution index (SPI) values indicate that most about 80% of the wells are seriously polluted. Langelier Saturation Index (LSI) indicates that most of data are either slightly scale forming or corrosive water or slightly corrosive but non-scale forming, and 75% of the wells are suitable for construction purposes (have SO₄ ^2? concentrations <300 mg/L). The groundwater reaches alarming situation, where almost chemically unsuitable for drinking purposes and the water to be used after proper treatment such as desalination.     

水力屏障和截渗墙在海水入侵防治中的数值模拟研究

基于SEAWAT-2000程序构建室内二维砂箱试验中咸水入侵数值模型,利用该模型分析了针对不同补给井井位、补给井流量、截渗墙位置及贯穿深度等多种情景下的咸淡水界面运移规律。二维砂箱实验模拟结果表明,当注水井位于盐水楔前锋附近,距咸水边界40 cm、砂箱顶部边界5 cm处时,注水井工程措施能达到最佳海水入侵驱退效果,回退系数达21.5%。当截渗墙布设于距咸水边界10 cm处,贯穿深度为35 cm时,截渗墙工程措施能达到最佳海水入侵驱退效果,回退系数达81.1%。在此基础上,结合实际场地条件,构建山东龙口地区滨海含水层中某典型二维剖面的海水入侵数值模型,探讨了不同截渗墙布设情景模式下海水入侵状况。模拟结果表明,当截渗墙布设于距海岸线600 m处,贯穿深度为18 m时,截渗墙工程措施达到最佳海水入侵驱退效果,回退系数达28.4%。研究结果揭示了补给井井位、补给井流量、截渗墙位置及贯穿深度等因素对咸淡水界面运移规律的影响,可为场地条件下滨海含水层海水入侵防治中的工程管理措施优化提供参考依据。     

Using radon-222 to study coastal groundwater/surface-water interaction in the Crau coastal aquifer (southeastern France)

Radon has been used to determine groundwater velocity and groundwater discharge into wetlands at the southern downstream boundary of the Crau aquifer, southeastern France. This aquifer constitutes an important high-quality freshwater resource exploited for agriculture, industry and human consumption. An increase in salinity occurs close to the sea, highlighting the need to investigate the water balance and groundwater behavior. Darcy velocity was estimated using radon activities in well waters according to the Hamada “single-well method” (involving comparison with radon in groundwater in the aquifer itself). Measurements done at three depths (7, 15 and 21 m) provided velocity ranging from a few mm/day to more than 20 cm/day, with highest velocities observed at the 15-m depth. Resulting hydraulic conductivities agree with the known geology. Waters showing high radon activity and high salinity were found near the presumed shoreline at 3,000 years BP, highlighting the presence of ancient saltwater. Radon activity has also been measured in canals, rivers and ponds, to trace groundwater discharges and evaluate water balance. A model of the radon spatial evolution explains the observed radon activities. Groundwater discharge to surface water is low in pond waters (4 % of total inputs) but significant in canals (55 l/m²/day).     

Groundwater beneath the urban area of Khan Younis City,southern Gaza Strip (Palestine): hydrochemistry and water quality

Arabian Journal of Geosciences - This paper presents the results of a groundwater hydrochemical assessment in the urban area of Khan Younis City in the southern Gaza Strip by measuring its...     

The deuterium and oxygen-18 isotopic composition of the groundwater in Khan Younis City,southern Gaza Strip (Palestine)

The environmental isotopes such as deuterium and oxygen-18 and the deuterium excess values have been used to assess groundwater recharge sources and their dynamics in Khan Younis City in the Gaza Strip in Palestine. Three isotopic lines for the relationship between δ²H and δ¹⁸O were used in the assessment. These lines are the global meteoric water line, the local meteoric water line and the groundwater evaporation line. The δ²H, δ¹⁸O and D-excess values indicate that deuterium and oxygen-18 isotopes originated in the groundwater from groundwater mixing with rainfall and other water sources; the groundwater in the area recharged from rainfall from a distant source that came from the Mediterranean Sea and from other sources such as wastewater, irrigation return flow and saline water.