1. Overview Saltwater Intrusion into Coastal Aquifers
Coastal aquifers serve as major sources for freshwater supply in many countries around the world, especially in the Mediterranean [
1]. The fact that coastal zones contain some of the most densely populated areas in the world makes the need for freshwater even more acute [
1]. The intensive extraction of groundwater from coastal aquifers reduces freshwater outflow to the sea and creates local water table depression, causing seawater to migrate inland and rising toward the wells [
2,
3,
4], resulting in deterioration in groundwater quality. This phenomenon, called seawater intrusion, has become one of the major constraints imposed on groundwater utilization in coastal areas.
Saltwater intrusion is one of the most widespread and important processes that degrade water quality to levels exceeding acceptable drinking and irrigation water standards, and endanger future water exploitation in coastal aquifers. Coupled with a continuing sea level rise due to global warming, coastal aquifers are even more under threat. This problem is intensified due to population growth, and the fact that about 70% of the world’s population occupies the coastal plain zones [
5,
6]. The intensity of the problem depends on the amount of the abstraction, in relation to the natural groundwater recharge, as well as on the well field location and design, the geometry, and the hydrogeological parameters of the pumped aquifer.
In recent years, there is an increasing interest in evaluating the extent of seawater intrusion in response to overexploitation and sea level rise [
7]. Seawater intrusion phenomena have been reported with different degree, in almost all coastal aquifers around the globe. In the United Sates, saltwater intrusion into coastal aquifers has been identified in the eastern Atlantic [
8,
9,
10,
11,
12,
13,
14], and the southern [
15] and western Pacific [
16,
17,
18,
19] coasts.
Reference [
20] provides an overview of saltwater intrusion in the 17 coastal states of Mexico, which is one of the most important cases around the globe. Seawater intrusion induced by groundwater development is also known in South America [
21] and Australia [
22]. In Africa, several cases of seawater intrusion into coastal aquifers have been reported [
23] with case studies in Morocco [
24,
25,
26], Tunisia [
27], Algeria [
28] and Dar es Salaam [
29,
30]. In Europe, seawater intrusion has been documented within many of the coastal aquifers, particularly along the North Sea [
31,
32] and the Mediterranean Sea and its eastern part in the archipelagos in the Aegean Sea [
33,
34,
35].
The arid and semi-arid areas are mostly chronically water-stressed. The problem of seawater intrusion is more severe in arid and semi-arid regions where the groundwater constitutes the main freshwater resource, which is mostly non-renewable. At present, developing countries of the Mediterranean Basin in North Africa and Middle East face environmental pressures induced by high population growth, rapid urbanization, and deficient water sector services reflecting on improper management of water resources [
36,
37]. The shortage of water in the Mediterranean region has been affected by the impact of climate change (increase of temperatures, variation of precipitations and high potential of evapotranspiration). Once again, the impacts have different effects in the Mediterranean region: the semi-arid and arid regions of the basin are exposed to desertification, increasing salinity of freshwater and exhaustion of water sources. Climatic change will also alter the marine environment, with an expected rise in sea level modifying several shores of the Mediterranean countries. To overcome the consequences of water scarcity and climate change, the aquifers and groundwater seem to be the solution. Many of the most important water projects in these regions focus on fossil water creating a sort of “pumping race” between the countries that share common aquifers, where overexploitation of groundwater in these regions is the major cause of seawater intrusion problems [
34,
35,
38,
39,
40,
41,
42,
43,
44,
45,
46]. Therefore, the main challenges in coastal areas in the semi-arid region are water conservation, management and planning of the water resources. This is further complicated with several complexities of the geological formations. With the semi-arid conditions, complex geological settings and over-shooting stresses, the aquifer system becomes extremely fragile and sensitive [
1]. Despite a good amount of research in this field, it is still needed to understand the behavior of such complex system precisely and apply the result in reasonably larger scales.
In the western Mediterranean, the situation of the groundwater in the Maghreb countries (Egypt, Libya, Tunisia, Algeria and Morocco) in North Africa has been marked by continuous decreases of water levels in coastal aquifers reaching alarming values. This decrease, caused by the synergistic effects of drought, flooding, changing land use, pollution from agriculture and industrialization, has intensified the problem of seawater intrusion [
1]. The Korba coastal aquifer situated in Cap-Bon, Tunisia, has been experiencing seawater intrusion since 1970 and currently the salt load in this unconfined aquifer has peak concentrations of 5–10 g/L [
27,
47]. The Algerian coastal aquifers have also not escaped overexploitation with the Mitidja aquifer suffering from seawater intrusion [
28], especially during the dry season. This aquifer system has a steady decline in water level in the order of 20–50 m per decade, which increases the rate of seawater intrusion on an annual basis [
48]. The origin of water salinity on the Annaba coast (North East Algeria) is attributed to several factors such as the geological features of the region, the climate and the salt deposits. The salinity increases steadily when approaching the sea, and indicates the influence of marine water [
49]. In Morocco, areas have been identified in which saltwater intrusion occurs (Temara-Rabat; Nador: [
24]; Saidia: [
25], however the aquifer system also contains marine deposits which contribute to the degradation of the groundwater quality. The rates of water abstraction in these areas have increased in the last 50 years, resulting in the lowering of the water table and eventually allowing seawater to intrude from coastal areas [
50].
In the Nile delta, seawater intrusion has been observed 60 km inland as a result of excessive pumping [
51]. An extensive saltwater body has developed from upper Egypt to eastern Libya in the past 50 years. The freshwater/saline water interface passes through the Qattara depression crossing the Libyan-Egyptian border and finally turning to the Southwest reaching the Tazerbo area, southeast Libya [
1]. The development of the Siwa oasis from the deep Nubian Sandstone Aquifer is close to the freshwater/saline water interface, and could cause the saline water to intrude into the freshwater aquifer (Internationally Shared Transboundary Aquifer Resources Management [
52]. The problem is further compounded since on the Libyan side large amounts of water are abstracted for urban development, causing saltwater intrusion along the Libyan coast. This overabstraction in combination with the sluggish flow of the Nubian Sandstone Aquifer causes the saline water body to encroach even further inland with considerable increases in salinity due to seawater intrusion and upconing of deep saline water [
46,
53].
Progressive seawater intrusion in the coastal aquifers of Libya has been experienced since 1930s because of its ever-increasing water demand from underground water resources. Since the 1960s, the risk of seawater intrusion is continuously threatening large coastal parts of the Jifarah Plain that forms one of the economically most significant areas in Libya, where TDS peaks up to 10 g/L are recorded [
54]. Numerous irrigated regions are located near the coast, principally in the northern part of Jifarah Plain including Tripoli region, where extensive irrigated areas have been established in the late 1970s and have evolved into advanced agricultural production zones; these activities are primarily dependent on groundwater extraction.
In Tripoli, the seawater intrusion has steadily increased from 1960 to 2007, a period during which potable water was available from the aquifer. Since 1999, a loss of 60% in well production in the upper aquifer has been observed [
55].
Because of the accelerated development of the coastal zone of Tripoli in the last decades, it is necessary to evaluate the saline water intrusion phenomena of the coastal aquifer to open different choices for the rational exploitation of the groundwater resources in this semi-arid zone avoiding the degradation of groundwater quality.
2. Introduction to the Study Area
Libya’s coastal area is one of the important cases in the arid and semi-arid regions, it is a south Mediterranean country and has a shoreline extent of about 1750 km. Groundwater is the main source for potable, industrial and irrigation water because of its semi-desert climate. Inevitably, groundwater extraction has been in excess of replenishment because of the rapid increase in agricultural and economic activities in the last 50 years. This has resulted in water level decline and deterioration in quality, including invasion of seawater along the coastal regions.
This situation has led to two significant problems linked to human activity: (1) salinization due to the formation of large piezometric drawdown cones, which have induced seawater intrusion by reversing the hydraulic gradients into aquifers; and (2) direct input of nitrate mainly from fertilizers and sewage. Agriculture is based on intensive irrigation and fertilization to improve the soils.
The Jifarah Plain in the northwest of the country, located between the Mediterranean coast and the Jabel Naffusah Mountain in the south, contains more than 60% of the country’s population and produces 50% of the total agricultural outputs. Tripoli area is a typical example showing the problems of coastal zones under high anthropogenic pressure in dryland regions. Tripoli city forms an almost rectangular area (763 km
2) between the Mediterranean Sea and the cities of Swani and Bin Gashir in the south (
Figure 1). This area extends for about 20 km along Tripoli coastal area and about 22 km inland. Topography is rising towards the south and east, a general trend in overall Jifarah Plain, which is bounded to the south and east by Jebal Nafusseh Mountains. The shortage of good quality water from surface sources has made groundwater to be very important in the study area. The scarcity of water in Tripoli is becoming more pronounced due to the increase of the population coupled with improvement of the standard of living over the last few decades, where the area of the study accommodates dense population with more than 1.5 million of inhabitants mostly concentrated in the coast.
The principal aquifer used by the population in Tripoli is the Upper Miocene-Pliocene-Quaternary aquifer system, called “first aquifer” or “upper aquifer”; intercalated thin clayey sand and marl series are dividing the aquifer into a number of horizons, all are considered as one unconfined unit [
56]. The Tripoli upper aquifer is affected by different sources of salinization, most serious is seawater intrusion [
54]. The aim of this study was to discover what processes have been responsible for variations in the chemical composition of groundwater in the upper aquifer of Tripoli and to recognize the different sources of pollution, and their relation to the intense water withdrawal.
The climate in the study area is arid to semi-arid and typically Mediterranean, with irregular annual rainfall. The average annual rainfall and evapotranspiration rates are 350 mm/year and 1520 mm/year, respectively [
54]. The estimation of groundwater exploitation (in the whole Jifarah Plain) from the main upper aquifer shows that the total amount of groundwater pumped in the Jifarah Plain for domestic, industrial and agricultural uses amounts to 1201.30 Mm
3/year [
54]. For drinking water supply and domestic wells, the overall amount pumped is 6% of the total amount of groundwater extraction. The yield of irrigation wells was estimated to be 1123 Mm
3/year, which is equal to 93% of the total amount of groundwater extraction. The industrial sector pumps only 1% of the total groundwater exploitation in the plain. Since 1996, the Great Man-Made River Project is supplying the plain with an amount of 149 million m
3/year of water. This amount has been considered in the total abstraction estimation [
54].
Geological and Hydrogeological Setting
The Jifarah Plain, including the study area, has been the subject of numerous geological studies [
57,
58,
59,
60].
Jifarah Plain is situated on the continental margin of Africa. Although the plain is thought to be underlain by Paleozoic rocks, the oldest encountered in boreholes are Triassic in age. These are continental, passing upwards into evaporites, a sequence thought to represent the progressive subsidence of the margin during major Mesozoic extension of the Tethyan Ocean. Continued subsidence through the Jurassic and Early Cretaceous led to the deposition of marine sequences.
Figure 2 shows the geological cross-section in Jifarah Plain crossing Tripoli region. The location of the cross-section is indicated in
Figure 1. The sediments of the Jifarah plain have been deposited since early Mesozoic times in a near shore lagoonal environment. The lithology of the upper aquifer varies widely and includes detrital limestone, dolomite, gravel, marl, clay, silt, sand, sandstone, gypsum/anhydrite and calcarenite. Middle Miocene clay separates the upper aquifer system in the area from the middle aquifer. The depth to the bottom of the upper aquifer varies between 30 and 200 m and depths of the wells that are utilizing this aquifer are between 10 and 180 m. Most of the wells tapping this aquifer give productivity of about 20–80 m
3/h [
56].
The geological deposits, that are playing a role in the hydrogeology of the area comprising the Upper-Miocene-Pliocene-Quaternary formations, are given in
Table 1.
The Pleistocene formations include terraces, which consist of cemented gravel and conglomerate. Al Kums Formation consists of limestone and dolomite. Qasr Al Haj Formation is mainly alluvial fans and cones consisting of clastic materials derived from the scarp. Jeffara Formation consists mainly of silt and sand, occasionally with gravel caliche bands; it covers extensive parts of the Jifarah Plain. Gergaresh Formation, which is known as Gergaresh Sandstone of Tyrrhennian age, occasionally contains silt lenses, conglomerate and sandy limestone.
The Holocene deposits include recent wadi deposits; these deposits consist of loose gravels and loam. Beach sands are represented by a narrow strip at the coast and are made up of shell fragments with a small ratio of silica sands. Eolian deposits are represented by sand dunes and sheets covering large parts of the coastal strip (coastal dunes). These coastal dunes consist of shell fragments with small amounts of silica sands. It is worth mentioning that the eolian material composing coastal dunes contains a large amount of grains of gypsum. In some places, it is composed of nearly pure gypsum (98%) especially in the immediate vicinity of the sebkhas, with a silty gypsum filling [
58]. Sebkha sediments are mainly gypsum deposits and are observed along the coastal area of the plain. They occupy the relatively low topographic areas and are separated from the sea by sea cliffs. Some of the sebkhas have occasional incursions of the sea and others may have subsurface connection with the seawater.
4. Results and Discussion
4.1. Water Level and Hydrodynamics
Figure 4 shows a piezometric map of Tripoli based on field measurements. In the coastal area, 64 shallow and deep wells located at different distances from the Mediterranean Sea were selected for water level measurement (see
Figure 1). Water depth was measured from the ground surface using water level meter and was converted into water level by subtracting from ground elevation.
The overall direction of groundwater flow in Jifarah Plain in general, is from the south from Jebal Naffusah Mountains to the coast. The important storage withdrawal by overexploitation from the upper aquifer in Tripoli is causing continuous drawdown of the water level, reducing the outflow rate to the sea, and the progressive degradation of the chemical quality of water. Groundwater level is mostly low, especially near the coast, where zero and negative heads are recorded for the majority of wells. The piezometric level in depression cones at the location of the public water supply well field of As Swani (123 wells) has dropped from 25 to 33 m below sea level (
Figure 4), which testifies the inversion of the hydraulic gradient and the intrusion of seawater.
From
Figure 4, it can be concluded that the general groundwater flow, from south to north following the topography, has altered, where, generally along the coast, flow is toward the reduced heads in the stressed areas around the depression cones.
To deal with the shortage of water in most coastal cities including Tripoli region, the Libyan government established the Great Man-Made River Project (GMMR) to transport millions of cubic meters of water a day from desert well fields to the coastal cities, where over 80% of the population lives. Since 1996, The GMMR is supplying Tripoli city with an amount of 149 million m3/year. The total planned supply by the project is 900 million m3/year for the whole Jifarrah Plain. The implementation of the project was interrupted since 2011, due to political situation and Tripoli is the only supplied part of Jifarrah coast.
Since the start of the supply in 1996, the wellfield of As Swani (123 public water supply wells) was stopped. The pumping from As Swani wellfield is the main cause of the depression cone in Tripoli; the depression cone location at the center of As Swani wellfield is the most affected part of the region. Now the GMMR is the main supplier for domestic use in the city and being used also in many farms. The drawdown value around the depression cone is being reduced from 80 m below sea level in 1996 to 33 m in 2008.
4.2. Major Hydrochemical Parameters
Major anions and cations were analyzed and pH, Eh, electrical conductivity (EC) as well as temperature were assessed on all samples. The results show that: temperature ranges 18–25 °C, pH range is 7.17–9.94, Eh range is 139–240 mV, EC range is 369–10,600 µS/cm (25 °C), TDS range is 340–6529 mg/L and chloride concentration ranges 39–3155 mg/L.
Table 6 shows analytical results for physico-chemical parameters of groundwater for selected representative samples in Tripoli. The high concentration of major ions such as Cl, Na, Ca, K, Mg, and a high EC indicate the presence of seawater in an aquifer [
68,
69].
Out of 64 samples analyzed, 38% have NO
3− higher than the highest desirable level of 45 mg/L [
70]. For sulfate 40% have SO
42− higher than the highest desirable level of 200 mg/L according to WHO (2008), with a maximum of 835 mg/L recorded southwards. Out of 64 analyzed samples, 58% exceed the recommended Cl
− value for standard drinking water (250 mg/L) and 26% have Cl
− greater than the highest admissible level of 600 mg/L [
70].
Levels of Cl
− and EC are the simplest indicators of seawater intrusion or salinization [
71,
72]. EC is positively correlated with the concentration of ions, mainly Cl
− concentration.
Figure 5 shows three zones on a plot of Cl
− vs. EC: freshwater zone, mixing zone and strong mixing (intrusion). It shows that groundwater samples with Cl
− exceeding 200 mg/L and EC exceeding ~1000 µS/cm are most likely influenced by seawater intrusion. Groundwater samples that are characterized by EC between 1000 and 5000 µS/cm represent a mixing between freshwater and saltwater. Samples with EC of more than 10,000 µS/cm represent strong seawater influence.
The spatial distribution of EC and Cl
− from analyzed groundwater samples across Tripoli city is presented in
Figure 6a,b. In general, the EC, which is tightly linked to TDS, is a measure of salinity, and therefore is generally closely related to the Cl
− content. Both EC and Cl
− show the same general decrease from the Mediterranean shoreline towards the south. Frequent local increases in both variables are observed at the depression cones as a result of the high pumping rate.
The high Cl− concentration is due to mixing with seawater. These high concentrations of chloride occur in most wells within a few kilometers of the coast, and are related to active seawater intrusion. However, high Cl− is also found far inland at and nearby the depression cones, whereby deeper saline groundwater is affected by upconing due to groundwater exploitation. The concentration of Cl− decreases gradually towards the south. However, in many farther inland wells, it is still the dominant anion. The higher concentration of Cl− than the 250 mg/L value for standard drinking water at the south of the region can be linked to the synsedimentary marine influence of the groundwater.
Sulfate concentration in Tripoli ranges from 29 to 835 mg/L.
Figure 6c shows the spatial distribution of SO
42− concentration in the study area. The main source for increasing SO
42− is mixing with seawater, which can add significant amounts of sulfate to freshwaters. High SO
42− is mostly linked to the high Cl
− concentration in the upper aquifer, both in the seawater intrusion zones and in the depression cones, where deep saline water upconing occurs.
Besides, more than 500 mg/L SO42− is observed towards the west of region, with much higher SO42−/Cl− compared to seawater. The main source of SO42− in this area is the dissolution of gypsum/anhydrite from the superficial sebkha deposits in those areas, as these wells are located near the vicinity of sebkhas. In this zone, lower Cl− is recorded, excluding seawater intrusion as the source.
4.3. Water Types and Piper Diagram
Classification of hydrochemical facies for groundwaters according to the Piper diagram is represented in
Figure 7.
In the Piper diagram, almost all water samples are plotted above the general seawater–freshwater mixing line [
64], comprising freshwater sample TJ17 and Mediterranean Seawater. Although various hydrochemical facies were observed (NaCl, CaCl, MgCl, CaMgHCO
3, NaHCO
3, NaSO
4 MgSO
4 and CaSO
4), CaCl and NaCl types are dominant. Large proportion of the groundwater shows NaCl type, which generally indicates a strong seawater influence [
50] or upconing of deep saltwater, while CaCl water type indicates salinization and cation exchange reaction [
73]. The region of the CaCl type water may be a leading edge of the seawater plume [
64,
74,
75]. Furthermore, sources of CaSO
4 water type are the dissolution of the scattered sebkha deposits.
4.4. Hydrochemical Profile and Water Classification According to Stuyfzand
Salinization is induced as the new saline end member is introduced into the freshwater aquifer. The main chemical reaction is cation exchange, resulting in deficit of Na
+ and surplus of Ca
2+:
where X represents the natural exchanger in the reaction. During cation exchange, the dominant Na
+ ions are adsorbed and Ca
2+ ions released, so that the resulting water moves from NaCl to CaCl water type, which is typical for salinization [
5]. The salinization process can be schematized as follows [
75] (water types according to classification of Stuyfzand, [
65]):
The chloride ion concentration is taken as a reference parameter [
64]. Therefore, as saltwater intrudes into coastal freshwater aquifers, the Na/Cl ratio decreases and the Ca/Cl ratio increases.
Upon the inflow of freshwater, a reverse process takes place:
Flushing of the saline aquifer by freshwater will thus result in uptake of Ca
2+ by the exchanger with concomitant release of Na
+. This is reflected in the increase of the Na/Cl ratio, and formation of the NaHCO
3 water type, which is typical for freshening. The anion HCO
3− is not affected because natural sediments behave as cation exchanger at the usual near-natural pH of groundwater [
64]. The freshening process can be schematized as follows [
73]:
The major hydrogeochemical processes occurring in the upper aquifer are: mixing with seawater end member, cation exchange during salinization, dissolution of gypsum from superficial sebkha sediments, carbonate dissolution and agricultural pollution.
The hydrogeochemical profile in Tripoli (Janzur) (
Figure 8) is selected as an example showing the distribution of water types along the flow path. Janzur profile is about 14 km long perpendicular to the sea (see
Figure 2). Seventeen wells were visited in this profile region. Their total depth is shallow (between 10 and 50 m) close to the coast, and reaches 120 m southward. They are pumping for irrigation in the private farms and for domestic use. The groundwater table in the profile is between 0 and 4 m a.s.l in the north, while it is down to −24 m a.s.l in the south at the depression cone. Groundwater is flowing from the south and north to the locally reduced heads, located at about 10 km from the coast.
The spatial distribution of water sub-types according to Stuyfzand in Tripoli is presented in
Figure 9 together with the pie plot of TDS distribution for selected representative samples in Tripoli region. The water classification scheme of Stuyfzand [
65] has the advantage that in brackish or saline groundwater, one can still identify many different water types even though the major anions and cations are the same and this may help to recognize processes such as upconing of more saline water in the aquifer [
67].
Figure 8 and
Figure 9 show water type is CaSO
4 in the west, CaHCO
3, NaHCO
3, CaMix and MgMix towards the south and NaCl, MgCl and CaCl in the north and at the depression cones. Close to the shoreline and mainly in the east in Tajura the water is NaCl type, due to the strong effect of seawater intrusion. CaCl results from cation exchange, due to mixing with seawater. Towards the south, CaMgMix(ClHCO
3) evolving further inland to CaMgMix(HCO
3Cl), indicates the location of the transition zone, where the groundwater changes from CaMgMix enriched with Cl
− ion to CaMix with HCO
3− as the dominant anion. The CaSO
4 water type observed in western Tripoli, up to about 14 km inland, shows the dissolution of the evaporitic rocks from the sebkha deposits.
The Mg
2+ content found in several wells, is mainly resulting from the freshwater end member coming from the recharge area, where Mg
2+-containing carbonate is dissolved [
54]. Thus, in this case, thee positive cation exchange code in the classification does not indicate freshening, as Mg
2+ is not supplied by the marine end member [
69]. At the depression cones and downstream, cation exchange equilibrium (cation exchange code “0”) exists for several wells, which in this case indicates the onset of the salinization process (the positive value of (Na
+ + K
+ + Mg
2+)
corrected is decreasing as the marine cations are adsorbed during salinization).
The increase of salinity in Tripoli is accompanied by an increase in NO
3− concentrations.
Figure 8 also shows the spatial distribution of NO
3− along Janzur profile. The average nitrate concentration of groundwater in the aquifer is about 38 mg/L, but contents as high as about 118 mg/L occur upstream in the south of the region. Irrigation with nitrogen fertilizers and domestic sewage and movement of contaminants in areas of high hydraulic gradients within the drawdown cones probably are responsible for localized peaks of the nitrate concentration for many wells in the region.
4.5. Ionic Ratio
Conservative seawater–freshwater mixing is expected to show a linear increase in Na
+ and Cl
− [
76], which is reflected by the high correlation coefficient (
r = 0.94) between both variables (
Figure 10a). Effects of seawater encroachment have been evaluated by studying the Na/Cl ionic ratio. Lower ratios of Na/Cl than seawater values (0.88) indicate seawater encroachment.
Figure 10b shows molar ratios of Cl
− versus Na/Cl concentrations. The Na/Cl ratios for the analyzed samples range from 0.23 to 1.68. Most of groundwater samples were less than or slightly higher than the Mediterranean seawater ratio (0.88). The lowered values with respect to the Mediterranean seawater ratio are resulting from cation exchange occurring when seawater intrudes freshwater aquifers, resulting in the deficit of Na
+ and surplus of Ca
2+. High ratio for several samples towards the recharge area (e.g., TJ17) indicates flushing the aquifer by freshwater from the south.
The ratio of SO4/Cl (meq/L/meq/L) for the Mediterranean seawater and the fresh recharge water from the study area are 0.54 and 0.25, respectively. The ratio of SO4/Cl for the analyzed samples ranges from 0.04 to 0.77, which indicates mixing between seawater and freshwater. Significantly higher values than Mediterranean seawater ratio (0.54) indicate dissolution of gypsum from Sebkha deposits.
The Na/K ratio is ranging in the area between 9.64 and 102.66, the largest values are observed in the area affected by the intrusion of seawater.
4.6. Saturation Indices
The saturation indices (SI) for calcite, dolomite, halite, aragonite, gypsum and anhydrite were calculated to verify precipitation and dissolution of these minerals. The selected minerals were based on the major ions in groundwater from the study area.
Figure 11 is a synthetic diagram showing SI values for calcite, dolomite, gypsum, anhydrite, halite and aragonite. The sample numbers are sorted according to their location from west to east and from north to south with each profile, but they are not plotted at spatial distances.
Out of 64 groundwater samples, 80% of groundwaters seem to be supersaturated (SI > 0) with respect to calcite (CaCO3), whereas 10% are undersaturated with respect to calcite (SI < 0) and 10% are at equilibrium (SI = 0). Dolomite (MgCa(CO3)2) seems to be oversaturated in 85% of groundwater samples analyzed, 5% are undersaturated and 10% are at equilibrium. Ninety-eight percent of groundwater samples in the study area are undersaturated with respect to gypsum (CaSO4·2H2O) and anhydrite (CaSO4).
In general, most of the analyzed samples have saturation indices close to saturation with respect to calcite (SI mostly 0–1) and dolomite (SI mostly 0–3). This slight supersaturation with respect to calcite and dolomite supersaturation rather points to groundwater in equilibrium with those minerals. During sampling, most often dissolved CO2 gas escapes, slightly raising pH and thus shifting carbonate equilibrium (more CO32−), such that SI > 0 is obtained, whereas in water in the aquifer SI with respect to calcite is close to zero. Thus, the water in the aquifer is not really oversaturated.
The majority of samples in the study area are undersaturated with respect to gypsum and anhydrite. Gypsum comes close to saturation (SI > −0.50) in several wells. The dissolution of gypsum from the superficial sebkha deposits for many wells in the coastal area raises the SI.
4.7. Deviation from Conservative Mixture of End Member Fraction
Figure 12 shows the ionic deltas calculated for Na
+, Ca
2+, Mg
2+, K
+, HCO
3− and SO
42− for all analyzed samples. The first thing to note is that the process of cation exchange due to salinization is very evident. For example, in
Figure 12, m
Na+,reaction (ΔNa
+) is plotted in the secondary axis; the ΔNa
+ is usually positive for freshwater, but a large number of samples have negative values particularly in the highly saline water, down to −23 mmol/L. The most logical explanation for this deficit of Na
+ is that a reverse cation exchange reaction is taking place during the salinization process, which releases Ca
2+ to the solution and captures Na
+. The reverse relationship between the two ions (Na
+ and Ca
2+) is noticed particularly in the highly saline groundwater, where samples with large negative values of ΔNa
+ generally show strong positive ΔCa
2+. Furthermore, also potassium shows negative (or very low positive) deltas characteristic for marine cations as a result of the salinization process.
ΔMg
2+ is mostly positive, due to more Mg
2+ added by dissolution of Mg
2+-rich carbonate than adsorbed at the clay exchange complex during salinization. Only very few samples show a deficit of Mg
2+.
Figure 12 also shows that the ionic delta of ∆HCO
3− is positive for most water samples. This is due to the dissolution of carbonate minerals in the aquifer deposits. In general, most samples show positive ∆SO
42−. The gypsum dissolution from sebkha deposits increases ∆SO
42− to high positive values for several wells.
5. Conclusions
Hydrochemistry of the coastal aquifers of arid and semi-arid regions is very complex. Extensive groundwater extraction in Tripoli region, mainly for agricultural development, has caused substantial seawater encroachment and upconing of the deep saline water into Tripoli shallow aquifer, with Cl−, SO42− and NO3− as the major pollutants. In this study, this potential problem is investigated. The dominant water types in the study area are NaCl, CaCl and CaMgMix(ClHCO3) except for several wells towards the recharge area, where CaHCO3-type prevails, and wells located near the superficial sebkha deposits, where the CaSO4 water type evolves.
Seawater intrusion is accompanied by chemical reactions, which modify the hydrochemistry of the coastal aquifer. The most remarkable reaction is that of the inverse cation exchange, characteristic of the changes of the theoretical mixture of seawater–freshwater, which is carried out between clays and the aquifer water. This exchange consists in the release of Ca2+ and the adsorption of Na+.
Great part of the observed high concentration of sulfate in Tripoli is coming from the effect of seawater intrusion. Furthermore, the scattered sebkha deposits in the north, containing large amounts of gypsum, produce high SO42− waters.
Another serious problem in the study area is the increased level of nitrate concentrations. It has been found that in large number of samples, nitrate contents range between 45 and 160 mg/L. The increased content originates from leaching of nitrates from the applied nitrogen fertilizers and from sewage. It is recommended that risk assessment of nitrate pollution is useful for a better management of groundwater resources, aiming at preventing soil salinization and minimizing nitrate pollution in groundwater.
The hydrochemical interpretation also indicates that the dissolution of calcite, dolomite and/or Mg2+ bearing calcite is an important process in most of the groundwaters. The saturation index shows mostly a slight tendency to precipitation of calcite and dolomite in the aquifer system, but this can be ascribed to lowering of CO2 pressure at sampling, while in the aquifer, there is equilibrium with respect to these minerals.
Although the Great Man Made River Project is supplying Tripoli with an amount of 149 million m
3/year of water since 1996, used mainly for domestic purposes, results show high degradation level of groundwater quality and most of water samples do not compare favorably with WHO standards [
70]; many samples exceed the maximum admissible concentrations, highlighting the degradation of groundwater quality. The recovery of groundwater quality is usually a very slow process as seawater intrusion is the result of a long-term negative mass balance in the aquifer. A balance between pumping demand and quality requirements is necessary. This balance is hard to maintain when the final goal is to reverse the qualitative status of the already contaminated aquifer. To protect the groundwater resource in the long-term, on which the future Tripoli residents depend, appropriate management against overexploitation from agricultural activity to control salinity is compelling, and especially urgent in the coastal fringe, where seawater intrusion is threatening. Artificial recharge of coastal aquifers, which are especially overexploited, may offer an efficient means of combating seawater intrusion and thus of preventing an inevitable degradation of the water quality which might prove irreversible.