Source Identiﬁcation of Trace Elements in PM 2.5 at a Rural Site in the North China Plain

: An intensive sampling of PM 2.5 was conducted at a rural site (Gucheng) in the North China Plain from 22 October to 23 November 2016. A total of 25 elements (Al, Na, Cl, Mg, P, S, K, Ca, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, As, Se, Br, Sr, Cd, Ba, Pb, and Sb) from PM 2.5 ﬁlter samples collected daily were measured using a wavelength dispersive X-ray ﬂuorescence spectrometer. Cl, S, and K were the most abundant elements, with average concentrations of 2077.66 ng m − 3 (range 118.88–4638.96 ng m − 3 ), 1748.78 ng m − 3 (range 276.67–4335.59 ng m − 3 ), and 1287.07 ng m − 3 (range 254.90–2748.63 ng m − 3 ), respectively. Among noncrustal trace metal elements, the concentration of Zn was the highest, with an average of 397.74 ng m − 3 (range 36.45–1602.96 ng m − 3 ), followed by Sb and Pb, on average, of 299.20 ng m − 3 and 184.52 ng m − 3 , respectively. The morphologies of PM 2.5 samples were observed using scanning electron microscopy. The shape of the particles was predominantly spherical, chain-like, and irregular. Positive matrix factorization analysis revealed that soil dust, following by industry, secondary formation, vehicle emissions, biomass and waste burning, and coal combustion, were the main sources of PM 2.5 . The results of cluster, potential source contribution function, and concentration weighted trajectory analyses suggested that local emissions from Hebei Province, as well as regional transport from Beijing, Tianjin, Shandong, and Shanxi Province, and long-range transport from Inner Mongolia, were the main contributors to PM 2.5 pollution.


Introduction
In the past few decades, China has experienced rapid economic growth, urbanization, industrialization, and population expansion. As a result of these factors, air pollution has become a primary environmental concern. The North China Plain (NCP) is one of the most polluted regions in China [1,2]. Severe air pollution frequently occurs in autumn and winter, which has become the focus and difficulty of improving the quality of the atmospheric environment. As an important pollutant in the atmosphere, particulate matter, especially that less than 2.5 µm in aerodynamic diameter (PM 2.5 ), can stay in the atmosphere for days to weeks and be transported over hundreds to thousands of kilometers, which has a profound effect on human health, visibility, and climate change.
Exposure to PM 2.5 can cause numerous adverse health effects. Increases in morbidity and mortality related to respiratory, cardiovascular, chronic pulmonary, and neurodegenerative diseases, as well a reduction in life expectancy, have been observed in exposed populations [3,4]. Because of its large surface area, small size, and the fact that it contains a considerable amount of harmful toxic characteristics, PM 2.5 can be breathed into the deepest recesses of the lung. It is deposited on the alveoli and even escapes into the blood stream to harm the health of humans [5]. PM 2.5 can also enter the human body through ingestion and dermal contact. Studies have indicated that trace elements in PM 2.5 play a key role in such adverse effects. Fe, Co, Ni, Cu, Zn, V, Cr, Mn, As, Pb, and Cd may boost the generation of reactive oxygen species (ROS). Excess ROS can overwhelm the antioxidant defense system in the body, causing oxidative stress, inflammation, and disease [6][7][8]. Therefore, trace elements pose a substantial threat to human health, even though they account for only a small fraction of the total mass of PM 2.5 .
Trace elements have proven to be useful tracers [9][10][11] and are widely used as input to aid in the interpretation of specific emission sources in receptor modeling studies. The positive matrix factorization (PMF) model developed by the Environmental Protection Agency (EPA) of USA is an advanced factor analysis technique, which employs a flexible modeling method to effectively use information in data and to identify the possible source contributions without the source profiles. This model has been used worldwide for the source apportionment of PM 2.5 [12,13]. The characterization and source identification of trace elements in PM 2.5 play an important role in the prevention and control of air pollution. Such approaches can improve air quality and reduce negative effects on human health, thus attracting extensive attention among researchers.
The NCP is a densely populated region with intensive agricultural and industrial activities. Several studies have measured atmosphere pollutants at Gucheng [14][15][16][17], a rural site in the core area of the NCP, and have shown that this site is ideally situated for monitoring the atmospheric pollution levels in the NCP and has favorable regional representativeness. The observed data indicate that Gucheng is a relatively polluted rural area with high levels of gaseous and particle pollutants. It was reported that the daily mean concentration of PM 2.5 reached a maximum of 394 µg m −3 , and only 8% of daily PM 2.5 concentrations were lower than the Chinese National Ambient Air Quality Standards (NAAQS, GB3095-2012) for daily PM 2.5 of 75 µg m −3 from December 2016 to January 2017 [16]. Studies related to source appointment of trace elements in PM 2.5 at Gucheng are still lacking. Thus, we undertook this study to characterize and identify potential sources of trace elements in PM 2.5 collected at Gucheng. The morphologies of PM 2.5 samples were observed using scanning electron microscopy (SEM). Enrichment factor (EF) and PMF were used for the source identification of trace elements in PM 2.5 . In addition, clusters of air mass back trajectories, the potential source contribution function (PSCF), and the concentration weighted trajectory (CWT) were utilized to explore the transport pathways and likely source regions.

