Nexus Between Carbon Emissions, Energy Consumption, Urbanization and Economic Growth in Asia: Evidence From Common Correlated Effects Mean Group Estimator (CCEMG)

This study investigates the dynamic relationships between carbon emission, urbanization, energy consumption, and economic growth in a panel of 42 Asian countries for the period 2000–2014 using dynamic common correlated effects panel data modeling. This study employs second generation cross-sectional Pesaran (J. Appl. Econom., 2007, 22(2), 265-312) panel unit root, Westerlund panel cointegration tests (Econom. Stat., 2007, 69(6), 709-748), and Pesaran’s (Econometrica, 2006, 74(4), 967-1012) common correlated effects mean group estimation technique. These approaches allow for cross-sectional dependence, and are robust to the presence of common factors, serial correlation, and slope heterogeneity. The Common Correlated Effect Mean Group test reveals a high average coefficient of 0.602 between carbon emission and energy consumption while low coefficients of 0.114 and 0.184 for the pairs of carbon emission-urbanization and carbon emission-GDP, respectively for the panel as a whole, suggesting a cointegration between carbon emission, urbanization, energy consumption, and economic growth. The results indicate that there is relatively high carbon emission especially for highly populated and geopolitical risk Asian countries in the short run. Findings reveal long run relationships between the variables, which is attributed to the on-going carbon taxation and energy prices. Our results are robust to PMG-ARDL estimator. Overall, these findings cast important implications on renewable energy policy and urban planning insights for the policymakers.


