International Journal of Geosciences
Vol.07 No.05(2016), Article ID:66882,10 pages
10.4236/ijg.2016.75055

Protective Capacity Assessment of Vadose Zone Material by Geo-Electrical Method: A Case Study of Pakistan

Yawar Hussain1, Sadia Fida Ullah2, Muhammad Babar Hussain3, Hernan Martinez-Carvajal1,4, Abdul Qayyum Aslam5

1Department of Civil and Environmental Engineering, University of Brasilia, Brasilia, Brazil

2Department of Biological Sciences, University of Brasilia, Brasilia, Brazil

3Department of Theoretical and Astrophysics, Xiamen University, Xiamen, China

4Department of Civil Engineering, National University of Colombia at Medellín, Medellín, Colombia

5College of Earth and Environmental Sciences, University of the Punjab, Lahore, Pakistan

Copyright © 2016 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

Received 22 March 2016; accepted 23 May 2016; published 27 May 2016

ABSTRACT

Rapid and poorly managed agricultural landuse in Punjab plain is a major contention in groundwater quality deterioration. This study is focused on the determination of groundwater vulnerability potential of vadose material using geo-electrical method. Total of eight electrical resistivity soundings with Schlumberger electrode configuration are taken near Khizarabad Dairy Form, district Sargodha. The instrument ABEM Terrameter SAS 300 is used for data acquisition. The quantitative interpretation is done by 1D modelling of resistivity soundings by using IX1D software. The results characterize the area into three geo-electrical layers as dry alluvium (15.7 - 32.4 Ωm), saturated alluvium with high resistivity (10 - 20 Ωm) and low resistivity (>10 Ωm). Aquifer vulnerability is assessed by calculating Dar Zarrouk parameter (Longitudinal Conductance) and thickness of vadose zone. Resultant vulnerability is divided into three zones as low (>3.8 siemens) moderate (1.6 - 2.5 siemens) and high (<1.6 siemens). It is concluded from this study that Vertical Resistivity Sounding (VES) can be used for vulnerability assessment of an unconfined sandy alluvial aquifer from agricultural contaminants.

Keywords:

Electrical Sounding, Vulnerability, Dar Zarrouk Parameter, Groundwater, Khizarabad Dairy, Sargodha

1. Introduction

The greater Punjab plain which is divided into five Doabs (land between two rivers) is suffering from agricultural stresses, because of rapid growth of agriculture and population. Groundwater is major source of water for growing agriculture and population in Pakistan. Groundwater is present in voids and it fully saturates it; porosity and permeability of that stratum have greater importance in groundwater recovery and protection [1] . Along with production, the groundwater protection is also very important, particularly in arid regions of the world. Agricultural Non-point Pollution Sources (NPS) and urban wastes are the major threat to the portability of groundwater in Punjab, Pakistan. An assessment of these NPSs on groundwater is an important environmental concern.

Several methods of aquifer vulnerability have been proposed; among these methods, index methods are the most popular. There are different such models but the most useful is DRASTIC model. Scientific history of the world has witnessed many examples of development and applications of index methods like DRASTIC, GOD, IRIH, pollutants [2] . Aquifer vulnerability assessment by S-model is a cheaper alternative to index models. In S- model the protective capacity of vadose zone is assessed by determining lithological composition and thickness through electrical resistivity sounding. Longitudinal conductance (S) values are grouped into different vulnerability zones. The resistivities of earth material depend on soil type, porosity and fluid present in pore spaces and can be determined by electrical resistivity surveying, a geophysical technique used for the determination of such resistivities variations [3] . The use of geophysical methods in groundwater studies has improved greatly because of development in computers [4] . Electrical resistivity sounding is economical and highly suitable for hydrological, geotechnical, mining and environmental studies. There are different geophysical methods used for groundwater development through determination of formation resistivity, properties of a fluid inside pores and protective capacities of material above the saturated zone [3] [5] . VES has proved a non-distractive and appropriate method for determination of various subsurface geologic strata and groundwater quality. This technique, VES with Schlumberger electrode configuration is very famous for the groundwater studies in arid regions [6] . ERS has been used in various groundwater studies for the achievement of multiple objectives in numerous hydrogeological conditions [7] . Some of important studies are, groundwater quality [8] , delineation of fresh and saline water boundary [9] , for imaging the industrial contamination plume [10] , determination of aquifer parameters [11] [12] and hydraulic conductivity and transmissivity [13] .