Sampling Site Description
Sampling was conducted at Gucheng, and Agrometeorological Field Experiment Station (39 • 08 N, 115 • 40 E) of the Chinese Academy of Meteorological Sciences, which is located in Dingxing County, Hebei Province. As shown in Figure 1, the site is approximately 110 km southwest of Beijing, 130 km west of Tianjin, and 160 km northeast of Shijiazhuang. It is surrounded by farmland, villages, and highway and railway networks.

Sampling Method and Analysis
PM2.5 samples were collected on the rooftop (about 10 m above the ground) of the observation building from 22 October to 23 November 2016. Each daily sampling was scheduled with two intervals of daytime (from 8:00 to 20:00) and nighttime (from 20:00 to 8:00 the following day). A highvolume sampler (GUV-15HBL1, Thermo Fisher Scientific, Waltham, MA, USA) equipped with a PM2.5 impactor (G1200-41) was employed to collect particles at a nominal flow rate of 1.13 m 3 min −1 using quartz fiber filters (20.32 × 25.40 cm, Pall, New York, NY, USA). The filters were baked at 550 °C for 6 h to remove organic species before use. After sampling, all filters were carefully wrapped in aluminum foil and stored at −20 °C until analysis. During the sampling period, the average temperature was 7.8 °C, the relative humidity was 67.5%, and the mean wind speed was 1.5 m s −1 . A total of 60 PM2.5 samples were taken after excluding damaged samples.
Twenty-five elements in PM2.5, namely Al, Na, Cl, Mg, P, S, K, Ca, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, As, Se, Br, Sr, Cd, Ba, Pb, and Sb, were analyzed using a wavelength dispersive X-ray fluorescence spectrometer (RIX3000, Rigaku, Tokyo, Japan), which provides a straightforward technique for determining elemental compositions without the pretreatment of filters. An X-ray tube (Rh anode, 50 kV, 50 mA), with a coarse slit and filed diaphragm diameter of 30 mm, was used for measurement. The quantitative calibration of the system was performed using standard reference samples . The calibration was verified using the NIST standards, and the relative error between measured value and standard value was within ± 10%.
The morphologies of PM2.5 samples were observed using SEM. In brief, a punch of the filter was fixed firmly on the sample holder with conductive tape. For higher conductivity and better observation, gold particles were sprayed onto the surface of the filter prior to SEM characterization. Magnification images were acquired using an FEI Quanta 650 FEG SEM operated at an accelerating voltage of 10.0 kV at 68.5 mPa.
EF has been widely employed to assess the anthropogenic influences on elements in PM2.5. This was calculated using the following equation:

