Document Type : Original Article


Department of Soil Sciences, Faculty of Agriculture, Payame Noor University, Kerman, Iran


The study of the spatial distribution of heavy metals is very important in land management and planning. The geostatistics theory is used to estimate spatial variables in unmeasured points. By using the CoKriging estimator and using information about the correlation between variables, a more accurate estimation of a variable (main variable) can be obtained using auxiliary information (secondary variables). Due to the importance of Isfahan Province because of the large population living in this province and the extent of agricultural and industrial activities, there is a need for sufficient information on the distribution of selenium in this area. The present study investigated the distribution of selenium and provided maps of its distribution in the surface soils of this province.
Material and methods:
The study region has an area of about 6800 km2 and covers a large part of the agricultural lands and the major industrial centers of Isfahan. Samples were taken from depths of 0 to 20 cm from the soil surface and 255 samples from the whole area in 4 × 4 km networks by using a stratified random sampling design. All soil samples were analyzed for pH, electrical conductivity (EC), sand, clay, silt content, and organic matter percentage. Total selenium concentrations were measured by XRF in only 72 samples. Statistical and geostatistical calculations were performed by SPSS and WinGslib, respectively, and maps were drawn by ILWIS and Surfer software.
Results and discussion:
The mean concentration of selenium in the study area was 0.63 mg/kg with a min and max amount of 0.5 and 1.5 mg/kg, respectively. There was a positive and significant correlation between the selenium concentration and soil’s EC. Therefore, EC was used to reduce the selenium estimation error by the CoKriging method. While toxicity level of selenium starts at 5 (mg/kg) for humans and animals, in this study, deficiency threshold for total soil selenium was 0.6 mg/kg and, therefore, the study area can be identified as selenium-deficient because more than 50% of the area had a concentration of less than 0.6 mg/kg and about 12% of the area has a high concentration of selenium in unpolluted soils. No toxicity of this element to humans and livestock was observed. Large parts of Isfahan Province have an arid climate with alkaline pH soils. Given the low levels of selenium in both agricultural and non-agricultural lands, it seems that the native material of this area, which is often the alluvial constituent of the Zayandehrood River, is poor in selenium. On the other hand, the deficiency of this element in agricultural lands indicates that farmers in these areas do not use selenium-containing fertilizers. The highest concentration of selenium was located in the city of Mobarake. There is a huge steel industry in this area. The dominant wind direction of the region is introduced as a threat factor in the transfer of selenium-containing atmospheric emissions from the industrial areas of Mobarka to more populated areas of the province
Industries are more concentrated in the city of Isfahan and its surrounding areas and thus selenium concentration is higher in these areas than the neighboring areas. The steel industry in Mobarakeh is an agent for the introduction of selenium into the atmosphere and finally the soil. Winds probably play a role in transmitting selenium to the northeastern part of the region. However, no toxicity of this element was observed in all the studied areas but even selenium deficiency was observed. Therefore, using selenium-rich fertilizers in the region is recommended.


  1. Amini, M., Afyuni, M., Khademi, H., Abbaspour, K.C. and Schulin, R., 2005. Mapping risk of cadmium and lead contamination to human health in soils of Central Iran. Science of Total Environment. 347, 64-77.
  2. Andronikov, S.V., Davidson, D.A. and Spiers, R.B., 2000. Variability in contamination by heavy metals: sampling implications. Water, Air and Soil Pollution. 120, 29-45.
  3. Askarian Amiri, A. and Sadegh Poor, M., 2018. Effect of Selenium Contaminant Transfer on Soil. In Proceedings 6th International Conference on Engineering Science and Technology (ICIEC 2017), 27th November, Mantissa College, Malaysia.
  4. Ayoubi, S.A., Mohamad Zamani, M. and Khormali, F., 2007. Prediction total N by organic matter content using some geostatistic approaches in part of farm land of Sorkhankalateh, Golestan province. Journal of Agricultural Science and Natural Researches. 14(4), 78-87.
  5. Bouyoucos, G.J., 1951. A recalibration of the hydrometer method for making mechanical analyses of soils. Agronomy Journal. 43, 434-438.
  6. Chen, T., Chang, Q., Liu, J., Clevers, J.G.P.W. and Kooistra, L., 2016. Identification of soil heavy metal sources and improvement in spatial mapping based on soil spectral information: a case study in northwest China. Science of the Total Environment. 565, 155–164.
  7. Clay, S.A., 2011. GIS Application in Agriculture. Tayler and Francis, CRC Press, USA.
  8. Combs, G.F., 2001. Selenium in global food systems. British Journal of Nutrition. 85, 517–47.
  9. Ersahin, S., 2003. Comparing ordinary kriging and Cokriging to estimate infiltration rate. Journal of Soil Science Society American. 67, 1848-1855.
  10. Fordyce, F., 2005. Selenium Deficiency and Toxicity in The Environment. In Essentials of Medical Geology. Springer, Dordrecht, Netherlands.
  11. Goovaerts, P., 1997. Geostatistics for Natural Resources Evaluation. Oxford University Press, UK.
  12. Gupta, U.C. and Gupta, S.C., 2010. Selenium deficiency in soils and crops and its impact on animal and human health. Current Nutrition & Food Science. 6, 268-280.
  13. Hasani Pak, A.A., 1997. Geostatistics. Tehran University Press, Tehran, Iran.
  14. Isaaks, E.H. and Srivastava, R.M., 1989. An introduction to applied geostatistics. Oxford University Press, UK.
  15. Kabata-Pendias, A., 2011. Trace Elements in Soils and Plants, 4th ed. Taylor and Francis Group, LLC., USA.
  16. Khosravi, Y., Zamani, A.A., Parizanganeh, A.H. and Yaftian, M.R., 2018. Assessment of spatial distribution pattern of heavy metals surrounding a lead and zinc production plant in Zanjan Province, Iran. Geoderma Regional. 12, 10-17.
  17. Lark, R.M., Ander, E.L., Cave, M.R. and Scanlon, R.P., 2014. Mapping trace element deficiency by Cokriging from regional geochemical soil data: a case study on cobalt for grazing sheep in Ireland. Geoderma. 226, 64–78.
  18. Lemly, A.D., 2004. Aquatic selenium pollution is a global environmental safety issue. Ecotoxicology Environmental Safety. 59, 44-56.
  19. Lindsay, W.L., 1979. Chemical Equation in Soils. John Wiley and Sons, INC., New York, USA.
  20. Magno, M.A., Luffman, I., Nandi, A. and Evanshen, B.G., 2018. Spatial interpolation of heavy metal concentrations in soils of bumpus cove. Available online at:
  21. Merian, E., Anke, M., Inant, M. and Stoeppler, M., 2004. Elements and Their Compounds in The Environment: Occurrence, Analysis and Biological Relevance. John Wiley & Sons, New York, USA.
  22. Mohammadi, J., 1999. Study of the spatial variability of soil salinity in Ramhormoz Area (Khuzestan) using Geostatistical Theory II. Cokriging. Journal of Science and Technology of Agriculture and Natural Resources. 3(1), 1-8.
  23. Moral, F.J., Terron, J.M. and Rebollo, F.J., 2011. Site-specific management zones based on the Rasch model and geostatistical techniques. Computers and Electronics in Agriculture. 75, 223-230.
  24. Mulla, D.J. and McBratney, A.B., 2002. Soil Spatial Variability. CRC press, LLC., USA.
  25. Murkers, J.P., 2000. Abiotic and Biotic Processes in the resource and control of selenium in the western phosphate resource area. MSc. Thesis, University of Idaho, Russia.
  26. Nazemi, L., Nazmara, S., Eshraghyan, M.R., Nasseri, S., Djafarian, K., Yunesian, M., Sereshti, H., Moameni, A. and Shahtaheri, S.J., 2012. Selenium status in soil, water and essential crops of Iran. Iranian Journal of Environmental Health Sciences and Engineering. 9, 1-8.
  27. Pais, I. and Jones, J.B., 2000. the Hand Book of Trace Elements. CRC press, USA.
  28. Reilly, C., 1996 Selenium in Food and Health. Springer, Germany.
  29. Renha, S.P., 1983. Clay dispersion in relation to changes in the electrolyte composition of dialyzed red- brown earths. Journal of Soil Science. 34, 723-732.
  30. Rhoades, J.D., 1996. Salinity: electrical conductivity and total dissolved solids. Methods of Soil Analysis: Part 3 Chemical Methods. 5, 417-435.
  31. Sadr, S. and Afyuni, M., 2017. Evaluation of vanadium distribution in agricultural and industrial land uses areas in Isfahan Region. Journal of Environmental Studies. 43(2), 207-218.
  32. Sadr, S., Afyuni, M. and Fathian Por, N., 2010. Spatial variability of arsenic under different land use in Isfahan Region. Journal of Science and Technology of Agriculture and Natural Resources. 13(50), 65-75.
  33. Sadr, S., Afyuni, M. and Movahedi Rad, Z., 2013. Mapping soil salinity in arid and semi-arid soil by Geostatistic Theory (Esfahan Province). Geography and Environmental Hazard. 2(7), 123-136.
  34. Thomas, G.W., 1996. Soil pH and soil acidity. Methods of Soil Analysis: Part 3, Chemical and Microbiological Properties. 5, 475-490.
  35. Walkley, A. and Black, I.A., 1934. An examination of the DEGTJAREFF method for determining soil organic matter and a proposed modification of the chromic acid titration method. Soil Science. 37, 29-38.