Characterization and Source Apportionment of PM in Handan—A Case Study during the COVID-19

: Handan is a typical city affected by regional particulate pollution. In order to investigate particulate matter (PM) characterization, source contributions and health risks for the general populations, we collected PM samples at two sites affected by a pollution event (12–18 May 2020) during the COVID-19 pandemic and analyzed the major components (SNA, OCEC, WSIIs, and metal elements). A PCA-MLR model was used for source apportionment. The carcinogenic and non-carcinogenic risks caused by metal elements in the PM were assessed. The results show that the renewal of old neighborhoods signiﬁcantly inﬂuences local PM, and primarily the PM 10 ; the average contribution to PM 10 was 27 µ g/m 3 . The source apportionment has indicated that all other elements came from dust, except Cd, Pb and Zn, and the contribution of the dust source to PM was 60.4%. As PM 2.5 grew to PM 10 , the PM changed from basic to acidic, resulting in a lower NH 4+ concentration in PM 10 than PM 2.5 . The carcinogenic risk of PM 10 was more than 1 × 10 − 6 for both children and adults, and the excess mortality caused by the renewal of the community increased by 23%. Authorities should pay more attention to the impact of renewal on air quality. The backward trajectory and PSCF calculations show that both local sources and short-distance transport contribute to PM—local sources for PM 10 , and short-distance transport in southern Hebei, northern Henan and northern Anhui for PM 2.5 , SO 2 and NO 2 .


Introduction
In recent years, combined particulate matter (PM) and ozone pollution has been detected in the air in Beijing, Tianjin, Hebei, and the surrounding areas [1].The PM concentration directly affects public physical and psychological health, and the annual economic loss associated with PM 2.5 (aerodynamic diameter of less than 2.5 µm) and PM 10 (aerodynamic diameter of less than 10 µm) health hazards in the Beijing-Tianjin-Hebei (BTH) region is 122.4 and CNY 118.34 billion, respectively [2].The main components of atmospheric PM include organic carbon (OC), elemental carbon (EC), sulfate (SO 4 2− ), nitrate (NO 3 − ), ammonium (NH 4 + ), and metal elements [3], which are hazardous to humans [4].PM 2.5 and its components, such as polycyclic aromatic hydrocarbon and hexavalent chromium, are carcinogens.Long-term exposure can increase the carcinogenic risk (CR) [5,6].Particulate aerosols not only affect public health, but also significantly reduce the downward shortwave flux and boundary layer height due to their radiation effect, and even affect surface temperature and relative humidity [7].In addition to national or provincial monitoring stations, sensors are being used more and more frequently for air pollutant monitoring [8].These sensors are favored by researchers because of their convenience and low cost.The use of sensors leads to higher spatial precision in monitoring, but the methods of component detection of particulates are still dominated by on-line aerosol mass spectrometry or off-line sampling plus laboratory testing [9][10][11][12].Satellite remote sensing has also been used in air quality monitoring in recent years, as it can provide long-term observations with the advantages of wide spatial coverage and multi-element synchronous acquisition.For example, moderate-resolution imaging spectral radiometry (MODIS) provides aerosol optical thickness (AOT) data, which has been widely used in the environment and other fields [13].In addition, microwave limb sounders (MLSs) [14], infrared atmospheric sounding interferometers (IASIs) [15] and total ozone mapping spectrometers (TOMSs) are being more commonly applied [16].Satellite remote sensing is mainly used for monitoring gas pollutants such as ozone, sulfur dioxide and formaldehyde, and it cannot achieve near-realtime monitoring because of its high altitude.Badr-Eddine Boudriki Semlali [17] developed a software architecture to combines complex event processing with remote sensing data from various satellite sensors to resolve the problems met in processing data in near-realtime.The "Ground-Satellite" method, which combines ground-based observations with satellite remote sensing, can obtain more pollutant characteristics at both temporal and spatial scales [18], and will be widely applied in the future.The prediction of PM concentration can provide sufficient information for environmental policy decision-makers to take control measures.The traditional prediction models include numerical prediction models (CMAQ [19], CHIMERE [20], AERMOD [21] et al.) and statistical prediction models [22][23][24].Recently, combinations of machine learning algorithms with numerical prediction models [25] or statistical prediction models [26,27] have achieved good results.For example, Dai H. et al. [28] built a hybrid model (XGBoost-GARCH-MLP) to predict PM 2.5 concentration and volatility, and obtained a better prediction result after using volatility as a benchmark for PM 2.5 .Based on the predicted results of PM concentration, some have used a haze risk assessment model [29] and health risk assessment model [30] to assess affected populations, transportation damage, crop damage area, direct economic loss and comprehensive disaster and health risk, then derived the optimal control measures to minimize losses.
The outbreak of COVID-19 provided a good research platform for people to study air pollution [31][32][33][34][35]. COVID-19 prevalence and the corresponding restrictions resulted in a significant reduction in anthropogenic emissions, but the reduction in emissions was offset by adverse meteorological factors, and the concentrations of PM in the BTH region remained high [36].At the same time, the ozone concentration increases with decreases in NO 2 [37].The researchers speculate that air pollution may increase the incidence, severity and mortality of COVID-19 [38].Ireri Hernandez Carballo's research indicated long-term exposure to air pollutants was positively associated with the incidence of COVID-19 [39].However, wind speed is also a significant factor increasing the number of people infected with COVID-19 compared to higher PM or ozone levels; high wind speed can clean the air of pollutants associated with COVID dynamics, thereby reducing the number of COVID-19 infections [40,41].As such, air quality should not be ignored during citywide shutdowns.In research into PM, we should not only pay attention to the mass concentration, but also to its components.The analysis for the components of PM (including water-soluble ions, metal elements, carbonaceous, etc.) can yield information on health risks and excess mortality assessments, as well as helping in source identification and apportionment.The source identification and apportionment of PM are usually performed using receptor models, including chemical mass balance (CMB) [42,43], principal component analysis (PCA) [44,45], and positive matrix factorization (PMF) [46,47].According to the characteristics of different PM components, the contributions of different pollution sources have been obtained to offer a clear approach to PM control.Currently, BTH, the Yangtze River Delta, and the Pearl River Delta, which are regions dominated by heavy industries, are still the most severely affected by air pollution, and studies on PM sources have also focused on these regions [48][49][50][51].
Handan City is located in southern Hebei province at the intersection of the BTH and central plains economic zones.It has high-emission heavy industries, such as those centered on thermal power, steel, and building materials.It experiences severe air pollution and was one of the 10 cities with the worst air pollution in China between 2005 and 2017 [52].With recent air quality management approaches, the PM concentration has decreased each year, and Handan city ranked first among 168 key cities in air improvement in 2021.However, its air quality still ranks in the bottom 20 [53].Many recent studies have focused on the characteristics, chemical composition and sources of PM 2.5 in Handan [54][55][56].Air quality improved significantly because of the city shutdown between January and April 2020, but rebounded in May 2020 [57], and PM pollution is still pronounced.Yang et al. [58] showed that direct or indirect emissions from the steel industry in eastern Tangshan and western Handan can impose a significant health burden, and concluded that government departments need to reduce emissions from steel enterprises in the BTH region.Soil dust is also an important source of PM emission in the BTH region, of which nearly 60% is from farmland soil [59].The neglected condensable PM (CPM) accounts for nearly half of organic aerosols (OAs), and is also an important component of PM [60].The spatial differences between the four national control stations within Handan are low, and differences in PM 2.5 mass concentrations are not evident, while air pollution in Handan is mainly regional [61].The main PM sources in Handan City are coal combustion, secondary inorganic aerosols, and industrial emissions.Moreover, the regional sources in southern Handan City may substantially contribute to haze pollution in Handan City [62].
Comprehensively promoting the renewal of old communities is one of the main national livelihood projects.Currently, several old communities are being renewed, which involves many residents.The construction dust generated during the renovating of old communities largely impacts PM-particularly TSP and PM 10 -concentrations [63].However, there are large differences between the results of domestic and foreign studies on the impact of construction dust on air quality, which mainly focus on emission factors [64] with few specific case studies of the impact on communities.The control of the COVID-19 epidemic has led to the cessation of most construction work.Due to the urgency of the renovation of old residential areas, some renewals are still under construction.However, the process of renovating residential areas is different from new construction involving only construction workers, as when a residential area is being renovated, most residents remain living in the building, and some residents do not even use respiratory protection.The air quality of a residential area will be affected by community renewal.Although construction processes are covered, and subjected to spraying and other measures, they still have a significant impact on the surrounding air quality [65,66].
This study intends to investigate the impact of community renewal on community air quality by collecting PM from different locations in Handan City and reconstructing the PM composition based on the mass concentration and chemical component characteristics of PM at different locations.We seek to expand the understanding of the health risks associated with PM components, including carcinogenic and non-carcinogenic risks, as well as excess mortality.PCA coupled with multiple linear regression (MLR) is used to estimate the contribution of pollution sources to PM, while the contribution of surrounding areas to Handan's air quality is analyzed according to the backward trajectory and potential source contribution factor (PSCF) method.

