key: cord-0793839-ci09xesw authors: Li, Li; Li, Qing; Huang, Ling; Wang, Qian; Zhu, Ansheng; Xu, Jian; Liu, Ziyi; Li, Hongli; Shi, Lishu; Li, Rui; Azari, Majid; Wang, Yangjun; Zhang, Xiaojuan; Liu, Zhiqiang; Zhu, Yonghui; Zhang, Kun; Xue, Shuhui; Ooi, Maggie Chel Gee; Zhang, Dongping; Chan, Andy title: Air quality changes during the COVID-19 lockdown over the Yangtze River Delta Region: An insight into the impact of human activity pattern changes on air pollution variation date: 2020-05-11 journal: Sci Total Environ DOI: 10.1016/j.scitotenv.2020.139282 sha: 74d33cb7e1dfbff481807a161b4963d8d2e3ba4e doc_id: 793839 cord_uid: ci09xesw Abstract The outbreak of COVID-19 has spreaded rapidly across the world. To control the rapid dispersion of the virus, China has imposed national lockdown policies to practise social distancing. This has led to reduced human activities and hence primary air pollutant emissions, which caused improvement of air quality as a side-product. To investigate the air quality changes during the COVID-19 lockdown over the YRD Region, we apply the WRF-CAMx modelling system together with monitoring data to investigate the impact of human activity pattern changes on air quality. Results show that human activities were lowered significantly during the period: industrial operations, VKT, constructions in operation, etc. were significantly reduced, leading to lowered SO2, NO x , PM2.5 and VOCs emissions by approximately 16–26%, 29–47%, 27–46% and 37–57% during the Level I and Level II response periods respectively. These emission reduction has played a significant role in the improvement of air quality. Concentrations of PM2.5, NO2 and SO2 decreased by 31.8%, 45.1% and 20.4% during the Level I period; and 33.2%, 27.2% and 7.6% during the Level II period compared with 2019. However, ozone did not show any reduction and increased greatly. Our results also show that even during the lockdown, with primary emissions reduction of 15%–61%, the daily average PM2.5 concentrations range between 15 and 79 μg m−3, which shows that background and residual pollutions are still high. Source apportionment results indicate that the residual pollution of PM2.5 comes from industry (32.2–61.1%), mobile (3.9–8.1%), dust (2.6–7.7%), residential sources (2.1–28.5%) in YRD and 14.0–28.6% contribution from long-range transport coming from northern China. This indicates that in spite of the extreme reductions in primary emissions, it cannot fully tackle the current air pollution. Re-organisation of the energy and industrial strategy together with trans-regional joint-control for a full long-term air pollution plan need to be further taken into account. At the end of the year 2019, the tragic coronavirus 2019 (COVID-19) pandemic occurred (Tian et al., 2020; Wang et al., 2020a) which has caused over 2.95 million global infections and 1931,000 deaths as of this writing (28 th April 2020). Moreover the pandemic has caused enormorous economic and social disruption (Wang et al., 2020a) . To control the rapid spread of the disease to protect people's health, China has enacted emergency response procedures, and required people to stay at home instead of going out beginning from late January 2020 (Tian et al., 2020; Wang et al., 2020b) , resulting in significant reduction in the number of vehicle kilometres travelled (VKT), industrial operations, constructions, and even restaurants in operation, etc. Human and industrial activities were reduced to essential or bare minimal only. As we know that these are general emission sources that cause air pollution, especially during this particular winter season (Zhang et al., 2019) . Thus the air quality during COVID-19 in China is much better than previous years in the same time period. The National Aeronautics and Space Administration (NASA) had published statellite imageries of the massive reduction in NO 2 over China resulting from the economic slow-down and the reduced human activities (NASA, 2020) . Other analyses have similarly found that ground-based concentrations of key pollutants fell substantially across much of J o u r n a l P r e -p r o o f 3 To understand the details of the air quality changes, we have divided the whole period into four stages: Pre-lockdown (1 st January to 23 rd January, 2020); Level I response (roughly 24 th January to 25 th February, 2020); Level II response (roughly 26 th February to 31 st March, 2020) and Level III response (31 st March onwards, 2020) . The China's Spring Festival (Chinese New Year) is covered by Level I response period. According to the 'National Emergency Response Plan for Public Emergencies' issued by China State Council, the level of early warning is based on the degree of harm that may be caused by public emergencies, the degree of urgency and development. It is generally divided into four levels: Level I (particularly serious), Level II (serious), Level III (heavier) and Level IV (general). On 24 th January 2020, the YRD region entered the first level of response, cities actively carried out epidemic prevention and control measures, including the adoption of compulsory measures in accordance with the law to stop all large-scale mass events such as bazaars and assemblies, etc. During Level I response period, the public significantly reduced the number of people and vehicles in public places, industrial enterprises, construction sites, catering enterprises and other large-scale work stoppage and closure. After two months of struggle, the number of new confirmed cases per day across the country has dropped dramatically and more and more people are being cured and discharged from hospital. As the situation of the epidemic improved, on 25 th February, 2020, the emergency response level in Anhui and Jiangsu provinces was downgraded from a primary response to a secondary response. Zhejiang and Shanghai were adjusted on 2 nd March and 24 th March, respectively. Secondary response period, under the premise of good protective measures, the YRD cities allow some industrial enterprises to resume work and resume production, construction sites can be built; blanket cancellation of the village closures, citizens can rely on the health code and wear a mask to travel locally, road traffic flow gradually increased. Since 31 st March, the YRD region enters into Level III period, with most activities gradually enetering into operation except schools, etc., of cause with strict protection measures. Here, we mainly foucus on the first three stages. J o u r n a l P r e -p r o o f 4 In this study, an integrated measurement-emission-modelling approach described in the next section including analysis of multi-pollutant observations, backward trajectory and potential source contribution analyses, estimates of pollutant emission reductions, and photochemical model simulations are adopted to conduct a comprehensive assessment of the impact of reduced human activity on air pollution reduction. We would like to achieve the following:  To study the correlation between the substantial change of human and industrial activites on the air pollution scenarios in YRD: pre-lockdown, Level I response during lockdown, Level II response during lockdown;  To study the changes of the source contributions to these air pollution scenarios which are related to local activity factors;  To figure out sources of the residual pollution and investigate policy implications for future air pollution control. To assess the effectiveness of the various emissions reductions linked with lowered human and industrial activities due to the COVID-19 restrictions, emission reductions associated with those limitations are calculated, and photochemical modelling is conducted to determine the changes in PM 2.5 attributed to specific emissions reductions. The hourly ambient mass concentrations of criteria air pollutants including SO 2 , NO 2 , CO, O 3 , PM 2.5 and J o u r n a l P r e -p r o o f 5 Agency (http://data.cma.cn/). Potential source contribution factor (PSCF) analysis is applied in this study to locate pollution sources using air mass trajectories (Duan et al., 2019; Liu et al., 2019) . It can be calculated for each 1°×1° cell by dividing the number of trajectory endpoints corresponding to samples with factor scores or pollutant concentrations greater than specified values by the number of total endpoints in the cell ( Hopke and Zeng, 1989; ) . Since the deviation of PSCF results can increase with the raise of distance between cell and receptor, a weight factor W ij is adopted in this study to lower the uncertainty of PSCF results, named weighted PSCF (WPSCF; Hopke and Zeng, 1989; Polissar et al., 1999; Zhang et al., 2019) . In this study, the TrajStat modelling system is used to analyse potential source contribution areas of PM 2.5 in Shanghai of YRD during different periods of the COVID-19 with the combination of Global Data Assimilation System (GDAS) meteorological data provided by the NCEP (National Center for Environmental Prediction). This shows us the air mass transport pathway during the three stages of lockdown. In this study, the WRF-CAMx air quality numerical modelling system is used to evaluate the improvement in air quality resulting from the reductions of human activities during COVID-19. For the mesoscale meteorological field, we adopt the WRF model Version 3.4 (https://www.mmm.ucar.edu/wrf-model-general), and J o u r n a l P r e -p r o o f 6 nested domains with two-way feedback between the parent and the nest grids. All the three domains utilise 27 vertical sigma layers with the top layer at 100 hPa. For the CAMx modelling domain shown in Figure 1 , we adopt a 36-12-4 km nested domain structure with 14 vertical layers, which are derived from the WRF 27 layers. The two outer domains cover much of Eastern Asia and Eastern China, respectively, while the innermost domain covers the YRD region. The simulation period is from 1 st January to 31 st March, 2020, during which 1 st to 5 th January was utilised for model spin-up. We focus on three stages: Pre-lockdown (Jan 6 to Jan 24, 2020); Level I response (24 th January to 25 th February, 2020) and Level II response (25 th February to 31 st March, 2020) periods. We select four major cities including Shanghai, Nanjing, Hangzhou, and Hefei to study (shown in Fig. 1 ), which are capital cities of the provinces in the YRD region. Initial and boundary conditions (IC/BCs) for the WRF modelling are based on 1°×1° grids FNL Operational Global Analysis data that are archived at the GDAS. Boundary conditions to WRF are updated at 6-hour intervals for D01. The Yonsei University (YSU) scheme (Hong et al., 2006) is applied to parameterise the boundary layer processes; the NOAH land surface scheme (Ek et al., 2003) is used to describe the land-atmosphere interactions; the Purdue-Lin microphysics scheme (Lin et al., 1983) is chosen to reproduce the cloud and precipitation processes; the RRTM long-wave and Goddard Short-wave radiation schemes (Chou et al., J o u r n a l P r e -p r o o f 7 1999; Mlawer et al., 1997 ) are adopted to refect the radiation. Anthropogenic source emission inventory for the YRD region is based on a most recent emission inventory developed by our group. Emissions for areas outside YRD in China is derived from the MEIC model (Multi-resolution Emission Inventory of China (http://www.meicmodel.org); anthropogenic emissions over other Asian regions are from the MIX emission inventory for 2010 (Li et al., 2017) . Biogenic emissions are calculated by an updated version of the Model of Emissions of Gases and Aerosols from Nature (MEGAN, v3.0, http://aqrp.ceer.utexas.edu/projects.cfm). Sea salt emissions are simulated using the OCEANIC pre-processor developed by Ramboll (http://www.camx.com/download/support-software.aspx, last access on 17 th April, 2020). We further develop a reduced emission inventory to account for the restricted human activities due to COVID-19. Estimation of emission reductions for each source sector are based on reported acitivity data and our best estimates (see details in Section 3.4). The Sparse Matrix Operator Kernel Emissions (SMOKE, https://www.cmascenter.org/smoke, last access on 17 th April, 2020) model is used to process emissions into model needed format. Prior to evaluating the impact of the restricted human activities and thus the reduced emissions, the performance of the modelling system is evaluated to ensure to reasonably reproduce the observed meteorological conditions and air quality levels. This is consistent with many other works where thermal fluxes are accurately captured. Simulations of the wind speed also achieve reasonably good correlation (0.53 -0.64). On the other hand, the agreement between observed and predicted humidity data are generally low, with negative mean bias. Overall, the model performances are consistent with observed data. In terms of CAMx model perforamcne evalution, we presented mean bias (MB), normalized mean bias (NMB), fractional bias (FB) and fractional error (FE) of PM 2.5 at 41 monitoring sites over the YRD region during our simulatin period (Table S1) Quantifying the air quality changes in response to lowered human activities and the resulting emission reductions was done using the so called Brute Force Method (BFM) (Burr and Zhang, 2011) , where a baseline scenario is simulated using unadjusted emissions (i.e. those emissions that would have occurred in absence of the COVID-19) and a reduction scenario is modelled based on the emission reductions estimation (See details in Section 3.4). For emission reductions outside the YRD region during lockdown, we apply the reduction ratio used by Wang et al. (2020b) . In both cases, the same meteorology are utilised to drive the photochemical model simulations. Through a comparative analysis of the scenarios, a relative improvement factor RF for a given atmospheric pollutant, resulting from emission reductions, can be calculated and combined with ground based observations to assess the changes in air quality associated with those emission reductions caused due to lowered human activities, J o u r n a l P r e -p r o o f where C b is the simulated pollutant concentration in the base case (µg m -3 ), C s is the pollutant concentration in the COVID-19 scenario (µg m -3 ), C o denotes the actual observed concentration at the site (µg m -3 ) and C d is the concentration reduction caused by the emissions reduction (µg m -3 ). We calculate the relative improvement factor (RF) for Level I and Level II periods, separately. These factors are applied at selected monitoring sites to reflect the changes of air quality associated with COVID-19 induced emission reductions. The Particulate Source Apportionment Technology (PSAT) coupled in the CAMx is utilised to quantify the sectoral contributions to PM 2.5 from eight source categories (aggregated from detailed source sectors listed in Table S2) We compare the hourly meteorological data of Shanghai, Hangzhou, Nanjing and Hefei during January to March from 2017 to 2020, results are shown in Table 2 . As indicated in Fig S2 and S3 , on the whole, there is no obvious change during COVID-19 in terms of air pressure, temperature and wind field from January to March, Figure S4 , it can be shown that NO 2 reduced sharply, especially during Level I period over the YRD reigon, decrease ratio of SO 2 is much lower than NO 2 , while O 3 rebuced almost all over the YRD region. In general, the level of PM 2.5 is fairly low over the entire YRD region, classified internationally at good air quality level, as is shown in (-29 .5% to -51.7%), while the amount drop of VOC is not as large as NO x , leading to a drop in titration effect towards ozone. Further policies regarding reducing regional atmospheric oxidation capacity is in urgent need to takle the increasing ozone issue in the region. To further study the air mass trajectories and source characteristics of PM 2.5 in the YRD during the three stages, PSCF is used to qualitatively analyse the source regions of PM 2.5 pollution. The PSCF threshold of PM 2.5 is set to be 75 μg m -3 , which is in line with the daily average standard. Figure Before the epidemic and during Level I response, the potential impact source areas of Shanghai mainly include three categories: 1) The local potential source areas are mainly from where the recipient city is adjacent to the neighbouring area, and are distributed in the North and Northwest of the Shanghai (Xuzhou-Suzhou-Wuxi-Changzhou-Nantong area), its WPSCF is above 0.6; 2). The regional potential sources are located in the southwest and southern cities of the recipient city, mainly including the impact of the northern cities of Zhejiang Province, Huzhou and Jiaxing, with the WPSCF is above 0.5; 3). The potential sources are located in Jiangsu. In the eastern coastal area, the polluted air mass mainly comes from the north. It is influenced by the northeast trade winds and the coast current from north to south. The polluted air mass landed in Shanghai after detouring at sea, and the WPSCF was 0.3-0.4. During the Level II period, the WPSCF results show that Shanghai is mainly affected by the contribution of cities in the region. The regional potential source is located in the southwestern city of the recipient city, mainly including Huzhou, a city in the north of Zhejiang Province, and the WPSCF is above 0.6. YRD is located in the monsoon area, affected by the monsoon climate since March, it is natural that Shanghai receives more southerly winds (southeast and southwesterly winds), which has caused polluted air masses. Industries have resumed production, and economic activities have brought more man-made sources after Level II response. Figures S6 and Table S5 shows the air mass trajectories of Pre-lockdown-C1, and indicate their sources are J o u r n a l P r e -p r o o f 16 similar. The air masses originate from Northern Mongolia and Siberia, and reaches the YRD region after passing through Beijing, Tianjin and Yellow Sea region. Through this air mass, it is shown that at all periods in 2019, the concentration of CO and PM 2.5 is higher. This shows that this air mass carried primary and secondary air pollutants through long distances. This indicates that in winter season, more large scale regional joint-control is necessary. It is important to realise the lockdown coincided with the Chinese New Year (CNY) holidays, which officially began on the 25 th January 2020. It is traditional that many industries and commercial activities will slow down two or three days before the CNY holidays, labour-intensive industries take holidays even earlier, which will last until a week after the CNY. Meanwhile, migrant workers return home and people start to visit friends and relatives. However, most people choose to stay at home during this CNY, which is the biggest difference from the holiday in previous years. Thus, the Spring Festival in 2020 is almost covered by the Level I response period. Officially, the YRD cities started its full lockdown on January 23 rd -25 th and remained in place until the end of February. In the meantime, the production activities of industries and pollutant emissions are not exactly at full capacity due to the seriousness of the COVID-19 epidemic. During the full lockdown, almost all medium and small industries except power plants and large-scale enterprises were closed. We attempt to analyse the activity levels of each respective sector. Power production and demand was subsequently significantly reduced as a result of decline in demand. According to published J o u r n a l P r e -p r o o f 17 respectively. The manufacturing sector did not actually show major slowdown: according to the industrial production data published by the bureau of statistics in the provinces of the YRD, the production of iron and other non-ferrous materials, medical and pharamaceutical remained roughly constant. The petrochemical industry, construction industry, facility manufacturing were strongly affected and hampered by both the upstream and downstream chain. All other industries were almost brought to a standstill as a result of these supply chain movements. The impact of the COVID-19 epidemic on industrial production can be seen in the sharp decline of 29% and 32% in industrial electricity consumption in January-February 2020 in Anhui and Zhejiang provinces. In the non-industrial sectors, citizens were at a restricted movement order and they could not go out to the streets, leading to a sharp decline in vehicle and public transport. Data published by the bureau of statistics of Anhui and Zhejiang provinces showed that passenger traffic in both provinces dropped significantly by 50% in January-February 2020 compared to the same period in 2019. While traffic flow monitoring data from Bengbu and Changzhou cities showed a 75% decline in the first-level response period and a 50% decline during the second-level response period, compared to pre-epidemic period. As to the total departures from China's 25 busiest airports (Flightradar 24, https://www.flightradar24.com/), the decline ratio is 80% and 60%, respectively. The construction and service industry came to almost a total halt: restaurants, construction sites, laundries, automobile service, renovation and refurbishment were all stopped proactively or involuntarily. It is interesting to know that during the epidemic, there is an observable increase in the use of solvents and volatile organic compound domestically and in the hospital, plus a small increase in domestic combustion and biomass burning during the first-level response. Based on changes of the activity data caused due to the lowered human activity, we estimated the emissions changes, as shown in Table 3 . Figure 6 shows the contribution from eight source sectors (see definition in Table S2 were also significantly lowered during Level I response period, with relative contribution of 2.6% to 7.7% and 3.9% to 8.1%, repectively, due to restricted construction and travelling. In contrast, relative contribution from residential sources (RES) are higher during lockdown periods (12.1% to 28.1% for Level I and 13.7% to 28.5% J o u r n a l P r e -p r o o f 20 for Level II), exceeding contributions from MOB and DST. Fig. 6 Source contribuitons to PM 2.5 at Shanghai, Hefei, Hangzhou and Nanjing during pre-lockdown, Level I and Level II response periods. Figure 7 shows the contributions from four major source sectors to PM 2.5 during the three stages over the YRD region. As shown in the figure, industry and residential contributions to PM 2.5 is significant, while mobile source contribution is most obvious during pre-lockdown. These results indicate that re-organisation of the industrial structure, transportation structure is urgently necessary to improve future air quality. During the COVID-19 control period, human activities have been lowered greatly, causing significant reductions in industrial operations, VKT, constructions and other activities which further leads to emissions reduction. During the most stringent Level I response period, the primary pollutants like SO 2 , NO x , PM 2.5 and VOCs have been reduced by 26%, 47%, 46% and 57%, this caused a decrease of PM 2.5 concentrations by 25.4% to 48.1% over the YRD region. However, the daily PM 2.5 still ranges between 15-79 μg m -3 during the lockdown period, and O 3 rebounded by 20.5% simultaneously. Source apportionment results show that PM 2.5 during lockdown periods mainly comes from industry and residential sources. This indicates that the co-control of PM 2.5 and O 3 is quite challenging, more stringent measures like adjustment of enrgy structure and industrial J o u r n a l P r e -p r o o f 22 structure, more stringent regional joint-control within YRD and even between YRD and northern China areas should be pushed forward to achive a better air quality. J o u r n a l P r e -p r o o f Ooi collected data and conducted the data analysis Wang performed modelling work Source apportionment of fine particulate matter over the Eastern U.S. Part I: source sensitivity simulations using CMAQ with the Brute Force method A 3-dimensional eulerian acid deposition model physical concepts and formulation Implementation of Noah land surface model advances in the National Centers for Environmental Prediction operational mesoscale Eta model A new vertical diffusion package with an explicit treatment of entrainment processes A study of the sources of acid precipitation in Ontario Sulfate formation during heavy winter haze events and the potential contribution from heterogeneous SO 2 + NO 2 reactions in the Yangtze River Delta region An Integrated Source Apportionment Methodology and Its Application over the Yangtze River Delta Region Evaluation of the effect of regional joint-control measures on changing photochemical transformation: a comprehensive study of the optimization scenario analysis. Atmospheric Chemistry and Physics MIX: a mosaic Asian anthropogenic emission inventory under the international collaboration framework of the MICS-Asia and HTAP. Atmospheric Chemistry and Physics Bulk parameterization of the snow field in a cloud model Pollution characteristics and source apportionment of ambient PM 2.5 during four seasons in Yantai City Radiative transfer for inhomogeneous atmospheres: RRTM, a validated correlated-k model for the longwave The aerosol at Barrow, Alaska: long-term trends and source locations An investigation of transmission control measures during the first 50 days of the COVID-19 epidemic in China A novel coronavirus outbreak of global health concern Severe air pollution events not avoided by reduced anthropogenic activities during COVID-19 outbreak Pollution Characteristics and Regional Transport of Atmospheric Particulate Matter in Beijing from This study was financially supported by the National Natural Science The authors declare that they have no conflict of interest.