Environmental Risk Assessment for PM 2.5 Pollution from Non-Point Sources in the Mining Area Based on Multi-Source Superposition and Diffusion

: To identify high-concentration contributing sources and highly dispersive pollution sources of ﬁne particulate matter, analyze the relationship between the location and distribution shape of emission sources and the concentration contribution and dispersion of particulate matter, and realize the atmospheric environment risk simulation and the differential control of non-point sources in the mining area, taking a large mining area in Inner Mongolia as an example, we classiﬁed the emission sources of PM 2.5 (particulate matter less than 2.5 µ m) and complied with the emission inventory. A CALPUFF model was used to simulate the contribution for the PM 2.5 concentration of six types of emission sources and a multi-source superposition. Through scenario simulation, we analyzed the relationship between the spatial distribution of emission sources and the emission concentration and diffusion in a large mining area. We analyzed the relative risks of six types of sources under the inﬂuence of other superimposed sources and the change of emission concentration during transmission. We established a three-dimensional evaluation model to assess the atmospheric environmental risk of PM 2.5 non-point sources in the mining area, considering the change rate of PM 2.5 concentration with migration, the relative contribution ratio of superimposed sources, and the equal contribution index of the standard concentration. The results show that the maximum equal contribution index of standard concentration of multi-source superposition was 4.40. Among them, non-paved roads, exposed surface sources of coal piles, and exposed surface sources of mine pits and dumps were the top three pollution contributors, and their maximum equal contribution indexes of standard concentration were 2.40, 2.21, and 2.10, respectively. The effect of superimposed pollution sources was affected by the wind ﬁeld and the spatial distribution density of emission sources, while the dispersion was affected by the wind direction and the distribution direction of pollution sources. In the case of the same source intensity and emission area, the denser the source distribution was, the higher the emission concentration was, the smaller the contribution ratio of superimposed sources was, and the greater the relative pollution risk was. When the angle between the direction of dispersed linear pollution research concludes that the dispersed distribution of pollution sources can reduce pollution risk, and the smaller the angle is between the linear distribution direction of pollution sources and the dominant wind direction, the smaller the diffusion risk is. Therefore, differentiated control can be carried out according to the characteristics of pollution sources. The conclusions can provide methods and theoretical support for the control of atmospheric environment risk, pollution prevention, and control planning in mining areas.


Introduction
There are many large-scale mining areas in northern China, and the emissions of various sources of particulate matter in open-pit mining areas has a complex impact on the formation of haze and dust in Beijing, Tianjin, Hebei, and other central regions [1,2]. In particular, non-point sources are difficult to control, with low control efficiency, high pollution risk, and great health hazard [3][4][5]. At present, a "one-size-fits-all" control approach is mostly adopted. The existing particulate matter risk assessment mostly takes the particulate matter concentration as the risk source and the population as the risk receptor to evaluate the exposure health risks [6][7][8][9][10]. In traditional risk assessments, all types of non-point sources in the study area have been used as a whole for simulation, which fails to consider the differences in the migration and diffusion rules of particles emitted by different non-point sources [11][12][13], or only the pollution sources within the study area are taken as the object, ignoring the overlapping risks that are caused by scattered and non-point sources around the research object. Therefore, calculating the relative contribution and dispersion of PM 2.5 (particulate matter less than 2.5 µm) among multiple types and quantities of emission sources, and analyzing the relationship between the spatial distribution of pollution sources and the emission concentration and dispersion of PM 2.5 will improve the accuracy of atmospheric environmental risk assessments, and provide a basis for differentiated management of pollution sources.
In addition, there has been no unified standard for emission inventory of particulate matter in mining areas. Currently, the AP-42 formula recommended by the US Environmental Protection Agency (EPA) is widely used [14,15], or the empirical formula obtained through field wind tunnel experiments is [16][17][18][19], the accuracies of which are varying [20][21][22]. In 2014, the Ministry of Environmental Protection of the People's Republic of China compiled the "Technical Guideline for the Preparation of Emission Inventory of Particulate Matter from Dust Sources (Trial) (Technical Guideline for short)", which is more suitable for China's reality [23][24][25]. Therefore, taking a large mining area in Inner Mongolia as the research object, and regarding the Technical Guideline, we classified the emission sources of PM 2.5 into paved roads, non-paved roads, loading and unloading, coal piles, quarries, and wind erosion on the exposed surface of mine pits and dumps, and then compiled PM 2.5 emission inventory. Based on the CALPUFF model [26][27][28][29], we performed a superposition of six sources in the simulation (multi-source superposition) and six types of single-source simulations, respectively. By constructing three indices, namely, the equal contribution index of standard concentration, the relative contribution ratio of superimposed sources, and the change rate of PM 2.5 concentration with migration, we calculated and evaluated the relative pollution risk and diffusion risk of the PM 2.5 of various pollution sources, as well as the overall environmental risk of multi-source superposition to identify the high-risk sources and get the spatial distribution of PM 2.5 environmental risk. Based on this, we divided the environmental risks into zones and proposed control schemes to avoid a one-size-fits-all control and achieve differentiated management.

