Next Article in Journal
How the Popularity of Short Videos Promotes Regional Endogeneity in Northwest China: A Qualitative Study
Next Article in Special Issue
Ecological Risk Evaluation and Source Identification of Heavy Metal Pollution in Urban Village Soil Based on XRF Technique
Previous Article in Journal
The Development of ARIMA Models for the Clear Sky Beam and Diffuse Optical Depths for HVAC Systems Design Using RTSM: A Case Study of the Umlazi Township Area, South Africa
Previous Article in Special Issue
Simulating Maize Productivity under Selected Climate Smart Agriculture Practices Using AquaCrop Model in a Sub-humid Environment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Priority Soil Pollution Management of Contaminated Site Based on Human Health Risk Assessment: A Case Study in Southwest China

1
Faculty of Architecture, Civil and Transportation Engineering, Beijing University of Technology, Beijing 100124, China
2
Appraisal Center for Environment and Engineering, Ministry of Ecology and Environment of China, Beijing 100012, China
3
Hebei Institute of Water Science, Shijiazhuang 050051, China
4
Key Laboratory of Environmental Change and Natural Disaster, Ministry of Education, Beijing Normal University, Beijing 100875, China
5
School of National Security and Emergency Management, Beijing Normal University, Beijing 100875, China
6
Faculty of Geographical Science, Academy of Disaster Reduction and Emergency Management, Beijing Normal University, Beijing 100875, China
*
Author to whom correspondence should be addressed.
Sustainability 2022, 14(6), 3663; https://doi.org/10.3390/su14063663
Submission received: 26 January 2022 / Revised: 11 March 2022 / Accepted: 14 March 2022 / Published: 21 March 2022
(This article belongs to the Special Issue Soil Management Practices to Promote Soil Health)

Abstract

:
The human risk assessment model can serve as a tool for regional contaminated site comprehensive management. However, site-specific risk assessment is still seldom applied in China as a basis for making decisions on risk management actions. In this study, a total of 112 soil samples were collected from ten polluted sites in Southwest China. The human health risk assessment method was used to assess the risk of Cd, As, Cu, Pb, Cr, Zn, and Hg pollution. According to the findings, the average concentrations of As, Zn, and Pb in contaminated sites are substantially higher than those of Cu, Cd, Cr, and Hg. Further studies conclusively showed the soil at contaminated site in its present condition pose risk to human health to potential future receptors. The main contribution of non-cancer and cancer risks was through incidental soil and dust ingestion. The priority control site remediation order is Region VI > Region V > Region IV > Region III > Region VII > Region X > Region IX > Region VIII > Region I > Region II. Finally, management recommendations are made, including reducing pollutant intake, implementing a stringent monitoring scheme, utilizing bioremediation, and strengthening the implementation of relevant laws. This study provides a case for the comprehensive evaluation of soil pollution at contaminated sites in China.

1. Introduction

Among all types of soil pollution, heavy metal(loid)s has become an emergent environmental concern worldwide [1]. The southwestern provinces in China (including Guizhou, Yunnan, and Guangxi) are well-known as karst areas for their richness in mineral resources which has resulted in the formation of a metal smelting industrial area [2,3]. Mining and metallurgical operations generate large amounts of wastes and derivatives, in which heavy metal(loid)s are readily dissolved under an oxidized atmosphere, causing contamination of surrounding soil [4]. Since 2017, China’s provincial capital cities have successively announced to the public the status of 174 contaminated sites in their jurisdiction [5]. In provincial capital cities alone, there are more than 144 contaminated sites are being remediated of which 109 are expected to redevelopment recently, and 25% of the contaminated sites have been sold, with a total amount of 104.96 billion RMB [6]. Therefore, the remediation and priority of heavy metal contaminated sites, which requires further study.
Contaminated site management is also a priority for China environmental policies [7,8,9]. The main Chinese regulation for the contaminated site management is represented by the regulation for contaminated site soil environmental [10], which require the environmental investigation, human health risk assessment, risk management and control, remediation, and supervision. Hechi city, located in the Southwest region of China, has been listed as the pilot area for comprehensive prevention and control of soil pollution. Following the selecting of a pilot area for contaminated sites management, the next issue that arises is the prioritizing the contaminate sites to determine which is the most urgent to be remediated [11,12]. Decision-makers must prioritize their management to maximize the benefit derived from limited funds [13]. This is a very complex task that relies on a variety of techniques and treatments, which incorporate multiple technical, social, economic, and environmental criteria required for any informed decision [14]. However, contaminated site prioritization was often founded on risk-based principles that assess the likelihood that a hazard will have an adverse im-pact on a receptor [15]. Health risk assessment (HRA) was based on a distinct assessment of the toxicity of the chemicals by exposure route (i.e., inhalation, ingestion, and dermal contact), as well as the evaluation system of pollution source-contact route-receptor-risk was constructed. Many contaminate site restoration projects may result either second pollution, and prioritization control area can provide remedies for remediation contaminated sites.
The heavy metals emitted from industrial and mining activities tend to initially accumulate in nearby soil [16] and may be transferred to environmental compartments, such as water, soil and plants, and eventually could be absorbed through ingestion, dermal contact and inhalation, increasing their risk to human health [17,18,19]. Human health risk assessment is the most rational approach for contaminated site risk management and has been the prevailing approach by international consensus [20,21,22,23]. By calculating and comparing these human health risk scores, decision-makers are able to determine which site is a priority. There are few published cases of site prioritization based on human health risk assessment in China. Moreover, frameworks that use risk assessment-based methodology to prioritize contaminated site risk management options, together with their application in two case studies in China, have recently been reported in the literature [24]. There is sufficient evidence to indicate the crucial role of human health risk assessment in the process of contaminated site redevelopment and subsequent risk management decision-making.
Therefore, it is essential to illustrate the crucial role of human health risk assessment (HHRA) in contaminated site prioritization and subsequent risk management decision-making. However, few attempts have been made to prioritize sites in industrial regions, particularly in regions with large numbers of abandoned industrial contaminated sites. We address this gap by developing a risk-based tool to prioritize contaminated sites in Southwest China. Herein, we report the following: (1) the soil contamination characteristics of contaminated sites in Southwest China; (2) the human health risk derived from contaminated sites; and (3) the priority control areas determined for contaminated sites. The findings of this study will be useful in providing decision-makers with a pragmatic solution for the determination of priority and, ultimately, regional contaminated site comprehensive management.

2. Materials and Methods

2.1. Study Area

The study area is located in the middle of Hechi city, Guangxi Province (106°34′–109°09′ E, 23°41′–25°37′ N), China (Figure 1). The study area, located in the met-allogenic belt of nonferrous metals, is a typical metalliferous industrial district with abundant mineral resources [25]. Due to lack of proper disposal and safe landfill, the waste residues caused serious pollution to the surrounding environment of the study area. Limestone soil is the major soil type in the karst area [26]. According to the China National Environmental Monitoring Center, the background value of heavy metals (loid) in limestone soil is higher than that in other soils types [27].
The study area has monsoon-influenced humid subtropical climate. The temperature difference between season is small, with the average temperature is about 16.9–21.5 °C. The frost-free period is long, the annual sunshine hours are 1447–1600 h in most areas [25]. On the whole, the terrain is high in the northwest and low in the southeast. According to monitoring data, as a result of long-term mining area, the quality of soil environment in the region declined, especially the contamination of arsenic, cadmium and lead that are in excess.

2.2. Sample Collection and Analysis

Using a regular grid, soil samples were collected from the 10 sites in Hechi city, Guangxi Province, from November 2016 to December 2016 (Figure 1). Before field work, the sampling coding system and sampling cell were created using the global positioning system on a digital map. A total of 112 samples were collected from the soil surface of the study area. An S-shaped sampling approach was adopted to sample flat locations with high soil heterogeneity. A quincunx approach was adopted to sample flat locations with homogeneous soil. Considering the identified primary pollution sources, the following contaminants were chosen for risk assessment purposes: cadmium (Cd), arsenic (As), copper (Cu), lead (Pb), chromium (Cr), zinc (Zn) and mercury (Hg). According to site characteristics and profiles length, the focus was on the surface soil layer up to the depth of 0.5 m, using a regular grid of 200 m × 200 m. Before sampling, the surface soil was levelled with a shovel to remove surface impurities. The samples were collected from the undisturbed soil and sealed in polyethylene bags with a wooden shovel. All samples were stored in a low-temperature incubator at 4 °C until all soil samples were sent to the laboratory.
During the pretreatment of soil samples, all samples were dried in a cool and ventilated place, and after natural air drying, animal residues, stones and plant roots were removed. Then, the soil samples were passed through a 100-mesh sieve and digested with HNO3-HF-HClO4 [1,18]. The concentrations of heavy metals (Cd, As, Cu, Pb, Cr, and Zn) were determined by inductively coupled plasma–mass spectrometry (ICP–MS, Bruker 820-MS) [28,29]. Atomic fluorescence spectrometry (AFS-9700) was used to determine the content of total mercury in the soil. Samples were analyzed using national standard soil samples (GSS-1, GSS-2, GSS-3, and GSS-4) as references, while quality assurance and quality control procedures were performed in accordance with USEPA 6020A [30]. The results indicated that the accuracy and bias of the analysis were generally less than 10%.

