Human Health Risk Distribution and Safety Threshold of Cadmium in Soil of Coal Chemical Industry Area

: Cadmium (Cd) is a highly carcinogenic metal that plays an important role in the risk management of soil pollution. In this study, 153 soil samples were collected from a coal chemical plant in northwest China, and the human health risks associated with Cd were assessed through multiple exposure pathways. Meanwhile, by the Kriging interpolation method, the spatial distribution and health risks of Cd were explored. The results showed that the average concentration of Cd in the soil was 0.540 mg/kg, which was 4.821 and 5.567 times that of the soil background value in Ningxia and China, respectively. In comparison, the concentration of Cd in the soil was below the national soil environmental quality three-level standard (1.0 mg/kg). In addition, health risk assessment results showed that the total carcinogenic risk of Cd was 1.269 × 10 − 6 –2.189 × 10 − 6 , both above the acceptable criteria (1 × 10 − 6 ), while the hazard quotient was within the acceptable level. Oral intake and ingestion of soil particles were the main routes of exposure, and the carcinogenic risk control value of oral intake was the lowest (0.392 mg/kg), which could be selected as the strict reference of the safety threshold for Cd in the coal chemical soil. From Kriging, a prediction map can be centrally predicted on heavy metal pollution in the area surrounding the coal entrance corridor and pedestrian entrance. This study can provide a theoretical basis for the determination of the heavy metal safety threshold of the coal chemical industry in China.


Introduction
The development of the modern coal chemical industry has accelerated the clean and efficient utilization of coal resources, especially in coal gasification [1][2][3]. In recent years, the development of the modern coal chemical industry is mainly concentrated in South Africa, the United States, and China [4]. However, the coal gasification process is often accompanied by the transfer and transformation of heavy metals and other harmful trace elements [5]. Because of the nonbiodegradability, toxicity, and cumulative properties of heavy metals, they can cause environmental pollution and threaten human health after being transferred to the soil, groundwater, air, and other environmental media [6,7]. From 2005 to 2013, the Ministry of Environmental Protection (MEP) and the Ministry of Land and Resources (MLR) of the People's Republic of China jointly carried out the first national soil pollution survey [8], but the overall pollution situation was relatively serious. According to the data of the National Soil Pollution Survey Bulletin in 2014, the standard exceeding ratios of soil heavy metal pollution reached 16.1% in total [9]. Among them, the standard exceeding ratio of cadmium (Cd) was 7.00%, and the proportion of heavy pollution points was 0.50% [10]. In China, soils in more than 11 provinces and 25 regions are rich in Cd [11]. As a worldwide environmental problem, Cd was listed in seventh place as a toxic substance of concern by the American Agency for Toxic Substances and Disease Registry (ATSDR). Moreover, it was listed as a highly toxic, hazardous, and carcinogenic substance The study area is a modern coal chemical plant in northwest China, which is located in Zone A of an energy chemical base. It has a typical continental monsoon climate, with an annual average temperature of 9.0 • C. Southerly wind prevails in summer and northerly wind prevails in winter. The surface water around the plant area is deficient. The annual average precipitation is 199.5 mm, and the annual evaporation is 1752.6 mm, much larger than the former [28]. The physical weathering and wind force are significant. In terms of vegetation, arid grassland is rare in this region. The main soil types are calcareous soil and aeolian sandy soil with poor development. Many days with strong winds exist in the winter and spring dry seasons. All of the above factors are likely to cause the migration of the surface soil polluted by heavy metals and then expand the hazardous area of soil heavy metals pollution.
The plant area occupies 400,000 square meters. It mainly consists of the following areas: a complete modification unit, a variable power distribution unit, an air separation and pressure unit, a water treatment unit, a methanol synthesis unit, a dimethyl ether synthesis unit, a gasification unit, a power unit, a product storage unit, a coal entry corridor, a railway loading/unloading area, two slag dumps, a chimney, a road for staff entry and exit, a road for logistics entry and exit, and main roads in the factory area. The geographical location of the study area and the regional distribution of the plant are shown in Figure 1.

