Seawater intrusion in Syrian coastal aquifers, past, present and future, case study

Overexploitation of shallow aquifers on the Syrian coast, north of Latakia (Damsarkho) for irrigation and tourism has caused an intrusion of seawater. The seawater intrusion into this aquifer has been presented by a three-dimensional finite element model using the FEFLOW numerical code. This conceptual model is based on field and laboratory data collected during the period 1966–2003. Meteoric infiltration and flows from the adjoining carbonate aquifer recharge the aquifer; natural outflow occurs through a diffuse flow into the sea; and artificial outflow occurs through intensive extraction of groundwater from wells. Water exchanges in the aquifer occur naturally (leakage) and artificially (multi-screened wells). The model was calibrated for transient conditions. The model helped establish that seawater intrusion is essentially due to withdrawals near the coast during the irrigation season and that it mainly occurs in the Damsarkho plain. The effects of hypothetical aquifer exploitation were assessed in terms of salt budget.     

Hydrochemical and statistical studies of the groundwater salinization in Mediterranean arid zones: case of the Jerba coastal aquifer in southeast Tunisia

Coastal aquifers are considered as major sources for freshwater supply worldwide, especially in arid zones. The weak rainfall as well as the intensive extraction of groundwater from coastal aquifers reduce freshwater budget and create local water aquifer depression, causing both seawater intrusion and a threat to groundwater. This phenomenon was observed in the Jerba Island which is located in southeast Tunisia. Jerba??s unconfined aquifer shows high values of groundwater salinity reaching, locally, 17?g/l and a strong contrast between some zones of the aquifer. High pumping rates and weak recharge disturb the natural equilibrium between fresh and saline water causing water salinization in most areas of the island. This study aims at establishing the salinity map of the aquifer and identifying the origin of groundwater salinization. The salinity map shows that zones characterized by low groundwater salinity are located in the center of the study area. High groundwater salinities are observed near the coast and in some parts having low topographic and piezometric levels. Groundwater geochemical characterization, and Br/Cl and Na/Cl ratios suggest that the origin of abnormal salinity is seawater intrusion. Considering groundwater salinity values and Br concentrations, a seawater intrusion map is established. It shows that many areas of the unconfined aquifer are contaminated by mixed groundwater and seawater. The statistical analysis demonstrates that high mineralization of the groundwater is due to gypsum and carbonate dissolution coupled with the mixed groundwater and seawater in many areas.     

Modeling of seawater intrusion in a coastal aquifer of Karaburun Peninsula,western Turkey

Seawater intrusion is a major problem to freshwater resources especially in coastal areas where fresh groundwater is surrounded and could be easily influenced by seawater. This study presents the development of a conceptual and numerical model for the coastal aquifer of Karareis region (Karaburun Peninsula) in the western part of Turkey. The study also presents the interpretation and the analysis of the time series data of groundwater levels recorded by data loggers. The SEAWAT model is used in this study to solve the density-dependent flow field and seawater intrusion in the coastal aquifer that is under excessive pumping particularly during summer months. The model was calibrated using the average values of a 1-year dataset and further verified by the average values of another year. Five potential scenarios were analyzed to understand the effects of pumping and climate change on groundwater levels and the extent of seawater intrusion in the next 10 years. The result of the analysis demonstrated high levels of electrical conductivity and chloride along the coastal part of the study area. As a result of the numerical model, seawater intrusion is simulated to move about 420 m toward the land in the next 10 years under “increased pumping” scenario, while a slight change in water level and TDS concentrations was observed in “climate change” scenario. Results also revealed that a reduction in the pumping rate from Karareis wells will be necessary to protect fresh groundwater from contamination by seawater.     

Groundwater resource sustainability in the Wadi Watir delta, Gulf of Aqaba, Sinai, Egypt

