Abstract: Background and Objective: The impact of pollution in the vicinity of overcrowded cities and from industrial effluents has reached a disturbing magnitude and is arousing public awareness. Therefore, the objective of this research was to determine the concentrations and characterize the spatial distribution of arsenic (As), cadmium (Cd), chromium (Cr), copper (Cu), nickel (Ni), lead (Pb) and zinc (Zn) heavy metals in soil of Inegol Plain by using geostatistical approaches combined with GIS based on 17 soil samples in both winter and summer periods. Methodology: Ordinary kriging interpolation technique was used to estimate a variable at an unmeasured location from observed values at nearer locations. Moran’s I and Getis-Ord Gi were used to investigate spatial autocorrelation of soil pollutants. Results: Very high heavy metal concentration levels were found throughout the study area. The local spatial autocorrelation analysis made clear that all heavy metals had important hotspots and coldspots and these areas should be targeted by best management practices to reduce heavy metal loads. Attention must be given to industrial activities performed within these high-risk areas. Conclusion: Results from this study can be used by ecologists, planners and managers to protect soil supplies.