Modeling the Effects of Agricultural Innovation and Biocapacity on Carbon Dioxide Emissions in an Agrarian-Based Economy: Evidence From the Dynamic ARDL Simulations

In this paper, we modeled the effects of income, agricultural innovation, energy utilization, and biocapacity on Carbon dioxide (CO2) emissions. We tested the validity of the environmental Kuznets curve (EKC) hypothesis for Nigeria from 1981 to 2014. We applied the novel dynamic autoregressive distributed lag (ARDL) simulations to develop conceptual tools for policy formulation. The empirical results confirmed the EKC hypothesis and found that agricultural innovation and energy utilization have an escalation effect on CO2 emissions whereas income and biocapacity have long-run emission-reduction effects. The causality results found agricultural innovation attributable to CO2 emissions and observed that income drives energy demand. Income, biocapacity, and energy utilization are found to predict changes in CO2 emissions. These results are validated by the innovation accounting techniques—wherein 22.79% of agricultural innovation corresponds to 49.43% CO2 emissions—5.95% of biocapacity has 35.78% attributable CO2 emissions—and 1.61% of energy spurs CO2 emissions by 16.27%. The policy implication for this study is that energy efficiency, clean energy utilization and sustainable ecosystem recovery and management are the surest ways to combat climate change and its impacts.


INTRODUCTION
Mitigation of climate change and its impacts on the environment and wellbeing are important global issues in recent times. Climate change has a traceable course to excessive use of "unclean" combustible energy, which disrupts the levels of carbon in the atmosphere, resulting to the preservation of heat in the atmosphere (See Kasman and Duman, 2015;Usman et al., 2019;Rafindadi and Usman, 2019;Agboola and Bekun, 2019;Usman et al., 2020a;Usman et al., 2020b). Research on energy utilization and economic outgrowth effects of CO 2 emissions has received significant attention in the literature of environmental management. Essentially, within the theoretical account of Environmental Kuznets Curve (EKC) hypothesis, it is reported that economic development initially triggers environmental pollution with increasing levels of income but declines afterward at specified threshold of income level where environmental awareness remains a priority (Grossman and Krueger, 1991). A significant number of the extant literature have tested the validity of the EKC hypothesis over the years lacking consensus. The empirical results from most studies are that economic growth trajectory heightens environmental pollution but declines thereafter following improvements in livelihood and environmental awareness, thereby validating the EKC hypothesis Rafindadi, 2016;Shahbaz et al., 2017;Mesagan et al., 2018;Rafindadi and Usman, 2019). On the contrary, some studies aptly posit that energyintensive economic outgrowth and environmental quality is not in line with the EKC hypothesis (Inglesi-Lots and Bohlmaann 2014; Nasr et al., 2015). Therefore, the EKC-based empirical findings are mixed and conflicting, hence, require further empirical validation. Despite mitigating efforts by world leaders geared toward CO 2 sequestration, a substantial rising of the contribution of CO 2 to greenhouse gas (GHG) emissions are reported over the years (IPCC, 2017). It is reported that CO 2 contributes 76.6% of GHG emissions generated mostly by developing economies in the quest to sustain economic productivity. Between 1961 and 2011, CO 2 emissions rose from ∼9.4 billion metric tons to ∼34.6 billion metric tons (IPCC, 2013). Equally, CO 2 emissions increased from ∼29.7 billion tons to ∼33.4 billion tons between 1999 and 2017 (BP, 2018). In Nigeria, CO 2 emissions remain a major threat to both human and ecosystem development. As reported by the World Bank (2015), as of 2014, Nigeria emitted 96,280.75 kilotons of CO 2 , which was lower than 106,067.98 kilotons in 2005.
A large body of literature has linked climate change to agricultural practices. As recently emphasized by Owusu and Asumadu (2016), Shabbir et al. (2020), Agboola and Bekun (2019), in addition to excessive consumption of energy from the fossil fuel sources, agricultural practices have a substantial effect on GHG emissions. Agriculture ranked is as the secondhighest contributor of GHG emissions and global warming, contributing roughly 21% of the global anthropogenic GHG emissions in the world (Blanco et al., 2014). This is because most agricultural practices require greater energy consumption, mostly sourced from fossil fuels (Blanco et al., 2014). Agriculture may affect the ability of land to absorb heat and light, which can lead to radioactive forcing. More so, deforestation and desertification resulting from land use and fossil fuels can exert upward pressure on anthropogenic carbon dioxide. Besides, raising livestock such as cattle, pigs and poultry may contribute to methane and nitrous oxide concentrations and emissions. On the other hand, agriculture can substantially reduce the level of carbon emissions as opined by the United Nations Food and Agricultural Organization (FAO), (2016). This is supported by Reynolds and Wenzlau (2012) who posit that agriculture innovation is reported to have a mitigation effect on CO 2 emissions. For example, some modern agricultural practices can be powered by clean energy to reduce the effects of the use of pesticides, irrigation, soil tillage, deforestation, and waste from the plastic mulch, stubble burning, and other channels of GHG emissions.
Our study, therefore, hypothesizes that the effects of agricultural innovation and biocapacity on CO 2 emissions have long-and short-term environmental consequences in Nigeria. Given that Nigeria is an agrarian nation blessed with natural resources, there are reports of its citizens engaging in crude methods of agricultural practices that hamper environmental sustainability. However, scientific literature on the scope is limited for policy formulation. More so, Nigeria is ranked among the top 10 countries with a dangerous precedent of ambient air pollution (HEI, 2018). Besides, a recent study ranked Nigeria as the sixth among 195 nations with the most approximate cases of disability-adjusted life years from exposure to air pollution (Owusu and Sarkodie, 2020a). Thus, justifies the need to investigate the effects of agricultural innovation and ecosystem dynamics on CO 2 emissions. This will have policy implication not only on carbon sequestration but mitigating mortality and morbidity rates. Therefore, insights from our study will provide supporting evidence for policymakers in designing appropriate energy and environmental policies for CO 2 sequestration that underpins the Sustainable Development Goals (SDGs). In terms of methodology, we use Lee-Strazicich (L-S) structural break, causality test and novel dynamic autoregressive distributed lag (ARDL) simulations approach-to estimate the out-sample parameters of counterfactual shocks in specific time periods and specified exogenous regressor useful for policy formulation. This is the first time such a novel out-sample, stochastic and simulation technique has been utilized in extant literature for the proposed theme.

