1 Introduction

Human history is riddled with infectious diseases that left thousands of people dead from the start. Global pandemics have always caused issues to human lives and the global environment. Researchers believe that pandemics' probabilities have been dramatically amplified over the past century, mainly due to greater exploitation of the natural environment, urbanization, improved global travel and integration, and land-use changes (Jones, Patel & Levy, 2008). These trends will keep on getting intensify with each passing day. In the early 2000s, when SARS (severe acute respiratory syndrome) was started spreading worldwide, the World Health Organization (WHO) issued a warning of massive damage to the whole world (Madhav & Wolfe, 2017).

No doubt, the world has suffered pandemics throughout history, and they keep on appearing with higher frequency each time. The world is currently dealing with another great pandemic of the century that is COVID 19. Coronavirus is a highly contagious virus in Wuhan, China, in December 2019 (Haleem & Javaid, 2020). From China, it has spread all across the world. Millions of people have lost the battle of their lives against this deadly virus. COVID-19 has rapidly affected our daily life and business activities. It has caused restricted movements and disrupted supply chain activities worldwide. Most businesses face lower or no sales, reducing manufacturing activities. COVID-19 has also greatly influenced the global environment. Figure 1 shows biweekly CPVID-19 cases per million people.

Fig. 1
figure 1

Biweekly COVID-19 daily-confirmed cases

Since the last few decades, the world has experienced a tremendous environmental pollution increase. Despite the enormous efforts, environmental pollution is a significant issue and poses a continuing risk to human lives. In the current global environment, where business activities are increasing, environmental pollution is also going upward at a higher pace. Most researchers believe that environmental pollution is a source of various diseases and affects human health by increasing exposure to different diseases (Kelishadi, 2012). Undoubtedly, this problem is most significant in developing countries. The traditional sources of pollution are still in processes, such as contaminated water supply, inadequate waste management, poor sanitation, and high exposure to indoor air pollution (Kumar, Taxak, et al., 2021). Pollution is directly correlated with the population. According to the United Nations (2019), the world population may rise to 11 billion by 2100, which will cause the environmental population.

Environmental pollution has reached to troublesome amount all over the world. Economic development, industrialization, and urbanization have caused an increase in greenhouse emissions (Kumar & Sharma, 2017; Kumar, Thanki, et al., 2021; Sharma et al., 2020; Tobías et al., 2020; Wang, Chen, et al., 2020; Wang, Zhang, et al., 2020). Thus, carbon dioxide and other toxic gases are released into the atmosphere polluting the air and increasing its temperature. This has become the reason for the depletion of the ozone layer, which significantly affects the daily lives of the people (Edwards, 2018).

Undoubtedly, the novel coronavirus has caused worldwide disruptions that affected the environment and the climate. A significant decrease in travel activities due to lockdowns has dropped air pollution considerably. According to recent research, China has experienced a 25% and 50% decrease in carbon emissions and nitrogen oxides, respectively, since the pandemic (Myllyvirta, 2020). A scientist estimated that it might save at least 77,000 lives over two months. Almost every country has implemented lockdowns to reduce coronavirus spread, which reduced industrial activities and transportation services. This has dramatically brought down the worldwide pollution level. The European Space Agency observed a significant decline in nitrous emissions from factories, power plants, and cars in Northern Italy, mainly due to strict curfew in the region (Green, 2020).

The European Space Agency (ESA) and the National Aeronautics and Space Administration (NASA) have been closely monitoring the ozone layer in the current pandemic. According to them, the decrease in nitrogen oxide started in Wuhan, China, and then slowly spread throughout the world drastically (Zhang & Zhang, 2020). No doubt, COVID-19 has caused havoc to the whole world. However, it is also considered a “Blessing in Disguise,” as nature is reclaiming itself. Poor quality is harmful to the entire population. Around 4.6 million people die every year due to poor quality worldwide. Air pollution is a huge problem, and its effects are visible to the developed nations. In 2012, around 193,000 people died in Europe due to air pollution (Brauer & Burnett, 2017).

1.1 COVID-19 in India