Sample and Data
Two sampling sites within Handan were selected for the study.Sampling site 1 is located approximately 100 m away from the national control station of the East Wastewater Treatment Plant, with geographic coordinates 36.61  1.The Hepingli community is undergoing renovation.The sampling site is 10-20 m from the ground, surrounded by residential areas and parks, and has no tall buildings or industrial sources.
Park in the Hepingli neighborhood, with geographic coordinates 36.62°N, 114.50°E. A schematic diagram of the sampling locations is shown in Figure 1.The Hepingli community is undergoing renovation.The sampling site is 10-20 m from the ground, surrounded by residential areas and parks, and has no tall buildings or industrial sources.
The samples were collected using a Laoying 2034 medium flow sampler adjusted to 100 L/min flow rate, while PM2.5 and PM10 were collected using quartz filter membranes, which were baked in a muffle furnace at 650 °C for 4 h before use and weighed after constant temperature and humidity treatment.The samples were collected from 12 to 18 May 2020, during the COVID-19 control period.Each sampling period was 23 h long and the membrane diameter was 77 mm.The collected samples were frozen at −20 °C before analysis and 28 effective sampling films were obtained.The samples were collected using a Laoying 2034 medium flow sampler adjusted to 100 L/min flow rate, while PM 2.5 and PM 10 were collected using quartz filter membranes, which were baked in a muffle furnace at 650 • C for 4 h before use and weighed after constant temperature and humidity treatment.The samples were collected from 12 to 18 May 2020, during the COVID-19 control period.Each sampling period was 23 h long and the membrane diameter was 77 mm.The collected samples were frozen at −20 • C before analysis and 28 effective sampling films were obtained.

Measures of Variables
The filter membranes were placed in a chamber at constant temperature and humidity for more than 24 h, then weighed after removing static electricity.The membranes were then cut into small pieces of 1.77 cm diameter with a cutting tool, and used for OC/EC, water-soluble ion, and inorganic element analyses.Blank filter membranes were analyzed simultaneously.

OC/EC
The OC/EC content was analyzed with a SUNSET RT-4 carbon analyzer.Briefly, a sample was taken in a quartz tube and the EC/OC content was determined using the protocol NIOSH 5040.The analysis process was as follows: OC (partly carbonized) was detected by continuous volatilization under a He atmosphere.Then, it was detected by the oxidative decomposition of the EC escaping under the He/O 2 environment, and the carbonized OC content was confirmed by the change in laser intensity [67].

Water-Soluble Ion
The contents of eight water-soluble ions (F − , Cl − , NO 3 − , SO 4 2− , NH 4 + , K + , Mg 2+ , and Ca 2+ ) were analyzed using a ion chromatograph (DIONEX ICS-1000, Thermo Scientific™, San Diego, CA, USA).Briefly, the sample was placed in a polypropylene vial, soaked in 12 mL ultrapure water for 10 min, and then ultrasonically extracted for 1 h.The sample was filtered through a 0.45 µm filter membrane and then analyzed.The cations and anions were detected on separation columns (IonPac AS23 and IonPac CS12A, Thermo Scientific™, San Diego, CA, USA), respectively.

Elemental
The contents of 19 metals (Ag, Al, Ba, Ca, Cd, Cr, Cu, Fe, K, Mg, Mn, Mo, Ni, Pb, Sr, Ti, V, and Zn) were analyzed by inductively coupled plasma emission spectroscopy (ICPE-9000; Shimadzu, Kyoto, Japan).According to the color depth of the filter membrane, 6-8 pieces were placed in a polytetrafluoroethylene digestion tank.The internal standard (0.010 mL 1000 µg/mL yttrium standard solution) and 5 mL digestion solution (nitric acid:perchloric acid:hydrofluoric acid = 3:1:1) were added and the reflux funnel was covered.The digestion tank was put into the digestion apparatus, and the temperature was raised to 170 • C for 3 h.The reflux funnel was removed and left for 1 h.The tank was then lifted and cooled to room temperature for 30 min.Then, the volume was fixed with 10% nitric acid to 10 mL and analyzed on the machine.This was follows by an analysis of the standard series before the samples were taken, the drawing of standard curves (the correlation coefficient must be more than 0.999), and analysis of the quality control samples (recovery should between 80 and 110%), after which samples could be analyzed.After the testing of 20 samples, a standard solution was analyzed to ensure no major fluctuations from the instrument.A sample blank and laboratory blank were required for every batch.