LITERATURE REVIEW
The EKC hypothesis from the pioneering work of Kuznets (1955) underpins the framework for this study. In its generic form, Kuznets observed a nexus between income per capita and inequality in such that income inequality would first rise and decline as income increases. This hypothesis led to what is known as EKC by Grossmann and Krueger (1991). The EKC hypothesis postulates a parallel increase of both income level and emissions until a threshold of income is achieved before a reduction in emissions can be noticed thereafter. This hypothesis explains the trade-off between sustained economic productivity and environmental sustainability.
The nexus between economic productivity and ecological degradation has gained prominence in extant literature since the mid-90s. For example, a study found an "inversed U-shaped" relationship where ecological pollution would increase at the early stages of economic development but after a specified threshold, economic outgrowth tends to mitigate ecological pollution (Selden and Song, 1994). Similarly, several studies have all reported an inverted U-shaped nexus between economic growth and CO 2 emissions (Galeotti et al., 2006;Shahbaz et al., 2013;Rafindadi and Usman, 2019;Ike et al., 2020a;Usman et al., 2020b). For example, Shahbaz et al. (2013) applied the ARDL cointegration approach to investigate the effect of energy intensity, economic growth, and globalization on CO 2 emissions in Turkey. The findings documented the presence of EKC and further revealed economic growth and energy intensity exert positive pressure on CO 2 emissions while globalization reduces CO 2 emissions. Similarly, a study by Rafindadi and Usman (2019) using ARDL modeling approach with controlled structural breaks validated the EKC hypothesis for South Africa. A recent paper by Ike et al. (2020a) using a novel quantile regression via quantile moments confirmed the EKC hypothesis by controlling for oil production in oil producing nations.
On the contrary, some studies reported that the EKC hypothesis might not hold always. For example, "N-shaped" relationship between productivity and emissions following a hike in CO 2 emission was observed for a small open economy and industrialized country (Fried and Getzner, 2003). Similarly, it is reported that the validity of the EKC is not certain in all circumstances, hence, there is no certainty that an inversed-U shaped link exists between economic productivity and pollution (See Spangenberg, 2001). In a study by Nasr et al. (2015) found no evidence to support the EKC in South Africa using a cosummability technique with a century of data.
In recent times, many studies have incorporated the role of energy utilization in testing the validity of the conventional EKC hypothesis. The EKC hypothesis was tested in Romania by incorporating energy utilization . The findings confirmed the EKC hypothesis and further revealed energy utilization attributable CO 2 emissions. Tiwari et al. (2013) found EKC and bi-directional causality between growth and CO 2 emissions from accounting for coal, growth, and trade in India. This means that economic growth first increases with environmental pollution but after reaching a turning point, increasing productivity improves environmental quality. Using the ARDL approach for Portuguese economy over the period 1971 to 2008, the EKC was validated in both short-and long-run in the presence of international trade, urbanization, and energy consumption. The effects of coal energy, industrial production and emissions were investigated in China and India (Shahbaz et al., 2014). The results identified an inversed U-shaped for India and U-shaped for China. It further showed that coal consumption causes CO 2 emissions in India while the feedback effect is observed in China. The impact of energy and democracy on CO 2 emissions was investigated in India using the ARDL methodology and found that, while energy increases CO 2 emissions, democracy perhaps mitigates emissions . Also, Usman et al. (2020b) incorporated globalization, democracy, and energy consumption in a standard EKC model for South Africa and confirmed an inversed U-shaped link between growth and emissions of CO 2 . Similarly, the EKC hypothesis was confirmed in Thailand, using heterogeneous fossil fuel sources (Ike et al., 2020a).
Based on panel data settings, the interaction of income and CO 2 emissions was assessed in 43 developing countries (Narayan and Narayan, 2010). The results revealed that CO 2 emissions significantly dropped with a rise in income, suggesting that the hypothesis of EKC fails to hold. Conversely, Apergis (2016) investigated the real GDP-CO 2 emissions nexus in 15 countries and showed evidence of the EKC in most of the countries. More recently, Ike et al. (2020b) reported EKC for 15 oil-producing countries while exogenizing crude oil, electricity, trade, and democracy. This understanding is supported by Ike et al. (2020c) who found EKC for a panel of G-7 both in country-specific and panel settings.
Unlike most studies, very few pieces of extant literature tested for the EKC by exogenizing agricultural production. For example, evidence of EKC with agriculture reducing the level of CO 2 emissions in Turkey was reported (Dogan, 2016). Gagnon et al. (2016) divulged that agriculture has no impact significantly on emissions of carbon dioxide in Canada. Gokmenoglu and Taspinar (2018) investigated the role of agriculture in inducing CO 2 emissions in Pakistan. The empirical results observed the existence of EKC and further discovered that agriculture increases CO 2 emissions. Furthermore, feedback causal relationships are noticed among GDP, energy, agriculture, and CO 2 emissions. The EKC position in Nigeria examined by controlling for agriculture and foreign direct investment (Agboola and Bekun, 2019). The results obtained echoed the EKC hypothesis and thus documented that agriculture deteriorates the environment in Nigeria.
A panel data methodology was used to analyze the effect of agriculture on CO 2 emissions for Southeast Asian countries (Liu et al., 2017). The finding failed to lend support for the EKC. The study revealed that agriculture reduces CO 2 emissions with causality from renewables to CO 2 emissions and from growth to agriculture. On the contrary, an increase in agriculture was found to reduce CO 2 emission in five MENA countries (Ben Jebli and Ben Youssef, 2017). Based on the causality, it was discovered that agriculture Granger-cause economic growth while energy causes agriculture. However, the bi-directional linkage was found for agriculture and CO 2 emissions. These findings of course are similar to Olanipekun et al. (2019) who found a positive effect of agricultural production on pollution in Africa.
The existing literature on agriculture-induced CO 2 emissions is very few and scanty, particularly for Nigeria. The only existing country-specific study on Agriculture-CO 2 emission linkage in Nigeria is a recent study by Agboola and Bekun (2019), which suffers from misspecification problems. For example, the authors used the log forms of agricultural value-added and trade which are in percentages and hence growth rates. Taking a log of growth rate is technically wrong and could lead to spurious regression. Another methodology problem suffered by the study is the application of a standard Granger causality test withoutf meeting its fundamental assumption. As noted in the literature, a traditional Granger causality is used only when the series are all in levels. The work by Olanipekun et al. (2019) is based on the panel of African countries, which have country-specific problems. Therefore, the findings may have limited policy implications for Nigeria. Also, the existing studies failed to capture structural breaks in the variables which could alter CO 2 emissions in the long run. Therefore, to properly model agriculture-induced CO 2 emissions and EKC in Nigeria, we incorporated the structural breaks into our model to examine their effects on the endogenous variable in the long run. Finally, since Nigeria is blessed with diverse natural resources, we control for biocapacity to capture the ability of the ecosystem to produce biological materials demand of the people.