Like any other country, India has also experienced a sudden decrease in carbon emission and environmental pollution due to lockdown implementation since March 2020 (Gautam et al., 2020). Delhi, India's capital is considered one of the most polluted megacities of the whole world. A report was issued by the World Health Organization in 2018, in which Delhi ranked high in the list of most polluted cities of the world (United Nations, 2018). Mumbai and Delhi, the megacities, are experiencing PM 2.5 concentrations for the last several years, far above tolerable limits, defined by NAAQS (National Ambient Air Quality Standards). The increased air pollution has become the source of various health issues, including acute asthma, pneumonia, respiratory disorders, etc. (Chakraborty & Maity, 2020; Rizou et al., 2020; SanJuan-Reyes et al., 2020). In 2017, owing to increased air pollution, the Indian Council of Medical Research (ICMR) has declared a community health emergency in Delhi (Verma & Prakash, 2020; Lokhandwala & Gautam, 2020; Wang & Su, 2020; Amdade et al., 2021; Gautam et al., 2020, 2021). India's bad air quality is a key reason for premature deaths (Dandona & Shukla, 2017).

To minimize coronavirus spread, the Indian government implemented a nationwide lockdown from March 24 for three weeks, and later, it further extended. The government has also banned all kinds of industrial activities and mass transportations all over the country. Resultantly, Sharma et al. (2020) stated pollution level had been significantly decreased in 88 cities of India, only after four days of commencing the lockdown event. The nitrogen oxide level has also dropped to 90 µmol/m2 in March 2020 in Delhi, 158 µmol/m2 in March 2019. Mumbai further dropped to 77 µmol/m2, 122 µmol/m2 last year (Poetzscher, 2020). The coronavirus has caused disruptions in India, but it has proved a “blessing in disguise” for minimizing the country’s environment and air pollution (Eroğlu, 2020; Shakil et al., 2020; Wilder-Smith & Freedman, 2020). Figure 2 represents daily COVID-19-confirmed cases and deaths in India from January 30, 2020, to February 21, 2021.

Fig. 2
figure 2

Daily confirmed cases and deaths in India

The current study aims to analyze the effect of carbon emission on sustainable development during COVID-19. As per the rules and regulations of lockdown, all types of human activities were banned, and India's people were restricted from living in their homes. The daily pausing activities reduce transportation emissions, and the businesses do not contribute to pollution and emission of toxic gases during the lockdown. Purposely, this research analyses the relationship between lockdown and carbon dioxide emission and the nature and direction of relation to its long-run relationship or short-run relationship.

2 Literature review

The coronavirus outbreak restricted almost all avoidable activities in China and worldwide since the first call of lockdown in Wuhan on January 23, 2020. The transmission rate of SARS-COVID-2 is comparatively high compared to its predecessors (Bherwani, Anjum, et al., 2020). By employing the Community Multi-Scale Air Quality (CMAQ) model with monthly anthropogenic emission data of China, it is observed that the concentration of PM2.5 decreased almost 20% in Guangzhou, Shanghai, Beijing, Shijiazhuang, Tianjin, Jinan, Taiyuan, Xi’an, Zhengzhou, and Wuhan. It is observed that the meteorology was unfavorable during the simulation episode. The PM2.5 emission’s concentration is reduced ~ 2 times in lockdown. A large emission reduction is observed in transportation and slight emission reduction in the industrial sector due to working being stopped in the industrial sector. Still, it does not decrease the pollution of China prominently. The pandemic and lockdown are imposed on all countries for a specific time and specific reason, so it is not applicable for reducing pollution. For avoiding severe air pollution, more efforts should be made to complete the task of environmental sustainability (Ambade et al., 2021; Bherwani, Nair, et al., 2020; Srivastava et al., 2020; Wang, Chen, et al., 2020; Wang, Zhang, et al., 2020). To control the rapidly increasing affected cases and death rate, the major strategies in clinical trials, medications, vaccinations, and social distancing for social and individual protection were designed by individual countries and the World Health Organization (WHO) (Gautam & Hens, 2020).

