Source Apportionment of Topsoil Heavy Metals and Associated Health and Ecological Risk Assessments in a Typical Hazy City of the North China Plain

: The North China Plain (NCP) is the most populous plain in China and forms the core of the Beijing-Tianjin-Hebei economic circle. With urbanization, anthropogenic heavy metals have increasingly dispersed and accumulated in urban topsoil, especially in hazy cities. To investigate the major sources of haze and their relation to topsoil heavy metals concentrations in cities of the NCP, 220 topsoil samples (0–10 cm) were collected from Shijiazhuang city (capital of Hebei Province). The concentrations of eight selected metals were determined. Statistical and spatial distribution analyses suggest that coal combustion and industrial processes (Ni, Cr, Cd and Hg) were the dominant anthropogenic sources of haze in Shijiazhuang city, followed by vehicle exhausts (Pb, Zn Cu and Hg). Contrastingly, As was derived from parent materials of the NCP. A health risk assessment showed that Pb, Cr and As pose signiﬁcant non-carcinogenic risks to children (hazard index > 1) via oral ingestion. A potential carcinogenic risk to children (CRs > 10 − 4 ) is also posed by As. While Cd and Hg do not pose health risks in Shijiazhuang city, they may pose important ecological risks as ecological risk factors > 40 were observed, resulting in ecological risk indexes of 150–600 (moderate to considerable ecological risks).


Introduction
The North China Plain (NCP) is a typical Cenozoic sedimentary basin and one of the largest alluvial plains in eastern Asia. Two municipalities (Beijing and Tianjin) and five provinces (Hebei, Shandong, Henan, Anhui, and Jiangsu) are located within the NCP, which has a total area of 300,000 km 2 and is located between 36 • 00 -40 • 00 N and 114 • 00 -118 • 00 E [1]. The NCP is the most populous plain in China with a population of 339 million, accounting for 24.2% of China's total population in 2019. Due to rapid urbanization and industrialization, large cities in the NCP, such as Beijing, Tianjin, Shijiazhuang, Zhengzhou, Qingdao and Xuzhou are facing severe environmental problems due to soil pollution by heavy metals, which are persistent and toxic [2][3][4][5]. Heavy metal pollution in urban soils influences the ecological functionality of cities and poses strong risks to human health [6][7][8][9].
The sources of heavy metals in urban soils are widespread and complex [9,10]. Soil is a natural buffer that controls the transport of heavy metals to the atmosphere, hydrosphere and biosphere [11]. In addition, soil acts as a geochemical sink; once contaminated, heavy metals adsorbed on the fine particles of urban topsoil can be easily transferred into the human body via hand-to-mouth and dermal contact pathways [4,12]. In particular, in the case of outdoor activities for children, the process of resuspension by wind erosion or the In this study, we determined concentrations of eight heavy metals (Cu, Pb, Zn, Cr As, Cd, Ni, and Hg) in topsoil sampled from Shijiazhuang. The specific objectives were to: (1) assess the enrichment level of heavy metals in the topsoil; (2) determine the dominant causes of haze in relation to topsoil heavy metals; and (3) evaluate the health and ecological risks posed by topsoil heavy metals to inform better management practices.

