The new method for the forecasting hourly concentrations of air pollutants is presented in the paper. The method was developed for a site in urban residential area in city of Zagreb, Croatia, for four air pollutants (NO 2, O 3, CO and PM 10). Meteorological variables and concentrations of the respective pollutant were taken as predictors. A novel approach, based on families of univariate regression models, was employed in selecting the averaging intervals for input variables. For each variable and each averaging period between 1 and 97 h, a separate model was built. By inspecting values of the coefficient of correlation between measured and modelled concentrations, optimal averaging periods for each variable were selected. A new dataset for building the forecasting model was then calculated as temporal moving averages (running means) of former variables. A multi-layer perceptron type of neural networks is used as the forecasting model. Index of agreement, calculated for the entire dataset including the data for model building, ranged from 0.91 to 0.97 for the respective pollutants. As suggested by the analysis of the relative importance of the input variables, different agreements for different pollutants are likely due to different sources and production mechanisms of investigated pollutants. A comparison of the new method with more traditional method, which takes hourly averages of the forecast hour as input variables, showed similar or better performance. The model was developed for the purpose of public-health-oriented air quality forecasting, aiming to use a numerical weather forecast model for the prediction of the part of input data yet unknown at the forecasting time. It is to expect that longer term averages used as inputs in the proposed method will contribute to smaller input errors and the greater accuracy of the model.