Abstarct: The purpose of this study is to estimate the spatial and geochemical distribution of heavy mextal pollution in soils around the Khatoon abad copper smelter located in east of Shahre babak, Iran. Elements were chosen baxsed on geology of region and also copper smelters potential for pollution. Total mextal concentration in samples along with soil physicochemical properties were determined. To assess the quantitative severity of pollution, geochemical indices such as EF, Igeo, Contamination factor and polluted load were calculated too. The results revealed high contamination for arsenic, considerable for antimony and silver and moderately for cadmium and copper. Sequential extraction and multivariate statistical analysis were done to determine the origin of elements baxsed on increasing distance from source of contamination (smelter). The results showed that the highest concentration was observed in samples taken from a distance of 1 kilometer of the smelter specially in the north east (prevailing wind direction) and the concentration will decrease when distance from the smelter and the depth of sampling increases. According to the analysis of the dust emitted from stack of furnace, the source of elements silver, arsenic, cadmium, copper, molybdenum, lead, sulfur, antimony and zinc is anthropogenic. Elements Al, Fe, Mn, Li and Sc showed a constant trend with increasing distance from the smelter and have low concentration in dust and concentrate, so they have geogenic origin. Since the study area is located in porphyritic Cu-Mo mineralization region, and also according to high concentration of these elements in dust and concentrate, Cu and Mo probably have two origin, anthropogenic and geogenic. These results have been proved through sequential extraction procedure, too. So that elements As, Cd. Zn and Pb were bounded to non-residual fractions while Cu and Mo were associated with residual fraction and primary phases. Finally to prepare spatial distribution map of elements, the discrete and ordinary kriging approach was used in GIS. Calculating the interpolation stimate error showed the most appropriate approach for spatial distribution map for Cu, Cd, Mo and Zn is ordinary kriging and for Pb and As is discrete kriging method. Spatial distribution maps showed a decreasing trend for concentration along with increasing distance.