Geofluids
Volume 2018 (2018), Article ID 7593430, 14 pages
https://doi.org/10.1155/2018/7593430

Research Article

Modern Recharge in a Transboundary Groundwater Basin Deduced from Hydrochemical and Isotopic Investigations: Al Buraimi, Oman

O. Abdalla,1,2 A. Izady,¹ T. Al-Hosni,² M. Chen,¹ H. Al-Mamari,¹ and K. Semhi²

¹Water Research Center, Sultan Qaboos University, Muscat, Oman
²Department of Earth Sciences, College of Science, Sultan Qaboos University, Muscat, Oman

Correspondence should be addressed to O. Abdalla; osman@squ.edu.om

Received 10 January 2018; Revised 5 June 2018; Accepted 30 July 2018; Published 12 September 2018

Academic Editor: Karsten Kroeger

Copyright © 2018 O. Abdalla et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Groundwater samples (54) collected from different geological units (alluvium, Tertiary, ophiolite, and Hawasina) located in the transboundary groundwater basin in north Oman at the United Arab Emirates (UAE) borders were analyzed for general hydrochemistry and water isotopes, and subsets thereof were analyzed for ¹⁴C and ³H and ⁸⁷Sr/⁸⁶Sr. The chemical composition, percentage of modern carbon (pmc), δ²H, δ¹⁸O, and ⁸⁷Sr/⁸⁶Sr of the groundwater in the study area progressively change from the recharge zone in the elevated area of the North Oman Mountains (NOM) to the flat plains at the UAE borders. While the water-rock interaction is the dominant process controlling the groundwater chemistry, evaporation and groundwater mixing affect the hydrochemistry at the UAE borders. Therefore, groundwater evolves from carbonate-dominant in the NOM into sodium chloride-dominant close to the UAE borders. It is also evident that groundwater lateral recharge from the ophiolites into the alluvium retains the chemical affinity of the ophiolites. Groundwater dating (high pmc), homogeneous ⁸⁷Sr/⁸⁶Sr ratios, and enriched δ²H and δ¹⁸O demonstrate the presence of modern recharge in the shallow zones of the ophiolites and alluvium. However, deep zones and areas at the UAE border contain older groundwater form during cooler and wetter climatic conditions as supported by the depleted δ²H and δ¹⁸O and lower ⁸⁷Sr/⁸⁶Sr ratios and pmc. Furthermore, the data clearly showed that modern groundwater mixes with older groundwater along the flow path from the NOM into the UAE border. Modern recharge occurs as lateral recharge from NOM and direct recharge in the plain area. The current findings support future development of aflaj system along NOM slopes and shallow wells in the plain areas.

1. Introduction

Transboundary groundwater basins, which cross countries’ borders, are normally of regional scale and thus represent important shared water resources. Proper management of these resources is fundamental for sustainability and development [1–3]. Groundwater in the border regions separating the Sultanate of Oman and the United Arab Emirates (UAE), near the Al Buraimi and Al Ain areas, represents an important resource for sustainable agricultural and urban development [4]. While North Oman Mountains (NOMs), located in the eastern part of this area, have been delineated as the principal recharge source [5], there is no prior hydrochemical and isotopic analyses in this region validating such argument. The hydrochemical study is an important component to understand the availability and nature of groundwater through identifying moisture sources and different geochemical processes that control the quality of water [6]. Hydrochemistry investigation may identify natural processes such as dissociations of salts in water [7, 8], redox processes [9], water-rock interaction, and mixing processes that strongly influence groundwater chemistry. Besides, environmental and some radioactive isotopes have become an important tool in identifying different patterns of precipitation, recognizing mixing processes, tracing the flow path, and dating groundwater. The current paper utilizes hydrochemistry and isotopes to understand recharge process and identify the processes controlling the groundwater chemistry in the transboundary aquifer of Al Buraimi-Al Ain areas.

2. Study Area and Hydrogeological Characteristics

The study area lies between 24°22 N to 24°38 N latitude and 55°44 E to 56°14 E longitude and covers around 1600 km² in the Al Buraimi northwest of Oman. It is bounded by border with United Arab Emirates (UAE) on the west and by the ridge of North Oman Mountain (NOM) on the east (Figure 1). The study area with an annual average rainfall of 82 mm (from 30 mm to 178 mm) is characterized by very low humidity and precipitation. Most of the rainfall is between December and April. In summer months, the daytime temperature may exceed 45°C. The long-term annual potential evapotranspiration is about 2700 mm [10–12].

Figure 1: Study area, geology, and sample location.

The geologic units of the study area are divided into three principal zones that include the Semail nappes (ophiolite), the Hawasina nappes, and the post-nappe strata units (Aruma Group and Tertiary bedrock and Quaternary alluvium) (Figure 1). The term “ophiolite” refers collectively to igneous rock that crops out in the study area with various dark, colored, crystalline, and microcrystalline characteristics. The Hawasina nappe exposures mainly occur as broken hills in the eastern piedmont zone. The post-nappe strata consist of the Cretaceous Aruma Group and Tertiary bedrock that were deposited in a foredeep basin downfolded along the frontal margin of the nappes. Folding associated with mountain building in the Late Tertiary turned over the nappes and Tertiary strata into their present structural configurations. Afterwards, erosive processes associated with flowing water led to the deposition of alluvium throughout the piedmont and alluvial fan zones to the west of the mountains. In terms of regional hydrogeological significance, surface waters running off the mountainous portion of the study area recharge the piedmont and alluvial fan zones. Groundwater recharged into the fractures of the ophiolite during rainfall and runoff events gradually drain laterally and vertically downhill through the fractures towards the plains of the lower basin. Groundwater flow from the ophiolite is generally to the west, and it is intercepted by aflaj and a few wells at the mountain front and continues to the alluvial aquifers underlying the upper basin. Compared to other geological units, the hydrological properties of the Hawasina nappes and Tertiary limestones are poor in terms of groundwater storage and transmission. However, they may provide significant groundwater supply through fractures. The alluvium in the piedmont and mountain front fan zones is the most important aquifer in the study area and is composed of an unconsolidated mixture of chert, limestone, and dolomite. The alluvium generally has good hydraulic conductivity (8 to 43 m/d) and water yields (0.1 to 31 l/s), which vary only with the coarseness of the alluvium and overall saturated thickness [10, 11, 13–17].