Enrichment Factor
The enrichment factor (EF) is commonly used to determine whether metal elements in PM are completely derived from crustal elements [68,69], as in Equation ( 1): where  [70].The reference values of crustal elements were adopted from a previous study [71].

Analysis of Secondary Conversion
The sulfur oxidation rate (SOR) and nitrogen oxidation rate (NOR) can be used to characterize the degree of conversion of SO 2 and NO 2 to SO 4 2− and NO 3 − .Calculations of SOR and NOR were performed according to Equations ( 2) and (3).When SOR > 0.25 and NOR > 0.10, the greater conversion of SO 2 and NO 2 into SO 4 2− and NO 3 − occurs in the PM [72].The SOR and NOR of the PM 2.5 and PM 10 at both sampling sites were >0.25, indicating that SO 2 and NO 2 were substantially converted to SO 4 2− and NO 3 − in the air of Handan City, and controlling the SO 2 and NO 2 emissions could effectively reduce the PM concentration.The formulas are shown in Equations ( 2) and ( 3):

PCA-MLR Model
PCA-MLR, with the input of the indicated inorganic and organic source tracers, can quantitatively generate outputs of PM source contributions [73].Compared with the PMF model, this method requires fewer samples and is more suitable for estimating samples in periods of heavy pollution period [74].The principle of PCA-MLR analysis is to reduce the dimension, and summarize different components of particulate matter into several specific factors.The calculation process includes the standardization of mass concentration, the calculation of main factors, and the contribution of identified sources; the formulas are shown in Equations ( 4)-( 7):

Analysis of Secondary Conversion
The sulfur oxidation rate (SOR) and nitrogen oxidation rate (NOR) can be used to characterize the degree of conversion of SO2 and NO2 to SO4 2− and NO3 − .Calculations of SOR and NOR were performed according to Equations ( 2) and (3).When SOR > 0.25 and NOR > 0.10, the greater conversion of SO2 and NO2 into SO4 2− and NO3 − occurs in the PM [72].The SOR and NOR of the PM2.5 and PM10 at both sampling sites were >0.25, indicating that SO2 and NO2 were substantially converted to SO4 2− and NO3 − in the air of Handan City, and controlling the SO2 and NO2 emissions could effectively reduce the PM concentration.The formulas are shown in Equations ( 2) and ( 3):

PCA-MLR Model
PCA-MLR, with the input of the indicated inorganic and organic source tracers, can quantitatively generate outputs of PM source contributions [73].Compared with the PMF model, this method requires fewer samples and is more suitable for estimating samples in periods of heavy pollution period [74].The principle of PCA-MLR analysis is to reduce the dimension, and summarize different components of particulate matter into several specific factors.The calculation process includes the standardization of mass concentration, the calculation of main factors, and the contribution of identified sources; the formulas are shown in Equations ( 4)-( 7): In Equation ( 4), Sij and Cij are the standardized value and the mass concentration of the jth composition species in the ith sample, respectively; Cj and σj refer to the average mass concentration and the standardized deviation of the jth composition species, respectively.In Equation ( 5), Fk is the factor score of the kth source; n represents the number of composition species; aij and Sij are the characteristic vector and the standardized value of the jth composition species in the ith sample, respectively.In Equation ( 6), CPM is the concentration of PM, m refers to the number of sources, βk means the regression coefficient of the kth source, and D is the constant value, while ɳ is the contribution ratio of the factor.

Inhalation Health Risk Assessment
Inhalation health risk assessment is an important tool for assessing the adverse health effects of the exposure of children and adults to air pollutants.Such health risks include carcinogenic risk (CR) resulting from well-defined carcinogenic substances and risk factors (THQ) resulting from non-carcinogenic substances.The carcinogenic risk can be calculated by Equations ( 8) and ( 9), as follows: In Equation ( 4), S ij and C ij are the standardized value and the mass concentration of the jth composition species in the ith sample, respectively; C j and σ j refer to the average mass concentration and the standardized deviation of the jth composition species, respectively.In Equation ( 5), F k is the factor score of the kth source; n represents the number of composition species; a ij and S ij are the characteristic vector and the standardized value of the jth composition species in the ith sample, respectively.In Equation ( 6), C PM is the concentration of PM, m refers to the number of sources, β k means the regression coefficient of the kth source, and D is the constant value, while Atmosphere 2023, 14, x FOR PEER REVIEW 6 of 24 40 indicates extremely high enrichment [70].The reference values of crustal elements were adopted from a previous study [71].

Analysis of Secondary Conversion
The sulfur oxidation rate (SOR) and nitrogen oxidation rate (NOR) can be used to characterize the degree of conversion of SO2 and NO2 to SO4 2− and NO3 − .Calculations of SOR and NOR were performed according to Equations ( 2) and (3).When SOR > 0.25 and NOR > 0.10, the greater conversion of SO2 and NO2 into SO4 2− and NO3 − occurs in the PM [72].The SOR and NOR of the PM2.5 and PM10 at both sampling sites were >0.25, indicating that SO2 and NO2 were substantially converted to SO4 2− and NO3 − in the air of Handan City, and controlling the SO2 and NO2 emissions could effectively reduce the PM concentration.The formulas are shown in Equations ( 2) and ( 3):

PCA-MLR Model
PCA-MLR, with the input of the indicated inorganic and organic source tracers, can quantitatively generate outputs of PM source contributions [73].Compared with the PMF model, this method requires fewer samples and is more suitable for estimating samples in periods of heavy pollution period [74].The principle of PCA-MLR analysis is to reduce the dimension, and summarize different components of particulate matter into several specific factors.The calculation process includes the standardization of mass concentration, the calculation of main factors, and the contribution of identified sources; the formulas are shown in Equations ( 4)- (7): = ∑   ×    =1 (5) In Equation ( 4), Sij and Cij are the standardized value and the mass concentration of the jth composition species in the ith sample, respectively; Cj and σj refer to the average mass concentration and the standardized deviation of the jth composition species, respectively.In Equation ( 5), Fk is the factor score of the kth source; n represents the number of composition species; aij and Sij are the characteristic vector and the standardized value of the jth composition species in the ith sample, respectively.In Equation ( 6), CPM is the concentration of PM, m refers to the number of sources, βk means the regression coefficient of the kth source, and D is the constant value, while ɳ is the contribution ratio of the factor.

