Projection of Greenhouse Gas Emissions for the Road Transport Sector Based on Multivariate Regression and the Double Exponential Smoothing Model

: The economic and health impacts resulting from the greenhouse e ﬀ ect is a major concern in many countries. The transportation sector is one of the major contributors to greenhouse gas (GHG) emissions worldwide. Almost 15 percent of the global GHG and over 20 percent of energy-related CO 2 emissions are produced by the transportation sector. Quantifying GHG emissions from the road transport sector assists in assessing the existing vehicles’ energy consumptions and in proposing technological interventions for enhancing vehicle e ﬃ ciency and reducing energy-supply greenhouse gas intensity. This paper aims to develop a model for the projection of GHG emissions from the road transport sector. We consider the Vehicle-Kilometre by Mode (VKM) to Number of Transportation Vehicles (NTV) ratio for the six di ﬀ erent modes of transportation. These modes include motorcycles, passenger cars, tractors, single-unit trucks, buses and light trucks data from the North American Transportation Statistics (NATS) online database over a period of 22 years. We use multivariate regression and double exponential approaches to model the projection of GHG emissions. The results indicate that the VKM to NTV ratio for the di ﬀ erent transportation modes has a signiﬁcant e ﬀ ect on GHG emissions, with the coe ﬃ cient of determination adjusted R 2 and R 2 values of 89.46% and 91.8%, respectively. This shows that VKM and NTV are the main factors inﬂuencing GHG emission growth. The developed model is used to examine various scenarios for introducing plug-in hybrid electric vehicles and battery electric vehicles in the future. If there will be a switch to battery electric vehicles, a 62.2 % reduction in CO 2 emissions would occur. The results of this paper will be useful in developing appropriate planning, policies, and strategies to reduce GHG emissions from the road transport sector. 2 emission to 27.16%. Use of both battery electric and plug-in vehicles on a 50% basis will result in 54.3% reduction of CO 2 emissions. In the present study, the data set from one region (North America region) was used to test the model. Data sets from other regions or countries may be used in the future to test the robustness of the proposed method and reduce any geographical bias. Because the model was based on the results of historical data, it may not include new fuel-e ﬃ cient technologies or processes that may evolve in the future. This means the model may need to be re-tested with newer sets of data which could further improve the model. Having said that, the generality of the approach should be applicable in the forecast of GHG emissions.


Introduction
The greenhouse effect is a major concern in many countries. This is due to the high degrees of GHGs emissions released into the atmosphere [1]. One of the contributors to GHG emissions is the emissions from the transportation sector. The growing number of cars and trucks has triggered an incredible rise in GHG emissions. In 2007, 23% of the worldwide GHG emissions came from the transportation sector with 73% of it coming from the road transportation sector [2].
The transport services sector has witnessed a rapid rise in the past four decades, with the continuous increase in carbon dioxide CO 2 emissions from this sector becoming an important worldwide issue. According to the International Transport Forum [3], CO 2 emissions from the road transport sector passenger cars, buses and single-unit trucks. The proposed approach will enable researchers and planners to calculate the GHG generated from the road sector and will assist in developing appropriate measures towards travel demand management and introducing clean vehicle technologies. A review of the related literature reveals that there is a lack of study that has examined the GHG emissions from the road transport sector considering the VKT and NTV ratio for different modes of road transport.
The structure of the paper is as follows. Part 2 provides a literature review, followed by the research methodology in Part 3. Part 4 presents the results obtained from a multivariate linear regression model and double exponential smoothing models to predict GHG emissions. The final part presents recommendations for future research and conclusions.

