Assessment of Fluxes and Ecological and Health Risks of Toxic Trace Elements in Atmospheric Deposition from the Baicheng-Songyuan Area, Jilin Province, Northeast China

: A total of 56 atmospheric deposition samples were collected on a yearly basis from the Baicheng-Songyuan areas, Jilin Province, Northeast China. Each sample was subdivided into wet (soluble) and dry (insoluble) fractions, and the concentrations of toxic trace elements including As, Cd, Co, Cr, Cu, Hg, Mn, Mo, Ni, Pb, Se, and Zn in both fractions were analyzed. The bulk fluxes and ecological and human health risk of these toxic trace elements in atmospheric deposition were evaluated. The bulk deposition fluxes of toxic trace elements decreased in the order of Mn > Zn > Cr > Pb > Cu > Ni > As > Co >Mo > Cd > Se > Hg. The fluxes of toxic trace elements in dry fraction accounted for 74.7–99.9% of their bulk deposition fluxes, indicating that the atmospheric deposition of toxic trace elements in the study area is predominantly dry atmospheric deposition. The mean values of the geo-accumulation index (Igeo) and enrichment factor showed moderately heavy contamination with Cd via dry atmospheric deposition, but no contamination from other toxic trace elements; Cd and Hg in dry atmospheric deposition were in the moderate to considerable ecological risk levels, while other trace elements were at low ecological risk levels. The health risk assessment showed that the effects of toxic trace elements in dry atmospheric deposition via three exposure pathways were in the order of ingestion > inhalation > dermal contact for adults and children. The mean values of hazard quotient (HQ) and hazard index (HI) of toxic trace elements via three pathways were less than one, indicating that their non-carcinogenic risks in dry atmospheric deposition may be low or negligible for adults and children. The mean values of carcinogenic risk (CR) and total carcinogenic risk (TCR) of As and Cr via the three pathways for adults and children were between 10 − 6 and 10 − 4 , indicating that the carcinogenic risk levels of As and Cr were tolerable or acceptable, and the mean TCR value of Cd through the three pathways for adult and children was less than 10 − 6 , implying that the carcinogenic risk level of Cd was negligible. Mn, Ni, Cr, and Co in dry atmospheric deposition were mainly contributed from the crustal sources, while As, Cd, Cu, Hg, Mo, Pb, Se, and Zn in dry atmospheric deposition were derived from both crustal and anthropogenic sources. The results obtained in this study advocate the necessity for monitoring atmospheric deposition in some rural areas, and also provide a scientific basis for controlling contamination posed by toxic trace elements in dry atmospheric deposition


Introduction
Atmospheric deposition in the forms of precipitation, particles, and aerosols is a significant pathway of pollutants entering natural ecosystems such as soils, wetlands, and lakes.Pollution by toxic trace elements, such as As, Cd, Cr, Cu, Hg, Ni, Pb, and Zn, in dry and wet atmospheric deposition processes into the ecosystems are of particular interest due to their negative effects on ecosystem safety and human health [1,2].Therefore, atmospheric with a short, windy, and dry spring, warm and rainy summer, sunshiny and cool autumn, and long and cold winter, with an annual mean temperature of about 3-6 • C, a yearly average rainfall of 350-550 mm and potential evaporation of 1600-2000 mm.The climate of the study area reflects a transitional zone between the semi-arid and arid climates in the west and the semi-humid climates in the east.The topography of the area is high in the northwest and low in the middle.The prevailing wind in the area is northwest to west in the winter, and southwest in other seasons, with a mean wind speed of 3.4-4.4m/s.The sand and dust storms mostly appear from March to June each year, lasting about 120-130 days, and a haze resulting from crop stalk incineration often occurs in the winter and spring.Both sandstorms and haze create serious atmospheric pollution in Changchun, the capital of Jilin Province, and other downwind cities. Soils in the study area comprise mainly chestnut soil, sub-chernozem soil, chernozem soil, saline-alkali soil, and aeolian sandy soil.The study area is covered by 40-75% grassland.The vegetation includes mainly terraced hills steppe, plain meadow steppe grassland, and plain meadow grassland [39].Some important wetlands, such as the Chagan wetland, the Yueliang wetland and the Xianghai wetland, are distributed in the area, and are considered as ecological barriers and freshwater fishing bases in Jilin Province.

Study Area
The study area is located in the Baicheng-Songyuan area, Jilin Province, Northeast China (Figure 1).The area is characterized by a temperate continental monsoon climate, with a short, windy, and dry spring, warm and rainy summer, sunshiny and cool autumn, and long and cold winter, with an annual mean temperature of about 3-6 °C, a yearly average rainfall of 350-550 mm and potential evaporation of 1600-2000 mm.The climate of the study area reflects a transitional zone between the semi-arid and arid climates in the west and the semi-humid climates in the east.The topography of the area is high in the northwest and low in the middle.The prevailing wind in the area is northwest to west in the winter, and southwest in other seasons, with a mean wind speed of 3.4-4.4m/s.The sand and dust storms mostly appear from March to June each year, lasting about 120-130 days, and a haze resulting from crop stalk incineration often occurs in the winter and spring.Both sandstorms and haze create serious atmospheric pollution in Changchun, the capital of Jilin Province, and other downwind cities. Soils in the study area comprise mainly chestnut soil, sub-chernozem soil, chernozem soil, saline-alkali soil, and aeolian sandy soil.The study area is covered by 40-75% grassland.The vegetation includes mainly terraced hills steppe, plain meadow steppe grassland, and plain meadow grassland [39].Some important wetlands, such as the Chagan wetland, the Yueliang wetland and the Xianghai wetland, are distributed in the area, and are considered as ecological barriers and freshwater fishing bases in Jilin Province.

