添加链接
link管理
链接快照平台
  • 输入网页链接,自动生成快照
  • 标签化管理网页链接
About us
Share on
Download Download PDF ReadCube XML (NLM)

ORIGINAL RESEARCH article

Front. Environ. Sci. , 13 June 2022
Sec. Soil Processes
This article is part of the Research Topic Analysis, Fate and Transformation of Emerging Contaminants in Soil and Sediment View all 7 articles

Pollution Characteristics, Sources, and Health Risk Assessment of Heavy Metals in the Surface Soil of Lushan Scenic Area, Jiangxi Province, China

  • 1 School of Marine Science and Engineering, Nanjing Normal University, Nanjing, China
  • 2 Jiangsu Center for Collaborative Innovation in Geographical Information Resource Development and Application, Nanjing, China
  • 3 School of Geography, Nanjing Normal University, Nanjing, China

Heavy metal pollution at tourist attractions centers has caused widespread concern. In this study, the concentration of seven heavy metals (Pb, Cu, Cr, Zn, Ni, As, and Mn) in the surface soil of the Lushan scenic area was measured, and their pollution characteristics were assessed using the Nemerow comprehensive pollution index and geo-accumulation index; further, the human non-carcinogenic and carcinogenic risks were evaluated by the U.S. Environmental Protection Agency health risk assessment model. Correlation analysis, cluster analysis, and a positive matrix factorisation model were used to analyse the heavy metal sources. The results indicated that the heavy metal concentrations did not exceed the pollution threshold levels. The hazard quotients for the six heavy metals (Pb, Cu, Cr, Zn, Ni, and As) are below the threshold for non-carcinogenic health risk, whereas the carcinogenic health risk in the region is at an “acceptable” level; however, because of the high contribution of As, the carcinogenic health risks to residents of this area require continuous monitoring. Analysis revealed six sources of the pollutants: soil parent material, altitude, surface accumulation of organic matter, industrial dustfall, local human activities (life and agriculture), and tourism. These findings provide a scientific basis for developing appropriate strategies for a sustainable development of the scenic area.

1 Introduction

With rapid industrialisation and urbanisation in the modern society, the increasing soil pollution has become a serious concern. Soil heavy metals are receiving increasing attention as the main soil pollutants and critical components of typical human social impacts ( Huang et al., 2019 ). Under natural conditions, the major sources of heavy metals in the soil are biomass residues and the host rock that forms the soil ( Abdu et al., 2017 ). The soil does have the potential to self-purify heavy metals and other pollutants ( Dashko and Shidlovskaya, 2016 ). However, large quantities of solid wastes, such as electronic products and plastics produced by human activities, have negatively impacted the soil’s ecological environment ( Bongoua-Devisme et al., 2018 ). Furthermore, with the rapid development of the tourism industry, many tourists and construction workers enter scenic spots, inevitably generating domestic garbage, construction waste, and other wastes that add considerable manufactured heavy metal pollution to the scenic spot ( Ciarkowska 2018 ; Brtnický et al., 2020 ). Heavy metal pollution is characterised by long-term concealment, easy accumulation, difficulty in degradation, irreversibility, and high toxicity. It poisons the ecological environment, affects the growth of plants and animals, and enters the human body through food chain enrichment, seriously threatening human life and health ( Ali et al., 2019 ). The travel behaviour of tourists’ affects the soil heavy metal content of the tourist destinations; further, there is a risk of accidental ingestion of polluted soils because the scenic spots are perceived to be natural and clean. Therefore, measuring the heavy metal content of surface soil in tourist attractions, evaluating its pollution characteristics, and assessing potential health risks are of great significance to the sustainable development of environmental resources in tourist destinations. In most recent studies on scenic spots, the focus has tended to be on air pollution or water pollution, and these studies confirm the need for soil pollution research. For example, Lushan Mountain’s surface water has been heavily polluted because of tourism activities ( Xie et al., 2020 ). Textile industry wastewater has polluted coastal scenic areas through watersheds ( Sun et al., 2020 ). Through regional transport, short-duration high-intensity human tourism activities have resulted in air pollution at the scenic spots ( Du et al., 2021 ).

