Next Article in Journal
Simulating the Flood Limits of Urban Rivers Embedded in the Populated City of Santa Clara, Cuba
Next Article in Special Issue
Hydrogeochemical Characteristics and Groundwater Quality in Phreatic and Confined Aquifers of the Hebei Plain, China
Previous Article in Journal
Maximizing Water Use Efficiency in Rice Farming: A Comprehensive Review of Innovative Irrigation Management Technologies
Previous Article in Special Issue
Hydrogeochemical Characteristics and Evolution of Karst Groundwater in Heilongdong Spring Basin, Northern China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessment of Hydrogeochemical Characteristics and Seawater Intrusion in Coastal Aquifers by Integrating Statistical and Graphical Techniques: Quaternary Aquifer, West Nile Delta, Egypt

1
Desert Research Center, El Matareya, Cairo 11753, Egypt
2
Geology Department, Faculty of Science, Tanta University, Tanta 31527, Egypt
3
National Water and Energy Center, United Arab Emirates University, Al Ain 15551, United Arab Emirates
*
Author to whom correspondence should be addressed.
Water 2023, 15(10), 1803; https://doi.org/10.3390/w15101803
Submission received: 29 March 2023 / Revised: 5 May 2023 / Accepted: 8 May 2023 / Published: 9 May 2023
(This article belongs to the Special Issue Groundwater Chemistry and Quality in Coastal Aquifers)

Abstract

:
The Quaternary aquifer in the western Nile Delta is threatened by seawater intrusion. Few studies have integrated diverse techniques for the assessment of seawater intrusion in this aquifer. The present study aims to determine the geochemical processes and impact of seawater intrusion on this aquifer. To accomplish this investigation, the integration of hydrogeochemical, statistical, multivariate statistical, and graphical tools were implemented on 75 groundwater samples and 5 soil samples. The physicochemical variables were analyzed using hierarchical cluster analysis (HCA), saturation index (SI), ionic ratios, ionic relationships, the seawater intrusion index (SWI) and the correlations among 16 hydrochemical parameters, to identify the influencing processes of groundwater quality in the study area. According to the statistical study, the groundwater is divided into four groups. Those are distributed, from north to south: Group1 (G1), Group2 (G2), Group4 (G4), and Group3 (G3). The samples of G1 and G2 are distinguished by Na–Cl chemical type. While G4 has two main ion associations, HCO3–Ca–Mg and Cl–SO4–Na, G3 is characterized by HCO3–Cl–SO4–Ca–Na type. The processes that affect the chemistry of the groundwater are the seawater intrusion, ion exchange, silicate and Ca-rich mineral weathering, and mineral deposition. G1 and G2 groups are primarily influenced by seawater incursion, evaporation, and the ion exchange mechanism. In addition, the weathering of silicate minerals has a substantial effect on G3 and G4 groups, resulting in the creation of carbonate minerals.

1. Introduction

Groundwater is a valuable natural resource that is essential for social and economic growth and human needs, especially in arid and semi-arid regions which suffer from shortage of precipitation and fresh surface water sources [1,2,3]. Aquifer rock material, climate conditions, and human activities in addition to seawater intrusion have a significant impact on groundwater quality and quantity in coastal aquifers [4,5,6,7].
Due to extreme use of the groundwater in coastal aquifers for industrial, agricultural, and urban expansion, seawater intrusion has come to be a global problem [8,9,10,11]. Climate change has raised atmospheric surface temperatures and melted ice caps, mountain glaciers, and polar ice sheets, which raises ocean and sea water levels, compounding the problem [12]. In addition, both the slope of the aquifer bed and the slope of the seaside influence the seawater intrusion; increasing the slope of the bed toward the sea will increase the interference [13]. Seawater intrusion hazard induces salinization of aquifers, reducing groundwater quality and causing a salinity hazard to the superimposed soils [14].
Globally, the coastal areas are occupied by about 20% of the world’s population [12]. However, due to groundwater overexploitation to face the increasing water demand associated with population expansion and fast urban development, these locations are vulnerable to seawater intrusion. Future climate change and the resulting sea level rise are both predicted to worsen this vulnerability. Nevertheless, groundwater management becomes more complicated when all these elements come together.
Groundwater salinization threatens coastal fresh water supplies such as north Kuwait [15], Kish Island, Iran [16], Thriassion Plain and Eleusis Gulf, Greece [17], China [18], Spain [19], Bangladesh [20], Mexico [21], and others. Contamination by seawater has been studied in many coastal aquifers. The assessment of this phenomenon involves determining the spatial variation of physicochemical properties such as electrical conductivity (EC), total dissolved solids (TDS), seawater mixing and groundwater level, groundwater hydrochemical analysis, and geoelectrical methods such as resistivity and electromagnetic techniques [22,23,24,25,26,27,28,29]. Several studies have been conducted to simulate the seawater intrusion using different numerical techniques, including Sefelnasr and Sherif [30], Mabrouk et al. [31], Mabrouk et al. [32], Abd-Elhamid et al. [12], Abd-Elhamid et al. [13], Sarker et al. [20], Masoud et al. [33]. The Nile Delta aquifer is one of the largest fresh groundwater sources in the world. Researchers found that the seawater intrusion front in the Nile Delta aquifer moved inland more than 100 km from the Mediterranean shoreline [34,35,36].
Eventually, the complex process of seawater intrusion is influenced by, among others, shoreline geomorphology, hydrogeochemical reactions, biological processes, and aquifer dynamics. Saltwater intrusion, mostly in unconfined coastal aquifers, is the main cause of salinization [19]. In confined or semi-confined aquifers, complicated hydrogeological processes such as water–rock interaction and mixing with different water sources may impact groundwater quality and geochemical evolution [30].
Understanding the hydrochemistry is essential for determining the groundwater quality and its suitability for various applications, and therefore, determining sound management options for this resource as well as efficient mitigation measures when required. Diverse statistical and graphical techniques, such as pattern diagrams and GIS, can also be used to characterize seawater incursion. However, evaluation of the groundwater chemistry in the region based on integration of hydrogeochemical methods, statistical analysis, and advanced geostatistical techniques is still weak. Such integration provides the means of investigating such complex coastal aquifers, and the results may function as the foundation for effective groundwater management.
The purpose of this study is to evaluate the anthropogenic impacts and saltwater intrusion on the geochemical constituents and processes governing groundwater in the study area. To achieve this objective, groundwater samples of soil and water were collected and analyzed. For organizing groundwater samples, managing this volume of data, and providing accurate interpretations, graphical representations of the data were accompanied by multiple statistical approaches and multiple indices.

2. Materials and Methods

2.1. Study Area

