key: cord-344501-kc0q6k6j authors: Zhu, Yongjian; Xie, Jingui; Huang, Fengming; Cao, Liqing title: The mediating effect of air quality on the association between human mobility and COVID-19 infection in China date: 2020-07-09 journal: Environ Res DOI: 10.1016/j.envres.2020.109911 sha: doc_id: 344501 cord_uid: kc0q6k6j BACKGROUND: Previous studies have found that human mobility restrictions could not only prevent the spread of COVID-19, but also improve the air quality because of the reduction of industrial production, transportation and traffic. It is noteworthy that air quality is also closely related to the risk of COVID-19 infection. Therefore, we aimed to assess the mediating role of air quality on the association between human mobility and the infection caused by this novel coronavirus. METHODS: We collected daily confirmed cases, human mobility data, air quality data and meteorological variables in 120 cities from China between January 23, 2020 and February 29, 2020. We applied the generalized additive model to examine the association of human mobility index with COVID-19 confirmed cases, and to assess the mediating effects of air quality index and each pollutant. RESULTS: We observed a significant positive relationship between human mobility index and the daily counts of COVID-19 confirmed cases. A unit increase in human mobility index (lag0–14) was associated with a 6.45% increase in daily COVID-19 confirmed cases, and air quality index significantly mediated 19.47% of this association. We also observed a positive relationship between human mobility index and air quality index. In the pollutant level analyses, we found significant mediating effects of PM(2.5), PM(10), and NO(2). CONCLUSIONS: Our study suggests that limiting human movements could reduce COVID-19 cases by improving air quality besides decreasing social contact. Coronavirus disease 2019 , caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), continues to spread rapidly around the world (Chen et al., 2020a; Dong et al., 2020; Gilbert et al., 2020; Zhou et al., 2020) . As of 28 May 2020, data from the World Health Organization have shown that more than 5,596,550 confirmed cases have been detected and 353,373 confirmed deaths worldwide. Most infected patients only have mild symptoms such as fever and dry cough, while some patients have fatal complications like Acute Respiratory Distress Syndrome (ARDS) (Chen et al., 2020c; Huang et al., 2020a; Sohrabi et al., 2020) . Thus, exploring risk factors affecting the transmission of SARS-CoV-2 is a meaningful study. Early studies have shown that human mobility plays an important role in the dispersal of infection since SARS-CoV-2 could spread by human-to-human transmission via direct contact or droplets (Chan et al., 2020; Kraemer et al., 2020; Li et al., 2020; Tian et al., 2020) . Therefore, to make people stay at home and decrease social contact, various public health interventions have been implemented in many countries, such as prohibiting travel in and out of cities, banning mass gatherings, and closing schools and entertainment venues (Bayham and Fenichel, 2020; Chen et al., 2020d; Chinazzi et al., 2020 ). An empirical study found that the lockdown of Wuhan city on 23 January 2020 reduced outflows from Wuhan by 56.35%, inflow into Wuhan by 76.64%, and within-Wuhan movements by 54.15% (Fang et al., 2020) . It is noteworthy that human mobility restrictions could not only prevent the spread of COVID-19, but also improve the air quality because of the reduction of industrial production, transportation and traffic (Kerimray et al., 2020; Mahato et al., 2020; Nakada and Urban, 2020; Tobías et al., 2020) . According to a study using data from 44 cities in China, the air quality index (AQI) decreased by 7.80% on average due to the implementation of travel restrictions (Bao and Zhang, 2020) . Importantly, air quality is also closely related to the risk of COVID-19 infection, since air pollutants could affect body's immunity and carry virus (Bianconi et al.; Fattorini and Regoli, 2020; van Doremalen et al., 2020; Zhu et al., 2020) . So, it is reasonable to hypothesize that air quality could partially mediate the association of human mobility with COVID-19 infection. In this study, we aimed to investigate the relationship between human mobility and COVID-19 2 infection, and to assess the mediating role of air quality on this association, based on data from 120 cities in China. 116 prefecture-level cities and 4 municipalities were included in this study, covering the majority of Chinese mainland (83.4 ° to 131.6 ° east longitude, and 20.0 ° to 51.4 ° north latitude). Fig. 1 plots the locations of these 120 cities and the cumulative COVID-19 confirmed cases in each city during our study period (i.e., January 23, 2020 to February 29, 2020). We focused our study on these cities because the air pollution data and the meteorological data were limited. We collected daily confirmed new cases for each city between January 23, 2020 (i.e., the date of Wuhan lockdown) and February 29, 2020 from the official websites of local health commissions. Human mobility was measured by the real-time intracity migration index calculated from the ratio of the number of people traveling in a city to the number of people living in that city (Bao and Zhang, 2020) . Human mobility index data were obtained from the Baidu Qianxi online platform (http://qianxi.baidu.com), which tracked human movements by mobile phone data from Baidu location-based services. Daily air quality data including air quality index (AQI), PM 2.5 , PM 10 , NO 2 , O 3 , SO 2 , and CO for each city were collected from an air quality monitor platform (https://www.aqistudy.cn/historydata/). The AQI is calculated by two steps, which could describe the level of air quality (larger AQI means more serious air pollution). The first step is to calculate the individual air quality index (IAQI) values for six pollutants (PM 2.5 , PM 10 , NO 2 , O 3 , SO 2 , and CO), and the second step is to choose the maximum of the six IAQI values as the AQI (Zhong, S., et al., 2019) . We obtained daily meteorological data on mean temperature, relative humidity, atmospheric pressure, and wind speed from the National Meteorological Information Center (http://data.cma.cn). We downloaded the dataset of daily climate data from Chinese surface stations for global exchange (http://data.cma.cn/data/cdcdetail/dataCode/SURF_CLI_CHN_MUL_DAY_CES_V3.0.html) on this platform. The surface stations could be matched with cities by the station number. Our statistical analysis included two parts. In the first part, we applied the generalized additive model (GAM) with a Gaussian distribution family to assess the association of human mobility index with daily COVID-19 confirmed cases (Hastie, 2017; Liu et al., 2020a; Xie and Zhu, 2020; Zhu et al., 2020) . Since the incubation period of COVID-19 is estimated to range from 1 to 14 days and there is also a delay (mean: 4.95 days) between the onset of illness and being confirmed according to the National Health Commission in China (National Health Commission of the PRC, 2020a; 2020b), a moving-average approach was used to capture the moving average lag effect (lag0-7, lag0-14, lag0-21) of human mobility index on daily confirmed cases of COVID-19 (Duan et al., 2019; Li et al., 2018; Zhu et al., 2020) . Our model was defined as follows: Here, log(Y it ) indicates the logarithm of the number of confirmed new cases in city i on day t (plus 1 to avoid taking the logarithm of 0). α is the intercept. HMI il denotes the (l+1)-day moving average of human mobility index (lag0-l) in city i. We controlled for meteorological variables (M il ) including mean temperature, relative humidity, atmospheric pressure, and wind speed during the same period according to previous studies, using s(•) (the thin plate spline function with the 3 degrees of freedom) to account for non-linear effect (Xie and Zhu, 2020; Zhu et al., 2020) . We also controlled the logarithm of the number of confirmed new cases in city i on day t-1 to account for potential serial correlation In the second part, we conducted the mediation analysis to test our hypothesis that air quality could partially mediate the association of human mobility with COVID-19 infection. That is, human mobility is not only directly associated with COVID-19 infection, but also indirectly related to COVID-19 infection through air quality. According to previous studies (Baron and Kenny, 1986; Chen 4 et al., 2020b) , Equation (2) an Equation (3) were constructed: Here, AQI il denotes the (l+1)-day moving average of air quality index (lag0-l) in city i. Therefore, the three Equations (1) -(3) constituted a complete mediation test based on three steps: the first step is to test whether the coefficient β in Equation (1) is statistically significant; the second step is to test whether the coefficient γ in Equation (2) is statistically significant; the third step is to test whether the coefficients δ and θ in Equation (3) were statistically significant. If β, γ, δ and θ are all significant, and the indirect effect (γ×θ) has the same sign as the direct effect (δ), our hypothesis is supported. The proportion of the association between human mobility and COVID-19 infection mediated was calculated by Equation (4) In the sensitivity analysis, we excluded Wuhan city from our data to examine the robustness of our results because the number of confirmed cases in Wuhan was much larger than that in others. The "mgcv" package (version 1.8-31) in R software (version 4.0.0) was used to conduct the GAM in this study. During the study period, our analysis included over 58,000 confirmed cases, covering 70% of confirmed cases in China. Table 2 shows the spearman correlation coefficients between human mobility index, air quality index, air pollutants, and meteorological variables. Human mobility index was positively correlated with air quality index (r = 0.05, p < 0.05), PM 2.5 (r = 0.06, p < 0.05), PM 10 (r = 0.08, p < 0.05), NO 2 (r = 0.05, p < 0.05), CO (r = 0.04, p < 0.05), mean temperature (r = 0.18, p < 0.05) and relative humidity (r = 0.07, p < 0.05), while had negative correlations with O 3 (r = -0.04, p < 0.05), SO 2 (r = -0.03, p > 0.05), atmospheric pressure (r = -0.10, p < 0.05) and wind speed (r = -0.05, p < 0.05). (1). We observed significant positive relationships between human mobility index and the daily counts of COVID-19 confirmed cases at lag0-14 and lag0-21. A unit increase in human mobility index was associated with a 6.45% (95% CI: 1.66 to 11.23) increase in daily COVID-19 confirmed cases at lag0-14, and a 18.95% (95% CI: 13.41 to 24.50) increase in confirmed cases at lag0-21. In the sensitivity analysis, the relationship between COVID-19 confirmed cases and human mobility index was robust after excluding Wuhan (Fig. 3 ). Table 3 shows the results of our mediation analysis. We found that air quality index partially mediated the association of human mobility index with COVID-19 confirmed cases, supporting our hypothesis. At lag 0-14, air quality index significantly mediated 19.47% of the relationship between human mobility and COVID-19 infection, and the direct effect of human mobility indicated that a unit increase in human mobility index (lag 0-14) was associated with a 5.15% (95% CI: 0.31 to 10.00) increase in daily COVID-19 confirmed cases. At lag 0-21, the proportion of the relationship between human mobility and COVID-19 infection mediated by air quality index was 9.98%. The direct effect showed that a unit increase in human mobility index (lag 0-21) was associated with a 17.16% (95% CI: 6 11.55 to 22.77) increase in daily COVID-19 confirmed cases. We also observed a positive relationship between human mobility index and air quality index. As shown in Table 4 , the results of the mediation analysis were still robust and the proportions of the association between human mobility and COVID-19 infection mediated by air quality index were also similar (19.14% at lag 0-14 and 8.87% at lag 0-21) after excluding Wuhan. The purpose of this paper was to determine the underlying mechanism of the association between human mobility and the spread of COVID-19 by assessing the mediating effect of air quality, different from our previous work (Zhu et al., 2020) which aimed to examine the relationship between air pollution and COVID-19 infection. In this study, we found that human mobility was positively associated with COVID-19 confirmed cases using a generalized additive model. The mediation analysis showed that air quality index, PM 2.5 , PM 10 , and NO 2 significantly mediated the association between human mobility and COVID-19 infection. However, we also found the suppressing effect of O 3 . These results providing evidence that human mobility restrictions could not only control the epidemic by decreasing social contact, but also affect the risk of COVID-19 infection by changing air quality. Our results are consistent with previous studies. First, for the positive relationship between human mobility and the spread of COVID-19, Tian et al. (2020) found that mobility restriction policies like suspending intracity public transport, banning public gatherings and closing entertainment venues were associated with reductions in case incidence of COVID-19. Additionally, a study also reported that these old-style public health measures could play an important role in the novel coronavirus epidemic (Wilder-Smith and Freedman, 2020). Second, for the association of human mobility with air quality, Bao and Zhang (2020) showed that the reductions in PM 10 and NO 2 during the city lockdown could be completely explained by variations in human mobility, and they also indicated that 44.9% and 9.3% of the variations in AQI and PM 2.5 were attributed to human mobility, respectively. The reason is that air pollution generated by consumption activities and industrial productions decreased due to the travel restrictions. Other studies reported that O 3 concentration increased during the COVID-19 lockdown, which may be caused by the increased garden and home activities (e.g. cleaning and biomass burning) infection, since ambient air pollutants could make pathogens more invasive to humans by carrying them and make people more susceptible by affecting body's immunity (Becker and Soukup, 1999; Horne et al., 2018; Xu et al., 2016) . Our study has several contributions. First, to the best of our knowledge, this is the first multi-city study to assess the mediating effect of air quality on the association between human mobility and COVID-19 infection, indicating that affecting air quality is a potential mechanism besides decreasing social contact. Second, this study contributes to the integration of literature investigating three associations (i.e., the association of human mobility with COVID-19 infection, the association of human mobility with air quality, and the association of air quality with COVID-19 infection). Third, our study provides new evidence that human mobility restrictions are useful for controlling the epidemic but could increase O 3 pollution. There are also some limitations. First, the human mobility index from Baidu does not represent the real number of people going outside. Second, subgroup analysis by age group or gender to explore the sensitive population could not be conduct due to the lack of related information. Third, only cities in China were included in this study, so our findings were not globally representative. For those countries that do not have serious pollution issues, the mediating effect of air quality may be different. Future studies are needed to overcome these limitations. 8 In conclusion, our study suggests that limiting human movements is related to the reduction in COVID-19 confirmed cases and indicates that affecting air quality is a potential mechanism besides decreasing social contact. Therefore, human mobility restrictions and other public health interventions encouraging people to stay at home are useful for controlling the epidemic. However, this is a preliminary study, and large datasets should be used to conduct global analyses in the future since the mediating effect of air quality may not exist in those countries that do not have serious pollution issues. 1 Locations of all 120 cities in our analyses and cumulative COVID-19 confirmed cases (log) in each city during our study period (i.e., January 23, 2020 to February 29, 2020). Cumulative comfirmed cases (log) Fig. 2 Percentage change (%) and 95% confidence interval of daily COVID-19 confirmed cases associated with a unit increase in human mobility index using Equation (1). Percentage change (%) and 95% confidence interval of daily COVID-19 confirmed cases associated with a unit increase in human mobility index using Equation (1) after excluding Wuhan. Climate variability and infectious diseases nexus: Evidence from Sweden Does lockdown reduce air pollution? Evidence from 44 cities in northern China. Science of The Total Environment The moderator-mediator variable distinction in social psychological research: Conceptual, strategic, and statistical considerations Impact of school closures for COVID-19 on the US health-care workforce and net mortality: a modelling study. The Lancet Public Health Exposure to urban air particulates alters the macrophage-mediated inflammatory response to respiratory viral infection Particulate matter pollution and the COVID-19 outbreak: results from Italian regions and provinces. Archives of Medical Science A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: a study of a family cluster Clinical characteristics and intrauterine vertical transmission potential of COVID-19 infection in nine pregnant women: a retrospective review of medical records. The Lancet The influence of increased population density in China on air pollution. Science of The Total Environment Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study COVID-19 control in China during mass population movements at New Year. The Lancet The effect of travel restrictions on the spread of the 2019 novel coronavirus (COVID-19) outbreak An interactive web-based dashboard to track COVID-19 in real time. The Lancet Infectious Diseases Effect of changes in season and temperature on cardiovascular mortality associated with nitrogen dioxide air pollution in Shenzhen, China. Science of The Total Environment Human mobility restrictions and the spread of the novel coronavirus (2019-ncov) in china Role of the chronic air pollution levels in the Covid-19 outbreak risk in Italy. Environmental Pollution. 114732 Preparedness and vulnerability of African countries against importations of This research was supported by the National Natural Science Foundation of China (grant no. The authors declare no competing interests.