Data Collection
We employed time-series data spanning 1981-2014, selected due to data availability. 1 The variables in the models include CO 2 emissions per capita as an endogenous variable while real GDP per capita, which represents "second order polynomial of real GDP per capita (GDP 2 ), agricultural value-added; biocapacity and energy per capita (EU) are exogenous variables. Generally, CO 2 emission per capita measures environmental quality. Real GDP per capita is used as a proxy for income or wealth, agricultural value-added per capita is used as a proxy for agricultural innovation since value is added to the raw materials of agriculture while Biocapacity per capital measures the ecosystem recovery. CO 2 emissions, real GDP, agricultural innovation measured by agricultural value-added, and Energy Use are obtained from the World Development Indicators (WDI) database, 2, while Biocapacity is retrieved from the Global Footprint Network (GFN) database. 3 The selections of these variables are guided by the United Nations' long-term plan for Sustainable Development Goals (SDGs) which emphasizes clean energy, growth, and environment. Particularly, we included energy use to tackle goal 6, which targets clean energy and water, energy use. Goal 7, which is centered on the affordability of clean energy, is facilitated by improvement in agriculture and biocapacity. We believe that once agriculture is stimulated coupled with biocapacity, people would be able to afford clean energy. We included GDP to capture goal 8, which is concerned with achieving decent work and growth without causing damage to the environment. Finally, goals 9 and 13, which are concerned with climate change and carbon sequestration, are represented by CO 2 emissions. The variables, measurements and source are described in Table 1.  (2019), Usman et al. (2020b), the standard EKC framework is expressed as:

