Abstract
India remained the third-largest energy consumer in the world, responsible for around 7% of global carbon emissions due to rising incomes and improving living standards. Although resource extraction has quadrupled since 1970 due to rising population and demand for natural resources, energy use and transformation, notably of fossil fuel energy, have increased by around 45%, thus increasing greenhouse gas (GHG) emissions. In this view, this study aims to explore energy intensity, renewable energy, natural resources, economic growth, and environmental degradation nexus in India. The novel dynamic simulated autoregressive distributed lag and kernel-based regularized least squares (KRLS) approaches are used to explore the effects of energy intensity, renewable energy, natural resources, and economic growth on carbon intensity for India from 1970 to 2020. The empirical results reveal that renewable energy and natural resources improve India's environmental quality via the mitigation of carbon emissions. It is also found that energy intensity and economic growth deteriorate the country's environmental quality by increasing carbon emissions in the short- and long run. A series of robustness estimation affirms the above evidence, thus providing requisite guideline for relevant policy recommendations for the country.
Keywords
Introduction
The concepts of energy intensity, renewable energy, natural resources and carbon intensity are central to environmental sustainability. Energy intensity is an indirect approach of measuring energy efficiency, suggesting that energy intensity is more concerned with energy leanness (quantity) than energy cleanliness. 1 Energy intensity tends to increase as populations increase due to the expected increase in the demand for energy product. According to the Organization for Economic Co-operation and Development (OECD), energy intensity has been examined and recognized as a significant factor in determining the direction of the energy transition towards establishing a low-carbon economy. 2 For example, India’s total primary energy consumption was 659 million tons of oil equivalent (Mtoe) in 2015; it is projected to increase to 1440 Mtoe by 2030. Therefore, affirming that energy efficiency remains a crucial option for advancing the development agenda and combating the effects of environmental degradation and climate change.
Renewable energy is a substitute for non-renewable energy sources, thus making renewable energy sources desirable for sustainable environment and being deliberately utilized to drive energy transition and mitigating environmental deterioration across the globe.3–5 It is noteworthy that as the demand for natural resources and the world's population grows, energy utilization and transformation, particularly of fossil fuel energy, have increased by about 45% even as resource extraction has tripled since 1970, leading to an increase in greenhouse gas (GHG) emissions. 6 Thus, natural resource rent (NRR) plays a vital role in determining the degree of carbon emissions and environmental sustainability. This is because natural resource-rich economies should ordinarily fare better and expand more quickly economically than those with limited natural capital and resources.7–9 However, one of the main limitations of the available literature is that there is almost lack of agreement on whether the NRR has a significant positive or negative impact on environmental sustainability. 10 Furthermore, the degradation of natural resources and pollution issues are global issues, as a result, the problems are intertwined and cannot be resolved separately. 11
The performance of a region's carbon emissions can be evaluated using carbon emissions per unit of GDP, commonly known as carbon intensity. Going by this, India is now the third-largest energy consumer in the world, responsible for around 7% of global carbon emissions due to rising incomes and improving living standards. Although it is still far behind China, the world's greatest emitter, and the United States, India is currently the third-highest producer of carbon dioxide (CO2) on the planet. Widespread wealth disparity (income inequality), unequal development, and other problems persist in India. 12 Specifically, International Energy Agency (IEA) noted that 700 million individuals in India gained access to electriicity since 2000 suggesting the determination of the government to tackle energy poverty while providing secure, affordable, and sustanable energy. 13 As a result, the government placed a higher priority in tackling the problem of poverty than on finding solutions to environmental issues. Therefore, a measure of carbon intensity which is associated with economic growth rather than just focusing on overall CO2 emissions will allow a more accurate evaluation of India's emissions reduction efforts. According to Wang et al., 14 the necessity for emerging countries to reduce their carbon intensity is currently being debated internationally due to the increasing industrialization and urbanization that is being fueled by excessive energy use. 15 Since 2000, energy consumption has increased, with coal, oil, and solid biomass still meeting 80% of the need. As a result, one-third of the world’s GHG emissions have come from coal, oil, and natural gas. Furthermore, with an annual urban population increase of the equivalent of a city the size of Los Angeles, India will soon surpass China as the world’s most populated nation. 16 As the population increases, it is expected that the country will make leaps in its development which will come at the cost of increased energy consumption.
Numerous studies have associated CO2 emissions and carbon intensity with environmental sustainability.17–22 Additionally, exploring carbon intensity alongside energy intensity, renewable energy, and natural resources is not well-reported in the literature especially for India's perspective. It is also essential to study India in light of recent events in climate change, population and developmental projections. While facing severe threats from the changing climate, including rising sea levels, vanishing glaciers, and extreme weather events, India is also committed to a robust energy transition plan that offers net zero emissions by 2070 and achieving 50% share of renewable electricity by 2030.23,24 By using the robust dynamic stimulated autoregressive distributed lag (DSARDL) and Kernel-based regularized least squares (KRLS), this paper makes a novel contribution to the literature by examining the effect of energy intensity, renewable energy, natural resources and economic growth on the carbon intensity of India. Given India’s growing relevance in achieving the sustainable development goals by 2030 and the potential environmental consequences it is anticipated to experience. The next section explains the model, variable description and the methodology employed. Section “Model, variables description, and methodology” contains the empirical result, while the conclusion and policy insight are explained in the last section.
Literature review: Indian perspective
The role of economic growth in environmental quality has been widely investigated in the literature for different cases. Specifically for the case of India, Udemba 25 implemented the autoregressive distributed lag (ARDL) to examine the environmental effect of economic growth (GDP), foreign direct investment (FDI), fossil fuels, and institutional quality over the period 1996Q1–2018Q4. The result from the investigation shows that institutional quality improves environmental quality by mitigating carbon emissions in the short- and long run. Additionally, GDP, FDI and fossil fuels are detrimental to environmental quality because they both causes surge in carbon emission in the country in short- and long-run. This study and its outcomes are also similar to the asymmetric evidence in Akadiri and Adebayo. 26 Specifically, with the implementation of non-linear ARDL for dataset that covers the period, the result reveals that a positive shock on GDP worsen environmental quality by increasing carbon emissions in India. In addition to the above studies, the environmental effect of economic growth and economic freedom has also been documented for India through its membership of the BRICS (Brazil, Russia, India, China, and South Africa) countries.27,28
The role of energy resources and by large natural resources in environmental quality in India have been covered in the literature.29,30 For instance, Itoo and Ali 30 employed the dataset that covers the period 1980–2018 to provide insight on the environmental effect of natural resources and energy aspects in India. By implementing ARDL approach alongside the robustness techniques of FMOLS, DOLS, and CCR (fully modified ordinary least squares, dynamic ordinary least squares, and canonical cointegrating regression, respectively), the study reveals that the depletion of natural resource alongside remittances inflow and industrial output shows a negatively long run but insignificant impact on carbon emissions in India. Also in the long run, GDP, population, and energy consumption negatively affect environmental quality by promoting carbon emission in the country and without validating the environmental Kuznets curve (EKC) hypothesis. In line with natural resources (land, fresh water) and GHG emission nexus, subnational environmental insight was investigated in Damerau et al.. 29 Specifically, the results show that regional usage of cropland and water resources could be up to 50% and 65%, respectively, thus accounting for total natural resources depletion of about 40%. This resource utilization is expected to help India achieve the national food self-sufficiency goals while also achieving GHG emissions reduction of about 26–34% especially through agroforestry emission sequestration.
Given the evidence from the literature on the drivers of environmental degradation in India, the current study extends the literature by implementing carbon intensity as environmental quality indicator while also considering energy intensity, renewable energy, and natural resources which are clearly missing approaches in previous studies. However, in related to other economies across the globe, studies have also captured large array of countries in similar perspectives.31,32
Model, variables description, and methodology
This study empirically analyzes energy intensity, renewable energy, natural resources, economic growth, and carbon intensity nexus in India by using the empirical model outlined in equation (1). This investigated model follows the STIRPAT (Stochastic Impacts by Regression on Population, Affluence and Technology) underlying framework in the earlier works.33,34
where t illustrates time, ln represents the logarithm,
Variables description
For empirical analysis, the paper uses the available annual time-series data ranging from 1970 to 2020. The predicted variable of this paper is carbon intensity, that is, environmental degradation. In contrast to CO2 emissions and ecological footprint, carbon intensity comprises both the environmental and economic consequences of human activities, which makes it an appropriate indicator of environmental degradation. 35 Therefore, the study employs carbon intensity as a proxy for environmental degradation rather than CO2 emissions or ecological footprint.
Additionally, the factor variables of this paper are energy intensity, renewable energy, natural resources, and economic growth. Following Goldemberg 36 and Ahmad and Wu 35 studies, carbon intensity is measured as per capita CO2 emissions per unit of GDP per capita and energy intensity is measured as per capita energy consumption per unit of GDP per capita, respectively. Following Khan et al., 37 total natural resources rents as a percentage of GDP are used as a proxy for natural resources. Finally, for renewable energy and economic growth, we use energy consumption from renewables as kilowatt per hour (kWh) and per capita gross domestic product (GDP) as constant 2015 US dollars, respectively. Detailed information and sources of the study variables are given in Table 1.38–42
Data information and sources.
Note: OWD and WDI are respectively Our wold in data and the World Bank's world development indicator.
By following Namahoro et al., 43 the study variables have been converted into the natural logarithm to perform a robust analysis and eliminate any potential heteroscedasticity problem. Summary statistics of the study variables are visualized in Table 2. As seen in Table 2, EG has the highest annual average and volatility, whereas CINT has the least.
Summary statistics.
Methodology
The autoregressive distributive lag (ARDL) technique of Pesaran et al. 44 has been extensively utilized in the studies examining the factors affecting environmental degradation (see e.g. Refs.45–52). However, as stated by Jordan and Philips, 53 this technique is not yet optimal due to the complicated models’ difficult interpretation of the short- and long-run coefficients. Therefore, Jordan and Philips 53 developed the dynamic simulated autoregressive distributive lag (DSARDL) method based on the ARDL(1,1) model to solve this problem with the traditional ARDL technique.
This innovative method uses stochastic simulations first to estimate, simulate, and provide the short- and long-run coefficients with p-values for statistical inference of the model used. It then automatically plots, while holding all other predictors constant, the impact of positive and negative shocks administered to a predictor at a point in the future time horizon defined by the practitioners on the predicted variable.
53
Overall, the DSARDL technique makes it simple to evaluate and examine the predictors’ present and prospective future impacts on the predicted variable, both in the short- and long-run. These advantages of DSARDL over the conventional ARDL technique have made the DSARDL approach more popular in recent research examining factors affecting environmental degradation (see also, e.g. Refs.54–59). It is possible to describe the DSARDL model of this study in error correction form as follows:
Empirical results
Unit root tests
The DSARDL approach needs to be performed with data stationarity at level, that is I(0), and/or first difference, that is I(1). 60 In this context, the empirical analysis of this study starts with assessing the stationary level of the study variables by utilizing two well-known unit root tests, namely augmented Dickey–Fuller (ADF) and Phillips–Perron (PP). The results of ADF and PP unit root tests are reported in Table 3. Both unit root test results divulge that the null hypothesis of the relevant variables is not stationary at level and cannot be rejected at I(0) for all study variables except NR. However, it can be rejected after getting the first difference of the variables. These findings imply that the stationary level of NR is I(0) and that of I(1) for other study variables.
Unit root tests results.
Cointegration test
Since the DSARDL approach also requires cointegration among the to-be-used variables, 63 we second employ the Pesaran, Shin, and Smith (PSS 44 ) bounds test with the lag length of 1, which is estimated as optimal by all criterions (see Table 4), to analyze the cointegration association among the study variables. The obtained F-statistic from the PSS bounds test and the Kripfganz–Schneider (KS 64 ) critical values for each confidence level are provided in Table 5. 1 Results show that the obtained F-statistic (5.653) is higher than the upper KS critical value (4.42) of the 5% confidence level, thereby rejecting the null hypothesis that there is no cointegration among the relevant variables at the 5% confidence level. These findings reveal a long-run statistically significant cointegration among the study variables.
Lag order selection process results.
Indicates the lag order selected by the criterion.
PSS bounds test results.
Notes: This table represents the obtained F-statistics for the estimated model of the study from the PSS bounds test and upper and lower KS critical values. Lower and upper bounds’ critical values are denoted as
Dynamic simulated autoregressive distributive lag (DSARDL) test
Following the verification of stationarity and cointegration requirements, we run the DSARDL approach to investigate energy intensity, renewable energy, natural resources, economic growth, and carbon intensity nexus in India. The short and long-run coefficients obtained from the DSARDL approach are provided in Table 6. At first glance, it can be seen that the ECT is statistically significant and negative, implying the presence of a long-term association. The ECT displays an imbalance in the short run; however, in the long run, it is adjusting at a pace of 72.9% each year towards equilibrium.
Short and long run coefficients of DSARDL.
Notes: This Table represents short- and long-run coefficients of the DSARDL approach. 1%, 5%, and 10% confidence levels are represented as a, b, and c, respectively while Δ represent the short run estimate.
Looking at the coefficient estimation results, it is revealed that there is a positive connection between energy intensity and carbon intensity. Specifically, a 1% increase in energy intensity in India surges carbon intensity by 0.9% in the short run and 0.8% in the long run. These results divulge that energy intensity which implies inefficient usage of primary energy sources raises India's environmental pollution both in the short and long run. This result would suggest that India must find alternative ways to reduce energy intensity by improving energy efficiency to mitigate the increase in carbon intensity. While our result is supported by the findings of Rahman et al., 66 the result of Zhang et al. 67 who found that energy intensity reduces carbon emissions in China, contradicts our result. Additionally, a positive interrelationship is also seen between economic growth and carbon intensity in the long run. A 1% increase in economic growth in India boosts carbon intensity by 0.08% in the long run, implying that economic growth in India has an adverse long-run impact on environmental quality. This implies that economic prosperity comes at the cost of increased carbon intensity, except measures are taken to mitigate this effect. This result confirms prior research especially for the developing economies such as revealed in these studies; Kalmaz and Kirikkaleli, 68 Adebayo et al., 69 Khan et al., 70 Ji et al., 71 and Pata et al. 72
On the other hand, we determine a negative association between natural resources and carbon intensity. A 1% rise in natural resources decreases carbon intensity by 0.02% in the short run and 0.01% in the long run, implicating that India's environmental quality is positively affected by natural resources both in the short- and long-run. More natural resources helping to reduce carbon intensity might explain India’s negative coefficient of NRR. Rich natural resources reduce the demand for imported fossil fuel energy; additionally, these outcomes are related to independent energy sources, which have lower emissions than fossil sources22,73. Similarly, Tufail et al. 74 reported that natural resources rent contributes to improving the atmosphere by lowering CO2 emissions and this claim was further amplified by Damerau et al. 29 that natural resources usage in form of cropland and water mitigates GHG emissions across geographical regions in India.
Moreover, and as expected, a negative connection between renewable energy and carbon intensity is also seen, but this connection is not statistically significant in the long run. 0.08% of the short-run mitigation in carbon intensity is caused by a 1% surge in renewable energy. The result demonstrates that renewable energy is an effective strategy against environmental degradation. Based on the negative association between renewable energy and carbon intensity, it further establishes the environmental desirability of renewable energy as well attesting to the richness of renewable energy sources in India. Moreover, Ji et al. 71 and Alola et al. 4 found similar results in their studies on the impact of renewable energy on environmental degradation.
Impulse response evidence
Following an examination of the DSARDL short- and long-run coefficients, we investigate the 15-year effect of 1% positive and negative shocks applied to any predictor at the fifth year on carbon intensity while holding the other predictors constant and the impulse response plots are displayed in Figures 1–4 automatically provided by the DSARDLS method. It is determined from the plots in Figure 1 that a positive (negative) shock to energy intensity raises (declines) carbon intensity both in the short and long run. This indicates that India's policymakers need to take strict measures to reduce energy intensity to prevent environmental degradation. In addition, we see from the plots in Figure 2 that a positive (negative) shock to renewable energy curbs (surges) carbon intensity both in the short and long run. This implies that more eco-friendly energy sources, such as renewable energy, are needed to promote India's environmental quality.