Inhalation Health Risk Assessment
Inhalation health risk assessment is an important tool for assessing the adverse health effects of the exposure of children and adults to air pollutants.Such health risks include carcinogenic risk (CR) resulting from well-defined carcinogenic substances and risk factors (THQ) resulting from non-carcinogenic substances.The carcinogenic risk can be calculated by Equations ( 8) and ( 9), as follows: is the contribution ratio of the factor.

Inhalation Health Risk Assessment
Inhalation health risk assessment is an important tool for assessing the adverse health effects of the exposure of children and adults to air pollutants.Such health risks include carcinogenic risk (CR) resulting from well-defined carcinogenic substances and risk factors (THQ) resulting from non-carcinogenic substances.The carcinogenic risk can be calculated by Equations ( 8) and ( 9), as follows: where CR is the carcinogenic risk, C is the components concentration in PM (µg/m 3 ), EF is the exposure frequency (250 day/year), ED is the exposure duration (6 years for children and 24 years for adults), ET is the exposure time (h/day) (8 h/day), AT is the average time of exposure (for non-carcinogens, AT = ED × 365 days × 4 h/day and for carcinogens AT = 70 year × 365 days/year × 24 h), and IUR is inhalation unit risk; the values for carcinogens were taken from USEPA [75].
In general, CR < 1 × 10 −6 indicates that the carcinogenic level caused by the substance poses a negligible risk, and there is no obvious carcinogenic risk.However, when CR > 1 × 10 −4 , this indicates that there is an obvious carcinogenic risk, which may lead to health problems.THQ > 1 and THQ < 1 indicate the presence and absence of noncarcinogenic risk, respectively.

Excess Mortality
The Generalized Additive Model (GAM), a traditional mode of time series analysis, was used to estimate the exposure-response relationship between regional air pollutant concentrations and daily mortality.The basic modeling strategy of time series analysis is the same as in the previous study [76].Based on a zero concentration of air pollutants, an exposure-response model was used to calculate the excess mortality caused by the daily PM 10 pollution level during the study period, as shown in Equation ( 10): where ER kt is excess mortality due to pollutant k on day t, β is the exposure-response relationship coefficient estimated by the regression model, that is, the daily increase in mortality due to each unit increase in the pollutant, and p kt is the average concentration of the k pollutant on day t.