The Research Area
The research area is located in Etoke Banner, Ordos City, Inner Mongolia, and has the characteristics of a large mining area in the west. The research scope is shown in Figure 1. The mining area is located in arid and semi-arid areas, which have a typical temperate continental climate, with the characteristics of sharp changes of cold and heat on the plateau, long sunshine hours, and intense solar radiation. Precipitation is low and concentrated, evaporation is extensive, the wind is strong, sand is abundant, natural vegetation is sparse, and ecological regulation ability is weak. Specific meteorological conditions are shown in Table S1. Meteorology contributes to the dust formation and diffusion of particulate matter and is not conducive to environmental governance. Moreover, frequent development activities can easily cause severe particulate pollution. Thus, it is necessary to control the emission sources of particulate matter in mining areas. The locations of the ground weather station (WS) and three monitoring points (MS1, MS2, and MS3) are shown in Figure 1. The three monitoring points observed the PM 2.5 concentration in the mining area.

Compilation of Particulate Matter Emission Inventory
For this paper, we selected the emission factor method, modified the calculation formulas and parameters according to the Technical Guideline, and combined it with the survey data to obtain the emission factors and activity levels of particulate matter from each pollution source in the mining area, and realize the localization of emission inventory. The calculation is shown in Equation (1).
In the formula, W i is the dust emission of emission source i; E i is the emission coefficient of the unit activity level corresponding to emission source i; A i is the activity level factor of dust source i; η is the removal efficiency of the pollution control technology on the dust; n is the number of emission source types; W t is the total particulate emissions from dust sources.
Based on various literature studies [6,14,16,30,31] and field investigations, there are many dust emission sources in mining areas. Their characteristics are generalized into 3 categories: road transport dust sources, exposed surface dust sources, and material loading and unloading dust sources (see Table 1). The short-term particulate matter emission sources, such as drilling and blasting, topsoil removal, and material transfer, were uniformly classified as the material loading and unloading dust sources, and the particulate matter emissions were converted into the corresponding material loading and unloading volumes. The parameter values [23,24,32] are shown in Tables S2 and S3. We allocated the inventory to a specific space according to land use types and field operations based on remote sensing interpretation and positioning [33][34][35].