Impulse response plots of energy intensity. Notes: The plot on the left shows the 15-year response of carbon intensity to a 1% positive shock applied to energy intensity at the fifth year, while the plot on the right shows the 15-year response of carbon intensity to a 1% negative shock applied to energy intensity at the fifth year. The average predicted carbon intensity values are depicted as ●. Also, 75%, 90%, and 95% confidence levels are represented in a vertical line using the colors from gray to light blue, respectively.

Impulse response plots of renewable energy. Notes: The plot on the left shows the 15-year response of carbon intensity to a 1% positive shock applied to renewable energy at the fifth year, while the plot on the right shows the 15-year response of carbon intensity to a 1% negative shock applied to renewable energy at the fifth year. The average predicted carbon intensity values are depicted as ●. Also, 75%, 90%, and 95% confidence levels are represented in a vertical line using the colors from gray to light blue, respectively.

Impulse response plots of natural resources. Notes: The plot on the left shows the 15-year response of carbon intensity to a 1% positive shock applied to natural resources at the fifth year, while the plot on the right shows the 15-year response of carbon intensity to a 1% negative shock applied to natural resources at the fifth year. The average predicted carbon intensity values are depicted as ●. Also, 75%, 90%, and 95% confidence levels are represented in a vertical line using the colors from gray to light blue, respectively.