Study Area
The NCP is characterized by flat terrain and numerous rivers and lakes. It has a convenient transportation system and a developed economy, and has developed into the political, economic and cultural centre of China. Shijiazhuang, the capital city of Hebei Province, is located in the hinterland of the NCP, bordering Beijing and Tianjin to the north, the Bohai Sea to the east and the Taihang Mountains to the west [28]. Geographically, Shijiazhuang is in the south-central part of Hebei Province within the NCP. It is mostly covered by cinnamon soils [29] and, from west to east, consists of piedmont alluvial, central alluvial-lacustrine, and littoral alluvial plains [12]. Shijiazhuang is also a large industrial city famous for raw material and energy production and its steel, power and cement industries [30]. In addition, Shijiazhuang is one of China's important commodity distribution centres, is a large northern commercial city, and is one of the main hosts of national trade conventions and exhibitions. It is also one of China's main railway transportation hubs, where four railways lines meet. In 2016, the number of vehicles was more than 2.0 million [30]. Shijiazhuang has obvious geographical advantages in terms of politics, economy and culture. The core urban area of Shijiazhuang is mainly a rectangular area (within one ring) bounded by Pingan Ave, Zhongshan Rd, Zhonghua Ave, Heping Rd, Tiyu Ave In this study, we determined concentrations of eight heavy metals (Cu, Pb, Zn, Cr As, Cd, Ni, and Hg) in topsoil sampled from Shijiazhuang. The specific objectives were to: (1) assess the enrichment level of heavy metals in the topsoil; (2) determine the dominant causes of haze in relation to topsoil heavy metals; and (3) evaluate the health and ecological risks posed by topsoil heavy metals to inform better management practices.

Study Area
The NCP is characterized by flat terrain and numerous rivers and lakes. It has a convenient transportation system and a developed economy, and has developed into the political, economic and cultural centre of China. Shijiazhuang, the capital city of Hebei Province, is located in the hinterland of the NCP, bordering Beijing and Tianjin to the north, the Bohai Sea to the east and the Taihang Mountains to the west [28]. Geographically, Shijiazhuang is in the south-central part of Hebei Province within the NCP. It is mostly covered by cinnamon soils [29] and, from west to east, consists of piedmont alluvial, central alluvial-lacustrine, and littoral alluvial plains [12]. Shijiazhuang is also a large industrial city famous for raw material and energy production and its steel, power and cement industries [30]. In addition, Shijiazhuang is one of China's important commodity distribution centres, is a large northern commercial city, and is one of the main hosts of national trade conventions and exhibitions. It is also one of China's main railway transportation hubs, where four railways lines meet. In 2016, the number of vehicles was more than 2.0 million [30]. Shijiazhuang has obvious geographical advantages in terms of politics, economy and culture. The core urban area of Shijiazhuang is mainly a rectangular area (within one ring) bounded by Pingan Ave, Zhongshan Rd, Zhonghua Ave, Heping Rd, Tiyu Ave and Yuhua Rd, and includes eight main roads; three horizontal and five vertical [31]. The region has high concentrations of workplaces, as well as social, economic, financial and trade activities. It is a typical centre of business and trade, economy and finance, culture and entertainment, and social and political activities. Shijiazhuang is characterized by a semi-arid continental monsoon climate with an annual mean temperature of 12-13 • C, annual precipitation of 400-800 mm and mean potential evaporation rate of 1100-1800 mm/year [24]. The average annual air humidity is 65% and about 63-70% of the annual precipitation fells from June to September.

Sample Collection and Chemical Analysis
The sampling methods in this study comply with the Geological Survey Technology Standard of the China Geological Survey. A total of 220 topsoil (0-10 cm) samples were collected from the central urban area of Shijiazhuang city and its surrounding suburban area according to the sampling density of 1 km 2 per sample ( Figure 1). During the sampling, we tried to avoid the interference of external soil and then removed impurities and gravel from the soil surface. Soil samples of about 1 kg were collected in a clean cloth bag with labels, and then dried naturally. The topsoil samples were gently tapped with a mallet before being sifted so that the soil samples remain in their natural size state. The sample was bottled with a nylon sieve with a grain size of 500 g less than 0.8 mm (20 mesh).
The Cr, Ni, Zn, Cu, and Pb concentrations were determined by X-ray fluorescence spectrometry (XRF) with detection limits of 2, 1, 2, 1, 1 mg/kg, respectively. The concentration of Cd in the samples were measured by ICP-MS (detection limits: 0.02 mg/kg) after the digestion (HCl-HNO 3 -HF-HClO 4 ). For the extraction of As and Hg metals, the topsoil samples were digested with aqua regia (HNO 3 : HCl = 1:3), and atomic fluorescence spectrometry (AFS) was used to determine their concentrations (with detection limits of 0.05 and 0.0003 mg/kg, respectively). The measure accuracy was verified using the Chinese standardized reference material (GSS-1). For quality control, reagent blanks, duplicate samples, reference standard and external monitor samples were inserted among the soil samples to evaluate the analysis methods. The recoveries were accepted within the range from 81.0 to 119%.