INTRODUCTION
As the debate surrounds anthropogenic carbon emissions and climate change typifies sustainable development dilemma, the global climate effort put forward by the Paris Agreement advocates urgent attention towards carbon mitigation and adaptation strategies. While economic development is crucial for the emerging markets, it is noteworthy to assess the nexus between energy consumption, urbanization, and carbon emission, especially from the Asian countries' perspective. Carbon emissions have been highly and widely researched in energy, finance, and economics literature. Carbon emissions have been linked to economic growth (Gao and Zhang, 2014;Arvin et al., 2015;Bekun et al., 2019;Mahembe and Odhiambo, 2019;Adedoyin et al., 2020), financial development Shahbaz et al., 2020), urbanization (Muhammad et al., 2020), trade openness (Farhani et al., 2013;Essandoh et al., 2020), energy consumption (Bekun et al., 2019;Yang et al., 2020), environmental pressures (Gingrich et al., 2011;, and renewable energy (O'Ryan et al., 2020).
Countries with high carbon emissions that are described to be highly technologically inclined and are among high-income countries, can reduce carbon emissions while developing regions and low-income countries with increasing volume and speed of emission do not have the sound financial system to curtail carbon emissions (Chen and Lee, 2020;Dong et al., 2020). The incessant and long-run crises rocking some countries, contribute to the rate of carbon emissions (Oreskes and Conway, 2008;Gingrich et al., 2011;Zou et al., 2016). Thus, the relationship between carbon emissions and other economic variables is constrained by per-capital income, economic development, technology advancement, geographical, and regional features Li et al., 2020;Tirkaso and Gren, 2020). Continental and regional issues are important to devise regional energy programs and improve strategies on rapid urbanization growth especially in Asia, where the continent receives the largest tourist arrivals with an increasing level of industrialization (Farzaneh, 2019;International Energy Agency, 2019). Such energy strategies including secure, affordable, and sustainable energy, reduction in the importation of fossil fuels, and energy efficiency are formulated and implemented to improve social welfare and quality of life for its citizens. The energyenvironment spectrum for regional energy efficiency in Asia is important for policymakers, researchers, and respective governments to understand which strategy would foster the reduction in environmental impacts of carbon emissions on economic activities, perhaps with the goals of low carbon emissions, and urbanization sustainability.
Our knowledge on energy consumption, urbanization, and economic growth and their effects on carbon emission in Asia is sorely lacking in terms of their interaction outcomes and possible policy implications both in the number of related studies and regional context (Maddock, 1995;Ng et al., 2016). Asian countries provide a unique phenomenon in which to better understand the long-run cointegration and causality between energy consumption, urbanization, carbon emission, and economic growth because of the increasing economic development and investment in green and renewable energy (Kannan et al., 2007;Balamurugan et al., 2012;Ali et al., 2017). Therefore, regional and geographical differences may call for heterogeneous investigation into the energyenvironment-growth phenomenon. Heterogeneous investigation into the relationship between carbon emissions, energy consumption, urbanization, and economic growth may be important for the formulation and implementation of dynamic carbon emission reduction policies across countries, suggesting that dynamic panel estimation modeling techniques may be appropriate in this context. Reliable and econometrically robust quantitative techniques are valid to analyze the heterogeneous individual effects and regional interconnectedness of energy-environment-growth nexus. Past findings on the relationship between carbon emissions, energy consumption, urbanization, and economic growth have been mixed. One major limitation is the used of first-generation panel unit roots and cointegration tests that assume that error terms are not correlated, therefore they fail to address the issue of cross-sectional dependence across countries. First-generation panel cointegration tests suffer from size distortions and fail to account for cross-sectional dependence, which often result in biased conclusions as to whether a long-run relationship exists between the variables (Bai and Kao, 2006). Andrews (2005) posit that common shocks (i.e. geopolitical risk and currency crisis) across countries may induce cross-sectional dependence among individual countries in the panel, in which the impacts vary across the cross-section units or countries. The increasing level of interconnectedness, financial integration, trade, and capital mobility accounts for cross section dependence across countries (Munir et al., 2020). In the presence of crosssectional dependence and slope heterogeneity in the series, the use of fixed effect, random effect, DOLS, and FMOLS estimators may give biased and inconsistent results (Murthy and Ketenci, 2020). Thus, it requires the application of second-generation panel unit root (i.e. cross-sectional panel unit root), cointegration, and dynamic common correlated effects panel modeling.
We contextualize our study in the Asian countries, longitudinally assessing the relationship between energy consumption, carbon emission, urbanization, and economic growth which presents an interesting setting whereby Asian countries face increasing economic development and investment in green and renewable energy. As such, the unique factors of geopolitical risk, war, economic growth, population, and economic growth programmes among the Asian countries may account for new evidence that differs from prior studies. Past studies find significant differences in the relationships between energy consumption, urbanization, carbon emissions, and economic growth in the MENA context (Al-Mulali et al., 2013;Farhani et al., 2013; compared to regions such as Sub-Saharan Africa (Gao and Zhang, 2014;Salahuddin et al., 2019), North America and G-20 countries Usman et al., 2020), Europe Kasman and Duman, 2015), and BRICS (Wang et al., 2016;Adedoyin et al., 2020). Regional context matters as a failure to address the regional interconnectedness that influences institutional and stakeholders' policies on urbanization, energy consumption, and carbon emissions may result in poor and ineffective energy solutions.
Accordingly, this study seeks to provide answers to the questions: Do long-run relationships exist between urbanization, energy consumption, carbon emissions, and economic growth in Asian countries? and to what extent do these variables interact in individual Asian countries?
Data from World Bank development indicators were employed. Findings reveal that a high average coefficient of 0.602 between carbon emission and energy consumption while low coefficients of 0.114 and 0.184 for the pairs of carbon emission-urbanization and carbon emission-GDP, respectively for the panel as a whole, suggesting a cointegration between carbon emission, urbanization, energy consumption, and economic growth. Our findings also reveal that there are considerable differences in the long-run relationships between urbanization, energy consumption, and carbon emissions across individual Asian countries. We link these relationship differences to the rise in rural-urban migration, concentration of industrial, commercial, and corporate industries and sectors in urban cities, population growth, and presence of geopolitical risks. These results depart from past studies on the relationships between urbanization, energy consumption, carbon emissions, and economic growth to income groups and distribution such as in the MENA context.
This study contributes to energy literature, first, by employing second generation cross-sectional Pesaran (2007), Westerlund (2007), and Pesaran's (2006) common correlated effects mean group estimation technique. These techniques are important to model the long-run relationships between variables and across panel observations that account for cross-sectional dependence which are robust to serial correlation and slope heterogeneity. Second, we compare our results with first generation estimation techniques i.e. DOLS and FMOLS to establish the differences between the firstand second-generation estimators. The positive and negative coefficients between carbon emission and GDP may suggest the phase or stage of EKC hypothesis across Asian countries, evidence shows that UAE, Bangladesh, China, India, Malaysia, Philippines and Singapore have negative carbon emissions-GDP coefficients, suggesting that these countries benefit from increased economic growth following carbon emissions reduction policies.
The remainder of the study is structured as follows. Section data methodology describes the model specification and panel data estimators. Section empirical results shows the estimation outputs of the panel data estimators. Discussion on findings and implications with conclusion are provided in Section discussions and conclusion.

