Estimation and forecasting of PM10 air pollution in Ankara via time series and harmonic regressions

In this study, monthly particulate matter (PM10) values in Ankara (39.9334° N, 32.8597° E) from January 1993 to December 2017 are examined. The PM10 are those thoracic particles whose aerodynamic diameter is less than 10 μm (micrometers), and it is of critical health importance due to the penetrability to the lower airways. As an alternative to classical unit root tests, a unit root test primarily based on periodograms is introduced owing to its advantages over alternatives. After examining the stationarity of the series through periodogram-based test as well as its standard rivals, periodic components in the series are examined and it is observed that the series has both periodic and seasonal components. These components are modeled, using the inherent dynamics of a time series alone, within a trigonometric harmonic regression setup, eventually yielding the forecast values for 2018 that turns out to be superior to those obtained by means of ARIMA (autoregressive integrated moving average). This is a striking result since the modeling framework requires no assumptions, no parameter estimations except for the variance of the white noise series, no simulations of the power of tests, no adjustments of test statistics with respect to sample size and no preliminary work as to independent variable which is simply time, i.e., the period of forecast.


Introduction
Air pollution can be defined as the phenomenon that various substances that need to be found in the air are out of the specified limits and the substances that should not be in the air are found at a hazardous rate to humans, plants, animals and environment (Cavkaytar et al. 2013). Along with various natural factors, air pollution increasing in parallel with human factors related to population growth, technological development and industrialization, has long been a global concern. Among all factors, increasing population density is one of the key human factors that cause air pollution without doubt. In that, the impacts of the natural pace of population increase were visibly augmented by the migration from rural to urban areas in Turkey especially since 1980. According to TurkStat, while Turkey's population was 43.9 million in 1980, this figure reached 63.2 million in 2000 and 80.8 million in 2017. Not surprisingly, the most dramatic change in Turkey's demographics has been the rate of urbanization. While the urban-rural divide was 24.2-75.8% in 1927, today this rate is reversed to 92.5-7.5% as of 2017.
A similar tendency is seen in metropolitan areas, as well. For instance, the population in Ankara reached from 3, 889,199 in 2000 to 4,771,716 in 2010 and 5,445,026 in 2017, where 94.3% of city population live in the centrum (TurkStat 2018). Such an increase in population and urbanization calls for a number of environmental issues. In Ankara, the second largest 1 3 city and capital of Turkey, air pollution has been a constant problem. The increase in the number of vehicles in traffic and their consumption of fossil fuels induced largely by population growth are its most important sources. Although the use of natural gas, removal of leaded gasoline consumption, emphasis on the use of natural gas-powered public transport vehicles initiated in 1988 slowed this increase, susceptibility of Ankara to air pollution in winters has not yet disappeared.
A brief look at the contaminants that threaten human health underlines carbon monoxide (CO), particulate matter (PM 2.5 and PM 10 ), sulfur dioxide (SO 2 ), nitrogen oxides (NO x ) and ozone (O 3 ) as noted by Li and Shi (2016). Particulate matter (PM) affects human health more than other pollutants (WHO 2018) and is defined as a mixture of solid particles and liquid droplets, hanging in the atmosphere. The PM 10 represents thoracic particles whose aerodynamic diameter is less than 10 μm and can penetrate to the lower airways. The PM 10 pollutants contain airborne particles due to the use of carbon-containing fuels, cigarette consumption and various industrial wastes, caused by humans, as well as particles formed due to the naturally occurring volcanic gases, seawater vapor, dust storms and forest fires which are mixed to atmosphere. The PM 10 causes many lung diseases, especially asthma, and cardiovascular diseases (Cavkaytar et al. 2013;Bayram and Dikensoy 2006). The World Health Organization (WHO) has set the annual average PM 10 amount standard as 20 μg/m 3 and the 24-hour average amount of PM 10 standard as 50 μg/m 3 in the 2005 Air Quality Directive WHO (2018). The European Parliament and Council has determined the annual average PM 10 limit as 40 μg/m 3 and the 24-hour average PM 10 limit as 50 μg/ m 3 under the directive 2008/50/EC of May 21, 2008, and that these values can be exceeded maximum 35 times per year European Commission (2018). In Turkey, according to Chamner of Environmental Engineers the amount of 24-hour average PM 10 limit is determined as 80 μg/m 3 in 2016 and this value has been updated as 60 μg/m 3 in 2018. According to the Air Pollution Report for 2017, the PM 10 concentrations recorded for Ankara in 2017 exceeded the upper bounds for more days than in 2016. As the Chamber of Environmental Engineers (2018) noted, such an incidence may be pointing at potential hazards to human health.
Note that, with regard to different particulate matters, PM 2.5 could also have been considered in our analysis. PM 2.5 are those particulate matters that have a diameter of less than 2.5 μm (micrometers), which is about 3% of the diameter of human hair. Since they are even smaller than their PM 10 counterparts, they are also called fine particles and have a higher likelihood to float longer in air. Nevertheless, our analysis could not consider PM 2.5 due to the inavailability of data on our side. Another important venue could be the consideration of the type of source of the PMs, i.e., a separate analysis of PMs from local sources (road traffic, industries, firewood burning, etc.,) and from distant sources (background PMs) could have shed more light to our understanding, basically since the variation in the concentration of PMs throughout the day and over the months of the year is very different and depends heavily on the source under study. However, we could not entail these in our analysis due to the inavailability of data with regard to sources of PM 10 .
Against this background, in this study, we propose an appropriate time series model for the monthly PM 10 amounts observed in Ankara from January 1993 to December 2017 and generate a set of forecasts for future. The monthly PM 10 measurements (μg/m 3 ) for Ankara from January 1993 to December 2017 are used. Considering the inherent features of seasonality and periodicity in data, we maintain a periodogram-based investigation. Our forecasting models are obtained through harmonic regression techniques which consider periodicity besides classical time series ingredients. Then, the results obtained for these two models are statistically compared and the forecast values obtained according to both models are compared with the observed values for the course of 2018. These comparisons reveal a great deal of forecasting superiority in our proposed estimation framework. What distinguishes the analysis of this study from the earlier work is that our proposed methodology makes use of a single time series and extracts the information embedded in the series by eloquently handling the repeating behavior, i.e., the seasonality and periodicity. So, even in the absence of other variables, which might be of concern under another analytical scheme, our work shows that a better forecasting performance has been viable.
In the remainder of the article, a brief review of the related literature is provided in Sect. 2 and our methodology is introduced in Sect. 3. Section 4 provides the implementation of method, prior to discussions in Sect. 5.

