Atmospheric Heavy Metals Biomonitoring Using a Local Pinus eldarica Tree

authors:

avatar Mohammad Miri 1 , 2 , avatar Mohammad Hassan Ehrampoush ORCID 1 , * , avatar Hamid Reza Ghaffari 3 , avatar Hamideh Ebrahimi Aval 2 , avatar Mehrab Rezai 1 , avatar Farzaneh Najafpour 1 , avatar Zeynab Abaszadeh Fathabadi 1 , avatar Mohsen Yazdani Aval 4 , avatar Asghar Ebrahimi 1

Environmental Science and Technology Research Center, Department of Environmental Health, School of Public Health, Shahid Sadoughi University of Medical Sciences, Yazd, IR Iran
Department of Environmental Health, School of Public Health, Sabzevar University of Medical Sciences, Sabzevar, IR Iran
Social Determinants in Health Promotion Research Center, Hormozgan University of Medical Sciences, Bandar Abbas, IR Iran
Faculty of Medical Sciences, Department of Occupational Health Engineering, Tarbiat Modares University, Tehran, IR Iran
Health Scope: Vol.6, issue 1; e39241
published online: August 9, 2016
article type: Research Article
received: May 20, 2016
revised: June 23, 2016
accepted: July 24, 2016
DOI: https://doi.org/10.17795/jhealthscope-39241

how to cite: Miri M, Ehrampoush M H, Reza Ghaffari H, Aval H E, Rezai M, et al. Atmospheric Heavy Metals Biomonitoring Using a Local Pinus eldarica Tree. Health Scope. 2017;6(1):e39241. https://doi.org/10.17795/jhealthscope-39241.

Abstract

Background:

Trees are one of the most commonly used bio-indicators and have been widely applied for heavy metals (HMs) monitoring in the ambient air because of their simple and efficient sampling and examination.

Objectives:

The main objective of this study was to evaluate local pine trees for their potential as a bio-indicator of HMs in urban areas. Besides, the distribution pattern of HMs in the study area, meteorological and land use factors affecting this pattern, and HMs source identification were investigated.

Methods:

Forty-eight sampling locations were selected in the Yazd city, capital of Yazd province in Iran, and concentrations of five HMs- iron (Fe), copper (Cu), zinc (Zn), lead (Pb), and cadmium (Cd)- were measured in the leaves of the most abundantly local pine trees, Pinus eldarica, using the atomic absorption spectroscopy method. The obtained data were analyzed for source identification using Spearman’s rank correlation coefficient and principal component analysis (PCA). The ordinary Kriging method was used to obtain a spatial map of HMs distribution in the study area.

Results:

Iron had the highest concentration (345.0053 mg/kg) in leaves when compared with other HMs. The concentration of Zn, Cu, Pb and Cd were 26.1384, 5.5905, 2.9963 and 0.7383 mg/kg, respectively. Spatial mapping indicated that industrial districts and areas with heavy traffic had the highest HMs pollution. Wind direction was recognized to be the most important factor in HMs transportation. Anthropogenic sources were known as the main responsible factor for HMs emission.

Conclusions:

Pine trees can be used as a worthy bio-indicator for HMs pollution monitoring in the ambient air in areas with extensive vegetation of pine trees.

1. Background

Air pollution is one of the main challenges in highly populated cities, which mainly originates from rapid and irregular urbanization and industrialization (1).

Heavy metals are one of the most important air pollutants facing humanity and other life forms in the ambient air. These are associated with some human health disorders and ecological problems. There are various anthropogenic sources and activities, which are responsible for releasing HMs into the atmosphere such as incomplete burning of multiple types of fuel in industrial activity, transportation and power plants. In contrast, transportation is recognized as the most important source of HMs in urban areas and has attracted increasing interest in this context (2-5).

Heavy metals are released into the urban atmosphere in the form of liquid and solid particulate matter (PM). The emitted PM has different sizes and is directly diffused into the ambient air (6). Moreover, some of the HMs is released in the fume phase related to the manufacturing processes and main fuel.

With respect to the above-mentioned points, it is important to measure and monitor HMs in the ambient air for control and management objectives. The trees that can be efficiently used as bio-monitors for atmospheric pollutants because of advantages such as simple sampling and examination are currently well established (7, 8). The other advantages of plants as a passive sampler for bio-monitoring aims are great spatial and temporal resolution due to better availability of vegetation and little sampling costs.

