Available online a t www.scholarsresearchlibrary.com

Scholars Research Library

Archives of Applied Science Research, 2011, 3 (2):414-421

(http://scholarsresearchlibrary.com/archive.html)

ISSN 0975-508X CODEN (USA) AASRC9

Identification of aquifer zones by VES method: A case study from Mangalore block, , S.

Ahilan J. and G. R. Senthil Kumar

Department of Earth Sciences, Annamalai University, , India ______

ABSTRACT

Aquifer zones are identified by the Geophysical Electrical Resistivity Survey method. This method is useful to delineate the subsurface formations, weathered zone, fracture pattern, etc. An attempt made to identify the subsurface lithology and aquifer zones by geoelectrical resistivity method in part of Mangalur block, District, Tamil Nadu, S.India. The study area consist of charnockite and charnockitic gneiss rocks. Electrical soundings have been conducted at every 2 Km intervals in grid pattern. Totally 36 VES were conducted, which have been qualitatively and quantitatively interpreted using software packages. Pseudosections have been generated by IPI 2 WIN ver.3.1 and isoresistivity maps were prepared. Geoelectrical parameter and sections have been prepared based on quantitative interpretation. Considering the hydrometeorological and hydrogeological conditions the VES interpretation was done. From the interpretation result the VES no. 21 (Alambadi) is a prospective zone for groundwater target.

Key words : Resistivity survey, quantitative and qualitative, interpretation, isoresistivity maps. ______

INTRODUCTION

Water is renewable resource occurs in three forms liquid, solid and gaseous. We known that the groundwater is essential for irrigation, industry and domestic purpose. Groundwater is the main source for potable water supply, domestic, industrial and agricultural uses. The scarcity of groundwater increases day by day due to rapid population, urbanization, industrial and agricultural related activities, natural calamities, and etc. the impact of trio on soil and groundwater is alarming with years of devastating effects on humans and the ecosystem [1]. Study of groundwater geology is much useful for all the activities of human life. Groundwater is more advantageous than the surface water. Water scarcity problem affects the human chain and other living things. To meet out the demand of water, people are depending more on aquifers. There are two end members in spectrum of types of aquifers; confined and unconfined (with

414 Scholar Research Library Ahilan J. et al Arch. Appl. Sci. Res., 2011, 3 (2):414-421 ______semi confined aquifer being in between them) [3]. For identifying the ground water potential in the hard rock aquifer terrain the main target is fractured zone. The present study focuses on identification of fracture zone and its thickness by using geoelectrical resistivity method.

Location of the study area The present study area has falls in Mangalore block of , Tamil Nadu, S.India (Location map shown in Fig.1). The study area lies between the latitudes North 11 °21’80” to 11 °30’11” and the longitudes East 78 °40 ’57” to 79 °03 ’11” in the Survey of India Toposheet no. M M 58 &58 . The present study area occupies an areal extent of 100 sq.km and the relief ranges 3 15 from 62 m to 121 m above MSL. The study area receives an average rainfall of 1100 mm with more than 80% of the rainfall received during the NE monsoon. The minimum and maximum temperature ranges between 20°C and 34°C in the month of January and May respectively. River Vellar flows in the southern part of the study area. Geomorphologically the area consists of old flood plains, pediments, duricrust and pediments covered by forest land [7].

Fig.1. Location map of the study area

Geology and Hydrogeology The study area rock types belong to early to mid Precambrian represented by charnockite and charnockitic gneiss, indicating the oldest and subjected to granulite facies of metamorphism. The charnockites are intermediate to acid in composition, coarse to medium grained and form the high land topography. The charnockitic rocks are massive to foliate and the foliations usually

415 Scholar Research Library

Ahilan J. et al Arch. Appl. Sci. Res., 2011, 3 (2):414-421 ______trending ENE – WSW with an average dip of 45° towards South. The charnockite shows different depth of weathered zone. In the study area groundwater occurs under water table conditions in the joints, fractures and weathered rocks. Generally, the charnockite of the study area is highly massive and compact and devoid of joints and fractures making it impervious, which in turn result in poor potential. The open wells give better yield than bore wells [8].

MATERIALS AND METHODS

Geophysical prospecting of groundwater comes under both surface and subsurface exploration. Under geophysical prospecting, one of the electrical methods is schlumberger array of electrical resistivity method. The schlumberger array was used to ensure deep penetration and for logistics of limited man power in the field [2]. Schlumberger configuration using Microprocessor based signal stacking digital resistivity meter of IGIS, Hyderabad make, Model – SSR-MP-AT-S. Both the survey procedures resistivity profiling and resistivity sounding (VES) have been carried out. Resistivity profiling has been conducted in a grid pattern. The total study area of 100 Sq.km. has been divided into square grids of 2.0 Sq.km and 6 resistivity profiles were conducted with a station interval of 2 Km. resistivity profiling with AB 10, 20, 30, 40 & 50 m have been carried 2 out. VES has been conducted at 36 locations with AB 50 m. The resistivity data have been 2 qualitatively and quantitatively interpreted and analyzed by software packages.

Resistivity Method The electrical resistivity is the resistance offered by the opposite faces of a unit cube of material to direct current is called as resistivity. In geophysical literature the unit of resistivity is taken as the ohm. m. The resistance (R) of the material having a resistivity ( ρ) over a length (L) and surface area of current flow (A) is given by L R = ρ A This is governed by ohms law. The inverse of resistance is termed as conductance [5]. The resistivity of the geological formation is generally very high under dry conditions and decreases in clayey rock. The presence of water containing salt even in minor amounts, geological formations makes them relatively conductive and as the moisture increases the resistivity falls deeply. As the salinity of water increases the resistivity of the rock formation decreases substanceously.

Measurement of Resistivity In general for measuring the resistivity of the sub surface formation four electrodes are required. The current of electrical intensity (I) is introduced between one pair of electrodes called current electrodes which can be identified as A & B. The potential difference produced as a result of current flow is measured with help of another pair of electrodes called potential electrodes represented as M & N. Let ∆r represent the potential difference. The apparent resistivity ∆v measure is K × , where K represents geometrical constant, which can be calculated if we I know the electrode arrangements. The basic needs for the resistivity survey are the power source, meter to measure current and potential, electrodes and cables. 416 Scholar Research Library

Ahilan J. et al Arch. Appl. Sci. Res., 2011, 3 (2):414-421 ______Schlumberger Configuration In schlumberger configuration all the four electrodes are kept in a line similar to that of Wenner but the outer electrode spacing is kept large compared to inner electrode spacing usually more than five times. For each measurement only the current electrodes are moved keeping the potential electrodes at the same locations. The potential electrodes are moved only when the signal become too week to be measured. The apparent resistivity for this configuration is computed with the formula;

 AB 2  MN 2  Π  −   R  2   2   ρa = MN

Vertical Electrical Sounding (VES) A series of measurements of resistivity are made by increasing the electrode spacing in successive steps about a fixed point. This method of vertical exploration is known as the expanding electrode method, “Resistivity sounding” or “Depth probing” or vertical electrical sounding (VES). The apparent resistivity values obtained with increasing values of electrode separation are used to estimate the thickness and resistivities of the subsurface formations. VES mainly employed in groundwater exploration to determine the disposition of the aquifers .

Geoelectrical Parameters The main objective of the quantitative interpretation of VES curve is to obtain the geoelectrical parameters. A geoelectric layer is called by its fundamental characters, resistivity ‘ ρ’ and thickness ‘h’. These two parameters are called the primary geoelectrical parameters. The secondary geoelctrical are also important to describe the geoelectric section consisting of several layers. The longitudinal (S) and the total transverse (T), transverse resistivity (T/h) an – isotropy (λ) are called the secondary geoelectrical parameters. The (Fig.2) Shows columnar prism used in defining geoelectric parameters.

P = RESISTIVITY h = THICKNESS

S = LONGITUDINAL CONDUCTANCE

T = TRANSVERSE RESISTANCE

Fig.2. Columnar prism defining geoelectrical parameters

417 Scholar Research Library

Ahilan J. et al Arch. Appl. Sci. Res., 2011, 3 (2):414-421 ______Fig.3.Pseudosections of the profiles I to VI

418 Scholar Research Library

Ahilan J. et al Arch. Appl. Sci. Res., 2011, 3 (2):414-421 ______

ISO RESISTIVITY MAP FOR 6 M ISO RESISTIVITY MAP FOR 2 M.

73.50 64.40 91.00 45.00 14.00 25.00 40.00 77.00 62.00 69.00 61.00 24.00 46.00 40.00

AR IN OHM /M 35.00 AR IN OHM/M 35.00 13.00 38.00 23.10 98.00 45.60 8.70 35.00 140.00 31.00 61.00 16.00 92.00 75.00 21.00 150.00 130.00 30.00 140.00 120.00 30.00 130.00 110.00 120.00 25.0015.00 43.00 38.00 27.00 74.00 7.40 100.00 110.00 25.0030.00 76.00 48.00 63.00 68.00 14.00 90.00 100.00 90.00 80.00 80.00 20.00 70.00 20.00 70.00 60.00 103.00 43.50 8.30 55.50 6.60 25.00 75.00 51.00 17.00 120.00 12.00 35.00 60.00 50.00 50.00 15.00 15.00 40.00 40.00 30.00 30.00 10.00 20.00 10.00 20.00 91.00 98.00 130.00 23.00 12.00 12.00 10.00 30.00 49.00 133.00 13.00 5.00 5.50 10.00 0.00 0.00 5.00 5.00

46.00 5.00 3.50 17.00 3.00 37.00 0.00 21.00 3.40 3.00 32.00 3.40 14.00 0.00 5.00 10.00 15.00 20.00 25.00 30.00 35.00 40.00 0.00 0.00 5.00 10.00 15.00 20.00 25.00 30.00 35.00 40.00 LEGEND

LEGEND 0 TO 50 OHM / M SCALE 1CM = 250M 0 TO 50 OHM /M 50 TO 100 OHM / M

SCALE 1CM = 250M 50 TO 100 OHM /M 100 TO 150 OHM / M

VES LOCATIONS 100 TO 140 OHM /M VES LOCATIONS

ISO RESISTIVITY MAP FOR 10 M ISO RESISTIVITY MAP FOR 15 M

87.00 64.00 60.00 64.00 44.00 70.00 70.00 61.00 78.00 56.00 71.00 74.00 40.00 40.00

35.00 AR IN OHM/M 35.00 AR IN OHM/M

44.00 69.00 121.00 86.00 74.00 23.00 150.00 65.00 101.00 258.00 106.00 96.00 35.00 260.00 140.00 30.00 30.00 130.00 240.00 120.00 220.00 110.00 200.00 25.0035.00 61.00 61.00 79.00 66.00 30.00 25.0058.50 94.00 61.00 79.00 81.00 35.00 100.00 180.00 90.00 160.00 20.00 80.00 20.00 20.00 140.00 70.00 120.00 78.00 60.00 19.00 130.00 18.00 48.00 60.00 72.00 81.00 26.00 145.00 26.00 63.00 100.00 15.00 50.00 15.00 40.00 80.00 30.00 60.00 10.00 20.00 40.00 135.00 133.00 85.00 31.00 18.00 18.00 10.00 10.00 235.00 41.00 72.00 39.00 28.00 25.00 20.00 0.00 0.00 5.00 5.00

71.00 6.00 5.00 24.00 5.00 56.00 0.00 0.00 5.00 10.00 15.00 20.00 25.00 30.00 35.00 40.00 85.00 10.00 8.00 26.00 7.00 64.00 0.00 5.00 10.00 15.00 20.00 25.00 30.00 35.00 40.00 0.00 0.00 5.00 10.00 15.00 20.00 25.00 30.00 35.00 40.00 LEGEND SCALE 1 CM = 250M LEGEND 0 TO 50 OHM / M 0 TO 100 OHM / M 50 TO 100 OHM / M SCALE : 1 CM = 250 M 100 TO 200 OHM / M 100 TO 150 OHM / M 200 TO 260 OHM / M

VES LOCATIONS VES LOCATIONS ISO RESISTIVITY MAP FOR 50M

ISO RESISTIVITY MAP FOR 20 M

144.00 126.00 147.00 77.00 273.00 184.00 65.00 73.00 96.00 52.00 85.00 89.00 40.00 40.00 40.00

AR IN OHM/M 35.00 35.00 AR IN OHM/M 35.00 178.00 278.00 840.00 140.00 301.00 92.00 85.00 132.00 335.00 115.00 126.00 43.00 820.00 340.00 780.00 30.00 30.00 320.00 30.00 740.00 300.00 700.00 280.00 660.00 260.00 620.00 25.00 25.00132.00 193.00 93.00 146.00 180.00 78.00 25.0072.50 116.00 68.00 89.00 99.00 44.00 240.00 580.00 540.00 220.00 500.00 200.00 460.00 20.00 20.00 180.00 420.00 160.00 380.00 79.00 99.00 31.00 148.00 32.00 78.00 140.00 112.00 175.00 76.00 228.00 65.00 134.00 340.00 300.00 120.00 15.00 15.00 260.00 100.00 220.00 80.00 180.00 60.00 10.00 140.00 10.00 40.00 100.00 238.00 45.00 62.00 48.00 35.00 31.50 200.00 81.00 87.00 95.00 72.00 74.00 20.00 60.00 0.00 20.00 5.00 5.00 5.00

82.00 13.50 11.00 34.50 10.00 79.00 197.00 33.00 28.00 75.00 32.40 159.00 0.00 0.00 0.00 5.00 10.00 15.00 20.00 25.00 30.00 35.00 40.00 0.00 5.00 10.00 15.00 20.00 25.00 30.00 35.00 40.00 0.00 5.00 10.00 15.00 20.00 25.00 30.00 35.00 40.00 LEGEND LEGEND 0 TO 300 SCALE 1CM = 250M 0 TO 150 SCALE 1CM = 250M 300 TO 600 150 TO 250

600 TO 850 250 TO 350

VES LOCATION VES LOCATION

Fig.4.Isoresistivity maps for the depth of 2, 6, 10, 15, 20 and 50 m

Quality and Quantitative method In the qualitative method the shape of the field curve is critically observed to get an idea qualitatively about the number of layers and the order of resistivities. The thirty six VES curves have been prepared by IPI2 Win software Program and studied in-depth VES curves of the study

419 Scholar Research Library

Ahilan J. et al Arch. Appl. Sci. Res., 2011, 3 (2):414-421 ______area. Majority of the curve types are ‘A’ and ‘AA’ types of the study area, the other curve types are ‘H’, combination of ‘HA’ and ‘KH’ types.

In the quantitative method geoelectrical parameter i.e., true resistivities and layer thickness are obtained. There are two types of quantitative analysis i) indirect method and ii) direct method. The field curve is compared with album of theoretical curves to get the geoelectrical parameters and this technique is known as curve matching technique.