2.3. Health Risk Assessment Model

According to the health risk assessment model recommended by the U.S. Environmental Protection Agency (US EPA), the carcinogenic and noncarcinogenic risks of soil heavy metal exposure in the study area were evaluated and predicted. The health risks of heavy metals depend on two aspects. One aspect is the level of environmental pollution, including the concentration, form and toxic effects of heavy metals [18]. The other aspect is human exposure behavior, including the behavior and characteristics of human exposure to heavy metals. After heavy metals enter the soil, the exposure of the population is estimated in three ways: skin contact, hand–oral ingestion, and inhalation of dust with attached heavy metals [31,32]. These exposure routes are represented as follows:
ADD ingest = c × IngR × EF × ED × EF BW × AT ,
ADD inh = c × SA × AF × ABS × EF × ED × CF BW × AT ,
ADD der = c × EF × ED PEF × AT ,
As adults and children have different physical characteristics and behavior patterns, Table 1 provides detailed values and descriptions of parameters in health risk assessment. The human body exposure parameters were mainly selected from CRNCDSR (2015), while the toxicity parameters for the heavy metals were derived from USEPA guidelines and international research findings. Various types of heavy metals enter the human body through different exposure routes that produce different toxic effects. The hazard index (HI) was defined as the relationship between the predicted exposure and the theoretical reference (RfD) for exposure. It is generally believed that when HI < 1.0, pollutants will not pose a threat to human health, so the risk can be ignored. However, when HI > 1.0, the pollutants will have adverse effects on human health, indicating that there is a risk of cancer. Cancer risk (CR) is used to assess the risk of cancer. CR was measured by multiplying the exposure by the carcinogenic slope factor (SF). The SF consists of oral intake, skin exposure and inhalation risk. HI and CR can be calculated using the following formula [18,31]:
HI = ADD RfD ,
CR = ADD × SF ,
Cancer risk is calculated only for some established heavy metals. When CR < 10−6, the cancer risk of a single pollutant is acceptable. When 10−6 < CR < 10−4, the combined cancer risk of contaminants is acceptable. When CR > 10−4, contaminated site cleanup is recommended [33,34].
Table 1. Calculation parameters used in human health risk assessment [18,27,35].
Table 1. Calculation parameters used in human health risk assessment [18,27,35].
ParameterSymbolUnitsValue
Body weightBWkg70 (adults); 15 (children)
Average timeATdED × 365 (noncarcinogenic), 25,550 (carcinogenic)
Exposure frequencyEFd/y350
Exposure durationEDy24 (adults), 6 (children)
Soil ingestion rateIngRmg/d100 (adults), 200 (children)
Conversion factorCFkg/mg1 × 10−6
Surface areaSAcm2/d5700 (adults), 2800 (children)
Adherence factor to skinAFsoilmg cm20.5 (trace metals), 0.2 (organic compounds)
Dermal absorption factorABSUnitless0.03 (As), 0.001 (trace metals except As)
Gastrointestinal absorption factorABSGIUnitless0.025 (Cd),1.0 (As), 1.0 (Cu), 1.0 (Pb), 0.025 (Cr), 1.0 (Zn), 0.07 (Hg)
Particle emission factorPEFm3/kg1.36 × 109
Ingestion reference valueOral RfDmg/kg d3.00 × 10−4 (As), 1.00 × 10−3 (Cd), 3.00 × 10−3 (Cr), 4.00 × 10−2 (Cu), 3.00 × 10−4 (Hg), 3.5 × 10−3 (Pb), 3.00 × 10−1 (Zn)
Dermal reference valueDermal RfDmg/kg d3.00 × 10−4 (As), 2.50 × 10−5 (Cd), 7.50 × 10−5 (Cr), 4.00 × 10−2 (Cu),
2.10 × 10−5 (Hg), 3.50 × 10−3 (Pb), 3.00 × 10−1 (Zn)
Inhalation reference valueInhal RfDmg/m31.50 × 10−5 (As), 1.00 × 10−5 (Cd), 1.00 × 10−4 (Cr), 4.02 × 10−2 (Cu), 3.00 × 10−4 (Hg), 3.25 × 10−3 (Pb), 3.00 × 10−1 (Zn)
Oral slope factorSFO(mg kg−1 d−1)−11.5 (As), ** (Cd), 5.0 (Cr), ** (Cu), ** (Hg), 8.50 × 10−3 (Pb), ** (Zn)
Dermal contact factorSFO×ABSGI(mg kg−1 d−1)−11.5 (As), ** (Cd), 2.0 (Cr), ** (Cu), ** (Hg), 8.50 × 10−3 (Pb), ** (Zn)
Inhalation unit riskIUR(mg/m−3)−14.3 × 10−3 (As), 1.80 × 10−3 (Cd), 8.40 × 10−2 (Cr), ** (Cu), ** (Hg),
1.2 × 10−5 (Pb), ** (Zn)
**: not available.

2.4. Development of Site-Specific Risk-Based Remedial Levels (RBRLs)

The second critical step in HHRA is the development of risk remediation levels at specific sites. When the total HI exceeds 1 or the total CR exceeds 1 × 10−4, the remedial level is calculated by the simplified method of the US EPA (2000) [36], and the formula is as follows:
RBRL = C × TR Calculated   Risk ,
where C is the concentration of heavy metals in the soil, TR represents the noncancer targeted risk (NTR) or cancer targeted risk (CTR) [37], the calculated risk is the value of the noncancer risk and cancer risk in the HHRA, and RBRL is the risk-based remediation level.
Based on the remediation recommendations of the US EPA (2002) [38], the site-specific RBRLs were calculated separately for carcinogenic and noncarcinogenic effects and separately for oral/dermal and inhalation exposures due to the potential for different health effects (target organs) via these routes. For noncarcinogens, an HI of 1 was established as the NTR since an HI above 1 indicates that toxic effects may appear [31]. If more than one noncarcinogen affects the same target organ system, the RBRLs should be adjusted to reflect the potential for additive risks [38]. The adjustment can be made by dividing the RBRL calculated for individual noncarcinogens by the number of chemicals with the same target organs/effects [37]. Concerning carcinogens, since multiple contaminants and multiple pathways of exposure occurred at contaminated sites, the CTR of 1 × 10−6 was initially considered sufficiently protective. Finally, the RBRLs could be modified so that the cumulative risk level would be at least at the 1 × 10−4 level corresponding to the upper end of the US EPA’s acceptable risk range of 1 × 10−6 to 1 × 10−4 [39].
The calculated soil remediation target values were compared, and the lowest of these values was suggested as a preliminary remedial level. The final repair target value needs to be compared with the local background value; if the lowest RBRL values were lower than the relevant background concentration, the background concentration could be considered the remedial level. However, this should be a site-specific decision and should not be carried out automatically.

2.5. Special Case: Lead

