Heavy Metal Pollution and Health Risk Assessment of Vegetable–Soil Systems of Facilities Irrigated with Wastewater in Northern China

Soil pollution by heavy metals is a major concern in China and has received much attention in recent years. Aiming to investigate the status of heavy metal pollution and the safety of vegetables in the soil of wastewater-irrigated facilities, this study investigated the distribution and migration characteristics of heavy metals in vegetable–soil systems of facilities in a typical sewage irrigation area of the Xi River, Shenyang City, northern China. Health risks due to the fact of exposure to heavy metals in the vegetable soil of facilities and ingrown vegetables through different exposure pathways were evaluated. Spatial interpolation and a potential ecological risk assessment were applied to evaluate the soil quality. Bioaccumulation factors (BCFs) were used to analyze the absorption and transportation capacity of Cd, Cu, Pb, and Zn by different parts of different vegetables. The results showed that the average concentration of Cd exceeded the standard values by 1.82 times and accumulated by 11 times, suggesting that Cd poses the most severe pollution among the four metals in the soil of facilities in the Xi River sewage irrigation area. In the city, a significant accumulation of Cd in the soil was identified with different spatial distributions. Cd also contributed the most in terms of the estimated potential ecological risk index, while the impacts of the other three metals were relatively small. The concentrations of heavy metals were mostly lower than the limit set by the corresponding Chinese standards. Various BCFs were observed for the four metals in the order Cd > Zn > Cu > Pb. Vegetables also demonstrated different BCFs in the order of leaf vegetables > Rhizome vegetable > Solanaceae vegetable. The magnitude of the noncarcinogenic risk for all four heavy metals was less than one for all three exposure routes and did not cause significant noncarcinogenic health effects in humans. However, the carcinogenic risk of Cd from some vegetables via dietary intake was considered higher. Protection measures should be taken to implement better pollution control and land use planning.


Introduction
China is the largest vegetable producing country in the world, accounting for 36% of the world's vegetable planting area and 66% of the world's vegetable production [1]. With China's economic and social development, the standard of quality for vegetables has increased. Hence, heavy metal pollution in vegetables has drawn great attention nationwide. Rapid development of facilities in the vegetable industry has occurred, and research on the accumulation of heavy metals in these lands has gradually increased. Several studies have suggested that the decline in soil quality was caused by the increase in the level of heavy metals, which was the result of the excessive usage of vegetable land, fertilizer, and a high cropping index [2][3][4][5]. Studies by several groups have also reported the 2 of 15 transfer of pollutants from the soil into the food chain, endangering residents' health [6][7][8][9]. The issue of declining soil quality is related to metal pollution in sewage irrigation areas, which has attracted many researchers' attention due to the accumulation of heavy metals, which is hazardous to food safety and human health [10][11][12][13][14][15]. Hence, it is necessary to perform investigations on heavy metal pollution, as monitoring, prevention, and controls are needed to ensure the quality and safety of the vegetables.
Differences in bioaccumulation ability can be found among different crops. In general, the abilities of absorption, conversion, and bioaccumulation of heavy metals in vegetables are higher than that in traditional plant crops. Some vegetables can not only absorb some heavy metals but also have special organs with the bioaccumulation ability to accumulate pollutants including As in the roots of carrot, Hg in bean pods, Pb and Cd in the roots of radish, and Se in the leaves of radish [16]. The concentration of heavy metals in soil generally increases over time in sewage irrigation areas. All of this results in the risk of the bioaccumulation of heavy metals by the vegetables planted in sewage irrigation areas. Due to the fact of accelerated urbanization and continuous increase in planting area and the production of vegetables in northern China, the issue regarding heavy metals in the soil in these areas has become prominent [1,17].
The Xi River sewage irrigation area in Shenyang, northern China, is an atypical irrigation area that was irrigated by industrial wastewater from the 1950s to the 1990s. It is one of the areas irrigated for the longest continuous period, with the most serious pollution in China [18][19][20]. Although there has been no continuous sewage irrigation in that area for more than twenty years, and the local government has implemented policies, including adjusting the planting crops, the soil and crops are still affected by the sustained effects. The potential ecological risks from the irrigation of accumulated sewage warrant continuous monitoring and investigations, as heavy metals are difficult to degrade and remove from the soil. In addition, whether the risk to human health has been altered due to the fact of changes in people's diet is unknown. Meanwhile, researchers have found that despite the heavy pollution, the Xi River water is still used as a source of agricultural irrigation due to the lack of clean water [20]. Xi River water has greatly accumulated heavy metals, as it has been polluted by sewage from industrial and agricultural use in addition to the surrounding cities. Different from a previous study on the farmland soil of the Xi River sewage irrigation area in Shenyang [21], this study took facilities' vegetable soil as the research object. There is a certain difference between the two farming and irrigation methods. The accumulation and spatial distribution characteristics of heavy metals in soil were investigated to evaluate the risks to the local communities' health [22][23][24]. Aimed at investigating the heavy metal pollution status and vegetable safety in the soil of wastewater-irrigated facilities, this study investigated the distribution and migration characteristics of heavy metals in a facility's vegetable soil-vegetable system in a typical sewage irrigation area of the Xi River, Shenyang city, northern China. Meanwhile, health risks due to the fact of exposure to the heavy metals in the facility's vegetable soil and the grown vegetables through different exposure pathways were evaluated to provide scientific data to support the adjustment of the spatial arrangement for vegetable food safety, the agricultural planting approach, and the land-use planning in that region.