The present study aimed at determining the pollution potential of the unsaturated layer present above the aquifer layer. Electrical resistivity method is utilized by adopting Schlumberger array. Dar Zarrouk parameter’s longitudinal conductance (S) values, which reflect the permeability of unsaturated layer, provide the vulnerability of aquifer to the surface contaminates.

2. Material and Methods

2.1. Study Area

Study area (Khezarabad Diary Form) lies in district Sargodha of Punjab at geographic coordinates Latitude 32˚14'3.01"N and Longitude 72˚45'30.76"E (Figure 1). District Sargodha is part of lower portion of Chaj Doabs, alluvium between Jhelum and Chenab River. There are two crop seasons as Rabbi and Kharif. The region is flat with well-developed canal system. Area is irrigated by tributaries from Jhelum River and motor driven pumps. This city is renowned for citrus production world-wide as climate conditions are highly suitable for its cultivation. Khizarabad livestock dairy farm has been selected as case study to assess protective capacity assessment of vados material zone. Case study research allows research to test hypothesis empirically. Availability of canal water is perennial study area consists of loamy soil (http://dlfpunjab.webs.com/khizerabad.htm).

2.2. Hydrometrology

The climate of the study area is arid as rainfall is less than evapotranspiration rate. The annual mean rainfall is 426 mm. Mean annual temperature is 24˚C with the mean maximum of 31˚C and mean minimum of 18˚C. Mean Relative Humidity at day time is 71% whereas the night time it is 79%. Study area is characterized by Monsoon showers during July and August contributing significantly in annual precipitation (Figure 2). In October and November, very little precipitation is recorded at this station.

Figure 1. Regional location of study area on Pakistan map with digital elevation model, major roads, rivers, canals, electrical sounding points and tube wells locations.

Figure 2. Average monthly rainfall recorded at Sargodha station for period 2001-2010.

2.3. Geology and Physiography

Geologically study area consists of thick alluvial deposits of Quaternary age over basement rocks. Major geological units are active flood plain and abandoned flood plain. There are also basement rocks exposed to surface [14] . The alluvium complex of Quaternary age deposited on metamorphic and igneous rocks of Precambrian age in Chaj Doab. The Alluvium consists of fine to medium sand with slit and minor clay lens. The strata lack both horizontal and vertical continuity. There are three important geological units as Potwar Losses, Piedment Deposits and Aeolian Deposits [15] . The underlain aquifer has good hydrogeological characteristics such as hydraulic conductivity and Transmissivity. Another important factor that influences hydrogeology of region is gradual exposure of bedrocks from central operation of the Doab to Kirana where bed rocks outcrop (Figure 3). This ridge influences greatly the groundwater flow everywhere in the Doab.

2.4. Concept of Vulnerability

It is dimensionless index which is considered as the function of hydrogeological factors, sources of contaminants and anthropogenic effects in a specific site [16] . The concept of groundwater vulnerability to contamination is based on the assumption that the physical environment may provide some degree of protection to ground- water against natural and human impacts with respect to contaminants in the groundwater. The vulnerability of a certain area can be described by the degree of susceptibility of that area to groundwater pollution [17] . In 1968, the French Margate was the first one who used the term vulnerability in Hydrogeology, thereafter, the concept was adopted worldwide [18] . Groundwater vulnerability to contamination is defined as “The sensitivity of groundwater quality to an imposed contaminant load, which is determined by the intrinsic characteristics of the aquifer” [19] . Groundwater vulnerability to contamination is the tendency or likelihood of contaminants to reach a specified position in the groundwater system after introduction at some location above the uppermost aquifer. As can be inferred from the above definition, groundwater vulnerability is not an absolute or measurable property, but an indication of the relative possibility with which contamination of groundwater resources will occur. This understanding implies a very basic vulnerability concept that all groundwater is vulnerable [2] . In the pre- sent study vulnerability is taken equal to protective capacity.

