Estimated Soil Ingestion Rates for Use in Risk Assessment

Assessing the risks to human health posed by contaminants present in soil requires an estimate of likely soil ingestion rates. In the past, direct measurements of soil ingestion were not available and risk assessors were forced to estimate soil ingestion rates based on observations of mouthing behavior and measurements of soil on hands. Recently, empirical data on soil ingestion rates have become available from two sources (Binder et al., 1986 and van Wijnen et al., 1986). Although preliminary, these data can be used to derive better estimates of soil ingestion rates for use in risk assessments. Estimates of average soil ingestion rates derived in this paper range from 25 to 100 mg/day, depending on the age of the individual at risk. Maximum soil ingestion rates that are unlikely to underestimate exposure range from 100 to 500 mg. A value of 5,000 mg/day is considered a reasonable estimate of a maximum single-day exposure for a child with habitual pica.     

Equation reliability of soil ingestion estimates in mass-balance soil ingestion studies

Exposure to chemicals from ingestion of contaminated soil may be an important pathway with potential health consequences for children. A key parameter used in assessing this exposure is the quantity of soil ingested, with estimates based on four short longitudinal mass-balance soil ingestion studies among children. The estimates use trace elements in the soil with low bioavailability that are minimally present in food. Soil ingestion corresponds to the excess trace element amounts excreted, after subtracting trace element amounts ingested from food and medications, expressed as an equivalent quantity of soil. The short duration of mass-balance studies, different concentrations of trace elements in food and soil, and potential for trace elements to be ingested from other nonsoil, nonfood sources contribute to variability and bias in the estimates. We develop a stochastic model for a soil ingestion estimator based on a trace element that accounts for critical features of the mass-balance equation. Using results from four mass-balance soil ingestion studies, we estimate the accuracy of soil ingestion estimators for different trace elements, and identify subjects where the difference between Al and Si estimates is larger (>3 RMSE) than expected. Such large differences occur in fewer than 12% of subjects in each of the four studies. We recommend the use of such criteria to flag and exclude subjects from soil ingestion analyses.     

Modeled Estimates of Soil and Dust Ingestion Rates for Children

Daily soil/dust ingestion rates typically used in exposure and risk assessments are based on tracer element studies, which have a number of limitations and do not separate contributions from soil and dust. This article presents an alternate approach of modeling soil and dust ingestion via hand and object mouthing of children, using EPA's SHEDS model. Results for children 3 to <6 years old show that mean and 95th percentile total ingestion of soil and dust values are 68 and 224 mg/day, respectively; mean from soil ingestion, hand‐to‐mouth dust ingestion, and object‐to‐mouth dust ingestion are 41 mg/day, 20 mg/day, and 7 mg/day, respectively. In general, hand‐to‐mouth soil ingestion was the most important pathway, followed by hand‐to‐mouth dust ingestion, then object‐to‐mouth dust ingestion. The variability results are most sensitive to inputs on surface loadings, soil‐skin adherence, hand mouthing frequency, and hand washing frequency. The predicted total soil and dust ingestion fits a lognormal distribution with geometric mean = 35.7 and geometric standard deviation = 3.3. There are two uncertainty distributions, one below the 20th percentile and the other above. Modeled uncertainties ranged within a factor of 3–30. Mean modeled estimates for soil and dust ingestion are consistent with past information but lower than the central values recommended in the 2008 EPA Child‐Specific Exposure Factors Handbook. This new modeling approach, which predicts soil and dust ingestion by pathway, source type, population group, geographic location, and other factors, offers a better characterization of exposures relevant to health risk assessments as compared to using a single value.     

A Probabilistic Arsenic Exposure Assessment for Children Who Contact CCA-Treated Playsets and Decks, Part 1: Model Methodology, Variability Results, and Model Evaluation