Inorganic lead is a special case with strong toxicity. There is no safety level of minimum allowable lead exposure, so HHRA is no longer used in the toxicity evaluation of lead pollution. For industrial and commercial land, the ALM model is usually used to assess lead exposure and formulate a RBRL [39,40,41]. Lead exposure pathways in the direct ingestion of soil and indoor dust are considered in the adult ALM health model. The blood lead concentration in women of childbearing age was calculated by using Formula (7):
PbB a d u l t , c e n t r a l = PbB a d u l t , 0 + PbS × BKSF × I R s × A F s × E F s A T s ,
where PbB a d u l t , c e n t r a l is the blood lead concentration of adults at the center of the site exposed to lead pollution. PbB a d u l t , 0 is the background lead level in blood from women of childbearing age not exposed to lead, with a reference value of 4.79 μg/dL. PbS: refers to the lead content in soil (mg/kg). BKSF refers to the slope coefficient of blood lead vs. daily intake of lead content in the body, and the reference value is 0.4 d/dl. IRs refers to the daily soil or indoor dust intake, with a reference value of 0.05 d/dL. AFs refers to the gastrointestinal absorption efficiency of lead intake in the body, with a reference value of 0.12 (dimensionless). EFs refers to the average number of days to lead pollution exposure every year, with a reference value of 240 d/a. Ats refers to the average time of long-term exposure, with a reference value of 365 days/y.
The risk-based criterion ensures that the probability of fetal PbB a d u l t , c e n t r a l being 10 μg/dL does not exceed 5%, and the risk-based target of adult blood lead concentration is calculated by Formula (8):
PbB a d u l t , c e n t r a l ,   g o a l = PbB f e t a l , 095 , g o a l G S D i ,   a d u l t 1.645 × R f e t a l / m a t e r n a l ,
where PbB a d u l t , c e n t r a l ,   g o a l is the target value of the average blood lead content of pregnant women exposed to lead-contaminated sites. PbB f e t a l , 0.95 , g o a l is the 95% probability target value of fetal blood lead content, with a reference value of 10 μg/dL. G S D i ,   a d u l t 1.645 is the geometric standard deviation of blood lead content of women of childbearing age, with a reference value of 1.48 (dimensionless). R f e t a l / m a t e r n a l is the correlation coefficient of blood lead content between a fetus and mother, with a reference value of 0.9 (dimensionless).
The soil lead concentration for a given exposure scenario and PbB a d u l t , c e n t r a l ,   g o a l can be calculated by Formula (9), and PbB a d u l t , c e n t r a l , g o a l can be replaced by PbB a d u l t , c e n t r a l .
RBPL = PbS = PbB a d u l t , c e n t r a l , g o a l PbB a d u l t , 0 × A T s BKSF × I R s × A F s × E F s ,
where RBPL represents the desired remediation value of soil lead concentration (PBs) for adult blood lead concentration ( PbB a d u l t , c e n t r a l ,   g o a l ) and the corresponding 95% foetal blood lead concentration ( PbB f e t a l , 095 , g o a l ) at a specific site.

3. Results and Discussion

3.1. Descriptive Statistics of Soil Heavy Metals

The descriptive statistics of heavy metals in soil at different sites are shown in Table 2. The elemental contents determined for each polluted site were compared with their mean background levels in Guangxi Province and China (Table A1). The mean background levels of all elements were higher in Guangxi Province than in China [26,42,43]. There were differences in the degrees of heavy metal pollution at the different polluted sites. The contents of Cu, Pb, Cr, Zn and Hg in region I are the lowest, and the contents of Cd and As in region II are the lowest. The contents of Cu and Pb in region V are the highest. The contents of As and Hg in region VI are the highest. The contents of Cu and Pb in area V are the highest. However, the contents of As and Hg in region VI are the highest. The content of Cr in region IX is the highest (Table 2). The contents of Cd and Pb in region X are the highest.
The coefficient of variation (CV) can reflect the data dispersion degree. According to relevant research on the degree of variation, the CV can be divided into 0 < CV < 0.16 (weak variation), 0.16 ≤ CV< 0.36 (medium variation), and CV ≥ 0.36 (high variation) [44,45]. Overall, from the CVs of heavy metals in regions I~X, Pb and Hg in region II show weak variation, Pb and Zn in region VIII show medium variation, and heavy metals in other regions show high variation (Table A1), indicating that the content of heavy metals in each region changes significantly in space and is subject to discrete inputs of certain natural or external factors.

3.2. Human Health Risk Assessment

3.2.1. Noncarcinogenic Risk

The noncarcinogenic risks of heavy metals by various exposure routes are shown in Table A2. The total noncarcinogenic risk was higher in children than in adults at 10 contaminated sites. The reason for this may be related to children’s lifestyle and physical characteristics (playing places and unhealthy eating habits), leading to children having more opportunities to contact soil particles polluted with heavy metals than adults [46]. In terms of individual elements, the THI values of adults for As in regions III, IV, V, VI, and VII exceeded the safety threshold (1.0). Except for regions II and X, the THI values of children for As exceeded the safety threshold (1.0). The potential mixture effect values of adults in regions III, IV, V, VI, and VII exceeded the safety threshold (1.0) and were 7.66, 11.09, 81.36, 210.79, and 1.86, respectively (Figure 2). Except for region II, the potential mixture effects value of children exceeded the safety threshold (1.0) (Figure 2). This is due to the high concentration of As in the soil, high noncarcinogenic risk to humans, and the low standard value of RfD.
Heavy metals in soil can enter the human body via three pathways: oral soil intake, skin contact and respiratory inhalation. Table A2 shows that for both adults and children, oral soil intake is the main exposure route in the study area, followed by dermal soil intake and soil intake by inhalation. Therefore, local residents and the government should consider protecting the health of people, limit people entering the vicinity of contaminated sites and encourage relevant protection research.

3.2.2. Cancer Risk

The total CRs basically exceeded the cancer risk value of 1 × 10−4, except for regions I, II and VIII (Figure 2). The total cancer risk resulted mainly from arsenic and comprised mostly oral and dermal risks, while the inhalation cancer risks were insignificantly similar to noncarcinogenic risks. Except for regions I, II and VII, the minimum total cancer risk value of As is one order of magnitude higher than the target value. Other research has found that arsenic is the primary contributor to total cancer and noncancer hazards in contaminated locations [47,48]. A previous study showed that health risks are mainly caused by direct oral intake of contaminated soils rather than dermal contact and inhalation. The inhalation cancer risks were insignificant and could be overlooked similarly to noncarcinogenic risks.

3.3. Remedial Analysis of Site-Specific Risks

When the health risk assessment value exceeds the baseline for the accepted risk level (noncancer target risk value and cancer target risk value), the remedial level repair target value is calculated (Table 3). The RBRLslowest for noncarcinogenic effects is equal to the NTR (HI) of 1. In terms of cancer risk, only heavy metals (arsenic) were the major carcinogenic risk substances in different study areas. Therefore, the TCR risk value of 1 × 10−6 was adjusted to 1 × 10−5 [35,49].
The results of Table 3 show that the RBRLslowest for children were significantly lower than that for adults in all regions. Since children are more susceptible to the health impacts of heavy metals, even at extremely low levels, they must be protected. For adults, the RBRLslowest values of Cu, Zn and Hg were based on noncarcinogenic risks and were higher than those of all analysed samples. The RBRLslowest values of Cd, As and Cr were based on carcinogenic risks and were lower than the maximum values. For children, the RBRLslowest values of Cd, As, Cu and Hg were lower than the maximum value for all analysed samples based on noncarcinogenic risk, while that of Zn was higher than the maximum value. The RBRLslowest of Cr was lower than that of all analysed samples based on carcinogenic risk. As shown in Table 3, the RBRLslowest is determined mainly by exposure by oral and respiratory routes.
This study took Pb as an example. The RBRLslowest values were assessed individually as a special case with a PbB of 10 μg/dL as the goal and a maximum value that was lower than the applicable maximum value. The industrial Pb RBRLs amounted to 331 mg/kg and were higher than the region I~X background levels for lead (Table 3).

3.4. Priority Control Regions and Recommendations

The health risk assessment results of different sites are shown in Table A2, and the urgency of the remediation of contaminated sites is analyzed according to the risk values. The results show that the total risk value of each region is ranked as follows: Region VI > Region V > Region IV > Region III > Region VII > Region X > Region IX > Region VIII > Region I > Region II.
Arsenic contributes the most to the total risk of each region and is the main metal causing the health risk of all regions, followed by the heavy metals Pb, Cr and Cd, and the contributions of the heavy metals Cu, Zn and Hg are almost negligible. Further analysis found that the order of total risk for As (adult or child, noncarcinogenic or cancer) was basically consistent with the order of total risk values for each region. Therefore, the urgency of remediation at each site is mainly related to the heavy metal arsenic. Reasonable management techniques should concentrate on how to regulate pollution cleanup as effectively as possible, some of which are detailed below.
(1) Given the enormous impact of mining operations in priority control regions, comprehensive environmental protection legislation enforcement is essential, particularly with regard to large emission sources such as mining smelting and other metal consuming sectors. Although all mining firms in China are obligated by law to engage in environmental restoration efforts, mining enterprises continue to cause substantial harm to the environment via pollution owing to the use of antiquated technology and inefficiency of mining operations.
(2) It is critical to establish strict monitoring in these areas, since this will allow for a realistic assessment of the existing contamination state and forecasting of future trends. Note that to acquire precise contamination characteristics of soil arsenic, more frequent routine sampling is needed. Only by obtaining real and exact concentrations and geographical dispersion findings will policy-makers be able to provide better guidelines for precision control.
(3) Physical, chemical, and biological cleanup of arsenic-contaminated soils are the most often used methods. For the remediation of polluted soil, bioremediation has been seen as a low-cost, ecologically benign method [50]. Plant-microbial joint remediation technology has been used to carry out research. The sizes and costs determine which methods may be employed in polluted site cleanup.

4. Conclusions

