Intersection of lineaments for groundwater prospect analysis using satellite remotely sensed and aeromagnetic dataset around Ibodi, Southwestern Nigeria

This research investigates groundwater potential of Ibodi and its environs using satellite imagery remotely sensed and aeromagnetic dataset. This hydrogeological investigation involved the extraction of lineaments from Landsat Thematic Mapper (TM) satellite imagery and aeromagnetic lineament of the area. The processed image revealed lineaments trending in approximately NE-SW direction; also, Landsat and aeromagnetic lineament trends tend to agree in the study area, suggesting real continuous fractures at depth. Hydrogeological maps, based on lineament and lineament intersection were produced from the generalized lineament trends in the area. The lineament intersection density map identified cluster zones in the range of 0 to 1.17, 1.17 to 2.34, 2.34 to 3.52, 3.52 to 4.69 and 4.69 to 5.86 kg/kg 2 , respectively. The map showed that the concentrations of lineament intersection nodes dominate the eastern and north-eastern parts of the study area while other areas are characterized by scanty or no lineament intersections nodes. The lineament analysis has been effectively done in a geographic information system (GIS) environment. Results from the vertical electrical sounding (VES) data reveals that the areas where the lineament intersected has tendency of very high groundwater prospect which reflects that the area will be good for groundwater development. The study thus displays that the remote sensing and aeromagnetic technique is capable of extracting lineament trends in an inaccessible tropical forest and has led to the delineation of areas where groundwater occurrences are most promising for sustainable supply.


INTRODUCTION
Groundwater is one of the most valuable natural resources on the earth surface and serves as main sources of drinking water.Basement complex have problem of potable groundwater supply due to the *Corresponding author.E-mail: bussytex4peace44@gmail.com.
Author(s) agree that this article remain permanently open access under the terms of the Creative Commons Attribution License 4.0 International License crystalline nature of the underlying rock which lack primary porosity (Olorunfemi and Fasuyi, 1993).
Groundwater storage capacity in those areas depends on the depth of weathering and the intensity of fracturing of the underlying rock.For basement rocks to become good aquifers, they must be highly fractured and highly weathered.Thickness of the weathered overburden and fractured zone determined the nature and intensity of hydrodynamic activities within the usually discrete bodies of aquifer in the terrain (Amudu et al., 2008;Amadi and Olasehinde, 2010).
Groundwater in crystalline rocks that have no intergranular porosity moves in a connected fracture network; but far from all fractures are permeable, and fracture permeability varies considerably (Banks et al., 1996).It is generally recognized that, in the prospect area of faults and fracture lineaments, the fault core and central zone have low permeability while the outer damage zone has enhanced permeability compared with the surroundings (Caine et al., 1996;Evans et al., 1997;Henriksen and Braathen, 2006).Lineament mapping was used long before this work in other geological applications and the first usage of the term lineament in geology is probably from a paper by Hobbs (1904Hobbs ( , 1912)), who defined lineaments as significant lines of landscape caused by joints and faults, revealing the architecture of the rock basement.This was later used by O' Leary et al. (1976) as a basis for developed definitions.Lineament is a mappable linear or curvilinear feature of a surface whose parts align in a straight or slightly curving relationship (Hung et al., 2005).They may be an expression of a fault, joints or other line weakness.These features are mappable at various scales, from local to regional, and can be utilized in mineral, hydrocarbons, and groundwater explorations.This interest has grown most rapidly in geological studies since the introduction of aerial photographs and satellite images.Generally, lineaments are underlain by zones of localized weathering and increased permeability and porosity.Meanwhile, some researchers studied relationships between groundwater productivity and the number of lineaments within specifically designated areas or lineament density rather than the lineament itself (Hardcastle, 1995).
This research focuses on intersecting landsat lineament and aeromagnetic lineament for regional groundwater prospect evaluation.Lineament intersections represent nodal point for two or more lineament lines (Bayowa et al., 2014).The nodes of intersections are the point of interest which indicate site for appreciable groundwater accumulation.