The morbidity and mortality are being increased due to air pollution, and the rate of it is rapidly increased year by year. Air pollution leads to health risks and environmental degradation by acid rains and reduced visibility. The rise in obstructive pulmonary disease (COPD), cardiovascular disease, and respiratory ailments like asthma and chronic bronchitis are just because of poor air quality management. As we know that the COVID-19 outbreak abruptly stopped all human activities, so in these situations and circumstances, the air pollution showed a significant reduction in it. If the rate of newly affected cases of COVID-19 was increased, the causes of morbidity and mortality with air pollution were decreased. The National Aeronautics and Space Administration (NASA) and European Space Agency (ESA) released the data of healing environmental sustainability during strict lockdown during an outbreak of COVID-19. The reduction rate of pollution in some of the epicenters like Wuhan, Italy, Spain, and the USA has 30%, which increases the sustainability of the environment and quality of air (Muhammad et al., 2020). Bherwani, Anjum, et al. (2020) monetarily quantified the declining trend in health impact due to air pollution reduction during the lockdown. Utilizing two types of data, i.e., environmental-related data and monetized health impact data from 2019 to 2020 and capturing four cities—Delhi, London, Paris, and Wuhan, analyzed air pollution reduction. Undoubtedly, these four cities suffered from economic damages, but this research noted the comparative preventing air pollution reduction benefits due to lockdown. This research also provides a pathway to compare the magnitude of actual damage environmentally and the damage-related lockdown.

Gautam et al. (2020) used seven environmental pollutants, i.e., PM2.5, PM10, NO2, NH3, SO2, CO2, and O3, to calculate the Air Quality Index (AQI) from February 17, 2020, to May 4, 2020. In this research, Gautam et al. (2021) measured the variations of air quality of Delhi (DTU, Okhla, and Patparganj), Haryana (Jind, Palwal, and Hisar), and Uttar Pradesh (Agra, Kanpur, and Greater Noida) and data have been retrieved from the Central Pollution Control Board website. The Air Quality Index has upgraded 30–47% approximately during the lockdown. The air mass contribution and healing layers of O3 have been observed by the air mass back trajectory (AMBT) analysis in India and Europe, Africa, and Gulf countries during the lockdown. Gautam et al. (2021) revised the research by including extra air quality factors and upgraded the research direction's accuracy. This research captured the eight most polluted cities of India, i.e., Mumbai, Delhi, Bangalore, Hyderabad, Lucknow, Chandigarh, Kolkata, and Ahmedabad. The environmental sustainability behavior of these cities is checked during the COVID-19 lockdown. The significant reduction of environmental pollutants over Major Dhyan Chand Stadium has been observed as PM2.5 (63%), PM10 (56%), NO2 (50%), SO2 (9%), and CO (59%). In comparison, the PM2.5 (44%), PM10 (50%) decreased over Chhatrapati Shivaji International Airport in India's initial few weeks of strict lockdown.

The overall Air Quality Index (AQI) of Chandigarh and Ballygunge was dramatically upgraded from 7 to 67%. Divergent places of India's major cities were observed, and individual cities of crowded places showed the same results. Long-range transportation has occupied a large portion of toxic emissions. During the lockdown, the trade activities between Iran, Afghanistan, Saudi Arabia, Indian Gangetic plains, and Delhi Non-capital region were stopped. The results of the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) analysis also confirmed the up-gradation of environmental quality during the lockdown. Bherwani, Nair, et al. (2020) analyzed India's data, especially the states of Maharashtra, Gujarat, Madhya Pradesh, Rajasthan, Tamil Nadu, West Bengal, Uttar Pradesh, and Delhi as a union territory on lockdown and social distancing policies by employing Bayesian probability model (BPM). The accuracy and positive resulting of lockdown were analyzed using the change point analysis (CPA), which confirmed the diminishing trend in newly reported cases and the death rate of COVID-19 patients. The delta values are between 9.6 and 4.7 for Gujarat and Madhya Pradesh and Maharashtra, respectively, At the same time, all states showed a strong significant correlation \((p<0.05,t-stat>t-critical)\) toward a different period of CPA and lockdown (cases per population (CPP) and cases per unit area (CPUA)). The weak correlation \((p<0.1, t-stat<t-critical)\) is observed by delta and cases per unit population density (CPD). Although the correlation is significant for CPP and CPUA, Pearson’s correlation confirmed a negative direction.

3 Methodology

3.1 Data description