Impulse response plots of economic growth. Notes: The plot on the left shows the 15-year response of carbon intensity to a 1% positive shock applied to economic growth at the fifth year, while the plot on the right shows the 15-year response of carbon intensity to a 1% negative shock applied to economic growth at the fifth year. The average predicted carbon intensity values are depicted as ●. Also, 75%, 90%, and 95% confidence levels are represented in a vertical line using the colors from gray to light blue, respectively.
Moreover, it is specified from the plots in Figure 3 that a positive (negative) shock to natural resources reduces (increases) carbon intensity both in the short and long run. These results imply that India should use its natural resources effectively to reduce its environmental pollution. Additionally, we understand from the plots in Figure 4 that a positive (negative) shock to economic growth escalates (mitigates) carbon intensity both in the short and long run. This result divulges that India's economic growth programs should be revised with more environmentally friendly policies.
Diagnostic tests
The outputs of the diagnostic tests to determine whether the findings of DSARDL mentioned above are valid or not are presented in Table 7. Results disclose that the model used has no issues with serial correlation, heteroscedasticity, ARCH effect, omitted variable, and non-normality. Also, plots of cumulative sum (CUSUM) and CUSUM of squares tests in Figure 5 reveal that the model is stable. The outputs of the diagnostic tests prove the validity of the obtained findings from the DSARDL method.