Several plant species, particularly evergreen trees, have been applied for atmospheric pollution monitoring and their benefits for this objective have been pointed out in the literature (9). In contrast to lower plants (fungi, algae, lichens or mosses), the trees that are not good indicators for bio-monitoring of atmospheric pollution are well-known (10-12). However, the main advantage of these trees, that has led to their selected as bio-indicators, is extensive dispersion in most cities as the main plant species. The other advantages of the tree are easy detection (13), and their long lifetime so that a replication of the examination is possible after a few years, while the time-trend distribution of pollutants can also be investigated (14).

In the recent years, there has been increasing interest in heavy metal monitoring in the atmospheric area using trees as bio-indicators.

Up to now, several trees have been used for HMs monitoring. The most popular trees used for this aim are Phoenix dactylifera (15, 16), Pinus sylvestris (17), Pinus pinea (18), Nerium oleander (19), Pittosporum tobira (20), Acacia retinodes, Eucalyptus torquata (21), Populus nigra (22), Tibouchina pulchra (23), and Cedrus libani (24). While some studies have been carried out on pine species as HMs bio-indicators, no studies have been reported for HMs monitoring in the urban atmosphere using Pinus eldarica specie.

2. Objectives

This study was therefore designed to assess the HMs concentration and their distribution pattern in the ambient air of Yazd city using Pinus eldarica as a local bio-indicator. For this aim, HMs concentration was measured in leaves of Pinus eldarica. In addition to the effect of meteorological conditions, traffic status and industrial activity on HMs pollution were investigated. Source identification was another objective of this study.

3. Methods

3.1. Study Area and Sampling Location

Yazd, the capital of Yazd province, is located in the center of Iran. The city has a total area of 2491 km2. It is located in the longitude of 54°36’ E, latitude of 31°88’ N, and altitude of about 1190 to 1300 meters above the sea level.

Yazd has dry periods during most months of the year and it is recognized as the driest city in Iran with an average annual rainfall of only 50 mm. On average, the warmest month is July with the maximum temperature of 33°C and the coolest month is January with the minimum temperature of 3°C. The annual mean temperature of the city is about 20.3°C. The annual mean relative humidity recorded in this city is 29%. The prevailing wind direction is west to north (25).

Forty-eight sampling locations were selected based on traffic status and land use in the study area. Selected locations were distributed uniformly throughout the city. Figure 1 shows the study area and sampling locations.

Study Area and Sampling Locations
Study Area and Sampling Locations
Figure 1.

Study Area and Sampling Locations

3.2. Bio-Indicator Selection, Sampling and Analyses

There are various tree species in Yazd city. However, one of the most abundant species is pine tree that can be seen throughout the city (26). Therefore, the Pinus eldarica, a species of pine tree, was selected as a bio-indicator to investigate HMs pollution of the urban atmosphere.

The trees in the main streets with 20 m distance from the traffic lights were selected for sampling. Samples were collected from different sides of same age trees (about 10 - 15 years old) at the heights of 1.5 to 2 m above ground level. The leaves with defects such as bird dropping, pesticide treatment, insect infestation, honeydew presence, chlorosis, necrosis, and unusual dust cover were avoided in the collection process. Further clean up steps were not conducted on taken leaves. The sampling process was conducted during late spring after a 10-day rainless period.

Sampling and analyzing were conducted in six steps, which can be listed as follows:

1- Sampling 50 g leaves from the tree in compliance with the conditions mentioned above.

2- Drying samples at laboratory temperature (about 25°C).

3- Extracting HMs based on acid digestion (27, 28), in which 1 g of dried leaves were added to 8 mL of concentrated HNO3 (Merck, 65%) and the mixture was left at laboratory temperature for one hour.

4- Placing the quartz tube containing mixture in the oven at 50°C for two hours and then at 160°C for four hours.

5- Filtering mixture through Whatman filters (0.48 µm) and diluting the final solutions to 50 mL with double de-ionized water.

