Residual aeromagnetic data of parts of Southeastern Nigerian sedimentary basin were reduced to the equator and subjected to magnetic vector inversion and spectral analysis. Average depths of source ensembles from spectral analysis were used to compute depth to magnetic tops (Z), base of the magnetic layer (Curie Point t Depth (CPD)), and estimate geothermal gradient and heat flow required for the evaluation of the geothermal resources of the study area. Results from spectral analysis showed depths to the top of the magnetic source ranging between 0.45 km and 1.90 km; centroid depths of 4 km - 7.87 km and CPD of between 6.15 km and 14.19 km. The CPD were used to estimate geothermal gradients which ranged from 20.3°C/km to 50.0°C/km 2 2 and corresponding heat flow values of 34.9 mW/m to 105 mW/m , utilizing an average thermal conductivity -1 -1 of 2.15 Wm k . Ezzagu (Ogboji), Amanator-Isu, Azuinyaba, Nkalagu, Amagunze, Nta-Nselle, Nnam, Akorfornor environs are situated within regions of high geothermal gradients (>38°C/Km) with models delineated beneath these regions using 3D Magnetic Vector Inversion, having dominant NW-SE and NE-SW trends at shallow and greater depths of <1km to >7 km bsl. Based on VES and 2D imaging models the geothermal system in Alok can be classified as Hot Dry Rock (HDR) type, which may likely have emanated from fracture systems. There is prospect for the development of geothermal energy in the study area. Keywords: Airborne Magnetics, Magnetic Vector Inversion, Geothermal Gradient, Heat Flow, Curie Point Depth, Geothermal Energy.