International Journal of
Physical Sciences

  • Abbreviation: Int. J. Phys. Sci.
  • Language: English
  • ISSN: 1992-1950
  • DOI: 10.5897/IJPS
  • Start Year: 2006
  • Published Articles: 2529

Full Length Research Paper

Spectral depth analysis for determining the depth to basement of magnetic source rocks over Nkalagu and Igumale areas of the Lower Benue Trough, Nigeria

Ikeh Joachims C.
  • Ikeh Joachims C.
  • Department of Industrial Physics, Enugu State University of Science and Technology, Enugu, Enugu State, Nigeria.
  • Google Scholar
Ugwu G. Z.
  • Ugwu G. Z.
  • Department of Industrial Physics, Enugu State University of Science and Technology, Enugu, Enugu State, Nigeria.
  • Google Scholar
Asielue K.
  • Asielue K.
  • School of Business Studies, Delta State Polytechnic, Ogwashi-Ukwu, Delta State, Nigeria.
  • Google Scholar

  •  Received: 05 August 2017
  •  Accepted: 07 September 2017
  •  Published: 16 October 2017


Structural interpretation of aeromagnetic data over Nkalagu-Igumale area of the Lower Benue Trough of Nigeria was carried out to determine the depth to magnetic basement and delineate the basement morphology and the structural features associated with the basin and their trends. The aeromagnetic data were subjected to series of computer based image and data enhancement techniques before spectral analysis. Results of the 2-D spectral analysis revealed two depths source models with the depth model (D1) for deep magnetic source bodies which are associated with intra-basement discontinuities and faults ranging from 2.15 to 5.25 km while the depth model (D2) of the shallow magnetic source bodies range from 0.35 to 0.99 km. From an economic viewpoint, the results indicate possible mineralization and existence of a reasonable Cretaceous sedimentary thickness in the area which is deep enough for hydrocarbon accumulation. The average sedimentary thickness obtained in the area is 3.75 km.

Key words: Benue Trough, depth to magnetic basement, aeromagnetic data, spectral analysis, intrusive bodies.


