The Influence of Seawater on a Coastal Aquifer in an International Protected Area, Göksu Delta Turkey ()
1. Introduction
The Göksu Delta is formed by the Goksu River near the Southern part of Silifke, a town in Turkey’s Mediterranean region. Göksu Delta is an important wetland (15000 ha) where the Göksu River reaches the sea in the eastern of the town Tasucu-Ice1 [1].
The Mediterranean coastline stretching from the city of Taşucu to the Susanoğlu region is heavily populated due to urban developments (villas, apartments, complexes, and multi-storey buildings) in the last fifteen years, though most of these structures are occupied only in the summer season for vacation purposes. Because of an increased population influx from the surrounding cities, especially during the peak season (May to September), the population of this region has increased several times.
The Göksu Delta is not only an urban area but it is also surrounded by densely cultivated orchards (mostly citrus), traditional vegetable farms and greenhouse cultivations, where farming activities continue all year long due to the favorable climate.
In the Göksu Delta area, urban and agricultural expansions have caused an ever-growing need for fresh water. In this region, water supply for most municipalities, domestic use water for urban developments and irrigated water for agricultural activities are almost exclusively provided through hand dug or drilled wells. Therefore, water resources in the Göksu Delta area are subjected to intensive demands, stresses and pollution risks [2]. The expansion of irrigated agriculture induces the risk of groundwater quality degradation through high groundwater pumping rates and overexploitation of the aquifers, leading to seawater intrusion in the coastal aquifers.
Göksu Delta is an internationally important wetland due to its location being on a bird migration route. The Environmental Protection Department of the Ministry of Environment has declared the Göksu Delta as a Special Environmental Protection Zone to protect the area against pollution and exploitation, and to ensure that natural resources and cultural assets have a future. The delta is classified as a Wetland of International Importance according to the Ramsar Convention on Wetlands of International Importance. The Goksu Delta has also a special significance for being one of the few remaining areas in the world where sea turtles (Caretta caretta, Chelonia mydas) and blue crabs (Callinectes sapidus) lay their eggs [1,3].
Under steady-state conditions, a state of equilibrium is established between seawater and freshwater in Göksu coastal alluvial aquifer. However, owing to increasing demand for water, groundwater may be subjected to over-exploitation and the natural equilibrium is thus disturbed. This results in aggressive seawater mixing, which may even reach the inland aquifer.
The purpose of this study is to understand the seawater intrusion into groundwater, to investigate the pollution of groundwater by sea water and pictorially represent it using the geographic information system (GIS).
2. Site Description
The Göksu Delta is situated in the Mediterranean Sea region of the southeastern part of Turkey and extends from 36°15-36°25 of latitude north to 33°55-34°05 of longitude west. The Göksu Delta area is bounded by the Taurus Mountains on the northern side and by the Mediterranean Sea on the southern and eastern side. The southern portion of the Göksu Delta area is a delta plain made up of sediments from Göksu River. The Goksu River flow regime is strongly dependent on the seasonal rains and temperature. The average flow of Goksu River is 130 m3/s where it reaches the highest value during May.
In the Göksu Delta area, climate is characterized by hot and dry periods in the summer and by warm and wet periods in the winter, which is typical for the coastal zones around the Mediterranean Sea. The mean annual temperature in this area is 19°C. Showers start in October, and continue until mid April and the maximum rainfall occurs in December. The Göksu Delta area receives slightly higher than 607 mm of precipitation annually, and extended periods (i.e., 3-4 months) without precipitation are common The study area was located in the southern part of Miocene carbonate rocks of the Taurus Mountains Belt.
The oldest rock unit of the Göksu Delta is the Akdere Formation from the Paleozoic Age, which consists of marble, schist and quartzite. The Akdere Formation (middleupper Devonian) is generally found in the northern part of the study area. Kusyuvası Formation (Middle Trias) consists of limestone. The Tokmar Formation (upper Jura-lower Cretaceous) is found in the western part of delta and contains dolomite and limestone.
Tertiary units are composed of lower-middle Miocene Karaisali formation and middle-upper Miocene Kuzgun formation (Figure 1). Tertiary rocks consist of a succession of marine, lacustrine and fluvial deposits, which display transitional characteristics both vertically and areally in the study area.
The Quaternary basin-fill deposits are a heterogeneous mixture of metamorphic and sedimentary rock detritus ranging from clay to boulder size. The mixture includes stream alluvium, stream-terrace deposits, fan deposits, delta deposits, shore deposits. The basin-fill deposits vary greatly in lithology and grain-size, both vertically and areally. Accordingly, the hydraulic properties of these deposits can differ greatly over short distances, both laterally and vertically.
The alluvial aquifer consists of a heterogeneous mixture of gravel, sand, silt, clay and sandy-clay. Conceptually the aquifer system in the delta is an unconfined aquifer and more than 500 m thick. Recharge occurs by means of precipitation and infiltration at the top of the delta. The recharge rate is determined as 35.23 mm by using Visual HELP model [4].
3. Materials and Methods
Salty water can be detected directly in observation wells by electrical electrodes or by sampling, or indirectly by geoelectrical methods [5]. In this study, the direct method by sampling is used. AquaChem software is used and it allows comparing a sample (Sea water sample) to another multiple samples (all groundwater samples). This comparation uses a linear regression algorithm to generate the correlation coefficient and the Euclidean distance between Sea water and all other samples.
For chemical analysis, a total of 16 water samples from the Göksu Delta (9 from alluvial aquifer, 4 from limestone aquifer, 1 from Paradeniz Lake, 1 from Akgöl and 1 from Mediterranean Sea) were obtained from 2006 to 2008 by four separate sampling campaigns at the sampling points shown in Figure 1. Table 1 summarizes the chemical analysis results for water samples collected from the Göksu Delta and the Table 2 shows the chemical characters of waters. The pH, temperature (T), electrical conductivity (EC) and salinity (sal) were measured at 25 sampling points in four different periods between 2006 and 2008 (Table 3).
Water samples obtained from the wells are from various depths because the wells in the area vary greatly in depth. Average well depth is 5 m for hand dug wells (showed with DSI) and 30-35 m for drilled wells (shown with ME). Electrical conductivity (EC), temperature (T), salinity (sal) and pH were monitored during pumping, and samples were collected only when values were stabilized or after at least three well volumes had been purged. Measurements of EC and pH were made in the field using a pH/Cond 340i WTW meter. For the pH measurements the electrode was calibrated against pH buffers at each location.
Cations were analyzed by inductively coupled plasma (ICP) and anions by ion chromatography (IC). SiO2 was analyzed mainly by visible spectrophotometer. Bicarbonates were determined by titration in the laboratory. Samples were analyzed in the laboratory of General Directorate of Mineral Research and Exploration (MTA) of Turkey in Ankara.
Followed by water quality analysis, a vector-based GIS software package MapInfo was used to map, query, and analyze the data in this study. GIS is an effective tool for storing large volumes of data that can be correlated spatially and retrieved for the spatial analysis and integration to produce the desirable output. GIS has been used by scientists of various disciplines for spatial queries, analysis and integration for the last three decades [6]. GIS is a powerful tool and has great promise for use in environmental problem solving. Most environmental problems have an obvious spatial dimension and spatially distributed models can interact with GIS [7]. Troge [8] reported that this computer-based tool has allowed successful integration of water quality variables into a comprehensible format.
4. Results and Discussion
The groundwater samples collected from the Göksu Delta are colorless, odorless and free from turbidity. In Table 1, the results of ion concentrations of groundwater for July 2006 are presented. The Table 2 shows the chemical types of groundwater from different aquifers and places.
Paradeniz Lake is a saltwater lagoon connected to the sea and the chemistry of Paradeniz water is similar to the Mediterranean Seawater. The EC values of sea water and lake water are 53000 µS/cm and 48000 µS/cm, respectively.
The types of groundwater from alluvial aquifer are Ca-Mg-HCO3, but in the region where sea water intrusion occurs it changes and the Na and Cl ions are added in the types of groundwater (Table 2). The similar occurrence can be observed in the groundwater from the limestone aquifer. The chemical types of groundwater are Ca-HCO3 or Ca-Mg-HCO3 and they change to CaNa-Mg-HCO3-Cl in the vicinity of seawater mixing regions.
Table 1. Results of the chemical analyses. Concentrations are in milligram per liter.
Table 2. The chemical characters of ground waters in the Göksu delta.
With very few exceptions, sample points fall within the mixing field of the Piper diagram (Figure 2). Three regions are distinguished in the diagram. First region where the samples ME3 and ME 28 are placed show the aquifer region where the seawater intrusion does not occurred. The second place, the central of the Piper diagram where the samples ME1 and ME 20 are shows the aquifer region where the sea water impact begins. The third place is very close to the Sea water sample, where the samples ME 12 and ME 23 are shows the aquifer region where the sea water mixing with sea water occurs. The samples ME12 and ME23 are very close to the Sea and Paradeniz waters.
The other indicators for the detection of sea water influence are T, EC and sal values. Table 3 contains the results of pH, T, EC and sal measurements for four periods between 2006 and 2008. By using MapInfo, many
thematic maps are produced and by overlaying these maps sea water affected areas are determined.
The extent of saltwater affect is reflected by the distribution of salinity and electrical conductivity as an equiconcentration map (Figures 4-9). The periodical sal and EC measurements are used to show the development of saltwater mixing in space for four selected times.
The thematic map of groundwater temperature shows that in the region the temperature is increased in the places where the sea water mixing occurs. These regions are the southern part of Delta and in the vicinity of Sultanoğlu and Altınkum towns. In these regions the groundwater temperature increased until to 25-26°C (Figure 3).
EC values vary between 434 and 3460 µS/cm and sal values between 0-1.7.
Table 4. The results of statistical (comparation) analyses.
Figure 3. Temperature of groundwater for 2006.
Figure 4. Thematic map for electrical conductivity (2006).
Figure 5. Thematic map for electrical conductivity (2007).
Figure 6. Thematic map for electrical conductivity (2008).
Figure 7. Thematic map for salinity (2006).
Figure 8. Thematic map for salinity (2007).
Figure 9. Thematic map for salinity (2008).
EC and sal contours prevailing in the delta, for 2006, 2007 and 2008, are presented in Figures 4-9, respectively. The EC and sal thematic maps give almost the same picture for all periods. The affect of seawater is reflected in a high value of EC (over 1500 µS/cm) and sal (over 0.6) of groundwater. It is obvious that groundwater quality is seriously degradated in the region. These thematic maps indicate that sea water is affected almost more than 2 km in inland.
The natural recharge of the aquifer seems not to be quantitative adequate to improve significantly the water quality.
Owing to increasing demand for water, groundwater has been pumped to over-exploitation and natural equilibrium has thus been disturbed. This has resulted in seawater mixing reaching inland aquifers.
5. Conclusions
Groundwater is an indispensable resource for water supply for the towns on Göksu Delta, for the urban and for sources has seriously impeded the development of industry, agriculture and the improvement of the people’s living standards in this internationally protected region.
The chemistry and the chemical types of groundwater indicate the sea water intrusion in Göksu delta. The chemical characters of waters in the region affected by the sea water are changed. The statistical comparation analysis supports also this event.
With the equi-concentration mapping of EC and sal the extent of the seawater pollution in the aquifer was determined. The seawater mixing reached more than 2 km inland in the north and east direction. The salinity and electrical conductivity measurements and the thematic maps indicate that the origin of groundwater salinity in Göksu Delta is due to seawater intrusion which resulted from intensive pumping.
6. Acknowledgements
The author would like to acknowledge the Scientific and Technological Research Council of Turkey, (TUBITAK) Ankara, for financing this project (105Y285).