Geochemistry and Origin of the Hyper Salinity of Groundwater in Wasia-Biyadh Aquifer, Saudi Arabia ()
1. Introduction
Salinity of groundwater is a function of many chemical and physical factors, such as the mineral composition, lithology of the aquifer, oxidation reduction reactions, natural or artificial recharge and discharge, weathering, exchange reactions (Elango & Kannan, 2007). Groundwater is classified generally according to the total dissolved solids (TDS) concentration as fresh (<1000 mg/L), brackish (1000 - 10,000 mg/L), saline (10,000 - 100,000 mg/L), and briny water (>100,000 mg/L) (Freeze & Cherry, 1979). Another classification based on chloride ion classifies groundwater into four categories as fresh (<300 mg/L), brackish (300 to 10,000 mg/L), salt (10,000 to 20,000 mg/L), and hypersaline (>20,000 mg/L) (Stuyfzand & Stuurman, 1994). Engelen et al. (2018) considered groundwater with salinity > 35,000 mg/L as “hypersaline”.
The origin of the hypersaline groundwater in aquifers all over the world is still unclear until today, yet several authors made attempts such as De Sitter (1947), White (1965), and Berry (1969). They introduced the so-called “salt-sieving”, also known as “natural reverse osmosis” which stated that shale and clays are considered as semi-permeable membranes, and are capable of preventing ions during the compaction of seawater-saturated sediments. Such mechanism could be applied to groundwater leakages through aquitards. However, certain conditions must exist to produce such groundwater flow patterns.
Avrahamov et al., 2014 indicated that the microbial evidence pointed to the anaerobic oxidation of methane (AOM) along with the bacterial sulfate reduction in the hypersaline groundwater of the Dead Sea (DS) Aquifer. Engelen et al., 2018 referred the origin of the hypersaline groundwater in the Nile Delta Aquifer (NDA) to the leakage of hypersaline groundwater from deeper deposits, or to the evapotranspiration of seawater. Sabarathinam et al., 2020 referred the origin of hyper salinity in the acidic hypersaline coastal groundwater in Kuwait to the bacterial oxidation of organic matter and the oxidation-reduction cycles of Fe.
In the present paper, we are seeking to explain the origin of the hyper salinity in deep seated groundwater of Wasia-Biyadh Aquifer.
2. Methods
The Ministry of Environment Water and Agriculture (MEWA) conducted the current research in cooperation with the Deutsche Gesellschaft für Internationale Zusammenarbeit (GIZ) as an internal report project, titled “Detailed Groundwater Resources Studies in the Rub Al-Khali Desert. Additional data was compiled from MOWE (2013a, 2013b, 2013c), BRGM (1976), and Al-Mashaikhi (2011). A total of 31 groundwater samples were collected by the GIZ from November 2010 to December 2011 for major cations and anions, trace elements and radionuclides. All Groundwater samples were collected in 100 mℓ bottles, filtered using a 0.45 μm filter and preserved in cooling boxes. Nitric acid (HNO3 70%) has been added additionally to samples assigned for cation analysis. All field measurements involved EC, T, pH dissolved oxygen, hydrogen sulfide and bicarbonate using multi-parameter pocket meters and a flow-through-cell (WTW, Multi 340i). Groundwater samples were also collected for trace elements in 250 ml bottles, then shipped to the laboratory in Darmstadt Technical University, Germany. The following analytical techniques have been applied to identify the chemical composition of the water samples:
- Ion chromatography (IC; major ions)
- Atomic Absorption Spectrometry (AAS; Fe2+)
- Photometric tests (B, SiO2)
- Gravimetry (TDS)
- Atomic Absorption Spectrometry, ICP/MS, and ICP-OES (trace elements)
- Isotope Ratio Infrared Spectrometry (IRIS) was used to for the 2H and 18O analysis
3. Hydrogeology
The Biyadh and the overlying Wasia Aquifer are of lower to mid Cretaceous age, composed of medium to coarse grained sandstone of a non-marine origin with some interbedded shale (MEWA, 2017d). In the western part of the study area (Figure 1), close to the outcrop, the Biyadh and the overlying Wasia Aquifer are in direct contact, forming the Wasia-Biyadh aquifer, which is considered a principal aquifer in Saudi Arabia, whereas in the eastern part they are separated by Shuaiba aquitard (Figure 2). The sandstone dips laterally in the subsurface under Rub’ Al Khali depression, grading into shaly facies that completely replaced by limestone and argillaceous limestone in the easten part of the study area (Kalbus et al., 2011). The Wasia-Biyad Aquifer system is unconformably topped with dolomite and limestone of the Aruma Formation as illustrated in the hydrogeologic cross section (Figure 2) which extends from west to east across Rub Al Khali area. Over most parts of the study area the Wasia-Biyadh aquifer is confined, while the unconfined conditions prevail only in the western part. The investigated area belongs commonly to an arid climate due to the low precipitation and high evaporation rate. The mean annual rainfall range between 100 - 220 mm/y, and the temperature reaches more than 50˚C in summers (Kalbus et al., 2011).
Figure 2. Hydrogeologic cross section for the investigated area from west to east.
4. Results and Discussion
4.1. Major Geochemical Processes
The groundwater samples in the western parts of the aquifer showed a broad range of TDS concentration, ranging between a few hundred to 3000 mg/ℓ, belonging to a Ca-HCO3 water type. A gradual increase in the ion concentrations was observed in groundwater along the flow path. However, some samples showed very high salinity with a TDS concentration of up to 230,000 mg/ℓ and a Na-Cl water type towards the eastern parts of the study area (Figure 3). A Na+-normalized Ca2+ VS
plot suggested by Gaillardet et al. (1999) was employed to identify the relative contributions of the major weathering/dissolution mechanisms (silicate, carbonate, and evaporite) to the ion concentrations in groundwater. Based on Figure 4, most of groundwater samples fell in the area of silicate weathering and evaporation/and or evaporite dissolutions, indicating that the major mechanisms which contribute ions and solutes are silicate mineral weathering and evaporite dissolution and/or evaporation.
4.2. Water-Rock Interaction
The main constituent of the sandstone in the western part of Wasia-Biyadh Aquifer is mainly Quartz (95%) with minor percentage of K-feldspars (1%) and clay (5%) according to the XRD (Figure 5) which has been carried out by the GIZ. The main geochemical reaction in such part of the aquifer is the incongruent dissolution of K-feldspars which results in releasing potassium, bicarbonate and silica into groundwater as indicated by Equation (1).
(1)
In the eastern part of the aquifer, The XRD for a core sample (RAK-02E) indicated that quartz is the main constituent of the sandstone aquifer (Figure 6) in addition to minor percentages of plagioclase, dolomite, pyrite, and clay. The precipitation of pyrite in this portion of the aquifer was driven by the Sulfate-reducing bacteria (SRB). It can perform anaerobic respiration utilizing sulfate as a terminal
Figure 4. Na+-normalized Ca2+ VS
plot, showing that silicate weathering and evaporite dissolution are the main sources of ions.
Figure 5. (a) Core sample RAK 04-E. (b) Pie chart showing the results from the X-Ray Diffraction (XRD) Analysis. (c) Low magnification view of PPL and XPL images showing abundant quartz (Qtz), minor K-feldspar, lithic fragment (Lith) and intergranular porosity (IgΦ). (d) High magnification view highlighting quartz (Qtz), and intergranular pores (IgΦ) and the dissociated K-feldspar (KFld) (MEWA, 2017b).
Figure 6. (a) Core sample of Wasia Aquifer at RAK 02-0 well. (b) Pie chart showing the result of XRD analysis. (c) Low magnifigation of wasia sandstone under the microscope. (d) High magnification of wasia sandstone under microscope showing Quartz (Qtz), trace plagioclase (Plg), intergranular pore (IgΦ) cementing ferroan dolomite (Fe-DolC) and kaolinite (Kao) observed. Minor pyrite (Py) noted (MEWA, 2017b).
electron acceptor (Equation (2)), reducing it to sulfur or hydrogen sulfide. Based on Figure 7, sulfate concentration decreases from west to north east due to the reduction process. The produced sulfur reacts with Fe2+ to precipitate pyrite (Equation (3)). Fe2+ is produced in the groundwater by the reduction of Fe3+ (Equation (4)) and the precipitation of ferroan dolomite (iron dolomite) that was also driven by the anaerobic bacteria in the aquifer (Equation (5)).
(2)
(3)
(4)
(5)
The incongruent dissociation of plagioclase, albite and anorthite introduced clays as secondary minerals to the aquifer in addition to sodium, calcium and silica to groundwater (Equations (6) & (7)). This was confirmed by the stability
Figure 7. The hydrochemical distribution of ions within Wasia-Biyadh aquifer, sulfate concentration decreases from west to north east due to reduction.
diagrams (Figure 8) for some groundwater samples that appeared to lie at/or close to equilibrium with kaolinite, illite, and montmorillonite. This also could be confirmed by applying the geochemical model phreeqc along the flow path in the eastern part of the aquifer to deduce the plausible reactions which account for solutes in the groundwater. The RAK-03E was selected as the initial well, and RAK-02E as the final. The result of the model is presented in Table 1.
(6)
(7)
The hyper salinity of some groundwaters in the investigated aquifer might be attributed to a different factor other than the typical known geochemical processes. The high temperature of some groundwater samples (Table 2) may suggest that the aquifer is subjected to an internal or external source of heat which might led to the concentration of salinity over time. A deep heating in a naturally occurring hydrothermal system might occur to groundwater due to the presence of a nearby magma chamber or high radionuclide content within an aquifer (Figure 9 & Figure 10). The significant content of the natural radionuclides within the investigated aquifer (Table 2) supports that the natural radioactivity is the source of heat within the aquifer. Based on the bivariate diagram of Na/Cl versus chloride relationship (Figure 11), the high-salinity samples (Na–Cl type) plot on the evaporation trend, which might support the proposed hypothesis.
Table 1. Result of the Phreeqc geochemical model along the flow path RAK-03 to RAK-02.
Table 2. Concentrations of the radionuclides and the environmental isotopes for some groundwater samples of Waia-Biyadh Aquifer (MEWA, 2017a).
Figure 8. Stability field of (a) Ca–Al silicate phases, (b) Na–Al silicate phases, and (c) K–Al silicate phases relative to some selected groundwater samples. Most of the groundwater samples appear to lie at/or close to equilibrium with montmorillonite, kaolinite and illite.
Figure 10. A Schematic diagram showing the effect of the accumulation of the radionuclides in Wasia-Biyad aquifer.
Figure 11. Na+/Cl− vs. Cl− relationship, indicating the sources of sodium ions.
5. Conclusion
The groundwater samples in the western parts of the aquifer showed a broad range of TDS concentration, ranging between a few hundred to 3000 mg/ℓ, belonging to a Ca-HCO3 water type. A gradual increase in the TDS concentrations was observed in groundwater along the flow path. However, some samples showed hyper salinity with a TDS concentration of up to 230,000 mg/ℓ and a Na-Cl water type towards the eastern parts of the study area. According to the hydrochemical and the natural radioactivity data, the presence of a naturally occurring hydrothermal system within the aquifer was considered the source of internal heat that led to heating up the groundwater and hence salinity concentrated over time. However, the geochemical reactions contributed salinity to groundwater. The geochemical reactions in the eastern part of the aquifer included the precipitation of pyrite, dolomite, clay, and the incongruent dissolution of plagioclase.