Literature Review
Emission variables in the transportation field have been considered in several studies in the literature. These variables are generally influenced by the transport routes, distance travelled, age and size of the vehicles [11].
In the literature, various approaches have been proposed by the researchers to forecast GHG emissions. These approaches can be classified into five main categories: time series study, decomposition analysis, regression analysis, bottom-up method and system optimization tool.
Saboori et al. [12] applied time series analysis to examine the nexus between the emissions of CO 2 from the transport sector and the economic growth of OECD countries. Forecasting the CO 2 emission trends from the vehicles, vehicular energy intake and population to improve management of CO 2 emissions have been a recent focus in Taiwan [13]. Meyer et al. [14], present a demand analysis of passenger vehicles in 11 regions and the associated CO 2 emissions.
Likewise, Tokunaga and Konan [15] and Konur [16] used panel data to provide estimations of CO 2 from the transportation sector. Tolón-Becerra et al. [17] focus on the emission reduction targets of EU countries. A study found a relationship of a co-combination of fuel price (FP) and pay per capita on transport fuel consumption (FC) [18]. Also, the relationship between CO 2 discharges and vitality utilization have been highlighted in the literature [18][19][20][21][22]. Begum et al. [23] examined the implications of GDP, FC and population on the emissions of CO 2 . Several recommendations were put forward by Ivy-Yap and Bekhet [24] that have the potential to reduce CO 2 emissions such as the use of low-carbon technologies. Additionally, Talbi [25] used the Vector Autoregressive (VAR) model to investigate the impact of energy consumption, economic growth, urbanization and fuel rate on CO 2 emissions. The influencing factors for the changes in CO 2 emissions from the transportation sector were also further analyzed. Magazzino et al. [26] used the time series analysis and a Machine Learning approach to develop a relationship between municipal waste generation, GDP, and GHG emissions.
The growing CO 2 emissions can influence changes in modal shift, fuel mix, economic growth, transport energy intensity, and emission coefficients. Using a decomposition method, Timilsina and Shrestha [27] found that transportation energy intensity and economic growth are the major variables for CO 2 emissions. Furthermore, economic growth roles in the contribution of CO 2 emissions in ASEAN-5 and EU27 economies have been established [28,29].
Lakshmanan and Han [30] applied decomposition analysis and found that population and GDP were the main causes of CO 2 emissions in the transport sector. Further, Logarithmic Mean Division Index (LMDI) and Kaya Identity approaches were used by Li [31] to estimate the CO 2 emissions of urban freight, identify key parameters of emissions change, and to use those parameters to predict CO 2 emissions from the transport sector. Likewise, a multivariate generalized Fisher index (GFI) decomposition model based on the extension of Kaya Identity was built by Fan and Lei [32] to estimate the influence of the energy structure, energy intensity, per unit traffic turnover, transportation intensity, economic growth and population size on carbon emissions in the transportation sector of Beijing.
By using regression analysis, Sadorsky [33] investigated the connection involving GDP, salary, urbanization, and vitality demand and discovered that CO 2 emissions could be reduced by shifting to renewable energies. While high CO 2 emissions are accumulated in dense road network (e.g., cities with high population density), low CO 2 emissions are distributed in low population density areas with scattered road networks (e.g., rural areas) [34]. Alhindawi et al [35,36] identified the main factors of GHG emissions for the roads sector using the ratio between vehicle-kilometers and number of transportation vehicles for six transportation modes. Likewise, Xu et al. [37] found the influence of vitality structure, population, vitality intensity and GDP on CO 2 emissions. GDP was found out to be the main factor for CO 2 emissions, followed by the vitality structure and population. Friedrich and Trois [38], through a combination of regression analysis and life cycle assessment approach, estimated the GHG emissions for solid waste management. They used specific emission factors for various waste management processes and included upstream, downstream and operational GHG emissions. To study the relationship between road transport energy consumption, transport CO 2 emissions, and economic activity in Saudi Arabia, the autoregressive distributed lag (ARDL) approach was used by Alshehry and Belloumi [39]. The impact of GDP, energy intensity, carbon intensity, and total population on carbon dioxide emissions in China's transport industry was investigated by Lin and Benjamin [40]. They used the quantile regression analysis and found that rather than urbanization, it is the GDP, energy intensity, and carbon intensity that have a greater impact on carbon emissions. Danish et al. [41] applied autoregressive distributive lag (ARDL) and vector error correction model (VECM) to examine the relationship between transport energy consumption, economic growth and CO 2 emissions to account for the effect of Foreign Direct Investment (FDI) and urbanization. The results also showed a significant impact of transport energy consumption on CO 2 emissions from the transportation sector. Van der Zwaan et al. [42] conducted a study on CO 2 emissions of transport divisions and the decarbonization strategies in Europe. In Italy, strategies have been developed through COPERT III methodology in the analysis of road transport emissions [43]. This method was also used to study vehicles emissions in urban areas [44][45][46]. A multi-model, a time series and bottom-up approaches, were used by Paladugula et al. [47] to estimate the energy consumption and emissions from India's transportation sector. In the literature, many optimization models have been proposed to forecast CO 2 emissions and recommend solutions to reduce CO 2 emissions. Several studies have utilized linear programming models to investigate the benefits of energy choices that can assist in the reduction of CO 2 emissions by industries [48][49][50]. Furthermore, a mixed-integer model has been recommended by Hashim et al. [51] that further sheds light on the effects of fuel switching and balancing on power generation. Their study shows that shift to fuel efficiency is one of the most effective choices to reduce CO 2 emissions. This was supported by other research that used mixed-integer linear programming for waste to vitality to ensure costs reductions were related to CO 2 emissions and electricity generation [52]. Güzel and Alp [53] examined the impacts of the transport sector on climate change in Istanbul.
Generally, after a period of time, the road transport sector's GHG outflows demonstrate a pattern [54]. Therefore, through the use of statistical forecasting techniques, researchers and planners can anticipate future outflows. Brown [55] and Brown and Meyer [56] pioneered the double exponential smoothing (DES) method to predict the variety of a noise, and a trend for forecasting. This technique has since been advanced by Goodman [57], Gardner [58] and Gijbels et al. [59] for various applications. For example, Goodman [57] made a remaining examination to upgrade the appraisal of the exactness of DES models, while Gardner [58] developed a general exponential smoothing technique that considers regularity. Likewise, Gijbels, Pope and Wand [59] added knowledge to the existing exponential smoothing hypothesis by using a DES model inside a nonparametric regression structure.
There are various studies which have demonstrated the use of DES models in different fields which include ecological contamination. For instance, a DES model was used in South Korea to anticipate the arrangement of ozone and how it contributes to air pollution [60]. In England, a DES model was used to gauge power request and help in minimizing the regular effects which are associated with the utilization of power [61]. In China, a DES model was used in the creation of an indicator for the indication of defilement of stream water [62]. Furthermore, a DES model was used in the assessment of CO 2 contamination by urban traffic roads [63]. In the US, a DES model was used to look to assess the diminishing pattern of CO 2 emissions as well as performing forecasts [54].
This section provided a comprehensive review of the tools, methods and factors used for the forecast of GHG emissions. Tables 1 and 2 provide a summary of the methodology adopted in the existing literature to predict GHG emissions in various industry sectors. Table 1 shows the pros and cons of different forecasting approaches. Table 2 shows the summary of models used in the transport sector along with the regions where those studies have been conducted. As can be seen, linear regression analysis and DES models have been less examined, particularly for predicting GHG emissions from the road transportation sector. Therefore, by considering linear regression analysis and DES models, the present study will contribute towards the modelling of GHG emissions estimation and thereby identify critical factors that impact GHG emissions from the road transport sector. It is a powerful technique for prediction as it considers the seasoned patterns. Utilizes past behaviour for predictions. It assists to examine the current functioning of the series by utilizing past data. Assists in comparing the performance of two different series of a different type for the same duration.
It is expensive considering that predictions are dependent on the past data designs that are used to anticipate the future market action. Results may not be perfect. The many variables that altered the Variances of a series cannot be perfectly adjusted.
The different variables which affect the timeline may not be static for a longer period, and so the prediction may not be accurate.

