Application of Magnetic Survey in the Investigation of Iron Ore Deposits and Shear Zone Delineation: Case Study of Mutomo-Ikutha Area, SE Kenya ()
1. Introduction
The use of magnetic data is a geophysical exploration technique based on mapping the subsurface distribution of magnetic minerals (oxides or sulphides). This technique can be a valid approach to mineral exploration either directly for mineral deposits containing magnetic minerals (e.g. magnetite, hematite, pyrhotite) or indirectly for the ones spatially associated with, for example, faults. Advances in spatial controls using global positioning systems (GPS) have not only reduced the survey costs but also increased data accuracy.
Total magnetization of rock is a vector sum of the two components: 1) induced magnetization (proportional in magnitude and generally parallel to the earth’s ambient field); and 2) remnant magnetization (which has direction and intensity dependent on the origin and geological history of the rock). The intensity of induced magnetization is related to the ambient field through the magnetic susceptibility of the rock considered (Hildenbrand et al. [1] ). This magnetic susceptibility constant is directly proportional to the modal mineral composition and chemistry of the rock.
Rocks have widely varying magnetic properties (Carmichael [2] ; Clark [3] ). At the regional scale, magnetic highs are commonly associated with major igneous provinces in crystalline basement. By contrast, magnetic lows often occur in areas dominated by thick sedimentary basins (Hildenbrand et al. [1] ) or where, for example, igneous rocks are altered and magnetite is replaced by hematite when hydrothermal fluids have migrated along faults. Telford et al. [4] emphasized that magnetic response observed in magnetic anomaly maps was influenced by the geometry and depth of the magnetic bodies, their azimuth with regard to magnetic north, and the inclination of the magnetic field at the latitude of the survey. In the neighbourhood of mineral deposits, mineralogical variety and thermal history related to mineral deposit can subsequently affect magnetic properties. The degree of oxidation, alteration, weathering, and metamorphism also affects the magnetic properties of the rock (Hildenbrand et al. [1] ). All these combined factors concur to generate complex sections of magnetic properties for ground magnetic survey in Mutomo-Ikutha area.
The main objective of this research was to investigate and determine the surface and subsurface distribution of the iron ore mineralization within the shear zone so as to establish its economic potential.
The study area is located between longitudes 38˚4'E to 38˚20'E and latitudes 1˚48'S to 2˚8'S in South Kitui within the Kitui County in the republic of Kenya (Figure 1).
Figure 1. The map showing the area of the project. Major towns within the area include Mutomo and Ikutha.
2. Geology of the Study Area
Mutomo-Ikutha area occurs within the Neoproterozoic Mozambique Belt. This belt is a major N-S trending metarmophic and lithotectonic domain that extends along the Eastern Coast of Africa, as well as in Saudi Arabia, Madagascar, India and Sri Lanka. It was first defined by Holmes [5] , mostly on metamorphic and structural criteria. The Mozambique is a Neoproterozoic Belt (900 to 550) which is affected by the Pan African tectonothermal event about 650 Ma, as suggested by (Kennedy [6] ; Kazmin [7] ; Kröner [8] ; Cartier [9] ) who considered it as a major phase of division of the African continent linked to the development of mobile belts. The Pan African event encompasses a wide time span from 650 to 500 million years. Mozambique Belt is considered as an orogen resulting from the oblique collision between two Gondwana fragments (Muhongo [10] ; Burke and Derwey [11] ; Gass [12] ; Key et al. [13] ; Kröner [8] ). Stern [14] proposed its name to be changed to East African Orogen. The geology, structures and tectonic the history of the Mozambique Belt in Kenya have been reviewed by Waswa et al. [15] . Other reviewers of this belt include Mathu and Tole [16] ; Mathu [17] ; Mathu et al. [18] ; Mathu et al. [19] ; Gaciri et al. [20] ; Biyajima et al. [21] ; Suwa et al. [22] ; Opiyo-Aketch and Nyambok [23] ; Hackman et al. [24] ; Mosley [25] ; Nyamai [26] -[32] ; Saggerson [33] ; Walsh [34] .
The Mutomo-Ikutha area is a medium- to high-grade gneiss terrane which comprises predominantly amphibolite-grade metasediments (pelitic, semi-pelitic, quartzo-felspathic, calc-silicate rocks, quartzite and marble migmatizedparagneisses), metavolcanics and intrusives (see Figure 2).
Figure 2. Geological map of Mutomo-Ikutha area (after Waswa et al., 2015).
The following mineral assemblages are represented by the lithologies within Mutomo-Ikutha;
・ Hornblende + Quartz + Biotite + Magnetite + Apatite ± Plagioclase, Microcline
・ Quartz + Plagioclase ± hornblende, biotite, muscovite, Magnetite and apatite
・ Hornblende + Quartz + Plagioclase + Biotite ± Epidote, apatite, Iron ore, Garnet, Zircon
・ Hornblende + Hypersthene + Quartz + Plagioclase + Diopside + Apatite + Biotite + Magnetite
・ Hornblende + Quartz + Plagioclase + Myrmekite + Biotite + Diopside ± Apatite, zircon, sericite, iron ore)
・ Plagioclase + Hornblende ± Biotite, Diopside, Hypersthene, Calcite, Sphene, Apatite, Magnetite.
3. Methodology
Geophysical investigation utilizes the measurements of physical quantities that are made at or above the ground surface and are used to investigate concealed geology. Proton magnetometer was used in this investigation. The subsurface structures in a part of the southern part of Ikutha area were estimated by the interpretation of ground magnetic data. Magnetic survey was used to delineate magnetic lineaments, analyze its relationship to tectonic fabric and estimate the depth of perturbing body source.
3.1. Magnetic Data Acquisition
Total intensity magnetic field measurements were carried out with MP2 SCINTREX proton precession magnetometer made by Scintrex Company of Canada with one nanoTesla (nT) sensitivity (Figure 3). The same magnetometer was used for base station recordings to apply the diurnal variation correction as well as taking field measurements. Magnetic surveys are based on the premise that a target is limited in space and has a different physical property (e.g. magnetic susceptibility), from the surrounding formation. Unlike gravity surveying, however, the variation in magnetic susceptibility for various rock types is orders of magnitude greater than the variation in density for the same rock types. Thus, even knowing the types of rocks in a specific area does not provide sufficient information to constrain susceptibilities. Like density contrast, variations in susceptibility tradeoff strongly with other model parameters. Therefore, if susceptibility, or other model parameters, cannot be constrained from different observations, it is difficult to make quantitative estimates of the geologic structure based on magnetic observations alone.
In a particular survey, additional constraints were applied that allowed the magnetic observations to be used in a quantitative fashion. This information was derived from other separate data sets, for example, formation layering from geological mapping. The general trend of the rock units was in the northwest direction. A magnetic survey was designed to estimate the spatial extent of the structure and its susceptibility by in-lab forward
Figure 3. Ground magnetic survey in Mutomo-Ikutha area using Proton precession magnetometer.
modeling. In planning the magnetic survey, the noise was predicted from sources not of interest in the survey, estimated the standard deviation of the random (operator and instrument) noise, calculated the shape of the signal (the theoretical anomaly produced by the assumed source), then decided whether the signal generated by the target of interest was above the noise level that allowed a meaningful interpretation to be conducted.
The Magnetic survey orientation in Mutomo-Ikutha area was in the northeast-southwest profile across the general strike of the rock units in the area. Survey was carried out at an interval of 20 m along each profile and 250 m from one profile to the next. MP2 SINTREX proton precession magnetometer was used to carry out the survey (Figure 4). Readings taken included position coordinates altitude, time and magnetic variations in nanoTesla (nT).
Magnetic Data Processing
1) Diurnal variation correction
Variation of earth’s magnetic field with time, due to the rotation of the earth and with respect to the solar wind, which may last for several hours to one day, is called diurnal variation (Riddihough [35] ). In order to correct for drift or diurnal effect in the magnetic readings, a base station within the area of the survey and free from magnetic interference was selected. Repeated readings were taken approximately after every one hour of the magnetic measurement for the diurnal correction. After the data collection, the diurnal effect was calculated and the magnetic data were filtered by using Equations (1) and (2):
(1)
where Drift is the drift of measurement in a loop; B base,f and B base,i are the final and initial total magnetic field at the base station; tbase,f and tbase,i are the final and initial time at the base station.
The drift corrected field at any station in the loop for a given day was calculated from:
(2)
where Bdrift,n is corrected magnetic field at the nth station in the loop; Bn is the total magnetic field at the nth station; tn is the measuring time at the nth station; tbase is the measuring time at the first station of the loop.
2) Removal of geomagnetic field
The regional magnetic field, often referred to as geomagnetic field, is needed to be subtracted from the acquired total magnetic field to obtain the magnetic field anomaly caused by the local source (Parasnis [36] ). After subtracting the diurnal effect from the original magnetic data observed, the geomagnetic field was calculated
Figure 4. The map showing the location where ground magnetic survey was taken.
by using the mathematical model of earth magnetic field called International Geomagnetic Reference Field (IGRF) model 2000-2015 in potent software, a product of Geosoft Oasis Montaj software. This model is calculated based on the dates, elevation and geographical locations (latitudes and longitudes) of the observed magnetic data with the generated average geomagnetic field of 33,000 nT, inclination of −25.3˚ and declination of 0.4˚. The IGRF values were subtracted from the observed magnetic values for each station to determine the residual magnetic field due to anomalous contribution from local magnetic sources in the area.
3.2. Geochemical Investigations
Geochemical investigations were carried out to establish the chemical composition of the geological material that was anomalous. Geochemical survey involved fieldwork and the samples collected were analyzed using a variety of analytical techniques. Amongst the instrumental methods used in analysis included atomic absorption spectrometry (AAS) and X-ray fluorescence (XRF). Detailed identification of individual minerals was undertaken by using a polarizing microscope.
4. Results and Discussion
The results obtained were used to delineate the iron mineralized zone. Both 3D Euler deconvolution and 2D Euler deconvolution were used to synthesize the results. 3D Euler deconvolution utilizes XYZ (X = length along X axis, Y = length along Y axis, Z = length along the Z axis) values to provide the subsurface geophysical information. 2D data interpretation gave information along designated profiles and the depth. The obtained values were corrected for the normal gradient of the earth’s magnetic field (IGRF); the corrected magnetic values were plotted using Euler 2D and Oasis Montaj.
4.1. Delineation of the Shear Zone Using Euler Convolution
The high anomalous values plot on the western part of the area of survey (Figures 5-7). These values display a
Figure 5. The figure showing the trend of Iron ore mineralization in Mutomo-Ikutha area. The boundary of the iron rich rocks is marked by the brown colour.
Figure 6. The contour map showing the trend of Iron ore mineralization in Mutomo-Ikutha area. The boundary of the iron rich rocks is marked by the tightly spaced contours on the western part of the research area.
Figure 7. The 3D figure showing the trend of Iron ore mineralization in Mutomo-Ikutha area in 3D. The boundary of the iron rich rocks is marked by the red colour.
linear trend which is in the north-south direction. This trend shows that the geological body/structure having anomalous values occurs closer to the surface trending in the north-south direction. The high anomalous values range from 550 nT - 1000 nT.
4.2. 2D Euler Data Analysis and Interpretation
Interpretation for Profile T1-T16
High magnetic anomaly is observed from 0 - 160 meters below the surface along all the profiles (Figures 8-11), indicating the presence of a geological body which is closer to the surface. The anomalous values indicate mineralization which varies along each profile. Geological material analyzed is found to be magnetite hosted in hornblende gneiss (Table 1). Iron ore occurs in form of veins pod form, dike, sill and disseminated form. The iron ore which is magnetic is hosted in non-magnetic body (hornblende gneiss). Discontinuities of the anomalies indicate that shearing and faulting has occurred along western part of all the profiles.
4.3. Field Proving with Laboratory Investigations
Extensive geological, field proving and geochemical investigations were carried out for the areas with anomalous values and delineated in the studied gneisses by the magnetic survey. During the field works, the geochemical investigations revealed that the Fe2O3 was relatively high in the sheared gneisses. The mineralized zones were found to be mainly connected or hosted by the hornblende gneisses that followed the direction of the strong and nearly vertical shear zone striking N320W-S140E in the western part of Ikutha town.
5. Conclusion
The results indicate that iron mineralization in Mutomo-Ikutha has been controlled by tectonics, especially in the shear zone on the western part of the area. Chemical analysis shows that this deposit is of magnetite variety, having between 35% - 96% Fe2O3, 2% - 19% P2O5 and 0.1% - 5% TiO2 and resembling the Apatite-Iron ores of Kiruna type. Magnetic survey has shown that Iron ore in Mutomo-Ikutha occurs from the surface to 160 meters below the surface in form of veins, dikes, sills, and pod form extending more than 10 kilometers from north- western part of the area towards the south. The iron ore deposit is structurally controlled by shear zone and hosted in
Figure 8. (a) Schematic section of profile T3-T′3 through a sheared and faulted showing the underlying mineralization; (b) structural index pseudo-section showing faulting and shearing.
Figure 9. (a) Schematic section of profile T6-T′6 through a sheared and faulted showing the underlying mineralization; (b) structural index pseudo-section showing faulting and shearing.
Figure 10. (a) Schematic section of profile T7-T′7 through a sheared and faulted showing the underlying mineralization; (b) Structural index pseudo-section showing faulting and shearing.
Figure 11. (a) Schematic section of profile T10-T′10 through a sheared and faulted showing the underlying mineralization; (b) Structural index pseudo-section showing faulting and shearing.
Table 1. Chemical analysis of the samples that proved to be Magnetite.
hornblende gneiss. Shearing along the western part of the area, especially within Tiva, gneisses creates room for iron ore deposition. Iron vein trending in the North-South direction on the western part of the area occurs within the shear zone.