Abstract

Upholding sustainability in the use of energies for the increasing global industrial activity has been one of the priority agendas of the global leaders of the West and East. The projection of different GHGs has thus been the important policy agenda of the economies to justify the positions of their own as well as of others. Methane is one of the important components of GHGs, and its main sources of generation are the agriculture and livestock activities. Global diplomacy regarding the curtailment of the GHGs has set the target of reducing the levels of GHGs time to time, but the ground reality regarding the reduction is far away from the targets. Sometimes, the targets are fixed without the application of scientific methods. The aim of the present study is to examine sustainability of energy systems through the forecasting of the methane emission and agricultural output of the world’s different income groups up to 2030 using the data for the period 1981–2012. The work is novel in two senses: the existing studies did not use both the Box–Jenkins and artificial neural network methods, and the present study covers all the major economic groups in the world which is unlike to any existing studies. Two methods are used for forecasting of the two. One is the Box–Jenkins method, where linear nature of the two variables is considered and the other is artificial neural network methods where nonlinear nature of the variables is also considered. The results show that, except the OECD group, all the remaining groups display increasing trends of methane emission, but unquestionably, all the groups display increasing trends of agricultural output, where middle- and upper middle-income groups hold the upper berths. The forecasted emission is justified to be sustainable in major groups under both methods of estimations since overall growth of agricultural output is greater than that of methane emission.

1. Introduction

From the last half of 19th century to till date, economic growth turns into the most important particle of almost all socioeconomic systems in our mother earth. To achieve the higher growth trajectory, each and every economy put all of their resources on the board without giving any potentiality to future generations. It is only in late 90s, when scarcity of resources and a relatively new term “global warming” knocking the door of the policy and law makers around the world, human beings push forward the agenda of sustainability. In the wake of the issues related to sustainability, researchers are often engaged themselves in a debate over the existence of whether substitutability or complementarity are working between the association of growth and environment [1, 2].

It has been historically evidenced that growth can revolutionize the structural changes in both production and consumption. Such changes may occur from either directions or both, that is, either from level or composition or from both of them [3]. Interestingly, both the level and the composition of production and consumption activities affect environmental degradation and raise the scope of greenhouse gas (GHG) emissions, owing to which the prospects of sustainable economic development may hamper in future [4, 5]. It is evidenced that GHG contributes global warming and, consequently, generates severe environmental matters. It is to be noted that, to control the global emissions of GHG, the Kyoto Protocol was proposed and signed by almost all the countries in the world. The Kyoto Protocol specified six GHGs, including methane (CH4), carbon dioxide (CO2), nitrous oxide (N2O), perfluorocarbons (PFCS), hydrofluorocarbons (HFCS), and sulphur hexafluoride (SF6) [6]. In 2014, the concentration of carbon dioxide (CO2), methane (CH4), and nitrous oxide (N2O) in the atmosphere was 397.7 ppm, 1833 ppb, and 325.9 ppb, respectively (World Meteorological Organization, 2015). On the average, the anthropogenic emissions grew 1.3% annually from 1970 to 2000 and 2.2% annually from 2000 to 2010 [7]. Moreover, after carbon dioxide, methane is the second most emitted GHG; its potential to catch heat in the atmosphere is 23 times higher than carbon dioxide [8] and so a clinical examination on increase in methane gas emission needs more attention.

Under 1996 IPCC revised guidelines, national GHG inventories includes energy, industrial process, solvent and other products, agriculture, land-use change and forestry, and waste, while the above-stated list is modified under 2006 IPCC guidelines [6, 9, 10]. Following the just-stated segregation of GHG, methane emissions are also generated from several production sectors. For instance, anthropogenic methane is emitted from sectors like cattle breeding, rice cultivation, extraction and transport of fossil fuels, and waste management [11]. These emissions result from very heterogeneous processes with several scopes for abatement. Accordingly, existing heterogeneity of production structures across countries introduces cross-country asymmetries broadly based on agriculture-based methane or industry-based methane emissions. Interestingly, methane emission and sectoral composition are rarely analyzed in the literature. However, such gap is widened enough in case of agricultural methane emissions. Methane is produced and emitted from the decomposition of livestock manure and the organic components in agro-industrial wastewater. These wastes are typically stored or treated in waste management systems that promote anaerobic conditions and produce biogas, a mixture of about 70 percent methane, 30 percent carbon dioxide, and less than 1 percent hydrogen sulfide. Globally, manure management added an approximated 237 million metric tons of carbon dioxide equivalent of methane emissions in 2010, roughly 4 percent of total anthropogenic methane emissions. Out of total emitted agro-based methane, almost 85 percent is accompanied by USA, China, and India together, followed by Brazil, Pakistan, and Vietnam [12]. It is to be noted that the agriculture methane may emit also from nations which use more capital-intensive production technique, and hence, a critical analysis between agriculture and methane emissions is needed abruptly.