. Model Setup and Simulation
The CALPUFF is an unsteady Lagrangian puff model that simulates pollutant migration, transformation, and removal in three-dimensional spatially and temporally variable meteorological conditions. The basic concentration equation of a single puff at a certain acceptance point is shown in Equation (2).
where C is ground concentration, g/m 2 ; Q is the source intensity; σ x , σ y , and σ z are diffusion coefficients; d a is the downwind distance; d c is the transverse distance; H e is the effective height; h is the height of the mixing layer; g is the vertical term of the Gaussian equation, which solves the problem of multiple reflections between the mixing layer and the ground. The model consists of a weather processor (CALMET), a dispersion model (CALPUFF), and a post-processor (CALPOST). The CALMET model is a weather model that generates a three-dimensional hourly wind field within a three-dimensional grid modeling domain. The CALPUFF model was used to simulate the PM 2.5 concentration emitted from pollution sources in 2018, from 1 January to 31 December. The data sources and parameter settings are shown in Table S4. For the CALMET simulations, the vertical discretization consisted of 11 layers that were 0 m, 20 m, 40 m, 80 m, 160 m, 320 m, 640 m, 1200 m, 2000 m, 3000 m, and 4000 m above the land surface. The model parameters BIAS, RMAX1, RMAX2, TERRAD, R1, and R2 were selected according to the specific site. The BIAS parameter assigns weights to the ground and upper air station data for each vertical layer. The surface data gives 100% weight (BIAS = 1) in the first layer and 0 weight (BIAS = 1) in the last 2 vertical layers. Equal weight was assigned to the fourth layer (BIAS = 0) with a gradation of weights between the lower and upper for the remaining layers. The TERRAD parameter defines the influence radius of the terrain features and is set to 10 km. The parameters RMAX1 and RMAX2 define the maximum influence radii for the surface and upper data over the land surfaces, and were set to 50 km and 100 km, respectively. The parameters R1 and R2 are the distances from the observation, where the weight of the observation and the initial prediction field to the surface and upper layer is equal. R1 and R2 were set to 18 km and 54 km, respectively. For the remaining parameters, we used the CALMET default values.
MS1, MS2, and MS3 monitor the average daily concentration of PM 2.5 . The monitoring periods are shown in Table 2. The model was simulated 5 times and was verified by the daily mean concentration data of PM 2.5 monitored by MS1, MS2, and MS3 (for a total of 12 days). Then, we screened out the meteorological conditions for the maximum probability of pollution in the mining area and used the CALPUFF to simulate the PM 2.5 concentration emitted by 6 types of dust sources and multi-source superposition.

Verification of Simulation
The accuracy of the model was evaluated by 6 indices [20,36,37], namely FB (fractional bias), MG (geometric mean bias), VG (geometric variance), R 2 (coefficient of determination), NMSE (normalized mean square error), and FAC2 (fraction of two). The calculation formulas and evaluation criteria are shown in Table 3. FB is a measure of the average deviation. A FB of 0.6 is equal to 2 times that of the under-prediction of the model, and a negative value indicates over-prediction of the model. The MG reflects the deviation degree of the geometric mean, and the values of 0.5 and 2.0 represent a factor of 2 over-prediction and under-prediction, respectively. The NMSE is a measure of the opposite difference, and a value of 1.0 indicates that the typical difference between the prediction and observation is approximately equal to the mean. The NMSE and FB are appropriate when the typical difference between the prediction and observation is about a factor of 2 and the range of prediction and observation in the data set is small (a factor less than 2). A VG value of 1.6 indicates a difference of 2 values between the predicted data and the observed data. The FB and MG are measures of mean relative deviation and represent only systematic errors, while the NMSE and VG are measures of mean relative scattering and reflect systematic and non-systematic (random) errors. The FAC2 is the most reliable measurement, reflecting the percentage of the predicted concentration within the range of the two observation factors.

Analysis of Dust Particle Distribution and Spatial Transport Laws from Pollution Sources in Multi-Scale Mining Areas
Based on the simulation, we first analyzed the spatial distribution law of PM 2.5 from each single source, and examined the spatial distribution law of PM 2.5 from multi-source superposition. Then, we analyzed the risk value changes of each single source before and after superimposed by multiple sources. We calculated the PM 2.5 concentration change rate from each single source with the increase of transport distance. There were 3 indices, namely the equal contribution index of standard concentration I ij , the relative contribution ratio of superimposed sources r ij , and the change rate of PM 2.5 migration concentration k ij ', that were used to characterize the contribution of PM 2.5 concentration, the increasing/reducing effect of superimposed sources on air pollution risk, and the spatial variation of PM 2.5 concentration. The calculation formulas are shown in Table 4. Table 4. The parameters and formulas of the particulate matter distribution and spatial transport law in the mining area.