Magnetic method is a measurement of the earth’s magnetic field intensity, typically involving the total magnetic field and/or vertical magnetic gradient, horizontal or vertical component or horizontal gradient of the magnetic field (Biswas et al., 2017). Anomalies in the earth’s magnetic field are caused by induced or remnant magnetism. Induced magnetic anomalies are the result of secondary magnetization induced  in  a  ferrous  body  by the earth’s magnetic field. The shape, dimensions and amplitude of an induced magnetic anomaly is a function of the orientation, geometry, size, depth and magnetic susceptibility of the body, as well as the intensity and inclination of the earth’s magnetic field in the survey area (Biswas, 2016; Biswas and Acharya, 2016). For exploration purposes, both ground and aeromagnetic data have been used to investigate the presence of mineral deposits in combination with gravity. In the mining industry, both gravity and magnetic methods are still widely used as exploration tools for mapping subsurface geology and to estimate ore reserves for some massive ore bodies (Mandal et al., 2015; Biswas and Sharma, 2016). Nowadays, aeromagnetic surveys are very useful for determining depths to magnetic source rocks of the anomalies over an area, and utilized the principle that the magnetic field measured at the surface can be considered an integral of magnetic signatures from all depths. These anomalies arise as a result of the interactions between the magnetic field and the rocks of the earth crust and also from secondary mineralization along the fault planes in many sedimentary basins. The average depth to sources of ensembles (sedimentary thickness) across the geological area can be obtained from the energy (power) spectrum of the surface field when plotted against frequency in a logarithmic scale. Interpretation of aeromagnetic data can be done quantitatively and qualitatively. Quantitative interpretation involves making numerical estimates of the depth and dimensions of the sources of the anomalies and this often takes the form of modelling of the sources which could  in theory, replicate the anomalies recorded in the survey (Biswas, 2015, 2016; Biswas et al., 2017). The computational efficiency and storage capacity of modern digital computers now enable the utilization of a wide range of sophisticated mathematical techniques for data processing and interpretation. Fourier spectral analysis as a mathematical technique has attracted increasing attention in recent years as a powerful tool for the interpretation of potential field data (magnetic anomalies).
The study area lies within the southern portion of the Lower Benue Trough and consists of Sheet 302 (Nkalagu sheet) and Sheet 288 (Igumale sheet) bounded within longitudes 7° 30’ and 8° 00’ East and latitudes 6° 30’ and 7° 00’ North and covering an area of about 6000 km2. Figure 1 shows the map of Nigeria and location of the study area.
The Lower Benue Trough as a rift system and aulacogen (Olade, 1976) is an area of importance in terms of economic mineral deposits. Consequently, this has inspired a lot of interest on geophysical investigations by many researchers (Ofoegbu and Onuoha, 1991; Obi et al., 2010; Ugwu and Ezema, 2012) to re-examine the economic potentials  of  this  mineral  belt  using  different geophysical methods. Previous magnetic anomaly studies over the Lower Benue Trough has shown that the anomalies in the area can best be explained in terms of the combined effects of deep lying basement and intermediate intrusions at both shallow and deep depths (Ofoegbu, 1984, 1985; Ofoegbu and Onuoha, 1991; Ugwu and Ezema, 2012). Ofoegbu (1985) estimated the sedimentary thickness in the Lower and Middle Benue Trough to vary between 0.5 and 7.0 km. Some of the anomalies were interpreted in terms of dykes and volcanic plugs (Ofoegbu, 1984; Ugwu and Ezema, 2012). Previous studies also suggested that Benue Trough has a typical network of fractures and lineaments with dominant trends of NE-SW direction. The intrusive bodies are of variable thickness and different directions of magnetization, which suggest that although they are derived from a common basic mantle material, they were probably emplaced at different polarity epochs (Ofoegbu, 1984).
To contribute to a better understanding of the depth and nature of mineralization, geology and development of the tectonic history in the Lower Benue Trough, we have considered the use of spectral analysis of the aeromagnetic data over Nkalagu and Igumale areas to compute for depths to magnetic sources.
Geology of the study area
The geology and evolution of the lower Benue Trough is now fairly well documented (Wright, 1968, 1976; Nwachukwu, 1972; Olade, 1975; Ofoegbu,  1985;  Obaje, 2009). The lower Benue Trough which is underlain by a thick sedimentary sequence was formed as a result of series of tectonics and repetitive sedimentation in the Cretaceous time. The depositional history of the Benue Trough is characterized by phases of marine regression and transgression (Reyment, 1965; Short and Stauble, 1967; Murat, 1972). The major component units of the lower Benue Trough include the Anambra Basin, the Abakaliki Anticlinorium and the Afikpo Syncline. The oldest sediment of the sequence belongs to the Asu River Group which unconformably overlies the Precambrian basement complex that is made up of granitic and magmatic rocks (Ofoegbu and Onuoha, 1991). The Asu River Group whose type outcrops in Abakaliki has an estimated thickness of about 2000 m (Ofoegbu, 1985) and is Albian to Cenomanian. It comprises of argillaceous sandy shales, laminated sandstone units and minor limestones with an interfingering of magnetic volcanics (Nwachukwu, 1972). The shales are fissile and highly fractured. Deposited on top of these Asu River Group sediments in the area are the Upper Cretaceous Eze-Aku shales. The Turonian Eze-Aku shales consist of nearly 1000 m of calcareous flaky shales and siltstones (Reyment, 1965). The geological map of the study area is shown in Figure 2.


Two aeromagnetic maps were acquired from the Nigerian Geological Survey Agency (NGSA). They are Sheets 302 (Nkalagu) and 288 (Igumale). The data were acquired along a flight line trend of 135° with a spacing of 500 m and an average  flight  elevation  of 75 m, while the tie lines trend and spacing are 225° and 500 m, respectively. The data were made available in digital form on scale of 1:100,000.
The geomagnetic gradient was removed from the data using the International Geomagnetic Reference Field (IGRF). The data sheets of Nkalagu (302) and Igumale (288) were combined using Microsoft Excel 2010 to form the study area covering a total area of about 6000 km2. The combined sheet was used to produce the total magnetic field intensity (TMI) map (Figure 3) of the study area employing WINGLINK Software and Surfer 10 binary grids. The regional-residual separation was done by polynomial fitting. Here, all the regional fields were calculated as two-dimensional first degree polynomial surfaces. Figure 4 shows the residual anomaly map of the study area after the regional-residual separation. The residual field was also gridded into horizontal and vertical derivatives (Figures 5 and 6). The 3-D surface anomaly map of the study area is also shown in Figure 7.
In order to calculate the depth to basement of the study area using spectral analysis, the study area was divided into six profiles with 18 cells containing 18 × 18 data points. In doing this, it was ensured that essential parts of each anomaly were not cut off in the cells as each cell was made to contain more than one maximum, as suggested by Hahn et al. (1976). To achieve this, the cells were made to overlap each other using Excel 2010 program  by  applying inequality formulas in which degree of overlap was insignificant. The Excel software which has the Fourier mathematical algorithm was used to compute the Fast Fourier Transform (FFT) magnitude and FFT frequency. The energy (magnitude) which is in logarithmic scale was plotted against logarithm frequency to produce a decay-curve segment which decreases in slope with increasing frequency. The depths to sources of ensembles were obtained by manually fitting a straight line to each linear interval of the logarithmic energy-decay curve which  yields  a  negative  slope.  The  line  of  the  first linear interval was fit from the deviation of the decay curve from the vertical axes of the logarithmic energy. The second line of the second linear interval was drawn from the deviation of the decay curve from the line of the first linear interval. This was done sequentially such that there is no other possible linear interval between the first and second linear intervals.
Smoothing was done to the signals to filter and determine a trend (linear curvature) on the decay curve by removing high frequency noise. The effect of this filter was to average the value at a given  grid node together with the values at its eight nearest neighbours to give a new value at the centre point. The smoothing operator was moved repeatedly row-by-row over the whole grid to produce new gridded values. The problems of aliasing effect and Gibb’s phenomenon which are normally encountered in the course of performing spectral analysis were taken care of by the use of small sampling intervals to reduce frequencies greater than Nyquist frequency.