Model Specification
Where Φ 0 is the constant, CO 2 is the carbon emissions, which measure environmental quality. Y t is Real GDP, which measures income while the squared term of real GDP (Y 2 t ) is added to determine whether the validity of the EKC hypothesis. μ t embodies the error term that is invariably presumed to be normally distributed. In this study, we incorporated agricultural and biocapacity variables into the standard EKC framework. This is because, agricultural activities and biocapacity of a country could contribute or mitigate the rate of carbon emissions as documented in the earlier studies by Dogan (2016), . Therefore, our model will be expressed as follows: Where CO 2 , Y, and Y 2 t remain as defined in Eq. 1. AGR represents agricultural value-added per capita, a measure of agricultural innovation; BCP t is the biocapacity per capita; EU t represents per capita energy consumption, t stands for time period while ε t denotes that the residual term is a white noise process with variance σ 2 , ε t ∼ iid(0, σ 2 ). The natural logarithmic regression of Eq. 2 is given as follows: Equation 3 is a log-log regression of Eq. 2 to explain the impacts growth in the long-run. To this extent, all the variables remain as defined in Eqs. 1 and 2. ln denotes the natural logarithm of the series. If the variables have a long-run relationship between them, it therefore, means that they will have a level relationship specified with long-run parameters so that they can follow the pattern of error correction model (ECM). The long-run and short-run parameters are obtained through a dynamic restricted ECM, resulting from the ARDL approach proposed by Pesaran et al. (2001) as given below: Where the variables remain as defined previously. Δ is a difference operator generically defined as Δy t y t − y t−1 . The long-run coefficients are obtained from the first part of Eq. 4. The errorcorrection term (ECT) can also be obtained as; The parameters θ 1 , θ 2 , θ 3 , θ 4 , and θ 5 are the long-run effects of all the explanatory variables on CO 2 emissions. Therefore, to capture the adjustment speed from short-run disequilibrium to long-run equilibrium, we estimate the conditional error correction model given as: where the speed of adjustment speed is captured by ect t−1 which is defined as the first lag of the residual of the short-run parameters are given by β i 's. To test for possible variable cointegration, we applied a level equation based on Eq. 4. As recommended by Pesaran et al. (2001), an F-test is used for testing the null hypothesis, which states that α 1 α 2 α 3 α 4 α 5 0 and the alternative hypothesis, which states that α 1 ≠ α 2 ≠ α 3 ≠ α 4 ≠ α 5 0. This methodology has some enviable advantages. First, it estimates both short-and long-run parameters of the model used in this study. Second, our model is suitable for mixed order of integration.
In other words, this model can be applied regardless of variables integrated of order zero, or order one, or mutually cointegrated. Third, the estimation approach yields robust and unbiased estimates irrespective of the sample size. This means that the model is more appropriate in our case-where the number of observations is thirty-four.

Lee-Strazicich Unit Root Test
The existing traditional unit root tests are found to be inadequate and as such provide false outcomes when structural breaks are present in the series. To avoid this, in addition to the Augmented Dickey-Fuller (ADF) and Phillips-Perron (PP) tests, we applied a minimum LM unit root test with one break (Lee and Strazicich, 2003). This test accommodates information concerning a single unknown break and tackles the inaccuracy problem of identified breakpoint under the null and alternative hypotheses. To this end, Lee-Strazicich unit root test is more superior to all other structural break unit root tests in the literature. In testing the unit root via this test, we applied a "crash" model which permits for a one-time change in intercept, under the alternative hypothesis with the optimal number of lag k determined by beginning the test from the general-to-specific method (Perron, 1989). To perform this test, we began with the maximum number of lagged first-differenced terms, k 8 and continue to reduce the lagged term if the model is insignificant. The null hypothesis H 0 : α 0 is checked against the alternative hypothesis H 1 : α < 0. These hypotheses also hold for typical unit root tests applied. We ruled that the series has a unit root in the presence of a break if the test value is less than the critical value at 1, 5, and 10% significant levels.