Study location and geology
Ibodi and its environs comprises sixty-eight villages which are located on the crystalline basement complex and falls within latitude 7° 32′ 20″ to 7° 38′ 20″ North and longitude 4° 36′ 20″ to 4° 44′ 20″ East (Figure 1).This area falls within the 1:50,000 topographic map of Ilesa sheet 234 SW.Major and minor road linkages characterize the study area linking both towns and villages together.In general, the location can be said to be fairly accessible with footpaths in areas where there are neither major nor minor roads.Regionally, the study area belong to the Southwestern Nigeria basement complex comprising migmatite-gneiss complex, metaigneous rock such as pelitic schist, quartzite, amphibolites, charnokitic rocks, older granite and unmetamorphorsed dolerite dykes.The rock sequence consists of basically weathered quartzite of older granite.The basement complex rocks of Nigeria are made up of heterogeneous assemblages (Rahaman, 1976).The geology of Ibodi and its environs is mainly Precambrian Basement rock (Figure 2) because of the tectonics and metamorphic changes that has occurred in the area (Olade, 1978).Schist and amphibolites, magnetitegneiss undifferentiated and amphibolites stands out in the Basement rock of the investigated area.

METHODOLOGY
The work utilized remote sensing, aeromagnetic and electrical resistivity data in the assessment of groundwater prospect around Ibodi, Osun State.
The remote sensing data (Landsat imagery) were acquired from the Global Land Cover Facility homepage and the automatic lineament extraction process was carried out using "ENVI" 4.5, "Geomatical'' pcl and "ArcGIS'' 10.2 software which was used for digital image processing and lineament extraction from the satellite imageries.The aeromagnetic data was acquired from Nigerian Geological Survey Agency and the data set gridded using minimum curvature gridding method with 50 m cell size to produce total magnetic intensity map.Data enhancement techniques involving Butterworth filter, reduction to equator and derivative filters were applied to the magnetic intensity map using Oasis Montaj 6.4.2 software package (Montaj TM Tutorial, 2004).3D Euler deconvolution data assisted in determining the locations and depths of the geologic sources.Structural index of 1.0 was found appropriate for delineation of structures that are favourable for groundwater accumulation.Lineaments were generated by georeferencing and digitizing 3D Euler deconvolution results.A total of nine Vertical Electrical Sounding involving Schlumberger arrays were used to acquire depth sounding data in the investigated area with the apparent resistivity values obtained from the area plotted on bi-log graph against the current electrode separation spacing.The initial quantitative interpretations were made using partial curve matching technique in which the field curves produced or generated were matched segment by segment with the appropriate master curves and auxiliary curves.The resistivities and thicknesses of the various layers were improved upon by employing an automatic iterative computer program following the main ideas of Zohdy and Martin (1993).The WINRESIST computer software was employed for carrying out the iteration and inversion processes.Intersection of landsat lineament and aeromagnetic lineament density were used to determine the groundwater prospect of the area.A total of nine (9) vertical electrical sounding obtained within the study area were used to validate lineament intersections in order to ascertain the significances in term of groundwater prospect.

Lineament analysis from landsat imagery
Results obtained from the satellite image interpretation (Figure 3) are discussed to demonstrate the usefulness of remote sensing in lineament mapping and analysis, that is, to delineate zones that are prone to groundwater development in the study area.The map shows that the northeastern region is characterized by very high concentration of lineaments, whereas the central parts of the study area are characterized by moderately

Landsat lineament on topographical and geological map
Figures 5 and 6 show landsat lineament on topographical and geological maps.Figure 5 shows that the northeastern parts, majority of the lineaments fall within the high hill with elevation between 380 m to 420 m and also at the southwestern part.Areas with high lineament density excluding (the residual hill environment) are good for groundwater development.Figure 6 show that there are more concentrations of lineaments on schist and amphibolites at the northeastern parts followed by amphibolites at the southwestern part while there exist scanty or low concentration of lineaments on gneiss and migmatite undifferentiated.This implies that there are more fracturing on schist and amphibolites, followed by amphibolites and gneisss and migmatite undifferentiated.Regarding groundwater exploration, these aforementioned areas may have high groundwater potentials due to their high concentration of lineaments, that is, since groundwater occurs within faults and fractures in the basement rock.

Landsat lineament density map
The lineament density of Ibodi and its environ was computed based on the number of lineaments per unit area (kg/kg 2 ) of grid, for quick graphical assessment of the lineament density values of the area.Areas with high lineament density excluding (the residual hill environment) are good for groundwater development.
Figure 7 shows five different hydrogeological potential zones distributed as patches in the study area.The zones are summarized in Table 1.The lineament density map shows that the lineament density is high around areas like Ibodi, Iloki, Iwori, Ijano, Ilotin, Aye, Igila, Agbao, Araromi and Oja Oko when compared to other areas in the study area.