3. Material and Methods

Several field trips were carried out to the study area between the years 2014 and 2015 that span different seasons. A total of 54 locations (Figure 1) were sampled during these field trips including hand-dug wells, boreholes, monitoring wells, precipitation, and aflaj (refers to an open channel that drains water from the slope of the mountain and transports water via gravity into demand areas). All samples’ locations were recorded using a global positioning system (GPS). Different parameters including temperature, electrical conductivity (EC), pH, and total dissolved solid (TDS) were measured immediately in the field using Aquaread Aquameter with multiple electrodes. Sampling procedure varies according to the location. For instance, the boreholes’ samples were collected directly from the pumping wells after a few hours of continuous pumping. The hand-dug wells were bailed using a Solinst stainless steel bailer after being purged with the aid of a purging pump. The monitoring wells were similarly sampled after purging while samples from aflaj and precipitation were collected directly. The hand-dug wells tap the shallow water system, and they comprise private wells mostly for irrigation purposes. The boreholes however tap deeper horizons and screened different intervals which may constitute inevitable groundwater mixing. Three samples (for cations, anions, and isotopes) at each location were collected in 500 ml plastic bottles for analyses. The bottles collected for cation analysis were acidified using 2-3 drops of diluted (5%) nitric acid. The bottles were transported in iced container and preserved in a refrigerator at a temperature of 5°C. All the samples were immediately analyzed for the major ions and alkalinity in the Central Analytical and Applied Research Unit in Sultan Qaboos University (CAARU-SQU) using ion chromatography (IC) and inductively coupled plasma mass spectrometry (ICP-MS). The stable isotopes deuterium (²H) and oxygen-18 (¹⁸O) were analyzed in the Hydrogeology Laboratory Unit in Sultan Qaboos University (SQU) using Triple Isotopic Water Analyzer (TIWA-DLT-EP) from Los Gatos Research (LGR). Selected samples were analyzed in the Environmental Isotopes Laboratory of the University of Waterloo for carbon isotopes (13C and 14C), tritium (3H) (11 samples), and strontium isotopes (⁸⁷Sr and ⁸⁸Sr) (17 samples).

4. Results and Discussion

4.1. Chemical Composition and Groundwater Classification

Chemical analyses of the collected samples from different aquifers are tabulated in Table 1. The hydrochemical characterization of the groundwater in the study area, aided by the Piper diagram (Figure 2), reveals three main groups: Mg-HCO₃, Mg-Cl-SO₄, and Na-Cl dominant groups. Except of one sample that may suggest a deeper source, all of the aflaj samples belong to the Mg-HCO₃, circulate within the ophiolite rocks, and are located along the slope of the NOM. The dominance of the carbonates indicates exchange with the atmospheric CO₂ and circulation of modern water at shallower depths. The Mg on the other hand is the weathering product of the ophiolites at shallower depth where the pH (around or less than 9) favors the dissolution of brucite which is a magnesium hydroxide. In addition, about half of the alluvium groundwater also falls within the Mg-HCO₃. This group of alluvium groundwater represents recent recharge from direct infiltration through streambeds. The dominance of Mg in these samples is attributed to the presence of ophiolitic fragments in the alluvium.

Table 1: Field parameters, physical characteristics, and isotopic composition of groundwater samples.

No.
Name
Depth to water (m)
Lithology
EC (μS/cm)
TDS (mg/l)
pH
δ ¹⁸O (per mil)
δ ²H (per mil)

1
B-9
9.8
Alluvium
1610
1046
7.27
−2.207
−5.498

2
WZDB1
34.12
Tertiary
1403
911
7.88
−1.075
−1.135

3
A-5
29.92
Hawasina
1059
686
7.99
−2.133
−5.747

4
WSKD1
23.81
Ophiolite
1234
818
7.88
−1.484
−2.472

5
ZG-3
27.14
Ophiolite
7017
4560
7.52
−1.024
−2.540

6
S-5
37.7
Hawasina
841
546
8.11
−1.359
−2.013

7
BU-EH11
33.25
Hawasina
988
641
8.85
−1.931
−5.241

8
S-28
46.75
Alluvium
708
460
8.4
−1.457
−2.374

9
S-6
43.7
Alluvium
1033
672
8.12
−1.484
−1.410

10
BU-EH 06
41.5
Alluvium
522
340
8.37
−1.247
−1.549

11
Davoodi falaj
—
Ophiolite
628
408
8.4
−1.397
−3.240

12
Mahdah falaj (Aini)
—
Ophiolite
599
390
8.39
−0.260
3.540

13
Al Khadra falaj
—
Ophiolite
748
487
8.87
−0.318
3.691

14
Wadi Sharm falaj
—
Ophiolite
841
546
8.23
−0.360
3.267

15
P-18
32.4
Alluvium
660
429
8.61
−0.017
7.739

16
Juwaif falaj
—
Ophiolite
714
464
8.45
−0.596
3.172

17
Juwaif well
>100
Ophiolite
2770
1809
11.8
−1.033
−0.061

18
Khatwa falaj
—
Ophiolite
624
406
8.6
0.918
7.159

19
Mehasen farm
>100
Ophiolite
605
395
7.91
−1.933
−4.985

20
S-25
25.78
Alluvium
749
486
8.65
−0.819
−0.531

21
P-8
35
Alluvium
639
415
8.81
−1.436
−3.279

22
WM-4
46.85
Alluvium
1555
1011
7.77
−0.933
−2.631

23
WM-3
46.75
Alluvium
1368
889
7.7
−0.965
−2.673

24
WM-6
47.45
Alluvium
1458
946
8.06
−1.060
−2.979

25
EW-4
25.35
Alluvium
662
430
8.23
−0.745
−1.774

26
EW-14
36.7
Alluvium
838
545
8.25
−1.773
−6.526

