Availability of Surface Water of Wadi Rajil as a Source of Groundwater Artificial Recharge: A Case Study of Eastern Badia /Jordan

Wadi Rajil catchment area is considered as one of the major wadis entering the Azraq Basin from the north. It is ungauged wadi and covers an area of about 3910km 2 . The annual average rainfall on Wadi Rajil catchment area is about 126.6mm. Heavy thunderstorms occur in April and May, causing significant floods covering the area. The flood waters are not utilized, and a small portion infiltrates into the ground, where the great portion of these waters remain over Qaa’ Azraqfew months before evaporation. Due to the absence of the hydrometric stream flow station, no data are available about surface water runoff in Wadi Rajil catchment area. Therefore, the first part of this study calculates the surface water potential of Wadi Rajil to be utilized for groundwater artificial recharge, applying the SCS curvilinear synthetic unit hydrograph method. The synthesis unit hydrograph of Wadi Rajil catchment is characterized by a peak value of 1146 m3/s (4047 cfs) per one inch of rainfall excess. Flood hydrographs for 10,25,50, and 100 years return periods were derived and their peak flow are found to be 10,8,186,412, and 680 m 3 /s, respectively and the corresponding flood volumes are 0.95, 16.53, 36.89, and 61.5 MCM, respectively.Groundwater artificial recharge conditions are suitably prevailing in the most northern and central part of the catchment area, whereas, geological, Hydrogeological, and water quality characteristics of the floodwater encourage artificial replenishment of the exploited aquifer in the study area.


INTRODUCTION
Azraq basin is located in the Northeastern part of Jordan and extends northwards into Syria and southwards into Saudi Arabia.The Azraq Oasis (called locally Sabkhah or Qa'a Azraq) which is located in the central part of the basin is at a distance of about 120 km northeast of Amman.Qa'a Azraq is a relative large mudflat located in the central part of the basin.Two villages are located on the western side of the Qa'a Azraq; these are Azraq Shishan and Azraq Druze.A well field called AWSA was established north of Azraq Druze Springs (northern springs) where about 15-20 million m 3 per year (MCM/a) of water is pumped to the capital Amman for drinking purposes since 1982.Farmers in the area are using around 45 MCM/a.Therefore, the total abstraction from the basin is about 65 MCM/a (El Naqa et al., 2007).
The over-pumping from the shallow groundwater aquifers, the water level dropped dramatically and signs of salinization and depletion are starting to occur.Therefore, it is necessary to determine the location of this interface, since saline water, intrusion has started and is predicted to accelerate in the near future.(Bajjali, and Hadidi, 2005).
A number of major wadis drain radially into the Azraq depression.These wadis are characterized by wadis shallow flow-beds with relatively low slopes.None of the wadis within the basin is being gauged.Wadi Rajilcatchment area is considered as one of the most significant wadi draining the Azraq basin from its northern part.The general shape of Wadi Rajilcatchment area is trapezoidal, with its longer axis in general oriented NE-SW direction.The slope of the area is from NW to SE and the average slope of stream channels varies from 1% to 3% Elevation ranges from 1550 m above mean sea level (a.m.s.1) at Salkhad in Syria to 500 m a.m.s.1 at Qaa' Azraq.Considerable part of the drainage is nearly flat.
The main objective of this study is to estimate the peak discharge and the flood volumes for different return period in Wadi Rajil catchment area, and the potentiality of using this water for groundwater artificial recharge.

Description of the study area
Wadi Rajil catchment area lies between 313 to 400 E and 138 to 228 N According to Palestine Grid, and it covers 3910 km 2  , Figure 1.Intermittent flow in the wadi occurs in winter and drains into Qaa' Azraq, where it evaporates within a few months without utilization.
The climate in the Azraq Basin belongs to the Mediterranean bio-climatic region, which is characterized by dry hot summers and wet, cold winter.Most of the study area is arid and small portion can be considered as semi-arid.As in most arid and semi-arid areas, temperatures exhibit large seasonal and diurnal variations with absolute daily temperatures ranging from a maximum of about 46 ºC (January of 1993), (JMD, 2010).
The thunderstorm rainfalls from the great part of the total rainfall in the study area, which is characterized by irregular intensity and duration.The heaviest rainfall over 24 hours is usually recorded between December and March and no significant normal rainfall is to be expected in October and May.The average annual rainfall over the sub-catchment area is 126.6 mm and ranges between less than 50mm in the south to more than 300 mm in the northern part of the sub-catchment.

Geology
Azraq Basin in general and Wadi Rajil catchment area in particular is a part of the limestone plateau in east Jordan.The northern part of this plateau is covered by six basalt flows, tuff and volcanic eruption originated mainly in Jabel Arab in Syria and many other relatively small volcanoes in the area occurred during the Miocene-Quaternary.These flows overlie Tertiary rocks in the northern part of the catchment area.
The youngest sediments in the Qaa' consists of gravels and sands of fluviatile origin (NRA, 1992).The outcropping formations are, the Rijam and wadi Shallala in the central and eastern part of the basin, both of Eocene age.
The northern part of the Azraq basin is dominated by Miocene to Pleistocene basalt whereas to the west and south by Rijam and Muwaqqar formations of late cretaceous-Early Tertiary age.Regional distribution of the outcropping geological formations is shown in Figure 2.
The area as the whole of Jordan was effected by the transgressions and regressions that occurred during the time from Cenomanian until upper Eocene (NRA, 1992).
These transgressions and regressions are represented in the Azraq Basin by the changes of facies in the Ajlun Group (limestone, dolomite, marl, shale and dolomitic limestone) and Belqa Group (marl,limestone ,chert, silicified, limestone, dolomite and chalk),that crop out in the central, western and southern part of the basin.During the period from upper Oligocene and into the lower Miocene, the Azraq Basin has gone under erosion and tectonic activity, (Bender, 1974).
In the Azraq playa (wetland reserve), the basalt is missing.Upper Tertiary sediments (B5) (Wadi Shallala) are located in the structural depression zones.The (B5) Formation consists of Marly Clayey layers in the area of AWSA well field and acts here as an aquitard between the B4 (Rijam) and the Basalt aquifer.Towards the southeast, the B5 Formation contains more sandy layers and it is classified as an aquifer in this area.South of the basalt areas Paleocene and Eocene,marly limestone, chalks, and chalky limestone with chert layers of the B4 formation, dominate the landscape.The B4 formation is underlain by the Maastrichtian B3 (Muwaqqar) formation.B3 formation reaches a thickness of about 300 m and consists of marl and marly limestone with some gypsum and evaporite.The underlying Campanian to Turonian B2/A7 formations (Amman/Wadi Sir) is mainly formed by chert and limestone, (El Naqa, 2010).

Methodolgy
There are several methods for determining the peak discharge in un-gauged watersheds.In this study, the SCS-curvilinear synthetic unit hydrograph concept was applied to estimate the peak discharge for the un-gauged Wadi Rajil catchment area.

Derivation of the Unit Hydrograph
The most common method of deriving unit hydrograph synthetically is the soil conversation services (SCS) curvilinear unit hydrograph method, Chow, et al, 1988, Wanielista, 1990).The time to peak or the time required to reach the peak discharge (Tp) is defined as: Where, Tp=is time to peak in hours Dr= is the standard duration of rainfall La= is the lag time The lag time (La) is calculated from the following formula; Where, C t , is aregional constant representing watershed slope and storage (Subramanya, 1984, Chow et al, 1988, Linsely, et al, 1988), L, is the hydraulic length, L c ,is the centroid length and S1 , is the slope.
The standard duration of rainfall (Dr) is calculated as: Dr=La/5.5 ..
The peak discharge (Qp) can be calculated using the SCS formula:

