Enrichment, Source Apportionment and Health Risk Assessment of Soil Potentially Harmful Elements Associated with Different Land Use in Coastal Tidelands Reclamation Area, Eastern China.

Coastal tidelands are important ecological habitat resources and valuable resources for agricultural land reclamation. Enrichment of potentially harmful elements (PHEs) in soil caused by anthropogenic activity is an important factor implicated in the ecological deterioration of soil in China. A total of 54 soil sample sites were selected from a 30-year reclaimed tideland and an adjoining coastal wetland. Descriptive and multivariate statistical analyses were performed to describe the enrichment, source, health risk status of eight PHEs (As, Co, Cr, Cu, Mn, Ni, Pb, and Zn) after long-term reclamation. Results indicated that after 30 years of reclamation, most soil PHEs are slightly enriched, whereas no serious threat of environmental pollution was observed. In the reclamation area, the enrichment of PHEs in the aquaculture land, industrial land, and cropland was relatively high compared with other land use types, such as tideland and halophyte land. The source analysis divided the PHEs into five categories: (1) Cu; (2) Co and Mn; (3) Cr; (4) As and Pb; (5) Zn and Ni. Cu was completely derived from natural parent materials and other elements were governed by both weathering of parent rock and human activities, including agricultural activities, industrial production, and transportation emissions. The health risk assessment showed that the soil PHEs potentially had no non-carcinogenic risk to the public, but there was an acceptable probability to have cancer due to Cr and As. Meanwhile, children are more susceptible to harm from the PHEs in soil than adults. According to the economic and social development situation in the coastal region, it is necessary to pay attention to the environmental threats of PHEs enrichment.


Introduction
The enrichment of potentially harmful elements (PHEs) in soil has been widely concerned as a significant ecological and environmental issue [1]. The bedrock geochemistry and human activities are in Jiangsu from 2009 to 2020. The high-intensity, rapid reclamation of coastal zones within a short time period has affected the physicochemical properties of soil in tidelands and even accelerated the rapid enrichment of pollutants [14,27]. There are over ten ports along the coast of Jiangsu Province, which undertake the transportation of numerous import and export commodities and energy (coal, petroleum and metal ore, etc.). Furthermore, port-centered industries such as metallurgy and petrochemicals, electricity, cement, and papermaking have also developed rapidly, which have posed huge threats to the soil ecosystem. Hence, the investigation of characteristics and trends of PHEs enrichment in the soil after coastal reclamation is important to clarify the influence of anthropogenic activities on tidelands and to understand and assess the ecological and environmental effects of coastal reclamation.
In this study, multivariate statistical methods and a US EPA soil health risk assessment model were applied to explore the enrichment, source, and health risk of PHEs in a long-term reclaimed coastal area. The objectives of this study are (1) to explore the variation of soil PHEs under different land use; (2) to identify the source of soil PHEs in a coastal reclamation area; and (3) to assess the level of health risk status after long-term reclamation. Our study is designated to provide a scientific basis for the monitoring and pollution control of the tideland.

Study Area
The study area is located in the coast of the Dongtai City (120 • 78 86.7"-120 • 90 29.2" E; 32 • 95 22.3"-33 • 00 38.9" N) with an area of 200 km 2 , which closely adjoins Yancheng National Nature Reserve, Jiangsu Province ( Figure 1). The regional climate is characterized as a northern humid subtropical monsoon climate with four distinct seasons, and the annual average rainfall and temperature are 1035 mm and 15.0 • C, respectively. The soils in the study area are mainly derived from the modern marine and fluvial sediments, which belong to sandy loam and are classified as coastal saline-sodic soils. Actually, the study area is considered one of the fastest silting regions in eastern China, and the tidal land is currently expanding seaward at a rate of 200 m per year [20]. The reclamation of the study area mainly depended on different land uses to carry out the process of desalination after the construction of seawalls. The land-use/cover types varied with increasing reclamation years and occurred in three stages [28]. In the early period of reclamation, newly reclaimed tidal land changed to aquaculture ponds and salt-tolerant crops, such as Sesbania cannabina. Over 10-30 years, the land-use pattern gradually turned into cropland with an increase of soil fertility and decrease of soil salinity, and some crops were planted in this period, such as cotton, wheat, and corn. After 30 years of reclamation, soil quality basically met the cultivation for most crops with consistent cultivation and fertilization, such as paddy, wheat, corn, rapeseed, and broad bean. In addition, driven by the demands of economic development, some reclaimed land has been converted into industrial land to support industries such as power generation, textiles, chemicals, and automobile infrastructure.

Sampling Design and Chemical Analysis
Soil sites were chosen in mid-August 2014 by using a uniform grid method in the study area. In total, 54 soil sampling sites were collected, including 5 sites respectively in tideland and halophyte land in coastal wetlands, and 44 sites (industrial land (10), forestland (6), residential land (7), aquaculture land (6), cropland (8), and vegetable land (7)) in the areas reclaimed for 30 years. Both surface (0-20 cm) and subsurface (20-40 cm) samples were collected from each sampling site and recorded using GPS ( Figure 2). Each soil sample consisted of 4-8 subsamples collected in a surrounding 10-m area. After air-drying at room temperature, the soil samples were grounded and sieved through a 2 mm nylon sieve to remove animal and plant residues, as well as other impurities, and then sieved through a 0.149 mm nylon sieve as described by Xie et al. [29]. From the sieved soil, 10-15 g preliminarily treated soil samples were used to form soil tablets. Eight PHEs (As, Co, Cr, Cu, Mn, Ni, Pb, and Zn) were analyzed according to X-ray fluorescence spectrometry methods, as described by Xia et al. [30].