Model Specification
This paper investigates the relationship between energy consumption, urbanization, carbon emissions and economic growth. This study utilizes the panel co-integration and dynamic common correlated effects mean group (CCEMG) techniques to achieve this objective. The panel data for this study consists of a sample period 2000-2014 of cross-country observations for 42 countries in Asia. The sample period is due to the availability of data. The data on the variables of this study are sourced from the World Development Indicators (WDI). The model controls for economic growth to avoid omitted variable bias (OVB) and to find evidence of an initial stage of the EKC hypothesis in the Asian context. The empirical model is as follows: Urbanization (URBAN) is expressed in terms of population. Carbon emissions (CO2EMS) are expressed in millions of metric tons. Energy consumption (ENERGYCON) is the kg of oil equivalent, measured in per $1,000 GDP (constant 2005 PPP). Economic growth (GDP) is expressed in real gross domestic product. Equations (1) to (3) were further written in their VECM and matrix structure as suggested in Pradhan and Arvin (2015) and Mahembe and Odhiambo (2019).
where CO2EMS, ENERGYCON, URBAN, and GDP which are individual dependent variables in the roles. Symbol Δ describes the first difference operator (I − L); i 1, 2, 3, . . . , N; t 1, 2, 3, . . . , T; α 1j to α 4j , c 11ik to c 43ik , and λ 1j to λ 4j are the system parameters to be estimated; ε 1it to ε 4it are the error terms; ECT it−1 represent the lagged values of the error correction terms from the system co-integration regressions. Symbols λ 1j to λ 4j also represent the speed of adjustment along the long-run equilibrium path. Short-run causality is usually inferred from the lagged dynamic variables of the independent variables (c 11ik to c 43ik ) by using partial χ 2 Wald test as suggested by Wald (1943). The long-run causality is tested through the lagged co-integrating vectors ECT t−j ( λ 1j to λ 4j ).

Panel Unit Roots
The first step in a panel cointegration test is to first examine the stationarity of the series. Several studies have used panel unit root to test for the stationarity of data series (Zhu et al., 2012;Al-Mmulali et al., 2013). The unit root examines the order of integration in the panel series. This study adopts four different first generation panel unit root tests, namely Levin et al. (2002) and Im et al. (2003), ADF-Fisher, and PP-Fisher. While LLC and IPS assume homogeneous unit root, ADF and PP-Fisher assume heterogeneous unit root (Baltagi, 2013). The IPS is advantageous because it caters for individual observation (i.e. country) heterogeneity (Mahembe and Odhiambo, 2019), however, we conduct three other tests as robustness unit root tests. The goal is to focus on unit root processes as a cointegrated process. We expect to see a stationary series at its first difference, d 1 i.e. I(1) process to enable us to perform the Pedroni cointegration test. For example, the principle of ADF-fisher is presented as follows: where d mt is the vector of the deterministic variables; ρ i is the lag-order and it is allowed to vary across observations which are determined by selecting a ρmax and thereafter using a t-stat of θ iL ; ε i, t is assumed not to correlate with the regressors and is permitted to be independently distributed across observations (i 1, 2, 3, . . . , N) and time-series (t 1, 2,3, . . . T). The results of the unit-roots are expected that the series must be stationary at their first difference [i.e. d I(1)], which will inform the panel causality test and meet the requirement of the panel VECM Granger causality tests. The null hypothesis is written as: H 0 : α i 0, for all the observations, i. The assumption is that if π i describes the p-value form, any unit test estimated for cross-section will be I. Choi (2001) shows that: where Φ −1 describes the inverse of the standard normal cumulative distribution function.