Sampling Method and Analysis
PM 2.5 samples were collected on the rooftop (about 10 m above the ground) of the observation building from 22 October to 23 November 2016. Each daily sampling was scheduled with two intervals of daytime (from 8:00 to 20:00) and nighttime (from 20:00 to 8:00 the following day). A high-volume sampler (GUV-15HBL1, Thermo Fisher Scientific, Waltham, MA, USA) equipped with a PM 2.5 impactor (G1200-41) was employed to collect particles at a nominal flow rate of 1.13 m 3 min −1 using quartz fiber filters (20.32 × 25.40 cm, Pall, New York, NY, USA). The filters were baked at 550 • C for 6 h to remove organic species before use. After sampling, all filters were carefully wrapped in aluminum foil and stored at −20 • C until analysis. During the sampling period, the average temperature was 7.8 • C, the relative humidity was 67.5%, and the mean wind speed was 1.5 m s −1 . A total of 60 PM 2.5 samples were taken after excluding damaged samples.
Twenty-five elements in PM 2.5 , namely Al, Na, Cl, Mg, P, S, K, Ca, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, As, Se, Br, Sr, Cd, Ba, Pb, and Sb, were analyzed using a wavelength dispersive X-ray fluorescence spectrometer (RIX3000, Rigaku, Tokyo, Japan), which provides a straightforward technique for determining elemental compositions without the pretreatment of filters. An X-ray tube (Rh anode, 50 kV, 50 mA), with a coarse slit and filed diaphragm diameter of 30 mm, was used for measurement. The quantitative calibration of the system was performed using standard reference samples . The calibration was verified using the NIST standards, and the relative error between measured value and standard value was within ± 10%.
The morphologies of PM 2.5 samples were observed using SEM. In brief, a punch of the filter was fixed firmly on the sample holder with conductive tape. For higher conductivity and better observation, gold particles were sprayed onto the surface of the filter prior to SEM characterization. Magnification images were acquired using an FEI Quanta 650 FEG SEM operated at an accelerating voltage of 10.0 kV at 68.5 mPa.
EF has been widely employed to assess the anthropogenic influences on elements in PM 2.5 . This was calculated using the following equation: where EF x is the EF of element x and r is a reference element for crustal material. (C x /C r ) PM and (C x /C r ) crust represent the concentration ratio of x to r in PM 2.5 and crust, respectively. Generally, an EF x value below 10 suggests that the element x primarily originates from a crustal source. By contrast, an EF x value greater than 10 implies that the element x is enriched and associated with human activities [18]. In this study, Al was chosen as a reference element, and the concentrations of elements in crust refer to their background concentrations in the surface soil of China [19].

Source Apportionment
PMF is a multivariate factor analysis technique based on the mathematical receptor model. EPA PMF version 5.0 was employed to identify the potential sources of elements in PM 2.5 . A detailed description of the PMF model can be found in related studies [20,21]. A speciated data set can be deemed as a data matrix X of i by j dimensions with uncertainties u, where i and j, respectively, represent the number of samples and chemical species. The goal of the PMF model is to identify the number of factors, p, the species profile, f, of each factor, and the amount of mass, g, contributed by each factor to individual samples. This can be expressed as: where e ij is the residual for each sample or species. Factor contributions and profiles were derived by the PMF model by minimizing the objective function Q: Two inputs are required to run PMF model, namely concentration and corresponding uncertainty. When the concentration is less than or equal to the method detection limit (MDL), it is replaced with half of the MDL. The uncertainty (Unc) is calculated as follows: When the concentration is higher than the MDL, the calculation is: In this study, 60 samples were included in the PMF analysis. Signal-to-noise (S/N) ratios were used to categorize species [22]. If the S/N ratio was greater than 0.5 but less than 1, the corresponding species was categorized as "weak". If the S/N ratio was less than 0.5, the species was categorized as "bad." Thus, V, Cr, Co, and Ni were set as "bad" and excluded from the rest of the analysis.