This study provides valuable information for health risk assessment and risk remediation at different contaminated sites. The results of health risk assessment shown that arsenic and chromium have higher non-carcinogenic and carcinogenic risks than other trace elements. The risk assessment results and the developed RBRLs could also be used to identify the priority area for remediation. Management recommendations such as minimizing pollutant input, creating tight monitoring systems, and implementing biological remediation were made in order to safeguard Chinese inhabitants from heavy metal contamination. Implementing these techniques necessitates not only technology advancement, but also socioeconomic evaluation and efficient implementation of environmental protection regulations by local governments.
The risk assessment employed in this study includes several potential uncertainties, such as model selection and localized parameter values. It should also be noted that when planning the remediation strategy of the site, apart from human health risk, several other aspects must be taken into account, such as economic cost, implementation, and environmental impact. Further studies are required to focus on areas with a high health risk from As combined with local specific conditions and application of more accurate spatial distribution expression approaches to make remediation actions conducive to implementation.

Author Contributions

Conceptualization, formal analysis, supervision, J.W.; methodology, data curation, software, R.J.; formal analysis, supervision, editing, H.X. and Y.X.; methodology, formal analysis, funding acquisition, D.Z. and G.Z.; All authors have read and agreed to the published version of the manuscript.

Funding

This study was financially supported by the investigation of arsenic contaminated sites in Southwest China (No. 201501).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data available on request due to privacy and ethical restrictions.

Acknowledgments

The authors thank the editors and anonymous experts for their comments, which have greatly improved the study.

Conflicts of Interest

The author states that there are no conflict of interest to declare.

Appendix A