Panel Cross-Sectional Dependency (CSD) Tests and Determination of Optimal Lags
One of the estimation procedures in panel granger causality tests is to test for cross-sectional dependency. A shock in an observation (i.e. a country) may affect other observations in a sample due to the level of interconnectedness through culture and increased globalization (Muye and Muye, 2017;Mahembe and Odhiambo, 2019). There is considerable evidence that a shock to a country would affect other countries in Asia. For instance, the Asian currency crises of 1997-1998 (Ito and Hashimoto, 2005;Azad, 2009;Chang et al., 2009), the transnational transmission of the global financial crisis (Liu and Arunkumar, 2019), and the ASEAN-Indochina and Rohingya crises (MacLean, 2019; Weatherbee, 2019).
Although the Breusch-Pagan (1980) test for the panel cross-sectional dependency is valid for a panel with small N and large T (Pesaran, 2004), this study uses four different tests to confirm the consistency of CSD results; Breusch-Pagan LM, Pesaran scaled LM, Bias-corrected scaled LM, and Pesaran CD. The null hypothesis of CSD is that CSD does not exist in residuals and the test statistics are standard normal i.e. asymptotically distributed (Tekin, 2012). Optimal lag length is important in panel granger causality estimations (Tekin, 2012). This study uses the Schwarz information criterion (SIC) (Schwarz, 1978) method of lag length selection. This is because SIC performs better than other lag length selection methods in the VECM framework (Wang and Liu, 2006). The SIC method is better than the Akaike information criterion (AIC) (Akaike, 1974) because the latter gives higher lag lengths.

Panel Cointegration
The significance of the panel cointegration can be seen in its use when examining the long-run and short-run causality among variables in economics and finance studies due to its estimation power. A lack of cointegration would suggest that a panel VAR should be estimated for casualty since it indicates that a long-run relation does not exist (Granger, 1988). Several panel cointegration tests are used to test for long-run relationships such as Kao (1999) test and Pedroni (1999), Pedroni (2004) panel cointegration tests. The Kao test assumes a common-integrating vector while the Pedroni test caters for heterogeneity because it combines the residual-based Lagrange multiplier (LM) tests, ADF, and PP principles (Mahembe and Odhiambo, 2019). Thus, this study utilizes the Pedroni cointegration test. This study uses the Pedroni cointegration test to examine the bi-directional longrun relationship between energy consumption, urbanization, carbon emission, and economic growth. Meanwhile, if the bidirectional relationship between the four variables is cointegrated, then the Granger causality test using the vector error correction model (VECM) would be performed. The goal of the granger causality-VECM is to investigate the bi-directional causal relationships (both short and long run) between energy consumption, urbanization, carbon emission, and economic growth. One of the benefits of the Pedroni cointegration is that it allows for heterogeneous intercepts and trend coefficients across sections. For robustness, the study also performs the Kao cointegration test. The regression pattern of Pedroni cointegration is presented as follows: where y is the dependent variable, x is the independent variable, and x N is the number of explanatory variables in the regression model. The intercept or constant across each observation or cross-section of the model is α i while c 1i . . . c Ni are the coefficients or slopes of the regression model. The residual term is represented as e i,t in the model. The rule of thumb is that the null hypothesis of no cointegration for the residual term is that it is e i,t I(1) while the alternative hypothesis of cointegration is that e i,t I(0). The equation below presents the regression of each observation to establish whether the residuals are I(1).