The Wadi Watir delta, in the arid Sinai Peninsula, Egypt, contains an alluvial aquifer underlain by impermeable Precambrian basement rock. The scarcity of rainfall during the last decade, combined with high pumping rates, resulted in degradation of water quality in the main supply wells along the mountain front, which has resulted in reduced groundwater pumping. Additionally, seawater intrusion along the coast has increased salinity in some wells. A three-dimensional (3D) groundwater flow model (MODFLOW) was calibrated using groundwater-level changes and pumping rates from 1982 to 2009; the groundwater recharge rate was estimated to be 1.58?×?10⁶ m³/year. A variable-density flow model (SEAWAT) was used to evaluate seawater intrusion for different pumping rates and well-field locations. Water chemistry and stable isotope data were used to calculate seawater mixing with groundwater along the coast. Geochemical modeling (NETPATH) determined the sources and mixing of different groundwaters from the mountainous recharge areas and within the delta aquifers; results showed that the groundwater salinity is controlled by dissolution of minerals and salts in the aquifers along flow paths and mixing of chemically different waters, including upwelling of saline groundwater and seawater intrusion. Future groundwater pumping must be closely monitored to limit these effects.     

Control of sea-water intrusion by salt-water pumping: Coast of Oman

A shallow alluvial coastal aquifer in the Batinah area of Oman, with sea-water intrusion that extends several kilometres inland, has been studied experimentally, analytically and numerically. The water table is proved to have a trough caused by intensive pumping from a fresh groundwater zone and evaporation from the saline phreatic surface. Resistivity traverses perpendicular to the shoreline indicated no fresh groundwater recharge into the sea. Using an analytical Dupuit-Forchheimer model, developed for the plain part of the catchment, explicit expressions for the water table, sharp interface location and stored volume of fresh water are obtained. It is shown that by the pumping of salt water from the intruded part of the aquifer, this intrusion can be mitigated. Different catchment sizes, intensities of fresh groundwater pumping, evaporation rates, water densities, sea level, incident fresh water level in the mountains and hydraulic conductivity are considered. SUTRA code is applied to a hypothetical case of a leaky aquifer with line sinks modeling fresh water withdrawal and evaporation. The numerical code also shows that pumping of saline water can pull the dispersion zone back to the shoreline.     

Geophysical characterization of the complex dynamics of groundwater and seawater exchange in a highly stressed aquifer system linked to a coastal lagoon (SE Spain)

Anthropogenic pressure influences the two-way interactions between shallow aquifers and coastal lagoons. Aquifer overexploitation may lead to seawater intrusion, and aquifer recharge from rainfall plus irrigation may, in turn, increase the groundwater discharge into the lagoon. We analyse the evolution, since the 1950s up to the present, of the interactions between the Campo de Cartagena Quaternary aquifer and the Mar Menor coastal lagoon (SE Spain). This is a very heterogeneous and anisotropic detrital aquifer, where aquifer–lagoon interface has a very irregular geometry. Using electrical resistivity tomography, we clearly identified the freshwater–saltwater transition zone and detected areas affected by seawater intrusion. Severity of the intrusion was spatially variable and significantly related to the density of irrigation wells in 1950s–1960s, suggesting the role of groundwater overexploitation. We distinguish two different mechanisms by which water from the sea invades the land: (a) horizontal advance of the interface due to a wide exploitation area and (b) vertical rise (upconing) caused by local intensive pumping. In general, shallow parts of the geophysical profiles show higher electrical resistivity associated with freshwater mainly coming from irrigation return flows, with water resources mostly from deep confined aquifers and imported from Tagus river, 400 km north. This indicates a likely reversal of the former seawater intrusion process.     

Salinisation des nappes côtières : cas de la nappe nord du Sahel de Sfax,Tunisie

The intensive agricultural and economic activities induce the increase of the risk of groundwater quality degradation through high groundwater pumping rates. The salinization and contamination are the main sources of this pollution, especially in coastal aquifers. The explanation of the origin of salinity for the shallow aquifer of Northern Sahel of Sfax was analysed by a chemical study of the groundwater main compounds. The partitioning of groundwaters into homogenous groups is undertaken by graphical techniques, including a Stiff pattern diagram, an expanded Durov diagram and several binary diagrams. The study indicates the presence of various salinization processes. In the recharge area, salinization is the result of dissolution/precipitation of the aquifer formation material (group I). The irrigation water return and the intensive pumping have been identified as major sources of salinization in the south by direct cation exchange and mixing reactions (groups II and III). The anomaly of high groundwater salinity observed near the Hazeg zone was explained by the presence of a seawater intrusion in this area. This hypothesis is related to the high chloride concentration, to the presence of inverse cation exchange reactions (group IV), and to the piezometric level inferior to sea level. To cite this article: R. Trabelsi et al., C. R. Geoscience 337 (2005).     