Lineament analysis using 3D Euler deconvolution
3D Euler solution using structural index (S.I = 1) was applied to the aeromagnetic data for lineament analysis (Figure 8).The delineated solutions were extracted and digitized to produce the lineaments from magnetic data (Figure 9).The rose (azimuth-frequency) diagram of the lineaments (Figure 10) was prepared from the extracted lineaments on the aeromagnetic structures which shows that there are two predominant sets of lineament which are closely related to tectonic activities such as fractures,

Aeromagnetic lineaments on topographical and geological map
Figures 11 and 12 displays aeromagnetic lineament on topographical and geological map. Figure 11 reveals that the majority of the lineaments fall within the plain land unlike surface lineaments (Landsat Lineament) which fall within the high hill.This implies that we cannot base our presumption on surface lineament without ascertaining it with subsurface lineaments in order to have clearer pictures of the groundwater prospect within the study area.Figure 12 shows that there are more concentration of lineaments on schist and amphibolites at the southeastern parts which go along with the landsat lineaments on geological map followed by gneiss and migmatite undifferentiated at the northwestern parts, while showing scanty or low concentration of lineaments on amphibolite.This implies that, there are more fracturing on schist and amphibolites.Regarding groundwater exploration, these aforementioned areas may have high groundwater potentials due to their high concentration of lineaments, that is, since groundwater occurs within faults and fractures in the basement rocks.

Aeromagnetic lineament density map
Lineament density map is one of the important maps prepared from the lineaments, which are critically used in groundwater studies related to hard rock terrain, Figure 13 displays the aeromagnetic lineament density map of the study area.The map reveals five different hydrogeological potential zones distributed as patches in the study area.The zones are summarized in Table 2.The lineament density map reflects that the lineament density is high around areas like Ibodi, Iworo Odo and Sefari, when compared to other areas in the study area.Regarding groundwater exploration, these aforementioned areas may have high groundwater potentials due to their high concentration of lineaments, that is, since groundwater occurs within faults and fractures in the basement rocks.

Intersection of aeromagnetic and landsat lineaments
Figure 14 indicates the lineament intersection map

Lineaments intersection on topographical and geological map
Figures 16 and 17 display lineament intersection on topographical and geological map. Figure 16 indicates that the lineament intersection nodes concentrated within northeastern part of the study area fall within the plain land while other areas are characterized by scanty or no lineaments intersection nodes.This implies that we cannot base our presumption on surface lineament without ascertain it with subsurface lineaments in order to have clearer pictures of the groundwater prospect within the study area.Figure 17 shows that there are more lineaments intersection on schist and amphibolites at the northeastern parts followed by gneiss and migmatite undifferentiated while having scanty or low lineaments intersection on amphibolite.This implies that there are more areas with high lineament intersections which are good for groundwater development on schist and amphibolites and gniess and migmatite undifferentiated while having low lineament intersection on amphibolites, which implies low groundwater development.Regarding groundwater exploration, these aforementioned areas may have high groundwater potentials due to their high concentration of lineament intersections, that is, since groundwater occurs within faults and fractures in the basement rocks.

Superposition of vertical electrical sounding (VES) on lineament intersections
Figure 18 indicates the superposition of vertical electrical sounding (VES) on lineament intersections.Nine VES data were obtained within the investigated area.Figure 19(a to i) shows typical iterated VES data curves and the estimated geoelectric parameters.Four (4) curve types were identified, viz; H, QH, KH and HKH, where curve type KH has the highest occurrences.The curve types were ground on the basis of the aquifer types.The curve   (confined).The fractured zone is concealed beneath a fresh basement rock, in which there is no hydraulic communication between the water contains in the weathered layer and the fractured zone, at least within the immediate finite of the abstraction point.Aquifer in KH corresponds to the weathered/fractured (unconfined)/fractured (confined) aquifer.This aquifer types produces very high groundwater yield particularly when the density of the fractures is high.However, area underlain by the weathered/fractured basement aquifer is more promising in terms of groundwater development than the weathered layer aquifer and also the points of intersection in lineament intersections are the point of interest for groundwater development.The VES interpretation results Figure 19a to i of the study area show a total of three to five subsurface layers namely: the topsoil, weathered layer, weathered/fractured basement and fresh basement.The first layer which is the topsoil has resistivity varying from 81 to 4333 Ωm and thickness ranging from 0.6 to 2.1 m.This layer is relatively thin and with high resistivity.The resistivity of the second layer which is the weathered zone varies from 38 to 700 Ωm while its thickness varies from 1.7 to 28.8 m.The third layer constitutes the weathered/fractured basement which is also being referred to as aquifer layer, and are interpreted as saturated of their relative lower resistivity ranging from 130 to 432 Ωm.The thickness ranges from 2.9 to 52.2 m while the fourth layer is the fresh basement with resistivity values ranging from 455 to 5286 Ωm.The result from the VES data reflects that the areas where the lineament intersected has tendency of very high groundwater prospect.This illustrates that the area will be good for groundwater development.of the study area.Lineament intersection density map is a measure of cluster of linear features in a particular area.The peaks in the lineament intersection density contour maps are believed to be potential zones for

