Characterization of lithological zones of the Isanlu sheet 225, North Central Nigeria, using aerogeophysical datasets

Gridded aeromagnetic and aeroradiometric datasets of Isanlu sheet 225 were processed and interpreted in order to map the geological structures and characteristic lithology, as well as potential mineral ore zones in the area. Recorded aeromagnetic data were reduced to the magnetic equator before filtering operations which enhanced the magnetic anomalies. Analytic signal map was created, and subsequently a model was generated along a selected profile for better interpretation of the subsurface materials. On the aeroradiometric data were performed both enhancement and statistical operations. The lithological zones were identified. They are metasedimentary rocks, migmatites and igneous rocks occupying varying proportions within the entire area. Some potential radioactive or precious metal ore zones were also localized. A major structure interpreted in the study area trends in the SE–NW direction. Other surface structures were also delineated and general orientation was evaluated. This information may serve as a valuable database for further prospective research and mineral exploration in the region.


Introduction
Conventional regional geological mapping and ground geophysical survey on a regional scale may be a Herculean task because of cost of execution, time of completion and inaccessibility at some instances. Therefore, appropriate knowledge and understanding of the general geology of an interesting area may be limited. However, the use of aerogeophysical datasets often enhances the success of the conventional geological mapping, especially for the regions of very large and/or mountainous, as well as allows quick regional reconnaissance (Boadi et al. 2013).
The magnetic method is capable especially of detecting and mapping of igneous extrusive rocks or locating the boundaries of the formations and determining the depth of crystalline basement in sedimentary basins. The measured parameter of the aeromagnetic surveys is total magnetic intensity (TMI) describing influences of stronger/weaker magnetic rocks on the geomagnetic field. The geological structures (lithology and tectonic) and their parameters (size, shape, depth to the surface, slope) as well as the magnetic properties of the formations can be investigated.
In radiometric surveys, gamma rays from rocks and soils are detected. Using the measured data, one can identify and locate the geological formations as well as occurrence of metal ore deposits. The terrestrial gamma radiation almost comes from the natural isotopes belonging to the decay series, i.e., uranium (U), thorium (Th) and potassium isotope ( 40 K).
Nowadays, airborne magnetic and radiometric surveys have been used extensively in the mineral exploration industry predominantly for the delineation of metalliferous deposits (Airo and Loukola-Ruskeeniemi 2004;Adepelumi and Falade 2017). Its application ranges from geological mapping, rock characterization and mineral exploration (Telford et al. 1990). Recently, the accuracy and resolution of these techniques have substantially increased, so that they can be used to provide useful enhanced information on lithology and geological structures. Additionally, advances in data processing and image analysis have improved the resolution of geophysical datasets to discover bodies causing very subtle variations in geophysical responses (Armstrong and Rodeghiero 2006).

3
Aero-survey measurements are usually made from lowflying airplanes along closely spaced, parallel flight profiles. The profiles are usually perpendicular to the main geological trending direction. These huge volumes of measurement data are then processed and presented as a digital aero-data map.
Potential gold deposits often occur in the hydrothermal altered zones rich in potassium (Ohioma et al. 2017a,b). The radioactive and heavy metals ores often occur in skarn or suture zones, accompanying with faults characterized by different petrophysical parameters including magnetic susceptibility. Based on the mentioned knowledge, an attempt has been made to interpret airborne radiometric and magnetic datasets to characterize the lithological zones and delineation of the potential metal ore places in the Isanlu sheet 225 North Nigeria.