Sample Collection and Analysis
The Xi River sewage irrigation area is located in west Shenyang (41 • 36 ~41 • 45 N, 123 • 02 ~123 • 18 E), and it is a typical sewage irrigation area in northern China. Wastewater and domestic sewage in the industrial zone of west Shenyang were used as irrigation in this area after being diluted with a small amount of Hun River water. As a tributary of the Hun River, the Xi River originated in the Tiexi District of Shenyang city, with a length of approximately 80 km and an average width of 8 m [19]. In this study, the sampling points were distributed along the river basin ( Figure 1). irrigation in this area after being diluted with a small amount of Hun River water. As a tributary of the Hun River, the Xi River originated in the Tiexi District of Shenyang city, with a length of approximately 80 km and an average width of 8 m [19]. In this study, the sampling points were distributed along the river basin ( Figure 1). In April 2019, 0-20 cm layer of arable soil and vegetable samples were collected. A total of 12 sampling points in the greenhouse area of typical facilities were set-up, and each sample consisted of five equiponderance subsamples taken from the surface of five sites. The GPS locations for the sampling sites are shown in Figure 1. After the soil samples were evenly mixed, each one was reduced to 100 g by the quartering method. Then, every sample was dried in an oven, with the tiny stones, animals, and plant residues in the samples removed. All samples were evenly mixed again after being ground using an agate mortar and passed a 0.149 mm nylon sieve (100 mesh). There were 35 vegetable samples, including 8 leafy vegetables (i.e., 1 cabbage, 4 coriander, 1 chlorella, 9 lettuce, 3 Chinese lettuce, 3 asparagus leaves, 2 celery, and 2 spinach),1 rootstock vegetable (i.e., 7 radishes), and 1 solanaceous vegetable (i.e., 3 tomatoes), collected from the study area. Vegetable samples were washed with tap water and purified water to remove the dirt and soil grains and then sterilized at 105 °C and dried at 65 °C. Fresh and dry weights were recorded before and after the drying process. Dried samples were further ground to pass through a2 mm sieve and stored in plastic bags at room temperature until analysis.
Prior to the measurement of the metal contents, the soil samples were digested with HNO3-HF-HClO4, and the vegetable samples were boiled with HNO3. Inductively coupled plasma was used to measure the concentrations of Cd, Cu, Pb, and Zn in the soil and vegetable samples. Parallel samples and blank samples were included for quality control purposes to evaluate the reliability and accuracy of the measurements.