Several methods for assessing soil heavy metal pollution have been developed ( Chonokhuu et al., 2019 ; Ahn et al., 2020 ; Xiang et al., 2021 ). The Nemerow comprehensive pollution index method examines the impact of multiple contaminants in the soil and compares the study area as a whole with other areas or historical measurements. This method is widely used for soil and water pollution assessment ( Wei et al., 2019 ) since it provides an objective evaluation based on certain regulations. The geo-accumulation index method considers the background value of each element in the soil and the influence of diagenesis on the background value, and it is a common indicator of heavy metal enrichment in the soil ( Muller 1969 ). The U.S. Environmental Protection Agency (USEPA) health risk assessment model is widely used to quantify the carcinogenic and non-carcinogenic risks of heavy metals to humans based on the number of pollutants that enter the human body as well as their toxicity and carcinogenicity ( Tepanosyan et al., 2017 ). This method is applicable to almost all conditions of exposure of humans to toxic substances ( Arshad et al., 2020 ; Bayati et al., 2021 ). There have been health risk studies for heavy metals in beach sediments of Red Sea and Gulf of Aqaba and confirmed the safety for tourists ( Nour et al., 2022 ). To obtain multidimensional data, this study assesses the soil pollution characteristics and health risks of the Lushan scenic area with the above three methods. Lushan is one of China’s first national-level scenic spots, a world cultural heritage, and one of the nation’s top 10 famous mountains, thus, an important tourism spot. Statistically, the local population of Lushan Mountain was 277,555 in 2018, and received 61,801,700 tourists, with a total tourism revenue of 40.017 billion yuan. The question arises as to whether the large number of residents, the staff serving tourists, the large number of tourists, and the large number of vehicles moving between various scenic spots in the Lushan scenic area have a significant impact on pollution at the tourist destinations as well as on the health of the residents and tourists of these scenic areas. Currently, these aspects are not clearly understood. This study aims to evaluate the pollution characteristics and health risks of heavy metal pollution in the Lushan scenic area by sampling the surface soils of tourist residences and various scenic spots in the area and analysing the sources of heavy metals to provide a scientific basis for the sustainable development of the area.

2 Study Area

Lushan (29°26–29°41N, 115°52–116°08E) is a typical horst block mountain with an elliptical range located in the north of Jiangxi Province. Its north–south length is 29 km, its east–west width is approximately 16 km, and the mountain area is approximately 300 km 2 . The rocks exposed in the study area are mainly feldspathic quartz sandstone, quartz sandstone and gravel-bearing quartz sandstone. Lushan Mountain is at the northern edge of China’s mid-subtropical zone, which has a humid subtropical monsoon climate and at high altitude it has a distinct monsoon and mountain climate. The climate varies with altitude from subtropical to warm temperate to temperate ( Feng et al., 2019 ). Lushan Mountain’s multi-year average temperature is 11.9°C, and the annual average precipitation is 2,024 mm, which is 500 mm more than that during April to July months, the annual average relative humidity is 78%, and the annual average sunshine hours are 1,675.2 h, the frost-free period is 260 d, the annual average foggy days are up to 197.5 d, and the atmosphere is cloudy throughout the year ( Feng et al., 2019 ). There are more than 3,000 types of plants on Lushan Mountain, which has a subtropical and subalpine landscape with a forest coverage rate of 67%, and the vegetation is a deciduous broad-leaved forest and coniferous and broad-leaved mixed forest. Conversely, 100 m below Lushan Mountain is subtropical, and the vegetation is evergreen broad-leaved forest. Due to the vertical climate variation, the distribution of organisms and soil also varies; for example, from bottom to top, the soil types are: red soil, mountain red soil, mountain yellow-red soil, mountain yellow soil, and mountain yellow-brown soil ( Hu 2007 ).