Conclusions
In this research, groundwater potential of Ibodi and its environs have been investigated using satellite imagery remotely sensed and aeromagnetic dataset.The hydrogeological investigation involved the extraction of lineaments from Landsat Thematic Mapper (TM) satellite imagery and aeromagnetic lineament from the investigated area.The processed image displays the lineaments trending NE-SW directions.Landsat and aeromagnetic lineament trends tend to agree in the study area, implying that these lineaments reflect real continuous fractures at depth.Hydrogeological maps based on lineament and lineament intersections were produced from the generalized lineament trends in the area.The lineament intersection density map identified cluster zones in the range of 0 to 1.17, 1.17 to 2.34, 2.34 to 3.52, 3.52 to 4.69 and 4.69 to 5.86 kg/kg 2 .The map revealed that the concentrations of lineament intersection nodes dominate the eastern and north-eastern parts of the study area while other areas are characterized by scanty or no lineament intersections nodes.Ibodi falls within the high lineament intersections density zones; however, localities like Olorombo, Iwaro-Odo reflected medium lineament intersection density.Agbao, Oke Agbe, Aye, Onigbodi, Ilosi, Sefari, Iwori, Oja Oko and Araromi fall within low lineament intersection density.On the other hand, Odo-Iju, Iyere, Iyemogun, Oloyin, Aye-Ile, Imosan, Shasha and Ijana areas have very low lineament intersection density values.Remote sensing and aeromagnetic dataset has proved to be a useful tool in lineament identification and mapping.The result from the depth sounding data interpretation indicates four curve types which are H, QH, KH and HKH, where curve type KH has the highest occurrence.The results from the vertical electrical sounding data reveal that the areas where the lineament intersected has tendency of very high groundwater prospect.This indicates that the area will be good for groundwater development.The study has led to the delineation of areas where groundwater occurrences are most promising for sustainable supply, suggesting that an area with high concentrations of lineament intersections has a high tendency for groundwater prospecting.This study demonstrates the application of remotely sensed data and aeromagnetic dataset for lineament interpretation in a hard rock hydrogeological environment.The lineament analysis has been effectively done in a GIS environment.The results from the study shows that the remote sensing and aeromagnetic technique is capable of extracting lineament trends in an inaccessible tropical forest.

Figure 1 .
Figure 1.Location map of Ibodi and its environ showing the study area (after Osun State Ministry of Lands and Physical Planning 1965.Topographic Map of Ilesha SW sheet 234).

Figure 2 .
Figure 2. Geological map of Ibodi and its environs (after Osun State Ministry of Lands and Physical Planning 1980.Geological Survey of Ilesha Iwo Shist 60).

Figure 3 .
Figure 3. Landsat lineaments map of the study area.

Figure 5 .
Figure 5. Landsat lineaments on topographical map of the study area.

Figure 6 .
Figure 6.Landsat lineaments on geological map of the study area.

Figure 7 .
Figure 7. Landsat lineament density map of the study area.

Figure 8 .
Figure 8. 3D Euler solution from the analysis the aeromagnetic data over the study area with structural index (S.I = 1).

Figure 9 .
Figure 9. Aeromagnetic lineaments of the study area.

Figure 11 .
Figure 11.Aeromagnetic lineament on topographical map of the study area.

Figure 12 .
Figure12.Aeromagnetic lineament on geological map of the study area.

Figure 13 .
Figure 13.Aeromagnetic lineament density map of the study area.

Figure 16 .
Figure 16.Lineament intersections on topographical map of the investigated area.

Figure 17 .
Figure 17.Lineament intersection on geological map of the investigated area.

Figure 20
Figure 20 exhibit the lineament intersection density map

Figure 20 .
Figure 20.Lineament intersection density map of the study area.

Table 1 .
Groundwater prospecting of the study area based on landsat lineament density.

Table 2 .
Groundwater prospecting of the study area based on aeromagnetic lineament density.
Figure 14.Intersection lineaments map of the study area.

Table 3 .
Groundwater prospect of the study area based on lineament intersection density.