The spectral plots of  logarithmic  energy  against frequency for the six representative cells of the profiles are shown in Figures 8 to 13. Two-depth models have been estimated for all the cells and these values are summarized in Table 1. The slope is given by the equation:
While the depth to slope relationship is given by D= -S/4  when frequency is in cycle/km (Bonde et al., 2014).


A visual inspection (subjective) of the total magnetic intensity (TMI) and residual maps show that the contour lines of the central and southern parts are widely spaced, suggesting that the depths to magnetic basement in these areas are relatively high. At the northern part of the maps, the contour lines are more closely together, suggesting that the depth to basement is shallow at this part. Also, there are spikes at the northeastern part which shares a geologic boundary with Ejekwe (very similar to Abakaliki, geologically) indicating the presence of magnetic mineral rocks or intrusive bodies within the sedimentary cover. So it is possible that the existence of intrusive bodies within the Abakaliki anticlinorium of the Benue Trough (Ofoegbu and Onuoha, 1991; Ugwu and Ezema, 2012; Abdulahi et al., 2014) extends to Igumale. The distribution of high magnetic intensities from just above the central part toward the southern part of the study area suggests the existence of deep penetrating fractures within the area while the low magnetic intensity indicated by closures and the linear sub parallel orientations of the contours at the uppermost northern part of the map suggest a shallow geologic structure. Structural trends within the study area as shown by delineation in the residual map are mainly NE-SW and NW-SE directions which agree  with  the  fault  orientation within the Benue Trough, with the NE-SW trends being dominant.
The 3-D surface map of the study area (Figure 7) shows high basement, relief indicating folded and undulating topography, which can be  interpreted  as  part of the Abakaliki folded belt. The results of the vertical derivative anomaly maps indicate the lithologic boundaries between different formations underlying the area. The horizontal derivative map shows the possibility of faults or local fractured zones passing through the north central and central part of the study area.
Results of the spectral analysis of the aeromagnetic data over the study area indicated a two-depth source model. The shallow  depth  model  is  between  0.35  and 0.99 km which probably indicates the presence of intrusive bodies. The deep source lies at a depth that varies between 2.15 and 5.25 km. These deep sources represented by the first segment of the spectrum in all the cells of the six profiles reflect the Precambrian basement of the study area. Profile two shows the thickest sedimentary cover in the Nkalagu area (cells 4, 5 and 6 of Table 1). This result closely agrees with the results from other  aeromagnetic   works   which   had   indicated depth (sedimentary thickness) of 1.5 to 4 km in Nkalagu (Nur et al., 1994; Obi et al., 2010). The range of the average profile depths implies a nearly evenly distributed sedimentary thickness within the study area and ranges between 3.5 and 4.2 km.


Quantitative interpretation of the aeromagnetic data over Nkalagu and Igumale areas of the Lower Benue Trough has been successfully carried out using spectral analysis to determine depth to the magnetic basement. Structural interpretation of the residual anomaly map was also used to delineate the basement morphology, relief, and the structural features associated with the basin and their trends. The results of the study have shown that the area is characterized with an average sedimentary thickness of 3.75 km for the deep source model. From economic point of view, the results indicate possible mineralization within Igumale area due to presence of a few intrusive bodies within the area. The results also show very high possibility for hydrocarbon occurrence due to the existence of folded basement, faults/fractures capable of trapping hydrocarbons. The occurrence of reasonable Cretaceous sedimentary thickness in the area supports the hydrocarbon potential of the area.


The authors have not declared any conflict of interests.


Abdulahi UA, Ugwu GZ, Ezema PO (2014). Magnetic Exploration of the Upper and Lower Benue Trough for Metallic Deposits and Hydrocarbons Using 2D/3D. J Nat. Sci. Res. 4(20):41-46.


Biswas A (2015). Interpretation of residual gravity anomaly caused by a simple shaped body using very fast simulated annealing global optimization. Geosci. Front. 6(6):875-893.