Amalgamation of methane emission with agriculture production creates doubt over the efficacy of sustainability issue. Massive agriculture production can emit more vulnerable methane along with other GHGs. Again, such methane emissions may affect weather variability and multiply climate change risks and the magnitude of global warming. This can affect dairy cattle feeding sector along with other agriculture-based activities more severely. As a consequence, the vulnerability of agriculture-based livelihoods may increase with induced disaster risks. However, there is no definite and robust model which can estimate social costs from such emissions [13]. Hence, by minimizing environmental degradation and pollution risks along with adaptation to climate and weather, variability risks should not only increase resilience of farmers’ production systems but also stabilize their output and income [14]. Identification and reduction of above-stated uncertainties and risk factors in terms of anticipatory adaptation may raise the potentiality of sustainability paradigm [15]. Therefore, climate change adaptation policies in the agricultural sector along with adaptation to control methane emission are to be implemented for getting sustainable development. Therefore, the question still remains in mind: “does complementarity between methane emission and agriculture production generate sustainability?” This paper also tries to locate, screen, and evaluate this issue for major income groups of the world.

This paper contributes original findings concerning methane emission and agriculture production with special emphasis on sustainability. First, it goes for forecasting of methane emission and agricultural output using the Box–Jenkins (BJ) and artificial neural network (ANN) methods. Second, it goes for testing the sustainability of methane emission vis-à-vis agricultural output.

The paper is organized as follows: literature review is presented first, followed by data, methodology, analysis of results, and conclusion.

2. Literature Review

Table 1 exhibits the brief information on the highly relevant works reviewed so far for the present study.

Analysis related to GHG emission and economic activities are not new in the literature. Study related to GHG emission and economic growth has been discussed in the literature quite rigorously [1619]. All these studies used the similar kinds of methodology to relate GHG emission with growth. In fact, these studies have used CO2 as a measure of GHG and per capita income for panel data to show the presence of EKC. Again, there are a few studies that have used several GHGs, and they have confirmed the existence of EKC for methane emissions [2022]. In this context, using a dataset for 22 OECD countries, it has shown a quadratic relationship between methane emission and GDP in the long run [20]. Such quadratic relationship between methane emission and GDP has also been established in the literature for different datasets [22]. Again, industrial methane emission of 39 countries explicitly claims N-shaped relationship between the methane emission and economic growth [21].

In a notable working series titled “OECD Environmental Outlook to 2030,” it studies the prediction of GHG emissions in 2030 if the present inaction on environment remains unchanged [24]. The report reveals that, by 2030, the world economy is expected to nearly double and world population to grow from 6.5 billion today to over 8.2 billion. Most of the growths in income and population will be in the emerging economies of Brazil, Russia, India, Indonesia, China, and South Africa (the BRIICS) and in other developing countries. Rising income and aspirations for better living standards will increase the pressure on the planet’s natural resources. In another series titled “OECD Environmental Outlook to 2050,” it envisages that, without more ambitious policies than those in force today, GHG emissions will increase by another 50% by 2050, primarily driven by a projected 70% growth in emissions from energy use [23]. This is primarily due to a projected 80% increase in global energy demand. Furthermore, it claims that, historically, although OECD economies have been responsible for most of the emissions, in the coming decades, increasing emissions will also be caused by high economic growth in some of the major emerging economies. Again, global energy-related carbon dioxide () emissions are projected to increase by one-third between 2012 and 2040. The continuing increase in total emissions occurs despite a moderate decrease in the carbon intensity ( per unit of energy) of the global energy supply [40].

