key: cord-0851755-ernebpj6 authors: Majeed, Muhammad Tariq; Tauqir, Aisha; Mazhar, Maria; Samreen, Isma title: Asymmetric effects of energy consumption and economic growth on ecological footprint: new evidence from Pakistan date: 2021-02-26 journal: Environ Sci Pollut Res Int DOI: 10.1007/s11356-021-13130-2 sha: cd4e3efaa000ffc2a32ba00956ffd99f19dab040 doc_id: 851755 cord_uid: ernebpj6 This study explores the asymmetric effects of both aggregate and disaggregate forms of energy consumption along with economic growth on environmental quality for Pakistan covering the period from 1971 to 2014. We have employed unit root test with breaks for stationary checks, BDS test for nonlinearity check and nonlinear autoregressive distributed lag (NARDL) approach for assessing the asymmetric co-integrating relationships among the variables by decomposing them into positive and negative shocks. The empirical findings for aggregate consumption reveal that only negative shocks have a significant impact on ecological footprint. Similarly, different sources of energy consumption have diverse asymmetric effects on ecological footprint. The positive (negative) shocks to oil and gas consumption increase (decrease) ecological footprint. Thus, an increase in oil consumption has a deteriorating impact on environmental quality while a decrease in gas consumption has a favorable impact on environmental quality. The asymmetric relationships also hold between coal consumption, electricity consumption, and ecological footprint. The positive shocks to coal and electricity consumption are negatively related with environmental quality while negative shocks are positively related with environmental quality. In the recent decades, human living standards and well-being have been improved to a greater extent because of high economic growth and prosperity. However, this growth has been compromising environmental quality as it heavily relies on energy-intensive inputs. Energy is a fundamental need for humans' survival because without heat, light, and power, firms and industries cannot be built, cities cannot survive, and resultantly goods cannot be produced. This ensures that energy being the lifeblood of the global economy acts as a main driving force behind all economic activities (Alam 2006) . According to IEA (2019), a tremendous rise in energy demand is observed across the globe. In 2018, the rate of global energy demand is recorded at about 2.3 % which was the ever fastest rise in the energy demand in the present decade. In particular, its major effects on the economies are twofold. First, it contributes to the economy through extraction, transformation, and distribution of energy goods and services by creating jobs and boosting economic growth. Second, it acts as an input and can affect the economy through energyrelated shocks (WEF 2012) and then to growth sustainability through a deteriorating environmental quality (Zhao et al. 2013; . Predominantly, energy comes from non-renewable resources (coal, oil, natural gas, nuclear) and has incredibly strong negative effects on the environment by producing carbon dioxide (CO 2 ), methane (CH 4 ) emissions, and other greenhouse gases (GHGs) in the atmosphere (Majeed and Luni 2019) . In Pakistan, like all other developing countries, the environmental issues have become the major concern of the Responsible Editor: Ilhan Ozturk government. The government of Pakistan (GoP) (2020) states that Pakistan is among the top ten countries which are highly vulnerable to climate change in the past 20 years. Eckstein et al. (2019) also inform that cost of environmental degradation is considerable in Pakistan as the country has lost 0.53 % per unit GDP, and economic losses of about 3792.53 million US dollars along with experiencing 152 extreme weather events from 1999 to 2018. Meanwhile, the country also witnessed a rise in energy demand from the last few decades as it is intently related to the development of the various sectors in the economy. According to the World Data Atlas (WDA 2020), Pakistan's primary energy consumption increased from 1.74 quadrillion btu in 1998 to 3.37 quadrillion btu in 2017, increasing at an average annual rate of 3.63%. The energy consumption in terms of oil equivalent stood at 460.23 million tons in 2014 as compared to 446.01 million tons in 2001. Along with this, the country also witnessed a severe energy crisis which has increased its reliance on imported oil. Therefore, the cost structure within the economy's power generation sector increased, putting the country into a severe shortage of gas and electricity. To address these issues, the government has formulated energy plans to provide indigenous, affordable, and sustainable energy for all. Consequently, Pakistan's dependence on some non-renewable (renewable) energy sources in the overall energy mix is slightly declining (increasing) from the past few years. In 2018, in terms of the energy mix, Pakistan's reliance on oil shrunk to 3.2 %, on gas 34.6 % (due to low consumption in the transport sector), and hydel energy to 7.7 % (due to short-sightedness of policy), respectively. The share of renewable energy remains 1.1 % in 2018 while the nuclear share increases about 2.7 % (GoP 2019). The overall petroleum product consumption remains about 19.68 million tons per annum. The consumption of oil and gas increased by 1.25952 million tons and 4.727285 million tons of oil equivalent in 2014 as compared to 18.461 million tons and 16.2979 million tons of oil equivalent in 2001, respectively. The significant rise in coal and electricity consumption has also been noticed in the last decade. In 2014 the coal and electricity consumption raised by 7.061729 million tons of oil equivalent and 420.5071167 kW h −1 per capita from 2.105 million tons of oil equivalent and 345.8106331 kW h −1 per capita in 2001, respectively (Bp global, 2020) . Along with this high energy use, the increase in carbon emissions is not surprising. The statistics show that total CO 2 emissions from energy consumption increased by 147.620 million metric tons in 2014 from 103.160 million metric tons in 2001. The links between energy consumption and environmental quality are also found in the existing literature Ozcan et al. 2020; Shujah-ur-Rahman et al. 2019 ). According to and Hassan et al. (2019) , Pakistan's economy is still in the development process where day by day increasing industrial activities are aggravating the existing energy demand. This demand is met by both renewable and non-renewable energy sources. The large share of energy comes from conventional energy sources such as coal, gas, and oil which contribute to high carbon emissions in the atmosphere. The increasing number of vehicles, transportation modes along with zero high-quality fuel, further exacerbates the environmental quality issues (Danish and Baloch 2017; Danish et al. 2018) . However, renewable energy has helped to control carbon emissions to a certain extent as claimed by Danish and Wang (2017) and Wang et al. (2018) in their recent studies for Pakistan. These studies have used linear econometric approaches and rely on a single measure of environmental CO 2 emission. Further, in the overall literature, there are only two studies that use the nonlinear econometric approach to explore the relationship between energy and the environment Baz et al. 2020) . The study by Baz et al. (2020) used an asymmetric approach but restricted to employ only one measure of energy consumption, while used disaggregated energy data but restricted to the use of standard econometric techniques (time series) and environmental indicator (CO 2 ). Hence, this study contributes to the existing literature in many ways. First, we provide the nonlinear impacts of energy sources on environmental quality by employing the nonlinear autoregressive distributed lag (NARDL) method. Second, the study incorporates both aggregated and disaggregated energy consumption data. The aggregated data comprises total energy consumption while disaggregated energy data includes oil, coal, gas, and electricity consumption. Third, we proxied the environmental quality by ecological footprint (EF) which is considered the best measure of environmental sustainability. It is a comprehensive environmental indicator capturing various dimensions of the ecosystem namely built-up land, carbon, cropland, fishing grounds, grazing land, and forest products (Bagliani et al. 2008; Majeed and Mazhar 2019; Danish et al. 2019) . Last, we use the maximum accessible available data covering the period from 1971 to 2014. In the methodological setting, the study intends to use various unit root tests (ADF, ADF breakpoint, and PP), for checking nonlinearity (BDS test) and recently popularized estimation technique (NARDL) to find out the nonlinear energy consumption impact on ecological footprint. The findings of this study will help in numerous ways. The overall findings will provide important insights on energy and the environment nexus. They will also help to recognize the type of relationship between various energy sources and ecological footprint. By providing the analysis of nonlinear relationships, the findings will assist policymakers and practitioners in formulating the energy-mix policies, because policymakers need literature, clear evidence, and quality research for formulating their policies. Further, the findings also provide assistance in policymaking regarding switching of conventional energy sources to clean energy sources by identifying the most harmful energy source. Through this, policymakers can get the idea of what type of policies will work for the gradual shift of nonrenewable energy sources to renewable energy sources to protect the environment. The structure of the remaining paper is as follows: Section 2 discusses the related literature, Section 3 describes the data set, model and methodology, Section 4 presents the empirical findings, and the last section concludes the study. The research on energy growth the environment has becomes the hot topic of today's world. The environment has a greater impact on the economy. In fact, economic system influences the environmental quality on a larger scale through excessive natural resource extraction and energy consumption. Therefore, the economic literature is abundant with "economic growth-environment" and "economic growth-energy consumption-environment" nexus. This section is segmented into the following sub-sections. The relationship between economic growth and environmental quality is puzzling and remains a matter of debate among various researchers and policymakers. Economic theory suggests that there is a positive relationship between economic growth and environmental degradation. As in the initial stages, an increase in economic activities increases the use of energy and fossil fuels in the absence of clean technology that in turn harms environmental health. However, this happens up to a certain point, and after that, economic development contributes to environmental sustainability by changing the economic structure of the economies. This renowned relationship between economic growth and the environment is known as the environmental Kuznets curve (EKC) hypothesis. Besides these theoretical predictions, the world academia is motivated to test this relationship empirically. For this purpose, the studies are conducted using cross-sectional, time series and panel data set and have used diverse methodologies and environmental indicators. Earlier studies, Shafik and Bandyopadhyay (1992) by using total and annual deforestation and Panayotou (1993) by using deforestation, sulfur dioxide (SO 2 ), nitrogenous oxides (NO x ), and solid particulate matter (SPM) for industrial-and energy-related pollutants tested the environment-national income relationship. They found that the income's linear term is positive and concluded a positive relationship between environmental indicators and income. According to them, an increase in economic development increases the use of energy and material inputs that causes more waste and emissions. Alongside this, growth fosters the process of structural changes in the economy, from agriculture to industrial and industrial to service, making people more aware of the environmental quality and requirement of green technology. Thus, after a threshold, growth promotes environmental sustainability. Selden and Song (1994) consider SO 2 , NO x , particulates, and carbon monoxide as environmental proxies while estimating the EKC hypothesis. They concluded the presence of an inverted U-shaped curve. Grossman and Krueger (1995) use more reliable data and several dimensions of water and air pollutants and validated the inverted U-shaped EKC. According to them, very poor countries are not in a position to manage environmental pressure with the increase in development level in the initial stages. However, after some time with the structural transformation, growth-environment relationship turns favorable. Akbostancı et al. (2009) by employing time series cointegration and panel estimated generalized least squares (EGLS) find a monotonically increasing relationship and Nshaped EKC between carbon emissions and GDP and particular matter, SO 2 and GDP, respectively. Fodha and Zaghdoud (2010) also obtained the monotonically increasing relationship between CO 2 and GDP for the economy of Tunisia. The EKC, however, was validated by the study of Fodha and Zaghdoud (2010) for SO 2 and GDP. Further, Narayan and Narayan (2010) and Jaunky (2011) use panel estimation techniques for 43 developing countries and 36 high-income countries, respectively. They find a positive relationship between income per capita and carbon emissions for most of the Middle Eastern and African countries and all of the South Asian, East Asian, and Latin American countries, and for the group of high-income countries, respectively. Piaggio et al. (2017) utilize the data from 1882 to 2010 for the Uruguay economy and argued that with the increase in economic activities, carbon emissions increase at a decreasing rate. Further, the recent studies by Uddin et al. (2017) and Ulucak and Bilgili (2018) used ecological footprint (EFP) as an environmental indicator. Uddin et al. (2017) employed the dynamic OLS on the panel data set of 27 highest emitting countries and concluded the positive impact of real income on EFP. Similarly, Ulucak and Bilgili (2018) confirm the EKC hypothesis for upper, lower, and middle-income countries while employing continuously updated fully modified and bias-corrected techniques. Recently, reinvestigate the EKC using a longer time horizon more than half century for different income groups and provide mixed evidence. They conclude that evidence on EKC is inconclusive based on methodology, sample data, and country heterogeneity. Economic growth, energy consumption, and environment nexus Economic development is a process of structural transformation that increases production activities and in turn the use of energy in the respective economy. Therefore, growth is highly linked with energy use which is exceedingly connected with the level of emissions. Numerous studies are conducted on this topic and present mixed results based on different econometric methods, sample size, and proxies used to measure environmental quality. Pao and Tsai (2010) and Pao and Tsai (2011) by employing panel co-integration techniques on BRIC economies data set have concluded the positive impact of energy consumption on CO 2 emissions and also validated the EKC hypothesis (except Russia). In another study for Brazil, Pao and Tsai (2011) used the gray prediction model and panel cointegration techniques to explore the relationship among GDP, energy consumption and carbon emissions. They argue that energy is the most important determinant of CO 2 emissions and have an inverted U-shaped relationship with income. As initially, energy consumption increases with the increase in income and then stabilizes and falls after some time. Wang et al. (2011) by taking the panel data of 1995 to 2007 for 28 Chinese's provinces and Shahbaz et al. (2012) by using the time series data of Pakistan's economy over the period 1971-2009 find the long-run bidirectional causal relationship among energy consumption, GDP and CO 2 emissions, and one-way causality from energy consumption to CO 2 emissions, respectively. In contrast, Saboori and Sulaiman (2013) use the Malaysian economic data of aggregate and disaggregate energy resources, namely, gas, oil, coal, and electricity. Their findings suggest the existence of an EKC hypothesis for the analysis based on disaggregated energy data. Ozturk and Acaravci (2013) also studied the relationship between energy consumption and carbon emissions along with economic growth. By taking the data from 1960 to 2007 and employing the ARDL technique, they confirm the positive relationship between energy consumption and CO 2 emission within the framework of the EKC hypothesis for Turkey. In a similar context, Ahmad et al. (2016) conducted the study by utilizing the data from 1971 to 2014 for the Indian economy. His results confirmed the positive impact of energy consumption (aggregate and disaggregate) on CO 2 emissions along with validating the EKC curve. According to him, energy consumption influence carbon emissions through the feedback effect. Ozturk et al. (2016) find different impacts depending on the income group of the countries. According to them, energy consumption is positively related to ecological footprint for most of the countries in all income groups. The EKC, however, is validated in most of the countries lying in the upper-middle and upper-income groups. Esso and Keho (2016) and Bekhet et al. (2017) elaborated this relationship for 12 sub-Saharan African countries and GCC countries and obtained similar findings. Ozturk (2017) proves fossil fuel energy consumption as a risk for growth sustainability. His findings assure that nuclear energy has the power to control pollutant emissions along with improving economic growth. Further, Charfeddine and Mrabet (2017) by using a fully modified least squares (FMOLS) estimation technique have also explored the unfavorable impact of energy consumption on EF. In another study, Charfeddine (2017) applied Markov switching equilibrium correction model and obtained mixed findings. His findings support the unfavorable impact of electricity consumption on EF and the favorable impact on carbon footprint and carbon emissions. Solarin et al. (2017) along with confirming the EKC hypothesis found out somehow surprising effects of renewable energy consumption. According to them, the impact of both fossil fuel and renewable energy on carbon emissions is positive for the economy of Ghana, while Danish and Wang (2017) found a negative (positive) association among renewable (non-renewable) energy consumption and carbon emissions. They also validate the EKC hypothesis for Pakistan's economy. Other studies such as Mrabet and Alsamara (2017) using ARDL, Atasoy (2017) Focusing on ecological management, He et al. (2018a) analyze ecological vulnerability by construing an index of vulnerability. For this purpose, they used spatial analysis of Geographic Information System (GIS) method and multicriteria decision analysis (MCDA). Their findings show that a larger population of China is living in an extremely ecologically valuable area where vulnerability is mainly caused by increasing economic activities (involving high energy consumption), natural resource consumption, inadequate environmental protection measures, and population pressure. In another study for Marcellus Shale in Pennsylvania and West Virginia, He et al. (2018b) emphasize on the integration of water management and greenhouse gas emission within the shale gas supply chain design and operations. Their analysis concludes that water discharge and water consumption both have a minute impact on economic efficiency, shale gas production, electric power generation, and overall greenhouse gas emissions. Emphasizing environmental issues from another perspective, Zhang et al. (2018) look the role of e-commerce in reducing electronic waste (e-waste) in China. Using the primary data of 896 Chinese (through a questionnaire) and employing the order logit regression, they conclude that e-commerce has not obtained universal acceptance yet because price is involved in recycling plus its acceptance is mostly limited to the customers involved in online shopping. However, ecommerce is an economic option for controlling e-waste that requires better policy options. Further, using the ARDL model by Danish et al. (2018) and the NARDL model by Sohail et al. (2021) explained that carbon emissions can also be controlled in Pakistan by limiting the energy consumption in the transport sector as their findings confirm the positive connection between transport energy consumption and CO 2 emissions. Lv et al. (2019) is of the same view as his analysis for China proves the emissions increasing impact of railway freight, waterway freight, road freight, and aviation freight. The impact of road freight is high compared to other modes. Moreover prove the EKC hypothesis and positive relationship among energy consumption and carbon emissions for Pakistan and ASEAN-5 countries, respectively. Zhao et al. (2020) consider the official driving cycle based on the fossil-fuel vehicle (FV) as a cause of environmental pollution in China comparative to the driving cycle based on electric vehicles (EV). They argue that EVs are more environmentally friendly as they don't release any pollutant emissions. Sarkodie et al. (2020) utilize the Prais-Winsten transformed regression with robust standard errors and dynamic ARDL simulations for China over the period 1961 to 2016. According to them, renewable energy consumption helps to control carbon emissions and EF while fossil fuel energy consumption exacerbates both environmental indicators. Supporting this, Wang et al. (2020) argue that replacing fossil fuel with biomass energy can help to control environmental issues. In another recent study, Shaari et al. (2020) find the positive impact of natural gas and energy consumption on carbon emissions of the Organization of Islamic Cooperation (OIC). Likewise, Majeed and Tauqir (2020) using extended STIRPAT modeling advocated a rise in carbon emissions as the result of energy consumption irrespective of different income categorizations of the sampled countries. In a similar setting, and Baz et al. (2020) used the methodology of asymmetric or nonlinear ARDL while exploring the relationship between energy consumption, economic growth, and the environment. By doing so, they confirmed the asymmetric relationship among the variables. So, following them, the present study employed nonlinear ARDL (NARDL) for the empirical analysis. The contribution of this study is that it employs NARDL for Pakistan's economy data using both aggregated and disaggregated energy sources taking EF as an environmental indicator. Although one similar study exists, its focus is on total energy consumption rather than types of energy (Baz et al. 2020) . Further, the ones who employed disaggregated data have used time-series data and CO 2 emissions ). In addition, other few studies focused on various types of energy but utilize standard econometric techniques (Charfeddine (2017) along with using carbon emissions as an environmental indicator (Saboori and Sulaiman 2013; Ahmad et al. 2016; Shaari et al. 2020 ). Data This study explores the asymmetric association between energy consumption, economic growth, gross fixed capital formation, and ecological footprint for Pakistan over the period 1971 to 2014. Energy is further divided into oil consumption, gas consumption, coal consumption, and electricity consumption. Ecological footprint is the dependent variable. Energy consumption, economic growth, gross fixed capital formation, oil consumption, gas consumption, coal consumption, and electricity consumption are taken as explanatory variables. The description and data source of all the variables are presented in Table 1 . We have employed the NARDL model proposed by Shin et al. (2014) to explore the nonlinear and asymmetric association among the selected variables. NARDL approach (2020) comprises of a dynamic error correction model that permits to capture the short-run and long-run asymmetries. According to Romilly et al. (2001) and Ahmad et al. (2017) , NARDL methodological framework allows co-integration and asymmetric nonlinearity in one equation jointly and does a better performance in the small samples as compared to other cointegration methods. Moreover, NARDL is more flexible to produce accurate and precise outcomes when variables are I (0), I (1), or a combination of both. Furthermore, when some of the regressors are endogenous, the NARDL will provide unbiased long-run approximations and test statistics. NARDL test takes numbers of lags and archives data generation procedure from general to a specific framework. Due to these features, the NARDL integrates short-run adjustment to long-run equilibrium without the loss of long-run information. Moreover, in NARDL model, we distinguish the shortrun and long-run impacts of independent variables on the dependent variable. According to Shin et al. (2014) problem of multicollinearity should be solved by selecting the appropriate lag order. The simple regression model in linear form can be written as: where EF is ecological footprint, GDP is economic growth, and GFCF is gross fixed capital formation. Considering the nonlinear behavior of variables, Shin et al. (2014) developed the method of NARDL by extending the Pesaran et al. (2001) bound test approach. The ARDL model can be presented as follow: In eq.2, the coefficients attached with difference operators (α 1 , α 2 , α 3 , α 4 ) measure short-run dynamics, whereas the terms with first lagged captures the long-run relationship. After this, to estimate the short-run dynamics, the error correction model (ECM) can be expressed in the following form: Where, ECT t-1 is the error correction term and η indicates the speed of adjustment. We expect a negative relationship between ECM and the dependent variable. By extending the ARDL approach the asymmetric NARDL model is presented as follow: Energy consumption is further disaggregated into oil consumption, gas consumption, coal consumption and electricity consumption. So, now the equations are as follows: where the α i presents the long run and γ i illustrates shortrun coefficient with i=0, 1….7. In Eqs. 4 to 8, the short-run analysis shows the immediate impact of independent variables on the dependent variable, whereas the long-run coefficients show the reaction time of the adjustment in the direction of the equilibrium level. In this study, we applied the time series method to examine how independent variables, namely, energy consumption, economic growth, and capital, describe the correspondence of the environmental quality as extended by EF. EF is presented by EF t , energy consumption by EC t , economic growth by GDP t and capital by GFCF t , oil consumption by OC t , gas consumption by GS t , coal consumption by CC t , and electricity consumption by ECT t . In the structural break date, a dummy variable was recognized through the unit root test (Kim and Perron 2009) and also incorporated in the equation. By applying the Wald test, we investigate the longrun (α = α + = α − ) and short-run (γ = γ + = γ − ) asymmetry for all the variables. p and q present the optimal lag for dependent (EF t ) and independent variables (EC t , GDP t , GFCF t , OC t , GS t , CC t and ECT t ) determined through Akaike information criteria. In the next step, the exogenous variables were decomposed into positive and negative partial sum as follows: The asymmetric co-integration has been done by applying the bound test presented by Shin et al. (2014) for all the regressors with lagged levels. In addition, F-statistic proposed by Pesaran et al. (2001) and t statistic (Banerjee et al. 1998) have been employed for testing the null hypothesis of no cointegration γ = γ + = γ − = 0 and γ =0, respectively. Granger's (1969) work on causality attracted the attention of numerous scholars across the globe. After it scholar applied symmetric approach for determining the direction of casual relationship among the variables. Recently, Hatemi-J (2012) introduced asymmetric causality test among the positive and negative shocks of variables. This test works according to theory in Toda and Yamamoto (1995) by including nonlinear impacts and differentiating between the positive and negative shocks. This test is known as (Hacker and Hatemi-J 2012) causality test and helps in examining the asymmetric relationship between the variables. Firstly, Hatemi-J (2012) supposed that the integrated variables y 1 and y 2 work in a following way: where y 1 and y 2 denote initial values, and t=1, 2…., T and two white noise disturbance term are represented by ε 1i and ε 2i . Onwards these error terms are decomposed into positive and negative shocks: ; 0 denotes the negative changes. In the next step, we decomposed the initial values into positive and negative shocks; the following was observed: The variables in asymmetric form can be presented as follow: After it, every variable is fractionated into positive and negative shocks and presented in cumulative form as Furthermore, for capturing the asymmetric causality between the variables, we applied vector autoregressive introduced by Hatemi-J (2012) within the order of p. For the optimal lag length, we employed the HJC model introduced by Hatemi-J (2003) and Hatemi-J (2008) as follow: The symbol jA J j represents the determinants of the variance-covariance matrix of the white noise disturbance term in the VAR model, where q is lag order, ln denotes the natural logarithm, and T represents the total observations. In the next step, we specified the null hypothesis of the k th elements of y2, which does not cause the w th element of y 1 in the causality test. The hypothesis is tested by applying the Wald test (Hatemi-j 2012) . H o : w th row, k th column in A r is equal to zero for r= 1…p. The null hypothesis of Wald test for no causality is Where β = vec (D) denotes the column-stacking operator, and C represents the p × n (1+ np) indicator matrix.⨂ denotes the Kronecker product, and S u shows the estimated VAR covariance matrix of the unrestricted VAR model as S u ¼ δ 0 This section presents and discusses the empirical findings obtained using EViews 9. Table 2 displays descriptive statistics of the variables used in this study. The statistics show that for Pakistan economy, the maximum and minimum value of ecological footprint is 0.9261 and 0.041, respectively. The country has maximally consumption energy at 6.215, whereas its minimum economy consumption is 5.652. Concerning GDP per capita constant 2010 US, Pakistan's per capita economic growth is 11.390, while the lowest is 10.545. Furthermore, the values of gross fixed capital formation lie between the maximum and minimum values of 10.158 and 8.984. The maximum value for oil, gas, coal, and electricity consumption is 3.032, 3.448, 2.391, and 6.499, respectively, in the country. Likewise, the minimum consumption of oil, gas, coil, and electricity in the country is 0.231, 0.838, −0.668, and 4.995, correspondingly. This study examines the stationarity properties of the variables before estimating asymmetries in the models. Table 3 reports the order of integration of the variables using Dickey and Fuller (1979) and Phillips and Perron (1988) unit root tests. The results of these tests show that all variables are stationary at first difference. Further, all variables at intercept and level have unit root process which later changes to stationary at first different in both ADF and PP unit root tests. Further, the existences of structure breaks in the series are also observed. Perron (1989) pointed out that the presence of structural break when ignored in traditional unit root tests produces biased empirical results. Therefore, structural breaks in unit root tests are important to accommodate for accurate and unbiased estimates. Perron (1989) resolved this problem. Unit root test with structural break dates is presented in Table 3 . The statistics show that all variables have a unit root process (non-stationary) in the presence of a structural break in them. For instance, the ecological footprint is nonstationary in structural break date 2011, energy consumption in 1986, economic growth in 2000, capital formation in 2012, oil consumption in 1988, gas consumption in 1986, coal consumption in 2002, and consumption of electricity in 1983. Therefore, the existence of unit root in these variables motivates us to investigate asymmetric co-integration among environmental quality, energy consumption (in aggregate and disaggregate forms), and economic growth for Pakistan using NARDL modelling. Table 4 presents both short-run and long-run estimates obtained using economic growth and energy consumption (in aggregate form) as focused variables. The long-run results indicate that asymmetric relationship between aggregate energy consumption and ecological footprint holds as coefficients size of positive shock (4.691) and negative shock (19.851) greatly differ. Besides, only the impact of the negative shock is significant, and the coefficient is carrying a positive sign which is inconsistent with the previous studies (Saboori and Sulaiman 2013; Charfeddine and Mrabet 2017; Destek and Sinha 2020) . Therefore, this surprising finding gives us room for moving towards the disaggregate energy consumption analysis for obtaining a clear understanding regarding the relationship between these two variables. The estimates for economic growth indicate that the effect of positive shocks to economic growth improves environmental quality in the long run. This long-run impact of economic growth is supported by the study by Zafar et al. (2019) and Tahir et al. (2020) . The coefficient suggests that 1% increase in economic growth enhances environmental sustainability by 8.886%. The increase in economic growth in the form of a rise in real incomes helps improve environmental quality. This positive impact of economic growth on ecological footprint is supported by studies such as York et al. (2003) , Bagliani et al. (2008) , and Moran et al. (2008) who concluded that after a certain level of economic growth, the rise in it does not damage the environment. As economic growth increases, real incomes of economic agents increase and the affordability and willingness of the masses to pay for mitigating the harmful pollutants increase, thereby enhancing environmental sustainability. Furthermore, high economic growth owing to improvement in technology also produces less pollution. As far as the negative shocks to economic growth are concerned, it has no significant impact on the ecological footprint in the long run. Furthermore, the impact of capital formation on ecological footprint is negative indicating that an increase in capital formation decreases environmental quality by 5.149%. This negative impact suggests that long-run fiscal investments hamper the environment of Pakistan. Our finding implies that as Pakistan's economy is in the developing stage, therefore, most of the investment is made in such technologies that help achieve economic growth targets. These attempts of achieving targeted economic growth at the cost of the environment create pollution and deteriorate environmental quality. Considering short run, the impact of both positive and negative shocks to energy consumption has no significant impact on ecological footprint. The same is the case with its lag term. These findings imply that there is no asymmetric impact of energy consumption on the environmental quality of Pakistan in the short run. Similarly, the energy consumption of the previous year also does not affect the environmental sustainability of the country in the short term. However, concerning the impact of economic growth of the country, it has a significant impact on the ecological footprint. The result indicates that both positive and negative shocks in economic growth affect ecological footprint positively. In other words, all shocks to economic growth lead to environmental sustainability in the short run. The coefficients suggest that a 1% rise in positive and negative shocks individually causes ecological footprint to increase by 3.588 and 9.768%, respectively. These findings are evident that Pakistan's economic growth, despite being small or large, helps sustaining environmental quality and consequently improves the environmental performance index of the country in the short run. Further, capital formation in the short run has the same negative impact on ecological footprint as in the long run. The finding suggests that an increase in capital formation in the country degrades the environment and contributes to global warming. This implies that an increase in capital investment (whether in the long or short run) in the country enhances technological advancement at the cost of environmental degradation and consequent ecological footprint decreases. Concerning the lag term of capital formation, ecological footprint increases by 0.918 percent as a result of 1% rise in the previous year's capital investment in the country. This implies that capital formation has a lag effect on the environment quality of the country. It is because after achieving a certain level of growth through particular capital formation, the country starts focusing on decreasing its environment deteriorating impact which in turn conserves environmental quality. Earlier studies such as Sodersten et al. (2018) also advocated that countries in the early phases of development invest in resource-intensive assets, basic infrastructure, and industrial machinery which lead to environmental deterioration. The error term is negative and Note: Positive and negative variation are denoted by "+" and "-" signs, respectively. *p < 0.1, **p < 0.05, ***p < 0.01 significant indicating that the error correction mechanism exists. Its efficiency suggests that the speed to adjustment is slow as only 40% of disequilibrium recovers to equilibrium in the short run. Table 5 reports the asymmetric effects of energy consumption in its several dis-integrated forms (oil, gas, coal, and electricity) along with the effect of economic growth. Here, findings reveal that different energy sources have different impacts on ecological footprint level of Pakistan which might depend on how much and how frequently each source is exploited. In the case of oil consumption, our study finds that a partial positive sum of oil consumption has a positive and significant impact on ecological footprint while the partial negative sum has a negative but insignificant impact on ecological footprint. The coefficient of positive shock (0.903) indicates that the increased amount of oil consumption is one of the main reasons for large ecological distortions in the country. In Pakistan, a steady increase in population, changes in people lifestyle, and changes in petroleum prices cause an increase in oil consumption, thereby enhancing the use and production of petroleum products. This in turn affects the environment through air and water pollution. The result is in line with the findings of Pao and Tsai (2010, 2011) , and Ozturk and Acaravci (2013) who also reported the environmental deteriorating impact of energy consumption. The positive shock of gas consumption is positively and insignificantly associated with ecological footprint, while the negative shock is negatively and significantly related with the ecological footprint. This means that a 1 % decline in gas consumption can decrease the ecological footprint by 6.976%. This finding is in line with (Munir and Riaz 2020) . In Pakistan, natural gas accounts for the largest share of Pakistan's energy use, representing about 44% of total energy consumption. However, from the past few years, Pakistan's reliance on natural gas consumption has declined due to the introduction and promotion of liquefied natural gas and compressed natural gas use and shifting of cement sector production activities from gas to coal fire system (GoP 2011). This helps the economy to control its ecological footprint resulting from high gas consumption. The asymmetry also holds in the case of coal consumption, electricity consumption, and ecological footprint. The positive shocks to coal consumption and electricity consumption are insignificant and negative while negative shocks are positive and significant. These findings are against our expectations; however, there can be various reasons behind this relationship. As in terms of energy-mix, Pakistan's dependence on thermal, local coal, and imported coal is decreasing from the last few years and also coal's share in total energy consumption is not much significant. Besides, the share of electricity consumption to overall energy is less relative to oil and gas consumption (GoP 2011). Note: Positive and negative variation are denoted by "+" and "-" signs, respectively. *p < 0.1, **p < 0.05, ***p < 0.01 Concerning the effect of shocks (positive) to economic growth on ecological footprint is positive and significant when coal consumption asymmetry is analyzed. Similar is the effect of both positive and negative shocks to economic growth when asymmetry in electricity consumption is estimated. This positive impact of economic growth irrespective of coal and electricity consumption implies that as the economic growth of Pakistan increases, it degraded the environmental quality of the country through increased consumption of nonrenewable resources. Previous studies such as Pao and Tsai (2010, 2011) and Ozturk and Acaravci (2013) also reported the rise in environmental degradation due to an increase in energy consumption. Considering the impacts in the short run, the results indicate that the impact of negative shocks to oil, gas, coal, and electricity consumption increases ecological footprint. Furthermore, the lag of negative shocks to gas consumption also raises environmental degradation. These findings are consistent with the findings of previous studies that advocated that less use of oil, gas, coal, and electricity helps sustain the environment. As we consume less oil, gas, and electricity, the emissions of toxic fumes in the atmosphere lower down. Moreover, less consumption also conserves natural resources and protects the ecosystem from degradation thereby increasing environmental sustainability. These findings of decreasing disaggregated forms of energy consumption and resultant betterment in environmental sustainability are in accordance with the results of Saboori and Sulaiman (2013) . On contrary, the positive shocks to the consumption of coal and electricity insignificantly impact ecological footprint. Furthermore, the lag effect of positive shocks to economic growth is positive and significant in the analysis of oil and gas consumption. Similarly, in the analysis of coal consumption, both positive and negative shocks and only positive shocks to economic growth in case of electricity consumption raise environmental degradation. This impact of economic growth on ecological footprint is in accordance with the results of Uddin et al. (2017) and Danish et al. (2019) Lastly, we derived multiplier dynamics adjustments graphs. Figure 1 shows the cumulative multipliers for energy consumption and environmental quality. The graph indicates that there exists a positive association between energy consumption and ecological footprint, and negative shock in energy consumption is initially dominant than positive shock. The cumulative multipliers for economic growth and ecological footprint are represented in Fig. 2 . The figure shows that a negative association exists between economic growth and ecological footprint (environmental quality). Similarly, the disintegrated forms of energy consumption such as oil, coal, gas, and electricity consumption and their association with ecological footprint are shown in Figs. 3, 4, 5, 6, 7, 8, 9, and 10 , respectively. These figures show that positive shocks in all disintegrated forms of energy consumption are dominant; the negative shocks expect in gas consumption where its negative shock dominates. In econometric analysis pre-estimation (diagnostic) tests are essentially required for accurate empirical estimates and sound policy implications. Table 6 illustrates estimates of several diagnostic tests employed in this study. The findings show that in all the models, there is no issue of serial Similarly, all models are correctly specified and do not suffer from omitted variables misspecification. However, the prevalence of the issue of heteroscedasticity varies across our models. There is the problem of heteroscedasticity when energy consumption in aggregate form is analyzed in model 1. On the other hand, the existence of homoscedasticity is observed when energy consumption in dis-aggregate forms (oil, gas, and electricity consumption) is analyzed except in coal consumption (model 4) where heteroscedasticity exists. The coefficients of R-square are reasonably high (0.9) in all estimated models. Furthermore, there is also no problem of autocorrelation in all the estimated models according to Durbin Watson statistics. Lastly in the diagnostic test, we have employed the nonlinear BDS test given by Broock et al. (1996) . Table 7 reports estimates obtained through the BDS test within VAR consisting of m dimension. The result reveals that nonlinearity exists as the null hypothesis of linearity gets rejected for all the variables. This existence of nonlinearity confirms the presence of disturbance and chaotic behavior in our data. Furthermore, all obtained diagnostic statistics confirm that the specified model of our study is well designed, appropriate for policy suggestion, and the results are consistent and reliable. Over the past few decades, the exploitation of natural resources and degradation of the overall environment have gone up at an alarming rate mainly because of high economic growth concerns. As growth increases production activities which are highly based on energy-intensive inputs, therefore, it has a worsening impact on environmental quality. This study investigated the impact of energy consumption in aggregate and disaggregate forms along with economic growth on ecological footprint over the period 1971-2014 for Pakistan economy. A neoclassical production function has been applied to investigate the asymmetric nexus using the nonlinear ARDL co-integration technique. Our findings confirm the presence of asymmetric associations among the concerned variables. The NARDL results suggest that in long-run equilibrium, negative shocks in energy consumption decrease ecological footprint significantly. This finding validates that as energy consumption decreases, it impedes environmental quality. The same is the case when energy consumption is taken disintegrated in the forms of coal and electricity consumption. However, positive shocks to oil and negative shocks to gas consumption increase ecological footprint. The economic growth of the country is also found to enhance ecological footprint. Furthermore, capital formation in the country decreases the ecological footprint. The empirical research on energy consumption, economic growth, and the environment has been extensively explored in the last few decades. However, these researchers broadly use emissions as a measure of environmental quality which is one component of ecological footprint (Pao and Tsai 2011; Wang et al. 2011; Alvarado et al. 2018; Majeed and Mazhar (2019) . Therefore, these existing studies capture an incomplete picture of the environment. This study contributes to the existing literature in the number of ways. First, it incorporates ecological footprint as a comprehensive indicator of the environment depicting the complete environmental impact of energy consumption and economic growth. Second, it analyses the dynamic effects of both aggregate and disaggregate forms of energy consumption on environmental quality. Third, this study investigated the asymmetric effects of energy consumption and economic growth on ecological footprint. To the best of our knowledge, this is the first study that analyzes the asymmetric effects of energy consumption and economic growth in Pakistan. Fourth, this study performs the analyses using overall and disaggregated forms of energy consumption. The study also covers the long time period of 1971 to 2014. Furthermore, asymmetry in energy consumption, economic growth, and ecological footprint is explored using an efficient technique NARDL for clear findings. On the policy front, this study suggests that policymakers may focus on negative shocks in energy consumption so that environmental quality can sustain. Similarly, disintegrated forms of energy consumption may also be analyzed while forming environmental protection policies. The government may encourage the use of clean energy (biomass) for the fulfillment of the country's energy requirements. Further, economic growth may be achieved through an eco-friendly growth stance. As Pakistan is the agro-based economy, the country may extensively focus on its agricultural sector growth which ultimately helps to achieve economic growth targets. The authorities may invest in pollution-free modern technologies and production equipment, thereby protecting environmental quality. The ongoing COVID-19 crisis has the negative effects on travel, industrialization, and other economic activities, alleviating the burden on the environment. However, the lower demand for energy may negatively affect efforts for transition to clean energy sources. In this scenario, government needs to support decarbonizing projects by strictly focusing on Wald test EC t 6.242*** GDP t 12.782*** GFCF t 5.011** OC t 424.614*** GC t 3.393*** CC t 11.202*** ELC t 3.003* Note: *p < 0.1, **p < 0.05, ***p < 0.01 EF t −0.001 −0.094*** −0.096*** −0.098*** −0.101** EC t 0.202*** 0.344*** 0.440*** 0.512*** 0.559*** GDP t 0.201*** 0.341*** 0.438*** 0.506*** 0.556*** GFCF t 0.116*** 0.169*** 0.267*** 0.351*** 0.415*** OC t 0.154*** 0.240*** 0.338*** 0.410*** 0.464*** GC t 0.177*** 0.294*** 0.386*** 0.448*** 0.491*** CC t 0.161*** 0.261*** 0.356*** 0.427*** 0.479*** ELC t 0.159*** 0.258*** 0.312*** 0.389*** 0.445*** Note: *p < 0.1, **p < 0.05, ***p < 0.01 environment-friendly capital formation projects. This can be done by using updated technology in capital formation to manage environmental sustainability. Availability of data Data sources are clearly mentioned. Interested person can access the data. Computational codes Computational codes are available on demand Author contribution This idea was given by Muhammad Tariq Majeed. Muhammad Tariq Majeed supported in all sections of this work and completed the final write up of the paper. Aisha Tauqir analyzed the data and discussed the results, Maria Mazhar helped in introduction and literature review sections, and Isma Samreen wrote the methodology and concluded the study. All authors have read and approved the manuscript. Ethical approval This article does not contain any studies with human participants or animals performed by any of the authors. Carbon emissions, energy consumption and economic growth: an aggregate and disaggregate analysis of the Indian economy Modelling the CO2 emissions and economic growth in Croatia: is there any environmental Kuznets curve The relationship between income and environment in Turkey: is there an environmental Kuznets curve? Environmental degradation and real per capita output: new evidence at the global level grouping countries by income levels Testing the environmental Kuznets curve hypothesis across the US: evidence from panel mean group estimators A consumption-based approach to environmental Kuznets curves using the ecological footprint indicator The effect of financial development on ecological footprint in BRI countries: evidence from panel data estimation Error-correction mechanism tests for cointegration in a single-equation framework Asymmetric impact of energy consumption and economic growth on ecological footprint: using asymmetric and nonlinear approach CO2 emissions, energy consumption, economic growth, and financial development in GCC countries: dynamic simultaneous equation models A test for independence based on the correlation dimension The impact of energy consumption and economic development on Ecological Footprint and CO2 emissions: evidence from a Markov Switching Equilibrium Correction Model The impact of economic development and social-political factors on ecological footprint: a panel data analysis for 15 MENA countries Dynamic linkages between road transport energy consumption, economic growth, and environmental quality: evidence from Pakistan Role of renewable energy and non-renewable energy consumption on EKC: evidence from Pakistan Energy production, economic growth and CO2 emission: Evidence from Pakistan Modeling the impact of transport energy consumption on CO2 emission in Pakistan: evidence from ARDL approach Linking economic growth and ecological footprint through human capital and biocapacity Renewable, non-renewable energy consumption, economic growth, trade openness and ecological footprint: evidence from organisation for economic Co-operation and development countries Distribution of the estimators for autoregressive time series with a unit root Global climate risk index 2020. Germanwatch Energy consumption, economic growth and carbon emissions: cointegration and causality evidence from selected African countries Economic growth and pollutant emissions in Tunisia: an empirical analysis of the environmental Kuznets curve Economic Advisor's Wing, Ministry of Finance Economic Advisor's Wing, Ministry of Finance Economic Advisor's Wing, Ministry of Finance Investigating causal relations by econometric models and cross-spectral methods Economic growth and the environment A bootstrap test for causality with endogenous lag length choice: theory and application in finance Role of institutions in correcting environmental pollution: an empirical investigation A new method to choose optimal lag order in stable and unstable VAR models Forecasting properties of a new method to determine optimal lag order in stable and unstable VAR models Asymmetric causality tests with an application Ecological vulnerability assessment for ecological conservation and environmental management A three-level framework for balancing the tradeoffs among the energy, water, and air-emission implications within the life-cycle shale gas supply chains Global energy demand rose by 2.3% in 2018, its fastest pace in the last decade The CO 2 emissions-income nexus: evidence from rich countries Unit root tests allowing for a break in the trend function at an unknown time under both the null and alternative hypotheses Effects of urbanization on freight transport carbon emissions in China: common characteristics and regional disparity Renewable energy, water, and environmental degradation: a global panel data approach Financial development and ecological footprint: a global panel data analysis Reexamination of environmental Kuznets curve for ecological footprint: the role of biocapacity, human capital, and trade Effects of urbanization, industrialization, economic growth, energy consumption, financial development on carbon emissions: an extended STIRPAT model for heterogeneous income groups The asymmetric relationship between financial development and CO 2 emissions: the case of Pakistan Measuring sustainable development-nation by nation Testing the Kuznets Curve hypothesis for Qatar: a comparison between carbon dioxide and ecological footprint Asymmetric impact of energy consumption on environmental degradation: evidence from Australia, China, and USA CO2 emissions, energy consumption and economic growth in the ASEAN-5 countries: a crosssectional dependence approach Carbon dioxide emissions and economic growth: panel data evidence from developing countries Ecological footprint, urbanization, and energy consumption in South Africa: including the excluded Energy consumption, economic growth and environmental degradation in OECD countries Measuring the impact of alternative and nuclear energy consumption, carbon dioxide emissions and oil rents on specific growth factors in the panel of Latin American countries The long-run and causal analysis of energy, growth, openness and financial development on carbon emissions in Turkey Investigating the environmental Kuznets curve hypothesis: the role of tourism and ecological footprint Empirical tests and policy analysis of environmental degradation at different stages of economic development (Working Paper No. 992927783402676). International Labour Organization, Geneva Pao HT, Tsai CM (2010) CO2 emissions, energy consumption and economic growth in BRIC countries Multivariate Granger causality between CO2 emissions, energy consumption, FDI (foreign direct investment) and GDP (gross domestic product): evidence from a panel of BRIC (Brazil The great crash, the oil price shock, and the unit root hypothesis Bounds testing approaches to the analysis of level relationships Testing for a unit root in time series regression The long-term relationship between CO2 emissions and economic activity in a small open economy: Uruguay 1882-2010 Car ownership and use in Britain: a comparison of the empirical results of alternative cointegration estimation methods and forecasts Environmental degradation, economic growth and energy consumption: evidence of the environmental Kuznets curve in Malaysia Mitigating degradation and emissions in China: the role of environmental sustainability, human capital and renewable energy Environmental quality and development: is there a Kuznets curve for air pollution emissions? The effects of energy consumption and national output on CO 2 emissions: new evidence from OIC countries using a panel ARDL analysis Economic growth and environmental quality: Time-series and cross-country evidence (Working Paper No. 904) Environmental Kuznets curve hypothesis in Pakistan: cointegration and Granger causality Modelling asymmetric cointegration and dynamic multipliers in a nonlinear ARDL framework Financial development and its moderating role in environmental Kuznets curve: evidence from Pakistan Environmental impacts of capital formation Pakistan management of green transportation and environmental pollution: a nonlinear ARDL analysis Investigating the pollution haven hypothesis in Ghana: an empirical investigation The impact of financial development and globalization on environmental quality: evidence from South Asian economies Statistical inference in vector autoregressions with possibly integrated processes Ecological footprint and real income: panel data evidence from the 27 highest emitting countries A reinvestigation of EKC model by ecological footprint measurement for high, middle and low-income countries CO 2 emissions, energy consumption and economic growth in China: a panel data analysis Renewable energy consumption, economic growth and human development index in Pakistan: evidence form simultaneous equation model The relationship between biomass energy consumption and human development: empirical evidence from BRICS countries WDA (2020) World Data Atlas WEF (2012) Industry agenda energy for economic growth Footprints on the earth: The environmental consequences of modernity The impact of natural resources, human capital, and foreign direct investment on the ecological footprint: the case of the United States Motivation and challenges for e-commerce in e-waste recycling under "Big data" context: a perspective from household willingness in China Regional differences and convergence analysis of energy efficiency in China: on stochastic frontier analysis and panel unit root Construction of electric vehicle driving cycle for studying electric vehicle energy consumption and equivalent emissions Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations Consent to participate I am free to contact any of the people involved in the research to seek further clarification and information. Conflict of interest The authors declare no competing interests.