The daily data of deaths, covid-19-confirmed cases, and lockdown are utilized to analyze the impact of COVID-19 on CO2 emission. This research comprises data from January 31, 2020, to December 1, 2020. The data series is divided into three parts: pre-lockdown, during the lockdown, and post-lockdown (Mohsin, Naiwen, et al., 2020; Mohsin, Naseem, et al., 2020). In India, the first confirmed COVID-19 case was reported on January 30, 2020. Under India's pandemic government, the global and domestic adverse spreading situation announced a strict lockdown on March 25, 2020, and it was prolonged till May 31, 2020. The visual presentation of CO2 emission data from January 30, 2020, to December 1, 2020, can be seen in Fig. 1. The daily data of confirmed cases and deaths are extracted from the World Health Organization (WHO). In contrast, the daily data of CO2 emission in India are obtained from carbon monitor (Asadi et al., 2020; Chakraborty & Ghosh, 2020; Guo et al., 2020).

The analysis is performed on the below-mentioned model:

$$\ln {\text{CO}}_{2t} = f\left( {{\text{NC}}_{t}^{ + } , {\text{NC}}_{t}^{ - } , {\text{ND}}_{t}^{ + } , {\text{ND}}_{t}^{ - } , {\text{LD}}_{t}^{ + } , {\text{LD}}_{t}^{ - } } \right)$$
(1)

The partial sum of positive and negative shocks in NC, ND, and LD is determined by employing the nonlinear autoregressive distributed lag (NARDL) method. Before exploring the primary model, the divergent test's data series have been checked for suitability of the NARDL method. The unsuitable data set can be misled due to fundamental analysis (Ouattara, 2004).

3.2 Unit root test

Unit root test is used to check the stationarity of data series. The augmented Dickey–Fuller (ADF) and Phillips–Perron (PP) tests are employed to inspect the unit root.

$$\Delta Y_{t} = \alpha + \beta_{t} + \rho Y_{t - 1} + \mathop \sum \limits_{i = 1}^{k} \gamma_{i} \Delta Y_{t} + e_{t}$$
(2)
$$\Delta Y_{t} = \left( {p - 1} \right)Y_{t - 1} + e_{t}$$
(3)

The equation-2 is ADF, and equation-2 PP test equation for unit root to check the data series's stationarity. These equations \(\Delta\) represent the first difference operator, and \({Y}_{t}\) shows a high order of autocorrelation concerning time. The endogenous variables are presented as \({Y}_{t-1}\). The \({e}_{t}\) is the error term in both tests of equation-2 and equation-3 (Fig. 3).

Fig. 3
figure 3

CO2 emission in India during pandemic (COVID-19)

3.3 Standard nonlinear autoregressive distributed lag (NARDL) model

The nonlinear autoregressive distributed lag (NARDL) model is introduced by Shin et al. (2014). The standard ARDL model is restricted to check the linear relationship between independent and dependent variables, while some are nonlinearly related (Hu et al., 2018; Shahzad et al., 2017). The nonlinear ARDL model is covered with negative and positive changes of variables. The \({{\varvec{x}}}_{{\varvec{t}}}\) in the following equation of the NARDL model shows the decomposition into a partial sum of positive and negative changes as:

$$x_{t} = x_{0} + x_{t}^{ + } + x_{t}^{ - }$$
(4)
$${\text{where}} x_{t}^{ + } = \mathop \sum \limits_{i = 1}^{t} \Delta x_{i}^{ + } = \mathop \sum \limits_{i = 1}^{t} \max (\Delta x_{i} ,0)$$
$$x_{t}^{ - } = \mathop \sum \limits_{i = 1}^{t} \Delta x_{i}^{ - } = \mathop \sum \limits_{i = 1}^{t} \min (\Delta x_{i} ,0)$$

The following equation can be used to examine the asymmetric long-run equilibrium relationship:

$$y_{t} = \beta^{ + } x_{t}^{ + } + \beta^{ - } x_{t}^{ - } + \mu_{t}$$
(5)

In equation-3, the parameter \({{\varvec{\beta}}}^{+}\) and \({{\varvec{\beta}}}^{-}\) are asymmetric long-run parameters which accordingly related to positive and negative variations in \({{\varvec{x}}}_{{\varvec{t}}}\). The combination of equation-4 and equation-5 with the ARDL model generates the following NARDL model:

$$\Delta y_{t} = \beta_{0} + \beta_{1i} \mathop \sum \limits_{i = 1}^{t} \Delta y_{t - i} + \beta_{n,i} \mathop \sum \limits_{i = 1}^{t} ( + \Delta x_{t - i}^{ + } + \ldots \Delta x_{t - i}^{ - } ) + \vartheta_{1} y_{t - i} + \vartheta_{n,i} x_{t - 1}^{ + } + \vartheta_{n,i} x_{t - 1}^{ - } + \mu_{t}$$
(6)

The equation-6 consists of two parts: the short-run asymmetry is tested by \(\beta _{{n,i}} \mathop \sum \limits_{{i = 1}}^{t} \Delta x_{{t - i}}^{ + } ,~\beta _{{n,i}} \mathop \sum \limits_{{i = 1}}^{t} \Delta x_{{t - i}}^{ - }\) while the long-run asymmetry determined by \({\beta }^{+}={\beta }^{-}\) (i.e., \(-\vartheta/\beta^{+}= -\vartheta/\beta^{-}).\)

The nonlinear ARDL model is formulated for this study as:

$$\begin{gathered} \Delta \ln {\text{CO}}_{2t} = \beta_{0} + \beta_{1i} \mathop \sum \limits_{i = 1}^{t} \Delta \ln {\text{CO}}_{2t - i} + \beta_{2a} \mathop \sum \limits_{i = 1}^{t} \Delta \ln {\text{NC}}_{t - i}^{ + } + \beta_{2b} \mathop \sum \limits_{i = 1}^{t} \Delta \ln {\text{NC}}_{t - i}^{ - } + \beta_{3a} \mathop \sum \limits_{i = 1}^{t} \Delta \ln {\text{ND}}_{t - i}^{ + } \hfill \\ \quad \quad \quad \quad \quad \quad + \beta_{3b} \mathop \sum \limits_{i = 1}^{t} \Delta \ln {\text{ND}}_{t - i}^{ - } + \beta_{4a} \mathop \sum \limits_{i = 1}^{t} \Delta \ln {\text{LD}}_{t - i}^{ + } + \beta_{4b} \mathop \sum \limits_{i = 0}^{t} \Delta \ln {\text{LD}}_{t - i}^{ - } + \vartheta_{1} \ln {\text{CO}}_{2t - i} + \vartheta_{2a} \Delta \ln {\text{NC}}_{t - 1}^{ + } \hfill \\ \quad \quad \quad \quad \quad \quad + \vartheta_{2b} \Delta \ln {\text{NC}}_{t - 1}^{ - } + \vartheta_{3a} \Delta \ln {\text{ND}}_{t - 1}^{ + } + \vartheta_{3b} \Delta \ln {\text{ND}}_{t - 1}^{ - } + \vartheta_{4a} \Delta \ln {\text{LD}}_{t - 1}^{ + } + \vartheta_{4b} \Delta \ln {\text{LD}}_{t - 1}^{ - } \hfill \\ \quad \quad \quad \quad \quad \quad + \mu_{t} \hfill \\ \end{gathered}$$
(7)

where \({CO}_{2t}\) represents the CO2 emission of India in the COVID-19 era at time \(t\), the partial sum of positive and negative changes in new confirmed cases of COVID-19 represented by \({NC}_{t-i}^{+}\) and \({NC}_{t-i}^{-},\) respectively, at time \(t\). The newly confirmed death is shown as \({ND}_{t-i}^{+}\) and \({ND}_{t-i}^{-},\) and the positive partial sum and negative partial sum of lockdown are presented as \({LD}_{t-i}^{+}\) and \({LD}_{t-i}^{-}\) with time \(t\). The lockdown is used as a dummy variable in this analysis (Barati & Fariditavana, 2020; Bahmani‐Oskooee & Saha, 2017). The dummy variable is justified in assigning the binary number; 1 is assigned to the lockdown period and 0 for the unlocked period. The error term is presented as \({\mu }_{t}\). The NARDL model can only be employed in variables integration at the order I(0) or I(1). The NARDL model is not suitable if the variables are integrated into two I(2) (Greenwood-Nimmo & Shin, 2013; Shin et al., 2014). Using the ADF and PP test, the PP test applies a nonparametric correction to the standard ADF test statistics. More general dependence in error and conditional heteroskedasticity is also checked. Both tests are used to check the stationarity of the series.

4 Results and discussion

The summary statistics of CO2 emission, newly confirmed cases, and newly confirmed deaths are presented in Table 1. The mean of CO2 emission (6.127) is lower than NC (30,924.210) and ND (449.745). The maximum and minimum values of variables present a clear picture of the relationship between CO2 emission and COVID-19. The maximum value of the CO2 emission is 8.580, and the minimum value is 3.630. The 4.950 metric tones reduction in CO2 emission is observed in India during COVID-19, while an increasing trend is observed in new confirmed cases and newly confirmed deaths. The drastic increase in the standard deviation of NC (30,921.330) and ND (405.705) is shown. The value of standard deviation variation of CO2 emission is 1.057. The latest confirmed cases and newly confirmed deaths show positive skewness, while CO2 emission is negatively skewed. The kurtosis value of all variables is significant at a 1% level of significance, and more than three, which means data are not normally distributed (Mohsin, Naiwen, et al., 2020; Mohsin, Naseem, et al., 2020; González et al., 2020; Shehzad et al., 2020; Hashmi et al., 2020).

Table 1 Descriptive statistics

The data series's stationarity is checked by performing two unit root tests, i.e., the augmented Dickey–Fuller (ADF) and Phillips–Perron (PP) test. The unit root test results are presented in Table 2, which shows that all variables are stationary at first difference. There are two terms to show cointegration I(0) for non-stationary series and I(1) for stationary data series. The data series's stationarity confirms the data fit for the cointegration process in the NARDL framework.

Table 2 Unit root test

The results of the bound test for cointegration are displayed in Table 3. The F-statistics value (3.993) is larger than the upper critical value (3.99) at a 1% level of significance. The automatic lag specification of the NARDL model is 1, 1, 1, 1, 0, 0, 1, which is selected based on the Akaike information criterion (AIC) (Mohsin, Naiwen, et al., 2020; Mohsin, Naseem, et al., 2020; Zhang et al., 2017; Shah et al., 2020). The model specification with AIC-based lag length is presented in Fig. 4. It vehemently rejects the null hypothesis of no cointegration (F-statistics < lower bound critical value). F-statistics significance at a 1% level confirms the cointegration occurrence for a long-run relationship between CO2 emission and COVID-19 with its determinants (Katrakilidis & Trachanas, 2012; Sarfraz et al., 2018, 2020; Zhang et al., 2017).

Table 3 The result of bound tests for cointegration test
Fig. 4
figure 4

AIC top 20 models

After confirming the long-run cointegration relationship between CO2 emission and COVID-19 (new confirmed cases, new confirmed deaths, and lockdown), the next step particularly determined the long-run or short-run asymmetry from the primary model (7). The long-run dynamic asymmetry estimation results are presented in Table 4. The long-run coefficient of CO2 emission (0.87225) is positively significant at a 1% level of significance. As per the results, we can say that a 1% increase in COVID-19 with its determinants becomes the reason for an 87.225% decrease in CO2 emission. The NC+ (−1), NC, LD+, LD, and LD (−1) are significant at 1% level, while NC+ and ND are significant at 10% and 5% level of significance. The NC (−1), ND+, and ND+ (−1) demonstrate the positivity with insignificance.

Table 4 Dynamic asymmetric estimation of CO2 emission