To obtain initial model parameters by indirect method the field curve is matched with standard curves and also by auxiliary point technique [4]. In the direct method the inverse slope technique which is a semi empirical method proposed by [6].

NANGUR, VES 8, BGL ALAMBADI, VES 21, BGL

1 1 0 0 0.29 FIRST LAYER - 217 Ω M 0.29 FIRST LAYER - 25.5 Ω M

0.51 SECOND LAYER - 22.3 Ω M 0.64 SECOND LAYER - 4.35 Ω M

1 1 THIRD LAYER - 408 Ω M 0.9 THIRD LAYER - 188 Ω M 0.97

FOURTH LAYER - 26.1 Ω M 2 2 FOURTH LAYER - 9.65 Ω M

FOUERTH LAYER FOURTH LAYER THIRD LAYER THIRD LAYER 3 3 SECOND LAYER SECOND LAYER 3.08 FIRST LAER FIRST LAYER 3.15 THICKNESS IN M THICKNESS IN M

4 4

5 5

6 6

Fig.5. Geoelectrical sections

RESULT AND DISCUSION

The pseudosection profiles (Fig.3.) shows that the southern part of the profiles (station 25 & 31) having high resistivity indicating poor weathering. The low resistivity anomaly extends down to a depth of 31 m. Another low resistivity anomaly is not prominently observed at the southern part i.e. at station 26. Along the station 27 a prominent high resistivity anomaly is observed from the depth of 7m, indicating the massive nature of the rock (Fig.4 & 5). Another low resistivity anomaly is exhibited at station.23, where the low anomaly is sandwitched. Another lower resistivity anomaly is observed between the stations. 23 and 29, the weathering thickness is restricted upto a depth of 5m. There are two high resistivity anomalies observed at station.06 (7.5m) and 18 (7.5m).