Table A1. A summary of soil heavy metals in different regions.
Table A1. A summary of soil heavy metals in different regions.
RegionConcentrationCdAsCuPbCrZnHg
ICV1.071.231.143.130.881.270.746
Over Cg13.20.6170.1310.3360.00050.3480.0008
Over Cc2.110.5670.1200.2430.00050.4920.0010
IICV0.5270.9350.92900.5380.2770
Over Cg0.4090.0530.0150.0030.00030.0110.029
Over Cc0.0660.0480.0140.0020.00030.0150.036
IIICV1.111.81-0.641-0.9521.12
Over Cg17986.9-1.62-4.1118.8
Over Cc28.779.9-1.17-5.8123.4
IVCV2.121.11-0.452-0.46816.3
Over Cg840125-6.06-21.219.2
Over Cc134115-4.38-30.123.9
VCV1.271.291.181.66-1.20-
Over Cg92292511.622.6-34.8-
Over Cc14885010.716.4-49.2-
VICV0.5252.14-1.25-0.5141.03
Over Cg7572402-5.33-8.4351.4
Over Cc1212209-3.86-11.963.9
VIICV0.8391.17-0.603-0.6331.52
Over Cg34.929.5-1.06-1.3629.4
Over Cc5.5827.1-0.764-1.9236.5
VIIICV0.5960.825-0.060-0.3440.373
Over Cg3.171.46-2.45-1.152.71
Over Cc0.5071.34-1.77-1.633.36
IXCV1.590.5021.532.210.3850.7880.710
Over Cg753.321.7411.51.583.544.25
Over Cc123.051.5998.321.525.015.29
XCV2.221.351.163.320.6061.321.97
Over Cg3567.264.0624.00.8073.786.62
Over Cc56.96.673.7317.30.7785.358.22
Note: CV/%, coefficient of variation; Cg, background value of metals in soils of Guangxi Province. Cc, background value of metals in soils of China.
Table A2. Summary of the mean values of noncancer and cancer risk.
Table A2. Summary of the mean values of noncancer and cancer risk.
Region IAdultsChildren
Noncancer Risk (Unitless)Cancer Risk (Unitless)Noncancer Risk (Unitless)Cancer Risk (Unitless)
ElementHIoralHIdermalHIinhalationTHICRoralCRdermalCRinhalationTCRHIoralHIdermalHIinhalationTHICRoralCRdermalCRinhalationTCR
Cd3.29 × 10−42.88 × 10−41.48 × 10−56.32 × 10−4NANA9.16 × 10−149.16 × 10−147.53 × 10−32.69 × 10−31.48 × 10−51.02 × 10−2NANA2.29 × 10−142.29 × 10−14
As2.90 × 10−22.48 × 10−22.99 × 10−45.41 × 10−24.48 × 10−63.83 × 10−66.61 × 10−128.31 × 10−69.48 × 10−12.71 × 10−12.99 × 10−41.22 × 101.04 × 10−57.31 × 10−81.65 × 10−121.05 × 10−5
Cu9.28 × 10−52.65 × 10−64.75 × 10−89.55 × 10−5NANANANA8.66× 10−46.07 × 10−54.75 × 10−89.27 × 10−4NANANANA
Pb2.47 × 10−37.05 × 10−51.37 × 10−62.55 × 10−32.52 × 10−8NA1.83 × 10−142.52 × 10−82.31 × 10−21.62 × 10−31.37 × 10−62.47 × 10−25.89 × 10−84.12 × 10−104.58 × 10−155.93 × 10−8
Cr1.45 × 10−51.27 × 10−51.96 × 10−72.74 × 10−56.54 × 10−9NA5.66 × 10−136.54 × 10−93.32 × 10−41.19 × 10−41.96 × 10−74.51 × 10−41.53 × 10−84.27 × 10−91.41 × 10−131.95 × 10−8
Zn7.37 × 10−52.10 × 10−63.79 × 10−87.58 × 10−5NANANANA6.88 × 10−44.82 × 10−53.79 × 10−87.36 × 10−4NANANANA
Hg3.15 × 10−71.28 × 10−71.62 × 10−104.44 × 10−7NANANANA2.94 × 10−62.94 × 10−61.62 × 10−105.89 × 10−6NANANANA
Total3.20 × 10−22.52 × 10−23.15 × 10−45.75 × 10−24.51 × 10−63.83 × 10−67.28 × 10−128.34 × 10−69.80 × 10−12.75 × 10−13.15 × 10−41.26 × 101.05 × 10−57.78 × 10−81.82 × 10−121.06 × 10−5
Region II
Cd1.02 × 10−58.97 × 10−64.62 × 10−71.97 × 10−5NANA2.85 × 10−152.85 × 10−152.34 × 10−48.37 × 10−54.62 × 10−73.19 × 10−4NANA7.13 × 10−167.13 × 10−16
As2.48 × 10−32.12 × 10−32.55 × 10−54.62 × 10−33.82 × 10−73.27 × 10−75.64 × 10−137.09 × 10−78.09 × 10−22.31 × 10−22.55 × 10−51.04 × 10−18.92 × 10−76.24 × 10−91.41 × 10−138.98 × 10−7
Cu1.08 × 10−53.08 × 10−75.54 × 10−91.11 × 10−5NANANANA1.01 × 10−47.06 × 10−65.54 × 10−91.08 × 10−4NANANANA
Pb1.92 × 10−55.47 × 10−71.06 × 10−81.97 × 10−51.96 × 10−10NA1.42 × 10−161.96 × 10−101.79 × 10−41.25 × 10−51.06 × 10−81.92 × 10−44.56 × 10−103.20 × 10−123.55 × 10−174.60 × 10−10
Cr1.02 × 10−58.90 × 10−61.37 × 10−71.92 × 10−54.58 × 10−9NA3.96 × 10−134.58 × 10−92.33 × 10−48.31 × 10−51.37 × 10−73.16 × 10−41.07 × 10−82.99 × 10−99.90 × 10−141.37 × 10−8
Zn2.28 × 10−66.50 × 10−81.17 × 10−92.35 × 10−6NANANANA2.13 × 10−51.49 × 10−61.17 × 10−92.28 × 10−5NANANANA
Hg1.14 × 10−54.65 × 10−65.88 × 10−91.61 × 10−5NANANANA1.07 × 10−41.07 × 10−45.88 × 10−92.13 × 10−4NANANANA
Total2.54 × 10−32.14 × 10−32.61 × 10−54.71 × 10−33.87 × 10−73.27 × 10−79.63 × 10−137.14 × 10−78.18 × 10−22.34 × 10−22.61 × 10−51.05 × 10−19.03 × 10−79.24 × 10−92.41 × 10−139.12 × 10−7
Region III
Cd4.48 × 10−33.93 × 10−32.02 × 10−48.60 × 10−3NANA1.25 × 10−121.25 × 10−121.03 × 10−13.66 × 10−22.02 × 10−41.39 × 10−1NANA3.12 × 10−133.12 × 10−13
As4.09E × 103.50E × 104.21 × 10−27.63E × 106.31 × 10−45.40 × 10−49.31 × 10−101.17 × 10−31.34E × 1023.82E × 1014.21 × 10−21.72E × 1021.47 × 10−31.03 × 10−52.33 × 10−101.48 × 10−3
CuNANANANANANANANANANANANANANANANA
Pb1.19 × 10−23.39 × 10−46.60 × 10−61.22 × 10−21.21 × 10−7NA8.82 × 10−141.21 × 10−71.11E× 10−17.78 × 10−36.60 × 10−61.19 × 10−12.83 × 10−71.98 × 10−92.21 × 10−142.85 × 10−7
CrNANANANANANANANANANANANANANANANA
Zn8.70 × 10−42.48 × 10−54.48 × 10−78.95 × 10−4NANANANA8.12 × 10−35.68 × 10−44.48 × 10−78.69 × 10−3NANANANA
Hg7.49 × 10−33.05 × 10−33.85 × 10−61.05 × 10−2NANANANA6.99 × 10−26.99 × 10−23.85 × 10−61.40 × 10−1NANANANA
Total4.12E × 103.51× 104.23 × 10−27.66 × 106.31 × 10−45.40 × 10−49.33 × 10−101.17 × 10−31.34 × 1023.83 × 1014.23 × 10−21.72 × 1021.47 × 10−31.03 × 10−52.33 × 10−101.48 × 10−3
Region VI
Cd2.10 × 10−21.84 × 10−29.47 × 10−44.03 × 10−2NANA5.85 × 10−125.85 × 10−124.81 × 10−11.72 × 10−19.47 × 10−46.54 × 10−1NANA1.46 × 10−121.46 × 10−12
As5.89 × 105.03 × 106.06 × 10−21.10E × 1019.09 × 10−47.77 × 10−41.34 × 10−91.69 × 10−31.92 × 1025.50 × 1016.06 × 10−22.47 × 1022.12 × 10−31.48 × 10−53.35 × 10−102.13 × 10−3
CuNANANANANANANANANANANANANANANANA
Pb4.46 × 10−21.27 × 10−32.47 × 10−54.59 × 10−24.55 × 10−7NA3.31 × 10−134.55 × 10−74.16 × 10−12.91 × 10−22.47 × 10−54.46 × 10−11.06 × 10−67.43 × 10−98.27 × 10−141.07 × 10−6
CrNANANANANANANANANANANANANANANANA
Zn4.50 × 10−31.28 × 10−42.31 × 10−64.63 × 10−3NANANANA4.20 × 10−22.94 × 10−32.31 × 10−64.49 × 10−2NANANANA
Hg7.62 × 10−33.10 × 10−33.92 × 10−61.07 × 10−2NANANANA7.12 × 10−27.12 × 10−23.92 × 10−61.42 × 10−1NANANANA
Total5.97 × 105.06 × 106.16 × 10−21.11 × 1019.09 × 10−47.77 × 10−41.35 × 10−91.69 × 10−31.93E × 1025.52 × 1016.16 × 10−22.49 × 1022.12 × 10−31.48 × 10−53.37 × 10−102.14 × 10−3
Region V
Cd2.30 × 10−22.02 × 10−21.04 × 10−34.43 × 10−2NANA6.42 × 10−126.42 × 10−125.28 × 10−11.89 × 10−11.04 × 10−37.18 × 10−1NANA1.61 × 10−121.61 × 10−12
As4.35 × 1013.72 × 1014.48 × 10−18.11 × 1016.71 × 10−35.74 × 10−39.90 × 10−91.24 × 10−21.42 × 1034.06 × 1024.48 × 10−11.83 × 1031.57 × 10−21.10 × 10−42.48 × 10−91.58 × 10−2
Cu8.26 × 10−32.35 × 10−44.23 × 10−68.50 × 10−3NANANANA7.71 × 10−25.40 × 10−34.23 × 10−68.25 × 10−2NANANANA
Pb1.67 × 10−14.75 × 10−39.23 × 10−51.71 × 10−11.70 × 10−6NA1.23 × 10−121.70 × 10−61.55 × 101.09 × 10−19.23 × 10−51.66 × 103.96 × 10−62.77 × 10−83.09 × 10−133.99 × 10−6
CrNANANANANANANANANANANANANANANANA
Zn7.38 × 10−32.10 × 10−43.80 × 10−67.59 × 10−3NANANANA6.88 × 10−24.82 × 10−33.80 × 10−67.37 × 10−2NANANANA
HgNANANANANANANANANANANANANANANANA
Total4.37 × 1013.72 × 1014.49 × 10−18.14 × 1016.71 × 10−35.74 × 10−39.91 × 10−91.24 × 10−21.42 × 1034.06 × 1024.49 × 10−11.83 × 1031.57 × 10−61.10 × 10−42.48 × 10−91.58 × 10−2
Region VI
Cd1.89 × 10−21.66 × 10−28.54 × 10−43.64 × 10−2NANA5.27 × 10−125.27 × 10−124.34 × 10−11.55 × 10−18.54 × 10−45.89 × 10−1NANA1.32 × 10−121.32 × 10−12
As1.13 × 1029.66 × 1011.16E × 102.11 × 1021.74 × 10−21.49 × 10−22.57 × 10−83.23 × 10−23.69 × 1031.05 × 1031.16 × 104.74 × 1034.07 × 10−22.85 × 10−46.43 × 10−94.09 × 10−2
CuNANANANANANANANANANANANANANANANA
Pb3.92 × 10−21.12 × 10−32.17 × 10−54.04 × 10−24.00 × 10−7NA2.91 × 10−134.00 × 10−73.66 × 10−12.56 × 10−22.17 × 10−53.92 × 10−19.34 × 10−76.54 × 10−97.27 × 10−149.40 × 10−7
CrNANANANANANANANANANANANANANANANA
Zn1.79 × 10−35.09 × 10−59.19 × 10−71.84 × 10−3NANANANA1.67 × 10−21.17 × 10−39.19 × 10−71.78 × 10−2NANANANA
Hg2.04 × 10−28.32 × 10−31.05 × 10−52.88 × 10−2NANANANA1.91 × 10−11.91 × 10−11.05 × 10−53.82 × 10−1NANANANA
Total1.13 × 1029.66 × 1011.16 × 102.11× 1021.74 × 10−21.49 × 10−22.57 × 10−83.23 × 10−23.69 × 1031.05 × 1031.16 × 104.75 × 1034.07 × 10−22.85 × 10−46.43 × 10−94.09 × 10−2
Region VII
Cd7.75 × 10−47.64 × 10−43.93 × 10−51.58 × 10−3NANA2.43 × 10−132.43 × 10−131.78 × 10−27.13 × 10−33.93 × 10−52.49 × 10−2NANA6.07 × 10−146.07 × 10−14
As1.39 × 104.40 × 10−11.43 × 10−21.84E × 102.14 × 10−46.79 × 10−53.15 × 10−102.82 × 10−41.68 × 1011.29 × 1011.43 × 10−22.98 × 1014.99 × 10−41.30 × 10−67.88 × 10−115.00 × 10−4
CuNANANANANANANANANANANANANANANANA
Pb3.63 × 10−38.56 × 10−52.01 × 10−63.71 × 10−33.70 × 10−8NA2.69 × 10−143.70 × 10−83.38 × 10−21.96 × 10−32.01 × 10−63.58 × 10−28.63 × 10−85.01 × 10−106.72 × 10−158.68 × 10−8
CrNANANANANANANANANANANANANANANANA
Zn2.87 × 10−48.03 × 10−61.48 × 10−72.95 × 10−4NANANANA2.68 × 10−31.84 × 10−41.48 × 10−72.87 × 10−3NANANANA
Hg1.17 × 10−23.88 × 10−36.01 × 10−61.55 × 10−2NANANANA1.09 × 10−18.88 × 10−26.01 × 10−61.98 × 10−1NANANANA
Total1.40 × 104.45 × 10−11.43 × 10−21.86 × 102.14 × 10−46.79 × 10−53.16 × 10−102.82 × 10−41.70 × 1011.30 × 1011.43 × 10−23.00 × 1014.99 × 10−41.30 × 10−67.89 × 10−115.00 × 10−4
Region VIII
Cd3.78 × 10−46.94 × 10−53.57 × 10−64.51 × 10−4NANA2.20 × 10−142.20 × 10−146.48 × 10−48.67 × 10−33.57 × 10−69.32 × 10−3NANA5.51 × 10−155.51 × 10−15
As2.89 × 10−16.87 × 10−27.07 × 10−43.58 × 10−11.06 × 10−54.46 × 10−51.56 × 10−115.52 × 10−56.41 × 10−11.10 × 1017.07 × 10−41.17 × 1012.47 × 10−58.52 × 10−73.91 × 10−122.56 × 10−5
CuNANANANANANANANANANANANANANANANA
Pb1.80 × 10−22.89 × 10−49.98 × 10−61.83 × 10−21.84 × 10−7NA1.33 × 10−131.84 × 10−71.68 × 10−16.62 × 10−39.98 × 10−61.75 × 10−14.28 × 10−71.69 × 10−93.34 × 10−144.30 × 10−7
CrNANANANANANANANANANANANANANANANA
Zn2.44 × 10−46.69 × 10−61.25 × 10−72.50 × 10−4NANANANA2.27 × 10−31.53 × 10−41.25 × 10−72.43 × 10−3NANANANA
Hg1.07 × 10−31.06 × 10−35.52 × 10−72.13 × 10−3NANANANA1.00 × 10−22.42 × 10−25.52 × 10−73.42 × 10−2NANANANA
Total3.09 × 10−17.01 × 10−27.21 × 10−43.80 × 10−11.08 × 10−54.46 × 10−51.58 × 10−115.54 × 10−58.22 × 10−11.11 × 1017.21 × 10−41.19 × 1012.51 × 10−58.54 × 10−73.95 × 10−122.60 × 10−5
Region IX
Cd1.87 × 10−31.64 × 10−38.46 × 10−53.60 × 10−3NANA5.22 × 10−135.22 × 10−134.30 × 10−21.53 × 10−28.46 × 10−55.84 × 10−2NANA1.31 × 10−131.31 × 10−13
As1.56 × 10−11.33 × 10−11.61 × 10−32.91 × 10−12.41 × 10−52.06 × 10−53.55 × 10−114.47 × 10−55.10 × 101.46 × 101.61 × 10−36.56 × 105.62 × 10−53.93 × 10−78.88 × 10−125.66 × 10−5
Cu1.24 × 10−33.53 × 10−56.34 × 10−71.27 × 10−3NANANANA1.16 × 10−28.09 × 10−46.34 × 10−71.24 × 10−2NANANANA
Pb8.47 × 10−22.41 × 10−34.70 × 10−58.72 × 10−28.64 × 10−7NA6.28 × 10−138.64 × 10−77.91 × 10−15.53 × 10−24.70 × 10−58.46 × 10−12.02 × 10−61.41 × 10−81.57 × 10−132.03 × 10−6
Cr4.23 × 10−24.83 × 10−26.54 × 10−49.13 × 10−22.18 × 10−5NA1.88 × 10−92.18 × 10−51.11 × 103.95 × 10−16.54 × 10−41.50 × 105.08 × 10−51.42 × 10−54.71 × 10−106.50 × 10−5
Zn7.50 × 10−42.14 × 10−53.86 × 10−77.72 × 10−4NANANANA7.00 × 10−34.90 × 10−43.86 × 10−77.49 × 10−3NANANANA
Hg1.69 × 10−36.88 × 10−48.70 × 10−72.38 × 10−3NANANANA1.58 × 10−21.58 × 10−28.70 × 10−73.15 × 10−2NANANANA
Total2.89 × 10−11.86 × 10−12.39 × 10−34.78 × 10−14.67 × 10−52.06 × 10−51.92 × 10−96.73 × 10−57.07 × 101.94 × 102.39 × 10−39.01 × 101.09 × 10−41.46 × 10−54.80 × 10−101.24 × 10−4
Region X
Cd8.88 × 10−37.79 × 10−34.01 × 10−41.71 × 10−2NANA2.48 × 10−122.48 × 10−122.04 × 10−17.27 × 10−24.01 × 10−42.77 × 10−1NANA6.19 × 10−136.19 × 10−13
As3.41 × 10−12.92 × 10−13.51 × 10−36.37 × 10−15.26 × 10−54.50 × 10−57.77 × 10−119.77 × 10−51.11 × 1013.18 × 103.51 × 10−31.43 × 1011.23 × 10−48.60 × 10−71.94 × 10−111.24 × 10−4
Cu2.89 × 10−38.24 × 10−51.48 × 10−62.97 × 10−3NANANANA2.70 × 10−21.89 × 10−31.48 × 10−62.89 × 10−2NANANANA
Pb1.76 × 10−15.03 × 10−39.78 × 10−51.82 × 10−11.80 × 10−6NA1.31 × 10−121.80 × 10−61.65 × 101.15 × 10−19.78 × 10−51.76 × 104.20 × 10−62.94 × 10−83.27 × 10−134.23 × 10−6
Cr2.47 × 10−22.17 × 10−23.35 × 10−44.67 × 10−21.11 × 10−5NA9.64 × 10−101.11 × 10−52.02 × 10−15.66 × 10−13.35 × 10−47.69 × 10−12.60 × 10−57.28 × 10−62.41 × 10−103.33 × 10−5
Zn8.02 × 10−42.29 × 10−54.13 × 10−78.25 × 10−4NANANANA7.48 × 10−35.24 × 10−44.13 × 10−78.01 × 10−3NANANANA
Hg2.63 × 10−31.07 × 10−31.35 × 10−63.70 × 10−3NANANANA2.45 × 10−22.45 × 10−21.35 × 10−64.91 × 10−2NANANANA
Total5.58 × 10−13.27 × 10−14.35 × 10−38.89 × 10−16.56 × 10−54.50 × 10−51.05 × 10−91.11 × 10−41.33 × 1013.97 × 104.35 × 10−31.72 × 1011.53 × 10−48.17 × 10−62.61 × 10−101.61 × 10−4