Number
Formula/Parameter Significance Equal contribution index of standard concentration of pollutant source i at the location j C ij The PM 2.5 concentration contribution of single pollutant source i at the location j Cs The standard concentration limit for PM 2.5 i The serial number of the pollution source n Total number of pollution sources j A place at some distance from the pollutant source The relative contribution ratio of superimposed sources for pollutant source i at the location j The PM 2.5 concentration contribution of the multi-source superposition at the location j The PM 2.5 concentration at a location that is h x meters away from the center of the source i C i0 The PM 2.5 concentration at the center of the source i The PM 2.5 concentration change at the location that is h meters away from the center of the source i C ij The PM 2.5 concentration at the location that is h meters away from the center of the source i The PM 2.5 concentration at the location that is h-∆h meters away from the center of the source i The change rate of PM 2.5 concentration emitted from the pollutant source i at the position j When I ij > 1, the higher the equal contribution index of standard concentration is, and the more the PM 2.5 concentration exceeds the standard. The larger the r ij is, the smaller the pollution risk of this single source i is compared with other pollution sources. When the r ij is smaller, the relative pollution risk of the single source i pollution source is larger. Considering the relative contribution ratio of superimposed sources can identify the sources with a high contribution ratio of concentration. The change rate of PM 2.5 concentration will be affected by the surrounding pollution source distribution, terrain, wind direction, and other factors. k ij ' represents the change rate of PM 2.5 concentration emitted from the pollutant source i with migration at position j. When k ij ' > 0, for the same pollution sources, the larger k ij ' is, the faster PM 2.5 concentration decays with diffusion distance, the smaller the transmission distance is, and the diffusibility decreases accordingly.

Air Pollution Risk Assessment and Management
Taking into account the concentration emitted from pollution sources, the relative contribution ratio of superimposed sources, and the change rate of PM 2.5 concentration emitted from the pollutant source, the PM 2.5 environmental risk assessment was carried out using the following methods: R(E ij ) represents the PM 2.5 environmental risk at the location j caused by pollution source i. R(P ij ) represents the relative PM 2.5 pollution risk at the location j caused by pollution source i. When the I ij is larger, the r ij is smaller, and the contribution proportion of concentration is larger, so the relative pollution risk is higher. R(D ij ) represents the PM 2.5 diffusion risk at the location j caused by pollution source i, and it depends on the emission concentration and the change rate of PM 2.5 concentration with migration. When I ij is fixed, k ij '(h) > 0, the smaller k ij '(h) is, the larger the migration and diffusion distance is, and the higher the diffusion risk is.
The risk values were graded, the risk areas were divided, and the corresponding control methods were proposed.

PM 2.5 Emission Inventory Calculation and Spatial Distribution
The PM 2.5 emission inventory is shown in Table 5. According to the proportion of PM 2.5 emission intensity, the emission sources were sorted as non-paved roads > coal piles > mine pits and dumps > paved roads > quarries > loading and unloading. Among them, the proportion of non-paved roads was 64.64%. According to the coal mine operation sites, the pollution sources were distributed in strips, concentrated in the west and southwest of Zhuozishan Mountain. There were 38 area sources and 6 linear sources ( Figure 2). The highest contribution came from non-paved roads, which is similar to previous studies. Huertas et al. [14,38,39] all showed that unpaved roads were the primary source of emissions, with PM 2.5 emissions accounting for 78.4%, 34%, and 60% of the total emissions, respectively, followed by transmission, loading and unloading, wind erosion, and other activities. The second major contribution source in this study was wind erosion on the surface of the coal piles. Due to different mining properties, geological conditions, management, mining technologies, and climatic conditions, the PM 2.5 emissions emitted from the same activity will vary; still, on the whole, unpaved transportation was the main emission source.

PM 2.5 Concentration Distribution Simulation
The linear relationship between the daily mean concentration of PM 2.5 at the monitoring points and the simulated values is shown in Figure 3. The calculations of the evaluation indexes are shown in Table 6. Since the average daily PM 2.5 concentration monitored occurred in discrete discontinuous days, the index calculations were aggregated. The data of the three monitoring points all showed good consistency, and the average FVC2 value reached 100%. The FB ranged from 0.08 to 0.31, with an average of 0.19; the MG ranged from 1.09 to 1.29, with an average of 1.2. The FB and MG showed that the difference between the predicted value and the monitored value was within the acceptable range, but the predicted value was underestimated. The slope of the scatter plot also proved this. This is most likely because the emission inventory did not include all pollution sources. Among them, MS2 had the largest deviation, probably because it was near the road and had a large traffic flow, leading to secondary dust emergence, which was not taken into account in the emission inventory. The VG and NMSE were also within the acceptable range. The VG of MS2 was 1.11, and the NMSE was 0.1, both of which were the highest values. This may be because the emission inventory did not consider the changes in a time series. However, in practice, the emission intensity of PM 2.5 sources was constantly changing in the heating season and non-heating season and in different periods of the day. MS2 was located near the road sources, and the emission intensity of the pollution source fluctuated at a high frequency, so the deviation was large. MS3 was located near the quarries, with fewer pollution sources and stable emission intensity, so there was little interference. MS1 was located in an area with dense pollution sources, and the intensity variation of surrounding pollution sources could cause deviation of the simulation, so the R 2 was low. In addition, errors in the model itself could also lead to bias. Overall, the inventory and the model apply to the PM 2.5 concentration simulation of the non-point sources in the study area, and the simulation can be used for a pollution risk assessment.
We selected the typical meteorological day to simulate the daily mean concentration and make an environmental risk assessment. The typical meteorological day is the one with the heaviest pollution under the most common weather conditions. The statistics of regional meteorological data in 2018 showed that the dominant wind direction was WNW and the average wind speed was 3.39 m/s (Table S5). The most serious PM 2.5 pollution period was concentrated in December, and the maximum concentration occurred in the condition of atmospheric stability D. We selected December 3 as the typical meteorological day (Table 7), and the atmospheric stability grade of this day is D.