From the international trade angle, a few studies have found positive relationship between trade openness and emissions [25, 29, 30]. Positive association between methane emission and trade openness is also acknowledged in the literature [26]. Moreover, economic growth and several socioeconomic activities are claimed responsible for methane emission [25]. In fact, rapid growth, population size, and foreign direct investment are made as the responsible factors behind methane emission for different cross-sections [30]. Again, through an interesting study, it is reported that the elasticity of methane emissions with respect to income per capita income is low and it may decrease over time [28]. In a recent study based on country specific efforts, it has been calculated by neural network method that the predicted methane emission from wastewater in China will be an increasing trend and a spatial transition of industrial wastewater emissions from eastern and southern regions to central and southwestern regions and from coastal regions to inland regions will occur [27].

Again, some studies are focused on the reduction potentiality of methane emission from agriculture [31, 41]. Usually, such studies have used static methodology and derived short-run estimates to locate the reduction possibility of methane emission from agricultural sector [31, 41]. However, long-run analysis has also been established, in which methane emissions from agriculture and reduction prospective under several marginal abatement costs, huge drop likely regions, and emission sources are claimed for long-run [32]. Again, establishing the significance of agricultural sector in the context of GHGs emission, it has been claimed that the anthropogenic methane emissions are mostly produced by a few economic sectors such as cattle breeding and rice cultivation [11].

Issues related to sustainable development in the context of agriculture production and methane gas emission have been discussed critically in the literature. Furthermore, it is argued that proper policy investigations, plans, programmes, and adaptation in terms of risks and opportunities can make GHG emissions as sustainability indicators to uphold sustainable development [34, 42]. To get sustainability, investigators are usually advocated for the attractive adaptation measures to pursue efficiently in long run [37]. Investigation in this aspect has revealed that, by controlling environmental degradation and pollution risks along with adaptation to climate and weather, variability risks may increase resilience of farmers’ production systems and also side by side stabilize their output and income [14]. In another series titled “Sendai framework for disaster risk reduction 2015–2030,” it has been argued that minimization of uncertainties and risk factors owing to climate change attached to agriculture can be optimized through anticipatory adaptation [15]. Sustainable development in terms of improvement of farmer’s livelihood is claimed and argued in favour of proper adaptation of changing policy regimes in the context of environmental degradation to opt sustainability [36]. Again, with inability to screen, evaluate, and treat risks augmented in dairy feeding, adaptation initiatives are declared as responsible factors to enhance risks embedded in climate change. Furthermore, studies claim that just-stated failure to adopt proper adaptation may aggravate small-scale farmers’ vulnerability to climate change and weather variability, and in return, economy will produce suboptimal outcomes [33, 35, 38]. In a more recent study, it has been claimed that methane gas emission along with other GHGs emissions from agriculture production and in dairy feeding strategies can be used as a measure and indicator of sustainability. It has been further argued that policy implementation to curb risks associated with agriculture production owing to methane emission must be embedded with the cognizance of poverty, maladaptation, and environmental degradation nexus [39].

3. Rationale of the Present Study

The review of literature highlights different aspects of GHG emission in general and methane emission in particular and their impacts in different sectors of different economies but does not cover studies related to forecasting of methane emission in world’s leading methane emitting countries. The present study has tried to fill the gap in the literature by means of forecasting methane emission for world’s leading economic groups up to the year 2030. Furthermore, the sustainability of the forecast values of methane emission has been analyzed by means of forecast values of agricultural output of the same economic groups. It is thus a novel work in our view.

3.1. Data