Geological setting of Nigeria
The landmass of Nigeria is situated within the Pan-African mobile belt. The belt is sandwiched between two super-crustal terrains, also referred to as African and Congo Craton (Burke and Dewey 1972). The process of the geology formation of the country begins in the Archean and Proterozoic eons in the Precambrian (Obaje 2009).
The geology of Nigeria is dominated by crystalline and sedimentary rocks both occurring approximately in equal proportions (Woakes et al. 1987). Essentially, the regional geology of Nigeria is made up of three major litho-petrological components: the Precambrian Basement Complex, Younger Granites and sedimentary basins. The Precambrian Basement Complex is made up of the migmatite-gneiss complex, the schist belts and the older granite. They are mostly found in the north-central, southwestern and the far eastern flank of Nigeria (Fig. 1). The Precambrian Basement rocks of South West Nigeria found in the Dahomeyan (Benin) Basin, consist of migmatites, banded and granite gneisses, with low-grade metasedimentary and metavolcanic schists, intruded by Pan-African-age granites and charnockites (Oyawoye 1972;Fitches et al. 1985). Within the migmatites and gneissic metasediments, the intrusive pegmatite veins and dykes are often observed (Oluyide et al. 1998). Older granites, granodiorites and syenites, with dolerite dykes also form part of the Precambrian Basement of SW Nigeria. The petrographic characteristic of the older granite is typical leucocratic with gray to greenish tint. The overall texture of the rock as seen in most outcrops is predominantly hypidiomorphic granoblastic, and its dominant mineral components are quartz, plagioclase, microcline and biotite.
The Younger Granites comprise several Jurassic magmatic ring complexes, locating around the Jos Plateau and other parts of North Central Nigeria. The acid nature is the most striking petrographic feature of the Younger Granite. Over 95% of the rocks can be classified as rhyolites, quartz  (Adelana et al. 2008) syenites or granites with basic rocks forming the remaining 5%.
The sedimentary sequence includes the Late Jurassic to Early Cretaceous Illo and Gundumi Formations, the Maastrichtian Rima Group, the Late Paleocene Sokoto Group and the Eocene-Miocene Gwandu Formation. Examples of rocks of these formations are shale with intercalations of limestone, siltstone and sandstone. These were deposited during a series of overlapping marine transgressions. The sedimentary basins occupy the central Y-shaped area in the country; they underlie all of southern and are also found in northwest in Sokoto, and in Borno states and the entire Lake Chad in the northeast. The sedimentary basins, containing sediments of Cretaceous to Paleogene, comprise the Niger Delta, the Anambra Basin, the Lower, Middle and Upper Benue Trough, the Chad Basin, the Sokoto Basin, the Mid-Niger (Bida-Nupe) Basin and the Dahomey Basin. The basal formation is made of sandstone. This is overlain by Awgu shale and Nkporo shale. The topmost formation is made of fine to medium sandstone. Solid mineral deposits of economic significance that include gold, iron ore, cassiterite, columbite, wolframite, barites, lead-zinc, pyrochlore, monazite, marble, coal, limestone, clays, etc., occur in the different geologic units of Nigeria (Basement, Younger Granites and sedimentary basins), and indeed each of the 36 federating states and the Federal Capital Territory has a fair share of the solid mineral inventory of the nation (Adelana et al. 2008).

Geology of the study area
The Isanlu sheet 225 is located in North Central Nigeria between latitude 114,100-169,600 m and longitude 885,600-941,300 m in coordinate system WGS84/UTM zone 32 N (Fig. 2). Major cities in this area include Isanlu, Egbe and Moro. The study area falls within the Precambrian Basement complex with rock classified as the schist belts. The schist belts are predominantly composed of metasedimentary and metavolcanic rocks trending N-S which are best developed in the western part of Nigeria (Obaje 2009).
The relief of the area is moderate to high with steep to gentle rolling terrain and flat plains. Massifs and round topped hills are common with the highest peak being 2175 m above sea level (a.s.l) in the southeastern parts of the study area. These parts are predominantly undulating, while in the north the area is relatively low with average height of 750 m a.s.l.

Methods applied in the study
The aerogeophysical surveys were performed between 2004 and 2009 by Nigerian Geological Survey Agency (NGSA) using a light aircraft with fixed wings Cessna Titan 404 (C-FYAU). The aircraft had a gamma ray spectrometer mounted on its wings and a magnetometer placed in a cage towed in a housing known as a "bird" behind the aircraft. In the study area, the total measured line amounted to 1,104,174 km with the used flight parameters summarized in Table 1.

