Investigation of groundwater potential using integrated geophysical methods in Moloko-Asipa, Ogun State, Nigeria

Water is essential for livelihood, development, and industrial growth. Its exploration in sufficient quantity is required where it does not freely occur on the surface. This research was aimed to delineate aquifer regions and provide information on the subsurface lithology of Moloko-Asipa Southwestern Nigeria. A combination of eight traverses investigated with very low frequency electromagnetic (VLF-EM) method at 5 m constant sampling interval and ten vertical electrical sounding (VES) were carried out in the survey. Measurements from the VLF-EM survey were processed with Karous and Hjelt filtering to give the resistivity contrast across the selected profiles. The VES data processing involved an automatic approximation of the initial resistivity and thickness of the geoelectric layers with IPI2Win and further filtering by WinResist iteration. Estimation of Dar-Zarrouk parameters was also employed to investigate the aquifer protective capacity of the area. The processed VLF-EM results showed the geology of the area to an average depth of 25 m. The geoelectric section of the VES data revealed minimum of 3 layers from sandy top soil to weathered layer and fresh basement with an average resistivity values of 1,816, 926 and 17,503 Ωm, respectively. The integration of VLF-EM and VES in the investigation revealed that the potential for groundwater exploration in the study area is poor due to the thin nature of the weathered layer and its shallow depth to basement. The aquifer protective capacity of the area was likewise inferred to be poor.