The study uses the time series data on methane emission (in kt equivalent) for the five groups of economies (high income, upper middle income, middle income, lower middle income, and low income) for the period 1981–2012. It also uses the time series data for the same period and same groups of economies on the total agricultural value added measured in current USD. Both the data series are borrowed from the World Bank (http://www.wordbank.org).

3.2. Methodology

Twin methods, not actually hybrid in usual sense, are used for forecasting of methane emission and agricultural value added. One is the Box–Jenkins method, where linear nature of the two variables is considered, and the other is artificial neural network methods, where nonlinear nature of the variables is also considered.

3.3. Box–Jenkins Method of Forecasting

Before going into the details of Box and Jenkins method of forecasting, we need to see how a time series data of a particular variable is generated.

There are three processes behind generation of a time series data:(1)AR process: past values of the variable and error term generate the data(2)MA process: only the errors or the disturbance term generate the data(3)ARMA process: data are generated by the combination of AR and MA processes

Sometimes, it is taken as ARIMA model, where “I” stands for integration of the series or how many differencing is done for making the time series of the variable to stationary.

In the AR (p) process, the current value of a variable “y” depends on only the past values plus an error term. If there are “p order in the process, i.e., the current value of y depends on the p order of past values and an error term of the current period,” then the AR(p) can be written aswhere ut is the white noise (WN) error term with zero mean, constant variance, and zero auto-covariance.

On the contrary, an MA(q) process is the linear combination of all the “q” terms of white noise terms depending on time. It is a white noise process in which the current value of yt depends on the current value of the WN error term and all past values of the error terms. Because all the errors are WN, an MA process is necessarily a stationary process further because it is the linear combination of all plus and minus values of the errors which hover around zero.

So, an MA (q) process can be written as

An AR process is stationary if the characteristic root lies outside the unit circle or having values >1, then φ becomes less than 1. This means the condition φ < 1 leads to the values lying inside the unit circle representing stationarity of the AR process, and the model will thus have stability property.

An ARIMA (p, q) process is the combination of AR and MA processes, “I” being the order of integration, which can be represented by “d,” number of differencing to convert the series from nonstationary to stationary. The model for ARMA (p, d, q) can then be written as

Using Lag operator, we have

This relation stands for invertibility between the AR and the MA process.

3.4. Forecasting in ARIMA Model: Box–Jenkins Method

The BJ methodology to determine which model is appropriate follows a four-step procedure:Step 1: identification: to determine the appropriate values of p, d, and q.(i)The main tools in this search are the correlogram and partial correlogram.Step 2: estimation: to estimate the parameters of the chosen model.Step 3: diagnostic checking: to check if the residuals from the fitted model are white noise.(i)If they are, accept the chosen model; if not, start afresh.(ii)That is why the BJ methodology is an iterative process.Step 4: forecasting. The ultimate test of a successful ARIMA model lies in its forecasting performance, within the sample period as well as outside the sample period. On the basis of the acceptable results obtained from steps 1 to 3, forecasting is made on the appropriate model of ARIMA. The forecasting results are accepted on the basis of the acceptable values of root mean square error (RMSE), bias proportion, variance proportions, and covariance proportions. The acceptable forecasted values will be those whose RMSE will be minimum possible, and covariance proportions will be greater than bias proportions and variance proportions.

4. Methodology of ANN-Based NAR

Real-world data always contains nonlinearity, and specifically, its behaviour is dynamic and depends on their current period. Under such circumstances, the nonlinear autoregressive (NAR) neural network structure is effective to make efficient prediction about future [43]. The first advantage of NAR networks is that they can accept dynamic inputs represented by time series sets. Time series forecasting using a neural network is a nonparametric method, which implies that knowledge of the process that causes the time series is not necessary. Moreover, the NAR model utilizes the past values of the time series to predict future values. This fact makes it hard to model time series using a linear model; therefore, a nonlinear approach should be preferred, and the present study has also attempted the method. A nonlinear autoregressive neural network, applied to time series forecasting, describes a discrete, nonlinear, autoregressive model that can be expressed in the following manner [44, 45]:where is data series of x variable at time t; is unknown in advance, and the training of the neural network aims to approximate the function by means of the optimization of the network weights and neuron bias; and is the error of the approximation of x at time t.

This training function is often operated efficiently with backpropagation-type algorithm, and to perform this with our stated, we use Levenberg–Marquardt backpropagation procedure (LMBP) [46, 47] to solve any specified NAR neural network. In Figure 1, we present the topology of a standard NAR network.

After getting the forecasted values of both the series for methane emission and agricultural value added up to the year 2030, we try to test whether the forecasted methane emission is sustainable by means of looking at the forecasted values in agricultural outputs for the selected five groups of economies. For this purpose, we have first calculated the growth of these two indicators over the forecast period and average values of these two growth rates for all the groups. After that, we have tested the mean difference of these two indicators, methane emission and agricultural value added, and tested their significance statistically. If the average growth of agricultural output is tested to be greater than that of methane emission, then it may be said that the emission is sustainable as it contributes positively and largely to agricultural output. The reverse results may say the unsustainable methane emission in the forecasted period.

4.1. Analysis of Results

As mentioned earlier, the study applies Box–Jenkins (BJ) and artificial neural network (ANN) for forecasting methane emission and agriculture output for the period 2013–2030 on the basis of data for the period 1981–2012. The results of both the methods are given one by one.

4.1.1. Forecasting by Box–Jenkins Method

For the BJ method, the following four steps are followed which are mentioned in Methodology:Step 1: identification: to determine the appropriate values of p, d, and q, we have done unit root test through ADF test and correlogram methods. The results for both the series are presented in Tables 2 and 3. Looking at the autocorrelation functions (ACFs) and partial autocorrelation functions (PACFs), we have identified the orders of AR and MA processes. There may be more than one alternative of the shapes of ACF and PACF, and we will have to determine the optimum structure of ARIMA. For this purpose, steps 2 and 3 are followed.Step 2: estimation: to estimate the parameters of the chosen model, we run equation (4).Step 3: diagnostic checking: to check if the residuals from the fitted model are white noise. The acceptable regression results are taken on the basis of where both AR and MA coefficients are significant, adjusted R2 is highest, and information criteria (AIC and SIC) are of lowest values. The results of steps 1 to 3 are given in Table 2 for methane emission and in Table 2 for agriculture output for all the groups of economies. The roots of the AR and MA should lie inside the unit circle, indicating stability of the models.Let us first discuss the results (for steps 1 to 3) on the methane emission with the help of Table 2. The results from the table show that, in all the groups of economies, the series are integrated of order 1. The optimum orders of the autoregressive and moving average terms are marked bold. They are (4, 4) for the OECD and lower middle group, (11, 11) for the upper middle group, (6, 6) for the middle group, and (1, 1) for the low group. And all of these terms are less than unity in values, indicating the stability of the models.Now come to the discussion on the results (for steps 1 to 3) of agriculture output with the help of Table 3. The results from the table show that, in all the groups of economies, the series are integrated of order 1. The optimum orders of the AR and MA terms are marked bold. They are (2, 12) for the OECD group, (4, 1) for the upper middle group, (3, 1) for the middle group, (3, 3) for the lower middle group, and (1, 1) for the low group. And the models in all the groups are stable.Step 4: forecasting: on the basis of the acceptable results obtained from steps 1 to 3, forecasting is made on the appropriate model of ARIMA. The forecasting results are accepted on the basis of the acceptable values of root mean square error (RMSE), bias proportion, variance proportions, and covariance proportions. Figures 2 and 3 present the graphical plots of forecasted values of methane emission and agricultural output, respectively. The numerical values of the two forecasted series are given in the Appendix (Tables 4 and 5).

It is observed from Figure 2 that, except the OECD group, all the remaining four groups of economies demonstrate rising trends of forecasted values of methane emission. Middle-income group leads the club followed by the upper middle-income group and low-income group. The lower middle-income group maintains a constant forecasted path for the entire period of prediction. The positive improvements are observed only for the countries in the OECD group. The results thus show that the agriculture activity in particular and all the economic activity in general is not going to put pressure on the pollution level measured by methane emission for the developed countries, whereas the countries in the remaining world are going to pollute the environment. The derived forecasted values of methane emission have maintained the desirable properties of forecasting as their RMSE is low and the covariance proportions are greater than the bias and variance proportions (the results are not shown to avoid crowding of figures and tables in the text).

On the contrary, the forecasted values of agricultural value added for all the groups under the BJ method, as depicted in Figure 3 and Table 5, show rising trends for all up to 2030. But the difference is observed in their relative positions. The middle-income group is at the top, followed by the upper middle-, lower middle-, OECD, and low-income group. The rate of growth is steeper for the middle-income group as well.

4.1.2. Forecasting by ANN-Based NAR

In the ANN method, only one hidden layer has been used while number of neurons in hidden layer has been varied at four levels (5, 10, 15, and 20 number of neurons). Our experiments suggest employing 2 feedback delays of the variables for model building. Here, we have used backpropagation algorithm proposed by Levenberg–Marquardt for training. Figures 4 and 5, respectively, present the predicted values of methane emission and agriculture output. The quantitative figures for both are presented in the Appendix (Tables 6 and 7).

It is observed from Figure 4 that the OECD group demonstrates falling trend of the predicted methane emission, while the other four groups from low- to upper middle-income countries produce rising trends of the said emission. Furthermore, it is to note that the results under the ANN method are similar to that under the BJ method.

Figure 5 depicts that all the groups’ predicted trends of agriculture output are upward rising over time which are alike to that under BJ method. But a little difference under ANN is observed for the OECD group as it turned downward trend after 2015.

Hence, the two methods of forecasting by and large produce the same results for both methane emission and agricultural output for all the groups of economies. Whatever differences observed are due to the differences in methodological structures. As mentioned in related literature [4850] that the ANN is applied for linear and nonlinear data and BJ only for linear data, the former one can be used as better predictor for a dynamic variable like methane emission. As having association between methane emission and agricultural output, it is now required to examine whether the predicted methane emission is sustainable for all the groups for the period 2013–2030. This is the second objective of the study.

One way of examining such sustainability is to see whether growth of the predicted agricultural output is greater than that of methane emission. In other words, whether good economic effect is greater than bad pollution effect. For the said purpose, we have calculated the average growth rates of predicted methane emission and agricultural output and took their difference and test statistically (by t-test) whether such mean difference is positive statistically. We have done these tests for all the groups of economies separately for BJ and ANN results. The test results are given in Table 8.

It is observed from both the two methods of forecasting that the average growth for predicted methane emission is negative for the OECD group and positive for all the remaining four groups. Furthermore, the average growth of agricultural output is greater than that of the methane emission for all the groups of economies. The correlation between the growth of methane emission and that of agriculture output is positive and significant for all in case of the BJ method, but the correlation result is not significant for all the groups in case of the ANN method.

The results for mean difference test are positive and significant under the BJ method for all the groups which mean that the forecasted values of agriculture output are significantly greater than that of methane emission. This further indicates that the methane emission is sustainable as it does not outweigh the agricultural output. But for the ANN-based results, the significant mean differences are observed for the OECD, upper middle-, and lower middle-income groups which further justify the sustainability of methane emission. The insignificant mean difference results for the middle-income and low-income groups may reveal unsustainable methane emission.

5. Discussion

As mentioned, we have attempted to make forecasting of methane emissions and agricultural value added by BJ and ANN methods and tested sustainability of such emissions vis-à-vis agricultural output for the major economic groups of the world for the time up to 2030. The results for methane emissions are seen to be declining for the OECD group but increasing for the remaining four groups. Referring to last row of Table 4 of the Appendix, the OECD group is expected to reduce the emission by 17.35 percent in 2030 in comparison to its value in 2012. The low-income group is expected to increase their emission levels by 27.12 percent, upper middle-income group by 20.44 percent, and middle income by 18.76 percent, and the lower middle-income group will face lowest emission of mere one percent. The predictions of OECDEO (2011) and USEIA (2016) are a little bit higher (30 percent) than that of the present study [24, 40].

Coming to the prediction of agricultural output, it is observed that all the groups have been showing increasing trends with the middle-income group at top of the list and the low-income group at the bottom with respect to the level values. Referring to Table 5 of the Appendix (last row), it is observed that the low-income group is expected to grow at the rate of 79 percent, middle group by 66.77 percent, upper middle by 61.45 percent, and OECD and lower middle group by 56.56 in 2030 with respect to 2012. All the results of forecasting are derived under the condition that all the associated indicators to methane emission will behave in the same manner in all the future period of prediction.

The sustainability of the methane emissions has been checked by the mean difference tests between the growth rates of the forecasted values of agricultural output and methane emissions. The results are positive and significant under the BJ method for all the groups which mean that the forecasted values of agriculture output are significantly greater than that of methane emission. This further indicates that the methane emission is sustainable as it does not outweigh the agricultural output. But for the ANN-based results, the significant mean differences are observed for the OECD, upper middle-, and lower middle-income groups which further justify the sustainability of methane emission. Hence, it is recommended that, considering all the other factors of forecasting to be unchanged for the forecasting period, sustained agricultural activities may be a better solution which will be viable in economic as well as environmental fronts.

6. Conclusion

In our journey to forecast methane emission and agricultural output of world’s leading groups by two methods, BJ and ANN, it is now to conclude the entire study. Both the methods of forecasting show that, except the OECD group, all the four remaining groups display increasing methane emission, but agricultural output is of increasing trends for all. Middle-income countries possess the top slot in both the methods. So, increase in methane emission is an alarming issue to the global leaders for the sake of environmental sustainability. Furthermore, testing for sustainability of such increasing emission vis-à-vis agricultural output, it is observed that the said emission is sustainable since the average growth rate of the latter is greater than that of the former. Hence, the environmental damage in true sense through methane emission may not be alarming as it boosts up the agricultural growth rate for all the groups. But the effect of methane emission upon other sectors of the economies for examining sustainability in a broader sense remains unverified. It may be kept as the agenda for future research.

Conflicts of Interest

The authors declare that they have no conflicts of interest.