Abstract
Emerging evidences have confirmed effects of meteorological factors on novel coronavirus disease 2019 (COVID-19). However, few studies verify the impact of air pollutants on this pandemic. This study aims to explore the association of ambient air pollutants, meteorological factors and their interactions effect confirmed case counts of COVID-19 in 120 Chinese cities. Here, we collected total confirmed cases of COVID-19 by combining with meteorological factors and air pollutants data from 15th January 2020 to 18th March 2020 in 120 Chinese cities. Spearman correlation analysis was employed to estimate the association between two variables; univariate and multivariate negative binomial regression analysis were applied to explore the effect of air pollutants and meteorological parameters on the COVID-19 confirmed cases. Positive associations were found between the confirmed cases of COVID-19 and carbon monoxide (CO), aerodynamic particulate matter with aerodynamic diameter ≤2.5 μm (PM2.5), relative humidity (RH) and air pressure (AP). And negative association was found for sulfur dioxide (SO2). In addition, multivariate negative binomial regression analysis suggested that confirmed cases of COVID-19 was positively correlated with ozone (O3) in lag 0 day while it was negatively associated with wind velocity (WV) in lag 14 days, and the pollutants-meteorological factors interactions also associate with COVID-19. In conclusions, air pollutants and meteorological factors and their interactions all associate with COVID-19.
1. Introduction
In December of 2019, a novel coronavirus disease (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) reported in Wuhan City, P. R. China (Wu et al., 2020; Zhou et al., 2020). COVID-19 outbreak has rapidly spread across China and many other countries (Holshue et al., 2020; Huang et al., 2017; Lai et al., 2020; Wang et al., 2020a), and emerging evidence indicates that the etiologic agent of COVID-19 may affect the physical and mental health of us like cough, myalgia, loss of appetite, have a fever, olfactory and gustatory dysfunctions, cause neurological complications (De Felice et al., 2020; Lechien et al., 2020) and anxiety (Wang et al., 2020b). Due to the great transmissibility and wide spreading, COVID-19 has officially been declared a pandemic by the World Health Organization (WHO) on March 11th, 2020 (https://www.who.int/emergencies/diseases/novel-coronavirus-2019). By March 18th, 2020, COVID-19 have become prevalent in more than 160 countries around the world, with a total of 191127 confirmed cases worldwide (https://www.who.int/emergencies/diseases/novel-coronavirus-2019/situation-reports/); and 80928 confirmed cases have been reported in China (http://www.nhc.gov.cn/xcs/xxgzbd/gzbd_index.shtml). This makes the world face a huge public health crisis.
Air pollutants can significantly affects health, and WHO announces 4.6 million individuals were died annually from diseases and illnesses directly related to poor air quality (Tobias and Molina, 2020). Particularly, increasing studies explore the relationship between air pollution and infectious diseases. A previous study reported that there was no relationship between air pollution and SARS outbreaks (Cai et al., 2007; Tobias and Molina, 2020; Yusuf et al., 2007). A positive association was observed between air pollutants (sulfur dioxide, nitrogen dioxide, carbon monoxide and ground-level ozone) and SARS case fatality (Cui et al., 2003). Another study in Brisbane suggested that high concentrations of O3 and PM10 were associated with more pediatric influenza cases (Xu et al., 2013).
Meteorological factors are regarded as an independent factor affecting the spreading of infectious diseases (Cai et al., 2007; Xie and Zhu, 2020). Although the association between virus that result in respiratory disease (SARS, influenza viral and SARS-CoV-2) and meteorological factors had been explored (Yusuf et al., 2007), the results of these studies were controversial or inconclusive. For example, Tan’s study showed sharply change of ambient temperature was correlated to increase the risk of SARS (Tan et al., 2005), while another study by Yao claimed ambient temperature had no significant effect on the transmissibility of SARS-CoV-2 (Yao et al., 2020a) and more studies indicated that the increase of average temperature contributed to an average relative risk (RR) and the decrease of daily SARS confirmed case counts (Cai et al., 2007). These inconsistent results implied that whether the interaction of pollutants-meteorological have influenced the transmission of these diseases.
The trends of the SARS outbreak may be influenced by air pollutants and meteorological factors which would have changed the host susceptibility and the survival time of SARS-CoV in vitro (Cai et al., 2007). COVID-19 infections can cause serious respiratory diseases, the common symptoms at the beginning includes fever, cough, muscle aches or fatigue, etc. which is similar to severe acute respiratory syndrome coronavirus (SARS-CoV) (Liu et al., 2020b; Siordia, 2020), and the similar effect of air pollutants or meteorological factors on confirmed COVID-19 cases was also observed (Jiang et al., 2020; Liu et al., 2020a; Oliveiros et al., 2020; Wang et al., 2020c; Zhu et al., 2020); however, the interaction effects of air pollutants and meteorological factors remained out of consideration. In this study, we investigate the association of ambient air pollutants and meteorological factors effects and, especially their interaction effects on COVID-19 confirmed cases, for extending and further understanding the relations of transmission of this disease and environmental factors.
2. Methods
2.1 Data collection
We observed the numbers of the confirmed cases in relation to air pollutants and meteorological parameters in 120 Chinese cities (to better fit the model, we delete the data of Wuhan since confirmed case counts are large in it and the data of Wuhan may reduce the robustness of the model) that covered the majority of the Chinese mainland with no less than 5 cases during 15 January 2020 to 18 March 2020. Daily confirmed case counts and daily air pollutants were collected from the official websites of Harvard University (https://dataverse.harvard.edu). Meteorological data were obtained from China Meteorological Data Sharing Service System (http://data.cma.cn/site/index.html), respectively. Air pollutants data included particulate matter with aerodynamic particulate matter with aerodynamic diameter ≤ 2.5 μm (PM2.5), nitrogen dioxide (NO2), and sulfur dioxide (SO2), carbon monoxide (CO), ozone (O3). Meteorological parameters including average temperature (AT), diurnal temperature range DTR, relative humidity (RH), wind velocity (WV), air pressure (AP), precipitation (PRE), and hours of sunshine (HS). The means of daily pollutants concentration and meteorological factors value in every city were calculated to represent the responding variable in that city during the whole period. Therefore, there is no need to have ethical review.
2.2 Statistical analysis
The descriptive analysis was performed for all the data to explore the characteristics of total confirmed case counts, the average concentration of O3, NO2, SO2, CO, and PM2.5, AT, DTR, RH, WV, AP, PRE and HS. Meanwhile, we account for the lag effect and the latent period of COVID-19, we evaluated the associations in lag 0, 3, 7 and 14 d between confirmed case counts and these air pollutants and meteorological factors (Liu et al., 2020a). Spearman correlation analysis was employed to explore the correlation of meteorological factors and air pollutants; associations between air pollutants, meteorological factors and confirmed case counts. We adopted generalized linear models with negative binomial regression to analyze the associations between air pollutants factors (O3, NO2, SO2, CO, PM2.5), meteorological factors (AT, DTR, RH, WV, AP, PRE and HS) and the total counts of confirmed case as this count data is clustered and may be over-dispersed. We first establish a univariate negative binomial regression models to analyze the air pollutants and meteorological factors and the confirmed case counts, then we conducted multivariate negative binomial regression models that including all above-mentioned factors since the independent variables are not serious related (all |rs|<0.9). The core negative binomial regression equation is: Log(λ) = β0 + β1X1 + β2X2 + … + βmXm, λ. here is the confirmed case counts, and X1, X2, Xm are the factors that affect λ, and β1, β2, βm are the correlation coefficient of X1, X2, Xm. To explore the interaction effects, we force Xi*Xj as an interaction term in the previous model, where Xi represents an air pollutant and Xj represents a meteorological factor. We put only one interaction term in one model fit for cautiousness.
All the statistical analyses were performed by employing SPSS version 22.0 (SPSS Inc., Chicago, IL, USA) and R software (version 3.6.3) with two-sided P value < 0.05 considered as statistically significant.
3. Results
3.1 Study area
Figure 1 showed the locations of these 120 cities in China mainland and the total confirmed cases in each city was collected in March 18, 2020. The latitude range of these 120 cities is 20°00’to 47°23’N, and the longitude range is 99°11’to 130°55’E.
3.2 Description of confirmed COVID-19 cases, meteorological factors and air pollutants
The maximum and minimum values, median, interquartile range, and other indicators number of confirmed cases, ambient air pollutants and meteorological factors of the assigned 120 cities were described in Table 1. Among these cities, during our study period (Jan 15 to Mar 18, 2020), the median of daily concentrations of CO, NO2, SO2, O3 and PM2.5, were 0.82 mg/m3, 20.89 (μg/m3, 8.23 (μg/m3, 56.38μg/m3, and 43.58 μg/m3, respectively. The ranges of RH, AT, PRE, AP, HS, WV, DTR were 35.42%~85.50%, -12.48~21.29 (°C), 0.00~6.82 mm, 758.41~1024.43 hPa, 1.60~9.20 (h), 1.00~5.05 (s/m), 5.13~18.37 (°C), respectively.
3.3 Correlation analysis of air pollutants and meteorological factors
Spearman correlation coefficient (rs) between air pollutants concentrations and meteorological factors were demonstrated in Figure 2 The larger circle indicates stronger relation of two variables. Figure 2 showed rs for any two variables are less than |0.9|. Therefore, in the multivariate negative binomial regression analysis, all the variables were integrated into the model.
3.4 Effects of air pollutants and meteorological parameters on confirmed case counts
The scatter plot displayed the relationship between pollutants, meteorological parameters and confirmed case counts for lag 14 d in Figure 3. There were positive associated between the confirmed case counts and PM2.5 (P=0.024), RH (P=0.001), PRE (P=0.015) and AP (P<0.001); while there were negative between the confirmed case counts and SO2 (P<0.001) and HS (P=0.018). Meanwhile, we haven’t found any correlations between the confirmed case counts and CO, NO2, O3, AT, WV, and DTR.
3.5 Negative binomial regression model analysis
Univariate negative binomial regression results were shown in Figure 4. CO, NO2, SO2, O3, PM2.5, RH, AT, AP, HS, WV were significantly associated with confirmed cases counts while PRE have no effect on it in lag 0, 3, 7 and 14 d (P<0.05); DTR have except an effect on confirmed cases counts in lag 7 and lag 14 d. Multivariate negative binomial regression model are shown in Figure 5, CO, PM2.5, RH and AP are positive associated with confirmed case counts but SO2 negative associated with it in lag 0, 3, 7, and 14 d. However, O3 have a positive effect on confirmed case counts in lag 0 d and WV have a positive effect on confirmed case counts in lag 14 d.
3.6 Interaction effects between air pollutants and meteorological factors
The significant interaction effects were showed in table 2. As it shows, most of the CO-meteorological factor interaction terms were listed. Except for CO*DTR and CO*DHS, other CO-meteorological factor interaction effect is positive to the total of COVID-19 conformed case. We also identify that PM2.5*DARH, PM2.5*DPRE, SO2*DAAP and NO2 * DAT are positive to the number of COVID-19 confirmed cases, and O3 * DAWV are negative to the number of COVID-19 confirmed cases.
4. Discussion
In this paper, we explored the relationship between air pollutants and meteorological factors and their interactions on confirmed cases counts in 120 Chinese cities, We focused on these cities due to the meteorological data we have obtained was limited (not including all the cities in China). Our results indicated that air pollutants and meteorological factors and their interactions all associate with COVID-19.
Many previous studies have showed that the effect of meteorological factors on infectious diseases (such as influenza, SARS, MERS, COVID-19) (Cai et al., 2007; Liu et al., 2020a; Prata et al., 2020; Tan et al., 2005; van Doremalen et al., 2013; Yusuf et al., 2007), however, the effect of air pollutants on infectious diseases especially COVID-19 is out of consider. Understanding the relationship between ambient air pollutants exposure and health can help us take protective measures to reduce the risk of exposure and disease. Some researchers have found that the higher concentration of PM2.5 the greater healthcare use for acute lower respiratory infection (Horne et al., 2018) and PM air pollutants was associated with acute respiratory illness (Pope, 1989). A Ecologic analysis carry out in Chinese demonstrated that air pollutants have a positive effect on SARS case fatality (Cui et al., 2003). Our research shows that, air pollutants of CO, PM2.5 was positively while SO2 inversely associated with confirmed COVID-19 cases; NO2 and O3 may positively associated with confirmed COVID-19 cases (Univariate negative binomial regression reach significance while multivariate negative binomial regression and Spearman correlation didn’t reach significance). Within our perception air pollutants should not affect survival of COVID-19 in vitro, though it may affect host local resistance. Therefore, we speculate that exposure to these pollutants makes the body more likely to infect SARS-CoV-2 and aggravate the severity of its infection. For example, exposure to ozone modulates human airway protease/ anti-protease balance contributing to increased influenza A infection, but epidemiological studies do not shed light on why air pollutants promotes these processes and it only reveals relationships, hence, further investigations are needed to explain it.
Weather variables like air temperature, diurnal temperature range relative humidity, wind velocity, air pressure, precipitation and hours of sunshine are important factors influencing survival of microorganisms in the external environment (Anwar et al., 2019; Bezirtzoglou et al., 2011; Casadevall, 2020; Chan et al., 2011; Liu et al., 2020a). For decades, experts have issued warnings that climate change may be related to changes in the epidemiology of infectious diseases (Casadevall, 2020). Global warming, rising sea levels, and climate change will be accompanied by new and unknown diseases (Bezirtzoglou et al., 2011; Casadevall, 2020). Weather variables can change the host’s behavior (such as spending more time indoors or outdoors), and it can change the host’s defense mechanism (inhalation of dry and cold air will cause damage to the cleaning function of cilia) or change the infectivity and stability of the virus, a study found that it seemed that the confirmed COVID-19 case number increase with temperature increasing in the range of −20°C~20°C. However, with DTR in the range of 5°C~15°C the confirmed case number decreased, and the relationship between it and absolute humidity shows a non-linear relationship (Liu et al., 2020a). In addition, other studies have shown that when the average temperature is below 3°C, it is positively correlated with confirmed COVID-19 case counts (Xie and Zhu, 2020). A Brazilian study shows that a negative linear relationship between temperatures and daily cumulative confirmed cases of COVID-19 in the range from 16.8 °C to 27.4 °C (Xie and Zhu, 2020). Besides, a study involving 244 cities in China shows no significant associations between relative humanity, maximum temperature, minimum temperature and cumulative incidence rate or R0 of COVID-19 (Yao et al., 2020b). The temperature range of this work is -12.48 ~ 21.29 °C, and our results show that RH is positively correlated with confirmed COVID-19 cases and AT, DTR may be related to confirmed case counts. This reveals that SARS-CoV-2 can adapt to higher temperatures. In addition, the experimental results of the University of Hong Kong show that the virus is very stable at 4 °C, but as the culture temperature increased to 70 °C, the virus inactivation time was shortened to 5 minutes (Chin et al., 2020). The relationship between coronavirus and RH is not a monotonic linear relationship, and the survival rate is greater at low RH (20%) and high RH (80%) than at moderate RH (50%).
Wenqi Hu’s (Hu et al., 2018) research suggests low wind velocity may increase mumps risk of Mumps. The research by Cai QC et al. shows that the larger the daily average wind velocity, the smaller the secondary attack rate (Cai et al., 2007). Our multivariate negative binomial regression (lag 14 d) showed with higher wind velocity may decrease confirmed case counts, consistent with previous research results. Wind can affect how long droplets are suspended in the air. In indoor environments, high wind speeds mean that the rooms are well ventilated and in outdoor environments, high wind speeds are beneficial to the dilution and removal of droplets, both of which can shorten their suspension time in the air (Cai et al., 2007). Our work found that air pressure may increase the confirmed case counts, but previous studies have reached the opposite conclusion that high air pressure may reduce the suspending time of SARS-CoV in the air (Cai et al., 2007). Our results indicate that hours of sunshine may reduce the confirmed case counts, which may be related to the sensitivity of SARS-CoV-2 to ultraviolet rays. About rainfall, some epidemiological data generally seem to support the relationship between influenza virus infection and precipitation (Pica and Bouvier, 2012). A study found a significant association between rainfall and influenza virus infection in India, and in the dry season there are few reports of infection (Agrawal et al., 2009). And association between precipitation and respiratory syncytial virus infection has also been seen in other studies. In our study, Spearman correlation analysis showed PRE was positive relate to confirmed case counts, however, univariate and multivariate negative binomial regression analysis did not find the relationship between precipitation and confirmed case counts as this finding results from investigations in only Chinese areas, further studies need to be conducted.
Although to our knowledge there was no report indicating the interaction effect between air pollutants and meteorological factors existing in the onset of respiratory infectious disease, the interaction effects of these factors were identified in noninfectious disease (Huang et al., 2017) or unclassified cases (Chen et al., 2019). In this study, we observed the interaction effects on total number of COVID-19 conformed case, especially the pollutant CO interacted with most other meteorological factors, which indicates that the effect of CO in increasing total number of COVID-19 conformed case varies depending on the value of a named meteorological factor. The interaction effects of PM2.5 by DARH or DPRE, NO2 by DAT were identified at significant level, even though these coefficients were much closed to zero. The interaction effects of CO*DTR is negative (−0.35), which indicates there was an antagonistic effect between these two factors. The exact mechanics behind all the interaction effects are still unclear for the limitation of this study.
This is the first time that we have comprehensively evaluated the effect of environmental pollution factors and meteorological factors on the confirmed COVID-19 case number, which gave us a further understanding of the risk factors of COVID-19. Of course, the limitation of our research is that we only analyzed the data of 120 cities in China. Further research should cross more countries and regions to clarify the effect of environmental air pollutants on COVID-19.
5. Conclusion
Our work implicates that air pollutants and meteorological factors and their interactions all associate with COVID-19, the increased number of confirmed COVID-19 case was accompanied by exposure to higher concentrations of CO, PM2.5 and higher RH, AP, and increased confirmed COVID-19 cases was accompanied by exposure to lower concentrations of SO2. Our research shows that it may be possible to reduce the transmission of COVID-19 by 1. reducing the concentrations of CO and PM2.5 in the air, 2. lower the RH and AP, but government public health precautions are essential.
Credit author statement
Jianli Zhou and Linyuan Qin: conceptualization, methodology, data curation, writing-original draft preparation, visualization, investigation. Nan Liu: supervision, validation, writing-reviewing and editing.
Data Availability
The [COVID-19 confirmed cases counts and air pollutants] data used to support the findings of this study were freely available from https://dataverse.harvard.edu. The [meteorological factors] data used to support the findings of this study were freely available from http://data.cma.cn/site/index.html if you registered users with real names in education and scientific research.
Declaration of Competing Interest
All the authors in the study declare no potential conflicts of interest.
Acknowledgements
The authors gratefully acknowledge the financial supports by the National Natural Science Foundation of China (Nos. 81872584 and 81273078), the National 863 Young Scientist Program (No. 2015AA020940), the Science and Technology Program of Guangzhou (No. 201704020056), and the Scientific Research Project for University of Education Bureau of Guangzhou (No. 201831841).