In the era of the Sustainable Development Goals (SDG’s) and the Paris Agreement on Climate Change, the world is struggling to find how to resolve the nexus between the need for economic growth and human happiness to be enhanced while living within environmental limits such as climate change. Three approaches to these issues are common in the past whether commentators come from three disciplinary bases: an economist’s perspective where economic growth is always optimised     , an environmental scientist’s perspective where environmental outcomes are always optimised   , and a human wellbeing perspective where that is always optimised   .
Despite the sustainable development discourse which tries to integrate all three, there is a strong reaction in the literature that keeps these three strands separate with little attempt to resolve how they can be integrated. The gap between these approaches seems to be widening; for example, a special volume on “Absolute Reductions” in the Journal of Cleaner Production  emphasised the need for a radical socio-technical transformation that can bring material, energy and emissions within ecological limits. In order to do this, economic activities were seen to be at the heart of the problem as consumption drives all the problems and this is increased by economic growth  . This is the degrowth movement and below this is expanded to show it is being described as the Treadmill of Production Theory.
The alternative approach to economic growth suggests that economic growth can also (eventually) be the driver of changes that lead to decreases in the flow of energy and emissions through technology that becomes more affordable and more politically acceptable as development occurs   . The tool used to describe this approach is the Environmental Kuznets Curve (EKC) hypothesis. This tool is widely used to help explain economic-environmental trade-offs, some studies supporting it while others contrasting it and some others finding mixed results (see Section 2). The EKC will be developed further below as the basis of the Ecological Modernization Theory. This literature does not include social issues in the triple bottom line nexus but now there is an emerging research strand that tests the EKC by investigating the nexus between economic growth and the environmental intensity of human well being (EIWB)   or carbon intensity of human well being (CIWB)  .
This study follows the emerging research strand by applying it to Bhutan. Using a regression analysis, the paper attempts to examine the case of Bhutan to see if the nexus between these three factors can be resolved and shed light on how an emerging nation may need to address the overlap behind these issues. Bhutan has pledged to remain carbon neutral while it pursues its famous Gross National Happiness (GNH) strategy for a more holistic development outcome whilst growing economically to be a middle range income country  . The case of Bhutan has potential to shed further light on this complex nexus because it has clear positions on all three areas of happiness (GNH), environment (carbon neutrality) and wealth (GDP).
Bhutan’s three goals are quite clearly specified:
1) Happiness, as defined by the Bhutan GNH, is about maintaining a symbiotic relation with nature and balancing material and non-material components of human well being; this is committed by Bhutan to be continually growing  .
2) Bhutan has pledged to remain carbon neutral in perpetuity which is to have net zero carbon emissions .
3) Bhutan wants as well to be a middle income country by 2023 whilst pursuing its GNH strategy.
This suggests that the goals of Bhutan are complex and challenging—requiring the integration of socio-economic and environmental dimensions which are at the heart of Bhutan’s GNH development philosophy and is increasingly associated with Sustainable Development   . However, until now no empirical studies exist for in-depth Bhutan-specific studies on the nexus between economic growth, human well being and environmental pollution, though a possibility of research on linking subjective well being and ecological footprint has been mentioned . Some studies   have included Bhutan in their cross-national level research on issues to do with development and footprint.
Thus, this paper attempts to fill this gap on Bhutan by using the concept of environmental intensity of human well being (EIWB)    to understand using data the various different conjectures about the nexus between happiness, wealth and the environment. It is hoped that the detailed study of an emerging nation can help the rather conflicted debates that are occurring about these three factors. This will have implications on the contested sustainability discourse especially in light of the Paris Agreement and the United Nations Sustainable Development Goals (SDG) for other emerging nations.
2. Environment and Economic Growth Theories
In the era of the SDGs and the Paris Agreement, the level of GHG emissions is a key environmental issue, which is perhaps the greatest environmental issue being faced by humanity in the 21st century. However, there are two strongly divergent theories about how to explain the relationship: Ecological Modernisation Theory which is based on a strongly positive acceptance of GDP-based economic growth as helping resolve issues of wealth and the environment and Treadmill of Production Theory which sees GDP-based economic growth as fundamentally causing environmental issues through its link to the IPAT formula.
2.1. Ecological Modernisation Theory
Ecological Modernisation Theory (EMT) has been using mainstream economic development theory to include how and when environmental factors are brought into the mainstream of development. EMT was developed by Joseph Huber  and the central tenet is to explore institutional reforms and technological advancement to mitigate ecological crises of modernity  . EMT uses the Environmental Kuznets Curve (EKC) as the main tool to explain the evolving relationship between environmental quality and the economic development for a country or countries. The EKC tool was first developed by Panayotou . The EKC suggests that an inverted U-shaped curve best explains how economic growth relates to the environment with a certain minimum in growth being required before environmental outcomes start to peak and then decline  . The EKC literature has been growing dramatically   suggesting the importance of and the interest to understand the environmental impact of economic development with which the world is still grappling.
Studies to test EKC theory were analysed using regression analysis on time series and cross-sectional data. For instance, Song et al.  used an Ordinary Least Square estimation for a cubic log-log model to assess EKC in the provinces of China between GDP per capita and the three environmental degradation parameters of waste water, solid waste and waste gas and found the EKC could explain all three. Similarly, Riti et al.  also found evidence of EKC for CO2 emissions and economic growth in China using various regression techniques. Al-mulali et al.  conducted multi-variant regression to examine EKC in 93 countries using ecological footprint as environmental variable and showed that the EKC holds for high incomes countries but not for low income countries. Similarly EKC was found for 10 Middle East and North African countries  and for a group of 132 countries . EKC was also confirmed by using ecological footprint in place of carbon emissions as the environmental impact parameter .
Other studies have demonstrated six different relationships besides the inverted U-shaped curve between economic growth and environmental pollution . Pérez-Suárez and López-Menéndez  examined the EKC puzzle and the forecasting ability of the EKC equation for 175 countries (including Bhutan). They added a third degree polynomial term of GDP per capita to extend the existing EKC equation and they found different patterns of EKC, while that for Bhutan along with other 10 countries were classified under “other cubic pattern” of EKC, a pattern other than the typical U-shaped, inverted U-shaped, N-shaped, or inverted N-shaped curves. The present study will explore whether Bhutan exhibits this pattern when examining the nexus between CIWB and economic growth.
The EKC has been used to support the need for economic growth in emerging nations as the basis not just of poverty reduction but also environmental improvement, especially if it is connected to urbanization . The evidence from Asia is very supportive of this, particularly China and has been used to generate support for such urbanization-related economic growth. The Asian Development Bank  has shown that Asian countries are significantly earlier in their inverted U-shaped curve and this can be attributed mostly to faster and denser urbanisation; in later years Asia also has higher levels of renewable energy. Similarly, Jamel and Derbali  have shown that urbanisation in selected Asian countries had a statistically significant negative relation to carbon emissions. Zhang et al.  found that the EKC could be established between urbanisation and carbon emissions for 141 countries; this was attributed to urban agglomeration economies that improve the economy and environmental outcomes through economies of scale and density. Similarly, Lin et al.  have shown that economic growth in non-high income countries will decrease carbon emissions and even accelerating urbanisation will lead to just small increases in carbon emissions. Wang et al.  have shown that urbanisation, income and energy efficiency can allow China and India to pursue economic development without growth in emissions.
2.2. Treadmill of Production Theory (TPT)
The TPT was first proposed by Schnaiberg in 1980  , based on political economy whereby a capitalistic system becomes like a treadmill creating social and environmental impacts without integrating ways to deal with such issues. Thus, TPT calls for fundamental changes to the economic system rather than just reform and mostly people in this tradition see GDP as a short hand for capitalism . The TPT also now incorporates a more biological approach developed by a number of scientists pioneered by Ehrlich and Holdren (1992); they began to try and explain environmental impact using the formula I = PAT. This is a combination of Population, Affluence measured as GDP per capita and Technology measured as impact per unit of GDP. GDP is used in this approach to be short for consumption. Quantitative discussions in this approach have been used to show that global impact is likely to continue to increase as long as population and GDP continue to grow  . This has been criticised as not understanding that GDP covers much more than consumption, not being able to distinguish the multiple interactions between PAT, not including social change or political change as factors in environmental impact and not doing justice to spatial factors and agglomeration economies as in the EKC and the EMT approach   .
The EMT suggests that rising economic opportunities will enhance technological advancement and their combination leads to reduced environmental impact once a certain minimum economic growth level is reached. The TPT, however, suggests that economic growth and technological advancement inevitably accelerate resource consumption—the treadmill—leading to further environmental degradation. The two theories are said to be opposing and orthogonal to each other .
2.3. Research Supporting EMT or TPT or Both
Research on the effect of population on carbon emissions found it was not simply population but was dependent on household size and age structure  ; others have found more than a proportional impact of population on carbon emissions   , though Dietz and Rosa  found that carbon emissions tended to decline at a per capita GDP of US $10,000 as would fit with the EKC approach and thus EMT. At national level, Zhang et al.  interpreted EMT for china and its relevance to the rapid socio-economic growth experienced by China in the past decade. At a household level, modest support for EMT was found but strong support for TPT was observed also . For Bhutan, it was found that income and education can lead to choosing cleaner energy  , which in turn can reduce emissions. This supports EMT especially from consumer behaviour and technological choices, though the study was not analysed under the EMT and TPT perspectives. Renewable energy sources make up cleaner forms of energy and it was found to suit the logics of EMT .
At the microeconomic level such as business and firms, the Porter hypothesis   argues that environmental regulation can trigger technological innovation to reduce impact; this approach has led to the phenomenon of decoupling. The decoupling agenda and green growth are closely associated with Ecological Modernisation Theory. Szigeti et al.  found decoupling between ecological footprint and GDP in 90% of the 131 countries examined. Evidence of economic growth decoupling from greenhouse gas emissions were observed globally  . However, Ewing  and Jackson  argue that green growth postulated by EMT should be rejected as a way to restore environmental degradation evident from climate change and biodiversity loss due to rebound effects.
York et al.  employed an Ordinary Least Square analysis under the STIRPAT framework (a flexible version of IPAT), and found a positive effect of both population and economic growth on environmental degradation—supporting TPT and contradicting the expectations of the EMT. They also provided data on the eco-efficiency for 142 nations ranging between 2.76 to 0.52 with smaller values being more efficient and that for Bhutan was reported as 0.75. Similarly, Dietz and Rosa  reported the technology multiplier for 111 nations ranging between 0.24 to 5.75 and that for Bhutan as 0.47, the smaller value indicating a more environmentally benign socio-economic system. Jorgenson and Clark  also found a positive effect of economic growth and population on total carbon emissions and per capita carbon emissions in their cross-national study of 86 countries. Interestingly, they found relative decoupling between economic growth and total carbon emissions but no relative decoupling was observed between economic growth and per capita carbon emissions. Thus, partly supporting both EMT and TPT. This suggests that the unit of measurement can affect the analysis and thereby the outcome of the analysis.
The analyses so far have all been about the nexus between economic growth and environmental outcomes. Now there is an emerging research strand that has applied the EKC hypothesis to incorporate variables that can capture human well being. These are discussed in the next section.
3. Environment, Economic Growth and Wellbeing
Applying a regression analysis for cross-sectional data for 107 countries, Jorgenson  found that economic development leads to increase in carbon intensity of human well being (CIWB) and hence greater environmental impact. Similarly, using a cross sectional data for 135 countries, Dietz et al.  showed that increasing environmental impact does not necessarily lead to greater human well being and that the effect of affluence on human well being characterises a diminishing return. A research at province level in China found that economic growth has a positive effect, while technological innovation has non-linear and negative effects on CIWB . This suggests that technological advancement can reduce the environmental stress generated to produce human well being. This is consistent with the expectations of EMT. Similarly, using cross-national data on GDP per capita, subjective well being index and ecological footprint for 105 countries, Knight and Rosa  supported the TPT but not the EMT.
Thus from Sections 2 and 3, the existing literature shows that the debate between the Treadmill of Production Theory and Ecological Modernisation Theory remains unresolved. For that matter, the EKC remains unsettled, while it stands as an economic tool for a progressive society. However, there is a clear link between the degree of urbanization and the rate at which the EKC turns around the inverted U-shaped curve  .
The question in this paper is therefore whether the EKC theories of EMT and TPT can explain the empirical relation for Bhutan that has a clear policy for the three dimensions of sustainability highlighted in Section 1.
In line with the existing literature, a multivariate regression technique is used to investigate the relationship between human well being, environmental impact and economic growth, arguably capturing the three dimensions of Sustainability. The regression analysis is conducted using EViews, a sophisticated econometrics tool developed with emphasis on time series analysis  . This paper used the method deployed by York et al.  , which was based on the STIRPAT framework, a stochastic version of the IPAT formulation. Robust least square MM estimation1 was also used since there were data outliers. The MM estimation was first used by Yohai in 1987 . The MM estimation and the Robust Efficient Weighted Least Square (REWLS) were found to outperform other robust estimation methods .
Considering the carbon neutral policy of Bhutan, the present study specifically uses the carbon intensity of human well being (CIWB) in place of the generic EIWB. The CIWB is the ratio between CO2eq per capita (denoted as CO2pc) and life expectancy (as LE) was used along with population (as pop), and GDP per capita (as GDPpc) representing economic growth. The details on the choice of these variables are discussed in the next section.
4.1. The Choice of Variables
The type of variables used for evaluating the environmental intensity of human wellbeing vary depending on the data availability and scope of the study. For instance, some studies used CO2 emission per capita    , others have used ecological footprint per capita  . Similarly, for human well being, life expectancy was used as a proxy variable in some studies    , while others have used a life satisfaction index . The choice of variables for the present study was limited by the availability of longer time series data. A longer life expectancy represents healthy living which arguably indicates physical fitness, health services, good food and all other social embedding that impacts longevity. While GDP per capita as a measure of economic growth is well known, the per capita carbon emission is of interest for a carbon-constrained world and for carbon neutral Bhutan. In the present study, per capita CO2eq emissions includes GHG emissions from the energy sector, cement production, and from farming and agriculture, thereby capturing most anthropogenic sources of CO2eq emissions compared to most of the existing research that has accounted for emissions from energy and cement production only   . The time series data used in the present study were sourced from the World Bank .
4.2. The Regression Model
The method based on the STIRPAT framework used in earlier studies    was employed to formulate a functional form specified in Equation (1). STIRPAT is flexible and it is being used for hypothesis testing and also the elasticity between the variables can take any value—not restricted to unity as in the IPAT.
where, CIWB is the ratio variable between CO2pc and LE adjusted for coefficient of variation to reduce the dominance of one variable driving the ratio  . C is the intercept. lnGDPpc is the per capita GDP in natural logarithmic form and the polynomial terms of GDPpc are centred by subtracting the mean of lnGDPpc to reduce collinearity   . It is also logical to add the lag term of the dependent variable as one of the regressors in a time series data analysis-lnCIWB(−1) for the present study. POP is the population and ε is the error term that captures those unobserved explanatory variables. Finally, β1 to β5 are the coefficients of the corresponding explanatory variables that are of primary interest in this study, which are to be estimated by regressing Equation (1) using EViews.
On regressing Equation (1) it was found to have serial correlation and very high variance inflation factor (VIF) and also the linear term and the quadratic term of the GDPpc were found to be less significant than the cubic term despite being jointly significant as per the Wald test. This led to constructing two reduced form models: without the quadratic term (Equation (2)) and without the linear term (Equation (3)).
Small models with fewer parameters were preferred over large models to reduce parameter uncertainty . Between 2) and 3), we chose model 3) which scores relatively better in most of the model selection criteria. It has higher adjusted R2 and lower Akaike info criterion (AIC), Hannan-Quinn criteria (HC); lower VIFs, and lower SE of regression and standard error for the estimated coefficients. In Equation (3), the quadratic- and the cubic-term were jointly significant at 5% significance level as per Wald test though the quadratic GDPpc was significant at 25% significance level only. The joint significance is important here as argued by Brown and McDonough . Moreover, this is within the significance level up to 25% argued by Agung  for cases where a theoretically important variable became statistically insignificant. Furthermore, statistical analysis is considered partly a science and partly an art  requiring trial and error methods. In keeping with this, interpretation of the results in the next section requires caution and it is valid within the data set used in this study.
4.3. Statistical Tests
For regression analysis, stationarity of data is essential to avoid any spurious regression. In this study, any presence of unit root for all the variables were tested by using the Augmented Dickey Fuller unit root test and Kwiatkowski-Phillips-Schmidt-Shin (KPSS) stationarity test. They were found to be non-stationary at levels but stationary in first difference. However, if the residual from their regression is stationary, then their combination can be considered stationary   . This condition was fulfilled in this study—that is the regression results were not spurious.
Statistical tests were conducted on Equation (3). The VIFs were below 16  and it is not reasonable to remove variables just to improve the VIF . The tests for the presence of serial correlation and heteroskedasticity were conducted using the Breusch-Godfrey (BG) LM test and Breusch-Pagan-Godfrey (BPG) test respectively  . The p-values failed to reject the null hypothesis of no serial correlation and homoskedasticity. The robustness of Equation (3) was also checked through a robust least squares MM estimation confirming the signs of the explanatory variables with small differences for the values of the estimated coefficients and their associated standard errors.
5. Results and Discussions
The estimated equation is provided in Equation (4) with standard errors in parenthesis. It has an overall goodness of fit (i.e. R2) is 0.72.
The regression result is shown in Table 1.
The estimated equation illustrates that the relation between the CIWB and GDPpc are positive and curvilinear and that with population is negative. The relationship between CIWB and GDPpc is shown in Figure 1.
The pattern shown by Figure 1 appears to concur with the ‘other cubic pattern’ observed in earlier study  considering that the cubic GDPpc term in Equation (4) is highly significant. Equation (4) suggests that Bhutan is both supporting the TPT (the GDP result) and the EMT (the population result). These findings are elaborated in the following sections.
5.1. Discussion on the GDP Result
The result of GDPpc causing greater carbon intensity per unit of well being suggests that GDP is not (yet) helping with the social and environmental issues associated with development in Bhutan as outlined in the TPT. The coefficient of the quadratic GDPpc was comparable with that reported in a cross national data . Thus GNH or at least the existing socio-economic pattern appears to be leading to greater carbon intensity in Bhutan because enhancing human well being follows a significantly positive and curvilinear relationship with economic growth. And the economic growth has not yet reached a level where the EKC curve would suggest that economic growth will start to reduce carbon intensity. Thus, at this juncture in Bhutan’s growth trajectory, the current per capita GDP for Bhutan is far below the various threshold levels reported in the literature suggested before the EKC turns around on the inverted U curve and begins to decline2. The way ahead for Bhutan which is seeking to resolve the three factors of wealth, happiness and the environment, is therefore quite complex.
The evidence from Bhutan is that the transition to reduced carbon intensity of well being is not yet under way. In other words, the result supports the Treadmill of Production Theory. The result reflects the transitioning of the Bhutanese economy from agriculture-based economy to manufacturing and industrial-based economy  . The latter being more polluting than the former. To avoid environmental degradation in the early stage of development, Panayotou  suggests that environmental protection should be an integral part of development plans especially in developing countries. For Bhutan with its GNH strategy this had always been the objective and it cannot be stressed more especially now with its carbon neutral pledge that Bhutan will need to bend the emission curve before other nations have managed to do this.
Table 1. EViews result.
Figure 1. CIWB plotted against GDP/capita.
One way to examine if the EKC is perhaps beginning in a developing country is to examine if relative decoupling of wealth and environmental impact has begun. Developing countries with low GDP have been found to exhibit relative decoupling . Bhutan’s GDP seems to exhibit a relative decoupling from CO2eq emissions as shown in Figure 2, although this was not the case for carbon intensity of well being. This suggests that there is a need to enhance life expectancy (the denominator of the CIWB) considering that the LE ranges from 45 years in 1980 to just 68 years in 2012. Moreover, Bhutan is still in its early stages of socio-economic development thus struggling to enhance liveability conditions for its citizenry.
Figure 2. Economic growth and CO2eq emissions (based on data source: World Bank .
This calls for greening Bhutan’s socio-economic activities in line with the aspiration of GNH and carbon neutral emphasised in their 12th five-year plan  and their economic development policy  . The need to further green growth in Bhutan was also emphasised by the World Bank . It can be expected that economic growth in middle-income countries will decrease carbon emissions  but it needs action sooner to be consistent with its global commitments to happiness and carbon neutrality.
5.2. Discussion of the Population Result
If the TPT approach had been the dominant story with Bhutan then the coefficient of population being negative in Equation (4) would not have been predicted. The result suggests that growing population will help carbon intensity per unit of well being. This finding may seem spurious to the IPAT identity that considers population as the main driver of environmental degradation   . Furthermore, the finding contradicts the adverse environmental impact of population growth speculated in the national strategy of Bhutan . One way to explain this in the TPT model is to suggest that it is consumption (GDP) that is over-riding the population factor as Toth and Szigeti  found in their studies. However, it is also possible to see within the EMT approach that population growth, if it is in cities, can be an important factor in more rapidly leading to the EKC turn-around that would help Bhutan achieve its multiple goals. Lamb et al.  found this to be the case for population growth in their global study, and for the case of China by Zheng et al.  and Meng and Huang . The plausibility of this is that in Bhutan farming and agriculture activities, which are concentrated in rural areas, are the major contributors of CO2eq emissions for the time period prior to 2000 . Rural Bhutan accounts for 66% of the total population .
Thus as wealth and population grow there is a faster move to the city and then a series of other factors follow that may indeed reduce consumption as economies of scale and density create smaller houses, better public transport and other potential impacts on consumption like greater recycling    . Bhutan expects its urbanisation rate to increase from 10% to 51% in just 55 years compared to 95 years for Asia and 150 years for Europe . This trend can be a positive phenomenon as Bhutan aspires to be a knowledge-based society  while also sustaining its carbon neutral goal. Wang  argues that rather than focusing on reducing population, the aim should be to educate and create awareness of how to reduce carbon emissions through behaviour changes. Similarly, Costantini and Monni  emphasised the need to invest in human and social capital to simultaneously achieve a higher living standard and environmental protection.
From Sections 5.1 and 5.2, it is clear that the case of Bhutan provides support for TPT with regard to GDP as the driver and EMT approach for its population growth that is showing positive outcomes as would be suggested with its rapid urbanization. Furthermore, urban areas were found to be happier  and wealthier   compared to their rural counterparts. In the interim as it moves towards the top of the inverted U-shaped curve there are appropriate policy interventions that can be put in place so well being and environmental quality are not degraded while the economy grows .
Using established econometric techniques, this study has attempted to interlink the relationships between the carbon intensity of human wellbeing and economic growth as a way of resolving how Bhutan is managing to balance happiness, environment and economic growth. It examined the two main approaches that are being used to explain these interactions: Ecological Modernisation Theory based on the Environmental Kuznets Curve, and Treadmill of Production Theory with its basis on the IPAT model. This study has used a 33-year set of data on Bhutan and through vigorous analysis of the time series data on some of the key variables, it has been able to show the following key conclusions.
1) The case of Bhutan appears to support the Treadmill of Production Theory against the Ecological Modernisation Theory as it is showing that GDP is causing the carbon intensity per unit of wellbeing to get worse. However, this is relatively easy to explain in terms of Ecological Modernisation Theory also as it is simply too early in the EKC curve to begin showing the potential reduction in carbon intensity found in more developed places, when examining the nexus between economic growth, human wellbeing and carbon emissions.
2) While the relationship between CIWB and GDP per capita is curvilinear, the relationship between GDP and carbon emissions does not appear to be in lockstep and it seems that relative decoupling is underway. However, there is not yet sufficient economic development for the wellbeing and environmental factors to be more optimised.
3) The case of Bhutan supports the Ecological Modernisation Theory when it examines population growth in that there is a statistically significant negative relation between population and carbon intensity of human wellbeing. It supports the theory of how the EKC turns around after some economic development because it is associated with population growth in cities rather than rural areas and Bhutan is showing all these characteristics.
4) Although Bhutan is still in its early stage of development and aspires to be a middle-income country, it has already pledged to grow in wellbeing and maintain its carbon neutrality. These can only be achieved by decoupling them from the growth in wealth through modifying economic growth to include energy efficiency and decarbonisation policies as well as continuing the commitment to GNH to enable growth in well being.
This research therefore formally places the case of Bhutan into the burgeoning literature on the Environmental Kuznets Curve proposition especially through the lens of the “environmental intensity of human well-being”. It suggests that its current growth trajectory is likely to turn around carbon intensity per unit of wellbeing as it moves more towards being a middle-income country. However, it should not wait to intervene on matters concerning the modification of economic growth in order to ensure its commitments to GNH and carbon neutrality are maintained. While this study is the first of its kind for Bhutan, it shows a promising line of research that may further inquire into the socio-technical structure of the Bhutanese economy that can complement the quantitative analysis undertaken in this study and provide a more nuanced picture of how it can balance the three goals. Similarly, other emerging nations that are yet to reach the top of the inverted U-shaped EKC may also seek direct ways to deal with social and environmental issues whilst maintaining their economic growth. Thus by Bhutan maintaining its carbon neutral goal and its happiness goal along with its economic goal, it is likely to continue to provide a model for the sustainable development discourse in general.
The lead author would like to thank the Government of Australia for providing the Research Training Programme (RTP) scholarship and Curtin University for the research top up scholarship (CRS). The authors would like to thank the anonymous reviewer for the comments.
1MM estimate is a combination of Maximum likelihood (M) estimate and Scale statistic (S) estimate  IHS Global (2014) EViews 9 User’s Guide I and II. in IHS Global Inc. (ed.): IHS Global Inc.
2The average threshold value was found to be United States Dollar (US$) 8910 as per the meta-analysis conducted by  Sarkodie, Samuel Asumadu and Vladimir Strezov (2019) A review on Environmental Kuznets Curve hypothesis using bibliometric and meta-analysis. Science of The Total Environment, 649: 128-145.
 Sachs, J. (2012) Introduction. In: Helliwell, J., Layard, R. and Sachs, J., Eds., World Happiness Report: The Earth Institute.
 Lyubomirsky, S., Sheldon, K.M. and Schkade, D. (2005) Pursuing Happiness: The Architecture of Sustainable Change. Review of General Psychology, 9, 111-131.
 Akenji, L., Bengtsson, M., Bleischwitz, R., Tukker, A. and Schandl, H. (2016) Ossified Materialism: Introduction to the Special Volume on Absolute Reductions in Materials Throughput and Emissions. Journal of Cleaner Production, 132, 1-12.
 Rockstrom, J., Steffen, W., Noone, K., Persson, A., Chapin III, F.S., Lambin, E., Lenton, T.M., Scheffer, M., Folke, C., Schellnhuber, H., Nykvist, B., De Wit, C.A., Hughes, T., van der Leeuw, S., Rodhe, H., Sorlin, S., Snyder, P.K., Costanza, R., Svedin, U., Falkenmark, M., Karlberg, L., Corell, R.W., Fabry, V.J., Hansen, J., Walker, B.H., Liverman, D., Richardson, K., Crutzen, C. and Foley, J. (2009) A Safe Operating Space for Humanity. Nature, 461, 472-475.
 Shuai, C., Chen, X., Wu, Y., Zhang, Y. and Tan, Y. (2019) A Three-Step Strategy for Decoupling Economic Growth from Carbon Emission: Empirical Evidences from 133 Countries. Science of the Total Environment, 646, 524-543.
 Steinberger, J.K. and Roberts, J.T. (2010) From Constraint to Sufficiency: The Decoupling of Energy and Carbon from Human Needs, 1975-2005. Ecological Economics, 70, 425-433.
 Sweidan, O.D. (2018) Economic Performance and Carbon Intensity of Human Well-Being: Empirical Evidence from the MENA Region. Journal of Environmental Planning and Management, 61, 699-723.
 Jorgenson, A.K. and Givens, J. (2015) The Changing Effect of Economic Development on the Consumption-Based Carbon Intensity of Well-Being, 1990-2008. PLoS ONE, 10, e0123920.
 Rogelj, J., Schaeffer, M., Meinshausen, M., Knutti, R., Alcamo, J., Riahi, K. and Hare, W. (2015) Zero Emission Targets as Long-Term Global Goals for Climate Protection. Environmental Research Letters, 10, Article ID: 105007.
 Allison, E. (2012) Gross National Happiness. In: Spellerberg, I., Ed., The Berkshire Encyclopedia of Sustainability: Measurements, Indicators, and Research Methods for Sustainability, Berkshire Publishing Group, Great Barrington, 180-184.
 Pérez-Suárez, R. and López-Menéndez, A.J. (2015) Growing Green? Forecasting CO2 Emissions with Environmental Kuznets Curves and Logistic Growth Models. Environmental Science & Policy, 54, 428-437.
 Dietz, T., Rosa, E.A. and York, R. (2009) Environmentally Efficient Well-Being: Rethinking Sustainability as the Relationship between Human Well-Being and Environmental Impacts. Human Ecology Review, 16, 114-123.
 Spaargaren, G. and Mol, A.P.J. (1992) Sociology, Environment, and Modernity: Ecological Modernization as a Theory of Social Change. Society & Natural Resources, 5, 323-344.
 Mardani, A., Streimikiene, D., Cavallaro, F., Loganathan, N. and Khoshnoudi, M. (2019) Carbon Dioxide (CO2) Emissions and Economic Growth: A Systematic Review of Two Decades of Research from 1995 to 2017. Science of the Total Environment, 649, 31-49.
 Sarkodie, S.A. and Strezov, V. (2019) A Review on Environmental Kuznets Curve Hypothesis Using Bibliometric and Meta-Analysis. Science of the Total Environment, 649, 128-145.
 Song, T., Zheng, T. and Tong, L. (2008) An Empirical Test of the Environmental Kuznets Curve in China: A Panel Cointegration Approach. China Economic Review, 19, 381-392.
 Riti, J.S., Song, D., Shu, Y. and Kamah, M. (2017) Decoupling CO2 Emission and Economic Growth in China: Is There Consistency in Estimation Results in Analyzing Environmental Kuznets Curve? Journal of Cleaner Production, 166, 1448-1461.
 Al-mulali, U., Choong, W.-W., Low, S.-T. and Mohammed, A.H. (2015) Investigating the Environmental Kuznets Curve (EKC) Hypothesis by Utilizing the Ecological Footprint as an Indicator of Environmental Degradation. Ecological Indicators, 48, 315-323.
 Farhani, S., Mrizak, S., Chaibi, A. and Rault, C. (2014) The Environmental Kuznets Curve and Sustainability: A Panel Data Analysis. Energy Policy, 71, 189-198.
 Ulucak, R. and Bilgili, F. (2018) A Reinvestigation of EKC Model by Ecological Footprint Measurement for High, Middle and Low Income Countries. Journal of Cleaner Production, 188, 144-157.
 Jamel, L. and Derbali, A. (2016) Do Energy Consumption and Economic Growth Lead to Environmental Degradation? Evidence from Asian Economies. Cogent Economics & Finance, 4, Article ID: 1170653.
 Zhang, N., Yu, K. and Chen, Z. (2017) How Does Urbanization Affect Carbon Dioxide Emissions? A Cross-Country Panel Data Analysis. Energy Policy, 107, 678-687.
 Lin, S., Wang, S., Marinova, D., Zhao, D. and Hong, J. (2017) Impacts of Urbanization and Real Economic Development on CO2 Emissions in Non-High Income Countries: Empirical Research Based on the Extended STIRPAT Model. Journal of Cleaner Production, 166, 952-966.
 Wang, Q., Su, M. and Li, R. (2018) Toward to Economic Growth without Emission Growth: The Role of Urbanization and Industrialization in China and India. Journal of Cleaner Production, 205, 499-511.
 Gould, K.A., Pellow, D.N. and Schnaiberg, A. (2004) Interrogating the Treadmill of Production: Everything You Wanted to Know about the Treadmill but Were Afraid to Ask. Organization & Environment, 17, 296-316.
 Waggoner, P.E. and Ausubel, J.H. (2002) A Framework for Sustainability Science: A Renovated IPAT Identity. Proceedings of the National Academy of Sciences, 99, 7860-7865.
 Newman, P. (2011) Sustaining Our Future: Resolving the Conflict over Population Models. In: 19th International Congress on Modelling and Simulation, Modelling and Simulation Society of Australia and New Zealand Inc., Perth, 25-37.
 Zhang, Z., Hao, Y., Lu, Z.-N. and Deng, Y. (2018) How Does Demographic Structure Affect Environmental Quality? Empirical Evidence from China. Resources, Conservation and Recycling, 133, 242-249.
 Adua, L., York, R. and Schuelke-Leech, B.-A. (2016) The Human Dimensions of Climate Change: A Micro-Level Assessment of Views from the Ecological Modernization, Political Economy and Human Ecology Perspectives. Social Science Research, 56, 26-43.
 Rahut, D.B., Behera, B. and Ali, A. (2016) Household Energy Choice and Consumption Intensity: Empirical Evidence from Bhutan. Renewable and Sustainable Energy Reviews, 53, 993-1009.
 Curran, G. (2018) Is Renewable Energy Still Green? Shaping Australia’s Renewable Energy Enterprise in an Age of Ecological Modernisation. Environmental Politics, 28, 950-969.
 Ramanathan, R., He, Q., Black, A., Ghobadian, A. and Gallear, D. (2017) Environmental Regulations, Innovation and Firm Performance: A Revisit of the Porter Hypothesis. Journal of Cleaner Production, 155, 79-92.
 Szigeti, C., Toth, G. and Szabo, D.R. (2017) Decoupling—Shifts in Ecological Footprint Intensity of Nations in the Last Decade. Ecological Indicators, 72, 111-117.
 Schandl, H., Hatfield-Dodds, S., Wiedmann, T., Geschke, A., Cai, Y., West, J., Newth, D., Baynes, T., Lenzen, M. and Owen, A. (2016) Decoupling Global Environmental Pressure and Economic Growth: Scenarios for Energy Use, Materials Use and Carbon Emissions. Journal of Cleaner Production, 132, 45-56.
 Jorgenson, A.K. and Clark, B. (2012) Are the Economy and the Environment Decoupling? A Comparative International Study, 1960-2005. American Journal of Sociology, 118, 1-44.
 Greene, W.H. (2001) Econometric Software. In: Wright, J.D., Ed., International Encyclopedia of the Social & Behavioral Sciences, Second Edition, Elsevier, Oxford, 1-7.
 York, R., Rosa, E.A. and Dietz, T. (2003) STIRPAT, IPAT and ImPACT: Analytic Tools for Unpacking the Driving Forces of Environmental Impacts. Ecological Economics, 46, 351-365.
 Lamb, W.F., Steinberger, J.K., Bows-Larkin, A., Peters, G.P., Roberts, J.T. and Wood, F.R. (2014) Transitions in Pathways of Human Development and Carbon Emissions. Environmental Research Letters, 9, Article ID: 014011.
 Brown, S.P.A. and McDonough, I.K. (2016) Using the Environmental Kuznets Curve to Evaluate Energy Policy: Some Practical Considerations. Energy Policy, 98, 453-458.
 Wu, Y., Zhu, Q. and Zhu, B. (2018) Comparisons of Decoupling Trends of Global Economic Growth and Energy Consumption between Developed and Developing Countries. Energy Policy, 116, 30-38.
 Zheng, X., Yu, Y., Wang, J. and Deng, H. (2014) Identifying the Determinants and Spatial Nexus of Provincial Carbon Intensity in China: A Dynamic Spatial Panel Approach. Regional Environmental Change, 14, 1651-1661.
 Meng, L. and Huang, B. (2018) Shaping the Relationship between Economic Development and Carbon Dioxide Emissions at the Local Level: Evidence from Spatial Econometric Models. Environmental and Resource Economics, 71, 127-156.
 IRP (2017) Assessing Global Resource Use: A Systems Approach to Resource Efficiency and Pollution Reduction. In: Bringezu, S., Ramaswami, A., Schandl, H., O’Brien, M., Pelton, R., Acquatella, J., Ayuk, E., Chiu, A.S.F., Flanegin, R., Fry, J., Giljum, S., Hashimoto, S., Hellweg, S., Hosking, K., Hu, Y., Lenzen, M., Lieber, M., Lutter, S., Miatto, M., Singh, A.N, Obersteiner, M., van Oers, L., Pfister, S., Pichler, P., Russell, A., Spini, L., Tanikawa, H., van der Voet, E., Weisz, West, H.W., Wijkman, A., Zhu, B. and Zivy, R., Eds., A Report by the International Resource Panel, United Nations Environment Programme, International Resource Panel, Nairobi, 100.
 Pearson, P. (1994) Energy, Externalities and Environmental Quality: Will Development Cure the Ills It Creates? Surrey Energy Economics Centre (SEEC), University of Surrey, Guildford.