Panel Causality and Panel Dynamic Estimation Tests
This study conducts a panel data causality after establishing a panel cointegration. Evidence has shown that the panel cointegration is a preconditioning test for dynamic panel causality test (Engle and Granger, 1987;Granger, 1988;Stock and Watson, 1993). The presence of panel cointegration infers the estimation of a panel VECM granger causality, otherwise, no cointegration will lead to the estimation of a panel VAR causality. While a panel VECM Granger causality test will estimate both the long-run and shortrun coefficients, a panel VAR causality test estimates only the short-run coefficients. This is because the panel VAR equation does not include the error correction component that addresses the long-run causality. A panel VECM can be estimated using the Granger causality test. The long run, short run, and error correction terms for the panel and individual countries can also be inferred from the pooled mean group estimation. This study estimates the first-generation panel dynamic Fully modified OLS (FMOLS) and dynamic OLS (DOLS) should the models be cointegrated. DOLS models are used to establish either negative or positive long-run relationships between variables in the cointegration system. DOLS models required the use of lags and leads that results in a cointegration equation error term (Pedroni, 2001). The cointegration equation error term should be orthogonal to the entire history of the stochastics regressor innovations (Al-Mulali et al., 2013). By adding the lags and leads of the different regressors, it will eliminate the long run correlations between the residuals. The limitation of the DOLS is that it does not produce efficient estimates like the Fully Modified Ordinary Least Square (FMOLS) (Phillips and Hansen, 1990). Phillips and Hansen (1990) initiated the FMOLS to deal with the problems of asymptotic bias and nuisance parameter dependency associated with cointegrating vector estimates in the conventional single equation case but it is not robust to endogeneity bias and serial correlation in dynamic panels. Pedroni (2001) developed the heterogeneous FMOLS estimator to adjust for endogeneity bias and serial correlation. FMOLS is appropriate and suitable for panel heterogeneous cointegration. Table 1 shows the descriptive statistics of the variables under review. It is to ascertain the basic measure of central tendency and dispersion of the variables for the period 2000-2014. As shown in Table 1, the minimum value of carbon dioxide is 19,546 kt, with a maximum value of 29,934 kt over the period 2000-2014. Real gross domestic product has a minimum value of $19,866 and a maximum value of $29,980. The variables are positively skewed. Energy consumption and GDP are not normally distributed as the null hypothesis of normality is rejected. Carbon emissions and urbanization are normally distributed, supporting the summary statistic reported on carbon dioxide by Bekun et al. (2019). The total sample size for the panel is 630 observations. Table 2 shows the correlation matrix analysis to examine the association between carbon emissions, urbanization, energy consumption and economic growth. Carbon emission has a positive association with economic growth, as expected. Urbanization also has positive association with economic growth. However, the correlation coefficient for CO2EMS-GDP is higher than the correlation coefficient for Urban-GDP. This supports the EKC hypothesis that at the first phase of the theory, countries benefit from increased real income through environmental degradation. Energy consumption has negative association on economic growth. The results also show that for every kt of carbon dioxide emitted, urbanization contributes to about 79.5% while energy consumption contributes about 5.5%. Since correlation matrix analysis cannot explain long-run relationships between the variables, Table 3 presents the panel unit root tests as the first step to cointegration analysis. Table 3 shows the results of the first-generation panel unit root tests. The series is required to be stationary at their first difference in the phases of panel Granger causality estimation. The absence of non-stationarity in the series may infer a spurious regression as suggested by Gujarati and Porter (2009) and Baltagi (2013). Four different panel unit root tests were tested, and results presented with the null hypothesis that there is no stationarity in the variable (i.e. the variable contains a unit root). Table 3 includes t-statistics at intercept and intercept and trend. This study considers the presence of stationarity at the 1% level of significance.  Estimating the first-generation panel unit root tests, the results demonstrate that CO2EMS is not significant at intercept for LLC, IPS, ADF, and PP at levels. At levels, ENERGYCON also depicts significance for all the four panel unit roots except IPS. URBAN is not significant at levels for LLC, IPS, and PP both at intercept and with intercept and trend. GDP is also insignificant at the level but significant at the 1% level and consistent at its first difference. However, at the first differences of the variables, there is considerable acceptance for the series being stationary at their first differences across the four panel unit-roots. As compared to their levels using IPS where all the variables are insignificantly mixed both at intercept and with intercept and trend, the results reveal that the variables are all significant at their first difference. IPS unit root test seems to have better results than other panel unit root tests. This is because IPS caters to individual observation (i.e. country) heterogeneity (Im et al., 2003). Thus, the results depict that the four panel variables are stationary at their first difference i.e. I(1).