Regression analysis
Simple to apply and understand The calculation process may be cumbersome.

Double exponential smoothing
It is easy to learn and implement. The accuracy level for forecasts is usually high It gives more significance to recent observations and weights them accordingly.
The spikes in the data aren't quite as detrimental to the forecast as previous methods.
It limits our ability to forecast demand using seasonality.

Model Development
Although there are various methods in the literature to forecast GHG emissions, there are differences in their scope and assumptions as they depend on data gathered from various countries. This may lead to variations in the results obtained. Although Alhindawi et al. [35] used a Multivariate Regression Model to forecast GHG emissions, its performance has not been compared with other similar methods such as the double exponential model. Further, the applicability of such forecasting tool in scenario analysis was not examined in Alhindawi et al. [35]. Therefore, in this study, the emphasis has been on a general robust approach which can be used in analysing and identifying various factors that influence the GHG emissions and developing a forecast of GHG emissions from the road transport sector. Further, we conduct scenario analysis by examining the effect of introducing plug-in hybrid electric vehicles and battery electric vehicles in the future in terms of greenhouse gas emissions reduction. The study process is comprised of two sections:
This model is constructed based on both NTV and VKT historical data which spans over a period of 20 years from the North America region.
To develop a good MRA model which can be used in GHG emissions forecasting, each variable needs to have a predicted value. In this paper, a forecasting tool which was based on the DES technique was used in generating these predicted values. The DES method has been highly recommended in the handling of historical data by Sullivan and Claycombe [64], Hyndman et al. [65] and is, therefore, a suitable tool in this study.