Concentration Contribution of Single Pollution Sources
The equal contribution indexes of standard concentration of six types of emission sources are shown in Figure 4. The PM 2.5 concentration distribution is significantly affected by the wind direction and spreads to the southeast, and the concentration contributions of six sources mainly depend on the emission intensities and the total emissions. The non-paved road source was the primary contribution source, followed by the exposed surface source of coal piles, and then the exposed surface source of mine pits and dumps, the maximum equal contribution indexes of which were 2.40, 2.21, and 2.10, respectively. The areas exceeding the standard were mainly concentrated near the emission sources, and the maximum equal contribution indexes of the other emission sources were all less than 1. The maximum equal contribution index of the multi-source superposition was 4.40.
Combined with the distribution of pollution sources, we found that the PM 2.5 concentration emitted by some dispersed high-intensity sources was lower than that of dense low-intensity sources. Ghannam [40] also mentioned that the difference in geometric shape and spatial distribution of pollution sources do not allow for correlation between emission magnitudes and contribution to ambient pollutant concentrations. In order to further clarify the relationship between the spatial distribution of pollution sources and emission concentration, we selected one of the coal piles and set three scenarios. Among them, the total area, emission intensity, altitude, emission height, and other conditions of pollution sources were all the same (Table S6), but the distribution of pollution sources ranged from dense to sparse (see Figure 5). The simulation showed that in terms of the maximum equal contribution index of standard concentration, Scenario 1, Scenario 2, and Scenario 3 were 3.68, 1.18, and 0.63 respectively, and the over standard areas of Scenario 1 and Scenario 2 were 8km 2 and 6km 2 , respectively. It shows that the more concentrated the pollution sources are, the higher the maximum ground concentration will be, and the larger the area exceeding the standard will be under the same source intensity and emissions. This is because the overlapping degree of the diffusion surface from each sub-pollution source becomes higher. Dispersed pollution sources can reduce pollution levels. In addition, we found that when the concentration contributions of the six types of sources were summed on the spatial grid chart, the result was basically equal to the emission concentration in the multi-source superposition scenario. Unlike CFD (Computational Fluid Dynamics) [41,42], CALPUFF is not a model based on particle simulations and does not take into account particle interference from different sources [43], so in a scenario with the same source intensity and source distribution, the simulations at the same point have a linear superposition law [40]. However, in the actual situation, fine particles emitted from multiple non-point sources will collide or agglomerate in the diffusion process due to turbulent diffusion, molecular diffusion, electrostatic attraction, chemical transformation, and other reasons, and interfere with each other, so the particle diffusion path is different from that of single-source emissions [44][45][46][47]. The emission concentration after the multi-source superposition may be smaller than the direct arithmetic summation of single sources. Therefore, when the equal contribution index of the standard concentration of the multi-source superposition is larger, the actual risk may be smaller than the calculation.