Potential Source Contribution Function (PSCF)
PSCF analysis identifies the main source areas of atmospheric pollutants based on the analysis of air mass trajectories.We used the meteoinfo software developed by Yaqiang Wang et al. (http://www.meteothink.org/,accessed on 1 December 2022).Meteoinfo can be used to calculate the backward trajectory, determine the spatial distribution of potential pollution sources, and combine this with the concentration of pollutants to calculate the PSCF [77,78].The PSCF model divides the study area into i × j grids with an accuracy of 0.5 • × 0.5 • (longitude × latitude), and the PSCF of each grid is calculated as in Equations ( 11)-(13): where Nij is the number of endpoints of trajectory segments on grid ij, and Mij is the number of endpoints of trajectory segments on grid ij with pollutant concentrations higher than the criterion.The threshold value is set to 75 µg/m 3 , the secondary average daily standard value according to the "China Ambient Air Quality Standard".To reduce the uncertainty caused by the small Nij in some grids, a weighting factor Wij was introduced as in Equations ( 12) and (13), and the weighted PSCF value was obtained by multiplying the PSCF value with Wij.

Monitoring Data
Handan's average daily concentrations of air pollutants such as SO 2 , NO 2 , O 3 and CO for 2020 were obtained from the China National Environmental Monitoring Station, the East Wastewater Treatment Plant and Congtai Park.Meteorological data were obtained from the China Meteorological Administration, including daily maximum temperature ( • C), daily minimum temperature ( • C), average temperature ( • C), and relative humidity (%).Handan residents' daily death data in 2020 were obtained from the Chinese Center for Disease Control and Prevention's cause of death registration and reporting information system, including the sex, age, location and underlying cause of death of the deceased.

Characteristics of PM
The time series of PM concentration, gaseous pollutants and meteorological conditions are shown in Figure 2. The average relative humidity is 49%, the average temperature is 22.6 • C, and the prevailing wind direction is southwest wind, with an average wind speed of 3.1 m/s.The average concentration of NO 2 is 26.1 µg/m 3 , the mean concentration of SO 2 is 19.7 µg/m 3 , and the average concentration of CO is 0.7 mg/m 3 , while the mean concentration of O 3-8H is 98.6 µg/m 3 .All the gaseous pollutants are present at values below the Ambient Air Quality Standard (GB3095-2012).The correlations between pollutants and meteorological favors are shown in Figure S1.The concentration of O 3-8h was significantly correlated with temperature (0.80, p < 0.05), and there is a strong negative correlation between PM 2.5 and wind speed (−0.89, p < 0.01), which means an increase in temperature will lead to much more O 3-8h , and a higher wind speed leads to lower PM concentrations.A lower PM concentration means the virus has no carrier, which will reduce the infection rate of COVID-19 [79].
as in Equations ( 12) and (13), and the weighted PSCF value was obtained by multiplying the PSCF value with Wij.

Monitoring Data
Handan's average daily concentrations of air pollutants such as SO2, NO2, O3 and CO for 2020 were obtained from the China National Environmental Monitoring Station, the East Wastewater Treatment Plant and Congtai Park.Meteorological data were obtained from the China Meteorological Administration, including daily maximum temperature (°C ), daily minimum temperature (°C ), average temperature (°C ), and relative humidity (%).Handan residents' daily death data in 2020 were obtained from the Chinese Center for Disease Control and Prevention's cause of death registration and reporting information system, including the sex, age, location and underlying cause of death of the deceased.

Characteristics of PM
The time series of PM concentration, gaseous pollutants and meteorological conditions are shown in Figure 2. The average relative humidity is 49%, the average temperature is 22.6 °C , and the prevailing wind direction is southwest wind, with an average wind speed of 3.1 m/s.The average concentration of NO2 is 26.1 μg/m 3 , the mean concentration of SO2 is 19.7 μg/m 3 , and the average concentration of CO is 0.7 mg/m 3 , while the mean concentration of O3-8H is 98.6 μg/m 3 .All the gaseous pollutants are present at values below the Ambient Air Quality Standard (GB3095-2012).The correlations between pollutants and meteorological favors are shown in Figure S1.The concentration of O3-8h was significantly correlated with temperature (0.80, p < 0.05), and there is a strong negative correlation between PM2.5 and wind speed (−0.89, p < 0.01), which means an increase in temperature will lead to much more O3-8h, and a higher wind speed leads to lower PM concentrations.A lower PM concentration means the virus has no carrier, which will reduce the infection rate of COVID-19 [79].The concentrations of PM 2.5 , PM 10 , and PM 2.5-10 (aerodynamic diameter between 2.5 and 10 µm, referring to coarse PM) and trends of PM 2.5-10 /PM 2.5 at the two sampling sites during the sampling period are shown in Figure 3.The particulate matter increased to the pollution level from a low value at the beginning, and then returned to a low value before the end of sampling, undergoing exactly one pollution process.The PM 2.5 , PM 10 , and PM 2.5-10 concentrations at sampling site 1 were 53 µg/m 3 (24-96 µg/m 3 ), 121 µg/m 3 (70-183 µg/m 3 ), and 68 µg/m 3 (28-134 µg/m 3 ), respectively, with the average PM 2.5-10 /PM 2.5 being 1.6 (0.4-2.7).The PM 2.5 , PM 10 , and PM 2.5-10 concentrations at sampling site 2 were 55 µg/m 3 (25-95 µg/m 3 ), 148 µg/m 3 (85-206 µg/m 3 ), and 93 µg/m 3 (59-153 µg/m 3 ), respectively, with the average PM/PM 2.5-10 being 2.0 (0.9-3.0).The PM 2.5-10 /PM 2.5 ratio is 1.6 (0.4-2.7) > 0.6 at both sampling sites, which indicates the PM was of the dust type [80].The PM 2.5 concentration at sampling site 2 was close to that at sampling site 1, but the PM 2.5-10 and PM 10 concentrations were significantly higher than those at sampling site 1.
Atmosphere 2023, 14, x FOR PEER REVIEW 10 of 24 Figure S2 shows the results of all measured PM mass concentrations, metal elements, water-soluble ions, and OC/EC concentrations during the sampling period.

Crustal Elements
The elements in the Earth's crust mainly include Si, Al, Ca, Fe, Mg, and K. Due to the limitations of the test method, the Si concentration in PM was not tested and was calculated based on the Si/Al ratio in the crust using [Si] = 3.41 [Al] [81].The concentrations of crustal elements in PM10 at both sampling sites were Si > Ca > Al > Fe > K > Mg > other, and six elements accounted for more than 20% of the PM10 mass concentration, while the other elements accounted for less than 0.5% of PM10 mass concentration.In contrast, the concentrations of crustal elements in PM2.5 were Si > K > Ca > Al > Fe > Mg > other, with the six elements accounting for about 10% PM2.5 mass concentration and the contents of the other elements being similar to those in PM10.It can be seen that the crustal elements in PM10 and PM2.5 had the same sources, except for K, and relatively more K was detected in the fine particles.
(  /  )  The paired t-test was used to analyze the significant difference in the levels of PM 2.5 , PM 10 and PM 2.5-10 concentrations between the two sampling sites (Table S1).The PM 2.5 concentrations at the two sampling sites were not significantly different (p = 0.26 > 0.05).However, the PM 2.5-10 and PM 10 concentrations at the two sampling sites were significantly different (p = 0.0002 < 0.01 and 0.0001 < 0.01, respectively).This shows that the difference in PM 10 concentration between two sampling sites was caused by the PM 2.5-10 concentration.The PM 2.5-10 and PM 10 concentrations in sampling site 2 were 25 and 27 µg/m 3 higher than those in sampling site 1, respectively, which shows that the PM 10 concentration around the old neighborhood increased by 27 µg/m 3 due to the renewal of old communities.
Figure S2 shows the results of all measured PM mass concentrations, metal elements, water-soluble ions, and OC/EC concentrations during the sampling period.

Crustal Elements
The elements in the Earth's crust mainly include Si, Al, Ca, Fe, Mg, and K. Due to the limitations of the test method, the Si concentration in PM was not tested and was calculated based on the Si/Al ratio in the crust using [Si] = 3.41 [Al] [81].The concentrations of crustal elements in PM 10 at both sampling sites were Si > Ca > Al > Fe > K > Mg > other, and six elements accounted for more than 20% of the PM 10 mass concentration, while the other elements accounted for less than 0.5% of PM 10 mass concentration.In contrast, the concentrations of crustal elements in PM 2.5 were Si > K > Ca > Al > Fe > Mg > other, with the six elements accounting for about 10% PM 2.5 mass concentration and the contents of the other elements being similar to those in PM 10 .It can be seen that the crustal elements in PM 10 and PM 2.5 had the same sources, except for K, and relatively more K was detected in the fine particles.
The enrichment factors of all determined elements EFx = (Cx/C re f ) PM (Cx/C re f ) crust are shown in Table S2.Influenced by the local metallurgical industry, the EF Cd values in PM 10 and PM 2.5 were 212 and 1422, respectively, indicating that the main sources were anthropogenic.In addition, the EF Pb and EF Zn were also >40, with a very high enrichment effect, while the enrichment factors of other elements were <5, reflecting no substantial contribution of anthropogenic activities.
The EF Mn was significantly linearly correlated with EF Fe (PM 10 : [EF Mn ] = 1.4646 × [EF Fe ] − 0.3587, R 2 = 0.8574; PM 2.5 : [EF Mn ] = 3.265 × [EF Fe ] − 1.0257, R 2 = 0.9383), which shows that the Fe and Mn in PM had the same source.However, the slopes of the regression lines in PM 10 and PM 2.5 are not consistent, reflecting that the sources in PM 2.5 and PM 10 were not consistent.The plot of the trends of enrichment factors and PM concentrations over time (Figure 4) indicates that the enrichment factors peaked two days earlier than the PM concentrations.The crust element is regarded as primary PM; it is assumed that the enrichment of crust element plays an important role in the process of haze formation.Perhaps crust elements can catalyze the generation of secondary PM.
Atmosphere 2023, 14, x FOR PEER REVIEW 11 of 24 addition, the EFPb and EFZn were also >40, with a very high enrichment effect, while the enrichment factors of other elements were <5, reflecting no substantial contribution of anthropogenic activities.
The EFMn was significantly linearly correlated with EFFe (PM10: [EFMn] = 1.4646 × [EFFe] − 0.3587, R 2 = 0.8574; PM2.5: [EFMn] = 3.265 × [EFFe] − 1.0257, R 2 = 0.9383), which shows that the Fe and Mn in PM had the same source.However, the slopes of the regression lines in PM10 and PM2.5 are not consistent, reflecting that the sources in PM2.5 and PM10 were not consistent.The plot of the trends of enrichment factors and PM concentrations over time (Figure 4) indicates that the enrichment factors peaked two days earlier than the PM concentrations.The crust element is regarded as primary PM; it is assumed that the enrichment of crust element plays an important role in the process of haze formation.Perhaps crust elements can catalyze the generation of secondary PM.

Water-Soluble Ions
Sulfur dioxide (SO 2 ) and nitrogen oxides (NO x ) produced from fossil fuel combustion are converted to disulfate (SO 4 2− ) and nitrate (NO 3 − ) through heterogeneous reactions under high relative humidity [82], and the presence of ammonia (NH 3 ) in the air generates secondary ammonium (NH 4 + ) aerosols [83].SO 4 2− , NO 3 − , and NH 4 + are collectively referred to as SNA.In addition, the water-soluble ions in PM include inorganic ions such as Ca 2+ , K + , Na + , Mg 2+ , Cl − , and F − , and organic ions such as formic acid and oxalate.Watersoluble ions in PM can affect the pH of atmospheric precipitation, reduce atmospheric visibility, and impact human health [84].
The concentrations of common water-soluble ions in PM at the two sampling sites are shown in Table S3.The total water-soluble ion (TWSI) concentrations in PM 2.5 at the two sampling sites were 24.8 and 24.0 µg/m 3 , respectively, while those in PM 10 at the two sampling sites were 32.1 and 34.8 µg/m 3 , respectively, and the TWSI concentrations as a percentage of the PM 2.5 and PM 10 concentrations at sampling sites 1 and 2 were 46.8%, 43.6%, 26.5%, and 23.5%, respectively.These proportions are consistent with those reported previously.The proportion of TWSI in PM 2.5 was substantially higher than that in PM 10 , indicating a low TWSI concentration in coarse PM.SNA is the main component of TWSI, accounting for more than 90% of the TWSI in PM 2.5 , while the SNA proportions in the TWSI in PM 10 were 80.1% and 75.9%, which shows that the secondary converted products in SNA were more enriched in fine PM.
The NO 3 − /SO 4 2− in PM is often used to evaluate the contribution of stationary and mobile sources to air quality [85].NO 3 − /SO 4 2− > 1 indicates that the contribution of mobile sources is higher than that of stationary sources; otherwise, the contribution of stationary sources is greater than that of mobile sources.The NO 3 − /SO 4 2− in the PM 2.5 and PM 10 of both sampling sites were >1, indicating that mobile sources contributed more to PM, and controlling mobile sources had a better effect on reducing PM concentrations.) were significantly correlated (R > 0.96; p < 0.01 Figure S3).Regardless, SNA is mainly present in the form of (NH 4 ) 2 SO 4 and NH 4 NO 3 in both PM 10 and PM 2.5 .However, the value of n(NH 4 + ) in PM 10 is smaller than that of n(SO 4 2− ) + n(NO 3 − ).Ca 2+ ions in coarse PM compete with NH 4 + for SO 4 2− to form CaSO 4 .As the Ca 2+ ion concentration in PM 10 increases, the CaSO 4 proportion increases, and more NH 4 + combines with NO 3 − and Cl − to form NH 4 NO 3 and NH 4 Cl, respectively.However, both NH 4 NO 3 and NH 4 Cl easily decompose, particularly in summer when the temperature increases; thus, the NH 4 + ion concentration in PM 10 at the same sampling site was even lower than that in PM 2.5 .On the contrary, PM 2.5 was found to be ammoniarich.In recent years, to reduce NO x emissions, businesses have been adding urea to the combustion process of diesel engines, increasing NH 3 emissions [87], particularly in nonagricultural cities where the original NH 3 emissions were mainly produced by fossil fuel combustion [88].As the NH 3 concentration in the air increases, the NH 4 + concentration in PM increases, which is more conducive to PM generation, particularly PM 10 growth.
The acidity and alkalinity of the PM were calculated according to Equations ( 14) and ( 15), and the regression analysis was performed with AE as the abscissa and CE as the ordinate.The correlation coefficients of cations and anions in PM 2.5 at both sampling sites were approximately 1.1 (Figure 5).The amounts of cations were greater than those of anions, causing the solution to be alkaline.The anion and cation concentrations in PM 10 at sampling site 1 were the same, thus the slope was 1.0 and the PM was neutral.The anion concentrations were higher than the cation concentrations in PM 10 at sampling site 2, thus the PM was acidic.As PM 2.5 grew into PM 10 , the PM gradually transformed from alkaline to acidic, thus the NH 4 + concentration in PM 10 was lower than that in PM