Literature
As mentioned earlier, people of today are exposed to air pollutants beyond standards, especially in areas with high population density, and their consequent hazards to health (Cavkaytar et al. 2013). When the literature on air pollution is examined, it can be seen that the studies are clustered around three main themes: (1) the studies intended to determine the relationship between pollutants and meteorological parameters and to determine the most important causes of air pollution, (2) the studies that demonstrate the relationship between pollutants and various health indicators and (3) the studies intended for forecasting various pollutants for the near future.
Among the studies intended to determine the relationship between pollutants and meteorological parameters and to determine the key causes of air pollution, Çiçek et al. (2004) use the stepwise regression to exhibit the relationship between SO 2 , PM 10 , NO, NO 2 , CO values and meteorological parameters such as the temperature, wind speed and relative humidity for the period of where ΔT is the daily temperature change, P is the daily mean atmospheric pressure, I d−1 is the air pollution index of the previous day and V is the daily average wind speed.
In order for such a regression model to be applied, the variables in the explanatory variable role must be non-random variables and the variables in the dependent variable role must be independent random variables. Also, since the observed values used are the values obtained in unit time intervals, it is more suitable to apply the time series models instead of the regression model. Nevertheless, it is still possible to use regression techniques, but firstly the data must meet requirements of the stationarity tests. Finally, Koutrakis et al. (2005) examine the relationships between PM 2.5 , PM 10 , PM 2.5-10 and meteorological variables using mixed regression models where they estimate the specific factors affecting the particle concentrations and their relative effects in Santiago, Chile, in 1989Chile, in -2001. Koutrakis et al. (2005) report significant relationships between meteorological variables and particulate matter amounts, like the reduction in the amount of particulate matter on Sundays as a result of reduced traffic and other polluting activities.
Silva (2015) is another important work in that the authors developed an index for air and noise quality via aggregation of relevant measurements and presented their results in comparison with the standardized legal limits for air pollution and noise. As an aid in decision making for urban planners and various policy makers, Silva's (2015) index is of high (1) value-added and widespread benefits. Owing to its intuitive methodology that involves the computation of a weighted linear combination of two base indices (with regard to noise and air separately), Silva (2015) provides a good description of the subject matter. The recent work by Ganguly et al. (2019) is also worthwhile with regard to measurement of air pollution along with its connection to other relevant factors. In that, the study provides not only an assessment of measurements at the urban stations in comparison with background monitoring facilities, but also the linkages of measured pollutant concentrations to seasonal factors and other pollutants' behavior. These two recent papers, in addition to those mentioned before, lay down a solid basis for understanding the problem at hand.
In the second strand, among the studies that investigate the relationship between pollutants and health indicators, to establish the relationship between the asthmatic response and the exposure of SO 2 and PM 10 , Berktaş and Bircan (2003) consider the number of patients who admitted to emergency room (Ankara Atatürk Chest Diseases and Chest Surgery Training and Research Hospital) with complaint of asthma between January 1, 1998, and December 31, 1998, meteorological conditions (daily mean rainfall, actual pressure, relative humidity, wind speed, duration and direction as well as minimum, average and maximum daily temperatures) and SO 2 and PM 10 concentrations in the Ankara region. Pearson and Spearman correlation tests and chi-square test were used to reveal the inherent statistical relationships. A significant correlation between the amount of PM 10 and the number of patients who applied was found, whereas a significant low level of association between the amount of SO 2 and the number of patients who applied to emergency service due to asthma was found. The results show that even short-term exposure to low-level air pollution in Ankara increases emergency room visits for suspected asthmatic reactions. Pope and Dockery (1992) study the acute health impact of respirable particulate pollution on symptomatic and asymptomatic children during 1990-1991 in the Utah Valley. Using logistic regression analysis, a positive association between PM 10 and respiratory symptoms was found. Ostro et al. (1999) reported a significant correlation between the PM 10 amount and the daily mortality in the Coachella Valley by means of Poisson regressions.
The last strand of the reviewed literature is devoted to forecasting of various pollutants for the near future. Turgut and Temiz (2015) apply the Box-Jenkins methodology to the weekly PM 10 concentration data obtained from the Ankara Sıhhiye station from January 1, 2010, to October 31, 2014, and forecast the future values of PM 10 . The study forecasts via an ARIMA (3,0,0) specification for November 2014, December 2014 and January 2015, yet omitting the seasonal effects which have been inherent to data. Another study for Ankara Province estimates an Air Quality Health Index comprising of concentration of pollutants such as PM 2.5 , O 3 and NO 2 (Bozkurt et al. 2015). Kurt et al. (2008) and Saral (2000) employ neural networks to predict air pollution in Istanbul, Kaplan et al. (2014) use Levenberg-Marquardt learning algorithm in artificial neural networks for PM 10 and SO 2 estimations for Kütahya Province and Yüksek et al. (2007) estimate the SO 2 pollution for Sivas centrum using a backpropagation artificial neural network framework.
All in all, it is salutary that the literature is deserted neither on air pollution nor on its determinants or consequences. Nevertheless, there is a common tendency to use regression rather than time series techniques to study PM 10 pollutant in Ankara, despite the genuine time series structure of the data sets. Indeed, this very time series structure might allow researchers to extract some rich information from data even in the absence of other explanatory variables or modeling peculiarities. Even when a pure times series approach is maintained, it is not rare to observe improper or no handling of seasonality. Among practitioners, it is not rare to observe the use of standard ADF 1 (augmented Dickey-Fuller) rather than DHF 2 (Dickey, Hasza and Fuller) or HEGY 3 (Hylleberg, Engle, Granger ve Yoo) tests even when seasonality is obvious. So, we try, in what follows to avoid these pitfalls through a genuine time series approach to modeling and forecasting.

