Groundwater is one of the most important sources of water in the world, since it represents 30% of the fresh water that consists (fresh water) in about 2.5% of the whole water on the earth   . This importance increases enormously especially for those countries that don’t have surface water. Iraq heavily relies (95%) on its surface water resources   . It was considered as one of the richest countries in its water resources where the allocation per capita was 6029 m3 in 1995 and dropped to 2100 m3 in 2015      . Recently, Iraq is suffering from water shortage problems due to the effect of climate change and building of dams upstream the Tigris and Euphrates Rivers and their tributaries    . With this shortage in surface water in Iraq, it becomes very important to look for another resource of fresh water, which can fulfill the requirements of both the municipal and agricultural practices. The most promising resource of fresh water in middle of Iraq is groundwater. This is due to many reasons among them are its quality and occurrence near the ground surface.
Karbala Governorate is one of the most important Iraq middle cities from its religious point of view. It covers an area of 5034 km2 and has a population that exceeds one million capita  . In this research the seepage velocity (magnitude and direction), it was evaluated for an area within Karbala Governorate, due to the importance of these variables in calculating many properties like the contaminants transporting velocity inside the aquifer, the well discharge rate, permeability of the aquifer, hydraulic conductivity, Darcy velocity and other geohydrology properties.
2. Study Area
The study area is a part of Karbala Governorate and it is located between Karbala and Najaf Governorates (in the middle of Iraq), between 412,000 - 432,000 UTM longitude and 3,574,000 - 3,594,000 UTM latitude. It covers an area of 400 km2. It is surrounded by Khan Al-Nukhaylah from the north, Khan Al-Hammad from the south, oil pipe from the west, and Euphrates River from the east (Figure 1). The study area is a part from Dibdibba Basin where Dibdibba formation is exposed. Dibdibba Basin (2485 km2) covers a part of Karbala Najaf region and can be represented by a triangle plateau surrounded by Tar-Al-Sayad from the west, Tar-AL-Najaf from the south, and by the flood plain from the east. This
Figure 1. Location of the study area.
basin has a dip that ranges from 0.5% to 0.1% from the west to the east (Figure 2)  .
Dibdibba formation is composed of poorly sorted sand, sand stone with gravel. The mineral composition of Dibdibba formation consists of 84.2% Quartz, 8.5% rock fragment, and 7.3% feldspar  . This (Pliocene-Pleistocene) formation is underlain by Injana (Upper Miocene) formation, which is in turn underlain by Nfayil (Middle Miocene) formation, then Euphrates (Early Miocene) formation as shown in Figure 2  . Injana formation in the study area represents a confining or semi-confining layer for Dibdibba formation  . This forces the ground water to flow horizontally in the dip direction of Dibdibba formation through its thickness which ranges, in general, from 1 meter or less to about 18 meters  . Within the studied area however, it ranges from 25 m to 70 m  .
Figure 2.(a) Boundaries of Dibdibba Basin; (b) cross section inside Dibdibba formation .
3. Review of Literature
Number of researchers realized the importance of the Dibdibba basin due to: Its arid or semi-arid climate, which leads to low rains and limited rain recharge, and low fluctuation in groundwater table. Beside that there are many reasons incentivize researchers to study this area, like:
1) It is located within a major city (Karbala).
2) Absence of surface water resources in the area where it is far from Euphrates and Tigris rivers.
3) Its arid or semi-arid area, having low rain, limited recharge, and low fluctuation in groundwater table.
4) Underlain by Injana formation and having a free water table, making it unconfined basin aquifer type, which facilitates the measuring process.
5) Shallow ground water, which produces an effective and quick recharge process.
6) Available hydrological and geo-hydrological information of the area.
7) Presence the natural boundaries around Dibdibba Basin, and having a dip from the west to east toward the flood plain lead to more accurate estimation process.
8) The increasing usage of the groundwater within the area.
Al-Ani  estimated the groundwater recharge of Dibdibba aquifer by using four techniques (groundwater fluctuation, water budget, numerical method/ Mod Flow/simulation model, and tracers). He found out that the tracer technique gave the closest results to the water exploitation in that area. Al-Mussawi  found the groundwater contour maps for the study area using Kriging and Inverse Distance Weighted (IDW) methods, and demonstrated the groundwater table decreasing in the direction from the west to the east. Ibrahim  studied the hydrochemistry of the area, and found high concentration of Ca2+ and SO42− in addition to heavy metal Zn and Al. Al-Saadi  classified the groundwater of Karbala Governorate according to WHO and Iraqi standards, US Salinity hazard laboratory classification  , Wilcox specification  , and permeability index (PI by Doneen  ) classifications. According to those classification he found that the groundwater in some places is suitable and in others not suitable for drinking and irrigation purposes. Ramadhan  evaluated the groundwater recharge in the study area using four methods: fluctuating water table, water balance of the basin, numerical modelling, and balance of chloride ion mass in unsaturated zone. He found that the values from the last method (chloride ion mass balance) are the most accurate values comparatively to the volume of water from the renewed inventories. Khalaf  located the best places to drill irrigation wells in the study area according to the topography of the region, sodium ratio, total dissolved salt, and groundwater elevation. In addition, he stated that there are two aquifers in the area: Dibdibba and Dammam. Thabit  generated 3-D image that simulate the resistivity for an area of 30 m * 30 m, then he used this image to demonstrate the seepage velocity-depth relation. He noticed that the maximum seepage velocity occurred in the depth ranged from 3 to 6 m below the ground surface, and the seepage velocity direction is from south west to north east of that area.
Despite the number of published papers, no article was found to evaluate the magnitude of seepage velocity of the study area (20 km * 20 km) shown in Figure 1. Therefore this article was performed to estimate this important geohydrological property.
Seepage velocity refers to the velocity of the flowing liquid through the pores of the porous media. It can be calculated by dividing Darcy’s velocity by the effective porosity, as follow  :
v: The seepage velocity (m/d),
q: Darcy’s velocity (m/d),
n: The effective porosity.
Darcy’s velocity can be found by the following equation:
K: The hydraulic conductivity (m/d),
: Head gradient (the change in the hydraulic head per unit length in the direction of the flow in an isotopic aquifer).
Since, hydraulic conductivity is equal to the transmissivity (T) divided by the effective thickness of the aquifer (b), therefore, Equation (1) can be rewritten as follow:
To generate seepage velocity maps (magnitude and direction maps) by GIS ArcMap software, four maps should be found first, they are:
1) Groundwater elevation head map.
2) Aquifer effective porosity map.
3) Aquifer saturated thickness map.
4) Aquifer transmissivity map.
All these maps should be in a raster type (with floating points) and have the same extent and cell size; otherwise the calculation process cannot be executed. Also all the dimensional parameters should be homogenous (i.e. the same unit for length and time). Ground water elevation (in the equations) is taken from a certain datum, which is usually the sea level.
The resulting maps after performing the commands (Darcy Velocity) in GIS/ArcMap are two maps: the first is seepage velocity magnitude map and the other is seepage velocity direction map. The direction map can be explained in an arrows map, while the magnitude map is explained as a raster map with floating values.
The process for finding Darcy’s maps includes the projection of the wells on the UTM coordinate system map either by interpolation using Kriging method which is more accurate or IDW which is less accurate. In this work Kriging method was used. Maps representing the four parameters (groundwater head elevation, effective porosity, saturated thickness, and transmissivity) were constructed. The same cell size for all maps was used. The input of these four maps was used in the command interference concerning Darcy velocity, to get the seepage velocity maps (magnitude and direction).
5. Results and Discussion
After collecting the necessary information to establish the four primary maps that are required for generate Darcy maps, GIS/ArcMap software was used to generate these maps using Kriging command (Figures 3-6). Figure 3, shows the groundwater head elevation above sea level in meter. From the figure it can be noticed that the water table dips toward the east. According to that, the ground water is flowing from the west to the east.
The effective porosity of aquifer in the study area is shown in Figure 4. Since there are no measured values for the porosity within this area, Todd  , and (GIS/ArcMap/Help) porosity ranges for different porous media were used. It had assumed that the average value for the porosity was 37%. Figure 5 shows the
Figure 3. Groundwater head elevation (m.a.s.l).
Figure 4. Effective porosity.
Figure 5. Saturated thickness (m) for the study area.
Figure 6. Transmissivity (m2/d).
saturated thickness of Dibdiba formation in the study area. The figure indicates that the saturated thickness of the western region is more than that the eastern region. The saturated thickness dips from 45 m in the west to 18.75 m in the east. Figure 6 demonstrates the transmissivity of the study area in (m2/day). It can be conducted that the transmissivity increase in the east and north east direction. And it ranges from 29 m2/d to 426 m2/d.
After producing the four maps (groundwater head elevation, effective porosity, saturated thickness, and transmissivity maps), they were used to generate seepage velocity maps (magnitude and direction maps), through inputting them in the command interface of Darcy Velocity. The resulting maps are shown in Figure 7 and Figure 8. Figure 7 shows that the seepage velocity range is 0 to 0.18 (m/d). Figure 8 reveals the seepage velocity direction. The overall flow direction is from the west to the east. The reasons beyond that are the dip of the Injana upper surface which represents an impervious bed underlying Dibdibba formation and the dip of water table toward the east direction. Figure 9 represent the combination of the seepage velocity direction map and the groundwater head elevation map. In Figure 9 the seepage velocity direction is perpendicular to the groundwater potential line which it is a reasonable result. Figure 10 represent the combination of the seepage velocity maps (magnitude and direction).
Figure 7. Seepage velocity magnitude.
Figure 8. Seepage velocity direction.
Figure 9. Seepage velocity direction and the groundwater head elevation map.
Figure 10. Seepage velocity maps (magnitude and direction).
Surface water is the main fresh water resource in Iraq. During the last years Iraq surface water was decreased dramatically. To solve or minimize the consequences of this shortage, other resources that can fulfill the municipality and agriculture requirements are to be used. One of the most important cities in the middle of Iraq is Karbala, due to its religious importance and population which exceeds one million. This city mainly depends on groundwater resources. The results indicated that the groundwater flows from the west to the east toward Mesopotamia Plain. This is due to the fact that Injana formation, which represents the bed of the unconfined Dibdibba aquifer, is dipping in that direction, and the water table is dipping in that direction also. The seepage velocity ranged from 0 m/d to 0.18 m/d. This variation in the seepage velocity is due to the changes in the geo-hydrological properties which govern the seepage velocity (i.e., groundwater elevation head, effective porosity, effective thickness, and transmissivity).
 Dhakras, P.A. (2011) Nanotechnology Applications in Water Purification and Waste Water Treatment: A Review. International Conference on Nanoscience, Engineering and Technology, Tamilnadu, 285-291.
 Issa, I.E., Al-Ansari, N.A., Sherwany, G. and Knutsson, S. (2014) Expected Future of Water Resources within Tigris—Euphrates Rivers Basin in Iraq. Journal of Water Resource and Protection, 6, 421-432.
 Al-Ansari, N.A., Ali, A. and Knutsson, S. (2014) Present Conditions and Future Challenges of Water Resources Problems in Iraq. Journal of Water Resources and Protection, 6, 1066-1098.
 Obaid, H.A., Shahid, S., Basim, K.N. and Shreeshivadasan, C. (2014) Modeling Sewerage Overflow in an Urban Residential Area Using Strom Water Management Model. Malaysian Journal of Civil Engineering, 26, 163-171.
 Ramadhan, A., Aljabbari, M. and Alkubaisi, R. (2013) Evaluation of Recharging Arid and Semi-Arid Regions (Case Study of Dibdibba Formation, Karbala Plateu-Najaf). Iraqi Journal of Science, 54, 902-910.
 Environmental Systems Research Institute (2017) How Darcy Flow and Darcy Velocity Work.