Empirical Results
The Pedroni (2004) cointegration test was adopted, with the null hypothesis that there is no cointegration. The results of the panel cointegration test are shown in Table 4. This is to confirm whether a long-run equilibrium relationship exists between the variables. In each cointegrating regression, one variable is the dependent variable while the other variables are the independent and vice versa.
The results for the first generation Pedroni cointegration test show that six out of the eleven Pedroni statistics reject the null hypothesis of a no cointegration for the four models (CO2EMS,  ENERGYCON, URBAN, and GDP). The results for the Pedroni cointegration test were also robust to the Kao (1999) cointegration test. Therefore, there is evidence of a long-run equilibrium between the four variables. Table 5 shows the results for dynamic OLS (DOLS) test. The significance of the relationships between urbanization, energy consumption, carbon emission, and economic growth varies across countries based on different reasons. This may be due to the level of development in the countries, the level of investments in renewable and non-renewable energies. It can be based on the country's law on urbanization, population, and energy consumption bills. There is a long-run relationship between energy consumption, urbanization, and carbon emissions. There is considerable evidence of a positive long-run relationship between energy consumption and carbon emissions in Afghanistan, Kuwait, Myanmar, Mongolia, Thailand, Tajikistan, and the Yemen Republic. This may be due to the long-term geopolitical risk, war, and other incessant regional crises in these countries. Evidence has it that incessant war contributes to carbon emission (Oreskes and Conway, 2008;Gingrich et al., 2011;Zou et al., 2016). However, energy consumption shows a negative long-run relationship with urbanization in countries like China, India, Bangladesh, Bahrain, Georgia, Indonesia, Saudi Arabia, Singapore, and Turkey. For example, In Singapore, the investment in renewable solar energy and other distributed renewable sources has reduced greenhouse gas emissions considerably, resulting in a growing urban-tourist destination (Kannan et al., 2007; Balamurugan et al., 2012;Ali et al., 2017). On the other hand, the majority of the Asian countries show a significant long-run relationship between urbanization and carbon emissions. For the results of the dependent variable, URBAN, carbon emissions show a significant positive long-run impact on urbanization across the majority of the Asian countries. At the 1% level of significance, the positive long-run impact of urbanization on carbon emissions is severe in some of the populous Asian countries including India, Indonesia, Pakistan, and Bangladesh with coefficients of 7.042, 10.634, 7.816, and 8.637, respectively. Concerning the positive long-run impact of carbon emissions on urbanization, the effect is highest in Qatar, UAE, Saudi Arabia, and the Maldives with coefficients of 1.552, 1.440, 0.829, and 0.616, respectively at the 1% level of significance. Although the DOLS estimator has been used in past studies, however, due to its weakness when compared with FMOLS estimator, the study also estimates the FMOLS and the results are presented in Table 6. Table 6 presents the results of the long-run relationships between carbon emission, energy consumption, urbanization, and economic growth, with each as a dependent variable based on individual countries and panel FMOLS tests. As shown in Table 6, the study finds positive long-run effects of energy consumption, urbanization, and GDP on carbon emissions for the 42 Asian countries-group (panel). Results also reveal that energy consumption reduces urbanization and GDP in the long run. Meanwhile, GDP and energy consumption have negative bi-directional long-run relationships. This may be due to carbon taxation and appropriate energy pricing. The negative effect of energy consumption in reducing urbanization may also be due to a feedback loop for urban development, and countries may not benefit from increased productivity, which is in support of Munir et al. (2020). It is an indication that more and better policies to improve environmental quality in such a way that countries can benefit from increasing urbanization without affecting is the amount of urban carbon emission is needed. The panel FMOLS test revealed that there is long-run elasticity of Urbanization to carbon emissions and from carbon emissions to Urbanization but the coefficient for Urbanization to carbon emissions is stronger than the coefficient for carbon emissions to Urbanization. The results of the FMOLS are more consistent compared with DOLS. The panel FMOLS results indicate that a 1% increase in urbanization results in a 0.025% decrease in energy consumption, suggesting that the alternative energy use, adoptions of solar panels in cities, transport planning, and compact urban form are yielding better results in rapidly urbanizing Asia.
The results show that energy consumption, urbanization, and GDP contribute positively to carbon emissions in the long run to highly populated countries in Asia including China, Bangladesh, India, Iran, Pakistan, The Philippines, and Turkey. However, Indonesia's economic growth (GDP) will benefit positively from the combined effect of urbanization, energy consumption, and carbon emissions. This may be due to Indonesia's green bond and green Sukuk initiative (a shariah-compliant bond) where 100% of the proceeds exclusive go to finance or refinance green projects that contribute to mitigating and adapting climate change as well as preserving biodiversity. In Singapore, while urbanization has a positive long-run relationship with carbon emissions at a 9.87% change, it, however, increases GDP at the 4.117% change, suggesting that the country should pursue and maintain a balance in urban-carbon emission ratio policy that will reduce carbon emissions and at the same time increase GDP. The results for the Panel showed that a percentage increase in economic growth (GDP) increases carbon emissions and urbanization by 0.217% and 0.043%, respectively, but it reduces energy consumption by 0.183%. Asian countries' economic growth benefits more from carbon emissions than GDP contributes to the amount of carbon emissions. Eight out of the 42 Asian countries (Armenia, Bangladesh, Bhutan, China, Iran, Kuwait, The Philippines, and Turkey) experience long-run negative effects of urbanization and energy consumption on economic growth while carbon emissions contribute to economic growth.
After estimating the DOLS and FMOLS, we proceed to estimate the panel cross-sectional dependence (CD) to ascertain whether there is cross-sectional dependence among the sampled countries and the need to adjust for regional interconnectedness. Table 7 presents the results for CD test.
Using four different methods to test for cross-sectional dependence, the results show that there is evidence of correlations in the variables at the statistics are all significant at the 1% level of significance. The results are expected since we expected that increasing Urban growth causes an increase in carbon emissions and Urban also increases energy consumption, and at the same time increasing energy consumption contributes to an increase in carbon emissions, and in turn impacts economic growth. Our panel CD results are in support of Munir et al. (2020) that also find the presence of cross-sectional dependence among ASEAN countries, which are among the Asian panel used in this paper. The increasing development in infrastructure and high level of tourists in Asia countries may also facilitate the relationships between the variables. The presence of integration and stationarity leads to the estimation of the second-generation panel cointegration test. First, we need to estimate a second-generation panel unit root using the Pesaran's CIPS panel unit root test. Second, estimate the  Table 8. Table 8 shows the Pesaran's panel unit root using the CIPS test. The results show that the null hypothesis of a unit root is not rejected at the 1% level of significance, suggesting that the variables/series are integrated at their first difference, i.e. I(1). Thus, it requires the computation of the Westerlund (2007) as shown in Table 9 below.
The study used the bootstrapped approach of Westerlund (2007) to allow for cross-sectional dependence and report bootstrapped p-values. As shown in Table 9, the Westerlund (2007) cointegration test revealed mixed findings, a few results fail to reject the null hypothesis of no cointegration between carbon emissions and energy consumption due to insignificant statistics. Only the statistic of Pa is significant for the cointegration between carbon emissions and GDP. Both the statistics of Gt and Pt are significant for the cointegration between carbon emissions and urbanization and between energy consumption and urbanization. We restricted the lags and leads to 1 due to our 15 years or 15 observations. Stata command 'xtwest' requires at least 21 observations for 3 lags and 3 leads. Notes: H 0 : homogeneous non-stationary; lag criterion decision: general to particular based on F joint test; critical values, CIPS with constant: 10% (-2.03), 5% (-2.11), 1% (-2.25); critical values CIPS with constant and trend: 10% (-2.54), 5% (2.62), 1% (-2.76). *indicates significance at the 1% level.  Meanwhile, the statistics of Gt, Pa, and Pt are significant for the cointegration between energy consumption and GDP. The null hypothesis of no cointegration between urbanization and GDP is rejected as the statistics for Gt, Ga, Pt, and Pa are all significant. Substantially, there is evidence of cointegration between carbon emissions, urbanization, energy consumption, and GDP, suggesting that there is a long-run equilibrium relationship between the variables when we consider cross-dependence in the panel except for the long-run equilibrium relationship between carbon emissions and energy consumption. In conclusion, the results of the bootstrapped Westerlund (2007) test provide consistent and strong evidence of cointegration relationship between the pairs of Urban-GDP and Energycon-GDP; moderate evidence of cointegration relationship between the pairs of CO2EMS-Urban and Energycon-Urban, and weak evidence of cointegration relationship between CO2EMS and GDP. However, it finds no evidence of a cointegration relationship between CO2EMS and Energycon. As stated, this study uses the CCEMG estimator to measure the relationship between carbon emissions, energy consumption, urbanization, and GDP in the panel and individual countries. Table 10 presents the empirical results of the panel estimators that take into consideration common factors, cross-sectional dependence, and slope heterogeneity. The carbon emissions-Urban coefficient is 0.114 and its significant difference from 1.000 indicates a relatively high degree of urbanization contributing to carbon emissions in the short-run among the panel countries. This suggests that in the long-run, the problem of urbanization is addressed, as the panel countries formulate and implement carbon reduction goals.
Concerning the EKC hypothesis of Grossman and Krueger (1995), the CCEMG results showed that there is evidence of the initial stage of the EKC hypothesis i.e. higher GDP does results in higher carbon emissions in Indonesia, Japan, Kuwait, Pakistan, Qatar, Saudi Arabia, Thailand, and Turkey. The formulation and implementation of energy and carbon emission reduction policies are still on-going in these countries, with slow improvement in economic development, a cleaner environment, and better life quality.