Carbon Fractions
The carbonaceous fractions of PM were divided into organic carbon (OC) and elemental carbon (EC), with OC divided into primary organic carbon emitted directly into the atmosphere and secondary organic carbon generated from precursors through photooxidation reactions.EC is mainly produced from incomplete fossil fuel or biomass combustion, as well as partially from transportation [89].The carbonaceous fraction can generally account for 10-50% PM [90], and the OC/EC ratio is often used to indicate the PM source.The OC/EC ratio for biomass combustion is lower than that for fossil fuel combustion products, and the similar OC and EC proportions in TC indicate consistent carbonaceous PM sources [91,92].An increase in OC concentration in PM implies a decrease in water-soluble ions, and because of the presence of polycyclic aromatic hydrocarbons in OC and the enrichment of elements, the health risk of people exposed to PM will increase [93].
The OC/EC values measured in this study were 1.95-3.23 and 2.58-4.09for PM10 and PM2.5, respectively, indicating that the main source of OCEC in PM is coal combustion [94].

PM Reconstruction
By reconstructing the PM mass, the relationship between the chemical component mass concentration and the total PM mass concentration can be understood.The formula for reconstructing the PM mass [95][96][97] has been revised to Equation (16).RCFM = OM + EC + 1.16 × GM + SNA + Salts + Trace element + others (16) where RCFM represents the reconstructed PM mass, OM represents organic matter, EC represents elemental carbon, GM represents mineral content, 1.16 was used to correct for unmeasured compounds, SNA represents secondary inorganic ions, Salts represents sea salt, Trace element represents the trace element, and others represents boundary moisture or loss.In addition to the carbon detected in organic matter, elements such as H, S, O, and N were also present.Therefore, the quantity of OM needs to be multiplied by a factor that varies regionally from 1.4 to 1.8 [98] in terms of organic composition.OM = 1.6 × OC was used in this study.Calculations of GM were performed using Equation (17).