Sampling Design and Chemical Analysis
Soil sites were chosen in mid-August 2014 by using a uniform grid method in the study area. In total, 54 soil sampling sites were collected, including 5 sites respectively in tideland and halophyte land in coastal wetlands, and 44 sites (industrial land (10), forestland (6), residential land (7), aquaculture land (6), cropland (8), and vegetable land (7)) in the areas reclaimed for 30 years. Both surface (0-20 cm) and subsurface (20-40 cm) samples were collected from each sampling site and recorded using GPS ( Figure 2). Each soil sample consisted of 4-8 subsamples collected in a surrounding 10-m area. After air-drying at room temperature, the soil samples were grounded and sieved through a 2 mm nylon sieve to remove animal and plant residues, as well as other impurities, and then sieved through a 0.149 mm nylon sieve as described by Xie et al. [29]. From the sieved soil, 10-15 g preliminarily treated soil samples were used to form soil tablets. Eight PHEs (As, Co, Cr, Cu, Mn, Ni, Pb, and Zn) were analyzed according to X-ray fluorescence spectrometry methods, as described by Xia et al. [30].

Descriptive and Multivariate Statistical Analysis
Descriptive analysis (range, mean, standard deviation, coefficient of variation) was applied to understand the distribution and variation of PHEs in the study area. The differences of PHEs among

Sampling Design and Chemical Analysis
Soil sites were chosen in mid-August 2014 by using a uniform grid method in the study area. In total, 54 soil sampling sites were collected, including 5 sites respectively in tideland and halophyte land in coastal wetlands, and 44 sites (industrial land (10), forestland (6), residential land (7), aquaculture land (6), cropland (8), and vegetable land (7)) in the areas reclaimed for 30 years. Both surface (0-20 cm) and subsurface (20-40 cm) samples were collected from each sampling site and recorded using GPS ( Figure 2). Each soil sample consisted of 4-8 subsamples collected in a surrounding 10-m area. After air-drying at room temperature, the soil samples were grounded and sieved through a 2 mm nylon sieve to remove animal and plant residues, as well as other impurities, and then sieved through a 0.149 mm nylon sieve as described by Xie et al. [29]. From the sieved soil, 10-15 g preliminarily treated soil samples were used to form soil tablets. Eight PHEs (As, Co, Cr, Cu, Mn, Ni, Pb, and Zn) were analyzed according to X-ray fluorescence spectrometry methods, as described by Xia et al. [30].

Descriptive and Multivariate Statistical Analysis
Descriptive analysis (range, mean, standard deviation, coefficient of variation) was applied to understand the distribution and variation of PHEs in the study area. The differences of PHEs among

Descriptive and Multivariate Statistical Analysis
Descriptive analysis (range, mean, standard deviation, coefficient of variation) was applied to understand the distribution and variation of PHEs in the study area. The differences of PHEs among different land-use types were compared using the analysis of variance (ANOVA) procedure and the mean comparisons were performed using the Fisher's least significant difference (LSD) test with a probability defined at 0.05 [24]. Pearson's correlation analysis and cluster analysis (CA) are effective to identify the source of PHEs in soil [15]. Pearson's correlation analysis could determine whether the sources of PHEs are the same through the correlation between PHEs, and CA could provide variable grouping results, which corroborates their correctness with correlation analysis results. Descriptive and multivariate statistical analysis was conducted in SPSS 22.0 for Windows (SPSS München, Germany).
The enrichment factor (EF) is applied to determine the degree of PHEs enrichment and to differentiate between elements originating from human activities and those from natural provenance by using Equation (1) [1,31].
where Mx and Nref are the concentration of PHEs and reference elements in soil, respectively. In this study, Fe was adopted as a reference because it is one of the abundant elements of soil [31]. The geochemical baseline of Fe (4.5 µg/g) in Jiangsu Province was chosen as reference background values [32]. Generally, EF is applied to differentiate whether the element originates from anthropogenic activities. The element is completely derived from natural processes under EF < 1, while it might be from man-made sources under EF > 1 and the proportion of human sources increased with the increasing of EF [33]. The EF is also applied to estimate polluted degree, which can be divided into smaller enrichment (EF < 2), moderate enrichment (2 ≤ EF < 5), significant enrichment (5 ≤ EF < 20), highly enrichment (20 ≤ EF < 40) and extremely high enrichment (EF ≥ 40) [31,33].