Materials and methods
The fundamental aim in the time series analysis is to forecast the future values of a series using its observed (past) values. Stationarity 4 is the most important assumption for the series to be forecastable. If the series is non-stationary, the forecasts obtained and the statistical inference about the model parameters will not be meaningful. MA series are always stationary, but AR series may not be stationary. Moreover, most of the economic series are non-stationary series. In order to make forecast by the non-stationary time series, you need to provide the stationary with the help of various transformations. There are many methods in the literature for testing the stationary of time series. But the two tests come into prominence in terms of both practicality and applicability: Dickey-Fuller test based on the distribution of the least squares estimators of the parameters and the Phillips-Perron test that used the critical values of this distribution. For these methods, the values of test statistics and p values are directly calculated by many package programs. However, the DHF method which is developed by depending on the distribution of the symmetric least squares estimator or the HEGY method is used to test the stationary of the seasonal series. In either case, an auxiliary regression model is utilized. If there is not a suspect about the periodicity 5 of the series, one of the above-mentioned tests can be used for the stationarity test. If the series contains a periodic component, the stationarity test based on periodograms can be used. Even if the series do not contain periodicity it can use, at the same time it can apply to the seasonal series (Akdi and Dickey 1999).
Periodograms are usually used to reveal hidden periodicities found in the series (Fuller 1996;Wei 2006;Brockwell and Davis 1987). Akdi and Dickey (1998) proposed a test based on periodograms. General explanations about periodograms are given below.
Periodic functions suggest trigonometric functions. So, to test whether the series contains a periodic component or not, any time series Y 1 , Y 2 , … , Y n can be defined by where , R , w and are referred to as expected value, amplitude, frequency and phase, respectively. It is necessary to estimate these parameters. Furthermore, when w k = 2 k∕n, the Fourier frequencies are obtained. Due to the characteristics of cosine functions, for = R cos ( ) and = R sin ( ) , this model can be written as (4) Y t = + R cos (wt + ) + e t , t = 1, 2, … , n, 1 Test procedure which proposed by Said and Dickey (1985) for unit root test in ARIMA(p, 1, q) models. 2 Seasonal unit root test which proposed by Dickey, Hasza, and Fuller (1984) and based on the symmetric least squares estimator of the parameters. 3 Seasonal unit root test which proposed by Hylleberg et al. (1990) and considered quarterly data. 4 A time series Y t ∶ t ∈ T is a collection of observations at unit time intervals, where T is a set of indices. Any time series can be handled as moving average (MA), autoregressive (AR), autoregressive moving average (ARMA) or seasonal processes based on how it fluctuates around its own mean. Despite the resemblance of especially autoregressive time series with regression models, they do differ in the basic assumptions. Omission of this fact may often result in misleading results while studying time series data. Assume that a given time series Y t , t = 1, 2, … , n , or Y 1 , Y 2 , … , Y n , fit to the AR(p) representation as in Eq. (2): where e t is a white noise series with expected value 0 and variance 2 , that is, e t ∼ WN 0, 2 . The characteristic equation for Eq. (2) is written as in Eq. (3): (2) , then f is said to be "periodic" with a period of p . If f is periodic with p , all multiples of p are periods of f since . The smallest of these periods is the fundamental period (or briefly period) of the function.
If all the roots, in absolute value, of Equation (3) are less than 1, the specification of Eq. (2) is said to be stationary. The case where the root is larger than 1 is not so common in practice, and the series with one of the roots has the absolute value of 1 are unit-rooted (nonstationary) series.