27
WS-15
38.4
Alluvium
889
577
8.24
−1.858
−6.498

28
P-6
43.51
Alluvium
668
433
8.72
−0.731
0.282

29
EW-10
41.48
Alluvium
595
387
8.61
−1.291
−2.398

30
E-5
26.35
Alluvium
2630
1710
7.72
−1.805
−3.182

31
EW-5
28.12
Alluvium
704
456
8.56
−1.504
0.425

32
WS-12
22.15
Alluvium
1368
889
8.13
−1.638
−5.307

33
WS-11
23.15
Alluvium
999
649
8.01
−1.543
−4.635

34
E-4
28.9
Alluvium
905
588
8.39
−1.608
−4.477

35
P-17
27.42
Alluvium
617
401
7.69
−1.745
−4.684

36
EW-3
19.98
Alluvium
1077
702
8.2
−1.289
−3.635

37
E-1
22.4
Alluvium
734
476
8.25
−1.551
−4.243

38
EW-2
25.5
Alluvium
978
634
8.19
−2.415
−9.344

39
EW-9
36.12
Alluvium
1147
745
8.02
−2.917
−11.299

40
P-7
41.45
Alluvium
899
584
8.61
−1.242
−1.073

41
P-1
42.7
Alluvium
4497
15917
6.95
−1.805
−8.589

42
EW-21
23.95
Alluvium
1567
1015
8.81
−1.416
−5.261

43
MF-1
43.85
Alluvium
1662
1082
7.96
−2.119
−8.415

44
S-15
42.55
Alluvium
1068
685
8.2
−2.192
−8.193

45
A-17
34.42
Alluvium
1477
962
8.25
−1.573
−4.563

46
SMD-1
32
Alluvium
15470
10054
7.33
−0.664
−2.159

47
SMD-5
30.95
Alluvium
1540
1000
7.94
−1.408
−3.162

48
AJ3-AJ3A
32.13
Tertiary
1973
1283
10.02
−2.206
−10.849

49
Farm 1
33.8
Alluvium
1797
1172
8.23
−1.465
−6.783

50
Farm 2 well 1
Deep
Tertiary
13570
8900
7.72
−2.212
−14.327

51
Farm 2 well 2
Shallow
Alluvium
3219
2063
8.11
−1.407
−4.925

52
Farm 3
Deep
Ophiolite
4760
3094
7.74
−0.919
−0.088

53
Farm 4
Deep
Ophiolite
2775
1803
9.48
−0.226
−0.371

54
Farm 5
—
Hawasina
2087
1343
8.16
−0.951
0.525

Figure 2: Piper diagram illustrating groundwater classification in different geological formations.

Contrary to the Mg-HCO₃ group, the Mg-Cl-SO₄ and Na-Cl groups may represent older groundwater (carbon-14 dating is discussed later). The deep boreholes penetrating the ophiolites (sample numbers 52 and 53) and the Tertiary (sample number 50) have shown the Mg-Cl-SO₄ and Na-Cl dominance, respectively. The hydrochemical dominance in ophiolite groundwater changes with depth from Mg-HCO₃ to Mg-Cl-SO₄ (Figure 2). While Cl is the dominant anion in both ophiolites and Tertiary deep groundwater, Na is the dominant cation in the former and Mg is the dominant cation in the latter. This cationic variation is attributed to the mineralogical composition of the hosting rock. Nonetheless, a couple of odd samples has shown a reversed dominance likely attributed to cationic exchange. Such cation variation supports the role of water-rock interaction and suggests that time circulation of the Mg-Cl-SO₄ and Na-Cl groundwater is larger compared to the Mg-HCO₃ groundwater type. The other half of the alluvium groundwater samples, mostly from relative deeper depths, belong to the Mg-Cl-SO₄ and Na-Cl groups, and similar to the deep borehole samples, their chemistry indicates no modern recharge, dissolution of evaporites and silicate minerals. Therefore, groundwater recharge processes seem to influence the hydrochemistry in the study area. Recharge takes place by two mechanisms and from two zones: higher elevated mountains and wadi systems such as Wadi Safwan. Recharge in the elevated high mountains is a direct recharge from rainfall into fractures, whereas recharge in the plain area of wadi systems is a direct percolation through streambed into the aquifer. Along its flow path from the NOM into the plains before the UAE borders, the groundwater hydrochemistry changes in the study area from carbonate-dominant into chloride/sulphate-dominant. The dissolution of halite, gypsum, and anhydrite present in the Hawasina and the Tertiary rocks may explain the dominance of the chloride and sulphate. Moreover, additional atmospheric source of Cl in such an arid area is possible where Cl from aerosols first accumulates into the soil profile and then leaches to the groundwater [18].

In general, the groundwater chemistry evolves with increasing depth and proximity to the NOM as indicated in the Piper diagram (Figure 2) signifying the role of water-rock interaction. Such role has been recognized in several arid areas across the world, for instance, in Chile [19], UAE [20] adjacent to the study area, and in Tunisia [21]. Hydrochemical investigations of groundwater in Oman (e.g., [22–28]) have made a consensus that water-rock interaction is the primary factor controlling the groundwater chemistry with a secondary contribution from evaporation. The NOM marks the region of fresh water and groundwater salinity increases away from NOM with increasing mineral dissolution. However, due to geological diversity of NOM, the anion/cation suit of groundwater varies.

Although nitrate is reported in most of the samples, it exceeds the limits for the Omani Standards (1998) in two samples number 49 and number 54 which were collected from farm numbers 1 and 5 (Tables 1 and 2). Thus, nitrate elevated concentration is due to agricultural activities and animal manure. The magnesium vs. calcium plot (Figure 3) shows that the Mg concentration exceeds the Ca concentration except for two samples collected from the alluvium aquifer. The highest ratio is observed in groundwater samples collected from the ophiolite aquifer and those collected from wells penetrating both alluvium and ophiolite (Mg/Ca ratio higher than 3). While the lowest ratio (less than 1) is observed in the groundwater samples collected from only the alluvium aquifer, the magnesium is likely to be driven from ophiolite weathering as well as from the dolomite in the Tertiary sediments. The limestone dissolution in the Tertiary sediments seems less effective; as in general, groundwater affected by limestone dissolution will reflect an average Mg/Ca ratio of 1 : 2 to 1 : 1.5 [29]. The excess of Ca + Mg over HCO₃ (Figure 4) indicates an excess dissolution of silicates over the carbonates, namely, the dolomite. If the main source was carbonate mineral dissolution, the Ca + Mg : HCO₃ should have shown a near to 1 value. The excess of the Ca + Mg indicates a source lacking HCO₃, and the available lithology will suggest the silicates found in the ophiolites.