Health Risk Assessment
The long-term accumulation of PHEs in soil may cause a threat to the health of animals and human beings, and increased risk of cancer [31,34]. The methodology used in this study to assess the exposure risks of PHEs to human health in the coastal reclaimed area based on risk assessment guidance for superfund (RAGS) health risk assessment as developed by the United States Environmental Protection Agency (US EPA) [35]. Direct exposure of an individual to PHEs in the soil through three main pathways: ingestion, inhalation, and dermal contact [36,37]. The average daily dose through three pathways can be respectively estimated by the following Equation [8].
Where ADD ing , ADD inh, and ADD derm are the average daily dose from ingestion, inhalation, and dermal contact, respectively (mg·kg −1 ·d −1 ). C soil is the concentration of PHEs in soil (mg·kg −1 ·d −1 ). Related definition and parameters used in the average daily dose of soil metals can be found in Table 1. Dermal absorption factor 0.001 [38] Non-carcinogenic risk assessment: Hazard quotient (HQ) is calculated by the ratio of average daily dose to the corresponding reference dose, which is used to assess the potential non-carcinogenic risk of PHEs under certain exposure pathways. The Hazard Index (HI) is the sum of the HQ and represents the overall potential non-carcinogenic effects posed by all exposure pathways. The HQ of a single element is determined by Equation (5) [35].
where n =8, i = As, Co, Cr, Cu, Mn, Ni, Pb, and Zn, m =3, j = ingestion, inhalation and dermal contact. RfD ij is the reference dose of i element in j exposure pathway (mg·kg −1 ·d −1 ). Carcinogenic risk assessment: Carcinogenic risk (CR) is the probability of cancer, which can be calculated by the product of daily average exposure multiplied by the corresponding slope factor. The CR of a single element is determined by Equation (6) [35]. SF ij is the slope factor i metal in j exposure pathway (mg·kg −1 ·d −1 ). The values of RfD ij and SF ij are shown in Table 2 [41][42][43].
If HQ or HI < 1, the non-carcinogenic risk did not exist to the exposed individual. However, if HQ or HI ≥ 1, the non-carcinogenic effect may occur with a probability which tends to increase with the increase of HQ or HI [8]. If CR < 10 −6 , it will neglect the carcinogenic risk to human health. 10 −6 < CR < 10 −4 is considered an acceptable range and CR > 10 −4 is viewed as unacceptable [44].
PHEs with carcinogenic risk are As, Cd, Cr, Ni and Pb; others are non-carcinogenic metals. "-" in the table means not available.

Descriptive Analysis
The average concentrations of soil PHEs (0-40 cm) showed that As, Co, Cr, Cu, Mn, Ni, Pb, and Zn were 3. 31 (Table 3). Compared with the Chinese Soil Environmental Quality Standard (GB 15618-1995), the average concentrations of PHEs were all within the range. In general, the concentrations of PHEs in the tidal flat reclamation area were relatively low, which is similar to the conclusions of other researches in Jiangsu coastal areas, and might be related to the local geochemical background [45,46]. Lv et al. [47] demonstrated that the background concentrations of Cr, Cu, Ni, Pb, and Zn in marine sediments were significantly lower than that in river alluvial deposits, lagoonal facie sediments, and delta sediments. However, the highest concentrations of Pb and Cr exceeded the limit values of the first and second-grade environmental quality standard for soils, respectively. This showed that partial land use had the opportunity to cause soil pollution, which might pose a threat to the regional natural background or agricultural production. Specifically, the average values of eight elements (As, Co, Cr, Cu, Mn, Ni, Pb, and Zn) did not exceed regional soil background values and the maximum values were close to or exceeded that value. The variation coefficient identified that As has a high coefficient of variation (60.06%), which suggests that the distribution of As is extremely non-uniform. The distribution of Cr, Cu, and Pb have been relatively moderate and non-uniform, owing to their coefficients of variation ranging from 20-30%. Although the average values of PHEs were lower than the risk control standard and regional soil background values, the human health risk of PHEs needs to be of concern due to the strong variation of PHEs (As, Cr, Cu, and Pb) and the difference between the exposure parameters and chronic reference dose of human health risk. Risk control standard China is Environmental quality standard for soils (GB 15618-1995) and has three levels of standards. (1) means the limit value of the soil environmental quality to protect the regional natural ecology and maintain the natural background. (2) means the limit value of the soil to ensure agricultural production and maintain human health; (3) means the limit value of the soil to ensure the agricultural and forestry production and normal plant growth. Mean 1 and Mean 2 represent the average concentration of PHEs in surface and subsurface. S.D 1 and S.D 2 represent standard deviations for the surface and subsurface concentrations.