Vehicle Kilometres by Mode (VKM)
The efficiency of a vehicle fleet in terms of its activity or the volume of traffic is expressed as vehicle kilometers by mode (VKM). In transport planning, VKM is commonly used for estimating vehicle emissions. Also, VKM can inform decision making and policies such as infrastructure investment decisions and road safety policy. Therefore, accurate VKM estimations are critical [66].

Number of Transportation Vehicles (NTV)
This study includes various categories of vehicles including light trucks, tractors, combination trucks as well as passenger cars [66]. An increase of NTV will result in higher GHG emissions.

Data Sources
The researchers collected the historical data used in this study on Vehicle-kilometers by Mode (VKM) and Number of Transportation Vehicle (NTV) of road transport, 1990-2012, from the official data source of the North American Transportation Statistics (NATS) online database [67]. These data are utilized to build the regression model and to forecast the GHG emissions using the estimation from a double exponential smoothing model. Table 3 shows the complete set of the data used in this study.

Data Limitations
NTV data are based on information collected by the Department of Transportation in the U.S. The vehicle categories include passenger cars, light trucks and combination truck tractors. Data for light trucks include pick-up trucks, vans and sport utility vehicles. Combination truck tractors were categorized under "tractors". Passenger cars include taxis. Buses include local motor buses [67].
Passenger cars, light trucks and motorcycles are included in the VKM data. In 1997, the Department of Transportation in the U.S. updated the VKM data on the highway mode for several years. The main change included the relocation of some vehicles from the passenger car category to the light truck category [67].
The discontinuities in the data set as stated above may lead to some errors in the regression model's estimation. However, when a detailed good data set is available, the generality of the proposed method and the results obtained in this study could be further examined in the future.