Table 2: The concentrations of different elements.

No.
Name
CO₃ (mg/l)
HCO₃ (mg/l)
F (mg/l)
Cl (mg/l)
NO₃ (mg/l)
SO₄ (mg/l)
Na (mg/l)
K (mg/l)
Ca (mg/l)
Mg (mg/l)
Sr (mg/l)
⁸⁷Sr/⁸⁶Sr

1
B-9
0
236
0.22
340.02
3.22
197.35
198.28
14.96
68.26
70.57
1.37
—

2
WZDB1
18
174
0.20
91.04
4.12
80.21
95.66
5.18
11.60
36.65
0.74
—

3
A-5
4
112
0.04
103.03
19.60
43.84
47.58
2.65
22.42
38.99
0.79
—

4
WSKD1
0
133
0.10
171.04
10.09
81.75
87.66
4.78
31.92
44.32
1.34
0.7086

5
ZG-3
4
101
0.04
32.02
6.45
30.36
20.41
1.93
14.89
24.79
0.50
—

6
S-5
0
212
0.06
148.64
39.74
94.66
74.28
4.62
42.16
66.07
1.54
0.7086

7
BU-EH11
6
171
0.28
38.42
13.81
59.35
70.34
3.65
23.57
21.15
0.87
—

8
S-28
0
172.2
0.00
90.91
12.73
91.01
59.14
3.25
17.83
30.88
0.12
—

9
S-6
0
115.62
0.12
346.68
1.07
2440.28
415.17
9.47
318.02
272.60
0.92
—

10
BU-EH 06
0
107
0.07
69.93
1.26
34.56
45.79
4.72
21.48
16.95
0.50
0.7086

11
Davoodi falaj
19.68
118.08
0.00
32.59
6.71
29.91
33.99
2.16
17.89
18.82
0.06
0.7086

12
Mahdah falaj (Aini)
0
114
0.02
45.21
5.94
51.17
15.28
1.15
15.07
42.84
0.65
0.7086

13
Al Khadra falaj
10
132
0.03
30.97
9.18
35.14
12.79
1.39
20.22
38.24
0.48
—

14
Wadi Sharm falaj
0
157.44
0.00
55.63
9.83
31.69
22.88
1.90
19.02
39.12
0.05
—

15
P-18
0
125.46
0.00
52.49
0.00
0.77
20.58
1.99
9.02
23.96
0.01
0.7087

16
Juwaif falaj
4.92
130.38
0.00
26.42
5.00
46.03
10.14
1.01
11.26
41.34
0.06
—

17
Juwaif well
13.2
39.6
0.00
260.87
0.00
0.52
134.26
7.06
4.11
11.26
0.05
0.7087

18
Khatwa falaj
0
22.14
0.00
26.56
7.53
36.25
16.08
1.40
16.91
31.61
0.04
—

19
Mehasen farm
4.92
137.76
0.00
17.81
5.93
19.35
45.94
1.68
7.73
10.35
0.02
0.7085

20
S-25
6
98
0.04
102.41
5.05
29.71
33.59
2.46
10.43
44.53
0.60
—

21
P-8
10
84
0.05
75.98
0.38
18.98
48.85
3.91
11.22
19.23
0.51
—

22
WM-4
9.84
137.76
0.00
185.55
20.97
82.05
91.54
4.94
37.66
43.61
0.22
—

23
WM-3
0
164.82
0.00
179.98
15.74
88.35
81.10
4.43
39.68
46.34
0.24
0.7085

24
WM-6
14.76
110.7
0.00
158.61
15.02
71.63
82.93
4.48
32.22
42.83
0.22
—

25
EW-4
6
111
0.21
47.33
7.10
37.04
34.70
2.96
21.58
22.58
1.17
0.7086

26
EW-14
0
125.46
0.00
104.28
3.39
44.42
56.97
3.07
22.27
30.49
0.14
—

27
WS-15
0
54.12
0.00
48.07
4.04
40.21
38.98
2.76
17.19
23.17
0.13
—

28
P-6
0
54.12
0.00
68.47
0.00
0.86
39.52
3.14
6.43
6.70
0.02
0.7086

29
EW-10
10
97
0.03
49.59
1.04
28.09
24.54
2.17
15.28
28.48
1.71
—

30
E-5
0
72.6
0.20
272.07
29.41
157.15
197.40
9.90
91.23
21.46
0.62
—

31
EW-5
0
97
0.56
82.64
6.49
61.94
58.20
4.84
29.19
37.95
0.26
—

32
WS-12
9.84
179.58
1.00
106.28
0.15
141.72
133.73
3.81
34.16
10.62
0.21
—

33
WS-11
0
167
0.56
82.13
6.45
61.63
57.88
4.83
29.67
37.23
1.43
—

34
E-4
0
130
0.18
82.81
5.52
52.76
75.57
3.91
17.69
17.77
1.23
0.7082

35
P-17
0
174.66
0.00
37.59
5.64
11.93
33.77
2.94
21.42
21.67
0.08
—

36
EW-3
0
120
0.08
62.04
3.30
37.65
37.64
2.99
21.25
27.09
0.98
—

37
E-1
0
116
0.07
69.35
4.35
35.71
38.33
3.66
22.68
27.59
1.20
—

38
EW-2
0
122
0.11
62.68
4.24
38.90
42.43
2.72
22.16
24.98
1.19
0.7086

39
EW-9
0
136
0.12
127.36
4.40
96.33
89.16
4.16
29.39
33.51
2.17
—