Biswas A (2016). Interpretation of gravity and magnetic anomaly over thin sheet-type structure using very fast simulated annealing global optimization technique. Model. Earth Syst. Environ. 2(1):1-30.


Biswas A, Acharya T (2016). A very fast simulated annealing method for inversion of magnetic anomaly over semi infinite vertical rod- type structure. Model. Earth Syst. Environ. 2(4):198.


Biswas A, Sharma SP (2016). Integrated geophysical studies to elicit the structure associated with Uranium mineralization around South Purulia Shear Zone, India: A Review. Ore Geol. Rev. 72:1307-1326.


Biswas A, Parija MP, Kumar S (2017). Global nonlinear optimization for the interpretation of source parameters from total gradient of gravity and magnetic anomalies caused by thin dyke. Ann. Geophys. 60(2):G0218, 1-17.


Bonde DS, Udensi EE, Rai JK (2014). Spectral Depth Analysis of Sokoto Basin. ISOR J. Appl. Phys. 6:15-21.


Hahn A, Kind EG, Mishra DC (1976). Depth estimation of magnetic sources by means of Fourier amplitude spectra. Geophys. Prospecting 24:287-308.


Jatau BS, Nandom A (2013). Morphology of parts of the middle Benue Trough of Nigeria from spectral analysis of aeromagnetic data. Int. J. Environ. Chem. Ecol. Geol. Geophys. Eng. 7(9).


Mandal A, Mohanty WK, Sharma SP, Bismas A, Sen J, Bhatt AK (2015). Geophysical signatures of uranium mineralization and its subsurface validation at Beldih, Purulia District, West Bengal, India: A case study. Geophys. Prospecting 63:713-726.


Murat C (1972). Stratigraphy and paleogeography of the Cretaceous and Lower Tertiary in South- Eastern Nigeria. In: Dessauvagie TFJ, Whiteman AJ (eds.), African Geology. Ibadan University Press. 251-266.


Nwachukwu SO (1972). The tectonic evolution of the southern portion of the Benue Trough, Nigeria. J. Min. Geol. 109:411-419.


Nur MA, Onuoha KM, Ofoegbu CO (1994). Spectral analysis of aeromagnetic data over the Benue Trough. J. Min. Geol. 30(2):211-217.


Obaje NG (2009). The Benue Trough. Geology and Mineral Resources of Nigeria, Spingerlink. 57-66.


Obi DA, Okereke CS, Obei BC, George AM (2010). Aeromagnetic modelling of subsurface intrusive and its implications on hydrocarbon evaluation of the Lower Benue Trough, Nigeria. Eur. J. Sci. Res. 47(3): 347-361.


Ofoegbu CO (1984). Interpretation of magnetic anomalies over the Lower and Middle Benue Trough of Nigeria. Geophys. J. R. Astron. Soc. 79:813-823.


Ofoegbu CO (1985). A review of the geology of the Benue Trough of Nigeria. J. Afr. Earth Sci. 3:293-296.


Ofoegbu CO, Onuoha KM (1991). Analysis of magnetic data over the Abakaliki anticlinorium of The Lower Benue Trough, Nigeria. Marine and Petroleum Geology. 8:174-183.


Olade MA (1975). Evolution of Nigeria's Benue Trough (Aulacogen); a tectonic model. Geol. Mag. 112:575-583.


Olade MA (1976). On the genesis of the lead-zinc deposits in Nigeria' Benue Rift (aulacogen). A re-interpretation. J. Min. Geol. 13(2):20-27.


Onuba LN, Onwemesi AG, Egboka BC, Anudu GK, Omali AA (2013). Review of Hydrocarbon prospects in the Lower Benue Trough, Nigeria: Adapted from extended abstract prepared in conjunction with oral presentation at AAPG Annual Convention and Exhibition, Pittsburgh, Pennsylvania.


Reyment RA (1965). Aspects of Geology of Nigeria, Ibadan University Press, Ibadan Nigeria.


Short KC, Stauble AJ (1967). Outline of Geology of the Niger Delta. Assoc. Petr. Geol. 5(6):761-779.


Ugwu GZ, Ezema PO (2012). Forward and Inverse Modelling of Aeromagnetic Anomalies over Abakaliki and Nkalagu areas of the lower Benue Trough, Nigeria. Int. Res. J. Geol. Min. 2(7):199-204.


Wright JB (1968). South Atlantic continental drift and the Benue Trough. Tectonophysics 6(4):301-310.


Wright JB (1976). Origins of the Benue Trough: A critical review. In: Geology of Nigeria, Kogbe CA (eds.). Elizabethan Publ. Co. Lagos. pp. 309-318.