Method for Contamination Assessment
Contamination factor (CF) has been used to assess the pollution of heavy metals in topsoil. It was determined by the ratio of the concentration of heavy metal n for the sampled soil (C n ) to the values for the background (B n ). The background values of the NCP (average regional values) were adopted to calculate the CF values for the selected metals [32]. The value of CF was interpreted using four classes: low pollution (CF < 1), moderate pollution (1 ≤ CF < 3), high pollution (3 ≤ CF < 6) and very high pollution (CF ≥ 6) [7,33].

Method for Human Health Risk Assessment
Human health risk assessment is the process of estimating the probability of adverse health effects of hazard chemicals to human health according to human exposure characters (see details in Supplementary Materials). Hazard quotient (HQ) is typically used to evaluate the non-carcinogenic risk level of human exposure for one metal [4]. While, for multiple chemicals, the sum of the HQ values of all chemicals is expressed as a hazard index (HI) to assess the overall non-carcinogenic risk effects. Whereas for carcinogens the dose is multiplied by the corresponding slope factor (SF) to calculate carcinogenic risks (CR), to indicate the incremental probability of an individual developing cancer over a lifetime as a result of exposure to the potential carcinogen [3]. Generally, the values of HI > 1 and CR > 10 −4 , respectively, represent adverse effects to human health and lifetime carcinogenic risks [4]. The reference values and data source of the abbreviation for the parameters were shown in Tables S1 and S2.

Method for Ecological Risk Assessment
When considering the different toxicities of the heavy metals to human and biota, the potential ecological risk index (RI) was developed by Hakanson (1980) to assess the pollution degree and ecological risk from heavy metals in soil [31]. It takes in to account the CF, potential ecological risk factor (E r ) and the toxicological response factor (T r ) (i.e., Hg = 40, Cd = 30, As = 10, Pb = Cu = Ni = 5; Cr = 2; Zn = 1) of each heavy metal (i) [6,34] and is calculated by the following equation: The ecological risk levels can indicate low risk (E i r < 40, RI < 150), moderate risk (40 ≤ E i r < 80, 150 ≤ RI < 300), considerable risk (80 ≤ E i r < 160, 300 ≤ RI < 600), very high risk (160 ≤ E i r < 320, 600 ≤ RI), and dangerous (E i r > 320) [23,33].

Data Analysis
Descriptive statistics and data analyses were conducted with Microsoft Excel 2019 and IBM SPSS Statistics Version 21. Pearson's correlation coefficient analysis and principal component analysis were used to determine the possible sources of eight heavy metals in the urban topsoils. Mapping that shows spatial distributions of heavy metals was visualized using ArcGIS version 9.0 (ESRI, Redlands, CA, USA). Table 1 shows the statistical characteristics of the eight heavy metals and Al 2 O 3 (n = 220) in Shijiazhuang city. The K-S test results show that none of the eight metals had normally distributed concentration data (significance > 0.1). The average concentrations of metals of topsoil (in mg/kg) were compared and can be ranked as follows: Zn (105) > Cr (71.9) > Pb (31), Ni (28.2) and Cu (27.4) > As (9.42) > Cd (0.27) > Hg (0.113) ( Table 1). This ranking is similar to those of other urban topsoils around the world [5]. The coefficient of variation (CV) of the metals in Shijiazhuang city decreased in the order Cd > Hg > 100% > Zn > Cr > Pb > As > Ni > Cu > Al 2 O 3 , indicating high variations of Cd, Hg and Zn. The contamination factors (CFs) of the eight metals were in descending order of Hg > Cd > Zn > Pb > Cu > Cr > 1 > Ni > As (Table 1). Of the sampling sites, 85.0%, 93.2%, 87.3%, 83.6% and 70.9% were moderately contaminated by Cd, Zn, Pb, Cu and Cr, respectively, while 18.2%, 53.2% and 28.2% were moderately, highly and very highly contaminated by Hg, respectively. In addition, 97.7% of sites showed pollution load indexes (PLIs) greater than 1. The average contamination level was 1.67 (unpolluted to moderately polluted) with only 18.2% of observations being greater than the moderate pollution level (PLI > 2). Some 79.5% of sites were unpolluted to moderately polluted and only 2.3% were unpolluted. These results suggest that the metals in the topsoil of Shijiazhuang city differed greatly with respect to pollution source and concentration, especially for Hg, Zn, Pb, Cu, Cd and Cr.