大陆岛地下水动力学特征—以湛江东海岛为例

东海岛是一具有独特水文地质条件的大陆岛,浅层含水层与大陆以浅海湾相隔,中深层承压水含水层与大陆地下水系统紧密相连。为了深入地认识大陆岛地下水水动力学特征,以湛江东海岛为例,阐述了其水文地质条件,并分析了东海岛浅、中、深层地下水的流场和动态特征。分析结果表明,东海岛为一个典型且独特的大陆岛,岛内和大陆的部分浅层含水层由湛江湾相隔,岛内中、深层含水层和大陆中、深层含水层通过湛江湾相连,且具有统一的水位分布,并保持着密切的水力联系,岛内中、深层地下水由南向北径流补给湛江市区的降落漏斗中心;滨海及海水区域浅层含水层及其下伏的粘土层构成了防止海水入侵中、深层地下水的保护层;浅层地下水流场基本保持天然状态,水位动态特征主要为入渗径流型,水位变化与降雨量相关;中、深层地下水流场以人工流场为主,地下水由南向北径流,水位动态类型主要为开采动态型,水位变化主要受到开采量变化的影响;在近海岸地区,地下水动态表现为潮汐效应型,在潮汐作用下,地下水位动态具有周期性。     

Hydrochemical changes in a small tropical island’s aquifer: Manukan Island, Sabah, Malaysia

Small islands groundwater are often exposed to heavy pumpings as a result of high demand for freshwater consumption. Intensive exploitation of groundwater from Manukan Island’s aquifer has disturbed the natural equilibrium between fresh and saline water, and has resulted increase the groundwater salinity and leap to the hydrochemical complexities of freshwater–seawater contact. An attempt was made to identify the hydrochemical processes that accompany current intrusion of seawater using ionic changes and saturation indices. It was observed that the mixing between freshwater–seawater created diversity in geochemical processes of the Manukan Island’s aquifer and altered the freshwater and seawater mixture away from the theoretical composition line. This explained the most visible processes taking place during the displacement.     

Monitoring and methods to analyse the groundwater quality degradation risk in coastal karstic aquifers (Apulia,Southern Italy)

A multi-methodological approach based on monitoring and spatio-temporal analysis of groundwater quality changes is proposed. The presented tools are simple, quick and cost-effective to give service to all sorts of users. The chief purpose of the monitoring network is the detection of the piezometric or potenziometric level in the aquifer. The spatial and multi-temporal analysis of usual chemical and physical data provides both an assessment of the spatial vulnerability of the aquifer to seawater intrusion, defining a salinity threshold between fresh groundwater and brackish groundwater and of the water quality trend in terms of salinity. The evaluation of the salinity trend or of salinity-correlated parameters highlights the effects of groundwater mismanagement. The multiparameter logging provides a rapid groundwater quality classification for each well. The whole approach allows evaluating the effects of current management criteria and designing more appropriate management targets. The Apulian karstic coastal aquifers have been selected as a case study (Southern Italy). Three types of aquifer zones can be distinguished: (1) areas with low vulnerability to seawater intrusion, (2) areas with high vulnerability and (3) areas with variable vulnerability in which the salt degradation largely depends on the ability to manage the well discharge. The water quality degradation caused by seawater intrusion appears to be a combined effect of an anomalous succession of drought periods observed from about 1980 onwards and increased groundwater pumping, particularly during drought periods. A management criterion based on aquifer zones is proposed.     

Numerical simulation of the influence of management alternatives of a projected reservoir on a small alluvial aquifer affected by seawater intrusion (Almuñécar, Spain)