6- Analyzing the final solution for HMs including zinc (Zn), copper (Cu), iron (Fe), lead (Pb) and cadmium (Cd) using graphite furnace varian spectra atomic absorption spectrometer (GF-AAS). For this aim, all standards and blanks were matrix-matched for the chemicals and reagents used in the sample preparation. For calibration of GF-AAS, standard solutions (0.001, 0.01, 0.1 and 1 mg/L for Zn, Cu, Pb, and Cd and 0.1, 1, 10, and 100 mg/L for Fe) were prepared from the stock solution by serial dilution. After this, the standard solutions were injected to GF-AAS through an auto sampler. Only the elements with given certified values and mean recoveries ranging from 80 to 120% were included in the data analysis. The detection limit for each element was considered as three times the standard deviation of eight replicate blank measurements.

3.3. Statistical Analysis and Spatial Mapping

The correlation between HMs concentrations in leaves were calculated using principal components analysis (PCA) and Spearman’s rank correlation coefficient.

In order to determine spatial distribution of heavy metal concentrations, the Kriging method in Arc GIS software version 10.1 was employed. Kriging is an advanced geostatistical method, by which the estimated surface is generated from a scattered set of points with z-values based on weighted moving average. The most important assumption with this method is that the distance or direction between sample points reflects a spatial correlation that can be used to explain variation in the surface. To determine the output value for each location, the Kriging tool fits a mathematical function to a specified number of points, or all points within a specified radius. The Kriging process involves the following steps: a) exploratory statistical analysis of the data, b) variogram modeling, c) creating the surface, and d) exploring a variance surface which is optional (29, 30). The estimator in the Kriging method is defined as follows (Equation 1):

Equation 1.

Where the Z × (xi), λi and z (xi) are the estimation grade, weight or sample quantity dependent on i, and measured variable value, respectively.

4. Results

4.1. Heavy Metals Concentrations

Table 1 shows the summary statistics for HMs concentration. As can be seen from the table, the highest and lowest average concentrations belonged to Fe (345 mg/kg) and Cd (0.73 mg/kg), respectively. The average value of 26.1 mg/kg for Zn, 5.59 for Cu, and 2.99 mg/kg for Pb were obtained. The maximum concentration of Zn (49.1 mg/kg) and Cd (4.85 mg/kg) were observed in the sampling point near the airport, industrial area and main bus terminal. However, the maximum concentration of Cu (161mg/kg), Pb (22 mg/kg) and Fe (717 mg/kg) were recorded in heavy traffic areas.

Table 1.

Descriptive Statistics of Measured Heavy Metals in Pine Leaves (mg/kg)

MinimumMaximumMeanStd. Deviation
Zn12.849.126.18.87
Cu0.041615.5922.9
Pb0.20222.993.19
Cd0.114.850.740.61
Fe15.4717345120

4.2. Source Identification

Table 2 shows the results obtained from the Spearman rank correlation analysis for HMs concentrations in different sampling locations. It is apparent from this table that correlations between HMs concentrations in sampling locations are very low (less than 0.475). These results suggest that sampling locations were impressed by different sources of contamination so that the tires and wheels of the automobiles are probably responsible for Zn, Cd and, Pb emissions while Fe and Cu are the main contaminants released from the corrosion process of automobiles and industry dumping areas.

Table 2.

Spearman Rank Correlations Between Measured Heavy Metal Concentrations in Pine Leaves (mg/kg)

ZnCuPbCdFe
Zn1.000
Cu0.475a1.000
Pb0.2040.1661.000
Cd0.331b0.417a0.1521.000
Fe0.2230.342b0.2600.337b1.000

aCorrelation is significant at the 0.01 level.

bCorrelation is significant at the 0.05 level.

The US environmental protection agency (EPA) has suggested the principal component analysis (PCA) technique for source identification of HMs pollution. In this study, PCA with the varimax rotation method was used for the apportionment of the probable HMs sources in pine leaves, according to previous studies (31, 32). Table 3 provides the results obtained from the PCA analysis including loading factors with varimax rotation and also eigenvalues.

Table 3.

Principal Component Analysis Results of Heavy Metals Concentrations in Pine Leaves

ElementRotated Component MatrixCommunalities
123
zn0.5120.6200.0750.652
pb0.847-0.026-0.1290.735
cu-0.1060.913-0.0250.845
cd-0.045-0.0090.9250.858
fe0.5380.0620.4250.474
Initial Eigenvalues1.2841.2221.059
% of Variance25.67124.43721.184
%Cumulative 25.67150.10871.292

4.3. Spatial Distribution