PHEs Concentrations under Different Land Use
As shown in Figure 3, the average concentrations of As, Cu, Mn, Ni and Zn in surface and subsurface soil were observed to be significantly different from different land-use types (p < 0.05), which indicates that PHEs are greatly affected by the land use. However, soil depth had less effect on the concentrations of determined elements. The concentrations of Co and Cr did not change significantly under different land-use types, whatsoever, in surface or subsurface soil. However, the average concentration of Pb changed significantly in subsurface soil with the highest concentration observed in aquaculture land, and the concentration of As varied significantly under different land-use types with the highest concentration observed in industrial land. The average concentration of Cu in industrial land and aquaculture land was obviously higher than that of tideland and halophyte land, and the differences of other land-use types were not significant. Besides, no significant differences of Mn were observed in surface soil under different land-use patterns, whereas in subsurface soil, it was significantly higher in aquaculture land than that of other land-use types. Compared with the industrial land and aquaculture land, the concentrations of Ni and Zn in tideland, halophyte land, forestland, and vegetable land were relatively lower and changed mildly. The marine sediment sources and the degree of disturbance of human activities revealed that the lower concentration of PHEs was distributed in uncultivated soil (Tideland and Halophyte land) [49,50]. As, Cu, Pb, and Zn of industrial land, residential land, aquaculture land, cropland, and vegetable land were relatively higher than that of forestland, which might be caused by the absorption of trees and the re-entry of toxic elements from industrial activities and agricultural practices [14]. land, and the differences of other land-use types were not significant. Besides, no significant differences of Mn were observed in surface soil under different land-use patterns, whereas in subsurface soil, it was significantly higher in aquaculture land than that of other land-use types. Compared with the industrial land and aquaculture land, the concentrations of Ni and Zn in tideland, halophyte land, forestland, and vegetable land were relatively lower and changed mildly. The marine sediment sources and the degree of disturbance of human activities revealed that the lower concentration of PHEs was distributed in uncultivated soil (Tideland and Halophyte land) [49,50]. As, Cu, Pb, and Zn of industrial land, residential land, aquaculture land, cropland, and vegetable land were relatively higher than that of forestland, which might be caused by the absorption of trees and the re-entry of toxic elements from industrial activities and agricultural practices [14].

Source Analysis of PHEs in Soil
As presented in Figure 4, the EF of measured elements were in the following order: Ni > Mn > Zn > Cr > Co > Pb > Cu > As. The average EF of measured elements ranged from 1.0 to 1.5, except for Cu and As, which indicated that Co, Cr, Mn, Ni, Pb, and Zn were at minimal enrichment level. However, the highest EF of Cr, Mn, Ni, Pb, and Zn were above 1.5 and that of Cr was above 2.0,

Source Analysis of PHEs in Soil
As presented in Figure 4, the EF of measured elements were in the following order: Ni > Mn > Zn > Cr > Co > Pb > Cu > As. The average EF of measured elements ranged from 1.0 to 1.5, except for Cu and As, which indicated that Co, Cr, Mn, Ni, Pb, and Zn were at minimal enrichment level. However, the highest EF of Cr, Mn, Ni, Pb, and Zn were above 1.5 and that of Cr was above 2.0, which was close to the medium enrichment level. This might be affected by anthropogenic activities (e.g., land use). As shown in Figure 5, the halophyte land, tideland, and forestland belong to the low EF value group whereas the cropland and aquaculture land belong to the high EF value group. It was worth noting that the average value of Cu was less than 1.0 in any kind of land-use types, which indicates that it was mainly controlled by the parent material in the study area. Only in industrial land did the EF value of As exceed 1.0, indicating that the As might be mainly affected by the parent material and local industrial activities such as steel smelting might be another reason [48]. For Pb, the mean EF of halophyte land, tideland and forestland were less than 1.0, while that of other land-use types was more than 1.0, especially in aquaculture land and industrial land. Similarly, the mean EF of Co ranged between 1.0-1.2 and was relatively high in aquaculture land and industrial land. Therefore, there was a slight enrichment of Pb and Co in the study area, which might be partly due to the natural processes, and local aquaculture and industrial activities. In addition, the average EF of Cr, Mn, Ni and Zn under different land-use types were more than 1.0, and the cropland and aquaculture land was much higher than that of other land use types, indicating that they were mainly controlled by agricultural and aquaculture activities. which was close to the medium enrichment level. This might be affected by anthropogenic activities (e.g., land use). As shown in Figure 5, the halophyte land, tideland, and forestland belong to the low EF value group whereas the cropland and aquaculture land belong to the high EF value group. It was worth noting that the average value of Cu was less than 1.0 in any kind of land-use types, which indicates that it was mainly controlled by the parent material in the study area. Only in industrial land did the EF value of As exceed 1.0, indicating that the As might be mainly affected by the parent material and local industrial activities such as steel smelting might be another reason [48]. For Pb, the mean EF of halophyte land, tideland and forestland were less than 1.0, while that of other land-use types was more than 1.0, especially in aquaculture land and industrial land. Similarly, the mean EF of Co ranged between 1.0-1.2 and was relatively high in aquaculture land and industrial land. Therefore, there was a slight enrichment of Pb and Co in the study area, which might be partly due to the natural processes, and local aquaculture and industrial activities. In addition, the average EF of Cr, Mn, Ni and Zn under different land-use types were more than 1.0, and the cropland and aquaculture land was much higher than that of other land use types, indicating that they were mainly controlled by agricultural and aquaculture activities.  Correlation analysis was performed to analyze the relationship between PHEs [36]. As depicted in Figure 6, Cr was weakly or uncorrelated with other elements (r < 0.4), which demonstrated the source of Cr was different from other measured elements. Most of the remaining PHEs were significantly correlated (r > 0.5). Especially, As-Cu-Pb had significant correlations (r > 0.5), Mn had significant correlation with Ni, Co (r = 0.62, 0.58), and Ni-Zn-Cu had extremely significant correlations (r > 0.75). As, Cu, and Pb had homology; part of the source of Mn might be the same as Ni and the other part might be the same as Co; Ni, Zn, and Cu had homology. The results of cluster which was close to the medium enrichment level. This might be affected by anthropogenic activities (e.g., land use). As shown in Figure 5, the halophyte land, tideland, and forestland belong to the low EF value group whereas the cropland and aquaculture land belong to the high EF value group. It was worth noting that the average value of Cu was less than 1.0 in any kind of land-use types, which indicates that it was mainly controlled by the parent material in the study area. Only in industrial land did the EF value of As exceed 1.0, indicating that the As might be mainly affected by the parent material and local industrial activities such as steel smelting might be another reason [48]. For Pb, the mean EF of halophyte land, tideland and forestland were less than 1.0, while that of other land-use types was more than 1.0, especially in aquaculture land and industrial land. Similarly, the mean EF of Co ranged between 1.0-1.2 and was relatively high in aquaculture land and industrial land. Therefore, there was a slight enrichment of Pb and Co in the study area, which might be partly due to the natural processes, and local aquaculture and industrial activities. In addition, the average EF of Cr, Mn, Ni and Zn under different land-use types were more than 1.0, and the cropland and aquaculture land was much higher than that of other land use types, indicating that they were mainly controlled by agricultural and aquaculture activities.  Correlation analysis was performed to analyze the relationship between PHEs [36]. As depicted in Figure 6, Cr was weakly or uncorrelated with other elements (r < 0.4), which demonstrated the source of Cr was different from other measured elements. Most of the remaining PHEs were significantly correlated (r > 0.5). Especially, As-Cu-Pb had significant correlations (r > 0.5), Mn had significant correlation with Ni, Co (r = 0.62, 0.58), and Ni-Zn-Cu had extremely significant correlations (r > 0.75). As, Cu, and Pb had homology; part of the source of Mn might be the same as Ni and the other part might be the same as Co; Ni, Zn, and Cu had homology. The results of cluster Correlation analysis was performed to analyze the relationship between PHEs [36]. As depicted in Figure 6, Cr was weakly or uncorrelated with other elements (r < 0.4), which demonstrated the source of Cr was different from other measured elements. Most of the remaining PHEs were significantly correlated (r > 0.5). Especially, As-Cu-Pb had significant correlations (r > 0.5), Mn had significant correlation with Ni, Co (r = 0.62, 0.58), and Ni-Zn-Cu had extremely significant correlations (r > 0.75).
As, Cu, and Pb had homology; part of the source of Mn might be the same as Ni and the other part might be the same as Co; Ni, Zn, and Cu had homology. The results of cluster analysis (Figure 7) indicated that all soil PHEs after reclamation were divided into four categories: (1) Co and Mn; (2) Cr; (3) As and Pb; (4) Cu, Zn and Ni. This was consistent with the results of the correlation analysis, indicating that PHEs in the study area were multi-sourced. analysis ( Figure 7) indicated that all soil PHEs after reclamation were divided into four categories: (1) Co and Mn; (2) Cr; (3) As and Pb; (4) Cu, Zn and Ni. This was consistent with the results of the correlation analysis, indicating that PHEs in the study area were multi-sourced.