Time Increment Rainfall
T(1) P(n) T( 2) Where, A, is the catchment area in (mi 2) and Qp in ft 3 /sec (cfs).

Critical Arrangement of the storm
The incremental values of the storm of certain return period (e.g.25-year), were arranged in descending order and re-arranged as shown in Table 1, below: From this arrangement, it is indicated that T (n) is the last time increment in 24 hrs andP (1) is the first and consequently the maximum incremental rainfall.This is the critical arrangement of the storm, which produces maximum peak discharge in its flood, (Matthai, 1969).

Derivation of the Flood Hydrographs
The flood hydrograph of each frequency storm is derived by combining of the hydrograph obtained for each incremental net rainfall or runoff (r1, r2,……., ri).In this combination procedure, the shifting of the hydrographs according to its time increment should be considered, (Chow et al., 1988).

Estimation of peak discharge
In this study, the unit hydrograph approach is applied to determine the peak discharge values.The hydrologic characteristic of the drainage area such as the area of the basin(A),hydraulic length (L), centroid length (Lc) and the elevation deference between the highest point of the main stream and the outlet(H) are calculated from the topographic maps related to the catchment.Whereas the Curve Number (CN) is calculated from the topographic maps, geologic maps and land use map available for Azraq Basin.The calculation of the unit hydrograph (UH) and the derivation of the flood hydrographs of 10, 25,50and 100years return period for Wadi Rajil catchment area were performed.

Unit Hydrograph of Wadi Rajil Catchment area
The calculation of UH of Wadi Rajil Catchment areawas done using the English unit, then the obtained peak discharge(Qp)values were converted to the metric units.The parameters taken from the topographic map were: The La value is corrected for Dr = 1 hr.instead of 3.34hr, using Snyder's Formula: Where, is the corrected lag time and is one hour duration.Then, theTp was calculated using formula (1): Finally, the estimated Qp was obtained using Equation ( 4 2, shows the calculation of Wadi Rajil UH.The UH of Wadi Rajil is also illustrated in Figure 3.