Plots obtained from the CUSUM and CUSUM of squares tests.
Model diagnostic tests results.
Kernel-based regularized least squares (KRLS) test
This section includes a robustness analysis to check the validity of the DSARDL findings. For robustness, we also analyze the impact of energy intensity, renewable energy, natural resources, and economic growth on carbon intensity in India by running the KRLS method established by Hainmueller and Hazlett. 75 This approach includes a machine learning algorithm with econometric properties. The KRLS method, in contrast to existing standard econometric methods, constructs pointwise derivatives and mean marginal effects, does hypothesis testing, and generates robust and reliable estimates. In a similar vein, the KRLS method outperforms existing machine learning algorithms with issues of misspecification bias over statistical judgments. In addition, the KRLS method provides flexible and interpretable parameters in regression and classification conundrums with undetermined functional forms. In other words, it determines the functional form of the data series under consideration and protects practitioners against misspecification bias. Furthermore, the KRLS approach is helpful for analysis incorporating learning about the data creation process, model-driven causal examination, forecast, and imputation of missing data.75,76
The results of the KRLS test are reported in Table 8. We determine that the KRLS model's predictive power (R2) is 0.989%, demonstrating that the factor variables (namely energy intensity, renewable energy, natural resources, and economic growth) explain 98.9% of the variance in the dependent variable, that is, carbon intensity. The estimated average marginal effects imply that energy intensity and economic growth increase carbon intensity by 0.832% and 0.084%, whereas renewable energy and natural resources decrease carbon intensity by 0.055% and 0.008%, respectively.
The KRLS results.
Notes: This table represents pointwise derivatives obtained from KRLS. 1%, 5%, and 10% confidence levels are represented as a, b, and c, respectively. Lambda is the value of regularization parameter, sigma specifies the bandwidth, and Looloss shows Leave-On-Out error. Also P25, P50, and P75 signify 1st quartile, median, and 3rd quartile, respectively.
Furthermore, we plot the pointwise marginal effects of energy intensity, renewable energy, natural resources, and economic growth on carbon intensity in Figures 6–9. Figures 6 and 9 illustrate that energy intensity and economic growth affect carbon intensity (environmental quality) positively (negatively) up to a certain point, after which the effect decreases slightly. After that, however, the positive (negative) effect on carbon intensity (environmental quality) continues.