Sample Collection and Pretreatment
A total of 56 atmospheric deposition sampling sites were chosen in the study area (Figure 1).Atmospheric deposition samples were taken using a self-designed collection container followed the protocol set by China Geological Survey [40] from July 2017 to June 2018.The sample collector was a cylindric porcelain container with an inner diameter of 20 cm, a height of 45 cm, and an available section area of 314 cm 2 .Prior to the sample collection, preparatory cleaning methods were used on the sample containers as follows: the containers were presoaked with 10% HCl solution for 24 h, then rinsed with deionized water, and air-dried at room temperature (20-25 °C).The cleaned and dried containers were transported to the sampling sites where they were positioned on the rooves of the farmers' houses at sampling sites (Figure 1), at a height of about 3-4 m above the ground,

Sample Collection and Pretreatment
A total of 56 atmospheric deposition sampling sites were chosen in the study area (Figure 1).Atmospheric deposition samples were taken using a self-designed collection container followed the protocol set by China Geological Survey [40] from July 2017 to June 2018.The sample collector was a cylindric porcelain container with an inner diameter of 20 cm, a height of 45 cm, and an available section area of 314 cm 2 .Prior to the sample collection, preparatory cleaning methods were used on the sample containers as follows: the containers were presoaked with 10% HCl solution for 24 h, then rinsed with deionized water, and air-dried at room temperature (20-25 • C).The cleaned and dried containers were transported to the sampling sites where they were positioned on the rooves of the farmers' houses at sampling sites (Figure 1), at a height of about 3-4 m above the ground, and filled with 3000 mL deionized water to avoid re-suspension of collected soil particles.Sample sites were not near major roads, high buildings, big trees, and industrial pollution sources.
The collected atmospheric deposition samples were transported into the laboratory, where the dry and wet fractions of samples were separated by pumping filtration with 0.45 µm cellulose polyester ester filters (Millipore, 47 mm i.d.).
The wet fraction of a sample was the filtrates passing through a 0.45 µm filter.The volume and weight of the wet fraction through the filter were measured; then, it was transferred to a 500 mL plastic container, and stirred uniformly after adding 5 mL HNO 3 (1:1, v/v), serving as a stock solution for analyzing multiple elements.
The dry fraction of a sample was the filtrate over the filter, which was air-dried or oven-dried at temperatures lower than 60 • C, then, weighed and ground with an agate mortar and passed through a 200-mesh sieve.The digested method for analyzing Cd, Co, Cr, Cu, Mn, Mo, Ni, Pb, and Zn in the dry fraction of a sample was as follows: weigh 0.1 g of the sieved portion of dry fraction sample with an accuracy of 0.0001 g as the test portion, and transfer it into a 50 mL polytetrafluoroethylene (PTFE) crucible.Moisten it with a small amount of water.Add 2.0 mL concentrated nitric acid, 2.0 mL hydrofluoric acid, 2.0 mL hydrochloric acid, and 1.0 mL perchloric acid separately into the crucible.The solution in the crucible was heated on a hot plate at a temperature of 200 • C until it was digested and evaporated into a white fume.The crucible was removed from the hot plate, and cooled at room temperature; then, 5.0 mL concentrated nitric acid was added to dissolve the residual, and the solution in the crucible was transferred into a 25 mL colorimetric tube with water and diluted with water to reach a 25 mL volume and mixed w.The final solution in the tube was used for testing Cd, Co, Cr, Cu, Mn, Mo, Ni, Pb, and Zn.
The digested method for analyzing As, Hg, and Se in the dry fraction of a sample was as follows: weigh 0.25 g of the sieved portion of the dry fraction sample with an accuracy of 0.0001 g as the test portion, then transfer it into a 50 mL colorimetric tube, and wet it with a small amount of water.Add 10.0 mL mixture of concentrated hydrochloric acid-nitric acid (3:1, v/v) into the tube.Cover the tube with a lid, then mix it.The tube was placed on a boiling water bath and digested for 2 h, and stirred for several times during the water bath period.After removing the tube from the water bath, it was cooled to room temperature, and diluted with water until the solution reached a volume of 50 mL.Then, 10.0 mL solution was transferred into a 25 mL colorimetric tube, 5.0 mL of concentrated hydrochloric acid and 2.5 mL thiourea solution (50 g/L) was added, and it was diluted with water to reach a 25 mL volume, and mixed.Afterwards, we let the tube stand for 30 min, and the final digested solution in the tube was used for testing As, Hg, and Se.