Footnote 4 (continued)
According to this model, if the null hypothesis H 0 ∶ = = 0 is rejected, it is decided that the data contain the periodic component. Standard F test can be used to test this hypothesis. However, the use of the F statistic is not significant in case w k frequencies are not known (Wei 2006). According to this model, the least squares estimators of , and parameters are, respectively: where a k and b k are called as Fourier coefficients. Due to the characteristics of cosine functions, since Fourier frequencies are invariant according to the mean. By means of these Fourier coefficients, the periodogram ordinate at w k frequency of the time series is calculated as Time series are usually examined under time domain and frequency domain. While the autocorrelation function is the most important point in the time domain, the spectral density function is important in the frequency domain. There is also a transition between autocorrelation and spectral density function of series as stated in Herglotz theorem. If f w k is the spectral density function of the stationary time series, then the asymptotic distribution of the statistic I n w k ∕f w k becomes a chi-square distribution with 2 degree of freedom (exponential distribution with expected value of 1). That is, the probability density function of the asymptotic distribution of normalized periodograms is Hence, periodograms can be taken as an estimator of the spectral density function. The periodograms are also used for the investigation of probable periodicities in the data (for testing the hypothesis H 0 above). For any stationary time series, the periodogram values I n w k at each k frequency are calculated. The statistic V is defined by (5) Y t = + cos w k t + sin w k t + e t , t = 1, 2, … , n.
where I n w (1) is the greatest periodogram value and m = n∕2 . If there is no any periodic component in the data (under H 0 ∶ = = 0 ), then for the V statistic can be written (Wei 2006). For any selected level of significance, the critical value c is calculated by If V > c , then the hypothesis H 0 ∶ = = 0 is rejected and it is concluded that the series include periodic component.
If the given time series is stationary, it was stated that the asymptotic distribution of the normalized periodogram is a chi-square distribution with 2 degree of freedom. In that case, it is written as (Fuller 1996;Wei 2006;Brockwell and Davis 1987). Under the assumption that the series is not stationary in other words that it is unit-rooted, for each constant w k , it is where Z 1 and Z 2 are the independent variables which have the standard normal distribution and ̂2 n is an estimator of the variance of the error term (Akdi and Dickey 1998). Briefly, the asymptotic distribution is Again, it has been shown by Akdi and Dickey (1999) that the method is also applicable for the seasonal time series (Akdi and Dickey 1999). That is, the statistic T n w k can also be used to test the stationary of the seasonal series (whether it is unit root or not). Although the asymptotic distribution is valid for each constant w k , the frequency w 1 is usually used in the hypothesis tests. The critical values of the distribution are given by the authors.
The structure elaborated in Eq. (4) through (15) has an array of advantages for the modeler with regard to assumptions, distributions of test statistics and accuracy of numerical outcomes. In that, first, no model assumption is needed and the method is invariant to model specifications as the periodograms can be calculated without reference to model assumptions. Second, no parameter estimation is required except for the variance of the white noise series, as opposed to the standard unit root tests which require estimated model parameters first. Third, as the distribution of the statistic T n w k is known under H 0 and H a , the analytical power T n w k D → 2 1 + 3 2 1 , n → ∞.
function exists for the test. Fourth, the critical values of the test statistic do not depend on the sample size. Finally, since the periodograms are calculated through trigonometric transformations of data, any periodic components of data are to be captured by the method, a clear strength of the framework to yield more meaningful and accurate results eventually.