Concerns have been raised regarding the safety of young children who may contact arsenic residues while playing on and around chromated copper arsenate (CCA)-treated wood playsets and decks. Although CCA registrants voluntarily canceled the production of treated wood for residential use in 2003, the potential for exposure from existing structures and surrounding soil still poses concerns. The EPA's Office of Research and Development developed and applied the probabilistic Stochastic Human Exposure and Dose Simulation model for wood preservatives (SHEDS-Wood) to estimate children's absorbed dose of arsenic from CCA. Skin contact with, and nondietary ingestion of, arsenic in soil and wood residues were considered for the population of children in the United States who frequently contact CCA-treated wood playsets and decks. Model analyses were conducted to assess the range in population estimates and the impact of potential mitigation strategies such as the use of sealants and hand washing after play events. The results show predicted central values for lifetime annual average daily dose values for arsenic ranging from 10(-6) to 10(-5) mg/kg/day, with predicted 95th percentiles on the order of 10(-5) mg/kg/day. There were several orders of magnitude between lower and upper percentiles. Residue ingestion via hand-to-mouth contact was determined to be the most significant exposure route for most scenarios. Results of several alternative scenarios were similar to baseline results, except for the scenario with greatly reduced residue concentrations through hypothetical wood sealant applications; in this scenario, exposures were lower, and the soil ingestion route dominated. SHEDS-Wood estimates are typically consistent with, or within the range of, other CCA exposure models.     

Probabilistic Health Risk Assessment for Exposures to Estuary Sediments and Biota Contaminated with Polychlorinated Biphenyls, Polychlorinated Terphenyls and Other Toxic Substances

A Latin Hypercube probabilistic risk assessment methodology was employed in the assessment of health risks associated with exposures to contaminated sediment and biota in an estuary in the Tidewater region of Virginia. The primary contaminants were polychlorinated biphenyls (PCBs), polychlorinated terphenyls (PCTs), polynuclear aromatic hydrocarbons (PAHs), and metals released into the estuary from a storm sewer system. The exposure pathways associated with the highest contaminant intake and risks were dermal contact with contaminated sediment and ingestion of contaminated aquatic and terrestrial biota from the contaminated area. As expected, all of the output probability distributions of risk were highly skewed, and the ratios of the expected value (mean) to median risk estimates ranged from 1.4 to 14.8 for the various exposed populations. The 99th percentile risk estimates were as much as two orders of magnitude above the mean risk estimates. For the sediment exposure pathways, the stability of the median risk estimates was found to be much greater than the stability of the expected value risk estimates. The interrun variability in the median risk estimate was found to be +/-1.9% at 3000 iterations. The interrun stability of the mean risk estimates was found to be approximately equal to that of the 95th percentile estimates at any number of iterations. The variation in neither contaminant concentrations nor any other single input variable contributed disproportionately to the overall simulation variance. The inclusion or exclusion of spatial correlations among contaminant concentrations in the simulation model did not significantly effect either the magnitude or the variance of the simulation risk estimates for sediment exposures.     

Health Risks of PCB Spills from Electrical Equipment

The Southern California Edison Company (SCE) has instituted a series of control strategies designed to minimize human exposure to polychlorinated biphenyls (PCBs) in electrical equipment used on its system. This paper describes a method of analyzing PCB risks using conservative estimates of human intake of PCBs originating from accidental spills from electrical equipment. The PCB releases from the Edison system were determined. The fate of these releases in soil, air, and water was analyzed to determine how much material reaches human receptors. The air and water pathways were determined to be the most likely candidates for the exposure and risk considerations. PCB intake via ingestion of soil at the spill site was neglected as an exposure pathway. Equipment spills without controls resulted in at the most 2 ng/day human intake of PCBs via the water exposure pathway. This was determined to be negligible in comparison with intake rates used in conjunction with the setting of food tolerance levels based on fish being the main dietary pathway of human exposure. The inhalation exposure of the hundred or so persons in the immediate vicinity of a spill was determined to equal the PCB intakes of the fish-eating subpopulation analyzed by the Food and Drug Administration for 2 ppm tolerance standard in the case of no controls or cleanup. Current cleanup procedures assure that even the persons in the immediate area are well below the intake of the subjects in the fish contamination analysis. All exposures were well below a "virtual safe dose" level estimated in the fish tolerance study.     

Ecological Risks to Fossorial Vertebrates from Volatile Organic Compounds in Soil

The past several years has seen an increased awareness of the need to conduct ecological risk assessments (ERAs) for hazardous waste sites. One technique used in ERAs involves estimating contaminant exposure to individual animals of selected species, which is then compared to a reference dose derived from the literature. Exposure estimates are conducted on those species which are representative of the different trophic levels found at the site. In many terrestrial systems, fossorial (burrowing) vertebrates are found in both lower and upper trophic levels. As part of the ERA conducted for Site 300, Lawrence Livermore National Laboratory's high-explosive test facility, contaminant exposures were estimated for fossorial and nonfossorial vertebrates spanning two trophic levels. The results of the evaluation revealed that a significant pathway by which fossorial vertebrates could be exposed to trichloroethylene in soil was through the inhalation of contaminated subsurface burrow air. This was the first time that the importance of this ecological exposure pathway has been recognized. The results of this analysis suggest that further research into the ecological significance of subsurface burrow air contaminated with volatile organic compounds is warranted.     

