A Gaussian plume model was modified to simulate the dispersion of non-reactive air pollutants under non-homogeneous wind conditions through a multi-puff approach. It was applied to the city of Lisbon and evaluated by comparison with measured sulphur dioxide data, showing a reasonable skill to estimate the transport and dispersion of pollutants under complex wind field and different atmospheric conditions. The modelling results were integrated with observed data, based on correlation functions determined from historical values, to obtain the improved analytical results by using optimal interpolation. A significant improvement over the predictions by the Gaussian puff model alone was achieved.