Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Scientific Reports volume 14, Article number: 21533 (2024 ) Cite this article Rb395 Max Rebar Tier
Soil heavy metals (HMs) pollution is a growing global concern, mainly in regions with rapid industrial growth. This study assessed the concentrations, potential sources, and health risks of HMs in agricultural soils near marble processing plants in Malakand, Pakistan. A total of 21 soil samples were analyzed for essential and toxic HMs via inductively coupled plasma‒optical emission spectrometry (ICP‒OES), and probabilistic health risks were evaluated via Monte Carlo simulation. The concentrations (mg/kg) of Ca (29,250), P (805.5) and Cd (4.5) exceeded the average shale limits of 22,100, 700, and 3.0 mg/kg, respectively, and indices such as Nemerow’s synthetic contamination index (NSCI) and the geoaccumulation index (Igeo) categorized the soil sites as moderately polluted. The potential ecological risk index (PERI) indicated considerable to high ecological risk for As and Cd. The deterministic analysis indicated non-carcinogenic risks for children (HI > 1), whereas the probabilistic analysis suggested no significant risk (HI < 1) for both adults and children. Both methods indicated that the total cancer risk for Cr, Ni, Cd, and As exceeded the USEPA safety limits of 1.0E-06 and 1.0E-04. Sensitivity analysis identified heavy metal concentration, exposure duration, and frequency as key risk factors. The study suggested that HM contamination is mainly anthropogenic, poses a threat to soil and human health, and highlights the need for management strategies and surveillance programs to mitigate these risks.
Soil is an indispensable natural resource, fundamental to food production and biodiversity maintenance, and plays an important role in national security and global ecological stability1,2. However, the integrity of this resource is increasingly threatened by the pervasive issue of soil pollution, which has become a pressing global environmental concern3. The contamination of soil with hazardous substances such as radioactive elements, nanomaterials, and heavy metals significantly undermines soil resilience and disrupts ecosystems4,5. These pollutants primarily originate from anthropogenic activities, including industrialization, urbanization, mining, and the excessive use of agrochemicals, as well as natural events such as volcanic eruptions and erosion of mineral deposits6,7. The global significance of soil pollution lies in its far-reaching impacts on agricultural productivity, human health, and environmental sustainability8,9. Contaminated soils can lead to the accumulation of toxic elements in the food chain, posing severe risks to food security and public health on a global scale10,11. Moreover, soil degradation due to pollution exacerbates the loss of arable land, which is already under pressure from a growing global population12,13. Addressing soil pollution is therefore critical to ensuring the long-term viability of global food systems, the health of populations, and the preservation of biodiversity14,15.
Millions of tons of marble waste, consisting of fine powder and fragmented rock, are generated worldwide in marble quarries and marble processing plants (MPPs)16,17. During the quarrying and refining processes, nearly 70% of the valuable mineral resources of marble are lost, with approximately 40% of marble waste being released into the surrounding atmosphere, leading to widespread environmental pollution18,19. MPPs, particularly in Khyber Pakhtunkhwa, Pakistan, contribute significantly to wastewater and environmental pollution, as large amounts of waste are released as sludge and sludge containing harmful marble particles20,21,22,23,24. This waste, which is often disposed of improperly in water bodies and land, poses severe ecological risks, including damage to aquatic ecosystems, soil, and human health20,25,26. Despite the introduction of pollution control programs by the government, poor enforcement has hindered their effectiveness27,28. Although the marble industry plays a crucial role in socioeconomic development worldwide, growing concern about its impact on the environment has led to an increasing focus on recycling waste and byproducts to align with sustainable development goals19,29. Recently, modified biochar composites (Mar-BC800 and Sep-BC800) were developed from marble waste, calcium-rich sepiolite, and agricultural waste, and have been shown to be highly effective, eco-friendly, and cost-efficient for phosphate removal from waste streams30. Moreover, their potential for repurposing as fertilizers presents a sustainable solution for enhancing soil fertility and improving soil properties31.
In developing countries without prior treatment, farmers use marble wastewater as an alternative to irrigating crops because it is a low-cost and readily available resource32,33,34. In the short term, marble wastewater can have a positive effect on plant growth as it contains high quality nutrients such as calcium, potassium, iron and copper24,25. However, farmers may not be fully aware of the potential long-term consequences of this practice, as marble waste also contains heavy metals35,36,37. Heavy metals such as Al, As, Cd, Cr, Co, Ni, Pb, and Hg are toxic to plants even at low concentrations38,39,40. The accumulation of these elements at relatively high concentrations in agroecosystems poses a major threat to soil microbes and disrupts the distribution and structure of the microbial community4,41,42. Similarly, their toxic behavior in soil disrupts the morphological and physiological processes of crops by reducing their growth rate, nutrient imbalance, stomatal movement, and inhibition of photosynthesis43. In human and animal body tissues, soil heavy metals accumulate through the food chain44,45 and cause serious effects such as nephrotoxicity, neurotoxicity, hypertension, infertility, and carcinogenicity20,46,47. Moreover, prolonged exposure to marble waste through inhalation and dermal contact poses serious health risks to local populations and workers, including conditions such as urolithiasis and abnormal lung function16,17,22,35,48. Therefore, it is crucial to assess the risk and extent of soil heavy metal contamination near MPPs49. Such assessments are essential for developing effective mitigation strategies to protect both the environment and public health6,50.
Geochemical factors in soil, water and sediments contain valuable information on element ratios and should be treated as compositional data51. Using classical statistical tools directly on the input concentrations can lead to spurious results because different geometric rules govern compositional data52. Traditional methods such as Euclidean distance and covariance are unsuitable for geochemical concentrations, resulting in spurious elements53,54. Geochemical data are positive vectors, represented by a constant sum, that are distinct from real Euclidean space55. Log-ratio transformations such as the centered log-ratio (CLR), isometric log-ratio (ILR) and additive log-ratio (ALR) can handle compositional element concentration data56,57,58. CLR provides an alternative to raw data but introduces collinearity, whereas ILR avoids collinearity57,59. The application of CLR and CoDa principles, such as variation matrices, biplots, and dendrograms, allows the establishment of linear relationships between CLR variables and the differentiation of geochemical sources42,60. Identifying unusual patterns related to heavy metal exposure in the area of interest is crucial26.
Similarly, the evaluation of soil pollution status, identification of potential sources of contamination and appraisal of the prospective risks to human populations are prerequisites for the successful implementation of any action plan61. Soils near the marble industry are most likely polluted by marble effluents, marble slurry dumping and the production of marble powder from marble processing plants25. Numerous studies have been conducted in Pakistan on heavy metal water pollution related to the marble industry20,37. These studies have shown that marble industry effluents strongly affect water and sediment quality and cause the bioaccumulation of heavy metals in aquatic fauna and flora. Similarly, with respect to health risk assessment, most studies have typically focused on deterministic risk quantification of exposure to soil contaminated with heavy metals, which could either underestimate or overestimate the risk62,63 In addition, owing to the uncertainty of heavy metal concentrations and the specific variability between individuals, identifying the most hazardous element for the population via a deterministic method is difficult64,65. Fortunately, as the most common method of probabilistic risk analysis, Monte Carlo simulation can determine the likelihood of a risk being exceeded and help identify the priority metal for risk control66,67,68. Therefore, a comprehensive investigation that considers the metal enrichment status and probabilistic risk assessment of heavy metals at the national level is valuable for understanding the contamination traits and contributions of various sources prior to taking efficient measurements to improve soil quality and human health61,69. In addition, for the production of healthy crops and ensuring the health security of farmers, researchers, soil managers, livestock and common people, it is essential to study agricultural soils contaminated with marble waste25,70,71,72. Therefore, the current study aimed to (i) evaluate the extent of heavy metals (HMs) contamination in agricultural soils withassociated human health risks using both deterministic and probabilistic health risk assessment modelsand (ii) identify potential sources of heavy metal contamination in soil via a compositional multivariate technique.
The present study was conducted between September and December 2019 on agricultural soil near MPPs in the district of Malakand (Fig. 1).
The study area map showing sampling sites was generated via ArcGIS software version 10.2.2.
Geographically, District Malakand is located between 34.3718° N latitude and 72.2420° E longitude and is the gateway to Dir, Chitral, Swat, Shangla, and Buner Districts and the currently merged, federally administered, tribal area of Mohmand and Bajaur. Malakand District covers an area of 952 km2 with a population density of 475 inhabitants per km273. The land of Malakand District is fertile and surrounded by mountains rich in biodiversity. A large part of the agricultural land in the study area is irrigated by the Swat River (Fig. 1). Winters are very cold, and the temperature can decrease to – 2 °C in January. Summers are hot, and the temperature reaches 41 °C from June to July74. The average annual rainfall is between 600 and 650 mm75. Rice is an important cash crop with a total irrigated area of 4991 hectares76. The valuable mineral reserves reported from several locations in the Malakand district are mica, quartz, marble, Malakand granite, and Dargai chromite77. The coordinates of the individual sampling points were recorded by a global positioning system (Garmin eTrex 10). The study area map (Fig. 1) was generated via ArcGIS version 10.2.2.
The agricultural soils in the study area are sandy loam. Soil samples were collected from 21 sites (Fig. 1) in triplicate at a depth of 0–20 cm via a hand auger following a quincunx sampling pattern78. Finally, one kilogram composite sample for each site was prepared by thoroughly mixing the five sub soil samples and sealed in clean, labeled polyethylene bags79,80. The soil samples were air dried, mechanically ground, sieved through 2 mm sieve and stored for further analysis81. 0.5 g of the dried soil sample was added in triplicate to a 50 mL Erlenmeyer flask with 26 mL of a triacid mixture of HNO3, HCl, and HF at a ratio of 9:3:1 and stored overnight82. The samples were completely digested on a hot plate by gently increasing the temperature to 155 °C until white dense fumes emerged83. After cooling, the digested sample was diluted to 50 mL with distilled water, filtered through Whatman filter paper (0.45 μm) in a clean volumetric flask and stored at room temperature84. The concentrations of essential elements and heavy metals in the extracted soil solutions were analyzed using inductively coupled plasma optical emission spectrometry (ICP‒OES; Thermo Scientific, iCAP 6000, UK).
A soil suspension was prepared using a soil to water ratio of 1:2.5 in double deionized water, then stirred and centrifuged at 3000 rpm for 4 h. The pH, electrical conductivity (EC), and total dissolved solids (TDS) of the soil samples were determined via a portable meter (Model PH-2012, WT01 China)85.
For quality assurance and quality control, special attention was given to sampling, preservation, transportation, storage and each experimental technique used in laboratory analysis86. Analytical-grade acids and reagents provided by Merck (Germany) were used in this study. The standard reference material (SRM) used to construct calibration curves of the tested elements was IAEA-Soil-7. The reference soil material (IAEA-Soil-7) was prepared and analyzed according to the certified analytical standard quality control procedure87. Eight different linear concentration standards were prepared. Each element with a good linear plot and a correlation coefficient greater than 0.999 was observed in the construction of standard curves. The traceability of the standards was verified via high-performance ICP‒OES against the corresponding SRM (ISO 2009). Each sample was analyzed in triplicate, and the mean value was used to interpret the results. The value of the relative standard deviation (RSD) for the analyzed elements was ≤ 2.0% (Table S1). The limit of detection (LOD) for each element was determined by digesting five blanks (Table S1). The recoveries were within the EPA limits of 92–123% (Table S1).
In the present study, the soil heavy metal pollution status was assessed through pollution indices developed by Muller87,88; Hakanson90 and Tomlinson et al.91.
The pollution index (Pi) reflects the pollutant level and primary exposure factors, serving as a tool and is used to assess heavy metal pollution in soil92. Pi can be calculated via the following equation.
where Ci is the measured concentration of the ith element and Bi is the background value of the ith element in the world average shale93. The Nemerow synthetic contamination index (NSCI) is used to represent the overall level of soil contamination by heavy metals. Soil HM contamination is a complicated process triggered by the synchronized action of numerous heavy metals94. The NSCI can be calculated via Eq. 2.
where Piave is the mean value of Pi for the analyzed HM and Pimax is the maximum Pi value of the HM in a sample. The soil risk categories based on the Pi and NSCI are listed in Table S2.
The index of geo-accumulation (Igeo), originally defined by Muller88 is used to estimate the natural and anthropogenic factors of soil contamination95. This method assesses the environmental impact while considering preindustrial and existing metal concentrations. The Igeo is calculated using the following equation, adapted from El Azhari et al.96.
where Cn signifies the measured concentration (mg/kg) of metals in the soil sample and Bn is the geochemical background (mg/kg) of the tested metal. A factor of 1.5 compensates for background values of the respective metal owing to minor environmental and anthropogenic influences. On the basis of Igeo, Muller88 proposed seven risk categories, as shown in Table S2.
Eri measures the potential ecological risk of a single element in soils90 and can be calculated as follows:
where \({\text{E}}_{\text{r}}^{\text{i}}\) is the potential ecological risk factor for the ith trace element and where \({T}_{r}^{i}\) is the coefficient of the ith element, which is 10 for As; 30 for Cd; 05 for Co, Cu and Ni; 01 for Mn and Zn; and 02 for Cr and V90,97,98,99,100. \({\text{P}}_{\text{i}}^{\text{i}}\) is the pollution index of trace elements taken from Eq. (1). The evaluations of Eri are listed in Table S2.
The comprehensive description of the ecological risk from a wide range of potentially toxic metals was determined through the application of the potential ecological risk index (PERI)90,101. The PERI results were calculated via the following equation.
where PERI designates the cumulative ecological risk of trace elements, \({\text{E}}_{\text{r}}^{\text{i}}\) is the potential ecological risk factor for the ith trace element and ni is the number of trace elements, which is 09 in this study. \({T}_{r}^{i}\) is the coefficient of the ith element, and \({\text{P}}_{\text{i}}^{\text{i}}\) is the pollution index of the trace elements taken from Eq. 1. The evaluations of PERI are listed in Table S2.
Human exposure to soil heavy metals (HMs) typically occurs through accidental ingestion, inhalation, and dermal contact with soil particles from water and air102,103. In the present study, the noncarcinogenic health risks from individual heavy metal across the three exposures routes were calculated via the following equations 104.
where ADDing, ADDinh, and ADDdermal are the average daily doses (mg/kg/day) of HM through ingestion, inhalation, and dermal contact, respectively. All other factors and values used in Eqs. 6–8 to estimate the exposure level and risks are listed in Table S3.
HQ is the ratio of the average daily dose (ADD; mg/kg/day) and the reference dose (RfD; mg/kg/day) for a given toxicant. HQ values > 1 signify the occurrence of a noncarcinogenic risk105 and can be calculated via Eq. (9).
where HQing, HQinh, and HQdermal are hazard quotients (HQs) for the exposure routes, and RFDing, RFDinh, and RFDdermal represent the corresponding reference doses listed in Table S4.
The hazard index (HI) estimates the overall noncarcinogenic health risks from exposure to multiple potentially toxic HMs. HI is the sum of all HQs computed for individual HMs for a given route of exposure106. The HI can be calculated via the following equation.
where 1, 2, …, n is the HQ calculated for individual toxic HMs. The population is considered safe if HI < 1 and vulnerable if HI > 1107.
The carcinogenic risk (CR) assessment estimates the likelihood of an increased risk of cancer for a person due to lifetime exposure to potentially toxic elements. The cancer risk from lifetime exposure to As, Cr, Cd, and Ni was calculated across the three exposure routes via the following Eqs. 105.
where CR is the cancer risk from the individual exposure pathway, and TCR is a cumulative cancer risk. CSF is the cancer susceptibility factor (mg/kg/day) of the HMs listed in Table S4. According to the USEPA108, cancer risk values below 1 × 10−6 are considered insignificant, while those higher than 1 × 10−4 are likely to be detrimental. According to the classification of the WHO109, As, Ni, Cd, and Cr (VI) compounds are being studied as Group 1 carcinogens; Pb compounds as Group 2 A; and Co and its compounds as Group 2B. In contrast, Cu, Zn, and Mn are not classified as carcinogenic, and their cancer risk has not been studied.
Uncertainty is inherent in risk assessments because of limited precise information, variability in environmental systems, and differences in individual human characteristics110,111. To address these uncertainties, Monte Carlo simulation, a random sampling method based on probability theory, is commonly used to quantify variability and uncertainty in risk assessment112. To achieve this, the probability distributions of the variables listed in Table S3 were used as input parameters to evaluate the probability functions for the health risks113. After performing 10,000 iterations, stable exposure distribution results were derived, and the values at quantiles 5th and 95th of the exposure distribution were used to assess the probabilistic risk114. A sensitivity analysis was conducted to evaluate the degrees of contribution of the exposure factors to the results. A positive value indicates a positive correlation between the exposure factor and health risks, while a negative value indicates a negative correlation115,116.
Descriptive statistical characteristics, including the mean, standard deviation (SD), standard error (SE), skewness, kurtosis, and coefficient of variation (CV), of the 19 elements were calculated via SPSS software (version 25 for Windows). Monte Carlo simulation, sensitivity analysis and the best-fitting distribution for each factor were performed via Oracle Crystal Ball version 11.1.3.0. Compositional data analysis approaches such as centered log-ratio (clr) transformation, variation matrix, clr-biplot, coda-dendrogram, and factor analysis were applied to establish a linear relationship between the variables and differentiate the source of contamination in agricultural soils42,59,60,117. The variation matrix (Table 4), clr-biplots (Fig. 14a,b), and Coda-dendrogram (Fig. 15) were generated via Coda-Pack software118 and the R package ‘‘compositions’’119. The factor analysis (Table 5) was conducted using the statistical software JMP version 16.0120.
The current research study was approved by the Institutional Research Committee of the Department of Botany at Islamia College Peshawar. The article does not include any humans or animals as research objects. All methods were conducted following the appropriate guidelines outlined in the methods section.
The assessment of physicochemical properties and heavy metal concentrations in soils near marble processing plants showed significant variation in the studied parameters (Table 1). The soil pH, ranging from 7.3 to 9.0, with a mean of 8.3, indicates slightly alkaline conditions. Similarly, EC and TDS values ranging from 213.3 to 556.6 µs/cm and 2231 to 1072.67 mg/L, with mean values of 374.60 and 635.80, respectively, indicate variability in soil salinity levels. These variations in soil pH and salinity may be attributed to differences in water quality, soil composition, and the influence of industrial activities, such as marble processing, on the surrounding environment23,121. High salinity levels could lead to soil degradation, impacting agricultural productivity122. Similarly, soil alkalinity impacts both physical and chemical properties, including texture, structure, color, infiltration, porosity, surface crusting, swelling, pH, and nutrient availability123. Managing soil alkalinity and salinity through proper irrigation practices, soil amendments, and crop selection is crucial for maintaining soil health and ensuring successful agricultural productivity124,125. Similar results were also reported by Eman et al.126, which is consistent with our current results.
The highest mean concentration (mg/kg) was observed for Ca (19881.50), followed by Al (10914.30), Mg (5420.71), Fe (4,699.71) and K (2105.70) (Table 1). Titanium (Ti) had a mean concentration of 600.50 (mg/kg), whereas phosphorous (P) recorded a mean concentration of 292.0 (mg/kg). The maximum concentrations (mg/kg) of Ca (29250) and P (805.5) exceeded the world average shale concentration (Table 1). These findings are consistent with those of Aukour and Al-Qinna33, Ozcelik127, Ahmad et al.16 and Noreen et al.128, who reported elevated levels of these elements in soils and water impacted by industrial activities such as marble processing and mining. The high Ca concentrations likely reflect the presence of calcareous materials, which are common in regions with intensive marble processing129. Elevated P levels could be linked to industrial discharge or the accumulation of phosphate-rich dust130. Such elevated concentrations may have ecological implications, including potential disruptions to soil nutrient balances and risks of eutrophication in nearby water bodies131. The findings underscore the need for monitoring and managing these elements to mitigate potential environmental and agricultural impacts.
The mean concentrations (mg/kg) of Mn (149.30), V (39.21), Sr (35.14), Cr (27.0), Cu (20, 10) and Ni (18.40) were found to be within the limits of world average shale (Table 1), indicating relatively lower contamination levels105. Similarly, the mean concentrations (mg/kg) of Si (12.50) Co (6.0), Zn (5.80), As (5.30), Zr (3.0), and Cd (2.0) were below the FAO permissible limits (Table 1), suggesting a lower risk to soil quality and ecosystem health132. However, the maximum concentration of Cd (4.50 mg/kg) exceeded the FAO limit of 3.0 mg/kg (Table 1), raising concerns about potential soil and human health risks133. Similar results were also reported by Khan et al.36 for Cd concentrations in different marble industries of Khyber Pakhtunkhwa, Pakistan, highlighting consistent patterns of cadmium accumulation in soils associated with marble processing activities in the region. Elevated cadmium levels have been associated with industrial emissions and waste, particularly in areas affected by mining and manufacturing activities134. Cadmium is a toxic metal that can accumulate in the food chain, posing long-term health hazards135.
The coefficient of variation (CV) quantifies the spatial variation in geochemical variables136. In the present study, the coefficient of variation (CV) values varied from 22.40 to 106.70% (Table 1). Low variations were obtained for Ni, Ca, Cr, Mg, V, and Zr; moderate variations were obtained for Si, Fe, Cu, Mn, K, Al, As, and Co; and high variations were obtained for Ti, P, Cd and Zn (Table 1). According to Zuo137, small CV values are associated with areas that have similar geochemical backgrounds, while large CV values characterize areas that have different geochemical backgrounds. The irregular distribution of elements in the current study area seems to be influenced by anthropogenic rather than natural influences138. The statistical skewness and kurtosis indices were used to ascertain the normality or abnormality of the heavy metal distribution139. The kurtosis values of Al, Fe, P, Ti and Zn were greater than zero, indicating steeper than normal distributions of these elements (Table 1). The skewness values of Al, K, Fe, P, Ti and Zn were greater than one, indicating that these elements were positively skewed to lower concentrations, highlighting potential environmental concerns140,141. Understanding these associations helps in identifying the sources of soil contamination and in developing strategies to mitigate their impact on soil health and agricultural productivity142.
Pollution indices are among the most effective tools for assessing heavy metal contamination in soils, offering more insight than simply measuring total metal concentrations143. In this study, single indices were calculated to evaluate the pollution levels of individual heavy metals in the soil. These indices are based on the total metal concentrations, along with geochemical background levels or preindustrial baselines, to provide a comprehensive assessment of soil contamination. Similarly, the Nemerow synthetic contamination index (NSCI) is employed to assess contamination levels across different sites144. The results of the pollution index (Pi) (Table S5) show that the Pi values (1 ≤ 2) for As are indicative of slight pollution at sites S1-S3 and S14; larger variations in Pi values (2 ≤ 3) for As indicate moderate pollution at soil sites S6-S9 and S15–S18; and the Pi values > 4 are indicative of heavy pollution at soil sites S4, S11-S14, and S19-S21. Similarly, Pi values between 1 and 2 for Cd indicate slight pollution at soil sites S1, S6, and S8-S10; moderate pollution (2 ≤ Pi < 3) at sites S2–S5, S7, and S17; and heavy pollution (Pi > 4) at soil sites S11–S16 and S18–21. The Nemerow synthetic contamination index (NSCI) values for heavy metals varied from 1.20 to 4.11, with a mean of 2.55 (Table S5). The results indicate that 38% of the soil samples were slightly contaminated, 29% were moderately contaminated, and 33% were heavily contaminated. Overall, the soil exhibited moderate pollution with a mean NSCI value of 2.55 (Table S5). Contamination of the soil sites may be associated with emissions from road traffic, marble processing plants (MPPs), and car repair workshops, which release pollutants into the surrounding atmosphere25,145,146. The results of the current study are consistent with the findings of Yüksel et al.147, Weissmannova et al.143, Wang et al.144 and Chen et al.148 who reported varying levels of soil and sediment contamination linked to industrial activities, traffic emissions, and other anthropogenic sources. These patterns align with established knowledge about the impact of such activities on environmental pollution, particularly in areas with similar sources of contamination.
The index of geoaccumulation (Igeo) originally proposed by Muller88 provides valuable insights into the severity, range, and abundance of heavy metal contents in the soils of the studied sites149. The Igeo results indicated that cadmium (Cd) had the highest value of 1.274, suggesting significant pollution, followed by arsenic (As) with a value of 0.663, indicating substantial contamination. Copper (Cu) and cobalt (Co), with Igeo values of 0.09 and 0.062, respectively, indicate moderate pollution. Vanadium (V), chromium (Cr), and nickel (Ni) exhibited relatively lower Igeo values of 0.061, 0.060, and 0.054, indicating comparatively lower contamination levels. Manganese (Mn), titanium (Ti), and strontium (Sr) had even lower Igeo values of 0.035, 0.026, and 0.024, suggesting minimal pollution. Likewise, Zinc (Zn) and zirconium (Zr) exhibited the lowest Igeo values of 0.012 and 0.004, indicating very low levels of contamination. The Igeo values of As, Cu, Co, V, Cr, Ni, Mn, Ti, Sr, Zn, and Zr categorized the soil sites as slightly to moderately polluted. The Igeo values between 1 and 2 for Cd are indicative of moderate pollution at soil sites S11 and S14–S16, whereas Igeo values between 2 and 3 at soil sites S12–S13 and S18–S21 are indicative of moderate to strong pollution. The Igeo values of Cd between 0 and 1 in the soil samples (S1–S10, S17) indicate unpolluted to moderate pollution. Similarly, the Igeo values of As (1 ≤ Igeo ≤ 2) indicate a moderate level of pollution at soil sites S13 and S19–S21. On the other hand, Igeo values of As between 0 and 1 classify the soil sites (S1–S12 and S14–S18) as uncontaminated to moderately contaminated (Table S6). The pollution pattern indicates that Cd and As are the primary contaminants of concern, with their levels varying from unpolluted to highly contaminated, depending on the site. Other metals, such as Cu, Co, V, Cr, and Ni, contribute to moderate pollution, while Mn, Ti, Sr, Zn, and Zr pose minimal environmental risk. The variation in Igeo values across different sites highlights the spatial heterogeneity in soil contamination, suggesting that targeted remediation efforts may be needed in specific areas with relatively high pollution levels144.
The potential ecological risk index provides an accurate assessment of heavy metal contamination in agricultural soil by considering metal toxicity relative to their reference concentrations in the Earth’s crust150. Table S7 presents the results of the ecological risk index (Eri) and the potential ecological risk index (PERI) assessments for heavy metals (HMs), showing substantial variation in the mean Eri values for individual heavy metals across the soil sites. The range of Eri values for Cd was (10–450) > As (2.5–65.5) > Cu (0.72–4.5) > Co (0.53–2.89) > Ni (0.92–1.99) > V (0.25–0.91) > Cr (0.29–0.89) > Mn (0.08–0.37) > Zn (0.01–0.2) (Table S7). The Eri values (< 15) for Cu, Co, Ni, V, Cr, Mn, and Zn at the soil sites indicated a low ecological risk associated with these metals. The Eri values for As, ranging from 15 to 30, indicate a moderate ecological risk at soil sites S2-S3, S6-S9, and S14-S18. Eri values between 30 and 60 signify considerable risk at sites S4-S5, S10-S13, and S19-S20. Site (S21) poses a high ecological risk for As (60 ≤ Eri < 120), whereas site (S1) is associated with a low ecological risk (Eri < 15). Furthermore, considerable ecological risk of Cd was characterized for S6 and S8-S10, while a high ecological risk was characterized for S4–S5, S7, and S17. Very high ecological risk was identified at sites S11–S16 and S18–S21, whereas sites S1–S3 exhibited lower ecological risk for Cd (Table S7). The current findings suggest that cadmium is the most concerning heavy metal, with several sites facing very high ecological risks, followed by arsenic, which presents moderate to high risks at various sites. Other metals such as Cu, Co, Ni, V, Cr, Mn, and Zn, generally pose low ecological risks across the studied areas. These results are strongly supported by the studies of Din et al.149, Chen et al.151 and Ajah et al.152, who also highlighted the significant ecological risks associated with cadmium and arsenic, while noting lower risks from other metals.
The PERI values reflect how different biological communities react to toxic substances, highlighting the potential ecological risks associated with hazardous elements153. The potential ecological risk values in the analyzed soil samples ranged from 18.52 to 511.35, with a mean value of 217.49 (Table S7). Almost 42.85% of the soil sites were at very high ecological risk, while 23.80% were at considerable ecological risk. In addition, approximately 14.28% of the sites were associated with moderate ecological risk, and 14.28% were classified as low ecological risk sites. The main contributor to the PERI appeared to be Cd (83.8%), followed by As (14.4%) (Fig. 2). In contrast, the cumulative contribution of Co, Mn, V, Cr, Zn, Cu and Ni to the PERI was 3.0%. The contribution rate of heavy metals to the PERI correlates not only with their concentration but also with their toxicity response factors138. In the current study, the pattern observed highlights Cd and As as the dominant hazardous elements in the soil, which pose the greatest ecological threats owing to their high toxicity and significant presence154. This finding is consistent with the results of Lu et al.155, Long et al.156 and Kumar et al.157, who also identified Cd and As as major contributors to potential ecological risks, especially in soils and sediments affected by human activities like mining. This study underscores the need for increased attention to Cd and As levels in agricultural soils because of their substantial ecological impact. The results also align with broader research, such as studies on sediment contamination in large freshwater lakes in China, where Cd was found to pose a moderate national-level risk due to significant emissions from human activities158.
Contribution percentage (%) of heavy metals to potential ecological risk.
A major global concern is human exposure to heavy metals due to environmental pollution. Heavy metals such as As, Cd, Cr, Pb, and Ni are highly toxic, classified as priority metals of significant public health concern159. Even at minimal exposure levels, they cause damage to multiple organs and are classified as human carcinogens160. Therefore, this study presents a health risk assessment of heavy metals for individuals working in marble processing plants (MPPs) and in surrounding agricultural soils, addressing both potential non-carcinogenic and carcinogenic health risks.
Table 2 presents the average daily dose (ADD, mg/kg/day), hazard quotient (HQ), and hazard index (HI) for non-carcinogenic risks of heavy metals (HMs) through ingestion, inhalation, and dermal contact for both children and adults. A comparatively greater trend in ADD values (mg/kg/day) was observed for children than for adults for both the ingestion and inhalation routes (Table 2). The ADDingestion values (mg/kg/day) for the children decreased in the order of Al (1.40E-01) > Fe (6.01E-02) > Mn (1.91E-03) > V (5.01E-04) > Cr (3.44E-04) > Cu (2.57E-04) > Ni (2.35E-04) > Co (7.49E-05) > Zn (7.37E-05) > As (6.76E-05) > Zr (3.59E-05) > Cd (2.44E-05). Compared to children, adults exhibited similar patterns of average daily dose (ADD) for ingestion and inhalation of the studied heavy metals, but with lower intake levels. Conversely, higher values of ADDdermal were obtained for adults than for children in the pattern of Al > Fe > As > Mn > V > Cr > Cu > Ni > Zn > Zr > Co > Cd (Table 2). According to this study, children appear to be more susceptible to health risks than adults are through the ingestion and inhalation routes. On the other hand, adults are more susceptible to such risks through the dermal route.
The hazard quotient (HQ) values of Al, As, Cd, Cr, Co, Ni, Cu, V, Fe, Mn, Zr, and Zn for adults and children were < 1.0 (Table 2). Therefore, no apparent health risk from exposure to a single HM is indicated. However, in children, the cumulative HI value from ingestion, inhalation and dermal exposure to heavy metals was > 1.0 (Table 2), indicating that prolonged exposure of children to soils near the marble industry is harmful to their health25. In addition, the probability distributions of total hazard index (HItotal) from ingestion, inhalation, and dermal exposure to heavy metals for adults (Fig. 3a) and children (Fig. 3b) were assessed via a probabilistic approach with Monte Carlo simulation. HI total scores were found to be less than 1.0 at the mean, 5th, and 95th percentiles for both children and adults, indicating no apparent risk for the target population. The deterministic approach estimated higher HItotal scores for both children and adults (Table 2) compared to the probabilistic approach shown in Fig. 3a,b. Similar results were also recorded by Sanaei et al.114, where the HI values calculated through the deterministic approach were higher than those from probabilistic values. This disparity can be attributed to the overestimation of risk caused by the inherent uncertainties of deterministic analysis115.
Probability distribution of total hazard index (HI) for soil heavy metals presented in Tabe 2 for (a) adults and (b) children from ingestion, inhalation and dermal contact route.
Cancer risk represents the probability of developing cancer from prolonged exposure to carcinogens161. The potential cancer risk from ingestion, inhalation, and dermal exposure to soil heavy metals (As, Cd, Co, Cr, and Ni) for adults and children is shown in Table 3. Commonly referenced acceptable risk levels range from 1.0E-06 to 1.0E-04 (i.e., one in 1,000,000 to one in 10,000)162. This study follows the USEPA162 guidelines, which define carcinogenic risk levels at 1.0E-04 and 1.0E-06, and considers these levels as potential indicators of adverse health risks for the exposed population. Through deterministic approach, the total carcinogenic risk (TCR) from soil heavy metal exposure ranged from 5.6E-10 to 6.9E-04 for children and from 1.2E-10 to 7.4E-05 for adults (Table 3). The carcinogenic risk from Cr, Ni, Cd and As ingestion was slightly higher in children than in adults and exceeded the USEPA limits of 1.0E-06 and 1.0E-04. The carcinogenic risk from inhalation and dermal exposure to these elements in children and adults was below the USEPA limits (Table 3). However, the TCR for Cr, Ni, Cd, and As was consistently higher in children than in adults across all exposure pathways (Table 3). Furthermore, the contribution of each element to the TCR decreased in the order of Cr > Ni > Cd > As, suggesting greater potential for the carcinogenic risk of Cr. The TCR for Co fell within the safety limits of 1.0E-04, indicating no apparent carcinogenic risk. Table 3 shows that Cr and Ni contribute significantly to the carcinogenic risk in children, and consequently, children are at higher risk than adults.
In addition to point estimation of carcinogenic risk, a probabilistic approach was applied to assess the carcinogenic risk in the exposed groups. This approach considers the proper distribution of parameters, including the metal concentration, ingestion rate (IR), exposure duration (ED), exposure frequency (EF), and body weight (BW) (Table S3). The simulated carcinogenic risk results for As, Cr, Cd, Ni, and Co from ingestion, inhalation, and dermal exposure for adults and children are shown in Figs. 4, 5, 6, 7, 8, 9, 10 and 11. In adults, the probabilistic cancer risk values from ingestion for As (Fig. 4a), Cr (Fig. 4b), and Ni (Fig. 4c) at the 5th, mean and 95th percentiles exceeded the minimum acceptable limit of 1.0E-06. The Cd level at the 95th percentile from the ingestion route also exceeded the minimum threshold of 1.0E-06 (Fig. 4d). The probabilistic cancer risk from dermal exposure to As (Fig. 5a) and Cr (Fig. 5b) in adults was slightly higher than from ingestion at the mean, 5th, and 95th percentiles. In adults, the potential cancer risk for As (Fig. 6a), Cd (Fig. 6b), Cr (Fig. 6c), Co (Fig. 6d), and Ni (Fig. 6e) via inhalation route was below the USEPA limits of 1.0E-04 and 1.0E-06. However, the total cancer risk for adults (Fig. 7d) from ingestion (Fig. 7a), inhalation (Fig. 7b), and dermal contact (Fig. 7c) route at the 95th percentile, exceeded the upper limit of 1.0E-04. This finding indicates a potential cancer risk for the adult group.
Shows the probability of cancer risk for (a) As (b) Cr (c) Ni (d) Cd in adults via ingestion route.
Probability of cancer risk for (a) As (b) Cr and (c) Ni in adults through dermal route.
Probability distribution of cancer risk for (a) As (b) Cr (c) Co (d) Cd and (e) Ni in adults through inhalation route.
Cumulative probability of cancer risk for adults from (a) ingestion, (b) inhalation (c) dermal and (d) total cancer risk from the three pathways.
Cumulative probability of cancer riskk for (a) As (b) Cr (c) Cd (d) Ni in children through ingestion route.
Probability distribution of cancer risk for (a) As (b) Cr (c) Cd (d) Co (e) Ni in children through inhalation route.
Cumulative probability of cancer risk for (a) As (b) Cr (c) Ni in chidren through dermal route.
Cumulative probability of cancer risk for children from (a) ingestion (b) inhalation (c) dermal route and (d) total cancer risk from the three exposure pathways.
For children, at the 95th percentile, the probabilities of cancer risk resulting from soil ingestion exposure to As (Fig. 8a), Cr (Fig. 8b) Cd (Fig. 8c), and Ni (Fig. 8d) were 1.83E-04, 2.52E-04, 1.99E-05, and 5.27E-04, respectively. These values exceeded the safety limit of 1.0E-04 and were relatively higher than the corresponding risks for the target group. Similarly, the probability of cancer risk due to inhalation exposure to soil Cr at the mean and 95th percentiles (Fig. 9b) was greater in children than in adults and exceeded the threshold of 1.0E-04. In addition, the probability of cancer risk in children due to inhalation exposure to soil As (Fig. 9a), Cd (Fig. 9c), Co (Fig. 9d) and Ni (Fig. 9e) at the mean and 95th percentiles was within the permissible limits of USEPA. In children, the potential cancer risks of As (Fig. 10a) and Cr (Fig. 10b) through the dermal pathway were lower than in adults but exceeded the lower safety limits of 1.0E-06. For Ni, the dermal exposure risk both in children (Fig. 10c) and adults (Fig. 5c) was within the permissible limits. In addition, the total cancer risk in children (Fig. 11d) from the combined pathways of ingestion (Fig. 11a), inhalation (Fig. 11b) and dermal contact (Fig. 11c) at the mean, 5th, and 95th percentiles were 6.95E-04, 5.49E-04, and 8.62E-04, respectively, and were above the upper safety limit of 1.0E-04, suggesting a greater carcinogenic hazard for children than adults. The slightly elevated risk in children can be attributed to their higher intake rates relative to body weight, increased absorption efficiency, and developmental susceptibility163. The exceedance of USEPA thresholds suggests that the levels of these metals pose a significant carcinogenic threat, necessitating urgent measures to mitigate exposure, particularly in vulnerable populations such as children164. Wu et al.153 also reported a high carcinogenic risk for children in the soil around electronics manufacturing sites.
In a previous study, the carcinogenic risk higher than 1.0E-04 through soil ingestion for adults and children was shown to predict health hazards for the target group165. In the current study, the slightly higher risk value in adults through the dermal pathway aligns with findings from a similar study conducted in the mining areas of central China165. Similarly, the cancer risk via the inhalation route was negligible for adults (Fig. 7b) and unacceptable at the 95th percentile for children (Fig. 11b). Cao et al.166 reported that consuming food grown on contaminated soil and direct skin contact with soil are significant human exposure pathways. Moreover, industrial and agricultural activities have led to the significant discharge of heavy metals into the surrounding soil. In accordance with Wang et al.168 and Yang et al.65, the use of pesticides and fertilizers significantly contributes to the occurrence of As and Cd in soil. Moreover, polluted water, smoke, and metal contaminated dust from smelting and metal processing can migrate into agricultural soils and increase the risk of cancer development65.
In the current study, the carcinogenic risks from As, Cr, Cd, and Ni are probably linked to the combustion of fossil fuel, inorganic fertilizers, battery workshops, and the emission of marble dust17,22,25,167. In addition, excessive accumulation of Cr in human tissues triggers gastric and lung cancer168. Similarly, exposure to nickel can produce a variety of human health side effects, such as pulmonary fibrosis, nasal and lung cancer, and kidney and cardiovascular diseases167. Similarly, soil arsenic exposure significantly associated with mortality rates from colon, stomach, kidney, lung and nasopharyngeal cancers169. According to Suwatvitayakorn et al.170, Cd has the potential to accumulate in the human body, resulting in acute and chronic effects on various organs and systems, with a particular impact on the kidneys and skeletal system. Therefore, it is equally important to carefully treat and monitor extremely noxious pollutants such as Cd, Cr, As, and Ni, even at very low concentrations in soil, air, and water.
Despite the carcinogenic or noncarcinogenic risk, children are more vulnerable to the potential health risk attributed to the existence of HMs in the soils surrounding MPPs. This suggests that children are more susceptible to adverse health effects from exposure to HMs, as they are more likely to ingest HMs orally from hand to mouth171. Among the different routes of exposure, the ingestion route contributes the most to the overall cancer risk in both children and adults, followed by the dermal route, and the inhalation route contributes very little to the overall cancer risk. Similar results were also reported by Saha et al.172, Kusin et al.171, Cai et al.173 and Hao et al.165.These findings strongly support our current results. Our study underscores the significant health risks from heavy metal exposure in agricultural soils near industries, emphasizing the need for focused risk control measures, as also suggested by Butt et al.48, Khan et al.25, Muhammad et al.174 and Noreen et al.128 in their respective studies.
A sensitivity analysis was performed to determine the importance of the variables involved in calculating noncarcinogenic and carcinogenic risk175. The total hazard index (HI) < 1 at the mean, 5th, and 95th percentiles for both children (Fig. 3b) and adults (Fig. 3a) suggests that there is no apparent noncancerous risk for these two groups. However, sensitivity analysis of the total HI in adults indicated that heavy metal concentration (60.4%), exposure duration (30.4%) and exposure frequency (8.4%) were the most predictive noncancer risk factors (Fig. 12b).These results are in agreement with those of the study conducted by Fallahzadeh et al.175. Similarly, for children, exposure duration (67.3%), heavy metal concentration (26.5%) and exposure frequency (6.4%) were important factors for noncancer risk (Fig. 12a). Sanaei et al.114 reported that the ingestion rate is the main factor influencing health risk, followed by exposure duration and metal concentration, while body weight is a less influential variable, partially supporting our study results. In addition, the assessment of the contributions of HMs to adult cancer risk through ingestion, inhalation, and dermal contact (Fig. 13b) indicated that As had the most significant impact, accounting for 83.6% of the total risk, followed by Cr (12.4%) and Ni (9%). For children, greatest contribution to the total cancer risk was simulated for Ni (59.5%), followed by Cr (21.6%) and As (18.3%) from the ingestion, inhalation and dermal exposure routes (Fig. 13a). The findings of this study align with those of Chen et al.169, Yang et al.65 and Mallongi et al.66, who identified As, Cd, Cr, and Ni as the main contributors to total risk variance.
Sensitivity analysis of HI total for (a) chidren (b) adults.
Sensitivity analysis of total cancer risk for (a) chidren (b) adults.
A normalized variation matrix used to show the linear association for the chemical components42. The analysis of the normalized variation matrix (Table 4) for the chemical parts (Table 1) reveals key associations between soil pH, EC, TDS, and the studied elements in agricultural soils near marble industries. These associations provide insights into the geochemical behavior of these elements and their interactions within the soil matrix176. According to Boente et al.42 variations less than 0.2 suggest a proportional or linear relationship, while variations greater than 1.0 indicate a lack of proportionality or linear connection. In the current study, strong linear relationships were observed between soil pH and elements such as Ca (0.04), Ni (0.04), Mg (0.07), Cr (0.10), Al (0.12), V (0.10), Sr (0.12), K (0.13), As (0.14) and Mn (0.18). These values suggest a proportional or linear relationship, indicating that changes in soil pH significantly affect the solubility and availability of these elements177. Therefore, soil pH plays a crucial role in determining the solubility and mobility of elements178. The very strong relationship of calcium with pH indicates that calcium carbonate plays a crucial role in buffering the pH of the soil and stabilizing pH levels179. However, acidic conditions often increase the solubility of these elements which can lead to higher concentrations in water and soil180. EC is an important indicator of soil salinity and nutrient availability181. The strong linear relationships of EC with Cr (0.07), Mg (0.09), V (0.09), Ni (0.11), Ca (0.12), Sr (0.14), K (0.15), As (0.15) and Al (0.17) (Table 4) suggest the impacts of nutrient availability and soil chemical processes182. The relationships of EC with Ca, K, Mg, and Sr contribute significantly to the ionic balance, reflecting their role in soil structure and maintaining soil health and fertility182,183. The presence of these elements in a balanced proportion ensures that the soil remains fertile and capable of supporting healthy plant growth184. However, an imbalance, often indicated by a high EC, could suggest salinity issues, which can hinder plant growth and lead to soil degradation over time185. The TDS is a measure of the total concentration of dissolved substances in soil water, which directly impacts soil water quality and its suitability for plant growth186. The linear proportions of TDS with Ni (0.15) and V (0.16) highlight how increases in TDS can affect soil and water quality187. Elevated TDS levels can indicate the presence of excessive salts or other dissolved substances, which may reduce the availability of water to plants, potentially leading to water stress and reduced plant growth188,189. The associations between pH, EC, TDS, and the studied elements provide valuable insights into the potential sources of these elements and their environmental implications12,26. For example, the strong relationship between pH and Ca suggests that calcium carbonate from marble industries could be a significant source, contributing to the observed soil chemistry23. The relationships involving EC and TDS with elements such as Ni, Cr, and V may indicate industrial or agricultural inputs, such as fertilizers or industrial emissions, as potential sources190. Similarly, values below 0.2 for Al versus K (0.07), Mn (0.12), Ni (0.17) and Cr (0.18) (Table 4) indicate that these element pairs are considered proportional and can be easily explained by the fact that where Al is present, it is likely to be accompanied by relatively high concentrations of K, Mn, Ni and Cr. Furthermore, the smaller contribution values for pairs such as Cr to V (0.09), Sr (0.14), Ni (0.16), and K (0.17), along with Sr to V (0.12), Ca to Ni (0.17), and K to Sr (0.18), underscore the existence of a linear binary relationship among these elements. This suggests that the concentration of one element in a pair can strongly influence the concentration of the other. For example, the linear relationship between Cr and V or Sr implies that fluctuations in chromium levels could directly impact vanadium or strontium concentrations. Similarly, the associations of Ca with Ni and K with Sr suggest a strong geochemical linkage between these pairs, likely due to co-deposition or similar environmental factors that influence their distribution. The observed associations may also indicate shared sources or pathways of these elements, which is consistent with the patterns identified in other studies by Boente et al.42, Khan et al.23; and Khorshidi et al.191.
In addition, appropriate measures to find sources of HM contamination in agricultural soils are of paramount importance for developing sustainable management approaches and protecting soil quality192. Therefore, a clr-biplots (Fig. 14a, b) based on a centered logarithmic transformation (clr) was constructed to interpret the inter-elemental associations and source identification of HMs in agricultural soils42,193. Figure 14a explains 62% of the variation between PC1 and PC2 and Fig. 14b explains 53% of the variation between PC1 and PC3 of the studied elements, indicating multiple element relationships related to natural geology and anthropogenic activities194. Additionally, strong element associations between Mn-P-Ti, Al, and K are explained by the vicinity of their vertices and their rays pointing in a similar direction, indicating a combination of geogenic and anthropogenic origins195. The vertices of Mn, P, and Ti are superimposed, and Ti was found to have the highest vector length, which could be associated with the co-precipitation of TiO2, MnO, and P2O5 found in marble powder and fertilizers25. Similarly, Cu, V, Sr, Cr, Mg, Ni, and Co form another association due to the close proximity of their vertices and the alignment of their rays in the same direction, as illustrated in Fig. 14a,b. This association is likely linked to vehicle emissions and industrial activities194. A close relationship between Zn, As, Cd, and Zr forms a separate group and could be associated with smelter activities and the regular use of inorganic fertilizers196,197. Fe and Si form another association because their rays also point in the same direction. However, in Fig. 15, Fe, Mn, Ti, and Ca shifted their positions and entered the elemental associations of Zn–Cd–As–Zr–Fe, Cu–Ti–Cr, Al–K–Mn, and Mg–Ca–Sr–Ni, indicating their affinity for these elements and sharing both natural and anthropogenic sources.
clr-biplot (a) PC1 vs PC2 (b) PC1 vs PC3 of the elements in study area.
CoDa-dendrogram shows association of the data set of the geochemical composition of the soil generated by the clustering algorithm (Wards method).
Factor analysis (Varimax rotation) of the centered log-ratio transform (clr) dataset was applied to discriminate potential geochemical source patterns and the associations and interrelationships between variables in the topsoil191. The three-factor model (Eigen > 1) explained 70.32% of the variance overall, with the contributions of each factor being F1 (30.76%), F2 (27.11%), and F3 (12.45%) (Table 5). A variable with factor loadings greater than 0.5 was included to explain the main components of each factor. The three-factor model captured a good correlation between the variables and their potential sources, showing a communality greater than 0.6 (60% variability). Factor one (F1) was dominated by Al, Mn, K, and Ti with strong positive loadings of 0.89, 0.80, 0.75, and 0.68, respectively, which could be attributed mainly to road transport and industrial activities198. For factor two (F2), the Cr, V, Cu, and Sr loadings were 0.92, 0.84, 0.80, and 0.78, respectively, indicating that the marble industry was the source of origin25. Factor three (F3) showed strong positive loadings, mainly with Ni (0.81), Mg (0.77), and Co (0.71), suggesting that river water could be another important source of pollution145.
Cluster analysis (CA) is commonly used with PCA to verify the results of associations and HMs sources193. For this purpose, the coda dendrogram (Fig. 15) is chosen as a visual method to identify possible relationships, pairwise or involving larger groups199. The short vertical bars in the coda dendrogram correspond to linear associations between the clustering variables. Therefore, three groups of elements can be identified: (i) Al, K, Mn, P, Mn, Si, Fe, and Ti, which are not highly toxic and are of geogenic origin, and (ii) Cr, Ni, V, Cu, Co, Mg, Ca, and Sr, which become toxic in excess quantity, indicating anthropogenic and lithogenic sources200. The third group consists of As, Cd, Zn, and Zr, which are potentially toxic and derived from the extensive application of chemical fertilizers197. According to Weissmannová et al.143, different elements grouped in the same cluster originate from a common source. The results of the coda dendrogram (Fig. 15) fully agreed with the results of the clr-biplots (Fig. 14a,b), variation matrix analysis (Table 4), and factor analysis (Table 5). These observations reinforce that industrial effluents, intensive fertilizer use, vehicle emissions, and irrigation with contaminated water are key sources of toxic elements, verifying findings from Jehan et al.145, Khan et al.25, and Sarwar et al.201.Therefore, these findings strongly support our current study results (Fig. 16).
Pollution sources of heavy metals and potential health risk assessment.
In the current study, the concentration of macro elements and heavy metals in soil near marble processing plants (MPPs) in district Malakand, was assessed to determine the extent of contamination and potential health risks. The concentrations of macro elements and heavy metals in the soil decreased in the following order: Al > Mg > Fe > K > Ti > P > Mn > V > Sr > Cr > Cu > Ni > Si > Co > Zn > As > Zr > Cd. The maximum concentrations of Ca, P and Cd exceeded the acceptable limits of the FAO and the world average shale. The Pollution Index (Pi) and Nemerow Synthetic Contamination Index (NSCI), indicated different levels of contamination from slight to heavy contamination. The Geoaccumulation Index (Igeo) indicated significant contamination of the soil sites by Cd and As, moderate contamination by Cu, Co, V, Cr, and Ni, and minimal contribution from Mn, Ti, Sr, Zn, and Zr. The potential ecological risk assessments revealed that Cd poses the highest ecological risk, followed by As. The pollution indices and compositional exploratory data analysis revealed that the spatial variability of soil contamination is primarily influenced by factors including industrial activities, ro ad traffic, and marble processing plants. The health risk assessments indicated elevated non-carcinogenic risks for children, with the total cancer risk fromthe potential toxic metals above safe levels, particularly for children. The key factors influencing health risks included metal concentration, exposure duration, and frequency, with As and Cr being significant contributors to cancer risk. The long-term monitoring of heavy metal contamination and its health impacts should be implemented, enforcing regulations on emissions from industrial and vehicular sources should be enforced, and soil remediation techniques such as phytoremediation and biochar application should be employed to reduce contamination levels. Additionally, public health programs should be developed to limit children’s exposure to contaminated soils, and alternative agricultural practices and land-use policies should be explored to minimize exposure and risks.
All the data are included in the article.
Wu, J., Lu, J., Li, L., Min, X. & Luo, Y. Pollution, ecological-health risks, and sources of heavy metals in soil of the northeastern Qinghai-Tibet Plateau. Chemosphere 201, 234–242 (2018).
Article ADS CAS PubMed Google Scholar
Abrahams, P. W. Soils: Their implications to human health. Sci. Total Environ. 291, 1–32 (2002).
Article ADS CAS PubMed Google Scholar
Yüksel, B. et al. Appraisal of metallic accumulation in the surface sediment of a fish breeding dam in Türkiye: a stochastical approach to ecotoxicological risk assessment. Mar. Pollut. Bull. 203, 116488 (2024).
Fazekašová, D. & Fazekaš, J. Soil quality and heavy metal pollution assessment of iron ore mines in Nizna Slana (Slovakia). Sustain. 12, (2020).
Hussain, C. M. & Keçili, R. Environmental pollution and environmental analysis. Mod. Environ. Anal. Tech. Pollutants https://doi.org/10.1016/b978-0-12-816934-6.00001-1 (2020).
Topaldemir, H., Taş, B., Yüksel, B. & Ustaoğlu, F. Potentially hazardous elements in sediments and Ceratophyllum demersum: An ecotoxicological risk assessment in Miliç Wetland, Samsun Turkey. Environ. Sci. Pollut. Pic. 30, 26397–26416 (2023).
Aradpour, S. et al. Metal contamination assessment in water column and surface sediments of a warm monomictic man-made lake: Sabalan Dam Reservoir, Iran. Hydrol. Res. 51, 799–814 (2020).
Brevik, E. C. et al. Soil and human health: current status and future needs. Air, Soil Water Res. 13, 117862212093444 (2020).
Nodefarahani, M. et al. Metal pollution assessment in surface sediments of Namak Lake, Iran. Environ. Sci. Pollut. Res. 27, 45639–45649 (2020).
Sarker, A. et al. Heavy metals contamination and associated health risks in food webs: A review focuses on food safety and environmental sustainability in Bangladesh. Environ. Sci. Pollut. Res. 29, 3230–3245 (2022).
Ustaoğlu, F., Yüksel, B., Tepe, Y., Aydın, H. & Topaldemir, H. Metal pollution assessment in the surface sediments of a river system in Türkiye: Integrating toxicological risk assessment and source identification. Mar. Pollut. Bull. 203, 116514 (2024).
Virto, I. et al. Soil Degradation and soil quality in Western Europe: Current situation and future perspectives. Sustainability 7, 313–365 (2014).
Prăvălie, R. Exploring the multiple land degradation pathways across the planet. Earth-Sci. Rev. 220, 103689 (2021).
Timmis, K. & Ramos, J. L. The soil crisis: the need to treat as a global health problem and the pivotal role of microbes in prophylaxis and therapy. Microb. Biotechnol. 14, 769–797 (2021).
Article PubMed PubMed Central Google Scholar
Rojas, R. V., Achouri, M., Maroulis, J. & Caon, L. Healthy soils: A prerequisite for sustainable food security. Environ. Earth Sci. 75, 180 (2016).
Ahmad, Z., Khan, S. M., Ali, M. I., Fatima, N. & Ali, S. Pollution indicandum and marble waste polluted ecosystem; role of selected indicator plants in phytoremediation and determination of pollution zones. J. Clean. Prod. 236, 117709 (2019).
Khoiroh, U. Marble dust exposure relationship to workers’ lung conditions in marble industries. J. Kesehat. Lingkung. 12, 285 (2020).
Shqiarat, M., Alnawafleh, H. & Abudanah, F. Characterization and provenance study of marble from Udhruh, Southern Jordan. Mediterr. Archaeol. Archaeom. 19, 1–8 (2019).
Tazzini, A., Gambino, F., Casale, M. & Dino, G. A. Managing marble quarry waste: Opportunities and challenges for circular economy implementation. Sustainability 16, 3056 (2024).
Mulk, S., Azizullah, A., Korai, A. L. & Khattak, M. N. K. Impact of marble industry effluents on water and sediment quality of Barandu River in Buner District, Pakistan. Environ. Monit. Assess. 187, 8 (2015).
Ali, S. et al. Effect of marble mines on physichochemical properties of soil in District Buner. 3, 312–318 (2017).
Khan, J., et al. Burden of marble factories and health risk assessment of kidney (renal) stones development in district Buner, Khyber Pakhtunkhwa, Pakistan. Expert Opin. Environ. Biol. (2014).
Khan, A. et al. Irrigation suitability, health risk assessment and source apportionment of heavy metals in surface water used for irrigation near marble industry in Malakand. Pakistan. PLoS One 17, 1–26 (2022).
Neelab, et al. Effects of marble industry effluents on soil quality, growth and productivity of tomato (Lycopersicon esculentum) in district Swat, Khyber Pakhtunkhwa, Pakistan. Polish J. Environ. Stud. 31, 3263–3270 (2022).
Khan, A., Khan, M. S., Hadi, F., Khan, A. N. & Saddiq, G. Energy-dispersive X-ray (EDX) Fluorescence Based Analysis of Heavy Metals in Marble Powder, Paddy Soil and Rice (Oryza sativa L.) With Potential Health Risks in District Malakand, Khyber Pakhtunkhwa, Pakistan. Environ. Pollut. Bioavailab. (2021) https://doi.org/10.1080/26395940.2021.1986427.
Khan, A., Khan, M. S., Shafique, M. A., Khan, Q. & Saddiq, G. Assessment of potentially toxic and mineral elements in paddy soils and their uptake by rice (Oryza sativa L.) with associated health hazards in district Malakand, Pakistan. Heliyon 10, e28043 (2024).
Azizullah, A., Khattak, M. N. K., Richter, P. & Häder, D.-P. Water pollution in Pakistan and its impact on public health: A review. Environ. Int. 37, 479–497 (2011).
Article CAS PubMed Google Scholar
Khan, M., Chaudhry, M. N., Ahmad, S. R., Saif, S. & Mehmood, A. Performance of EIA authority and effectiveness of EIA system in Pakistan. Environ. Impact Assess. Rev. 81, 106357 (2020).
Migliore, M., Lavagna, M. & Talamo, C. Strategies for the environmental management of the marble districts through the integrated contribution of life cycle assessment and Gis. Ital. J. Eng. Geol. Environ. 2, 79–87 (2017).
Deng, W. et al. Adsorption recovery of phosphate from waste streams by Ca/Mg-biochar synthesis from marble waste, calcium-rich sepiolite and bagasse. J. Clean. Prod. 288, 125638 (2021).
Rostocki, A., Unyay, H., Ławińska, K. & Obraniak, A. Granulates based on bio and industrial waste and biochar in a sustainable economy. Energies 16, 56 (2022).
Nafees , M. , Saeed , T. , Zuljajal , F. , Akbar , H. & Rehman , A. Marble production and environmental pressure : A case study of Hayatabad industrial estate methodology (2023) https://doi.org/10.21203 /rs.3.rs-3244303/v1.
Aukour , F. & Al-Qinna , M. Marble production and environmental constraints: Case study from Zarqa Governorate, Jordan. In (2008).
Salem, H. S., Yihdego, Y. & Muhammed, H. H. The status of freshwater and reused treated wastewater for agricultural irrigation in the occupied palestinian territories. J. Water Health 19, 120–158 (2021).
Noreen, U. et al. Water pollution and occupational health hazards caused by the marble industries in district Mardan, Pakistan. Environ. Technol. Innov. 16, 100470 (2019).
Said, A. K., & Muhammad Suleman, M. A. Assessment of pollution load in marble waste water in Khairabad, District. 8, 35–39 (2017).
Mulk, S., Korai, A. L., Azizullah, A., Shahi, L. & Khattak, M. N. K. Marble industry effluents cause an increased bioaccumulation of heavy metals in Mahaseer (Tor putitora) in Barandu River, district Buner, Pakistan. Environ. Sci. Pollut. Res. Int. 24, 23039–23056 (2017).
Article CAS PubMed Google Scholar
Zulfiqar, U. et al. Cadmium toxicity in plants: Recent progress on morpho-physiological Effects and Remediation Strategies. J. Soil Sci. Plant Nutr. 22, 212–269 (2022).
Nagajyoti, P. C., Lee, K. D. & Sreekanth, T. V. M. Heavy metals, occurrence and toxicity for plants: A review. Environ. Chem. Lett. 8, 199–216 (2010).
Rahman, Z. & Singh, V. P. The relative impact of toxic heavy metals (THMs) (arsenic (As), cadmium (Cd), chromium (Cr)(VI), mercury (Hg), and lead (Pb)) on the total environment: an overview. Environ. Monit. Assess. 191, 419 (2019).
Keshavarzi, A. & Kumar, V. Ecological risk assessment and source apportionment of heavy metal contamination in agricultural soils of Northeastern Iran. Int. J. Environ. Health Res. 29, 544–560 (2019).
Article CAS PubMed Google Scholar
Boente, C., Albuquerque, M. T. D., Gallego, J. R., Pawlowsky-Glahn, V. & Egozcue, J. J. Compositional baseline assessments to address soil pollution: An application in Langreo. Spain. Sci. Total Environ. 812, 152383 (2022).
Article CAS PubMed Google Scholar
Rai, P. K., Lee, S. S., Zhang, M., Tsang, Y. F. & Kim, K. H. Heavy metals in food crops: Health risks, fate, mechanisms, and management. Environ. Int. 125, 365–385 (2019).
Article CAS PubMed Google Scholar
Gabarrón, M., Faz, A. & Acosta, J. A. Use of multivariable and redundancy analysis to assess the behavior of metals and arsenic in urban soil and road dust affected by metallic mining as a base for risk assessment. J. Environ. Manage. 206, 192–201 (2018).
Wang, S. et al. Source and health risk assessment of heavy metals in soil–ginger system in the jing river basin of Shandong province, North China. Int. J. Environ. Res. Public Health 18, (2021).
Nduka, J. K., Kelle, H. I. & Amuka, J. O. Health risk assessment of cadmium, chromium and nickel from car paint dust from used automobiles at auto-panel workshops in Nigeria. Toxicol. Rep. 6, 449–456 (2019).
Article CAS PubMed PubMed Central Google Scholar
Anyimah, E., Isaac, A., Ofosu, W., Erick, H. & Godfred, L. Exposures and health risks associated with elements in diets from a gold mining area. Biol. Trace Elem. Res. https://doi.org/10.1007/s12011-021-02777-0 (2021).
Butt, I. M., Mustafa, T., Rauf, S., Razzaq, A. & Anwer, J. Pulmonary function parameters among marble industry workers in Lahore, Pakistan. F1000Research 10, 938 (2021).
Fawad, M. et al. Pollution hotspots and potential impacts on land use in the Mohmand Marble Zone, Pakistan. Environ. Earth Sci. 80, 1–12 (2021).
Alizadeh, M., Noori, R., Omidvar, B., Nohegar, A. & Pistre, S. Human health risk of nitrate in groundwater of Tehran-Karaj plain, Iran. Sci. Rep. 14, 7830 (2024).
Article ADS CAS PubMed PubMed Central Google Scholar
Pawlowsky-Glahn, V. & Buccianti, A. Compositional Data Analysis: Theory and Applications. Compositional Data Analysis: Theory and Applications (John Wiley & Sons, Ltd, 2011). https://doi.org/10.1002/9781119976462.
Fačevicová, K., Bábek, O., Hron, K. & Kumpan, T. Element chemostratigraphy of the Devonian/Carboniferous boundary: A compositional approach. Appl. Geochem. 75, 211–221 (2016).
Aitchison, J. The Statistical Analysis of Compositional Data. Monographs on Statistics and Applied Probability. (Chapman and Hall Ltd, London, UK, 1986).
Pawlowsky-Glahn, V. & Egozcue, J. J. Compositional data and their analysis: an introduction. Geol. Soc. London Spec. Publ. 264, 1–10 (2006).
Pawlowsky-Glahn, V., Egozcue, J. J. & Tolosana-Delgado, R. Modelling and Analysis of Compositional Data. (John Wiley & Sons, Ltd, 2015). https://doi.org/10.1002/9781119003144.
Coenders, G., Hlebec, V. & Kogovsšek, T. Measurement quality in indicators of compositions. A compositional multitrait-multimethod approach. Surv. Res. Methods 5, 63–74 (2011).
Egozcue, JJ, Pawlowsky-Glahn, V., Mateu-Figueras, G. & Barceló-Vidal, C. Isometric logratio transformations for compositional data analysis. Math. Geol. 35, 279–300 (2003).
Pawlowsky-Glahn, V. & Egozcue, J. J. Geometric approach to statistical analysis on the simplex. Stoch. Environ. Res. Risk Assess. 15, 384–398 (2001).
Templ, M., Filzmoser, P. & Reimann, C. Cluster analysis applied to regional geochemical data: Problems and possibilities. Appl. Geochem. 23, 2198–2213 (2008).
Article ADS CAS Google Scholar
Thiombane, M. et al. Exploratory analysis of multi-element geochemical patterns in soil from the Sarno River Basin (Campania region, southern Italy) through compositional data analysis (CODA). J. Geochem. Explor. 195, 110–120 (2018).
Xiao, R. et al. Accumulation, ecological-health risks assessment, and source apportionment of heavy metals in paddy soils: A case study in Hanzhong, Shaanxi, China. Environ. Pollut. 248, 349–357 (2019).
Article CAS PubMed Google Scholar
Ren, W. et al. Reconsidering brownfield redevelopment strategy in China’s old industrial zone: A health risk assessment of heavy metal contamination. Environ. Sci. Pollut. Res. 22, 2765–2775 (2015).
Islam, M. S., Ahmed, M. K. & Habibullah-Al-Mamun, M. Heavy metals in cereals and pulses: Health implications in Bangladesh. J. Agric. Food Chem. 62, 10828–10835 (2014).
Article CAS PubMed Google Scholar
Albanese, S. et al. Potentially toxic elements in the soils of Campi Flegrei (south Italy) and the immediate surroundings: Spatial distribution, origin and probabilistic human health risk. Chemosphere 313, 137297 (2023).
Article CAS PubMed Google Scholar
Yang, S. et al. Status assessment and probabilistic health risk modeling of metals accumulation in agriculture soils across China: A synthesis. Environ. Int. 128, 165–174 (2019).
Article CAS PubMed Google Scholar
Mallongi, A. et al. Health risk assessment of potentially toxic elements in Maros karst groundwater: a Monte Carlo simulation approach. Geom. Nat. Hazards Risk 13, 338–363 (2022).
Wang, C. Monte carlo simulation. Springer Ser. Reliab. Eng. https://doi.org/10.1007/978-3-030-62505-4_3 (2021).
Moazamnia, M. et al. Probabilistic human health risk assessment for arsenic, nickel and lead exposures based on two-dimensional Monte Carlo simulation. Groundw. Sustain. Dev. 27, 101312 (2024).
Yang, Z. et al. Spatial distribution and sources of seven available heavy metals in the paddy soil of red region in Hunan Province of China. Environ. Monit. Assess. (2018).
White, C. M. & Barbercheck, M. E. Managing Soil Health: Concepts and Practices. https://extension.psu.edu/managing-soil-health-concepts-and-practices (2012).
Aqeel, M., Jamil, M. & Yusoff, I. Soil Contamination, Risk Assessment and Remediation. in Environmental Risk Assessment of Soil Contamination vol. i 13 (InTech, 2014).
Hou, D., Bolan, N. S., Tsang, D. C. W., Kirkham, M. B. & O’Connor, D. Sustainable soil use and management: An interdisciplinary and systematic approach. Sci. Total Environ. 729, 138961 (2020).
Article CAS PubMed PubMed Central Google Scholar
Hussain, S. Species diversity of basidiomycetes of district Malakand. (PhD Dessertation, Department of Botany, Hazara University. pp 1–249.http://prr.hec.gov.pk/jspui/handle/123456789/7580, (2016).
Jan, S., Khan, K., Hameed, I. & Ahmad, N. Ethnobotanical studies of the medicinal plants of Malakand agency, Khyber Pakhtunkhwa. Pakistan J. Plant Sci. 18, (2012).
Khan, A., Khan, N. & Ali, K. An assessment of the floristic diversity, life-forms and biological spectrum of vegetation in Swat Ranizai, district Malakand, Khyber Pakhtunkhwa, Pakistan. Technol. Dev. 36, 61–78 (2017).
Iqbal, M., Khalid, N., Sajjad, M., Pervaiz, U. & Ullah, S. Farmers field schools and rice productivity: An empirical analysis of District Malakand. Sarhad J. Agric. 28, 143–149 (2012).
Malkani, M. S., Mahmood, Z., Alyani, M. I. & Siraj, M. Mineral Resources of Khyber Pakhtunkhwa and FATA, Pakistan. Geol. Surv. Pakistan, Inf. Release 996, 1–61 (2017).
Wroble, J., Frederick, T., Frame, A. & Vallero, D. Comparison of soil sampling and analytical methods for asbestos at the Sumas mountain asbestos site, working towards a toolbox for better assessment. PLoS One 12, 1–26 (2017).
Rafique, N. & Tariq, S. R. Distribution and source apportionment studies of heavy metals in soil of cotton/wheat fields. Environ. Monit. Assess. 188, (2016).
Khan, M. U., Malik, R. N. & Muhammad, S. Human health risk from Heavy metal via food crops consumption with wastewater irrigation practices in Pakistan. Chemosphere 93, 2230–2238 (2013).
Article ADS CAS PubMed Google Scholar
Ali, F., Ullah, H. & Khan, I. Heavy Metals Accumulation in Vegetables Irrigated with Industrial Influents and Possible Impact of such Vegetables on Human Health. Sarhad J. Agric. (2017).
Shi, C., Ding, H., Zan, Q. & Li, R. Spatial variation and ecological risk assessment of heavy metals in mangrove sediments across China. Mar. Pollut. Bull. 143, 115–124 (2019).
Article CAS PubMed Google Scholar
Usman, K., Al-Ghouti, M. A. & Abu-Dieyeh, M. H. The assessment of cadmium, chromium, copper, and nickel tolerance and bioaccumulation by shrub plant Tetraena qataranse. Sci. Rep. 9, 1–11 (2019).
Doabi, S. A., Karami, M., Afyuni, M. & Yeganeh, M. Pollution and health risk assessment of heavy metals in agricultural soil, atmospheric dust and major food crops in Kermanshah province, Iran. Ecotoxicol. Environ. Saf. 163, 153–164 (2018).
Article CAS PubMed Google Scholar
Liang, Y. et al. Heavy metal contamination and health risk assessment in the vicinity of a tailing pond in Guangdong, China. Int. J. Environ. Res. Public Health (2017).
Paetz, A. & Crmann, G. Problems and results in the development of international standards for sampling and pretreatment of soils. In Environmental Sampling for Trace Analysis 321–334 (Wiley-VCH Verlag GmbH, 1994). https://doi.org/10.1002/9783527615872.ch16.
Adebiyi, F. M. & Ore, O. T. Assessing highway-related metal pollution using surface soil and tree bark as indicators. Environ. Forensics (2021) https://doi.org/10.1080/15275922.2021.1976314.
Muller, G. Index of geoaccumulation in sediments of the Rhine River. GeoJournal 2, 108–118 (1969).
Muller, G. The heavy metal pollution of the sediments of Neckars and its tributary: A stocktaking. Chem Zeit 105, 157–164 (1981).
Hakanson, L. An ecological risk index for aquatic pollution control.a sedimentological approach. Water Res. 14, 975–1001 (1980).
Tomlinson, D. L., Wilson, J. G., Harris, C. R. & Jeffrey, D. W. Problems in the assessment of heavy-metal levels in estuaries and the formation of a pollution index. Helgoländer Meeresuntersuchungen 33, 566–575 (1980).
Lu, X., Wang, L., Lei, K., Huang, J. & Zhai, Y. Contamination assessment of copper, lead, zinc, manganese and nickel in street dust of Baoji, NW China. J. Hazard. Mater. 161, 1058–1062 (2009).
Article CAS PubMed Google Scholar
Turekian, K. K. & Wedepohl, K. H. Distribution of the elements in some Major Units of the Earth’s Crust. GSA Bull. 72, 175–192 (1961).
Lu, X., Pan, H. & Wang, Y. Pollution evaluation and source analysis of heavy metal in roadway dust from a resource-typed industrial city in Northwest China. Atmos. Pollut. Res. 8, 587–595 (2017).
Karimi, A. et al. Assessment of human health risks and pollution index for heavy metals in farmlands irrigated by effluents of stabilization ponds. Environ. Sci. Pollut. Res. 27, 10317–10327 (2020).
El Azhari, A., Rhoujjati, A., El Hachimi, M. L. & Ambrosi, J. Pollution and ecological risk assessment of heavy metals in the soil-plant system and the sediment-water column around a former Pb/Zn-mining area in NE Morocco. Ecotoxicol. Environ. Saf. 144, 464–474 (2017).
Weissmannová, H. D. & Pavlovský, J. Indices of soil contamination by heavy metals—methodology of calculation for pollution assessment (minireview). Environ. Monit. Assess. 189, 616 (2017).
Emenike, P. G. C. et al. An integrated assessment of land-use change impact, seasonal variation of pollution indices and human health risk of selected toxic elements in sediments of River Atuwara, Nigeria. Environ. Pollut. 265, (2020).
Lu, X. et al. Risk assessment of toxic metals in street dust from a medium-sized industrial city of China. Ecotoxicol. Environ. Saf. 106, 154–163 (2014).
Article CAS PubMed Google Scholar
Lu, X., Shi, D., Yin, N., Pan, H. & Smith, P. Risk assessment of heavy metals in finer than 63-μm dust particles from various functional areas in Xi’an, China. Air Qual. Atmos. Heal. 10, 907–915 (2017).
Ke, X. et al. Ecological risk assessment and source identification for heavy metals in surface sediment from the Liaohe River protected area, China. Chemosphere 175, 473–481 (2017).
Article ADS CAS PubMed Google Scholar
Adimalla, N. Heavy metals contamination in urban surface soils of Medak province, India, and its risk assessment and spatial distribution. Environ. Geochem. Health 42, 59–75 (2020).
Article CAS PubMed Google Scholar
Lian, M. et al. Profiles and potential health risks of heavy metals in soil and crops from the watershed of Xi River in Northeast China. Ecotoxicol. Environ. Saf. 169, 442–448 (2019).
Article CAS PubMed Google Scholar
USEPA. Risk Assessment Guidance for Superfund Volume I : Human Health Evaluation Manual ( PartE, Supplemental Guidance for Dermal Risk Assessment ) Final. (2004).
Ahmad, W. et al. Toxic and heavy metals contamination assessment in soil and water to evaluate human health risk. Sci. Rep. 11, 1–12 (2021).
Dessie , BK , Gari , SR , Mihret , A. , Desta , AF & Mehari , B. Determination and health risk assessment of trace elements in the tap water of two sub-cities of Addis Ababa . Ethiopia. Helium 7, e06988 (2021).
Article CAS PubMed Google Scholar
Sultana, M. S., Rana, S., Yamazaki, S., Aono, T. & Yoshida, S. Health risk assessment for carcinogenic and non-carcinogenic heavy metal exposures from vegetables and fruits of Bangladesh. Cogent Environ. Sci. 3, 1–17 (2017).
USEPA. Risk Assessment Guidance for Superfund. Human Health Evaluation Manual (Part A). vol. I https://rais.ornl.gov/documents/HHEMA.pdf (1989).
WHO. Agents Classified by the International Agency for Research on Cancer (IARC) Monographs, Volumes 1–23. https://monographs.iarc.who.int/wp-content/uploads/2018/09/ClassificationsAlphaOrder.pdf (2016).
Mesa-Frias, M., Chalabi, Z., Vanni, T. & Foss, A. M. Uncertainty in environmental health impact assessment: Quantitative methods and perspectives. Int. J. Environ. Health Res. 23, 16–30 (2013).
Mari, M., Nadal, M., Schuhmacher, M. & Domingo, J. L. Exposure to heavy metals and PCDD/Fs by the population living in the vicinity of a hazardous waste landfill in Catalonia, Spain: Health risk assessment. Environ. Int. 35, 1034–1039 (2009).
Article CAS PubMed Google Scholar
Li, X. et al. Status of copper accumulation in agricultural soils across China (1985–2016). Chemosphere 244, 125516 (2020).
Article CAS PubMed Google Scholar
Wang, Z., Yin, Z., Caers, J. & Zuo, R. A Monte Carlo-based framework for risk-return analysis in mineral prospectivity mapping. Geosci. Front. 11, 2297–2308 (2020).
Sanaei, F. et al. Health risk assessment of potentially toxic elements intake via food crops consumption: Monte Carlo simulation-based probabilistic and heavy metal pollution index. Environ. Sci. Pollut. Res. 28, 1479–1490 (2021).
Guo, G., Zhang, D. & Wang, Y. Probabilistic human health risk assessment of heavy metal intake via vegetable consumption around Pb/Zn Smelters in Southwest China. Int. J. Environ. Res. Public Health 16, 3267 (2019).
Article CAS PubMed PubMed Central Google Scholar
Jiang, C. et al. Distribution, source and health risk assessment based on the Monte Carlo method of heavy metals in shallow groundwater in an area affected by mining activities, China. Ecotoxicol. Environ. Saf. 224, 112679 (2021).
Article CAS PubMed Google Scholar
Egozcue, J. J. & Pawlowsky-Glahn, V. Compositional data: The sample space and its structure. TEST 28, 599–638 (2019).
Thio-Henestrosa, S. & Comas, S. CoDaPack v2 User’s Guide (University of Girona, 2016).
van den Boogaart, K. G. & Tolosana-Delgado, R. Analyzing compositional data with R. Analyz. Compos. Data R https://doi.org/10.1007/978-3-642-36809-7 (2013).
Sall, J., Stephens, M., Loring, S. & Lehman, A. JMP start statistics: A guide to statistics and data analysis using JMP (SAS Institute Inc, 2017).
Zornoza, R., Faz, A., Carmona, D. M., Martínez-Martínez, S. & Acosta, J. A. Plant cover and soil biochemical properties in a mine tailing pond five years after application of marble wastes and organic amendments. Pedosphere 22, 22–32 (2012).
Mohanavelu, A., Naganna, S. R. & Al-Ansari, N. Irrigation induced salinity and sodicity hazards on soil and groundwater: An overview of its causes, impacts and mitigation strategies. Agriculture 11, 983 (2021).
Zewd, I. & Siban, M. The effects of alkalinity on physical and chemical properties of soil. J. Plant Biol. Agric. Sci. 3, 1–5 (2021).
Saha, L. & Bauddh, K. Sustainable agricultural approaches for enhanced crop productivity, better soil health, and improved ecosystem services BT—ecological and practical applications for sustainable agriculture. In Bauddh, K., Kumar, S., Singh, R. P. & Korstad, J. (eds) 1–23 (Springer Singapore, 2020). https://doi.org/10.1007/978-981-15-3372-3_1.
Mishra, A. K. et al. Promising management strategies to improve crop sustainability and to amend soil salinity. Front. Environ. Sci. (2023).
Shakir, E., Zahraw, Z. & Al-Obaidy, A. H. M. J. Environmental and health risks associated with reuse of wastewater for irrigation. Egypt. J. Pet. 26, 95–102 (2017).
Ozcelik, M. Environmental pollution and its effect on water sources from marble quarries in western Turkey. Environ. Earth Sci. 75, 796 (2016).
Noreen, U. et al. Water pollution and occupational health hazards caused by the marble industries in district Mardan, Pakistan. Environ. Technol. Innov. 16, 100470 (2019).
Singare, P., Meshram, N. & Kausadikar, H. Calcareus soil and their management: A review. Just Agric. 3, 1–6 (2022).
Carrillo, V., Fuentes, B., Gómez, G. & Vidal, G. Characterization and recovery of phosphorus from wastewater by combined technologies. Rev. Environ. Sci. Bio/Technol. 19, 389–418 (2020).
Feng, W. et al. Chemical composition, sources, and ecological effect of organic phosphorus in water ecosystems: A review. Carbon Res. 2, 12 (2023).
Krishna, A. K. & Govil, P. K. Heavy metal distribution and contamination in soils of Thane-Belapur industrial development area, Mumbai Western India. Environ. Geol. 47, 1054–1061 (2005).
Bhuiyan, M. A. H., Parvez, L., Islam, M. A., Dampare, S. B. & Suzuki, S. Heavy metal pollution of coal mine-affected agricultural soils in the northern part of Bangladesh. J. Hazard. Mater. 173, 384–392 (2010).
Article CAS PubMed Google Scholar
Charkiewicz, A. E., Omeljaniuk, W. J., Nowak, K., Garley, M. & Nikliński, J. Cadmium toxicity and health effects—a brief summary. Molecules 28, 6620 (2023).
Article CAS PubMed PubMed Central Google Scholar
Genchi, G., Sinicropi, M. S., Lauria, G., Carocci, A. & Catalano, A. The effects of cadmium toxicity. Int. J. Environ. Res. Public Health 17, 3782 (2020).
Article CAS PubMed PubMed Central Google Scholar
Timber Floor Shi, D. & Lu, X. Accumulation degree and source apportionment of tr