40
P-7
0
120
0.03
40.10
0.17
15.98
23.17
2.25
14.80
20.87
0.48
—

41
P-1
0
29.52
0.00
627.72
0.00
162.16
414.65
6.02
102.11
91.18
0.23
—

42
EW-21
0
206.64
0.26
135.07
0.37
123.37
154.55
4.54
31.15
28.77
0.11
—

43
MF-1
0
167.28
0.20
187.38
5.37
80.54
107.67
4.83
33.51
37.97
0.13
0.7085

44
S-15
0
105.78
0.00
93.20
43.52
31.80
37.72
2.51
24.25
30.94
0.08
—

45
A-17
0
123
0.00
81.13
41.14
39.05
30.26
2.58
25.37
39.73
0.09
0.7086

46
SMD-1
0
236.16
3.42
1594.06
1.68
2766.70
808.35
26.03
618.91
643.27
1.07
—

47
SMD-5
0
201.72
0.24
115.60
6.17
90.98
80.24
4.02
35.51
31.55
0.08
—

48
AJ3-AJ3A
9.84
132.84
0.22
239.57
0.00
2.61
162.08
9.75
14.47
7.06
0.07
—

49
Farm 1
19.68
147.6
0.00
124.98
107.40
85.26
116.30
4.42
27.23
35.12
0.13
—

50
Farm 2 well 1
0
226.32
2.13
1874.77
13.66
1041.92
1307.70
15.46
161.37
174.41
0.98
0.7078

51
Farm 2 well 2
14.76
167.28
1.22
404.08
1.58
220.81
326.71
7.37
61.01
43.51
0.53
—

52
Farm 3
9.84
169.74
0.00
475.62
20.61
424.11
200.74
7.12
151.63
168.98
0.47
—

53
Farm 4
0
34.44
0.00
233.25
0.52
48.64
84.08
6.69
21.51
98.42
0.19
—

54
Farm 5
0
73.8
0.00
133.92
239.38
33.33
43.82
3.30
34.11
73.04
0.13
0.7086

Figure 3: Mg vs. Ca ratio. The solid line indicates 1 : 1 ratio.

Figure 4: Variations of Ca⁺² + Mg⁺² vs. HCO₃⁻ + SO₄⁻² concentrations.

The plot of different samples in the Gibbs diagram (Figure 5), a plot of TDS versus weight ratio of Na/(Na + Ca), suggests that the chemical composition of most of the groundwater samples collected from the different aquifers is controlled by weathering rather than evaporation except for a few samples that show evaporation dominance. Evaporation can be atmospheric which takes place during precipitation or may occur during runoff as direct evaporation of surface water. The atmospheric evaporation effect is reflected in the isotopic signature of groundwater (will be discussed in the following section) and is seen in the highlands. The surface water evaporation is observed in the plain areas, particularly in the Safwan area. The surface water evaporation in the highlands is unlikely due to the steeper slope and thus the residence time of water on the surface is short which allows no direct surface water evaporation. Thus, only the effect of atmospheric evaporation could be seen. However, in the plain areas where the runoff along the wadis get ample time on the surface, direct evaporation is significant, and therefore, the TDS of the groundwater progressively increases. Further discussion of the evaporation effects is elaborated in the next section.

Figure 5: Gibbs diagram compares the role of evaporation and water-rock interaction processes on groundwater.

4.2. Isotopic Characteristics

4.2.1. Stable Isotopes: δ²H and δ¹⁸O

The stable isotopes δ²H and δ¹⁸O data analyzed during this study (Table 1) are presented in Figure 6 with reference to the global meteoric water line (GMWL) after Craig [30]. The δ¹⁸O ranges from −2.917 (EW-9, sample number 39) to +0.918 (Khatwa falaj, sample no. 18) while that of δ²H ranges from −4.327 (farm 2 well 1, sample number 50) to +7.739 (P-18, sample number 15). The regression line for all groundwater and aflaj data yields a slope of 5.3 indicating an evaporation effect. The most isotopic enriched group of groundwater samples is the aflaj discharging from the ophiolite aquifer (shown as orange diamonds in the Figure 6). This group of groundwater samples is located at the slopes of the highlands, and they receive direct recharge from precipitation and thus reflects the effect of atmospheric evaporation. Secondary isotopic enrichment due to evaporation from aflaj channels is minimum as samples were collected directly from the mother well (first water discharge point). The short time circulation of groundwater within this shallow system does not allow for mixing with deeper older water that is likely depleted in its isotopic signature. Comparing the isotopic composition of the fossil stalagmites deposited during the early to middle Holocene with the modern stalagmites in Oman has revealed that the fossil stalagmites are much more isotopically depleted [31]. Reconstructed groundwater values for these time intervals indicate that groundwater was predominately recharged by monsoon precipitation from the Indian Ocean remarkably characterized by depleted isotopic signature. Hence, the most isotopically depleted samples in the study area, which are the Tertiary and a few alluvium samples close to the UAE border, are likely to form during the early to middle Holocene time when the climatic conditions were different from todays. The groundwater sample collected from farm 2 well 1 (sample number 50) penetrating the Tertiary rocks shows the lowest percentage of modern carbon (pmc) indicating older ages (pmc is discussed in the next section). In addition, two groundwater samples from the alluvium located closer to the UAE border (E4, sample number 34, and MF1, sample number 43) show low pmc. The general trend of the isotopic composition when different aquifers are considered is that it gets more depleted in the direction of the groundwater flow from the elevated mountains of the ophiolites in the east to the plain areas of the alluvium in the west crossing the piedmonts of the Hawasina and Tertiary. This phenomenon can be attributed to mixing between the current isotopically heavy current precipitation and depleted older groundwater already present in the system. Therefore, the groundwater becomes progressively depleted as one goes away from the recharge zone at the NOM (Figure 7). It is also suggested that climatic conditions in the early to middle Holocene were cooler and wetter than todays and the evaporation was less.

Figure 6: δ ²H and δ¹⁸O are plotted with reference to the global meteoric water line (GMWL).

