The onset of the Messinian salinity crisis in the deep Eastern Mediterranean basin

Astronomical tuning of the Messinian pre‐salt succession in the Levant Basin allows for the first time the reconstruction of a detailed chronology of the Messinian salinity crisis (MSC) events in deep setting and their correlation with marginal records that supports the CIESM ( 2008 ) 3‐stage model. Our main conclusions are (1) MSC events were synchronous across marginal and deep basins, (2) MSC onset in deep basins occurred at 5.97 Ma, (3) only foraminifera‐barren, evaporite‐free shales accumulated in deep settings between 5.97 and 5.60 Ma, (4) deep evaporites (anhydrite and halite) deposition started later, at 5.60 Ma and (5) new and published ⁸⁷Sr/⁸⁶Sr data indicate that during all stages, evaporites precipitated from the same water body in all the Mediterranean sub‐basins. The wide synchrony of events and ⁸⁷Sr/⁸⁶Sr homogeneity implies inter‐sub‐basin connection during the whole MSC and is not compatible with large sea‐level fall and desiccation of the Mediterranean.     

The Mediterranean Messinian salinity crisis: an Apennine foredeep perspective

Owing to its expanded stratigraphic sections, the Apennine thrust belt offers the opportunity to better understand the evaporitic and post-evaporitic Messinian events. A physical stratigraphic framework of Messinian deposits, based on facies analysis and basin-wide correlation of key surfaces and sedimentary cycles, is presented. It is shown that the Messinian Apennine foredeep had marginal basins with shallow-water primary evaporites and deeper basins where resedimented evaporites accumulated under relatively deep-water conditions. Like many other Mediterranean examples, primary shallow-water evaporites of Apenninic marginal basins show evidence for subaerial exposure and erosion. However, the development of such an erosional surface does not correspond to the deposition of primary evaporites in the deepest part of the basin(s); here, the unconformity can be traced towards the base of resedimented evaporites or to a level within them, implying that the deeper basins of the Apennine foredeep never underwent desiccation during the Messinian salinity crisis, but rather received the eroded marginal evaporites. This fact, usually overlooked, raises important questions about the deep desiccation model of the Mediterranean.     

The Messinian Sicilian stratigraphy revisited: new insights for the Messinian salinity crisis

Controversies around the Messinian salinity crisis (MSC) are because of the difficulties in establishing genetic and stratigraphic relationships between its deep and shallow‐water record. Actually, the Sicilian foreland basin shows both shallow and deep‐water Messinian records, thus offering the chance to reconstruct comprehensive MSC scenarios. The Lower Gypsum of Sicily comprises primary and resedimented evaporites separated in space and time by the intra‐Messinian unconformity. A composite unit including halite, resedimented gypsum and Calcare di Base accumulated between 5.6 and 5.55 Ma in the main depocentres; it records the acme of the Messinian Salinity Crisis during a tectonic phase coupled with sea‐level falls at glacials TG14‐TG12. These deposits fully post‐date primary gypsum, which precipitated in shallow‐water wedge‐top and foreland ramp basins between 5.96 and 5.6 Ma. This new stratigraphic framework results in a three‐stage MSC scenario characterized by different primary evaporite associations: selenite in the first and third stages, carbonate, halite and potash salt in the second one associated with hybrid resedimented evaporites.     

A new facies model for the Upper Gypsum of Sicily (Italy): chronological and palaeoenvironmental constraints for the Messinian salinity crisis in the Mediterranean