Children''s Exposure to Arsenic from CCA-Treated Wooden Decks and Playground Structures

CCA-treated wood is widely used in the fabrication of outdoor decks and playground equipment. Because arsenic can be removed from the surface of CCA-treated wood both by physical contact and by leaching, it is important to determine whether children who play on such structures may ingest arsenic in quantities sufficient to be of public health concern. Based on a review of existing studies, it is estimated that arsenic doses in amounts of tens of micrograms per day may be incurred by children having realistic levels of exposure to CCA-treated decks and playground structures. The most important exposure pathway appears to be oral ingestion of arsenic that is first dislodged from the wood by direct hand contact, then transferred to the mouth by children's hand-to-mouth activity. The next most important pathway appears to be dermal absorption of arsenic, while ingestion of soil that has become contaminated by leaching from CCA-treated structures appears to be of lesser importance, except possibly in the case of children with pica. Considerable uncertainty, however, is associated with quantitative estimates of children's arsenic exposure from CCA-treated wood. Priorities for refining estimates of arsenic dose include detailed studies of the hand-to-mouth transfer of arsenic, studies of the dermal and gastrointestinal absorption of dislodgeable arsenic, and studies in which doses of arsenic to children playing in contact with CCA-treated wood are directly determined by measurement of arsenic in their urine, hair, and nails.     

Urinary Excretion of Chromium Following Ingestion of Chromite-Ore Processing Residues in Humans: Implications for Biomonitoring

Biomonitoring programs for urinary chromium (Cr) typically attempt to evaluate occupational exposure via the inhalation route. This study investigated whether Cr can be detected in the urine of people following the ingestion of soils that contain relatively high concentrations of chromium in chromite ore processing residue (COPR). To evaluate the reasonableness of using urinary monitoring to assess environmental exposure, six volunteers ingested 400 mg of soil/day (low-dose group), two others ingested 2.0 g of soil/day (high-dose group) for 3 consecutive days, and one person ingested a placebo on each of 3 days. The soil and COPR mixture contained concentrations of total chromium (Cr) and hexavalent chromium [Cr(VI)] of 103 ± 20 and 9.3 ± 3.8 mg/kg, respectively. Therefore, the low-dose group ingested 41 μg Cr/day [including 3.7 μg Cr(VI)] and the high-dose group ingested 206 μg Cr/day [including 18.6 μg Cr(VI)] on each of 3 consecutive days. All urine samples were collected and analyzed individually for total Cr on the day prior to dosing, during the 3 days of dosing, and up to the first void 48 h after the last dose. No significant increases in urinary Cr excretion were found when background excretion data were compared with data following each of the 3 days of dosing or in daily mean urine concentrations of the high- vs the low-dose groups. It appears that Cr present in a soil and COPR mixture at Cr doses up to 200 μg/day is not sufficiently bioavailable for biomonitoring of urine to be informative. These results are consistent with previously published findings suggesting that incidental exposure to dusts and soils containing comparable levels of Cr will not result in increased concentrations of Cr in urine.     

Current Views on the Oral Bioavailability of Inorganic Mercury in Soil: Implications for Health Risk Assessments

Due to the presence of mercury at a number of major contaminated sites in the United States, the bioavailability of inorganic mercury in soil following ingestion has emerged as an important public health issue. Studies of the leachability/solubility of inorganic mercury in soil have shown that it is largely immobile, thereby suggesting that it will not be readily available for absorption in the gastrointestinal tract. Ignoring the effect of the soil matrix on decreasing bioavailability may result in a substantial overprediction of risks due to ingestion of contaminated soil. This paper discusses current knowledge about the oral bioavailability of inorganic mercury in soil and offers suggestions about how these data may be applied in human health risk assessment. Though precise estimates are not available, in vivo and in vitro estimates of the bioavailability of different inorganic mercury species in different matrices suggest that the bioavailability of mercury in soil is likely to be significantly less, on the order of at least three- to tenfold, than the bioavailability of mercuric chloride, the species used to derive the toxicity criteria for inorganic mercury. Because bioavailability can vary significantly with soil type, soil aging, the presence of co-contaminants and other factors, it is suggested that whenever the fiscal aspects justify a more precise estimate of bioavailability, site-specific estimates be developed. To develop a database for identifying a less expensive and more efficient method for estimating bioavailability, it is suggested that in vivo studies be conducted concurrently with in vitro studies. However, due to the lack of precision associated with the derivation of the most widely-used health guidance value for inorganic mercury (the US EPA RfD), additional work to address the uncertainties in the RfD is recommended.     