Heavy Metal Accumulation in Topsoil
The economic development level and industrial structure vary among cities throughout China; therefore, the impacts of human activities on topsoil are diverse [4]. In general, the concentrations of heavy metals in the topsoil of different types of cities usually follow the order: industry-based cities (such as Baoji, Qingdao and Dongguan) > more developed cities (Hangzhou and Lishui) > metropoles (such as Beijing, Tianjin, Shanghai and Guangzhou) > western underdeveloped cities (Urumqi, Xining and Lhasa) ≈ county and agricultural areas (such as Xiangfen, Wulian and Xinglonggang) as shown in Table 2. The concentrations of most heavy metals (except As and Hg, but especially Zn and Cr) were higher in Shijiazhuang city than in county and agricultural areas and underdeveloped cities, and were generally lower than in some metropoles such as Shanghai and Tianjin. In such urbanization -leading metropoles, the potential emission intensity of anthropogenic heavy metal pollutants, including Pb, Cu and Zn in Shijiazhuang city were much weaker than those in industry-based and more developed cities with a large population.

Spatial Distribution of Heavy Metals
Spatially, the areas defined by Zhonghua Ave, Heping Rd, Tiyu Ave and Yuhua Rd ( Figure 1B) represent the core areas of Shijiazhuang city with the highest traffic flows, greatest workplace densities and highly concentrated social and economic activities. These areas usually showed the highest concentrations of Cu, Pb and Zn. The topsoil adjacent to railway stations and bus stations in northern Shijiazhuang also showed high concentrations of Cu, Pb and Zn ( Figures 1B and 2). In addition, the areas with high values corresponded to the industrial areas in the middle, east and northeast of Shijiazhuang ( Figure 1B). In contrast, high concentrations of Cr, Cd and Ni were distributed in the industrial area in the northeast of the urban area, which includes steel, coking, cement and coal machinery plants. Areas of high As concentrations occurred near the pharmaceutical plants in the north-central part of Shijiazhuang city, while, other areas had As concentrations below the background value. In contrast, there was a wide distribution of high-Hg-concentration areas suggesting that Hg pollution in the Shijiazhuang urban area is widespread and has multiple sources.

Source Identification
Correlation analysis (CA) and principal component analysis (PCA) were used to identify the sources of heavy metals in urban topsoil in Shijiazhuang city (Tables S3 and S4, respectively). A total of four principal components (PC; eigenvalues > 1) were extracted ( Figure 3).
For F1, significant positive loading values (>0.9) of Cr, Ni and Cd usually reflect contamination by atmospheric deposits in urban soils [35]. Coal-burning power plants and smelters could cause high Cr and Ni accumulation in nearby topsoil [5,36]. Previous studies show that high concentrations of Cr, Ni and Cu are discharged during metal electroplating and rolling processes, and have accumulated as dust in streets near the urban industrial area of Delhi city, India [37]. Yuswir et al. (2015) showed that Cd and Cr in urban topsoil may have originated from the electronics and semiconductor industries in Klang district, Malaysia [38]. Therefore, we infer that PC1 is primarily controlled by industrial activities and coal burning.