Cluster and Source Regions
Air mass back trajectories over 48 h arriving at the Gucheng site during the sampling period were calculated and clustered using the geographic information system-based software, TrajStat [23]. The meteorological data used for the calculation was obtained from the Global Data Assimilation System of the NOAA National Center for Environmental Prediction with a resolution of 1 • × 1 • . The model was run four times (00:00, 06:00, 12:00, and 18:00 UTC) a day at 100 m above ground level.
The potential source contribution function (PSCF) and contribution weighted factor CWT methods were implemented using TrajStat to identify the potential source regions of PM 2.5 [24]. The study field of 60 • E to 125 • E and 30 • N to 65 • N with a resolution of 1 • × 1 • was divided into small equal grid cells (i × j). The PSCF value in the ij cell was defined as: where m ij and n ij represent the number of "polluted" trajectory endpoints and trajectory endpoints that fall in the ij cell, respectively. A PM 2.5 concentration of 75 µg m −3 was chosen as the pollution threshold. In order to determine the relative contribution of different potential source region, the CWT method was used. In the CWT method, each grid cell was assigned a weighted concentration by averaging the sample concentrations that have associated trajectories crossing that grid cell, as follows: where C ij is the average weighted concentration in the ij cell, l is the index of the trajectory, and M is the total number of trajectories. C l and τ ijl represent the concentration observed on arrival of trajectory l and the time spends in the ij cell by trajectory l, respectively. Grid cells with high PSCF ij and C ij values indicated the potential source regions.
To reduce the uncertainty in cells with a small n ij value, the PSCF and CWT value was multiplied by a weight function, W ij , which was defined as:

General Characteristics
The concentrations of PM 2.5 and 25 trace elements measured at Gucheng site are summarized in Table 1 ) were significantly higher than those in the guidelines of the WHO (6.6 ng m −3 and 5 ng m −3 for As and Cd, respectively), the EU Directive and the NAAQS of China (6 ng m −3 and 5 ng m −3 for As and Cd, respectively), indicating potential threats to human health.    The concentration of S (1748.78 ng m −3 ) was lower than that measured in Beijing (3700 ng m −3 ) and Baoding (7800 ng m −3 ), comparable with that that recorded in Ulaanbaatar (1740 ng m −3 ), and higher than that observed in Athens (1409.9 ng m −3 ), Košetice (1011.6 ng m −3 ), and Treviso (938 ng m −3 ) ( Table 2). The concentration of K (1287.07 ng m −3 ) was lower than that in Baoding (1920 ng m −3 ) and higher than that in Terviso (688 ng m −3 ). The concentrations of Mn, Cu, Zn, and Pb were lower than that in Baoding, but higher than that in Beijing and other sites. In addition, Gucheng, Beijing, Baoding, Ulaanbaatar, and Treviso had high levels of crustal elements, such as Al Mg, Ca, and Fe.
A morphological analysis of PM 2.5 was conducted using SEM. According to SEM images (Figure 3), the particles collected on the quartz fiber filter mainly had a spherical, chain-like, and irregular shape, and such particles varied in size (no larger than 2.5 µm). Spherical particles are easily recognized. They generally originate from fluid melts due to high-temperature combustion, and over 95% of fly ash particles are smoothly spherical [31]. Coal combustion is the major source of this type of particle [32]. Chain-like particles are typical soot aggregates consisting of many ultrafine particles, which are mainly produced by combustion emissions, such as from gasoline and diesel exhaust emissions and coal and biomass burning [33]. Irregular particles are mainly composed of mineral matter from windblown dust, resuspended road dust, and agricultural and construction dust [34].
An EF analysis was performed to evaluate the contributions of crustal and noncrustal sources of elements in PM 2.5 . As shown in Figure 4, the EF values for Ti, Mg, Fe, and Na were lower than 10, demonstrating that these elements were mainly from crustal sources. Ca, K, as well as Mn, Sr, V, and Ba had values ranging from 10 to 100, representing slight enrichment due to an increase in anthropogenic activities, such as agriculture and construction (Ca) and biomass burning (K). The EF values for Cr, Co, Ni, Cu, Sn, and As ranged from 100 to 1000, indicating that these elements were moderately enriched due to human activities. Zn, Pb, Br, Se, Cd, and Sb had values larger than 1000, providing evidence of high enrichment from anthropogenic sources. These elements may be derived from emission-related sources.
Atmosphere 2020, 11, 179 8 of 15 moderately enriched due to human activities. Zn, Pb, Br, Se, Cd, and Sb had values larger than 1000, providing evidence of high enrichment from anthropogenic sources. These elements may be derived from emission-related sources.