Relative Pollution Risk after the Multi-Source Superposition
The relative contribution ratios of superimposed sources for six types of sources are shown in Figure 6. The primary relative pollution risk contribution source was non-paved road sources, with an average value of 2.34 × 10 −2 and a maximum value of 1.46. The least was the loading and unloading source, with an average value of 2.86 × 10 −6 and a maximum value of 1.92 × 10 −4 . Concentrated sources with the same source intensity and emissions are more likely to have high pollutant concentrations, so the pollution risk is higher than that of dispersed sources. On the whole, the order of the relative risk of the six sources was consistent with the order of concentration contribution. However, for the proportion of the relative pollution risk of the single sources, unpaved roads and quarries increased by 15.5% and 1.0%, while the other sources decreased, compared with the absolute pollution risks without superimposed sources. This is because the source intensity and emission area of non-paved roads were the largest among the six types of pollution sources, so the concentration contribution ratio was the highest. Although the overall concentration level of the quarries was low, it was in the upwind direction, and there were few other disturbing pollution sources in the region. Therefore, the relative contribution ratio of the quarries to PM 2.5 concentration in the region was higher than that of other pollution sources, and the risk increases accordingly. For multiple sources, it is more practical to calculate the relative pollution risk in a certain meteorological condition than the absolute pollution risk of a single source. When the influence of other superimposed sources is considered in the relative pollution risk assessment, the risk of high contribution sources in the low-concentration pollution area is increased, and the risk of low-contribution sources in the high-concentration pollution area is reduced.

Spatial Transport Law of Dust Particles in Multi-Scale Mining Areas
The diffusion of pollutants is affected by wind direction, wind speed, temperature, terrain, source intensity, and other factors [48,49], among which the meteorological and terrain factors have the most significant influence, but they are difficult to control artificially. Feng Liu's study [50] mentioned that the concentration distribution of pollutants is not only related to meteorological factors but is also largely affected by the relative position between sources, and this effect is caused by the change of wind direction. In order to clarify the relationship between the distribution shape of pollution sources and emission concentration and diffusion distance, we set up three comparative scenarios. In the three scenarios, the emission intensity, area, altitude, effective height, and other conditions of pollution sources were all the same (Table S6), but the distribution direction of pollution sources was different (Figure 7). The pollution source in Scenario 1 was most concentrated. In Scenario 4, the pollution source was distributed laterally, and the angle with the dominant wind direction was the largest, and the windward side was the largest. In Scenario 5, the pollution source was distributed longitudinally, and the angle with the dominant wind direction was the smallest, and the windward side was the smallest. The simulation showed that for the pollution sources with the same emission intensity and source area, in terms of the maximum ground concentration, Scenario 1 > Scenario 5 > Scenario 4; in terms of the area exceeding the pollution standard, Scenario 4 > Scenario 1 > Scenario 5, which were 9 km 2 , 8 km 2 , and 6 km 2 respectively. When the angle between the distribution direction of the linear pollution sources and the dominant wind direction is smaller, the emission concentration is higher, but the diffusion surface is smaller. This is because the particles diffuse along the wind direction; the particles emitted by the upstream pollution sources and the particles emitted by the downstream pollution sources accumulate and superposition continuously, resulting in the increase of the maximum ground concentration, and the diffusion surface is closer to the arrangement direction of the pollution sources. However, when the angle between the direction of the linear pollution source distribution and the dominant wind direction is larger, the maximum ground concentration is lower, but the diffusion surface is larger because the diffusion directions of particles emitted by each sub-pollution source are almost parallel to each other and the superposition interference is weaker. A concentrated distribution of pollution sources spread in the middle, but the emission concentration is the highest, so the area of excessive pollution is more than that of Scenario 5. On the whole, the influence of pollution source distribution direction on concentration and diffusion is not as obvious as that of the meteorological factors. Figure 8 shows the change rates of PM 2.5 concentration with migration in six types of sources and multi-source superposition scenarios under the meteorological conditions of a typical day. On the whole, the change rate of PM 2.5 concentration with migration was the highest in the multi-source superposition scenario, with an average value of 6.62 × 10 −2 and a maximum value of 3.92. Among the six types of sources, the change rate of PM 2.5 concentration of unpaved roads with migration was the largest, and the change rate of PM 2.5 concentration of the quarries with migration was the smallest. The transport distance of particulate matter is determined by the emission concentration and the change rate of PM 2.5 concentration with migration. Therefore, the calculations of diffusion risk showed that among the six types of sources, non-paved roads had the largest diffusion risk, with an average risk value of 3.28 × 10 −2 and a maximum value of 2.37; second, coal piles, with an average of 1.42 × 10 −2 and maximum of 1.85; the diffusion risk of the mine pits and dumps was close to that of the coal pile sources, with an average value of 1.27 × 10 −2 and a maximum value of 1.65; loading and unloading sources had the lowest diffusion risk, and the maximum value was 0.02. Combined with the emission intensity, emission area, and distribution, we found that the six types of pollution sources were in a north-south linear arrangement macroscopically. For the non-paved road sources, a large number of large-scale sources with moderate emission intensity were concentrated in the north-south direction, while a small number of small-scale sources with high emission intensity were scattered in the southeast corner. The non-paved road sources were distributed in a wide range, and the particulate matter emitted by the pollution sources in the upwind direction overlapped with the particulate matter emitted by the pollution sources in the downwind direction, increasing the pollution concentration and the diffusion area, so the diffusion risk was the highest. The distribution of loading and unloading sources was the same as that of the non-paved roads, but the source intensity was the smallest, so the transmission distance of particles was also the smallest, and the diffusion risk was the lowest. The distribution of coal piles was similar to that of non-paved roads. There were both concentrated linear sources and scattered sources of high emission intensity. However, the source intensity and total pollution source area were smaller than those of the non-paved roads, so the diffusion risk was the second. For mine pits and dumps, the distribution was relatively dispersed macroscopically, while, locally, it was densely distributed. Although the diffusion surface was large, the concentration exceeded the standard only in the vicinity of the area with dense pollution sources, so there was a greater diffusion risk only in the downwind direction of the area with dense pollution sources.

