Strong link between coronavirus count and bad air: a case study of India.

The present study aims to highlight the contrast relationship between COVID-19 (Coronavirus Disease-2019) infections and air pollutants for the Indian region. The COVID-19 data (cumulative, confirmed cases and deaths), air pollutants (PM10, PM2.5, NO2 and SO2) and meteorological data (temperature and relative humidity) were collected from January 2020 to August 2020 for all 28 states and the union territory of India during the pandemic. Now, to understand the relationship between air pollutant concentration, meteorological factor, and COVID-19 cases, the nonparametric Spearman's and Kendall's rank correlation were used. The COVID-19 shows a favourable temperature (0.55-0.79) and humidity (0.14-0.52) over the Indian region. The PM2.5 and PM10 gave a strong and negative correlation with COVID-19 cases in the range of 0.64-0.98. Similarly, the NO2 shows a strong and negative correlation in the range of 0.64-0.98. Before the lockdown, the concentration of pollution parameters is high due to the shallow boundary layer height. But after lockdown, the overall reduction was reported up to 33.67% in air quality index (AQI). The background metrological parameters showed a crucial role in the variation of pollutant parameters (SO2, NO2, PM10 and PM2.5) and the COVID-19 infection with the economic aspects. The European Centre for Medium-Range Weather Forecasts derived monthly average wind speed was also plotted. It can see that January and February of 2020 show the least variation of air mass in the range of 1-2 m/s. The highest wind speed was reported during July and August 2020. India's western and southern parts experienced an air mass in the range of 4-8 m/s. The precipitation/wet deposition of atmospheric aerosols further improves the AQI over India. According to a study, the impact of relative humidity among all other metrological parameters is positively correlated with Cases and death. Outcomes of the proposed work had the aim of supporting national and state governance for healthcare policymakers.

Introduction

In late December 2019, many cases related to pneumonia were reported due to infection to novel coronavirus in Wuhan City, China (Gautam, Dilwaliya, et al., 2020; Gautam, Joshi, et al., 2020; WHO, 2020a; Wu et al., 2020). According to Xu et al. (2020), novel coronavirus was identified as SARS-CoV-2 for its likeness with SARS (severe acute respiratory syndrome) in 2002–2003. Many studies (Bherwani et al., 2021; Chen et al., 2017; Sohrabi et al., 2020) reported that the possible effects (i.e. respiratory problems, organ failure, pneumonia, etc.) had been identified in people affected by SARS-CoV-2 infections. In this regard, older individuals having severe symptoms are the primary cases with previous comorbidities (i.e. respiratory diseases, cardiovascular, etc.) (Changotra et al., 2020; Wang et al., 2020). The worst scenario was observed by the World Health Organization (WHO), and they declared it a pandemic disease with the acronym COVID-19 for coronavirus disease 2019 (WHO, 2020b). The Chinese government has taken initial actions, unfortunately, these steps were botched to control the diffusion of coronavirus, moreover have spread globally within a short time.

Figure 1 presents the transmission classification with different categories (i.e. community transmission, a cluster of cases, sporadic cases, no cases, and pending). In the present study, we are highlighting the Indian region with few points, these are: (i) India is growing very fast with positive cases; (ii) India is now in the second stage (cluster cases); (iii) India is one of the most polluted countries in the world; (iv) India is the only country affected by COVID-19 in south Asia region even more extensive than that initially identified in China. On the other hand, different countries like the USA, Brazil, Russia, etc., are rapidly registered with the exponential growth of positive cases and deaths. However, all possible measures (i.e. lockdowns, social distancing, etc.) have been applied to prevent or break the chain of COVID-19 (ECDC, 2020). The recent update from WHO has released confirmed cases (25,327,098), including deaths (848,255) worldwide.

Fig. 1

A schematic diagram of transmission classification of COVID-19 (Source: WHO, 2020b)