Regression Analysis and Double Exponential Model
Regression analysis has been generally used in the examination of multifaceted information which is achieved through the creation of mathematical statements that shows how a response is related to a set of independent variables or predictors. GHG emissions are used as the response variable in this model. The VKM to NTV ratios for the six transportation modes are used as the causal variables in this study. These modes are depicted as follows; Light Trucks (LT), Passengers Cars (PC), Tractors (T), Motorcycles (C), Single Unit Truck (SUT) and Buses (B). It should be noted that there could be other several causal variables, but due to data unavailability, only the above parameters were considered in this study. Additionally, aggregated data analysis has been conducted by the researchers rather than carrying out a disaggregation of data into multiple modules of transport. In order to avoid the average squared residual which is associated with a recent one-step-ahead during forecasting, exponential smoothing is done by specifying a smoothing factor. According to Hyndman et al. [65] and Sullivan Sustainability 2020, 12, 9152 9 of 18 and Claycombe [66], the equation and formula which are used in double exponential forecasting are as shown below(Equations (1)-(6)): where: • F t+m is the forecast's representation after the performance of m of periods. • a t depicts the forecasted intercept while b t is the forecasted slope.

•
The intercept a t and the slope b t are approximated as shown in Equations (2)-(4): α is the smoothing constant which is used when weighing the past and current observations and it lies in this range; 0 < α < 1. When the value of α is less than one it is considered to have less of a smoothing effect, and this gives a greater weight the latest changes that have taken place in the data. When the value of α is closer to 0 it is concluded that it has a greater smoothing effect hence it becomes less responsive to changes that have been conducted recently. The selection has no specific formal procedure which is correct for DES and single values for time X t , and t respectively. The values S t and S t are computed using the formula shown below in Equations (5) and (6).
On the other hand, when α is higher it means that more weight has been given to the current observations which are most recent. Before one runs a model, one needs to choose α first. This is followed by a computation of forecasts which is achieved by the use of a variety of α values. The value which results into a small mean squared error during the calculations is concluded to have shown the expected future growth. Besides the selection of α value approximates, S t−1 and S t−1 values should be assumed with the values of t = 1 since there are no values that exist at that period. This is achieved after consideration of the equality of both values when compared with the original set of historical data [64,65].

Results and Discussions
This section brings out a presentation and discussion of the main results which have been obtained after carrying out an analysis of a multivariate regression of the GHG emissions model.

Regression Analysis Results
Minitab software was used in the analysis of this regression. It was also used together with the application of an Analysis of Variance (ANOVA) test to check the significance of the multivariate linear regression model.
Results which were obtained after conducting this analysis have been presented well in Table 4 below. All coefficients together with their expected signs have been well indicated. The negative sign shows that the predicted value of the dependent variable will be a value which is less than zero with a condition where the predictor values used are set to 0. This should always be the expectation in cases where the dependent variable has a negative mean value. It is worth noting that it is the overall and general relationship among variables which forms one of the most significant parts of the regression model. This section utilizes some of the aforementioned six modes of transportation which are abbreviated as passenger cars (PC), light trucks (LT), single unit trucks (SUT), buses (B), motorcycles (M) and tractors (T). After standardization of all coefficients, PC has been established to have the largest absolute value and then it is followed by LT, B and SUT respectively. Hence, PC is established as the most important variable in this study. To establish the importance of the variable in a given model, an estimation of the p-value is carried out. In this study, the p-values have been established and presented in the table as shown below. Consideration was given to the variables that had a p-value of less than 0.05. The output below shows predictor values for LT, B, PC and SUT since their p-values fall within the range of 0.0-0.05. The p-values for the Tractors and Motorcycles were 0.317 and 0.693 respectively. Hence, they were excluded from the model.