The Upper Gypsum unit of the Caltanissetta Basin (Sicily) records the last phase of the Messinian salinity crisis comprising the so‐called ‘Lago Mare’ event. A new facies analysis study recognizes nine to ten depositional cycles consisting of seven rhythmically interbedded primary gypsum bodies, and two to three sandstone bodies separated by marly terrigenous horizons showing laterally persistent vertical organization. A basal thin gypsum bed is overlain by a cluster of five thicker gypsum bodies. A marly interval containing two distinct sandstone horizons separates this cluster from the overlying uppermost (seventh) gypsum body. The terrigenous Arenazzolo Formation, in turn followed by the lower Pliocene Trubi Formation, is considered here to form the uppermost part of the Upper Gypsum unit. The rhythmic alternation in the sandy marls and gypsum/sandstone bodies records the response of sediments from shelfal to deltaic systems to precession‐driven arid‐wet climate fluctuations causing cyclical changes of both base‐level and water concentration. During wet climate phases (at insolation maxima) marl and sandstone were deposited in a hypohaline environment as suggested by: (i) the typical Lago Mare faunal assemblage and (ii) the negative δ¹⁸O values. During arid phases (at insolation minima) the reduced meteoric supply, recorded by higher δ¹⁸O values in the carbonate, caused the development of a negative hydrological budget leading to evaporite precipitation. At a basinal scale the Upper Gypsum unit unconformably overlies a mainly clastic evaporite unit containing carbonate breccia (the so‐called ‘Calcare di Base’) and/or clastic gypsum. Towards the basin centres, where the basal contact becomes conformable, a primary gypsum cumulate horizon is present. This layer is interpreted as a possible lateral equivalent of the Halite unit present only in the deepest depocentres. Based on astronomical calibration of the depositional cycles, the Upper Gypsum unit, including the Arenazzolo Formation, spans the interval between 5·33 and 5·53 Ma. This new age calibration allows the deposition of the Halite unit to be dated between 5·6 Ma (top of the Lower Evaporites) and 5·55 Ma (base of the Upper Evaporites) corresponding to isotopic stages TG12 and/or TG14.     

Quantitative evaluation of the depth of the western Mediterranean before, during and after the Late Miocene salinity crisis

Quantitative geophysical calculations which take into consideration the isostatic loading of sediment overburden, the overlying water cover, and the thermal cooling history of the continental edge, and its adjacent oceanic lithosphere, demonstrate the foundering of the margins of the western Mediterranean had already commenced in the Aquitanian stage of the early Miocene. The calculations are based on magnitudes and rates of sediment accumulation observed along a profile of three commercial boreholes into the subsurface of the continental shelf of southern France. By the time of the late Miocene (Messinian) salinity crisis, the depth of the seafloor within the Balearic basin exceeded 2.5 km. Sea-level fluctuations induced by evaporitic draw-down permitted the exposure of large tracts of the former submerged continental margins to subaerial processes. The measured magnitude of sediment removal by erosion and channel incision near the outer shelf of the modern Gulf of Lion surpasses 1 km. The subsidence history of this shelf platform south of France provides new evidence that the continental lithosphere behaves as if it is rigidly coupled to its oceanic counterpart commencing with the initial phase of the pull-apart. No major vertical fault displacements have subsequently offset their overlying crustal layers. The sedimentary shaping and construction of the margin seaward of the Rhône delta resulted in a pronounced shelf edge migration and slope progradation during the pre-salinity crisis Miocene. It has taken 5 million years of predominant upbuilding to establish a new equilibrium profile similar in cross-section to the precrisis depositional surfaces created by outbuilding.     

The Messinian Guadalhorce corridor: the last northern, Atlantic–Mediterranean gateway