Analysis
As mentioned earlier, a chief capability of the methodology maintained in this paper is to reveal the information embedded in a time series by focusing solely on the time series itself, i.e., not explicitly referring to other variables. Such a capability on the side of modeling allows us to keep some factors outside the analysis. For instance, the prevailing meteorological conditions like wind direction, wind speed, etc., are not under our direct consideration. The bright side is that such an omission does not reduce the information content of our findings, since the impacts of these factors have already been established in the evolution of the time series over time. Equivalently, the numerous effects that seem to have been omitted are already captured as an array of periodic components in our periodogram-based analysis. Consequently, in this section the monthly PM 10 measurements (μg/ m 3 ) for Ankara from January 1993  Although there is a significant reduction in the PM 10 amounts especially after substitution of coals with natural gas in the late 1990s, the values are still above the European Union, World Health Organization and United Nations standards. When the monthly behavior of PM 10 values is considered (Table 1), it is seen that the averages for winter months are higher than those for the others. The highest values are observed in November, December and January. The lowest value is observed in June, July and August, though not falling below international standards.
To see whether the PM 10 values differ according to the months or not, the following one-way ANOVA (analysis of variance) model is considered: where e ij 's are the independent variables which have the normal distribution. As a result of the analysis (Table 2), the null hypothesis of H 0 ∶ 1 = 2 = ⋯ = 12 = 0 is rejected (16) y ij = + i + e ij , i = 1, 2, … , 12, j = 1, 2, … , 25, ( F value = 34, 33 , p value < 0.0001 ), i.e., the air pollution is said not to differ across months.

