Source Apportionment of Heavy Metals in Wet Deposition in a Typical Industry City Based on Multiple Models

: Handan city as a transportation hub in North China, air quality ranks at the bottom all year round, causing environmental pollution that has aroused widespread concern. In order to explore the pollution characteristics and main sources of heavy metals in atmospheric wet deposition in the city, and comparison of the applicability of multiple traceability models, a total of 76 wet deposition samples were collected in the three functional areas from December 2017 to November 2019 by a dry and wet deposition automatic sampler. Concentrations of Cu, Zn, Cr, Ni, Pb, and As were determined and enrichment factors were calculated. Sources of these heavy metals were apportioned by PMF, Unmix, and APCS-MLR models, and analyzed using a backward trajectory analysis model. The results showed that the concentrations of these heavy metals in the atmospheric wet deposition were in order of Zn, Cr, Pb, Cu, Ni, and As, and their mean concentrations were 29.53, 14.11, 9.18, 7.03, 6.41, and 1.21 µ g · L − 1 , respectively. According to the results of EF, the studied heavy metals were mainly affected by anthropogenic activities. The source apportionment results showed that heavy metal pollution in the wet deposition was mainly affected by trafﬁc sources, industrial sources, and coal combustion sources, and PMF identiﬁed an additional source factor: metal smelting source. By comparing the relevant parameters of the source apportionment results of the three models, the APCS-MLR model has better accuracy results than PMF and Unmix models. The analysis of the backward trajectory of the air mass showed that the wet deposition of Handan in the study time was mainly from the southwest direction, accounting for 54.35%. In the future, more evaluation methods and models will be used to compare and analyze the different application scenarios and parameter selection requirements in order to contribute to urban atmospheric environmental pollution prevention and control work.


Introduction
In recent years, heavy metal pollution has received widespread attention all over the world, due to its toxicity, difficult degradation, and strong bioaccumulation [1,2]. Heavy metals in the environment seriously threaten the stability and balance of ecosystems, so their spatial distribution characteristics and main sources of pollution have also become a hot topic at present. Studies have shown that the primary source and carrier of pollutants are atmospheric deposition [3][4][5]. Wet deposition refers to the process in which particulate and gaseous pollutants in the atmosphere settle to the surface after weather processes such as wet deposition, snowfall, hail, and fog, and is one of the main ways these air pollutants are moved [6]. Wet deposition carries heavy metal pollutants into the surface environment, which will have a serious impact on the ecosystem health of land and water bodies [7]. Heavy metals are difficult to degrade and accumulate continuously after deposition on the The major objectives of this research were to: (a) the concentrations of Cr, Cu, N Pb and As in wet deposition of Handan were determined, calculate enrichment fac investigate their pollution enrichment degree; (b) identify the sources of heavy met Unmix, PMF, and APCS-MLR models, and the suitability of the three models was pared; (c) further analyze the region and source of influencing pollutant concentrat backward trajectory model.

Study Area and Sampling
Handan is located in Southern Hebei Province, China, bordering Shanxi Pro and Henan Province, with Taihang Mountain in the west and North China Plain east. The total urban area is 12,073.8 km 2 , and the resident population is 9,549,700. T into account the urban function division of Handan City, three sampling stations laid in the main urban area City (KY, DW, and HG). KY was located in Hebei Univ of Engineering (N 36.57°, E 114.50°), representing the high parish area; DW was lo in the East Sewage Treatment Plant in Handan (N 36.63°, E 114.53°), representing t ing area; HG was located in Handan Iron and Steel Co, Ltd. (N 36.60°, E 114.44°) of Group Hansteel Company, representing the industrial area. The specific distribut sampling points was shown in Figure 1.  All sampling points were set up on an open building roof, with a height of about 12 to 15 m. There were no tall shielding objects and local pollution sources around. Dry and wet deposition samples were monitored and collected monthly from December 2017 to November 2019 by dry-wet deposition automatic collector (TE-78-100X, Tisch. Cleveland, OH, USA). The sampling barrel was soaked in 10% HCl solution for 24 h, washed with deionized water, dried, and placed at the sampling point for collection [8,12]. In this study, 76 wet deposition samples were collected after wet deposition.