Figure 7: Contour lines of the oxygen-18 isotope. The isotopic concentration of the groundwater is progressively depleted in the direction of the flow.

4.2.2. Strontium Isotopes

Strontium in groundwater is mainly released from hosting rocks because of the water-rock interaction. Aquifers with high amount of potassium-bearing minerals contain more Rb and therefore more ⁸⁷Sr isotopes than carbonate aquifers. However, the concentration of strontium may change due to evaporation or mixing of different water’s sources. Furthermore, the radioactivity of the strontium increases with increasing groundwater residence time in the aquifer. However, the ratio of ⁸⁷Sr/⁸⁶Sr does not change with latitude or altitude. The correlation between the ratio of ⁸⁷Sr/⁸⁶Sr and 1/Sr is a useful tool to identify the source rocks. Two mixing groundwater sources will show a strong positive correlation, whereas a poor correlation indicates that there are more than two groundwater sources [27].

The Sr concentration in the groundwater samples (17 samples) collected during this study varies from 0.02 to 2.17 mg/l, whereas the ⁸⁷Sr/⁸⁶Sr ratio remains unchanged except for 2 samples (E4, sample number 34, and farm 2 well 1, sample number 50). For the 15 samples, the ⁸⁷Sr/⁸⁶Sr ratio is either 0.7085, 0.7086, or 0.7087 while a slightly lower ratio is reported for E4 (0.7082) and farm 2 well 1 (0.7078). These two samples are those which have shown extra depleted stable isotopes and low pmc, thus suggesting older groundwater form under different climatic conditions. The homogeneity of the ⁸⁷Sr/⁸⁶Sr ratio, except for the two samples E4 (sample number 34) and farm 2 well 1 (sample number 50), suggests a unified Sr isotopes source, likely the ophiolites. Such source can either be the main ophiolite rock or the fragments present in the alluvium. The ⁸⁷Sr/⁸⁶Sr homogeneity of the groundwater in all aquifers also suggests that these aquifers are hydraulically connected regardless of their lithology and hydraulic property variations.

Contrary to the homogenous ⁸⁷Sr/⁸⁶Sr isotopic ratio, the Sr concentration shows considerable variation. The highest Sr concentration is observed in the sample from well EW-9 (sample number 39), whereas the lowest concentration characterizes the groundwater from Mehasen farm (Table 2). The concentration of Sr in the groundwater samples collected from the ophiolite aquifer in the study area ranges between 0.02 and 0.65 mg/l with an ⁸⁷Sr/⁸⁶Sr ratio about 0.7085–0.7087 (Figure 8). For samples collected from the Tertiary aquifer, the Sr concentration varies from 0.06 to 1.54 mg/l. In groundwater samples from ophiolite aquifer, the ⁸⁷Sr/⁸⁶Sr ratio is higher than the values reported by Weyhenmeyer et al. [32]. This is attributed to the absence of the carbonate rocks of the Hajar Super Group (HSG) in the study area, which according to Semhi et al. [27] and Weyhenmeyer et al. [32] marks the main source of recharge. In the study area, only the ophiolites are present, which when interact with the groundwater will release more Sr as opposed to the HSG. There is a little variation of the ⁸⁷Sr/⁸⁶Sr ratio, for all aquifers as opposed to Sr concentration, which exhibit a large variation. Large variation of Sr concentration may indicate the effect of groundwater residence time. This is well demonstrated by reasonable correlation between the Sr concentration and δ¹⁸O ‰ (Figure 9) of groundwater. Similar to δ¹⁸O ‰ that gets depleted with the flow direction, the Sr concentration increases which signifies the increasing groundwater residence time.

Figure 8: ⁸⁷Sr/⁸⁶Sr ratio vs. 1/Sr for groundwater from ophiolite, Tertiary, and alluvium aquifers.

Figure 9: Correlation between the Sr concentration (mg/l) and δ18O ‰. The diagram shows a fair correlation, which indicates the change due to groundwater increasing residence time.

4.2.3. Carbon 14 and Tritium

In parallel to tracing, isotopes can also be used to estimate the age of groundwater. During infiltration of groundwater through the soil, carbon 14 is introduced to the groundwater from CO₂ produced in the soil and from the atmosphere. The dissolved inorganic carbon decays with time. The term percentage of modern carbon “pmc” is used to express the content of modern carbon in the sample. The higher pmc, the younger will be the groundwater sample. The highest age, which can be determined using this method, cannot exceed 40000 years. If the water is older than 40000-year, the result will be almost 0% pmc [33]. Data of carbon 14 in the samples of the present study show that most of the groundwater in the ophiolite have high pmc (Table 3). The water collected from the ophiolite aquifer have a depth less than 32 m and the area itself is elevated and sloped. One sample shows 1% modern carbon, which is “farm 2 well 1” sample located close to the UAE borders. The well “farm 2 well 1” penetrates the Tertiary, and its actual depth is reported as exceeding 100 m (the real depth is unknown). Relatively lower pmc values are also reported in 2 groundwater samples collected from the alluvium: E-4 (sample number 34) and MF1 (sample number 43) (42.51 and 54.55, respectively). Both samples are also located close to the UAE boarder, i.e., the lower reaches of the groundwater flow path. Also, sample number 54 (farm 5) shows a 64 pmc which may also suggest a relatively older groundwater.

Table 3: Tritium and percentage of modern carbon results for selected groundwater samples.

No.
Sample
³H
±1
¹⁴C

±8.0 T.U.
pmc

6
S-5
<6.0
3.0
92.67

11
Davoodi falaj
<6.0
3.1
83.14

12
Mahdah falaj (Aini)
<6.0
3.1
75.25

15
P-18
<6.0
3.1
91.64

19
Mehasen farm
<6.0
3.0
86.09

28
P-6
<6.0
4.0
83.17

34
E-4
<6.0
4.1
42.51

43
MF-1
<6.0
4.1
54.55

48
AJ3-AJ3A
<6.0
3.1
84.91

50
Farm 2 well 1
<6.0
4.1
1.517

54
Farm 5
<6.0
4.1
64.00