Time-domain approach
To determine a suitable time series model for the data, different candidates are considered and the model with the smallest AIC 6 (Akaike information criterion) statistic value (Table 3) is chosen.
Considering the results in Table 3, it can be said that the AR (13) is the most suitable model (Model I) for the data. Accordingly, the model considered is where e t ∼ WN 0, 2 . The parameter estimates for this model are given in Table 4.
According to the results obtained with both the PROC ARIMA and the PROC REG, the most suitable model (Model II) for the data is On the other hand, the value of the AIC statistic obtained for this model is very close to the value obtained for the AR(13) model. The parameter estimates are presented in Table 5. The results show that the model parameters are significant.

Stationarity
If the sum of the estimated values of the parameters in the model is close to 1, then this case indicates that the model is stationary. When the estimation results obtained for the Model II are considered (Table 5), it is seen that the results are different from each other such that ̂1 +̂1 2 +̂1 3 ≅ 0.9717 for PROC ARIMA and ̂1 +̂1 2 +̂1 3 ≅ 0.8361 for PROC REG. In order to specify exactly whether the assumption of stationarity holds, the unit root test has been performed. The results of the ADF and PP 7 (Phillips-Perron) unit root test with 13 the maximum delay length are presented in Table 6. The results show that the series is stationary.
After determining that the series is stationary, 12 monthly forecast values for 2018 year are calculated by taking into account the Model II considered appropriate for the data. The forecasted values calculated, standard errors of these forecasted values and 95% confidence limits are given 6 A criterion of model goodness which judges a model by how close its fitted values tend to be to the true values. 7 The unit root test which is proposed by Phillips and Perron (1988). Table 7. When the forecast values are examined, it is expected that the highest value will be observed in November. It is followed by December and February.
Due to the advantages offered by the unit root test based on the periodograms, the stationarity of the series is also tested with periodograms. For Model II, the periodogram value and estimate of variance are calculated as I n w 1 = 6133.66 and ̂2 n = 322.2444 , respectively. So, the value of the test statistic computed as (19) T n w 1 = 2 1 − cos w 1 * I n w 1 ∕̂2 n = 0.008349.
Following Akdi and Dickey (1998), the critical values of the test statistic are given as So, the null hypothesis that the series is not stationary has been rejected at the significant levels = 0.01 and = 0.05 ( T n w 1 = 0.008349 < 0.178 and T n w 1 = 0.008349 < 0.0348).

Periodicity
We then investigate whether there is hidden periodicity in the series by means of periodograms. Figure 2 displays the periodogram against frequencies from 0 to . Figure 2 suggests that one of the periodograms differs significantly from the others. This can be regarded as a sign of a possible periodicity in the data (Wei 2006). The largest 5-periodogram values obtained from the data, the frequencies that corresponded to those and their periods are given in Table 8. Also, the total value of the periodograms in all frequencies is calculated as The statistic V is and the critical values that corresponded to the statistic V , with the approximate value Wei (2006), are calculated as c = 1 − ( ∕m) 1∕(m−1) , where m is n∕2 and is defined by In order to look whether the other frequencies have the periodicity or not, the following formula is used: where I n w (i) for i = 1, 2, 3, 4, 5 are the periodogram values from big to small. Herefrom, V 1 = 0.518 , V 2 = 0.088 , If V > c , it will be rejected the null hypothesis that the dataset does not contain the periodic component and it will be decided that the series has a periodic component. According to the obtained results, since V = 0.5180 > c 0.01 = 0.0625 , it is concluded that there is the periodic component in the data. The period of the series is obtained as P = 2 ∕w k = 2 ∕0.52360 = 12 . This is an expected result because the data are monthly.
On the other hand, although it is This shows that there is not a periodicity at the frequency corresponding to the fifth largest periodogram value. Therefore, the periodic components corresponding to the first four frequencies must be considered.