Causality Test
fWe ascertained the direction of causality by applying a Granger causality test within the Toda-Yamamoto framework (Toda and Yamamoto, 1995) which applies a modified Wald statistic. The method involves estimating a vector autoregressive VAR (p) with extra lag d. this generally denotes (p + d max ), where p denotes the VAR order and d is the extra lag (d max ) which is the maximum order of integration in the VAR system. To apply this method, we augmented the correct VAR order p with d extra lag and used the asymptotic χ 2 distribution of the Wald statistic to assess the existence of a causal relationship. This method is widely accepted in the literature to be superior and richer than the standard Granger causality test or VECM causality test. Particularly, the test is suitable and provides robust results regardless of the integration order of the series and their co-integration. Therefore, the VAR (p + d max ) is expressed as follows: From Eq. 6, the Granger causality running from ln EU t to lnCO 2t implies that ξ 16i ≠ 0∀ i similarly Granger causality running from lnCO 2t to ln EU t implies that ξ 16j ≠ 0∀ j .
The framework for our model is shown in Figure 1, which begins with ARDL specification and estimation as well as residual and stability diagnostic tests. The second stage is the estimation of the structural model based on impulse-response and variance decomposition analyses.

Statistical Analysis
The mean of the variables showed that incomes have the highest meanwhile CO 2 emissions and biocapacity have low and negative mean scores. The standard deviations are also low with energy use having the lowest. This suggests that all the variables are less volatile over the study period. The skewness of the variables indicates that CO 2 emissions and biocapacity are negatively skewed while income, the squared term of income, agriculture, and energy use are positively skewed with the values tending toward zero. More so, the kurtosis of the variables indicated that all the series have a positive kurtosis with Jarque-Bera values exceeding the region of normal distribution as can be seen by the probability values.
The graphical plots of the variables in described in Figure 2. This is necessitated by the presence of drift, trend, and seasonality as well as structural breaks. As shown by the Figure, all the variables seem to have structural breaks. These breaks are more evident in CO 2 emissions, biocapacity, and energy use with no precise evidence of a trend. For income, squared income, and agriculture, it is observed that the variables begin to trend upward.

Stationary Test Results
Before estimating the model for this study, we first, applied the usual unit root tests via ADF and PP as earlier stated. The results given in Table 2, Panel A show that all the series (CO 2 emissions, Income, the square of income, agricultural innovation, biocapacity, and energy use) are not stationary in their levels. However, after we took their first differences, they all turn out to be stationary. This means that the variables are classified as I (1) process. To circumvent the inadequacy of conventional unit root tests, we applied the minimum LM unit root test with one break. The results as displayed in Table 2, Panel B validated the earlier results that all the series are integrated of order one, i.e., I (1) process. Also, the identified breakpoint for CO 2 emissions is 1999, income and its squared term is 2006; agricultural value added is 2001, biocapacity is 2010, and energy use is 2002. The break in 1999 could be attributed to the effect of general elections which lowers the pressure on stimulating growth and hence CO 2 emissions. The break in 2002 may be caused by the effect of pre-2003 general elections. The 2006 break in income and its squared term can be attributed to exchange rate volatility, which significantly affected income levels. Finally, the 2010 break was caused by 2008 worldwide financial disaster which affected the agricultural sector significantly.

