key: cord-0960976-vup1ufti authors: Liu, Yue; Wang, Yinan; Cao, Yang; Yang, Xi; Zhang, Tianle; Luan, Mengxiao; Lyu, Daren; Hansen, Anthony D. A.; Liu, Baoxian; Zheng, Mei title: Impacts of COVID‐19 on Black Carbon in Two Representative Regions in China: Insights Based on Online Measurement in Beijing and Tibet date: 2021-06-03 journal: Geophys Res Lett DOI: 10.1029/2021gl092770 sha: dcb668148b8b5893463f15e56d96dfe3c7cecb35 doc_id: 960976 cord_uid: vup1ufti Under the influence of Coronavirus Disease 2019 (COVID‐19), China conducted a nationwide lockdown (LD) which significantly reduced anthropogenic emissions. To analyze the different impacts of COVID‐19 on black carbon (BC) in the two representative regions in China, one‐year continuous online measurements of BC were conducted simultaneously in Beijing and Tibet. The average concentration in the LD period was 20% higher than that in the pre‐LD period in Beijing, which could be attributed to the increase of transport from southwestern neighboring areas and enhanced aged BC. In contrast to megacity, the average concentration of BC in Tibet decreased over 70% in the LD period, suggesting high sensitivity of plateau background areas to the anthropogenic emission reduction in South Asia. Our study clearly showed that BC responded very differently in megacity and background areas to the change of anthropogenic emission under the lockdown intervention. Disease 2019 (COVID-19) on black carbon (BC) was significantly different in Beijing and Tibet • Transport from southwestern areas and enhanced aged BC could contribute to higher BC in Beijing during the lockdown period • Large decrease of BC in Tibet under showed high sensitivity to the anthropogenic emission reduction in South Asia Supporting Information may be found in the online version of this article. Correspondence to: Y. Wang and B. Liu, wangyinan@mail.iap.ac.cn; liubaoxian28@163.com Liu, Y., Wang, Y., Cao, Y., Yang, X., Zhang, T., Luan, M., et al. (2021) . Impacts of COVID-19 on black carbon in two representative regions in China: Insights based on online measurement in Beijing and Tibet. Geophysical Research Letters, 48, e2021GL092770. https://doi.org/10.1029/2021GL092770 Different from Beijing, a typical urban site and a megacity, Tibet is a representative background area which is largely influenced by the cross-border regional transport from South Asia, especially during non-monsoon seasons (November-May) (Cong et al., 2015; Kang et al., 2019) . Previous studies reported that primary emission in India decreased significantly during the LD period from March to May, resulting in reduction of BC concentration (Kalluri et al., 2021; Panda et al., 2021) . Under the nationwide LD in South Asia in 2020, it is essential to investigate the sensitivity of BC in Tibet to the emission reduction, which has not been investigated to our best knowledge. In this study, 1-year continuous online measurement of BC was conducted simultaneously in Beijing and Tibet. The aethalometer model was used for two sites while the receptor model (Positive Matrix Factorization [PMF] ) was applied for Beijing. This is mainly because more online measurements for different chemical components could be carried out in Beijing, while it is very hard to deploy multiple instruments in Tibet. The aim of the study is to investigate the different response of BC in Beijing and Tibet under the impact of COVID-19 as well as the major factors impacting BC concentration and temporal variation during the LD period, which could help to understand the sensitivity of urban and background areas to the change of anthropogenic emission under the LD intervention. To analyze the different influence of COVID-19 on BC in megacity and remote areas in China, Beijing and Tibet were selected as two representative sites in this study. Figure S1 ). The BC mass concentration by AE33 is calculated based on the change of optical attenuation at 880 nm with a time resolution of 1 min (Drinovec et al., 2015) . The detailed information of AE33 measurement can be found in Drinovec et al. (2015) . As other online measurements of multiple chemical species can be carried out at the Beijing site, the receptor model could be applied for BC source apportionment in Beijing. The online measurements of other PM 2.5 components include multiple elements by the Xact™ 625 (Cooper Environmental Services LLC), water soluble ions by the URG-9000 ambient ion monitor, and organic carbon (OC) by the Semi-continuous OCEC Carbon Aerosol Analyzer (Sunset Laboratory Inc.). These online measurements were conducted simultaneously with BC at the BJMEMC site. The detailed information of instruments could be found in and Manigrasso et al. (2010) . The quality control was discussed in the supporting information. The operation and maintenance of the above multiple online instruments is very challenging at the YBJ site, especially during the LD period. Therefore, only AE33 measurement was made in Tibet. To compare the diurnal variation of BC and the correlation between BC and secondary inorganic aerosol (including sulfate, nitrate, ammonium; hereafter SNA) in 2020 with those in previous years, the long-term online measurements of BC by AE33 and SNA by the In-situ Gas and Aerosol Composition monitoring system (Model S-611, Fortelice International Co., Ltd.) at the Peking University site (the PKU site, 39°59′N, 116°18′E; 52 m MSL) from January 2016 to November 2019 were used for this analysis. The PKU and BJMEMC sites are both representative urban site of Beijing, and two sites are 7 km apart. The BC data from these two urban sites are comparable (see Figure S2 ). The aethalometer model was applied for source apportionment of BC in both Beijing and Tibet. It is assumed that BC absorption can be attributed to fossil and non-fossil combustion. Therefore, abs, 880 abs, 880 fossil abs, 880 non fossil The wavelength dependence of absorption coefficient is proportional to λ −α , where λ is the wavelength and α is the angstrom exponent (AAE). The contribution of BC from fossil fuel and non-fossil fuel can be derived as follows: where α fossil and α non-fossil represents the AAE for BC from fossil fuel (α fossil = 1.0) and BC from non-fossil fuel (α non-fossil = 2.0). How to determine α values can be found in the supporting information. The PMF model requires input of source tracers, which were not measured in Tibet in this study. Therefore, PMF was only applied to investigate BC sources in Beijing. The principle and detailed information of this model could be found in Paterson et al. (1999) . Source profiles of PMF are shown in Figure S3 . The identification of each factor can be found in the supporting information. Both two models were applied in Beijing. The fossil fuel sources of BC mainly includes traffic emission and coal combustion in China (Q. ). However, the AAE of BC from coal combustion was larger than 1.0 (Liu et al., 2018) , thus the fossil fuel contribution apportioned by the aethalometer model mainly includes traffic source . Therefore, the contribution of traffic source by the PMF model and fossil fuel combustion source by the aethalometer model were compared (seen in Figure S4 and Table S1 ). It could be seen that the two models agreed well with each other with R 2 over 0.9. The total potential source contribution function (TPSCF) method was conducted for potential source region analysis and the detailed information of TPSCF is provided in the supporting information. To investigate the influence of COVID-19 on BC in Beijing and Tibet, different periods were defined according to the timeline of government regulations. On December 31, 2019, the new coronavirus was found and reported in Wuhan, China. The first LD of city started on January 23 and the first level response to major public health emergencies was issued in all 31 provinces from January 29, 2020. Different from China, which experienced a recovery period in March 2020, the lockdown of India started on March 25 and lasted until the end of May 2020 (Chatterjee et al., 2021) . Therefore, five special period was defined as follows The temporal variations of BC concentration in Beijing and Tibet from November 2019 to October 2020 are shown in Figure 1a . Figure 1b shows the average BC concentration in different periods and the log-normal distributions are shown in Figure S5 . It could be seen that the influence of COVID-19 on BC were significantly different in megacity and plateau background site. In Beijing, the average concentration of BC in the LD period was 2.22 ± 1.66 μg/m 3 , which was higher than that in the pre-LD period by about 20%. The log-normal distribution also showed similar variation trend, with the mean value increased from 0.09 in pre-LD period to 0.20 in LD period ( Figure S5 ). Unlike the decrease of BC in other cities (Chen et al., 2020; Xu et al., 2020) , the increase of BC average concentration in the LD period in Beijing could be mainly attributed to high concentrations during two haze episodes. Without the two haze episodes, BC in the LD period was actually lower than that in the pre-LD period by about 25% ( Figure S6 ). Three haze episodes (HP) and two non-haze periods (NHP) were selected with daily PM 2.5 concentration larger than 75 μg/m 3 for HP and less than 75 μg/m 3 for NHP ( Figure S7 ). As shown in Figure 1c , HP1 and NHP1 were selected for the pre-LD period while HP2, HP3, and NHP2 were selected for the LD period. The date and the average, maximum, and minimum concentration of BC during different episodes are listed in Table S2 . It could be seen that BC in NHP2 during the LD period was lower than that in NHP1 during the pre-LD period. Compared with the pre-LD period, the haze episodes in the LD period were characterized by longer duration with lower maximum value of BC. In contrast to Beijing, the average concentration of BC in Tibet decreased over 70% from pre-LD (0.81 ± 1.16 μg/m 3 ) to LD period (0.23 ± 0.17 μg/m 3 ). Moreover, BC concentration showed an increasing trend after the LD-IN period, with an average concentration two times higher than that in the LD period LIU ET AL. (0.51 ± 0.46 μg/m 3 ). The monthly variation of BC in our study was different from previous studies in Tibet. Figure 1d showed the monthly variation of BC at the YBJ site in this study and at Lulang in 2008 and Namco site in 2013 in Tibet (X. Zhang et al., 2017; Zhao et al., 2013) . In 2013, BC concentration increased significantly from November to May and decreased from June to September (X. Zhang et al., 2017) . Similar temporal variation could also be found in Lulang in 2008 (Zhao et al., 2013) . The relatively high BC from January to March in previous years could be attributed to the regional transport from highly polluted South Asia. However, instead of an increasing trend seen in previous years, BC in our study decreased from December to March (abnormal low in 2020) while slightly increased after June. This was probably due to low BC concentration in South Asia in 2020, providing supporting evidence that Tibet was very sensitive to the impact of COVID-19 and emission reduction in the source region in 2020. The variation of meteorological conditions was different during the LD period in Beijing and Tibet. As reported in several studies, the unfavorable meteorological conditions during the LD period in Beijing could enhance the accumulation of atmospheric pollutants (Su et al., 2020) . found that the relative humidity (RH) increased during the LD period, especially in the two haze episodes with highest RH exceeding 80%, while the wind speed and the highest boundary layer height decreased to around 1.0 m/s and 1 km, respectively. In contrast to Beijing, the variation of meteorological conditions was not significant during the LD period in Tibet. As shown in Figure S8 , the RH in Tibet was the lowest in January and February, and only increased in the end of February due to snow. During the LD period, the highest RH was around 40%, which was only half of that in Beijing. Moreover, the wind speed in Tibet remained rather constant than that in Beijing, with slight fluctuation around 2.0 m/s. Seven sources of BC were identified by the PMF model, including industrial, dust, traffic, biomass burning, coal combustion, firework, and aged BC. It has been reported that BC could be coated with secondary species in winter in Beijing (D. Liu et al., 2015) , thus the aged BC factor in our study was defined with high loadings of BC and SNA. The temporal variation of BC sources in Beijing from November 2019 to March 2020 is shown in Figure 2a . It could be seen that before the LD, traffic and industrial source were predominant. After January 23, the contribution of primary source decreased due to the strict regulation. Instead, the increase of BC concentration in the LD period could be mainly attributed to aged BC, especially during haze episodes. Figures 2b-2d showed the temporal variation of sources in HP1 in the pre-LD period, HP2 and HP3 in the LD period. HP1 was mainly influenced by traffic source with average contribution of 60%. In HP2, which was close to the Spring Festival, traffic source contribution decreased while the influence of firework was important. As for HP3, the contribution of aged BC increased significantly, with the average contribution two times higher than that in HP1. The absolute concentration of aged BC increased from 0.77 ± 0.62 μg/m 3 in HP1 to 1.61 ± 0.59 μg/m 3 in HP3. The high BC concentration in HP3 during the LD period could be influenced by regional transport from neighboring province. As shown in Figures 2e and 2f, HP1 in the pre-LD period was mainly influenced by local emission, while the contribution of regional transport to BC in HP3 during the LD period was more significant and was mainly originated from Hebei province. As reported by , 80% of SNA in Beijing in rapid increasing stage of a haze episode was from emissions in Hebei province. Therefore, the regional transport from the southwestern areas in the LD period could be the reason for the enhanced contribution of aged BC to high BC concentration in HP3. The less impact of primary sources to BC in the LD period could be found in the diurnal variation. As shown in Figure 3a , the diurnal variation of BC in 2016, 2017, and 2019 showed more significant peaks in the morning rush-hour and in the evening, suggesting the impact of traffic source, especially heavy-duty trucks in Beijing. However, the diurnal variation in 2020 was not as significant as previous years, indicating lower contribution of primary traffic source to BC. The diurnal variation in 2018 was minor as well due to the strict policy control and favorable meteorological conditions in that year (Cheng et al., 2019) . The associations of BC with SNA in the LD periods in 5 years are presented in Figure 3b . It could be seen that SNA and BC were best correlated in 2020 (R 2 = 0.77) with the highest slope of 19.6, suggesting that BC was more associated or coated with SNA in 2020. As shown in Figure 2 , aged BC contribution increased during the LD period based on the PMF model. A previous study used a single particle soot photometer (SP2) to analyze the relative coating thickness (Dp/Dc) of aged BC and BC from primary sources by the PMF model (B. Zhang, 2019) . The results showed that aged BC was characterized by larger Dp/Dc compared to primary sources, indicating the aging and mixing of BC with secondary species. LIU ET AL. As shown in Figure 4a , BC in Tibet exhibited a clear reduction from January to March within the LD periods, which was abnormally low compared to previous years. It is important to understand the reason for BC reduction in Tibet based on source apportionment analysis. As only AE33 was able to be deployed continuously for 1 year in Tibet in our study, BC sources were investigated by the aethalometer model. Figures 4b and 4c showed the monthly variation of fossil fuel contribution and the AAE at the YBJ site. It could be seen that although BC concentration decreased abnormally, the relative contribution of non-fossil fuel (mainly biomass burning) and AAE increased during the LD and LD-IN periods. That increasing trend was similar with previous years in Tibet, in which higher biomass burning contribution was found in February and April (Cong et al., 2015; J. Huang et al., 2010) . The increased biomass burning contribution indicated that BC in Tibet during the LD periods was still quite influenced by regional transport, which was similar to previous years (Cong et al., 2015; Kang et al., 2019; Wan et al., 2017; Yang et al., 2021) . It was reported that biomass burning activities were significant in South Asia with the highest intensity during February-May (Bhardwaj et al., 2016) and polluted air mass from biomass burning in South Asia could be transported over the Himalayas to Tibet (Cong et al., 2015) . The regional transport from South Asia influenced the spatial distribution of BC sources in Tibet, with non-fossil fuel contribution higher in Himalayas close to South Asia and lower in Northern Tibet (66% ± 16%) (Kang et al., 2019) . The TPSCF method was applied in this study to analyze the potential source region of BC in the pre-LD and LD period in 2020. As shown in Figure S9 , the potential source region during the LD period was more concentrated in the south of Tibet, including Nepal and Northern India, where biomass was burned extensively (Bhardwaj et al., 2016) . Overall, the impact of regional transport from South Asia was significant during the LD periods in Tibet. However, as the main source region of BC in Tibet, South Asia was also influenced by COVID-19 and its anthropogenic emission significantly decreased, resulting in lower concentration of atmospheric pollutants (Panda et al., 2021) . Figure S9 showed the BC surface mass concentration derived by the NASA GIOVANNI tool (https://giovanni.gsfc.nasa.gov/giovanni/) during the pre-LD and LD periods. It could be seen that the BC concentration level significantly decreased during the LD period in South Asia, with the highest concentration reduced by half. Based on the satellite data, it can be seen that there were fewer fire spots in South Asia during the LD-IN period in 2020 than the same period in 2019 ( Figure S10) . Therefore, the abnormal LIU ET AL. reduction of BC in Tibet during the LD periods could be attributed to high sensitivity of plateau background areas to the emission reduction in South Asia. The results suggested that the investigation and control of cross-border transport was essential to reduce BC concentration in Tibet. In this study, 1-year continuous online measurement of BC was conducted simultaneously in Beijing and Tibet from November 2019 to October 2020. The results showed that the influence of COVID-19 on BC was significantly different in Beijing and Tibet. Due to the two haze episodes in the LD period in Beijing, the average concentration in the LD period was 20% higher than that in pre-LD period. In contrast to Beijing, the average concentration of BC in Tibet decreased over 70% from pre-LD to LD period, which exhibited very different seasonal variation pattern compared to those in previous years. For Beijing, the PMF results showed higher aged BC contribution during haze episodes in the LD period. Moreover, BC and SNA were best correlated in 2020 compared with 2016-2019. The results indicated that although the local emission from primary sources was reduced, higher BC during the LD period could be attributed to the contribution of aged BC and transport from southwestern neighboring areas. LIU ET AL. . Monthly variation of (a) BC concentration, (b) the percentage of fossil fuel contribution, and (c) angstrom exponent (AAE) at the Yangbajain (YBJ) site. Similar to previous years, the contribution from biomass burning source increased during spring (pre-monsoon season) although BC concentration was abnormally low during lockdown (LD) and LD-IN periods in 2020. As for Tibet, although BC concentration decreased, the increased relative contribution of biomass burning in spring indicated impact of regional transport from South Asia during the LD periods in Tibet. Due to COVID-19, the anthropogenic emissions and fire spots reduced significantly in South Asia. Therefore, the decrease of BC in Tibet confirmed high sensitivity of plateau background areas to the anthropogenic emission reduction in South Asia. Overall, this study helps to better understand different sensitivity of BC in megacity and background areas in China to change of anthropogenic activities due to COVID-19. Nevertheless, there is still space for improvement in the study. In the future, the mixing state of BC could be measured simultaneously, which could help better understand aging of BC. In addition, as BC has significant effects on radiative forcing, especially in Tibet, the radiative forcing model could be combined with source apportionment results to further study the impacts on radiative forcing due to anthropogenic emission reduction during COVID-19. Seasonal, interannual, and long-term variabilities in biomass burning activity over South Asia High rise in carbonaceous aerosols under very low anthropogenic emissions over eastern Himalaya, India: Impact of lockdown for COVID-19 outbreak Impact of the COVID-19 pandemic and control measures on air quality and aerosol light absorption in southwestern China Dominant role of emission reduction in PM 2.5 air quality improvement in Beijing during 2013-2017: A model-based decomposition analysis Carbonaceous aerosols on the south edge of the Tibetan Plateau: Concentrations, seasonality and sources The "dual-spot" Aethalometer: An improved measurement of aerosol black carbon with real-time loading compensation Seasonal variations and sources of ambient fossil and biogenic-derived carbonaceous aerosols based on 14 C measurements in Lhasa Enhanced secondary pollution offset reduction of primary emissions during COVID-19 lockdown in China Strong radiative heating due to the mixing state of black carbon in atmospheric aerosols Impact of aerosols on surface ozone during COVID-19 pandemic in southern India: A multi-instrumental approach from ground and satellite observations, and model simulations Linking atmospheric pollution to cryospheric change in the third pole region: Current progress and future prospects Unexpected air pollution with marked emission reductions during the COVID-19 outbreak in China Rapid formation of a severe regional winter haze episode over a mega-city cluster on the north China plain Contrasting physical properties of black carbon in urban Beijing between winter and summer Source apportionment of black carbon during winter in Beijing High-time-resolution source apportionment of PM 2.5 in Beijing with multiple models Source-Receptor relationship revealed by the halted traffic and aggravated haze in Beijing during the COVID-19 lockdown Time-resolved measurement of the ionic fraction of atmospheric fine particulate matter A study on variation of atmospheric pollutants over Bhubaneswar during imposition of nationwide lockdown in India for the COVID-19 pandemic Analysis of air quality data using positive matrix factorization Using aerosol light absorption measurements for the quantitative determination of wood burning and traffic emission contributions to particulate matter Abnormally shallow boundary layer associated with severe air pollution during the COVID-19 lockdown in China Chemical characteristics of ice residual nuclei in anvil cirrus clouds: Evidence for homogeneous and heterogeneous ice formation Organic molecular tracers in the atmospheric aerosols from Lumbini, Nepal, in the northern Indo-Gangetic plain: Influence of biomass burning Variation in concentration and sources of black carbon in a megacity of China during the COVID-19 pandemic Contribution of South Asian biomass burning to black carbon over the Tibetan Plateau and its climatic impact Effects of black carbon with different sources and mixing states on respiratory inflammatory response (Master's thesis) Asian emissions in 2006 for the NASA INTEX-B mission The online measured black carbon aerosol and source orientations in the Nam Co region Aerosol particles at a high-altitude site on the Southeast Tibetan Plateau, China: Implications for pollution transport from South Asia Aethalometer® model AE33 user manual Short-term variability of mineral dust, metals and carbon emission from road dust resuspension Qualitative determination of source regions of aerosol in Canadian high Arctic Identification and estimate of biomass burning contribution to the urban aerosol organic carbon concentrations in Beijing Contribution of wood burning to PM 10 in London. Atmospheric Environment Insights into the chemical characterization and sources of PM 2.5 in Beijing at a 1-h time resolution The identification of source regions of black carbon at a receptor site off the eastern coast of China Identification of the typical metal particles among haze, fog, and clear episodes in the Beijing atmosphere A novel technique for quantifying the regional component of urban aerosol solely from its sawtooth cycles Aerosol absorption over Bay of Bengal during winter: Variability and sources Monitoring and source apportionment of trace elements in PM 2.5 : Implications for local air quality management High-time-resolution source apportionment of PM 2.5 in Beijing with multiple models Source-Receptor relationship revealed by the halted traffic and aggravated haze in Beijing during the COVID-19 lockdown EPA Positive Matrix Factorization (PMF) 5.0 Fundamentals and User Guide Investigation of black and brown carbon multiple-wavelength-dependent light absorption from biomass and fossil fuel combustion source emissions Methods for estimating uncertainty in factor analytic solutions Estimation of the contribution of road traffic emissions to particulate matter concentrations from field measurements: A review Black carbon and wavelength-dependent aerosol absorption in the north China plain based on two-year aethalometer measurements Using aerosol light absorption measurements for the quantitative determination of wood burning and traffic emission contributions to particulate matter Chemical profiles of urban fugitive dust PM 2.5 samples in Northern Chinese cities pH of aerosols in a polluted atmosphere: Source contributions to highly acidic aerosol PM 2.5 pollution in a megacity of southwest China: Source apportionment and implication Influences of fireworks on chemical characteristics of atmospheric fine and coarse particles during Taiwan's Lantern Festival Trace elements in coal: Associations with coal and minerals and their behavior during coal utilization-A review Estimating middle-, neighborhood-, and urban-scale contributions to elemental carbon in Mexico City with a rapid response aethalometer Aerosol optical properties during firework, biomass burning and dust episodes in Beijing Comparison of water-soluble inorganic ions and trace metals in PM 2.5 between online and offline measurements in Beijing during winter Asian emissions in 2006 for the NASA INTEX-B mission The authors declare no conflicts of interest relevant to this study. The associated data can be downloaded online (