Messinian marine deposits of the Guadalhorce River valley in southern Spain record evidence of the last northern gateway that existed between the Mediterranean and the Atlantic. They comprise sandstones and conglomerates with unidirectional cross-bed sets up to nearly 1 km long in their down-sedimentary-dip direction. These cross-bed sets relate to extremely fast (1.0–1.5 m s⁻¹) bottom currents flowing from the Mediterranean into the Atlantic. The Guadalhorce gateway (which had a maximum width of 5 km and a maximum water depth of 120 m) was an important element controlling the Messinian pre-evaporitic oceanic circulation in the Mediterranean Sea, as it acted as a major outflow channel. Its closure limited the exchange of water between the Atlantic and the Mediterranean to the Rifian corridors of Morocco, inducing water-mass restriction and stratification in the western Mediterranean immediately prior to the `Messinian Salinity Crisis'.     

Management scenarios for the Jordan River salinity crisis

Recent geochemical and hydrological findings show that the water quality of the base flow of the Lower Jordan River, between the Sea of Galilee and the Dead Sea, is dependent upon the ratio between surface water flow and groundwater discharge. Using water quality data, mass-balance calculations, and actual flow-rate measurements, possible management scenarios for the Lower Jordan River and their potential affects on its salinity are investigated. The predicted scenarios reveal that implementation of some elements of the Israel–Jordan peace treaty will have negative effects on the Jordan River water salinity. It is predicted that removal of sewage effluents dumped into the river (∼13 MCM/a) will significantly reduce the river water’s flow and increase the relative proportion of the saline groundwater flux into the river. Under this scenario, the Cl content of the river at its southern point (Abdalla Bridge) will rise to almost 7000 mg/L during the summer. In contrast, removal of all the saline water (16.5 MCM/a) that is artificially discharged into the Lower Jordan River will significantly reduce its Cl concentration, to levels of 650–2600 and 3000–3500 mg/L in the northern and southern areas of the Lower Jordan River, respectively. However, because the removal of either the sewage effluents or the saline water will decrease the river’s discharge to a level that could potentially cause river desiccation during the summer months, other water sources must be allocated to preserve in-stream flow needs and hence the river’s ecosystem.     

Correlations and sequence stratigraphic model for Messinian carbonate platforms of the western and central Mediterranean

Distance correlations of Late Tortonian–Messinian littoral carbonate complexes are proposed from the study of eight platforms in the western and central Mediterranean. Correlations are based on the identification of two major biological sedimentary cycles and of two index surfaces. Surface A is a maximum flooding surface during cycle 1 at around 6.7 Ma. Surface B is a regional marine planation surface at around 5.95 Ma, at the base of cycle 2 (Terminal Carbonate Complex). A general sedimentary model is proposed for the 7–5.6-Ma time-span. The boundary between cycles 1 and 2 is coincident with the onset of the Messinian Salinity Crisis, and appears to be related to major environmental–paleo-oceanographic changes in the Mediterranean, rather than to a major sea-level drop or to climatic change.     

Messinian events in the Sorbas Basin in southeastern Spain and their implications in the recent history of the Mediterranean

The Messinian (Late Miocene) marine stratigraphic record of the Sorbas Basin (S.E. Spain) is well preserved and can be considered as being representative of the entire western Mediterranean. It exhibits a series of features relating to: (1) the composition, characteristics and evolution of coral reefs; (2) changes between temperate and subtropical climates; and (3) the extensive development of microbial carbonates (stromatolites and thrombolites) at the end of the Messinian. Each of these features has global significance.

Porites, which is the major and almost only coral component in reefs, is heavily encrusted with stromatolites. These reefs grew at the edge of the subtropical belt and were totally eliminated at the end of the Messinian because of global cooling.

Lowermost-Messinian carbonate sediments in the Sorbas Basin reflect a temperate climate, whereas those immediately above, which contain bioherms and coastal reefs, are subtropical. The shift from temperate to subtropical conditions during the early Messinian was accompanied by an important change in water circulation within the western Mediterranean. Temperate times were marked by cold surface Atlantic waters entering the Mediterranean, whereas subtropical times coincided with warm surface waters entering the western Mediterranean from the east. The subtropical waters were thermally stratified, which favoured the deposition of euxinic marls and diatomites at the centre of the basin. The upwelling of nutrient-rich water promoted stromatolite development within reefs and Halimeda growth on adjacent slopes.

Lastly, microbial carbonates (stromatolites and thrombolites) attained giant dimensions during the late Messinian, which can be regarded as a measure of their success in occupying a variety of ecological niches. This abundance of available habitats is believed to have resulted from the Messinian “salinity crisis”, which was followed by a re-colonization of the western Mediterranean. In this context stromatolite proliferation was due to opportunism of microbial communities in colonizing the new environments, rather than to a complete absence of other competitive biota. We do not believe that hypersaline conditions were a causal factor in stromatolite development because of the normal-marine biota associated with them.     

The end of the Messinian Salinity Crisis in the western Mediterranean: Insights from the carbonate platforms of south-eastern Spain

How the Messinian Salinity Crisis (MSC) ended is still a matter of intense debate. The Terminal Carbonate Complex (TCC) is a late Messinian carbonate platform system that recorded western Mediterranean hydrological changes from the final stages of evaporite deposition till the advent of Lago-Mare fresh- to brackish water conditions at the very end of Messinian times. A multidisciplinary study has been carried out in three localities in south-eastern Spain to reconstruct the history of TCC platforms and elucidate their significance in the MSC. Overall, this study provides evidence that the TCC formed following a regional 4th order water level rise and fall concomitant with an opening-restriction trend. It can be subdivided into four 5th order depositional sequences (DS1 to DS4) recording two phases: (1) from DS1 to DS3, a tide-dominated ooidic to oobioclastic system with stenohaline faunas developed as a result of a 70 m water level rise. During this period, the TCC developed in a shallow sea with close to normal marine salinity; (2) in depositional sequence 4, a microbialite-dominated platform system developed. This is indicative of a significant environmental change and is attributed to a 30 to 40 m water level fall in the basins under study. These restricted conditions were coeval with intense evaporite deformation and brine recycling. The syn-sedimentary deformation of evaporites had a major impact on platform architecture and carbonate production, affecting the Messinian series throughout south-eastern Spain at the end of the TCC history. At that time, the TCC developed in a lake with fluctuating, brackish- to hypersaline water. These findings suggest a temporary restoration of marine conditions in the western Mediterranean marginal basins due to Atlantic water influxes prompted by a global sea level rise around 5.6 Ma. Whether marine conditions extended to the entire western Mediterranean still needs to be investigated.     

The effect of the Messinian Deep Stage on karst development around the Mediterranean Sea. Examples from Southern France

It is difficult to explain the position and behaviour of the main karst springs of southern France without calling on a drop in the water table below those encountered at the lowest levels of Pleistocene glacio-eustatic fluctuations. The principal karst features around the Mediterranean are probably inherited from the Messinian period (“Salinity crisis”) when sea level dropped dramatically due to the closing of the Straight of Gibraltar and desiccation of the Mediterranean Sea. Important deep karst systems were formed because the regional ground water dropped and the main valleys were entrenched as canyons. Sea level rise during the Pliocene caused sedimentation in the Messinian canyons and water, under a low hydraulic head, entered the upper cave levels.

The powerful submarine spring of Port-Miou is located south of Marseille in a drowned canyon of the Calanques massif. The main water flow comes from a vertical shaft that extends to a depth of more than 147 m bsl. The close shelf margin comprises a submarine karst plateau cut by a deep canyon whose bottom reaches 1,000 m bsl. The canyon ends upstream in a pocket valley without relation to any important continental valley. This canyon was probably excavated by the underground paleoriver of Port-Miou during the Messinian Salinity Crisis. Currently, seawater mixes with karst water at depth. The crisis also affected inland karst aquifers. The famous spring of Fontaine de Vaucluse was explored by a ROV (remote observation vehicle) to a depth of 308 m, 224 m below current sea level. Flutes observed on the wall of the shaft indicate the spring was formerly an air-filled shaft connected to a deep underground river flowing towards a deep valley. Outcroppings and seismic data confirm the presence of deep paleo-valleys filled with Pliocene sediments in the current Rhône and Durance valleys. In the Ardèche, several vauclusian springs may also be related to the Messinian Rhône canyon, located at about 200 m below present sea level. A Pliocene base level rise resulted in horizontal dry cave levels. In the hinterland of Gulf of Lion, the Cévennes karst margin was drained toward the hydrologic window opened by the Messinian erosional surface on the continental shelf.     

Messinian events in the Black Sea

Past hydrological interactions between the Mediterranean Sea and Black Sea are poorly resolved due to complications in establishing a high‐resolution time frame for the Black Sea. We present a new greigite‐based magnetostratigraphic age model for the Mio‐Pliocene deposits of DSDP Hole 380/380A, drilled in the southwestern Black Sea. This age model is complemented by ⁴⁰Ar/³⁹Ar dating of a volcanic ash layer, allowing a direct correlation of Black Sea deposits to the Messinian salinity crisis (MSC) interval of the Mediterranean Sea. Proxy records divide these DSDP deposits into four intervals: (i) Pre‐MSC marine conditions (6.1–6.0 Ma); (ii) highstand, fresh to brackish water conditions (~6.0–5.6 Ma); (iii) lowstand, fresh‐water environment (5.6–5.4 Ma) and (iv) highstand, fresh‐water conditions (5.4–post 5.0 Ma). Our results indicate the Black Sea was a major fresh‐water source during gypsum precipitation in the Mediterranean Sea. The introduction of Lago Mare fauna during the final stage of the MSC coincides with a sea‐level rise in the Black Sea. Across the Mio‐Pliocene boundary, sea‐level and salinity in the Black Sea did not change significantly.     

Late Miocene sedimentary architecture of the Ebro Continental Margin (Western Mediterranean): implications to the Messinian Salinity Crisis

The Messinian Salinity Crisis (MSC) resulted from a significant multi-phase drop and subsequent reflooding of the Mediterranean Sea from 5.96 to 5.33 Ma. Well-developed drainage networks, characterized by step-like profiles and abrasion platforms, are associated to this event. The Ebro Continental Margin (Western Mediterranean) presents an additional complexity since the capture of the drainage of the adjacent subaerial Ebro Basin took place sometime prior to the Messinian stage. Using 3D seismic reflection data, this work provides new insights into the origin of the step-like profile of the Messinian erosional surface (MES) and timing of the capture of the subaerial Ebro Basin. The results obtained indicate a sedimentary-active continental slope and delta progradation during Middle-Late Miocene, in a normal regressive context associated to a pre-Messinian proto-Ebro River. The mature development attained by the Messinian Ebro River network during the MSC corroborates that the capture of the Ebro Basin occurred prior to the MSC. The configuration of the clinoforms below the MES suggests that deltaic sediments of the Messinian Paleo-Ebro River deposited during the Tortonian and initial Messinian sea-level drawdown. The MES formed at the top of the Tortonian Highstand, where a fluvial network was deeply carved, and in the topset region of the Messinian Falling Stage Systems Tract, where minor erosion occurred. Fluvial deposits are outstandingly preserved on the main valleys of the MES. Therefore, the step-like profile of the MES was not created during Zanclean inundation, but during the latest stages of the main Messinian sea-level fall and lowstand.     

Messinian deposits and erosion in northern Tunisia: inferences on Strait of Sicily during the Messinian Salinity Crisis

Outcrops, offshore wells, electric logs and seismic profiles from northern Tunisia provide an opportunity to decipher the Messinian Salinity Crisis in the Strait of Sicily. Messinian deposits (including gypsum beds) near the Tellian Range reveal two successive subaerial erosional surfaces overlain by breccias and marine Zanclean clays, respectively. In the Gulf of Tunis, Messinian thick evaporites (mostly halite) are strongly eroded by a fluvial canyon infilled with Zanclean clays. The first erosional phase is referred to the intra-Messinian tectonic phase and is analogous to that found in Sicily. The second phase corresponds to the Messinian Erosional Surface that postdates the marginal evaporites, to which the entire Sicilian evaporitic series must refer. The Western and Eastern Mediterranean basins were separated during deposition of the central evaporites.     

Messinian palaeoclimate and palaeo-environment in the western Mediterranean realm: insights from the geochemistry of continental deposits of NW Sardinia (Italy)

To evaluate the palaeo-environmental parameters of a portion of the Sardinia–Corsica microplate during the Messinian drop in sea level, we examined the chemistry and mineralogy of upper Tortonian–lower Messinian (late Miocene) clayey continental deposits from NW Sardinia. Differences exist between the uppermost part of the succession, which is devoid of carbonate phases, and the lower part, reflecting changes in provenance and climate. The carbonate-free samples were probably derived from quartzite of the metamorphic basement and were deposited under a climate characterized by alternating dry and relatively wet periods. The other samples were derived from basement phyllite and were deposited under a warm, dry climate that promoted the capillary rise of Ca²⁺ and bicarbonate from a shallow water table, and therefore, the precipitation of carbonate. This part of the succession contains both calcite and dolomite. The presence of barite indicates an important concentration of SO₄ ^2? in the solution from which the CaMg(CO₃)₂ precipitated. The formation of dolomite under hypersaline conditions may be explained by bacterial degradation of organic matter, which produced CO₂ and ammonia, thereby increasing the solution alkalinity. The succession formed in an oxic environment, except for a calcite-rich level that formed under relatively reducing conditions. For this level, the large amount of calcite and the lack of dolomite indicate an alkaline environment and a very low Mg²⁺/Ca²⁺ ratio in the soil solution. These observations, coupled with the reducing conditions, indicate the availability of large amounts of degraded organic matter, probably related to a period typified by a wetter climate.     

Evolution of groundwater salinity in the unconfined aquifer of Bou-Areg, Northeastern Mediterranean coast, Morocco

The Bou-Areg plain in the Mediterranean coast at the North-eastern of Morocco is characterized by a semi-arid climate. The aquifer consists of two sedimentary formations of Plio-quaternary age: the upper formation of fine silts and the lower one of coarse silts with sand and gravels. The aquifer is underlain by marly bedrock of Miocene age that dips toward the coastal lagoon of Bou-Areg. The hydrodynamic characteristics vary between 10⁻⁴ and 10⁻³ m/s; and transmissivities range between 10⁻⁴ and 10⁻¹ m²/s. The general direction of flow is SW to NE, toward the lagoon. The aquifer is crossed by the river Selouane, which also ends in the lagoon. The groundwater is characterized by a high salinity that can reach 7.5 g/l. The highest values are observed in the upstream and in the downstream sectors of the aquifer. The temporal evolution of the physico-chemical parameters depends on the climatic conditions and piezometric variations. The analysis of the spatio-temporal distribution of the physico-chemical parameters suggests different sources of groundwater salinization: the seawater intrusion, the influence of marly gypsum-bearing terrains, and the influence of anthropogenic products as the agricultural fertilizers, which cause great nitrate concentrations that vary between 80 and 140 mg/l.     

Deep-water clastic evaporites deposition in the Messinian Adriatic foredeep (northern Apennines,Italy): did the Mediterranean ever dry out?

A new genetic facies model for deep-water clastic evaporites is presented, based on work carried out on the Messinian Gessoso-solfifera Formation of the northern Apennines during the last 15 years. This model is derived from the most recent siliciclastic turbidite models and describes the downcurrent transformations of a parent flow mainly composed of gypsum clasts. The model allows clearer comprehension of processes controlling the production and deposition of clastic evaporites, representing the most common evaporite facies of the northern Apennines, and the definition of the genetic and stratigraphic relationship with primary shallow-water evaporites formed and preserved in marginal settings. Due to the severe recrystallization processes usually affecting these deposits, petrographic and geochemical analyses are needed for a more accurate interpretation of the large spectrum of recognized gravity-driven deposits ranging from debrisflow to low-density turbidites. Almost all the laminar ‘balatino’ gypsum, previously considered a deep-water primary deposit, is here reinterpreted as the fine-grained product of high to low-density gravity flows. Facies associations permit the framing of the distribution of clastic evaporites into the complex tectonically controlled depositional settings of the Apennine foredeep basin. The Messinian Salinity Crisis occurred during an intense phase of geodynamic reorganization of the Mediterranean area that also produced the fragmentation of the former Miocene Apennine foredeep basin. In this area, primary shallow-water evaporites equivalent to the Mediterranean Lower Evaporites, apparently only formed in semi-closed thrust-top basins like the Vena del Gesso Basin. The subsequent uplift and subaerial exposure of such basins ended the evaporite precipitation and promoted a widespread phase of collapse leading to the resedimentation of the evaporites into deeper basins. Vertical facies sequences of clastic evaporites can be interpreted in terms of the complex interplay between the Messinian tectonic evolution of the Apennine thrust belt and related exhumation–erosional processes. The facies model here proposed could be helpful also for better comprehension of other different depositional and geodynamic contexts; the importance of clastic evaporites deposits has been overlooked in the study of other Mediterranean areas. Based on the Apennine basins experience, it is suggested here that evaporites diffused into the deeper portions of the Mediterranean basin may consist mainly of deep-water resedimented deposits rather than shallow-water to supratidal primary evaporites indicative of a complete basin desiccation.     

Assessment of groundwater salinity using GIS and multivariate statistics in a coastal Mediterranean aquifer

The integration of the statistical approaches and GIS tools with the hydrogeological and geological contexts allowed the assessment of the processes that cause groundwater quality deterioration in the great important deltaic aquifer in the northeastern Tunisia (Medjerda Lower Valley Aquifer). The spatial variation of the groundwater parameters and the molar ratio (Cl^?/Br^?) were also used to determine the possible impacts from seawater intrusion and from the septic tank leachate. Sixty shallow groundwater samples were collected in 2014 and analyzed for major and trace ions over an area of about 1090 km² to determine the suitability for drinking or agricultural purposes. The total dissolved solids (TDS) content ranges from 1005 to 19,254 mgl^?1 with a mean value of 3477.18 mgl^?1. The chemistry is dominated by the sodium–chloride waters (55%). Mapping of TDS, Cl^?, Na⁺, SO₄^2? and NO₃^? using kriging method shows a clear increase in salinity toward the coastline accompanied by Na⁺ and Cl^? increase which may be related to seawater intrusion and halite dissolution. Locally, higher nitrate concentration is related to the agricultural activities inducing contribution of chemical fertilizers and irrigation with treated wastewater. The saturation indices indicate that all carbonate minerals tend to reach saturation equilibrium confirming water–rock interactions, while evaporitic minerals are still in sub-saturation state and may increase the salinity of the groundwater. The principal component analysis proves the occurrence of groundwater contamination principally by seawater intrusion in the factor I (74.15%) and secondary by an anthropogenic source in the factor II (10.35%).