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.     

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.     

Management options for a multipurpose coastal aquifer in Oman

Using MODFLOW 2005, this study numerically evaluated the effects of managed aquifer recharge (MAR) using treated wastewater (TWW) in managing the Al-Khawd coastal aquifer northeast of Oman. Our primary objective is to increase the urban water supply and to sustain the aquifer service with the lowest possible damage to the aquifer. A number of managerial scenarios were simulated and progressively developed to reduce seawater intrusion and outflow of the groundwater to the sea. An economic analysis was conducted to characterize the trade-off between the benefits of MAR and seawater inflow to the aquifer under increased abstraction for domestic supply. The results show that by managing irrigation wells and relocating public wells in conjunction with MAR practices, the abstracted volume for drinking purposes could be doubled. Even though injection of TWW is more expensive (due to the injection cost), it was observed to result in greater benefits. The results indicate that managing the aquifer would produce a net benefit ranging from $8.22 million (scenario 7) to $15.21 million (scenario 4) compared to $1.57 million with the current practice. In conclusion, MAR using TWW is a feasible solution to develop water resources in arid regions, and the best scenario depends on the decision maker’s preference when weighing the benefits of MAR and the level of damage to the aquifer. MAR could help manage stressed aquifer systems in arid zones to maximize the benefit of using the water for domestic purposes while minimizing the damage to the aquifer.     

Effect of irrigation pumpage during drought on karst aquifer systems in highly agricultural watersheds: example of the Apalachicola-Chattahoochee-Flint river basin,southeastern USA

In the Apalachicola-Chattahoochee-Flint (ACF) river basin in Alabama, Georgia, and Florida (USA), population growth in the city of Atlanta and increased groundwater withdrawal for irrigation in southwest Georgia are greatly affecting the supply of freshwater to downstream regions. This study was conducted to understand and quantify the effect of irrigation pumpage on the karst Upper Floridan Aquifer and river–aquifer interactions in the lower ACF river basin in southwest Georgia. The groundwater MODular Finite-Element model (MODFE) was used for this study. The effect of two drought years, a moderate and a severe drought year, were simulated. Comparison of the results of the irrigated and non-irrigated scenarios showed that groundwater discharge to streams is a major outflow from the aquifer, and irrigation can cause as much as 10 % change in river–aquifer flux. The results also show that during months with high irrigation (e.g., June 2011), storage loss (34 %), the recharge and discharge from the upper semi-confining unit (30 %), and the river–aquifer flux (31 %) are the major water components contributing towards the impact of irrigation pumpage in the study area. A similar scenario plays out in many river basins throughout the world, especially in basins in which underlying karst aquifers are directly connected to a nearby stream. The study suggests that improved groundwater withdrawal strategies using climate forecasts needs to be developed in such a way that excessive withdrawals during droughts can be reduced to protect streams and river flows.     

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.     

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.     

Numerical simulation of sea water intrusion near Beihai, China

 A leaky aquifer system occurs in the coastal plain near Beihai, China. Seawater intrusion into the confined aquifer took place along the northern coast. Chloride concentrations at some observation wells increased steadily from 1988 and were at their peak in 1993. A quasi-three-dimensional element model has been developed to simulate the spatial and temporal evolution of hydraulic heads and chloride concentrations of the groundwater near the northern coast. The simulation model was based on the transition zone approach, which requires simultaneous solution of the governing water flow and solute transport equations. An irregular grid of a quadrangle was used to discretize the flow domain. Various aquifer parameters were verified with the numerical model in order to obtain satisfactory matches between computed values and observed data from an investigation. Three pumpage schemes were designed to use the calibrated model for prediction of future changes in water levels and chloride concentrations in groundwater in the study area. Results show that seawater intrusion would worsen in the confined aquifer if the current rates of groundwater pumpage continue. The alternative, to eliminate pumpage in the intruded area and to moderate pumpage rates from water supply wells far from the seashore, may limit seawater intrusion significantly and is considered attractive in the area. Received: 27 September 1999 · Accepted: 27 December 1999     

Intervention scenarios to manage seawater intrusion in a coastal agricultural area in Oman

