In this study, three different sensors of satellites including the Moderate Resolution Imaging Spectroradiometer (MODIS), Multi-angle Imaging SpectroRadiometer (MISR), and Total Ozone Mapping Spectrometer (TOMS) were used to study spatial and temporal variations of aerosols over ten populated cities in Iran. Also, the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model was used for analyzing the origins of air masses and their trajectory in the area. An increasing trend in aerosol concentration was observed in the most studied cities in Iran during 1979–2016. The cities in the western part of Iran had the highest annual mean of aerosol concentration. The highest aerosol optical depth (AOD) value (0.76 ± 0.51) was recorded in May 2012 over Ahvaz, and the lowest value (0.035 ± 0.27) was recorded in December 2013 over Tabriz. After Ahvaz, the highest AOD value was found over Tehran (annual mean 0.11 ± 0.20). The results show that AOD increases with increasing industrial activities, but the increased frequency of aerosols due to land degradation and desertification is more powerful in Iran. The trajectory analysis by the HYSPLIT model showed that the air masses come from Egypt, Syria, and Lebanon and passed over the Iraq and then reached to Iran during summer. Aerosol radiative forcing (ARF) has been analyzed for Zanjan (Aerosol Robotic Network site) during 2010–2013. The ARF at surface and top of the atmosphere was found to be ranging from − 79 to − 10W m−2 (average − 33.45 W m−2) and from − 25 to 6 W m−2 (average − 12.80 W m−2), respectively.