Introduction
The study of groundwater exploration is increasing in recent times due to its inevitable needs for livelihood, development and industrial growth. Groundwater plays a vital role as it is a dependable primary source of freshwater for all needs in areas where surface water is inadequate or not available (WWDR 2015). Though groundwater exploration in some areas like hard rock terrain is very challenging and sometimes tricky as its potential depends mainly on the thickness of the weathered or fractured layer overlying the basement, a thorough hydrogeological investigation dwindles these challenges (Al-Garni 2009). New developments have been carried out to characterize the nature and the pattern of flow of subsurface water which generally categorized groundwater into groundwater found in the unsaturated zone and groundwater located in the saturated site (Delleur 2006). This categorization is due to the differences in the physics of water flow in saturated and unsaturated zones. In the study of groundwater potential, the knowledge of the geology of the area under investigation is very crucial as the understanding of the geology of the site helps to delineate possible saturated zones from unsaturated zones beneath the crust (Nimmo 2005;Hassan et al. 2017). Apart from having a good knowledge of the geology of the study area, correct interpretations of hydraulic test's result and perfect interpretations of resulting predictions are needed to have a successful groundwater exploration (Badmus and Olatinsu 2012;Akintayo et al. 2014;Lukuman et al. 2015). The application of different principles of physics and earth science is generally related to providing a good understanding of the subsurface structure, which are very important for both scientific and economic reasons (Kamil 2008). The field source that is required to generate information about an object of interest in the subsurface can either be artificial or natural (Henderson 1992). The observed natural fields can originate from a combination of sources in the environment of the target that is being investigated. Hence a filtering process is needed to eliminate the field signals from all other sources to get the specific signal that will indeed characterize the target (Janvier et al. 2015). One advantage of natural field source is that it does not require any additional deployment to generate the field, unlike the artificial source where a further deployment is needed. However, the artificial fields have an advantage such that the source parameters can be predetermined and can be used to eliminate noise from the obtained signal (James 2009;Ram and Rao 2018;Odipe et al. 2018;Nuraddeen et al. 2019). This study was aimed to delineate the regions of aquifer and provide information on the geology of the study area. With the view of characterizing the subsurface lithology, identifying structures that would reveal information on the potential for groundwater exploration in the area is also a goal. Olayanju (2011) combined very low frequency electromagnetic (VLF-EM) with offset Wenner resistivity methods of geophysical survey to evaluate the groundwater spring source in Akure, Southwestern Nigeria. It was concluded that the adoption of an integrated method of hydrogeological survey provides a better result for a successful groundwater exploration process. A geophysical investigation for groundwater development involving the electrical resistivity and electromagnetic VLF methods was likewise carried out by Afolayan et al. (2004) in the premises of Obafemi Awolowo University, Ile-Ife, southwestern Nigeria. Three VLF-EM profiling were carried out and used as a reconnaissance survey to locate the conductive zones where six vertical electrical soundings were done. The VLF normal and filtered real component anomalies from the survey revealed five major geological interfaces while the geoelectric section prepared from VES interpretation delineates four subsurface layers. The researchers used a lithological log from a test borehole at one of the VES stations to corroborate the predicted sequence and its structural disposition of the subsurface by the VLF-EM and VES investigation. Falade et al. (2019) integrated very low frequency electromagnetic (VLF-EM) with vertical electrical sounding (VES) electrical resistivity methods to determine groundwater prospect in Ilara-Mokin, Akure, Nigeria. The VLF-EM profiling was carried out on eight traverses and the results revealed the points which are relatively conductive with depth. Eighteen VES stations with maximum electrode separation of 65 m along the six traverses trending East-West were then positioned based on the results of the VLF-EM survey. The researchers concluded that the integration of VLF-EM and VES aids groundwater development and recommend that borehole logging can be adapted to characterize the aquifer for better integration results. The integration of VLF-EM and VES geophysical methods for groundwater exploration in hard rock areas was also studied by Adelusi et al. (2014). The problems associated with VES in hard rock terrains with surface or near-surface expression of outcrops were identified. The researchers submitted that the use of VLF-EM becomes imperative in mapping the basement structures relevant to groundwater development in such areas. The VLF-EM method yields a higher depth of penetration in hard rock areas because of their high resistivity (McNeill and Labson 1991). One of the advantages of integrating VLF-EM with VES is that its collection of data is relatively fast and inexpensive compared to many other geophysical methods (Olafisoye et al. 2012). VLF-EM on its own is very effective in locating zones of high electrical conductivity, though all electrical conductivity affects it, and it needs to be filtered during interpretation (Yusuf 2016). VLF-EM method has probably been the most popular tool for quick mapping of near-surface geologic structures in mineral exploration. It is, however, increasingly used for shallow groundwater exploration as a reconnaissance tool for weathered layer investigations (Palacky et al. 1981). Also, VLF-EM has proved to be successful in identifying deep water-bearing fractures in bedrock (Sundararajan et al. 2007).
Vertical electrical soundings (VES) and geophysical logs were employed by Oladele and Odubote (2017) to map and characterize aquifer units in a typical complex transition zone of Ijebu Ode. The study area is a zone between the Precambrian Basement rocks and Cretaceous sediments of Dahomey Basin, Southwestern Nigeria (Osinowo and Olayinka 2012). Sixteen VES were acquired using Schlumberger array method with maximum current electrode separation of 900 m. The interpretation of the VES data revealed four geoelectric layers with the weathered layer at 80 m depth in the subsurface. The researchers concluded that the study area has a low potential for groundwater exploration. Badmus and Olatinsu (2010) likewise carried out aquifer characterization and groundwater recharge pattern in a typical basement complex Formation at Federal College of Education, Osiele Abeokuta southwestern Nigeria. The survey involved thirty-four (34) vertical electrical sounding (VES) using Schlumberger electrode array at 200 m electrode separation. Seven geologic formations were revealed in their studies, viz topsoil, sandy clay, clayey sand, shale/clay, sandstone, fractured basement and fresh basement. It was concluded that the weathered and fractured basement in the study area constitutes the main aquiferous unit. Mosuro et al. (2011) carried out geophysical investigation for groundwater exploration at Kobape, via Abeokuta southwestern Nigeria using Schlumberger array for the vertical electrical sounding (VES) data acquisition. It was submitted that Kobape area lies within the sedimentary terrain of Dahomey Basin, which consist Abeokuta group and according to Agagu (1985). It was concluded based on the results of the VES data that there exist four to five geoelectric layers that their lithology is essentially sandy and highly resistive. The researchers reported that the promising aquiferous layer in the study area was observed at a depth of 270 m below the surface, and recommended that a borehole drilling should be carried out in the peak of the dry season in the area.

