Next Article in Journal
Ultrasound-Assisted Extraction of Antioxidants from Baccharis dracunculifolia and Green Propolis
Next Article in Special Issue
Non-Intrusive Load Monitoring Based on Deep Pairwise-Supervised Hashing to Detect Unidentified Appliances
Previous Article in Journal
Enhancing the Thermal Performance of Slender Packed Beds through Internal Heat Fins
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nexus between Energy Usability, Economic Indicators and Environmental Sustainability in Four ASEAN Countries: A Non-Linear Autoregressive Exogenous Neural Network Modelling Approach

by
Siti Indati Mustapa
1,
Freida Ozavize Ayodele
2,
Bamidele Victor Ayodele
1,* and
Norsyahida Mohammad
1
1
Institute of Energy Policy and Research, Universiti Tenaga Nasional, Jalan Ikram-Uniten, Kajang 43000, Malaysia
2
Department of Accounting and Finance, UCSI University Kuala Lumpur, 1 Jalan Menara Gading, Taman Connaught, Cheras Kuala Lumpur 56000, Malaysia
*
Author to whom correspondence should be addressed.
Processes 2020, 8(12), 1529; https://doi.org/10.3390/pr8121529
Submission received: 27 August 2020 / Revised: 1 October 2020 / Accepted: 3 October 2020 / Published: 25 November 2020

Abstract

:
This study investigates the use of a non-linear autoregressive exogenous neural network (NARX) model to investigate the nexus between energy usability, economic indicators, and carbon dioxide (CO2) emissions in four Association of South East Asian Nations (ASEAN), namely Malaysia, Thailand, Indonesia, and the Philippines. Optimized NARX model architectures of 5-29-1, 5-19-1, 5-17-1, 5-13-1 representing the input nodes, hidden neurons and the output units were obtained from the series of models configured. Based on the relationship between the input variables, CO2 emissions were predicted with a high correlation coefficient (R) > 0.9. and low mean square errors (MSE) of 3.92 × 10−21, 4.15 × 10−23, 2.02 × 10−19, 1.32 × 10−20 for Malaysia, Thailand, Indonesia, and the Philippines, respectively. Coal consumption has the highest level of influence on CO2 emissions in the four ASEAN countries based on the sensitivity analysis. These findings suggest that government policies in the four ASEAN countries should be more intensified on strategies to reduce CO2 emissions in relationship with the energy and economic indicators.

1. Introduction

In the past three decades, there has been rapid urbanization and industrialization among the Association of Southeast Asian Nations (ASEAN), which invariably has translated to economic growths, and high energy demand and consumption [1]. The increasing energy demand amongst the ASEAN countries is often meet with energy derived from fossil sources most especially coal, natural gas, and oil [1,2]. The share of coal in the energy mix of ASEAN has been on the rise and has been projected to continue to increase in a long time [3,4]. Recent energy outlook by the international energy agency revealed that the rising fuel demands among Southeast Asian countries have surpassed the regional production hence, the possibility of becoming a net importer of fuel [5]. Although, the high energy demand is seen as the powerhouse that drives industrial and economic activities within the region, and it often comes with the cost of high carbon dioxide (CO2) emissions [6]. As shown in Figure 1, the increase in economic indicators and energy demand also translated into energy-related CO2 emissions. To this effect, several energy economists have investigated the relationship between CO2 emissions, economic growth and energy consumption, as summarized in Table 1. Proper understanding of this relationship would help in finding lasting solutions to the age-long dilemma of sustainable economic growth and environmental pollution.
In the last one-decade various modelling techniques have been employed to investigate the causal relationship between CO2 emissions, economic indicators, and energy consumption in various countries [7,8,9,10]. Modelling techniques, such as panel quantile regression approach, Simultaneous Equation Modelling, system of simultaneous equations using seemingly unrelated regression, and Cointegration-Vector Error Correction Modelling have been widely applied. Zhu et al. [11] and Heidari et al. [12] employed panel quantile regression model to investigate the effect of foreign direct investment (FDI), economic growth, and energy consumption on CO2 emissions in Malaysia, Philippines, Indonesia, Malaysia, Singapore, and Malaysia. The study revealed that FDI, economic growth, and energy consumption have a heterogeneous effect on CO2 emissions. However, the level of importance of the input variables, as well as the prediction of the CO2 emissions, based on the variables, could not be ascertained using the panel quantile regression model. In a recent study, Salman et al. [13] also employed panel quantile regression for modeling the impact of export and import on carbon emissions in seven ASEAN countries namely Brunei, Indonesia, Malaysia, Philippines, Singapore, Thailand, and Vietnam. Using datasets between 1990–2017, the study revealed that CO2 emissions in the seven ASEAN countries were significantly influenced by export and import. Besides export and import, the authors reported that CO2 emissions increased by population size and energy intensity. However, the extent of influence and the level of importance of these variables were not reported. Also, the technique could not predict CO2 emissions based on the non-linear relationship between the variables.
Beside panel quantile regression, other modelling techniques that have been employed in investigating relationship between CO2 emissions, economic indicator, and energy consumption include Cointegration and Vector Error Correction Model [14], Auto-Regressive Integrated Moving Average (ARIMA) and Simple Exponential Smoothing Models [15], Autoregressive Distributive Lag (ARDL) model [16], Linear regression [17], Backpropagation neural network [17], non-linear dynamic neural network [17], and Simultaneous Equation Models [18]. These techniques were effective in explaining the causal relationship between CO2 emissions and the various input variables investigated. However, one major shortcoming of these techniques is their inability to accurately predict CO2 emissions based on the various input variables and their level of importance on the prediction, due to the non-linear relationship that exists between them. This shortcoming can be overcome using the nonlinear autoregressive networks with exogenous input (NARX) [17]. The NARX model is a time series recurrent dynamic feedback neural networks that have the capability to model the interrelation between input and output variable for prediction purposes [19]. Moreover, the level of importance of these variables could also be determined using the vector weights that connect the input variables with the output [20]. Xu et al. [17] employed NARX technique to model the prediction of CO2 emissions peak in China. Using datasets from 1978–2016, the study revealed that the NARX model was robust in predicting the CO2 emissions in China with GDP having the highest level of influence based on the sensitivity analysis. To the best of the authors’ knowledge, the use of NARX model for predictive modeling and to investigate the nonlinear relationship between economic indicators, energy consumptions and CO2 emissions in ASEAN countries such as Malaysia, Thailand, Indonesia and the Philippines have not reported in the literature. This study is aimed at employing the NARX techniques for modelling the prediction of CO2 emissions in Malaysia, Thailand, Indonesia and Philippines based on the nonlinear relationship that exist between them and to determine the level of importance of the prediction.

2. Materials and Methods

2.1. Data Description

The variables used in this study for the NARX modelling were carefully selected based on previous literature summarized in Table 1. The CO2 emissions as a function of kg per $US of GDP, energy consumption per capita measured in kWh, GDP per capita (current $US), population, oil consumption, coal consumption, were annual data from 1965–2018 obtained from World Bank Open database and British Petroleum energy database. The CO2 emission (kg per 2010 $US of GDP) represents the CO2 emissions from industrial and human activities in the utilization of fossil fuel. The energy consumption per capita, which is one of the fundamental markers of economic development of any country, explains energy utilization by an individual per year. The causal relationship between energy consumption, CO2 emissions and economic growth has been investigated by several authors using various modelling techniques as shown in Table 1. The annual coal consumption is the amount of coal required for electricity generation. The rapid industrialization of the four ASIAN countries requires high energy demand which heavily dependent on coal. Liu et al. [27] employed coal consumption as one of the variables in the scenario analysis of energy consumption and CO2 emission in China using system dynamic model. Similarly, Chandran Govindaraju and Tang [28] also investigated the dynamic link between CO2 emissions, coal consumption and economic growth in China and India. Oil consumption is the amount of crude oil utilized for refining purposes. There has been a steady increase in oil consumption in the four ASEAN countries over the years. The GDP per capita is an indicator of a country’s economic growth as a function of the population. Over the years, the GDP per capita of the four ASEAN countries has increased steadily. This could be attributed to the rapid industrialization and economic growth of the region.