Sample Analysis and Quality Controls
Both the dry and wet fractions of an atmospheric deposition sample were analyzed on the basis of the final digested solutions by pretreatment of the sample.For the wet fraction of a sample, As, Hg, and Se levels were determined with an atomic fluorescence spectrometer (AFS-230E, Langfang, China), and Cd, Co, Cr, Cu, Mn, Mo, Ni, Pb, and Zn were measured with an inductively coupled plasma mass spectrometer (ICP-MS, VG PQ ExCell, Thermo Elemental, Franklin, MA, USA).The detection limits of toxic trace elements were As 0.01 µg/L, Cd 0.01 µg/L, Co 1.0 µg/L, Cr 1.0 µg/L, Cu 1.0 µg/L, Hg 0.001 µg/L, Mn 1.0 µg/L, Mo 1.0 µg/L, Ni 1.0 µg/L, Pb 0.01 µg/L, Se 0.01 µg/L, and Zn 1.0 µg/L, respectively.For the dry fraction of a sample, Co, Cr, Cu, Mn, Ni, Pb, and Zn were measured by inductively coupled plasma atomic emission spectrometer (ICP-AES, IRIS Intrepid II XPS, Thermo Fisher, Waltham, MA, USA).As, Hg, and Se were measured with an atomic fluorescence spectrometer (AFS-230E, Langfang, China).Cd was analyzed by a graphite furnaced atomic absorption spectrometer (GF-AAS, M6, Thermo Elemental, Franklin, MA, USA), and Mo was determined by an inductively coupled plasma mass spectrometer (ICP-MS, VG PQ ExCell, Thermo Elemental, Franklin, MA, USA).The detection limits of toxic trace elements were As 0.6 mg/kg, Cd 0.01 mg/kg, Co 2 mg/kg, Cr 3 mg/kg, Cu 1 mg/kg, Hg 0.005 mg/kg, Mn 10 mg/kg, Mo 0.3 mg/kg, Ni 2 mg/kg, Pb 2 mg/kg, Se 0.01 mg/kg, and Zn 2 mg/kg, respectively.In order to facilitate the evaluation of contamination and sources of toxic trace elements, Al and Fe contents in both fractions of a sample were also analyzed.For the wet fraction of a sample, Al and Fe were measured with ICP-MS (VG PQ ExCell, Thermo Elemental, Franklin, MA, USA), and detection limits were 0.01 mg/L for Al, and 0.05 mg/L for Fe, respectively.For the dry fraction of a sample, Al and Fe were determined with ICP-AES (IRIS Intrepid II XPS, Thermo Fisher, Waltham, MA, USA); detection limits were 0.05% and 0.05%, respectively.
For quality control (QC) and quality assurance (QA), the reagents and chemicals used in the study were of the super grades, and blank control, duplicate samples, and standard reference materials from the Center for Certified Reference Materials, China (GSS-1 and GSS-3) were used to ascertain the accuracy and precision of analytical methods.Blank samples were used throughout all of the experiments as well.The analytical precision for duplicate samples was within ±5%.The precision (RSD) and accuracy (RE) for analyzing standard substances with different analytical methods were −5.01-4.95% and 1.07-5.75%,respectively.The recovery rate of toxic trace elements was in the range of 95.3-105.9%.The concentrations of toxic trace elements in both wet and dry fractions in atmospheric deposition samples were all above their corresponding detection limits.Therefore, the analytical results of toxic trace elements were reliable and could be used for this study.

Geo-Accumulation Index
A geo-accumulation index (Igeo) was used to characterize contamination of toxic trace elements in dry atmospheric deposition, as defined in Equation (1): where C i is the measured concentration of toxic trace element i, B i is a background or reference concentration of same metal i, and 1.5 is a correction constant.In general, the concentration of toxic trace elements in shale average value was often used for B i ; however, the averaged contents of toxic trace elements in the surface soil data from mid-west Jilin Province [41] were applied for B i in this study to facilitate the interpretation of EF of toxic trace elements.Seven contamination categories were classified based on Igeo [42]: Igeo < 0 practically uncontaminated; 0 < Igeo < 1 uncontaminated to moderately contaminated; 1 < Igeo < 2 moderately contaminated; 2 < Igeo < 3 moderately to heavily contaminated; 3 < Igeo < 4 heavily contaminated; 4 < Igeo < 5 heavily to very heavily contaminated; Igeo > 5 very heavily contaminated.

Enrichment Factor
The enrichment factor (EF) was used to assess the degree of the enrichment, depletion, or possible specific sources of toxic trace elements in dry atmospheric deposition [15], as defined in Equation (2): where C i and C ref are the measured concentrations of toxic trace element i and reference element for normalization, respectively; B i and B ref are the contents of the toxic trace element i and reference element, respectively.In this study, Al was chosen as the reference element.Here, the surface soil data from the mid-west Jilin Province [41] rather than the upper continental crust data were used for EF calculations in order to reduce the primary uncertainty resulting from the differences between chemical composition of local soil and reference crustal composition.Five enrichment categories were recognized on the basis of EF [43,44]: EF < 2, deficiency to minimal enrichment; EF = 2-5, moderate enrichment; EF = 5-20, significant enrichment; EF = 20-40, very high enrichment; EF > 40, extremely high enrichment.