Non-Carcinogenic Risk Assessment from Eight PHEs
According to the risk assessment model, the non-carcinogenic risk exposure dose (ADD), hazard quotient (HQi), and hazard index (HI) of PHEs for adults and children through three exposure pathways (ingestion, inhalation, and dermal contact) were calculated. It could be seen from Table 4 that both the HQi and the HI were less than 1, indicating that the PHEs in soil potentially had no noncarcinogenic risk to the public. Compared with the national average of non-carcinogenic risks, the non-carcinogenic risks of soil Cr, Cu, and Ni in the study area were lower than the national average, while the non-carcinogenic risks of As, Pb and Zn were slightly higher than the national average [31]. In general, the daily average exposure dose and non-carcinogenic risk for different people varied analysis ( Figure 7) indicated that all soil PHEs after reclamation were divided into four categories: (1) Co and Mn; (2) Cr; (3) As and Pb; (4) Cu, Zn and Ni. This was consistent with the results of the correlation analysis, indicating that PHEs in the study area were multi-sourced.

Non-Carcinogenic Risk Assessment from Eight PHEs
According to the risk assessment model, the non-carcinogenic risk exposure dose (ADD), hazard quotient (HQi), and hazard index (HI) of PHEs for adults and children through three exposure pathways (ingestion, inhalation, and dermal contact) were calculated. It could be seen from Table 4 that both the HQi and the HI were less than 1, indicating that the PHEs in soil potentially had no noncarcinogenic risk to the public. Compared with the national average of non-carcinogenic risks, the non-carcinogenic risks of soil Cr, Cu, and Ni in the study area were lower than the national average, while the non-carcinogenic risks of As, Pb and Zn were slightly higher than the national average [31]. In general, the daily average exposure dose and non-carcinogenic risk for different people varied