The study area is a part of the Northwest Rosetta branch. It lies between 30°42′ to 31°28′ N and 29°45′ to 30°48′ E (Figure 1). It is bounded to the north by the Mediterranean Sea, to the east by the branch of Rosetta, and to the west by El Nubaryia Canal. The study area has an arid to semi-arid climate with an average annual temperature of 30.4 °C and an average annual humidity of 59% [37]. The rainfall ranges from 50 mm/a in the south to about 200 mm/a in the north. The evapotranspiration ranges from 1648 in the north to 1680 mm/y in the south [38]. Geomorphologically and hydrogeologically, the investigated area was subjected to different studies such as Allam et al. [39], Dawoud et al. [40], Morsy [41], El Fakharany and Hagran [42], and Salem et al. [43]. Young alluvial plains characterize the research region. It features a nearly elevated ground surface that ranges from −3 m (amsl) in the northwest to +6 m (amsl) in the southeast [41]. The area is characterized by sand dunes, aeolian deposits that occasionally occupy the surface, brackish lakes, and water-logged areas that occur in several places [44]. The slope is very smooth, representing one meter per 10 km to the north direction [44,45]. The Quaternary aquifer is comprised of Mit Ghamr Formation on the bottom and Bilqas Formation on the top [43]. The Mit Ghamr Formation is made up of sand and gravel with clay intercalations of continental origin. The aquifer layer is semiconfined with a thickness reaches its maximum near the Rosetta branch (about 850 m) and decreases rapidly towards the west, where it notably decreases near the El Nubaryia canal (<200 m) [36,46]. The Holocene-aged cap layer of the Bilqas Formation consists of clay and silt and is less than 20 m thick in the south and up to 50 m thick in the north. The Nile Delta aquifer system is recharged in the flood plain by infiltration from surface water (especially from irrigation systems), rainfall, and by leakage from subsurface drainage water in the cultivated lowlands [40]. Groundwater discharge takes place either naturally through the outflow into the drainage system, via evapotranspiration and groundwater baseflow, or artificially through direct well extraction, which is primary discharge component of the aquifer [47]. The groundwater level within the study area ranges from 1 m in the north to 4 m close to El Nubaryia canal. The general groundwater flow direction is from south to north and from southwest to northeast.

2.2. Sampling and Analytical Techniques

The study of the hydrogeochemical aspects of groundwater in the study area was mainly based on the results of the chemical analyses of 75 groundwater samples that were collected from the drilled wells during March, June, and October 2018 (Figure 1). In addition, five soil water samples were also collected from the subsoil drainage system and analyzed to estimate the origin of dissolved ions. The sampled wells have depths varying from 15 to 200 m. The chemical analyses were performed in the laboratory of the Geology Department, Faculty of Science, Tanta University, and in the Center of Scientific Research and Measurements, Tanta University. The analysis techniques are shown in Table 1. While water samples were being analyzed, a number of quality assurance and quality control procedures were carried out. The validation of analytical methodologies was accomplished by calibrating devices and assessing the precision of the samples being analyzed. Charge error balance (CBE) was calculated to ensure that the analytical error of the analyzed ions concentration in meq/L−1 was accurate [48,49]. The CBE of every examined sample was found to be within the permissible range of ±5%. The obtained chemical data are listed in Table 2.

2.3. Cluster Analysis and Saturation Index (SI)

Multivariate statistical analyses were used to provide a quantitative measure of water quality parameters relative to each other and to estimate correlations between chemical parameters and groundwater samples [50,51,52,53,54,55,56,57]. Hierarchical cluster analysis HCA using Ward’s linking method and squared Euclidean distance as a dissimilarity metric was used to group parameters by multivariate similarities.
The saturation percentage, specified as SI = log (IAP/Ks), defines the degree of mineral saturation in water. IAP is the ionic activity product of the proper ions, and Ks is the constant of saturation under a specific temperature. The equilibrium distribution of minerals in an aqueous solution and the state of saturation is determined using the hydrochemical program PHREEQC Interactive Program [58]. An SI < 0 means that the mineral’s groundwater is undersaturated. An SI > 0 indicates that the groundwater is supersaturated with this mineral phase and thus incapable of dissolving more of the mineral.

2.4. Seawater Intrusion Quality Index (SWI)

To identify the threat of seawater intrusion, Tomaszkiewicz et al. [59] created the groundwater quality index for seawater intrusion (GQISWI) (Equation (1)). This model translates the Piper diagram information GQIpiper (mix) (freshwater seawater mixing index of Piper diagram) and the seawater fraction index GQIfsea to create a new two-stage numerical indicator for seawater intrusion [59].
G Q I S W I = G Q I p i p e r ( m i x ) + G Q I f s e a 2
The first stage included the calculation of GQIpiper (mix) as follows:
G Q I P i p e r m i x ( m e q / l ) = C a + 2 + M g + 2 T o t a l   c a t i o n s + H C O 3 T o t a l   a n i o n s × 50
The second stage involves computing the seawater fraction index GQIfsea (Equation (3)) from the seawater fraction (fsea) (Equation (4)). GQIfsea was then used to classify the SWI. These values vary from 0 to 100, the fresher waters of which have lower values.
G Q I f s e a = 1 f s e a × 100
The following formula can be used to measure the seawater fraction (fsea) [60]:
f s e a = m C l s a m p l e m C l ( f r e s h w a t e r ) m C l s e a w a t e r m C l ( f r e s h w a t e r )
where mCl (sample) is the Cl concentration of the sample, mCl (seawater) is the Cl concentration of the Mediterranean Sea (603 meq/L), and mCl (freshwater) represents the Cl concentration of the freshwater (0.48 meq/L).
The difference in hydrochemistry between seawater and freshwater was classified by concentrations of anions and cations. Indeed, the freshwater is defined by an abundance of calcium and bicarbonate, while the seawater is classified by its dominance of chloride and sodium [61]. The concentration of chloride is the base which determines the mixing of these two types of water, which is clarified by the conservative presence of anions [62].

3. Results and Discussions

3.1. Hydrochemical Characteristics and Spatial Distribution

The hydrochemical characteristics of groundwater are discussed with respect to the study area and according to the results of water samples chemical analysis. As shown in Table 2, the groundwater pH values range from slightly acidic (5.05) to slightly alkaline (8.7) averaging 7.1. The differences in pH values are primarily due to the aquifer–rock interaction and the distance from the sea [63]. Based on the classification of Winslow and Kister [64], 53%, 21%, 15%, and 11% of the groundwater samples are fresh, slightly, moderately, and very saline water, respectively (Table 3). The EC of groundwater in the study area ranges from 0.38 to 55.3 mS/cm with an average of 7.1 mS/cm (Table 2). TDS levels range from 190 to 27680 mg/L, averaging 3806.2 (Table 2). The Rosetta branch and many irrigation canals nearby produce fresh to slightly salty water in the southern and central study area. Near irrigation canals, Nile water seepage dilutes groundwater. Figure 2 shows that seawater incursion increases salinity northward, where the shore is 30 km away.
As demonstrated in Figure 3a–c, Ca, Mg, and Na concentrations were lowest in the south and center regions and rose northward. It is largely seawater interference. The southern and central parts may have lower cations due to Nile water seepage and irrigation.
Figure 3d shows that the concentration of bicarbonate increases from south and northwest towards the central parts. Bicarbonates are introduced to natural water through recharging with meteoric water and interaction with sediments rich in either carbonate minerals or organic matter. The dissolution of CaCO3 or the decay of organic matter increase in the presence of CO2 leading to the formation of more soluble bicarbonate [65]. SO4 and Cl spatial distributions showed that they have almost the same pattern. The lower values are located in the south and increase towards the north (Figure 3e,f). Such distribution patterns could be related to seawater intrusion effect.
Figure 4a–f shows Al, Mn, Ba, Fe, Sr, and Si spatial distribution maps. Al concentration increases southward and decreases northward (Figure 4a). Mn is highest in the northeast and declines westward to its minimum at El Nubaryia Canal (Figure 4b). Ba, Fe, and Sr have similar distributions. The concentration increases northeastward from the south (Figure 4c–e). Si levels are higher in central areas from northeastern to southwestern directions and lower along the western boundary (Figure 4f).

3.2. Hydrochemical Classification

3.2.1. Cluster Analysis