Isoresistivity contour map showing the high resistivity anomalies at station the prominent low resistivity anomalies observed at stations 21, 15 and 09 indicating a weaker zone along N – S direction (Fig.4- 6m map). The prominent low resistivity anomaly exhibited in 2 m and 6 m contour maps indicating a prominent weaker zone. From the low resistivity anomalies at a station 21, 35, 36, 24, 12, 13 and 05 indicating a NW – SE oriented weaker zone, the high resistivity anomalies at a station 25, 26, 22, 09, 01 and 06 indicates a massive rock formation (Fig. 4- 10m map). The prominent low resistivity anomalies at stations are 21, 36, 12 and 18 indicating a weaker zone along N – S direction. High resistivity anomalies at a station 09, 22, 420 Scholar Research Library

Ahilan J. et al Arch. Appl. Sci. Res., 2011, 3 (2):414-421 ______and 25 indicate a poorly weathered rock (Fig.4- 15m map). The prominent low resistivity anomalies at stations 21, 22, 36, 12, 18 and 05 indicate N – S trending weaker zone and E - W trending weaker zones along station. 21 and 22 (Fig.4- 20m map). Figure 4 i.e. 50m map shows the prominent low resistivity anomalies at stations 36, 22, 18, 12, 11, 05, and 04 indicting a fractured zone along N – S direction.