2.5. Resistivity Method

2.5.1. Electrical Resistivity Survey

Apparent resistivities in the field are measured by using SAS-300 equipment by GeoWatAssociates, Islamabad (PVT), a consulting company in July 2011 [20] . Vertical Electrical Soundings provide information of different geological section and vertical variations of different water quality zones [21] . In this survey Schlumberger array is used because of its data acquisition and interpretations are less labor intensive and economical [22] . In this

Figure 3. Physiographic units of Chaj Doab of which study area a small portion modified from (Hussain et al., 1990) [15] .

array the central potential electrodes are fixed and the other peripheral current electrodes are moved horizontally in order to target different depths. The two outer electrodes A and B are used for the current, and the resulting potential difference is measured across the two inner electrodes M and N. MN/2 is always kept sufficiently small relative to AB/2. The distance from center of current electrodes is high comparative to potential electrodes (Figure 4). Apparent or field resistivity is calculated by using Equation (1).

(1)

where ρa is apparent resistivity of earth material measured in ohm-meter, I is current applied in milliamp ere (mA) and K is geometric constant, its value depends on the electrode configurations. For Schlumberger electrode configuration its value can be calculated as

(2)

where AB is current electrode separation and MN is potential electrode separation. Resistivity values of dry rocks are very high while presence of clay and salted water in the pores reduces the resistivity significantly [12] .

Resistivity data can be interpreted both qualitatively and quantitatively. Qualitative interpretation can be performed by apparent resistivity curves plotted between filed resistivity values and different AB/2 values. This will provide an initial guess about the variation of resistivity in field. Qualitative interpretation of resistivity data is based on two basic assumptions that earth is horizontal and resistivity of last layer has infinite thickness [23] . Each layer in 1D model is assumed as uniform and isotropic.

Primary data about VES locations was collected from report [20] . Resistivity values are grouped into three zones as upper unsaturated, lower high resistivity and low resistivity based on the modeling results of resistivity models

IX1D is used for the interpretation of resistivity curves in terms of layered resistivity model by matching measured curves with theoretical master curves using master curve matching techniques [24] . These models provide resistivity and thickness of each layer (Figure 5). Geoelectrical section from the results of ID models is obtained along profile AB (Figure 6). These layered parameters are used for the calculation of Longitudinal Conductance (S) and resistivity of unsaturated zone. The Dar Zarrouk (DZ) Parameter (Longitudinal Conductance S) is calculated by layered parameters obtained from ID resistivity curve modeling. Following equation is used for the calculation of DZ parameter.

(3)

where S is used in for the assessment of protective capacity of aquifer, hi is the thickness of unsaturated layer and is the resistivity of that layer

2.6. Preparation and Generation Thematic Map

2.6.1. Unsaturated Thickness Map

Isopack thematic map is created by using IWD interpolation algorithm of ArcGIS 10 software. IWD works on the principle that value at unmeasured point can be predicted from neighboring measured points by assigning a weight factor, which varies inversely with the distance from measured point [25] .

Figure 4. Sketch diagram of Schlumberger configuration array, AB are current electrodes and MN are potential electrodes.

Figure 5. ID resistivity model of VES-1 created in IX1D resistivity software.

Figure 6. Sketchof geoelectrical section along profile AB.

2.6.2. Layer Resistivity Map