3 Materials and Methods

3.1 Sample Collection and Chemical Analysis

In August 2018, 49 surface soil samples were collected in Guling Town, in the Lushan scenic area along the main tourist routes where tourists and residents gather. The distribution of the sampling points is shown in Figure 1 .

FIGURE 1

FIGURE 1 . Distribution of samples in Lushan scenic area.

After collection, all the samples were brought to the laboratory, dried in a constant temperature drying oven at 40°C, and the plant rhizomes and other sundries were removed, and approximately 10 g of each sample was ground in an agate mortar and then manually filtered through a 200-mesh sieve. Approximately 5 g of the ground sample was weighed for preparing tablets for analysis ( Shu et al., 2021 ). The sample tablets were placed into an X-ray fluorescence spectrometer (PANalytical Co., Almelo, Netherlands) for measurement; the GSS + GSD mode was selected, and the concentration of 27 elements was measured. The instrument selects Hongze Lake sediment GSS-9 (GBW07423) as the standard material for quality control, and the instrumental analysis error is less than 5%. The six heavy metals, Pb, Cu, Cr, Zn, Ni, and As, with evident biological toxicity that are also listed as priority control pollutants by USEPA, were identified. Furthermore, Pb, Cr, and As are listed as priority control metals by the Chinese government because of their high toxicity ( He et al., 2013 ). Mn is an essential trace element in organisms, but its excessive intake can cause central nervous system damage ( Sule et al., 2020 ). This study, therefore, primarily assessed the characteristics of heavy metal pollution and the associated health risks for the seven heavy metals Pb, Cu, Cr, Zn, Ni, As, and Mn.

3.2 Assessment of Heavy Metal Contamination

3.2.1 Pollution Characteristics

The Nemerow pollution index and geo-accumulation index were used to assess the characteristics of heavy metal pollution. Geo-accumulation index was calculated as follows:

I g e o = l o g 2 ( C n 1.5 B n ) ( 1 )

where C n is the measured value of the heavy metal n in the sampled soil (mg·kg −1 ); B n represents the geochemical background value of the heavy metal at the sampling site at Lushan scenic spot mg·kg −1 ; k represents the correction coefficient, considering the influence of diagenesis on the background value, k = 1.5 in this study ( Kong and Zhang, 2021 ). The sample with I geo >0, is concluded to have certain degree of contamination.

Nemerow comprehensive pollution index was calculated using Equations 2 , 3 :

P i = C i S i ( 2 )
P t = ( C i / S i ) a v 2 + ( C i / S i ) m a x 2 2 ( 3 )

where C i is the actual concentration of heavy metal pollutant i in the soil (mg∙kg −1 ); S i is the evaluation standard value of the heavy metal element i in mg∙kg −1 , which refers to the National Technical Regulations on the Evaluation of Soil Pollution Status ( MEPC, 2008 ) No. 39; P i is the single pollution index of heavy metal pollutant i in the soil; (C i /S i )av and (C i /S i )max are the average and maximum values of the single pollution index of heavy metals in the soil, respectively. Generally, P i ≤1 indicates the heavy metal is non-polluting; when 1 < P i ≤2, the heavy metal is in a light polluting state; when 2 < P i ≤3, the heavy metal pollution is moderate ( Yan et al., 2022 ). P t represents the pollution index of the seven heavy metal pollutants in the surface soil of the 49 sampling locations in the Lushan scenic area. When Pt ≤ 0.7, the sampling site is in a safe state; when 0.7 < Pt ≤ 1, it means that the sampling site is in a cordon state; when 1 < Pt ≤ 2, the sampling site is in a lightly polluted state ( Wang et al., 2022 ).

Interpolation analysis was used to estimate the above indices. The interpolation results of the Nemerow comprehensive pollution index represent the pollution status in the whole area, whereas those of the geo-accumulation index are used to analyse the spatial variations of each heavy metal element and provides a reference for source apportionment.