Source Identification
Correlation analysis (CA) and principal component analysis (PCA) were used to identify the sources of heavy metals in urban topsoil in Shijiazhuang city (Tables S3 and  S4, respectively). A total of four principal components (PC; eigenvalues > 1) were extracted ( Figure 3). For F1, significant positive loading values (>0.9) of Cr, Ni and Cd usually reflect contamination by atmospheric deposits in urban soils [35]. Coal-burning power plants and

Source Identification
Correlation analysis (CA) and principal component analysis (PCA) were used to identify the sources of heavy metals in urban topsoil in Shijiazhuang city (Tables S3 and  S4, respectively). A total of four principal components (PC; eigenvalues > 1) were extracted ( Figure 3). For F1, significant positive loading values (>0.9) of Cr, Ni and Cd usually reflect contamination by atmospheric deposits in urban soils [35]. Coal-burning power plants and Principal component F2 also had high loadings on Pb, Cu and Zn (0.7−0.9). Several studies have proven that vehicle emissions are the main source of Pb in urban topsoil [39][40][41]. Zn and Cu may come from mechanical wear, brake pads, oil leak pools and cylinder head washers on vehicles [42,43]. Engine and tire oil production requires Zn as an additive [17]. Cu, Pb and Zn pollution in urban topsoil where industrial factories and plants once operated remains an unresolved problem, especially in large cities of China [18,44]. Thus, F2 represents traffic sources coupled with industrial input.
For PC F3, only Al 2 O 3 and As were grouped in a positive correlation (r = 0.254, p < 0.01) (Tables S3 and S4). There were low concentrations of Al 2 O 3 and As (comparable to background levels) and low CF (<1) and CV% values (4.43% and 23.5%, respectively). Similar results have been reported in the eastern Hebei plain (Hebei Province) [45]. Therefore, F3 primarily indicates natural sources of heavy metal input. Principal component F4 was only loaded by Hg (Table S3). However, a 3D-loading plot shows that F4 was distributed in the space constrained by F1, F2 and F3 (Figure 3). This suggests that Hg has multiple-sources. Some 99.5% of the topsoil samples exceeded the background concentration of 0.02 mg/L. This can explain the relatively high average Hg concentrations in comparison to those of underdeveloped cities and agricultural based counties ( Table 2). Table 3 shows that the topsoil Cr concentration in Shijiazhuang was higher than that of Beijing (natural source) and Zhengzhou (non-industrial source). The Cd concentration in Shijiazhuang was higher than in Beijing (anthropogenic and industrial sources) and Tianjin (industrial source). This indicates that coal combustion and industrial processes are the dominant factors causing haze in Shijiazhuang. In contrast, the topsoil Pb concentration in Shijiazhuang was much lower than that in Zhengzhou, which seems to indicate that although there are much greater traffic volumes in central Shijiazhuang city, the influence of traffic emissions on haze is less than those of coal combustion and industrial processes.