The first layer resistivity is interpolated. The resultant map gives resistivity of first layer at different locations in the study area which is then used in aquifer protective capacity assessment. According to Braga et al. (2014), lithology of unsaturated zone is difficult to map with ERS, however, litho-facies of saturated zone can be determined with this method. Study is underlain by unconfined alluvial sandy aquifer and in this case lithology of saturated and unsaturated can be considered as same [26] . The lithologies of unsaturated zone are built by extrapolation of lithology of saturated zone [27] [28] .

2.6.3. Longitudinal Conductance Map

Earth martial above water table acts as natural filter against contaminants passing through it. The vadose zone material lithology in term of permeability and porosity is very important in determining pollution potential of underlying aquifer. Such vulnerability potential is related directly with longitudinal conductance S for clayey aquifer, that can be used to determine protective capacity which is taken as infiltration time of contaminates. Here PC is related to S which is measured in units of meter ohm [5] .

Above layer parameters are used for the calculation of protective capacity as given in Equation (3). The input parameters as well as the protective capacity is evaluated in GIS environment by using ArcGIS 10.1. The algorithms utilized by current investigation is IWD present in Spatial Analyst Tools Module of ArcGIS 10.1 software. Raster Calculator is used for the protective capacity evaluation (Equation (3)). Other major tools are Raster to Polygon conversion, Geo-referencing, reclassify and Geometry calculation.

3. Results and Discussions

3.1. Thickness and Resistivity of Unsaturated Zone

The thickness of dry strata over aquifer varies from 14 (minimum value) to 57 meters depth (maximum value). The north western portion has shallow groundwater depth increases vulnerability of surface contaminates to reach aquifer by reducing their travel times. The south eastern and south western parts have very deep groundwater and provide greater protection to underlying aquifer (Figure 7). Higher the resistivity of material of vadose zone lower is the protective capacity (Equation (3)). Resistivity of first geoelectrical layer varies from less than 12 to greater than 15 Ωm (Figure 8). Major portion of area is occupied with resistivity range of 13 - 14 Ωm.

3.2. Protective Capacity Mapping

The aquifer protective capacity or vulnerability map is classified into three classes using quartile classification of ArcGIS 10.1. The resultant vulnerability is divided into three classes as low (>3.8 siemens), moderate (1.6 - 2.5 siemens) and high (<1.6 siemens). Major portion of study area is under moderate vulnerability zone (Figure 9). The high vulnerability is the reflection of shallow water table and low apparent resistivity of material overlie

Figure 7. Spatial distribution of thickness of unsaturated zone.

Figure 8. Spatial distribution of resistivity of 1st geoelectrical layer.

Figure 9. Protective capacity of unsaturated zone divided into different vulnerability zones.

the aquifer while reverse is true for low vulnerability zones.

4. Conclusion

This study was carried out to provide vulnerability assessment of unsaturated layer by using electrical resistivity surveying technique of geophysical method. Longitudinal conductance (S) is used for the calculation of relative aquifer protection from the surface contaminants. The resultant vulnerability capacity is divided into three zones high, moderate and low. The results indicate that the most part of study is under low to moderate vulnerability while a small portion is under high vulnerability. The areas having high vulnerability should be taken under consideration of Khezarabad Diary Form administration.

Acknowledgements

Authors would like to thank GeoWat Associates, Islamabad (PVT) for providing necessary data for this study. We would also like to thank Mis. Farzana Ahmad Deen, who made English corrections.

Cite this paper

Yawar Hussain,Sadia Fida Ullah,Muhammad Babar Hussain,Hernan Martinez-Carvajal,Abdul Qayyum Aslam,1 1, (2016) Protective Capacity Assessment of Vadose Zone Material by Geo-Electrical Method: A Case Study of Pakistan. International Journal of Geosciences,07,716-725. doi: 10.4236/ijg.2016.75055