Notes: ³H is tritium, T.U. is tritium unit, and PDB and pmc are, respectively, Pee Dee Belemnite and percentage of modern carbon.

Tritium was also analyzed during this study to complement carbon-14 dating. The results of the tritium analysis during this study do not tell much. The tritium results (direct method that has a precision of 6 T.U.) indicate for all the analyzed samples a concentration less than 6 T.U. (Table 3). In the view of the banned nuclear activities, tritium is no longer introduced to the groundwater systems as during the 1950s, and therefore, its efficacy as an indirect measurement of groundwater age is questioned [34].

5. Discussion Summary

In the view of the findings described in the above sections, a hydrogeological/hydrogeochemical cross section is shown in Figure 10 about the hydrochemistry of the area. The cross section illustrates the evolution of the groundwater chemistry from the recharge zone in NOM dominated by ophiolites to the discharge zone in the plain area at the UAE borders dominated by alluvium cropping at the surface. Induced recharge in the shallow zones of the ophiolite and alluvium increases Mg and HCO₃ concentration, enriches δ²H and δ¹⁸O, provides the source for the strontium isotopes, and introduces fresh recent groundwater (high pmc) to the system. Flowing away from the NOM, the groundwater hydrochemistry gradually changes to dissolve more Cl, Na, and SO₄ ions, gets depleted in δ²H and δ¹⁸O, maintains the same strontium isotopes but increases Sr concentration, and decreases pmc which means becomes older. Similar lateral changes in groundwater are also observed vertically as deeper groundwater contains more Cl, Na, and SO₄ ions, is depleted in δ²H and δ¹⁸O, and has low pmc. The current qualitative recharge assessment in the study area constrains the estimated groundwater recharge in the same region conducted by Izady et al. [4]. The present study added the hydrochemical aspect and identified processes controlling water chemistry. The findings of the three studies encourage further developments in the region along the slope of the NOM and in the shallow zone of the plain area. However, this development should not exceed the estimated sustainable recharge values. Izady et al. [4] suggested 18.09 Mm³/year as long-term regional groundwater recharge and 5.67 Mm³/year as long-term lateral groundwater flux from the ophiolites in NOM into the alluvial zone. Exceeding these values will cause hydraulic imbalance that will encourage upconing of deeper saline groundwater. Conducting a proper solute transport modeling is recommended to identify possible zones of upward groundwater flow to prevent salinization.

Figure 10: Hydrogeological/hydrogeochemical schematic cross section summarizes the groundwater chemistry and the prevailing physiochemical processes (see Figure 1 for the location).

6. Conclusion

The hydrochemical characterization and isotopic investigation of groundwater in the study area have revealed different water types (Mg-HCO₃, Mg-Cl-SO₄, and Na-Cl). The variation in the chemical composition of the groundwater is attributed to recharge processes, and the mineralogical composition of the hosting rocks as the mineral dissolution was found to be the main factor affecting the chemistry of the groundwater. Lateral groundwater recharge from the ophiolite induces an increase of CO₂ from atmospheric circulation along with Mg concentration from the dissolution of silicates. Increasing Cl, SO₄, and Na in the downstream is due to lack of recharge and the dissolution of evaporites. Radiogenic isotopes of carbon-14 indicate modern groundwater in the shallow zones of the ophiolites and the alluvium and old water in the deeper horizons. Such modern recharge and its mechanism are indicated from the stable isotopes of deuterium and O-18 that shows evaporation effects on the shallow groundwater and a completely different signature in the deep groundwater formed under different climatic conditions. The analysis of strontium isotopes reveals considerable homogeneity in the groundwater, which demonstrates a single source of Sr isotopes and hydraulic connectivity of the aquifers. The future development of groundwater resources in the investigated area can be planned in the view of the current findings. Aflaj systems along the slope of the NOM and the downstream shallow alluvium aquifer can provide additional water supply. Such future exploitation must be supported by confirmed presence of modern groundwater evident from application of recent techniques (e.g., helium/tritium).

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Acknowledgments

We would like to thank the Sultan Qaboos University for providing financial support under the project # CL/SQU-UAEU/14/04.

References

C. Bourne, “The International Law Commission’s draft articles on the law of international watercourses: principles and planned measure,” Colorado Journal of International Environmental Law and Policy, vol. 3, pp. 65–92, 1992.

M. E. Campana, “Foreword: transboundary ground water,” Groundwater, vol. 43, no. 5, pp. 646–650, 2005.

A. T. Wolf, S. B. Yoffe, and M. Giordano, “International waters: identifying basins at risk,” Water Policy, vol. 5, no. 1, pp. 29–60, 2003.

A. Izady, O. Abdalla, A. Joodavi, and M. Chen, “Groundwater modeling and sustainability of a transboundary hardrock-alluvium aquifer in North Oman Mountains,” Water, vol. 9, no. 3, p. 161, 2017.

A. Izady, O. Abdalla, A. Joodavi, A. Karimi, M. Chen, and A. F. B. Tompson, “Groundwater recharge estimation in arid hardrock-alluvium aquifers using combined water-table fluctuation and groundwater balance approaches,” Hydrological Processes, vol. 31, no. 19, pp. 3437–3451, 2017.

R. A. Freeze and J. A. Cherry, Groundwater, Prentice-Hall, Inc., 1979.

A. M. Al Dahaan, N. Al-Ansari, and S. Knutsson, “Influence of groundwater hypothetical salts on electrical conductivity total dissolved solids,” Engineering, vol. 8, no. 11, article 72338, 2016.

A. L. Herczeg, S. Dogramaci, and F. W. J. Leaney, “Origin of dissolved salts in a large, semi-arid groundwater system: Murray Basin, Australia,” Marine & Freshwater Research, vol. 52, no. 1, pp. 41–52, 2001.

R. D. Lindberg and D. D. Runnells, “Ground water redox reactions: an analysis of equilibrium state applied to Eh measurements and geochemical modeling,” Science, vol. 225, no. 4665, pp. 925–927, 1984.

M. B. Kaczmarek, M. Brook, T. Haig, R. E. Read, and D. George, Water Resources Assessment and Hydrogeology of the Mahdah Watershed, Northern Oman, Ministry of Water Resources, Sultanate of Oman, 1993.

