J. Nig. Soc. Phys. Sci. 1 (2019) 1–5 Journal of the Nigerian Society of Physical Sciences Original Research The Delineation of Groundwater and Geotechnical Parameters within Marmara Area of Chikun Local Government of Kaduna State, Nigeria N. Kurea,∗, I. H Daniela, C. G. Afuwaia, E. J. Adoyia, I. A. Bellob aDepartment of Physics, Kaduna State University, Kaduna, Nigeria bDivision of Agricultural Colleges, Ahmadu Bello University, Zaria, Nigeria Abstract Geophysical investigation for engineering studies was carried out around Faith Academy, Marmara village, Kaduna within the basement complex of North-Western Nigeria. Six (6) Vertical Electrical Sounding (VES) were established using Schlumberger configuration. The geoelectric section revealed five subsurface layers defined by the topsoil, which comprises humus, clay and laterite; followed by weathered layer which comprises of sandy soil, fractured layer which constitutes of coarse grain sands and gravels and fresh basement which is porphyritic granite. The resistivities and thicknesses of the topsoil range from 47 − 4212Ωm and 2 − 5m respectively. The weathered/fractured basement has average thickness of 42m with resistivity ranging from 350 − 774Ωm. The Areas found to meet the conditions of high basement resistivity and shallow depth to basement, are found to be at VES A1, A2, A3, A4 and A6. Hence, VES A1 have been found to be the competent zones within the study area, and are good for the construction of high-rise buildings, roads and bridges. VES A3, A4 and A6 are areas that do not meet the conditions for construction of a high-rise buildings. However, VES A3, A4 and A6 appears to be viable for groundwater exploration with aquifer thickness ranges between 18 − 42m and depth ranges of 6 − 10m, while the average depth to fresh basement is 43m. Keywords: Geotechnical, Resistivity, Chikun, Topsoil/Laterite, Aquifer Article History : Received: 13-03-2019 Received in revised form: 26-03-2019 Accepted for publication: 27-03-2019 Published: 30-03-2019 c©2019 Journal of the Nigerian Society of Physical Sciences All Rights Reserved. Communicated by: O. J. Abimbola 1. Introduction The development of technology has made the quest for wa- ter for all purpose in life to drift from ordinary search to prospect- ing to steady and reliable subsurface or groundwater from bore- holes. In Nigeria, presently, borehole has rescued the citizenry from acute shortage of water. The geoelectric method has been found to be very reliable for groundwater studies over the years. The geotechnical evalu- ation of subsurface condition of a site is necessary in generating ∗Corresponding Author Tel. No: 09036041268 Email address: nicodemus.kure@kasu.edu.ng (N. Kure ) relevant data inputs for the design, and construction of struc- tures and boreholes [1]. The search for groundwater dates back to history. This is because, of all the natural resources, water permeates perhaps must deeply into the all aspect of our life. Thus, water is essential to life. Its presence or lack of it deter- mines, to great extent, the nature of the natural environment in which we live and majority of our economic activities depend on it. The development and utilization of groundwater resource therefore becomes very imperative [2, 3]. Due to the collapse of bridges, buildings and other struc- tures in Nigeria, there are needs to provide information on the subsurface sequence and structural disposition necessary for the 1 Kure et al. / J. Nig. Soc. Phys. Sci. 1 (2019) 1–5 2 buildings. Detailed geophysical and geotechnical investigation provides such information to provide adequate solution to pre- vention against future occurrence. Furthermore, lack of proper soil investigation and interpretation of sites conditions amongst the causes that may pose danger to the existing structures in Nigeria. The importance of geological exploration methods, as means of exploring the groundwater and zones for engineering struc- tures such as high rising building and bridges cannot be over emphasized. Methods including Vertical Resistivity Sounding (VES), Seismic Refraction, Horizontal Profiling (Wenner) and Very Low Frequency (VLF) have been successfully used to ex- plore for groundwater [4, 5]. The Vertical Electrical Sounding is a Direct Current Resis- tivity tool of geophysics for the determination of the subsurface of a place. It has been used extensively for the determination of the aquifer potential for the drilling of boreholes. These days, it is used for in-depth geotechnical studies to determine the suit- ability of a site for the building of heavy structures including high rise buildings, bridges, and stadia. Hence, the method can be used in the study of subsurface layering of a place, to pro- vide information that will help in determining the nature of soil at the site and its stratification, selecting the type and depth of foundation problems (e.g. expansive soil and collapsible soil) and determining the depth and nature of bedrock, if and when encountered. In this work, D.C resistivity method-considered to be the quickest and the most economic technique, using VES- will be used to provide information on the zone. 2. Theory The theoretical study of the earth resistivity methods is to consider the case of completely homogenous isotropic medium. The equation which gives the potential due to a single point source of current at surface can be deduced from two basic equations; ohm’s law and divergence condition. The ohm’s law is given as: ∆ · J = 0 (1) The potential due to a single point source of current at earth surface is given as: U = ρI 2πr (2) where ρ denotes the earth density, I is the Current, 2π is a Con- stant, and r is Resistance. The current I, is passed through cur- rent electrode C1, C2 as shown in figure below (Figure 1): The potential difference at M is Um = Iρ 2π [ 1 r1 − 1 r2 ] . (3) The potential difference at N is Un = Iρ 2π [ 1 r3 − 1 r4 ] . (4) Figure 1. General four electrode configuration for resistivity measurement. The difference in potential at M and N is given as ∆U = Um−Un Un = Iρ 2π [ 1 r3 − 1 r4 − 1 r3 + 1 r4 ] . (5) Hence, the resistivity is given by ρa = ∆U I [ 2π { 1 r3 − 1 r4 − 1 r3 + 1 r4 }] . (6) Let K = 2π { 1 r3 − 1 r4 − 1 r3 + 1r4 } , therefore, ρa = KR (7) where, K is a Geometric factor and R = ∆U/I. 3. Methodology 3.1. Geology of the study area Marmara is a settlement of Chikun Local Government Area of Kaduna State, Nigeria. It is located geographically between latitude 10◦3′ and 10◦50′N and longitude 6◦4′ and 7◦5′E as shown in Figure 2. Chikun Local Government Area of Kaduna State geographically is in the North-West part of Nigeria which is part of West Africa Zone. These unconsolidated materials are known to reflect some dominant hydrologic properties, and the highest groundwater yields in Basement complex area found in areas of thick overburden overlying fractured zones and are characterized by relatively low resistivity [6]. Groundwater oc- currence in the area could be grouped into three; the weath- ered/fractured basement complex, newer basalts and river allu- vium. Climate and vegetation of study area has a typical sa- vannah climate with distinct wet and dry seasons, the wet sea- son ranges from April to October, while the dry season between end October to March, the average annual rainfall for Kaduna is 300mm. Rainfall generally reaches its peak in August and it has a characteristic mean annual temperature of 29◦C in March/April. The dominant River Kaduna controls the course of most of the rivers and streams in the area [7] 2 Kure et al. / J. Nig. Soc. Phys. Sci. 1 (2019) 1–5 3 Figure 2. Map of the Study Area showing all the VES Stations 3.2. Materials and method The basic requirement used for the study is Ohmega Resis- tivity Meter which displays resistance value digitally as com- puted from Ohm’s law. The Ohmega Resistivity Meter is pow- ered by a 12.5 Volts, D.C. power source. Other materials that accompany the equipment are, two measuring tapes, four ham- mers, thirty steel electrodes and four connecting wires (cables) for current and potential electrodes and a Global Positioning System (GPS). Generally four-electrode array are used at the surface, one pair for introducing current into the earth and the potential dif- ference these established are measured in the vicinity with the second electrode pair. Large numbers of electrode arrangement have been used for resistivity exploration[5]. Any of these can be used for lateral variation of resistivity or vertical variation of resistivity as a function of the current depth. In resistivity method, current are driven into the ground. Any variation of subsurface resistivity (ρ) alter the current flow which in turn af- fects the distribution of electric potentials. The potential, these established, are measured at the surface. Deviation of these pattern of potential difference from uniform earth represent the geological features of resistivity exploration. 3.3. Choice of electrode configuration The choice of a configuration for prospecting is dependent on a number of factors amongst which are the type of inves- tigation required, the terrain of the area of the investigation, and the position of the suspected geological structure. In ver- tical electrical sounding (VES) measurement, the centre of the electrode spread remains fixed but the separation of electrodes is progressively increased, hence choice of Schlumberger con- figuration (Figure 1), this configuration is therefore more cost effective since it saves time and manpower[1]. Applying equations 5 and 6 to the above array gives K = 2π/ { 1 a − b − 1 a + b − 1 a + b + 1 a − b } K = π/ { 1 (a − b) − 1 a + b } K = π/ { a2 2b − b 2 } , (8) therefore, ρa = π/ { a2 2b − b 2 } ∆/I (9) where, • a = separation between current electrode at the centre of the configuration, • b = separation between potential electrode at the centre of the configuration, and • ρa= apparent resistivity of the earth for a Schlumberger array. This last expression gives the approximate relationship between the apparent resistivity ρa, and the approximate resistivity value term usually given by terrameter during field measurement. The coefficient of ∆U/I is the geometric factor which is character- istics of the spread used. The electrical resistivity data was ac- quired with ABEM terrameter SAS 300 along with other equip- ment cables and electrodes. Table 1. Resistivity Range of some Rock Types (Compiled from [2, 4, 9, 10]) Rock type Resistivity (Ω − M) Topsoil/clay/silt 65 - 200 Laterite/Indurated Laterite 45 - 800 Weathered Basement 2 – 220 Fractured basement 218 - 520 Fresh basement < 1000 3.4. Data Interpretation The interpretation of the data was done using a computer program (RES1D version 1.00.07 Beta). To interpret the data, some initial value parameters were found from the observed field curve and fed into the program. The initial value parame- ters were taken such that each points of maxima, minima and in- flection indicate the existence of layer boundaries. The program uses these values to design a model curve, comparing it with the observed field resistivity curve. Where the modelled curve and the observed field curve are not in agreement, the initial value parameters were altered until the best agreement between the modelled curve and the observed field curve is obtained. At the point when the two curves were in best agreement and mini- mum error, the layer resistivity and thickness were recorded[8]. The geoelectric sections show the average variations of re- sistivity and thickness values of layers within the depth pene- trated in the study area at the indicated sections. Generally, the sections revealed two to four subsurface layers: topsoil/laterite, indurated laterite, weathered and fresh basement layers. 3 Kure et al. / J. Nig. Soc. Phys. Sci. 1 (2019) 1–5 4 Figure 3. Typical resistivity curve. 3.5. Data Processing Vertical Electrical Sounding (VES) data was collected on six (6) VES points using the Schlumberger configuration. The first stage in any interpretation of apparent resistivity sounding curves is to note the curve shape [10, 11]. Before applying a more complicated method of interpretation, it is useful to con- sider a rough idea. The data collected in the study area were reduced and computed. The computed apparent resistivity val- ues were then fed into a computer Geosoftware (ResID ver- sion 1.00.07 Beta) for resistivity curve interpretation. The in- terpreted data was used to generate the final model geoelectric parameters presented in Figure 3 which shows a typical exam- ple of resistivity curve. 3.6. Resistivity Values Table 2 shows the resistivity range of geoelectrical layer compiled from some previous works around the study area. Table 2. Resistivity values of some rock types basement area (compiled from [1, 2, 3, 9, 10]) Soil and Rock Types Resistivity Value (Ωm) Fadama Loam 30 - 90 Sandy Clay and Sandy Silty 100 - 200 Sand and Gravel Laterite 10 - 150 Weathered Basement 20 - 200 Fractured Basement 500 - 1000 Fresh Basement > 1000 4. Results and Discussion Results obtained indicates that the terrain is highly variable in thickness and resistivity. According to Alao[11] and Fadele et al. [12], the engineering competence of the subsurface can be qualitatively evaluated from the layer resistivity, since the higher the layer resistivity value, the higher the competence of the Layer. The final model geoelectric parameters along the six VES stations were used for the preparation of the geoelectric/geologic Figure 4. Shows the Geoelectric and Geologic section of the Study Area. sections for one of the profiles (Figure 4). VES A1 is recom- mended for the construction of high-rise buildings, this is be- cause the depth to basement is approximately 25m. Hence, the subsurface Layers at A1 are not favorable for groundwater ex- ploitation. VES A3, A4 and A6 are not recommended for high- rise buildings. However, surface Layers appear to be highly favorable for groundwater exploration which has the highest aquifer thickness of 46m at a depth of 48m. This is because the thickness of the aquifer at this point ranges from 18m to 42m with a depth of approximately 6m to 10m respectively. According to Aboh [1], the areas underlain with low bed rocks resistivity (< 1000Ωm) in Figure 4 could indicate frac- tures or faults. Hence, are less competent. Based on the resistiv- ity values of the different geoelectric layers and the various ge- ologic units obtained, the study has identified VES stations A1, A2 and A5 as competent zones. These regions were chosen due the fact that, there was no indication of geologic trend such as fracture or fault that could cause structure failure. Fadele, et al. [12], Dogara, et al. [13] and Alao [11], noted that the seasonal variation in the saturation of clay which causes ground move- ment could cause by clay swells and shrinkages; these factors could contribute to structural defect of building but the presence of consolidated laterite in the study area beneath the clayey top- soil which extends beyond 6.5m, could reduce the danger posed by clay formation to large buildings. However, the topsoil resis- tivity of the chosen zones is characterized with low resistivity and thin in thickness, which may not be good enough to sup- port massive buildings. Hence, this superficial cover is loose humus/loam and clay formation at the topsoil in some part of the study area; thus need to be totally excavated in any future structural and engineering works. In summary, the interpreted data shows that there are no indications of any major trends such as fracture or faults in the area chosen as competent zones (A1, A2 and A5), that could aid building subsidence. 4 Kure et al. / J. Nig. Soc. Phys. Sci. 1 (2019) 1–5 5 5. Conclusion The investigation shows the geo-electric and geologic sec- tions derived from the study area which has revealed five sub- surface layers, namely: lateritic topsoil followed by clayey/ sandy/ silty, layer, weathered basement and fractured/ fresh base- ment rocks. A broad section of the area was classified as com- petent for civil works. There are no indications of any major linear structure such as fracture, faults or voids that could aid building subsidence in the study area. However, the seasonal variation in the saturation of clay which causes ground movement could cause havoc in building construction due to the sensitiveness of soils to moisture loss or gain (swells and shrinkages) in VES A3 and A5. Consider- ing the presence of consolidated laterite which extends beyond 3.0 m could reduce the danger posed by clay formation to large buildings as well as tarred roads. In summary, the interpreted data shows that there are no indications of any major trend such as fracture, cavities, voids or faults that could aid building sub- sidence in the area identified as competent zones. Other prob- able factors such tree roots, organic deposits, superficial cover that could cause foundation defects can therefore be eliminated by total evacuation of the topsoil. The lithologies Thickness and Depth are as follows: Top- soil ranges from 1.0m − 1.8m; Indurated Laterite ranges from 1.9m−6.8m; Weathered/Fractured Basement has infinite thick- ness. The Areas found to meet the conditions of high Basement resistivity and shallow Depth to Basement, are found to be at VES A1, A2, A3, A4 and A5. Hence, A1 have been found to be the competent zones within the Study Area, and are good for the construction of highrises, roads and bridges. VES A3, A4 and A6 those not meet the conditions for construction of high- rises building. However, VES A3, A4 and A6 Areas appears to be very viable for groundwater exploration with aquifer thick- ness ranges from 18m−42m with depth ranges 6m−10m. And the average Depth to Fresh Basement is 43m. References [1] H. O. Aboh, “Geotechnical characterization of surface materials in Kaduna Area, Kaduna State”, Zuma Journal of Pure and Applied Science 4 (2002) 23. [2] C. O. Ajayi & M. Hassan, “The Delineation of Aquifer Overlying the Basement Complex in the Western Part of Kubani Basin of Zaria Nige- ria”, Journal of Mining and Geology 26 (1990) 117. [3] N. Kure, H. O. Aboh, R. Jimoh, O. A. Joseph & I. H. Daniel, “The Delin- eation of Potential Groundwater Aquifers within Basement Complex in ABU Zaria, Nigeria”, British Journal of Applied Science and Technology 19 (2017) 1. [4] E. M. Shemang, Electrical Depth Sounding of Selected Well Sites within the Kubani River Basin, Zaria, Nigeria, Unpublished MSc Thesis, Physics Department, ABU Zaria, 1990. [5] O. Koefoed, “An appropriate method of resistivity sounding interpreta- tion”, Geophysical Prospecting 24 (1979) 617. [6] M. O. Olorumfemi, S. A. Fasuyi, “Aquifer types and geoelec- tric/hydrogeologic characteristics of part of central basement terrain of Nigerian (Niger State)”, Journal of Earth Science 16 (1993) 309. [7] P. O. Udo, Physical geography of Nigerian, Heineman Education Pub- lishers, 1982. [8] K. M. Lawal, The application of reduction to the pole technique for the interpretation of aeromagnetic data of Zaria area, North Central Nigeria, 1998. [9] N. K. Abdullahi, E. E. Udensi, A. Iheakanwa & B. E. Eletta, “Geo- electrical method applied to evaluation of groundwater potential and aquifer protective capacity of overburden units”, British Journal of Ap- plied Science and Technology 4 (2014) 2024. [10] M. W. Telford, L. P. Geldart, R. E. Sheriff & D. A. Keys, Applied geo- physics, Cambridge University Press, London, 1970. [11] J. O. Alao, Delineation of the Interfacial Configuration in a Section of Millennium City, Kaduna, Unpublished MSc. thesis, Physics Dept. Kaduna State University, Kaduna, 2017. [12] S. I. Fadele, S. B. Jatau, & N. O. Patrick, “Geophysical engineering inves- tigation around Makiyaye Village, shika area within the basement com- plex of North-Western Nigeria”, International Journal of Engineering Re- search and Applications (2012) 1143. [13] M. D. Dogara, J. Alao, H. Abdullahi, J. Ezekiel, J. George, & R. A. Ahammed, “Delineation of the geotechnical parameters within the Kaduna refining and petrochemical corporation layout”, World Journal of Applied Physics 2 (2017) 36. 5