Based on the HCA procedures, four groups were identified. Figure 5a shows that G2, G3, and G4 are related, but G1 samples are distinct. This is primarily due to their occurrences in the south and center of the study area, which are relatively located far from the sea. On the other hand, G1 is located in the northern region near to the sea, thus, its chemical composition was significantly affected by the seawater.
Table 4 shows the mean hydrochemical data values for each sample group to explain their features. Figure 6 depicts group dispersion by location.
Scholler’s diagram is commonly used to correlate many of the chemical water compositions. The average hydrochemical data of the four groups (Table 4) are plotted on the semi-logarithmic diagrams (Figure 5b). The samples of G1, G2 and G4 have Na–Cl water type (Na > Ca > Mg > and Cl > SO4 > HCO3) with an average TDS values 22,686.7 mg/L, 6486.4 mg/L and 1829.1 mg/L, respectively. The spatial distribution of G1, G2 and G4 (Figure 6) indicated that these groups are influenced by seawater intrusion, while G3 samples have Na–Cl (Na > Ca > Mg and Cl > HCO3 > SO4) water type and lower concentrations of TDS (average 597.9 mg/L) and ions compared to the other groups. G3 is distributed in the southern part of the study area and its origin might be related to direct seepage from the irrigation canals.

3.2.2. Piper’s (Trilinear) Diagram

The collected samples were plotted in the trilinear diagram [66] (Figure 7) according to the HCA classification (Figure 5a). In most groundwater samples, Na is the major cation and Cl is the main anion. Generally, the samples of G1, G2 and G4 were plotted on sub-area 7 (Figure 7) indicating primary salinity character, where NaCl and Na2SO4 salts dominate the chemical properties. This is largely due to the impact of seawater intrusion and soil salinity. Most samples of G4, which occupied the southern parts, are located in sub-area 9 (Figure 7) indicating that none of the cation–anion pairs exceeds 50%. In the right triangle, samples of G3 and 4 show an increase in concentrations of either SO4 or HCO3 and decrease of Cl. In the lift triangle, due to ion exchange process, concentrations of Ca and Mg increase in groups 3 and 4 and Na decrease [67].

3.2.3. Statistical Ions Classification

Figure 8 shows sample group ion-association dendrograms. Na, Cl, and TDS in groundwater samples from groups G1 and G2 (Figure 8a,b) show seawater intrusion in the northern section of the research area (Figure 6). Assuming G3, HCO3–Cl–SO4–Ca–Na is the major ion association in the southern region, indicating freshwater origin and evaporation (Figure 8c). Figure 8d demonstrates ion-association in group four samples. HCO3–Ca–Mg-minors and Cl–SO4–Na showed that freshwater impacted by seawater intrusion and exchange processes is widespread. Al, Mn, Ba, Fe, Sr, Si, and NO3 are linked in a short distance, suggesting a shared origin that may affect fertilizer or water–rock interactions (Figure 8).

3.3. Saturation Index (SI)

Table 5 shows the statistical summary of dissolved mineral saturation indicators. The dissolved minerals were classified into five different groups: carbonate, sulphate, chlorides, oxides, and silicate minerals. Detected carbonate minerals include aragonite, calcite, dolomite, and witherite. Unlike groups 1, 2, and 4, group 3 is saturated with aragonite, calcite, and dolomite. All groups are undersaturation with respect to witherite. As the investigated aquifer is composed of clastic sediments, production of bicarbonate ion is due to the formation of CO2 from decay of the organic matter (Equation (5)) [60]. Such conditions could lead to saturation conditions in the southern part where group 3 samples are located but the unsaturation conditions in groups 1, 2, and 4 might be related to carbonate minerals deposition where the increase in sulfate and chloride minerals due to seawater intrusion in the northern parts is accompanied by deposition of the carbonate minerals [36,68].
CO2 + H2O → H2CO3 + H → HCO3 + H
SO4−2 + 2C + 2 H2O → H2S +2 HCO3
The studied groundwater is mostly unsaturated with respect to anhydrite and gypsum minerals, whereas most of G3, which is located in the southern part of the area, was saturated with respect to barite. Hydrochemical studies of major and trace elements made by Drevaliene et al. [69] showed that the probable barium source in groundwater is the dissolution of witherite. According to Table 5, all four groups of groundwater samples showed an undersaturation state with the dissolved chloride minerals (Halite and Sylvite).
G1, G2, G3, and G4 are undersaturated with respect to pyrolusite but most of them are saturated with hematite (Table 5) due to the reductive dissolution of Fe and Mn. The dissolution of Fe and Mn oxides in the aquifer is at the cost of the oxidation of organic matter [70,71], and the bicarbonate will increase. The reaction Equations (7) and (8) are as follows:
CH2O + 2MnO2 + 3H+ = 2Mn2+ + HCO3 + 2H2O
CH2O + 4Fe(OH)3 + 7H+ = 4Fe2+ + HCO3 + 10H2O
The encountered dissolved silicate minerals are albite, anorthite, K-feldspar, kaolinite, and quartz. Groundwater samples are undersaturated with albite and anorthite. Most of G1, G2, and G4 are in the undersaturation state with K-feldspar but G3 is mostly saturated with this mineral. G2, G3, G4, and most of G1 are saturated with quartz and kaolinite (Table 5). Following Hosono et al. [72], incongruent dissolution reactions of silicate weathering in the study area can be hypothesized by following Equations (9) and (10):
2NaAlSi3O8 (Albite) + 2CO2 + 11H2O → Al2Si2O5(OH)4 (Kaolinite) + 2Na+ + 2HCO3 + 4H4SiO4
CaAl2Si2O8 (Anorthite) + 2CO2 + 3H2O → Al2Si2O5(OH)4 (Kaolinite) + 2Ca2+ + 2HCO3

3.4. Hydrogeochemical Processes

Ions-TDS Relationships

  • Gibbs diagram
Gibbs diagrams show the groundwater chemistry’s main natural process. It shows how rock weathering, mineral precipitation, and evaporation affect groundwater chemistry [73]. High concentration of Ca2+ and HCO3 reveals rock–water interaction. On the other hand, the presence of high concentration of Na+ and Cl indicates evaporation processes or seawater intrusion [74,75]. TDS concentrations are plotted against the weight ratios of Na/(Na + Ca) for cations (Figure 9a) and the weight ratios of Cl/(Cl + HCO3) for anions (Figure 9b). The major processes governing groundwater chemistry are evaporation/seawater intrusion and rock interaction. Most G3 samples have dissolved HCO3 and Ca via water-rock interaction. Evaporation and seawater intrusion supply Na and Cl to Groups 1, 2, 4, and soil water.
  • TDS vs. Ca, Mg, Na + K, Cl, HCO3, and SO4
TDS vs. Ca, Mg and Na + K relationship (Figure 10a–c) are very strong (r2 = 0.9, 0.9, and 0.9, respectively, Table 6). For TDS relationships with Ca and Mg, samples were plotted under the mixing line which indicate ion-exchange process due to seawater intrusion, where Na replaces either Ca or/and Mg as shown in Equations (11) and (12). TDS vs. Cl relationship shows a very strong correlation (r2 = 0.9, Table 6, Figure 10d). For the relationships between TDS and Na + K and Cl, most plotted groundwater samples are clustered along the irrigation water–seawater mixing line suggesting that the variance in concentrations of theses ions is due to the mixing between irrigation canal water and seawater. TDS vs. HCO3 (Figure 10e) showed weak inverse relationship (r2 = −0.1, Table 6), which means that freshwater seepage is not only the primary source of dissolved carbonates but may also be the reason for the weathering of silicate minerals. as shown in Equations (11) and (12) [76].
2Na+ + Ca − X2 → 2Na − X + Ca2+
2Na+ + Mg − X2 → 2Na − X + Mg2+
where x is an ion exchanger.
On the other hand, TDS–SO4 (Figure 10f) have a strong relationship (r2 = 0.7, Table 6). The origin of sulfates is not linked to the trend of surface water–seawater mixing but possibly related to soil origin, especially in low salinity water. The relationships between TDS vs. Ba, Fe and Sr (r2 equals 0.5, 0.4, and 0.5, respectively, Table 6) showed intermediate correlation, on the other hand it showed weak relationships with Al, Mn, and Si (r2 equals −0.2, 0.1, and 0.1 respectively, Table 6). Figure 10a–f shows that plotted soil water samples and indicates that the soil hydrochemical processes are very effective in the hydrochemical characteristics of the studied low salinity groundwater (G3 and G4) where it acts as an important origin of the dissolved ions.

