Nwo 1.
Nwo 1.
Analysis of aeromagnetic data over Igabi, Kaduna state, which is part of the Basement Complex of
Nigeria, was carried out to highlight linear structures, and to infer the effects of such features on regional
groundwater, solid mineral availability, as well as hydrocarbon potential. The study area lies within
latitudes 100 33’ - 100 35’ N and longitudes 7 0 26’ - 7 0 46’ E, and it is located in the Younger Granite
Province, known to be characterized by granitic plutonism, igneous intrusion and magmatism. Data was
processed using Oasis Montaj software and ArcGIS. Structural trend characteristics of the study area, and
depth to basement of magnetic anomalies were highlighted using the Euler deconvolution methods.
Deconvolution of aeromagnetic profiles, were utilized to determine the depth to magnetic sources. These
depths were found to range from (0 – 3000m). The study area was found to be characterized by major and
minor faults, with a dextral movement in the N-E, N-W and S-E directions. Predominant lineament trends
were also observed to be in the N-E, and S-W directions. Areas overlain by a major crossover of
lineaments implying the presence of a weak zone was found to correspond to areas with low magnetic
intensities on the residual aeromagnetic field map of the study area. These were suggested as regions for
possible groundwater targeting. Findings suggested that study area is also good for mineral prospecting,
but not favorable for hydrocarbon prospecting.
CHAPTER ONE
INTRODUCTION
1.1 Preamble
Obtaining data about the subsurface either from the extrapolation of surface geological
information or from direct measurement of physical properties of rocks and earth materials
by geophysical methods is still the most potent fundamental means of deriving knowledge
about structural disposition of an area (Reynolds, 2011; Parsekian et al., 2015).
Various geophysical methods are often applied to solve specific geological problems
especially when the nature of the study requires rapidity and less disturbance of the
subsurface (Lowe and Zaccheo, 1991; O’Neil, 2001; Holden et al., 2002; Rana, 2019).
Geological problems often requiring geophysical solutions include exploration and mining
studies, groundwater/hydrological studies, environmental and nuclear studies, and
engineering site investigations and developmental planning (Fairbridge et al., 1998; National
Research Council, 2004; Gore and Davies, 2013).
Engineering site investigations and developmental planning is an aspect of geophysics, which
though have been practiced for decades, is just gaining recognition and patronage in the
developing economies of the world (Ribes and Finholt, 2009; Curtis and Low, 2016). It is
gradually becoming obvious that modern civil engineering practice needs information about
the conditions of the subsurface (geologic information) as much as it needs accurate
structural detailing which in time past has usually been their singular focus in order to
maximise land resources (Walters, 2008; Griffiths and King, 2013; Dearman, 2013).
Some of the most effective geophysical methods for engineering site investigations and
developmental planning include seismic refraction method (Cardarelli et al., 2010; Adamo et
al., 2021), electrical methods (Oyedele and Olorede, 2010; Cardarelli et al., 2010; Hossain et
al., 2010), magnetic methods (Adagunodo et al., 2015; Balogun, 2019; Oyeniyi et al., 2016)
and gravity method (McDowell et al., 2002; Davis et al., 2008) e.t.c.
Magnetic and electrical (electrical resistivity) methods remain two of the most patronised
methods in engineering site investigation and developmental planning because of its
affordable instrumentation and ease of use (Bogoslovsky and Ogilvy, 1997; Robinson et al.,
2008; Cardarelli et al., 2010). Magnetic survey data often depict differences in the earth's
magnetic field caused by the magnetic characteristics of the underlying rocks (e.g. magnetic
susceptibilities). Magnetic susceptibility is lowest in sedimentary rocks, whereas it is highest
in metamorphic and acidic igneous rocks, intermediate and basic igneous rocks (Telford et
al., 1990; Kearey et al., 2002). Magnetic survey can be land borne (ground magnetic survey),
ship borne (ocean magnetic survey) or air borne (aeromagnetic survey) (Kearey et al., 2002;
Telford et al., 1990; Reynolds, 1997).
Geologic zones of structural weaknesses (fracture and shear zones) can be delineated from
aeromagnetic survey data (Udensi and Osazuwa, 2004; Balogun, 2019). Such zones of
weakness have been known to be prospective hosts for a range of minerals. They are also
indicators for exploring epigenetic, stress-related mineralization in rocks (Paterson and
Reeves, 1985). They could as well serve as potential productive regions for ground water
abstraction (Oni et al., 2020), however, their occurrence in large numbers and diverse
orientations may render a place not suitable for construction purposes (Isaksson et al, 2017;
Momoh et al., 2008).
One of the most significant pieces of information needed for planning in pre-construction site
study, where geophysics has grown into a reliable, fast, cost-effective, and non-invasive tool,
is the knowledge of the orientation, distribution, and deep extent of linear geologic features
which include fractures and rock contacts (Isaksson et al, 2017). In most situations where pre-
construction site research studies were deemed superfluous, avoidable structural flaws were
discovered, and a large amount of money would have been squandered at this point (Mattsson
and Wahlgren, 2010).
As a result, major construction companies have begin to spend in pre-construction studies
that include geophysical mapping of linear geologic features, many of which could be
harmful to engineering constructions. The benefit of such pre-construction mapping is that
resources are conserved, land is effectively utilized, engineering buildings are guaranteed to
be durable, and human lives are not put at risk.
In line with this evolving good practice of pre-construction site investigation consisting of
various geophysical methods for mapping linear geologic structures that could pose threats to
civil engineering structures, the Study area was investigated for potential linear geologic
structures that could be inimical to engineering structures using the aeromagnetic data.
The Igabi Local government area of Kaduna state is considered to be an area of high
mineralization due to its tectonic formation. Aeromagnetic survey is best suitable because of
accessibility, the data was interpreted to delineate tectonic fault In order to indicate specific
zones of mineralization.
1.3 Study Area
Igabi Local Governmnet Area lies within latitudes 100 25′ 28″ N and 110 35′ 53″ N and Longitudes 70
21′ 49″ E and 70 50′ 00″ E (Figure 1.1). The area covers an area 9 approximately 3,727 square
kilometers and shares boundaries with Kaduna North, Kaduna South, Zaria, Kajuru, Kauru, Igabi and
Birnin-Gwari Local Government Councils. Turunku is the headquarters of Igabi LGA which was the
seat of power of the famous Queen Amina of Zazzau. The seat of Queen Amina and many historical
evidences are still located on Turunku Hill as well as tomb of Bakwa Turunku the father of Queen
Amina andher sister Zaria where Zazzau Emirate derived its name.
1.3.2 Relief, Climate and Vegetation
Geologically the study area is covered by sedimentary rock formation consisting of sand stones and
alluvial deposited especially around the site of the park. The area lies on a flat gently rolling land
with few hills. In and around the Park, there are important land marks which are of potential interest
to tourist. Some of the exquisite barren rocks make the landscape alluring and this create attractive
scenery. Some of these lands are Tsunin Buffa, Kashiga Mountains, Marmara Hills, etc. the beauty of
the under field land marks inhabited 13 by both animals and birds sharing the luxury of tropical
sunshine alongside meandering rivers makes the park attractive and beautiful. Most of the hills are
rounded granite domes and flat topped mensas or residual hills with reddish-brown sand stones are
characteristics relief features of the study area (Jatau et al., 2013). Generally the topography of the
study area is that of an undulating Plateau that forms part of the rich tourist attractions in areas like
Kufena in Zaria, Kagoro, Kwoi, Gwantu etc. These areas have protruding hard resistant granite rocks
that are so attractive for sightseeing. Variable elevations of land topography ranging from the lowest
height of 380m above sea level in river valleys and as high as 450m occurring mostly in upland areas
leads to the spring up of streams arounds the study area. The Ukinkina, Kuzomani and Kuwaimbana
rivers are some of the major streams that drain around the study area. Other streams include the
Dagara stream, Bogoma stream and Numashi stream. Kwiambana Rivers forms quite an extensive
length of the northern boundary of the study area, draining from west to east. The two streams
Kuzomani and Ukinkina both drains from the centre parts of the park. Kuzomani drains to the west
into the Kuwaimbana while the Ukinkina drains out from the park in the South West (Shehu, 2011).
Its main rivers are River Kaduna, River Gurara, River Kogum and River Kubani (Shehu, 2011).
The study area has an average annual rainfall of 1250mm. The rains occurs between months of April
– October when the South Westerly humid winds brings in rain. The dry seasons last between
November and March when the prevailing North Easterly winds Study Area 11 (Harmattan) brings
with it dusty, dry and cool air of the Sahara desert ushering in the dry season. The mean annual
rainfall in the study area ranged from 1000mm to 1500mm. The month of August-September
recording highest rains of 300mm (Yamusa, et al., 2015). The rainy season starts between 10th of
The temperature of the study area resembles that of the North Central Zone of the country.
Temperature ranges between 250 C – 350C during the dry season. The temperature may rise to
about 420C in March/April which is the hottest period. The coldest month is December/January.
During the harmattan popularly referred to as the West Africa Doctor, temperature sometimes
reaches freezing point (Record from the nearest meteorological station over a period of 15years)
1.3.2.3 Vegetation
The vegetation cover of the study area is uniform and monotonous typical of the Northern Guinea
Savanna Zone. Green (2007) further described the vegetation as containing three vegetation sub-
types viz: 1. Isoberlinia doka woodland grassland 2. Gallery or riparian vegetation 3. Shrub land
composed of trees – shrub association. Other areas of the forest are covered by Uapacamonotes
Afornosia association, Vtogoensis, M. Kekrstingi afzelia Africana, etc. In the low ground near
Numanshi stream which stretches into the park serving as a source of water for animals during the
dry season, 12 one finds a terminalia Pseudocederela Mitrogya association with; Terminallia
macrotera, Mitragyna inermis, philiostigma inermis. The presence of these vegetation types signifies
clay composition in the soil. The soil types of the study area corresponds fairly closely with the
vegetation belts of the guinea savanna which extends from east to west of Kaduna metropolis.
Excepts for highland areas, Kaduna metropolis has high temperature which promotes active
chemical and biological changes in the soil. There are few outcrops and some fossilized literate hills
in the study area. Soil types includes; sandy loams and a little clay in valleys (Oyedele, 2011). The
climate, soil and hydrology supports the cultivation of most of the study area staple crops and still
leaves ample scope for grazing, fresh water, fishing and forestry. Rich fertile alluvial soils are visible
around the stream banks. Estimate shows that 80% of the people in the local government area
engage directly or indirectly in agricultural activities (Oyedele, 2011). Oyedele (2011) opined that the
soils are yellowish to reddish, deep, well drained and fine or medium textured. Bleached layers of
ferruginous tropical soils associated with ferralitic soils are also found along the River Kaduna flood
plain.
1.3.3 Geology of the Study Area
The area understudy is underlain by Precambrain rocks of the Nigerian Basement Complex. The
weathering of the crystalline Basement Complex rocks under tropical condition is well known to
produce a sequence of unconsolidated material whose thickness and lateral extent vary extensively.
Groundwater localization within the Basement Complex occurs either in 14 the weathered mantle or
in the fracturing, fissuring and jointing systems of the bedrock (Ako and Olorunfemi, 1999; Olayinka
and Olorunfemi, 1992). These unconsolidated materials are known to reflect some dominant
hydrologic properties, and the highest groundwater yield in Basement Complex area are found in
areas of thick overburden overlying fractured zones and are characterized by relatively low
resistivity. The Basement Complex rocks in the areas are mostly migmatite, granite gneiss,
This study is aimed at delineating tectonic faults that could be linear geological structures
from aeromagnetic data acquired over Igabi and its environ, so as to make the information
(i) reduce the aeromagnetic data to the magnetic equator and separate the regional
(ii) perform basic structural enhancement operations which include Total Horizontal
(iii) identify the magnetic lineaments delineated within the study area.
(iv) present the magnetic lineament map of magnetic lineaments found within the
study area.
Figure 1.3: Geologic Map of the Study Area (digitised from Geological Map of
Nigeria 2004).
1.5 Literature Review
technique) has proved to become reliable as a means of rapidly checking the subsurface for
geologic structural discontinuity (Oyeniyi et al., 2014; Balogun, 2019; Oni et al., 2020). This
practice has continued to gain more popularity and patronage because it is effective and
reliable, easily applicable to vast potential field data and the mathematical procedure are not
The enhancement methods commonly applied include the Total Horizontal Derivative
(THD), Second Vertical Derivative (SVD), Analytical Signal (AS), Tilt Derivative (TDR),
Euler Deconvolution among others (Reid et al., 1990; Oyeniyi et al., 2014; Oladele and
Ayolabi, 2014). These methods rely on simple applications of calculus and Fourier
transformations to enhance edges of spatial field signals (Reid et al., 1990; Balogun, 2019).
Linear geologic structures are important surface/subsurface features with a lot of significance
(Abdullah et al., 2013). Being potential host for the accumulation of groundwater, they are
the major structure of interest in ground water prospecting (Ejepu et al., 2017). However, in
construction works, they are avoided as much as possible (Amadi et al., 2012).
One of the most significant pieces of information needed for planning in pre-construction site
studies where geophysics is gradually growing to become a reliable, fast, cost-effective, and
non-invasive tool, is the knowledge of orientation, distribution, and deep extent of linear
geologic features (Isaksson et al, 2017; Adewumi and Olorunfemi, 2005). In most situations
where pre-construction site research studies were deemed unnecessary, avoidable structural
flaws are often discovered and a large amount of money wasted (Mattsson and Wahlgren,
2010). To avoid this wastage and risk pose to lives, major construction companies have begin
Magnetic method has been successfully employed in the mapping of geological structures for
the purpose of construction and ground water prospecting both in sedimentary and basement
complex terrains. One of such studies where magnetic method was used is the work of
Oladejo et al., (2020), where high resolution magnetic data of the sedimentary terrain and
part of the basement complex of Southwestern Nigeria were processed and interpreted to
provide fault architecture of the area. It represents a good case study for the applicability of
magnetic method in mapping geologic structure both within the sedimentary terrain and
Another case study where aeromagnetic mapping was used in the delineation of lithological
structures is the work of Okpoli and Oludeyi (2019). The total magnetic intensity field data
was reduced to the magnetic equator, regional and residual fields were separated, before
suites of enhancement techniques which include upward and downward continuation, tilt
derivative, second vertical derivative, analytical signal, and Euler deconvolution were used to
High resolution aeromagnetic data of the region (Sheet 279 – Ilaro Sheet) where the study
area is located was procured from Nigerian Geological Survey Agency (NGSA) and the study
Near-surface (cultural) noises were first filtered using the upward continuation approach. The
data were subsequently reduced to the magnetic equator using the “Reduction to the
Magnetic Equator (RTE)” filter operator. This was done in order to reposition anomalies at
magnetic low latitude to appear as if there were actually acquired at zero magnetic inclination
The “Reduced to the Magnetic Equator” TMI data were decomposed into residual and
regional field components and the residual field was enhanced using the Total Horizontal
Derivative (THD), Euler Deconvolution and Tilt Derivative (TDR) enhancement operations
for proper delineation of the magnetic linear structures. The delineated lineaments were
Tomography (ERT) data were subsequently acquired within the study area and the result
The study will reveal the locations of plausible linear geologic structures and their density
Within the vicinity of a bar magnet, magnetic flux is developed which flows from one end of
the magnet to the other (Figure 2.1). This flux can be mapped from the directions assumed by
a small compass needle suspended within it. The points within the magnet where the flux
converges are known as the poles of the magnet (Telford et al., 1990; Kearey et al., 2002). A
freely-suspended bar magnet similarly aligns in the flux of the earth’s magnetic field. The
pole of the magnet which tends to point in the direction of the earth’s north pole is called the
north-seeking or positive pole, and this is balanced by a south seeking or negative pole of
identical strength at the opposite end of the magnet (Lowrie, 2007). The force, F, between
μ 0 m1 m2
F= 2 2.1
4 πμ R r
where μ0 and μ R are constants corresponding to the magnetic permeability of vacuum and the
The force is attractive if the poles are of different sign and repulsive if they are of like sign.
The magnetic field B due to a pole of strength m at a distance, r, from the pole is defined as
μ0 m
B= 2 2.2
4 π μR r
Figure 2.1: Magnetic Flux Surrounding a Bar Magnet (Telford et al., 1990)
Magnetic fields can be defined in terms of magnetic potentials in a similar manner to
gravitational fields. For a single pole of strength, m, the magnetic potential, V, at a distance, r,
μ0 m
V= 2.3
4 π μR r
The magnetic field component in any direction is then given by the partial derivative of the
potential in that direction. In the SI system of units, magnetic parameters are defined in terms
of the flow of electrical current. If a current is passed through a coil consisting of several
turns of wire, a magnetic flux flow through and around the coil annulus which arises from a
The magnitude of H is proportional to the number of turns in the coil and the strength of the
current, and inversely proportional to the number of turns in the coil and the strength of the
current, and inversely proportional to the length of the wire, so that H is expressed in Am−1 .
The density of the magnetic flux, measured over an area perpendicular to the direction of
flow, is known as the magnetic induction or magnetic field B of the coil. B is proportional to
H and the constant of proportionality μ is known as the magnetic permeability. Lenz’s law of
induction relates the rate of change of magnetic flux in a circuit to the voltage developed
within it, so that B is expressed in Vsm−2 also known as Weber (Wb) (Telford et al., 1990;
The unit Wbm−2 also designated as the tesla (T) is the unit of measurement in magnetic
survey. Permeability is consequently expressed in Wb A−1 m−1 or Henry per meter, Hm−1. The
c.g.s unit of magnetic field strength is the gauss (G), numerical equivalent to 10−4 T. The
Tesla is too large a unit in which to express the small magnetic anomalies caused by rocks,
and a subunit, the nanotesla (nT), is employed (1 nT = 10−9 T ¿ . The c.g.s system employs the
numerically equivalent gamma (γ ), equal to 10−5 G. Common magnets exhibit a pair of poles
and are therefore referred to as dipoles. The magnetic moment M of a dipole of strength, m, at
M = ml 2.4
The magnetic moment of a current-carrying coil is proportional to the number of turns in the
coil, its cross-sectional area and the magnitude of the current, so that magnetic moment is
magnetization in the direction of the field which is lost when the material is removed from
the field (Telford et al., 1990; Reynold, 1997; Lowrie, 2007). This phenomenon is referred to
as induced magnetization or magnet polarization, and results from the alignment the material
has magnetic poles distributed over its surface which correspond to the end of the dipoles.
According to Telford et al. (1990) and Lowrie (2007), the intensity of induced magnetization
M
J i= 2.5
LA
J i=κ H 2.6
Where κ is the magnetic susceptibility of the material, since J I and H are both measured in A
−1
m . Susceptibility is dimensionless in the SI system. In the c.g.s system susceptibility is
similarly dimensionless.
The magnetic induction, B, of a magnetic material such as a rock with positive magnetic
susceptibility is given by the total field within the rock material which includes:
(a) the strength of the earth’s magnetic field, Hi, that the rock material is exposed to
(b) the strength of the magnetic field in the rock materials itself, H, arising from its being
According to Blakely (1996), the field Hi is a product of the intensity of magnetization J iand
a solid angle 4 π
H i=4 π J i 2.7
B=H + 4 πκ H 2.9
The ratio of the magnetic induction B to the magnetizing field H is referred to as the
magnetic permeability.
B = μH 2.12
“ μ is the magnetic permeability of the medium between the magnetizing field H and the rock
material.
2.1.2 The Geomagnetic Field
The Earth’s magnetic field is believed to originate in the core. Although the mechanism is
dynamo theory) (Lowrie, 2007). It is currently believed that the geomagnetic field is
generated by shearing and twisting fluid motions within the liquid, electrically conducting
outer core, the result being a dynamo, which maintains and regenerate the earth’s magnetic
field (Reynold, 1997). It is believed that only 1% of the earth’s magnetic field is due to
sources outside the earth (atmosphere). The geomagnetic field is usually described by a bar
magnet at the core, with the lines of magnetic force looping from the south magnetic pole to
In practice, it is a little more complex with some multi-polar patterns, but to a first
approximation the dipolar model is convenient. While the N-S magnetic polar axis is
currently inclined at an angle of about 100o to the earth’s axis of rotation, it is thought that the
magnetic and geographic poles would coincide when averaged over a long period of time.
The earth’s magnetic field is a vector quantity; that is, it has both magnitude and direction.
The magnitude of this field is relatively weak, with maximum intensity near the magnetic
poles. The geomagnetic elements (Figure 2.2) are taken to be components parallel to the
geographic north and east directions. From a geological point of view, there are two
(b) The angle of inclination (I) which is the angle of dip of a magnetized needle relative
to the horizontal. This angle varies from 90 o (vertical) at the magnetic poles to 0 o
And consequently;
2 2 2 2
X =Y + Z =F 2.14
Majority of magnetic surveys are conducted in the air, with the sensor either towed in a 'bird'
to shield it from the aircraft's magnetic field or placed in a 'stinger' in the aircraft's tail, in
which case the inboard coil installations compensate for the aircraft's magnetic field. Aircraft
usually employed for airborne magnetic surveys can either be a fixed-wing or rotary-wing
(helicopter) aircraft. The rotary-wing (helicopter) is most suitable for rough terrain. Fixed-
wing aircraft are less expensive, cover more ground faster, and produce higher-quality data in
Two basic sensor configurations are usually employed in airborne magnetic survey and are
explained below;
1. The towed-bird Configuration
The sensor is placed far enough from the airplane in this configuration so that the magnetic
field generated by the aircraft is minimal at the sensor. The towed-bird system produces
excellent data at the low speeds typical of helicopter surveys, and it is the most popular
(bird) motion and tow cable vibrations at the comparatively high speeds of fixed-wing
aircraft (about 150 mi/h, or 240 km/h). The towed-bird installations are sometimes employed
with fixed-wing aircraft. Figure 2.3 is a typical line plan for an airborne magnetic survey.
The sensor is mounted on a stiff rod at the aircraft's back. The sensor is sometimes installed
2.1.4.2 Instrumentation
The following are the essential instruments required for airborne magnetic survey
1. Stinger or towed-bird
and GPS.
4. Analog Recorder: Used to record specific parameters, such as magnetics and altimeter
system that provides actual, visible track information in addition to the Doppler
navigation system.
6. Magnetic Compensation Device (fixed-wing only): used for accounting for the aircraft's
induced (both electrical and platform motion) and permanent magnetic fields.
Auxiliary equipment used in airborne magnetic survey include: inertial frame of reference
support equipment includes a magnetometer and recording device at the base station and a
1. Aeromagnetic surveying is quick and inexpensive, costing about 40% less per line
2. Large areas can be quickly surveyed without incurring the expense of deploying a field
Position fixing used to be the most challenging difficulty in airborne surveys, but the
Either aeromagnetic, ship borne or ground magnetic survey, magnetic readings are usually
affected by stray fields from iron and steel (cultural noise) most of which we cannot fully
steer clear off during magnetic data acquisition. To eliminate this kind of noise, the noise is
filtered from the aeromagnetic data by upward continuing the data to a moderate height of 80
m so that the cultural noises (which are known to be dominated by high wavenumber
spectrums) can decay rapidly as we move higher away from the source.
CHAPTER THREE
The elevation data adopted for the study is the high-resolution, Shuttle Radar Topographic
Mission (SRTM) global elevation data (Digital Terrain Elevation Data – DTED®) having a
spatial resolution of 30 m × 30 m and a near world-wide coverage. The data was acquired by
radar interferometry and made use of the Earth’s Gravity Model 1996 vertical datum. The
unit of height measurement is meter(m) and the C-band wavelength is 5.6 cm.
The magnetic data used for the study is the high-resolution 100 m × 100 m grid aeromagnetic
(Total Magnetic Field Intensity (TMI)) data, having a mean terrain clearance of 80 m and
recording interval of 0.1 s, acquired by the Nigerian Geological Survey Agency between
2003 and 2010. For this study, the data was reduced to the magnetic equator to remove
asymmetries associated with low magnetic latitude anomalies. The regional field was
removed and the resulting residual field was processed and analyzed.
The aeromagnetic data were processed using the Geosoft’s Oasis Montaj TM software version
6.4.2 (HJ). The suites of basic data enhancement operations which include the least square
This technique was first suggested for use for regional-residual field separation in potential
field data processing by DeLury (1950) and was used by Oldham and Sutherland (1955). It
Let the observed magnetic field T (x, y) be composed of a regional field T R (x, y) and a
Then
Let T Lbe a random variable with zero mean and constant variable σ 2. Let the observed
T (x, y) = a 00 f 0 ( x ) f 0 ( y ) +a 10 f 1 ( x ) f 0 ( y )
a 01 f 0 ( x ) f 1 ( y ) +…
a pq f p ( x ) f q ( y ) 3.2
where constant a rsare determined by the least square techniques using the equation
a rs =¿ ¿ 3.3
accentuate edges of structures. The principle is based on the fact that the gradient of a
field function such as the gravity and magnetic fields, where the derivative is peak is usually
THD = √ ¿ ¿ 3.4
δT δT
Where is the derivative of the magnetic field in the “x: direction, and is the derivative
δX δy
field data processing. It is widely used because it requires only little prior information about
anomaly source geometry and requires no information about anomaly source vector in the
case of magnetic anomalies (Reid et al., 1990; 1980; Balogun et al., 2016; Balogun, 2019).
One major important advantage of Euler deconvolution is that it gives the estimate of the
depth of magnetic sources in addition to defining their edges. The Euler Deconvolution
operation is based on solving Euler’s homogeneity equation given in equations 3.5 and 3.6 as
δT δT δT
(x- x 0 ¿ + (y- y 0 ¿ + (z- z 0 ¿ = η ( β−T ) 3.5
δχ δy δz
δT δT δT δT δT δT
x +y +z + ηT = x 0 + y0 + z0 + ηβ 3.6
δχ δy δz δχ δy δz
where β is the regional value of the potential field and x 0 , y 0 , z o defines the position of the
source, which produce the total potential field T measured at (x, y, z), η is the structural
index.
3.2.1.4 Tilt Derivative
A generalized definition for the local phase of a potential field function is known as the tilt
derivative (or angle). Due to the efficient Automatic Gain Control (AGC) enforced by the arc
tangent function, which restricts the tilt angle to within the range -90 to +90 independent of
the amplitude or wavelength of the magnetic field, mapping the magnetic tilt angle offers the
Tilt derivative is used to estimate the location and depth of magnetic sources. It is based on a
model of a buried 2D vertical contact and it is a very straightforward way to estimate the
location and strike of geological contact/faults, as well as the depth to basement of magnetic
anomalies (Fairhead et al., 2009). Depths to a magnetic source are estimated by taking the
average of the distance between +45 o and -45o contour lines (Verduzco et al., 2004; Salem
VDR
TDR = tan-1
THD
3.7
where VDR is the first vertical derivative and THD is the total horizontal derivative of the
TMI respectively.
REFERENCES
Arora, K., Cazenave, A., Engdahl, E. R., Kind, R., Manglik, A., Roy, S., and Uyeda, S.
(2011). Encyclopedia of solid earth geophysics. Springer Science & Business Media.
complex of Ilorin area from aeromagnetic data. NRIAG Journal of Astronomy and
university press.
Ejepu, J. S., Olasehinde, P., Okhimamhe, A. A., and Okunlola, I. (2017). Investigation of
Fairbridge, R.W., Askew, A.J., Herschy, R.W., Herschy, R.W., Herschy, R.W., Hordon,
R.M.,
Muller, R.A., Grymes, J.M., Fairbridge, R.W., Herschy, R.W. and Brassington, F.C.,
Resources, pp.699-706.
Fairhead, J. D., Salem, A., and Williams, S. E. (2009). PS Tilt-Depth: A Simple Depth-
Estimation Method Using First Order Magnetic Derivatives.
magnetic anomalies using tilt derivatives: Part II — Application to real data," SEG
Google Maps (2020), Mountain Top University Permanent Site, Makogi-Oba, Ogun State,
Gore, D. and Davies, P. (2013). Background paper on produced water and solids in relation to
Grant, F.S (1957). A problem in the analysis of geophysical data. Geophysics 22, 309-344.
Griffiths, D.H. and King, R.F., 2013. Applied geophysics for geologists and engineers: the
Kearey, P., Brooks, M., and Hill, I. (2002). An introduction to geophysical exploration. John
Southwestern
suspected fault zones around Iwaraja-Ijebu Ijesa area of Osun state, Nigeria. Ife Journal
Oladejo, O. P., Adagunodo, T. A., Sunmonu, L. A., Adabanija, M. A., Enemuwe, C. A., and
Oladele, S., and Ayolabi, E. A. (2014). Geopotential imaging of the benin basin for
Paterson, N. R., and Reeves, C. V. (1985). Applications of gravity and magnetic surveys: The
Conditions. ISRM India Journal-Half Yearly Technical Journal of Indian National Group
Reid, A. B., Allsop, J. M., Granser, H., Millett, A. T., and Somerton, I. W. (1990). Magnetic
John
Wiley &Sons.
Ribes, D. and Finholt, T. A. (2009). The long now of infrastructure: Articulating tensions in
development. Paul Edwards, Geoffrey C. Bowker, Steven Jackson, and Robin Williams
Robinson, D. A., Binley, A., Crook, N., Day‐Lewis, F. D., Ferré, T. P. A., Grauch, V. J. S.,
and
near‐surface geophysics: A vision for, and review of, electrical and magnetic
Salem, A., Williams, S., Fairhead, J. D., Ravat, D. and Smith, R. (2007). Tilt-depth method:
simple depth estimation method using first-order magnetic derivatives. The leading
Telford, W. M., Telford, W. M., Geldart, L. P., and Sheriff, R. E. (1990). Applied
geophysics.
Udensi, E. E., and Osazuwa, I. B. (2004). Spectral determination of depths to magnetic rocks
Schewe, O., Tao, H., Gemmer, M., Bai, Y., Su, B. and Mao, W. (2011), ―Trends of stream
flow in the Tarim River Basin during the past 50 years: Human impact or climate change?‖,
Journal of Hydrology, Vol. 400, pp.1-9. Shehu, S. (2011). The Growth and Development of
Nigeria Since the 20th Century. Kaduna State University, Kaduna-Nigeria. Environmental
Reuben, F. and Barau J. (2012), ―Assessing climate change impacts and adaptation
strategies for smallholder agricultural systems in Uganda‖, African Crop Science Journal,
Vol. 20, No. 2, pp. 303 – 316. Selvaraju, H., Benessaiah, K., Barrera, J. F., Cruz-Bello, G. M.
and Morales, H. (2012), ―Livelihoods and landscapes at the threshold of change: Disaster
and resilience in a Chiapas coffee community‖, Regional. Environmental Change, Vol. 22,
Kaduna South and environs using Wenner Offset Method of Electrical Resistivity Sounding.
Journal of Earth Sciences and Geotechnical Engineering, 3(1): 41-54 Kaduna State Ministry
of Land and Geoinformatics, 2018. Kandlinkar, M & Risbey, J. (2010). Agricultural impacts
of climate change: If adaptation is the answer, what is the question? Climatic change, 45,
529-539.