In order to confront pollution events concerning the city of Elefsis, in the environmentally aggravated area of Thriassion Plain, an effort is undertaken to create a model forecasting maximal daily concentrations of NO(x) (NO(2)+NO), NO(2) and O(3). The data analyzed were obtained from the Bureau of Pollution Control and Environments Quality based in Elefsis. The model in question uses hourly values of the pollutants as well as meteorological data recorded at the center of the city of Elefsis from 1993 to 1999. Three fitting methods are utilized, namely ordinary least squares, piecewise, and quantile regression. The verification and reliability of the forecasting models are based on the measurements of the year 2000. The results are considered to be satisfactory, with the forecasted values following the general tendencies.