Health Risk Assessment
The mean HQ and HI values of each metal were calculated to evaluate their noncarcinogenic risks. As shown in Figure 4, due to their high topsoil concentrations or low RfD values, As, Cr and Pb showed high HQ values (>1), which suggests that they present higher non-carcinogenic risks to humans than Cd, Cu, Zn and Hg (<0.1). Oral ingestion is the main exposure pathway for metals in Shijiazhuang, followed by dermal contact. It is noteworthy that children always had higher HI values for each metal than adults, demonstrating that children are at greater risk of metal exposure from topsoil than adults. For example, the HI values of Cr, Pb and As in topsoil were >1 for children. The HI value of Ni was >0.1, indicating a potential health hazard to children [47]. These results indicate that As, Cr, Pb and Ni in topsoil may pose non-carcinogenic risks to children in Shijiazhuang city, while there were acceptable HI values (<1) for adults.
In this study, three exposure pathways of carcinogenic risk (CR) from As and Cd were considered due to their relatively high concentrations in different types of cities (Table 2). However, only CRs due to inhalation were estimated for Cr and Ni [4]. The results show that oral ingestion was the main exposure pathway of As in both children and adults, because the CRs of~10 −4 were much greater than the acceptable limit of 10 −6 ( Figure 5). In addition, for children, Cd had a CR value slightly exceeding 10 −6 , while As had a CR value close to 10 −6 for the oral ingestion and dermal contact pathways. In contrast, the inhalation CRs of these four metals were much lower than 10 −6 , suggesting that the inhalation exposure pathway from topsoils in Shijiazhuang city has an acceptable risk.
Given the relatively high concentration of naturally sourced As in Shijiazhuang (11.0 mg/kg), mitigation measures should be taken during earth-moving projects to minimize the exposure of people, especially children. Children near transportation hubs and industrial areas should be protected against the non-carcinogenic risks of Cr, Pb and Ni and the carcinogenic risks of Cd through oral ingestion. It is worth paying additional attention to reducing the emissions of industries and enterprises near urban residential areas. In this study, three exposure pathways of carcinogenic risk (CR) from As and Cd were considered due to their relatively high concentrations in different types of cities (Table 2). However, only CRs due to inhalation were estimated for Cr and Ni [4]. The results show that oral ingestion was the main exposure pathway of As in both children and adults, because the CRs of ~10 −4 were much greater than the acceptable limit of 10 −6 ( Figure  5). In addition, for children, Cd had a CR value slightly exceeding 10 −6 , while As had a CR value close to 10 −6 for the oral ingestion and dermal contact pathways. In contrast, the inhalation CRs of these four metals were much lower than 10 −6 , suggesting that the inhalation exposure pathway from topsoils in Shijiazhuang city has an acceptable risk.
Given the relatively high concentration of naturally sourced As in Shijiazhuang (11.0 mg/kg), mitigation measures should be taken during earth-moving projects to minimize the exposure of people, especially children. Children near transportation hubs and industrial areas should be protected against the non-carcinogenic risks of Cr, Pb and Ni and the carcinogenic risks of Cd through oral ingestion. It is worth paying additional attention to reducing the emissions of industries and enterprises near urban residential areas.

Potential Ecological Risks
The accumulation of heavy metals in urban topsoil influences urban ecosystems chronically through direct exposure or indirect trophic transfer [34]. Unlike Ni and As,

Potential Ecological Risks
The accumulation of heavy metals in urban topsoil influences urban ecosystems chronically through direct exposure or indirect trophic transfer [34]. Unlike Ni and As, the metals Cr, Cd, Pb, Cu and Zn were enriched in the topsoil of Shijiazhuang city to various degrees (69.5% of the topsoil had concentrations greater than the background). In general, the soil suffered moderate pollution (70.9-93.2% of topsoil with a CF between 1-3 and 79.5% of topsoil with a PLI between 1-2). However, for Cd and Hg in topsoil, 87.7% and 100% of the samples had Er values above 40, while other heavy metals had Er values below 40 ( Figure 6). As a result, Shijiazhuang soil is faced with potential ecological risks that are moderate (52.7% of topsoil with RI of 150-300) to considerable (35.5% of topsoil with RI = 300-600).  (Figure 6). As a result, Shijiazhuang soil is faced with potential ecological risks that are moderate (52.7% of topsoil with RI of 150-300) to considerable (35.5% of topsoil with RI = 300-600).
Mapping can provide useful information for pollution risk assessment (Figure 7). Heavy metals with a PLI > 1 were observed in 97.7% of topsoil samples, while 94.5% had RI > 150. The highest PLI and RI values were found in the east near the industrial area (with steel, coking, cement and coal machinery plants) in the northeast of the urban area. It is worth noting that the spatial distribution of RI in the topsoil was similar to that of Hg concentrations, while PLI was similar to those of Cr and Cd. Although Cd and Hg do not pose health risks to human, they do pose important ecological risks to the urban area of Shijiazhuang. Figure 6. Contamination factor (CF) and potential ecological risk factor (Er) for heavy metals in topsoil from Shijiazhuang city. Figure 6. Contamination factor (CF) and potential ecological risk factor (Er) for heavy metals in topsoil from Shijiazhuang city.
Mapping can provide useful information for pollution risk assessment (Figure 7). Heavy metals with a PLI > 1 were observed in 97.7% of topsoil samples, while 94.5% had RI > 150. The highest PLI and RI values were found in the east near the industrial area (with steel, coking, cement and coal machinery plants) in the northeast of the urban area. It is worth noting that the spatial distribution of RI in the topsoil was similar to that of Hg concentrations, while PLI was similar to those of Cr and Cd. Although Cd and Hg do not pose health risks to human, they do pose important ecological risks to the urban area of Shijiazhuang.

