Generally, the significant sources of contamination in groundwater are anthropogenic. However, groundwater contamination is not exclusively due to the direct input of contaminants. Further it’s related to geogenic contamination as a result of the dissolution of rock-forming minerals enriched with trace elements in the aquifer  . The leaching of metals from these rocks poses a threat to the groundwater resources and may persist for a long time in groundwater  . This type of contamination is a public health issue worldwide with respect to molybdenum (Mo).
The concentration of molybdenum is generally present in small proportions in the natural water system. However, in the study area at Wadi Al Arab (WA) located in northwest Jordan; the carbonate aquifer is affected by geogenic contamination of groundwater by Mo. Whereas, the concentration of Mo reached 1.62 mg/L    that exceeded the WHO Guidelines for Drinking-Water Quality of 0.07 mg/L, and also exceeded the Jordan Institution for Standards and Metrology (JISM) standards of 0.09 mg/L  .
The geogenic potential sources of Mo in the study area are the Cretaceous-Paleogene limestone sediments of the oil shale (B3) and the phosphorites deposits (B2/A7) which contained the highest level of molybdenum  . These rocks are in direct contact with the groundwater resources. The rock-water interaction processes play an important role in the evolution of groundwater .
Based on the previous studies most researchers believed that the source of Mo is geogenic, and refers to naturally occurring elevated concentration of certain elements in groundwater having negative health effects. But the source of molybdenum is not deeply studied; it is from the organic matter in the oil shale deposits or phosphate deposits rather from pyrite or other oxides. Therefore, further detailed research will be conducted to describe the spatial variations in the mineralogical composition of the aquifer matrix lithology and its effect on the hydrochemical characteristics of groundwater resources. The two potential sources of Mo in the study area are the oil shale and the phosphorites deposits, which are in a direct contact with the groundwater resources. Their spatial distribution, dissolution processes and the prevailing water-rock interactions first will be evaluated.
2. Study Area
2.1. Site Description
Wadi Al Arab is located at the northwestern part of Irbid Governorate, about 81 km from the capital city of Amman-Jordan (Figure 1), it is situated at the transition between the mountainous highland in the southeast, the high plateaus and hilly areas in the central part, the rift escarpment, and the Jordan Valley along the western part.
Wadi Al Arab well field was commenced in 1982 with five wells. Currently the well field is composed of 18 operational wells used for the water supply in Irbid. Initially the early drilled five wells were artesian flowing wells and the discharge was ranging between 640 m3/h and more than 1500 m3/h  . The annual abstraction from these wells was 26 Million Cubic Meters (MCM). This high abstraction amount caused the water level to decline significantly and consequently leading all wells to become non-artesian flowing wells. Due to this; the abstraction rate from these wells reduced to approximately 16 - 18 MCM/Year . Since 1982, the measured decline in the water level is 118.64 m/35year. The average annual decline rate of groundwater level is estimated as 3.39 m/year .
The change in groundwater level leads to a dynamic groundwater system characterized by changing aquifer conditions, accompanied by changes in water quality   . As a result, the quality of produced water in the study area shows increased concentrations of heavy metals, in particular Molybdenum, Nickel, and Arsenic, which is in many wells proves a rising trend of contamination.
Figure 1. The hydrogeological map of Wadi Al Arab Aquifer modified after Basem (2000).
2.2. Geology and Hydrogeology
The geological succession in the study area is hydrogeologicaly subdivided into lithostratigraphic units, which form systems of aquifers and aquitards. The aquifer systems in the study area include two major (Kurnub, and Amman Wadi As Sir Aquifer System), and three minor (Na‘ur (A l-2); Hummer (A4); and Rijam (A4))  .
Amman-Wadi As Sir Aquifer (Silicified Limestone and Phosphatic Limestone) (A7/B2) is the main aquifer in the study area. The Muwaqqar Formation (B3) conformably overlies the Amman Formation (B2) in the study area. The B3 is a regional aquitard, which confines the underlying A7/B2 Aquifer system   .
Hydro-geochemical studies were carried out to identify the geochemical processes and their relation to groundwater quality of high Mo concentration. The water quality data collected from Wadi Al Arab boreholes during the period from 1986 until 2018 The collected hydrochemical dataset were tested using different hydrochemical software such as RockWare AQ-QA and Phreeqc Interactive Software in order to calculate the saturation index (SI) value for various related minerals such as calcite, aragonite, and dolomite; in order to evaluate the chemical equilibrium between minerals and water in natural environment. On the other hand, multivariate statistical analysis was conducted using SigmaPlot 14 software. Whereas; the correlation based on the Pearson’s correlation coefficient, principle component analysis (PCA) and the relationship between the different variables of the analyzed hydrochemical parameters. This had been done to evaluate the factors affecting hydrochemical components of the groundwater resources in the study area.
4. Results and Discussion
4.1. Major Cations and Anions
Most of the water samples anions are dominated by with relative abundance of > > Cl− > NO3, whereas; the main cation is Ca2+, with relative abundance of Ca2+ > Mg2+ > Na+ > K+. The trace element content of As, Mo, and Ni exhibits increase concentration in the groundwater of Wadi Al Arab, which is beyond the concentration threshold of the Jordan Institution for Standards and Metrology (JISM), and the World Health Organization (WHO) guidelines.
4.2. Hydro-Chemical Facies of Groundwater
Scatter plots of the cations (Na+, K+, Ca2+, and Mg2+) and anions ( ,Cl−
Table 1. The major ions composition of Wadi Al Arab groundwater samples for the 1986 year.
Table 2. The major ions composition of Wadi Al Arab (WA) groundwater samples for the 2017 year.
and , ) were shown in Piper’s diagram ,classifying the hydrochemical characteristics. One main water type has been identified; all water samples were concentrated on the field of Ca-HCO3, The reason is groundwater passing through carbonate rock, for cation concentration Calcium type and for anion concentration; HCO3-Type of water Predominated (Figure 2). Total hardness of the sample is classified as very hard water, which exceeds 300 mg/L.
The Schöeller diagram (Figure 3) further demonstrates the uniformity and similarity of composition among the groundwater samples by the cluster of the lines of concentration from cations to anions.
Figure 2. Piper diagram classification water type.
Figure 3. Schöeller diagram for different Wadi Al Arab groundwater samples for 2017 and 1986 present the variation of different ion concentrations.
4.3. Correlation Matrix for Analyzed Parameters of Groundwater
The relationship between dependent and independent variables is known as Correlation Matrix, and is statistically measured. It assists in predicting the mutual relationship between two variables. The correlation between the parameters differentiates between strong (the range is +0.8 to 1.0 and −0.8 to −1.0), and weak (the range is +0.0 to 0.5 and −0.0 to −0.5). The correlation values obtained in the present study in Wadi Al Arab are presented respectively in (Table 3 and Table 4). As per the findings of the correlation matrix, it can be inferred that the correlation matrix between parameters is changed during the period between 2017 and 1986, which marks the change of the geochemical processes.
4.4. Hydrochemical Process of Groundwater-Correlation of Major Ions
Groundwater composition is controlled by three main hydrogeochemical processes; including water rock interaction, evaporation and precipitation. Gibbs diagram  is used to give an indication of the main process affecting the water quality. The distribution of the sampling points in the Gibbs diagram showed that the majority of samples are plotted within the field of rock weathering (Figure 4), which indicates that the main factor controlling dissolved hydrochemical components of the water samples is the dissolution of the rock-forming minerals.
Hydrogeochemical data for groundwater in the study area suggested the occurrence of water type as a result of the dissolution of the carbonate minerals, as shown in the following Equation (1) and Equation (2).
Table 3. Pearson correlation matrix for Wadi Al Arab Groundwater Samples in 1986.
Table 4. Pearson correlation matrix for Wadi Al Arab Groundwater Samples in 2017.
Figure 4. Gibbs diagram showing factors controlling the groundwater chemistry .
To further know the functional relationship of the dissolved ions and their contribution to the groundwater the data are plotted in the scatter diagram in Figure 5; Based on the Equations (1) and (2) above the dissolution of calcite and dolomite would produce the and with a molar ratio of (0.5). In this study; the ratios of most samples are above (1) as a result of the dissolution of calcite, while ratios are below (1).
Figure 5. Distributions diagram of ionic ratios in water samples.
Additionally when the Ca2+/Mg2+ molar ratio below (1) the groundwater is undersaturated with respect to Mg2+ which triggers the dissolution of dolomite, most of the plotted samples had the Ca2+/Mg2+ molar ratios largely less than (1) which indicate the dissolution of the aquifer matrix dolomites. Moreover, if the Ca2+/Mg2+ molar ratio greater than (2) can be the result of the dissolution of silicate minerals (silicified limestone and chert).
The relation of Na+/Cl− molar ratio indicates whether the Na+ was derived from the dissolution of halite. The Na+/Cl− molar ratio of most water samples is less than 1 in the analyzed samples collected during 2017. While the 1986 ratio is equal or more than 1. This is attributed to the weathering action of silicate minerals.
The Ca2+ and concentrations of groundwater are controlled by gypsum dissolution and precipitation processes (Equation (3)), the ratio between Ca2+ and showed enriched Ca2+ anddeficient concentration.
The relation between (Ca2+ + Mg2+) versus identifies whether ion exchange or reverse ion exchange is the dominant process based on Equations (4) and (5).
The result of the scatter plot showed that most of the samples falling along the equiline of ratio. Can assume that the Ca2+ and Mg2+ are in equilibrium and signify the result of carbonate weathering.
4.5. Minerals Saturations
The equilibrium state of the water with respect to a mineral phase can be determined by calculating a saturation index (SI) which gives insight of minerals reactivity and can reflect the thermodynamic process of natural water systems, SI values (Equation (6)) of water samples were calculated using Phreeqc Interactive Software.
SI = log (IAP)/K (6)
The calculated results for groundwater samples are shown in (Table 5) Most of the water samples yielded SI calculated results as follows: calcite is supersaturated, dolomite is approximately saturated to supersaturated; which released from the dissolution of carbonate rock (A7/B2) aquifer of Wadi Al Arab. All of the groundwater samples were under-saturated to gypsum and Halite because the major lithological units of the Wadi Al Arab area are insufficient in respect to these minerals.
4.6. Critical Parameter
In order to consider the groundwater safe and portable, its contents must concurrently meet the set standards. Some of these constituents occur in significant amounts with the main ions and others are present as trace elements. The water in the study area shows enrichment of trace elements largely with Mo followed by less concentration amounts of Ni, As, and U. Those categorized as critical parameter relevant from a health prospective based on the JSMO standard (2015) and WHO (2017) and can be harmful to the human health.
The increased concentration of Molybdenum in Wadi Al Arab wells starting from 2008 exceeding the threshold in standards with time (Figure 6), which clearly present the spatial heterogeneity of Mo Groundwater concentration in Wadi Al Arab Aquifer.
The incremental increase of molybdenum concentration in the groundwater
Table 5. Saturation index value for Calcite, Dolomite, Gypsum, and Halite.
with maximum concentration of 1.62 mg/L in well 11a. A correlation of the heavy metal Mo, As, Ni, and U are conducted. In the scattered plot of the Mo with other critical heavy metals (Figure 7), it clearly shows a different trend line. Such concentrations suggest that the Mo was originated and contributed to the groundwater from two different sources which could be the oil shale and the phosphate host rocks.
The wells are divided into two groups based on the correlation behavior of the heavy metal;
· Group 1: includes well name WA11a, WA15, WA5, WA14 which are located in the southern part of the catchment area, at the unconfined aquifer (Figure 7) and have the same behavior.
· Group 2: includes well name WA02, WA3, WA3a, WA13a, which are located in the northern part of the catchment area, at the confined aquifer (Figure 7) and have the same behavior too.
Figure 6. Spatial Variability maps of Mo concentration.
4.7. Redox Condition of Groundwater in the Carbonate Aquifer
Reduction/oxidation (Redox) processes affect the chemical quality of the groundwater in the aquifer system . The redox conditions of groundwater reflect those of their reservoir rocks, even though chemical equilibrium is not usually reached between all minerals present and the groundwater.
The effects of pH and redox potential on Mo mobility in water are reviewed using the paurbaix diagram (Figure 8). As the Mo is highly soluble in water, it is usually present in the form of an oxyanion, and relatively weak in terms of sorption onto geological media with a low organic content .
The paurbaix diagram present (electro-potential/pH (pE/pH) diagram) for Mo. Oxidizing and reducing conditions are described by the y-axis of the diagram (higher positive electrode potential = more oxidizing environment). The acidic and alkaline conditions of the environment are described by x-axis (pH conditions).
The common redox states of Mo in the environment ranges from +4 to +6 (Mo (IV) to Mo (VI)). Under the oxic or nitrate reducing condition the most oxidizing Mo species molybdate ( ) and hydrogen molybdate ( )
Figure 7. The correlation behavior of the heavy metal.
are expected to be the dominant species under neutral pH condition (pH = 6 - 8). Both and are very mobile.
Under iron reducing conditions, Mo(VI) species can be reduced to Mo(V) species ( ). Under anaerobic highly sulfidic conditions, Mo can be further reduced to a solid form as Mo disulfide (or Molybdenite, MoS2), which is considered immobile in the aquifer .
The Mo species in the twelve sampled wells are confined to the field of ; Anoxic/non sulfidic condition (Fe reduction) which is located in the unconfined part of the aquifer, and which is Sulfidic condition (low and high) in the confined aquifer.
The MoS2 is not soluble in water; oxidation of this mineral will result in the production of molybdate oxyanion, which is very mobile in water. The stoichiometry of the reaction can be expressed as in Equation (7).
Figure 8. Paurbaix diagram for molybdenum adapted from Wang (2007).
Similarly, oxidation of other Mo-Bearing sulfides, such us Pyrite, may also result in production of molybdate. Iron or Sulfide oxidizing bacteria have the potential to facilitate metal leaching but in anoxic or anaerobic condition bioleaching reaction is expected to be very slow  .
The highest value of molybdenum measured in Wadi Al Arab is confined to the field of : anoxic/non-sulfidic condition (Fe reduction), which is located in the unconfined part of the aquifer, and which is Sulfidic condition (low and high) in the confined aquifer.
The result of the correlation of heavy metal Mo, As, Ni, and U shows two different trend lines, suggesting that Mo originated and contributed to groundwater from two different sources are mainly the oil shale and the phosphate deposits.
The various concentrations and the different behaviors of molybdenum in the wells located in the northern and southern well fields are still questionable. This could be related to the hydrogeochemical evolution along different separate groundwater flow paths. This evolution represents the chemical weathering and oxidation-reduction reactions, which reflect the degree of Mo mobilizations, as well as mineral co-precipitation/re-adsorption reactions that dictate the degree of Mo immobilization that needs to demonstrate further explanations.
More tests should be conducted to set and to evaluate the complex of Mo behaviors regarding the factors that may affect the transport of Mo under varying redox, pH, and geological media condition, and to assess the fractionation of Mo in the two sources of aquifer matrix.
The authors gratefully acknowledged the financial support provided by the University of Jordan, Deanship of Scientific Research at for the presented work in this paper.
 Pichler, T. and Mozaffari, A. (2015) Occurrence, Distribution and Mobility of Geogenic Molybdenum and Arsenic in A Limestone Aquifer Matrix. Applied Geochemistry, 63, 623-633.
 Pichler, T., Renshaw, C.E. and Sültenfuβ, J. (2016) Geogenic As and Mo Groundwater Contamination Caused by an Abundance of Domestic Supply Wells. Applied Geochemistry, 77, 68-79.
 Amini, M., Abbaspour, K.C., Berg, M., Winkel, L., Hug, S.J., Hoehn, E., Yang, H. and Johnson, C.A. (2008) Statistical Modeling of Global Geogenic Arsenic Contamination in Groundwater. Environmental Science & Technology, 42, 3669-3675.
 Al Kuisi, M., Al-Hwaiti, M. and Mashal, K. (2015) Spatial Distribution Patterns of Molybdenum (Mo) Concentrations in Potable Groundwater in Northern Jordan. Environmental Monitoring and Assessment, 187, 148.
 El-Naqa, A., Ibrahim, K. and Abu Alrous, M. (2018) Hydrogeochemical Investigation of the Origin of Molybdenum (Mo) in Groundwater Resources in Wadi Al Arab Area, Northern Jordan. Acque Sotterranee-Italian Journal of Groundwater, 7, 19-30.
 Rimawi, O. (1985) Hydrogeochemistry and Isotope Hydrology of the Ground and Surface Water in North Jordan (North-Northeast of Mafraq, Dhuleil-Hallabat, Azraq-Basin). PhD Thesis, TU Muenchen, Germany.
 Abed, A.M. and Amireh, B. (1983) Petrography and Geochemistry of Some Jordanian Oilshale from North Jordan. Journal of Petroleum Geology, 5, 261-274.
 Abed, A.M., Sadaqah, R. and Al Kuisi, M. (2008) Uranium and Potentially Toxic Metals During the Mining, Beneficiation and Processing of Phosphorite and Their Effect on Groundwater in Jordan. Mine Water Environment, 27, 171-182.
 Smedley, P.L. and Kinniburgh, D.G. (2017) Molybdenum in Natural Waters: A Review of Occurrence, Distributions and Controls. Applied Geochemistry, 84, 387-432.
 Margana, A. (2006) Hydrogeological Proposal for the Delineation of a Groundwater Protection Area for the Wadi Al Arab Well Field. Technical Cooperation Project WI-BGR, Ministry of Water and Irrigation, Amman.
 Margane, A. (1995) Groundwater Resources of Northern Jordan, Vol.2, Part 2: Monitoring of Groundwater Levels in Northern Jordan. Technical Cooperation Project “Advisory Services to the Water Authority of Jordan. BGR & WAJ, BGR Archive No. 118701, Amman.
 Margane, A., Al-Qadi, M. and Al-Kurdi, O. (2015) Updating the Groundwater Contour Map of the A7/B2 Aquifer in North Jordan. Technical Cooperation Project ‘Syrian Refugee Response’, Technical Report No. 1, BGR & MWI, Amman, 129 p.
 Salameh, E. and Rimawi, O. (1984) Isotope Analyses and Hydrochemistry of the Thermal Springs along the Eastern Side of the Jordan Dead Sea-Wadi Araba Rift Valley. Journal of Hydrology, 73, 129-145.
 Basem, K.M. (2000) The Geology of Irbid and Ash Shuna Ash Shamaliyya (Waqqas) Map Sheets No. 3154-II and 3154-III. Geology Directorate, Natural Resources Authority, Amman, Jordan, Bulletin 46, 9 p.
 Smedley, P.L., Cooper, D.M., Ander, E.L., Milne, C.J. and Lapworth, D.J. (2014) Occurrence of Molybdenum in British Surface Water and Groundwater: Distributions, Controls and Implications for Water Supply. Applied Geochemistry, 40, 144-154.
 Vlek, P.L.G. and Lindsay, W.L. (1977) Thermodynamic Stability and Solubility of Molybdenum Minerals in Soils. Soil Science Society of America Journal, 41, 42-46.
 Erickson, B. and Helz, G. (2000) Molybdenum (VI) Speciation in Sulfidic Waters: Stability and Lability of Thiomolybdates. Geochimica et Cosmochimica Acta, 64, 1149-1158.