Model Adequacy Check
To ensure verification of the multivariate linear regression model, it is important to ensure that adequacy of the model has been checked. Regression analysis involves the use of an ANOVA tool that tests the importance and validity of the model. This is carried out based on some few assumptions. These include the residuals having an even distribution and a constant variance. These assumptions are validated using a graph which analyzes the residuals. To determine whether a small set of data are part of a given normal distribution, a plot showing normal probability is made. Since the straight line represents data points, it is concluded that the distribution is normal. This is well depicted in Figure 1a. The assumptions are checked further by plotting a residual vs fitted values as shown in Figure 1b. According to these plots, a constant variance has been associated with errors. This has led to the dispersion of the residuals being around zero. For instance, an increase or decrease of residuals in a pattern with the fitted values leads to errors which do not have a constant variance. This plot has points which seem to have a random dispersion around zero and this brings forward an idea suggesting the error has a mean of 0 which is reasonable. The fitted values do not appear to increase or decrease as a result of the vertical width of the scatter plot. Hence, the variance of the error is assumed to be constant [68]. Therefore, this analysis demonstrates satisfactory results which fall within the horizontal band. Secondly, in this model, there was no detection of influence or leverage points. Thirdly, multicollinearity was also not detected in this model. When large Variation Inflation Factors (VIFs), usually larger than 10, are detected, this means that there was a poor estimation of the regression coefficients. This is normally caused by near-linear and multicollinearity dependencies among the regression variables which leads to outcomes which are misleading. In this study, variables which have VIF that is less than 10 has been demonstrated in Table 4. These results can be translated to indicate the non-existence of multicollinearity in this model. Fourthly, the model points represent the behavior of its data adequately as the coefficient of multiple determinations (R 2 ), the adjusted R 2 and the predicted R 2 statistics are 91.8%, 89.46%, and 84.88% respectively.
Finally, the contribution of each vehicle category towards emissions of GHG is depicted by the last model (Equation (7)). Based on tests conducted previously, it can be concluded that the proposed model represents data accurately and has no violations towards the main assumptions. The equation for this model is as shown below in Equation (7). By application of Equation (7), there has been calculation of percentage error in order to draw a comparison between GHG real value and expected value. This has given a result value of 1.48% percentage error (Table 5).  By application of Equation (7), there has been calculation of percentage error in order to draw a comparison between GHG real value and expected value. This has given a result value of 1.48% percentage error (Table 5). Application of the DES model, which was described in Equation (1), can be used to calculate GHG emissions forecasted variables over a long period of time which could even be from 1990-2060. This is made possible by the use of smoothing constants (α's) as depicted in Table 6. By using a predicted period of 1990 to 2012, the resulting variables have been used in the calculation of percentage error of projection, which is depicted in Table 7. A summary of GHG emissions predictions has been made using the DES technique as shown in Figure 2. It has been predicted that by the year 2060, GHG emissions could rise to a figure around 2,957,400 K Metric Tons (Scenario 1). This increase is expected due to the increased modes of transportation. According to the Union of Concerned Scientist nationwide [69], GHG emissions for battery electric vehicles, the plug-in hybrids and gasoline vehicles are 144, 204 and 381 as depicted in Table 8. These figures are measured in terms of grams of CO 2 per mile. Based on this data, Table 9 presents five possible scenarios that could result from shifting from one type of vehicles to another as described below: • Gasoline-Plug-in vehicles-This would lead to about 46.5% reduction of national CO 2 emissions (Scenario 2).

Conclusions
GHG emissions from the road transport sector pose a significant challenge. To cut these emissions, planners and decision makers need to formulate appropriate policies and strategies to reduce emissions from the transport sector.
At present, models that forecast GHG emissions have some constraints. In particular, the unavailability of detailed data (e.g., fuel consumption per vehicle, average annual distance travelled) poses a constraint to the development of the forecasting tool. Further, the models are complicated with numerous equations to satisfy many assumptions. Also, the methodology and assumptions made are sometimes too rigid, which makes such models and equations difficult to be used in practice [9]. By use of the multivariate linear regression model and double exponential smoothing analysis,  Table 8. Average Emissions Nationwide [69].

Type of Vehicles GHG Emission (Grams of CO 2 Equivalents per Mile)
Gasoline vehicles 381 plug-in electric hybrids 204 battery electric vehicles 144

Conclusions
GHG emissions from the road transport sector pose a significant challenge. To cut these emissions, planners and decision makers need to formulate appropriate policies and strategies to reduce emissions from the transport sector.
At present, models that forecast GHG emissions have some constraints. In particular, the unavailability of detailed data (e.g., fuel consumption per vehicle, average annual distance travelled) poses a constraint to the development of the forecasting tool. Further, the models are complicated with numerous equations to satisfy many assumptions. Also, the methodology and assumptions made are sometimes too rigid, which makes such models and equations difficult to be used in practice [9]. By use of the multivariate linear regression model and double exponential smoothing analysis, this paper was able to examine GHG emissions factors which are associated with the road transportation sector. It made use of data collected in a span of 20 years (1990-2012) in the North America region. From the six modes of transportation, the VKM to NTV ratios were identified. The results indicate that the VKM to NTV ratio for the different transportation modes has a significant effect on GHG emissions, with the coefficient of determination adjusted R 2 and R 2 values of 89.46% and 91.8%, respectively. The model has 1.48% error which is associated with the difference between the actual value and the value which had been predicted. It is found that the most crucial regression model variable is the use of a personal vehicle.
This study made use of a DES model in order to carry out predictions of variables which affect the emission of GHG. This gave a 1.94% error between the real and the predicted value. From these results, it is predicted that by 2060, the emission of CO 2 will rise to a figure around 2,957,400 K metric tons. This will happen if we continue using gasoline vehicle types only. A switch to battery electric would lead to a potential 62.2% decrease in emission of CO 2 . If we shift from 50% use gasoline vehicles and 25% use of plug-in vehicles as well as 25% of battery vehicles, this will yield a decrease of CO 2 emission to 27.16%. Use of both battery electric and plug-in vehicles on a 50% basis will result in 54.3% reduction of CO 2 emissions. In the present study, the data set from one region (North America region) was used to test the model. Data sets from other regions or countries may be used in the future to test the robustness of the proposed method and reduce any geographical bias. Because the model was based on the results of historical data, it may not include new fuel-efficient technologies or processes that may evolve in the future. This means the model may need to be re-tested with newer sets of data which could further improve the model. Having said that, the generality of the approach should be applicable in the forecast of GHG emissions.
Many countries are aiming to reduce fuel/energy consumption and tailpipe GHG emissions of vehicles as road vehicles are a noticeable source of GHG emissions. For example, the European Commission set out a 60% reduction of transport emissions from 1990 levels by 2050. The GHG emissions forecasting tool as developed in this study may be used by planners and policymakers to estimate CO 2 emissions resulting from different modes of road transport and the vehicle kilometers travelled. Further, the estimation of CO 2 emissions can assist decisionmakers in developing a variety of policies or action plans to reduce the GHG emissions from the road transport sectors. For example, strict registration policies to promote the use of newer fuel-efficient vehicles, providing subsidies for fuel-efficient vehicles or implementing additional levies for older and polluting road vehicles could be implemented. This may help to reduce the greenhouse effect. In the future, researchers should conduct more analysis which is more detailed occurring at a disaggregated level and involving different and various transportation modes. While carrying out this study, a number of suggestions arose for further work on the topic of modelling GHG emissions from the transportation systems. First, an analysis which is more detailed occurring at a disaggregated level and involving different and various transportation modes including freights in fuel efficiencies needs to be conducted. Second, a more thorough analysis of the emissions factors for different kinds of travel and weekday-to-annual expansion factors may be insightful in terms of evaluation of GHG emissions in a community. Furthermore, it is recommended the five different scenarios as examined in this study are evaluated from a cost and economic point of view. Thirdly, more variables could be considered, such as the age of vehicles and new technology class of vehicles, that have a much lower carbon footprint as compared to older vehicles.
Author Contributions: This article is a joint work of the four authors. R.A, Y.A.N, A.K and N.S. R.A. conceptualization, collected the data, wrote the original draft and prepared the table. Y.A.N. analyzed the data and prepared figures. A.K. and N.S. reviewed the manuscript. All authors have read and agreed to the published version of the manuscript.
Funding: This research received no external funding.