Chemistry Analysis
After the wet deposition samples were collected and sent to the laboratory, the falling leaves, insects, and other foreign bodies were removed, the particulate matter in the samples was filtrated with a 0.45 µm microporous membrane, the filtrate was loaded into a dry and clean polyethylene plastic bottle, a part of the samples were taken to determine the pH (PHS-3C) and conductivity (DDS-11A) on the spot, and the remaining samples were placed in a refrigerator at 4 • C for preservation. Wet deposition samples adopted water quality-digestion of total metals-nitric acid digestion method (HJ677-2013) standard. Take 100 mL water sample into a beaker, add 5 mL concentrated HNO 3 , and mix it with sample. Cover the beaker and put it on the electric heating plate (EH20A plus, LabTech. Beijing, China) and kept at 95 ± 5 • C until the sample had completely reacted with HNO 3 , then removed the cover and evaporated the solution to about 5 mL. After the reaction, the resulting solution was cooled down to room temperature, added 3 mL H 2 O 2 , and continued heating until the general appearance of the solution did not change. Transfer the solution to a volumetric bottle of 100 mL, add distilled water to the scale mark, and filter with 0.45 µm filter film. Cu, Zn, Cr, Ni, and Pb were analyzed using an inductively coupled plasma atomic emission spectrometer (ICP-OES) (Agilent 5100, Agilent Technologies. Santa Clara, CA, USA) [16]. As was determined by atomic fluorescence spectrometer (AFS-230E. Runqee, Beijing, China). In the wet deposition water-soluble blank ion test, rinsed sample bottle with 50 mL deionized water, preserved and determined according to the above test method. Prepare 2 blank tests for each batch of digested samples and took the average [8].

Quality Control and Quality Assurance
All data were subjected to strict quality assurance and control procedures. The detection limits (3× the standard deviation of the blanks) were as follow: Cr (0.0002 µg·L −1 ), Ni (0.004 µg·L −1 ), Cu (0.0002 µg·L −1 ), Zn (0.0003 µg·L −1 ), As (0.0006 µg·L −1 ) and Pb (0.008 µg·L −1 ). Calibrating the standard curve for every 10 samples analyzed, and the analytical error was <5%, the recovery rates for the elements were between 85 and 105%, the correlation coefficients (R 2 ) of the calibration curves were ≥0.999 for all elements, each sample was operated in parallel three times. LODs of heavy metals was determined by triple standard deviations from seven replicates at the lowest level of the standard solution (1 to 20 ug/L), and it showed the linearity16. When parallel samples were collected, the number of parallel samples was controlled at 10% of the total number of samples, and the deviation between the determination results of parallel samples and samples was less than 10% [10]. The concentration of each metal in the blank sample determined was lower than the detection limit.

Enrichment Factor
The enrichment factor (EF) was a method used to express the enrichment degree of elements, and judge and evaluate the source of elements in particulate matter (natural source and anthropogenic source). We generally choose elements such as Ti and Al which have little anthropogenic pollution impact and stable geochemical properties as reference elements, Al was chosen as the reference element in this study. The calculation formula is as follows [31]: where C i is the concentration of element i (µg·L −1 ), C Al is the concentration of Al. The element concentration average value of layer A soil in Hebei Province was used as the background value [32]. The enrichment factor was calculated for each sample collected and averaged separately based on the number of collections within each month. Pollution degrees can be divided into three categories according to the EF values, and their specific grading criteria and pollution source are presented in Table 1.