Potential Ecological Risk Index
The potential ecological risk index (RI) proposed by Hakanson [45] was used to assess the potential ecological risk posed by toxic trace elements in dry atmospheric deposition, which was calculated using Equation (3): where E i r is the potential ecological risk factor (Er) of toxic trace element i, RI is the potential ecological risk index of multiple toxic trace elements, C i f is the pollution coefficient of toxic trace element i, Ti is the toxic response factor of toxic trace element i.The values for toxic trace elements were Cu = Ni = Pb = Co = 5, Mn = Zn = 1, As = 10, Cd = 30, Hg = 40, and Mo = Se = 15, respectively [46-49].C i s is the concentration of toxic trace element i, C i n is the background level of toxic trace element i.The averaged concentrations of toxic trace elements in the surface soil data from mid-west Jilin Province [43] were used as the background level of toxic trace elements.The following criteria values were used to assess the ecological risk levels [45]: E i r < 40, low risk; 40 ≤ E i r < 80, moderate risk; 80 ≤ E i r < 160, considerable risk; 160 ≤ E i r < 320, high risk; E i r ≥ 320, serious risk.RI <150, low risk; 150 ≤ RI < 300, moderate risk; 300 ≤ RI < 600, considerable risk; RI ≥ 600, very high risk.

Health Risk Assessment
The health risks posed by toxic trace elements in dry atmospheric deposition were evaluated with the method recommended by the US Environmental Protection Agency [33][34][35].
For non-carcinogenic risk assessment of As, Cd, Co, Cr, Cu, Hg, Mn, Mo, Ni, Pb, Se, and Zn in dry atmospheric deposition, their average daily dose (ADD), hazard quotient (HQ), and the hazard index (HI) through ingestion (hand-mouth route), inhalation (mouthnose route), and dermal absorption (skin contact) were estimated using Equations ( 4)-( 8), respectively: where ADD (mg/kg•d) represents the dose through ingestion (ADD ing ), inhalation (ADD inh ), and dermal absorption (ADD derm ).Hazard quotient (HQ) is calculated as ADD divided by reference dose (RfD).The hazard index (HI) is the sum of HQ.If HQ or HI values of toxic trace elements is below one, the non-carcinogenic risk or the adverse health effects may be low or negligible.If HQ or HI is over one, non-carcinogenic risk or adverse health effects may be indicated.
For carcinogenic risk assessment of As, Cd, and Cr in dry atmospheric deposition, their carcinogenic risks (CR) and total carcinogenic risk (TCR) were calculated as Equations ( 9)- (10), respectively: where carcinogenic risk (CR i ) for toxic trace element i via three exposure pathways is calculated by the chronic daily dose (ADD) with the corresponding cancer slope factor (SF), respectively.TCR is the sum of CR i .The carcinogenic risks can be divided into three levels: CR < 10 −6 , negligible; 10 −6 < CR < 10 −4 , tolerable or acceptable; CR > 10 −4 , considerable or high.
The parameters and reference values used in the Equations ( 4)-( 9) are listed in Tables 1 and 2.

Multivariate Statistical Analysis
PCA with varimax rotation and Kaiser normalization was carried out to minimize the number of variables with high loading on each component, interpreted in accordance with the possible sources of toxic trace elements in dry atmospheric deposition.
CA was performed with the standardized data by Z score using Ward's method with squared Euclidean distances to group toxic trace elements with similar origin or sources, and to check and supplement PCA results.
Geostatistical analysis was carried out employing MapGis 6.7 software with Kriging interpolation to characterize the spatial distribution of elements in atmospheric deposition.
Geostatistical analysis was carried out employing MapGis 6.7 software with Kriging interpolation to characterize the spatial distribution of toxic trace elements in dry atmospheric deposition, which was helpful for finding hot spots or assessing pollution sources of toxic trace elements in dry atmospheric deposition.

Dry Atmospheric Deposition
The concentrations of toxic trace elements in dry and wet fractions of atmospheric depositions from the study area are listed in Table 3.The toxic trace elements in dry atmospheric deposition followed a decreasing abundance order of Mn > Zn > Cr > Pb > Cu > Ni > As > Co > Mo > Cd > Se > Hg.Mn, Zn, and Cr had mean concentrations of 415.36 mg/kg, 105.69 mg/kg, and 46.53 mg/kg, respectively (Table 3), while Hg had the lowest mean concentration of 0.061 mg/kg.The mean concentrations of Cd, Zn, Hg, Cu, Mo, Se, and Pb in dry atmospheric deposition were higher than their corresponding values in soils from the mid-west Jilin Province [43].The variation coefficients (VC, mean/standard derivation) of toxic trace elements in dry atmospheric deposition were less than one, with a decreasing order of Cd > Cu > Se > Hg > Zn > Mo > Ni > Pb > As > Co > Cr = Mn.Lower CV values for Co, Cr, and Mn indicate that their concentrations were less variable, implying that they were mainly from natural sources, while higher CV values for Cd, Cu, Hg, Mo, Se, and Zn suggested a wide variation in their concentrations, resulting from the effects of anthropogenic activities.Al had the lowest mean VC value of 0.12, reflecting that Al is a less variable and conservative element in dry atmospheric deposition, so Al was chosen as the reference element to calculate the enrichment factor in this study.The mean concentrations of Mn and Zn were 495.02 µg/L, 318.75 µg/L, respectively (Table 3), while Hg had a minimum mean concentration of 0.041 µg/L.The concentrations of toxic trace elements in wet atmospheric deposition decreased in the order of Mn > Zn > As > Pb > Cr > Ni > Cu > Co > Mo > Se > Cd > Hg.In wet atmospheric deposition, Cd and Pb had VC values of more than one, while other toxic trace elements had values less than one.The concentrations of toxic trace elements in wet atmospheric deposition actually included the dissolved materials in precipitation and dry atmospheric deposition.