Non-Carcinogenic Risk Assessment from Eight PHEs
According to the risk assessment model, the non-carcinogenic risk exposure dose (ADD), hazard quotient (HQ i ), and hazard index (HI) of PHEs for adults and children through three exposure pathways (ingestion, inhalation, and dermal contact) were calculated. It could be seen from Table 4 that both the HQ i and the HI were less than 1, indicating that the PHEs in soil potentially had no non-carcinogenic risk to the public. Compared with the national average of non-carcinogenic risks, the non-carcinogenic risks of soil Cr, Cu, and Ni in the study area were lower than the national average, while the non-carcinogenic risks of As, Pb and Zn were slightly higher than the national average [31]. In general, the daily average exposure dose and non-carcinogenic risk for different people varied slightly, showing that children were higher than adults were. The average hazard quotients of PHEs decreased in the order of HQ Mn > HQ Cr > HQ As > HQ Pb > HQ Co > HQ Ni > HQ Cu > HQ Zn , and the contribution rates to HI were 48.87%-73.30%, 13.65%-25.48%, 7.52%-16.09%, 1.29%-1.67%, 0.70%-1.49%, 0.13%-0.28%, 0.10%-0.21% and 0.15%-0.16%, respectively. For all PHEs, the average daily exposure dose of ingestion was almost 10 3 times that of inhalation and dermal contact. Similarly, the non-carcinogenic risk of most PHEs through ingestion obviously exceeded the other pathways, and risk caused by dermal contact of Mn also accounted for a large proportion. Therefore, the non-carcinogenic risk was mainly contributed by the exposure pathway of ingestion in Cr and dermal contact in Mn. Although there was no non-carcinogenic risk in the study area, the HQ of children were nearly three times that of adults, indicating that children were more vulnerable to PHEs (Figure 8).

The Influence of Different Land Use on PHEs Enrichment in Coastal Tideland Region
High-intensity reclamations have resulted in drastic changes in land use of tidal flats, which have made ecosystems more fragile and sensitive, including the increased risk of soil pollution [19]. Table 3 indicated that the concentrations of PHEs in surface and subsurface soil layers were almost equal, but the maximum and minimum concentrations of whole sites differed greatly. Furthermore,

Carcinogenic Risk Assessment from As, Co, Ni, Gr and Pb
Due to only As, Co, Ni, Cr, and Pb have carcinogenic slope factor parameters, the ADD ing , RI i, and RI of adults and children through three exposure routes (oral intake, respiratory inhalation, and dermal contact) were evaluated. As shown in Table 5, the daily average exposure of children is higher than that of adults under different exposure routes, and the RI i and RI of children are higher than that of adults. The average RI of PHEs is following the order of RI Cr > RI As > RI Pb > RI Co > RI Ni , with the contribution rate of 83.04%-83.05%, 16.49%-16.52%, 0.42%-0.43%, 0.01%-0.04%, 0.0016%-0.008%, respectively. This indicated that the RI of PHEs (As, Co, Ni, Cr, and Pb) in the soil of the study area is at an acceptable risk level, which is mainly contributed by the oral exposure of Cr and As. From different exposure routes, the daily average exposure of oral intake is 2-3 orders of magnitude higher than that of inhalation and dermal contact. The RI ing of As, Cr, and Pb in oral intake was much higher than that of RI inh and RI derm .
1 "-" in the table means not available.