2.2. The NARX Neural Network Modelling

Artificial Neural Networks are biological system inspired learning models comprised of inbuilt algorithms with a number of interlinked neurons between the input and the output [29]. The network usually receives input signals which are transformed to output signals. The NARX is a time series feedback artificial neural networks represented in Equation (1) [30]. It comprises interconnected feedback synaptic and delays that provides a flow of signals between the neurons [31]. The major advantage of the NARX as a robust time-series modeling technique is the significance of the delay that supply the precise historical information of a set of data at the current moment and the feedback loop that help to filter the historical data thereby allowing for the prediction of more accurate output [32]. The model architecture of the NARX neural network is depicted in Figure 2. It consists of the input, hidden and output layers connected by artificial neurons. The artificial neurons execute certain actions on the input signals as shown in Equation (2). In Equation (2), the hidden nodes are estimated in terms of the weights of the paths associated with the input nodes and the bias path which is a component of the back-propagation algorithm [33]. The configuration makes use of the past and present values of u(t) as well as the actual past values of y(t) for the prediction of y(t + 1). The multilayer perceptron topology which is a feedforward backpropagation network is often adopted in the training of the NARX [34]. The training entails the adjustment of the weights in such as a manner that the errors between the actual and the predicted values are minimized [35]. These errors are usually measured as mean square error (MSE) defined in Equation (2) [36]. The robustness of the NARX predictive models is measured using the coefficient of correlation (R) represented in Equation (3). The model development was performed using the Neural Network toolbox in MATLAB version 2019 a (MathWorks Inc.),
y ( t ) = f ( y ( t 1 ) , y ( t 2 ) , , y ( t n y ) ,   u ( t 2 ) , , u ( t d u ) )
where y ( t ) = depicts the output of the NARX model and u ( t ) depicts the input of the NARX model.
M S E =   1 n i i n ( y i y ^ i ) 2
R =   n   y i y ^ i (   y i ) (   y ^ i ) ( n   y i 2 (   y i ) 2 ) ( ( n   y ^ i 2 (   y ^ i ) 2 )
where n, y i , . y ^ i , are the number of observed variables, the observed outputs, and the predicted output, respectively.
In Equation (1) each of the successive output signals ( y ( t ) ) is regressed on the preceding value of the ( y ( t ) ) and the values of an exogenous input signal as shown in Equation (4).
ϑ =     k w k x k + b
where ϑ , w k , x k , b , θ are the neurons’ output, the weight factor, k t h input, the activation function and the bias. The sigmoid activation function represented in Equation (5) is often employed in predictive modeling due to it smoothen effect:
θ =   1 1 + e x p ϑ

2.3. Sensitivity Analysis of the Input Variables

The sensitivity analysis to determine the level importance ( ρ i z ) . of the input variables on the CO2 emission per capita was determined using modified Garson’s algorithm represented in Equation (6) [33,37]. The algorithm involves the partitioning of the connected weights of the hidden-output neurons within the NARX network into elements associated with each input neurons. The relative predictive importance of the input variables is revealed by the connection weights in the paths of each of the input nodes the corresponding output nodes [33],
ρ i z =   k = 1 i | ω i k v j z | / x = 1 i | ω x j | i = 1 n k = 1 i | ω i k v j z | / x = 1 i | ω x j |
where ω i k , depicts the weight between i th input and j th hidden unit; v j z signifies the weight between the j th hidden unit and the z th output while x is the number of hidden neurons, ω i k is the weight of the xth input and jth hidden unit, n is the number of input neurons.

3. Results and Discussion

3.1. Optimization of the Hidden Neurons

Optimization of the hidden neurons was performed to determine the number of hidden neurons that will minimize the MSE in the NARX neural network architecture. The hidden neurons were varied from 1 to 29 using different delays ranges from 1 to 15 (See Table A1, Table A2, Table A3 and Table A4 in the Appendix A for detail of the optimization). Table 2, Table 3, Table 4 and Table 5 shows the optimum neuron obtained for each of the countries investigated. In Table 2, optimized hidden neurons of 29 with MSE values of 3.92 × 1021 was obtained for NARX neural network architecture using 13 delays. This resulted in a NARX model architecture of 5-29-1 representing the input layer, hidden neurons, and the output layer. This optimized NARX model architecture was subsequently obtained for modelling the prediction of CO2 in Malaysia. Table 3 shows the variation of the hidden neurons for each iteration of the delay using the dataset for Thailand. Optimized hidden neurons of 19 with MSE value of 4.15 × 10−23 were obtained for the NARX neural network architecture using 13 delays resulting in an optimized network with the architecture of 5-19-1 denoting the input, hidden, and output layers, respectively. The subsequent modelling of the CO2 prediction in Thailand was modelled using the optimized NARX model architecture. Table 4 and Table 5 show the optimized neurons obtained using the datasets for Indonesia and the Philippines. An optimized 17 and 13 hidden neurons were obtained for the NARX neural network used for training the datasets for Indonesia, and the Philippines, respectively. This was achieved using 15 delays each for the optimized hidden neuron. Based on the optimized hidden neurons, NARX neural network model architectures of 5-17-1, and 5-13-1 were obtained for modeling the prediction of CO2 emissions in Indonesia, and the Philippines, respectively. Besides the low MSE obtained, the suitability of the four optimized NARX models as robust predictive techniques were further ascertained from the high values of the R and R2 which are greater than 0.9 as shown in Table 2, Table 3, Table 4 and Table 5. The optimized architecture of the NARX models used for the predictive modelling of the CO2 emissions in Malaysia, Thailand, Indonesia and the Philippines is generally represented in Figure 3 while the detailed architectures are depicted in Table 6. The optimized NARX models consist of the five input nodes, the hidden layers (which is made up of the optimized hidden neurons and delays) and the output node. The validation of each of the optimized NARX model is depicted in Figure 4. The best validation performance with MSE values of 6.93 × 10−4, 9.26 × 10−4, 2.03 × 10−3, and 7.99 × 10−4 was obtained at the epoch of 6, 2, 7, and 11, respectively for Malaysia, Thailand, Indonesia and Philippines (Figure 4a–d). High epoch values indicate that it takes longer iterations for the best performance to be validated. Zounemat-Kermani et al. [38] have reported optimized NARX model architectures of 10-12-1, 10-10-1, 8-9-1, and 6-8-1 for modelling the prediction of gaseous emissions within the influent chamber of four different wastewater treatment plant. The variation in the NARX model architectures could be attributed to the nature of the input variables, the datasets and the non-linear relationship between the input and the output variables.

3.2. Performance Evaluation of the Optimized NARX Model

The performance of the optimized NARX models based on their predictability is depicted in Figure 5. It is obvious from the dispersion plots in Figure 5a–d that the actual values of the CO2 emissions from Malaysia, Thailand, Indonesia, and the Philippines are in proximity with those predicted by the NARX models. The actual values of the CO2 emissions are consistent with the predicted with high R2 > 0.9. This can further be confirmed by the residuals shown in Figure 6. The residual range of −0.0005 to +0.0094 obtained for the predicted CO2 emissions in the four ASEAN countries signifies that the NARX models possessed the capability to predict with minimal residuals. Besides very low MSE values (Table 2, Table 3, Table 4 and Table 5) obtained for the optimized NARX model architecture is evidence of its robustness as a predictive modelling tool. The training, validation and testing of the datasets using the optimized NARX models are depicted in Figure 7. In each stage of the training, validation, and testing, the actual and the predicted CO2 emissions of the four ASEAN countries understudied were obtained with minimized errors. The autocorrelation analyses of the NARX models used for the prediction of CO2 emissions in the four ASEAN countries investigated are in Figure 8a–d. Most of the prediction errors are within the 95% confidence limits, which is an indication of the reliability of the estimations of the weight, bias and the network parameters [39]. The robustness of the NARX models employed in this study for the prediction of CO2 emissions in four ASEAN countries is consistent with that reported by several authors reported in Literature and summarized in Table 7. Alcan et al. [40] for the prediction of NOx emissions from diesel engines. Using eight different input variables, the authors employed an optimized NARX model architecture predictive modelling of NOx emissions. The findings revealed that the optimized NARX model architecture predicted the NOx emissions with a high degree of accuracy. Similarly, the work of Xu et al. [17] also confirms the applicability of NARX as a predictive modeling tool. The authors employed an optimized NARX model for the prediction of CO2 emission in China using input variables which include gross economic output, urbanization, industrialization, urbanization, economic structure, total population, energy consumption and energy productivity. The findings revealed that the NARX model presented a more accurate prediction of CO2 emissions in China compared to backpropagation neural networks and linear regression models. In a comparative study between NARX models and support vector machine for prediction of energy consumption in non-residential buildings, Koschwitz et al. [41] reported that NARX models displayed a better prediction of the energy consumption in the non-commercial buildings compared to support vector machine model. One major advantage of the NARX neural networks over the other modelling techniques is the ability to employed historical time series data for predictive modelling. Similarly, a multilayer perceptron type of Neural Network has been reported to accurately predict CO2 absorption and solubility with low MSE values. [42,43].

3.3. Sensitivity Analysis to Determine the Level of Importance of the Input Variables

The level of importance of the five variables namely GDP per capita, Population, Energy consumption per capita, coal consumption, and oil consumption are depicted in Figure 9. As shown in Figure 9a–d, the five variables have a significant level of importance on the CO2 emissions in Malaysia, Thailand, Indonesia, and Thailand based on the sensitivity analyses. This is evidence from the values of the normalized importance which is greater than 50% for each of the variables. However, coal consumption with 100% normalized importance can be said to have the most significant influence on CO2 emissions in Malaysia, Thailand, Indonesia, and the Philippines. According to the recent report by the International Energy Agency, Coal consumption used for electricity generation in the four ASEAN countries has been on the increase in the last decade [5]. For instance, coal consumption increased from 19.312 TOE December 2017 to 21.122 TOE in December 2018. Based on the Energy Commission Malaysia 2017 report, over 90% of electricity, generated in Peninsular Malaysia was from fossil fuel in which coal accounted for about 50% [44]. Similarly, coal consumption used for electricity generation in Thailand has also increased over the years [45]. Most of the electricity generated in Thailand are from fossil fuel, such as coal, natural gas and fuel oil [46]. The consumption of coal for electricity generation has also increased in Indonesia over the years [5]. Studies have shown that coal consumptions in Indonesia is expected to rise up to 157 million tons by 2027 thereby increasing the share of electricity generation from coal by 33.6% [47]. Hence, upwards in the amount of CO2 emissions in Indonesia is also expected if there no commitment to diversify the electricity generation mix. In the Philippines, 70% of electricity generation is from fossil fuel mainly natural gas and coal out of which coal accounted for 42.8% as of 2018 [48]. There has been an annual increase in coal consumptions over the year which directly translated to an increase in CO2 emissions. The causal relationship between coal consumption, CO2 emission and economic growth has been reported by Al-Mulali and Che Sab [49]. The study which investigated ten countries revealed that coal consumption and CO2 emissions have a short-run positive bidirectional relationship. This implies that an increase in coal consumption directly translated into more CO2 emissions. Also, in conformity with these study authors such as Chang (2010) [50] and Chandran Govindaraju and Tang [28] have also reported the direct relationship between coal consumption and CO2 emissions in China. Besides the significant influence of coal consumption on CO2 emissions, other factors such as oil consumptions, energy consumption per capita, population, and GDP per capita has been reported to have various level of relationship with CO2 emissions. Apergis and Payne [51,52], Ozturk and Acaravci [26] and Chang [50] reported that GDP and CO2 emissions in six America countries and China have short-run unidirectional relationship. Also, Al-mulali et al. (2013) [53], Chandran Govindaraju and Tang [28], Saboori and Sulaiman [54] also reported a unidirectional relationship between energy consumption and CO2 emissions in NEBA and five ASEAN countries. However, the modelling techniques did not capture the level of importance of these variables on the CO2 emissions which we have reported in this study.

4. Conclusions and Policy Implications

This study demonstrates the application of a non-linear autoregressive exogenous neural network for modelling the prediction of CO2 emissions in Malaysia, Thailand, Indonesia, and the Philippines. Based on the non-linear relationships between the CO2 emissions and the input variables (energy and economic indicators), fifteen NARX model architectures were developed and trained using historical data obtained for the four ASEAN countries. Among the various NARX model trained, optimized configurations of 5-29-1, 5-19-1, 5-17-1, 5-13-1 were obtained using datasets for Malaysia, Thailand, Indonesia, and the Philippines, respectively. Robust predictions of CO2 emissions in Malaysia, Thailand, Indonesia, and the Philippines were achieved using the optimized NARX model configurations. The predicted CO2 emissions in the four ASEAN countries were in close agreement with the actual values from the historical data, as indicated by the high values of R2. The findings also show that the NARX models have high predictability with minimal residual and MSE due to its excellent learning ability, rapid convergence, and superior generalization. The level of importance analyses of the various factors revealed that the economic and energy indicators significantly influenced CO2 emissions in the four ASEAN countries. Nevertheless, coal consumption displayed the highest level of importance on the CO2 emissions in the four ASEAN countries. This calls for more commitment by the government and stakeholders in each of the four ASEAN countries for more stringent legislation and policies on the modalities to drastically reduce CO2 emissions. Although, the governments of each of the four countries investigated have made various commitments to reduce energy-related CO2 emissions. These include the mandatory use of smart street lighting systems and sustainable urban transport policies in Indonesia; the national biofuel policy, national renewable energy policy and the low carbon cities framework in Malaysia; the National Renewable Energy program in the Philippines; Climate Change Master plan in Thailand. However, more commitment is needed to reduce the CO2 emissions in each of the four countries investigated. Governments of the four ASEAN countries should give it a serious thought of domesticating the climate-neutral now policy proposed by the United Nation Climate Change which entails the measurements of the GHGs emissions, proposed action plans to reduce the GHGs emissions as much as possible, and the implementation of carbon credit, which involve compensating the GHGs emissions that cannot be avoided. Moreover, the ongoing efforts by the four countries in incorporating renewable energy as part of the primary energy mix should be intensified. The regional efforts should be strengthened to address the issues of CO2 emissions in the ASEAN region. Considering the peculiarity of each country, collective policies targeting the reduction of CO2 emissions could be formulated and each member committed to implementation. Based on the level of importance analysis in this study, direct target and strategies on the gradual replacement of coal as the main sources of energy in ASEAN countries should be put in place. Although, the NARX model used in this study has been demonstrated to be robust in predicting the CO2 emissions in the four ASEAN countries investigated, the findings cannot be generalized and adapted for other ASEAN countries dues to the uniqueness of each country’s data. However, the NARX model could also be employed in a broader perspective by incorporating more data from other ASEAN countries. This could provide a robust CO2 emission prediction for the whole ASEAN region and hence provide opportunity for a regional effort to develop modalities for reduction of CO2 emissions in the region.

Author Contributions

S.I.M. and B.V.A. conceptualized the idea. F.O.A., B.V.A. and N.M. wrote the draft manuscript. SIM reviewed the draft manuscript for technicality. S.I.M. acquired the funding for the research. All authors have read and agreed to the published version of the manuscript.

Funding

This research and the APC was funded by Universiti Tenaga Nasional by grant number BOLD2025 (10436494/B2019141).

Acknowledgments

B.V.A., S.I.M., and S.M. acknowledge the financial support of Universiti Tenaga Nasional, Malaysia.

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Abbreviations

ASEANAssociation of Southeast Asian Nations
ARIMAAuto-Regressive Integrated Moving Average
ARDLAutoregressive Distributive Lag
CO2Carbon dioxide
EKCEnvironmental Kuznets Curve
FDIForeign Direct Investment
MSEMean Square Error
NARXNon-linear Autoregressive neural network with Exogenous input
MAPEMean Absolute Percentage Error
TOETonne of Oil Equivalent

Appendix A

Table A1. Detail results of the optimization of the hidden neurons at varying delay using training data for Malaysia.
Table A1. Detail results of the optimization of the hidden neurons at varying delay using training data for Malaysia.
Number of Delay13579111315
Hidden NeuronsMSERMSERMSERMSERMSERMSERMSERMSER
12.54 × 10−30.8711.56 × 10−30.9071.85 × 10−30.8992.32 × 10−30.8823.24 × 10−30.8132.80 × 10−30.8441.77 × 10−30.8183.42 × 10−30.892
32.13 × 10−30.8929.46 × 10−40.9392.01 × 10−30.9005.75 × 10−30.6413.84 × 0−30.7911.74 × 10−30.9124.24 × 10−40.9791.35 × 10−30.936
51.78 × 10−30.9142.47 × 10−40.8751.24 × 10−30.9213.68 × 10−40.9785.72 × 10−30.9051.42 × 10−30.9425.07 × 10−30.8781.21 × 10−20.516
71.80 × 10−30.8962.88 × 10−30.8811.43 × 10−30.9171.43 × 10−30.9531.42 × 10−40.9942.15 × 10−201.002.58 × 10−30.9418.34 × 10−170.999
92.07 × 10−30.9094.59 × 10−30.8101.14 × 10−30.9241.28 × 10−50.9992.66 × 10−30.9891.43 × 10−30.9649.05 × 10−40.9584.08 × 10−30.897
111.68 × 10−30.9121.08 × 10−30.9321.24 × 10−30.9288.17 × 10−40.9611.46 × 10−30.9642.78 × 10−50.9981.06 × 10−40.9957.46 × 10−50.996
131.17 × 10−30.9311.55 × 10−30.9011.74 × 10−40.9911.83 × 10−40.926.57 × 10−50.9984.13 × 10−40.9792.79 × 10−30.93342.05 × 10−90.999
151.75 × 10−30.8921.22 × 10−30.9331.05 × 10−30.9588.39 × 10−40.9744.16 × 10−50.9988.37 × 10−30.7022.51 × 10−40.9972.08 × 10−30.949
175.59 × 10−40.9642.47 × 10−30.9862.03 × 10−40.9891.63 × 10−40.9521.32 × 10−30.9373.16 × 10−50.9981.32 × 10−40.9983.71 × 10−50.999
191.08 × 10−30.9364.36 × 10−30.9761.22 × 10−30.9466.96 × 10−40.968.69 × 10−30.9634.20 × 10−201.008.27 × 10−50.9971.32 × 10−30.974
211.650.8981.10 × 10−30.9453.46 × 10−40.9891.45 × 10−30.9311.55 × 10−30.9922.54 × 10−30.8838.38 × 10−90.9994.24 × 10−40.977
232.38 × 10−30.8877.88 × 10−40.9578.78 × 10−40.9753.25 × 10−30.9833.15 × 10−40.9883.37 × 10−170.9993.25 × 10−160.9997.37 × 10−40.986
254.54 × 10−40.9741.24 × 10−30.9033.28 × 10−30.8857.19 × 10−40.9691.17 × 10−30.9782.99 × 10−40.9874.41 × 10−90.9995.86 × 10−40.977
279.17 × 10−40.9578.12 × 10−50.9968.48 × 10−40.9532.35 × 10−40.9941.03 × 10−40.9951.11 × 10−40.9863.63 × 10−40.9841.61 × 10−20.699
291.37 × 10−30.9258.71 × 10−30.9651.15 × 10−40.9984.28 × 10−30.9244.84 × 10−30.9152.09 × 10−80.9993.92 × 10−210.9992.16 × 10−40.998
Table A2. Detail results of the optimization of the hidden neurons at varying delay using training data for Thailand.
Table A2. Detail results of the optimization of the hidden neurons at varying delay using training data for Thailand.
Number of Delays13579111315
Hidden NeuronMSERMSERMSERMSERMSERMSERMSERMSER
12.3 × 10−30.3817.07 × 10−40.9337.37 × 10−40.9448.96 × 10−30.5487.69 × 10−30.1171.44 × 10−30.8981.12 × 10−30.9672.4 × 10−30.805
31.31 × 10−30.9125.52 × 10−30.5672.49 × 10−30.8062.23 × 10−30.8272.77 × 10−30.7411.560.9053.69 × 10−40.9978.2 × 10−40.957
51.1010.9378.00 × 10−30.6131.93 × 10−30.8691.30 × 10−30.9002.55 × 10−50.9982.85 × 10−30.7861.44 × 10−40.999.2 × 10−190.999
78.14 × 10−40.9539.76 × 10−40.9481.5 × 10−30.8991.14 × 100.9311.21 × 10−50.9997.27 × 10−40.9364.06 × 10−30.9958.8 × 10−120.999
92.37 × 10−30.8634.19 × 10−30.9731.65 × 10−30.8511.59 × 10−30.8751.77 × 10−30.8048.56 × 10−50.9977.53 × 10−40.9586.2 × 10−80.999
113.64 × 10−40.9791.48 × 10−30.8981.86 × 10−70.9991.51 × 10−80.9999.25 × 10−50.9941.87 × 10−30.9234.19 × 10−30.7725.9 × 10−50.997
137.88 × 10−30.6271.04 × 10−40.9925.19 × 10−40.9665.02 × 10−40.9471.81 × 10−30.9095.88 × 10−40.9647.93 × 10−90.9994.7 × 10−30.813
151.88 × 10−30.8921.71 × 10−30.8814.86 × 10−30.6131.05 × 10−30.9688.54 × 10−40.9233.02 × 10−30.9981.14 × 10−30.9142.5 × 10−30.914
178.71 × 10−30.9516.62 × 10−40.9545.56 × 10−40.953.94 × 10−30.8175.30 × 10−50.9971.11 × 10−30.9375.08 × 10−60.9996.9 × 10−180.999
195.13 × 10−30.9633.7 × 10−100.9999.38 × 10−50.9932.86 × 10−40.9848.52 × 10−50.9956.90 × 10−40.9564.2 × 10−230.9991.9 × 10−50.879
218.86 × 10−40.9492.75 × 10−40.9593.03 × 10−30.8472.45 × 10−60.9994.44 × 10−30.9638.09 × 10−40.9545.54 × 10−50.9977.9 × 10−50.994
233.26 × 10−30.9845.54 × 10−30.7756.32 × 10−190.9991.32 × 10−40.9981.52 × 10−30.9074.58 × 10−80.9991.21 × 10−90.9993.1 × 10−70.999
252.62 × 10−30.8748.83 × 10−40.9497.36 × 10−40.9454.36 × 10−40.9717.21 × 10−60.9996.86 × 10−210.9998.92 × 10−40.9319.5 × 10−50.992
272.53 × 10−30.8562.81 × 10−30.8414.28 × 10−130.9991.52 × 10−40.9931.11 × 10−50.9996.12 × 10−110.9991.08 × 10−40.9961.7 × 10−30.994
291.98 × 10−40.9877.06 × 10−40.9746.58 × 10−30.5462.45 × 10−50.9981.84 × 10−40.9929.78 × 10−40.9992.72 × 10−50.9992.4 × 10−30.999
Table A3. Detail results of the optimization of the hidden neurons at varying delay using training data for Indonesia.
Table A3. Detail results of the optimization of the hidden neurons at varying delay using training data for Indonesia.
Delay13579111315
Hidden NeuronMSERMSERMSERMSERMSERMSERMSERMSER
12.71 × 10−30.8522.75 × 10−30.8116.84 × 10−40.9461.47 × 10−30.8372.14 × 10−30.7536.96 × 10−50.9932.14 × 10−30.6971.25 × 10−30.932
39.91 × 10−40.9358.39 × 10−40.8855.21 × 10−40.9358.73 × 10−40.9293.44 × 10−40.9527.33 × 10−40.8821.82 × 10−30.8658.37 × 10−40.904
52.17 × 10−30.8554.34 × 10−40.9651.74 × 10−30.8322.85 × 10−30.8691.18 × 10−30.9113.98 × 10−30.6846.76 × 10−30.9915.15 × 10−40.939
71.52 × 10−30.9091.71 × 10−40.991.12 × 10−30.9233.44 × 10−30.5631.28 × 10−30.9842.59 × 10−30.9811.17 × 10−30.8984.89 × 10−30.739
92.07 × 10−30.8681.54 × 10−30.8575.68 × 10−50.9963.46 × 10−30.8151.19 × 10−50.9994.81 × 10−30.4199.86 × 10−40.9012.54 × 10−40.972
111.44 × 10−30.9018.28 × 10−50.9922.09 × 10−30.9919.12 × 10−30.9918.07 × 10−50.9919.27 × 10−60.9957.8 × 10−100.9996.69 × 10−30.959
131.50 × 10−30.9193.23 × 10−40.9621.08 × 10−30.93261.16 × 10−30.9311.41 × 10−30.8318.68 × 10−40.861.82 × 10−40.9654.73 × 10−50.998
155.54 × 10−40.9511.39 × 10−40.9896.39 × 10−60.9992.01 × 10−30.8917.08 × 10−50.9935.01 × 10−30.972.54 × 10−30.9027.45 × 10−30.995
174.51 × 10−30.5711.04 × 10−30.9111.86 × 10−40.9862.09 × 10−60.9991.38 × 10−30.9221.35 × 10−30.8048.04 × 10−70.9999.28 × 10−50.994
194.63 × 10−30.8472.70 × 10−40.9851.83 × 10−30.9451.19 × 10−30.891.36 × 10−130.9994.65 × 10−30.4076.11 × 10−40.9531.4 × 10−130.999
212.13 × 10−30.8127.18 × 10−60.9991.23 × 10−40.9911.31 × 10−30.8642.87 × 10−40.9632.17 × 10−30.7341.94 × 10−30.8671.29 × 10−40.996
231.12 × 10−30.9387.83 × 10−40.9474.01 × 10−40.97281.98 × 10−30.9811.06 × 10−30.8191.8 × 10−110.9995.14 × 10−40.9299.79 × 10−50.987
255.940.9576.06 × 10−40.9568.70 × 10−40.9923.03 × 10−30.973.35 × 10−30.9973.86 × 10−40.9695.35 × 10−80.9991.75 × 10−30.928
271.75 × 10−30.8746.69 × 10−40.9521.57 × 10−30.8571.83 × 10−30.9861.41 × 10−30.8472.46 × 10−30.8292.32 × 10−30.7592.12 × 10−60.999
294.37 × 10−40.9842.05 × 10−40.9861.75 × 10−170.9991.25 × 10−210.9995.82 × 10−40.9543.05 × 10−30.9673.02 × 10−30.9922.0 × 10−190.999
Table A4. Detail results of the optimization of the hidden neurons at varying delay using training data for Philippines.
Table A4. Detail results of the optimization of the hidden neurons at varying delay using training data for Philippines.
Delay13579111315
Hidden NeuronMSERMSERMSERMSERMSERMSERMSERMSER
13.55 × 10−40.9554.27 × 10−40.9552.12 × 10−30.9784.06 × 10−30.9661.14 × 10−60.9998.8 × 10−50.9921.21 × 10−40.9859.64 × 10−50.989
33.73 × 10−40.966.95 × 10−40.9372.08 × 10−30.9816.64 × 10−30.9376.92 × 10−50.9922.35 × 10−40.9731.05 × 10−40.9952.11 × 10−30.976
53.13 × 10−40.9672.32 × 10−40.9731.13 × 10−30.9042.56 × 10−40.9075.61 × 10−50.9946.12 × 10−40.9921.41 × 10−190.9995.59 × 10−50.995
78.76 × 10−40.9143.41 × 10−40.9691.64 × 10−40.9872.71 × 10−40.9731.94 × 10−40.9797.09 × 10−50.9944.08 × 10−110.9993.88 × 10−50.996
92.61 × 10−40.9732.40 × 10−80.9991.45 × 10−30.9067.52 × 10−50.9939.48 × 10−50.9934.98 × 10−40.9574.56 × 10−60.9995.95 × 10−50.994
114.24 × 10−40.9482.16 × 10−40.9772.56 × 10−40.9812.9 × 10−130.9992.04 × 10−40.9836.07 × 10−50.9951.21 × 10−40.9925.51 × 10−40.969
131.18 × 10−40.9797.92 × 10−50.9911.95 × 10−50.9982.14 × 10−50.9987.89 × 10−50.9938.48 × 10−50.9922.46 × 10−40.9781.32 × 10−200.999
151.14 × 10−30.8272.48 × 10−50.9973.31 × 10−40.9751.73 × 10−40.9854.36 × 10−50.9952.27 × 10−50.9983.66 × 10−50.9969.91 × 10−80.999
173.17 × 10−40.9669.55 × 10−40.9894.99 × 10−60.9998.8 × 10−100.9994.86 × 10−50.9966.65 × 10−50.9944.32 × 10−50.9971.21 × 10−50.999
198.43 × 10−50.9895.94 × 10−50.9934.56 × 10−50.9967.15 × 10−40.9364.18 × 10−100.9996.22 × 10−90.9992.82 × 10−50.9982.22 × 10−40.979
213.59 × 10−40.9662.09 × 10−40.9813.56 × 10−40.9741.64 × 10−30.9081.09 × 10−60.9992.11 × 10170.9993.42 × 10−50.9961.49 × 10−80.999
237.05 × 10−40.9921.47 × 10−40.9831.81 × 10−40.9821.45 × 10−30.9011.62 × 10−40.9876.85 × 10−50.9938.38 × 10−60.9995.22 × 10−40.972
252.38 × 10−40.9711.76 × 10−50.9985.59 × 10−40.9317.29 × 10−50.9956.97 × 10−190.9992.23 × 10−50.9983.45 × 10−40.9592.89 × 10−50.998
275.80 × 10−40.9294.94 × 10−50.9941.26 × 10−40.9877.90 × 10−50.9932.61 × 10−−50.9974.36 × 10170.9994.91 × 10−140.9997.96 × 10−60.999
291.24 × 10−40.9886.29 × 10−50.9941.78 × 10−40.9812.18 × 10−40.9764.88 × 10−40.9356.25 × 10−70.9998.82 × 10−200.9998.92 × 10−60.999

References

  1. Wang, Y.; Chen, L.; Kubota, J. The relationship between urbanization, energy use and carbon emissions: Evidence from a panel of Association of Southeast Asian Nations (ASEAN) countries. J. Clean. Prod. 2016, 112, 1368–1374. [Google Scholar] [CrossRef]
  2. Tongsopit, S.; Kittner, N.; Chang, Y.; Aksornkij, A.; Wangjiraniran, W. Energy security in ASEAN: A quantitative approach for sustainable energy policy. Energy Policy 2016, 90, 60–72. [Google Scholar] [CrossRef]
  3. Shi, X. The future of ASEAN energy mix: A SWOT analysis. Renew. Sustain. Energy Rev. 2016, 53, 672–680. [Google Scholar] [CrossRef]
  4. Silberglitt, R.; Kimmel, S. Energy scenarios for Southeast Asia. Technol. Forecast. Soc. Chang. 2015, 101, 251–262. [Google Scholar] [CrossRef]
  5. International Energy Agency Southeast Asia Energy Outlook 2017; Southeast Asia Energy Outlook: Paris, France, 2017.
  6. Mensah, I.A.; Sun, M.; Gao, C.; Omari-Sasu, A.Y.; Zhu, D.; Ampimah, B.C.; Quarcoo, A. Analysis on the nexus of economic growth, fossil fuel energy consumption, CO2 emissions and oil price in Africa based on a PMG panel ARDL approach. J. Clean. Prod. 2019, 228, 161–174. [Google Scholar] [CrossRef]
  7. Mardani, A.; Štreimikienė, D.; Cavallaro, F.; Loganathan, N.; Khoshnoudi, M. Carbon dioxide (CO2) emissions and economic growth: A systematic review of two decades of research from 1995 to 2017. Sci. Total Environ. 2019, 649, 31–49. [Google Scholar] [CrossRef]
  8. Kasman, A.; Duman, Y.S. CO2 emissions, economic growth, energy consumption, trade and urbanization in new EU member and candidate countries: A panel data analysis. Econ. Model. 2015, 44, 97–103. [Google Scholar] [CrossRef]
  9. Tzeremes, P. Time-varying causality between energy consumption, CO2 emissions, and economic growth: Evidence from US states. Envrion. Sci. Pollut. Res. 2017, 25, 6044–6060. [Google Scholar] [CrossRef]
  10. Ozcan, B.; Tzeremes, P.G.; Tzeremes, N.G. Energy consumption, economic growth and environmental degradation in OECD countries. Econ. Model. 2020, 84, 203–213. [Google Scholar] [CrossRef]
  11. Zhu, H.; Duan, L.; Guo, Y.; Yu, K. The effects of FDI, economic growth and energy consumption on carbon emissions in ASEAN-5: Evidence from panel quantile regression. Econ. Model. 2016, 58, 237–248. [Google Scholar] [CrossRef] [Green Version]
  12. Heidari, H.; Katircioğlu, S.T.; Saeidpour, L. Economic growth, CO2 emissions, and energy consumption in the five ASEAN countries. Int. J. Electr. Power Energy Syst. 2015, 64, 785–791. [Google Scholar] [CrossRef]
  13. Salman, M.; Long, X.; Dauda, L.; Mensah, C.N.; Muhammad, S. Different impacts of export and import on carbon emissions across 7 ASEAN countries: A panel quantile regression approach. Sci. Total Environ. 2019, 686, 1019–1029. [Google Scholar] [CrossRef] [PubMed]
  14. Nugraha, A.T.; Osman, N.H. CO2 emissions, economic growth, energy consumption, and household expenditure for Indonesia: Evidence from cointegration and vector error correction model. Int. J. Energy Econ. Policy 2019, 9, 291–298. [Google Scholar] [CrossRef]
  15. Fatima, S.; Ali, S.S.; Zia, S.S.; Hussain, E.; Fraz, T.R.; Khan, M.S. Forecasting Carbon Dioxide Emission of Asian Countries Using ARIMA and Simple Exponential Smoothing Models. Int. J. Econ. Environ. Geol. 2019, 10, 64–69. [Google Scholar] [CrossRef]
  16. Hanif, I.; Raza, S.M.F.; Gago-De-Santos, P.; Abbas, Q. Fossil fuels, foreign direct investment, and economic growth have triggered CO2 emissions in emerging Asian economies: Some empirical evidence. Energy 2019, 171, 493–501. [Google Scholar] [CrossRef]
  17. Xu, G.; Schwarz, P.; Yang, H. Determining China’s CO2 emissions peak with a dynamic nonlinear artificial neural network approach and scenario analysis. Energy Policy 2019, 128, 752–762. [Google Scholar] [CrossRef]
  18. Kahouli, B. The causality link between energy electricity consumption, CO2 emissions, R&D stocks and economic growth in Mediterranean countries (MCs). Energy 2018, 145, 388–399. [Google Scholar] [CrossRef]
  19. Marcjasz, G.; Uniejewski, B.; Weron, R. On the importance of the long-term seasonal component in day-ahead electricity price forecasting with NARX neural networks. Int. J. Forecast. 2019, 35, 1520–1532. [Google Scholar] [CrossRef]
  20. Alsaffar, M.A.; Ayodele, B.V.; Mustapa, S.I. Scavenging carbon deposition on alumina supported cobalt catalyst during renewable hydrogen-rich syngas production by methane dry reforming using artificial intelligence modeling technique. J. Clean. Prod. 2020, 247, 119168. [Google Scholar] [CrossRef]
  21. Yang, Z.; Zhao, Y. Energy consumption, carbon emissions, and economic growth in India: Evidence from directed acyclic graphs. Econ. Model. 2014, 38, 533–540. [Google Scholar] [CrossRef]
  22. Omri, A.; Nguyen, D.K.; Rault, C. Causal interactions between CO2 emissions, FDI, and economic growth: Evidence from dynamic simultaneous-equation models. Econ. Model. 2014, 42, 382–389. [Google Scholar] [CrossRef] [Green Version]
  23. Feng, Y.; Chen, S.; Zhang, L. System dynamics modeling for urban energy consumption and CO2 emissions: A case study of Beijing, China. Ecol. Model. 2013, 252, 44–52. [Google Scholar] [CrossRef]
  24. Pao, H.-T.; Yu, H.-C.; Yang, Y.-H. Modeling the CO2 emissions, energy use, and economic growth in Russia. Energy 2011, 36, 5094–5100. [Google Scholar] [CrossRef]
  25. Alam, M.J.; Begum, I.A.; Buysse, J.; Rahman, S.; Van Huylenbroeck, G. Dynamic modeling of causal relationship between energy consumption, CO2 emissions and economic growth in India. Renew. Sustain. Energy Rev. 2011, 15, 3243–3251. [Google Scholar] [CrossRef]
  26. Maziarz, M. A review of the Granger-causality fallacy. J. Philos. Econ. Reflect. Econ. Soc. Issues 2015, 8, 86–105. [Google Scholar]
  27. Liu, X.; Mao, G.; Ren, J.; Li, R.Y.M.; Guo, J.; Zhang, H.-L. How might China achieve its 2020 emissions target? A scenario analysis of energy consumption and CO2 emissions using the system dynamics model. J. Clean. Prod. 2015, 103, 401–410. [Google Scholar] [CrossRef]
  28. Govindaraju, V.C.; Tang, C.F. The dynamic links between CO2 emissions, economic growth and coal consumption in China and India. Appl. Energy 2013, 104, 310–318. [Google Scholar] [CrossRef]
  29. Hossain, M.A.; Ayodele, B.V.; Cheng, C.K.; Khan, M.R. Artificial neural network modeling of hydrogen-rich syngas production from methane dry reforming over novel Ni/CaFe2O4 catalysts. Int. J. Hydrog. Energy 2016, 41, 11119–11130. [Google Scholar] [CrossRef] [Green Version]
  30. Saghafi, M.; Ghofrani, M.B. Real-time estimation of break sizes during LOCA in nuclear power plants using NARX neural network. Nucl. Eng. Technol. 2019, 51, 702–708. [Google Scholar] [CrossRef]
  31. Rahmani, A.; Deihimi, A. Reduction of harmonic monitors and estimation of voltage harmonics in distribution networks using wavelet analysis and NARX. Electr. Power Syst. Res. 2020, 178, 106046. [Google Scholar] [CrossRef]
  32. Ezzeldin, R.; Hatata, A.Y. Application of NARX neural network model for discharge prediction through lateral orifices. Alex. Eng. J. 2018, 57, 2991–2998. [Google Scholar] [CrossRef]
  33. Garson, G.D. Comparison of Neural Network Analysis of Social Science Data. Soc. Sci. Comput. Rev. 1991, 9, 399–434. [Google Scholar] [CrossRef]
  34. Cheng, J.; Ji, Z.; Li, M.; Dai, J. Study of a noninvasive blood glucose detection model using the near-infrared light based on SA-NARX. Biomed. Signal Process. Control 2020, 56, 101694. [Google Scholar] [CrossRef]
  35. Pazikadin, A.R.; Rifai, D.; Ali, K.; Malik, M.Z.; Abdalla, A.N.; Faraj, M.A. Solar irradiance measurement instrumentation and power solar generation forecasting based on Artificial Neural Networks (ANN): A review of five years research trend. Sci. Total Envrion. 2020, 715, 136848. [Google Scholar] [CrossRef]
  36. Giannopoulos, A.; Aider, J.-L. Prediction of the dynamics of a backward-facing step flow using focused time-delay neural networks and particle image velocimetry data-sets. Int. J. Heat Fluid Flow 2020, 82, 108533. [Google Scholar] [CrossRef]
  37. Das, S.K.; Basudhar, P.K. Prediction of residual friction angle of clays using artificial neural network. Eng. Geol. 2008, 100, 142–145. [Google Scholar] [CrossRef]
  38. Zounemat-Kermani, M.; Stephan, D.; Hinkelmann, R. Multivariate NARX neural network in prediction gaseous emissions within the influent chamber of wastewater treatment plants. Atmos. Pollut. Res. 2019, 10, 1812–1822. [Google Scholar] [CrossRef]
  39. Oko, E.; Wang, M.; Zhang, J. Neural network approach for predicting drum pressure and level in coal-fired subcritical power plant. Fuel 2015, 151, 139–145. [Google Scholar] [CrossRef]
  40. Alcan, G.; Unel, M.; Aran, V.; Yilmaz, M.; Gurel, C.; Koprubasi, K. Predicting NOx emissions in diesel engines via sigmoid NARX models using a new experiment design for combustion identification. Measurements 2019, 137, 71–81. [Google Scholar] [CrossRef] [Green Version]
  41. Koschwitz, D.; Frisch, J.; Van Treeck, C. Data-driven heating and cooling load predictions for non-residential buildings based on support vector machine regression and NARX Recurrent Neural Network: A comparative study on district scale. Energy 2018, 165, 134–142. [Google Scholar] [CrossRef]
  42. Pakzad, P.; Mofarahi, M.; Izadpanah, A.A.; Afkhamipour, M. Experimental data, thermodynamic and neural network modeling of CO2 absorption capacity for 2-amino-2-methyl-1-propanol (AMP) + Methanol (MeOH) + H2O system. J. Nat. Gas Sci. Eng. 2020, 73, 103060. [Google Scholar] [CrossRef]
  43. Hamzehie, M.; Mazinani, S.; Davardoost, F.; Mokhtare, A.; Najibi, H.; Van Der Bruggen, B.; Darvishmanesh, S. Developing a feed forward multilayer neural network model for prediction of CO2 solubility in blended aqueous amine solutions. J. Nat. Gas Sci. Eng. 2014, 21, 19–25. [Google Scholar] [CrossRef]
  44. Energy in Malaysia. Toward a World-Class Energy Sector; Suruhanjaya Tenaga: Putrajaya, Malaysia, 2017; Volume 12, pp. 1–45. [Google Scholar]
  45. Kenney, W. Energy Conservation in the Process Industries; Elsevier: Amsterdam, The Netherlands, 1984; pp. 1–17. [Google Scholar]
  46. Gross Energy Generation and Purchase; EGAT: Bangkok, Thailand, 2019.
  47. Handbook of Energy & Economic Statistics of Indonesia 2018; Ministry of Energy and Mineral Resources, Republic of Indonesia: Jakarta, Indonesia, 2018.
  48. Philippines Power Statistics; Department of Energy: Taguig, Philippines, 2019.
  49. Al-Mulali, U.; Sab, C.N.B.C. Electricity consumption, CO2 emission, and economic growth in the Middle East. Energy Sources Part B Econ. Plan. Policy 2018, 13, 257–263. [Google Scholar] [CrossRef]
  50. Chang, C.-C. A multivariate causality test of carbon dioxide emissions, energy consumption and economic growth in China. Appl. Energy 2010, 87, 3533–3537. [Google Scholar] [CrossRef]
  51. Apergis, N.; Payne, J.E. CO2 emissions, energy usage, and output in Central America. Energy Policy 2009, 37, 3282–3286. [Google Scholar] [CrossRef]
  52. Apergis, N.; Payne, J.E. The emissions, energy consumption, and growth nexus: Evidence from the commonwealth of independent states. Energy Policy 2010, 38, 650–655. [Google Scholar] [CrossRef]
  53. Al-Mulali, U.; Fereidouni, H.G.; Lee, J.Y.; Sab, C.N.B.C. Exploring the relationship between urbanization, energy consumption, and CO2 emission in MENA countries. Renew. Sustain. Energy Rev. 2013, 23, 107–112. [Google Scholar] [CrossRef]
  54. Saboori, B.; Sulaiman, J. Environmental degradation, economic growth and energy consumption: Evidence of the environmental Kuznets curve in Malaysia. Energy Policy 2013, 60, 892–905. [Google Scholar] [CrossRef]
Figure 1. Share of Southeast Asia in global economic indicators, energy consumption and CO2 emissions [5].
Figure 1. Share of Southeast Asia in global economic indicators, energy consumption and CO2 emissions [5].
Processes 08 01529 g001
Figure 2. The model architecture of the NARX neural network.
Figure 2. The model architecture of the NARX neural network.
Processes 08 01529 g002
Figure 3. General representation of the optimized architecture for the NARX model used for prediction of CO2 emissions in Malaysia, Thailand, Indonesia, and Philippines.
Figure 3. General representation of the optimized architecture for the NARX model used for prediction of CO2 emissions in Malaysia, Thailand, Indonesia, and Philippines.
Processes 08 01529 g003
Figure 4. The validation performance of the optimized neuron at various epoch for (a) Malaysia; (b) Thailandl (c) Indonesia; and (d) the Philippines.
Figure 4. The validation performance of the optimized neuron at various epoch for (a) Malaysia; (b) Thailandl (c) Indonesia; and (d) the Philippines.
Processes 08 01529 g004
Figure 5. Dispersion plot showing the actual and predicted CO2 emissions from (a) Malaysia, (b) Thailand, (c) Indonesia, (d) Philippines.
Figure 5. Dispersion plot showing the actual and predicted CO2 emissions from (a) Malaysia, (b) Thailand, (c) Indonesia, (d) Philippines.
Processes 08 01529 g005
Figure 6. Residuals obtained from the NARX model prediction of CO2 emissions.
Figure 6. Residuals obtained from the NARX model prediction of CO2 emissions.
Processes 08 01529 g006
Figure 7. Comparison between the training, validation and testing responses of the NARX model. (a) Malaysia (b) Thailand, (c) Indonesia, (d) Philippines.
Figure 7. Comparison between the training, validation and testing responses of the NARX model. (a) Malaysia (b) Thailand, (c) Indonesia, (d) Philippines.
Processes 08 01529 g007
Figure 8. Autocorrelation analysis for the prediction of CO2 emissions in (a) Malaysia, (b) Thailand, (c) Indonesia, (d) the Philippines.
Figure 8. Autocorrelation analysis for the prediction of CO2 emissions in (a) Malaysia, (b) Thailand, (c) Indonesia, (d) the Philippines.
Processes 08 01529 g008aProcesses 08 01529 g008b
Figure 9. The level of importance of the input variables on the predicted CO2 emissions in (a) Malaysia, (b) Thailand, (c) Indonesia, (d) Philippines.
Figure 9. The level of importance of the input variables on the predicted CO2 emissions in (a) Malaysia, (b) Thailand, (c) Indonesia, (d) Philippines.
Processes 08 01529 g009
Table 1. Summary of literature on modelling of CO2 emissions.
Table 1. Summary of literature on modelling of CO2 emissions.
Countries InvestigatedMethodsVariablesReference
ASEAN (Malaysia, Thailand, Indonesia and the Philippines)Autoregressive Exogenous Neural Network ModelingEnergy Consumption per capita, Population, GDP per capitaThis study
ASEAN (Brunei, Indonesia, Malaysia, Philippines, Singapore, Thailand, Vietnam)Panel quantile regression approachExport and Import[13]
IndonesiaCointegration and Vector Error Correction ModelEnergy Consumption, GDP, House Expenditure[14]
Asian countries (Japan, Bangladesh, China, Pakistan, India, Sri Lanka, Iran, Singapore, and Nepal)Auto-Regressive Integrated Moving Average (ARIMA) and Simple Exponential Smoothing ModelsHeat and electricity, manufacturing industries, residential and commercial buildings, transport[15]
Asia countriesAutoregressive Distributive Lag model.Fossil fuel, FDI, GDP[16]
ChinaLinear regression, Backpropagation neural network, non-linear dynamic neural networkGDP, total population, urbanization rate, total energy consumption, percentage of coal consumption, percentage of non-fossil consumption [17]
Mediterranean countriesSystem of simultaneous equations using seemingly unrelated regression.Research and Development Stocks, GDP. Electricity Consumption, [18]
CroatiaEnvironmental Kuznets Curve (EKC) modelGDP[18]
ASEAN (Malaysia, Indonesia, Philippines, Singapore, and Thailand)Panel quantile regression modelForeign Direct Investment, GDP, Energy Consumption[11]
ASEAN (Malaysia, Indonesia, Philippines, Singapore, and Thailand)Panel smooth transition regression modelEnergy Consumption and GDP[12]
IndiaDirected acyclic graphsEnergy Consumption, GDP, fixed capita formation, and Trade openness[21]
Middle East, North Africa and Sub-Sahara AfricaSimultaneous Equation ModelsForeign direct investment, [22]
ChinaSystem Dynamic ModellingEnergy Consumptions[23]
RussiaEnvironmental Kuznets Curve (EKC) modelEnergy Consumptions and GDP[24]
IndiaSystem Dynamic ModellingEnergy Consumption and GDP[25]
Turkeyautoregressive distributed lag bounds testing approach of cointegrationEnergy consumption and economic indicator[26]
Table 2. Optimized hidden neuron obtained using data for Malaysia.
Table 2. Optimized hidden neuron obtained using data for Malaysia.
Model ArchitectureDelayHidden Neurons RR2MSE
5-25-11250.9740.9494.54 × 10−4
5-27-13270.9960.9928.12 × 10−5
5-29-15290.9980.9961.15 × 10−4
5-9-1790.9990.9981.28 × 10−5
5-15-19150.9980.9964.16 × 10−5
5-7-11170.9990.9982.15 × 10−20
5-29-1 13290.9990.9983.92 × 10−21
5-7-11570.9990.9988.34 × 10−17
Table 3. Optimized hidden neuron obtained using data for Thailand.
Table 3. Optimized hidden neuron obtained using data for Thailand.
Model ArchitectureDelayHidden NeuronsRR2MSE
5-29-11290.9870.9741.98 × 10−4
5-19-13190.9990.9983.72 × 10−10
5-23-15230.9990.9986.32 × 10−19
5-11-17110.9990.9981.51 × 10−8
5-25-19250.9990.9987.21 × 10−6
5-25-111 250.9990.9986.86 × 10−21
5-19-113190.9990.9984.15 × 10−23
5-5-1 1550.9990.9989.18 × 10−19
Table 4. Optimized hidden neuron obtained using data for Indonesia.
Table 4. Optimized hidden neuron obtained using data for Indonesia.
Model ArchitectureDelayHidden Neurons RR2MSE
5-29-11290.9840.9684.37 × 10−4
5-21-13210.9990.9987.18 × 10−6
5-29-15290.9990.9981.75 × 10−17
5-29-17290.9990.9981.25 × 10−21
5-19-19190.9990.9981.36 × 10−13
5-23-111230.9990.9981.83 × 10−11
5-25-113250.9990.9985.35 × 10−8
5-17-115170.9990.9982.02 × 10−19
Table 5. Optimized hidden neuron obtained using data for Philippines.
Table 5. Optimized hidden neuron obtained using data for Philippines.
Model ArchitectureDelayHidden NeuronsRR2MSE
5-23-11230.9920.9847.05 × 10−5
5-9-1390.9990.9982.40 × 10−8
5-17-15170.9990.9984.99 × 10−6
5-11-17110.9990.9982.94 × 10−13
5-25-19250.9990.9986.97 × 10−19
5-21-111210.9990.9982.11 × 10−17
5-29-113290.9990.9988.82 × 10−20
5-13-115130.9990.9981.32 × 10−20
Table 6. Details of the optimized NARX architecture used for the modeling.
Table 6. Details of the optimized NARX architecture used for the modeling.
MalaysiaThailandIndonesiaPhilippines
Input units5555
Optimized hidden neurons13191713
Optimized delay13131515
Output units1111
Table 7. Comparison of the present work with literature.
Table 7. Comparison of the present work with literature.
ReferenceModelling TechniquesObjectivesInput VariablesMeasurement of Accuracy Conclusions
This workNARX neural networkPredictive modelling of CO2 emissions in four ASEAN countries Energy consumption per capita, GDP per capita, population, oil consumption, coal consumptionMSE values of 6.93 × −4, 9.26 × 10−4, 2.03 × 10−3, and 7.99 × 10−4 obtained for CO2 emissions in Malaysia, Thailand, Indonesia, and the PhilippinesThe optimized NARX models efficiently predicted CO2 emissions I the four ASEAN countries with coal consumptions having the highest level of importance.
Alcan et al. [40]NARX neural networkPredictive modelling of NOx emissions from diesel enginesEngine speed, Manifold Absolute Pressure, Mass Air Flow, rail pressure, main and pilot injection fuel quantities, Main and pilot start of injectionsNot reportedThe optimized NARX model architecture predicted the NOx emissions with high degree of accuracy
Xu et al. [17]NARX neural networkPredict CO2 emission in ChinaGDP, total population, industrialization, urbanization, secondary sector, tertiary sector, total energy consumption, coal consumption, non-fossil consumption, energy productivity, investment in environment governanceRMSE = 0.0311Economic scale, secondary sector
(manufacturing and construction), coal consumption, industrialization,
and energy productivity significantly influenced CO2 emissions
Koschwitz et al. [41]NARX neural network and Support Vector MachinePrediction of energy consumption in commercial buildingdew point temperature, mean wind direction, mean wind velocity, outdoor temperature, precipitation intensity, precipitation quantity, relative humidity, school holiday time, working time scheduleMSE = 2.35–2.89NARX models displayed a better prediction of the energy consumption in the non-commercial buildings compared to support vector machine
Pakzad et al. [42]Multilayer Perceptron Neural Network Modelling the CO2 absorption CO2 partial pressure,
temperature, AMP concentration, and MeOH concentration
Average absolute relative deviation (AARD%) = 1.95The findings revealed that the ANN models predicted
CO2 absorption was in proximity with the observed.
Hamzehie et al. [43]Multilayer Perceptron Neural NetworkModelling the CO2 solubilityTemperature, pressure, overall concentration, apparent molecular weight of the mixture.MSE = 2.300 × 10−4The solubility of CO2 in the mixed aqueous solution was accurately predicted by the Neural Network.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mustapa, S.I.; Ayodele, F.O.; Ayodele, B.V.; Mohammad, N. Nexus between Energy Usability, Economic Indicators and Environmental Sustainability in Four ASEAN Countries: A Non-Linear Autoregressive Exogenous Neural Network Modelling Approach. Processes 2020, 8, 1529. https://doi.org/10.3390/pr8121529

AMA Style

Mustapa SI, Ayodele FO, Ayodele BV, Mohammad N. Nexus between Energy Usability, Economic Indicators and Environmental Sustainability in Four ASEAN Countries: A Non-Linear Autoregressive Exogenous Neural Network Modelling Approach. Processes. 2020; 8(12):1529. https://doi.org/10.3390/pr8121529

Chicago/Turabian Style

Mustapa, Siti Indati, Freida Ozavize Ayodele, Bamidele Victor Ayodele, and Norsyahida Mohammad. 2020. "Nexus between Energy Usability, Economic Indicators and Environmental Sustainability in Four ASEAN Countries: A Non-Linear Autoregressive Exogenous Neural Network Modelling Approach" Processes 8, no. 12: 1529. https://doi.org/10.3390/pr8121529

APA Style

Mustapa, S. I., Ayodele, F. O., Ayodele, B. V., & Mohammad, N. (2020). Nexus between Energy Usability, Economic Indicators and Environmental Sustainability in Four ASEAN Countries: A Non-Linear Autoregressive Exogenous Neural Network Modelling Approach. Processes, 8(12), 1529. https://doi.org/10.3390/pr8121529

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop