Modeling and Forecasting of Carbon Dioxide Emissions in Bangladesh Using Autoregressive Integrated Moving Average ( ARIMA ) Models

In the present paper, different Autoregressive Integrated Moving Average (ARIMA) models were developed to model the carbon dioxide emission by using time series data of forty-four years from 1972-2015. The performance of these developed models was assessed with the help of different selection measure criteria and the model having minimum value of these criteria considered as the best forecasting model. Based on findings, it has been observed that out of different ARIMA models, ARIMA (0, 2, 1) is the best fitted model in predicting the emission of carbon dioxide in Bangladesh. Using this best fitted model, the forecasted value of carbon dioxide emission in Bangladesh, for the year 2016, 2017 and 2018 as obtained from ARIMA (0, 2, 1) was obtained as 83.94657 Metric Tons, 89.90464 Metric Tons and 96.28557 Metric Tons respectively.


Introduction
Generally, emission of carbon dioxide (CO 2 ) from living animals, humans, wetlands, volcanoes, and other sources is nearly balanced by the same amount being removed from the atmosphere by plant photosynthesis and by the oceans.Human activity is disturbing this equilibrium by generating increased CO 2 from fossil fuels like as coal, gas, and petroleum products; and combustion via electricity generation, transportation, industry, and domestic use.The results of these imbalances are believed to be greenhouse effects: global warming, melting of polar ice sheets and caps, a rise in sea levels and subsequent coastal inunda-A.Rahman [1].IPCC (2007) [2] shows that carbon dioxide is the most dominant GHGs (greenhouse gases) which accounted 77% of the total global GHG emission where CH 4 , N 2 O and other gases contributed 14%, 8% and 1% respectively.
Climate change as a result of global warming has become one of the most important issues in the recent years.Reddy et al. (1995) shows that global mean temperature will also rise to 3˚C -4˚C with doubling of the CO 2 concentration [3].In 2007, the Intergovernmental Panel on Climate Change reported that there would be an estimated rise in the average global temperature between 1.1˚C and 6.4˚C within the next 100 years [2].Lau et al. (2009) shows that a meager 2˚C increase in temperature would greatly hamper many ecosystems and would cause an increase in the sea-level that would adversely impact the lives of people living in coastal zones [4].
The developed countries have a much higher share in global emissions than the developing ones.Nakicenovic (1994) [5] studied the perspectives of GHG emission on a regional basis.He stated that the developing countries are responsible for less than 16% of the CO 2 concentration due to their past consumption of fossil energy.The developed countries have a much higher share in global emissions than the developing ones.In Bangladesh, the power sector alone contributes 40% of the total CO 2 emissions [6] [7].Sarkar et al. (2015) has shown that the growth of CO 2 emission found 6.7% per year which is higher than the growth of GDP and energy consumption in Bangladesh [8].According to the third Assessment Report of IPCC (2007) [2], Bangladesh is the most vulnerable region of the world to climate change impacts [6] [9] [10] [11] [12].Global warming has caused fundamental changes to our climate.According to the IPCC (2001) [13], a 45cm sea-level rise could cause a potential land loss of 10.9% and a one meter sea-level rise a loss of 20.7%.Recently, the IPCC, in its Fifth Assessment Report, has suggested a 1.5 mm sea-level rise per year for Bangladesh.According to the IPCC (2007) Bangladesh will loose about 8% of its rice and 32% wheat production by the year 2050.The UNDP (2007) [14] predicts that 11% of the population will be directly threatened by a one meter sea-level rise.Rahman et al. ( 2015) [15] shows that in Bangladesh the maximum temperature increased significantly by 0.021 degree Celsius per year in Cox's Bazar and Sylhet region which is alarming for Bangladesh.

Materials
In the present study, time series secondary data on carbon dioxide emissions in Bangladesh were considered for the period 1972 to 2015 from Boden et al.

Methods
ARIMA is one of the most traditional methods of non-stationary time series analysis.In contrast to the regression models, the ARIMA model allows time series to be explained by its past or lagged values and stochastic error terms.The models developed by this approach are usually called ARIMA models because they use a combination of autoregressive (AR), integration (I) -referring to the reverse process of differencing to produce the forecast and moving average (MA) operations [19].
The ARIMA model is denoted by ARIMA ( ) , , , p d q where "p" stands for the order of the auto regressive process, "d" is the order of the data stationary and "q" is the order of the moving average process.The general form of the ARIMA ( ) , , p d q can be written as [20], (i.e.) MA (q) can be written as,

Result and Discussion
In Figure 1, from the autocorrelation (ACF) and partial autocorrelation (PACF), it is clear that there is no significant spike in the original series which also indicates that there are no significant effects of Auto-Regressive and Moving Average in the original series, that is, the carbon dioxide emission series is stationary without any difference.
After making the series stationary, different parametric combinations of ARIMA (p, d, q) model were tried to analyze the forty-four-year data (1972 to 2015) of carbon dioxide emission and the best fitted model is accepted on the basis of minimum value of all selection criteria as mentioned above in methodology.The results of performance of developed ARIMA (p, d, q) model is presented in Table 1.Table 1 exposed the performance of twenty ARIMA Models out of which ARIMA (0, 2, 1) was best out of all.ARIMA (0, 2, 2) and ARIMA (1, 2, 1) ranks second and third respectively while remaining ARIMA Models are not as good as these three.
Therefore, it was concluded that the appropriate model for forecasting the carbon dioxide emission in Bangladesh during 2015 was ARIMA (0, 2, 1) having minimum value of all selection criteria as compared to remaining nineteen models.
Table 2 shows the forecasted value of carbon dioxide emission using the best

Conclusion
This paper aimed to model the emission of carbon dioxide during 2015 in Bangladesh, by Autoregressive Integrated Moving Average (ARIMA) approach.On basis of results obtained, it is concluded that ARIMA (0, 2, 1) model having minimum value of all measures of selection criteria was found to be the appropriate model amongst all for predicting the carbon dioxide emission in Bangladesh.The model showed a good performance in case of explaining variability in Seasonal ARIMA model is to denoted by ARIMA (p, d, q) (P, D, Q), where P denotes the number of seasonal autoregressive components, Q denotes the number of seasonal moving average terms and D denotes the number of seasonal differences required to induce stationarity[21].The steps which are followed in order to define an ARIMA model as stated by Box & Jenkins: a) Identifying a model; b) Estimating the parameters of the model; c) Diagnostic checking.In the present paper, time series yearly data on emissions of carbon dioxide in Bangladesh were considered so there is no seasonal variation in the data which means non-seasonal ARIMA (p, d, q) models are applicable only.Comparison among family of different parametric combination of ARIMA (p, d, q) was done on the basis of minimum value of selection criteria which are Root mean squared error (RMSE), Mean percentage error (MPE), Mean absolute percentage error (MAPE), Mean absolute error (MAE), Maximum absolute percentage error (MAPE), Maximum absolute standard error (MASE) and Bayesian information criteria (BIC) [22] [23].
, M. M. Hasan DOI: 10.4236/ojs.2017.74038561 Open Journal of Statistics tions, and damage to agriculture and natural ecosystems, among others.There is an increasing trend of greenhouse gas (GHG) emissions worldwide due to human activities, which indicates a substantial increase in atmospheric concentrations of carbon dioxide (CO 2 ), methane (CH 4 ), nitrous oxide (N 2 O), hydrofluoro-

Table 1 .
Performances of different ARIMA (p, d, q) models of carbon dioxide in Bangladesh.