Source Apportionment
On the basis of the PMF modeling results, six factors were identified, namely soil dust, vehicle emissions, industry, biomass and waste burning, secondary formation, and coal combustion. The source profiles are shown in Figure 5, and the relative contributions from each source to PM2.5 are illustrated in Figure 6.
The first factor, soil dust, is characterized by high loadings of Ca (85.09%), Al (68.24%), Mg (68.00%), Ti (58.35%), and Fe (46.50%). These typical crustal elements may be derived from natural and anthropogenic sources, including natural mineral dust from desert, arid and semiarid regions, the resuspension of road dust by vehicles, and dust generated by agricultural and construction activities [35]. This factor comprised a major proportion (26.1%) of PM2.5.
The second factor, vehicle emissions, results in high percentages of Cu (55.13%), Zn (53.03%), Pb (37.59%), and Br (36.11%). The contribution of this source to PM2.5 was 14.0%. Cu is a key tracer of brake-wear particle emissions [36]. Zn is associated with tire wear. It is employed in tires as an activator for vulcanization, and it is also released from additives in lubricating oils for engines [37]. moderately enriched due to human activities. Zn, Pb, Br, Se, Cd, and Sb had values larger than 1000, providing evidence of high enrichment from anthropogenic sources. These elements may be derived from emission-related sources.

Source Apportionment
On the basis of the PMF modeling results, six factors were identified, namely soil dust, vehicle emissions, industry, biomass and waste burning, secondary formation, and coal combustion. The source profiles are shown in Figure 5, and the relative contributions from each source to PM2.5 are illustrated in Figure 6.
The first factor, soil dust, is characterized by high loadings of Ca (85.09%), Al (68.24%), Mg (68.00%), Ti (58.35%), and Fe (46.50%). These typical crustal elements may be derived from natural and anthropogenic sources, including natural mineral dust from desert, arid and semiarid regions, the resuspension of road dust by vehicles, and dust generated by agricultural and construction activities [35]. This factor comprised a major proportion (26.1%) of PM2.5.
The second factor, vehicle emissions, results in high percentages of Cu (55.13%), Zn (53.03%), Pb (37.59%), and Br (36.11%). The contribution of this source to PM2.5 was 14.0%. Cu is a key tracer of brake-wear particle emissions [36]. Zn is associated with tire wear. It is employed in tires as an activator for vulcanization, and it is also released from additives in lubricating oils for engines [37].