The present study concerns the application of a numerical approach to describe the influence of anthropogenic modifications in surface flows (operation of a projected reservoir) on the freshwater-seawater relationships in a downstream coastal aquifer which has seasonal seawater intrusion problems (River Verde alluvial aquifer, Almuñécar, southern Spain). A steady-state finite element solution to the partial differential equation governing the regional motion of a phreatic surface and the resulting sharp interface between fresh water and salt water was used to predict the regional behavior of the River Verde aquifer under actual surface flow conditions. The present model approximates, with simple triangular elements, the regional behavior of a coastal aquifer under appropriate sinks, sources, Neumann and open boundary conditions. A steady-state solution to this numerical approach has been shown to precisely calculate freshwater heads, saltwater thicknesses, and freshwater discharges along steeply sloping coasts. Hence, the adequate treatment and interpretation of the hydrogeological data which are available for the River Verde aquifer have been of main concern in satisfactorily applying the proposed numerical model. Present simulated conditions consider steady-state yearly averaged amounts of external supplies of fresh water in order to determine the influences of the projected Otívar reservoir on the further behavior of the River Verde coastal aquifer. When recharges occur at the coastline, essentially because of freshwater deficits due to groundwater overexploitation, a hypothesis of mixing for the freshwater-saltwater transition zone is made in order to still allow the model to continue calculating groundwater heads under the sea level, and, as a consequence, the resulting seawater intrusion and recharges of saltwater from the sea. Simulations show that a considerable advance in seawater intrusion would be expected in the coastal aquifer if current rates of groundwater pumping continue and a significant part of the runoff from the River Verde is channeled from the Otívar reservoir for irrigation purposes.     

Groundwater quality of porous aquifers in Greece: a synoptic review

Greece is dependent on groundwater resources for its water supply. The main aquifers are within carbonate rocks (karstic aquifers) and coarse grained Neogene and Quaternary deposits (porous aquifers). The use of groundwater resources has become particularly intensive in coastal areas during the last decades with the intense urbanization, tourist development and irrigated land expansion. Sources of groundwater pollution are the seawater intrusion due to over-exploitation of coastal aquifers, the fertilizers from agricultural activities and the disposal of untreated wastewater in torrents or in old pumping wells. In the last decades the total abstractions from coastal aquifers exceed the natural recharge; so the aquifer systems are not used safely. Over-exploitation causes a negative water balance, triggering seawater intrusion. Seawater intrusion phenomena are recorded in coastal aquifer systems. Nitrate pollution is the second major source of groundwater degradation in many areas in Greece. The high levels of nitrate are probably the result of over-fertilization and the lack of sewage systems in some urban areas.     

Salinity mapping of coastal groundwater aquifers using hydrogeochemical and geophysical methods: a case study from north Kelantan,Malaysia

Integrated hydrogeochemical and geophysical methods were used to study the salinity of groundwater aquifers along the coastal area of north Kelantan. For the hydrogeochemical investigation, analysis of major ion contents of the groundwater was conducted, and other chemical parameters such as pH and total dissolved solids were also determined. For the geophysical study, both geoelectrical resistivity soundings and reflection seismic surveys were conducted to determine the characteristics of the subsurface and groundwater contained within the aquifers. The pH values range from 6.2 to 6.8, indicating that the groundwater in the study area is slightly acidic. Low content of chloride suggests that the groundwater in the first aquifer is fresh, with an average concentration of about 15.8 mg/l and high geoelectrical resistivity (>45 ohm m). On the other hand, the groundwater in the second aquifer is brackish, with chloride concentration ranging from 500 mg/l to 3,600 mg/l and very low geoelectrical resistivity (<45 ohm m) as well as high concentration of total dissolved solids (>1,000 mg/l). The groundwater in the third aquifer is fresh, with chloride concentrations generally ranging from 2 mg/l to 210 mg/l and geoelectrical resistivity of greater than 45 ohm m. Fresh and saltwater interface in the first aquifer is generally located directly in the area of the coast, but, for the second aquifer, both hydrogeochemical and geoelectrical resistivity results indicate that the fresh water and saltwater interface is located as far as 6 km from the beach. The considerable chloride ion content initially suggests that the salinity of the groundwater in the second aquifer is probably caused by the intrusion of seawater. However, continuous monitoring of the chloride content of the second aquifer indicated no significant changes with time, from which it can be inferred that the salinity of the groundwater is not affected by seasonal seawater intrusion. Schoeller diagrams illustrate that sulphate concentrations of the groundwater of the second aquifer are relatively low compared to those of the recent seawater. Therefore, this result suggests that the brackish water in the second aquifer is probably from ancient seawater that was trapped within the sediments for a long period of time, rather than due to direct seawater intrusion.     