Derivation of Wadi Rajil Flood Hydrographs
The intensity duration frequency curves (IDF) of Azraq rainfall station and the CN method mentioned previously, were utilized to calculate the effective rainfall (runoff).The IDF curves are shown in Figure 4.From these curves ,the 10 ,25,50 and 100 years return period rainfalls for duration of 24 hours were obtained hourly, the results of calculation are tabulated in Tables (3,4,5,and 6).These tables also show the critical arrangement of the incremental rainfall to obtain the maximum flood of 10,25,50 and 100 years return period storm.The incremental runoff values were applied to the UH with the time lag, and the individual hydrographs were obtained for each incremental runoff.The addition of these hydrographs give the total storms hydrographs .The application of the procedure to 10,25,50,and 100 years storms for Wadi Rajil catchment are shown in Table 7and the resulted flood hydrographs are illustrated in Figure 5.The calculated peak discharges were 10.8,186,412and 680 m 3 /sfor the 10, 25, 50, and 100 year return period respectively.The calculated flood volumes were 0.95,16.5,36.9and61.5 million cubic meter (MCM)for the 10,25,50,and 100 year return period respectively.

Possibilities of Artificial Recharge
As noticed from the previous section, the study area in particular as well as desert area in Jordan in general receive a low precipitation amount, where most of rainfall occurs in few storms of high intensity and short duration.Flush floods occur in desert areas, due to the low infiltration index.Most of the floodwater evaporates and very little amount for the time being has been utilized.Due to this, the ideal solution is the better management of such water by using it in the replenishment of the over pumped aquifer by artificial recharge techniques.
Artificial recharge is essential on enhancement of natural recharge ; therefore an understanding of the hydrologic cycle precedes developments of sites selection criteria .Five facets of the hydrologic cycle that need to be well understood (Chidly,1981).These are precipitation,  The water quality of surface runoff of Wadi Rajil is considered to be of good quality,Table 8.whereas, the electrical conductivity (EC µS/cm) is considered low and highly appropriate since the water flows mostly on the basalt floes or the weathering products of the basalts or the volcanic ashes in most of the drainage areaand in general will not affect negatively the groundwater quality in the area after its percolation through the highly fractured basalt.
Talking into consideration the hydrological conditions prevailing in the Upper Aquifer Complex in the area, there is an enough space for ground water aquifer since it is all over the basin is over pumped as indicated from the groundwater withdrawals in the area of study, Figure 6.In addition, the hydraulic characteristics of the aquifer as indicated from the pump test analyses of the groundwater wells operating in the areaare as follows: the discharge, specifics capacities and transmissivities are ranging between 70-300m 3 / hr,0.2-35200m 3 /hr/m and 15-72300m 2 /d, respectively, (El-Naqa and Al-Shayeb 2008).Therefore, the aquifer has a good potential to accept recharged water.Furthermore, the depth of the water table in the area is ranging between 20-70m from the ground surface.
As indicated from the above discussion on the situation of Wadi Rajil catchment, artificial recharge is considered of high potential to replenish the Upper aquifer.Concerning the technique of artificial recharge, which should be developed in the study area.It is suggested to use deep trench technique, Figure 7, whereas the course of Wadi Rajil itself as well as its tributaries are of wide type and under these circumstances check or desert dams are not feasible (Rimawi et al., 1995).The proposed technique is usually preferred by constructing successive trenches cross the wadi course.Trenches are filled with boulder, gravels and sands (volcanic ashes) in sequence from bottom to top.Using these techniques, the suspended mater will be removed and good water quality will reach the aquifer and no negative impacts will be expected in the groundwaterquality.

CONCLUSIONS
Generally, the northern part of the Azraq Basin receives higher rainfall amounts than the other parts of the basin and accordingly most of the groundwater resources are originated from north of the Azraq basin.The investigation of the availability of surface water and the prevailing conditions for groundwater artificial recharge for Wadi Rajil catchment area has led to the following results: ´The thunderstorm rainfalls from the great part of the total rainfall in the study area, which is characterized by irregular intensity and duration.The heaviest rainfall over 24 hours is usually recorded between December and March.´The calculated peak discharges for the studied catchment area were 10.8,186,412and 680m 3 /s for the 10, 25, 50, and 100 return period, respectively.´The estimated flood volumes from the resulted flood hydrographs of Wadi Rajil catchment area were 0.95, 16.5, 36.9 and 61.5 MCM for the 10, 25, 50, and 100-year return periods, respectively.´Geological, hydrogeological and waterquality offloodwater are found to be highly suitable for artificial recharge of the Upper Aquifer System.´Deep long narrow trenches refilled with the available basalt boulders, gravels as well as volcanic ashes will be the appropriate method to be used for artificial recharge in Wadi Rajil catchment area .
): = 40471.6or 1146 m 3 / sec (cms) The T/Tp and Q/Qp values of the generalized dimensionless UH of the SCS were used to derive the synthetic UH of Wadi Rajil catchment area.Table