This study employs geospatial statistical technique to assess the spatial distribution of heavy metals along the Nun River. Core sediment samples were collected from relatively undisturbed areas (twenty-five different stations) using Uwitec Triple sediment cutter. The rectangular coordinates of the sediment sample location were determined with the aid of Germin handheld GPS receiver. The concentrations of cadmium, lead chromium and zinc present in the sediments was determined with the aid of an atomic absorption spectrophotometer. For geospatial analysis, five semi-variogram models (stable, circular, spherical, exponential and K-Bessel) were fitted for each of the four critical parameters (heavy metals). In addition, four goodness-of-fit statistics (mean square error, root mean square error, root mean square standardized error and average standard error) were utilized to decide the most suitable model used to develop the final prediction map for each parameter. From the results obtained, it was observed that; regions with red color codes signify higher concentrations of cadmium, lead, chromium and zinc. Further assessment of the results showed that Otuan, Obeleli, Angiama, Odobio, Kasama, Akedda and Akele experienced high concentration of cadmium while Tombia, Ewoi, Abilabio, Agudama and Yenikpa experienced high concentration of lead.