Physiographic and geological setting of the study area
The study area falls along Abeokuta-Sagamu Expressway in which Sagamu is a sedimentary terrain and Abeokuta town is Basement complex. The site lies between latitude 7 0 3 ′ 18 ′′ N to 7 0 3 ′ 25 ′′ N, and longitude 3  with Lagos state to the south (Fig. 1). The rejuvenated crystalline rocks found in the study area are dominated by granites, schist, gneiss and quartzite and are of Precambrian age to early Paleozoic age. The rocks extend from the north-eastern part of the Ogun state running southwest and dipping towards the coast (Ako 1979). Selby (1985) reported that rock often breaks down quickly to produce a zone of weathered materials of saprolite or laterite. The surface soils are usually underlain by red-brown silty clay which does not serve as a good aquifer, and according to , basement aquifers occur within the weathered residue overburden (the regolith) and the fractured bedrock. Development of the regolith components is by shallow boreholes and wells which are liable to be drilled by lightweight percussion rigs as viable aquifers occur within fractured bedrock. Development of bedrock components requires interaction with available storage in overlying or adjacent saturated regolith, or other suitable formations such as alluvium for effectiveness. Typical metamorphic rocks found in basement complex are characterized by various folds, structures of different degrees of complexity, faults and foliation. The rock types that characterize the study area range from granite, schist, gneiss and pegmatite. The individual rocks have various hydrogeologic characteristics and belong to the stable plate, which was not subjected to intense tectonics in the past (Usman and Ibrahim 2017). Hence, the underground faulting system is minimal, and this has contributed to the problem of underground water occurrence in the area. The geology of Kobape shows that it lies within the sedimentary terrain of Dahomey Basin, which consists of Abeokuta group (Agagu 1985). This formation reveals sandstones, grits, siltstones, limestones, shales and sands lithologies (Ariyo and Adeyemi 2012). The extracted geology of the study area is shown in Fig. 1.