The Influence of Different Land Use on PHEs Enrichment in Coastal Tideland Region
High-intensity reclamations have resulted in drastic changes in land use of tidal flats, which have made ecosystems more fragile and sensitive, including the increased risk of soil pollution [19]. Table 3 indicated that the concentrations of PHEs in surface and subsurface soil layers were almost equal, but the maximum and minimum concentrations of whole sites differed greatly. Furthermore, soil depth was observed to have less influence on the concentration of PHEs than land use. Soil environmental status of the study area was generally in good condition due to the average concentrations of all PHEs lower than the environmental quality standard for soils in China (Table 3). However, the maximum concentration of some PHEs was close to or exceeded the limits of first grade (As and Pb) and second grade (Cr) of the environmental quality standard, which might be affected by reclamation activities, specifically manifested in the changes of land-use intensity and soil properties [24,26,45]. The traditional evolution pathway of land use in the reclaimed coastal areas in eastern China was mudflat→halophyte land→aquaculture land→arable land→construction land [51], which showed that the land-use intensity gradually increased with the increase of human activity interference. In this study, the concentrations of most PHEs (except Co) were in the following order: aquaculture land > industrial land/cropland > forestland/residential land > vegetable land > tideland/halophyte land (Figure 3). Yao et al. [52] also reported similar conclusions after soil investigation on the reclaimed tideland of Dongtai City, Jiangsu Province, which found there was no significant excess of soil PHEs, while the total ecological risk and polluted level showed the trend following the order: industrial park > residential area > suburb farmland > newly-reclaimed tidal flats. ANOVA (Figure 3) illustrated that the concentrations of As, Cu, Pb and Zn in soil after 30 years of reclamation (industrial land, aquaculture land, farmland, residential land, and vegetable land) were significantly higher than those in uncultivated coastal wetlands (tideland and halophyte land), which might be caused by changes in soil physical and chemical properties under reclamation activities [19]. It is reported that the sand content, salinity, and pH gradually decreased, while the silt and clay content and organic matter content gradually increased following tidal land reclamation, which finally reached a stable state after 30 years reclamation [14,19]. The dynamics of soil physicochemical properties have been proven as the key factors affecting the mobility of PHEs in soil [14]. Specifically, fine-grained sediments are easy to solidify PHEs in soil [53], the increase of SOM promotes PHEs accumulation [1,24], and the soluble PHEs in pore water decreases with the decreasing of soil salinity after reclamation [27].
The sources of PHEs in the study area were complex based on the results of enrichment analysis, correlation analysis, and cluster analysis. The EFs were observed less than 2.0 in the study area, which indicated the sources of soil PHEs were largely contributed by parent materials and partly determined by human activities, such as land use. The enrichment analysis showed Co (0.86-1.39), Cr (0.96-2.01), Mn (1.21-1.75), Ni (1.11-1.84), Pb (0.78-1.85), and Zn (0.97-1.61) were generally in slight enrichment, while Cr, Ni, and Pb in some sampling sites were close to moderate enrichment ( Figure 4). Therefore, it is necessary to pay attention to the pollution of Cr, Ni, and Pb in the future reclamation and utilization of tidal land. The EFs of the determined PHEs in halophyte land, tideland, and forestland were generally classified as low concentrations group, and that in cropland and aquaculture land were classified as high concentrations group. Similarly, Sowan et al. [26] reported that the PHEs contamination in municipality areas, industrial zones, and dockyard were more serious, while shrimp farming and traditional land uses such as salt marshes, paddy fields, orchards, and mangrove forests showed low levels of metals in Pattani Bay. These conclusions proved the intensity of human activity further influenced the soil PHEs enrichment [7,54].
It could be seen in Table 3, the highest concentration of Cu (14.10 mg·kg −1 ) was much smaller than the background of Jiangsu (22.3 mg·kg −1 ), and the EFs of Cu did not exceed 1.0 in all land-use types. This indicates that no significant Cu enrichment was observed in the study area. Actually, the average concentration of Cu (8.10 mg·kg −1 ) was close to the deep sediment of coastal mudflat in Jiangsu (13.0 mg·kg −1 ) [49], which indicated that Cu most likely originated in natural factors and was rarely disturbed by anthropogenic activities [55].
Co and Mn might have the same sources and be slightly accumulated by the influence of human activities, as their correlation coefficient was 0.58 and the average EF ranged from 1 to 1.5. Table 3 showed that the average concentrations of Co and Mn were less than the background parameters, and their concentrations in the surface (8.00 mg·kg −1 , 457.49 mg·kg −1 ) were lower than that in the subsurface layer (8.01 mg·kg −1 , 460.03 mg·kg −1 ). This phenomenon was primarily attributed to the natural leaching, which indicated that Co and Mn were controlled by natural geological processes [15]. Li et al. [7] and Franco-Uría et al. [56] also presented that Co and Mn were originated from a lithogenic source [24]. The EF of Co and Mn in tideland (EF = 1.14, 1.43) and halophyte land (EF = 1.16, 1.39) were observed to be slightly lower than that in aquaculture land (EF = 1.17, 1.47), which also revealed that the Co and Mn background content in the tidal flats under the natural sedimentary process was accounted for high proportions [24]. Co and Mn highly accumulated in aquaculture land compared with other land-use types. The marine aquaculture fish feed usually contains Zn, Cu, Cd, Fe, Mn, Co, Ni, P [55], and it may be a major source of Co and Mn. In addition, the water submerged environment promotes metals and S 2associate and precipitate to insoluble metal sulfides under redox status [57]. Thereby, natural sources mostly contribute to the distribution of Co and Mn, and aquaculture is another source.
Cr seemed to be an isolated element, having low correlations with other elements (r < 0.4) [58]. Even if Cr originates from a natural source which has been extensively mentioned in many pieces of literature [15,59]. In this study, the distribution of Cr was mainly affected by agricultural input as the average EF of Cr in cropland (EF = 1.36) and aquaculture land (EF = 1.30) was significantly higher than that of others. The average concentration of determined Cr (52.50 mg·kg −1 ) in cropland was higher than that in phosphate fertilizers (47.3 mg·kg −1 ), whereas much less than that in sludge (86.73 mg·kg −1 ) [24]. Consequently, the enrichment of Cr in the study area may be the result of the application of fertilizer using sludge as the raw material. Field investigation found printing and dyeing mills and mining enterprises using Cr compounds around the reclamation area [60]. In summary, sewage irrigation and fish feed fish have become important sources of soil Cr in this area [55,61].
The slight accumulation of soil As appeared locally rather than the whole area (EF average = 0.58, EF max = 1.19). As, Cu and Pb had homology as the correlation coefficients among them all exceeded 0.5. It had been concluded that Cu originated from natural factors, which indirectly revealed that As and Pb were partially derived from natural parent materials. The average concentration and EF of As in industrial land were significantly higher than those in other land use types, indicating that the industrial activities might be another source of As. Huang et al. [62] and Khairy et al. [63] reported that As was related to the emissions of industrial flue gas, sewage, and sludge. The average EF of Pb ranged from 0.93 to 1.21 under different land-use types and tended to accumulate in aquaculture land, industrial land, vegetable land, and cropland. The cropland and fishponds were adjacent to each other, which might cause mutual pollution. It is widely reported that pesticides and fertilizers (e.g., phosphate fertilizers, nitrate fertilizers, and organic fertilizers) applied in agricultural soils caused high concentrations of Cd, Pb, and Zn [15,64].
Zn-Ni-Cu were significantly correlated (r > 0.75), which revealed that partly, sources of Zn and Ni might derive from natural parent materials as well as Cu. This was consistent with the result of Sun et al. [59]. Slight enrichment of Zn and Ni occurred in aquaculture land, industrial land, and cropland (1.30 < EF < 1.60, Figure 4). Zn accumulated in fishery sediments may be attributed to the input of fish feed or fish manure [55]. In addition, the application of poultry manure in agricultural land increased Zn accumulation [15,56]. Furthermore, the emission of automobile exhaust and wear of the tires significantly affected the Zn and Ni concentrations in road dust, and then entered the surrounding soil through atmospheric sedimentation [15,63]. Ni might be directly related to the "three wastes" (including some domestic waste) emissions during the production of batteries near the study area as the battery slag and wastewater contains high concentrations of Pb, Cd, Ni, etc. [49].