Aeromagnetic data
The total magnetic intensity (TMI) was acquired using an optically pumped magnetometer Scintrex Cesium SC-2™. The magnetic data were sampled at a rate of 0.1 s. The data were made available in the form of grid on a scale of 1:100,000. The geomagnetic global field was removed using The acquired aero-data were processed using Oasis montaj computer software (version 8.0) and Surfer version 10 developed by Geosoft and Golden software incorporation, respectively.
In the area close to the magnetic equator, the magnetization (M) of the body is consistent with the horizontal component of the Earth's magnetic field (H). Therefore, the TMI anomaly of a single magnetic body surrounded by nonmagnetic medium has the same value and direction as the horizontal component (ΔH) and is visible on maps as negative magnetic anomaly.
The measured magnetic intensities are presented as a map (Fig. 3a). Reduction to the magnetic equator (RTE) was done using the IGRF model and taking into account the inclination and declination average values of the studied area I = − 9.08°, D = − 2.17°, respectively (Fig. 3b). The RTE method is based on theoretical moving of the magnetization of the investigated area to the magnetization of the "hypothetic region" located in the equator. It causes the magnetization of rocks and magnetic anomaly positioned directly on the center of magnetic body (Leu 1982).
The data transformation was performed using the MAG-MAP™ extension of the Oasis montaj software. The operations in MAGMAP are a range of linear and nonlinear mathematical algorithms which selectively enhance the anomalies related to the groups of the interesting geological source (Milligan and Gunn 1997). They also eliminate the influence of the inclined magnetization occurring in the study region and make the interpretations to be easier. The two-dimensional fast Fourier transform (2D-FFT) is particularly useful in the transformation from the frequency domain to the wave number domain and also for the derivative calculation (Telford et al. 1990). Subsequently, the RTE TMI anomaly was divided into regional and residual anomaly maps.
The regional magnetic anomaly ( Fig. 4a) was calculated from the RTE magnetic anomaly using low-pass Gauss and fifth-order Butterworth filters. The residual anomaly map ( Fig. 4b) was obtained by the subtraction of the regional magnetic anomaly from the TMI anomaly reduced to the magnetic equator (RTE). In the residual map, some enhancements and transformations were performed to extract the near-surface geology.
Various vertical derivatives of the magnetic field can be computed by multiplying the amplitude spectra of the field by a factor of the form: where n is the order of the vertical derivative and (u, v) is the wave number corresponding to the (x, y) directions. The first vertical derivative of the RTE anomaly map for n = 1 is presented in Fig. 5a, which is showing significant enhancement of the near-surface structures and attenuation of the deeper objects.
The influence of the deeper magnetic responses was also examined using the upward continuation algorithm of the MAGMAP extension. The upward continued ∆T (the total magnetic field anomaly) at higher level (z = h) is given by (Everett 2013) The RTE anomaly map resulted from the upward continuation to 1000 m (Fig. 5b) clearly reflects only the large structures and is consistent with the regional anomaly map (Fig. 4a).
For the delineation of edges of similar magnetic responses which may coincide with lithological boundary, the Magnetic anomaly maps; a the regional anomaly map, b the residual anomaly map, cutoff wavelength 5000 in both cases Fig. 5 a First vertical derivative of the RTE anomaly map, b upward continuation to 1000 m of the RTE anomaly map analytical signal filtering was performed. Analytical signal is a process of combining the three directional gradients of the magnetic field to obtain the total gradient and to remove the complexities of derivatives response. The analytical signal A(x, y) at point (x, y) is given by (Dentith and Mudge 2014): In addition, a simple magnetic model was generated, constituting a significant part of the quantitative magnetic interpretation (Foss 2011;Petecki and Rosowiecka 2017). The original TMI anomaly data set (with no reduction or transformations of the magnetic data) was used to create a magnetic model along the profile AB with 78,000 m length.

Aeroradiometric data
The radiometric data were measured at a rate of 1 s using a gamma 512 channel spectrometer "Exploranium GR-820" coupled with a solid scintillation detector. The detector is composed of eight crystals with a total volume of 1280 cubic ins. Each crystal has individual pole-zero cancellation, semi-Gaussian shaping and advanced baseline restoration circuitry. In the spectrometer, there is a special channel to record the cosmic rays with high-energy, special analyzer for pulse pileup of signals from single events and reduce the Compton effect. Detector resolution is automatically computed for each crystal during peak analysis and is displayed for operator monitoring when required. The summed down resolution is also output on the data stream. The detailed technical parameters of this tool can be found in the Web site www.gpxsu rveys .com.au and the monograph published by Foote and Humphrey (1976). The individual channels show the absolute concentrations of the K (%), eTh (ppm) and eU (ppm). Based on the obtained data, the lithological variations were characterized. Moreover, the enhanced aerodatasets were analyzed statistically and presented in various formats such as tables, graphs, rose diagram, pseudo-cross section and maps for interpretations.

Results and discussion
In this section, the results of the processing of the datasets are presented. Deductions made are based on responses influenced by variations in lithology, alterations, mineralogy and microstructures which form major control of geophysical methods of investigation (Dentith and Mudge 2014). The important considerations to interpretations made in this study are the effects of near-surface covers, regolith and deep responses which may be distorted or attenuated.

Digital topography map
Digital topography variation of the study area was generated from digital elevation model (DEM) of NASA short radar topography mission (SRTM) data and displayed as image which indicates variation in elevation in meters. Digital elevation data are useful for identifying areas of erosion and sedimentation, where dissected cover is likely as a guide to where bedrock responses are most likely (Dentith and Mudge 2014). The southern half of the area is occupied by high-elevation terrain, ranging from 280 to 481 m, while the northern half has an elevation range of 180 to 260 m above sea level (Fig. 6). Inferences from elevation variation coincided with some of the responses obtained from the magnetic data (Fig. 3b). The southern zones, particularly SE part, are of very high elevation which may indicate the region of volcanic igneous activities and this may have caused the lifting of the surface.

Aeromagnetic results
In general, the positive magnetic anomalies with maximum intensity value about 5.0 nT dominate in the south and negative anomalies with the lowest intensities about − 115 nT are observed in the north (Fig. 3b). The difference between the highest and lowest magnetic intensities (amplitude) is equal to 120 nT, suggesting that in the study area there are predominantly rocks of weak or even very weak magnetic properties (magnetic susceptibility, magnetization). This type of magnetic property is characteristic of occurrence of acid igneous rock or sedimentary rocks. Sources of the positive magnetic anomalies, occurring near the southern boundary of the studied sheet, may be mafic metavolcanic rocks (e.g., amphibolites) and ultramafic rocks (e.g., peridotite), while the sources of the negative magnetic anomalies characterized of the northern parts suggest the signature of metasedimentary rocks (e.g., phyllites, schist). However, in the northwestern part there is the area of the high enough magnetic intensities, within the negative magnetic intensities regions which may suggest an occurrence of gabbro dykes or faults. The main magnetic anomalies are trending in the SW-NE direction indicating the main geological structures of the area. In the center of the area, there is a magnetic anomaly with ~ 30 km of length trending SW-NE direction; this anomaly could be considered as a fault zone (Fig. 3b, with arrow). Moreover, the map resulted from the upward continuation processing using the Gauss filter (Fig. 4a) indicates the following observations: 1. The positive weak magnetic anomalies on the regional magnetic map are dominant in the southern part of the sheet, and in the NW part several single positive anoma-lies are also observed, but their amplitudes are significantly lower in comparison with that in the southern part. 2. The magnetic lineament zone occurring in the central part probably corresponds to a near-surface tectonic unit (fault). The SW-NE trending magnetic lineament almost disappears after performing the upward continuation with the value of 3000 m (figure not shown in this article).
The first vertical derivative map (Fig. 5a) consists mostly of the local anomalies trending in the SW-NE directions and the near-surface structures which may include lineaments such as fractures and dykes crosscutting the country rocks (Fig. 2). The near-surface structural map of the entire area was generated from the vertical derivative analysis as shown in Fig. 7a. Plot of the orientations of the linear structures on rose diagram (Fig. 7b) shows a general orientation in the ENE-WSW direction.
From the analytical signal (AS) filtering of the magnetic data (Fig. 3b), edges of similar magnetic responses are discernable and are thought to delineate lithological boundaries that produce the responses (Fig. 8a). The edges have enabled the characterization of the entire study area into three major lithology zones (Fig. 8b).
The three major lithological mappable zones are: acid igneous rocks that occupied almost all the four corners of the entire study area (pink), the metasedimentary rocks that formed the dominant lithology extending from the central part to the northern border (yellow), and the migmatite rocks that formed large pockets within the whole area (grey). In Fig. 8b, the igneous rocks include the porphyritic granite, gabbro and granite; the metasediment group is composed of intervolcanic metasediments and metasediment rocks. The rocks in each group are characterized by the similarity in magnetic properties and are not easily distinguished on the analytical signal map. What is more, migmatites are often characterized by a wide range of magnetic properties and occur on the large part of the sheet (Fig. 9).
In order to confirm the assumptions and recognize the magnetic susceptibility (S) of the rocks as well as to view the depth to the top of the basement, a cross section of the profile AB (Fig. 10) was drawn as an example. The profile AB is selected because it cuts diagonally the sheet with the NW-SE direction and crosses most of the main geological structures of the area. The TMI values measured along this profile were taken to generate the two-dimensional (2D) magnetic model using the GM-SYS Profile Modelling extension of Oasis montaj. Based on the geological data as found in the literature, the magnetic model built on the profile AB consists of a sedimentary formation and a basement, which relative to AS map is divided into seven blocks with different susceptibility (Fig. 8a). Since there was no magnetic susceptibility (S) data from this   Magnetic susceptibility values of gabbro dykes and acid igneous rocks amount to about 0.1 and 0.00001, respectively. It is in agreement with the statement that mafic intrusive igneous rocks (gabbros) are definitely magnetically stronger in comparison with acid igneous rocks (e.g., granite without magnetite grains) (Parasnis 1986). Migmatites in the model have strong responses with high negative magnetic anomalies at the SE end part of the AB profile (Fig. 10), and their magnetic susceptibility varies from 0.021 to 0.087. The magnetic susceptibility of metasedimentary rocks ranges from 0.00001 to 0.003. In the southern part of the studied area, the basement occurs relatively at shallow depth in comparison with other parts. This assertion is corroborated with information extracted from the digital topography map (Fig. 6). In the central and NW parts, the basement could be zones of fracture/ faults which exhibit undulating surface and much thicker sedimentary cover (Fig. 10).
The magnetic model generally is in agreement with the results of qualitative interpretation regarding the basement morphology and the formations crossed along the profile AB (Figs. 8b, 10).

Aeroradiometric results
The aeroradiometric data distribution of potassium (K), equivalent uranium (eU) and equivalent thorium (eTh) on the study area are presented in Fig. 11a-c, respectively.
The measured concentrations range from 0.2 to 5% for potassium, 0.5-12 ppm corresponds to the 6-150 Bq/kg for uranium and 3-100 ppm correspond to the range of 12 to near 400 Bq/kg for thorium. Generally, the zone of high concentration of all measured natural radionuclides is located in the south of the study area. In this zone, the eTh concentration varies from 30 to 49 ppm, eU from 5 to 10 ppm and K from 2.7 to 4.1% and suggests areas connected with  outcrops of crystalline basement and the presence of the intrusive formation. The high radiometric anomalies are also present on the southeast corner and near center as well as some local points in the north. The local places with high concentrations of K, eU and eTh may be connected with the intrusive suits (granite, pegmatite). The low concentrations of K (0.4-1.2%), eU (0.5-3.5 ppm) and eTh (3-20 ppm) are prevalent; such low radioactivity is more than 60% of the total area, suggesting the presence of the carbonate, silicate formation and quaternary sediments. For better recognition of the lithology, the relations between the concentrations of the radioactive elements measured on some selected profiles (as shown in Fig. 12) are presented in Fig. 13a-f. Generally, the relationship between U and Th concentration is linear, suggesting that the U and Th are accompanying each other. The high concentrations of U and Th are characteristic of the widespread intrusive rocks (granite, syenite) in the studied region as well as in other parts of Nigeria (Adepelumi and Falade 2017). The high concentration of Th and U can be related also to the Fe or another metal element deposits formed by the magmatic activities.
The relations between K, Th or U are more complicated in comparisons with that of U and Th. Since U and Th belong to the actinide elements, they have similar geochemical properties, while K is the alkaline element, so its geochemical characteristics differ from that of U and Th. Generally, potassium concentration ranges from a few tenth percent up to five percent, suggesting no potassium salt deposit. The high concentration of potassium together with high U and Th can indicate the places of the outcrop of the intrusive formation.
Based on the measured concentrations of eU, eTh and K, the rocks in the studied area can be divided into four types of lithology. The intervals of the variation of the U, Th and K for each lithology are summarized in Table 2.
The A type can be characteristic for the most of common rocks, such as metasedimentary rocks and some quaternary sediments. The B type is composed of the rocks or deposit with the highest concentrations of U, Th and K. The rocks could be the igneous, metavolcanic rocks or deposits rich in radioactive elements. The C type has a high concentration of U and Th and moderate concentration of K suggesting the occurrence of the fine-grained clastic, syenite or pegmatites.  Table 2) on the background of eTh concentration map (see Fig. 11b) The D type is characterized by moderate U and Th concentration but relatively high concentration of potassium suggesting the occurrences of metamorphic rocks such as granite-gneiss, migmatite and schists (Dentith and Mudge 2014). The rocks of this type are present majorly in the southern part of the area and some places, where there are outcrops of crystalline basement. All the lithology types are shown in color on the eTh map (Fig. 14). Generally, the formations A, B, C and D are trending in the SW-NE direction conscious with the area tectonics (compare Figs. 7a and 14). The boundaries between the identified rock types presented in Fig. 14 generally are trending in the SE-NW. Figure 15a-c illustrates the maps of eTh/eU, eTh/K and eU/K, respectively.
The average mass concentration of thorium, uranium and potassium in the Earth crust is 10 ppm, 3.5 ppm and 2.5%, respectively, so the average ratio amounts to ~ 3.5 for eTh/eU, ~ 4.0 for eTh/K and ~ 1.4 for eU/K. Reference to the mentioned ratio values, Fig. 15a-c indicates that the gamma radiation of the study region is mostly generated from the isotopes of thorium and uranium decay series. The zones are rich in K (eU/K < 1.4 and eTh/K < 3.5), contribute above 15% of the region area and mostly occur in SE and SW (Fig. 15b, c) and partly cover with the D zones (Fig. 14). According to Dickson and Scott (1977), Ohioma et al. (2017a, b) and Airo (2002), the K-rich zones could indicate the former hydrothermal alteration and can be potential gold deposits. The zones B in Fig. 14 have the similar concentrations of eTh, eU and K, but the ratio of eTh/eU at the Fig. 15 Maps of radionuclides mass ratios a eTh/eU, b eTh/K and c eU/K center place signed as B2 is near 3.0 and significantly lower than that in the south places signed as B1 (eTh/eU: 5.0-7.0). The fact is probably related to reduction in Th in the zone B2 (eTh ~ 36 ppm, eU ~ 9.5 ppm, K ~ 4.1%) in comparison with the zones B1 (eTh ~ 45 ppm, eU ~ 9.5 ppm, K ~ 3%) (Fig. 15a-c). The observation is confirmed by the map of U-Th/3.5 (Fig. 16). The B2 zone is more homogeneous and suggests that it is the magmatic rock without controlled by the geological processes; on the other hand, the B1 zones also are dominated by intrusive formations, but they were inhomogeneous and may be controlled by different geological processes; in consequence, the rocks of the B1 zones could be metamorphosed or fractured.

Conclusions
The aeromagnetic and aeroradiometric data obtained for the Isanlu sheet 225, North Central Nigeria, were interpreted with a view to delineate geological structures and the lithology. The aeromagnetic data delineate main structures trending principally in the SW-NE direction, which is in general consistent with that of structures of similar terrain (schist belts) in Nigeria. In the western part, there are zones with many deeper subsurface structures which are thought to be faults/fractures or possibly veins/dykes trending in the WNW-ESE direction. In the center of the study area, there is a long fault (occurring on ~ 30 km (Figs. 3a,b,4a,b,9)site with arrow) crosscutting different subsurface lithologies.
The analytic signal map shows the zones occurrence of igneous, metasedimentary and migmatite rocks (Fig. 8b).
Aeroradiometric results in the SE corner of the Isanlu sheet 225 depict crystallized hardness formation with very high concentration of eTh and moderate eU but relatively low K. Generally, the concentration of eU is linearly proportional to that of eTh, suggesting the elements accompanying themselves in most rocks present in the studied area. The zones of eU with the concentration of 2-8 ppm and Th from 5 to 35 ppm contribute near 90% of the whole area. The zone with high concentration of eU (8-12 ppm) and eTh (36-95 ppm) dominates the south-center, SE corner and some local places distributed in the NE corner and a few points in the upper parts. These zones can be connected with the outcrop of basement, igneous or crystalline metamorphose suites or the relatively skarn districts connected with precious metals (Ohioma et al. 2017a,b). The zones of low eU, eTh and K covers above 60% of the area and mainly predominate in the north part, indicating the presence of the migmatite-gneiss-quartzite complex of the pelitic schists, phyllites and carbonate rocks and quaternary sediments on the surface. Based on the radiometric data, the rocks can be classified into four types A, B, C and D. The A type is composed of the rocks with relatively low concentrations of the natural radionuclides, the B type with the highest concentration of the natural radioactive elements, the C type with high concentration of both eU and eTh but relatively low potassium, and the type D with moderate eU and eTh but relatively high potassium. From Fig. 16 Map of U-eTh/3.5 the isotope ratio point of view, the zones B can be divided into B1 and B2. The rocks in the B1 zones probably are inhomogeneous and fractured, while the rocks in the B2 zone are homogeneous. The rocks in the D zones are most probably hydrothermal-altered suggesting gold potential elevated regions.