3.2.2 Health Risk Assessment

The health risk assessment of heavy metals in soil is widely used to quantify the carcinogenic and non-carcinogenic risks to humans during exposure. The basic equations for human exposure and health risk assessment (including the carcinogenic and non-carcinogenic risks of heavy metals) applied in this study are based on the recommendations and methodology for risk assessment introduced by the USEPA, including hazard identification, dose-response assessment, exposure assessment, and risk characterisation.

Step 1—Hazard identification: The aim is to determine whether exposure to these elements at the scenic spot can cause harm to human health and the probability of harm. Six heavy metals with chronic non-carcinogenic risks selected in this study included: Cr, Ni, Cu, Zn, As, and Pb. Concurrently, some heavy metals have carcinogenic risks. Therefore, these six heavy metals can harm human health after they accumulate in the human body beyond a certain limit, and the impact of these pollutants on human health is often irreversible ( Huang et al., 2007 ).

Step 2—Dose-response assessment: This primarily includes integrating the toxicity data of Cr, Ni, Cu, Zn, As, and Pb and their carcinogenic slope. The reference coefficients used in this study are primarily derived from the carcinogenic slope factor (SF) and reference dose (RfD) of pollutants published by the USEPA. In addition, the parameters from other health risk studies were referred. The specific values are presented in Table 1 .

TABLE 1

TABLE 1 . Values of reference dose (RfD) [mg·(kgd) −1 ] and slope factor (SF) [per mg·(kgd) −1 ] ( Chen et al., 2016 ; Liu et al., 2016 ; Jiang et al., 2017 ; Zeng et al., 2019 ; Yu et al., 2021 ).

Step 3—Exposure assessment: This determines the potential sensitive population in the Lushan scenic area as adults, and three exposure pathways were assumed: direct hand-mouth ingestion, respiratory system inhalation, and skin contact exposure ( EPA 1989b ; Ferreira-Baptista and De Miguel, 2005 ). The Formula for the assessment is based on the exposure parameter values such as the average weight and exposure time of the exposed population. The equations for the three exposure pathways are as follows:

A D D i n g = C × I n g R × CF × EF × ED B W × AT ( 4 )
A D D i n h = C × I n h R × EF × ED P E F × BW × AT ( 5 )
A D D d e r m = C × S A × CF × SSAR × ABS × EF × ED B W × AT ( 6 )

where ADD ing , ADD inh , and ADD derm , indicate the daily average exposure mg·(kgd) −1 of the human body exposed to the Lushan scenic area via ingestion, inhalation, and skin contact routes, respectively. The specific meanings and values of the related parameters are listed in Table 2 .

TABLE 2

TABLE 2 . Definition and the reference values of parameters for health risk assessment of heavy metals in soils ( EPA 2011 ).

Step 4—Risk characterisation: Since tourism is the major business in the area, tourists were considered as objects for risk characterisation. This step calculates the potential health risks caused by the six heavy metal pollutants in the soil of the study area to tourists visiting the area by integrating the access results of the first three steps. Hazard quotient (HQ) and hazard index (HI) quantify non-carcinogenic health risks. The non-carcinogenic health risks are calculated using the following equations:

HQ i j = A D D i j R f D ( 7 )
HQ i = j = 1 n H Q i j ( 8 )
HI = i = 1 n H Q i ( 9 )

where HQ ij represents the non-carcinogenic health risk of a certain heavy metal i through a single route j to the body of tourists exposed to the study area; HQ i is the total non-carcinogenic health risk probability of a certain heavy metal i through the three exposure pathways assumed in this study to the tourists exposed to the Lushan scenic area ( Chen et al., 2017 ); HI represents the total non-carcinogenic health risk probability of the six heavy metal pollutants through the three routes assumed in this study for the tourists exposed to the Lushan scenic area. When HQ i or HI < 1, indicates that the non-carcinogenic health risk can be ignored; when either of these is > 1, it implies a non-carcinogenic health risk ( EPA 1989a ).