Co-Integration Tests
Having established the integrating properties of the variables in our model, the next is to check whether co-integration exists among the variables. To do this, we applied the ARDL bounds testing approach. The robustness of this test is carried out based on the combined cointegration test (Bayer and Hanck, 2013). The lag length selection of three based on the Akaike information criterion (AIC) is shown in Table 3 while Table 4 provided the reports of the bounds-testing co-integration. According to the reports, we found that when each of the variables is treated as endogenous, we confirmed five cointegrating vectors, which by implication means that a longrun relationship exists between the sampled series. These findings are validated by the combined co-integration test of Bayer-Hanck Table 5, which found a co-integration in all the six equations, implying that there is a long-run nexus between the investigated series.  Table 6 reports the long-run and short-run parameters of the ARDL model estimator. Based on the parameters of the model, we find evidence that real income and its squared term have a positive and negative relationship with CO 2 in the long run and short run, respectively. The negative effect of squared term of income indicates a breakaway of CO 2 emissions and real income at higher income level. This result, therefore, suggests the validation of the EKC hypothesis in Nigeria both in the long run and short run. The plausible reason for the findings is that Nigeria being an oilexporting country mostly engages in excessive use of fossil fuels and cement manufacturing. Furthermore, a larger carbon is emitted during the utilization of liquid and gas fuels as well as gas flaring. Therefore, the validity of the EKC hypothesis in this study is consistent with previous studies such as Galeotti et al., (2006 Ike (2020aIke ( , 2020bIke ( , 2020c, Iorember et al. (2020). The effect of agricultural innovation on CO 2 emissions is positive, inelastic, and statistically significant both in the long run and short run. This implies that a 1% increase in agricultural innovation would cause CO 2 emissions to rise by 0.5145% in the long run and 0.5329% in the short run. The economic reason supporting this result is that agricultural practices such as bush burning, tillage, fertilization, deforestation, and desertification as well as raising livestock like cattle, pigs, fish and poultry could accelerate the level of anthropogenic carbon emissions. This finding agrees with a study that found a positive relationship between agriculture and CO 2 emissions in Nigeria (Agboola and Bekun, 2019). Our result also corroborates with Olanipekun et al. (2019) who found a similar result for African countries and for Tunisia (Ben Jebli and Ben Youssef, 2017). Moreover, we found that after taking the first lag of agricultural value-added, its effect on CO 2 emissions was negative, indicating that the historical effects of agricultural value-added underpin CO 2 emissions mitigation. The negative relationship between agriculture and CO 2 emissions is supported by a finding documented for 53 countries in the world (Rafiq et al., 2016); five MENA countries (Ben Jebli and Ben Youssef, 2017), and Turkey (Dogan, 2016). Furthermore, the influence of biocapacity on CO 2 emissions is negative, inelastic and substantial in the long run while in the short run, it is negative, elastic and significant. Particularly, a 1% increase in biocapacity would reduce CO 2 emissions by 0.1853% in the long run, while in the short run, it reduces CO 2 emissions by 1.1179%. This is because biocapacity is a non-carbon measurement of the ability of the ecosystem to renew the biological materials demand by the people from the earth's surfaces. Therefore, it is consistent with the Sarkodie and Strezov (2018) who found a negative relationship between biocapacity and CO 2 emissions for the US, Australia, China and Ghana. Finally, the influence of energy use is positive, elastic and statistically significant with CO 2 emissions. This means that a 1% increase in energy use would increase CO 2 emissions by 4.2620% in the long run and 1.1211% in the short run. The results further showed that from one lag period afterward, the effect of energy use on CO 2 emissions turns negative. The implication for this result is that most of the Nigerian energy sources are stemming traditional biomass and waste, which could explain about 83% of the total primary production, while 16% is accounted for by the fossil fuels and 1% by hydropower. These energy sources are renewables, which emit low carbon and GHGs. This reason is also attributed to the negative effects of agriculture from the first lag afterward.
The speed of adjustment (ECT t-1 ) is negative and significant with a value −0.9938. This implies that the speed of convergence from short-run variation toward equilibrium long run is about 99% yearly. We also tested the diagnostics of the model estimated. The results showed that there is no case of serial correlation and conditional heteroscedasticity problems. Similarly, the functional form of the model is correctly constructed with evidence that the error term is normally distributed. Furthermore, apart from the RAMSEY RESET test, we applied the cumulative sum (CUSUM) and CUSUM squares (CUSUM Sq.) to test the stability of the model. As shown in Figure 3, both tests revealed that the model is stable and adequate both in the long and short run.

Causality and Innovation Accounting
Theoretically, if a co-integration is found, there must be at least causality between the variable. As displayed in Table 7, we found evidence that a uni-directional causality runs from agriculture to CO 2 emissions, which contradicts the earlier finding by Agboola and Bekun (2019). The plausible reason could be attributed to the fact that the study applied a standard Granger causality which tends to produce a spurious result if the variables are not all integrated at levels. However, our finding agrees with Ben Jebli and Ben Youssef (2017) who found  Notations: ***, **, and * denote statistical significance level at p-value < 0.01, <0.05, and <0.10. ADF, Augmented Dickey-Fuller Test; P-P, Phillips-Perron Test; L-S, Lee-Strazicich test.  Notations: ** refers to the rejection of no level relationship at 5% significance level. The critical value is determined with unrestricted intercept and no trend. The maximum lag order is three and the optimal lag order is selected by the Akaike Information Criterion (AIC).
Frontiers in Energy Research | www.frontiersin.org February 2021 | Volume 8 | Article 592061 8 agriculture and CO 2 emissions to have a causal link in the long run for five MENA countries. We also found that CO 2 emission could predict income, biocapacity, and energy use. Furthermore, our results provide evidence that a bidirectional causal relationship exists between agriculture and biocapacity as well as agriculture and energy use. These results imply that agriculture causes biocapacity and energy use and vice versa. The results that agriculture has predictability for energy use are consistent with Agboola and Bekun (2019). This is also consistent with Ben Jebli and Ben Youssef (2017) who found a long-run causality running from renewable energy to agriculture. There is also evidence that income level and its squared term have predictability for energy use. This result also corroborates a similar reported case in Ben Jebli and Ben Youssef (2017).
We step forward to validate our findings via the innovation accounting test of variance decomposition and impulse-response function analyses based on 10-year forecast horizons. From Table 8, we found that except for energy use, CO 2 emissions have the highest contribution to the variance decomposition of all the variables in the model. Similarly, energy use has the lowest contribution to the variance decomposition of all the variables. Starting from the variance decomposition of CO 2 emissions, we observed that own shock contributed about 65.4%, followed by the contribution from agriculture which accounted for about 11.65%. Energy use has the lowest contribution of 1.46%, which confirms the earlier results that about 83% of total energy consumption in Nigeria stems from the renewables which emit low carbon dioxide. More so, from the variance decomposition of income and its squared term, we found that CO 2 emissions contributed about 56.01 and 56.4%. This is followed by the contribution of agriculture, which accounted for about 22.8 and 22.7%, respectively. The contribution of energy use is about 1.60%. We further found that while agriculture contributed about 32.01% due to own shock, the contribution of CO 2 emissions is about 49.42% while energy use is about 2.41%. The results further suggested that for variance decomposition of bio-capacity, own shock contributed just 13.23% while CO 2 emissions contributed about 35.78% with 1.71% contribution from energy use. Additionally, the highest contributor to the variance decomposition of energy use is squared term of income with about 31.24%, apparently followed by agriculture with about 23.09%. The contribution from its own shock is about 3.83%. Therefore, from the results of the forecast error variance decomposition, we observed that 22.79% of agriculture corresponded to 49.43% CO 2 emissions. We also found that 5.95%   Notations: ***, ** and * denote rejection of the null hypothesis at 1, 5, and 10% significant levels. p-values are presented in parenthesis (. biocapacity caused 35.78% CO 2 emissions, while 1.61% of energy use led to just 16.27% CO 2 emissions. Figure 4 presents the impulse responses of all the variables to an innovation shock. As shown, CO 2 emissions responded positively to the innovation shocks up to the eighth horizons and consequently turned negative. This implies that about 8th horizons, CO 2 emissions responds negatively to innovation shocks. For income, we found that the response of income to innovation shocks is positive until sixth horizon. The response became negative between sixth and eighth horizons, after which it became positive. The same is not observed in the case of income squared. The response of the square of income is positive with no visible evidence of a trend (i.e., response moves ups and downs) until it became negative after the sixth horizons. The response of agriculture to innovation shocks is positive over the periods of horizons, while that of biocapacity is characterized by upward and downward movements over the entire horizons. Finally, the response of energy use to innovation is positive up to the fourth horizon. However, between fourth and sixth horizons, the response turned negative and consequently crossed to the positive region in the mid-sixth horizons. The results have validated the causality we have found between the variables.

Counterfactual Change
The traditional ARDL estimation procedure produces in-sample parameters that often complicate for statistical inferences. The novel dynamic ARDL simulations technique was developed by Jordan and Philips (2018) and utilized in the seminal work of . The versatility and policy usefulness of the estimation method has been applied in several disciplines (Owusu and Sarkodie, 2020b;Sarkodie et al., 2020;Shabbir et al., 2020). Thus, we utilized the novel dynamic ARDL simulations to examine the outsample effects of counterfactual shocks in exogeneous independent variable at a given time period. This is appropriate to examine how CO 2 emissions will respond to future shocks from a specified exogeneous regressor. The counterfactual shocks observed in Figures 5A,B reveals that −1% change in predicted income has no potential effect in the first 9 years but a 1.4% positive rebound effect of CO 2 emissions is observed in the 10th year and stabilizes from the 13th year and thereafter. Contrary, a 1.4% negative rebound  effect is observed after 10 years of 1% shock in predicted income but turns steady after 13 years. This implies that wealth has a long-term mitigating effect on CO 2 emissions-corroborating the notion of pollute in poverty, clean when wealthy. A similar trend is observed in Figures 5C,D, however, −1% shock in predicted energy utilization leads to over 5% decline in CO 2 emissions from the 10th year but reaches a steady-state in the 14th year and afterward. In contrast, over 5% increase in CO 2 emissions is observed after 9 years of no impact at a 1% change in energy consumption. Thus, energy intensity has an escalation effect on CO 2 emissions, which calls for energy efficiency, management, and conservation options to decarbonize energy utilization. Unlike income and energy utilization, there is evidence of very little impact (∼0.001%) observed in CO 2 emissions at −1% change in predicted agricultural innovation in the first 10 years. However, a sharp decline of CO 2 emissions by 0.48% is noticed in the 11th year but a "noisy" effect of reduction in CO 2 is observed thereafter. A 1% shock in predicted agricultural innovation escalates CO 2 emissions by 0.47% from the 11th year and afterward ( Figures  6A,B). The relatively low impact of agricultural innovation on CO 2 emissions compared to energy utilization-can be attributed to vintage agricultural practices, mechanization, and technologies for value addition. This implies that the implementation of modernized and sustainable agricultural process may have a long-term effect, leading to a decarbonized agrarian economy. Regeneration of the ecosystem plays an essential role in reducing carbon footprint. We observe in Figures 6C,D that −1% shock in predicted biocapacity increases CO 2 emissions by almost 1% after 9 years whereas 1% change in predicted biocapacity declines CO 2 emissions by ∼1%. Implying that a reduction in ecological footprint will improve the regenerative capacity of the ecosystem, hence, reducing emissions in the long-term.

CONCLUSION AND POLICY IMPLICATIONS
Agrarian-based economies are often characterized by natural resource exploitation and ecological degradation. In this study, we assessed the impact of agricultural innovation and biocapacity on carbon-based emissions and tested the validity of the EKC hypothesis in Nigeria. We applied several estimation techniques including the novel dynamic ARDL simulations, with data from 1981 to 2014. The empirical evidence based on the ARDL procedure confirmed the long-and short-run validity of the EKC hypothesis for Nigeria. We found that agricultural innovation and energy utilization escalate the levels of anthropogenic CO 2 emissions. In contrast, an expansion of the regenerative capacity of the ecosystem is found to decline the outgrowth of CO 2 emissions. The causality revealed that agricultural innovation has strong predictive power on CO 2 emissions. Similarly, we found income level to predict long-term energy utilization. The results indicated that CO 2 emissions predict income, biocapacity and energy use while feedback causality occurs between agriculture and biocapacity and again agriculture and energy use. These findings were validated by the variance decomposition and impulse-response function analyses. Particularly, we found that 22.79% of agricultural innovation corresponds to 49.43% CO 2 emissions. We also found that 5.95% biocapacity caused 35.78% CO 2 emissions, while 1.61% of energy use led to just 16.27% CO 2 emissions. In contrast to the insample estimation techniques, the counterfactual shocks from the novel dynamic ARDL simulation showed favourable mitigation effects of income, and biocapacity on CO 2 emissions whereas escalation effects of agricultural innovation and energy utilization were also noticed. These findings demonstrate that improvement in livelihoods, environmental awareness creation, and prioritization of ecosystem management and restoration will have a long-term effect on environmental sustainability. Therefore, based on these findings, to achieve carbon sequestration in Nigeria, there is a need for sustainable and clean energy policies that have low environmental damaging effects. Such policies would discourage the excessive use of traditional biomass and fossil fuels. Second, decarbonizing agricultural innovations will include sustainable agricultural practices that encourage clean and renewable energy utilization for agricultural activities. For example, solar energy can be extended to greenhouse heating and cooling, product drying and lighting, in addition to irrigation in the farm field. More so, for improvement in soil, greenhouse, and barns, as well as heating the soil and drying agricultural products, geothermal can be applied. Furthermore, bioenergy can be used to power machinery whereas wind and hydro can be used to generate electricity, irrigate, and process crops. In addition to the above, we suggest a number of policy instruments to mitigate carbon dioxide emissions in Nigeria. These instruments include the use of fiscal instruments such as taxes and fees. We also suggest that financial instruments like subsidies can also be used to regulate the behaviors of the polluters. The coercive power of the state can be applied on polluters who go beyond the avoidable levels of pollution as prescribed by industrial emission standards.
To this end, future studies can shift from theoretical investigation to undertake a randomized controlled trial that examines the effect of sustainable agricultural practices on CO 2 emissions. Such experimental studies will improve the global debate on emissions.

DATA AVAILABILITY STATEMENT
Publicly available datasets were analyzed in this study. This data can be found here: https://databank.worldbank.org/ reports.aspx?source=World-Development-Indicators. All scripts from the estimation method are available upon reasonable request.   The filled-black circle denotes the predicted shock whereas sunflower lime, emerald and cranberry colored spikes refer to statistical significance at 25, 10, 5% level.