Study Area
The study area is a modern coal chemical plant in northwest China, which is located in Zone A of an energy chemical base. It has a typical continental monsoon climate, with an annual average temperature of 9.0 °C . Southerly wind prevails in summer and northerly wind prevails in winter. The surface water around the plant area is deficient. The annual average precipitation is 199.5 mm, and the annual evaporation is 1752.6 mm, much larger than the former [28]. The physical weathering and wind force are significant. In terms of vegetation, arid grassland is rare in this region. The main soil types are calcareous soil and aeolian sandy soil with poor development. Many days with strong winds exist in the winter and spring dry seasons. All of the above factors are likely to cause the migration of the surface soil polluted by heavy metals and then expand the hazardous area of soil heavy metals pollution.
The plant area occupies 400,000 square meters. It mainly consists of the following areas: a complete modification unit, a variable power distribution unit, an air separation and pressure unit, a water treatment unit, a methanol synthesis unit, a dimethyl ether synthesis unit, a gasification unit, a power unit, a product storage unit, a coal entry corridor, a railway loading/unloading area, two slag dumps, a chimney, a road for staff entry and exit, a road for logistics entry and exit, and main roads in the factory area. The geographical location of the study area and the regional distribution of the plant are shown in Figure 1.

Soil Sampling and Analysis
To comprehensively and accurately evaluate the pollution degree of soil heavy metals at different locations in the plant, regional soil samples were collected through a 50 m × 50 m checkerboard distribution method in this study. A total of 153 soil samples from the surface layer (0-20 cm) were collected. In the process of sampling, some sample points in pools, workshops, slag piles, and other areas, were detrimental to soil collection. Therefore, it was necessary to make appropriate adjustments according to the actual environment around the preset sampling points. Sampling information is shown in Figure 2. The

