key: cord-0738624-gkz2c0z1 authors: Anil, Ismail; Alagha, Omar title: Source Apportionment of Ambient Black Carbon during the COVID-19 Lockdown date: 2020-12-03 journal: Int J Environ Res Public Health DOI: 10.3390/ijerph17239021 sha: 3d4212d257bd9d93d61efece0d3f23616834dc84 doc_id: 738624 cord_uid: gkz2c0z1 Black carbon (BC) particles being emitted from mobile and stationary emission sources as a result of combustion activities have significant impacts on human health and climate change. A lot of social activities have been halted during the COVID-19 lockdowns, which has evidently enhanced the ambient and indoor air quality. This paper investigates the possible emission sources and evaluates the meteorological conditions that may affect the dispersion and transport of BC locally and regionally. Ground-level equivalent BC (eBC) measurements were performed between January 2020 and July 2020 at a university campus located in Dammam city of the Kingdom of Saudi Arabia (KSA). The fossil fuel (eBC(ff)) and biomass burning (eBC(bb)) fractions of total eBC (eBC(t)) concentrations were estimated as 84% and 16%, respectively, during the entire study period. The mean eBC(bb), eBC(ff), and eBC(t) concentrations during the lockdown reduced by 14%, 24%, and 23%, respectively. The results of statistical analyses indicated that local fossil fuel burning emissions and atmospheric conditions apparently affected the observed eBC levels. Long-range potential source locations, including Iraq, Kuwait, Iran, distributed zones in the Arabian Gulf, and United Arab Emirates and regional source areas, such as the Arabian Gulf coastline of the KSA, Bahrain, and Qatar, were associated with moderate to high concentrations observed at the receptor site as a result of cluster analysis and concentration-weighted trajectory analysis methods. The most valuable resource and future asset of any country is its people. Exposing them to low air quality is a considerable risk for their health and welfare [1] [2] [3] [4] . Nowadays, an enormous consumption of non-renewable energy resources by industrial expansion, urbanization, and energy demand has been leading to the deterioration of air quality locally and globally. Black carbon (BC) particles are health-threatening pollutants to human health, particularly in urban cities and industrial areas [1] [2] [3] [5] [6] [7] [8] . It is known that BC particles have both anthropogenic and natural sources, and they are mainly linked to the incomplete combustion of fossil fuels and biomass burning [9] [10] [11] . Fine and coarse particulate matter fractions can suspend in the atmosphere for a long time and be transported over long distances [12, 13] . BC particles act as a host to many chemical and physical species during their lifetime (4-12 days) in the atmosphere [14] . Thus, the SARS-CoV-2 virus, having a diameter of around 0.1 µm, might be carried by BC particles to remote areas [15] . Cluster analysis (CA) of backward air mass trajectories, potential source contribution function (PSCF), and concentration weighted trajectory (CWT) analyses have been widely practiced to investigate the potential source regions of atmospheric pollutants [3, 12, [16] [17] [18] [19] . The PSCF method's constraint is that the grid cells might indicate identical PSCF values when the concentrations of the targeted air pollutant are slightly or metropolitan area has been classified as "arid", since the annual total rainfall is 90 mm. The Air Quality Monitoring Station (AQMS) used in this study is located at the top of the engineering college building on the south campus of Imam Abdulrahman Bin Faisal University (26.3955 N, 50 .1978 E) ( Figure 1 ). An aethalometer for measuring BC mass concentrations is installed in the AQMS, which is 15 m above the sea level and 3 km away from the Arabian Gulf. Local anthropogenic emission spots, such as King Abdulaziz Sea Port (the second busiest and largest seaport of the KSA), industrial area #1, and two main arteries with high traffic density (King Abdulaziz and King Fahd roads), are present within the five-km perimeter of the AQMS. Jubail city, known as the largest industrial city in the world, and Qatif city, hosting the oil production field of Saudi Aramco (the Saudi national oil company) and Saudi Basic Industries Corporation (SABIC) petrochemical and chemical company, are the most significant regional industrial emission sources located in the northwestern part of the study area. Furthermore, Al Hofuf city in the southwestern section of the monitoring station accommodates the largest oil fields in the KSA and can be considered as a regional emission source. In this study, the measurements of BC mass concentrations and meteorological parameters were performed between 8 January and 2 July 2020. The obtained dataset was categorized into three phases: before lockdown (Pre-LD) (8 January-22 March), during lockdown (LD) (23 March-20 June), and after lockdown (Post-LD) (21 June-2 July) to find out any potential effects of COVID-19 lockdown countermeasures on ambient BC emissions. The meteorological parameters, including temperature (°C), humidity (%), pressure (mbar), wind speed (m/s), rainfall (mm), wind direction (0-360°), and solar radiation (W/m 2 ), used during this study were collected via a wireless meteorology station (Davis ® Vantage Pro2+, Hayward, CA, USA). The meteorology station is located next to the AQMS. The meteorology station is mounted 4 m above the roof to ensure that no obstacle impairs the measurements. The wind data statistical analysis revealed that the dominant wind directions during the study period originated from northwesterly sectors (Figure 2 ). The frequencies of the winds blowing from the WNW, NW, and NNW directions through the whole study, Pre-LD, LD, and Post-LD periods were found to be 47.1%, 49 .5%, 43.9%, and 55.6%, respectively. The mean wind speeds for the NW prevailing direction during the Pre-LD, LD, and Post-LD periods were 3.3, 3.7, and 3.3 m/s, respectively. The meteorological parameters, including temperature ( • C), humidity (%), pressure (mbar), wind speed (m/s), rainfall (mm), wind direction (0-360 • ), and solar radiation (W/m 2 ), used during this study were collected via a wireless meteorology station (Davis ® Vantage Pro2+, Hayward, CA, USA). The meteorology station is located next to the AQMS. The meteorology station is mounted 4 m above the roof to ensure that no obstacle impairs the measurements. The wind data statistical analysis revealed that the dominant wind directions during the study period originated from northwesterly sectors ( Figure 2 ). The frequencies of the winds blowing from the WNW, NW, and NNW directions The equivalent BC (eBC) mass concentrations were derived using a rack mount aethalometer (Magee Scientific ® Model AE-31), which is an optical absorption method [71, 72] . The device installed in the AQMS has an 'extended range' inlet capable of measuring seven different wavelengths (370, 470, 520, 590, 660, 880, and 950 nm), as it is recommended for urban areas of moderate to high aerosol concentrations. A weatherproof PM10 inlet (M4121, Magee Scientific ® , Berkeley, CA, USA) was installed on the top of the inlet tube. The PM10 inlet tube line with a length of two meters was covered with a temperature-controlled jacket heater to prevent water condensation. The possible loading effect due to the aethalometer's filter-based optical method needs to be compensated to obtain accurate ambient eBC concentrations [73, 74] . The procedure developed by Virkkula et al. [73] was employed to correct the reported raw eBC concentrations by the instrument using Equation (1): where k is the compensation parameter and ATN is the optical attenuation value recorded by the instrument before and after each eBC reading. The k-parameter may indicate seasonal variations depending on the aerosol properties [73, 74] . In this study, the k-parameter for the nth filter spot was calculated by use of Equation (2) , since the ATN values of the first measurement data of the new filter spots were all higher than zero. During this study, the mean k values were found to be 0.003 ± 0.003 at 370 nm and 0.005 ± 0.004 at 880 nm. where tn+1,first is the time of the first measurement data for the next filter spot and tn,last is the time of the last measurement data for filter spot n. The derivation of fossil fuel BC (BCff) or biomass burning BC (BCbb) fractions is only an indicative datum when it is not supported by additional chemical composition measurements. BCff and BCbb concentrations cannot precisely be calculated from the aethalometer measurements alone [71] . For The equivalent BC (eBC) mass concentrations were derived using a rack mount aethalometer (Magee Scientific ® Model AE-31), which is an optical absorption method [71, 72] . The device installed in the AQMS has an 'extended range' inlet capable of measuring seven different wavelengths (370, 470, 520, 590, 660, 880, and 950 nm), as it is recommended for urban areas of moderate to high aerosol concentrations. A weatherproof PM 10 inlet (M4121, Magee Scientific ® , Berkeley, CA, USA) was installed on the top of the inlet tube. The PM 10 inlet tube line with a length of two meters was covered with a temperature-controlled jacket heater to prevent water condensation. The possible loading effect due to the aethalometer's filter-based optical method needs to be compensated to obtain accurate ambient eBC concentrations [73, 74] . The procedure developed by Virkkula et al. [73] was employed to correct the reported raw eBC concentrations by the instrument using Equation (1): where k is the compensation parameter and ATN is the optical attenuation value recorded by the instrument before and after each eBC reading. The k-parameter may indicate seasonal variations depending on the aerosol properties [73, 74] . In this study, the k-parameter for the nth filter spot was calculated by use of Equation (2) ATN values of the first measurement data of the new filter spots were all higher than zero. During this study, the mean k values were found to be 0.003 ± 0.003 at 370 nm and 0.005 ± 0.004 at 880 nm. k n = eBC raw t n+1, f irst − eBC raw t n,last ATN t n,last × eBC raw t n,last − ATN t n+1, f irst × eBC raw t n+1, f irst (2) where t n+1,first is the time of the first measurement data for the next filter spot and t n,last is the time of the last measurement data for filter spot n. The derivation of fossil fuel BC (BC ff ) or biomass burning BC (BC bb ) fractions is only an indicative datum when it is not supported by additional chemical composition measurements. BC ff and BC bb concentrations cannot precisely be calculated from the aethalometer measurements alone [71] . For this reason, BC ff and BC bb concentrations reported in this study indicate estimated concentrations. The estimation of biomass burning and fossil fuel fractions of the ambient eBC is dependent on the absorption coefficients (b abs ) at the wavelengths of 370 and 880 nm, respectively. Initially, the attenuation coefficient of the analyzed particles (b ATN ) should be calculated. In this study, the calculation of b ATN and b abs values was performed according to the methods developed by Hansen [75] . Equations (3) and (4) were used for the calculation of b ATN and b abs , respectively: where A is the aerosol collecting spot area of filter (1.67 cm 2 ); Q is the sampling flow rate (5 L/min); ∆ATN is the difference between the ATN values recorded by the instrument before and after each eBC reading, and ∆t is the sampling time (5 min). The multiple scattering parameter (C = 2.14) and a linear function of ln(ATN) (R ATN ) in Equation (4) are defined as calibration factors, which are applied to b ATN values to calculate b abs values [37] . The shadowing effect of aerosols due to the filter overload is corrected by the R ATN parameter, which is computed using Equation (5) , where f is the free parameter and is found to be 1.14 and 1.08 for the wavelengths of 370 and 880 nm, respectively. The concentrations of equivalent biomass burning fraction (eBC bb ) and equivalent fossil fuel fraction (eBC ff ) were predicted according to the model developed by Hansen [75] , using the Beer-Lambert Law. The model assumes that b abs is the sum of optical absorption coefficients of biomass burning (b abs,bb ) and fossil fuel (b abs,ff ) fractions, as given in Equations (6) and (7) . b abs (370 nm) = b abs (370 nm) ff + b abs (370 nm) bb (6) b abs (880 nm) = b abs (880 nm) ff + b abs (880 nm) bb The b abs,ff and b abs,bb are computed for the wavelengths of 370 and 880 nm using Equations (8) and (9) , respectively, by assuming that b abs,bb is negligible at 880 nm [68, 76] . The absorption exponents for fossil fuel (α ff ) and for biomass (α bb ) were used as 1 and 2, respectively. The eBC ff and eBC bb concentrations can be calculated by using Equations (11) and (12) , respectively, after obtaining the fossil fuel fraction (FF fraction ) via Equation (10) . eBC ff = FF fraction × eBC (11) eBC bb = eBC − eBC ff (12) Multivariate linear correlation analysis is a useful method to statistically determine the relationships between pollutants or other atmospheric factors influencing the air quality and to reveal the most significant parameters on the concentrations of atmospheric pollutants [12, 43, 77] . In this work, IBM ® SPSS ® statistics software (Ver. 24) was used to perform descriptive statistics, sample comparison, and multivariate linear correlation analyses for the hourly mean values of the measured eBC concentrations and meteorological parameters. Bivariate polar plots of concentrations have recently been used to reveal the collective impact of wind velocity and wind direction on air pollutants' observed concentrations [78] [79] [80] . In this research, bivariate polar plots were drawn to illustrate the impact of wind components on hourly mean values of eBC concentrations through Pre-LD, LD, and Post-LD phases. Cluster analysis is a statistical method to classify trajectories that are close to each other into particular groups. In a cluster analysis, similar trajectories are merged until the total variance is minimized and distinct groups are formed. The TrajStat air mass trajectory plugin (Ver. 1.4.9) of the MeteoInfo GIS (Geographic Information System) meteorology-dispersion model analysis software (Ver. 2.2.7) was employed to calculate the backward air mass trajectories. Furthermore, the previous model was used to perform cluster analysis, cluster statistics, and concentration weighted trajectory analysis [81] . The three-day backward trajectories arriving at 500-m height above ground level (AGL) to the receptor study site were computed for every 6 h throughout the whole studied period using the global data assimilation system (GDAS) meteorological data with half-degree resolution and 10,000 m (AGL) of model top height. Following these, the daily mean eBC concentrations were added to the respective trajectories. In this study, the targeted percentage change criterion of 20% of the total spatial variance was selected to determine the optimum number of clusters [82] . Besides, the hourly time-averaged surface mass BC concentrations (within the domain of 10-45 • N and 10-60 • E and with a spatial resolution of 0.5 • × 0.625 • ) were obtained from the Goddard Earth Observing System Model (GEOS-5) satellite observation database in netCDF data format, which was further processed by using the GIS interface of the MeteoInfo software. The satellite observation datasets were used to identify possible contributions from BC emission points to the clustered trajectories [2, 83, 84] . The CWT method has been widely practiced in investigating the relative importance of potential source areas. In a CWT analysis, each grid cell has a weighted average concentration for each pollutant related to the backward trajectories traveling over that grid cell, as given in Equation (13) [20, 21] . Like the PSCF analysis, a point filter was employed as the last stage of CWT analysis to eliminate grid cells with limited endpoints [20] . As a result of the CWT analysis, weighted concentration zones indicate concentration gradients throughout potential source locations. c l t ijl (13) where C ij is the average weighted concentration in the grid cell (i,j); l is the trajectory index; n is the total number of trajectories; C l is the measured concentration on the arrival of trajectory l; and t ijl is the time spent by trajectory l in the grid cell (i,j). Types of emission sources, emission rates, meteorological conditions, and atmospheric stability are the main factors that significantly impact the hourly variations of ambient BC concentrations [33, 47, 72, 85, 86] . Descriptive statistical parameters, such as mean, median, mode, standard deviation (S.D.), minimum (Min.), maximum (Max.), and percentiles, were calculated for both hourly mean values of meteorological parameters and eBC concentrations ( Table 1 ). The relative differences between the mean and median values of temperature, humidity, wind speed, and pressure parameters during the Pre-LD, LD, and Post-LD periods were found to be less than 6%. The mean values of temperature, humidity, pressure, and solar radiation parameters along the LD changed by 59%, −24%, −1.1%, and 31%, respectively, as compared to the Pre-LD period. The changes in mean wind speed values during the LD and Post-LD periods in comparison with the Pre-LD stage were 3.1% and −12.5%, respectively. The comparison of wind speed and wind direction ( Figure 2 ) suggests that the wind data were quite similar throughout Pre-LD and LD periods, while lower wind speeds during the Post-LD might result in higher eBC concentrations due to the atmospheric stability. During the whole study period, eBC bb , eBC ff , and eBC t concentrations varied between 0.01 and 3.6, 0.12 and 13.4, and 0.25 and 14.2 µg/m 3 , respectively. The eBC concentrations measured within LD phase exhibited the lowest percentile (25th, 50th, and 75th) values. The independent sample comparison t-test was applied to hourly means, medians, and standard deviations of the eBC datasets of Pre-LD, LD, and Post-LD periods. The comparison test results revealed that the measured eBC concentrations (eBC bb , eBC ff , and eBC t ) in each period were statistically different at the 95.0% confidence level. Hourly changes in eBC concentrations throughout the Pre-LD, LD, and Post-LD stages are depicted in Figure 3 . Hourly variations in eBC ff concentrations during Pre-LD, LD, and Post-LD indicated distinct bimodal distributions. Within the Pre-LD period, a significant and sharp peak between 5:00 am and 10:00 am was observed for the mean eBC ff concentrations, which reached the maximum value of 3.2 ± 2.2 µg/m 3 at 8:00 am and then gradually decreased to 1.4 ± 0.9 µg/m 3 by 2:00 pm. Thereafter, eBC ff concentrations demonstrated a minor peak between 2:00 pm and 7:00 pm and a gradual build-up from 7:00 pm to 5:00 am. The highest eBC ff concentrations within the morning peak are attributed to the emissions resulting from intense vehicular activities during rush hours. The decrease in the eBC ff concentration between 10:00 am and 2:00 pm could be explained by the enhanced unstable atmospheric conditions as a consequence of stronger solar radiation and winds. The afternoon peak and increasing eBC ff concentrations until 5:00 am might be ascribed to two main factors: (i) emissions from heavy-duty diesel trucks (HDDTs), which are allowed to enter main arterial roads of the Dammam metropolitan area within the following time intervals: 3:00-5:00 pm and 10:00 pm-5:00 am, and (ii) stable nocturnal atmospheric conditions due to the decreased vertical thermal convection and weakened surface wind speeds. The hourly cycle of eBC ff during the LD phase showed lower concentrations within a narrower distribution in comparison with the Pre-LD and Post-LD stages. The significant and sharp peak between 5:00 am and 10:00 am observed during the Pre-LD period dramatically weakened, which might be ascribed to the imposed lockdown measures. The lowest mean eBC ff concentrations ranging within 1.2 ± 0.7 and 1.1 ± 0.7 µg/m 3 were noticed between 10:00 am and 2:00 pm. The eBC ff concentrations gradually increased starting from 3:00 pm, reached the maximum mean value of 2.1 ± 1.8 µg/m 3 at 3:00 am, and displayed a stepwise decline trend until 10:00 am. Shipping activities and HDDT movements continued within the allowed time intervals to supply essential products to the market and public in compliance with COVID-19 safety measures and regulations forced by the KSA. These emission sources and the shallower nocturnal atmospheric boundary layer were the main reasons for the gradual eBC ff concentration build-up between 3:00 pm and 5:00 am during the LD stage. The daily mean eBC ff concentrations within the Post-LD period were 1.5-and 2.0-times higher than those measured at Pre-LD and LD phases, respectively. The highest eBC ff concentrations varied between 3.1 ± 1.9 and 4.4 ± 2.7 µg/m 3 within a distinct broad peak period from 5:00 pm to 5:00 am. In contrast, the time interval between 8:00 am and 3:00 pm indicated the lowest concentrations in the range of 2.2 ± 1.3 and 2.7 ± 2.3 µg/m 3 . The sharp morning peak observed during the Pre-LD phase was less pronounced within the Post-LD period, which could possibly be explained by the suspension of all schools and educational institutions, which drastically reduced the number of vehicles in rush-hour traffic. Boosted traffic activities from the supply chain and logistic disruptions because of the COVID-19 pandemic resulted in a surge in fossil fuel consumption. The lower mean wind speed and long-range transport might be linked to the increased eBC ff concentrations during the Post-LD period compared to the Pre-LD period. The afternoon peak and increasing eBCff concentrations until 5:00 am might be ascribed to two main factors: (i) emissions from heavy-duty diesel trucks (HDDTs), which are allowed to enter main arterial roads of the Dammam metropolitan area within the following time intervals: 3:00-5:00 pm and 10:00 pm-5:00 am, and (ii) stable nocturnal atmospheric conditions due to the decreased vertical thermal convection and weakened surface wind speeds. The hourly cycle of eBCff during the LD phase showed lower concentrations within a narrower distribution in comparison with the Pre-LD and Post-LD stages. The significant and sharp peak between 5:00 am and 10:00 am observed during the Pre-LD period dramatically weakened, which Monthly variations in eBC concentrations and their changes during Pre-LD, LD, and Post-LD phases are illustrated in Figure 4 . Monthly changes in eBC ff concentrations are more apparent than those of eBC bb concentrations. The eBC ff concentrations observed in January slightly changed in February by increasing from 2.1 ± 1.1 to 2.3 ± 0.9 µg/m 3 . Afterward, the mean eBC ff concentration sharply decreased in March by 29.4%, which might presumably be attributed to the imposed measures to combat the COVID-19 pandemic. Compared to February, the change in eBC ff concentrations was most evident in April as a 24-h LD was set nationwide between 6 April and 28 May. The lowest mean eBC ff and eBC bb concentrations were recorded as 1.3 ± 0.6 and 0.26 ± 0.07 µg/m 3 , respectively, in April. During the entire LD period, the public's movement was restricted, except when accessing essential nearby facilities, such as food stores, pharmacies, and hospitals. The previous activities were allowable through obtaining an emergency permit between 6:00 am and 3:00 pm using the official Tawakkalna mobile application. Despite the full LD stringent measures, both eBC ff and eBC bb concentrations observed during May increased by 50.5 and 15.4%, respectively, compared to those measured in April. Besides, the mean eBC t concentration in May was 11% lower than that recorded in February. The mean eBC ff and eBC bb concentrations during June were 2.1 ± 1.2 and 0.29 ± 0.12 µg/m 3 , respectively, and reverted the levels observed in January as the 24-h LD was partially lifted on 28 May. Domestic heating, burning of crop residues, and forest fires have been reported as primary sources of BCbb emissions, contributing significantly to the BC budget in other regions, such as Brazil [47] , Canada [46] , China [87] , Finland [88] , Greece [89] , India [90] , Italy [91] , Spain [77] , Turkey [41] , and the USA [92] . However, wood-burning stoves in restaurants and open-air barbecues are the only eBCbb sources within the study area and most parts of the KSA, since domestic heating is provided by air conditioners/electric heaters and farms/forested lands are very scarce. The study area's biomass-burning characteristic mentioned herein is the main reason why the ratios of eBCbb concentrations to eBCt concentrations (BBfraction) were estimated as low as 0.16 ± 0.12, 0.17 ± 0.09, and 0.13 ± 0.06 throughout the Pre-LD, LD, and Post-LD stages, respectively. These results suggested that fossil fuel emissions dominated the ambient eBCt budget over the study area by representing 84% of the eBCt concentrations during the entire study period. The mean eBCbb and eBCff concentrations during the LD stage dropped by 14.2% and 24.2% and then increased by 31.7% and 51.4%, respectively, after the release of LD measures in reference to the Pre-LD phase. This unexpected increase, especially in eBCff concentrations, after lifting of the LD was most probably derived from (i) exertions for restoring supply-chain losses, such as elevated fossil fuel consumption in transportation and industrial activities, (ii) escalated construction activities in the vicinity of the study area, (iii) more stable atmospheric conditions due to lower wind speed compared to the Post-LD phase, and (iv) long-range transport of BC particles from other regions to the receptor site. Domestic heating, burning of crop residues, and forest fires have been reported as primary sources of BC bb emissions, contributing significantly to the BC budget in other regions, such as Brazil [47] , Canada [46] , China [87] , Finland [88] , Greece [89] , India [90] , Italy [91] , Spain [77] , Turkey [41] , and the USA [92] . However, wood-burning stoves in restaurants and open-air barbecues are the only eBC bb sources within the study area and most parts of the KSA, since domestic heating is provided by air conditioners/electric heaters and farms/forested lands are very scarce. The study area's biomass-burning characteristic mentioned herein is the main reason why the ratios of eBC bb concentrations to eBC t concentrations (BB fraction ) were estimated as low as 0.16 ± 0.12, 0.17 ± 0.09, and 0.13 ± 0.06 throughout the Pre-LD, LD, and Post-LD stages, respectively. These results suggested that fossil fuel emissions dominated the ambient eBC t budget over the study area by representing 84% of the eBC t concentrations during the entire study period. The mean eBC bb and eBC ff concentrations during the LD stage dropped by 14.2% and 24.2% and then increased by 31.7% and 51.4%, respectively, after the release of LD measures in reference to the Pre-LD phase. This unexpected increase, especially in eBC ff concentrations, after lifting of the LD was most probably derived from (i) exertions for restoring supply-chain losses, such as elevated fossil fuel consumption in transportation and industrial activities, (ii) escalated construction activities in the vicinity of the study area, (iii) more stable atmospheric conditions due to lower wind speed compared to the Post-LD phase, and (iv) long-range transport of BC particles from other regions to the receptor site. A bivariate correlation analysis was applied to hourly mean values of both meteorological parameters and eBC fractions for the whole study period in order to reveal the effects of local meteorology on ambient eBC ff , eBC bb , and eBC t concentrations. The Pearson correlation coefficients between each parameter are shown in Table 2 . Moderate negative correlations were found between wind speed and eBC ff (r = −0.445), eBC bb (−0.375), and eBC t (−0.435) concentrations, suggesting that BC particles were more dispersed under higher wind speeds or that they accumulated at lower wind speeds, as also reported by other studies [33, 88, 93] . The effect of wind direction on eBC t concentrations was further evaluated using bivariate polar plots because of the polar characteristic of wind direction data (between 0 and 360 • ). Other meteorological parameters did not yield any statistically significant correlations with the measured eBC concentrations. eBC t concentrations resulted in strong positive correlations with eBC ff (0.996) and eBC bb concentrations (0.634), implying that the eBC t composition was mostly represented by the eBC ff fraction during the study period. Bivariate polar plots of concentrations were graphed to find out the conjunct effect of wind velocity and wind direction on hourly mean eBC t concentrations and potential source locations of eBC t emissions affecting the receptor site ( Figure 5 ). It can be inferred from Figure 5a that low wind speeds (<2 m/s) were mostly accompanied by high eBC t concentrations (>4 µg/m 3 ), revealing that local sources mainly dominated the observed high concentrations during the entire study period. On the other hand, moderate (2-4 µg/m 3 ) and high (>4 µg/m 3 ) concentrations were also brought on by high-speed winds blowing from varying directions to the monitoring site, implying the potential contributions of regional and long-range sources to the measured concentrations. The Pearson correlation coefficients between the wind speed and eBC t concentrations were found to be −0.674, −0.271, and −0.515 for the Pre-LD, LD, and Post-LD stages, respectively. These correlation levels could suggest that local sources' effects were more evident in the measured concentrations at the Pre-LD and Post-LD stages. The low-speed winds (<4 m/s) blowing from westerly sectors most likely transported the fossil-fuel-combustion plumes originating from nearby emission spots, such as the King Abdulaziz Sea Port, King Abdulaziz and King Fahd roads, and industrial area #1, which resulted in moderate-to-high concentrations at the monitoring site throughout the Pre-LD, LD, and Post-LD phases, as seen in Figure 5b -d, respectively. During the Pre-LD period, the easterly sectors demonstrated moderate concentrations at variable wind speeds. Contrary to the Pre-LD stage, the effect of westerly winds on eBC t concentrations was less pronounced, which could possibly be due to dramatically reduced traffic volume during the LD. The highest concentrations at low wind speeds and moderate concentrations at variable wind speeds were, instead, associated with the LD's easterly directions. The emissions from cargo ships docking to and hailing from the King Abdulaziz Sea Port, land reclamation activities, and transboundary pollution transport might be the main reasons giving rise to the elevated eBC t concentrations at the receptor site since there are no other known local emission sources on the east section of the monitoring site. Throughout the Post-LD phase, eBC t concentrations showed a distinct distribution on the bivariate polar plot, depicted in Figure 5d . It was observed that intense land reclamation construction activities had started at the sector between NE and SE directions with respect to the monitoring location during the Post-LD period. Diesel-fuel-combustion plumes of heavy-duty excavators and trucks operating on the nearby construction site were carried to the monitoring site by low-speed winds (<3 m/s) and, thus, led to the highest concentrations (>8 µg/m 3 ) in the presence of atmospheric inversion conditions. The low-speed winds (<4 m/s) blowing from westerly sectors most likely transported the fossilfuel-combustion plumes originating from nearby emission spots, such as the King Abdulaziz Sea Port, King Abdulaziz and King Fahd roads, and industrial area #1, which resulted in moderate-tohigh concentrations at the monitoring site throughout the Pre-LD, LD, and Post-LD phases, as seen in Figure 5b -d, respectively. During the Pre-LD period, the easterly sectors demonstrated moderate concentrations at variable wind speeds. Contrary to the Pre-LD stage, the effect of westerly winds on eBCt concentrations was less pronounced, which could possibly be due to dramatically reduced traffic volume during the LD. The highest concentrations at low wind speeds and moderate The cluster analysis model grouped backward air mass trajectories to the whole study, Pre-LD, LD, and Post-LD periods into five distinct clusters as an optimum number of clusters according to the objective percentage change criterion applied for the total spatial variance of each cluster solution. The frequency of backward trajectories of each cluster in all clusters and the mean value of measured eBC t concentrations assigned to backward trajectories in each cluster are given in Table 3 . The mean trajectories of clusters were superimposed on the satellite-based BC surface mass concentration data to visualize the BC emissions spots that might enrich the air masses traveling over them and eventually reaching the receptor site ( Figure 6 ). to the monitoring site by low-speed winds (<3 m/s) and, thus, led to the highest concentrations (>8 μg/m 3 ) in the presence of atmospheric inversion conditions. The cluster analysis model grouped backward air mass trajectories to the whole study, Pre-LD, LD, and Post-LD periods into five distinct clusters as an optimum number of clusters according to the objective percentage change criterion applied for the total spatial variance of each cluster solution. The frequency of backward trajectories of each cluster in all clusters and the mean value of measured eBCt concentrations assigned to backward trajectories in each cluster are given in Table 3 . The mean trajectories of clusters were superimposed on the satellite-based BC surface mass concentration data to visualize the BC emissions spots that might enrich the air masses traveling over them and eventually reaching the receptor site ( Figure 6 ). Figure 6a demonstrates that air mass trajectories traversing the NNE zone of Egypt, the United Arab Emirates, and central Iraq and the SE section of Iraq with borders with Kuwait and Iran were presumably highly enriched with BC since these areas yielded the highest BC emissions according to the satellite observations. The mean trajectories of clusters reaching the study area from the NNW section have different origins, but they followed a familiar route during the last 12 h of their travel time. The zone, including the SE region of Iraq and Kuwait and the NE section of the Arabian Gulf, commonly passed over by the corresponding air masses, has a significant potential to contribute to the elevated eBC concentrations measured at the monitoring site. Anil et al. [82] also classified this zone as a hot spot emitting high amounts of SO 2 , NO 2 , and PM 10 pollutants in their recent study investigating transport pathways on the chemical composition of wet deposition over the Dammam area. Satellite observations show higher mean BC concentration distributions over the mapped domain for the Pre-LD period as compared with the LD and Post-LD periods, except over the Mediterranean Sea, Red Sea, and most parts of Africa (Figure 6b-d) . During the Pre-LD period, clusters 1 and 2 classifying 56% of all air mass trajectories represented the highest mean eBC concentrations of 2.7 and 2.8 µg/m 3 , respectively, among other clusters. Clusters 1 and 2, including slow-moving air masses originated from central Iraq and the NE part of the KSA, were enhanced with long-range and regional BC emissions and arrived at the receptor site from NW and S directions, respectively. The trajectories in cluster 2 could be ascribed to regionally polluted air masses by fossil fuel burning emissions as they crossed over heavily industrialized areas of Jubail, Ras Tanura, and Qatif cities, leading to high eBC concentrations in the study area. The lower mean BC concentration distributions had been experienced over most of the countries within the studied domain, which might be explained by the LD countermeasures imposed in those countries (Figure 6c ). The air mass trajectories with the highest mean eBC concentrations of 2.3 and 2.0 µg/m 3 were associated with clusters 1 and 2, describing 31.1% and 27.8% of all trajectories, respectively, throughout the LD phase. Cluster 1, deriving from the Red Sea coast of Egypt, traveled over the northern part of the KSA, enriched with BC emissions over the SE region of Iraq, Kuwait, and the eastern coast of the KSA, and contributed to the receptor's BC budget. Cluster 2, representing slow-moving and low-altitude (850 m AGL) air masses, was more effectively polluted by the industrial and urban emissions along the pathway through the United Arab Emirates, Buqayq, and Al Hofuf cities of the KSA, and eventually arrived at the study area from the SSW section. Satellite observations indicated slightly higher mean BC concentration dispersions during the LD period compared to the Pre-LD period (Figure 6d ). Cluster 1, accounting for 50% of the back trajectories, originated from Lebanon, traversed Syria, Iraq, and Kuwait, and eventually brought moderately enriched air masses to the receptor site, resulting in the measured mean eBC concentration of 3.2 µg/m 3 . Clusters 2 and 5, representing the lowest number of air mass trajectories, induced the highest observed mean concentrations of 6.2 and 7.0 µg/m 3 , respectively. These extreme concentrations could be attributed to slow and low-altitude (820 m AGL) flowing characteristics of the clusters that are intensely loaded with anthropogenic emissions arising from Iran's western border, the eastern region of Iraq, Kuwait, United Arab Emirates, and Bahrain during their residence times. The CWT analysis was further applied in order to quantitatively determine the relative contributions of regional and long-range source regions to the ambient eBC t concentrations in the study area, since the capability of the cluster analysis is quite limited in delineating the effects of potential source regions at the receptor site in detail. The CWT analysis results of the entire study, Pre-LD, LD, and Post-LD periods are visualized in Figure 7 , where the weighted concentration (WCWT, µg/m 3 ) calculated for each cell (0.5 • × 0.5 • ) demonstrates the contribution level of a potential source location at the receptor site. During the entire study period, long-range source locations with WCWT values > 1 μg/m 3 , including the central and southeastern section of Iraq, Iran's western border with Iraq, Kuwait, distributed zones in the Arabian Gulf, and the United Arab Emirates, were associated with moderate to high eBC concentrations observed at the receptor site (Figure 7a ). The impacts of long-range source zones mentioned herein were also recently reported for the BC levels over Riyadh and Dammam's rainwater chemical compositions [52, 83] . In addition, it is evident from Figure 7a that the effects of regional source locations, such as the Arabian Gulf coastline of the KSA, Bahrain, and Qatar, with moderate to high WCWT values on the observed eBC levels were substantial. These medium-and long-range potential source zones were, thus, coherent with the idea that the advections of emissions from fossil fuel combustion at factories and power plants, gas flaring at oil refineries, transportation activities at densely populated urban cores, and heavy maritime traffic of oil tankers and cargo ships on the Arabian Gulf contributed to the BC budget over Dammam's atmosphere. The advection of air masses from long-range potential source zones, such as Iraq, Iran, Kuwait, and the United Arab Emirates, with WCWT values between 1 and 2 μg/m 3 , moderately contributed to the observed eBC concentrations in the study area through the Pre-LD phase (Figure 7b ). The highest concentrations measured during the Pre-LD period were accompanied by regional source emissions from heavy industrial activities in Jubail city, oil production fields and petrochemical industries in Qatif city, and intense transportation activities on the Arabian Gulf, further adding to high local BC emissions within the Dammam area. The WCWT values remarkably decreased during the LD period. This decrease could be possibly due to the reduced traffic and industrial emissions in response to the applied LD measures by many countries (Figure 7c ). Within this period, contributions from regional source locations, including Qatar, Bahrain, and Jubail city, to the observed eBC levels were significant. The WCWT values indicated a relatively narrow distribution in the post-LD stage During the entire study period, long-range source locations with WCWT values > 1 µg/m 3 , including the central and southeastern section of Iraq, Iran's western border with Iraq, Kuwait, distributed zones in the Arabian Gulf, and the United Arab Emirates, were associated with moderate to high eBC concentrations observed at the receptor site (Figure 7a ). The impacts of long-range source zones mentioned herein were also recently reported for the BC levels over Riyadh and Dammam's rainwater chemical compositions [52, 83] . In addition, it is evident from Figure 7a that the effects of regional source locations, such as the Arabian Gulf coastline of the KSA, Bahrain, and Qatar, with moderate to high WCWT values on the observed eBC levels were substantial. These medium-and long-range potential source zones were, thus, coherent with the idea that the advections of emissions from fossil fuel combustion at factories and power plants, gas flaring at oil refineries, transportation activities at densely populated urban cores, and heavy maritime traffic of oil tankers and cargo ships on the Arabian Gulf contributed to the BC budget over Dammam's atmosphere. The advection of air masses from long-range potential source zones, such as Iraq, Iran, Kuwait, and the United Arab Emirates, with WCWT values between 1 and 2 µg/m 3 , moderately contributed to the observed eBC concentrations in the study area through the Pre-LD phase (Figure 7b ). The highest concentrations measured during the Pre-LD period were accompanied by regional source emissions from heavy industrial activities in Jubail city, oil production fields and petrochemical industries in Qatif city, and intense transportation activities on the Arabian Gulf, further adding to high local BC emissions within the Dammam area. The WCWT values remarkably decreased during the LD period. This decrease could be possibly due to the reduced traffic and industrial emissions in response to the applied LD measures by many countries (Figure 7c ). Within this period, contributions from regional source locations, including Qatar, Bahrain, and Jubail city, to the observed eBC levels were significant. The WCWT values indicated a relatively narrow distribution in the post-LD stage (Figure 7d ). The high eBC concentrations measured at the monitoring station were mainly of local origin. They could mostly be attributed to the industrial emissions from Ras Tanura city, located 40 km to the NNW direction, and local emissions in the vicinity of the study area, as discussed in Section "3.3. Effect of meteorology and local sources". Worldwide ground-level BC concentrations were compared with the BC data observed in the study area before and during the LD, as given in Table 4 . Before the LD period, the highest BC levels were reported in populated urban cores, such as Suzhou and Chongqing of China and Milan of Italy, where local emissions from traffic, industrial, household heating, and cooking activities have been significantly contributing the ambient BC budget over those regions. The mean BC concentration measured in Dammam was higher than those in Bhubaneswar/India and Massachusetts/USA, most likely due to more extensive industrial and transportation activities in the study area. The BC levels indicated decreasing trends with varying rates for all compared sites in line with COVID-19 LD measures. The potential effect of LD countermeasures on BC emissions was more pronounced in Milan/Italy, where the highest decrease ratio of 71% was noticed. The mean BC reduction rate observed in Dammam is one of the lowest among others. This result could be attributed to a wide range of emission sources with varying speeds, intrinsic atmospheric stability, specific arid climate conditions of the Dammam area, and long-range sources. This research presents the first attempt to report ground-level measurements of eBC mass concentrations for the Eastern Province of Saudi Arabia, particularly for the Dammam metropolitan area. Continuous monitoring of hourly eBC concentrations was initiated at a university campus' premises between 8 January and 2 July 2020. The eBC ff and eBC bb fractions in total eBC concentrations were estimated as 84% and 16%, respectively, by using the aethalometer data during the entire study period. Hourly variations of eBC ff and eBC bb levels during Pre-LD, LD, and Post-LD periods indicated distinct bimodal distributions. Before the implementation of COVID-19 LD measures, the highest eBC concentrations within the significant and sharp morning peak (5:00-10:00 am), the lowest concentrations between 10:00 am and 2:00 pm, and the minor peak (2:00 pm and 7:00 pm) along with the gradual concentration build-up from 7:00 pm to 5:00 am were attributed to (i) intense vehicular activities during morning rush hours, (ii) enhanced atmospheric instability around midday, and (iii) emissions from HDDTs under stable nocturnal atmospheric conditions, respectively. The mean eBC bb , eBC ff , and eBC t concentrations during the LD decreased by 14%, 24%, and 23%, respectively. The reductions in eBC levels within the LD period could be possibly attributed to the strict restrictions on public movement and activities for preventing the spread of COVID-19, further adding to changes in atmospheric conditions and effects of long-range transport. The results of the bivariate correlation analysis and polar plots of concentrations revealed obvious influences of local fossil fuel burning emissions and atmospheric conditions on the observed eBC levels. Besides, distinct regional and long-range source regions contributing to the measured concentrations at the receptor site were identified for the Pre-LD, LD, and Post-LD periods by using a cluster analysis of backward air mass trajectories and CWT analysis. Long-range potential source locations, including Iraq, Kuwait, Iran, distributed zones in the Arabian Gulf, and the United Arab Emirates, and regional source areas, such as the Arabian Gulf coastline of KSA, Bahrain, and Qatar, were associated with moderate to high concentrations observed at the receptor site. The nationwide COVID-19 LD has procured a unique chance for the atmospheric researchers, policymakers, and administrative bodies to assess the effects of emission reductions on the ambient air quality and to rethink the current and future air pollution mitigation strategies. Author Contributions: All authors contributed to the conceptualization and design of the study. Experimental preparations, data collection, and analysis were performed by I.A. The first draft of the manuscript was written by O.A., and all authors commented on the manuscript. Project administration and funding acquisition were provided by O.A. All authors have read and agreed to the published version of the manuscript. Distribution and health risks of aerosol black carbon in a representative city of the Qinghai-Tibet Plateau Assessing public health burden associated with exposure to ambient black carbon in the United States Health impact of PM10, PM2.5 and black carbon exposure due to different source sectors in Stockholm, Gothenburg and Umea, Sweden Particulate matter (PM2.5, PM10-2.5, and PM 10) and children's hospital admissions for asthma and respiratory diseases: A bidirectional case-crossover study Respiratory health effects from exposure to carbon black: Results of the phase 2 and 3 cross sectional studies in the European carbon black manufacturing industry Assessing outdoor air quality and public health impact attributable to residential black carbon emissions in rural China Respiratory health effects of carbon black: A survey of European carbon black workers Ambient black carbon particles reach the fetal side of human placenta Effects of Personal Exposures to Micro-and Nano-Particulate Matter, Black Carbon, Particle-Bound Polycyclic Aromatic Hydrocarbons, and Carbon Monoxide on Heart Rate Variability in a Panel of Healthy Older Subjects Characteristics of PM2.5 and Black Carbon Exposure Among Subway Workers Association between Mortality and Short-Term Exposure to Particles, Ozone and Nitrogen Dioxide in Long-range potential source contributions of episodic aerosol events to PM10 profile of a megacity ATR-FTIR Spectroscopic Study of Functional Groups in Aerosols: The Contribution of a Saharan Dust Transport to Urban Atmosphere in Istanbul The atmospheric lifetime of black carbon Is there an association between the level of ambient air pollution and COVID-19? Ambient black carbon particulate matter in the coal region of Dhanbad Chemical characteristics of trace metals in PM10 and their concentrated weighted trajectory analysis at Central Delhi, India Characteristics of Black Carbon aerosols over Patiala Northwestern part of the IGP: Source apportionment using cluster and CWT analysis Transport pathways of PM10 during the spring in northwest China and its characteristics of potential dust sources Comparison of hybrid receptor models to locate PCB sources in Chicago Transport of polluted boundary layer air from the PO valley to high-alpine sites Spatial and temporal assessment of organic and black carbon at four sites in the interior of South Africa Real time black carbon measurements in West and Central Africa urban sites The anthropogenic contribution to atmospheric black carbon concentrations in southern Africa: A WRF-Chem modeling study The vertical variability of black carbon observed in the atmospheric boundary layer during DACCIWA Local Emissions and Regional Wildfires Influence Refractory Black Carbon Observations Near Palmer Station Changes in black carbon deposition to Antarctica from two high-resolution ice core records Holocene black carbon in Antarctica paralleled Southern Hemisphere climate Black carbon aerosols over coastal Antarctica and its scavenging by snow during the Southern Hemispheric summer Black Carbon Aerosol at McMurdo Station High altitude (∼4520 m amsl) measurements of black carbon aerosols over western trans-Himalayas: Seasonal heterogeneity and source apportionment Black carbon and carbon monoxide over Bay of Bengal during W_ICARB: Source characteristics Diurnal and seasonal variations of black carbon and PM2.5 over New Delhi, India: Influence of meteorology The characteristics of brown carbon aerosol during winter in Beijing Spatial-Temporal Pattern of Black Carbon (BC) Emission from Biomass Burning and Anthropogenic Sources in New South Wales and the Greater Metropolitan Region of Sydney Estimates of greenhouse gas and black carbon emissions from a major Australian wildfire with high spatiotemporal resolution Absorption of light by soot particles: Determination of the absorption coefficient by means of aethalometers Estimates of Black Carbon and Size-Resolved Particle Number Emission Factors from Residential Wood Burning Based on Ambient Monitoring and Model Simulations Vertical profiles of aerosol absorption coefficient from micro-Aethalometer data and Mie calculation over Milan A European aerosol phenomenology-5: Climatology of black carbon optical properties at 9 regional background sites across Europe Black carbon and size-segregated elemental carbon, organic carbon compositions in a megacity: A case study for Istanbul The relationship between optical attenuation and black carbon concentration for ambient and source particles A trajectory analysis of atmospheric transport of black carbon aerosols to Canadian high Arctic in winter and spring Long term trends in Black Carbon Concentrations in the Northeastern United States Contribution of particulate brown carbon to light absorption in the rural and urban Southeast US Ambient measurements and source apportionment of fossil fuel and biomass burning black carbon in Ontario Local and remote black carbon sources in the Metropolitan Area of Buenos Aires Relationship between black carbon (BC) and heavy traffic in São Paulo Exposure to fine particulate, black carbon, and particle number concentration in transportation microenvironments Aerosols physical properties at Hada Al Sham, western Saudi Arabia Ambient air quality in the holy city of Makkah: A source apportionment with elemental enrichment factors (EFs) and factor analysis (PMF) Sources of PM2.5 carbonaceous aerosol in Riyadh, Saudi Arabia WHO Coronavirus Disease (COVID-19). Available online Changes in Air Quality during the First-Level Response to the Covid-19 Pandemic in Shanghai Municipality Impact on Air Quality of the COVID-19 Lockdown in the Urban Area of Palermo (Italy) Impact of COVID-19 Induced Lockdown on Environmental Quality in Four Indian Megacities Using Landsat 8 OLI and TIRS-Derived Data and Mamdani Fuzzy Logic Modelling Approach Observed and Potential Impacts of the COVID-19 Pandemic on the Environment The impact of COVID-19 lockdown on the air quality of Eastern Province, Saudi Arabia. Air Qual. Atmos. Heal. 2020, 78 Impact of COVID-19 Event on the Air Quality in Iran Morales Suárez-Varela, M.M. The Effect of the Covid-19 Lockdown on Air Quality in Three Italian Medium-Sized Cities. Atmosphere 2020 Impact of the COVID-19 Outbreak on Air Quality in Korea The Effects of COVID-19 Measures on Air Pollutant Concentrations at Urban and Traffic Sites in Istanbul Effect of Road Traffic on Air Pollution. Experimental Evidence from COVID-19 Lockdown situ Measurements and City Scale Modelling to Assess the COVID-19 Lockdown Effects on Emissions and Air Quality in Assessing air quality changes in large cities during COVID-19 lockdowns: The impacts of traffic-free urban conditions in Almaty Reductions in traffic-related black carbon and ultrafine particle number concentrations in an urban neighborhood during the COVID-19 pandemic Lockdown for CoViD-2019 in Milan: What are the effects on air quality? Impact of the COVID-19 pandemic and control measures on air quality and aerosol light absorption in Southwestern China Characterization of the aerosol chemical composition during the COVID-19 lockdown period in Suzhou in the Yangtze River Delta A study on variation of atmospheric pollutants over Bhubaneswar during imposition of nationwide lockdown in India for the COVID-19 pandemic Recommendations for reporting black carbon measurements Properties of particulate pollution in the port city of A simple procedure for correcting loading effects of aethalometer data The "dual-spot" Aethalometer: An improved measurement of aerosol black carbon with real-time loading compensation Scientific Aethalometer Handbook Source apportionment of black carbon during winter in Beijing Characterization of atmospheric black carbon and co-pollutants in urban and rural areas of Spain Analysis of the variability of airborne particulate matter with prevailing meteorological conditions across a semi-urban environment using a network of low-cost air quality sensors The Interaction between Air Quality and Meteorological Factors in an Arid Environment of Makkah, Saudi Arabia The impact of movement control order (MCO) during pandemic COVID-19 on local air quality in an urban area of Klang valley GIS-based software that uses various trajectory statistical analysis methods to identify potential sources from long-term air pollution measurement data Effects of transport patterns on chemical composition of sequential rain samples: Trajectory clustering and principal component analysis approach Source Identification of Episodic Rain Pollutants by New Approach: Combining Satellite Observations and Backward Air Mass Trajectories A novel approach for determining source-receptor relationships in model simulations: A case study of black carbon transport in northern hemisphere winter Evaluation of spatial and temporal heterogeneity of black carbon aerosol mass concentration over India using three year measurements from IMD BC observation network Continuous measurement of black carbon aerosol in urban Nanjing of Yangtze River Delta Characterization of black carbon in an urban-rural fringe area of Beijing Characteristics and source apportionment of black carbon in the Helsinki metropolitan area Long-term brown carbon spectral characteristics in a Mediterranean city (Athens) -A study during agriculture crop residue burning period using ground measurements and satellite data Source apportionment of black carbon (BC) from fossil fuel and biomass burning in metropolitan Light absorption of organic carbon and its sources at a southeastern U.S. location in summer Measurement and analysis of black carbon aerosols over a tropical semi-arid station in Kadapa The authors gratefully thank Yaqiang Wang at the Chinese Academy of Meteorological Sciences for the provision of the MeteoInfo desktop application (Ver 2.2.7) (http://meteothink.org/) used in this publication. The authors acknowledge the NOAA Air Resources Laboratory (ARL) for the access to Global Data Assimilation System (GDAS) meteorological data through the READY website (https://www.ready.noaa.gov). Satellite-based black carbon surface mass concentrations used in this study were produced with the Giovanni online data system, developed and maintained by the NASA GES DISC (https://giovanni.gsfc.nasa.gov). The authors declare no conflict of interest.