In the present study, we assessed the concentration of airborne HMs (Zn, Cu, Pb, and Cd) and their probable sources using the bark of Pinus eldarica as a bio-indicator. Hence, 47 tree bark samples were harvested according to the land uses and biomonitoring techniques in the city of Yazd, Iran. The potential health risks in 13 age groups, ecological risk, as well as the possible relationship between HM concentrations and traffic indicators, were evaluated. The order of average HM concentrations in the P. eldarica bark samples was as Zn > Pb > Cu > Cd. The mean values of non-carcinogenic risks of all HMs in entire age groups were within secure range (HQ < 1); however, the carcinogenic risk of Cd was higher than the allowed level (TCR > 1 × 10). About Pb, it was in the safe level. The main element causing potential ecological risks was Cd, indicating moderate to very high ecological risk in most of the study areas. There was an inverse significant association between distance from major roads and Pb concentration (β = -0.011 95% confidence interval (CI): 0.022, -0.0001). All HMs in bark samples render the negative Moran's index, representing a random spatial distribution pattern. Besides, according to principal component analysis (PCA), the first component accounted for 36.55% of the total variance, dominated by Cd, Pb, Cu, and Zn, respectively, and characterized by vehicle and industrial emissions. Our results infer that industrial activities and traffic are the main sources of HMs pollution in urban environments that should be considered by decision-makers.