Soil Sampling and Analysis
To comprehensively and accurately evaluate the pollution degree of soil heavy metals at different locations in the plant, regional soil samples were collected through a 50 m × 50 m checkerboard distribution method in this study. A total of 153 soil samples from the surface layer (0-20 cm) were collected. In the process of sampling, some sample points in pools, workshops, slag piles, and other areas, were detrimental to soil collection. Therefore, it was necessary to make appropriate adjustments according to the actual environment around the preset sampling points. Sampling information is shown in Figure 2. The collected soil samples were packed in sealed polyethylene bags, labeled, sent to the laboratory, and dried in a natural environment. After picking out the stone, animal, and plant debris, the samples were ground in an agate mortar and passed through a 100-mesh nylon sieve.
collected soil samples were packed in sealed polyethylene bags, labeled, sent to th atory, and dried in a natural environment. After picking out the stone, animal, an debris, the samples were ground in an agate mortar and passed through a 100-mes sieve.
The content of Cd in soil was determined by graphite furnace atomic absorpti trometry with an AAS-ZEEnit 700 Atomic Absorption Spectrometer (Jena, Germa detection limit was 0.006 mg/kg. A 0.300 g sample was accurately weighed and pu beaker. Ten milliliters of concentrated HNO3, 4ml 3% (mass fraction) H2O2, and 10 were added in turn for microwave digestion. After digestion, 1% (mass fraction was added to a 50 mL colorimetric tube for standby (GB/T 17141-1997).

Quality Control
To control quality, each batch of samples used the same reagents and follo same steps to make 2 reagent blanks. Into each batch of samples was added a na certified reference material (GSS-10 and GSS-11 from the Ministry of Environmen tection Standard Sample Institute). The relative standard deviation was usually  The detection limit was 0.006 mg/kg. A 0.300 g sample was accurately weighed and put into a beaker. Ten milliliters of concentrated HNO 3 , 4ml 3% (mass fraction) H 2 O 2 , and 10 mL HF were added in turn for microwave digestion. After digestion, 1% (mass fraction) HNO 3 was added to a 50 mL colorimetric tube for standby (GB/T 17141-1997).

Quality Control
To control quality, each batch of samples used the same reagents and followed the same steps to make 2 reagent blanks. Into each batch of samples was added a nationally certified reference material (GSS-10 and GSS-11 from the Ministry of Environmental Protection Standard Sample Institute). The relative standard deviation was usually within ±2%; otherwise, the sample was measured again. The measurement results of the blank samples were all less than the detection limit of the method, and the recovery rate of the matrix spiked was in the range of 80-130% [29].

Human Health Risk Assessment
The Ministry of Environmental Protection of the PRC has released the Technical Guidelines for Risk Assessment of Soil Contamination of Land for Construction [27]. The human health risk assessment models recommended in the Technical Guidelines are based on the EPA models. However, in order to reflect the actual contamination situations in China, a suitable parameter is given according to the environment and living habits of Chinese residents. Based on the pollution characteristics of coal chemical sites, the site type was verified, the exposure route selected, and the exposure amount, carcinogenic risk, hazard quotient, and risk control value calculated. On this basis, the contribution rate of Cd carcinogenic risk under different exposure routes was analyzed, and the environmental safety threshold of Cd in the soil of the study area was determined.

Exposure Pathways
The study site is a coal chemical industry area, a non-insensitive industrial land, with no surface water around and no groundwater as drinking water. In this paper, oral intake, skin contact, and inhalation of soil particles were selected as the main exposure routes to conduct a human health risk assessment of Cd. The soil exposure model (Table 1) and exposure factor parameters ( Table 2) corresponding to carcinogenic and noncarcinogenic effects of a single pollutant were selected [25]. There are no children in the coal chemical plant area, so children's exposure was not considered, and only adults were selected as recipients. Table 1. Calculating models of soil exposure dose in three soil exposure routes.

Exposure Routes Explanation Formula Expression of Exposure
Mouth intake Carcinogenic Inhalation of soil particles Note: In Table 1, OISER ca soil exposure dose in oral intake (carcinogenic) in milligrams per kilogram per day, OISER nc soil exposure dose in oral intake (noncarcinogenic) in milligrams per kilogram per day, DCSER ca soil exposure dose in skin contact (carcinogenic) in milligrams per kilogram per day, DCSER nc soil exposure dose in skin contact (noncarcinogenic) in milligrams per kilogram per day, PISER ca soil exposure dose in inhalation (carcinogenic) in milligrams per kilogram per day, PISER nc soil exposure dose in inhalation (noncarcinogenic) in milligrams per kilogram per day. Reference dose distribution coefficient exposed to soil Note: In Table 2, the exposure parameters only targeted adults.

Risk Characterization and Contribution Rate
Human health risk brought about by Cd pollution was characterized by a single pollutant risk. The risk characterization and contribution rate can be calculated by the formula of carcinogenic risk and hazard quotient of different soil exposure routes (Table 3) [30].  Table 3, CR ois -carcinogenic risk of mouth-intake soil; CR dcs -carcinogenic risk of skin-contact soil; CR pis -carcinogenic risk of soil particle inhalation; HQ ois -hazard quotient of mouth-intake soil; HQ dcs -hazard quotient of skin-contact soil; HQ pis -hazard quotient of soil particle inhalation.
Based on the Cd pollution risk assessment, the contribution rates to the risk of the three different routes were calculated (Formula (1)), and the main route was analyzed to obtain evidence for prevention and control of the subsequent risk. In the above formula: CR i means the carcinogenic risk contribution rate or hazard quotient level of a certain exposure route, and the dimension is 1; ΣCR i means the total carcinogenic risk or total hazard quotient.

Risk Control Thresholds
The acceptable risk value of carcinogens defined by the USEPA is that the lifetime risk of cancer exceeds the normal value of 1 × 10 −4 -1 × 10 −6 . When pollutant hazard exceeds the acceptable level of human carcinogenic risk (1 × 10 −6 ) or hazard quotient (1), the risk control value should be calculated corresponding to its soil exposure route. The calculation formula of the safety threshold of three kinds of soil exposure paths can be shown in Table 4 [30].  Table 4, RCVSois-soil safety threshold based on carcinogenic effect of mouth-intake soil pathways, mg/kg; ACR-acceptable carcinogenic risk, which is dimensionless and takes the value of 10 −6 ; AHQ-acceptable hazard quotient, which is dimensionless and takes the value of 1; RCVSdcs-soil safety threshold based on carcinogenic effect of skin-contact soil pathways, mg/kg; RCVSpis-soil safety threshold based on carcinogenic effect of soil particle inhalation pathways, mg/kg.

Kriging Interpolation Method
The Kriging interpolation of the GIS spatial model is a geostatistical method that is used in smoothing surfaces and predicting the values of unsampled locations. It provides the best linear unbiased estimation and estimation error distribution information, which shows a strong statistical advantage [31,32]. Its formula is expressed as [33]: In the Formula (2), Z (x0) represents the concentration of heavy metal in the estimation points, and Z (xi) represents the concentration of heavy metal in the sampling point i. n represents the total number of sampling points, and λ i is a set of weight coefficients, the value of which is to ensure minimum variance of the result and an unbiased estimate.

The Descriptive Statistics of Cd
The statistical characteristic values of Cd content in the soil of the coal chemical plant area are listed in Table 5. It was measured that the concentration of Cd in the soil of the coal chemical factory was 0.400~0.690 mg/kg, with an average concentration of 0.540 mg/kg. According to the total Cd concentration standard of 20 mg/kg (industrial land) in the soil environmental quality standard GB15618-2008 [34], the Cd concentration did not exceed the standard. However, compared to the national background value of 0.097 mg/kg and the Ningxia background value of 0.112 mg/kg, the rate of exceeding the standard in all production units reached 100%. From the spatial distribution of Cd content in the soil, the coefficient of variation was small, with an average value of 11.043%, indicating that the spatial distribution of Cd in the soil was uniform and the degree of dispersion was small.

Spatial Autocorrelation Test
The spatial autocorrelation (Moran's I) tool of ArcGIS spatial statistics was used to analyze the carcinogenic risk and hazard quotient of Cd pollution. The results are shown in Table 6. It can be seen from Table 6 that the value of Moran's I was relatively close to 1, so the hazard quotient of Cd was a significant cluster with a positive correlation [35]. By the criterion of spatial autocorrelation, compared with the z value and p value in Table 6, it was indicated that the spatial autocorrelation was very significant. Therefore, simulation analysis could be used in the spatial interpolation method of Kriging.

Analysis of Spatial Structure Variation
Semivariance analysis is the key to spatial variation structure analysis and plays an important role in reasonably controlling the number of sampling points and improving the accuracy of interpolation. When choosing the semivariance function model, the closer the error standard average of the cross-validation result is to 0, the closer the root mean square of error standard is to 1, which means the higher precision of the simulation of the model [36]. Table 7 showed the fitting results of the semivariance model of total carcinogenic risk and total hazard quotient of Cd in Kriging interpolation points. From the ratio of nugget value C 0 to the sill value (C 0 /Sill), the interpolated value of Cd was ranged from 33.3 to 34.64% and between 25 and 75% (Table 7), which was a medium correlation, indicating that the influence of regional factors on the variation of soil Cd content in the chemical plant was significantly greater than that of nonregional factors (i.e., random factors).

Human Health Risk Assessment of Soil Cd in Coal Chemical Area
The risk values of soil Cd in the coal chemical industry area were calculated according to oral intake, skin contact, and inhalation of soil particles. The total carcinogenic risk value and total risk quotient were calculated by summing. The results are shown in Table 8. The Kriging interpolation method was used to simulate the multipath total carcinogenic risk value of Cd in the soil in Table 8, and the total carcinogenic risk assessment map of the whole plant area was obtained (Figure 3). According to the results in Table 8 and Figure 3, the multipath combined carcinogenic risk value of Cd in the soil area of the chemical unit in the plant area was in the range of 1.269 × 10 −6 -2.189 × 10 −6 , with an average value of 1.754 × 10 −6 , which exceeded the acceptable level of human health (1 × 10 −6 ).

Human Health Risk Assessment of Soil Cd in Coal Chemical Area
The risk values of soil Cd in the coal chemical industry area were calculated according to oral intake, skin contact, and inhalation of soil particles. The total carcinogenic risk value and total risk quotient were calculated by summing. The results are shown in Table  8.

Carcinogenic Risk
The Kriging interpolation method was used to simulate the multipath total carcinogenic risk value of Cd in the soil in Table 8, and the total carcinogenic risk assessment map of the whole plant area was obtained (Figure 3). According to the results in Table 8 and Figure 3, the multipath combined carcinogenic risk value of Cd in the soil area of the chemical unit in the plant area was in the range of 1.269 × 10 −6 -2.189 × 10 −6 , with an average value of 1.754 × 10 −6 , which exceeded the acceptable level of human health (1 × 10 −6 ).  It can also be seen from the figure that the spatial distribution of the pollution area was uneven. The high pollution area was mainly divided into two parts. One part was near the coal entry corridor. Mainly heavy-duty transport vehicles pass through the road system in the mining area, producing a large amount of dust, and the transported mineral materials may be dispersed to a certain extent, constituting an important source of heavy metals, especially lead, nickel, and Cd [37]. The content of heavy metals in soils located along roads is strongly related to traffic and decreases with increasing distance from the road [38].
Cd pollution in the soil around the coal mine was caused by coal combustion, coal mining, the coal transportation process, and the migration and settlement of coal gangue and coal dust. This is consistent with our results [39,40].
Another part of the high pollution area was near the entrance and exit for people. The accumulation of Cd in the soil was related to the common effects of coal mineralization and accumulation and human disturbance. Anthropogenic sources account for more than 90% of the Cd released from the environment [41]. Research by Xu et al. showed that Cd pollution levels were relatively high in areas with intensive human activities in the coal mine area while exploring Cd's harm to human health [42]. The risk value decreases from the northwest and southeast to the center and surroundings. There is a chimney in the power unit, so Cd settles down with the smoke and dust, which leads to a higher cancer risk of Cd around the chimney and in the southwest in a phenomenon consistent with plume dispersion patterns as reflected in the regulatory understanding of atmospheric emissions [43]. Regulatory guidance suggests carrying out harmless treatment of the slag of the total transformation unit and the methanol synthesis unit and controlling the emission of flue gas to reduce or block the potential carcinogenic risk.

Hazard Quotient
The Kriging interpolation method was used to simulate the multipath total hazard quotient of soil Cd in Table 8, and the hazard quotient evaluation map of the entire plant area was obtained ( Figure 4). It can be seen from Table 8 and Figure 4 that the Cd multipath hazard quotients in the soil area of the chemical units were 0.014-0.024 (with an average of 0.020), all of which were within the acceptable level of 1. The main polluted areas were those near the coal transportation corridor, power unit, complete modification unit, and the entry/exit for people.
It can also be seen from the figure that the spatial distribution of the pollution area was uneven. The high pollution area was mainly divided into two parts. One part was near the coal entry corridor. Mainly heavy-duty transport vehicles pass through the road system in the mining area, producing a large amount of dust, and the transported mineral materials may be dispersed to a certain extent, constituting an important source of heavy metals, especially lead, nickel, and Cd [37]. The content of heavy metals in soils located along roads is strongly related to traffic and decreases with increasing distance from the road [38]. Cd pollution in the soil around the coal mine was caused by coal combustion, coal mining, the coal transportation process, and the migration and settlement of coal gangue and coal dust. This is consistent with our results [39,40].
Another part of the high pollution area was near the entrance and exit for people. The accumulation of Cd in the soil was related to the common effects of coal mineralization and accumulation and human disturbance. Anthropogenic sources account for more than 90% of the Cd released from the environment [41]. Research by Xu et al. showed that Cd pollution levels were relatively high in areas with intensive human activities in the coal mine area while exploring Cd's harm to human health [42]. The risk value decreases from the northwest and southeast to the center and surroundings. There is a chimney in the power unit, so Cd settles down with the smoke and dust, which leads to a higher cancer risk of Cd around the chimney and in the southwest in a phenomenon consistent with plume dispersion patterns as reflected in the regulatory understanding of atmospheric emissions [43]. Regulatory guidance suggests carrying out harmless treatment of the slag of the total transformation unit and the methanol synthesis unit and controlling the emission of flue gas to reduce or block the potential carcinogenic risk.

Hazard Quotient
The Kriging interpolation method was used to simulate the multipath total hazard quotient of soil Cd in Table 8, and the hazard quotient evaluation map of the entire plant area was obtained ( Figure 4). It can be seen from Table 8 and Figure 4 that the Cd multipath hazard quotients in the soil area of the chemical units were 0.014-0.024 (with an average of 0.020), all of which were within the acceptable level of 1. The main polluted areas were those near the coal transportation corridor, power unit, complete modification unit, and the entry/exit for people. On the whole, the multipath total hazard quotient of Cd in the soil of the main process unit of coal gasification was relatively low. Therefore, its fluctuation under an acceptable level can be ignored. However, Cd is a common environmental contaminant with a long biological half-life in the soil. Moreover, its microbial or chemical loss is small, so it can exist for a long time in the soil. With the long-term development of the factory, Cd accumulates continuously in the soil [44][45][46]. Gorospe [47] analyzed 16 different heavy metals in 91 soil samples from a vegetable garden in San Francisco. The results showed that soils from most (>75%) gardens contained higher levels of heavy metals than those of the California human health screening, and the level of Cd was the highest, reaching 84%. Cd has high plant-soil mobility and enters the human body through the food chain. When it enters the human body, the half-life of Cd can last for 10-30 years, which causes chronic toxicity [48]. Therefore, some attention should be paid to the situation of the entry and exit of people to avoid harm to human health caused by Cd accumulation pollution.

Contribution Ratios of Human Health Risk for Different Exposure Routes
The risk contribution rates of carcinogenic risk and hazard quotient under different exposure routes ( Figure 5) were calculated to provide the basis for targeted prevention and control of human health risk in coal chemical industry areas. It can be seen from Figure 5 that the carcinogenic risk of Cd in the soil of the coal chemical industry areas was mainly caused by oral intake of soil dust, accounting for 80.51% of the total carcinogenic risk, which was 4 and 73 times of the contribution rates of skin exposure and inhalation of soil particles, respectively. However, for the noncarcinogenic risk of Cd in the coal chemical industry area, the contribution rates of the three exposure pathways of Cd in the soil of the coal chemical industry area were respiratory inhalation (79%) > oral intake (17%) > skin contact (4%). Inhalation of soil particles was the main noncarcinogenic exposure route of Cd. Heavy metals in street dust mainly come from vehicle emissions and coal combustion. Some studies [49,50] have showed that coal transportation has an impact on the environment. Dust is rich in heavy metal pollutants and could be easily inhaled or ingested. Also, the regional evaporation is much greater than the precipitation, so the vegetation is sparse and scattered, and the physical weathering and wind effects are obvious. These factors are likely to lead to a higher carcinogenic and noncarcinogenic risk from inhalation and oral intake.

Safety Threshold of Cd in the Coalification Zone Soil
Based on the human health risk assessment of Cd in the coal chemical industry, the risk control values for different exposure routes were calculated (Table 9). When the hazard quotient does not exceed the standard, the control value can be ignored. As is shown in Table 9, the risk control values of Cd under different exposure routes were different within the acceptable carcinogenic risk level of 1E-06. The risk control value of Cd for oral intake was the minimum (0.392 mg/kg), slightly higher than the soil background value, and far less than the industrial land standard. Therefore, the human health risk assessment based on different exposure pathways can be used as a strict standard for the safety  In conclusion, considering the human health risks of Cd in the coal chemical industry area, oral intake and inhalation of soil particles were the main exposure routes, so targeted control and prevention should be carried out during the production process of the coal chemical area. Yang et al. [51] evaluated the health risks of heavy metals in the soil of Changchun and found that Cd was mainly exposed through inhalation and oral intake. This is consistent with our findings. Therefore, in the production process of the coal chemical industry, it is necessary to prevent the exposure of soil particles by oral intake and inhalation. For workers in the plant area, protective measures should be strengthened to reduce exposure of Cd in the soil due to their high exposure. For example, water can be sprayed regularly in the exposed area of soil in the factory districts to suppress dust, and protective masks can be used, to reduce potential risks.

Safety Threshold of Cd in the Coalification Zone Soil
Based on the human health risk assessment of Cd in the coal chemical industry, the risk control values for different exposure routes were calculated (Table 9). When the hazard quotient does not exceed the standard, the control value can be ignored. As is shown in Table 9, the risk control values of Cd under different exposure routes were different within the acceptable carcinogenic risk level of 1E-06. The risk control value of Cd for oral intake was the minimum (0.392 mg/kg), slightly higher than the soil background value, and far less than the industrial land standard. Therefore, the human health risk assessment based on different exposure pathways can be used as a strict standard for the safety threshold of Cd in the soil environment. Different choices of risk levels will also lead to different safety thresholds. For example, the USEPA recommends 1 × 10 −6 -1 × 10 −4 , the United Kingdom generally adopts 1 × 10 −5 in practice, and the Netherlands suggests a looser 1 × 10 −4 [52]. In addition, the background value, geological conditions, and other data should be comprehensively considered to determine the safety threshold of Cd in the soil environment of certain sections of coal chemical areas, which can provide a comprehensive safety guarantee and guidance for human health and regional development of Cd in China's coal chemical industry field [53].

Conclusions
In this study, the distribution of human health risks and the safety threshold of Cd in the soil of the coal chemical plant area were researched. The Cd content in the soil of coal chemical plants was higher than the regional and national background values. Human health risk assessment revealed that the total carcinogenic risk values of Cd exceeded the acceptable level (1 × 10 −6 ). It is necessary to regularly spray water on the soil exposure areas of the plant to suppress the dust in the plant and to wear protective masks to reduce the potential risk.
The spatial distribution of Cd pollution in the plant area was uneven, mainly distributed in the outflow entrance and coal entrance corridor, which is highly dangerous for health. The hazard quotient of Cd was within the acceptable level, but Cd is liable to accumulate in the soil, so it is necessary to pay close attention to avoid its harmful effects on the human body.
The safety threshold of the plant was taken as 0.392 mg/kg, higher than the soil background value. However, compared with the industrial land exceeding standard, there is still a long way. This standard may be too strict for the determination of the safety threshold of the coal chemical industry area, accounting for the difference in carcinogenic risk level and the selection of parameters in different counties and districts as well as the gap between the high toxicity of Cd to human health. The concentration threshold that poses a risk or hazard to human health is low. Therefore, a series of factors such as regional soil background values, geological environments, and different carcinogenic risk levels should be comprehensively considered to determine the safety threshold of heavy metals in the coal chemical industry area.