References

  1. Wcisło, E.; Bronder, J.; Bubak, A.; Rodríguez-Valdés, E.; Gallego, J.L.R. Human health risk assessment in restoring safe and productive use of abandoned con-taminated sites. Environ. Int. 2016, 94, 436–448. [Google Scholar] [CrossRef] [PubMed]
  2. Xia, Y.; Gao, T.; Liu, Y.; Wang, Z.; Liu, C.; Wu, Q.; Qi, M.; Lv, Y.; Li, F. Zinc isotope revealing zinc’s sources and transport processes in karst region. Sci. Total Environ. 2020, 724, 138191. [Google Scholar] [CrossRef]
  3. Chen, H.; Teng, Y.; Lu, S.; Wang, Y.; Wang, J. Contamination features and health risk of soil heavy metals in China. Sci. Total Environ. 2015, 512, 143–153. [Google Scholar] [CrossRef] [PubMed]
  4. Pu, W.; Sun, J.; Zhang, F.; Wen, X.; Liu, W.; Huang, C. Effects of copper mining on heavy metal contamination in a rice agrosystem in the Xiaojiang River Basin, southwest China. Acta Geochim. 2019, 38, 753–773. [Google Scholar] [CrossRef]
  5. Wang, Y.; Guo, G.; Zhang, D.; Lei, M. An integrated method for source apportionment of heavymetal(loid)s in agricultural soils andmodel uncertainty analysis. Environ. Pollut. 2021, 276, 116666. [Google Scholar] [CrossRef] [PubMed]
  6. Greenpeace; Lishui Institute of Ecology and Environment of Nanjing University (LIEEN). Problems and Countermeasures in the Development and Utilization of Polluted Urban Land in China. Report. 2019. Available online: https://www.greenpeace.org.cn/redeveloping-the-polluted-land-under-chinas-cities-problems-and-solutions-report (accessed on 10 September 2021).
  7. Wu, J.; Xiong, Y.; Ge, Y.; Yuan, W. A sustainability assessment-based methodology for the prioritization of contaminated site risk management options. Environ. Sci. Pollut. Res. 2021, 29, 7503–7513. [Google Scholar] [CrossRef]
  8. Zhang, Y.; Guo, Z.; Peng, C.; Deng, H.; Xiao, X. A questionnaire based probabilistic risk assessment (PRA) of heavy metals in urban and suburban soils under different land uses and receptor populations. Sci. Total Environ. 2021, 793, 148525. [Google Scholar] [CrossRef]
  9. Wang, X.; Zhang, C.; Wang, C.; Zhu, Y.; Cui, Y. Probabilistic-fuzzy risk assessment and source analysis of heavy metals in soil considering uncertainty: A case study of Jinling Reservoir in China. Ecotoxicol. Environ. Saf. 2021, 222, 112537. [Google Scholar] [CrossRef]
  10. State Council of the People’s Republic of China (SCC). The National Soil Pollution Prevention and Treatment Action Plan of China. 2016. Available online: http://www.gov.cn/zhengce/content/2016-05/31/content_5078377.htm (accessed on 27 August 2017). (In Chinese)
  11. Ministry of Environmental Protection (MEP) and Ministry of Land and Resources (MLR) of P. R. China. The National Soil Pollution Survey Bulletin; 2014. Available online: http://www.sdpc.gov.cn/fzgggz/ncjj/zhdt/201404/t20140418_607888.html (accessed on 17 April 2021). (In Chinese)
  12. Ministry of Environmental Protection (MEP). Announcements on the Release of Five National Environmental Protection Standards. 2014. Available online: http://www.mee.gov.cn/gkml/hbb/bgg/201402/t20140224_268193.htm (accessed on 25 January 2022). (In Chinese)
  13. Teng, Y.; Wu, J.; Lu, S.; Wang, Y.; Jiao, X.; Song, L. Soil and soil environmental quality monitoring in China: A review. Environ. Int. 2014, 69, 177–199. [Google Scholar] [CrossRef]
  14. Wu, J.; Li, J.; Teng, Y.; Chen, H.; Wang, Y. A partition computing-based positive matrix factorization (PC-PMF) approach for the source apportionment of agricultural soil heavy metal contents and associated health risks. J. Hazard. Mater. 2020, 388, 121766. [Google Scholar] [CrossRef]
  15. Alvarez-Guerra, M.; Viguri, J.R.; Voulvoulis, N. A multicriteria-based methodology for site prioritisation in sediment management. Environ. Int. 2009, 35, 920–930. [Google Scholar] [CrossRef] [PubMed]
  16. Wu, J.; Teng, Y.; Lu, S.; Wang, Y.; Jiao, X. Evaluation of soil contamination indices in a mining area of Jiangxi, China. PLoS ONE 2014, 9, e112917. [Google Scholar] [CrossRef] [PubMed]
  17. Li, P.; Lin, C.; Cheng, H.; Duan, X.; Lei, K. Contamination and health risks of soil heavy metals around a lead/zinc smelter in southwestern China. Ecotoxicol. Environ. Saf. 2015, 113, 391–399. [Google Scholar] [CrossRef] [PubMed]
  18. Teng, Y.; Li, J.; Wu, J.; Lu, S.; Wang, Y.; Chen, H. Environmental distribution and associated human health risk due to trace elements and organic compounds in soil in Jiangxi province, China. Ecotoxicol. Environ. Saf. 2015, 122, 406–416. [Google Scholar] [CrossRef]
  19. Wu, Q.; Leung, J.Y.S.; Geng, X.; Chen, S.; Huang, X.; Li, H.; Huang, Z.; Zhu, L.; Chen, J.; Lu, Y. Heavy metal contamination of soil and water in the vicinity of an abandoned e-waste recycling site: Implications for dissemination of heavy metals. Sci. Total Environ. 2015, 506, 217–225. [Google Scholar] [CrossRef]
  20. Song, Y.; Hou, D.; Zhang, J.; O’Connor, D.; Li, G.; Gu, Q.; Li, S.; Liu, P. Environmental and socio-economic sustainability appraisal of contaminated land remediation strategies: A case study at a mega-site in China. Sci. Total Environ. 2018, 610, 391–401. [Google Scholar] [CrossRef]
  21. Huang, J.; Wu, Y.; Sun, J.; Li, X.; Geng, X.; Zhao, M.; Sun, T.; Fan, Z. Health risk assessment of heavy metal(loid)s in park soils of the largest megacity in China by using Monte Carlo simulation coupled with Positive matrix factorization model. J. Hazard. Mater. 2021, 415, 125629. [Google Scholar] [CrossRef]
  22. Świdwa-Urbańska, J.; Batlle-Sales, J. Data quality oriented procedure, for detailed mapping of heavy metals in urban topsoil as an approach to human health risk assessment. J. Environ. Manag. 2021, 295, 113019. [Google Scholar] [CrossRef]
  23. Wang, S.; Gao, Z.; Zhang, Y.; Zhang, H.; Wu, Z.; Jiang, B.; Liu, Y.; Dong, H. Source and Health Risk Assessment of Heavy Metals in Soil–Ginger System in the Jing River Basin of Shandong Province, North China. Int. J. Environ. Res. Public Health 2021, 18, 6749. [Google Scholar] [CrossRef]
  24. Ekosse, G.I.; Nkeng, G.E.; Bukalo, N.; Oyebanjo, O. Geophagic clays from Cameroon: Provenance, metal contamination and health risk assessment. Int. J. Environ. Res. Public Health 2021, 18, 8315. [Google Scholar] [CrossRef]
  25. Zhang, B.; Liu, L.; Huang, Z.; Hou, H.; Zhao, L.; Sun, Z. Application of stochastic model to assessment of heavy metal(loid)s source apportionment and bio-availability in rice fields of karst area. Sci. Total Environ. 2021, 793, 148614. [Google Scholar] [CrossRef] [PubMed]
  26. China National Environmental Monitoring Centre (CNEMC). The Background Values of Chinese Soil; Environmental Science Press of China: Beijing, China, 1990. [Google Scholar]
  27. Zhao, Y.; Deng, Q.; Lin, Q.; Zeng, C.; Zhong, C. Cadmium source identification in soils and high-risk regions predicted by geographical detector method. Environ. Pollut. 2020, 263, 114338. [Google Scholar] [CrossRef] [PubMed]
  28. Cao, S.; Duan, X.; Ma, Y.; Zhao, X.; Qin, Y.; Liu, Y.; Li, S.; Zheng, B.; Wei, F. Health benefit from decreasing exposure to heavy metals and metalloid after strict pollution control measures near a typical river basin area in China. Chemosphere 2017, 184, 866–878. [Google Scholar] [CrossRef] [PubMed]
  29. Liang, J.; Feng, C.; Zeng, G.; Gao, X.; Zhong, M.; Li, X.; Li, X.; He, X.; Fang, Y. Spatial distribution and source identification of heavy metals in surface soils in a typical coal mine city, Lianyuan, China. Environ. Pollut. 2017, 225, 681–690. [Google Scholar] [CrossRef] [PubMed]
  30. USEPA. Inductively Coupled Plasma-Mass Spectrometry; USEPA: Washington, DC, USA, 1998.
  31. USEPA. Risk Assessment Guidance for Superfund; Human Health Evaluation Manual Part A, Volume I; Office of Emergency and Remedial Response, US Environmental Protection Agency: Washington, DC, USA, 1989; p. 1.
  32. USEPA. Risk Assessment Guidance for Superfund; Volume I: Human Health Evaluation Manual, Part E, Supplemental Guidance for Dermal Risk Assessment; Office of Superfund Remediation and Technology Innovation, U.S. Environmental Protection Agency: Washington, DC, USA, 2004. Available online: https://nepis.epa.gov/Exe/ZyNET.exe/P1002UOM.TXT?ZyActionD=ZyDocument&Client=EPA&Index=2006+Thru+2010&Docs=&Query=&Time=&EndTime=&SearchMethod=1&TocRestrict=n&Toc=&TocEntry=&QField=&QFieldYear=&QFieldMonth=&QFieldDay=&IntQFieldOp=0&ExtQFieldOp=0&XmlQuery=&File=D%3A%5Czyfiles%5CIndex%20Data%5C06thru10%5CTxt%5C00000006%5CP1002UOM.txt&User=ANONYMOUS&Password=anonymous&SortMethod=h%7C-&MaximumDocuments=1&FuzzyDegree=0&ImageQuality=r75g8/r75g8/x150y150g16/i425&Display=hpfr&DefSeekPage=x&SearchBack=ZyActionL&Back=ZyActionS&BackDesc=Results%20page&MaximumPages=1&ZyEntry=1&SeekPage=x&ZyPURL (accessed on 25 January 2022).
  33. USEPA. Risk Assessment Guidance for Superfund; Volume I: Human Health Evaluation Manual, Part F, Supplemental Guidance for Inhalation Risk Assessment; Office of Superfund Remediation and Technology Innovation, U.S. Environmental Protection Agency: Washington, DC, USA, 2009; p. 1.
  34. USEPA. Risk Assessment Guidance for Superfund (RAGS); Volume I: Human Health Evaluation Manual, Part E, Supplemental Guidance for Dermal Risk Assessment; Office of Superfund Remediation and Technology Innovation, U.S. Environmental Protection Agency: Washington, DC, USA, 2009; p. 1. Available online: https://www.epa.gov/risk/risk-assessment-guidance-superfund-rags-part-e (accessed on 25 January 2022).
  35. Li, J.; Fan, J.; Jiang, J.; Zhai, Y.; Luo, Z.; Zhang, Z.; Wu, J. Human health risk assessment of soil in an abandoned arsenic plant site: Implications for contaminated site remediation. Environ. Earth Sci. 2019, 78, 673. [Google Scholar] [CrossRef]
  36. USEPA. Superfund Supplemental Guidance to RAGS: Region 4 Bulletins, Human Health Risk Assessment Bulletins; USEPA: Washington, DC, USA, 2000.
  37. CEHT. Technical Report: Development of Cleanup Target Levels (CTLs) for Chapter 62–777, F.A.C; Division of Waste Management Florida Department of Environmental Protection, Center for Environmental & Human Toxicology, University of Florida: Gainesville, FL, USA, 2005. [Google Scholar]
  38. USEPA. Supplemental Guidance for Developing Soil Screening Levels for Superfund Sites; OSWER 9355.4-24; Office of Solid Waste and Emergency Response: Washington DC, USA, 2002.
  39. USEPA (United States Environmental Protection Agency). Regional Screening Levels Frequent Questions. 2015. Available online: https://www.epa.gov/risk/regional-screening-levelsfrequent-questions-november-2015#FQ13 (accessed on 11 November 2015).
  40. USEPA. Recommendations of the Technical Review Workgroup for Lead for an Approach to Assessing Risks Associated with Adult Exposure to Lead in Soil; EPA-540-R-03-001; Technical Review Workgroup for Lead: Washington, DC, USA, 2003.
  41. USEPA. Update of the Adult Lead Methodology’s Default Baseline Blood Lead Concentration and Geometric Standard Deviation Parameters; OSWER 9200.2-82; Office of Waste and Emergency Response: Washington, DC, USA, 2009.
  42. Chinese Ministry of Environmental Protection (CMEP). Environmental Quality Standard for Soils; (GB15618–1995); CMEP: Beijing, China, 1995. (In Chinese)
  43. Chinese Ministry of Environmental Protection (CMEP). Technical Guidelines for Risk Assessment of Contaminated Sites; HJ 25.3—2014; CMEP: Beijing, China, 2014. (In Chinese)
  44. Jacinthe, P.-A.; Lal, R. Spatial variability of soil properties and trace gas fluxes in reclaimed mine land of southeastern Ohio. Geoderma 2006, 136, 598–608. [Google Scholar] [CrossRef]
  45. Guan, X.; Yang, P.; Lv, Y. Analysis on spatial variability of soil properties based on multifractal theory. J. Basic Sci. Eng. 2011, 19, 712–720. [Google Scholar]
  46. Yang, Q.; Li, Z.; Lu, X.; Duan, Q.; Huang, L.; Bi, J. A review of soil heavy metal pollution from industrial and agricultural regions in China: Pollution and risk assessment. Sci. Total Environ. 2018, 642, 690–700. [Google Scholar] [CrossRef]
  47. Fernández-Caliani, J.C. Risk-based assessment of multimetallic soil pollution in the industrialized peri-urban area of Huelva, Spain. Environ. Geochem. Health 2012, 34, 123–139. [Google Scholar] [CrossRef]
  48. Ordóñez, A.; Álvarez, R.; Loredo, J. Asturian mercury mining district (Spain) and the environment: A review. Environ. Sci. Pollut. Res. 2013, 20, 7490–7508. [Google Scholar] [CrossRef]
  49. Li, Z.; Ma, Z.; van der Kuijp, T.J.; Yuan, Z.; Huang, L. A review of soil heavy metal pollution from mines in China: Pollution and health risk assessment. Sci. Total Environ. 2014, 468–469, 843–853. [Google Scholar] [CrossRef] [PubMed]
  50. Km, M.L.K.; Tu, C.; Zhang, W.H.; Cai, Y.; Kennelley, E.D. A fern that hyperaccumulates arsenic: A hardy, versatile, fast-growing plant helps to remove arsenic contaminated soil. Nature 2001, 409, 579. [Google Scholar]