Materials and methods
Very low frequency electromagnetic (VLF-EM) profiling and vertical electrical sounding methods were employed in the survey. The VLF-EM survey was carried out using the ABEM WADI VLF instrument, while electrical resistivity meter was utilized for the resistivity data collection. The VLF-EM enables surveying for electrical conductors beneath the surface without any contact with the ground, and eight different traverses were carried out in the study area (Fig. 2). Electromagnetic field signals that are being generated worldwide from a high powered transmitter are essentially uniform in the atmosphere and can penetrate the ground to a depth of about 30 m. The VLF signal transmitters generate a plane electromagnetic (EM) waves which consist of both the vertical electric field and horizontal magnetic field components; each perpendicular to the direction of propagation. The vertical electric fields induce electric currents, particularly on electrically conductive zones in the subsurface to produce secondary magnetic fields that are detectable at the surface through the deviation of the normal radiated field. The magnetic component of the VLF signal is mainly used for the field measurement. According to the basic electromagnetic (EM) theory, the primary EM field is shifted in phase when encountering a conductive body which then becomes the source of a secondary field. The VLF instrument detects the primary and secondary fields. It then separates the secondary field into real (in-phase) and imaginary (quadrature) components based on the phase lag of the secondary field. These two components of the secondary fields are sometimes referred to as the tilt ( ) and ellipticity (e).
A VLF receiver is tuned to a specific frequency within the range of 15-30 kHz of the high powered transmitter and directs the signal towards the Earth's surface. The VLF receiver then measures the field tilts while the real and the imaginary components of the field data are being recorded to processing (Ndatuwong et al. 2013). In VLF-EM survey, the in-phase and quadrature are the most important field measurements and can be expressed as the normalized real and quadrature components of the vertical magnetic field. Equations 1 and 2 can be used to derive the apparent conductivity of the ground (Reynolds 1997).
where Hx, Hy and Hz are the magnetic field components in x, y and z directions.
Linear inversion of the acquired real component of the field data involved the application of the Fraser and Karous-Hjelt filtering and processing. The filtering is done to remove geologic noise and the noise originating from the transmitter (Benson et al. 1997). A qualitative interpretation of VLF-EM data is based on the crossover points between the real and imaginary lines, which appears as positive peaks in the Fraser-filtered real curve. These regions constitute anomalous zones which can be attributed to the presence of vertical conductor or lateral contacts of different resistivity values beneath the surface (Srigutomo et al. 2005).
The process of Karous-Hjelt filtering yields a pseudosection of relative current density variation with depth. A higher value of relative current density corresponds to resistive subsurface structures and a lower value tending towards y negative corresponds to a conductive subsurface structure. The application of Karous-Hjelt filtering process to the real component of the field signal reveals the distribution of the current density, which is displayed and interpreted as a map of current density (Benson et al. 1997). The filtering process, therefore, provides a visual indication of the depth of the various current concentrations and hence the spatial dispositions of subsurface geological features, such as mineral veins, faults, shear zones and stratigraphic conductors (Khalil et al. 2009). Fraser Filter is another smoothening technique similar to the Karous-Hjelt filter which transforms the VLF anomalies to contours in such a way that proper crossovers are transformed to positive peak amplitude. Reverse crossovers, on the other hand, become negative values, thereby enhancing the signals of the conductive structures (Telford et al. 1990). The centre of the anomalous structure may fall directly under the peak of the Fraser filtered data. The advantages of the Fraser filter include (i) complete removal of DC bias and great attenuation of long wavelength signals; (ii) complete removal of Nyquist frequency related noise; (iii) phase shifts in all frequencies by 90° and (iv) having the band-pass centered at a wavelength of five times the station spacing (Sundararajan et al. 2006;Khalil et al. 2009). Two-dimension image of Karous-Hjelt filter shows a typical inverted current density sections using a software; KHFfilt ver. 1.1a (Pirttijärvi 2004). The VLF-EM data are useful to obtain a qualitative view of the sub-structure, particularly after filtering the data and analysing the current density across the section. The results obtained from the VLF-EM survey with KHFfilt package give the Fraser outputs and 2-D subsurface conductivity images. The Fraser filter uses four consecutive data points in a situation where the data have been acquired at a regular interval. The sum of the first and second data point is subtracted from the sum of the third and fourth values and plotted at the midpoint between the second and third tilt-angle stations. The point of inflection or crossover point on the Fraser filtering curve indicates the presence of a target which the depth of occurrence can be estimated as half of the distance between the peak negative and peak positive amplitudes (Sriramulu et al. 2017).
A vertical electrical sounding (VES) using Schlumberger array method was employed to determine the depth of conductive zones and estimate the layers associated with the point under investigation. The array is capable of isolating successive hydrogeologic layers beneath the surface using their resistivity contrast (Mosuro et al. 2011). The observed field data, which is the ratio of the resulting voltage to the injected current, was obtained from the resistivity meter, and it is only a measure of the resistivity of the subsurface material. This measured resistance (R) was multiplied by the corresponding Geometric Factor (a function of the electrode spacing and configuration) to obtain the apparent resistivity measured in Ohmmeter (Ωm) needed for the interpretation. The sounding curves for each point were obtained by plotting the computed apparent resistivity against the half current electrode spacing (AB/2), and IPI2Win software was first used to perform an automatic approximation of the initial resistivity and thickness of the various geoelectric layers. The generated results were further processed with WinResist iteration to filter possible error from the equipment and smooth the sounding curves.

Dar-Zarrouk parameter and its application
The Dar-Zarrouk parameters (DZP) are the most essential parameters in electrical sounding prospecting, which can be used directly in aquifer protection studies and evaluation of hydrologic properties of the aquifer (Adeniji et al. 2013). The parameters are longitudinal conductance (S) and transverse resistance (T). The longitudinal unit conductance in Ω −1 is the ratio of the layer thickness (h), to resistivity (ρ) while the transverse unit resistance in Ωm 2 is the product of the layer thickness (h) and resistivity (ρ) shown in Eqs. 3 and 4 (Ahamefula et al. 2012). Mathematically: A geoelectric layer is always characterized by these two basic parameters; layer resistivity i and layer thickness h i for every ith layer (where i = 1 for the surface layer and n for number of layers) The sum of the entire longitudinal conductance for a layered ground is referred to as Dar-Zarrouk function and the sum of the corresponding transverse resistance is called Dar-Zarrouk variable as shown in Eqs. 5 and 6 (Isaac 2014).