Carcinogenic health risk equations:

CR i = ADD i × SF ( 10 )
TCR = i = 1 n CR i ( 11 )

where SF represents the slope factor of a certain heavy metal pollutant in the bodies of tourists exposed, CR i represents the probability of cancer patients in a certain number of tourists, and TCR is the total carcinogenic risk to tourists exposed to the Lushan scenic area. According to the Technical Guidelines for Risk Assessment of Contaminated Sites published by the Chinese government, generally, when TCR or CR i <10 –6 , it indicates no carcinogenic risk to tourists exposed to the Lushan scenic area; when 10 –6 <TCR or CR i <10 –4 , it indicates acceptable carcinogenic risk for tourists exposed to the Lushan scenic area. Pollutants within the acceptable carcinogenic risk range do not cause any harmful or undesirable health hazards to the exposed humans ( EPA 1989b ). Among the six heavy metal pollutants investigated in this study, Cr, Ni, As, and Pb have carcinogenic risks. Due to the lack of carcinogenic SF data for Ni and Pb exposure pathways, this study assessed the carcinogenic risk of Ni inhalation and Pb ingestion routes only. It must be noted that the carcinogenicity of Cr in the human body is due to the reduction of Cr 6+ and production of the highly reactive intermediates Cr 4+ and Cr 5+ ( Dayan and Paine, 2001 ). Cr 3+ is an essential trace element in the human body. The anthropogenic source of Cr 6+ is primarily industrial production ( Liang et al., 2021 ) and Cr exists in this hexavalent form mostly in alkaline soil ( Pradas Del Real et al., 2020 ) rather than the acidic mountainous soil as in Lushan Mountain (the simultaneous traceability work also confirmed the small industrial impact here). Therefore, it is believed that the carcinogenic risk of Cr in the soil of Lushan Mountain to tourists can be ignored and cannot be calculated.

3.3 Source Apportionment of Heavy Metals

Pearson correlation analysis was used to determine the correlation between heavy metals, and cluster analysis (CA) was used to predict possible sources. Data is normalized using Z-Score. In addition, the positive matrix factorisation (PMF) model was used for source analysis for better assessment of the sources and their contribution, and this method does not have unexplainable negative values in the matrix factorisation process compared with that of the PCA method ( Sofowote et al., 2008 ).

PMF was performed using receptor models ( EPA 2014 ):

x i j = k = 1 p g i k f k j + e i j ( 12 )

where x i j is the concentration of pollutant j measured in sample i , g i k is the contribution of source k to sample i , and f k j is the concentration of pollutant j in source k , and e i j is the residual for each sample/pollutant.

3.4 Data Processing and Spatial Analysis

Geostatistical methods are commonly used to study and characterise the spatial distribution of pollutants in soil. This method was used to describe the spatial distribution of the pollutant concentrations in this study. The data in this study were processed using SPSS, and interpolation processing was performed using empirical Bayesian kriging (EBK) method because it has better data adaptability and predictability than other interpolation methods ( Krivoruchko and Gribov, 2019 ) (although it also requires longer computing time). EBK interpolation was completed in ArcGIS 10.5, a box plot was constructed using R, and EPAPMF 5.0, provided by USEPA, was used for source analysis.

4 Results and Discussion

4.1 Metal Element Content and Pollution Characteristics

The concentration of the seven heavy metals in the soil samples from Lushan Mountain are listed in Table 3 . The background value of heavy metals in the soil samples in the table is taken from the “Seventh Five-Year” National Science and Technology Research Project National Soil Environmental Background Value Research Subtopic Soil Environmental Background Value Research in Jiangxi Province ( He et al., 2006 ). The national standard value is from the National Technical Regulations on the Evaluation of Soil Pollution Status .

TABLE 3

TABLE 3 . Variations in the concentration of heavy metal elements in the surface soil samples of Lushan Mountain (mg·kg −1 ).

