Model outputs were mapped using ArcView GIS 8.2 [19] on an "intelligent" grid based on 100 × 100 fixed nodes (rectangle of 7518 × 7618 m.) of the study area to develop concentration grids for each tracer pollutant. We constructed map layers by means of surface interpolations (simple kriging) and obtained the concentration maps for heavy metals and NO2.