Robustness Test Using Pooled Mean Group-ARDL (PMG-ARDL)
This study further uses the Panel Pooled Mean Group-Autoregressive Distributive Lag (PMG-ARDL) estimator to establish the magnitude of cointegration as shown in Table 11.
As shown in Table 11, long-run relationship is inferred from the ECT coefficient values, confirming the presence of convergence speeds. The results depict that there is a long-run relationship between carbon emissions, energy consumption, urbanization, and economic growth. Our results are in line with Al-Mulali et al. (2013) who found long run relationships between carbon emissions, energy consumption, and urbanization. The long-run coefficient is the same across all the groups that make up the panel. In the long run, urbanization, energy consumption, and economic growth (GDP) granger cause carbon emissions at the 1% level. For the urbanization model, carbon emissions has significant relationship with economic growth positively while energy consumption and urbanization are negatively related in the long run at the 1% level. The error correction term (ECT) shows that there is cointegration among the variables across the different models in the panel and long run relationship. Thus, there is a long-run relationship at the 1% level for the carbon emissions model, energy consumption    Table 12.
Using the PMG-ARDL, the results for short run relationship in group panel show that energy consumption and carbon emission positively significantly related in the short run with a coefficient of 0.5014 at the 1% level. In relation to Urban model, energy consumption has a significant negative relationship with urbanization in the short run at the rate of 37.81%. For energy consumption model, both carbon emission and GDP have positive and negative relationships respectively, on energy consumption in the short run. Urbanization reduces GDP in the short run at the panel level. Although PMG estimator assumes common long run relationships between variables, however, the short run dynamics vary across countries. So, relying on grouped or panel short run estimates as shown in Table 12 may be misleading. Table 13 addresses individual countries' short run relationships. Table 13 shows the short-run relationship between carbon emissions, urbanization, energy consumption, and economic growth for individual countries using the PMG-ARDL estimator. It also shows that each country has a different error correction term, indicating the extent to which deviations from the long-run is corrected at the individual country level. The results show unidirectional relationship from GDP to energy consumption for China, UAE, and Russia. The unidirectional relationship is found from GDP to carbon emissions for Pakistan, Qatar, Singapore, Israel, Lebanon, and the Maldives. Unidirectional relationship is found from energy consumption to carbon emission in India, Kazakhstan, Myanmar, and the Philippines. However, bi-directional relationship (i.e. feedback hypothesis) is found from carbon emission to GDP in Russia, Georgia, Turkmenistan, Afghanistan, Armenia, and Uzbekistan in the short run.