On the other side, many studies (Ambade et al., 2021; Bashir, 2020; Dutheil et al., 2020; Frontera et al., 2020; Zhu et al., 2020) have been done on the potential correlation between COVID-19 outbreaks and air pollutants. Moreover, every country has prepared the documentary on SARS-CoV-2 diffusion and pandemic has been stated; however, a discussion is pending about the number of death cases or death rate and severity of pneumonia encountered in few particular countries (i.e. China, Italy, India, USA, etc.) (WHO, 2020a). In India, COVID-19 is in the second phase (clusters of cases) where few states (i.e. Maharashtra, Andhra Pradesh, Tamil Nadu, Karnataka, Uttar Pradesh, New Delhi) have been reported a majority of positive cases. The total number of confirmed cases and cumulative deaths makes approximately 50% of COVID-19 instances in the Indian region. This information indicates to a precise point of the selected states that may allow for the spreading of the novel coronavirus. From an air quality point of view, the mentioned states are well known to be Indian states having the highest polluted cities or highest pollution levels. Here, we investigate the link between air pollutants and the severity of novel coronavirus, which could have vital indications towards managing the current situation and preventing the pandemic.

Materials and methods

Sample and data

The COVID-19 related data (i.e. positive cases, cumulative cases, new death cases, and cumulative cases) have been collected from the WHO (https://covid19.who.int/). Air pollutants (i.e. PM2.5, PM10, SO2, NO2) data and meteorological parameters (MPs) (i.e. temperature, humidity) have been collected from the Indian government organization called Central Pollution Control Board (CPCB) (https://cpcb.nic.in/). The data for monthly wind patterns were extracted from European Centre for Medium-Range Weather Forecasts (ECMWF) platform (Patel & Kumar, 2015) and plotted by using GRADS software. We have collected the data from the 21 states and union territories (which have at least one monitoring station) to obtain the average values of India for the measured variables in the particular month(Also explained in Fig. 2).

Fig. 2

Flow chart of data extraction and processing

Measures of variables: We have used air quality index (AQI), PM2.5, PM10, SO2, NO2, T, RH, wind speed and wind direction as a measured variables in the presented manuscript.

Techniques and data analysis

Spearman's and Kendall's rank correlation was calculated to determine the nature of the association between the parameters. The nonparametric Spearman's rank correlation can be evaluated by using the following mathematical relation:where = Difference between the rank of two parameters, = Number of alternatives.

The correlation coefficients lie in the typical range −1 (Positive correlation) to + 1 (Negative correlation).

Whereas, Kendall’s rank correlation is also a nonparametric test that estimates the ordinal association between variables. The mathematical equation for Kendall’s rank correlation is given by the following relation:where concor = Number of concordant, discor = Discordant pairs and n = number of pairs. The value of also ranges from −1 to + 1 similar to Spearman's coefficient (rs). The calculation of correlation coefficients and p values was performed by using HIMSC and MASS packages in R as suggested by the Gautam, Dilwaliya, et al. (2020)), Gautam, Joshi, et al. (2020)), Kumar, Hama, et al. (2020) and Kumar, Lin, et al. (2020).

For the assessment of the impact of lockdown enforced due to COVID-19, the air quality index (AQI) is also calculated using the methodology suggested by the (Gautam, Dilwaliya, et al., 2020; Gautam, Joshi, et al., 2020; Mahato et al., 2020) in their study. The Indian AQI system is defined as impact of individual pollutants on human health, which is combinations of two steps involvements (a) formation of individual sub-indices (SIs) of pollutants and (b) combination of sub-indices to obtain AQI. The SIs for N pollutant is calculated by using sub-index functions as follows:where

Now Aggregation of sub-indices (I) is calculated by using the following mathematical function (given below)

The mathematical expression for sub-index (I) of known pollutant concentration (C) is given by

BHI = Breakpoint concentration greater or equal to given concentration,

BLO = Breakpoint concentration smaller or equal to given concentration,

IHI = AQI value corresponding to BHI,

ILO = AQI value corresponding to BLO.

Furthermore;where .

Mahato et al. (2020) (National Air Quality Index. Central Pollution Control Board,1–44).

Result and discussion

The present trend in the Indian region straightaway revealed potential differences at the regional level for reporting cases or infections of a virus. The probably massive number of positive instances concentrated, especially in the northern part of India (Bherwani et al., 2020; Gupta et al., 2020). The number of COVID-19 positive (new and cumulative) and deaths (new and incremental) cases reported from January 2020 to August 2020 is presented in Fig. 3 (data re-elaborated from the official daily reports of WHO 2020; Data assessed on September 3, 2020). Figure 3b and c shows that the approximate number of COVID-19 positive cases and death cases due to COVID-19 in India are > 300,000 and > 10,000, respectively.

Fig. 3

Distribution of COVID-19 outbreak in India (January–August 2020). a An absolute number of positive cases (new and cumulative cases); b the Absolute number of deaths due to COVID-19 (recent deaths and cumulative deaths); c Two-point segment comparison between new COVID-19 positive cases and recent Deaths due to COVID-19 and d Two-point segment comparison between new deaths cases and cumulative deaths cases due to COVID-19

The two-point segment graph is presented in Fig. 3c, where it can be easily found that positive cases have been increased from April and its very high in August and so on. Similarly, the rate of deaths is growing with time; however, the death ratio, according to the number of positive cases, 1.78% (new cases and new deaths), and 2.19% (cumulative positive cases and cumulative deaths) only (Fig. 3c, d). Many studies (Bukhari & Jameel, 2020; Liu et al., 2020; Wang et al., 2020) indicated that confirmed positive cases were observed in lower temperatures. In contrast to the reported all studies, it can be simply observed the confirmed cases are higher in number, especially in summer or dry season (Fig. 2c, d).

In this regard, a comprehensive analysis should be done to understand the relationship between meteorological factors and air pollutants to transmit a novel coronavirus in the Indian region. We collected the air pollutants and MPs data sets from the Central Pollution Control Board (CPCB), Indian government organization, and summarized in Fig. 4a and b.

Fig. 4

Variation of the concentration profile of air pollutant a and meteorological parameters b during pandemic (January–August 2020)

From Fig. 4a, as many Indian studies (Bherwani et al., 2020; Gautam, Dilwaliya, et al., 2020; Gautam, Joshi, et al., 2020; Gupta et al., 2020) suggested on air quality and lockdown, it can be seen that the sudden fall of air pollutants due to lockdown applied during the pandemic, especially in between late March and April. On the other side, the ion and mass concentration of air pollutants is affected by ambient temperature (T) and relative humidity (RH) (Liu et al., 2020; Chow et al., 2015). Therefore, we have collected the data sets of T and RH along with air pollutant data to understand the role of air pollution in COVID-19 outbreak risk, especially in the Indian region. The ambient T and RH showed the ambient atmosphere's general profile, where they are offering a reverse relationship with each other (Fig. 4b). Some studies have shown that air pollutant’s concentration changes smoothly as atmospheric T and RH increases or decreases (Chan et al., 1992; Pathak et al., 2004). That means variation in T and RH might play a significant role in changing air pollutant concentrations. (Xie & Zhu, 2020) obtained the threshold value of T at 3° C by using the response curve COVID-19. MPs such as wind direction, wind speed, T, rain fall, pressure, solar radiation, etc. have the impacts on the removal, dispersion and chemical formation of the particulate matters (PMs). Few studies have reported a significant linkage of PM concentration and MPs. PM10 mass concentrations decreased due to an increase of T, wind speed and precipitation (Keary et al., 1998). Negative correlation of PM10 and temperature is due to traffic density and domestic heating during winter (Keary et al., 1998). Relationship of five MPs and total suspended particles (TSP) showed the low correlations during winter at Elaziğ city (Turkey) as reported by (Akpinar et al., 2008). Yadav et al. (2015) also performed a model based study to understand the variation of PM and MPs over Pune (India).

The monthly AQI shows a remarkable variation in due lockdown across India (Fig. 5). Before lockdown, the AQI varies in the range of 91.45 ± 45.21 to 285.40 ± 75.30 with an average (mean + standard deviation) of 174.19 ± 99.42, but after lockdown, the AQI levels has been drastically changed in the range of 41.87 ± 20 to 75.87 ± 42.2 with an average value 58.66 ± 8.5. The overall reduction in the pandemic periods (January–August 2020) was reported by up to 33.67% over India.

Fig. 5

Monthly variation of the air quality index (AQI) during pandemic (January–August 2020)

Maji et al. (2015) estimate the trends of air quality from 2001 to 2010 and also explained the seasonal and monthly variation of AQI over Delhi based on CPCB datasets. (Garg & Gupta, 2020) also, study the monthly variation of AQI during different months and applied the statistical approach to determine the responsible sources to degrade the AQI during 2017 over Delhi. The lower values of AQI in June, July and August 2020 are also explained in the upcoming sections.

The monthly average wind pattern is also plotted in Fig. 6 by using the ECMWF data set over India. In January (a) and February (b), the variation of air mass over the Indian subcontinents is relatively low. The air mass is blowing in the range of 1–2 m/s over India excepts some parts of Tamilnadu, Gujrat, Rajasthan and Delhi non-capital region in the typical range of 2–3 m/s. During this period (January and February 2020), the concentration of pollution parameters is high due to boundary layer dynamic (Allabakash & Lim, 2020). In March (c) and April (d), the average wind speed is slightly increased mainly over the western part of Gujrat, Tamilnadu and Indian Gangetic Plain (IGP) in the range of 2–3 m/s, but the dominantly the air mass is coming from the westerly-north westerly. But in May (e) and June (f) 2020, most of the air mass is blowing range of 2–7 m/s over the western and southern part of India.

Fig. 6

Variation of monthly wind patterns in months of January (a), February (b), March (c), April (d), May (e), June (f), July (g) and August 2020 (h) over Indian subcontinent during pandemic (January–August 2020)

Now in July (g) and August 2020 (h), the wind speed was high as compared to previous months. The western and southern parts of India experienced an air mass in the range of 4–8 m/s. The air mass from south-southwest usually dominants during June, August and September 2020. The South- south-westerly wind picks moisture from the Indian Ocean and leads to heavy precipitation over India. The process of precipitation/wet deposition of atmospheric aerosols helps to improve the AQI over India (Fig. 5) during monsoon (Garg & Gupta, 2020; Mamta & Bassin, 2010; Patel & Kumar, 2015). Coccia (2020) has examined the spread of COVID-19 with concentration of pollutants and wind speed over Italy and showed the higher risk of COVID-19 infections at low wind speed and high concentration of pollutants. The high wind speed and temperature play important role in the reduction of COVID-19 infections in Brazil. It is clear that wind speed playing a crucial role in the spread process of COVID-19 in the India.

Babu et al. (2020) obtained a significant positive correlation of wind speed and COVID-19 daily cases over Delhi, which indicate that the higher wind speed is helpful to spreading over Delhi.

The COVID-19 (New Cases, Cumulative Cases, New Deaths and Cumulative Deaths) cases show a strong and positive with maximum, average and minimum values of temperature (Table 1). The average temperature shows a stronger correlation (r values = 0.57–0.79) than maximum (r = 0.50–0.71) and the minimum (r = 0.55–0.78) temperature, which means at the lower and higher T the COVID-19 cases are less as compare average T range. A similar pattern of variation of T was also reported by (Gautam, Dilwaliya, et al., 2020; Gautam, Joshi, et al., 2020) over India. The RH shows a positive but weak correlation with COVID-19 cases in the range of 0.14–0.52 during observation. The fine particulate (PM2.5) gives a strong and negative correlation with COVID-19 cases in the range of −0.90 to −0.98 for the Spearman and −0.64 to −0.91for the Kendall correlations during calculations (Kumar, Hama, et al., 2020; Kumar, Lin, et al., 2020). The average PM2.5 (PM2.5 ave) shows a relatively high correlation than minimum PM2.5 (PM2.5 min) and maximum PM2.5 (PM2.5 Max) with COVID-19 cases (Table 1), whereas in the case of PM10, show a strong and negative correlation at the higher (PM10 max) and average (PM10 ave) levels. At the lower (PM10 min) levels, the correlation is relatively weaker. Similarly, the NO2 also, shows a strong and negative variation of correlations in the range of 0.64–0.98, it means at more than average levels of PM (PM10 and PM2.5) and NO2 the COVID-19 cases are less (Fig. 3c, d) as compared to lower levels of PM10 and PM2.5. But SO2 shows negative and non-significant correlations with COVID-19 cases in the range of 0.10–0.69 during the study, however, at lower SO2 (SO2 min) levels some significant correlation was also observed (Table 1). Kumar, Hama, et al. (2020), Kumar, Lin, et al. (2020) estimate health and economic impacts due to reduction in the PM2.5 concentration during lockdown period over Delhi, Mumbai, Bangalore, Kolkata and Chennai and found that this reduction in PM2.5 can prevent 630 premature deaths worth of 0.63 billion US Dollar (USD). Tanvir et al. also observed a negative correlation between COVID-19 cases and PM over Italy.

Table 1

Spearman's rank and Kendall's rank correlation with different variables during study

Variables	Spearman's rank correlation				Kendall's rank correlation
	New cases r values (p values)	Cumulative cases r values (p values)	New deaths r values (p values)	Cumulative deaths r values (p values)	New cases r values (p values)	Cumulative cases r values (p values)	New deaths r values (p values)	Cumulative deaths r values (p values)
Tmax	0.69 (< 0.05)	0.71 (< 0.05)	0.71 (< 0.05)	0.71 (< 0.05)	0.50 (< 0.05)	0.57 (< 0.05)	0.55 (< 0.05)	0.55 (< 0.05)
Tave	0.76 (< 0.05)	0.79 (< 0.05)	0.78 (< 0.05)	0.78 (< 0.05)	0.57 (< 0.05)	0.64 (< 0.05)	0.62 (< 0.05)	0.62 (< 0.05)
Tmin	0.75 (< 0.05)	0.78 (< 0.05)	0.77 (< 0.05)	0.77 (< 0.05)	0.55 (< 0.05)	0.62 (< 0.05)	0.59 (< 0.05)	0.59 (< 0.05)
RH max	0.52 (0.18)	0.48 (0.23)	0.50 (0.20)	0.50 (0.20)	0.43 (0.17)	0.36 (0.27)	0.40 (0.17)	0.40 (0.17)
RH ave	0.62 (0.10)	0.57 (0.13)	0.60 (0.11)	0.60 (0.11)	0.50 (0.10)	0.43 (0.17)	0.47 (0.10)	0.47 (0.10)
RH min	0.24 (0.57)	0.29 (0.49)	0.26 (0.53)	0.26 (0.53)	0.14 (0.71)	0.21 (0.54)	0.18 (0.54)	0.18 (0.54)
PM2.5 max	−0.90 (< 0.05)	−0.93 (< 0.05)	−0.92 (< 0.05)	−0.92 (< 0.05)	−0.79 (< 0.05)	−0.86 (< 0.05)	−0.84 (< 0.05)	−0.84 (< 0.05)
PM2.5 ave	−0.95 (< 0.05)	−0.98 (< 0.05)	−0.97 (< 0.05)	−0.97 (< 0.05)	−0.86 (< 0.05)	−0.93 (< 0.05)	−0.91 (< 0.05)	−0.91 (< 0.05)
PM2.5 min	−0.76 (< 0.05)	−0.74 (< 0.05)	−0.75 (< 0.05)	−0.75 (< 0.05)	−0.71 (< 0.05)	−0.64 (< 0.05)	−0.69 (< 0.05)	−0.69 (< 0.05)
PM10 max	−0.98 (< 0.05)	−0.95 (< 0.05)	−0.97 (< 0.05)	−0.97 (< 0.05)	−0.93 (< 0.05)	−0.86 (< 0.05)	−0.91 (< 0.05)	−0.91 (< 0.05)
PM10 ave	−0.98 (< 0.05)	−0.95 (< 0.05)	−0.97 (< 0.05)	−0.97 (< 0.05)	−0.93 (< 0.05)	−0.86 (< 0.05)	−0.91 (< 0.05)	−0.91 (< 0.05)
PM10 min	−0.71 (< 0.05)	−0.74 (< 0.05)	−0.73 (< 0.05)	−0.73 (< 0.05)	−0.57 (< 0.05)	−0.64 (< 0.05)	−0.62 (< 0.05)	−0.62 (< 0.05)
SO2 max	−0.12 (0.78)	−0.10 (0.82)	−0.11 (0.80)	−0.11 (0.80)	−0.14 (0.71)	−0.07 (0.54)	−0.11 (0.53)	−0.11 (0.53)
SO2 ave	−0.67 ( 0.07)	−0.62 (0.10)	−0.65 (< 0.08)	−0.65 (< 0.08)	−0.50 (0.10)	−0.43 (0.17)	−0.47 (0.17)	−0.47 (0.17)
SO2 min	−0.71 (< 0.05)	−0.67 (< 0.07)	−0.69 (< 0.05)	−0.69 (< 0.05)	−0.57 (0.06)	−0.50 (0.10)	−0.55 (0.60)	−0.55 (0.60)
NO2 max	−0.93 (< 0.05)	−0.88 (< 0.05)	−0.91 (< 0.05)	−0.91 (< 0.05)	−0.79 (< 0.05)	−0.71 (< 0.05)	−0.76 (< 0.05)	−0.76 (< 0.05)
NO2 ave	−0.98 (< 0.05)	−0.95 (< 0.05)	−0.97 (< 0.05)	−0.97 (< 0.05)	−0.93 (< 0.05)	−0.86 (< 0.05)	−0.91 (< 0.05)	−0.91 (< 0.05)
NO2 min	−0.88 (< 0.05)	−0.86 (< 0.05)	−0.87 (< 0.05)	−0.87 (< 0.05)	−0.71 (< 0.05)	−0.64 (< 0.05)	−0.69 (< 0.05)	−0.69 (< 0.05)

p values are written in the bracket bellow the correlation coefficient of r

Where, max maximum, min minimum, ave average levels of variables, RH relative humidity, T temperature

Conclusion

An association between air pollution, COVID-19 confirmed cases, and meteorological factors (T and RH) is plausible. This may have a potential impact on the transmission of the virus and the high rate of infection and mortality. In conclusion, it is observed that infectious diseases come to society with one of the key challenges during a pandemic that has a strong correlation between the periodic emergence of viral agents and socio-economic and environmental factors.

To understand the relationship between air pollutant concentration, meteorological factor, and COVID-19 cases, the nonparametric Spearman's and Kendall's rank correlation were calculated. The COVID-19 shows a strong and positive with T in the range of 0.55–0.79 over the Indian region. The outcomes indicate that the cumulative death and cumulative cases are increasing rapidly with T and RH. The RH shows a positive but weak correlation with COVID-19 cases in the range of 0.14–0.52. The PM2.5 gives a strong and negative correlation with COVID-19 cases in the range of −0.64 to −0.98 during calculations Similarly, the NO2 shows a strong and negative variation of correlations in the range of 0.64–0.98, which indicates that before lockdown (January–March 2020) the levels of PM (PM10 and PM2.5) and NO2 were high. Before the lockdown, the concentration in the pandemic periods (January–August 2020) was reported by up to 33.67% over India. The background MPs (wind Speed, wind direction, T and RH) play a very crucial role in the variation of pollution parameter and strongly influence the COVID-19 cases in India. The improvements in the levels of pollution parameters influenced the health and economic aspects due to lockdown.

The monthly ECMWF derived average wind speed is relatively low in January and February in the range of 1–2 m/s over India, except some parts of Tamilnadu, Gujrat, Rajasthan and Delhi non-capital region in the typical range of 2–3 m/s. But in July and August 2020, the western and southern parts of India experienced an air mass in the range of 4–8 m/s. The south-southwest during June August and September 2020 helps to reduce the AQI due to heavy precipitation. Policy and Managerial Implications: The proposed work's primary purpose has only supported state, national governance for future healthcare policy. The MPs and pollution parameters have a significant role in the spread of COVID-19 cases so policy-maker should consider impact of both parameter during making the policy.