Bulk Atmospheric Deposition Fluxes
Bulk atmospheric deposition fluxes of toxic trace elements in the study area were obtained as the sum of both dry and wet deposition fluxes, which were calculated using Equation (11): where Q is the annual atmospheric deposition fluxes (mg/m 2 /year) of toxic trace elements of interest, C d and W d represent the concentration (mg/kg) and the weight (g) of toxic trace elements of interest in the dry fraction of an atmospheric deposition sample, C w and V w denote the concentration (mg/L) and the volume of toxic trace elements of interest in the wet fraction of an atmospheric deposition sample, and S is the available section area of sampling container (314 cm 2 ).The annual mean atmospheric deposition fluxes of toxic trace elements are shown in Table 4.The yearly atmospheric depositional bulk fluxes of t toxic trace elements in the study area decreased in the order of Mn > Zn > Cr > Pb > Cu > Ni > As > Co > Mo > Cd > Se > Hg.Mn had bulk atmospheric depositions of 237.41 mg/m 2 /year, while Hg had the lowest annual mean atmospheric depositions of 0.02 mg/m 2 /year.The annual deposition fluxes of toxic trace elements in dry atmospheric depositions were much greater than those in wet atmospheric depositions, implying that dry atmospheric deposition predominated in the study area.Therefore, the main concern regarding toxic trace elements in dry atmospheric deposition in the study area should be evaluating their source, and ecological and health risks in atmospheric deposition.
The fluxes for As, Co, Cr, Mn, and Ni in dry atmospheric deposition exceeded those reported atmospheric deposition for other regions in China (Table 4), while the depositional fluxes for Cd, Cu, Hg, Pb, and Zn were comparatively lower than those reported for other regions in China [12][13][14]17,19,53], implying that the atmospheric emissions of As, Co, Cr, Mn, and Ni in the study area are higher than those in other areas.

Spatial Variations of Bulk Atmospheric Deposition Fluxes
The annual bulk accessions of toxic trace elements through atmospheric deposition showed a clear spatial heterogeneity in different counties in the study area (Figure 2).Comparatively, the yearly depositional fluxes for most toxic trace elements in Tongyu and Taonan counties were substantially higher than those in other counties.The calculated bulk depositional fluxes for As, Co, Cr, Mn, and Ni in the study area exceeded those reported for other regions in China (Table 4), while the bulk depositional fluxes for Cd, Cu, Hg, Pb, and Zn were comparatively lower than those reported for other regions in China [12,13,17], implying that the atmospheric emissions of toxic trace elements such as As, Co, Cr, Mn, and Ni in the study area are higher than those in other regions.
Taonan counties were substantially higher than those in other counties.The calculated bulk depositional fluxes for As, Co, Cr, Mn, and Ni in the study area exceeded those reported for other regions in China (Table 4), while the bulk depositional fluxes for Cd, Cu, Hg, Pb, and Zn were comparatively lower than those reported for other regions in China [12,13,17], implying that the atmospheric emissions of toxic trace elements such as As, Co, Cr, Mn, and Ni in the study area are higher than those in other regions.1% exceeded one, implying that these elements did not contaminate most samples.Igeo values of Cd of more than one were present in 70% of samples; 14%, 45%, and 11% of samples had Igeo values of Cd in the range of 1-2, 2-3, and 3-4, respectively, indicating that most samples were moderately to heavily contaminated with Cd.

Ecological Risk Evaluation
The potential ecological risk coefficient (E ) and ecological risk index (RI) of toxic trace elements in dry atmospheric deposition are shown in Figure 3.The E values of Mn, Cu, Co, Ni, Zn, As, Pb, and Cr in all samples were less than 40, indicating these toxic trace elements had low ecological risk.The E values of Se, Mo, Hg, and Cd were in the ranges

Ecological Risk Evaluation
The potential ecological risk coefficient (E i r ) and ecological risk index (RI) of toxic trace elements in dry atmospheric deposition are shown in Figure 3.The

Non-Carcinogenic Risk Evaluation
The results of health risk for non-carcinogenic and carcinogenic effects on children and adults are presented in Table 3.The average values of HQ via the ingestion pathway were generally higher than those through pathways of inhalation and dermal contact for adults and children (Figure 4).Moreover, HQ values of toxic trace elements calculated for children were generally greater than those for adults, meaning that adults have more tolerance than children.The HI values of As, Cd, Co, Cr, Cu, Hg, Mn, Mo, Pb, Se, Ni, and Zn for adults were in the range of 1.04 × 10  5, Figure 4) values of toxic trace elements for all three pathways were less than one, indicating that no non-carcinogenic effects to adults and children were posed by toxic trace elements in dry atmospheric deposition the study area.

Carcinogenic Risk Evaluation
The averaged values of CR for As, Cd, and Cr via the ingestion pathway were generally higher than those through pathways of inhalation and dermal contact for both adults and children (Figure 5), implying that carcinogenic risk of As, Cd, and Cr through ingestion are higher than those through inhalation and dermal contact.The carcinogenic risks of As and Cr in dry atmospheric deposition through three exposure pathways decreased in the order of ingestion > inhalation > dermal contact for both adults and children, while the carcinogenic risks of Cd in dry atmospheric deposition through three exposure pathways decreased in the order of ingestion > dermal contact > inhalation for both adults and children.The mean TCR of As, Cd, and Cr for adults was 4.18 × 10 6 , 6.34 × 10 −8 , and 5.89 × 10 −5 , with a range of 9.60 × 10