Pointwise marginal effects of energy intensity on carbon intensity.

Pointwise marginal effects of renewable energy on carbon intensity.

Pointwise marginal effects of natural resources on carbon intensity.

Pointwise marginal effects of economic growth on carbon intensity.
On the other hand, Figures 7 and 8 show that renewable energy and natural resources negatively (positively) affect carbon intensity (environmental quality) up to a certain point, after which the effect decreases slightly. And then, the negative (positive) effect on carbon intensity (environmental quality) increases again. The findings obtained from the DSARDL method are robust by all these findings. Finally, the findings of the study are summarized in Figure 10.

Visualization of the results.
Conclusion and policy recommendations
For most nations, a sustainable and pollution-free environment has been a priority. However, it should be emphasized that it is practically hard to arouse a shared interest due to the glaring inequalities in economic and social standards between nations. Several industrialized and developing nations intend to reduce their carbon intensity to realistic targets set by national policies, both independently and as participants in multilateral agreements. 18 This investigated energy intensity, renewable energy, natural resources, economic growth, and carbon intensity nexus in India using data from 1970 to 2020. The DSARDL method, employed in this study, uses stochastic simulations to estimate, simulate, and provide the short and long-run coefficients with p-values for statistical inference of the model used, which typical approaches may neglect. According to the DSARDL, energy intensity and carbon intensity are positively correlated. Economic growth and carbon intensity are also observed to have a positive relationship over time. However, we find a negative relationship between natural resources and carbon intensity. Additionally, a negative association between renewable energy and carbon intensity is observed, but this correlation is not statistically significant over the long term. The KRLS approach was used for robustness since it beats current machine learning algorithms in terms of bias caused by misspecification over statistical judgments.
Although this study profers interesting result, yet it suggests that salient environmental-related questions could further answered in the context of India if regional-level and cross industry examination is considered. Moreover, this study could further be explored by employing more macroeconomic and energy indicators alongside the GHGs or ecological footprint components considering that India state is ecological deficit.
Overall, this study strengthens the idea that as the demand for energy increases and a country advances economically, it prompts a surge in carbon intensity. Therefore, it is pertinent for India to introduce policies to help abate the environmental impact of this demand as they prosper economically. For instance, investment in clean energy development could be advanced through provision of credit and energy subsidies as a way of depleting share of fossil fuels in total primary energy. Additionally, to further discourage carbon emission activities, environmental tax especially for energy-intensive industries could be increased while providing subsidies that encourage energy transition. In addition to the important implications of this study, India has to keep a close watch on how household energy use varies by income level in light of the country’s significant population of people that lack access to clean cooking. Additionally, even though increasing natural resources rent reduces carbon intensity, the government should take urgent steps to ensure sustainable natural resources along with rising revenue. This is because increased natural resource extraction through agriculture, deforestation, and mining which are prominent economic activities in the country could continue to hurt it environment if not effectively checked.
Footnotes
Abbreviations
Author note
Andrew Adewale Alola is currently affiliated with Faculty of Economics, Administrative and Social Sciences, Nisantasi University, Istanbul, Turkey.
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
