48,06,2005misc 937 ANNALS OF GEOPHYSICS, VOL. 48, N. 6, December 2005 Key words electrical resistivity – sounding – ground- water condition – Shooro – Iran 1. Introduction The Shooro Basin that lies in Southeast Iran is located in Sistan and Baluchestan Province. The elevation of the basin is in the range of a minimum of 1360 m above the mean sea level in north to a maximum of 2960 m above the mean sea level in the southeast. The study area lies in the north of the basin about 80 km in the southwest of Zahedan city, the capital of Sistan and Baluchestan Province (29º00lN-29º12lN and 60º05lE-60º55l; fig. 1), and covers about 210 km2. A nearly flat plain with bedrock outcrops characterizes this area by a mild relief, with slops between 30 to 40 m topo- graphic contours in the south and central parts. The purpose of this paper is to use the resis- tivity data and interpreting geoelectrical sound- ings to study the aquifer conditions. The use of electrical resistivity for both groundwater re- source mapping and for water quality evalua- tions has increased dramatically over recent decades in much of the world due to the rapid ad- vances in microprocessors and associated nu- merical modeling solutions. The VES has proved very popular with groundwater studies due to simplicity of the technique. Water resources of the study area are threat- ened by increasing population trend with a re- sulting increase in water demand and the stress- Geoelectrical investigation for the assessment of groundwater conditions: a case study Gholam R. Lashkaripour (1) and Mohammad Nakhaei (2) (1) Department of Geology, University of Sistan and Baluchestan, Zahedan, Iran (2) Department of Geology, University of Tarbiat Moallem, Tehran, Iran Abstract An electrical resistivity survey involving Vertical Electrical Soundings (VES) was carried out in the Shooro Basin in Southeast Iran in order to study groundwater conditions such as depth, thickness and aquifer bound- aries. Vertical electrical soundings by Schlumberger array were conducted in this area. The resistivity Schlum- berger soundings which have a maximum current electrode spacing (AB) ranging from 200 m to 600 m were carried out at 207 positions in 19 profiles. Interpretation of these soundings indicates the presence of an alluvial aquifer. This aquifer is divided into eastern and western parts by the Shooro River, which comprises a variable thickness and resistivity of deposits. The average permeability coefficient and resistivity in the western part, es- pecially southwest is higher than the eastern part of the aquifer. Therefore, it seems that Shooro River follows a fault zone in the region. The high resistivity of west part is due to the water quality and the existence of alluvial fan with coarse grain materials. Low aquifer resistivities in the east are associated with finer materials and also brackish water infiltration from the adjacent basin mainly in the central part of the aquifer. Furthermore, zones with high yield potential have been determined in this research based on the resistivity data. Mailing address: Dr. G.R. Lashkaripour, Department of Geology, University of Sistan and Baluchestan, P.O. Box 161-98135, Khash Road, Zahedan, Iran; e-mail: La- shkarg@hamoon.usb.ac.ir 938 Gholam R. Lashkaripour and Mohammad Nakhaei es of water use for various activities. The need for this research is to study groundwater condi- tions and to protect groundwater supplies as a unique source of water for this area. 2. Geological background From the geological point of view, Shooro Basin is situated in the Flysch zone of Eastern Iran. The major portions of the Flysch zone con- sist of shale, sandstone and limestone. In the Fly- sch zone, sediments of older than Cretaceous age are absent (Tirrul et al., 1983; McCall, 1997). Approximately the total bedrock of the aquifer consists of shale and slats mainly green of poor permeability, in the age of Upper Creta- ceous and partly Eocene. In the west part of the basin, major outcrop consists of sandstone, shale and philit. Quaternary sediments cover the study Fig. 1. Location of the studied area. Geoelectrical investigation for the assessment of groundwater conditions: a case study area. The Quaternary deposits which cover the plain area are composed mainly of fine grain al- luvium of Shooro River, gravel to pebble of allu- vial fan, sand to gravel and rock debris of the nearby mountains (Geological Survey of Iran, 1995). 3. Hydrogeological conditions The study area is characterized by an arid climate with extremely low rainfall (average 84 mm/yr) and with high evaporation rate (Lash- karipour, 2000). The scanty rainfall is confined to the winter season and rain usually occurs as thunderstorms and showers. The fresh ground- water in the region is originated from infiltra- tion of rainwater, which has created the uncon- fined Shooro aquifer in the study area. The Shooro aquifer is the major exploited aquifer in this region and it is an alluvial deposit aquifer. This aquifer is divided into east and west parts by the Shooro River, which comprises a variable thickness of deposit. The average per- meability coefficient in the west part is higher than that in the east. Therefore, it seems that Shooro River follows a fault zone in the region. The aquifer is mainly recharged by the an- nual rainfall, and discharges naturally by under- flow to the Shooro River. It is also discharged artificially through a number of dug wells and galleries called qanat. The newly drilled wells in north and central parts of the aquifer are now used to supply water for agricultural purposes and drinking water to the villages. Rapid agri- cultural development in the study area has led to increased on demand for water supply. The monitoring of the groundwater level exhibits a decreasing trend of water level. The main rea- son responsible for this lowering of the ground- water table is that the wells pumping from groundwater resources have exceeded natural recharge in the recent years. The present study is based on data from the vertical electrical soundings that carried out dur- ing April to June 2001; lithological and hydroge- ological information from 13 wells; chemical analyses of groundwater samples from the avail- able boreholes. These data are used to evaluate the subsurface hydrogeological and structural conditions to a depth of about 150 m. In addition, estimation of the groundwater quality and recom- mendation for possible site-selections for drilling productive wells will be made. The aquifer hy- draulic parameters are shown in table I. 4. Geoelectrical resistivity survey Geoelectrical resistivity techniques are pop- ular and successful geophysical exploration for study groundwater conditions in the world. The resistivity of material depends on many factors such as groundwater, salinity, saturation, aquifer lithology and porosity. For example, the resis- tivity of an aquifer is related to Electrical Con- ductivity (EC) of its water. When the groundwa- ter EC is high, the resistivity of the aquifer could reach the same range as a clayey medium and the resistivity parameter is no longer useful to determine the aquifer (Vouillamoz et al., 2002). However, this method has been carried out suc- cessfully for exploration of groundwater. This technique is widely used to determine depth and nature of an alluvium, boundaries and location of an aquifer. The resistivity method was used to solve more problems of groundwater in the types allu- vium, karstic and another hard formation aquifer as an inexpensive and useful method. Some uses of this method in groundwater are: determination of depth, thickness and boundary of an aquifer (Zohdy, 1969; and Young et al. 1998), determi- nation of interface saline water and fresh water (El-Waheidi, 1992; Yechieli, 2000; and Choud- hury et al., 2001), porosity of aquifer (Jackson et al., 1978), water content in aquifer (Kessels 939 Table I. Hydraulic parameters of the aquifer. No. Parameters East part West part 1 Aquifer thickness (m) 24 40 2 Average groundwater 23 41 depth from ground level (m) 3 Aquifer resistivity (µm) 12 45 4 EC range (µs cm−1) 2800-9600 760-2500 5 Storage coefficient (%) 25 10 940 Gholam R. Lashkaripour and Mohammad Nakhaei et al., 1985), hydraulic conductivity of aquifer (Yadav and Abolfazli, 1998; Troisi et al. 2000), transmissivity of aquifer (Kosinski and Kelly, 1981), specific yield of aquifer (Frohlich and Kelly, 1987), contamination of groundwater (Kelly, 1976; and Kaya, 2001). Contamination usually reduces the electrical resistivity of pore water due to increase of the ion concentration (Frohlich and Urish, 2002). However, when re- sistivity methods are used, limitation can be ex- pected if ground inhomogeneties and anisotropy are present (Matias, 2002). This research project deals with detecting the aquifer conditions and groundwater quality in the study area. The resistivity survey was completed with 207 Vertical Electrical Sound- ings (VES) by the Schlumberger array in 19 profiles, with a maximum current electrode spacing (AB) ranging from 200 to 600 m. The profile spacing was 1 km and sounding spacing was 750 m. Five of these VES have been con- ducted adjacent to the dug wells for subsequent calibration process. Position and extension of all the VES’s are indicated in fig. 2. Twenty nine geoelectrical sections were drawn along the profiles. The apparent resistivity measure- ments were made with a KD Sound 7800B Ter- rameter Unit. This equipment is light and pow- erful for deep penetration. 5. Discussion The VES curves were obtained by plotting the apparent resistivity against electrode spacing. The resistivity for each of the vertical electrical sounding was drawn on transparent double log graph paper and a smooth field curve. Computer programs for reducing geoelectri- cal sounding curves into values of thickness and Fig. 2. The location of VES and aal profile. Geoelectrical investigation for the assessment of groundwater conditions: a case study resistivity of individual layers are described by Koefoed (1979) and Zohdy and Bisdorf (1989). The field curves were interpreted by the well- known method of curve matching with the aid of Russian software IPI7.63. However, the thick- ness and characteristics of the aquifer are fairly well known due to the number of wells dug in the center of the aquifer. A number of geoelectric sta- tions were purposely located near about 13 wells. The key to success of any geophysical survey is the calibration of the geophysical data with hy- drogeological and geological ground truth infor- mation. VES success must rely on the careful inter- pretation and integration of the results with the other geologic and hydrogeologic data for the site. Therefore, lithologic information obtained from log could be used to calibrate the VES field curves. Where test hole-log information was available, the solution to automatic interpretation procedure was constrained by keeping known layer thicknesses constant during the program computations. Lastly, the results of the Schlumberger elec- trical soundings were compared with the geo- logical sections obtained from 13 pizometer wells. These results are in a good agreement with the geological sections. Measuring the current penetration depth was calculated by a correlation between depth and AB/2 as shown in fig. 3. From this figure eq. (5.1) was obtained to calcutate depth of current penetration (5.1) From eq. (5.1), the current penetration depth is equal to approximately AB/4. The boundary of the aquifer, thickness and resistivity of subsurface layers was also deter- mined by the electrical survey in this research. From the interpretation of the resistivity curves a four-layer resistivities and thicknesses indi- cated four subsurface layers. These layers con- sisting of surface layer (topsoil), alluvium, sat- urated layer, and bedrock. In some cases more than one layer was evident in the saturated zone but these cases were also treated in this analy- sis as single layers. Depth and thickness of sub- surface layers were identified and dimension of the aquifer and type of bedrock were also indi- cated. Bedrock of the area is generally slat but in some parts has appeared as shale. For exam- ple, the geoelectrical section of profile aal in the central west of aquifer is shown in fig. 4. The aquifer consists of two main east and west parts that is separated along the Shooro River. The thickness of the aquifer and depth of bedrock is different in the two parts. Average resistivity of the thin top layer, alluvium, aquifer and bedrock calculated in the east as 70, 74, 12 and 113 Ω-m and in the west as 175, 116, 46 and 106 Ω-m respectively. These data indi- cate that bedrock in both sides of the aquifer is the same, but permeability and storage coeffi- cient in the eastern part is higher. Depth and . /Depth AB0 57 2 .0 97= .] g 941 Fig. 3. Correlation between depth of current penetration and length of current electrodes. 942 Fig. 4. Geoelectrical section of profile aal. Fig. 5. Subsurface layers in the eastern and western parts of the aquifer. T – thickness; R – resistivity. thickness of the aquifer were measured in the east as 23 and 24 m and in the west as 40 and 41 m respectively (fig. 5). The higher resistivi- ty in the west part is due to the existence of an alluvial fan that consists of a mixture of gravel, sand and silt. The lower resistivity in the east- ern part is due to intering of brackish water (EC in the range of 2800 to 9600 µs cm−1) from an Gholam R. Lashkaripour and Mohammad Nakhaei 943 Geoelectrical investigation for the assessment of groundwater conditions: a case study adjacent basin and finer materials of sand and silt mixed with clay. Yield potential in the western part of the aquifer is more than the eastern part. RT map (Aquifer Resistivity multiply by Aquifer Thick- ness) that shows the potential of the aquifer for water resources and quality in the whole study area is shown in fig. 6. The best part of the aquifer for future development is the central portion of the west part of the aquifer. In this part, the best water quality and quantity of groundwater sources are found with respect to high thickness and resistivity. Furthermore, the boundary of the aquifer and type of the bedrock were indicated based on the geoelectrical data. 6. Conclusions This study reveals that surface electrical measurements are effective to study groundwa- ter conditions. Based on the interpretation of geoelectrical data, the following conclusions are drawn: 1) The alluvial aquifer mainly consists of gravel and sand in the west part and sand and silt mixed with clay materials in the east part. The re- sistivity of the aquifer increases towards the west due to decreasing salinity of water and/or clay content. 2) VES tests revealed four subsurface geo- electric layers; thin top layer, the alluvium, the aquifer and the bedrock respectively. 3) The aquifer thickness increases towards the west, the regional direction of increasing deposi- tion in the basin. The average thickness of the sat- urated alluvial aquifer in the west part has been estimated about 40 m and in the east about 24 m. 4) The bedrock of the aquifer shows different resistivity values due to type of bedrock and re- spect to degree of saturation and values of frac- tures. 5) The boundary of the aquifer has been esti- mated and zones with high yield potential have Fig. 6. RT map of the aquifer. 944 Gholam R. Lashkaripour and Mohammad Nakhaei been determined for future development in the basin and for choosing drilling sites. 6) The resistivity data of the aquifer shows a decrease in the eastern part due to the intru- sion of low quality (saline) water from the ad- jacent basin. Acknowledgements The research presented in this paper has been supported by a grant from the Sistan Baluchestan Province Water Organization of the Ministry of Power, for which we express our sincere thanks. REFERENCES APPARAO, A. and T.G. RAO (1974) Depth of investigation in resistivity methods using linear electrodes, Geophys. Prospect., 22, 211-223. CHOUDHURY, K., D.K. SAHA and P. CHAKRABORTY (2001): Geophysical study for saline water intrusion in a coastal alluvial terrain, J. Appl. Geophys., 46, 189-200. DARVISHZADEH, A. (1981): Geology of Iran (Nasher Danesh Emrooz Publisher), (in Persian). EL-WAHEIDI, M.M., F. MERLANTI and M. PAVAN (1992): Geoelectrical resistivity survey of the central part of Azraq Basin (Jordan) for identifying saltwater/fresh- water interface, J. Appl. Geophys., 29, 125-133. FROHLICH, R.K. and W.E. KELLY (1987): Estimates of spe- cific yield with the geoelectric resistivity method in glacial aquifers, J. Hydrol., 97, 33-44. FROHLICH, R.K. and D. URISH (2002): The use of geo- electrics and test wells for the assessment of groundwa- ter quality of a coastal industrial site, J. Appl. Geo- phys., 50, 261-278. GEOLOGICAL SURVEY OF IRAN (1995): Geological Map of Khash, Scale 1:250 000. JACKSON, P.N., S.D. TAYLOR and P.N. STANFORD (1978): Re- sistivity- porosity-particle shape relationships for ma- rine sands, Geophysics, 43, 1250-1268. KAYA, G.K. (2001): Investigation of groundwater contami- nation using electric and electromagnetic methods at an open waste-disposal site: a case study from Isparta, Turkey, Environ. Geol., 40, 725-731. KELLY, E.W. (1976): Geoelectric sounding for delineating ground water contamination, Ground Water, 14, 6-11. KESSELS, W., I. FLENTGE and H. KOLDITZ (1985): DC geo- electric sounding to determine water content in the salt mine asse (FRG), Geophys. Prospect., 33, 456-446. KOEFOED, O. (1979): Geosounding Principles, 1. Resistivi- ty Sounding Measurements (New York, NY, Elsevier Scientific Pub. Co.), pp. 276. KOSSINSKI, W.K. and W.E. KELLY (1981): Geoelectric sounding for predicting Aquifer Properties, Ground Water, 19, 163-171. LASHKARIPOUR, G.R. (2000): Groundwater pollution of Za- hedan city in the East of Iran, J. Nepal Geol. Soc., 21, 99-102. MATIAS, M.J.S. (2002): Squary array anisotropy measure- ments and resistivity sounding interpretation, J. Appl. Geophys., 49, 185-194. MCCALL, G.J.H. (1997): The geotectonic history of Makran and adjacent area of Southern Iran, J. Asian Earth Sci., 15, 517-531. TIRRUL, R., I.R. BELL, R.J. GRIFFIS and V.E. CAMP (1983): The Sistan suture zone of Eastern Iran, Geol. Soc. Am. Bull., 94, 134-150. TROISI, S., C. FALLICOS, S. STRAFACE and E. MIGLIARI (2000): Application of kriging with external drift to es- timate hydraulic conductivity from electrical resistive- ly data in unconsolidated deposits near Montato Uffu- go, Italy, Hydrogeol. J., 8, 356-367. VINCENZ, S.A. (1968): Resistivity investigations of lime- stone aquifers in Jamaica, Geophysics, 33, 980-994. VOUILLAMOZ, J.M., M. DESCLOITRES, J. BERNARD, P. FOUR- CASSIER and L. ROMAGNY (2002): Application of inte- grated magnetic resonance sounding and resistivity methods for borehole implementation. A case study in Cambodia, J. Appl. Geophys., 50, 67-81. YADAV, G.S. and H. ABOLFAZLI (1998): Geoelectric sound- ings and their relationship to hydraulic parameters in semiarid regions of Jalore, Northwestern India, J. Ap- pl. Geophys., 39, 35-51. YECHIELI, Y. (2000): Fresh-saline ground water interface in the Western Dead Sea area, Ground Water, 38, 615-623. YOUNG, M.E., R.G.M. DE BRUIJIN and A. SALIM AL-ISMAI- LY (1998): Reports: exploration of an alluvial aquifer in Oman by time-domain electromagnetic sounding, Hy- drogeol. J., 6, 383-393. ZOHDY, A.A.R. (1969): Application of deep electrical sound- ings for groundwater exploration in Hawaii, Geophysics, 34, 584-600. ZOHDY, A.A.R. and R.J. BISDORF (1989): Programs for the automatic processing and interpretation of Schlum- berger sounding curves in Quick Basic, U.S. Geol. Surv. Open File Rep. 89-137-2, p. 64. (received January 24, 2005; accepted November 28, 2005)