Results and discussions
The processing of the VLF data generates a graph for real and imaginary component for each traverse and corresponding pseudo-sections which are the plots of depth versus station interval in 2-D inversion. Colour scales are usually shown beneath the 2-D inversion image which represents conductivity band. The conductivity band, generated from the processing of the obtained VLF measurements, ranges from −60 to 10% with a colour spectrum ranging from deep blue, turquoise blue, yellowish-green, yellow and orange to red. The red colour indicates the high current density (conductive body), while the dark blue colour indicates low current density (resistivity body) and the intermediate green colour reveals a moderate resistive bodies such as shear zone (Srigutomo et al. 2005;Sundararajan et al. 2006).

VLF interpretation for traverse 1
The Fraser filtering curve of the VLF responses obtained from traverse 1 shows relative distribution of conductive and resistive bodies of various amplitudes across the east-west direction. The maximum conductivity response of 40% was observed at position 108 m while a minimum conductivity response of -55% occurs at lateral position of 116 m. Other relatively high and low responses along the traverse indicate weakly conductive and resistive bodies at various depths. The 2-D inversion by Karous-Hjelt filtering revealed moderately resistive bodies from the base of the profile at position 82 m to 100 m close to the surface and same bodies are longitudinally distributed from 100 to 132 m along the traverse. Pocket of conductive material was observed close to the surface between 8 and 22 m, 36 to 48 m and also from the base of the profile at position 62 m towards 88 m close to the surface (Fig. 3).

VLF interpretation for traverse 2
Traverse 2 Fraser filtering curve revealed maximum positive amplitude of 20% and minimum negative amplitude of −15% at positions 68 and 92 m along the traverse, respectively. Another high negative amplitude of −14% occurs at position 108 m which revealed the presence of another resistive body along the east-west direction of the profile. The 2-D inversion by Karous-Hjelt filtering revealed that no anomaly is present from the beginning of the traverse to around 38 m where a conductive zone is observed extending to 58 m along the traverse. Another conductive body is observed between 70 m from the base of the profile to 98 m towards the surface along the traverse (Fig. 4).

VLF interpretation for traverse 3
The Frazer filtering curve of traverse 3 has conductivity range of −25 to 16%. The peak positive and negative responses were observed at 62 m and 142 m in the east-west direction of the profile. Another high conductivity response of + 15% was also observed at position 72-78 m showing the presence of a conductive anomaly in the region. The 2-D inversion by Karous -Hjelt filtering revealed the presence of a resistive anomaly between 44-80, 92-98 and 132-148 m along the traverse. Scattered conductive bodies were also observed to vary laterally between positions 38-132 m along the traverse (Fig. 5).

VLF interpretation for traverse 4
Traverse 4 revealed a peak negative anomalous response of −28% at position 70 m and the highest conductivity response is 23% at position 135 m. Other relatively dwarf peaks in the order of 18% correspond to the weakly conductive material distribution across the traverse. The 2D-Inversion by Karous-Hjelt filtering revealed the material distribution along traverse to be a complex varying mixture of resistive to conductive bodies at closely varying intervals. The investigation carried out on this traverse covers up to 200 m thereby revealing the subsurface structure up to a depth of 30 m (Fig. 6).

VLF interpretation for traverse 5
The Fraser filtering curve from the measurement taken from traverse 3 has a maximum positive amplitude response of 23% at position 112 m and minimum negative amplitude response of -15% at position 122 and 132 m along the east-west direction of the traverse. The 2-D inversion by Karous-Hjelt filtering revealed moderately resistive anomaly of about 12 m thickness between 42-62, 90-100 and 130-148 m along the traverse. The high conductivity responses from the Fraser filtering curve correspond to the bodies observed from the base of the profile at 52 m to 78 m on the surface. Other isolated conductive bodies were observed between 102 and 132 m in the North-South direction of the profile (Fig. 7).

VLF interpretation for traverse 6
The Frazer filtering curve of the measurements from traverse 6 revealed has conductivity range of −16 to 18%. The peak positive and negative amplitudes were observed at   (Fig. 9).

VLF interpretation for traverse 8
The Frazer filtering curve of the measurements taken from traverse 7 has peak negative and positive amplitudes -22% and 22% at 72 m and 88 m in the East-West direction of the profile. Relatively weak amplitude of an average of ± 15% on the traverse indicates weakly conductive and resistive bodies of various depths along the traverse. The 2-D inversion by Karous-Hjelt filtering revealed moderately resistive bodies at 50-90 m and 96-132 m in the North-South direction of the profile (Fig. 10).