Figure 1. Distribution of soil sampling sites and contaminated site locations.
Figure 1. Distribution of soil sampling sites and contaminated site locations.
Sustainability 14 03663 g001
Figure 2. HHRA values for different contact pathways for children and adults.
Figure 2. HHRA values for different contact pathways for children and adults.
Sustainability 14 03663 g002
Table 2. Descriptive statistics of soil heavy metals in different regions.
Table 2. Descriptive statistics of soil heavy metals in different regions.
RegionnConcentrationCdAsCuPbCrZnHg
I11Max0.7825.208.7868.900.1062.202.30 × 10−4
Min0.020.250.010.050.010.095.00 × 10−5
Median0.142.801.020.050.028.735.00 × 10−5
33II2Max0.011.050.610.050.030.643.00 × 10−3
Min3.00 × 10−30.040.020.050.010.362.50 × 10−3
Median0.010.540.320.050.020.503.00 × 10−3
3III11Max10.805810.00-74.40-653.006.46
Min0.1012.50-10.80-49.500.39
Median1.17151.00-22.90-131.000.71
IV8Max17.303915.00-165.00-1554.005.02
Min4.4952.50-73.30-682.000.73
Median13.401290.00-114.00-985.001.67
V6Max14.60530.00-116.00-882.001.19
Min0.90527.0032.4023.80-113.00-
Median4.861526.0063.9079.60-445.00-
VI8Max22.50164,000.00-429.00-833.0015.50
Min3.66377.00-24.30-129.000.43
Median12.403895.00-53.90-368.003.29
VII15Max1.981350.00-43.60-192.0012.70
Min0.0811.50-8.81-22.900.14
Median0.38184.00-14.30-49.501.02
VIII3Max0.0932.40-49.50-70.200.33
Min0.014.39-42.80-27.800.12
Median0.058.28-45.70-62.100.25
IX20Max9.1078.00269.002030.00212.00611.001.40
Min0.1014.703.102.7050.3043.800.10
Median0.6525.4018.1018.4085.00121.500.30
X28Max56.00463.00484.007600.00151.00978.006.27
Min0.1710.504.002.8020.0018.800.10
Median1.6738.0076.5030.2039.00105.500.29
Total112Max56.00164,000.00776.007600.00212.005078.0015.50
Min3.00 × 10−30.040.010.050.010.095.00 × 10−5
Note: The units of heavy metals are mg/kg, n: the sample size.
Table 3. Summary of soil risk-based remedial levels (RBRLs: mg/kg) for regions I–X.
Table 3. Summary of soil risk-based remedial levels (RBRLs: mg/kg) for regions I–X.
Regions I–XAdultsChildren
RBRLs for Noncancer RiskRBRLs for Cancer RiskRBRLs for Noncancer RiskRBRLs for Cancer Risk
ElementOralDermalInhalationOralDermalInhalationOralDermalInhalationOralDermalInhalation
Cd1.75 × 1041.60 × 1031.42 × 1043.38 × 101NA6.57 × 1047.82 × 1013.43 × 101.42 × 1041.45 × 101NA2.63 × 105
As5.26 × 1036.40 × 1022.13 × 1041.42 × 1021.70 × 1022.74 × 1042.35 × 1013.65 × 101~8.61 × 1012.13 × 1046.08 × 1015.10 × 1031.10 × 105
Cu7.01 × 1057.68 × 1055.70 × 107NANANA3.13 × 1034.12 × 1035.70 × 107NANANA
Cr5.26 × 1043.84 × 1034.06 × 1045.07 × 10NA9.85 × 1032.35 × 1022.06 × 1014.06 × 1042.17 × 10NA3.94 × 104
Zn5.26 × 1061.92 × 1074.25 × 108NANANA2.35 × 1042.06 × 1044.25 × 108NANANA
Hg5.26 × 1031.34 × 1034.25 × 105NANANA2.35E × 1017.20 × 104.25 × 105NANANA
Note: The background values of arsenic in regions I~X are 36.5 mg/kg, 39.1 mg/kg, 58.1 mg/kg, 86.1 mg/kg, 62.3 mg/kg, 43.2 mg/kg, 76.1 mg/kg, 66.4 mg/kg, 70.1 mg/kg, 61.2 mg/kg, respectively. The background values of lead in regions I~X are 65.1 mg/kg, 79.1 mg/kg, 67.1 mg/kg, 68.1 mg/kg, 76.3 mg/kg, 54.2 mg/kg, 81.1 mg/kg, 46.4 mg/kg, 80.1 mg/kg, and 81.2 mg/kg, respectively.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wu, J.; Jia, R.; Xuan, H.; Zhang, D.; Zhang, G.; Xiao, Y. Priority Soil Pollution Management of Contaminated Site Based on Human Health Risk Assessment: A Case Study in Southwest China. Sustainability 2022, 14, 3663. https://doi.org/10.3390/su14063663

AMA Style

Wu J, Jia R, Xuan H, Zhang D, Zhang G, Xiao Y. Priority Soil Pollution Management of Contaminated Site Based on Human Health Risk Assessment: A Case Study in Southwest China. Sustainability. 2022; 14(6):3663. https://doi.org/10.3390/su14063663

Chicago/Turabian Style

Wu, Jin, Ruitao Jia, Hao Xuan, Dasheng Zhang, Guoming Zhang, and Yuting Xiao. 2022. "Priority Soil Pollution Management of Contaminated Site Based on Human Health Risk Assessment: A Case Study in Southwest China" Sustainability 14, no. 6: 3663. https://doi.org/10.3390/su14063663

APA Style

Wu, J., Jia, R., Xuan, H., Zhang, D., Zhang, G., & Xiao, Y. (2022). Priority Soil Pollution Management of Contaminated Site Based on Human Health Risk Assessment: A Case Study in Southwest China. Sustainability, 14(6), 3663. https://doi.org/10.3390/su14063663

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop