Regional estimation of Curie-point depths and succeeding geothermal parameters from recently acquired high-resolution aeromagnetic data of the entire Bida Basin , north-central Nigeria

A regional estimation of Curie-point depths (CPDs) and succeeding geothermal gradients and subsurface crustal heat flow has been carried out from the spectral centroid analysis of the recently acquired highresolution aeromagnetic (HRAM) data of the entire Bida Basin in north-central Nigeria. The HRAM data were divided into 28 overlapping blocks, and each block was analysed to obtain depths to the top, centroid, and bottom of the magnetic sources. The depth values were then used to assess the CPD, geothermal gradient, and subsurface crustal heat flow in the basin. The result shows that the CPD varies between 15.57 and 29.62 km with an average of 21.65 km, the geothermal gradient varies between 19.58 and 37.25 C km−1 with an average of 27.25 C km−1, and the crustal heat flow varies between 48.41 and 93.12 mW m−2 with an average of 68.80 mW m−2. Geodynamic processes are mainly controlled by the thermal structure of the Earth’s crust; therefore this study is important for appraisal of the geo-processes, rheology, and understanding of the heat flow variations in the Bida Basin, north-central Nigeria.


Introduction
This study aims at quantitative estimation of regional Curiepoint depth (CPD) and succeeding geothermal structures, namely geothermal gradients and subsurface crustal heat flow anomalies, in the whole of Bida Basin in north-central Nigeria using the spectral centroid analysis of the recently acquired high-resolution airborne magnetic (HRAM) data.The HRAM surveys were carried out by Fugro Airborne Survey Limited for the Nigerian Geological Survey Agency (NGSA) between 2004 and 2009.Acquisition, processing, and compilation of the HRAM data were jointly financed by the Federal Government of Nigeria and the World Bank as part of the Sustainable Management for Mineral Resources Project (SMMRP) in Nigeria.
Several studies have shown that regional magnetic data can be used extensively to determine the thermal structure of the Earth's crust in various geologic environments (Spec-tor and Grant, 1970;Bhattacharyya andLeu, 1975, 1977;Byerly and Stolt, 1977;Blakely and Hassanzadeh, 1981;Okubo et al., 1985Okubo et al., , 2003;;Blakely, 1988Blakely, , 1995;;Maus et al., 1997;Tanaka et al., 1999;Chiozzi et al., 2005;Eppelbaum and Pilchin, 2006;Ross et al., 2006;Ravat et al., 2007;Trifonova et al., 2009;Gabriel et al., 2011Gabriel et al., , 2012;;Bansal et al., 2011Bansal et al., , 2013Bansal et al., , 2016;;Nabi, 2012;Hsieh et al., 2014;Nwankwo and Shehu, 2015;etc.).For example, dominant magnetic minerals in the Earth's crust pass from ferromagnetic to paramagnetic state at temperature, commonly called Curie-point temperature (CPT).Magnetite (Fe 3 O 4 ) is the most common magnetic material in igneous rocks and has an approximate CPT value of 580 • C (Stacey, 1977).At temperature above CPT, the thermal agitation causes the spontaneous alignment of the various domains in the mineral to be destroyed (or randomized) to the extent that the ferromagnetic minerals become totally paramagnetic (Langel and Hinze, 1998).Therefore CPD, which is defined as the depth at which CPT is reached within the subsurface, can be considered as an index of depth to the bottom of magnetic sources (DBMS) and can consequently be calculated from geomagnetic anomalies (Bansal et al., 2011(Bansal et al., , 2013;;Hsieh et al., 2014).However, in some circumstances DBMS can be caused by contrasts in lithology instead of CPT and may not necessarily coincide with CPT in detail (Bansal et al., 2011;Trifonova et al., 2009).For instance, Trifonova et al. (2009) opined that even if the spectral method provides a good estimate of DBMS there is no assurance that it represents the CPD.They reasoned that a variety of geologic reasons exist for truncated magnetic sources that are unrelated to crustal temperatures; for example, a sequence of relatively non-magnetic sediments below young volcanic material may limit the depth of magnetic sources regardless of the CPT, and another reason is the variety of magnetic minerals like titanomagnetite (Fe 3−x Ti x O 4 ).Titanomagnetite is the most important iron oxide in crustal magnetic sources; it has a CPT that is strongly influenced by the amount of titanium and ranges from 150 to 580 • C. In some geologic environments, alloys of iron with CPTs in excess of 620 • C may be significant contributors to magnetic anomalies.In spite of these limitations, many studies (Tanaka et al., 1999;Trifonova et al., 2009;Bansal et al., 2011;Hsieh et al., 2014;etc.)have reasonably used DBMS as an estimate of CPD and therefore serve as a proxy for temperature at depth.Again, Trifonova et al. (2009) pointed out that several studies have identified low-titanium titanomagnetite as the dominant magnetic phase, and CPTs at these depths are estimated to be between 575 and 600 • C.This confirms the estimated value of 580 • C by Stacey (1977) as the case in this study.Another important justification is that DBMS/CPD estimations can similarly be used to complement geothermal data in re-gions where deep boreholes are unavailable (Chapman and Furlong, 1992;Ross et al., 2006;Bansal et al., 2011Bansal et al., , 2013)).
The Bida Basin is the least studied of all Nigeria's inland frontier basins.To date, the basin has no information on seismicity, no exploratory wells have penetrated its sequences, and deep crustal data are limited (Obaje et al., 2015).Therefore, this present work is expected to contribute immensely to a better understanding of the geothermal structures and geodynamic processes in the entire Bida Basin in north-central Nigeria.

Location and geology of the study area
The Bida Basin (also known as the Middle Niger Basin or Nupe Basin) is an elongated NW-SE-trending depression perpendicular to the main axis of the Benue Basin of Nigeria.The entire basin (Fig. 1) is bounded by latitudes 8 • 00 N and 10 • 30 N and longitudes 4 • 30 E and 7 • 30 E, and it covers an area of approximately 90 750 km 2 .
The geology of the Bida Basin is believed to be a gentle down-warped shallow trough filled with Campanian-Maastrichtian marine to fluviatile strata.The strata are also believed to be more than 300 m thick (Adeleye, 1976).Those with marine affinity, the limestones, often form cappings (under variable thickness of laterites) to the means of the basin.Some form prominent intermediate breaks of slope along the mesa walls.The buried basement complex probably has a high relief, and the sedimentary formations have been shown to be about 2 km thick (Obaje et al., 2015) with a constituted post-tectonic molasse facies and thin marine strata, which are all unfolded.
The basin might also be regarded as the north-western extension of the Anambra Basin, which is found in the southeast, both of which were major depocentres during the sec-  ond major sedimentary cycle of southern Nigeria in the Upper Cretaceous (Obaje, 2009).
The stratigraphy of the basin consists of mainly Patti, Lokoja, and Agbaja formations (Ikumbur et al., 2013).According to Akande et al. (2005), Patti Formation is the only stratigraphic unit containing carbonaceous shale in the basin and is sandwiched between the older Campanian-Maastrichtian Lokoja Formation -which contains conglomerates, sandstones, and claystones -and the younger Agbaja Formation, which comprises mostly ironstones.
Tectonically, rift hypothesis has been proposed as the possible origin and evolution of the basin (Kogbe et al., 1983).Some researchers have inferred that the rifting in the basin started in the Upper Cretaceous and that the process started initially from the Benue Trough in the early Cretaceous and eventually spread to other rifted neighbouring basins in the Cretaceous period (Kogbe et al., 1983).
The mathematical models of the centroid method are based on the examination of the shape of isolated magnetic anomalies introduced by Bhattacharyya andLeu (1975, 1977) and the study of the statistical properties of magnetic ensembles by Spector and Grant (1970).Blakely (1995) subsequently introduced power spectral density of total magnetic field, φ T (k x , k y ), as where k x and k y are wave numbers in x and y direction, φ M (k x , k y ) is the power spectra of the magnetization, C M is a constant, θ M and θ f are factors for magnetization direction and geomagnetic field direction, and Z b and Z t are depths to bottom and top of magnetic layer respectively.
If the layer's magnetization, M(x, y), is a random function of x, y, it implies that φ M (k x , k y ) is a constant, and therefore the azimuthally averaged power spectrum, φ(|k|), would be given as (2) The depth to the top of the magnetic source is therefore derived from the slope of the high-wave-number portion of the power spectrum as where P (k) is the azimuthally averaged power spectrum, k is the wave number (2π km −1 ), A is a constant, and Z t is the depth to the top of magnetic sources.The centroid depth of magnetic sources can also be calculated from the low-wave-number portion of the wavenumber-scaled power spectrum as (Tanaka et al., 1999) where B is a constant and Z o is the centroid depth of magnetic sources.
The depth to the bottom of the magnetic source (Z b ) can subsequently be obtained from the relation (Okubo et al., 1985) (5) Using the depth to the bottom of magnetic sources (Z b ), the geothermal gradient (dT /dz) can be estimated as (Tanaka et al., 1999;Ross et al., 2006) where θ c is the Curie temperature.
Next, using Z b and dT /dz, the heat flow (q z ) can similarly be estimated as (Okubo et al., 1985) where σ is thermal conductivity.Thermal conductivity of 2.5 W m −1 • C −1 as the average for igneous rocks and a Curie temperature of 580 • C (Stacey, 1977;Trifonova et al., 2009) are used as standard in this work.