Unmix Model
The Unmix model is one of the source apportionment models recommended by the U.S. Environmental Protection Agency (EPA), which is an acceptor model developed on the basis of PCA. The Unmix model is simple to perform, does not require complex adjustments to parameters, and does not require knowing the source component spectrum [33]. The basic formula is as follows [34,35]: where C ij is the concentration of the j th specie in the i th sample; F jk is the concentration of the j th specie in source k (1 . . . . m) mass fraction; S ik is the total amount of source k in the i th sample, representing the contribution rate of the source; E is the standard deviation of each source composition. It is generally believed that the use of this model must meet the following two conditions simultaneously: 1 Min R 2 > 0.8, R 2 is the correlation coefficient of the calculated value and the measured value of specie concentration; 2 Min S/N > 2, S/N represents the signal to noise ratio, the model analysis results are credible.

PMF Model
The PMF model was first proposed by Paatero in 1994 [27]. The sample concentration in the model is expressed by factor contribution, factor distribution and residual error. The principle is based on the weighted least squares method for bounded and iterative calculation. The matrix X is continuously decomposed to obtain the matrices G and F, and the new matrix is optimized and solved. When the objective function Q is the smallest, it is the optimal solution [26,28]. In the operation process, the model needs to repeatedly calculate and select the optimal result according to the signal and noise ratio (S/N) and Q value. The calculation formula is as follows: where X ij is the measurement matrix of the j th heavy metal in sample i; G ik is the contribution matrix of the k th source factor in sample i; F kj is the distribution matrix of the j th heavy metal in the k th source factor; E ij is the residual value of the j th heavy metal in sample i; Q is the objective function defined by the PMF model; U ij is the uncertainty value of the j th element in the i th sample. Running the PMF model requires the import of uncertainty data, and if the heavy metal concentration is less than the species-specific method detection limit (MDL), then the uncertainty value is equal to 5 6 of the MDL; if heavy metal concentration exceeds MDL value, the uncertainty value is calculated by the Equation (5): where C is measured heavy metal concentration and RSD is the relative standard deviation of the heavy metal concentration.

APCS-MLR Model
Absolute principal component score-multiple linear regression (APCS-MLR) receptor model was proposed by Thurston and Spengler (1985) [26]. The purpose of multiple linear regression is to further determine the contribution rate of different pollution sources on the basis of source identification. To standardize the concentration of data processing, the basic formula is as follows: Z ij is the standardization concentration value (dimensionless); C ij is the concentration of element i; C i and σ i are the mean concentration and standard deviation of element i, respectively.
Introduced the artificial sample with concentration of 0, calculate the factor fraction of this sample, and calculate the formula as: The standardized factor scores of heavy metal concentration samples can be obtained by principal component analysis (PCA), and the absolute principal component scores (APCS) of each heavy metal element can be obtained by subtracting the factor score of the sample from the factor score of the zero-concentration human sample corresponding to heavy metals. Then taking APCS as independent variable and the concentration of heavy metal element as dependent variable, multiple linear regression was performed. The calculation formula is as follows: where C i is the element concentration; b i0 is the constant term of multiple linear regression; n is the total source number; b pi is the regression coefficient of multiple linear regression; APCS p is the absolute principal component factor score; b pi × APCS p is the contribution rate of source p to C i concentration, and the average value of b pi × APCS p of all samples is the contribution rate of this source.