Heavy Metal Ecological Risk Index
The Potential Ecological Risk Index method was developed by the Swedish scientist Hakanson to quantitatively evaluate the extent of heavy metal contamination in sediments [25]. The method takes into account the heavy metal content, ecological effects and environmental toxicological characteristics at the same time and is calculated by the formula: In April 2019, 0-20 cm layer of arable soil and vegetable samples were collected. A total of 12 sampling points in the greenhouse area of typical facilities were set-up, and each sample consisted of five equiponderance subsamples taken from the surface of five sites. The GPS locations for the sampling sites are shown in Figure 1. After the soil samples were evenly mixed, each one was reduced to 100 g by the quartering method. Then, every sample was dried in an oven, with the tiny stones, animals, and plant residues in the samples removed. All samples were evenly mixed again after being ground using an agate mortar and passed a 0.149 mm nylon sieve (100 mesh). There were 35 vegetable samples, including 8 leafy vegetables (i.e., 1 cabbage, 4 coriander, 1 chlorella, 9 lettuce, 3 Chinese lettuce, 3 asparagus leaves, 2 celery, and 2 spinach),1 rootstock vegetable (i.e., 7 radishes), and 1 solanaceous vegetable (i.e., 3 tomatoes), collected from the study area. Vegetable samples were washed with tap water and purified water to remove the dirt and soil grains and then sterilized at 105 • C and dried at 65 • C. Fresh and dry weights were recorded before and after the drying process. Dried samples were further ground to pass through a2 mm sieve and stored in plastic bags at room temperature until analysis.
Prior to the measurement of the metal contents, the soil samples were digested with HNO 3 -HF-HClO 4 , and the vegetable samples were boiled with HNO 3 . Inductively coupled plasma was used to measure the concentrations of Cd, Cu, Pb, and Zn in the soil and vegetable samples. Parallel samples and blank samples were included for quality control purposes to evaluate the reliability and accuracy of the measurements.

Heavy Metal Ecological Risk Index
The Potential Ecological Risk Index method was developed by the Swedish scientist Hakanson to quantitatively evaluate the extent of heavy metal contamination in sediments [25]. The method takes into account the heavy metal content, ecological effects and environmental toxicological characteristics at the same time and is calculated by the formula: where i represents the names of heavy metals; RI refers to comprehensive potential ecological risk factor; Er i is the potential ecological risk factor of substance i; Tr i is the toxic reaction factor of element i in the reference value from the toxic factors of Hakanson's research (the values for Cd, Cu, Pb, and Zn are 30, 5, 5, and 1, respectively) [25]; C i f is the pollution factor of a single substance i; C i is the tested concentration of substance i in mg·kg −1 ; m is the number of heavy metals; C i r is the reference value of element i for comparison, where the background values of heavy metals in Liaoning Province were used. The results were grouped into five classes, namely, low risk (Er i < 40, RI < 150), moderate risk (40 ≤ Er i < 80, 150 ≤ RI < 300), considerable risk (80 ≤ Er i < 160, 300 ≤ RI < 600), high risk (160 ≤ Er i < 320, RI ≥ 600), and very high risk (Er i ≥ 320) [25].

Vegetables' Bioaccumulation Index for Heavy Metals
The vegetables' bioaccumulation index for heavy metals was used to measure the ability of vegetables to absorb and accumulate heavy metals [26]. BCFs were calculated using the following equation: where C vegetable is the concentration of heavy metals in the edible parts of the vegetables (dry weight, mg·kg −1 ); C soil is the total concentration of heavy metals in the soil.

Health Risk Assessment
There are four human exposure pathways for heavy metals in soil. The first one is dietary ingestion through the food chain, where heavy metals are transferred into human bodies through fruits, vegetables, and crops grown in polluted soil. The second one is the direct ingestion of soil grains. The third one is the inhalation of polluted soil dust via the mouth or nose. The fourth one is dermal intake, where skin directly contacts polluted soil particles. Formulas for the calculation of the health risk, proposed by the US Environmental Protection Agency in combination with other studies [3,21,27], were used to evaluate the risk levels of heavy metals in the soil through dietary ingestion. Health risks can be divided into carcinogenic and noncarcinogenic risks. Noncarcinogenic elements in this study included Cd, Pb, Cu, and Zn, with Cd also being a carcinogenic element [28]. Cd is classified as a confirmed carcinogen for human lung and prostate cancers by the International Agency for Research on Cancer, and it can also cause multiorgan damage, mainly in kidney and bone, cardiovascular and cerebrovascular diseases, and metabolic diseases [29][30][31].
The daily exposure doses of contaminants via the various exposure pathways were calculated using the following equations: where ADI is the average daily intake of heavy metals in vegetables (mg·kg −1 ·day −1 ); C is the tested concentration of element i in the sampling area (mg·kg −1 ). All parameters and their definition are summarized in Table 1.  (7) where HQ is the ratio of exposure to a single noncarcinogenic substance; HI is the total risk index to the noncarcinogenic composition of multiple elements; RFD is the noncarcinogenic reference dose of a contaminant (mg·kg −1 ·day −1 where CRI is the probability of health risks from a carcinogen to an individual; SF is the carcinogenic slope factor of a carcinogenic contaminant (mg·kg −1 ·day −1 ). We only included the risk of cancer from dietary intake (SF = 6.1) and inhalation (SF = 6.4) in this study. RT is the cumulative cancer risk of the composition of multiple heavy metals. The only carcinogenic element was Cd in this study.

Data Analysis
Descriptive statistics were calculated using SPSS 19. Origin 8 was employed for the graph plotting. Spatial distribution mapping of heavy metal concentrations using ordinary Kriging was performed with ArcGIS 10.6 (ESRI Co., Redlands, CA, USA).

Concentrations of Heavy Metals in Soil
The concentrations of heavy metals, in general, exceeded the standards in the greenhouse soil of the Zhangshi sewage irrigation area ( The heavy metal cadmium (Cd) was found to be one of the important pollutants in the Xi River sewage irrigation area in this study, which is consistent with previous studies [18,21]. The coefficient of variation (CV) is regarded as an important indicator of the degree of interference of soil heavy metal content distribution by human activities. The degree of variability can be classified as weak variability (CV < 10%), moderate variability (10% ≤ CV < 100%), and strong variability (CV > 100%) [38,39]. The coefficients of variation (CV) for Cd and Cu were 40.84% and 46.65%, respectively, and they were highly variable (CV > 36%); the CV of Pb and Zn were 20.59% and 27.00%, respectively, and they were moderately variable (15% < CV < 36%). These data show that the concentrations of cadmium and copper varied between samples from different locations, while the distributions of Zn and Pb were relatively even with smaller differences. Compared to existing studies over the last 20 years [40][41][42], there was a significant decline in the levels of Cu and Zn in the soil, while no decrease in Cd and Pb in the soil was observed, even after an extended period of recovery from sewage irrigation. In comparison with similar sewage irrigation areas in other cities (i.e., Tianjin and Beijing), the contents of Cd and Pb in the sewage irrigation soil in Shenyang were higher than that in Tianjin [43]. In addition, the content of Cu was higher than that in the sewage irrigation areas in Tongzhou District, Beijing [44]. The difference in heavy metal contents between these sewage irrigation areas may be attributed to the distinctive sources of sewage, irrigation frequency, and soil background values. In the studied area, there were many manufacturing facilities, including metal casting, nonferrous metal, chemical, and leather-making industries, which were responsible for the discharge of heavy metal pollutants into the water and soil. These local industries are likely the determinant of the type and level of heavy metal pollutants in the study area. Although previous reports suggest that heavy metal pollution may also be affected by atmospheric sedimentation and traffic [18,45], the actual impact of these factors on the pollution in the study area requires further investigation.

Spatial Distribution of the Heavy Metals in the Soil
The spatial distribution of the sampling points in the study area is shown in Figure 2. Cd, Cu, Pb, and Zn demonstrated different spatial distribution patterns. The overall contamination level of Cd was moderate, with the maximum values mainly distributed in the northwest part, decreasing from the northwest to the southwest along the river. The locations with the maximum levels of Cu were relatively scattered, with higher values in the central and west parts and the lower values in the north. Moreover, the concentration of Cu demonstrated a large variation, ranging from15.23 to 88.58 mg·kg −1 . The spatial distribution of Pb was also scattered, but the pattern was opposite, with the high-value 7 of 15 points mainly distributed in the northwest and central parts. The spatial distribution of Zn was similar to Cd, decreasing from the southeast to the northwest. The maximum value area for Zn was located in ZS2, with the minimum value in the north. In general, the concentration distribution of Zn decreased from the ZS2 point to other areas, where the concentration of Zn was mainly 150 mg·kg −1 . Lian et al. [21] compared the concentration of heavy metals in the upstream, middle, and downstream of the Xi River sewage irrigation area and found a declining trend in the concentrations of Cd, Pb, and Zn along the river. The high level of Cd pollution in the upstream was near the Tiexi District, which was the main receiver of industrial and domestic sewage. The major sources of metal contamination in the soil were sewage irrigation and domestic and solid wastes. Although the middle and downstream of the Xi River were far away from the main pollution source, pollutant concentrations were still relatively high. This phenomenon suggests that self-purification capacity could not reduce the concentration of metals to a safe level. Finally, the use of fertilizer and livestock manure may have also led to the accumulation of heavy metals in the soil due to the fact that the middle area of the Xi River was the main agricultural area.
The locations with the maximum levels of Cu were relatively scattered, with higher va ues in the central and west parts and the lower values in the north. Moreover, the con centration of Cu demonstrated a large variation, ranging from15.23 to 88.58 mg·kg −1 . Th spatial distribution of Pb was also scattered, but the pattern was opposite, with th high-value points mainly distributed in the northwest and central parts. The spatial dis tribution of Zn was similar to Cd, decreasing from the southeast to the northwest. Th maximum value area for Zn was located in ZS2, with the minimum value in the north. I general, the concentration distribution of Zn decreased from the ZS2 point to other area where the concentration of Zn was mainly 150 mg·kg −1 . Lian et al. [21] compared th concentration of heavy metals in the upstream, middle, and downstream of the Xi Rive sewage irrigation area and found a declining trend in the concentrations of Cd, Pb, an Zn along the river. The high level of Cd pollution in the upstream was near the Tiex District, which was the main receiver of industrial and domestic sewage. The majo sources of metal contamination in the soil were sewage irrigation and domestic and soli wastes. Although the middle and downstream of the Xi River were far away from th main pollution source, pollutant concentrations were still relatively high. This phenom enon suggests that self-purification capacity could not reduce the concentration of metal to a safe level. Finally, the use of fertilizer and livestock manure may have also led to th accumulation of heavy metals in the soil due to the fact that the middle area of the X River was the main agricultural area.

Ecological Risk of Soil Heavy Metal Pollution
The potential ecological risk of heavy metals in soil irrigated by the Xi River sewag is shown in Figure 3. The potential ecological risk factors for the four heavy metals in th soil samples demonstrated a descending order of Cd > Cu > Pb > Zn, and the potentia

Ecological Risk of Soil Heavy Metal Pollution
The potential ecological risk of heavy metals in soil irrigated by the Xi River sewage is shown in Figure 3. The potential ecological risk factors for the four heavy metals in the soil samples demonstrated a descending order of Cd > Cu > Pb > Zn, and the potential risk factors for Cu, Pb, and Zn were all lower than the mild risk level [25]. In addition, the factors for Cd were higher than the intense risk level. Secondly, among the four metals, the contribution of Cd was the highest, and the combined potential ecological risk level varied among each sampling point. Since the study area mostly has a high ecological risk level, the resulting ecological risk issue cannot be ignored. Overall, our results showed that Cd contributed the most to the potential ecological risk, while Cu, Pb, and Zn contributed much less. factors for Cd were higher than the intense risk level. Secondly, among the four metals, the contribution of Cd was the highest, and the combined potential ecological risk level varied among each sampling point. Since the study area mostly has a high ecological risk level, the resulting ecological risk issue cannot be ignored. Overall, our results showed that Cd contributed the most to the potential ecological risk, while Cu, Pb, and Zn contributed much less.

Concentrations of Heavy Metals in Vegetables
The collected vegetables were eight leaf vegetables (i.e., cabbage, coriander, chlorella, lettuce, Chinese lettuce, asparagus leaves, celery, and spinach), one Rhizome vegetable (i.e., radish), and one Solanaceae vegetable (i.e., tomato). The moisture content of all vegetables was above 90%. The dry weights of the heavy metals in the vegetables were converted into fresh weights and compared with the safe concentration limit for Cu in "Tolerance limit of copper in foods (GB 15199-1994)" [46], the safe concentration limit for Zn in "Tolerance limit of copper in foods (GB 13106-1991)" [47], and the safe concentration limits for Cd and Pb in "Standard of tolerance limit of pollution in food security (GB 2762-2017)" [48], as shown in Table 3.

Concentrations of Heavy Metals in Vegetables
The collected vegetables were eight leaf vegetables (i.e., cabbage, coriander, chlorella, lettuce, Chinese lettuce, asparagus leaves, celery, and spinach), one Rhizome vegetable (i.e., radish), and one Solanaceae vegetable (i.e., tomato). The moisture content of all vegetables was above 90%. The dry weights of the heavy metals in the vegetables were converted into fresh weights and compared with the safe concentration limit for Cu in "Tolerance limit of copper in foods (GB 15199-1994)" [46], the safe concentration limit for Zn in "Tolerance limit of copper in foods (GB 13106-1991)" [47], and the safe concentration limits for Cd and Pb in "Standard of tolerance limit of pollution in food security (GB 2762-2017)" [48], as shown in Table 3. The concentrations of Cd, Cu, and Zn were lower than the national standard, while the level of Pb in only one sample exceeded the national standard. The concentrations of metals in Rhizome and Solanaceae vegetables were all lower than the standard values. The concentrations of Cd, Cu, and Zn in all vegetables were lower than the standard values. Except for leaf vegetables, the average concentrations of Pb in other vegetables were lower than the standard values. Considering that most vegetables are consumed in their fresh form, the dry weights of heavy metals in the vegetable samples were converted to fresh weights by the moisture coefficient.

Accumulation of Heavy Metals
To further investigate the accumulation of heavy metals in the soil-vegetable system, the bioaccumulation factors (BCFs) for Cd, Cu, Pb, and Zn were calculated ( Table 4). The average bioaccumulation factors for Cd, Cu, Pb, and Zn were in the order of Cd (0.4969) > Zn (0.1532) > Cu (0.0529) > Pb (0.0056) on the ground and Cd (0.4937) > Zn (0.1547) > Cu (0.0768) > Pb (0.0219) underground. The BCFs for Cd were higher than that of the other metals, with higher migration rates in the soil crops. The order of the bioconcentration factors (BCFs) of the heavy metals was Cd > Zn > Cu > Pb in the vegetable parts on the ground and the parts underground. In addition, the BCFs of the four metals in different vegetables varied and demonstrated the order: leaf vegetables > root vegetable > solanaceous vegetable, which is similar to the conclusions from studies of the sewage irrigation areas in Tianjin, China [43,49]. Among the vegetables, the BCFs of leaf vegetables were the highest, indicating that this type of vegetable has a strong ability to accumulate metals in soil. This finding is similar to a report by Li et al. [50] in which leaf vegetables exhibited higher accumulation levels of heavy metals than non-leaf vegetables.

Health Risk Assessment
The estimated daily exposure doses are presented in Table 5, i.e., the heavy metals in males, females, and children by the four exposure pathways (i.e., direct ingestion of soil grains, inhalation, dermal contact, and dietary ingestion). Table 5 also demonstrates the differences in exposure doses for the four metals in different groups. The total daily exposure doses were in the following order: Zn > Cu > Pb > Cd for all groups. Meanwhile, by comparing the four exposure pathways, dietary ingestion was the dominant route for the total exposure doses of heavy metals (Table 5). Our data also suggest that the doses of hazardous metals in these four pathways for children are higher than that for male adults, similar to findings by Li et al. [50], Zeng et al. [51], and Chen et al. [23].The exposure doses of Cd, Pb, and Zn for female adults were higher than that of children in this study, in contrast to the results of Lian et al. [21]. Overall, our results suggest the necessity of considering the differences in exposure among genders and ages. Since no published reference dose (RFD)for the direct ingestion of soil grains was available, this study only included the noncarcinogenic risk index of dietary intake, inhalation ingestion, and dermal contact. The noncarcinogenic risks caused by heavy metals to children and adults in the soil and vegetable samples are presented in Figure 4. The HQ of heavy metal elements Pb and Zn demonstrated a descending order of dietary ingestion > inhalation ingestion > dermal contact, while the HQ for Cd and Cu demonstrated a descending order of dietary ingestion> dermal contact > inhalation ingestion. Similar to the daily ingestion dose, the daily ingestion dose estimated by HQs was 5-7 magnitudes higher than that of nondietary ingestion. Comparing the exposure risk by dietary ingestion of the four heavy metals, the ingestion exposure risk of Cd was the highest, which was the dominant contaminant. The HQ for Cd for male adults, female adults, and children was 0.117, 0.135, and 0.208, respectively. When comparing the exposure risks for the three routes of intake for the four elements, the highest dietary exposure risk was found for cadmium. On the one hand, cadmium was the main contaminant of noncarcinogenic risk and, on the other hand, dietary intake was the main route. Although the overall noncarcinogenic risk (HI) for all four heavy metals was less than 1, there may still be some potential health risks to the local population from dietary exposure to cadmium. According to the HI values for each heavy metal, vegetable consumption by adults and children would not cause any obvious hazard to health. In comparison, Lian et al. [21] found that Cd was the most important contaminant, and its HQ for male adults, female adults, and children was higher than the threshold value of 1, indicating harmful effects to the health of residents.
At present, the slope factor (SF) for the carcinogenic risk of soil intake and skin contact of Cd has not been published, which led to this study only including dietary intake and respiratory intake in the carcinogenic risk assessment, as shown in Figure 5. For male and female adults and children, the carcinogenic risk (CRI) of Cd by inhalation ingestion was 2.90 × 10 −10 , 3.35 × 10 −10 , and 1.30 × 10 −10 , respectively. The carcinogenic risk of Cd by dietary ingestion for them was 3.06 × 10 −4 , 3.53 × 10 −4 , and 1.09 × 10 −4 , respectively, which was 1.09-3.06 times higher than the suggested value of 10 −6 -10 −4 by the US Environmental Protection Agency [27], and higher than the suggested value of 5 × 10 −5 by the International Commission on Radiological Project (ICRP). Similar to the noncarcinogenic risk, dietary ingestion was the dominant pathway causing the estimated carcinogenic risk. In this study, dietary ingestion was the main pathway for heavy metal exposure, and the potential risk to health for all ages would increase once heavy metal pollutants were transferred into daily edible vegetables in the food chain. Further studies are needed to improve the accuracy of the model by focusing on the heavy metal exposure by dietary ingestion, to provide more detailed data on the health risk for the local community, and to provide a scientific and theoretical basis for local governments to consider adjusting agricultural production activities in the area. All this information can support the government to adjust agricultural production activities in the polluted area. The exposure risk of heavy metals to children is higher than that of adults due to the fact of several characteristics of children including eating, drinking, breathing, weight parameters lower than that of adults, and relatively weak metabolism and detoxification abilities. Hence, children are more easily affected by environmental hazards than adults [51][52][53]. As a result, the same doses of heavy metals may have negative effects on children's development while being nonhazardous to adult health. adults, and children was 0.117, 0.135, and 0.208, respectively. When comparing the exposure risks for the three routes of intake for the four elements, the highest dietary exposure risk was found for cadmium. On the one hand, cadmium was the main contaminant of noncarcinogenic risk and, on the other hand, dietary intake was the main route. Although the overall noncarcinogenic risk (HI) for all four heavy metals was less than 1, there may still be some potential health risks to the local population from dietary exposure to cadmium. According to the HI values for each heavy metal, vegetable consumption by adults and children would not cause any obvious hazard to health. In comparison, Lian et al. [21] found that Cd was the most important contaminant, and its HQ for male adults, female adults, and children was higher than the threshold value of 1, indicating harmful effects to the health of residents. At present, the slope factor (SF) for the carcinogenic risk of soil intake and skin contact of Cd has not been published, which led to this study only including dietary intake and respiratory intake in the carcinogenic risk assessment, as shown in Figure 5. For male and female adults and children, the carcinogenic risk (CRI) of Cd by inhalation ingestion was 2.90 × 10 −10 , 3.35 × 10 −10 , and 1.30 × 10 −10 , respectively. The carcinogenic risk of Cd by dietary ingestion for them was 3.06 × 10 −4 , 3.53 × 10 −4 , and 1.09 × 10 −4 , respectively, which was 1.09-3.06 times higher than the suggested value of 10 −6 -10 −4 by the US Environmental Protection Agency [27], and higher than the suggested value of 5 × 10 −5 by the International Commission on Radiological Project (ICRP). Similar to the noncarcinogenic risk, dietary ingestion was the dominant pathway causing the estimated carcinogenic risk. In this study, dietary ingestion was the main pathway for heavy metal exposure, and the potential risk to health for all ages would increase once heavy metal pollutants were transferred into daily edible vegetables in the food chain. Further studies are needed to improve the accuracy of the model by focusing on the heavy metal exposure by dietary ingestion, to provide more detailed data on the health risk for the local community, and to provide a scientific and theoretical basis for local governments to consider adjusting agricultural production activities in the area. All this information can support the government to adjust agricultural production activities in the polluted area. The exposure risk of heavy metals to children is higher than that of adults due to the fact of several characteristics of children including eating, drinking, breathing, weight parameters lower than that of adults, and relatively weak metabolism and detoxification abilities. Hence, children are more easily affected by environmental hazards than adults [51][52][53]. As a result, the same doses of heavy metals may have negative effects on children's development while being nonhazardous to adult health.

Discussion
The accumulation of four heavy metal elements in the soil of the sewage irrigation area of the Xi River occurred in varying degrees, and the content of three metals (i.e., Cd,

Discussion
The accumulation of four heavy metal elements in the soil of the sewage irrigation area of the Xi River occurred in varying degrees, and the content of three metals (i.e., Cd, Pb, and Zn) had a clear downward trend from the upstream to the middle and lower reaches, where the upper reaches were significantly higher than the middle and lower reaches. For Cd, the high level of upstream heavy metal pollutants was due to the proximity to the Tiexi District, which is the largest recipient of industrial wastewater and urban sewage in the city and an important source of heavy metal pollution from the same sewage irrigation, domestic waste, and solid waste soil. The area around the midstream is also the main agricultural base of Shenyang; thus, the application of fertilizers and livestock manure may also lead to the accumulation of heavy metals in the surrounding soils. Relevant scholars have also confirmed that large and long-term application of organic fertilizers and phosphorus-containing fertilizers from livestock and poultry manure has led to the serious accumulation of Cd [54,55]. In addition, greenhouse vegetables are commonly cultivated with mulch, and heat stabilizers containing Cd are added during the production of mulch, which may also contribute to the accumulation of Cd [56,57]. At the same time, studies [58] showed that vegetable plots are transformed from farmland with a history of frequent use of heavy metal-containing pesticides, which can leave a large amount of heavy metal contaminants.
The contents of all four elements in the soil in this study were higher than the background values and accumulated significantly, but the contents in the vegetables were relatively low, and only one leafy vegetable samples had Pb contents exceeding the national standard, indicating that the contents of heavy metals in vegetables were influenced by factors other than the total amount of heavy metals in the soil. Studies have shown that there are many factors affecting the uptake of heavy metals in vegetables such as soil physicochemical properties, the heavy metal content in soil, morphological distribution, biological effectiveness, vegetable varieties, planting management conditions, and spatial differences [59][60][61][62][63]. Therefore, when studying the heavy metal content in vegetables, attention should be paid to the possible effects of factors other than the heavy metal content in the soil. In addition, in most studies, the weight of heavy metals was calculated on a dry weight basis, while some researchers used fresh weight, and whether the standard should be unified is an issue that deserves in-depth study.
Meanwhile, due to the lack of specific ingestion data for soil particles and local vegetables, the calculation of the risk assessment model could lead to uncertainties in this research. For instance, seriously polluted areas may be far away from residents, which may overestimate the daily exposure doses via dietary ingestion, inhalation ingestion, and dermal contact. In addition, residents tend to reduce outside activities in the cold and hot seasons to reduce the opportunities of exposure to heavy metals, causing overestimation. As this study mainly referenced overseas exposure parameters [64,65], there may be uncertainties caused by factors related to regional and community differences, economic and social development levels, and different lifestyles. Parameters published in the past may also need to be updated to more accurately represent the current situation [28][29][30][31]. Many exposure parameters suitable for Western countries may not be suitable for Chinese communities [21], which would inevitably increase the uncertainty of the results. Hence, considering the uncertainty in the model's parameters is necessary when interpreting the risk assessment in this study. Meanwhile, there may be other exposure pathways. For example, local residents may be exposed to heavy metals by drinking polluted underground water or ingesting Hg evaporated from the soil [51].
Although the abovementioned factors may bring some uncertainties, this study can still provide valuable scientific data for better pollution control and the planning of land use by local governments and policymakers.

Conclusions
This study investigated the accumulation of four heavy metals in the Xi River sewage irrigation area in Shenyang. The spatial distributions of the four metals were different, with the high-value area of Cu being scattered, while the concentrations of Cd, Pb, and Zn showed a declining trend from upstream to downstream of the Xi River. The levels of all heavy metals, except Pb, in the vegetable samples were lower than the limits in the national standard, with a slight difference. Considering the accumulation of heavy metals in the soil-vegetable system, the BCFs were in the order of Cd > Zn > Cu > Pb for all of the vegetables. The BCFs of the four heavy metals indifferent vegetables varied: leaf vegetables > Rhizome vegetable > Solanaceae vegetable. When total soil heavy metals were used for the risk assessment for children and adults by different routes of intake, the noncarcinogenic risks of the four elements were in safe ranges, while the carcinogenic risks of Cd was higher than the limit. Considering the strong ecological risk of Cd, protection measures should be taken to implement better pollution control and land use planning.
Author Contributions: All authors contributed to the study conception and design. Z.X., investigation, formal analysis, and writing-original draft; M.S., co-writing; X.Y., co-writing and review and editing; M.L., conceptualization and supervision. All authors have read and agreed to the published version of the manuscript.