Data acquisition and analysis
The regional airborne magnetic surveys over the entire Bida Basin and adjoining areas were carried out using 3 Scintrex   138-143, 159-166, 180-187, 203-208, and 224-229) on a scale of 1 : 100 000, covering a total area of 130 075 km 2 , were used in this work.All data, covering the entire Bida Basin and adjoining areas, were procured from NGSA as a composite residual total magnetic field intensity (RTMI) map (Fig. 2).Regional correction, which was based on the International Geomagnetic Reference Field (IGRF -11) derived to spherical harmonic degree 13, was carried out by the NGSA prior to the publication of the map.
The composite residual map was then divided into 28 overlapping square blocks (Fig. 3), for the purpose of spectral analysis.Each block covers a square area of 200 km × 200 km.Since magnetic source bodies having bases deeper than window length/2π may not be properly resolved by the spectral method (Shuey et al., 1977), a window length of 200 km is found to be suitable in this study.
www.geoth-energ-sci.net/5/1/2017/Geoth.Energ.Sci., 5, 1-9, 2017  The estimated results are shown in Table 1.The table shows that the estimated CPD varies from 15.57 to 29.62 km with an average of 21.65 km.The CPD isotherm map for the Bida Basin is consequently shown in Fig. 5. CPD varies greatly with different geological settings (Tanaka et al., 1999;Salk et al., 2005).Tanaka et al. (1999), after a compilation of CPD results from several researchers across the globe, inferred that volcanic, tectonic, and associated geodynamic environments have CPD shallower than 10 km, while CPDs ranging between 15 and 25 km are as a result of island arcs and ridges, and deeper than 25 km in plateaus and trenches.Figure 5 also shows that the CPD values in the Bida Basin trend mostly in the NE-SW direction, with the shallowest portion (less than 15 km) in the north-western part of the basin; the CPD extends and deepens both north-eastward and south-eastward in the basin, with its deepest depth (about 32 km) in the north-eastern parts.
Previous studies have shown that two major regional fault lines (namely St Paul and Romanche) are likely to have traversed the basin; these are believed to be extensions of the  onshore lineaments in West Africa, which are part of the major weakness in the crust that predates the opening of the Atlantic Ocean, and were reactivated in the early stages of continental rifting (Fig. 6) (Ajakaiye et al., 1991;Buser, 1966).Thus, the relatively low CPD values over the central portion of the basin may be fairly consistent with probable positions of St Paul and Romanche palaeofracture zones.Table 1 similarly discloses that the geothermal gradients in the basin vary between 19.58 and 37.25 • C km −1 with an average of 27.25 • C km −1 , while the crustal heat flow varies between 48.41 and 93.12 mW m −2 with an average of 68.80 mW m −2 .Contour maps for the geothermal gradient and heat flow are shown in Figs.7 and 8 respectively.The geothermal gradient map also exhibits mostly NE-SW trending.The observed major trends are similar to the regional trending faults in the basin.The lowest values for the geothermal gradient were found in the south-western portion of the basin.The north-westward trend of gradient increase was found to result in a maximum value of 42 • C km −1 in the north-western part.The minimum heat flow value required for considerable generation of geothermal energy is approximately 60 mW m −2 , whereas values ranging from 80 to 100 mW m −2 and above indicate anomalous geothermal conditions (Jessop et al., 1976).Crustal heat flow in the basin also exhibits NE-SW trending, while the derived amounts increase from the central portion towards the north-west, with maximum values above 90 mW m −2 observed in the northcentral portion.This portion signifies an anomalous crustal thermal state and, therefore, is recommended for further investigations.

Conclusion
The newly acquired high-resolution aeromagnetic anomaly data over the Bida Basin, north-central Nigeria, have been analysed to estimate the Curie-point depths, geothermal gradients, and near-surface crustal heat flow.The result shows that the CPD varies between 15.57 and 29.62 km with an average of 21.65 km, the geothermal gradient varies between 19.58 and 37.25 • C km −1 with an average of 27.25 • C km −1 , and the crustal heat flow varies between 48.41 and 93.12 mW m −2 with an average of 68.80 mW m −2 .
Regions are observed in the basin with shallow Curie-point depths (below 15 km) and corresponding high heat flows (above 80 mW m −2 ), thus suggesting anomalous geothermal conditions (Jessop et al., 1976).Hence, further detailed studies are recommended in such regions.Finally, oftentimes, direct crustal temperature measurements may not be too feasible for regional studies; hence, the derived geothermal gradients suffice for the entire basin.Moreover, geodynamic processes are mainly controlled by the thermal structure of the Earth's crust; therefore this study is anticipated to contribute significantly to the quantitative appraisal of the geoprocesses, rheology, and understanding of the heat flux variations in the Bida Basin in north-central Nigeria.

Figure 1 .
Figure 1.Geological map of Nigeria showing location of the Bida Basin (after Nwankwo and Shehu, 2009).

Figure 2 .
Figure 2. Residual total magnetic intensity map of the entire Bida Basin with superimposed federal survey half-degree sheets and showing major towns flown over.A constant TMI value of 33 000 nT was removed.

Figure 3 .
Figure 3. Sketch of procedure for achieving overlapping blocks.

Figure 5 .
Figure 5. CPD map of the study area.

Figure 6 .
Figure 6.West African Rift System (WARS) and Central African Rift System (CARS) (after Heine et al., 2013).Red circle indicates position of the Bida Basin with extensions of inferred major fracture zones -St Paul, Romanche, and Chain drawn with broken blue lines.These extensions had been suggested by Ajakaiye et al. (1991).

Figure 7 .
Figure 7. Geothermal gradient map of the study area.

Figure 8 .
Figure 8.Heat flow map of the study area.

Table 1 .
Estimated Curie-point depths and succeeding geothermal parameters.