VES.08 is ‘KH’ type curve and VES.21 is ‘AA’ type curve both are five layered nature, through VES. 21 have very least resistivity, indicating more weathered nature of formation (Fig.5).

CONCLUSION

Electrical resistivity survey has been carried out at every 2Km interval in a grid pattern and totally 36 VES were conducted. Pseudosections and isoresistivity maps were generated. Both low and high resistivity anomalies have been demarcated in pseudosection and isoresistivity map. Majority of the VES curves are ‘A’ and ‘AA’ types. With multi layered geoelectrical sections. From the pseudosections it may be revealed that the N, NW and NE are comparatively more weathered than the other parts. A prominent low resistivity anomaly has been exhibited at VES.21 in all the six maps indicating more weathered nature of the formation. Based on the hydrogeological condition and geoelectrical resistivity survey findings VES.21 (Alambadi) may be the prospective zone for groundwater development.

REFERENCES

[1] Ehirim C. N. and C. N. Nwankwo. Archives of Physics Research , 2010 , (2): 396-403. [2] Emmanuel C. Okolie. (2010), Archives of Physics Research , 2010 , 1 (3):62-71. [3] Ogungbe1 A.S, J.A Olowofela2, O.O Oresanya1, A.A Alabi., Archives of Physics Research, (2010) , 2 (2): 24-34. [4] Orellana E, and H. M. Mooney; Master tables and curves for electrical sounding over layered structure interciencia, (Madrid, 1966 ). [5] Ramanuja chary K. R., and T. Gangadhra Rao; A monograph on groundwater exploration by resistivity, S.P and well logging techniques, (IGIS, 2005) Hyderabad. [6]Sanker Narayana P. K. and K. R. Ramanuja Chary; Short notes – An inverse slope method of determining absolute resistivity, (Geophysics, 1967 ) XXXII (6). [7] Senthil Kumar, G.R; Hydrogeological investigations in Tittagudi Taluk, Cuddalore, Tamil Nadu, S. India, (Unpublished thesis of Annamalai University, 2006 ) pp. 23-27. [8] TWAD; Hand book for TWAD Hydrogeologist, (TWAD Board, 2002 ) Chennai.

421 Scholar Research Library