References

  1. 1. Obiora, D.N., Ajala, A.E. and Ibuot, J.C. (2015) Investigation of Groundwater Flow Potential in Makurdi, North Central Nigeria, Using Surficial Electrical Resistivity Method. African Journal of Environmental Science and Technology, 9, 723-733.
    http://dx.doi.org/10.5897/AJEST2015.1983

  2. 2. Hussain, Y. (2014) Hydrogeophysical Investigations and GIS Vulnerability Mapping of Kot Addu, M.Phil. Dissertation, Quaid-i-Azam University, Islamabad, Pakistan

  3. 3. Daniel, A., Louis, O., Emmanuel, C. and Orji, N. (2015) Delineation of Potential Groundwater Zones Using Geo- Electrical Sounding Data at Awka in Anambra State, South-Eastern Nigeria. European Journal of Biotechnology and Bioscience, 3, 1-5.

  4. 4. Fadele, S.I., et al. (2013) The Use of Vertical Electrical Sounding (VES) for Groundwater Exploration around Nigerian College of Aviation Technology (NCAT), Zaria, Kaduna State, Nigeria. The Pacific Journal of Science and Technology, 13, 549-555.

  5. 5. Atakpo, E.A. and Ayolabi, E.A. (2009) Evaluation of Aquifer Vulnerability and the Protective Capacity in Some Oil Producing Communities of Western Niger Delta. The Environmentalist, 29, 310-317.
    http://dx.doi.org/10.1007/s10669-008-9191-3

  6. 6. Kelly, W.E. and Stanislav, M. (1993) Applied Geophysics in Hydrogeological and Engineering Practice. Elsevier, Amsterdam, London, New York, Tokyo.

  7. 7. Mosuro, G.O., et al. (2016) Assessment of Groundwater Vulnerability to Leachate Infiltration Using Electrical Resistivity Method. Applied Water Science, 1-13.

  8. 8. Sonkamble, S. (2014) Electrical Resistivity and Hydrochemical Indicators Distinguishing Chemical Characteristics of Subsurface Pollution at Cuddalore Coast, Tamil Nadu. Journal Geological Society of India, 83, 535-548.
    http://dx.doi.org/10.1007/s12594-014-0081-7

  9. 9. Swartx, J.H. (1937) Resistivity Studies of Some Salt-Water Boundaries in the Hawaiian Islands. Transactions of American Geophysical Union, 18, 387- 393.
    http://dx.doi.org/10.1029/TR018i002p00387

  10. 10. Lima, O.A.L., Sato, H.K. and Porsani, M.J. (1995) Imaging Industrial Contaminant Plumes with Resistivity Techniques. Journal of Applied Geophysics, 34, 93-108.
    http://dx.doi.org/10.1016/0926-9851(95)00014-3

  11. 11. Sikandar, P., Bakhsh, A., Arshad, M. and Rana, T. (2010) The Use of Vertical Electrical Sounding Resistivity Method for the Location of Low Salinity Groundwater for Irrigation in Chaj and Rachna Doabs. Environmental Earth Sciences, 60, 1113-1129.
    http://dx.doi.org/10.1007/s12665-009-0255-6

  12. 12. Srinivasan, K., Poongothai, S. and Chidambaram, S. (2013) Identification of Groundwater Potential Zone by Using GIS and Electrical Resistivity Techniques in and around the Wellington Reservoir, Cuddalore District, Tamilnadu, India. European Scientific Journal, 9, 312-331.

  13. 13. Amos-Uhegbu, C., Igboekwe, M.U., Chukwu, G.U., Okengwu, K.O. and Eke, T.K. (2014) Geoelectrical Delineation and Geochemical 2 Characteristics of Aquifer Systems in Kwa-Ibo River 3 Watershed, Abia State, Nigeria. Journal of Scientific Research and Reports, 3, 818-843.
    http://dx.doi.org/10.9734/JSRR/2014/6215

  14. 14. Greenman, D.W., Swarzenski, W.V. and Bennett, G.D. (1967) Groundwater Hydrology of the Punjab, West Pakistan with Emphasis on Problems Caused by Canal Irrigation. US Geological Survey Water-Supply Paper.

  15. 15. Hussain, S.D., Sajjad, M.I., Akram, W., Ahmad, M., Rafiq, M. and Tariq, J.A. (1990) Surface Water/Groundwater Relationship in Chaj Doab. Pakistan Institute of Nuclear Science and Technology, Pakistan.

  16. 16. Plymale, C.L. and Angle, M.P. (2002) Groundwater Pollution Potential of Fulton County, Ohio. Ohio Department of Natural Resources Division of Water, Water Resources Section, Groundwater Pollution Potential, Report No. 45.

  17. 17. Al Hallaq, A.H. and Elaish, B.S.A. (2012) Assessment of Aquifer Vulnerability to Contamination in Khanyounis Governorate, Gaza Strip-Palestine, Using the DRASTIC Model within GIS Environment. Arabian Journal of Geosciences, 5, 833-847.
    http://dx.doi.org/10.1007/s12517-011-0284-9

  18. 18. Albinet, M. and Margat, J. (1970) Groundwater Pollution Vulnerability Mapping. Bulletin du Bureau de Researches Geologicques et Minieres Bull BRGM 2nd Series, 3, 13-22.

  19. 19. Duijvenbooden, W.V. and Waegeningh, H.G. (1987) Vulnerability of Soil and Groundwater to Pollutants. Proceedings and Information No. 38 of International Conference, Noordwijk aan Zee, 30 March-3 April 1987.

  20. 20. Geo Wat (2011) Technical Report—Measurement of Sub-Surface Resistivity.

  21. 21. Al Dulaymi, A.S., Al-Heety, E.A.R. and Hussien, B.M. (2012) Geo-Electrical Investigation of Mullusi Aquifer, Rutba, Iraq. International Journal of Geosciences, 3, 549-564.

  22. 22. Selvam, S. (2016) 1D Geoelectrical Resistivity Survey for Groundwater Studies in Coastal Area: A Case Study from Pearl City, Tamil Nadu. Journal of the Geological Society of India, 87, 169-178.
    http://dx.doi.org/10.1007/s12594-016-0385-x

  23. 23. Anomohanran, O. (2013) Geophysical Investigation of Groundwater Potential in Ukelegbe, Nigeria. Journal of Applied Sciences, 13, 119-125.
    http://dx.doi.org/10.3923/jas.2013.119.125

  24. 24. Patil, S., Kachate, N., Marathe, N., Ingle, S. and Golekar, R. (2015) Electrical Resistivity Studies for Groundwater Exploration in the Some Parts Chopda Block of Jalgaon District, Maharashtra India. International Research Journal of Earth Sciences, 3, 2321-2527.

  25. 25. Hussain, Y., Dilawar, A., Ullah, S., Akhter, G., Martinez-Carvajal, H., Hussain, M. and Aslam, A. (2016) Modelling the Spatial Distribution of Arsenic in Water and Its Correlation with Public Health, Central Indus Basin, Pakistan. Journal of Geoscience and Environment Protection, 4, 18-25.
    http://dx.doi.org/10.4236/gep.2016.42003

  26. 26. Aller, L., Bennett, T., Lehr, J.H., Petty, R.J. and Hackett, G. (1987) DRASTIC: A Standardized System for Evaluating Groundwater Pollution Potential Using Hydrogeologic Settings. Environmental Protection Agency NWWA/EPA Series EPA-600/2-87-035, National Water Well Association, Dublin.

  27. 27. de Oliveira Braga, A.C. and Francisco, R.F. (2014) Natural Vulnerability Assessment to Contamination of Unconfined Aquifers by Longitudinal Conductance—(S) Method. Journal of Geography and Geology, 6, 68-79.

  28. 28. Braga, A.C.O. (2008) Estimation of the Natural Vulnerability of Aquifers: A Contribution from the Resistivity and Longitudinal Conductance. Brazilian Journal of Geophysics, 26, 61-68.

  29. 29. Obiora, D.N., Ajala, A.E. and Ibuot, J.C. (2015) Investigation of Groundwater Flow Potential in Makurdi, North Central Nigeria, Using Surficial Electrical Resistivity Method. African Journal of Environmental Science and Technology, 9, 723-733.
    http://dx.doi.org/10.5897/AJEST2015.1983

  30. 30. Hussain, Y. (2014) Hydrogeophysical Investigations and GIS Vulnerability Mapping of Kot Addu, M.Phil. Dissertation, Quaid-i-Azam University, Islamabad, Pakistan

  31. 31. Daniel, A., Louis, O., Emmanuel, C. and Orji, N. (2015) Delineation of Potential Groundwater Zones Using Geo- Electrical Sounding Data at Awka in Anambra State, South-Eastern Nigeria. European Journal of Biotechnology and Bioscience, 3, 1-5.

  32. 32. Fadele, S.I., et al. (2013) The Use of Vertical Electrical Sounding (VES) for Groundwater Exploration around Nigerian College of Aviation Technology (NCAT), Zaria, Kaduna State, Nigeria. The Pacific Journal of Science and Technology, 13, 549-555.

  33. 33. Atakpo, E.A. and Ayolabi, E.A. (2009) Evaluation of Aquifer Vulnerability and the Protective Capacity in Some Oil Producing Communities of Western Niger Delta. The Environmentalist, 29, 310-317.
    http://dx.doi.org/10.1007/s10669-008-9191-3

  34. 34. Kelly, W.E. and Stanislav, M. (1993) Applied Geophysics in Hydrogeological and Engineering Practice. Elsevier, Amsterdam, London, New York, Tokyo.

  35. 35. Mosuro, G.O., et al. (2016) Assessment of Groundwater Vulnerability to Leachate Infiltration Using Electrical Resistivity Method. Applied Water Science, 1-13.

  36. 36. Sonkamble, S. (2014) Electrical Resistivity and Hydrochemical Indicators Distinguishing Chemical Characteristics of Subsurface Pollution at Cuddalore Coast, Tamil Nadu. Journal Geological Society of India, 83, 535-548.
    http://dx.doi.org/10.1007/s12594-014-0081-7

  37. 37. Swartx, J.H. (1937) Resistivity Studies of Some Salt-Water Boundaries in the Hawaiian Islands. Transactions of American Geophysical Union, 18, 387- 393.
    http://dx.doi.org/10.1029/TR018i002p00387

  38. 38. Lima, O.A.L., Sato, H.K. and Porsani, M.J. (1995) Imaging Industrial Contaminant Plumes with Resistivity Techniques. Journal of Applied Geophysics, 34, 93-108.
    http://dx.doi.org/10.1016/0926-9851(95)00014-3

  39. 39. Sikandar, P., Bakhsh, A., Arshad, M. and Rana, T. (2010) The Use of Vertical Electrical Sounding Resistivity Method for the Location of Low Salinity Groundwater for Irrigation in Chaj and Rachna Doabs. Environmental Earth Sciences, 60, 1113-1129.
    http://dx.doi.org/10.1007/s12665-009-0255-6

  40. 40. Srinivasan, K., Poongothai, S. and Chidambaram, S. (2013) Identification of Groundwater Potential Zone by Using GIS and Electrical Resistivity Techniques in and around the Wellington Reservoir, Cuddalore District, Tamilnadu, India. European Scientific Journal, 9, 312-331.

  41. 41. Amos-Uhegbu, C., Igboekwe, M.U., Chukwu, G.U., Okengwu, K.O. and Eke, T.K. (2014) Geoelectrical Delineation and Geochemical 2 Characteristics of Aquifer Systems in Kwa-Ibo River 3 Watershed, Abia State, Nigeria. Journal of Scientific Research and Reports, 3, 818-843.
    http://dx.doi.org/10.9734/JSRR/2014/6215

  42. 42. Greenman, D.W., Swarzenski, W.V. and Bennett, G.D. (1967) Groundwater Hydrology of the Punjab, West Pakistan with Emphasis on Problems Caused by Canal Irrigation. US Geological Survey Water-Supply Paper.

  43. 43. Hussain, S.D., Sajjad, M.I., Akram, W., Ahmad, M., Rafiq, M. and Tariq, J.A. (1990) Surface Water/Groundwater Relationship in Chaj Doab. Pakistan Institute of Nuclear Science and Technology, Pakistan.

  44. 44. Plymale, C.L. and Angle, M.P. (2002) Groundwater Pollution Potential of Fulton County, Ohio. Ohio Department of Natural Resources Division of Water, Water Resources Section, Groundwater Pollution Potential, Report No. 45.

  45. 45. Al Hallaq, A.H. and Elaish, B.S.A. (2012) Assessment of Aquifer Vulnerability to Contamination in Khanyounis Governorate, Gaza Strip-Palestine, Using the DRASTIC Model within GIS Environment. Arabian Journal of Geosciences, 5, 833-847.
    http://dx.doi.org/10.1007/s12517-011-0284-9

  46. 46. Albinet, M. and Margat, J. (1970) Groundwater Pollution Vulnerability Mapping. Bulletin du Bureau de Researches Geologicques et Minieres Bull BRGM 2nd Series, 3, 13-22.

  47. 47. Duijvenbooden, W.V. and Waegeningh, H.G. (1987) Vulnerability of Soil and Groundwater to Pollutants. Proceedings and Information No. 38 of International Conference, Noordwijk aan Zee, 30 March-3 April 1987.

  48. 48. Geo Wat (2011) Technical Report—Measurement of Sub-Surface Resistivity.

  49. 49. Al Dulaymi, A.S., Al-Heety, E.A.R. and Hussien, B.M. (2012) Geo-Electrical Investigation of Mullusi Aquifer, Rutba, Iraq. International Journal of Geosciences, 3, 549-564.

  50. 50. Selvam, S. (2016) 1D Geoelectrical Resistivity Survey for Groundwater Studies in Coastal Area: A Case Study from Pearl City, Tamil Nadu. Journal of the Geological Society of India, 87, 169-178.
    http://dx.doi.org/10.1007/s12594-016-0385-x

  51. 51. Anomohanran, O. (2013) Geophysical Investigation of Groundwater Potential in Ukelegbe, Nigeria. Journal of Applied Sciences, 13, 119-125.
    http://dx.doi.org/10.3923/jas.2013.119.125

  52. 52. Patil, S., Kachate, N., Marathe, N., Ingle, S. and Golekar, R. (2015) Electrical Resistivity Studies for Groundwater Exploration in the Some Parts Chopda Block of Jalgaon District, Maharashtra India. International Research Journal of Earth Sciences, 3, 2321-2527.

  53. 53. Hussain, Y., Dilawar, A., Ullah, S., Akhter, G., Martinez-Carvajal, H., Hussain, M. and Aslam, A. (2016) Modelling the Spatial Distribution of Arsenic in Water and Its Correlation with Public Health, Central Indus Basin, Pakistan. Journal of Geoscience and Environment Protection, 4, 18-25.
    http://dx.doi.org/10.4236/gep.2016.42003

  54. 54. Aller, L., Bennett, T., Lehr, J.H., Petty, R.J. and Hackett, G. (1987) DRASTIC: A Standardized System for Evaluating Groundwater Pollution Potential Using Hydrogeologic Settings. Environmental Protection Agency NWWA/EPA Series EPA-600/2-87-035, National Water Well Association, Dublin.

  55. 55. de Oliveira Braga, A.C. and Francisco, R.F. (2014) Natural Vulnerability Assessment to Contamination of Unconfined Aquifers by Longitudinal Conductance—(S) Method. Journal of Geography and Geology, 6, 68-79.

  56. 56. Braga, A.C.O. (2008) Estimation of the Natural Vulnerability of Aquifers: A Contribution from the Resistivity and Longitudinal Conductance. Brazilian Journal of Geophysics, 26, 61-68.