Spatial distribution of HMs concentration is shown in Figure 2. From the figure, the maximum concentrations of Zn were observed in the center and north of the city, which are areas with heavy traffic. The maximum concentration of Cd and Pb were detected in areas near the airport and heavy traffic locations, respectively. The eastern areas of the city had maximum concentration of Fe, and also maximum concentration of Cu was observed in the south of the study area.

Spatial Mapping of Heavy Metals Concentrations in Pine Leaves (mg/kg)
Spatial Mapping of Heavy Metals Concentrations in Pine Leaves (mg/kg)
Figure 2.

Spatial Mapping of Heavy Metals Concentrations in Pine Leaves (mg/kg)

The spatial map for the total concentrations of HMs is shown in Figure 3. The maximum concentrations of HMs were observed in southern areas of the city. A possible explanation for these results may be the wind direction pattern, which is often from north to south and moves pollutants from northern to southern areas. The other significant result indicated by the figure is that the central areas of the city with heavy traffic have higher concentration of HMs in contrast to those with low traffic.

Spatial Mapping of Total Heavy Metals Concentrations in Pine Leaves (mg/kg)
Spatial Mapping of Total Heavy Metals Concentrations in Pine Leaves (mg/kg)
Figure 3.

Spatial Mapping of Total Heavy Metals Concentrations in Pine Leaves (mg/kg)

5. Discussion

The literature has emphasized on the importance of trees as bio-indicators for air pollution monitoring. The HMs concentrations in the tissue of the tree display contamination levels of the soil and air. Soil and air pollution with HMs is a major concern, particularly in highly populated cities. To date, there has been a dramatic increase in bio-monitoring of air quality using trees for detection purposes as well as for monitoring the effects of the pollutants (33, 34). Also, trees can efficiently trap air-born particles and have been identified as major contributing agents for the declining of respirable PM, which is considered as the potential cause of major health complications for the exposed population.

The most important criteria for selecting suitable tree species as bio-indicator for HMs can be listed as follows: a) it should be able to differentiate between soil borne and airborne HMs, b) be simple and cheap to sample, c) be present in large numbers throughout the monitoring area, and d) be accessible in a wide geographical range. Pine trees due to the greater ability to survive in Yazd city represent all these criteria and our research completely supports this view (35, 36).

The higher Fe and Zn concentrations in the study area may be due to the local geochemistry and industrial activity. Several studies have also reported that Fe and Zn concentrations are much more than other heavy metals in leafy trees located in the urban area (12, 26).

The normal range of the Cu concentration in trees is 5 to 30 mg/kg (37) and the results of this study showed that measured Cu concentrations in all sampling points, except one, were in the normal range.

Plant roots can uptake only very small amounts of Pb (12). It has been described that Pb cannot be absorbed from the soil by the leaves of the tree since Pb is a HM with little plant motion (38). It can be concluded that measured Pb in the leaf of the trees is mainly airborne Pb. The toxic range of Pb is different in various literature reports. For example, according to Kabata-Pendias (2010), this range is 30 to 300 mg/kg (37), while Markert (1994) considered 3 to 20 mg/kg as a toxic range (39). The normal range of Pb is 5 to 10 mg/kg (37). The concentration of Pb in the study area was higher than the normal range. The most likely causes of high Pb concentrations in the study area are old automobiles driven with leaded petrol, the combustion of the diesel oil, and the manufacturing and metallurgical processes.

Ferrous-steel industry is the major source of environmental Cd contamination. Furthermore, mineral oils, vehicle wheels and waste sludge usage may emit Cd into the air and soil, which can increase the level of Cd in trees (40). A small amount of Cd also originates from synthetic fertilizers. The results of the present study showed high Cd concentrations in industrial areas. These results further support this idea.

Table 4 shows the HMs concentration obtained in this study and other studies in which leafy trees have been used as bio-indicators for investigation of HMs pollution in urban areas. The concentrations obtained in this study are inconsistent with data obtained in previous studies. A possible explanation for this may be the difference in studied areas standpoint of meteorological conditions, traffic status, industrial activity and consumed fuel.

Table 4.

Concentration of Heavy Metals in Leaves of Different Species From Different Studies

LocationTree SpeciesHeavy Metals, mg/kgReference
ZnCuPbCdFe
IranPinus eldarica26.15.592.990.74345Present study
IranPlatanus orientalis L3515.11.127-404Norouzi et al. (2015) (26)
Serbia and MontenegroAesculus hippocastanum9.952.961.07--Tomasevic et al. (2005) (41)
Turkey Robinia pseudo-acacia13920.872.61.33414Celik et al. (2005) (12)
SpainPopulus alba5427.635.003.82336Madejon et al. (2004) (42)
JapanR. pulchrum63.47.150.61.24-Suzuki et al. (2009) (43)