DISCUSSIONS AND CONCLUSION
The upsurge in economic development, income level, and perhaps the geopolitical risk and war situations have influenced the link between energy consumption, urbanization, carbon emissions, and economic growth. The relationships between these variables informed the objective of this study. This study utilizes the second generation cross-sectional Pesaran (2007), Westerlund (2007), and Pesaran's (2006) common correlated effects mean group estimation technique to examine the relationships between energy consumption, urbanization, carbon emissions, and economic growth in 42 Asian countries for the period 2000-2014. The Pedroni cointegration estimation results showed that urbanization and carbon emissions were cointegrated but the no cointegration when the energy consumption model is considered. Results for CCEMG estimator showed that there are long run relationships between urbanization, carbon emissions, energy consumption, and economic growth. The results show that carbon emissions play a significant role in increasing urbanization for countries experiencing civil and political war such as Afghanistan, Iraq, Turkey, and Vietnam. It suggests that war in these countries are mostly in the rural, boundaries and resources areas, influencing rural-urban migration, and thus, contributing to the high rate of urbanization. The level of energy consumption also increases the rate of urbanization in countries like the United Arab Emirates, the Russian Federation, and Saudi Arabia. It indicates that most industrial, commercial, and household use of energy are highly concentrated in the urban areas because as governments of these countries reduce unemployment, seekers of job opportunities focus on industrial, commercial, and corporate sectors. Thus, increasing urban employment cycle and destination change. China along with Bahrain, Indian, Malaysia, The Philippines, Saudi Arabia, Singapore, Thailand, and Vietnam contributes to carbon emissions through their number of urban destinations. Policy-wise, the government should shift investments to rural destinations to reduce rural-urban migration. Although for some high-income countries, like China, Malaysia, and Saudi Arabia where a large percentage of their investments and revenue generation is concentrated in commercial cities, shifting investment to local communities may discourage foreign investors. However, continued investments into the conservation and efficiency of energy may reduce future carbon emissions but the ruraltourism destination programs must be fostered to reduce further urban migration to sustain the energy efficiency programs for reduced carbon emissions while at the same time not influencing a further increase in urbanization.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.