The Batinah coastal plain in northern Oman has experienced a severe deterioration of groundwater quality due to seawater intrusion as a result of excessive groundwater abstraction for agricultural irrigation. Upgrading all farms to fully automated irrigation technology based on soil moisture sensors may significantly reduce the water demand and lead to recovering groundwater levels. This study compares the effects of smart irrigation technology, recharge dams, and a combination of both on seawater intrusion in the coastal aquifer of the Batinah. A groundwater flow and transport model is used to simulate the effect of reduced pumping rates on seawater intrusion for various intervention scenarios over a simulation period of 30 years, and an economic analysis based on cost-benefit analysis is conducted to estimate the potential benefits. Results indicate that a combination of smart irrigation and recharge dams may prevent further deterioration of groundwater quality over the next 30 years. In conjunction with increased efficiency, this combination also generates the highest gross profit. This outcome shows that the problem of seawater intrusion needs to be tackled by a comprehensive, integrated intervention strategy.     

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.     

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.     

Managed aquifer recharge (MAR) by the construction of subsurface dams in the semi-arid regions: A case study of the Kalangi river basin, Andhra Pradesh

Overuse of groundwater in coastal areas, due to high population and agricultural activity results in seawater intrusion into the coastal aquifer. This paper presents the control measures taken to manage aquifer recharge (MAR) and also to overcome the problem of seawater intrusion into the coastal aquifer along the Kalangi river, Nellore district of Andhra Pradesh, India having connectivity with Pulicat (saltwater) lake estuary. Due to overexploitation of groundwater and less rainfall in past years, adjacent seawater has started intruding in the Kalangi river sub-surface and deteriorating groundwater quality up to 11.6 km from the confluence of the river with Pulicat lake. To prevent this situtation, subsurface dams were constructed in traditional manner using local earth material in three different places across the Kalangi river near Sullurpet town. The water storage capacities calculated after the sub-surface dams’ construction are 1.28 mcft at GK Engineering College, 6.23 mcft at Challamagudi and 3.143 mcft at Holy Cross School sites. The Holy Cross School sub-surface dam is the first full scale dam-cum-check dam constructed to prevent salt water intrusion in the Kalangi river at Sullurpet, Nellore district, Andhra Pradesh. At the Kalangi river estuary portion (at the mouth of sea) a groyne was reconstructed over old groyne site with the introduction of clay bed and wooden sheet piles at down stream. Apart from prevention of sea water entry into Kalangi river sub-surface (during seasons) the groyne top level was raised to prevent mixing of high sea water tides with fresh water and ensuring additional storage of fresh water at upstream side. The reconstructed groyne was serving the purpose of obstructing the surface seawater entry in the Kalangi river and water quality has improved in the river as well as in the wells. After construction of sub-surface dam, as per the Simpson ratio classification, there is substantial improvement of water quality in the SHAR infiltration well situated near the Holy Cross School sub-surface dam.     

Groundwater management by riverbank filtration and an infiltration channel: the case of Obrenovac, Serbia

The Vi? Bare groundwater source is used to supply water to the population and industry of Obrenovac, one of the municipalities of Belgrade (the capital of Serbia). It is a typical riverbank filtration site; exploitation is performed through 30 drilled wells and two radial wells located in the meander of the Sava River. The established hydraulic connection between the river and tapped aquifer is so great that the river regime has a dominant influence on the aquifer. As a consequence of this, water-delivery reduction occurs in the dry months (summer–autumn), when the population needs water the most. Based on the data associated with the river’s gauges, precipitation, quantity of pumped water and groundwater-level fluctuation, a simulation of the groundwater regime for non-steady-state flow conditions has been undertaken through a numerical model. To overcome problems of water shortage during the dry season, the possibility of artificial recharge using an infiltration channel, made up of two connected parts, was analyzed. During the dry months, 80 % of the wells receive water partly from the infiltration channel. In this way, possibilities for extracting additional water are created. The application of this concept is discussed.     

Evaluating freshwater lens morphology affected by seawater intrusion using chemistry-resistivity integrated technique: a case study of two different land covers in Carey Island, Malaysia

Freshwater lenses are vital to small island communities but are susceptible to seawater intrusion due to the physical changes in the shoreline land cover. The effect of seawater intrusion and irrigation water on a coastal unconfined aquifer beneath naturally preserved mangrove and deforested mangrove-barren belt was investigated in Carey Island. Analysis of the total dissolved solids (TDS) and earth resistivity (ER) using a geochemistry-electrical integrated technique gave a TDS–ER relationship capable of predicting freshwater lens morphology affected by sea-irrigation water. The study result shows freshwater was fourfold thicker in close proximity of the mangrove forest than the mangrove barren area; the further the shoreline from the mangrove thickest section, the less vulnerable was the seawater intrusion and the more fresh the irrigation water, and hence the greater the freshwater availability potential.     

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.     

山东烟台夹河中、下游地区海水入侵三维水质数值模拟研究

建立了三维变密度对流弥散水质数学模型来研究山东省烟台夹河中、下游地区咸淡水界面的运移规律。以四面体为基本离散单元 ,推导出三维海水入侵变密度水质模型求解的数值方法 ,其中水流方程求解时运用了迦辽金有限单元法。溶质运移方程求解时运用了欧拉拉格朗日混合方法 ,将对流项与弥散项分离 ,用传统迦辽金有限元方法求解弥散项 ;采用自适应MOC MMOC法求解对流项 ,以消除人工过量和数值弥散。根据地下水的潮汐效应观测信息 ,确定了含水系统的海底延伸边界 ;利用该地区地下水水头及水质长观资料识别了模型的水文地质参数 ,探讨了夹河地区海水入侵的原因 :认为夹河下游地区滨海地带地下水过量开采是造成烟台地区海水入侵的主要原因。此外 ,海水随潮定期地倒灌进入夹河 ,通过局部岩性天窗侵入淡水含水层加剧了沿夹河河床两侧地下水的咸化。同时还预测了几种情况下地下水的水质演化趋势 ,为防止和减轻夹河地区海水入侵提供合理、科学的依据。     

Characterization of salt-water intrusion in the lower Esino Valley, Italy using a three-dimensional numerical model

A seawater-intrusion study was conducted at an oil-refinery site located on the coast in the lower Esino Valley, Italy. A steady-state density-dependent flow model was used in order to understand the position of the freshwater/salt-water interface, as influenced by the hydrogeologic structure and the presence of industrial activities and a river. Collected data and model results showed that in a large part of the area, the salt-water interface is steep and can penetrate only a few meters inland. On the other hand, close to the river mouth, seawater represents the main saline source for the aquifer. The river, in connection with the sea, can enhance seawater encroachment into the coastal aquifer; a long-term survey of river level and chloride concentrations in groundwater is recommended to further improve the physical model and to obtain a better calibration. At the refinery site, two “secondary” sources of saline water were identified and were demonstrated to have had a great influence on the presence of brackish waters in the unconfined aquifer: leakage from the fire-extinguishing system (network of pipes containing seawater) and rough sea events. This confirmed that groundwater contamination by chloride can result from means other than 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.     

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.     

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 resource degradation in coastal plains: The example of the Cecina area (Tuscany – Central Italy)

The paper describes the degradation of the groundwater resources in the Cecina area, where seawater intrusion, B contamination and NO₃ pollution are all affecting the heavily exploited Pleistocene aquifer. Over-pumping has brought water levels to about 0 m.a.s.l. as far as about 7 km from the shore line, thereby promoting the seawater intrusion. The intrusion, which is characterized by cation exchange phenomena and Ca–Cl type waters, enters the plain mostly through the shallower horizons. The saline front, which advanced from 0.5 to 1 km in 4 a, has by now reached the foot of the hills to the east of the town, where it is also affecting wells of the local aqueduct. Boron contamination, linked to past discharge of industrial waste transported downstream by the river, reached concentrations as high as 3.5 mg/L in the mid-1980s. Although a decreasing trend is now under way, B content is still close to 1 mg/L. The presence of high NO₃, which, together with the seawater intrusion, represents a major issue for groundwater management in the area, is linked to the widespread utilization of fertilizers. Nitrate concentration, which reaches a maximum of about 300 mg/L in the shallow aquifer horizons and then decreases rather regularly with depth, is strongly influenced by precipitation. However, irrigation also contributes significantly to transporting the NO₃ contamination to depth, as clearly shown by δ¹⁸O data. The severe decline in the quality of the groundwater resource in the Cecina area is further compounded by an overall decrease in water availability in the region of Tuscany, as evidenced by long-term monitoring of precipitation and fluvial discharge.