Sheshtamad area is located in the western part of Khorasan Razavi province and includes Sheshtamad city and its surrounding villages. A variety of sedimentary and igneous rock units, mainly belonging to the Cretaceous or younger age, are found in the area. A series of orogenic movements, lithofacies changes and magmatic activities led to the distribution of heavy elements that can cause environmental problems in the area. This study aimed to assess the heavy metal contamination of the soil in the extent of the study area. It has been attempted to investigate the origin of pollution and its impact on the environment using geochemical data.

In this study distribution and possibility of contamination from arsenic (As), lead (Pb), zinc (Zn), copper (Cu), chromium (Cr), nickel (Ni), cobalt (Co), antimony (Sb), strontium (Sr) and vanadium (V) were investigated. These variables were determined using chemical analysis of the sediments. For this purpose, network sampling was carried out at 1.5×1.5 km and 3 to 4 samples were taken from each cell. Throughout the study area (geological sheet of 1:100,000 Sheshtamad), 1248 samples of sediments were collected by Jiangxi Company and analyzed by ICP-MS method. The analysis results were summarized in ArcGIS software and elemental concentration zoning was performed by Kriging and Inverse Distance Weighting methods. Some indices have been proposed to evaluate heavy metal contamination in sediments and soils. In order to determine soil contamination with heavy metals, some parameters such as Contamination Factor and Pollution Index were calculated. Multivariate statistical analyses such as correlation analysis have been applied to identify the geochemical behaviors of different geochemical groups.

Surface distribution of As with values above 270 mg kg-1 was observed in the drainages of the southern part of the study area. Sr and Sb were measured with values above 2900 mg kg-1 in the drainages from south to eastern part, and above 4 mg kg-1 in south and central drainages, respectively. Surface dispersion of Zn above 210 mg kg-1 was present in the drainages from central, south, southwest and northwestern parts of the region. Cu and Pb were measured with values greater than 240 mg kg-1 in North and Northwestern parts, and greater than 110 mg kg-1 in East, Central and southwest, respectively. Surface distribution of Co element with values above 58 mg kg-1 from drainages in central, north and southwest and V with values higher than 85 mg kg-1 in central and southwest part were also observed. Finally, distribution of Cr element with values above 7400 mg kg-1 in east central and southwestern part and Ni with values higher than 890 mg kg-1 in center and north section were found. Pb had significant contamination only in a small area of the central part. Cu also showed little pollution only in a small area of the center of the region on the fringes of mining activities. According to the Pollution Index, As had a high pollution in the south and west parts. The central and southwestern parts of the study area show moderate Cr pollution, based on the Contamination Factor. Based on the correlation matrix used to determine the relationship of heavy elements with each other, some elements have a positive correlation with each other, which may indicate their common origin. For example, Co had a positive correlation with Ni and Cr indicating a common origin associated with the distribution of mafic rock units in the region. Pb and Zn have a positive correlation with each other as well as with Co and V, whereas the absence or negative …