Source Apportionment
On the basis of the PMF modeling results, six factors were identified, namely soil dust, vehicle emissions, industry, biomass and waste burning, secondary formation, and coal combustion. The source profiles are shown in Figure 5, and the relative contributions from each source to PM 2.5 are illustrated in Figure 6.
The first factor, soil dust, is characterized by high loadings of Ca (85.09%), Al (68.24%), Mg (68.00%), Ti (58.35%), and Fe (46.50%). These typical crustal elements may be derived from natural and anthropogenic sources, including natural mineral dust from desert, arid and semiarid regions, the resuspension of road dust by vehicles, and dust generated by agricultural and construction activities [35]. This factor comprised a major proportion (26.1%) of PM 2.5 .
The second factor, vehicle emissions, results in high percentages of Cu (55.13%), Zn (53.03%), Pb (37.59%), and Br (36.11%). The contribution of this source to PM 2.5 was 14.0%. Cu is a key tracer of brake-wear particle emissions [36]. Zn is associated with tire wear. It is employed in tires as an activator for vulcanization, and it is also released from additives in lubricating oils for engines [37]. Motor vehicle exhaust was previously the largest source of Pb and Br emissions [38] until leaded gasoline use was banned and the Pb concentration in unleaded gasoline was limited to below 0.005 g L −1 in China in 2000 (GB17930-1999). Since then, Pb emissions from motor vehicle gasoline have substantially declined, but vehicular exhaust emissions remain an important source of Pb emissions [15], owing to the sustained increase in the number of automobiles.
Atmosphere 2020, 11, 179 9 of 15 Motor vehicle exhaust was previously the largest source of Pb and Br emissions [38] until leaded gasoline use was banned and the Pb concentration in unleaded gasoline was limited to below 0.005 g L −1 in China in 2000 (GB17930-1999). Since then, Pb emissions from motor vehicle gasoline have substantially declined, but vehicular exhaust emissions remain an important source of Pb emissions [15], owing to the sustained increase in the number of automobiles.  The third factor, industry sources, results in an abundance of Cd (62.09%), Se (52.37%), Ba (46.30%), Sr (40.92%), and Mn (38.35%), which originate from steel metallurgy and nonferrous metal production [9,[39][40][41]. The contribution of this factor to PM2.5 was 22.9%. Motor vehicle exhaust was previously the largest source of Pb and Br emissions [38] until leaded gasoline use was banned and the Pb concentration in unleaded gasoline was limited to below 0.005 g L −1 in China in 2000 (GB17930-1999). Since then, Pb emissions from motor vehicle gasoline have substantially declined, but vehicular exhaust emissions remain an important source of Pb emissions [15], owing to the sustained increase in the number of automobiles.  The third factor, industry sources, results in an abundance of Cd (62.09%), Se (52.37%), Ba (46.30%), Sr (40.92%), and Mn (38.35%), which originate from steel metallurgy and nonferrous metal production [9,[39][40][41]. The contribution of this factor to PM2.5 was 22.9%. The third factor, industry sources, results in an abundance of Cd (62.09%), Se (52.37%), Ba (46.30%), Sr (40.92%), and Mn (38.35%), which originate from steel metallurgy and nonferrous metal production [9,[39][40][41]. The contribution of this factor to PM 2.5 was 22.9%.
The fourth factor, biomass and waste burning, is represented by high levels of Cl (63.12%), K (34.61%), and Sb (31.45%). Cl is released during the combustion of high-chlorine biomass, such as straw, grass, wood barks, and grains [42,43]. K, another element with a high concentration in natural biomass, has been extensively used as a tracer element for biomass burning [44]. In the NCP, some rural households still use agriculture and forestry residue, including crop residues, weeds, branches, and leaves, as biomass fuel for cooking or heating [45]. Furthermore, studies have demonstrated that Cl is related to waste incineration [46,47]. Chlorinated plastics, such as polyvinylchloride, chlorinated polyethylene, and polyvinyl dichloride, and chloride salts in food residue, are the main sources of Cl in waste burning [48]. Refuse incinerators are often major sources of Sb [49]. This factor accounted for 10.9% of PM 2.5 .
The fifth factor, secondary formation, is dominated by a high concentration of S (69.96%) [50,51], accounting for 15.2% of PM 2.5 emissions. S is related to secondary sulfates, which are likely generated by the homogeneous or heterogeneous reaction of the precursor gas SO 2 emitted from coal burning [13]. Secondary formation is affected by meteorological conditions, such as temperature, humidity, and solar radiation.
The sixth factor, coal combustion, is typified by an abundance of As (80.79%) and Pb (23.07%); the contribution of this factor to PM 2.5 was 10.9%. Coal is the major energy source in China [52]. The main consumers of coal include industry and power plants, and it is also used in small amounts for heating residences and for other purposes [53]. As is recognized as a tracer of coal combustion [54]. Coal combustion was found to be the predominant source of Pb emissions after the phasing out of leaded gasoline [55]. Exposure to these toxic elements may cause serious harm to human health, especially to children, due to their larger absorption rate and intake and higher frequency of engaging in hand-to-mouth activities [56].
The PMF modeling results were compared with previous studies. In this work, factors could be interpreted as a mixture of soil dust, vehicle emissions, industry, biomass burning, waste burning, secondary formation, and coal combustion, which were similar to previous studies in Beijing [25,57] and Baoding [30], and the contributions were found to be reasonable.

Cluster and Source Regions
Air mass back trajectories were grouped into five clusters, with clustering pathways illustrated in Figure 7a. Cluster 1, 2, 3, 4, and 5 accounted for 9.85%, 46.21%, 17.42%, 14.39%, and 12.12% of the total trajectories, respectively. Cluster 1 gathered air mass trajectories from Inner Mongolia, across the northeast of Shanxi Province, the northwest of Hebei Province, the southwest of Beijing, and the west of Tianjin before arriving at Gucheng. Cluster 2 was characterized by trajectories originating from the west of Mongolia, passing through Inner Mongolia and the north of Shanxi Province, and then reaching Gucheng. This represented fast-moving trajectories and long-range transport. Cluster 3 began in Mongolia and traveled through Inner Mongolia, the northeast of Hebei Province, and Tianjin before arriving at Gucheng. Cluster 4 originated from the west of Inner Mongolia, passed though Gansu, Ningxia, Shaanxi, and Shanxi Province, and finally reached Gucheng. Cluster 5 originated from the east of Hebei Province and passed through the northwest of Shandong Province to Gucheng, indicating low speed and short-range transport. Air masses from Mongolia and Inner Mongolia may deliver dust from the Gobi Desert, then mix with industrial emissions on their pathways though industrial areas in Tianjin and Hebei Province, or carry combustion pollutants from Shanxi Province, which is an important coal consumption and power generation region in China.
The weighted PSCF (WPSCF) map and the weighted CWT (WCWT) map of PM 2.5 is shown in Figure 7b,c. The central and southern Hebei Province were identified as the primary source regions, with the highest WPSCF and WCWT values. The southwest of Beijing, the west of Tianjin, the northwest of Shandong Province, the northeast of Shanxi Province, and Inner Mongolia were potentially the main source areas, with WPSCF values greater than 0.5 and WCWT values greater than 110 µg m −3 .
Therefore, local emissions, as well as regional and long-range transport processes, had a substantial role in PM 2.5 pollution in the studied area.

Conclusions
The NCP is one of the most polluted regions in China. This study investigated the concentrations of trace elements in PM 2.5 collected at Gucheng, a rural site in the NCP, from 22 October to 23 November 2016. Cl, S, and K were the most abundant elements. Among noncrustal trace metal elements, Zn had the highest concentration, followed by Sb and Pb. SEM was employed for the morphological observation of PM 2.5 . The shape of the particles was predominantly spherical, chain-like, and irregular. An EF analysis indicated that most of the trace elements in PM 2.5 were related to anthropogenic sources. The PMF model identified six sources, namely soil dust, industry, secondary formation, vehicle emissions, biomass and waste burning, and coal combustion (in descending order). Cluster, PSCF, and CWT analyses revealed that local emissions from Hebei Province, as well as regional transport from Beijing, Tianjin, Shandong, and Shanxi Province, and long-range transport from Inner Mongolia, contributed mainly to PM 2.5 . Therefore, strict and coordinated regional air pollution prevention and control measures should be adopted, such as strengthening the regulation of dust, industrial emissions, high-emission vehicles, bulk coal and waste combustion, adjusting the industrial structure, and promoting the use of clean energy. The results of this study will help researchers better understand the characteristics and sources of trace elements in PM 2.5 in the studied area and provide a reference for emission reduction measures, which are important for improving air quality and reducing the effects of air pollution on human health.