3.5. Seawater Intrusion Quality Index (SWI)

In ArcGIS 10 framework, the calculated GQISWI results are transformed into a map using a kriging interpolation function. The GQISWI value varies from 17.8 to 82.3 as shown in Figure 11. According to Table 7, the groundwater of the study area is classified into saline, mixed, and fresh water. The coastal area in the northern parts has low values and reflects the salinization of groundwater.

4. Conclusions

In this work, the seawater intrusion and the processes governing the groundwater of the Quaternary aquifer of the Nile Delta were assessed using hydrochemical analysis in conjunction with several multivariate statistical methods and graphical approaches. It was concluded that groundwater chemistry is influenced by the evaporation process, seawater intrusion, ion exchange process, dissolution and weathering of Ca-rich minerals, and weathering of silicate minerals. The groundwater samples are divided into four groups in accordance with HCA, which were identified as G1, G2, G3, and G4. G1 and G2, which are found in the northern regions near to the Mediterranean Sea, are primarily controlled by the evaporation process and seawater intrusion. High Cl and Na concentrations and low HCO3 concentrations mostly serve to demonstrate this. Additionally, seawater intrusion is joined by another process such as cation exchange where Na replaces Ca and/or Mg. On the other hand, in the case of G3 and G4, the weathering of silicate minerals led to the development of carbonate minerals. Most of the groundwater samples from the four groups included more calcium than HCO3 and SO4; this is mostly because calcium-rich minerals dissolve and weather.
Climate change and its parameters, sea-level rise, extensive groundwater pumping to meet various needs, land-use changes due to population growth and urbanization, changes in rainfall patterns, and reduced groundwater recharge all threaten saltwater intrusion coastal aquifers.
Since the coastal aquifer in the Nile Delta is responsive to climate change and groundwater extraction, sustainable water resource management is needed to protect it from saltwater intrusion. This research provides a useful assessment tool for seawater intrusion variability in coastal aquifers. The results can improve understanding of present and future saltwater intrusion in the investigated area and provide specific adaptive keys. As a response to the aquifer stresses, adaptation styles to any of the given hazards or challenges may be designed. Thus, the adaptation in the case of the Quaternary aquifer in the Nile Delta should be transformational, where changes in the fundamental attributes of a system (especially pumping) must be adopted. This may not preclude combining different adaptation strategies or using biodiversity and environmental services as part of an adaptation plan. Constant monitoring of the saline interface is essential for maintaining the aquifer’s equilibrium and determining the most effective aquifer management strategy. An assessment study for the saltwater intrusion under different climate and anthropogenic activities is still needed for the aquifer.

Author Contributions