The asymmetrically short- and long-run results are displayed in Tables 5 and 6. The recent confirmed case is significant at 5% with a positive sign and the negative new confirmed cases also significant at a 1% level. The new confirmed death positive is displayed in a table, which is positively insignificant. The lockdown coefficient shows negative significance at a 1% level for the short run. The ECT term confirms the long-run relationship between CO2 emission and COVID-19 with its determinants with its significance at 1% and the negative sign's assumption. The second half of Table 5 presents long-run asymmetric coefficients with their confirmation of significance or insignificance. The NC+ and NC are significant at a 5% level of significance, while ND + and ND are significant at a 10% level of significance. The LD+ is significant at the 1% level, and LD is insignificant. The significance of LD+ at 1% with value −2.5040 is explained by the nonlinear relationship between CO2 emission and lockdown. The CO2 emission decreases due to strict lockdown in India during pandemic spread (Azam et al., 2020; Mohsin, Naiwen, et al., 2020; Mohsin, Naseem, et al., 2020; Omoke et al., 2020). As per the graphical presentation, the CO2 emission started decreasing from the starting point of COVID-19 breakout because people stayed in their homes to avoid the risk of viruses. A rapid decrease in CO2 emission is observed after strict lockdown. All types of business and the flow of transport were stopped in strict lockdown, which becomes the reason for CO2 emission reduction in India. This research has covered three phases: pre-lockdown, during-lockdown, and post-lockdown. The adjusted assessment of asymmetry in the presented long-run equilibrium with positive and negative shocks is displayed in a dynamic multiplier graph (Figs. 5, 6, 7). The negative and positive CO2 emission shocks generate a linear mixture of dynamic multiplier asymmetric curves.

Table 5 Result of asymmetric short run
Table 6 Result of asymmetric long-run coefficients
Fig. 5
figure 5

CO2 and NC dynamic multiplier

Fig. 6
figure 6

CO2 and ND dynamic multiplier

Fig. 7
figure 7

CO2 and LD dynamic multiplier

Table 7 contains the results of diagnostic tests. The value of Breusch–Godfrey is 1.301, which confirms the model is free from serial correlation. The insignificance of Jarque–Bera with 1.524, ARCH test with 1.734, and Ramsey RESET test with F-statistics value of 0.511 confirmed the normality, absence of heteroskedasticity, and correct specification of the model, respectively. The CUSUM and CUSUM square test shows the stability of the model. The CUSUM and CUSUM2 are visually presented in Figs. 8, 9, which shows the cumulative sum of recursive residuals and its square lies between critical bound 5% level of significance. The employed model is stable and trustworthy to determined short- and long-run relationship coefficients (Van & Bao, 2019).

Table 7 Diagnostic test
Fig. 8
figure 8

CUSUM for CO2 emission

Fig. 9
figure 9

CUSUM2 for CO2 emission

5 Conclusion

Every field of life has been adversely shaken by spreading pandemic COVID-19, but it left some positive effect on environmental sustainability. The holes of ozone layers were healing in this pandemic era because of lockdowns worldwide. People were stuck in their homes and all types of businesses and transportations were closed, which reduced the emission of toxic gases in the air. This research focused on the relationship between CO2 emission and COVID-19 with its determinants. The asymmetric relationship among CO2 emission, newly confirmed cases, newly confirmed deaths, and lockdown in the COVID-19 era is observed by utilizing daily data from January 30, 2020, to December 1, 2020, for India. The nonlinear autoregressive distributed lag (NARDL) or asymmetric ARDL model is employed to check the long- and short-run relationship between variables. This model's constructive contribution is to simultaneously capture both short- and long-run relationships of variables with negative and positive nature of the relationship. The results of this research confirm the long-run relationship between CO2 emission and COVID-19 in India.

The CO2 emission is negatively influenced by strict lockdown in India. In India, the pandemic starts spreading on January 30, 2020. At the start of COVID-19 spreading, people individually took care to avoid the risk of being infected. Later on, the rapid increase in confirmed cases of COVID-19 patients and increasing death toll compel India's government for strict lockdown. On March 25, 2020, government of India announced a strict lockdown, which extended until May 31, 2020. From January 30, 2020 to March 25, 2020, a minor decrease in CO2 emission has been observed, but a sudden downward trend with a high rate observed after one day of strict lockdown in India. The unlocked situation again contributed to India's CO2 emission, but it is still below the starting value (30 January 2020) of CO2 emission. This research suggested that government, policymakers, and businesses reduce the toxic emission of gases for a sustainable environment. No doubt, the lockdown became the cause of CO2 emission reduction in India’s environment, but a short-period strategy under special circumstances. This lockdown declared one thing that the human activities and industrial activities in India are unhealthy and environmental unfriendly. The government of India, policymakers, and environmentalists should make the people know about environmentally friendly activities and restricted industries in the sense of toxic gas emissions to enhance environmental purity and prolong its sustainability.