Carbon Fractions
The carbonaceous fractions of PM were divided into organic carbon (OC) and elemental carbon (EC), with OC divided into primary organic carbon emitted directly into the atmosphere and secondary organic carbon generated from precursors through photooxidation reactions.EC is mainly produced from incomplete fossil fuel or biomass combustion, as well as partially from transportation [89].The carbonaceous fraction can generally account for 10-50% PM [90], and the OC/EC ratio is often used to indicate the PM source.The OC/EC ratio for biomass combustion is lower than that for fossil fuel combustion products, and the similar OC and EC proportions in TC indicate consistent carbonaceous PM sources [91,92].An increase in OC concentration in PM implies a decrease in watersoluble ions, and because of the presence of polycyclic aromatic hydrocarbons in OC and the enrichment of elements, the health risk of people exposed to PM will increase [93].
The OC/EC values measured in this study were 1.95-3.23 and 2.58-4.09for PM 10 and PM 2.5 , respectively, indicating that the main source of OCEC in PM is coal combustion [94].

PM Reconstruction
By reconstructing the PM mass, the relationship between the chemical component mass concentration and the total PM mass concentration can be understood.The formula for reconstructing the PM mass [95][96][97] has been revised to Equation (16).RCFM = OM + EC + 1.16 × GM + SNA + Salts + Trace element + others (16) where RCFM represents the reconstructed PM mass, OM represents organic matter, EC represents elemental carbon, GM represents mineral content, 1.16 was used to correct for unmeasured compounds, SNA represents secondary inorganic ions, Salts represents sea salt, Trace element represents the trace element, and others represents boundary moisture or loss.In addition to the carbon detected in organic matter, elements such as H, S, O, and N were also present.Therefore, the quantity of OM needs to be multiplied by a factor that varies regionally from 1.4 to 1.8 [98] in terms of organic composition.OM = 1.6 × OC was used in this study.Calculations of GM were performed using Equation (17). 2− and NO 3 − do not combine completely with NH 4 + in the PM, so we added three concentrations for reconstruction.Since Ca is present in mineral dust as CaO and CaCO 3 , the coefficient of Ca was calculated using 1.94 [100].
The reconstructed and measured PM values were significantly correlated (R 2 > 0.96), and the slope was greater than 0.81, except for in the PM 10 of sampling site 2 (Figure S4).This indicates that the reconstructed PM accounts for more than 81% of the measured values, and its chemical composition can represent the PM of the measured values.The reconstructed results are consistent with the results reported earlier [101,102].The reconstructed PM 10 of sampling site 2 accounts for 78.2% of the measured values, which shows that under the influence of community renewal, the reconstructed results deviate from the measured values.

Source Apportionment of PM
To some extent, the correlation between different components in particulate matter can reflect the similarity of their sources [103].The correlations between the components of PM 2.5 and PM 10 are shown in Figure S5.In the figure, red indicates positive correlation, blue indicates negative correlation, larger circles indicate better correlation, and asterisks indicate significant correlation (p < 0.05).The correlations between the components of PM 10 in the two sampling sites are better than those of PM 2.5.The correlation between crustal elements such as Al, Fe, K, Mg, Ti, etc., is significant in PM 10 , but not in PM 2.5 .The correlation between secondary pollutants NO 3 − , SO 4 2− , NH 4 + and OC is good, except for the poor correlation between NH 4 + and OC in the PM10 of sampling site 2. A good correlation between Cl − and K indicates that they share common sources.
The results of the PCA of PM show that 85.9% of the total variance can be expressed by three factors.Table 1 shows the load matrix of the PC rotation factor for the components in PM.Al, Ca, Fe, K, Mg, Ti, Fe and V showed higher loads in factor 1 (F1), with the values all higher than 0.82.F1 was considered as a source of dust.The loads of OC, Cd, Pb, Zn, NO 3 − , SO 4 2− and NH 4 + were high in factor 2 (F2).F2 was found to mainly contain secondary pollutants, which are oxidized by their precursors of SO 2 , NO 2 , VOC and NH 3 , and Pb is the indicator of motor vehicle source, so F2 was considered as mixed source including secondary transformation, vehicle exhaust and fossil fuels.Factor 3 (F3) contained high levels of Cl − , but no K + , and this was considered to be a result of COVID-19.The disinfectant used to disinfect COVID-19 is chlorine-rich [104].We predict that F3 was the source of the high presence of disinfectants, but more research is needed to confirm this.
In order to further quantitatively analyze the main pollution sources and their relative contributions to PM, MLR analysis was performed with the normalized principal factor score as independent variable and PM concentration as the dependent variable, and we derived the regression equation: Z = 0.831F1 + 0.513F2 + 0.031F3.This equation indicates that during the sampling period, the later stage of the COVID-19 epidemic, 60.4% of the PM was contributed by dust, 37.3% by secondary sources, including motor vehicles and industrial coal-fired sources, and the remaining 2.3% was assumed to be derived from disinfectants.

PSCF Analysis
In this study, the 48 h backward trajectory of the air mass from Handan City (36.61 • N, 114.19 • E, 500 m above ground) from 12 May to 18 May 2020 (during the COVID-19 pandemic) was calculated using meteorological data downloaded from the National Center for Environmental Prediction (NCEP) Global Data Assimilation System (GDAS).The air mass trajectories were clustered into five types using the Euler method (Figure S6).Cluster 1 was short-distance transportation from areas such as Anyang and Kaifeng, contributing 33.93% of the trajectory; cluster 2 was from Inner Mongolia via Shaanxi and southern Henan and Hebei, accounting for 22.2% of the trajectory.Both clusters 1 and 2 reflect the south of Handan City, accounting for 56% of the trajectory; the three trajectories from the northwest of Handan account for 20.8%, 16.7%, and 6.6% of the trajectory, respectively.
The potential source contributions of PM 10 , PM 2.5 , SO 2 , and NO 2 in Handan were calculated using PSCF (Figure S6).From Figure S6, it can be seen that the potential PM 10 sources were mainly local (WPSCF > 0.7), and the spatial patterns of the potential PM 2.5 , SO 2 , and NO 2 sources were similar.The main source areas were essentially distributed in southern Hebei, while most of Henan and northern Anhui constituted the potential PM 2.5 , SO 2 , and NO 2 sources, with relatively high WPSCF.The results of the PSCF calculation are similar to those of previous studies [105,106]; in general, local sources of PM 10 in Handan City were significant contributors, while the short-distance transportation of PM 2.5 , SO 2 , and NO 2 in the southern region was obvious.The results are consistent with the weak southwest wind during the sampling period.