In order to identify the probable source of HMs, Spearman rank correlation coefficient and PCA were employed (Tables 2 and 3). The results of these analyses showed that there were no strong correlations between HMs. The observed correlations are due to the difference between emission sources of HMs. In PCA, factor 1 was dominated by Pb, accounting for 25.67% of the total variance. As generally identified, the main release source of Pb is the leaded fuel incineration (44). Therefore, this factor emphasizes on anthropogenic activities such as industrial process and transportation as a probable emission source and is the same as the results reported in previous researches (45, 46).

The second factor contributes to 24.437% of the total variance and has high loadings for Cu and Zn. Zinc may be produced by mechanical abrasion of automobile (47), the motor vehicle tires, and lubricating oils (48). Brake dust is previously known as an important carrier for Cu in PM composition (49, 50). Copper is commonly used in brakes to control heat transfer (51). Therefore, based on factor 2, the atmospheric concentrations of HMs are also from road traffic (50, 52).

In factor 3, Cd had a high loading of 0.925 and 21.174% of the total variance in contrast to other HMs. The concentration of Cd may potentially be related to coal combustion and traffic sources in the urban atmosphere (53). Therefore, the results of factor 3 show that in addition to transportation, industrial activity can be a probable source of HMs.

Spatial distributions showed that HMs concentrations in the center (heavy traffic area) and south of the city are more than northern areas, while a large number of industries exist in the north of the city. Therefore, wind direction is one of the major parameters in transporting PM and HMs from north (industrial area) to the south of the study area.

5.1. Conclusion

The leaves of pine trees were found to be a good bio-indicator for Zn, Cu, Pb, Cd and Fe contamination in the ambient air of the study area. The samples obtained from heavy traffic and industrial areas had maximum concentration of HMs due to greater human activity. This study showed that manufacturing and transportation are the main sources of HMs pollution. Also, wind direction was identified as the most important meteorological factor affecting HMs transport. More studies are required in order to use other local trees as bio-indictors for HMs and other air born pollutants in the studied and other regions.

Acknowledgements