Conclusions
This study analyzed the sources of heavy metals in topsoil in Shijiazhuang and the associated health and ecological risks. The CA and PCA methods indicated that the major pollution sources were industrial processes, coal-burning emissions and vehicle exhausts. The spatial distributions of traffic and industry have an important influence on the enrichment of anthropogenically sourced heavy metals. High Ni, Cr and Cd levels were concentrated near industrial districts but were also related to coal combustion activities, while high concentrations of Pb, Zn and Cu were mainly attributed to urban transport in central Shijiazhuang city and industrial production around the city. Hg pollution was widespread and showed multiple sources. Coal combustion and industrial production may be the most important factors causing haze in Shijiazhuang.
The health risk assessment indicated that Cr, Cd and Pb concentrations in urban topsoil are generally acceptable, while that of As is not. Arsenic may pose high carcinogenic (CR > 10 −4 ) and non-carcinogenic risks (HI > 1) to children through oral ingestion, even though the source of As was mainly parent materials. This shows that contact between local residents and soil As should be avoided during urban development. Therefore, special attention should be paid to As contamination in topsoil and its potential health risks.

Conclusions
This study analyzed the sources of heavy metals in topsoil in Shijiazhuang and the associated health and ecological risks. The CA and PCA methods indicated that the major pollution sources were industrial processes, coal-burning emissions and vehicle exhausts. The spatial distributions of traffic and industry have an important influence on the enrichment of anthropogenically sourced heavy metals. High Ni, Cr and Cd levels were concentrated near industrial districts but were also related to coal combustion activities, while high concentrations of Pb, Zn and Cu were mainly attributed to urban transport in central Shijiazhuang city and industrial production around the city. Hg pollution was widespread and showed multiple sources. Coal combustion and industrial production may be the most important factors causing haze in Shijiazhuang.
The health risk assessment indicated that Cr, Cd and Pb concentrations in urban topsoil are generally acceptable, while that of As is not. Arsenic may pose high carcinogenic (CR > 10 −4 ) and non-carcinogenic risks (HI > 1) to children through oral ingestion, even though the source of As was mainly parent materials. This shows that contact between local residents and soil As should be avoided during urban development. Therefore, special attention should be paid to As contamination in topsoil and its potential health risks. In addition, Cr and Pb may pose non-carcinogenic risks to children, while Cd may pose carcinogenic risks. Hg can pose significant ecological risks to residents of Shijiazhuang city (35.5% of topsoil had RI = 300-600), even though it is considered to be safe for human health (HI < 1). Studies related to heavy metals in urban topsoil should focus on potential health and ecological risks, rather than simply conducting source analysis or singular risk assessments.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/su131810046/s1, Table S1: Definition and reference value of some parameters for health risk assessment of heavy metal, Table S2: Summary of reference doses (RfD) and slope factors (SF) of heavy metals, Table S3: Rotated component matrix for heavy metals data of urban soils from Shijiazhuang city, Table S4: Correlation matrix for the heavy metals of urban soil in Shijiazhuang city.