The harmonic regression model
Considering that the model is stationary and the data contain the periodic component, it is assumed that the following regression model (Trigonometric Model I) is appropriate for the data: Estimates of model parameters are given in Table 9.
Y t = + A 1 cos 2 t 12 + B 1 sin 2 t 12 + A 2 cos 2 t 12.5 + B 2 sin 2 t 12.5 + A 3 cos 2 t 300 + B 3 sin 2 t 300 + A 4 cos 2 t 6 + B 4 sin 2 t 6 + e t , t = 1, 2, … , 300.  Based on the p values in Table 9, it is decided that the appropriate model (Trigonometric Model II) is The parameter estimates obtained according to this model are as given in Table 10.
Based on the Trigonometric Model II, the following prediction model is established for the monthly PM 10 data:  The graph of the PM 10 values observed and of the prediction values which were calculated according to this model is shown in Fig. 3a. Also, the graph of the PM 10 values observed and of the prediction values which were calculated according to the Model II is shown in Fig. 3b.
The monthly average PM 10 values observed for the last 23 years and the monthly average prediction values obtained from the ARIMA (Model II) and the harmonic regression (Trigonometric Model II) are presented in Fig. 3c. Note that the observations in the first 13 months are not taken into consideration in predictions, so the averages of the last 23 years are covered in the graphs.
As shown in Fig. 3c, the monthly average prediction values obtained by the harmonic regression are the closer to the real average values in comparison with the monthly average prediction values obtained by the ARIMA. The sum of the squares of the distances between these averages and the realized averages is calculated by where y i 's are the monthly realized averages, y a,i 's are the monthly average predictions obtained by the ARIMA and y h,i 's are the monthly average predictions obtained by the harmonic regression. Again, when the sums of squares are examined, it is concluded that the prediction values obtained by the harmonic regression are closer to the realized values in comparison with the prediction values obtained by ARIMA. This result is more apparent in Table 11, in which are present the realized values for the first two months of 2018 and the forecast values which were calculated by taking into account the prediction values obtained by both methods.
When Fig. 3a and b is examined, it is seen that the predicted values obtained for both models are very close to the observed values.
When the forecast values given above are examined, it is seen that the forecasts obtained from the both models are lower in the summer months (May, June, July and August). The result that the air pollution values observed in the summer months are lower than the other months is an expected result. The forecast values obtained from the harmonic model are much closer to the monthly average, whereas the forecast values obtained   with ARIMA are well above the monthly average. Accordingly, it can be said that the harmonic regression model is more consistent and reliable than the ARIMA. In fact, this result can be also seen in Fig. 3c. The monthly average values and the graph of prediction equation are presented in Fig. 4.

Results and discussion
The model estimates and its yielding forecasts we presented in Sect. 4 were said to be superior to those produced by rivalling frameworks as evidenced by the comparison of alternative forecasts with actuals. Equivalently, the analytical approach of this paper comes with sizable benefits for the analysts and policy makers. It is salutary that these benefits are further augmented by the reduced costs of modeling, i.e., modeling as a professional effort. Among those, absence of the need for model assumptions seems to be an appealing one for any researcher or analyst. The ability of one, especially of the field specialist, to compute periodograms without reference to model assumptions provides a certain improvement to reliability of projection practices of institutions. In a similar fashion, as opposed to the standard unit root tests that mandate the preliminary estimation of model parameters, the approach maintained here removes that need except for the variance of the white noise series. The third advantage, from a purely academic rather than field work perspective, is the knowledge of the distribution of the statistic T n w k under both H 0 and H a , implying the existence of the analytical power function for the test. Independence of the critical values of the test statistic from the sample size, in addition, provides another strength.
Finally, the periodograms, being calculated through trigonometric transformations of data, can practically capture any periodic components, so promising more meaningful and accurate results. Over a practical domain of work, analysts or policy makers can easily obtain forecasts via our framework 1 3 simply by substituting a particular value of t , the period of forecast, into the estimated function. In a nutshell, the harmonic regression approach as presented in this paper seems to be a source of ease and flexibility on the side of the field practitioners while providing the modelers with desired statistical properties.
With regard to the maintenance of a univariate modeling scheme, we again develop a benefit-cost perspective. Following the widespread principle of parsimony, also called the Occam's razor, among the models that yield similar explanatory or forecasting performances, the one(s) with lesser explanatory factors are to be preferred. This simple-looking yet strong principle ensures that a researcher could avoid redundant factors in a well-guided manner. In the current case as presented in this study, our pure time series approach based on periodograms departs from the knowledge of (1) embodiment of the periodic and non-periodic effects of all related factors in a time series of interest and (2) ability to express any time series as a sum of properly defined sinusoidal. As a matter of fact, these allowed us to produce quality forecasts even in the absence of explicit reference to potentially relevant factors like weather conditions or the behavior of other pollutants.