References

  • 1.

    Kumar S. Appraisal of heavy metal concentration in selected vegetables exposed to different degrees of pollution in Agra, India. Environ Monit Assess. 2013;185(3):2683-90. [PubMed ID: 22773079]. https://doi.org/10.1007/s10661-012-2739-1.

  • 2.

    Mertens J, Luyssaert S, Verheyen K. Use and abuse of trace metal concentrations in plant tissue for biomonitoring and phytoextraction. Environ Pollut. 2005;138(1):1-4. [PubMed ID: 16023913]. https://doi.org/10.1016/j.envpol.2005.01.002.

  • 3.

    Naghipour D, Amouei A, Nazmara S. A Comparative Evaluation of Heavy Metals in the Different Breads in Iran: A Case Study of Rasht City. Health Scope. 2014;3(4):18175.

  • 4.

    Sayadi MH, Shabani M, Ahmadpour N. Pollution Index and Ecological Risk of Heavy Metals in the Surface Soils of Amir-Abad Area in Birjand City, Iran. Health Scope. 2015;4(1):21137.

  • 5.

    Oliva SR, Espinosa AJF. Monitoring of heavy metals in topsoils, atmospheric particles and plant leaves to identify possible contamination sources. Microchemical J. 2007;86(1):131-9. https://doi.org/10.1016/j.microc.2007.01.003.

  • 6.

    Anagnostatou VA. Assessment of heavy metals in central Athens and suburbs using plant material. Surrey: University of Surrey; 2008. [PubMed ID: 18523992]. https://doi.org/10.1002/pbc.21633.

  • 7.

    Markert B. Sample preparation (cleaning, drying, homogenization) for trace element analysis in plant matrices. Sci Total Env. 1995;176(1-3):45-61. https://doi.org/10.1016/0048-9697(95)04829-4.

  • 8.

    Aksoy A, Sahin U, Duman F. Robinia pseudo-acacia L. as a posssible biomonitor of heavy metal pollution in Kayseri. Turk J Botany. 2000;24(5):279-84.

  • 9.

    Murin A. Basic criteria for selection of plant bioindicators from regional flora for monitoring environmental pollution. Biologia (Slovakia). 1995;50:37-40.

  • 10.

    Sawidis T, Marnasidis A, Zachariadis G, Stratis J. A study of air pollution with heavy metals in Thessaloniki city (Greece) using trees as biological indicators. Arch Env Contam Toxicol. 1995;28(1):118-24. https://doi.org/10.1007/bf00213976.

  • 11.

    Sawidis T, Chettri MK, Papaioannou A, Zachariadis G, Stratis J. A study of metal distribution from lignite fuels using trees as biological monitors. Ecotoxicol Environ Saf. 2001;48(1):27-35. [PubMed ID: 11161674]. https://doi.org/10.1006/eesa.2000.2001.

  • 12.

    Celik A, Kartal AA, Akdogan A, Kaska Y. Determining the heavy metal pollution in Denizli (Turkey) by using Robinio pseudo-acacia L. Environ Int. 2005;31(1):105-12. [PubMed ID: 15607784]. https://doi.org/10.1016/j.envint.2004.07.004.

  • 13.

    Berlizov AN, Blum OB, Filby RH, Malyuk IA, Tryshyn VV. Testing applicability of black poplar (Populus nigra L.) bark to heavy metal air pollution monitoring in urban and industrial regions. Sci Total Environ. 2007;372(2-3):693-706. [PubMed ID: 17140640]. https://doi.org/10.1016/j.scitotenv.2006.10.029.

  • 14.

    Ballach HJ, Wittig R, Wulff S. Twenty-five years of biomonitoring lead in the Frankfurt/Main area. Environ Sci Pollut Res Int. 2002;9(2):136-42. [PubMed ID: 12008294].

  • 15.

    Aksoy A, Ozturk M. Phoenix dactylifera L. as a biomonitor of heavy metal pollution in Turkey. J Trace Microprobe T. 1996;14(3):605-14.

  • 16.

    Al-Shayeb SM, Al-Rajhi MA, Seaward MRD. The date palm (Phoenix dactylifera L.) as a biomonitor of lead and other elements in arid environments. Sci Total Environment. 1995;168(1):1-10. https://doi.org/10.1016/0048-9697(95)04556-g.

  • 17.

    Dmuchowski W, Bytnerowicz A. Monitoring environmental pollution in Poland by chemical analysis of Scots pine (Pinus sylvestris L.) needles. Environ Pollut. 1995;87(1):87-104. [PubMed ID: 15091612].

  • 18.

    Dongarra G, Varrica D, Sabatino G. Occurrence of platinum, palladium and gold in pine needles of Pinus pinea L. from the city of Palermo (Italy). Appl Geochemistry. 2003;18(1):109-16. https://doi.org/10.1016/s0883-2927(02)00051-3.

  • 19.

    Dongarra G, Sabatino G, Triscari M, Varrica D. The effects of anthropogenic particulate emissions on roadway dust and Nerium oleander leaves in Messina (Sicily, Italy). J Environ Monit. 2003;5(5):766-73. [PubMed ID: 14587847].

  • 20.

    Lorenzini G, Grassi C, Nali C, Petiti A, Loppi S, Tognotti L. Leaves of Pittosporum tobira as indicators of airborne trace element and PM 10 distribution in central Italy. Atmos Environ. 2006;40(22):4025-36.

  • 21.

    Sawidis T, Krystallidis P, Veros D, Chettri M. A study of air pollution with heavy metals in Athens city and Attica basin using evergreen trees as biological indicators. Biol Trace Elem Res. 2012;148(3):396-408. [PubMed ID: 22410948]. https://doi.org/10.1007/s12011-012-9378-9.

  • 22.

    Djingova R, Ivanova JU, Wagner G, Korhammer S, Markert B. Distribution of lanthanoids, Be, Bi, Ga, Te, Tl, Th and U on the territory of Bulgaria using Populus nigra 'Italica' as an indicator. Sci Total Environ. 2001;280(1-3):85-91. [PubMed ID: 11763275].

  • 23.

    Moraes RM, Delitti WB, Moraes JA. Gas exchange, growth, and chemical parameters in a native Atlantic forest tree species in polluted areas of Cubatao, Brazil. Ecotoxicol Environ Saf. 2003;54(3):339-45. [PubMed ID: 12651190].

  • 24.

    Onder S, Dursun S. Air borne heavy metal pollution of Cedrus libani (A. Rich.) in the city centre of Konya (Turkey). Atmos Environ. 2006;40(6):1122-33. https://doi.org/10.1016/j.atmosenv.2005.11.006.

  • 25.

    Meteorological Organization of Yazd province. Iran: Meteorological Organization of Yazd province; 2016. Available from: https://yazdmet.ir/SCphp?type=static&id=397.

  • 26.

    Norouzi S, Khademi H, Faz Cano A, Acosta JA. Using plane tree leaves for biomonitoring of dust borne heavy metals: A case study from Isfahan, Central Iran. Ecological Indicators. 2015;57:64-73. https://doi.org/10.1016/j.ecolind.2015.04.011.

  • 27.

    Sawidis T, Breuste J, Mitrovic M, Pavlovic P, Tsigaridas K. Trees as bioindicator of heavy metal pollution in three European cities. Environ Pollut. 2011;159(12):3560-70. [PubMed ID: 21907471]. https://doi.org/10.1016/j.envpol.2011.08.008.

  • 28.

    Turan D, Kocahakimoglu C, Kavcar P, Gaygisiz H, Atatanir L, Turgut C, et al. The use of olive tree (Olea europaea L.) leaves as a bioindicator for environmental pollution in the Province of Aydin, Turkey. Environ Sci Pollut Res Int. 2011;18(3):355-64. [PubMed ID: 20686863]. https://doi.org/10.1007/s11356-010-0378-z.

  • 29.

    Dai F, Zhou Q, Lv Z, Wang X, Liu G. Spatial prediction of soil organic matter content integrating artificial neural network and ordinary kriging in Tibetan Plateau. Ecological Indicators. 2014;45:184-94. https://doi.org/10.1016/j.ecolind.2014.04.003.

  • 30.

    Peng X, Wang K, Li Q. A new power mapping method based on ordinary kriging and determination of optimal detector location strategy. ‎Ann Nucl Energy. 2014;68:118-23. https://doi.org/10.1016/j.anucene.2014.01.002.

  • 31.

    Ragosta M, Caggiano R, Macchiato M, Sabia S, Trippetta S. Trace elements in daily collected aerosol: Level characterization and source identification in a four-year study. Atmos Res. 2008;89(1-2):206-17. https://doi.org/10.1016/j.atmosres.2008.01.009.

  • 32.

    Fang GC, Chang CN, Chu CC, Wu YS, Fu PP, Yang IL, et al. Characterization of particulate, metallic elements of TSP, PM(2.5) and PM(2.5-10) aerosols at a farm sampling site in Taiwan, Taichung. Sci Total Environ. 2003;308(1-3):157-66. [PubMed ID: 12738209]. https://doi.org/10.1016/S0048-9697(02)00648-4.

  • 33.

    Mingorance MD, Oliva SR. Heavy metals content in N. oleander leaves as urban pollution assessment. Environ Monit Assess. 2006;119(1-3):57-68. [PubMed ID: 16738782]. https://doi.org/10.1007/s10661-005-9004-9.

  • 34.

    Gajic G, Mitrovic M, Pavlovic P, Stevanovic B, Djurdjevic L, Kostic O. An assessment of the tolerance of Ligustrum ovalifolium Hassk. to traffic-generated Pb using physiological and biochemical markers. Ecotoxicol Environ Saf. 2009;72(4):1090-101. [PubMed ID: 19285727]. https://doi.org/10.1016/j.ecoenv.2009.01.010.

  • 35.

    Markert B. Plants as biomonitors: indicators for heavy metals in the terrestrial environment. New York: VCH; 1993.

  • 36.

    Wittig R. General aspects of biomonitoring heavy metals by plants. In: Market BA, editor. Plants as biomonitors. Germany: VCH; 1993. p. 3-27.

  • 37.

    Kabata-Pendias A. Trace elements in soils and plants. 4 ed. US: CRC Press; 2010.

  • 38.

    Gunthardt-Goerg MS, Vollenweider P. Linking stress with macroscopic and microscopic leaf response in trees: new diagnostic perspectives. Environ Pollut. 2007;147(3):467-88. [PubMed ID: 17050053]. https://doi.org/10.1016/j.envpol.2006.08.033.

  • 39.

    Markert B. Element concentration cadasters in ecosystems, Progress report on the element concentration cadaster project ECCE of INTECOL/IUBS. Paris: General Assembly of IUBS; 1994.

  • 40.

    Gunes A, Alpaslan M, Inal A. Plant nutrition and fertilizer. Agr Publication. 2004;(1539).

  • 41.

    Tomasevic M, Vukmirovic Z, Rajsic S, Tasic M, Stevanovic B. Characterization of trace metal particles deposited on some deciduous tree leaves in an urban area. Chemosphere. 2005;61(6):753-60. [PubMed ID: 15893796]. https://doi.org/10.1016/j.chemosphere.2005.03.077.

  • 42.

    Madejon P, Maranon T, Murillo JM, Robinson B. White poplar (Populus alba) as a biomonitor of trace elements in contaminated riparian forests. Environ Pollut. 2004;132(1):145-55. [PubMed ID: 15276282]. https://doi.org/10.1016/j.envpol.2004.03.015.

  • 43.

    Suzuki K, Yabuki T, Ono Y. Roadside Rhododendron pulchrum leaves as bioindicators of heavy metal pollution in traffic areas of Okayama, Japan. Environ Monit Assess. 2009;149(1-4):133-41. [PubMed ID: 18347924]. https://doi.org/10.1007/s10661-008-0188-7.

  • 44.

    Manno E, Varrica D, Dongarra G. Metal distribution in road dust samples collected in an urban area close to a petrochemical plant at Gela, Sicily. Atmos Envron. 2006;40(30):5929-41. https://doi.org/10.1016/j.atmosenv.2006.05.020.

  • 45.

    Kong S, Lu B, Bai Z, Zhao X, Chen L, Han B, et al. Potential threat of heavy metals in re-suspended dusts on building surfaces in oilfield city. Atmos Environ. 2011;45(25):4192-204. https://doi.org/10.1016/j.atmosenv.2011.05.011.

  • 46.

    Miri M, Allahabadi A, Ghaffari HR, Fathabadi ZA, Raisi Z, Rezai M, et al. Ecological risk assessment of heavy metal (HM) pollution in the ambient air using a new bio-indicator. Environ Sci Pollut Res Int. 2016;23(14):14210-20. [PubMed ID: 27053049]. https://doi.org/10.1007/s11356-016-6476-9.

  • 47.

    Jiries AG, Hussein HH, Halaseh Z. The quality of water and sediments of street runoff in Amman, Jordan. Hydro Proc. 2001;15(5):815-24. https://doi.org/10.1002/hyp.186.

  • 48.

    Soltani N, Keshavarzi B, Moore F, Tavakol T, Lahijanzadeh AR, Jaafarzadeh N, et al. Ecological and human health hazards of heavy metals and polycyclic aromatic hydrocarbons (PAHs) in road dust of Isfahan metropolis, Iran. Sci Total Environ. 2015;505:712-23. [PubMed ID: 25461074]. https://doi.org/10.1016/j.scitotenv.2014.09.097.

  • 49.

    Adachi K, Tainosho Y. Characterization of heavy metal particles embedded in tire dust. Environ Int. 2004;30(8):1009-17. [PubMed ID: 15337346]. https://doi.org/10.1016/j.envint.2004.04.004.

  • 50.

    Sternbeck J, Sjodin A, Andreasson K. Metal emissions from road traffic and the influence of resuspension—results from two tunnel studies. Atmos Environ. 2002;36(30):4735-44. https://doi.org/10.1016/s1352-2310(02)00561-7.

  • 51.

    Huang S, Tu J, Liu H, Hua M, Liao Q, Feng J, et al. Multivariate analysis of trace element concentrations in atmospheric deposition in the Yangtze River Delta, East China. Atmos Environ. 2009;43(36):5781-90. https://doi.org/10.1016/j.atmosenv.2009.07.055.

  • 52.

    Bem H, Gallorini M, Rizzio E, Krzeminska M. Comparative studies on the concentrations of some elements in the urban air particulate matter in Lodz City of Poland and in Milan, Italy. Environ Int. 2003;29(4):423-8. [PubMed ID: 12705939]. https://doi.org/10.1016/S0160-4120(02)00190-3.

  • 53.

    Ji YQ, Zhu T, Feng YC, Bai ZP. Concentration and Enrichment Characteristics of Elements in PM_ (10) of Tianjin [J]. Environ Sci Technol. 2006;7:19.