MRMWR: Ministry of Regional Municipalities and Water Resources, Consultancy for Drilling and Aquifer testing program of Buraimi at Ad Dhairah Region. Ministry of Regional Municipalities, Environment and Water Resources, Sultanate of Oman, 2004.

W. R. Turner, A. S. Forbes, and J. R. Rapp, Groundwater Resources of the Wadi Safwan Area, Public Authority for Water Resources (PAWR), Sultanate of Oman, 1986, Report: PAWR I-86-9.

W. D. Davison, Results of Test Drilling in the Buraimi Area, Sultanate of Oman, Public Authority for Water Resources (PAWR), Water Supply Paper 1, Sultanate of Oman, 1982.

M. B. Kaczmarek, Hydrogeology and Groundwater Availability Buraimi Production Well-Field, Part I – Final Report. Director General Water, Water Projects, Hydrogeology Section, Sultanate of Oman, 1988.

M. B. Kaczmarek, M. Brook, T. Haig, R. E. Read, and D. George, Water Resources Assessment and Hydrogeology of the Zarub Gap Watershed, Wadi Musayliq and Wadi Al Ayn, Northern Oman, Ministry of Water Resources, Sultanate of Oman, 1993.

M. B. Kaczmarek, M. Brook, T. Haig, R. E. Read, and D. George, Water Resources Assessment and Hydrogeology of Wadi Safwan and Wadi Sharm, Northern Oman, Ministry of Water Resources, Sultanate of Oman, 1993.

M. Onanda, V. Price, T. Jolley, A. Stuck, and R. Hall, Water Balance Computation for the Sultanate of Oman. Final Report, Ministry of Regional Municipalities and Water Resources, Sultanate of Oman, 2013.

A. J. Love, A. L. Herczeg, L. Sampson, R. G. Cresswell, and L. K. Fifield, “Sources of chloride and implications for ³⁶Cl dating of old groundwater, Southwestern Great Artesian Basin, Australia,” Water Resources Research, vol. 36, no. 6, pp. 1561–1574, 2000.

C. Zang, J. Dame, and M. Nüsser, “Hydrochemical and environmental isotope analysis of groundwater and surface water in a dry mountain region in Northern Chile,” Environmental Monitoring and Assessment, vol. 190, no. 6, p. 334, 2018.

M. M. Mohamed and A. B. Hassane, “Hydrochemistry assessment of groundwater quality in Al-Ain city, UAE,” Environmental Earth Sciences, vol. 75, no. 4, p. 353, 2016.

J. Hammami Abidi, B. Farhat, A. Ben Mammou, and N. Oueslati, “Characterization of recharge mechanisms and sources of groundwater salinization in Ras Jbel coastal aquifer (Northeast Tunisia) using hydrogeochemical tools, environmental isotopes, GIS, and statistics,” Journal of Chemistry, vol. 2017, Article ID 8610894, 20 pages, 2017.

O. Abdalla, R. Al Abri, K. Semhi, T. Al Hosni, M. Amerjeed, and I. Clark, “Groundwater recharge to ophiolite aquifer in North Oman: constrained by stable isotopes and geochemistry,” Environmental Earth Sciences, vol. 75, no. 15, article 1117, 2016.

O. Abdalla, R. Al Abri, K. Semhi et al., “Hydro-chemical evolution of groundwater in a sequence of Tertiary Formations in Northwest Oman,” Environmental Earth Sciences, vol. 75, no. 21, article 1410, 2016.

O. Abdalla and R. Al-Abri, “Factors affecting groundwater chemistry in regional arid basins of variable lithology: example of Wadi Umairy, Oman,” Arabian Journal of Geosciences, vol. 7, no. 7, pp. 2861–2870, 2014.

B. Askri, “Hydrochemical processes regulating groundwater quality in the coastal plain of Al Musanaah, Sultanate of Oman,” Journal of African Earth Sciences, vol. 106, pp. 87–98, 2015.

J. M. Matter, H. N. Waber, S. Loew, and A. Matter, “Recharge areas and geochemical evolution of groundwater in an alluvial aquifer system in the Sultanate of Oman,” Hydrogeology Journal, vol. 14, no. 1-2, pp. 203–224, 2006.

K. Semhi, O. Abdalla, R. Al-Abri, T. Al-Hosni, and I. D. Clark, “Strontium isotopes as a tool for estimation of groundwater recharge and aquifer connectivity,” Groundwater for Sustainable Development, vol. 4, pp. 1–11, 2017.

C. E. Weyhenmeyer, S. J. Burns, H. N. Waber, P. G. Macumber, and A. Matter, “Isotope study of moisture sources, recharge areas and groundwater flowpaths within the eastern Batinah coastal plain, Sultanate of Oman,” Water Resources Research, vol. 38, article 1184, 2002.

S. Mandel and Z. L. Shiftan, Groundwater Resources: Investigation and Development, Academic Press, New York, NY, USA, 1981.

H. Craig, “Isotopic variations in meteoric waters,” Science, vol. 133, no. 3465, pp. 1702-1703, 1961.

D. Fleitmann, S. J. Burns, U. Neff, A. Mangini, and A. Matter, “Changing moisture sources over the last 330,000 years in Northern Oman from fluid-inclusion evidence in speleothems,” Quaternary Research, vol. 60, no. 2, pp. 223–232, 2003.

C. E. Weyhenmeyer, H. N. Waber, J. Kramers, S. J. Burns, and A. Matter, Strontium Isotope Ratios (87Sr/86Sr) as Tracers for Recharge Areas, Groundwater Movement and Mixing in an Arid Coastal Region of the Sultanate of Oman, American Geophysical Union, Fall Meeting, 2001.

C. A. J. Appelo and D. Postma, Geochemistry, Groundwater and Pollution, A. A. Balkema Publishers, Leiden, The Netherlands, 2nd edition, 2005.

I. Clark, Groundwater Geochemistry and Isotopes, CRC Press, Taylor and Francis group, Boca Raton, FL, USA, 2015.