Air Pollution Risk Assessment and Management
Taking the relative pollution risk and the diffusion risk of emission sources into consideration, based on Equation (14), the environmental risks of the six types of emission sources are shown in Figure 9. The order of the environmental risk contribution of the six types of emission sources from large to small is the following: non-paved road > coal piles > coal pits and dump > paved road > quarries > load and unloading. The average environmental risks of the top three pollution sources are 5.61 × 10 −2 , 2.06 × 10 −2 , and 1.89 × 10 −2 , respectively, and the average environmental risks of loading and unloading sources are 3.66 × 10 −4 . In the multi-source superposition scenario, emission sources in the red frame area (see Figure 10) were densely distributed and arranged in a north-south linear manner. Due to the high PM 2.5 concentration and the angle between the distribution direction of sources and the dominant wind direction, both pollution risk and diffusion risk were higher, so the environmental risk was higher. In particular, the downstream pollution sources along the linear alignment direction, due to the PM 2.5 concentration accumulation of upstream pollution sources, the pollution risk would continue to increase. In the blue frame area (see Figure 10), the emission sources were distributed sparsely and irregularly. The pollution risk was lower, but the diffusion risk was higher, which decreased with the increase of the diffusion distance, and the overall environmental risk is lower. Regarding previous studies [51,52], we formulated the risk classification standard of the mining area (Table 8) and divided the mining area into three categories (Figure 10). High-risk areas are mainly concentrated areas of emission sources, which are heavily polluted areas. Therefore, it is necessary to significantly reduce the intensity of high-emission sources in this area, reduce the density of pollution sources and adjust the layout of pollution sources. The medium-risk areas belong to the moderate and light pollution areas where the pollution sources are dispersed with high diffusion risk. Appropriate reductions in the intensity of emissions sources are needed, and the diffusion risk can be reduced by installing wind-proof and dust-suppressing nets upstream and downstream of the dominant wind direction (west, northwest, and north), or by distributing the easily mobile pollution sources, such as unpaved roads, along the dominant wind direction. Low-risk areas are mainly those with low pollution contribution and low dispersion, where PM 2.5 concentration is within the safety standard and its impact on the environment is acceptable. Conventional management can be implemented. According to the concentration distribution of the multi-source superposition, it can be seen that the middle and high-risk areas can completely cover the areas exceeding the PM 2.5 standard, and we can avoid the one-size-fits-all control. It is not necessary to reduce the intensity of the source altogether. Instead, we can reduce the risk and the management cost by reducing the distribution density of pollution sources and adjusting the layout of pollution sources.