Carcinogenic Risk Evaluation
The averaged values of CR for As, Cd, and Cr via the ingestion pathway were generally higher than those through pathways of inhalation and dermal contact for both adults and children (Figure 5), implying that carcinogenic risk of As, Cd, and Cr through ingestion are higher than those through inhalation and dermal contact.The carcinogenic risks of As and Cr in dry atmospheric deposition through three exposure pathways decreased in the order of ingestion > inhalation > dermal contact for both adults and children, while the carcinogenic risks of Cd in dry atmospheric deposition through three exposure pathways decreased in the order of ingestion > dermal contact > inhalation for both adults and children.The mean TCR of As, Cd, and Cr for adults was 4.18 × 10 6 , 6.34 × 10 −8 , and 5.89 × 10 , with a range of 9.60 × 10

The Uncertainty of Human Risk Assessment
It should be noted that the health risk assessment results have some uncertainties in this study.First, some exposure parameters used in the study come from the citizens of the USA, and may not be suitable for Chinese citizens.Second, some parameter values may be overestimated using the UAEPA risk assessment model.For example, the mean weight of dry atmospheric deposition was 21.87 mg/day with a range of 1.83-176.80mg/day in this study area, and only 5.4% samples had weights of dry atmospheric deposition of more than 100 mg/day; while the intake via ingestion (IngR) values used in the study were 100 mg/day, and 200 mg/day for adults and children, respectively.Therefore, the overestimated daily dose of toxic trace elements through ingestion may have enhanced their corresponding non-carcinogenic risks.Third, the particle sizes of atmospheric deposition and available forms of toxic elements in atmospheric deposition may also have contributed some uncertainties in assessment of their human health risks.In order to deduce the uncertainty of risk assessment and obtain reliable results, the exposure parameter values suitable for China should be adopted and the particle sizes of atmospheric deposition and available forms of toxic elements in atmospheric deposition should also be considered.