VES results
The nature of the VES curves obtained from the analysis of each sounding results indicates the interplay between low and high resistivity layers in the study area. The output is the inverse resistivity model providing layer-wise distribution of resistivity value and thickness of the corresponding layer with curve types are shown in Table 1. Most of the resulting models produced a low RMS relative error of the order of 2%. The resistivity values extracted from Loke (2000) for common rocks and soil materials in the survey area were used to categorize the lithology of the VES data into geoelectric layers (Table 2).
Four geoelectric layers ranging from sandy topsoil, sandy horizon, and weathered layer to fresh basement were inferred at VES 1 with resistivity value of 564, 3592, 726 and 17,204 Ωm, respectively. The respective depths of the first three layers are 0.5, 1.3 and 10.3 m while that of the last layer cannot be determined. The curve type generated from WinResist iteration is KH type ( 1 < 2 > 3 < 4 ) as shown in Fig. 11. H type curves in a typical basement complex environment that contains a low resistivity intermediate layer underlain and overlain by more resistant materials are commonly water saturated. Such regions are often characterized wit high porosity, low specific yield and low permeability. The main There are four geoelectric layers which are sandy topsoil, sandy horizon, weathered layer and fresh basement at VES 2 with corresponding resistivity values of 3596, 1406, 725 and 27,995 Ωm. The relative depths of the layers are 0.5, 5.8, 5.7 m and that of the last layer cannot be determined. QH-type curve where 1 > 2 > 3 < 4 is revealed by the iteration carried out for this second sounding point (Fig. 12).
The geoelectric layers of VES 3 show that it consists of three layers, namely sandy topsoil, weathered layer and fresh basement with corresponding thickness of 1.6 and 19.8 m while the thickness of the third layer cannot be determined. The resistivity values of the layers are 1448.3, 859.5 and 56,382.1 Ωm from the topsoil to fresh basement. The curve type for this station is H ( 1 > 2 < 3 ) given that this curve type has the main aquifers found at the base of the weathered profile (Fig. 13).  Three geoelectric layers inferred as sandy topsoil, weathered layer and fresh basement with resistivity values of 1179, 1753 and 10,285 Ωm are shown at VES 4 with respective depths of 0.71 and 26.1 m for the first two layers while the depth of the third layer cannot be determined (Fig. 14). The WinResist iteration for this sounding point reveals A-type curve ( 1 < 2 < 3 ).
The layer model interpretation for VES 5 reveals five layers ranging from sandy topsoil, sandy horizon, weathered layer, sandstone to fresh basement with corresponding resistivity values of 1241, 4333, 635, 4829 and 12,679 Ωm. The relative thicknesses of the layers are 0.5, 0.84, 2.23, 64.9 m where that of the last layer cannot be determined. The curve type from the iteration for this fifth sounding point is KHKA-type showing that 1 < 2 > 3 < 4 > 5 (Fig. 15).
The sixth sounding also reveals four geoelectric layers: sandy topsoil, sandstone, weathered layer and fresh basement, with resistivity values of 825, 1378, 972 and 6608 Ωm and respective thicknesses of 0.5, 0.8 and 8.4 m as the thickness of the fresh basement cannot be determined. KH-type curve ( 1 < 2 > 3 < 4 ) is thus revealed from the iteration carried out on the sounding data (Fig. 16).
Three geoelectric layers ranging from sandy Topsoil, Weathered Layer to Fresh Basement with A-type curve ( 1 < 2 < 3 ) were revealed from the iteration of VES 7 data. The resistivity values of the layers are 341, 882 and 30,043 Ωm with corresponding thickness of 0.6, 13.2 m and undetermined thickness of the last layer (Fig. 17).
Three geoelectric layers were also revealed on VES 8 section with resistivity values of 1735, 907 and 4576 Ωm and respective thickness of 0.5, 9.5 m and that of the last layer cannot be determined. The layers are sandy Topsoil, Weathered Layer and Fresh Basement with H-type curve ( 1 > 2 < 3 ) as revealed from the resistivity iteration in Fig. 18.
The ninth sounding reveals four geoelectric layers: sandy topsoil, sandstone, weathered layer and fresh basement, with resistivity values of 6139, 2021, 1000 and 5539 Ωm with corresponding thickness of 0.5, 1.7 and 9.8 m as the thickness of the fresh basement cannot be determined. QH-type of curve is thus revealed from the iteration carried out with WinResist ( Fig. 19) revealing that 1 > 2 > 3 < 4 .
Four geoelectric layers were also recorded at VES 10 ranging from sandy topsoil, sandy horizon and weathered layer to fresh basement with corresponding resistivity values of 1092, 472, 804 and 3721 Ωm. The relative thicknesses of the layers are 0.5, 2.4 and 14.5 m while that of the last layer cannot be determined. The curve type of the iteration for this sounding point is HA-type revealing that the 1 > 2 < 3 < 4 (Fig. 20).
The characteristic longitudinal unit conductance prepared with Eq. 5 provides an estimate of the Dar-Zarrouk parameters for all the VES locations (Table 4) and was used to evaluate the protective capacity rating of the study area. The aquifer protective capacity characterization is based on the values of the longitudinal layer conductance of the overburden rock units at each sounding stations. The standard values for Longitudinal Conductance and aquifer protective capacity rating are shown in Table 3. The application of these parameters is aimed at unraveling the potential for viable groundwater exploration in the study area. The ability of the earth medium to acts as a natural filter to percolating fluid and retard or filter percolating ground surface polluting fluid is a measure of its protective capacity (Olorunfemi and Dan-Hassan 1999). The highly impervious sandy or clayey overburden, which is characterized by relatively high longitudinal conductance, offers a better protection to the underlying aquifer; hence, the protective capacity inferred from the study area is poor as the values of the longitudinal unit conductance (S) range between 0.010 and 0.024 mhos ( Table 4).
The potential for groundwater exploration in the study area is therefore evaluated based on the weathered layer thickness and the aquifer protective capacity. The inferred weathered layers constitute the water saturated zones which are aquifer units. The interpretation carried out is useful to provide information on the curve type and determine the geoelectric properties of the subsurface layers. The geoelectric section was drawn for the closest sounding in the study area to provide information on the subsurface geology (Fig. 21).