The PHEs Health Risk in Coastal Tideland Reclamation Area
It is found that the total non-carcinogenic risk of soil PHEs in the study area is at the no non-carcinogenic risk level, while the total carcinogenic risk is at the acceptable risk level. This is consistent with the result that the average concentration of PHEs in Table 3, indicating that the concentration of soil PHEs in the study area has no obvious harm to human health. It is worth noting that the carcinogenic risk index of Cr (1.50 × 10 −5 -2.30 × 10 −5 ) and As (2.99 × 10 −6 -4.58 × 10 −6 ) exceeds the maximum acceptable carcinogenic risk of the human body. This might be attributed to the high slope coefficient of the carcinogenic risk of Cr and As [31,43]. It is reported that As is a carcinogen for lung cancer and skin cancer, and excessive intake of Cr will cause respiratory and digestive diseases [3]. Accordingly, As and Cr should receive more attention in future soil environmental monitoring of coastal reclamation areas. Tideland reclamation accelerated the transformation from bare flat to cultivated land, aquaculture land and construction land, which had an important impact on the health risk of PHEs in soil. In this study, the highest non-carcinogenic risk and the carcinogenic risk appeared in industrial land, followed by farmland and vegetable land (Figure 8). Yang et.al [60] found that the health risks of industrial areas were higher than that of agricultural areas in China. Tables 3 and 4 showed that the ratio of average daily exposure dose between ingestion and inhalation, dermal contact was all close to 10 3 . However, the non-carcinogenic risk was mainly contributed by ingestion and dermal contact of Mn and ingestion of Cr (Table 4). Carcinogenic risk assessment (Table 5) showed that ingestion was the main exposure route for As, Cr, and Pb, while the inhalation was the unique exposure pathway for Co and Ni. This might be due to the increased risk exposure pathway by the reclamation of tideland. The smoke discharged by industry and transportation enters the human body through inhalation, while the rest enters the soil through natural sedimentation and rain. PHEs from industrial activities and agricultural fertilization management accumulate continuously after entering the soil and migrate in different systems (such as water systems and crops), thereby entering the human body through the food chain, resulting in harm to human health [3,65]. Besides, compared with daily average exposure, non-carcinogenic risk index and carcinogenic risk index, it can be found that children are more vulnerable to PHEs in soil than adults (Figure 8), which is similar to previous research results [44,60,65]. It is mainly because the children's immune resistance is much lower than that of adults, which is more sensitive to environmental pollution. In addition, children's exposure to soil due to outdoor activities is higher than that of adults, and children's physical behavior (such as finger sucking) makes their health risk higher than adults [3,65].

Conclusions
The enrichment, sources apportionment, and health risk assessment of eight PHEs (As, Co, Cr, Cu, Mn, Ni, Pb, and Zn) were identified in Jiangsu coastal reclamation area, eastern China. In general, none of the PHEs we concerned about exceeding the environmental quality standard for soils in China. Enrichment analysis showed that As and Cu were not enriched, while other elements were slightly enriched. Land use is observed to have an obvious enrichment effect on most of the PHEs, which resulted in a relatively high level of PHEs enrichment in aquaculture land, industrial land, and cropland compared with other land use types, such as tideland and halophyte land. The correlation analysis and hierarchically cluster divided PHEs into five categories: (1) Cu; (2) Co and Mn; (3) Cr; (4) As and Pb; (5) Zn and Ni. Additionally, it was observed that the Cu was completely derived from natural parent materials; the Mn and Co were mainly controlled by natural sources and aquaculture is another source, and; the Cr was an independent element and was mainly derived from agricultural fertilization and aquaculture activities. However, the As, Pb, Ni, and Zn are governed by both weathering of parent rock and human activities, including agricultural activities, industrial production, and transportation emissions. The health risk assessment showed that the soil PHEs potentially had no non-carcinogenic risk to the public, while there was an acceptable probability to have cancer due to Cr and As. Meanwhile, children were more susceptible to harm of soil PHEs than that of adults. Although the soil PHEs in the study area had not caused significant negative effects on the environment and individuals, the uncertainty of some PHEs (Pb, As, and Cr) was still high. That is to say, it may increase PHEs pollution and further increase the probability of human body damage if land management could not be rationally planned in the future.