Evaluation of the impacts of climate changes on the coastal Chaouia aquifer, Morocco, using numerical modeling

The aquifer of the Chaouia Coast, Morocco constitutes an example of groundwater resources subjected to intensive and uncontrolled withdrawals in a semi-arid region. The analysis of the trends of precipitation and piezometric levels of the Chaouia coastal aquifer, with the use of moving averages, emphasized the impact of the climate on the groundwater resources of the system. The results showed that the periods 1977–1993 and 1996–2000 are characterized by a deficit in precipitation, although the precipitation increased slightly during the periods 1973–1977 and 1993–1996. Numerical modeling of the Chaouia aquifer showed that the groundwater resources of this system are less sensitive to the variations in precipitation. Severe degradation of the resource is related to intensive pumping during the periods of drought, which has forced abandonment of wells due to seawater intrusion.     

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     

Evaluation of groundwater quality in coastal areas: implications for sustainable agriculture

Seawater intrusion is a problem in the coastal areas of Korea. Most productive agricultural fields are in the western and southern coastal areas of the country where irrigation predominantly relies on groundwater. Seawater intrusion has affected agricultural productivity. To evaluate progressive encroachment of saline water, the Korean government established a seawater intrusion monitoring well network, especially in the western and southern part of the peninsula. Automatic water levels and EC monitoring and periodic chemical analysis of groundwater help track salinization. Salinization of fresh groundwater is highly associated with groundwater withdrawal. A large proportion of the groundwaters are classified as Na–Cl and Ca–Cl types. The Na–Cl types represent effects of seawater intrusion. The highest EC level was over 1.6 km inland and high Cl values were observed up to 1.2 km inland. Lower ratios of Na/Cl and SO₄/Cl than seawater values indicate the seawater encroachment. A linear relation between Na and Cl represents simple mixing of the fresh groundwater with the seawater. The saline Na–Cl typed groundwaters showed Br/Cl ratios similar to or less than seawater values. The Ca–HCO₃ type groundwaters had the highest Br/Cl ratios. Substantial proportions of the groundwaters showed potential for salinity and should be better managed for sustainable agriculture.     

Groundwater salinization processes in shallow coastal aquifer of Djeffara plain of Medenine, Southeastern Tunisia

Urban and industrial development and the expansion of irrigated agriculture have led to a drastic increase in the exploitation of groundwater resources. The over-exploitation of coastal aquifers has caused a seawater intrusion and has seriously degraded groundwater quality. The shallow coastal aquifer of the Djeffara plain, southeastern Tunisia constitutes an example of water resource suffering an intensive and uncontrolled pumping for irrigation. Intensive exploitation of the aquifer and climate aridity caused a decrease in piezometric level and an increase in salinity. According to the hydrochemical data (Cl⁻, SO₄ ²⁻, NO₃ ⁻, HCO₃ ⁻, Br⁻, Ca²⁺, Mg²⁺, Na⁺, K⁺) and the stable isotope composition (oxygen-18 and deuterium content), groundwater salinization in the investigated system is caused by three main processes: (i) salts dissolution especially in the central part of Jerba and around Medenine plain; (ii) evaporation process; and (iii) seawater intrusion which caused the increase in salinity in the peninsula of El Jorf, in Jerba and in the North of Ben Gardane.     

环境同位素特征对滨海岩溶地区海水入侵过程的指示意义