The Uncertainty of Human Risk Assessment
It should be noted that the health risk assessment results have some uncertainties in this study.First, some exposure parameters used in the study come from the citizens of the USA, and may not be suitable for Chinese citizens.Second, some parameter values may be overestimated using the UAEPA risk assessment model.For example, the mean weight of dry atmospheric deposition was 21.87 mg/day with a range of 1.83-176.80mg/day in this study area, and only 5.4% samples had weights of dry atmospheric deposition of more than 100 mg/day; while the intake via ingestion (IngR) values used in the study were 100 mg/day, and 200 mg/day for adults and children, respectively.Therefore, the overestimated daily dose of toxic trace elements through ingestion may have enhanced their corresponding non-carcinogenic risks.Third, the particle sizes of atmospheric deposition and available forms of toxic elements in atmospheric deposition may also have contributed some uncertainties in assessment of their human health risks.In order to deduce the uncertainty of risk assessment and obtain reliable results, the exposure parameter values suitable for China should be adopted and the particle sizes of atmospheric deposition and available forms of toxic elements in atmospheric deposition should also be considered.PCA is commonly used to distinguish the possible sources of toxic trace elements in dry atmospheric deposition.The concentration data of toxic trace elements in dry atmospheric deposition were suitable for PCA analysis according to the KMO value (0.907) and Bartlett spherical test (0.000).Two principal components were extracted with eigenvalues more than one, accounting for 76.223% of the total variance.PC1 accounted for 68.864% of the total variance, which was dominated by Cu, Zn, Hg, Pb, Se, Cd, Mo, As, and Ni, with positive loadings of more than 0.5, indicating their strong correlations with a common source.The anthropogenic sources of As, Cd, Hg, and Se in atmospheric emission in the study area are possibly coal combustion and biomass burning [10,19,21,[53][54][55][56][57][58].The anthropogenic sources of Cu, Pb, and Zn are mainly coal combustion, traffic emission, and waste incineration [21,59,60].The anthropogenic Mo more likely originated from coal combustion, waste incineration, and biomass burning [60].Therefore, toxic trace elements in the first component are considered as the impacts of anthropogenic activities.PC2 with a total variance of 7.359% was heavily loaded by Al, Co, Fe, Cr, Mn, Ni, Pb, and As, with positive loadings over 0.5, suggesting that they have a good correlation and a similar source.Al and Fe in this component are naturally distributed as major elements in soil minerals.Al is mainly associated with aluminosilicate minerals such as plagioclase (14.91%), alkali feldspar (13.14%), and clay minerals such as illite (3.41%), mixed-layer illite-smectite (2.82%) and kaolinite (1.59%), while Fe and Mn mainly occur as oxide minerals, and Co, Cr, and Ni are commonly adsorbed onto clay and oxide mineral surfaces or substituted for Fe and Mn into amphiboles or other silicate minerals [61].The toxic trace elements in this component are those typical of soil-derived dust particles (crustal source).Ni, As, and Pb occurred in two components with loadings of more than 0.5, implying that they may have been from natural and anthropogenic sources.
In Figure 6, the toxic trace elements, together with Al an Fe, in dry atmospheric deposition can be divided into two groups: group one includes Co, Cr, Mn, Ni, Al, and Fe, and another group is composed of As, Mo, Se, Pb, Hg, Zn, Cu, and Cd.[61].The toxic trace elements in this component are those typical of soil-derived dust particles (crustal source).Ni, As, and Pb occurred in two components with loadings of more than 0.5, implying that they may have been from natural and anthropogenic sources.
In Figure 6, the toxic trace elements, together with Al an Fe, in dry atmospheric deposition can be divided into two groups: group one includes Co, Cr, Mn, Ni, Al, and Fe, and another group is composed of As, Mo, Se, Pb, Hg, Zn, Cu, and Cd.

Cluster Analysis (CA)
The CA results for toxic trace element contents in dry atmospheric deposition are illustrated in Figure 7. Two distinct clusters are distinguished from the dendrogram of CA.The first cluster includes Fe, Co, Mn, Al, Ni, and Cr.The second cluster consists of Zn, Pb, Cu, Se, Cd, As, Hg, and Mo.The results of CA fit those of PCA well.Moreover, toxic trace elements in cluster 1 and cluster 2 join together at a relatively higher correlation level, implying toxic trace elements in both clusters are correlated.The high correlation coefficients may imply that the dry atmospheric deposition is partly contaminated to some extent by Mo, As, Hg, Zn, Pb, Cu, Se, and Cd.The CA result suggests that toxic trace elements in cluster 2 may have originated from both crustal and anthropogenic sources.

Spatial Distribution of Toxic Trace Elements
The spatial distribution of toxic trace elements can be used as an aid to identify hot spots and their possible sources.The spatial distributions of toxic trace elements in dry atmospheric deposition in the study area are shown in Figure 8. As, Cd, Cu, Hg, Mo, Pb, Se, and Zn show a similar spatial distribution pattern, with a relatively higher concentration or hot spots in the densely populated urban area of the Baicheng city (sampling location 37) and its surrounding areas (sampling locations 31, 32, and 43), where the mean concentrations of As, Cd, Cu, Hg, Mo, Pb, Se, and Zn are over 12.1, 1,688, 67.88, 0.095, 1.459, 50.47, and 0.509,170 mg/kg, respectively, implying that their sources may be anthropogenic activities, such as coal combustion, traffic, and industrial emissions.Some hot spots of As, Cd, Hg, Mo, Pb, Se, and Zn appear in the sparsely populated rural areas, suggesting their crustal sources.Co, Cr, Mn, and Ni have quite different distribution patterns from other elements, without the hot spots in the urban area of the Baicheng city, indicating that these four elements may have primarily originated from crustal sources.

Spatial Distribution of Toxic Trace Elements
The spatial distribution of toxic trace elements can be used as an aid to identify hot spots and their possible sources.The spatial distributions of toxic trace elements in dry atmospheric deposition in the study area are shown in Figure 8. As, Cd, Cu, Hg, Mo, Pb, Se, and Zn show a similar spatial distribution pattern, with a relatively higher concentration or hot spots in the densely populated urban area of the Baicheng city (sampling location 37) and its surrounding areas (sampling locations 31, 32, and 43), where the mean concentrations of As, Cd, Cu, Hg, Mo, Pb, Se, and Zn are over 12.1, 1.688, 67.88, 0.095, 1.459, 50.47, and 0.509, 170 mg/kg, respectively, implying that their sources may be anthropogenic activities, such as coal combustion, traffic, and industrial emissions.Some hot spots of As, Cd, Hg, Mo, Pb, Se, and Zn appear in the sparsely populated rural areas, suggesting their crustal sources.Co, Cr, Mn, and Ni have quite different distribution patterns from other elements, without the hot spots in the urban area of the Baicheng city, indicating that these four elements may have primarily originated from crustal sources.
67.88, 0.095, 1.459, 50.47, and 0.509,170 mg/kg, respectively, implying that their sources may be anthropogenic activities, such as coal combustion, traffic, and industrial emissions.Some hot spots of As, Cd, Hg, Mo, Pb, Se, and Zn appear in the sparsely populated rural areas, suggesting their crustal sources.Co, Cr, Mn, and Ni have quite different distribution patterns from other elements, without the hot spots in the urban area of the Baicheng city, indicating that these four elements may have primarily originated from crustal sources.

Conclusions
This study was carried out to measure the concentrations, to assess the risks, and to identify the sources of toxic trace elements in atmospheric deposition in the Baicheng-Songyuan area, Jilin province, China.The bulk fluxes of toxic trace elements decreased in the order of Mn > Zn > Cr > Pb > Cu > Ni > As > Co > Mo > Cd > Se > Hg.The dry deposition fluxes of toxic trace elements accounted for 74.7-99.9% of their bulk deposition fluxes, indicating that dry atmospheric deposition predominates in the atmospheric deposition in the study area.Cd showed moderate to heavy contamination in dry atmospheric deposition, while other toxic trace elements showed no contamination; Cd and Hg in dry atmospheric deposition showed moderate to considerable ecological risk, while other toxic trace elements were at low ecological risk levels.The non-carcinogenic risks of toxic trace elements in dry atmospheric deposition through three exposure pathways decreased in the order of ingestion > inhalation > dermal contact for adults and children.HQ and HI of toxic trace elements via three pathways in dry atmospheric deposition were less than one, suggesting a low non-carcinogenic risk.The averaged values of CR for As, Cd, and Cr via the ingestion pathway were generally higher than those through pathways of inhalation and dermal contact for both adults and children.The mean TCR values of As and Cr for adult and children were between 10 −6 and 10 −4 , sugesting that the carcinogenic risk levels of As and Cr were tolerable or acceptable, and the mean TCR value of Cd for adult and children was less than 10 −6 , implying that the carcinogenic risk level of Cd was negligible.Mn, Ni, Cr, and Co in dry atmospheric deposition were mainly originated from the crustal sources, while As, Cd, Cu, Hg, Mo, Pb, Se, and Zn in dry atmospheric deposition were derived from both crustal and anthropogenic sources.
It should be emphasized that health risk assessment of toxic elements in atmospheric deposition may be uncertain due to the exposure parameters being unsuitable to Chinese, and consideration only of the total concentration of toxic elements rather than their speciation in atmospheric deposition, and the particle sizes of atmospheric deposition.
The results obtained in this study advocate the necessity for monitoring atmospheric deposition in the vast rural area, and also provide a scientific basis for controlling contamination posed by toxic trace elements in atmospheric deposition.The approach by integration of multiple analysis and spatial distribution to identify the source of toxic trace elements in dry atmospheric deposition is available for other researchers.

Figure 1 .
Figure 1.Sample locations of atmospheric deposition in the study area.

Figure 1 .
Figure 1.Sample locations of atmospheric deposition in the study area.

Figure 2 .
Figure 2. Comparison of yearly depositional fluxes of toxic trace elements in different counties of the study area.

Figure 2 .
Figure 2. Comparison of yearly depositional fluxes of toxic trace elements in different counties of the study area.

Figure 3 .
Figure 3. Igeo, EF, and potential risk index of toxic trace elements in dry fraction of atmospheric deposition.

Figure 3 .
Figure 3. Igeo, EF, and potential risk index of toxic trace elements in dry fraction of atmospheric deposition.
E i r values of Mn, Cu, Co, Ni, Zn, As, Pb, and Cr in all samples were less than 40, indicating these toxic trace elements had low ecological risk.The E i r values of Se, Mo, Hg, and Cd were in the ranges 1.67-61.58,7.38-48.93,5.58-228.54,and 13.77-637.41,with their average values being 21.44, 22.27, 71.76, and 198.25, respectively.The E i r values of Se in 10.7% samples, and Mo in 5.4% samples exceeded 40, implying that both elements in most samples were of low ecological risk.The E i r values of Cd in 87.5% samples, and Hg in 76.8% were over 40, indicating that Cd and Hg had the highest potential ecological risk in the study area.The samples with E i r values of Cd in the ranges 40-80, 80-160, 160-320, and over 320 accounted for 14.3%, 14.3%, 42.9%, and 16.1%, respectively, indicating Cd was a moderate to serious risk.The samples with E i r values of Hg in the ranges of 40-80, 80-160, and 160-320 accounted for 39.3%, 32.1%, and 19.4%, respectively, indicating Hg was a moderate to high risk.RI values ranged from 45.66 to 1010.1, with a mean value of 352.7.The samples with RI values of 150-300, 300-600, and over 600 accounted for 23.2%, 55.4%, and 3.6%, respectively, reflecting that toxic trace elements in dry atmospheric deposition in most samples posed a moderate to considerable potential ecological risk.

Figure 4 .
Figure 4. Comparison of HQ values of toxic trace elements via three pathways for adults and children.

Figure 4 .
Figure 4. Comparison of HQ values of toxic trace elements via three pathways for adults and children.

Figure 5 .
Figure 5. CR and TCR values of As, Cd, and Cr via three pathways for adults and children.

Figure 5 .
Figure 5. CR and TCR values of As, Cd, and Cr via three pathways for adults and children.

Figure 6 .
Figure 6.Plot PC1 against PC2 for toxic trace elements in dry atmospheric deposition.

Figure 6 .
Figure 6.Plot PC1 against PC2 for toxic trace elements in dry atmospheric deposition.

Atmosphere 2024 , 20 Figure 7 .
Figure 7. Dendogram from Ward's method of hierarchical cluster analysis of toxic trace elements in dry atmospheric deposition.

Figure 7 .
Figure 7. Dendogram from Ward's method of hierarchical cluster analysis of toxic trace elements in dry atmospheric deposition.

Figure 8 .
Figure 8.The spatial distribution of toxic trace elements in dry atmospheric deposition.

Table 3 .
Descriptive statistics of toxic trace elements in dry and wet atmospheric depositions.

Table 4 .
Flux comparisons of toxic trace elements from atmospheric deposition in the study area and others (mg/m 2 /year).
1.67-61.58,7.38-48.93,5.58-228.54,and 13.77-637.41,with their average values being 21.44, 22.27, 71.76, and 198.25, respectively.The E values of Se in 10.7% samples, and Mo in 5.4% samples exceeded 40, implying that both elements in most samples were of low ecological risk.The E values of Cd in 87.5% samples, and Hg in 76.8% were over 40, indicating that Cd and Hg had the highest potential ecological risk in the study area.The samples with E values of Cd in the ranges 40-80, 80-160, 160-320, and over 320

Table 5 .
Non-carcinogenic and carcinogenic risks for adults and children.