In the last decades, big developments took place in all fields of knowledge and industry with intensive attention on investment of all available resources to get the maximum benefit with minimum cost and minimum side effects on the environment.
Groundwater is one of the main water resources on the earth   , especially for arid and semi-arid countries for two reasons. The first is the scarcity of both surface water, and the second is the relatively low susceptibility of groundwater to pollution  . After the industrial revolution the world witnessed many changes due to the global warming as a result of the large development in industry. Consequently, many countries were affected and became arid or semi-arid. Middle East and North Africa countries are on the top of the affected countries. Some studies highlighted the size of the problem by comparing the annual water share for each person. For Iraq, it was found that the share reached about 6029 (m3/capita・year) during 1995, and it was expected to be 2100 (m3/capita・year) in 2015  .
Recent work indicates that Iraq is facing water crises problems and its water resources are declining due to climate change and building of dams upstream the Tigris and Euphrates rivers  -  . Due to this fact, Iraq will be depending more on its groundwater resources. Since groundwater is threatened to be contaminated due to natural reasons or different human activities (e.g. residential, industrial, agricultural activities municipal and commercial)  , in addition for the fact that water quality is an effective measure for human welfare  , then it becomes very important to study the groundwater potential to be contaminated and try to avoid contamination hazards.
There are two types of vulnerability intrinsic and specific vulnerability. Intrinsic vulnerability may be defined as how much easy for the contaminant that is found on the surface ground to be carried by the infiltrated water and diffused in to the groundwater   , while the specific vulnerability could be defined as the ability of a specific contaminant or a group of contaminants to be transported by the infiltrated water and diffuses in groundwater   .
Some researchers studied the vulnerability of groundwater in various parts of Iraq   , but none of them studied the groundwater conditions in Babylon Governorate. The groundwater in this area is relatively shallow and can be easily contaminated. For this reason this research was performed to see the vulnerability of groundwater to be polluted in the northern part of Babylon Governorate. This can help decision makers for putting prudent management planning for the groundwater resources in this governorate.
2. Study Area
The studied area represents the north part of Babylon governorate, the famous ancient city which is located in the middle part of Iraq about 100 km south of the capital Baghdad (Figure 1). It is bounded between longitudes (44˚2'42.245''E and 45˚2'2.964''E) and latitudes (32˚25'55.287''N and 33˚7'34.229''N) (Figure 1). The study area includes the two northern administrative units from the five units of Babylon governorate, which are Al-Musayyab Qadhaa’ and Al-Mahawil Qadhaa’ (Figure 1) covering an area of about 2794 km2. Its population is about 710832 capita  . It is bordered by Baghdad from the north, Wasit from the east, Karbala and Anbar from the west, and the southern part of Babylon from the south (Figure 1). The Euphrates River flows through this area for about 52 km (with 34 km before the Hindyia barrage which controls Euphrates flow into the downstream water courses: Euphrates, Shatt Al-Hilla, Kifil, Hussainina, Beni Hassan), and Tigris River flows for a distance of about 112 km about 20 km to the east of this area.
Figure 1. Study area.
The study area forms 2.4% of the Mesopotamia plain (Figure 2) which covers an area of about 116,000 km2  . The general topography of this area is characterized by its imperceptible gradient from northwest to the southeast. Where the highest point in this area is about 62 meter above sea level (m.a.s.l.) in its northern part, and the lowest point is about 21 (m.a.s.l.) in its southern part. This coincide with the slope of Mesopotamia Plain, where its highest point is 140 (m.a.s.l.), in Fatha vicinity in the north, whereas its lowest point is about 1 (m.a.s.l.), in the extreme southeastern margin along the Gulf  .
Generally, the soil in this area is alluvial fluvial silty clayey loam similar to the whole Mesopotamia Plain soil, where its origin is the fluvial deposition from Euphrates and Tigris Rivers. The climate of this area is arid or semiarid. It is extremely hot in summer, mild winter and short transition seasons of spring and autumn. The mean annual rainfall ranges between 200 mm to 125 mm from the north to the south, and the average annual evaporation from this area coincide with the general evaporation range for the Mesopotamia Plain which is 2600 mm in the northern parts to 3250 mm in the southern parts (Figure 2)  .
Figure 2. Location of the study area within the Mesopotamia plain  .
The stratigraphy of the Mesopotamian region was summarized by Yacoub  . The sequence of deposits can be summarized as follow (Figure 3).
1) Pliocene-Early Pleistocene Rock Units: These are represented by Bai Hassan Formation, Dibdibba Formation, Mahmudiya Formation. Bai Hassan is the formation that is found in the studied area, where the other formations are located in other locations in the Mesopotamia Plain.
2) Pleistocene Sediments (Rock Units): These are sediments of the Mesopotamia Fluvial Basin, alluvial fan sediments, and river terraces (Figure 4). The study area contains fluvial sediments and some of river terraces.
Figure 3. Subsurface geological cross section between Iskandariyah?Mandali, Northern Mesopotamia plain  .
Figure 4. Euphrates river terraces in the vicinity of Iskandariyah after  .
3) Pleistocene-Holocene Sequence: These are Sheet Run-off sediments, slope sediments, and gypcrete. This type of sediment is not observed in the study area.
4) Holocene Sequence: This is the existing sedimentary environment. This type of sediments covers the majority of Mesopotamia plain (Figure 5). They make the upper part (about 15 - 20 m) of the Quaternary sediments of the Mesopotamia Basin. The Holocene sequence contains different types like fluvial, deltaic, lacustrine and estuarine/marine units. But the types that cover the study area are mostly flood plain of Euphrates and Tigris Rivers, some Aeolian sediment exists in the north and south of the studied area. There are some areas in the north of the studied area having exposures of Bai Hassan formation, which represents Pliocene-Early Pleistocene sediments (Figure 5). The depth of Quaternary deposits (Pleistocene deposits and Holocene deposits) reaches in some sites up to 100 m (see Figure 3).
Figure 5. Geological map of the Mesopotamia plain  .
The aquifer-aquitard in Mesopotamia Plain is considered lithological complex aquifer system  . This is due to the fact that there are abrupt lithological changes in the Mesopotamian plan which affects the hydrogeological properties. These effects repre- sented by the heterogeneous lithification (laterally and vertically) on the sediments which lead to heterogeneous hydraulic conductivity in the plain, consequently difficulty in aquifer-aquitard delineation.
2.2. Population and Water Use
The population of the study area is 710,832. The major cities cover an area of 32.12 km2 with a population density greater than 5000 (capita/km2), and the remainder of the area (2691 km2) has population density less than that (see Figure 6).
The main share of the used water in this area is consumed in the agriculture purposes with the fact of presence about 2961 km2 arable area.
Figure 6. North Babylon population distribution.
2.3. Hydrology of the Mesopotamia Plain
The groundwater piezometric level for the plain is based on the assumption of the hydraulic continuity within the entire Quaternary aquifer system. This means that all the aquifers in the plain are in a hydraulic continuity, which depends on the degree of lithification. Furthermore, there are two more facts related to the hydrogeological conditions in the sedimentary plain. The first is the hydraulic continuity between the aquifer complex system and the rivers in the plain and the second is the hydraulic continuity between the aquifer system and the deep formation underlying the fluvial deposits. Thus, depending on these assumptions the water table can be drawn in Mesopotamia Plain  .
In this study the DRASTIC model with the aid of GIS program was used to evaluate the intrinsic vulnerability of the aquifer of the study area   . DRASTIC model was first developed by the United State Environmental Protection Agency (EPA) for evaluating the groundwater pollution potential    , and during the last years it was widely used  , because it gives a good indication for the groundwater sensitivity for the contaminants, and it is easy to get the required information for application of this model on groundwater  .
The concept of DRASTIC model is that groundwater is affected by the contaminants that are carried by infiltrated water from the surface of the ground. In this model the contaminants are assumed to be on the ground surface, they are flushed into the groundwater by precipitation, and they have the mobility of the water  .
According to DRASTIC model, the intrinsic vulnerability of the aquifer is affected by seven geological, hydrological parameters. These can cause changes of the groundwater quality of the aquifer  . The seven parameters are collected in the acronym DRASTIC which refers to: (D) depth to groundwater (which means depth from the ground surface to the water table), (R) net recharge, (A) aquifer media, (S) soil media, (T) topography, (I) impact of vadose zone media, and (C) hydraulic conductivity of the aquifer   .
In DRASTIC model each parameter has two values the first value is the rating value and the second is the weighting value. The rating value for each parameter ranges from 1 to 10 depending on the different effects of the variable types of that parameter on the vulnerability of the aquifer. The weighting value of each parameter is a constant value in the model, i.e. it will be constant for different rating values of that parameter.
The parameter weight value depends on the potential effect of that parameter on the groundwater vulnerability. Table 1 and Table 2 show the weight values with a brief description for each parameter and how it affects the vulnerability of the aquifer and the rating values for different cases for each parameter can be seen respectively.
DRASTIC index is the output of the DRASTIC model, which represents the vulnerability index of the groundwater, i.e. the higher DRASTIC index the greater chance for the groundwater to be affected by the pollution of the contaminants. The DRASTIC
Table 1. Assigned weight for DRASTIC parameters  .
Table 2. DRASTIC parameters and their rating scale  .
model can be written as    :
where D, R, A, S, T, I, and C are the seven parameters, r and w refer to rating and weight values for the corresponding parameter.
After collecting all the required data to calculate the DRASTIC Index for a specific region, ArcMap10.3 of GIS-Software is needed to perform DRASTIC model on that region. The ArcMap10.3 gives the capabilities to construct a digital map (raster map) for each parameter of the seven parameters. The constructed maps work as layers for the same region. Then by using Raster Calculator Tool in Map Algebra Tool Box from Spatial Analyst Tools, an integrated map for the DRASTIC Index can be generated. The resultant map represents the Vulnerability index map of the region. From the literature        the vulnerability can be divided into seven classes depending on RASTIC index, these classes are summarized in Table 3. From Table 3, it can be seen that high DRASTIC index corresponds to high vulnerability, i.e. groundwater is easily affected by the contaminants.
4. Results and Discussion
After collecting and inputting all the required data to find the aquifer vulnerability for north Babylon region, Arc Map software was used to construct seven rated layers for that region representing seven parameters’ maps, with one layer for each parameter. The resulting maps were two for each variable (before the rating and after the rating).
4.1. Depth to Groundwater Table Map
The results from the Arc Map program showed that the area of north Babylon is 2794.42 km2. The unrated map for the depth to the groundwater is shown in Figure 7(a). The results indicate that most of the region (about 87% from the area) has a depth to groundwater ranging from 1.5 to 4.5 m, and a small area (about 11%) has groundwater depth with a range of 4.5 to 9.1 m and a very small area with depth to ground-
Table 3. Vulnerability classification.
Figure 7. Depth to groundwater map ((a)-before the rating; (b)-after the rating).
water less than 1.5 m. By using the Arc Map toolboxes the rated map for the depth to the groundwater was produced, and it is shown in Figure 7. The rated map (Figure 7(b)) shows that there are three zones (depending on the rating) in the studied area which are 7, 9 and 10 zones.
4.2. Net Recharge Map
Figure 8(a) shows the unrated recharge map of the studied area. It is obvious from the map that all the area has a net recharge less than 50 mm/year. Consequently the rated map for the recharge is one zone map which has rating equal to 1 (Figure 8(b)).
4.3. Aquifer Media Map
The aquifer media for the studied area is from one type which is silty clay, as shown in Figure 9(a). Depending on that classification the rated map of the aquifer media has also one type rating which is 6 (Figure 9(b)).
4.4. Soil Map
According to Table 1 and the soil map for the study area (Figure 10(a)), there are three type of the soil rating in the studied area. These are clayey loam, silty loam and sandy loam. The results of ArcMap indicate that the area is covered by clayey loam soil (20%), silty loam soil (74%), and sandy loam soil (5%). Consequently the rated map of the soil is shown in Figure 10(b).
Figure 8. Recharge map ((a)-before the rating; (b)-after the rating).
Figure 9. Aquifer media map ((a)-before the rating; (b)-after the rating).
Figure 10. Soil map ((a)-before the rating; (b)-after the rating).
4.5. Topographic Map (Slope Map)
The topographic map can be produced by using the Digital Elevation Map (DEM) and the Arc toolboxes in ArcMap. The topographic map of the area is shown in Figure 11(a). The results show that most of the area (which is about 93%) has a slope ranging between 0% and 2%, while an area about 6.5% has a slope ranging from 2% to 6%. Finally, a very small area (0.026%) has a slope ranging between 6% and 10%. According to the unrated map of the slope the topographic rated map can be drawn as in Figure 11(b), in which three types of zones can be distinguished according to the slope value which is 5, 9 and 10 rated zones.
4.6. Impact of Vadose Zone Map
Figure 12(a) shows the unrated map of impact of vadose zone, in this figure the area is classified as one type, because the entire region has one type of vadose zone which is silty clay. Due to that classification the resultant rated map for impact of vadose zone in Figure 12(b) is also one type and has the rating 3 according to the classification in Table 2.
4.7. Hydraulic Conductivity Map
Figure 13(a) represents the unrated map of hydraulic conductivity, according to the results, the region is divided into 3 zones due to the difference in the hydraulic
Figure 11. Topographic map ((a)-before the rating; (b)-after the rating).
Figure 12. Impact of Vadose zone map ((a)-before the rating; (b)-after the rating).
Figure 13. Hydraulic conductivity map ((a)-before the rating; (b)-after the rating).
conductivity of the aquifer. These zones are less than (4.1), (4.1 - 12.2) and (12.2 - 28.5) m/d hydraulic conductivity. The results indicate that about 87% of the total area has hydraulic conductivity of less than 4.1 m/d. The remainder area (10%) has hydraulic conductivity between 4.1 and 12.2 m/d, and about 1.5% has hydraulic conductivity between 12.2 and 28.5 m/d. Depending on these values the rated map of hydraulic conductivity was produced (Figure 13(b)). The resultant rated map shows three zones with rating 1, 2 and 4.
Table 4 shows some statistical values for the rated maps for all the seven parameters. The integrated vulnerability map is shown in Figure 14. The map indicates that the region is divided into two zones according to the vulnerability. The first and second zones have vulnerability of (86 - 100) and (100 - 120) respectively. So the region in general is either has very low vulnerability or low vulnerability. This is due to the low rating values of the factors or because of the low weight values of these factors, especially recharge factor and impact of vadose zone, where they have high weight values but low rating values.
Due to water shortage problems in Iraq, it is expected that the demand of using groundwater resources will increase with time. Pollution of groundwater will create severe consequences. DRASTIC model had been applied to look at the vulnerability of
Figure 14. Vulnerability map.
Table 4. The statistical values for the rated maps.
groundwater to pollution in the northern part of Babylon Governorate in Iraq. The results indicated that the area is divided into two regions according to its vulnerability. It either has very low vulnerability or low vulnerability. This leads to the possibility to utilize the groundwater in this region with small risks on its chemical or physical properties. It is noteworthy to mention that further checking is required due to the fact that the low rating values of the factors or because of the low weight values of these factors.
The authors would like to thank the Iraqi Government for awarding a scholarship for the first author. Mr. Twana Abdullah, Mr. Maher Zainy and Mr. Ali Chabuk gratefully helped the authors.
 Villeneuve, J., Banton, O. and Lafrance, P. (1990) A Probabilistic Approach for the Groundwater Vulnerability to Contamination by Pesticides: The Vulpest Model. Ecological Modelling, 51, 47-58.
 Sinha, M.K., Verma, M.K., Ahmad, I., Baier, K. Jha, R. and Azzam, R. (2016) Assessment of Groundwater Vulnerability Using Modified DRASTIC Model in Kharun Basin, Chhattisgarh, India. Arabian Journal of Geosciences, 9, 1-22.
 Jamrah, A., Al-Futaisi, A., Rajmohan, N. and Al-Yaroubi, S. (2008) Assessment of Groundwater Vulnerability in the Coastal Region of Oman Using DRASTIC Index Method in GIS Environment. Environmental Monitoring and Assessment, 147, 125-138.
 Hussain, H.M., Al-Haidarey, M., Al-Ansari, N. and Knutsson, S. (2014) Evaluation and Mapping Groundwater Suitability for Irrigation Using GIS in Najaf Governorate, Iraq. Journal of Environmental Hydrology, 22, 1-16.
 Al-Ansari, N., Abdellatif, M., Ali, S. and Knutsson, S. (2014) Long Term Effect of Climate Change on Rainfall in Northwest Iraq. Central European Journal of Engineering, 4, 250-263.
 Al-Ansari, N., Abdellatif, M., Al-Khateeb, M. and Knutsson, S. (2014) Desertification and Future Rainfall Trends North Iraq. 8th Edition of the International Scientific Congress of GIS and Geospace Applications Geotunis, Tunis/Hammamet, 2-6 April 2014.
 Al-Ansari, N.A., Abdellatif, M., Zakaria, S., Mustafa, Y. and Knutsson, S. (2014) Future Prospects for Macro Rainwater Harvesting (RWH) Technique in North East Iraq. Journal of Water Resource and Protection, 6, 403-420.
 Al-Ansari, N.A., Abdellatif, M.,Ezeelden, M.,Ali, S. and Knutsson, S. (2014) Climate Change and Future Long Term Trends of Rainfall at North-Eastern Part of Iraq. Journal of Civil Engineering and Architecture, 8, 790-805.
 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.
 Al-Ansari, N.A., Ali, A.A. and Knutsson, S. (2015) Iraq Water Resources Planning: Perspectives and Prognoses. 13th International Conference on Civil and Construction Engineering, Jeddah, 26-27 January 2015, 2097-2108.
 Baier, K., Schmitz, K.S., Azzam, R. and Strohschon, R. (2014) Management Tools for Sustainable Ground Water Protection in Mega Urban Areas-Small Scale Land Use and Ground Water Vulnerability Analyses in Guangzhou, China. International Journal of Environmental Research, 8, 249-262.
 Balakrishnan, P., Saleem, A. and Mallikarjun,N. (2011) Groundwater Quality Mapping Using Geographic Information System (GIS): A Case Study of Gulbarga City, Karnataka, India. African Journal of Environmental Science and Technology, 5, 1069-1084.
 Barzegar, R., Moghaddam, A.A. and Baghban, H. (2016) A Supervised Committee Machine Artificial Intelligent for Improving DRASTIC Method to Assess Groundwater Contamination Risk: A Case Study from Tabriz Plain Aquifer, Iran. Stochastic Environmental Research and Risk Assessment, 30, 883-899.
 Vrba, J. and Zaporozzec, A. (1994) Guidebook on Mapping Groundwater Vulnerability. (IAH)International Association of Hydrogeologists, IAH International Contributions to Hydrogeology, 16,131pp. Heise Publ., Hannover.
 Gogu, R.C. and Dassargues, A. (2000) Current Trends and Future Challenges in Groundwater Vulnerability Assessment Using Overlay and Index Methods. Environmental Geology, 39, 549-559.
 Abdullah, T.O., Ali, S.S., Al-Ansari, N.A. and Knutsson, S. (2015) Groundwater Vulnerability Mapping Using Lineament Density on Standard DRASTIC Model: Case Study in Halabja Saidsadiq Basin, Kurdistan Region, Iraq. Engineering, 7, 644-667.
 Lindblom, J., Al-Ansari, N. and Al-Madhlom, Q. (2016) Possibilities of Reducing Energy Consumption by Optimization of Ground Source Heat Pump Systems in Babylon, Iraq. Engineering, 8, 130-139.
 Al-Adamat, R., Foster, I.D.L. and Baban, S.M.J. (2003) Groundwater Vulnerability and Risk Mapping for the Basaltic Aquifer of the Azraq Basin of Jordan Using GIS, Remote Sensing and DRASTIC. Applied Geography, 23, 303-324.
 Babiker, I.S., Mohamed, M.A.A., Hiyama, T. and Kato, K. (2005) A GIS-Based DRASTIC Model for Assessing Aquifer Vulnerability in Kakamigahara Heights, Gifu Prefecture, Central Japan. Science of the Total Environment, 345, 127-140.
 Aller, L., Lehr, J.H. and Petty, R. (1987) Drastic: A Standardized System to Evaluate Ground Water Pollution Potential Using Hydrogeologic Settings. National Water Well Association Worthington.
 Al-Zabet, T. (2002) Evaluation of Aquifer Vulnerability to Contamination Potential Using the DRASTIC Method. Environmental Geology, 43, 203-208.
 Shrestha, S., Semkuyu, D.J. and Pandey, V.P. (2016) Assessment of Groundwater Vulnerability and Risk to Pollution in Kathmandu Valley, Nepal. Science of the Total Environment, 556, 23-35.
 Zhang, B., Li, G., Cheng, P., Yeh, T.C. and Hong, M. (2016) Landfill Risk Assessment on Groundwater Based on Vulnerability and Pollution Index. Water Resources Management, 30, 1465-1480.
 Al-Rawabdeh, A.M., Al-Ansari, N., Knutsson, S., Al-Taani, A. and Al-Khateeb, F. (2014) Modeling the Risk of Groundwater Contamination Using Modified DRASTIC and GIS in Amman-Zerqa Basin, Jordan. Central European Journal of Engineering, 4, 264-280.