大连大魏家水源地位于中国北方典型滨海岩溶地区。近30年来,地下淡水的不合理开采造成的地下水位降落漏斗引发了严重的海水入侵。以大魏家水源地为研究对象,通过大量的水文地质调查和水化学及同位素采样测试分析,探讨海水入侵形成的水动力条件,通过分析滨海岩溶含水层中地下水主要水化学和多种同位素(δ2H-δ18O,δ34S,δ13C)组成特征,识别了海水入侵过程中发生的主要水文地球化学作用,并对其进行了定量模拟,从而阐明了岩溶含水层中的海水入侵机理。研究结果表明:大连大魏家海水入侵主要通道为大魏家地区存在的导水断裂、岩溶裂隙以及第四系松散地层。对δ2H-δ18O同位素的组成分析表明,研究区地下水主要来自大气降水补给,结合Cl-浓度分布,认为除海水入侵淡水含水层后增加了地下水中的盐分外,浅层地下水的蒸发也对地下水中盐分的累积起到了重要作用。根据不同水体中δ34SSO4,δ13CHCO3等同位素特征,结合水化学成分(如SO2-4,Cl-)分析认为,研究区微咸水和咸水并不是地下水淡水和海水简单混合而成。利用反向水文地球化学模拟揭示了控制滨海岩溶含水层中水化学演化的主要水文地球化学反应有方解石、蒙脱石和石膏的溶解作用,伊利石的沉淀作用以及Ca-Na离子交换作用,伴随着CO2的释放。     

Effects of the planned ephesus recreational canal on freshwater–seawater interface in the Selçuk sub-basin, İzmir-Turkey

An artificial water canal opening is planned between the Agean Sea and the historical Ephesus site for the sake of tourism in the Selçuk sub-basin. In order to predict the effects of the planned canal on freshwater–seawater interface and related contamination in the aquifer, 3-D numerical density dependent flow and solute transport simulations were carried out. The simulations included the pre-pumping and pumping periods without a canal and the prediction period in the presence of the canal. Chloride concentration comparisons of the results obtained from the pre-pumping period and the pumping period indicate that the freshwater-seawater interface in the aquifer has progressed inland due to artificial discharge in the sub-basin. Drawdown during the pumping period is about 15 cm. The planned canal opening could further lower the groundwater levels in the area and would change the groundwater flow directions in the first 4 years. Then the levels and flow directions will nearly recover. However, the canal opening could cause further seawater intrusion into the aquifer to the extent that groundwater would be unfit to use for irrigation after the seventh year of the canal opening in the irrigation cooperative II wells area and would be unfit to use for drinking purposes after the tenth year in the municipality wells area located at the south of the cooperative II wells. On the other hand, the cooperative I wells would not be effected by the opening of the canal.     

Regional hydrostratigraphy and groundwater flow modeling in the arsenic-affected areas of the western Bengal basin, West Bengal, India

The first documented interpretation of the regional-scale hydrostratigraphy and groundwater flow is presented for a ～21,000-km² area of the arsenic-affected districts of West Bengal [Murshidabad, Nadia, North 24 Parganas and South 24 Parganas (including Calcutta)], India. A hydrostratigraphic model demonstrates the presence of a continuous, semi-confined sand aquifer underlain by a thick clay aquitard. The aquifer thickens toward the east and south. In the south, discontinuous clay layers locally divide the near-surface aquifer into several deeper, laterally connected, confined aquifers. Eight 22-layer model scenarios of regional groundwater flow were developed based on the observed topography, seasonal conditions, and inferred hydrostratigraphy. The models suggest the existence of seasonally variable, regional, north–south flow across the basin prior to the onset of extensive pumping in the 1970s. Pumping has severely distorted the flow pattern, inducing high vertical hydraulic gradients across wide cones of depression. Pumping has also increased total recharge (including irrigational return flow), inflow from rivers, and sea water intrusion. Consequently, downward flow of arsenic contaminated shallow groundwater appears to have resulted in contamination of previously safe aquifers by a combination of mechanical mixing and changes in chemical equilibrium.