In ensuring livelihood and sustainability across the world, groundwater resources play a major and concrete role. Consumption of groundwater as a feasible source of drinking, domestic, industrial and agricultural needs has proven to be not only safer but also more economical than surface water, as it is commonly unpolluted and available. In recent years, investigation of groundwater sources has become a burning issue and a major concern as groundwater basins are being rapidly stressed due to population explosion, high level of urbanisation, industrialisation and other human activities. Presently, the percentage increase in water usage on a global scale has exceeded twice that of the population [1].
Pores and fractured rock formations in the sub-surface are usually hosts of groundwater. In the basement terrain, groundwater occurs within the overlying unconsolidated material derived directly from weathering of rocks and fractured/faulted bedrock, while in the sedimentary terrain, groundwater occurs within the porous and permeable layer of the saturated zone in the sub-surface [2, 3]. Over the years, groundwater exploration has been carried out using geophysical methods, which include electrical resistivity surveying, electromagnetic techniques and seismic methods, to obtain accurate information about the sub-surface settings, such as aquifer’s nature, type and depth of materials (consolidated or unconsolidated), depth of weathered or fractured zone, depth to groundwater, depth to bedrock and salt intrusions into groundwater [4].
Aquifers in the Precambarian basement complex are vulnerable to surface or near-surface contaminants as they commonly occur at shallow depths. Hence, successful exploration of groundwater in a basement terrain requires proper understanding of the hydrogeological characteristics of the aquifer units in relation to their susceptibility to environmental pollution and assessment of their protective capacity [3, 5]. One of the most effective ways of evaluating an environment without interfering with the hydrogeological system is through geophysical studies [6]. Over the years, geophysical survey using the vertical electrical sounding (VES) method has been applied in groundwater exploration within the basement complex rocks in Nigeria [7, 8, 9, 10].
VES using the Schlumberger array method was carried out at 30 different stations in the study area with the aim to determine the geoelectric parameters, such as resistivities and thicknesses of the sub-surface layers and their hydrogeological properties. This study was also aimed at evaluating the groundwater potential of the area, establishing the aquifer protective capacity of the overlying formations, especially its isolation from contamination, and recommending suitable points for groundwater positioning.
The study area, as shown in Figure 1, is located between latitudes 7°10’N and 7°12’N and between longitudes 3°23’E and 3°28’E. The study area is characterised by tropical climate with distinct wet and dry seasons. The annual rainfall ranges from 1400 mm to 1500 mm; the mean temperature is 30°C and varies from 25.7°C in July to 30.2°C in February [11]. The study area is underlain by Precambarian basement rocks (Figure 2), which are innately characterised by near-negligible permeability and low porosity. These rocks, according to a previous paper [12], were acknowledged to belong to the youngest of the three major provinces of the West African Craton. These rocks are of Precambarian age to early Palaeozoic age, which extends from the Northeastern part of Ogun State and dips towards the coast [13].
The electrical resistivity survey method using the VES method was carried out in the study area. The resistivity data were acquired using Campus Ohmega Terrameter. Thirty VES points were positioned in the study area using Schlumberger electrode configuration with half-current electrode separation (AB/2) ranging from 1 m to 132 m. The apparent resistivity values were obtained as the product of the resistance read from the resistivity meter and its corresponding geometric factor (K) for each electrode separation. The apparent resistivity data were then plotted against AB/2 on a bi-logarithm graph as sounding curves. The plotted sounding curves were interpreted manually by partial curve matching using different master curves [15]. The geoelectric parameters from the partial curve matching served as the input model for computer-assisted iteration using WINRESIST.
The values of the longitudinal unit conductance of the overburden rock units in the study area serve as the basis for the characterisation of its aquifer protective capacity. The longitudinal unit conductance gives a measure of the impermeability of the confining clay layer, which has low resistivity and low hydraulic conductivity. The protective capacity of the overburden layers in a particular area is proportional to the longitudinal unit conductance [16]. The longitudinal layer conductance (S) of the overburden at each VES station was obtained as shown in Equation (1) [17].
where
Longitudinal unit conductance/protective capacity rating (source: Oladapo and Akintorinwa [18]).
Total longitudinal unit conductance (mhos) | Rating of overburden’s aquifer protective capacity |
---|---|
<0.10 | Poor |
0.1–0.19 | Weak |
0.2–0.79 | Moderate |
0.8–4.90 | Good |
5.0–10.0 | Very good |
>10.0 | Excellent |
Olayinka [19] opined that in identifying areas of favourable aquifers within a basement terrain, the resistivity of the basement cannot be exclusively relied upon; hence, one has to consider the basement’s reflection coefficient in effectively evaluating groundwater potential in the study area. The degree of fracturing of the underlying basement is shown by the reflection coefficient [3]. The reflection coefficients (
where
In a basement terrain, groundwater yield can be grouped into high, medium and low depending on the overburden thickness and/or reflection coefficient (Table 2), as stated by Bayewu et al. [3]. The highest groundwater yield is found in areas where thick overburden overlies the fractured zone [18].
Groundwater potential yield (modified after Bayewu et al. [3]).
Overburden thickness (m) | Reflection coefficient | Groundwater yield |
---|---|---|
>13 | <0.8 | High |
>13 | >0.8 | Medium |
<13 | >0.8 | Low |
<13 | <0.8 | Very low |
The summary of the geoelectric parameters and inferred lithologies in the study area is presented in Table 3. The curve types obtained after partial curve matching ranges from the three-layer H type (66.7%), A type (3.3%) and K type (3.3%), through the four-layer KH type (23.3%), to the five-layer curve HKH type (3.3%). The predominant H type curve recorded in the study area further affirmed the findings of Oloruntola and Adeyemi [9], who recorded >72% of H curve type in the basement geological terrain at Abeokuta. Figures 3 and 4 show the typical iterated curves generated from the field measurements. The geoelectric interpretations revealed three to five layers, as follows: topsoil (23–700 Ωm); the weathered layer, which is composed of clay/sandy clay and clayey sand/sand (4–790 Ωm); underlying this layer are the fractured layer (93–437 Ωm) and the fresh basement (532–2106 Ωm). The aquifer unit in the study area is basically found in the fractured layer, the yield being dependent on the amount of clay content.
Summary of geoelectric parameters and inferred lithologies.
VES no. | No. of layers | Curve types | Resistivity (Ωm) | Thickness (m) | Depth (m) | Reflection coefficient | Inferred lithology |
---|---|---|---|---|---|---|---|
1 | 117 | 2.3 | 2.3 | Topsoil | |||
1 | 2 | H | 29 | 25.5 | 27.8 | 0.876 | Clay |
3 | 437 | Fractured basement | |||||
1 | 211 | 0.9 | 0.9 | Topsoil | |||
2 | 2 | H | 110 | 5.0 | 5.9 | 0.421 | Sandy clay |
3 | 270 | Fractured basement | |||||
1 | 324 | 1.3 | 1.3 | Topsoil | |||
3 | 2 | H | 162 | 13.1 | 14.4 | 0.437 | Clayey sand |
3 | 413 | Fractured basement | |||||
1 | 75 | 2.1 | 2.1 | Topsoil | |||
4 | 2 | H | 30 | 24.1 | 26.2 | 0.843 | Clay |
3 | 352 | Fractured basement | |||||
1 | 100 | 1.2 | 1.2 | Topsoil | |||
5 | 2 | H | 5 | 7.2 | 8.4 | 0.977 | Clay |
3 | 427 | Fractured basement | |||||
1 | 134 | 0.6 | 0.6 | Topsoil | |||
2 | 35 | 1.8 | 2.4 | Clay | |||
6 | 3 | HKH | 79 | 3.4 | 5.8 | 0.899 | Sandy clay |
4 | 17 | 8.8 | 14.6 | Clay | |||
5 | 318 | Fractured basement | |||||
1 | 94 | 0.8 | 0.8 | Topsoil | |||
7 | 2 | H | 24 | 4.9 | 5.7 | 0.778 | Clay |
3 | 192 | Fractured basement | |||||
1 | 30 | 1.0 | 1.0 | Topsoil | |||
8 | 2 | H | 4 | 2.0 | 3.0 | 0.996 | Clay |
3 | 2106 | Fresh basement | |||||
1 | 94 | 1.2 | 1.2 | Topsoil | |||
2 | 120 | 6.6 | 7.8 | Sandy clay | |||
9 | 3 | KH | 60 | 11.8 | 19.6 | 0.701 | Clay |
4 | 342 | Fractured basement | |||||
1 | 109 | 0.8 | 0.8 | Topsoil | |||
10 | 2 | 32 | 7.2 | 9 | 0.489 | Clay | |
3 | 93 | Fractured basement | |||||
1 | 26 | 0.6 | 0.6 | Topsoil | |||
11 | 2 | 655 | 5.1 | 5.7 | Sand | ||
3 | KH | 51 | 13.1 | 18.8 | 0.841 | Clay | |
4 | 592 | Fresh basement | |||||
1 | 161 | 0.9 | 0.9 | Topsoil | |||
2 | 191 | 9.1 | 10 | Clayey sand | |||
12 | 3 | KH | 114 | 8.6 | 18.6 | 0.315 | Sandy clay |
4 | 219 | Fractured basement | |||||
1 | 369 | 2.6 | 2.6 | Topsoil | |||
13 | 2 | H | 62 | 19.7 | 22.3 | 0.625 | Clay |
3 | 269 | Fractured basement | |||||
1 | 269 | 0.7 | 0.7 | Topsoil | |||
2 | 327 | 3.2 | 3.9 | Sand | |||
14 | 3 | KH | 198 | 8.8 | 12.7 | 0.534 | Clayey sand |
4 | 652 | Fresh basement | |||||
1 | 167 | 1.2 | 1.2 | Topsoil | |||
15 | 2 | H | 66 | 12.7 | 13.9 | 0.576 | Clay |
3 | 245 | Fractured basement | |||||
1 | 700 | 1.5 | 1.5 | Topsoil | |||
16 | 2 | H | 38 | 8.3 | 9.8 | 0.887 | Clay |
3 | 637 | Fresh basement | |||||
1 | 417 | 1.0 | 1.0 | Topsoil | |||
17 | 2 | 61 | 8.7 | 9.7 | 0.744 | Clay | |
3 | 415 | Fractured basement | |||||
1 | 81 | 2.3 | 2.3 | Topsoil | |||
18 | 2 | H | 27 | 10 | 12.3 | 0.880 | Clay |
3 | 422 | Fractured basement | |||||
1 | 34 | 0.5 | 0.5 | Topsoil | |||
2 | 790 | 2.5 | 3.0 | Sand | |||
19 | 3 | KH | 59 | 12.6 | 15.6 | 0.761 | Clay |
4 | 435 | Fractured basement | |||||
1 | 103 | 1.0 | 1.0 | Topsoil | |||
2 | 146 | 2.1 | 3.1 | Clayey sand | |||
20 | 3 | KH | 40 | 16.8 | 19.9 | 0.615 | Clay |
4 | 168 | Fractured basement | |||||
1 | 273 | 1.7 | 1.7 | Topsoil | |||
21 | 2 | H | 35 | 11.1 | 12.8 | 0.702 | Clay |
3 | 200 | Fractured basement | |||||
1 | 329 | 3.5 | 3.5 | Topsoil | |||
22 | 2 | H | 135 | 35.5 | 39 | 0.237 | Sandy clay |
3 | 219 | Fresh basement | |||||
1 | 226 | 2.0 | 2.0 | Topsoil | |||
23 | 2 | H | 42 | 19.4 | 21.4 | 0.921 | Clay |
3 | 1019 | Fresh basement | |||||
1 | 285 | 0.8 | 0.8 | Topsoil | |||
24 | 2 | H | 99 | 7.8 | 8.6 | 0.759 | Clay |
3 | 722 | Fresh basement | |||||
1 | 23 | 0.8 | 0.8 | Topsoil | |||
25 | 2 | A | 49 | 16.6 | 17.4 | 0.777 | Clay |
3 | 391 | Fractured basement | |||||
1 | 63 | 1.0 | 1.0 | Topsoil | |||
26 | 2 | K | 153 | 16.0 | 17.0 | 0.553 | Clayey sand |
3 | 532 | Fresh basement | |||||
1 | 149 | 1.0 | 1.0 | Topsoil | |||
27 | 2 | H | 81 | 10.4 | 11.4 | 0.643 | Clay |
3 | 373 | Fractured basement | |||||
1 | 108 | 1,6 | 1.6 | Topsoil | |||
2 | 296 | 6.7 | 8.3 | Sand | |||
28 | 3 | KH | 82 | 9.2 | 17.5 | 0.889 | Clay |
4 | 1391 | Fresh basement | |||||
1 | 161 | 1.0 | 1.0 | Topsoil | |||
29 | 2 | H | 65 | 10.0 | 11.0 | 0.522 | Clay |
3 | 207 | Fractured basement | |||||
1 | 274 | 2.6 | 2.6 | Topsoil | |||
30 | 2 | H | 40 | 19.6 | 22.2 | 0.830 | Clay |
3 | 431 | Fresh basement |
The resistivities and thicknesses of the underground layers were used to compute the longitudinal unit conductance (S) of the layers in the study area. Table 4 shows the calculated longitudinal unit conductance in mhos and the protective capacity rating for the study area. The longitudinal unit conductance values of the overburden materials in the study area ranged from 0.049720 to 1.452000 mhos. It can be observed that the protective capacity in the study area reveals poor, weak, moderate and good capacity rating. Four VES stations have poor protective capacity, six shows weak protective capacity, 17 show moderate protective capacity, while three show good protective capacity rating, with 33% of the study area falling within the poor/weak overburden protective capacity. About 57% falls within the moderate range, while 10% falls within good overburden protective capacity. The longitudinal unit conductance map of the study area in Figure 5 shows that the northeastern and northwestern parts of the study area are characterised by moderate-to-good protective capacity, and this signifies that there is a little or no infiltration due to precipitation.
Longitudinal unit conductance and aquifer protective capacity of the study area.
VES no. | Longitudinal unit conductance (mhos) | Overburden’s aquifer protective capacity rating |
---|---|---|
1 | 0.898968 | Good |
2 | 0.049720 | Poor |
3 | 0.084877 | Poor |
4 | 0.824333 | Good |
5 | 1.452000 | Good |
6 | 0.616591 | Moderate |
7 | 0.212677 | Moderate |
8 | 0.533333 | Moderate |
9 | 0.264433 | Moderate |
10 | 0.232339 | Moderate |
11 | 0.287726 | Moderate |
12 | 0.128673 | Weak |
13 | 0.324788 | Moderate |
14 | 0.056833 | Poor |
15 | 0.199610 | Moderate |
16 | 0.220564 | Moderate |
17 | 0.145021 | Weak |
18 | 0.398765 | Moderate |
19 | 0.231430 | Moderate |
20 | 0.444092 | Moderate |
21 | 0.323370 | Moderate |
22 | 0.273601 | Moderate |
23 | 0.470754 | Moderate |
24 | 0.081595 | Poor |
25 | 0.373558 | Moderate |
26 | 0.120448 | Weak |
27 | 0.135106 | Weak |
28 | 0.149645 | Weak |
29 | 0.160057 | Weak |
30 | 0.499489 | Moderate |
Figures 6 and 7 show the reflection coefficient and overburden thickness map of the study area. The reflection coefficient in the study area varies from 0.24 to 1.00. The groundwater prospects in the study area are categorised into high, medium and low potentials. In this study, zones where the overburden thickness is >13 m and the reflection coefficient is <0.8 are considered as zones with high groundwater potential, while zones with overburden thickness <13 m and reflection coefficient <0.8 are considered as zones having very low groundwater potential.
Generally, about 33% of the study area has high groundwater potential, which is restricted mostly to areas underlain by porphyritic granite and porphyroblastic gneiss, as established by a previous paper [22], while 23% of the study area has medium groundwater potential; moreover, 43% of the area has low groundwater potential. This result invariably indicates the significance of detailed groundwater survey and exploration in the study area for locating areas where successful boreholes can be sited.
Groundwater potential across the VES locations.
VES no. | Overburden thickness (m) | Reflection coefficient | Groundwater yield |
---|---|---|---|
1 | 27.8 | 0.876 | Medium |
2 | 5.9 | 0.421 | Very low |
3 | 14.4 | 0.437 | High |
4 | 26.2 | 0.843 | Medium |
5 | 8.4 | 0.977 | Low |
6 | 14.6 | 0.899 | Medium |
7 | 5.7 | 0.778 | Very low |
8 | 3.0 | 0.996 | Low |
9 | 19.6 | 0.701 | High |
10 | 9.0 | 0.489 | Very low |
11 | 18.8 | 0.841 | Medium |
12 | 18.6 | 0.315 | High |
13 | 22.3 | 0.625 | High |
14 | 12.7 | 0.534 | Very low |
15 | 13.9 | 0.576 | High |
16 | 9.8 | 0.887 | Low |
17 | 9.7 | 0.744 | Very low |
18 | 12.3 | 0.880 | Low |
19 | 15.6 | 0.761 | High |
20 | 19.9 | 0.615 | High |
21 | 12.8 | 0.702 | Very low |
22 | 39.0 | 0.237 | High |
23 | 21.4 | 0.921 | Medium |
24 | 8.6 | 0.759 | Very low |
25 | 17.4 | 0.777 | High |
26 | 17.0 | 0.553 | High |
27 | 11.4 | 0.643 | Very low |
28 | 17.5 | 0.889 | Medium |
29 | 11.0 | 0.522 | Very low |
30 | 22.2 | 0.830 | Medium |
The geoelectric investigation of the study area has revealed three to five subsurface geoelectric layers: top soil, weathered basement and fresh basement rocks. The fractured layer constitutes the sole aquifer unit in the study area. The protective capacity in the study area is more of the moderate type and is therefore not exposed to pollution. About 33% of the study area falls within the high rated groundwater potential zone, while the remaining 67% constituted the medium/low groundwater potential zone. Hence, the groundwater potential rating of the area is considered generally as medium/low. Therefore, areas for locating groundwater should be narrowed to zones of moderate/good groundwater protective capacity.