Conceptualization, S.S.H., Z.E.S. and A.S.; methodology, S.S.H. and Z.E.S.; software, S.S.H.; validation, S.S.H., Z.E.S. and A.S.; formal analysis, S.S.H., Z.E.S. and A.S.; investigation, S.S.H., Z.E.S. and A.S.; resources, S.S.H. and Z.E.S.; data curation, S.S.H., Z.E.S. and A.S.; writing—original draft preparation, S.S.H.; writing—review and editing, Z.E.S. and A.S.; visualization, S.S.H.; supervision, Z.E.S. and A.S.; project administration, Z.E.S. and A.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, X.; Miao, J.; Hu, B.X.; Liu, H.; Zhang, H.; Ma, Z. Hydrogeochemical characterization and groundwater quality assessment in intruded coastal brine aquifers (Laizhou Bay, China). Environ. Sci. Pollut. Res. 2017, 24, 21073–21090. [Google Scholar] [CrossRef] [PubMed]
  2. Makri, P.; Stathopoulou, E.; Hermides, D.; Kontakiotis, G.; Zarkogiannis, S.D.; Skilodimou, H.D.; Bathrellos, G.D.; Antonarakou, A.; Scoullos, M. The Environmental Impact of a Complex Hydrogeological System on Hydrocarbon-Pollutants’ Natural Attenuation: The Case of the Coastal Aquifers in Eleusis, West Attica, Greece. J. Mar. Sci. Eng. 2020, 8, 1018. [Google Scholar] [CrossRef]
  3. El Osta, M.; Masoud, M.; Alqarawy, A.; Elsayed, S.; Gad, M. Groundwater Suitability for Drinking and Irrigation Using Water Quality Indices and Multivariate Modeling in Makkah Al-Mukarramah Province, Saudi Arabia. Water 2022, 14, 483. [Google Scholar] [CrossRef]
  4. Salem, Z.E.-S.; Abdelrahman, k.; Kováčiková, S.; Badran, O.M. Use of various statistical techniques to assess the vertical and lateral change in the groundwater chemistry of Quaternary aquifer in an irrigated highly populated area. J. King Saud Univ. -Sci. 2021, 33, 101556. [Google Scholar] [CrossRef]
  5. Abu Salem, H.S.; Gemail, K.S.; Junakova, N.; Ibrahim, A.; Nosair, A.M. An Integrated Approach for Deciphering Hydrogeochemical Processes during Seawater Intrusion in Coastal Aquifers. Water 2022, 14, 1165. [Google Scholar] [CrossRef]
  6. El Yousfi, Y.; Himi, M.; El Ouarghi, H.; Elgettafi, M.; Benyoussef, S.; Gueddari, H.; Aqnouy, M.; Salhi, A.; Alitane, A. Hydrogeochemical and statistical approach to characterize groundwater salinity in the Ghiss-Nekkor coastal aquifers in the Al Hoceima province, Morocco. Groundw. Sustain. Dev. 2022, 19, 100818. [Google Scholar] [CrossRef]
  7. Makri, P.; Hermides, D.; Kontakiotis, G.; Zarkogiannis, S.D.; Besiou, E.; Janjuhah, H.T.; Antonarakou, A. Integrated Ecological Assessment of Heavily Polluted Sedimentary Basin within the Broader Industrialized Area of Thriassion Plain (Western Attica, Greece). Water 2022, 14, 382. [Google Scholar] [CrossRef]
  8. Dibaj, M.; Javadi, A.A.; Akrami, M.; Ke, K.-Y.; Farmani, R.; Tan, Y.-C.; Chen, A.S. Modelling seawater intrusion in the Pingtung coastal aquifer in Taiwan, under the influence of sea-level rise and changing abstraction regime. Hydrogeol. J. 2020, 28, 2085–2103. [Google Scholar] [CrossRef]
  9. Prusty, P.; Farooq, S.H. Seawater intrusion in the coastal aquifers of India—A review. HydroResearch 2020, 3, 61–74. [Google Scholar] [CrossRef]
  10. Rajendiran, T.; Sabarathinam, C.; Chandrasekar, T.; Panda, B.; Mathivanan, M.; Nagappan, G.; Natesan, D.; Ghai, M.; Kumar Singh, D.; Alagappan, R. Geochemical variations due to salinization in groundwater along the southeast coast of India. SN Appl. Sci. 2021, 3, 581. [Google Scholar] [CrossRef]
  11. Hajji, S.; Allouche, N.; Bouri, S.; Aljuaid, A.M.; Hachicha, W. Assessment of Seawater Intrusion in Coastal Aquifers Using Multivariate Statistical Analyses and Hydrochemical Facies Evolution-Based Model. Int. J. Environ. Res. Public Health 2021, 19, 155. [Google Scholar] [CrossRef] [PubMed]
  12. Abd-Elhamid, H.F.; Abd-Elaty, I.; Hussain, M.S. Mitigation of seawater intrusion in coastal aquifers using coastal earth fill considering future sea level rise. Environ. Sci. Pollut. Res. 2020, 27, 23234–23245. [Google Scholar] [CrossRef] [PubMed]
  13. Abd-Elhamid, H.F.; Abd-Elaty, I.; Sherif, M.M. Effects of Aquifer Bed Slope and Sea Level on Saltwater Intrusion in Coastal Aquifers. Hydrology 2019, 7, 5. [Google Scholar] [CrossRef]
  14. Ding, Z.; Koriem, M.A.; Ibrahim, S.M.; Antar, A.S.; Ewis, M.A.; He, Z.; Kheir, A.M.S. Seawater intrusion impacts on groundwater and soil quality in the northern part of the Nile Delta, Egypt. Environ. Earth Sci. 2020, 79, 313. [Google Scholar] [CrossRef]
  15. Yihdego, Y.; Al-Weshah, R.A. Assessment and Prediction of Saline Sea Water Transport in Groundwater Using 3-D Numerical Modelling. Environ. Process. 2017, 4, 49–73. [Google Scholar] [CrossRef]
  16. Mahmoodzadeh, D.; Ketabchi, H.; Ataie-Ashtiani, B.; Simmons, C.T. Conceptualization of a fresh groundwater lens influenced by climate change: A modeling study of an arid-region island in the Persian Gulf, Iran. J. Hydrol. 2014, 519, 399–413. [Google Scholar] [CrossRef]
  17. Hermides, D.; Makri, P.; Kontakiotis, G.; Antonarakou, A. Advances in the Coastal and Submarine Groundwater Processes: Controls and Environmental Impact on the Thriassion Plain and Eleusis Gulf (Attica, Greece). J. Mar. Sci. Eng. 2020, 8, 944. [Google Scholar] [CrossRef]
  18. Han, G.; Liang, C.-Z.; Chung, T.-S.; Weber, M.; Staudt, C.; Maletzko, C. Combination of forward osmosis (FO) process with coagulation/flocculation (CF) for potential treatment of textile wastewater. Water Res. 2016, 91, 361–370. [Google Scholar] [CrossRef]
  19. Llamas, M.R.; Martínez-Santos, P. Intensive Groundwater Use: Silent Revolution and Potential Source of Social Conflicts. J. Water Resour. Plan. Manag. 2005, 131, 337–341. [Google Scholar] [CrossRef]
  20. Sarker, M.M.R.; Hermans, T.; Van Camp, M.; Hossain, D.; Islam, M.; Ahmed, N.; Bhuiyan, M.A.Q.; Karim, M.M.; Walraevens, K. Identifying the Major Hydrogeochemical Factors Governing Groundwater Chemistry in the Coastal Aquifers of Southwest Bangladesh Using Statistical Analysis. Hydrology 2022, 9, 20. [Google Scholar] [CrossRef]
  21. Escolero, O.; Marin, L.E.; Domínguez-Mariani, E.; Torres-Onofre, S. Dynamic of the freshwater–saltwater interface in a karstic aquifer under extraordinary recharge action: The Merida Yucatan case study. Environ. Geol. 2007, 51, 719–723. [Google Scholar] [CrossRef]
  22. Vengadesan, M.; Lakshmanan, E. Management of Coastal Groundwater Resources. In Coastal Management; Elsevier: Amsterdam, The Netherlands, 2019; pp. 383–397. ISBN 978-0-12-810473-6. [Google Scholar]
  23. Torres Martínez, J.A.; Mora, A.; Ramos Leal, J.A.; Morán Ramírez, J.; Arango Galván, C.; Mahlknecht, J. Constraining a density-dependent flow model with the transient electromagnetic method in a coastal aquifer in Mexico to assess seawater intrusion. Hydrogeol. J. 2019, 27, 2955–2972. [Google Scholar] [CrossRef]
  24. Senthilkumar, S.; Vinodh, K.; Johnson Babu, G.; Gowtham, B.; Arulprakasam, V. Integrated seawater intrusion study of coastal region of Thiruvallur district, Tamil Nadu, South India. Appl. Water Sci. 2019, 9, 124. [Google Scholar] [CrossRef]
  25. Vann, S.; Puttiwongrak, A.; Suteerasak, T.; Koedsin, W. Delineation of Seawater Intrusion Using Geo-Electrical Survey in a Coastal Aquifer of Kamala Beach, Phuket, Thailand. Water 2020, 12, 506. [Google Scholar] [CrossRef]
  26. Shin, J.; Hwang, S. A Borehole-Based Approach for Seawater Intrusion in Heterogeneous Coastal Aquifers, Eastern Part of Jeju Island, Korea. Water 2020, 12, 609. [Google Scholar] [CrossRef]
  27. Hermans, T.; Paepen, M. Combined Inversion of Land and Marine Electrical Resistivity Tomography for Submarine Groundwater Discharge and Saltwater Intrusion Characterization. Geophys. Res. Lett. 2020, 47. [Google Scholar] [CrossRef]
  28. Sae-Ju, J.; Chotpantarat, S.; Thitimakorn, T. Hydrochemical, geophysical and multivariate statistical investigation of the seawater intrusion in the coastal aquifer at Phetchaburi Province, Thailand. J. Asian Earth Sci. 2020, 191, 104165. [Google Scholar] [CrossRef]
  29. Palacios, A.; Ledo, J.J.; Linde, N.; Luquot, L.; Bellmunt, F.; Folch, A.; Marcuello, A.; Queralt, P.; Pezard, P.A.; Martínez, L.; et al. Time-lapse cross-hole electrical resistivity tomography (CHERT) for monitoring seawater intrusion dynamics in a Mediterranean aquifer. Hydrol. Earth Syst. Sci. 2020, 24, 2121–2139. [Google Scholar] [CrossRef]
  30. Sefelnasr, A.; Sherif, M. Impacts of Seawater Rise on Seawater Intrusion in the Nile Delta Aquifer, Egypt. Groundwater 2014, 52, 264–276. [Google Scholar] [CrossRef]
  31. Mabrouk, M.; Jonoski, A.; HP Oude Essink, G.; Uhlenbrook, S. Impacts of Sea Level Rise and Groundwater Extraction Scenarios on Fresh Groundwater Resources in the Nile Delta Governorates, Egypt. Water 2018, 10, 1690. [Google Scholar] [CrossRef]
  32. Mabrouk, M.; Jonoski, A.; Oude Essink, G.; Uhlenbrook, S. Assessing the Fresh–Saline Groundwater Distribution in the Nile Delta Aquifer Using a 3D Variable-Density Groundwater Flow Model. Water 2019, 11, 1946. [Google Scholar] [CrossRef]
  33. Masoud, M.; El Osta, M.; Alqarawy, A.; Elsayed, S.; Gad, M. Evaluation of groundwater quality for agricultural under different conditions using water quality indices, partial least squares regression models, and GIS approaches. Appl. Water Sci. 2022, 12, 244. [Google Scholar] [CrossRef]
  34. Sherif, M.; Sefelnasr, A.; Javadi, A. Incorporating the concept of equivalent freshwater head in successive horizontal simulations of seawater intrusion in the Nile Delta aquifer, Egypt. J. Hydrol. 2012, 464–465, 186–198. [Google Scholar] [CrossRef]
  35. Salem, Z.E.; Al Temamy, A.M.; Salah, M.K.; Kassab, M. Origin and characteristics of brackish groundwater in Abu Madi coastal area, Northern Nile Delta, Egypt. Estuar. Coast. Shelf Sci. 2016, 178, 21–35. [Google Scholar] [CrossRef]
  36. Salem, Z.E.-S.; Osman, O.M. Use of major ions to evaluate the hydrogeochemistry of groundwater influenced by reclamation and seawater intrusion, West Nile Delta, Egypt. Environ. Sci. Pollut. Res. 2017, 24, 3675–3704. [Google Scholar] [CrossRef]
  37. El Deek, M.; Nessim, R. Organic pollutants in Abu Qir Bay. J. High Inst. Public Health 1995, 25, 221–226. [Google Scholar]
  38. Armanuos, A.M.; Negm, A. Integrated Groundwater Modeling for Simulation of Saltwater Intrusion in the Nile Delta Aquifer, Egypt. In Groundwater in the Nile Delta; Negm, A.M., Ed.; The Handbook of Environmental Chemistry; Springer International Publishing: Cham, Switzerland, 2019; pp. 489–544. ISBN 978-3-319-94283-4. [Google Scholar]
  39. Allam, A.R.; Saaf, E.-J.; Dawoud, M.A. Desalination of brackish groundwater in Egypt. Desalination 2003, 152, 19–26. [Google Scholar] [CrossRef]
  40. Dawoud, M.A.; Darwish, M.M.; El-Kady, M.M. GIS-Based Groundwater Management Model for Western Nile Delta. Water Resour. Manag. 2005, 19, 585–604. [Google Scholar] [CrossRef]
  41. Morsy, W. Environmental Management of Groundwater Resources in the Nile Delta Region. Ph.D. Thesis, Cairo University, Cairo, Egypt, 2009. [Google Scholar]
  42. Awad, S.R.; El Fakharany, M.A.; Hagran, N.M. Environmental Impact on Water Resources at the Northwestern Part of the Nile Delta, Egypt. J. Am. Sci. 2015, 11, 1–11. [Google Scholar]
  43. Salem, Z.E.; Sefelnasr, A.M.; Hasan, S.S. Assessment of groundwater vulnerability for pollution using DRASTIC Index, young alluvial plain, Western Nile Delta, Egypt. Arab. J. Geosci. 2019, 12, 727. [Google Scholar] [CrossRef]
  44. USGS. Earth Resources Observation And Science (EROS) Center Shuttle Radar Topography Mission (SRTM) 1 Arc-Second Global; USGS: Reston, VA, USA, 2017. [CrossRef]
  45. Salem, Z.E.; Hasan, S.S. Use of GALDIT model and HFE-Diagram to assess seawater intrusion vulnerability in West Nile Delta, Egypt. Arab. J. Geosci. 2021, 14, 1318. [Google Scholar] [CrossRef]
  46. El-Arabi, N.E.; Morsy, W.S. Applying integrated ground- and surface-water management (case study: Nubaryia Basin, West Delta, Egypt). J. Am. Sci. 2013, 9, 43–53. [Google Scholar]
  47. REGWA. The Groundwater in a Study Area Southern of Wadi El Farigh; REGWA: Cairo, Egypt, 1993; p. 74. [Google Scholar]
  48. Domenico, P.A.; Schwartz, F.W. Physical and Chemical Hydrogeology, 2nd ed.; Wiley: Hoboken, NJ, USA, 1998; ISBN 978-0-471-59762-9. [Google Scholar]
  49. Deutsch, W.J.; Siegel, R. Groundwater Geochemistry: Fundamentals and Applications to Contamination; CRC Press: Boca Raton, FL, USA, 1997; ISBN 978-1-00-306994-2. [Google Scholar]
  50. HACH. Chemical Procedures Explained; Hach Technical Center for Applied Analytical Chemistry: Loveland, CO, USA, 1990. [Google Scholar]
  51. US EPA. Method 200.7: Determination of Metals and Trace Elements in Water and Wastes by Inductively Coupled Plasma-Atomic Emission Spectrometry. Available online: https://www.epa.gov/esam/method-2007-determination-metals-and-trace-elements-water-and-wastes-inductively-coupled (accessed on 7 March 2023).
  52. Okiongbo, K.S.; Douglas, R.K. Evaluation of major factors influencing the geochemistry of groundwater using graphical and multivariate statistical methods in Yenagoa city, Southern Nigeria. Appl. Water Sci. 2015, 5, 27–37. [Google Scholar] [CrossRef]
  53. Alfy, M.E.; Lashin, A.; Faraj, T.; Alataway, A.; Tarawneh, Q.; Al-Bassam, A. Quantitative hydro-geophysical analysis of a complex structural karst aquifer in Eastern Saudi Arabia. Sci. Rep. 2019, 9, 2825. [Google Scholar] [CrossRef] [PubMed]
  54. Dişli, E.; Gülyüz, N. Hydrogeochemical investigation of an epithermal mineralization bearing basin using multivariate statistical techniques and isotopic evidence of groundwater: Kestanelik Sub-Basin, Lapseki, Turkey. Geochemistry 2020, 80, 125661. [Google Scholar] [CrossRef]
  55. Ren, X.; Li, P.; He, X.; Su, F.; Elumalai, V. Hydrogeochemical Processes Affecting Groundwater Chemistry in the Central Part of the Guanzhong Basin, China. Arch. Environ. Contam. Toxicol. 2021, 80, 74–91. [Google Scholar] [CrossRef]
  56. Li, P.; Tian, R.; Liu, R. Solute Geochemistry and Multivariate Analysis of Water Quality in the Guohua Phosphorite Mine, Guizhou Province, China. Expo Health 2019, 11, 81–94. [Google Scholar] [CrossRef]
  57. Wu, J.; Li, P.; Wang, D.; Ren, X.; Wei, M. Statistical and multivariate statistical techniques to trace the sources and affecting factors of groundwater pollution in a rapidly growing city on the Chinese Loess Plateau. Hum. Ecol. Risk Assess. Int. J. 2020, 26, 1603–1621. [Google Scholar] [CrossRef]
  58. Charlton, S.R.; Parkhurst, D.L. PHREEQCI—A Graphical User Interface to the Geochemical Model PHREEQC; Fact Sheet; U.S. Geological Survey: Reston, VA, USA, 2002; Volume 031–032.
  59. Tomaszkiewicz, M.; Abou Najm, M.; El-Fadel, M. Development of a groundwater quality index for seawater intrusion in coastal aquifers. Environ. Model. Softw. 2014, 57, 13–26. [Google Scholar] [CrossRef]
  60. Appelo, C.A.J.; Postma, D. (Eds.) Geochemistry, Groundwater and Pollution, 2nd ed.; CRC Press: London, UK, 2005; ISBN 978-0-429-15232-0. [Google Scholar]
  61. Kreitler, C.W. Geochemical Techniques for Identifying Sources of Ground-Water Salinization; CRC Press: Boca Raton, FL, USA, 1993; ISBN 978-1-56670-000-9. [Google Scholar]
  62. Panteleit, B.; Hamer, K.; Kringel, R.; Kessels, W.; Schulz, H.D. Geochemical processes in the saltwater–freshwater transition zone: Comparing results of a sand tank experiment with field data. Environ. Earth Sci 2011, 62, 77–91. [Google Scholar] [CrossRef]
  63. Sharaky, A.; Atta, S.; El Hassanein, A.; Khallaf, K. Hydrogeochemistry of Groundwater in the Western Nile Delta Aquifers, Egypt; Cairo University: Giza, Egypt, 2007; Volume 19, pp. 1–23. [Google Scholar]
  64. Winslow, A.G.; Kister, L.R. Saline-Water Resources of Texas; U.S. Government Printing Office: Washington, DC, USA, 1956.
  65. Abdel Baki, A. Hydrogeological and Hydrochemical Studies on the Area West of Rosetta Branch and South El Nasr Canal. Ph.D. Thesis, Ain Shams University, Cairo, Egypt, 1983. [Google Scholar]
  66. Piper, A. A graphic procedure in the geochemical interpretation of water-analyses. Eos Trans. Am. Geophys. Union 1944, 25, 914–928. [Google Scholar] [CrossRef]
  67. Sharaky, A.M.; El Hasanein, A.S.; Atta, S.A.; Khallaf, K.M. Nile and Groundwater Interaction in the Western Nile Delta, Egypt. In The Nile Delta; Negm, A.M., Ed.; The Handbook of Environmental Chemistry; Springer International Publishing: Cham, Switzerland, 2016; Volume 55, pp. 33–62. ISBN 978-3-319-56122-6. [Google Scholar]
  68. Berkowitz, B.; Singurindy, O.; Lowell, R.P. Mixing-driven diagenesis and mineral deposition: CaCO3 precipitation in salt water–fresh water mixing zones: Mixing-Driven Diagenesis. Geophys. Res. Lett. 2003, 30. [Google Scholar] [CrossRef]
  69. Drevaliene, G.; Karro, E.; Mokrik, R.; Savitskaja, L. The origin of barium in the Cambrian–Vendian aquifer system, North Estonia. Est. J. Earth Sci. 2009, 58, 193–208. [Google Scholar] [CrossRef]
  70. Rawson, J.; Siade, A.; Sun, J.; Neidhardt, H.; Berg, M.; Prommer, H. Quantifying Reactive Transport Processes Governing Arsenic Mobility after Injection of Reactive Organic Carbon into a Bengal Delta Aquifer. Environ. Sci. Technol. 2017, 51, 8471–8480. [Google Scholar] [CrossRef] [PubMed]
  71. Biswas, A.; Gustafsson, J.P.; Neidhardt, H.; Halder, D.; Kundu, A.K.; Chatterjee, D.; Berner, Z.; Bhattacharya, P. Role of competing ions in the mobilization of arsenic in groundwater of Bengal Basin: Insight from surface complexation modeling. Water Res. 2014, 55, 30–39. [Google Scholar] [CrossRef]
  72. Hosono, T.; Hossain, S.; Shimada, J. Hydrobiogeochemical evolution along the regional groundwater flow systems in volcanic aquifers in Kumamoto, Japan. Environ. Earth Sci. 2020, 79, 410. [Google Scholar] [CrossRef]
  73. Marandi, A.; Shand, P. Groundwater chemistry and the Gibbs Diagram. Appl. Geochem. 2018, 97, 209–212. [Google Scholar] [CrossRef]
  74. Tiwari, A.K.; De Maio, M.; Singh, P.K.; Singh, A.K. Hydrogeochemical characterization and groundwater quality assessment in a coal mining area, India. Arab. J. Geosci. 2016, 9, 177. [Google Scholar] [CrossRef]
  75. Bouderbala, A.; Gharbi, B.Y. Hydrogeochemical characterization and groundwater quality assessment in the intensive agricultural zone of the Upper Cheliff plain, Algeria. Environ. Earth Sci. 2017, 76, 744. [Google Scholar] [CrossRef]
  76. Shin, K.; Koh, D.-C.; Jung, H.; Lee, J. The Hydrogeochemical Characteristics of Groundwater Subjected to Seawater Intrusion in the Archipelago, Korea. Water 2020, 12, 1542. [Google Scholar] [CrossRef]
  77. Hussien, R.A. Groundwater Quality Index Studies for Seawater Intrusion in Coastal Aquifer Ras Sudr, Egypt Using Geographic Information System. Sciences 2015, 5, 209–222. [Google Scholar]