Significance and Limitations
In the current research on the environmental risk assessment of particulate matter, when mathematical models are used to simulate the diffusion of particulate matter of pollution sources, the various assumed ideal conditions will bring some practical deviations. Therefore, in this study, with the concentration of simulation model for environmental risk assessment, we constructed an index system based on different scenarios, highlighted the environmental risks of sources with high pollution concentration contribution and high diffusibility to more easily and accurately assess and control the risks of non-point sources in mining areas.
However, in the compilation of the fine particulate matter emission inventory in the mining area, due to the limitation of conditions, the pollution sources such as coal spontaneous combustion and secondary dust emissions were not counted. Moreover, meteorological conditions were not fixed. The emission calculation formulae were adjusted for different meteorological conditions, and the emission inventory had not been allocated in a more detailed time series. These factors could lead to the uncertainty of the inventory. In addition, the particle interference from different sources is not considered in the model CALPUFF simulation, so when we use this model for environmental risk assessment, we need to further consider the particle interference of multiple sources in the future and incorporate more comprehensive risk factors, which is expected to improve a more scientific, standardized, and accurate risk assessment system.

Conclusions
The PM 2.5 emission inventory was obtained from the study, which indicated that among the key control types of non-point sources in the mining area, unpaved road transport was the primary emission source, followed by the wind erosion on the exposed surface of coal piles, and the third was the wind erosion on the exposed surface of mine pits and dumps. The simulation and calculations showed that the maximum equal contribution index of the standard concentration of these three pollution sources were 2.40, 2.21, and 2.10, respectively. The maximum equal contribution index of the standard concentration of the multi-source superposition was 4.40. The simulation showed that the wind field and the spatial distribution of emission sources could affect the effect of superimposed pollution sources, thus affecting the risk. The denser the distribution of pollution sources, the higher the emission concentration, the smaller the relative contribution ratio of superimposed sources for a single source, and the greater the relative pollution risk. Dispersed pollution sources can reduce the pollution risk by reducing the concentration level. The relative pollution risk ranking of the six types of sources was consistent with the ranking of the concentration contribution level. Among them, non-paved road sources had the largest relative pollution risk, with an average risk value of 2.34 × 10 −2 and a maximum risk value of 1.46. However, for the proportion of the relative pollution risk of single sources, unpaved roads and quarries increased by 15.5% and 1.0%, while other sources decreased, compared with the absolute pollution risks without superimposed sources. When the influence of superimposed sources is considered in the relative pollution risk, the risk of high contribution sources in the low-concentration pollution area is increased, and the risk of low-contribution sources in the high-concentration pollution area is reduced.
In the case of the same source intensity and emission area, the smaller the angle between the linear pollution source arrangement direction and the dominant wind direction, the higher the cumulative emission concentration, but the smaller the diffusion surface. The transport distance of particulate matter is determined by the emission concentration and the change rates of PM 2.5 concentration with migration. In the multi-source superposition scenario, the change rate of PM 2.5 concentration with migration was the highest, with an average value of 6.62 × 10 −2 and a maximum value of 3.92. Among the six types of sources, the change rate of PM 2.5 concentration with migration of the unpaved roads was the least, and that of the quarries was the largest. Considering the influence of concentration, the diffusion risk of unpaved roads was the highest, coal piles were the second, mine pits and dumps were the third, and loading and unloading sources were the lowest; the average risks are 3.28 × 10 −2 , 1.42 × 10 −2 , 1.27 × 10 −2 , and 3.63 × 10 −4 , respectively. Diffusion risk could better evaluate the potential of pollution risk and realize a more accurate classified control.
Combined with the relative pollution risk and diffusion risk of emission sources, the environmental risk of non-paved roads was the highest, followed by coal piles; the third is mine pits and dumps, and that of loading and unloading sources is the lowest, and their average environmental risks are 5.61 × 10 −2 , 2.06 × 10 −2 , 1.89 × 10 −2 and 3.66 × 10 −4 , respectively. At the same time, different control measures were proposed to realize the risk control by zoning and classification.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/su13126619/s1, Table S1: Multiyear meteorological conditions in Qipanjing town (1988-2018); Table S2: The formulas, parameter values, and significances of emission factors and total emissions of particulate matter emission sources in the mining area; Table S3: Data of production and operation in the mining area; Table S4: The sources of the model data and Settings; Table S5: Regional wind speed and wind frequency statistics in 2018; Table S6: Parameters of pollution sources in different scenarios for simulation.