A Case Study Comparing Static and Spatially Explicit Ecological Exposure Analysis Methods

Exposure to chemical contaminants must be estimated when performing ecological risk assessments. A previous article proposed a habitat area and quality conditioned population exposure estimator, E[HQ]P, and described an individual-based, random walk, Monte Carlo model (SE3M) to facilitate calculation of E[HQ]P. In this article, E[HQ]P was compared with exposure estimates from a baseline risk assessment that evaluated mink and great blue heron exposure to fluoride at a federal Superfund site. Calculation of E[HQ]P took into consideration a receptor's forage area, movement behavior, population size, and the areal extent and quality of suitable habitat. The baseline assessment used four methods that did (total and unit Tier 2) and did not (total and unit Tier 1) consider habitat area or quality; where "total" included all exposure units on site and "unit" only a given exposure unit. Total Tier 1 estimates were consistently higher than E[HQ]P (e.g., 169.1 mg/kg x d versus 21.6 mg/kg x d). Risk managers using total Tier 1 results for decision making would be unlikely to underestimate exposure; however, implementability of correspondingly lower remedial objectives could be challenging. Unit Tier 1 estimates were higher (e.g., 96.5 mg/kg x d versus 61.6 mg/kg x d) or lower (e.g., 3.5 mg/kg x d versus 51.1 mg/kg x d) than E[HQ]P depending on variations in landscape features. Total Tier 2 and E[HQ]P estimates were similar (e.g., 20.7 mg/kg x d versus 21.6 mg/kg x d) when an ecologically questionable average exposure was assumed. Unit Tier 2 estimates were consistently well below E[HQ]P (e.g., 17.8 mg/kg x d versus 61.6 mg/kg x d) when an average exposure was not assumed. Risk managers using unit Tier 1 or 2 results could be basing their decisions on potentially large underestimates of exposure. By forgoing average exposure assumptions, and explicitly addressing landscape heterogeneity, SE3M appears capable of yielding exposure estimates that are not as potentially misleading to risk managers as those produced with traditional averaging methods.     

Bioaccessibility of Fukushima‐Accident‐Derived Cs in Soils and the Contribution of Soil Ingestion to Radiation Doses in Children

Ingestion of contaminated soil is one potential internal exposure pathway in areas contaminated by the Fukushima Daiichi Nuclear Power Plant accident. Doses from this pathway can be overestimated if the availability of radioactive nuclides in soils for the gastrointestinal tract is not considered. The concept of bioaccessibility has been adopted to evaluate this availability based on in vitro tests. This study evaluated the bioaccessibility of radioactive cesium from soils via the physiologically‐based extraction test (PBET) and the extractability of those via an extraction test with 1 mol/L of hydrochloric acid (HCl). The bioaccessibility obtained in the PBET was 5.3% ± 1%, and the extractability in the tests with HCl was 16% ± 3%. The bioaccessibility was strongly correlated with the extractability. This result indicates the possibility that the extractability in HCl can be used as a good predictor of the bioaccessibility with PBET. In addition, we assessed the doses to children from the ingestion of soil via hand‐to‐mouth activity based on our PBET results using a probabilistic approach considering the spatial distribution of radioactive cesium in Date City in Fukushima Prefecture and the interindividual differences in the surveyed amounts of soil ingestion in Japan. The results of this assessment indicate that even if children were to routinely ingest a large amount of soil with relatively high contamination, the radiation doses from this pathway are negligible compared with doses from external exposure owing to deposited radionuclides in Fukushima Prefecture.     

Application of Geostatistics to Risk Assessment