Conclusion
The geophysical investigation for groundwater exploration in Moloko-Asipa was carried out using VLF-EM and VES methods of survey. The interpreted results from the VLF-EM survey revealed the nature of the subsurface in the area to an average depth of 25 m and were inferred to be dominated with highly resistive material. These results were verifiable from the VES interpretation that revealed the subsurface structure of the area investigated to be composed of 3 to 4 geoelectric layers: sandy topsoil, sandstone, weathered layer and fresh basement. VES1, VES 2, VES 6, VES 9 and VES 10 revealed four layers, while VES 3, VES 4, VES 7 and VES 8 revealed three layers. VES 5, on the other hand, revealed 5 layers with the weathered layer found in between two sandstone layers. The resistivity range of the topsoil was 341-3596 Ωm with thickness ranging from 0.5 to 1.6 m. The resistivity of the sandstone is considerably high, ranging from 1378 to 4333 Ωm with thickness ranging from 0.8 to 25.4 m. The weathered layer had a resistivity range of 635-1000 Ωm with thickness ranging from 2.2 to 19.8 m. The fresh basement has resistivity range of 3721-56,382 Ωm which extends to infinity. The highly impervious sandy overburden is characterized by relatively low longitudinal conductance, and it offers a weak protection to the underlying aquifer. The protective capacity inferred from the study area is poor as the values of the longitudinal unit conductance range from 0.010 to 0.024 mhos. The geophysical investigation for groundwater potential in the study area with the integration of VLF-EM and VES methods revealed a thin and shallow weathered layer that has an average depth to bedrock of 12 m. Generally, the investigation inferred a low groundwater yield with poor protective capacity of sandy overburden. Manual well digging that is targeted on a suspected aquifer at a peak of dry season is advised to economically manage the quest for groundwater exploration in the area. The application of other verifiable geophysical techniques like well logging is, however, recommended for further studies.

Conflict of interest
The authors declare that they have no conflict of interest in this research work. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.