Heavy Metal Concentration Level in Wet Deposition
In this study, the precipitation during the wet deposition samples collection period is shown in Figure 2  According to the results of chemical analysis, the average concentrations of six metal elements in wet deposition from December 2017 to November 2019 in three functional areas of Handan were statistically analyzed. It can be seen from Table 2 that the average concentrations of elements from high to low were Zn, Cr, Pb, Cu, Ni, and As. Among them, The average concentration of Zn was 29.52 µg·L −1 , which was much higher than other elements; the mean value of As was only 1.21 µg·L −1 . Elements also vary in concentration in different functional regions. The concentrations of Cr, Ni, and Pb were highest in high parishes, As was highest in living areas, and Cu and Zn were highest in industrial areas. In the three functional regions, Zn and Cr were the highest. Among them, Zn is a metallurgical dust characteristic element, and Zn is also produced by the wear of rubber tires [36,37], while Cr is closely related to industrial production [38], indicating that industrial and traffic sources contribute the most to heavy metal particulate matter in atmospheric wet Atmosphere 2022, 13, 1716 7 of 16 deposition in Handan City. This city is also dominated by heavy industry, and there are industrial enterprises such as machinery factories and smelters around [39]; due to its special geographical location of the city, it is located at the junction of the four provinces of Hebei, Shandong, and Henan, and its developed transportation may also be the reason for the high concentration of elements in wet deposition in this area [40]. According to the results of chemical analysis, the average concentrations of six metal elements in wet deposition from December 2017 to November 2019 in three functional areas of Handan were statistically analyzed. It can be seen from Table 2 that the average concentrations of elements from high to low were Zn, Cr, Pb, Cu, Ni, and As. Among them, The average concentration of Zn was 29.52 μg·L −1 , which was much higher than other elements; the mean value of As was only 1.21 μg·L −1 . Elements also vary in concentration in different functional regions. The concentrations of Cr, Ni, and Pb were highest in high parishes, As was highest in living areas, and Cu and Zn were highest in industrial areas. In the three functional regions, Zn and Cr were the highest. Among them, Zn is a metallurgical dust characteristic element, and Zn is also produced by the wear of rubber tires [36,37], while Cr is closely related to industrial production [38], indicating that industrial and traffic sources contribute the most to heavy metal particulate matter in atmospheric wet deposition in Handan City. This city is also dominated by heavy industry, and there are industrial enterprises such as machinery factories and smelters around [39]; due to its special geographical location of the city, it is located at the junction of the four provinces of Hebei, Shandong, and Henan, and its developed transportation may also be the reason for the high concentration of elements in wet deposition in this area [40].

Enrichment Factor
Because there had been few rainfalls and snowfalls from October to next February in Handan, no effective wet deposition samples were collected, so only samples from March to September were selected for comparative analysis. The EF value of each heavy metal was calculated according to Equation (2), as shown in Figure 3. There were differences in the enrichment levels between different heavy metals during different months in the study area, while their EF values were all above 10. Most of the EF values of Cr, Cu, Pb, Zn, and Ni were above 100, and the EF values of As were mostly between 10 and 100, indicating that the study area was mainly affected by anthropogenic sources. Some of the EF values of Cr, Cu, Pb, Zn, and Ni were greater than 1000, showing a high enrichment, indicating that they were seriously affected by anthropogenic activities. The above results were consistent with other scholars' research reports, indicating that the concentration of six elements in wet deposition in this study was seriously affected by human activities [6,41].