After calculating the pollutant concentration in the soil sample using Formula 1 , the geo-accumulation index of the heavy metals at each sample point was obtained. Figure 2 shows the relationship between the distribution of the geo-accumulation index and heavy metal elements in box plots.

FIGURE 2

FIGURE 2 . Box plot of the relationship between the geo-accumulation index and heavy metal elements.

The average and median values of the accumulation index of heavy metals in the Lushan scenic area are less than 0, indicating that the area is generally pollution-free (see Table 3 ). Among the seven elements, the average values for the accumulation index of Cr, Ni, and Zn were relatively higher, and the Q3 (upper quartile) of Ni and Zn in the box plot exceeded the pollution-free warning line ( Figure 2 ). Chromium and its salts are widely used as raw materials in the chemical, electroplating, tanning, and several other industries, which are the sources of chromium-containing liquid and gaseous wastes ( GracePavithra et al., 2019 ). Ni pollution primarily comes from industrial wastes, gas, and fuel combustion ( Chaharlang et al., 2016 ). There are multiple sources of Zn, major being manufacturing activities such as smelting ( Li et al., 2015 ) and traffic exhaust gases ( Johansson et al., 2009 ).

The coefficient of variation (CV) statistically shows the degree of change in the observed values of the sample. Classical statistics classifies CV as follows: CV ≤ 0.1 indicates weak variation, 0.1 < CV < 1 indicates medium variation, and CV > 1 indicates high variation ( Nielsen and Wendroth, 2003 ). In this study, the CV of the observed values of the concentration of each element is in the medium variation range, indicating that each element in the soil has a medium degree of spatial variation. Among these, Mn had the highest CV because Mn is primarily derived from soil-forming parent material in the soil, and it is inferred that the spatial variation of heavy metals in the study area is affected by the manufactured and natural sources.

For a better understanding of each element’s spatial variation, I geo values were used, and the search neighbourhoods of I geo Zn and I geo Pb were found to be smooth circles, whereas the others were standard circles. Compared with the concentration value, the I geo value not only reflects the spatial variation of the element’s concentration but also standardises the concentration of each element, which is more relevant for discovering the similarities and differences of the changes caused by human activities. As shown in Figure 3 , some commonalities in the spatial distribution of certain elements can be found: Ni, Cu, and Mn have high-value areas in almost the same location (north of the town), and their impact has spread to the town; As and Pb have similar high-value areas to the east of the town and the second high-value area to the west of the town; the distribution of Zn is similar to that of Pb, but Zn has a high-value area in the town; Cr has two high-value areas to the north (similar to the Cu group) and west of the town which encircle the town. These commonalities imply that some of the pollutants have similar sources, but there is no doubt that the town is more affected by the surrounding heavy metal sources than internal sources. This phenomenon implies that in areas with complex geological environments such as the Lushan Mountain, the effect of soil parent material on heavy metals cannot be ignored ( Zinn et al., 2020 ). This study focuses more on the phenomenon of accumulation of the heavy metal pollutants. These artificially constructed groups of elements can be used as motivations or references in subsequent studies on processes and accumulation mechanisms.

FIGURE 3

FIGURE 3 . The spatial distribution of I geo values after empirical Bayesian kriging (EBK) interpolation.

Equation 2 and 3 were used to calculate the Nemerow comprehensive pollution index (Pt) for each sample location, and the data are shown in Table 4 . The Pt value is interpolated by the EBK method ( Figure 4 ) to show a spatial change in the image map. The data transformation type is empirical, the semi-mutation model type is K_BESSEL, and the search neighbourhood is a standard circle. The layer transparency was maintained at 20%.

TABLE 4

TABLE 4 . Lushan scenic area soil heavy metal Nemerow comprehensive pollution index.

FIGURE 4

FIGURE 4 . The spatial distribution of Pt values after empirical Bayesian kriging (EBK) interpolation.

The highest Pt value in the study area is 0.92 ( Table 4 ), which is below the pollution warning line. The highest value appeared at sample point No. 32 near Longgong Road. Notably, although the colour bar used in the interpolation layer expresses high values in red, the highest value in the interpolation result is only 0.65. After the interpolation process, the Pt value of this area is in a low and stable state (0.567–0.652), two high-value areas appear northeast and southwest of the town, and the urban area is a low-value area. The high-value areas in the northeast are primarily distributed in the tea plantations and scenic spots, and the high-value area in the southwest is centred on the famous scenic spot Datianchi in Lushan Mountain. This phenomenon shows that the heavy metal content in the region tends to be affected by tourism activities and natural factors.

4.2 Health Risk Assessment Results

4.2.1 Exposure Assessment

The major source of activity in this area is tourism, hence, this study considered the tourist group as the object to conduct health risk assessment. According to the recommended playtime in the Lushan scenic area, the duration of play was 2–3 days, the exposure duration was ED = 2.5d = 0.08a, and the exposure frequency EF was negligible. Furthermore, Report on the Nutrition and Chronic Disease Status of Chinese Residents shows that the average weight of the Chinese population is BW = 61.9 kg. Furthermore, it is assumed that each tourist visits the Lushan scenic area only once in his lifetime; therefore, the carcinogenic and non-carcinogenic substance exposure time AT is considered to be 0.08. Table 5 shows the average daily exposure doses from different exposure routes in the Lushan scenic area.

TABLE 5

TABLE 5 . Average daily exposure dose from different exposure routes in Lushan scenic area.

The data in Table 5 shows that tourists visiting Lushan scenic area are likely to be exposed to a maximum daily average pollutant dose through hand-mouth ingestion, followed by the skin contact route, and the least exposure is via inhalation. This suggests that ingestion is a significant route of human exposure to heavy metals, which is consistent with the findings of previous studies ( Chonokhuu et al., 2019 ; Yu et al., 2021 ). Therefore, developing good outdoor hygiene habits: washing hands frequently and avoiding hand-mouth contact are effective ways to reduce heavy metal exposure (and also infectious diseases).

4.2.2 Non-carcinogenic Health Risks

Similar to the data for exposure assessment ( Table 5 ), the non-carcinogenic health risks in the Lushan scenic area are also dominated by hand and mouth ingestion ( Table 6 ). The HQ contributed by Cr and As is higher (85% in total); however, they are still far below the safety threshold and not likely to cause any non-carcinogenic health risks to tourists exposed to the Lushan scenic area. Due to the variability of environment, it is important to take conservative estimates when assessing health risks. In the non-carcinogenic health risk assessment, the different valency states of Cr are not distinguished, and the RfD of the more toxic Cr 6+ is uniformly used. Likewise, the RfD for inorganic form of As was used to calculate the HI for total As. Therefore, the HQ for these two heavy metals is overestimated.

TABLE 6

TABLE 6 . Non-carcinogenic health risks under different exposure routes.

4.2.3 Carcinogenic Health Risks

The carcinogenic health risks calculated using Eq. 10 and Eq. 11 according to the carcinogenic SF data are shown in Table 7 .

TABLE 7

TABLE 7 . Carcinogenic health risks of heavy metals in Lushan scenic area.

According to the USEPA soil treatment standards, since the calculated carcinogenic health risk values for Ni and Pb are lower than 10 –6 , their carcinogenic risk is nil, and the risk values for As are lower than 10 –4 , which fall within the acceptable range. The highest carcinogenic health risk is from As, and oral intake is the most probable route of exposure. The TCR was lower than 10 –4 , which is still acceptable; thus, we can conclude that in the Lushan scenic area, concentration of none of the heavy metals is high enough to cause cancer. However, the accumulation of heavy metals in the soil surface under the influence of numerous factors may pose a threat of low carcinogenic health risk in the Lushan scenic area. Hence, there is a need for continuous attention in this regard. Carcinogenesis being a long-term process, such a minimal risk may pose a threat to the residents living in towns within the scenic spot.

Although the non-carcinogenic and carcinogenic health risks in the study area are at safe levels, high HIs from Cr and As and high CRs from As indicate that Cr and As are potential health threats in the area. This trend is similar to the health risk assessment data of a rural-industrial town in Jiangsu Province from a previous study ( Jiang et al., 2017 ). To protect the health of tourists and residents in Lushan scenic area, the anthropogenic sources of Cr and As need to be carefully monitored.

4.3 Analysis of Possible Sources of Heavy Metals

Source identification is a challenge due to lack of historical data, since most of the pollution monitoring studies have focused on plains due to the poor continuity and regularity of mountain soils. In this study, several analytical methods and reference factors were used in combination for traceability and to obtain reasonable and more realistic results. Because the concentration of total phosphorus (TP) in the soil is related to agricultural activities and natural biological processes ( Katayama et al., 2021 ; Zhang et al., 2021 ), it is a better indicator for the source. Lushan Mountain has complex geomorphic conditions, and altitude is considered a potential factor affecting the spatial distribution of heavy metals ( Ding et al., 2017 ). Therefore, in the traceability analysis, the TP and the altitude of the sample point were used as references.

The Pearson correlation coefficient used to determine the potential relationship between the concentration of heavy metals ( Table 8 ), found a significant correlation between most heavy metal elements. TP was significantly positively correlated with Cu, Zn, and Pb, and significantly negatively correlated with Ni, and altitude was significantly negatively correlated with Cr. Such a complicated relationship shows that various factors are affecting the distribution of heavy metals.

TABLE 8

TABLE 8 . Pearson correlation coefficient of the concentration of each element.

In CA, the elements are divided into four categories at 15 distances ( Figure 5 ): Zn, Pb, Cu, As, and TP are considered in one category; Cr and Ni are placed together; Mn and altitude are considered as one. At a distance of 20, Mn, Cr, and Ni were classified into one category. This shows that altitude has a limited impact on the distribution of heavy metals, and as an influencing factor on mountain soil and parent material. Mn is primarily affected by natural sources in the soil ( Dong et al., 2019 ), whereas TP tends to be affected by human activity; therefore, it can be temporarily placed in a group that is more affected by human activities, such as Zn, Pb, Cu, and As, and a lower group, Cr, Ni, and Mn.

FIGURE 5

FIGURE 5 . Heavy metal cluster analysis (CA) dendrogram.

In running the PMF model, the residual error (95.2%) of most data was within ±3% when six factors were selected. Therefore, it is believed that the PMF model identified six possible sources ( Figure 6 ). Figure 7 shows the proportion of each element’s contribution to this factor. Factor 1 contributes most of the Mn, which is considered to be the parent material for soil formation; Factor 2 primarily comes from the altitude and has little effect on each element, which is considered to be the effect of altitude on the concentration of heavy metals; Factor 3 primarily comes from TP, which contributes significantly to As. According to the Soil Environmental Background Value Research in Jiangxi Province , surface accumulation of As in the Lushan Mountain soil is evident. Other elements are more leached, whereas the surface accumulation is considered to be caused by organic matter. Factor 4 primarily comes from altitude, which exerts a certain impact on Cr, Ni, Mn, and Zn and is considered to be an industrial source (local construction and industrial dustfall); Factor 5 primarily comes from TP, which has an impact on Pb, Zn, and Cu. It is considered a local human influence (life and agriculture) in Lushan Mountain. Factor 6 has almost no effect on Mn and TP and indicates more effect on Cr, Ni, and Cu, which is considered the traffic impact caused by tourism. After removing the contribution of each element to altitude and TP, the impact on heavy metals was analysed, and the proportions of each factor was 43.24, 16.53, 7.14, 20.53, 5.59, and 6.96%, for factors 1–6, respectively, indicating that the parent material plays a major role in the heavy metal concentration, followed by industrial sources and tourism related activities.