Geostatistics offers two fundamental contributions to environmental contaminant exposure assessment: (1) a group of methods to quantitatively describe the spatial distribution of a pollutant and (2) the ability to improve estimates of the exposure point concentration by exploiting the geospatial information present in the data. The second contribution is particularly valuable when exposure estimates must be derived from small data sets, which is often the case in environmental risk assessment. This article addresses two topics related to the use of geostatistics in human and ecological risk assessments performed at hazardous waste sites: (1) the importance of assessing model assumptions when using geostatistics and (2) the use of geostatistics to improve estimates of the exposure point concentration (EPC) in the limited data scenario. The latter topic is approached here by comparing design-based estimators that are familiar to environmental risk assessors (e.g., Land's method) with geostatistics, a model-based estimator. In this report, we summarize the basics of spatial weighting of sample data, kriging, and geostatistical simulation. We then explore the two topics identified above in a case study, using soil lead concentration data from a Superfund site (a skeet and trap range). We also describe several areas where research is needed to advance the use of geostatistics in environmental risk assessment.     

Development of a Standard Soil-to-Skin Adherence Probability Density Function for Use in Monte Carlo Analyses of Dermal Exposure

It has recently been suggested that "standard" data distributions for key exposure variables should be developed wherever appropriate for use in probabilistic or "Monte Carlo" exposure analyses. Soil-on-skin adherence estimates represent an ideal candidate for development of a standard data distribution: There are several readily available studies which offer a consistent pattern of reported results, and more importantly, soil adherence to skin is likely to vary little from site-to-site. In this paper, we thoroughly review each of the published soil adherence studies with respect to study design, sampling, and analytical methods, and level of confidence in the reported results. Based on these studies, probability density functions (PDF) of soil adherence values were examined for different age groups and different sampling techniques. The soil adherence PDF developed from adult data was found to resemble closely the soil adherence PDF based on child data in terms of both central tendency (mean = 0.49 and 0.63 mg-soil/cm²-skin, respectively) and 95th percentile values (1.6 and 2.4 mg-soil/cm²-skin, respectively). Accordingly, a single, "standard" PDF is presented based on all data collected for all age groups. This standard PDF is lognormally distributed; the arithmetic mean and standard deviation are 0.52 ± 0.9 mg-soil/cm²-skin. Since our review of the literature indicates that soil adherence under environmental conditions will be minimally influenced by age, sex, soil type, or particle size, this PDF should be considered applicable to all settings. The 50th and 95th percentile values of the standard PDF (0.25 and 1.7 mg-soil/cm²-skin, respectively) are very similar to recent U.S. EPA estimates of "average" and "upper-bound" soil adherence (0.2 and 1.0 mg-soil/cm²-skin, respectively).     

Use of Monte Carlo Simulation for Human Exposure Assessment at a Superfund Site

This work presents a comparison of probabilistic and deterministic health risk estimates based on data from an industrial site in the northeastern United States. The risk assessment considered exposures to volatile solvents by drinking water ingestion and showering. Probability densities used as inputs included concentrations, contact rates, and exposure frequencies; dose-response inputs were single values. Deterministic risk estimates were calculated by the "reasonable maximum exposure" (RME) approach recommended by the EPA Superfund program. The RME non-carcinogenic risk fell between the 90th and the 95th percentile of the probability density; the RME cancer risk fell between the 95th percentile and the maximum. These results suggest that in this case (1) EPA's deterministic RME risk was reasonably protective, (2) results of probabilistic and deterministic calculations were consistent, and (3) commercially available software Monte Carlo software effectively provided multiple risk estimates recommended by recent EPA guidance.     

Assessing Risk‐Based Upper Limits of Melamine Migration from Food Containers

Melamine contamination of food has become a major food safety issue because of incidents of infant disease caused by exposure to this chemical. This study was aimed at establishing a safety limit in Taiwan for the degree of melamine migration from food containers. Health risk assessment was performed for three exposure groups (preschool children, individuals who dine out, and elderly residents of nursing homes). Selected values of tolerable daily intake (TDI) for melamine were used to calculate the reference migration concentration limit (RMCL) or reference specific migration limit (RSML) for melamine food containers. The only existing values of these limits for international standards today are 1.2 mg/L (0.2 mg/dm²) in China and 30 mg/L (5 mg/dm²) in the European Union. The factors used in the calculations included the specific surface area of food containers, daily food consumption rate, body weight, TDI, and the percentile of the population protected at a given migration concentration limit (MCL). The results indicate that children are indeed at higher risk of melamine exposure at toxic levels than are other groups and that the 95th percentile of MCL (specific surface area = 5) for children aged 1–6 years should be the RMCL (0.07 mg/dm²) for protecting the sensitive and general population.     

Explosion Probability of Unexploded Ordnance: Expert Beliefs

This article reports on a study to quantify expert beliefs about the explosion probability of unexploded ordnance (UXO). Some 1,976 sites at closed military bases in the United States are contaminated with UXO and are slated for cleanup, at an estimated cost of $15–140 billion. Because no available technology can guarantee 100% removal of UXO, information about explosion probability is needed to assess the residual risks of civilian reuse of closed military bases and to make decisions about how much to invest in cleanup. This study elicited probability distributions for the chance of UXO explosion from 25 experts in explosive ordnance disposal, all of whom have had field experience in UXO identification and deactivation. The study considered six different scenarios: three different types of UXO handled in two different ways (one involving children and the other involving construction workers). We also asked the experts to rank by sensitivity to explosion 20 different kinds of UXO found at a case study site at Fort Ord, California. We found that the experts do not agree about the probability of UXO explosion, with significant differences among experts in their mean estimates of explosion probabilities and in the amount of uncertainty that they express in their estimates. In three of the six scenarios, the divergence was so great that the average of all the expert probability distributions was statistically indistinguishable from a uniform (0, 1) distribution—suggesting that the sum of expert opinion provides no information at all about the explosion risk. The experts' opinions on the relative sensitivity to explosion of the 20 UXO items also diverged. The average correlation between rankings of any pair of experts was 0.41, which, statistically, is barely significant (p= 0.049) at the 95% confidence level. Thus, one expert's rankings provide little predictive information about another's rankings. The lack of consensus among experts suggests that empirical studies are needed to better understand the explosion risks of UXO.     

Estimating Consumer Exposure to PFOS and PFOA

Perfluorinated compounds have been used for more than 50 years as process aids, surfactants, and for surface protection. This study is a comprehensive assessment of consumer exposure to perfluorooctane sulfonate (PFOS) and perfluorooctanoic acid (PFOA) from a variety of environmental and product-related sources. To identify relevant pathways leading to consumer exposure to PFOS and PFOA a scenario-based approach has been applied. Scenarios represent realistic situations where age- and gender-specific exposure occurs in the everyday life of consumers. We find that North American and European consumers are likely to experience ubiquitous and long-term uptake doses of PFOS and PFOA in the range of 3 to 220 ng per kg body weight per day (ng/kg_bw/day) and 1 to 130 ng/kg_bw/day, respectively. The greatest portion of the chronic exposure to PFOS and PFOA is likely to result from the intake of contaminated foods, including drinking water. Consumer products cause a minor portion of the consumer exposure to PFOS and PFOA. Of these, it is mainly impregnation sprays, treated carpets in homes, and coated food contact materials that may lead to consumer exposure to PFOS and PFOA. Children tend to experience higher total uptake doses (on a body weight basis) than teenagers and adults because of higher relative uptake via food consumption and hand-to-mouth transfer of chemical from treated carpets and ingestion of dust. The uptake estimates based on scenarios are within the range of values derived from blood serum data by applying a one-compartment pharmacokinetic model.     

Improving Risk Assessment Calculations for Traditional Foods Through Collaborative Research with First Nations Communities

As industrial development is increasing near northern Canadian communities, human health risk assessments (HHRA) are conducted to assess the predicted magnitude of impacts of chemical emissions on human health. One exposure pathway assessed for First Nations communities is the consumption of traditional plants, such as muskeg tea (Labrador tea) (Ledum/Rhododendron groenlandicum) and mint (Mentha arvensis). These plants are used to make tea and are not typically consumed in their raw form. Traditional practices were used to harvest muskeg tea leaves and mint leaves by two First Nations communities in northern Alberta, Canada. Under the direction of community elders, community youth collected and dried plants to make tea. Soil, plant, and tea decoction samples were analyzed for inorganic elements using inductively coupled plasma‐mass spectrometry. Concentrations of inorganic elements in the tea decoctions were orders of magnitude lower than in the vegetation (e.g., manganese 0.107 mg/L in tea, 753 mg/kg in leaves). For barium, the practice of assessing ingestion of raw vegetation would have resulted in a hazard quotient (HQ) greater than the benchmark of 0.2. Using measured tea concentrations it was determined that exposure would result in risk estimates orders of magnitude below the HQ benchmark of 0.2 (HQ = 0.0049 and 0.017 for muskeg and mint tea, respectively). An HHRA calculating exposure to tea vegetation through direct ingestion of the leaves may overestimate risk. The results emphasize that food preparation methods must be considered when conducting an HHRA. This study illustrates how collaboration between Western scientists and First Nations communities can add greater clarity to risk assessments.