Health Risk Assessment
Health assessments need to include the concentrations of chromium in trivalent CR (III) and hexavalent CR (VI), whereas the present study measured total concentrations of Cr.We took as reference [107], which assumed a concentration ratio of CR (VI) to CR (III) of 1:6.The carcinogenic and non-carcinogenic risks related to human exposure through inhalation are shown in Table 2.The CR of PM 10 in two sampling sites exceeded the acceptable limit of 1 × 10 −6 , and the main contributor was CR (VI), as shown in Figure S7.The results are concerning, especially in sampling site 2, where the renewal of old community construction has led to higher values of CR than in sampling site 1.The THQ values of the two sampling sites are both below 1, which implies the absence of non-carcinogenic risk.Hongya Niu's [108] findings are similar, but Xing Li [109] found that the non-carcinogenic and carcinogenic risks of PM pollution to children exceeded the acceptable limits in almost all cities in Hebei province.The carcinogenic risk of the PAH in particulate matter is generally higher than the acceptable limit in the North China Plain [110].The health risks associated with PM-especially the carcinogenic risks-are significantly higher than the acceptable limit, and reducing the carcinogens in PM should be the priority when developing control measures.

Excess Mortality
Excess mortality due to PM 10 at both sites is shown in Figure 6, which shows that excess mortality was higher at sampling site 2 than at sampling site 1, with an average ∆ER 0.11 and a relative increase of 23%.The particulate pollution caused by the renovation of old residential areas not only increases the carcinogenic risk, but also significantly increases the excess mortality.The CR of PM10 in two sampling sites exceeded the acceptable limit of 1 × 10 −6 , and the main contributor was CR (VI), as shown in Figure S7.The results are concerning, especially in sampling site 2, where the renewal of old community construction has led to higher values of CR than in sampling site 1.The THQ values of the two sampling sites are both below 1, which implies the absence of non-carcinogenic risk.Hongya Niu's [108] findings are similar, but Xing Li [109] found that the non-carcinogenic and carcinogenic risks of PM pollution to children exceeded the acceptable limits in almost all cities in Hebei province.The carcinogenic risk of the PAH in particulate matter is generally higher than the acceptable limit in the North China Plain [110].The health risks associated with PMespecially the carcinogenic risks-are significantly higher than the acceptable limit, and reducing the carcinogens in PM should be the priority when developing control measures.

Excess Mortality
Excess mortality due to PM10 at both sites is shown in Figure 6, which shows that excess mortality was higher at sampling site 2 than at sampling site 1, with an average ΔER 0.11 and a relative increase of 23%.The particulate pollution caused by the renovation of old residential areas not only increases the carcinogenic risk, but also significantly increases the excess mortality.

Comparison with Other Studies
A comparison of the characteristics, source apportionments and health risk assessments of PM in Handan city in this study with other similar studies is shown in Tables 3 and 4. It can be seen from the comparison that the content of this study is more comprehensive.Studies on PM in Handan generally focus on winter, when air quality is obviously worse than it is in summer.Research on the PM during winter seems to be more valuable.However, during the period of COVID-19, it was found that PM could carry the virus and aggravate its spread, so research on summer PM cannot be neglected.
As regards source apportionment, the main pollution source identified in most of the literature was secondary conversion, or mixed factors such as secondary conversion and motor vehicle and/or coal burning, but dust was the main source in our paper.There are two possible reasons for this: on the one hand, the times of analysis are different-winter is more beneficial to the transformation of secondary pollutants than summer, leading to an increase in the proportional contribution; on the other hand, the contribution of dust to PM may be increased because sampling site 2 in this study is affected by the renewal of old communities.The F3 discovered in this study also differs significantly from other studies.

Conclusions
In the latter stage of COVID-19, we chose two sites in Handan to collect PM samples; one of the sites is undergoing community renovation.The characteristics of PM 10 and PM 2.5 were to be learned.It was found that the concentration of PM 10 increased by 27 µg/m 3 during the renovation of the old residential area, which should be of great concern to construction units and the ecological and environmental protection departments.At the same time, the concentration of ammonium in PM 10 was lower than that in PM 2.5 at the same sampling site, the mechanism of which requires further study.The source apportionment shows that dust is the most important contributor, accounting for more than 60% of the total, which provides a means of controlling PM.Although the carcinogenic risk caused by PM is no more than 1 × 10 −4 , this value is higher than the acceptable limit of 1 × 10 −6 , which needs to be taken seriously, because the transformation of old residential areas leads to an increase in PM 10 concentration, causing a 23% increase in the excess mortality rate.Meanwhile, PM may be a carrier for viruses; under conditions of low wind speed, the spread of the virus will be accelerated, so it is necessary to take protective measures during periods of pollution.
As sampling took place during the COVID-19 lockdown period, the number of samples is small; more sampling sites and a longer study period are thus needed to improve the accuracy of the results.

Figure 1 .
Figure 1.Schematic diagram of sampling locations.Figure 1.Schematic diagram of sampling locations.

Figure 1 .
Figure 1.Schematic diagram of sampling locations.Figure 1.Schematic diagram of sampling locations.

Figure 2 .
Figure 2. Time series of PM concentration, gaseous pollutants and meteorological conditions.

Figure 2 .
Figure 2. Time series of PM concentration, gaseous pollutants and meteorological conditions.

Figure 4 .
Figure 4. Relationship between Cd, Pb and Zn enrichment factors and particle concentration over time (A: the PM10 of sample site 1; B: the PM10 of sample site 2; C: the PM2.5 of sample site 1; D: the PM2.5 of sample site 2).

Figure 4 .
Figure 4. Relationship between Cd, Pb and Zn enrichment factors and particle concentration over time (A: the PM 10 of sample site 1; B: the PM 10 of sample site 2; C: the PM 2.5 of sample site 1; D: the PM 2.5 of sample site 2).

Figure 5 .
Figure 5. Relationship between AE and CE in PM.

Figure 5 .
Figure 5. Relationship between AE and CE in PM.

Figure 6 .
Figure 6.The excess mortality caused by PM10 between the two sampling sites.

Figure 6 .
Figure 6.The excess mortality caused by PM 10 between the two sampling sites.
• N, 114.53 • E. Sampling site 2 is located approximately 1.1 km away from the national control station of Congtai Park in the Hepingli neighborhood, with geographic coordinates 36.62 • N, 114.50 • E. A schematic diagram of the sampling locations is shown in Figure 2.5 .

Table 1 .
Rotated principal component analysis for components in PM.
Note: The boldface is the load of the component with larger load in this factor.

Table 2 .
The carcinogenic risk and non-carcinogenic risks of human exposure to PM.

Table 2 .
The carcinogenic risk and non-carcinogenic risks of human exposure to PM.

Table 3 .
Comparison with other studies on characteristics of PM 2.5 .

Table 4 .
Comparison with other studies on source apportionment.