Results of the PMF Model
The contribution load and rate of each element source analyzed by the PMF was shown in Figure 4. The factor probability value obtained by PMF model deco tion was an important basis for the quantitative allocation of heavy metals to sour In this study, the PMF model completed 20 probability matrix decomposition u factors number was 4, and the QRobust/QTrue value tended to converge and reached mum value (QTrue was the goodness-of-fit parameter calculated including all poin The contribution load and rate of each element source analyzed by the PMF model was shown in Figure 4. The factor probability value obtained by PMF model decomposition was an important basis for the quantitative allocation of heavy metals to sources [6]. In this study, the PMF model completed 20 probability matrix decomposition until the factors number was 4, and the Q Robust /Q True value tended to converge and reached a minimum value (Q True was the goodness-of-fit parameter calculated including all points, Q Roubust was the goodness-of-fit parameter calculated excluding points not fit by the model), while the residual size was basically located between −3 and 3, S/N ≥ 8.9, R 2 ≥ 0.54, which met the optimal result requirements, and the model analysis effect was good.  The concentration load of Zn in factor 1 was the highest, and its contribution rate was as high as 83.83%, and the contribution rate of Ni was 39.33%. Zn was mainly derived from smelting and mining, was a marker element of industrial pollution, and was also widely used in the production of tires in the rubber products industry; the corrosion-resistant alloy that constitutes the shell of motor vehicles contains Ni, and it was also associated with coal burning [42,43]. Therefore, factor 1 can be identified as a mixed source of industry and transportation.
As had the highest concentration load in factor 2, and its contribution rate was 96%. Ni and Pb also had a certain concentration load on factor 2, but the contribution rates were below 10%. Therefore, As was the main element of factor 2, and it was also a characteristic and tracer element of coal-burning emissions [44], Handan is a typical heavy industrial city, with coal consumption of 4.01 × 10 7 and 4.10 × 10 7 t in 2018 and 2019, respectively [45,46]. Therefore, it was considered that factor 2 was mainly the source of coal combustion.
Pb had the highest concentration load in factor 3, with a contribution rate of 79.03%, and the other elements had a low concentration load. Pb was generally derived from traffic emissions and industrial activities, but with the use of unleaded gasoline, the contribution of traffic emissions to Pb in the atmosphere had decreased, and industrial emissions had gradually become the main source of atmospheric dust Pb [45]. Therefore, factor 3 was considered to be of industrial origin.
The concentration loading of Cr and Cu in factor 4 was much higher than those of other elements, and the contribution rates were as high as 56.87% and 64.45%. Cu was a marker element in the nonferrous metallurgy industry, while Cr had a close association with metal smelting in the atmospheric environment [46,47]. Factor 4 was therefore considered to be a metal smelting source.

Results of the Unmix Model
The contribution rate of different sources by Unmix model is shown in Figure 5. The The concentration load of Zn in factor 1 was the highest, and its contribution rate was as high as 83.83%, and the contribution rate of Ni was 39.33%. Zn was mainly derived from smelting and mining, was a marker element of industrial pollution, and was also widely used in the production of tires in the rubber products industry; the corrosion-resistant alloy that constitutes the shell of motor vehicles contains Ni, and it was also associated with coal burning [42,43]. Therefore, factor 1 can be identified as a mixed source of industry and transportation.
As had the highest concentration load in factor 2, and its contribution rate was 96%. Ni and Pb also had a certain concentration load on factor 2, but the contribution rates were below 10%. Therefore, As was the main element of factor 2, and it was also a characteristic and tracer element of coal-burning emissions [44], Handan is a typical heavy industrial city, with coal consumption of 4.01 × 10 7 and 4.10 × 10 7 t in 2018 and 2019, respectively [45,46]. Therefore, it was considered that factor 2 was mainly the source of coal combustion.
Pb had the highest concentration load in factor 3, with a contribution rate of 79.03%, and the other elements had a low concentration load. Pb was generally derived from traffic emissions and industrial activities, but with the use of unleaded gasoline, the contribution of traffic emissions to Pb in the atmosphere had decreased, and industrial emissions had gradually become the main source of atmospheric dust Pb [45]. Therefore, factor 3 was considered to be of industrial origin.
The concentration loading of Cr and Cu in factor 4 was much higher than those of other elements, and the contribution rates were as high as 56.87% and 64.45%. Cu was a marker element in the nonferrous metallurgy industry, while Cr had a close association with metal smelting in the atmospheric environment [46,47]. Factor 4 was therefore considered to be a metal smelting source.

Results of the Unmix Model
The contribution rate of different sources by Unmix model is shown in Figure 5. The numerical value for the solution's diagnostic indicators (Min R 2 = 0.85 and Min S/N = 2.14) is consistent with the recommendations, and the analytical results obtained from these three sources were credible.
1, which were typical industrial sources and were also closely related to metal smelting, and large state-owned smelters were located in the study area [48,49], so source 1 was defined as the industrial source. Source 2 accounted for the largest proportion of interpretation of heavy metals, Cu, Ni, and Zn were representative elements, accounting for 51.65%, 57.14%, and 87.93%, respectively. Ni and Zn were a significant correlation with transportation, and Cu was used as an additive for automotive braking and heat dissipation systems [50], so source 2 was defined as a traffic source. The representative element in source 3 was As, accounting for 80.56% and consisting of of the factor 2 elucidation result of PMF. Cr and Cu also had a certain concentration load in source 3, but the contribution rate was below 50%, much lower than that of As, so source 3 was judged as coalfired.

Results of the APCS-MLR Model
In this study, six heavy metals were tested for suitability (Kaiser-Meyer-Olkin, KMO), their KMO was 0.568 > 0.5, and Bartlett's sphericity test p value was ≤0.001, indicating that there was a good correlation between the elements, and they were suitable for factor analysis. According to the APCS-MLR receptor model, the scores of the three principal factors in factor analysis were converted into absolute principal factor scores, and then analyzed by multiple linear regression, respectively, to obtain the multiple linear regression equations on the six heavy metals, and the fitted values of the concentrations of each element could be obtained from the equations and compared with the measured values, the closer this value was to 1 indicating that the better the fitting degree of multiple linear regression [49,50]. The factor analysis results of heavy metals are shown in Table 3. After Varimax rotation, selecting the factor with characteristic value greater than 1, the three factors cumulatively explained 78.22%. It was inferred that there were three main sources of heavy metal pollution.  The heavy metals accounted for 30.90%, 55.80%, and 13.30% in source 1, source 2, and source 3, respectively. Pb and Cr accounted for high concentrations of load on source 1, which were typical industrial sources and were also closely related to metal smelting, and large state-owned smelters were located in the study area [48,49], so source 1 was defined as the industrial source. Source 2 accounted for the largest proportion of interpretation of heavy metals, Cu, Ni, and Zn were representative elements, accounting for 51.65%, 57.14%, and 87.93%, respectively. Ni and Zn were a significant correlation with transportation, and Cu was used as an additive for automotive braking and heat dissipation systems [50], so source 2 was defined as a traffic source. The representative element in source 3 was As, accounting for 80.56% and consisting of of the factor 2 elucidation result of PMF. Cr and Cu also had a certain concentration load in source 3, but the contribution rate was below 50%, much lower than that of As, so source 3 was judged as coal-fired.

Results of the APCS-MLR Model
In this study, six heavy metals were tested for suitability (Kaiser-Meyer-Olkin, KMO), their KMO was 0.568 > 0.5, and Bartlett's sphericity test p value was ≤0.001, indicating that there was a good correlation between the elements, and they were suitable for factor analysis. According to the APCS-MLR receptor model, the scores of the three principal factors in factor analysis were converted into absolute principal factor scores, and then analyzed by multiple linear regression, respectively, to obtain the multiple linear regression equations on the six heavy metals, and the fitted values of the concentrations of each element could be obtained from the equations and compared with the measured values, the closer this value was to 1 indicating that the better the fitting degree of multiple linear regression [49,50]. The factor analysis results of heavy metals are shown in Table 3. After Varimax rotation, selecting the factor with characteristic value greater than 1, the three factors cumulatively explained 78.22%. It was inferred that there were three main sources of heavy metal pollution. According to the results of reference PMF and Unmix models, these three factors represent traffic source, coal-burning source, and industrial source, respectively. The proportion of different pollution factors could be seen in Figure 6. Factor 3 (44.18%) accounts for the highest proportion, followed by factor 1 (41.75%) and factor 2 (14.07%). After the Varimax rotation method, the coefficients of the factor matrix (Table 3), and factor 1 (traffic source) explained 30.01% of the total factor score, with a large load on Cr, Ni, and Pb, and the load decreased in turn; factor 2 (coal-fired source) explained 27.60% of the whole, with high loading on As and Cu, and the load decreased in turn; factor 3 (industrial source) explained 20.61% of the whole, with a large load on Zn, and the load was much higher than other elements.  According to the results of reference PMF and Unmix models, these three factors represent traffic source, coal-burning source, and industrial source, respectively. The proportion of different pollution factors could be seen in Figure 6. Factor 3 (44.18%) accounts for the highest proportion, followed by factor 1 (41.75%) and factor 2 (14.07%). After the Varimax rotation method, the coefficients of the factor matrix (Table 3), and factor 1 (traffic source) explained 30.01% of the total factor score, with a large load on Cr, Ni, and Pb, and the load decreased in turn; factor 2 (coal-fired source) explained 27.60% of the whole, with high loading on As and Cu, and the load decreased in turn; factor 3 (industrial source) explained 20.61% of the whole, with a large load on Zn, and the load was much higher than other elements.

Comparison of Three Receptor Models
In this study, we compared the fitting degree of the three models as well as R 2 to determine their advantages and disadvantages. The fitting degree of the predictive value and the measured value (P/M) were considered the fitting degree of the models. The regression coefficient and residual value were the common parameters used to estimate the fitting goodness of the model. The closer the R 2 value and slope of the former were to 1, the closer the intercept was to 0, and the better the fitting effect; while the closer the mean value of the latter was to 0, the smaller the range was, the better the fitting effect was. The relevant results of the model were shown in the following Table 4.

Comparison of Three Receptor Models
In this study, we compared the fitting degree of the three models as well as R 2 to determine their advantages and disadvantages. The fitting degree of the predictive value and the measured value (P/M) were considered the fitting degree of the models. The regression coefficient and residual value were the common parameters used to estimate the fitting goodness of the model. The closer the R 2 value and slope of the former were to 1, the closer the intercept was to 0, and the better the fitting effect; while the closer the mean value of the latter was to 0, the smaller the range was, the better the fitting effect was. The relevant results of the model were shown in the following Table 4.
The fitting effect of Pb in the Unmix model was far lower than in other models (R 2 = 0.35 and P/M = 0.60), the R 2 of Pb was as high as 0.77 and 0.78 in PMF and APCS-MLR, respectively. The PMF and APCS-MLR models were therefore more persuasive for the source apportionment results of Pb. In the PMF model, the fitting effect of each heavy metal was good, the R 2 of the remaining elements except Ni was more than 0.60, and the P/M of all elements was more than 0.70, of which the R 2 value of As was as high as 0.99. In the APCS-MLR model, the R 2 of the six heavy metals was above 0.70, both P/M values were 1.00, indicating that the operation range error of the model was small, the predictive value was basically consistent with the actual value, and the analysis results of each heavy metal in the model have good accuracy. In summary, although the PMF, Unmix, and APCS-MLR receptor models were operated with different algorithms and mechanisms, the results were more consistent in the study area, and the analysis results could validate each other. The three source apportionment models were used in this study concluded that metal elements in wet deposition in Handan affected by anthropogenic impacts such as industry, traffic, and coal combustion, and were basically consistent with the soil, PM 2.5 , and dry deposition results studied by Cai et al. [51][52][53][54]. Compared with the results of source apportionment in this study, the APCS-MLR model had a better fitting degree, the fitting effect was more prominent when explaining individual elements, which could more accurately explain the source category and contribution of heavy metals in the study area, this model is more suitable for the study of wet deposition.

Air Pellet Backward Trace
Air mass transport could bring pollutants or clean air to the local area, affecting the local wet deposition and wet deposition chemical composition [55]. In this study, a backward trajectory model (Meteoinfo) was used to analyze the trajectory of wet deposition clouds at 500 m height for 48 h in 25 wet depositions fields in Handan during 2018 to 2019 ( Figure 7). Combined with the backward trajectory of the air mass, it could be seen that the wet deposition air mass in Handan in this year mainly came from the southwest direction, the wet deposition times reached 14 times, and the wet deposition reached 1207.57 mm, accounting for 54.35% of the total; the wet deposition brought about by the northwest and northeast air masses for Handan was relatively small, 627.89 mm and 386.38 mm, respectively, accounting for 28.26% and 17.39% of the total wet deposition. The clouds in the northeast direction mainly came from Inner Mongolia, through Beijing, Langfang, and other places; the two clouds in the northwest direction pass through Shanxi Province, and pass through Taiyuan, Jinzhong and other areas; the southwest clouds mainly came from Changzhi City, Shanxi Province. Through the analysis of the air mass trajectory in the southwest direction, it was found that the air mass operation trajectory operation path was short, indicating that Handan was mainly controlled by the short-distance transport air mass in the inner and southern regions of Shanxi Province. Combined with the topographic analysis of Shanxi Province, the Taihang Mountains are distributed in the east, most of which were above 1200 m above sea levels, and there was a block of Luliang Mountains in the west, most of which were above 1000 m above sea level, and the city was an energy heavy chemical base in China, with coal power, coking, metallurgy as the main industries, which contributes greatly to the metal Zn, Cr and As in the air mass; there were large thermal power plants in Jinzhong and Changzhi, such as Yangcheng Power Plant and Zhangze Power Plant, which contribute greatly to the metal Ni, Cu and Pb in the air mass [56].
Mountains in the west, most of which were above 1000 m above sea level, and the city was an energy heavy chemical base in China, with coal power, coking, metallurgy as the main industries, which contributes greatly to the metal Zn, Cr and As in the air mass; there were large thermal power plants in Jinzhong and Changzhi, such as Yangcheng Power Plant and Zhangze Power Plant, which contribute greatly to the metal Ni, Cu and Pb in the air mass [56].

Conclusions
(1) The mean concentrations of Zn, Cr, Pb, Cu, Ni, and As in wet deposition were 29. 53, 14.11, 9.18, 7.03, 6.41, and 1.21 μg·L −1 , respectively. The concentration of elements in different functional areas varied, indicating that human activity was one of the main reasons affecting the concentration of different heavy metals in wet deposition. (2) In this study, the EFs of the six heavy metals were above 10, and Pb and Zn were above 100, and the highest values were even above 1000, which were greatly affected by anthropogenic activities and showed moderate enrichment and high enrichment characteristics. In addition, according to the backward trajectory analysis, the wet deposition in Handan was mainly affected by the gas mass from the southwest direction. The coking and power industrial areas through which the gas mass trajectory passed had a greater impact on the wet deposition of heavy metals, the joint prevention and control of the area should be the main improvement measure for the heavy metal pollution of wet deposition in the future, while requiring all regions to strictly implement the emission limit requirements of air pollutants. (3) By comparing the analytical results of PMF, Unmix, and APCS-MLR multiple receptor models, the analytical results of the three models confirm each other and have a good consistency. In particular, the APCS-MLR model had better analysis results, and its fitting degree and applicability were better than other models, which was more suitable for source analysis of heavy metals in wet deposition in this study.
In this study, multiple models were used for comparison, and more suitable models for the case area were selected to ensure the scientificity and rationality of the results, with specific characteristics and applications. In the future study, we will select more elements and models to analyze and compare, so that the source assignment is more objective.

Conclusions
(1) The mean concentrations of Zn, Cr, Pb, Cu, Ni, and As in wet deposition were 29. 53, 14.11, 9.18, 7.03, 6.41, and 1.21 µg·L −1 , respectively. The concentration of elements in different functional areas varied, indicating that human activity was one of the main reasons affecting the concentration of different heavy metals in wet deposition. (2) In this study, the EFs of the six heavy metals were above 10, and Pb and Zn were above 100, and the highest values were even above 1000, which were greatly affected by anthropogenic activities and showed moderate enrichment and high enrichment characteristics. In addition, according to the backward trajectory analysis, the wet deposition in Handan was mainly affected by the gas mass from the southwest direction. The coking and power industrial areas through which the gas mass trajectory passed had a greater impact on the wet deposition of heavy metals, the joint prevention and control of the area should be the main improvement measure for the heavy metal pollution of wet deposition in the future, while requiring all regions to strictly implement the emission limit requirements of air pollutants. (3) By comparing the analytical results of PMF, Unmix, and APCS-MLR multiple receptor models, the analytical results of the three models confirm each other and have a good consistency. In particular, the APCS-MLR model had better analysis results, and its fitting degree and applicability were better than other models, which was more suitable for source analysis of heavy metals in wet deposition in this study.
In this study, multiple models were used for comparison, and more suitable models for the case area were selected to ensure the scientificity and rationality of the results, with specific characteristics and applications. In the future study, we will select more elements and models to analyze and compare, so that the source assignment is more objective.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to participants not having consented to the data being shared.