Figure 1. Location map showing the extent of the study area and the distribution of the collected water samples.
Figure 1. Location map showing the extent of the study area and the distribution of the collected water samples.
Water 15 01803 g001
Figure 2. Spatial distribution of TDS across the study area.
Figure 2. Spatial distribution of TDS across the study area.
Water 15 01803 g002
Figure 3. The spatial distribution maps of major ions.
Figure 3. The spatial distribution maps of major ions.
Water 15 01803 g003
Figure 4. The spatial distribution maps of minor ions.
Figure 4. The spatial distribution maps of minor ions.
Water 15 01803 g004
Figure 5. (a) Dendrogram shows the grouping of the groundwater samples into four groups depending on the major ion concentrations in milliequivalents per lite and (b) Scholler diagrams characterizing the chemical type for samples groups.
Figure 5. (a) Dendrogram shows the grouping of the groundwater samples into four groups depending on the major ion concentrations in milliequivalents per lite and (b) Scholler diagrams characterizing the chemical type for samples groups.
Water 15 01803 g005
Figure 6. Spatial distribution map of the groundwater samples groups.
Figure 6. Spatial distribution map of the groundwater samples groups.
Water 15 01803 g006
Figure 7. Piper diagrams show the hydrogeochemical characters of the studied groundwater. Samples are labeled according to their HCA groups shown in Figure 5a.
Figure 7. Piper diagrams show the hydrogeochemical characters of the studied groundwater. Samples are labeled according to their HCA groups shown in Figure 5a.
Water 15 01803 g007
Figure 8. Dendrograms show the ion-associations in sample groups.
Figure 8. Dendrograms show the ion-associations in sample groups.
Water 15 01803 g008
Figure 9. Gibbs diagrams illustrate (a) TDS concentrations are plotted against the weight ratios of Na/(Na + Ca) for cations and (b) the weight ratios of Cl/(Cl + HCO3) for anions. Concentrations of TDS and ion are in mg/L.
Figure 9. Gibbs diagrams illustrate (a) TDS concentrations are plotted against the weight ratios of Na/(Na + Ca) for cations and (b) the weight ratios of Cl/(Cl + HCO3) for anions. Concentrations of TDS and ion are in mg/L.
Water 15 01803 g009
Figure 10. The relationships between Total Dissolved Solids (TDS) and the concentrations of (a) Ca, (b), Mg, (c) Na + Ca, (d) Cl, (e) HCO3, and (f) SO4.
Figure 10. The relationships between Total Dissolved Solids (TDS) and the concentrations of (a) Ca, (b), Mg, (c) Na + Ca, (d) Cl, (e) HCO3, and (f) SO4.
Water 15 01803 g010
Figure 11. GQISWI in the study area.
Figure 11. GQISWI in the study area.
Water 15 01803 g011
Table 1. The analyzing techniques of the measured parameters.
Table 1. The analyzing techniques of the measured parameters.
ParameterUsed InstrumentAnalyzing TechniquesAnalyzing Place
Electric conductance
Ec (mS/cm)
Hach’s Portable conductivity/ Total dissolved solids (TDS)The typical analytical methods explained by HACH [50]Field
Total dissolved solids
TDS (mg/L)
Temperature (°C)
Hydrogen ion activity (pH)Portable Consort pH Meter (Model p 314)
SO4 and NO3Hach’s Direct Reading (DR/2000) SpectrophotometerDepartment of Geology, Faculty of Science, Tanta University
Cl and HCO3Hach’s Digital Titrator Model 16900-01
Na, K, Ca, Mg, Al, Ba, Fe, Mn, Sr, and SiICP (Inductive Coupled Plasma Optima 7000 DV)EPA method 200.7 by USEPA [51]Center of Scientific Research and Measurements, Tanta University
Table 2. Statistical analysis of the obtained hydrochemical data.
Table 2. Statistical analysis of the obtained hydrochemical data.
ParameterMinMaxAvSD
TDS (mg/L)19027,6803806.27001.4
EC (ms/cm)0.3855.37.112.3
pH5.058.77.10.8
TH (mg/L)25.933315762.8839.4
Ca (mg/L)4.74670153.6159.5
Mg (mg/L)3.4343092.3110.9
Na (mg/L)37.9892001156.72339.8
K (mg/L)2.34337.734.868.4
HCO3 (mg/L)13.02394222.674.1
SO4 (mg/L)02300365.7517.7
Cl (mg/L)4814,3001844.13647.6
NO3 (mg/L)020222.936.6
Al (mg/L)9 × 10−40.360.10.1
Ba (mg/L)12 × 10−515.390.31.7
Fe (mg/L)12 × 10−534.561.14.1
Mn (mg/L)9 × 10−42.8530.50.5
Si (mg/L)4.5723.4212.13.9
Sr (mg/L)0.0462.54.611.1
Table 3. Classification of water according to Total Dissolved Solids [64].
Table 3. Classification of water according to Total Dissolved Solids [64].
Water ClassSalt Concentration (mg/L)No. of Wells%
Fresh>10004053
Slightly saline1000–30001621
Moderately saline3000–10,0001115
Very saline10,000–35,000811
Brine<35,000--
Table 4. The average composition of the hydrochemical parameters of the four groundwater sample groups.
Table 4. The average composition of the hydrochemical parameters of the four groundwater sample groups.
ParameterG1G2G3G4
TDS mg/L22,686.756486.4597.941829.1
Ca meq/L26.0914.274.245.37
Mg meq/L30.1214.972.575.94
Na meq/L329.7882.245.5119.17
K meq/L5.180.720.170.56
HCO3 meq/L3.421873.493.95
SO4 meq/L29.7113.942.925.93
Cl meq/L332.5992.925.620.91
Si mg/L12.6312.8911.3110.59
Al mg/L0.090.030.010.00
Ba mg/L0.123.110.130.24
Fe mg/L0.2510.010.551.47
Mn mg/L0.491.550.490.88
Sr mg/L0.8028.011.976.14
Table 5. Dissolved minerals saturation indices.
Table 5. Dissolved minerals saturation indices.
No.AragoniteCalciteDolomiteWitheriteAnhydriteGypsumBariteHalite
G1Min−1.7−1.6−2.7−8.4−1.3−0.9−2.5−3.1
Max0.20.31.0−3.6−0.40.02.5−2.2
Av−0.7−0.6−0.8−5.5−0.7−0.4−0.1−2.7
SI > 1%12.5252500057.20
SI < 1%87.5757510010010042.8100
G2Min−1.9−1.7−3.1−8.6−2.2−1.9−2.7−4.2
Max1.11.22.5−3.7−0.8−0.6−0.7−3.8
Av−0.5−0.4−0.5−6.0−1.4−1.1−1.5−3.9
SI > 1%40404000000
SI < 1%606060100100100100100
G3Min−1.1−1.0−1.9−6.4−3.8−3.5−4.2−7.2
Max1.21.32.6−1.8−1.1−0.70.7−5.4
Av0.40.50.9−2.9−2.1−1.8−0.2−6.3
SI > 1%84.9087.9087.900063.70
SI < 1%15.1012.1012.110010010036.3100
G4Min−1.9−1.7−3.4−8.0−3.8−3.5−3.3−7.1
Max0.80.91.7−1.8−1.2−0.81.2−4.4
Av−0.8−0.6−1.1−5.2−2.1−1.8−1.2−5.4
SI > 1%22.625.922.600025.90
SI < 1%77.474.177.410010010074.1100
No.SylviteHematitePyrolusiteAlbiteAnorthiteK-feldsparKaoliniteQuartz
G1Min−4.6−2.7−1.8−4.4−12.1−4.3−1.80.1
Max−3.716.52.90.1−4.30.64.80.7
Av−4.010.60.7−2.2−7.7−1.61.30.4
SI > 1%087.502502550100
SI < 1%10012.51007510075500
G2Min−5.71.7−0.4−2.9−8.4−2.91.10.1
Max−4.920.30.4−1.7−5.1−1.11.70.5
Av−5.49.70.1−2.5−6.9−2.21.40.3
SI > 1%01000000100100
SI < 1%100010010010010000
G3Min−7.95.5−0.2−3.8−6.6−2.60.70.1
Max−6.420.32.10.1−0.71.35.80.5
Av−7.115.01.1−1.0−2.20.13.60.3
SI > 1%010005.8079.5100100
SI < 1%100010094.210020.500
G4Min−7.61.1−0.6−4.6−9.6−3.80.5−0.1
Max−4.921.22.7−0.6−1.60.06.40.5
Av−6.411.40.5−2.9−6.5−2.01.70.2
SI > 1%0100000010088
SI < 1%1000100100100100012
Table 6. Correlation matrix between physico-chemical parameters in groundwater.
Table 6. Correlation matrix between physico-chemical parameters in groundwater.
TDSECCaMgNaKHCO3SO4ClNO3AlBaFeMnSiSr
TDS1
EC0.91
Ca0.90.91
Mg0.90.90.91
Na0.90.90.80.91
K0.60.60.50.60.71
HCO30.10.10.20.10.10.11
SO40.60.60.50.70.60.70.01
Cl0.90.90.90.90.90.70.10.61
NO3−0.1−0.1−0.10.0−0.10.10.00.1−0.11
Al−0.2−0.2−0.2−0.3−0.2−0.1−0.1−0.1−0.2−0.11
Ba0.50.50.40.40.50.4−0.10.50.4−0.1−0.11
Fe0.40.40.40.40.50.4−0.10.50.4−0.1−0.10.91
Mn0.10.10.20.1−0.1−0.10.2−0.10.1−0.2−0.10.10.11
Si0.10.10.10.10.10.0−0.10.10.1−0.20.30.10.10.21
Sr0.50.50.60.60.50.6−0.10.60.50.1−0.20.90.90.10.11
Note: values listed in red font have strong positive correlation, blue font have intermediate positive correlation, values listed in green font have weak correlation.
Table 7. GQISWI ranges [77].
Table 7. GQISWI ranges [77].
Water TypeGQISWI Based on Worldwide LiteratureTypical GQISWI
MinMaxMeanMinMax
Fresh water73.590.182.775100
Mixed groundwater47.879.963.45075
Saline groundwater4.858.827.51050
Seawater3.19.25.8010
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hasan, S.S.; Salem, Z.E.; Sefelnasr, A. Assessment of Hydrogeochemical Characteristics and Seawater Intrusion in Coastal Aquifers by Integrating Statistical and Graphical Techniques: Quaternary Aquifer, West Nile Delta, Egypt. Water 2023, 15, 1803. https://doi.org/10.3390/w15101803

AMA Style

Hasan SS, Salem ZE, Sefelnasr A. Assessment of Hydrogeochemical Characteristics and Seawater Intrusion in Coastal Aquifers by Integrating Statistical and Graphical Techniques: Quaternary Aquifer, West Nile Delta, Egypt. Water. 2023; 15(10):1803. https://doi.org/10.3390/w15101803

Chicago/Turabian Style

Hasan, Samia S., Zenhom E. Salem, and Ahmed Sefelnasr. 2023. "Assessment of Hydrogeochemical Characteristics and Seawater Intrusion in Coastal Aquifers by Integrating Statistical and Graphical Techniques: Quaternary Aquifer, West Nile Delta, Egypt" Water 15, no. 10: 1803. https://doi.org/10.3390/w15101803

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop