key: cord-0067916-s0nypmp3 authors: Golui, Debasis; Datta, S. P.; Dwivedi, B. S.; Meena, M. C.; Ray, P.; Trivedi, V. K. title: A new approach to establish safe levels of available metals in soil with respect to potential health hazard of human date: 2021-09-25 journal: Environ Earth Sci DOI: 10.1007/s12665-021-09988-7 sha: 55140ec4087f3f73fb066f0a9046e83e4b06a593 doc_id: 67916 cord_uid: s0nypmp3 Safe levels of extractable pollutant elements in soil have not been universally established. Prediction of metal solubility in polluted soils and the subsequent transfer of these metals from soil pore water to the human food supply via crops are required for effective risk assessment from polluted soils. Thus an attempt has been made to develop a novel approach to protect human health from exposure to toxic metals through assessing risk from metal polluted soils utilised for agriculture. In this study, we assess the relative efficacy of various forms of ‘free ion activity model’ (FIAM) for predicting the concentration of cadmium (Cd), lead (Pb), nickel (Ni), zinc (Zn) and copper (Cu) in spinach and wheat as example crops, thereby providing an assessment of risk to human health from consumption of these crops. Free metal ion activity in soil solution was estimated using the Windermere Humic Aqueous Model VII (WHAM-VII) and the Baker soil test. Approximately 91, 81, 75, 94 and 70% of the variability in Cd, Pb, Ni, Zn and Cu content, respectively, of spinach could be described by a FIAM using an estimate of the free ion activity of the metals provided by WHAM-VII. Owing to the different concentration of ethylenediamine tetraacetic acid (EDTA) and diethylenetriamine pentaacetic acid (DTPA) used in the present experiment, higher prediction coefficients were obtained using EDTA (0.05 M), rather than DTPA (0.005 M), as the metal extractant in an integrated solubility-FIAM model. Out of three formulations, the FIAM, based on free ion activity of metals in soil pore water, determined from solution extracted with Rhizon samplers, was distinctly superior to the other formulations in predicting metal uptake by spinach and wheat. A safe level of extractable metal in soil was prescribed using a hazard quotient derived from predicted plant metal content and estimated dietary intake of wheat and spinach by a human population. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s12665-021-09988-7. Use of waste waters (industrial effluents and domestic sewage) for irrigating crop lands is expanding steadily owing to the paucity of fresh water reserves (Golui et al. 2020) . The presence of substantial amounts of essential plant nutrients in waste water may enhance the agricultural productivity (Chen et al. 2005) . Application of sewage effluents has been reported to improve soil chemical and physical characteristics including organic carbon as well as major and micronutrients (Meena et al. 2016) . Besides, enrichment of the edible portions of rice and wheat has also been reported due to irrigation with sewage (Meena et al. 2016) . Conversely, the application of such effluents on farming land frequently results in build-up of toxic metals in soils in the long term (Rattan et al. 2005; Deshmukh et al. 2015; Meena et al. 2016 ). It has been established beyond doubt that crops grown on such metal-loaded soils may surpass the legislative limits for toxic metals in the palatable portions of crops (Nabulo et al. 2011) . Excessive intake of metals through the human diet may lead to several health complications. For example, excess cadmium (Cd) may result in aminoaciduria, namely itai-itai disease. Encephalopathy, failure in reproduction and metabolic disorders have been associated with lead toxicity in humans. Zinc toxicity may lead to a reduction in the functioning of Fe in human, thereby causing anaemia (Rattan et al. 2009 ). Nickel and Cu are tumourdeveloping agents, whose mutagenic effect has attracted worldwide concerns. Exposure to nickel-enriched dust may induce nasopharyngeal carcinoma in humans (Chen 2011) . Therefore, an assessment of the suitability of agricultural land for crop production in respect of its metal status should be an integral part of preventing food-chain contamination. Over the years, total soil metal content has been used as the simplest index of metal hazard, but this does not take into consideration the effect of soil properties on metal solubility and bioavailability. Hence, such indices are not always meaningful in either protecting human health from metal hazards or judging the suitability of agricultural lands for crop production. Use of too stringent permissible limits is not desirable, considering the current continuing loss of agricultural land. Furthermore, such permissible limits should also be crop-specific. The free ion activity in the soil pore water, described as an 'intensity factor', is a key indicator of metal-related risk in terms of leaching and bioavailability (Oorts et al. 2006) . Most studies suggest that measurement of 'intensity' of metals in soil is preferred to the determination of metal 'quantity factors' for prediction of bioavailability of metals and associated risk (Datta and Young 2005; Meena et al. 2016; Golui et al. 2017; Mandal et al. 2019) . Thus the intensity of metal in pore water has a direct effect on the metal bioavailability to crops and ecotoxicity to microbial communities (Vulkan et al. 2000; Groenenberg et al. 2010) . One way of measuring the intensity of metal in soils is in-situ withdrawal of soil pore water through Rhizon sampler from the root zone of growing crops, followed by speciation using geochemical models such as WHAM VII (Tipping et al. 2011; Marzouk et al. 2013; Mao et al. 2017; Mishra et al. 2019) . Free ion activity may be estimated by the simple Baker soil test programme of Baker and Amacher (1981) . Empirical and semiempirical procedures have also been applied to replace more tedious soil solution extraction and speciation procedures (Hough et al. 2003; Tye et al. 2003; Datta and Young 2005) . Recently, the efficacy of the simpler approaches (Baker soil test, empirical and semi empirical) for determining free metal ion activity in soil solution was evaluated against the more direct technique utilizing Rhizon sampler for extraction of soil solution and following speciation by WHAM VII model (Golui et al. 2020) . However, the usefulness of these approaches have not been evaluated as predictors of plant metal uptake. Precise prediction of plant uptake of metal is of supreme consequence for prescribing the permissible limit of metals in soil in respect of the knowledge of the dietary uptake of metals via ingestion of foodstuff cultivated on metal polluted soils. In this regard, there is limited information on the application of free ion activity models (FIAM) for predicting transport of metal from soil to crop, particularly in tropical soils with low organic carbon content. The FIAM for prediction of transfer of metal to plant was conceptualized in the early 1980's and has been adapted following the principles of the biotic ligand model (Campbell 1995; Datta and Young 2005) . The FIAM is conceptually built on the interchange and binding of free metal ions with the cellular binding sites on plant roots built on a common chemical equilibrium theory. Several formulations of the FIAM have been implemented using both measured free ion activity (M 2+ ) and modelled (M 2+ ) data (Hough et al. 2004; Datta and Young 2005) . In an extended procedure, uptake of metal by plants can be predicted reasonably well using an integrated solubility and FIAM approach, based on the predicted ion activity of metal in the soil solution (Hough et al. 2003 (Hough et al. , 2004 Golui et al. 2017) . Plant metal content as predicted by the most accurate free ion activity model could be used for assessing risk from dietary metal intake by humans and thereby prescribing permissible limits for metal concentrations in soil. The present study attempted to develop a comprehensive tool for risk assessment of polluted soil. Objectives of the present study are 1) assessing the relative efficacy of several formulations of the FIAM for predicting the transport of zinc, copper, nickel, lead and cadmium using spinach and wheat as test crops, (2) to develop risk assessment protocols for polluted soil pertaining to health hazard of human, and (3) to prescribe safe limits of extractable metals in soil. Contaminated soil samples (25) in bulk were brought from various sites ( Fig. 1 and Table S1 ). Four soil samples in bulk were also brought from nearby agricultural sites which did not have the history of receiving industrial effluents or solid waste. After processing (air-drying followed by sieving with 2-mm sieve), samples were analysed for important chemical properties. Samples were analysed in triplicate for pH, organic carbon and texture following the protocol of Jackson (1973) , Walkley and Black (1934) and Bouyoucos (1962) , respectively. Total concentrations of zinc, copper, nickel, lead and cadmium were measured in inductively coupled plasma-mass spectrometry (ICP-MS) (Perkin Elmer NexIon 300) following digestion with mixture of HNO 3 and HCl @ 1:3 (Quevauviller, 1998) . Concentrations of metals extractable with EDTA (0.05 M) and DTPA (0.005 M) solution were also measured by ICP-MS following the protocol of Quevauviller (1998) and Lindsay and Norvell (1978) , respectively. An experiment was carried out in greenhouse, using processed soil from each location, with spinach and wheat (details in supplementary information, section A1). Soil pore water was collected from the pot experiment using Rhizon samplers. Prior to sampling of soil solution, pots were irrigated. After draining of excess moisture, soil solution was collected in situ according to Golui et al. (2020) . A portion of the filtered soil solution was used for analysis of Na, Mg, K, Ca, Zn, Cu, Fe, Mn, Ni, Pb, Cd; the remainder was used for determination of Cl, NO 3 , SO 4 , PO 4 , dissolved organic-and inorganic carbon as reported by Golui et al. (2020) . Speciation analysis was carried out using WHAM-VII for determination of the free metal ion activity following Tipping et al. (2011) . The free metal ion activities in soils were also estimated following the standard procedure described by Baker and Amacher (1981) (details in supplementary information, section A2). The following forms of the FIAM were used to predict the metal uptake by spinach and wheat. An 'integrated solubility-FIAM' was used for prediction of metal uptake by spinach and wheat grain. The FIAM suggests that metal uptake by plant is governed by metal ion activity in soil solution. A transfer factor (TF) is given as the ratio of [M plant ] to (M 2+ ) (Eq. 1). where [M plant ] is the plant metal content; (M 2+ ) is free metal ion activity. The free ion activity of metal (M 2+ ) was predicted from a pH-dependent Freundlich equation according to Golui et al. (2020) . Metal uptake by spinach and wheat can be predicted by combining Eq. (1) with (M 2+ ) as follows: where C, β 1 and β 2 are empirical metal and plant-specific coefficients. The 'Solver' facility in Microsoft Excel 2020 was used to parameterize Eq. (2) through non-linear error minimization. For calculation of error sum of squares, numerical data on plant metal content were used rather than logarithmic data. Prediction of plant metal content was made using a biotic ligand model (BLM) formulation. The underlying principle of this approach is that metal sorption (1) occurs on assumed root surface and (2) is in competition with positive ions (H + and cations) for the root surface. Total root surface can be given as: R T = total root surface; RM = Metal absorbed sites; RH = H + absorbed sites; R = Free sites. The procedure of absorption of metal and H + can be given as Eqs. 4 and 5: (3) R T = RM + RH + R K M and K H do not consider possible effects of variation in pH or positive ion absorption on root membrane properties. Adding Eqs. 3-5 and assuming that metal uptake by spinach /wheat grain, namely [M spinach/wheat ], is proportional to RM, leads to Eq. (6). where K t is a proportionality constant presuming that the content of metal ions on root sorption sites during the entire growth period of crop is reflected by the concentration of metal in the shoot and grain of the plant. Parameterization of model equation (Eq. 6) was done in 'Solver' in MS Excel 2020. For Model III, the formulation of the FIAM was identical to Model II except that the free ion activity of metal (M 2+ ) in soil pore water was estimated using the Baker soil test programme. Hazard quotients (HQ) were computed following US-EPA protocols (IRIS 2020) to appraise health hazard of human due to dietary ingestion of metals through spinach and wheat grain grown on polluted soils. where ADD is average daily dose (mg kg −1 day −1 ) and R f D is reference dose (mg kg −1 day −1 ) of metals. The R f D values are 0.3 for Zn, 0.5 for Cu, 0.02 for Ni, 0.0035 for Pb and 0.001 for Cd (IRIS 2020; WHO 1982) . Recommended dietary consumption of spinach of 0.2 kg day −1 (fresh weight) was considered for computation of HQ, while in case of wheat it is 0.4 kg day −1 . The HQ for metal intake with the consumption of spinach and wheat was calculated according to Eq. (8) for an adult (average body weight is 70 kg). where M plant is the concentration of metal in spinach and wheat; W is the consumption of vegetables and wheat (kg); F is the conversion factor i.e. 0.082 for spinach and 1 for wheat. Values of hazard quotient up to unity are considered as safe, when all probable routes of entry of metal to human body are calculated i.e. a complete risk assessment. But in the present case, assessment of risk is not complete, since metal intake by humans can occur by other possible routes as well. Hence, Meena et al. (2016) proposed a safe limit of HQ for rice and wheat as 0.5. Owing to the fact that dietary intake of green vegetable by human being is far less than staple food like rice or wheat, a safe limit of HQ for green vegetable is proposed here as 0.25. Hence, critical values of HQ in case of spinach and wheat were considered as 0.25 and 0.50, respectively for ascertaining the permissible level of available metal at a certain organic carbon and soil pH under modelling framework. The reason behind fixation of relatively lower critical value of HQ for spinach is based on the presumption that a small amount of the daily human diet is constituted by green vegetables. The mean value of soil pH was 7.61 ± 0.31; the extent of soil pH in the twenty nine experimental soils varied from neutral to alkaline (Table S2 ). The mean soil organic carbon content was 1.12 ± 0.24%; the Dhapa soil had the highest value (2.68 ± 0.40%), where solid wastes of municipalities have been deposited for the last forty years. High organic matter contents may be associated with long-term application of municipal solid wastes, sewage-sludge, and industrial effluents. The mean clay content was 22.5 ± 6.71%. The soils covered six textural classes: silty clay loam, clay, clay loam, sandy loam, loam, and sandy clay loam. The range of total metal concentration across the studied soils were 28. 2-28,662, 23.6-2305, 11.4-1513, 7 .90-3793 and 0.22-352 mg kg −1 for Zn, Cu, Ni, Pb and Cd, respectively with the corresponding median value of 381, 226, 25.6, 22.2 and 2.13 mg kg −1 (Table S3 ). An apparent elevation of Cd, Pb and Zn in Debari soils may have resulted from the application of Zn-smelter effluents from smelter plants on a long-term basis. The Dhapa site receiving municipal solid was located adjacent to Kolkata city roads with high traffic volumes and exhibited the highest level of Pb, thereby (probably) reflecting emissions from automobiles (USEPA 2017). The total Ni content was higher in the soils of Sonepat, reflecting irrigation using cycle industry effluents where Ni is included in the production of stainless steel (Cempel and Nikel 2006 (Table S4 ). The concentrations of extractable metals were used as indicator of the labile pool of metal for prediction of metal uptake by plants using Model I. The Baker extracts showed significantly higher concentrations of all metals (5.57, 2.38, 0.33, 0.73 and 0.41 mg L −1 for Zn, Cu, Ni, Pb and Cd, respectively) compared to the Rhizon extracts (Table S5 ). Soil solution extracted through in-situ Rhizon samplers directly represents the pore water of the root zone. Whereas, the Baker extracts use DTPA for soil extraction to estimate free metal ion activity in soil pore water. The mean free metal ion activities in soil pore water extracted from the rhizosphere of spinach using Rhizon samplers, and estimated using WHAM VII, were 6.93 ± 0.32 for pZn 2+ , 10.1 ± 0.68 for pCu 2+ , 7.70 ± 0.28 for pNi 2+ , 10.3 ± 0.32 for pPb 2+ and 9.08 ± 0.45 for pCd 2+ (Table S6) . Corresponding values under wheat rhizospheres were 7.03 ± 0.38, 10.2 ± 0.59, 7.65 ± 0.31, 10.7 ± 0.42 and 9.28 ± 0.54, respectively (Table S7 ). The Baker soil test also provided similar location-specific changes in free metal ion activities. The average free metal ion activities as estimated with the Baker soil test programme were 10.1 ± 1.12 for pZn 2+ , 13.4 ± 1.23 for pCu 2+ , 12.9 ± 0.85 for pNi 2+ , 11.6 ± 0.74 for pPb 2+ and 12.6 ± 2.26 for pCd 2+ (Table S8) . Total metal concentration in spinach and wheat varied considerably among the experimental soils. In edible portion of spinach, metal concentration ranged from 57.2 to 1245 mg kg −1 for Zn, 3.85 to 40.2 mg kg −1 for Cu, 0.55 to 15.5 mg kg −1 for Ni, 0.40 to 82.0 mg kg −1 for Pb and 0.15 to 48.3 mg kg −1 for Cd (Table S9) . The corresponding median values were 164, 8.31, 1.04, 1.93 and 0.48 mg kg −1 for Zn, Cu, Ni, Pb and Cd, respectively. By contrast, total metal concentrations in wheat grain varied from 12.6 to 98.5 mg kg −1 for Zn, 1.53 to 10.3 mg kg −1 for Cu, 0.18 to 2.80 mg kg −1 for Ni, 0.08 to 0.74 mg kg −1 for Pb and 0.02 to 10.4 mg kg −1 for Cd. The median values of metal concentrations in wheat grain were 54.5, 5.26, 0.38, 0.21 and 0.17 mg kg −1 for Zn, Cu, Ni, Pb and Cd, respectively. Metal concentrations in dicotyledonous spinach (high root CEC) were considerably higher than in monocotyledonous wheat grain (low root CEC). Higher rates of transpiration in spinach lead to greater absorption of metals by plants (Zhou et al. 2016 ). These two contrasting crops in respect of metal uptake efficiency and dietary intake by human were used as test crops to ensure wider applicability of permissible limits. The prediction coefficients (R 2 ) and model parameters (C, β 1 and β 2 ) of the integrated solubility and free ion activity model (Model I) are presented in Table 1 . It has been found that 78, 57, 64, 67 and 93% of the change in Zn, Cu, Ni, Pb and Cd concentrations of spinach, respectively, could be described by soil reaction and M C (labile pool of metals assumed to be adsorbed on organic matter as extracted by EDTA) ( Figure S1 ). By contrast, the prediction coefficients of the model as obtained for Zn, Cu, Ni, Pb and Cd were 0.51, 0.55, 0.63, 0.66 and 0.87, respectively using labile pool of metals as extracted with DTPA ( Figure S2 ). The Model I (EDTA) based on EDTA extractable metal could describe the changes in metal concentration in wheat grain to the extent of 53% for Zn, 42% for Cu, 85% for Ni, 52% for Pb and 88% for Cd (Table 1) . Prediction coefficients of 0.71 for Zn, 0.19 for Cu, 0.81 for Ni, 0.45 for Pb and 0.87 for Cd were obtained for wheat grain in Model I (DTPA). This model involves determination of soil characteristics like soil reaction, soil organic carbon and available metals. Meena et al. (2016) reported that organic carbon and pH are among the key soil chemical characteristics which govern the metal solubility in polluted soils. Considerable variation in the crop and the metal-specific constants indicated the uniqueness of model parameters for each metal and crop. Values of β 2 in Model I were positive for all metals in both the crops. Thus, as expected, elevated concentrations of available metals in experimental soil will magnify the concentration of studied metals in the edible portion of spinach and wheat. In Model I, the free ion activity of metal in soil solution was predicted through Freundlich equation (solubility model). Hence, it is likely that the use of a more directly estimated free metal ion activity would enhance the predictive power of this formulation of the FIAM, as presented in the following section of the paper. However, values of prediction coefficients (R 2 ) were more than 0.5 in all cases. Overall, for both crops, the use of EDTA-extractable metal concentration in the model yielded higher prediction coefficients compared to the DTPA-extractable metals. This may be due to the higher efficiency of EDTA in extracting metals from highly polluted soils due to the higher concentration of the extractant (ten times) in comparison to DTPA (0.005 M). The DTPA extractant may be capacity limited and may not reflect the full labile metal pool in soils with a large available metal concentration. Findings of the present research work is in concurrence with the findings of (Zan et al. 2013) , where it was reported that prediction of free ion activity by solubility model (pH dependent Freundlich equation) based on EDTA-extractable metal yielded higher prediction coefficient in comparison to model based on labile pool of metal as extracted with DTPA. The present study attempted to predict the metal uptake by plant using FIAM based on estimated free ion activity of metal in soil solution (WHAM VII and Baker soil test extract). The R 2 values of model II (Rhizon) and model III (Baker) for different metals are presented in expected to be closer to that of the pore water as compared to that of Baker soil test. Further, WHAM VII is a robust ion speciation model based on a large experimental data set, whereas the Baker test algorithm is based on more limited data and requires that soil is extracted with a chelating agent and salt solution. On the other hand, in theory, Rhizon samplers extract the soil solution which is actually bathing the surface of the plant roots and therefore provides a more realistic measure of metal ion intensity. The relative efficacy of different formulations of FIAM in predicting uptake of metal by spinach and wheat was compared in terms of the mean prediction coefficient ( Table S10 ). The highest mean prediction coefficient (R 2 = 0.82) for uptake of metals by spinach was recorded for Model II (Rhizon) based on the integrated use of Rhizon-WHAM VII. On the other hand, Model III (Baker), based on Baker soil test data was less effective in predicting uptake of metals by spinach (mean R 2 = 0.62). Model I, based on EDTA-extractable metal yielded a higher mean prediction coefficient (mean R 2 = 0.72) in comparison to the use of DTPA-extractable metal (mean R 2 = 0.64). The efficacy of different formulations of FIAM for predicting metal uptake by plants follows the order of Model II (Rhizon) > Model I (EDTA) > Model I (DTPA) > Model III (Baker) . The value of HQ below 1 for complete risk assessment, i.e. where HQ is based on total metal intake from various exposure paths (food materials, drinking water ingestion of soil and inhalation of dust) has been considered safe as per USEPA. However, in the present study a partial Table 2 Model parameters for predicting uptake of Zn, Cu, Ni, Pb and Cd by spinach and wheat as a function of estimated free metal ion activity (M 2+ ) and soil pH The values of R 2 > 0.18 are significant at 5% probability level a As estimated by WHAM VII-Model II risk assessment was carried out, i.e. only metal intake by consumption of palatable portion of spinach and wheat was considered. Being a staple food, wheat constitutes the major portion of the diet in India, whereas the proportion of spinach in the diet is small. Considering this fact, the upper safe limits of HQ for staple food (wheat) and green leafy vegetables (spinach) have been fixed at 0.50 and 0.25, respectively. The values of HQ as calculated for spinach ranged from 0.006 to 0.790 for Ni, 0.027 to 5.492 for Pb and 0.035-11.32 for Cd (Table 3) . The values of HQ in respect of Ni were less than 0.25 except for Sonepat soils (Soil No. 26) . Debari and Sonepat soils showed high values of HQ for Pb and Cd, thereby indicating that leafy vegetables raised on these soils may not be suitable for human intake. The HQ values of Pb exceeded the safe limit of 0.25 for spinach grown on Dhapa soil. The HQ values in connection with the consumption of wheat grain are presented in Table 4 . Average values of HQ for Ni, Pb and Cd were 0.22, 0.44 and 5.97, respectively. For Ni, the HQ values were lower than 0.50 in all the studied soils except Sonepat soil (all soil samples from Sonepat) and Keshopur soil (Soil No. 13) . Higher HQ for Cd in twenty four soil samples indicates that wheat crop raised on these soils is not suitable for human intake. There is strongest evidence that Zn is required for immune-support of human being (Gombart et al. 2020) . Zinc is also required to check the incidence of respiratory tract infection, pneumonia and diarrhea, thereby increasing the immunity in Figure 2 shows that values of (M 2+ ) and M spinach were either, respectively, (i) both calculated, (ii) calculated and predicted through solubility and FIAM (iii) calculated and predicted through the FIAM. There was reasonable agreement among all the four plots ( Fig. 1) , which may explain the utility of soil chemical properties, namely soil pH, organic carbon and extractable metal to measure risk from edibles consumption. Permissible limits for soil metal concentrations were also ascertained on the basis of predicted HQ using the plant metal content modelled by model II i.e. FIAM based on measured free ion activity involving Rhizon-WHAM, for intake of Cd in humans due to consumption of spinach and wheat grain (Figs. 3, 4) . The level of extractable Cd in soil corresponding to the HQ > 0.25 for Cd intake via consumption of spinach grown thereon was considered as the permissible limit. The DTPA-extractable Cd in soil corresponding to the HQ > 0.50 was considered as safe limit for wheat crop. Figure 2 shows that the safe limit of DTPA-extractable Cd in soil will be 0.05 mg kg −1 for spinach at pH of 6 and organic carbon content of 0.25%; the corresponding safe limit of Cd will be 0.30 mg kg −1 at pH of 8.0 and organic carbon of 0.5%. Similarly, the safe limit of DTPA-extractable Cd in soil will be 0.03 mg kg −1 for wheat at pH 6.0 with 0.25% organic carbon; the corresponding safe limit of Cd will be 0.11 mg kg −1 , if pH and organic carbon remains at 8.0 and 0.5%, respectively (Fig. 3) . It is evident that HQ with respect to Cd for spinach and wheat grown in the experimental soils exceeded in several cases. A ready reckoner was generated for calculating crop specific safe level of extractable Cd in soil in connection with human health hazard. In the experimental soil, the safe level of extractable Cd in soil ranged broadly with the variation in soil organic carbon, while such variation was not noted with soil pH. This can be ascribed to the narrow scale of pH (alkaline pH scale) of the experimental soils. Such research findings have practical implication for fixing the safe limit of extractable metal considering Fig. 2 Comparison of observed and predicted hazard quotients for Cd in spinach; prediction was made from the free ion activity model using values of (Cd 2+ ) ions either from speciation of the extracted pore water or using the generic solubility model Fig. 3 Permissible limit of DTPA extractable Cd in soils in relation to solubility of metal for intake of Cd through spinach by human. *Values in parentheses indicate the toxic limit of extractable Cd in soil Fig. 4 Permissible limit of DTPA extractable Cd in soils in relation to solubility of metal for intake of Cd through wheat grain by human. *Values in parentheses indicate the toxic limit of extractable Cd in soil the important soil properties, since total metal content is not a good metal hazard index to human health vis-à-vis phytoavailability. The novelty of this study is a comparative evaluation of efficacy of different formulations of FIAM in predicting the metal uptake by crops. Furthermore, prescription of safe level of metal in soils in connection with health hazard of human for intake of metal by food materials is a new concept, and should constitute priority area of research in the field of metal pollution. There is also scarcity of free metal ion activity data as speciated through most recent and robust speciation model, WHAM VII. The FIAM, based on measured free ion activity of metals in soil solution extracted by Rhizon sampler, was distinctly superior to other formulations in predicting the metal uptake by spinach and wheat. The efficacy of Model I was more or less at par in predicting the metal uptake by plants. The protocol used in the present investigation for assessing the safe level of extractable metal may be practical for the appraisal of risk related to contaminated soils on routine basis, as well as in devising the management options for specific contaminated sites. We expect that this novel approach of risk assessment will prove to be useful and promising. Availability of data and material Not applicable. Code availability Not applicable. The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. Development and interpretation of a diagnostic soil testing program. Pennsylvania Agricultural Experiment Station Bulletin No. 826. Pennsylvania Agricultural Experiment Station Interaction between trace metals and aquatic organisms: a critique of the free-ion activity model Nickel, A review of its sources and environmental toxicology Review of the research on heavy metal contamination of China's city soil and its treatment method Residues and source identification of persistent organic pollutants in farmland soils irrigated by effluents from biological treatment plants Predicting metal uptake and risk to human food chain from leafy vegetables grown on soils amended by long-term application of sewage sludge Safe limit of arsenic in soil in relation to dietary exposure of arsenicosis patients from Malda district, West Bengal-a case study Prediction of free metal ion activity in contaminated soils using WHAM VII, baker soil test and solubility model A review of micronutrients and the immune system-working in harmony to reduce the risk of infection Transfer functions for solid-solution partitioning of cadmium, copper, nickel, lead and zinc in soils: derivation of relationship for free metal ion activities and validation with independent data Assessing potential risk of heavy metal exposure from consumption of home produced vegetable by urban populations Dietary reference intakes for Vitamin A, Vitamin K, Arsenic Integrated risk information system database. US Environmental Protection Agency Soil chemical analysis. Prentice Hall of India Private Limited Development of a DTPA soil test for zinc, iron, manganese and copper Assessing equilibria of organoarsenic complexes and predicting uptake of arsenic by wheat grain from organic matter amended soils Predicting trace metal solubility and fractionation in Urban soils from isotopic exchangeability Predicting the solubility and lability of Zn, Cd, and Pb in soils from a minespoil-contaminated catchment by stable isotopic exchange Longterm impact of sewageirrigation on soil properties and assessing risk in relation to transfer of metals to humanfood chain. Environ Enhancing the effectiveness of zinc, cadmium and lead phytoextraction in polluted soils by using amendments and microorganisms Trace metal uptake by tropical vegetables grown on soil amended with sewage sludge Soil properties affecting the toxicity of CuCl 2 and NiCl 2 for soil microbial processes in freshly spiked soils Operationally defined extraction procedures for soil and sediment analysis Long-term impact of irrigation with sewage effluents on heavy metal content in soils, crops and groundwater, a case study Micronutrients in soil, plant, animal and human health Humic Ion-Binding Model VII: a revised parameterisation of cation-binding by humic substances Predicting the activity of Cd 2+ and Zn 2+ in soil pore water from the radio-labile metal fraction USEPA (2017) Sources of lead in soil: a literature review Copper speciation and impacts on bacterial biosensors An examination of the degtjareff method for determining soil organic matter, and a proposed modification of the chromic acid titration method Joint FAO/WHO expert committee on food additives. WHO Food Additive Series No. 683. World Health Organization Accumulation of heavy metals in vegetable species planted in contaminated soils and the health risk assessment Animal research Not applicable. Consent to publish Not applicable.