Exposure to Persistent Organochlorines in Canadian Breast Milk: A Probabilistic Assessment

Exposure to persistent organochlorines in breast milk was estimated probabilistically for Canadian infants. Noncancer health effects were evaluated by comparing the predicted exposure distributions to published guidance values. For chemicals identified as potential human carcinogens, cancer risks were evaluated using standard methodology typically applied in Canada, as well as an alternative method developed under the Canadian Environmental Protection Act. Potential health risks associated with exposure to persistent organochlorines were quantitatively and qualitatively weighed against the benefits of breast-feeding. Current levels of the majority of contaminants identified in Canadian breast milk do not pose unacceptable risks to infants. Benefits of breast-feeding are well documented and qualitatively appear to outweigh potential health concerns associated with organochlorine exposure. Furthermore, the risks of mortality from not breast-feeding estimated by Rogan and colleagues exceed the theoretical cancer risks estimated for infant exposure to potential carcinogens in Canadian breast milk. Although levels of persistent compounds have been declining in Canadian breast milk, potentially significant risks were estimated for exposure to polychlorinated biphenyls, dibenzo-p-dioxins, and dibenzofurans. Follow-up work is suggested that would involve the use of a physiologically based toxicokinetic model with probabilistic inputs to predict dioxin exposure to the infant. A more detailed risk analysis could be carried out by coupling the exposure estimates with a dose–response analysis that accounts for uncertainty.     

Bounding Analysis of Drinking Water Health Risks from a Spill of Hydraulic Fracturing Flowback Water

A bounding risk assessment is presented that evaluates possible human health risk from a hypothetical scenario involving a 10,000‐gallon release of flowback water from horizontal fracturing of Marcellus Shale. The water is assumed to be spilled on the ground, infiltrates into groundwater that is a source of drinking water, and an adult and child located downgradient drink the groundwater. Key uncertainties in estimating risk are given explicit quantitative treatment using Monte Carlo analysis. Chemicals that contribute significantly to estimated health risks are identified, as are key uncertainties and variables to which risk estimates are sensitive. The results show that hypothetical exposure via drinking water impacted by chemicals in Marcellus Shale flowback water, assumed to be spilled onto the ground surface, results in predicted bounds between 10⁻¹⁰ and 10⁻⁶ (for both adult and child receptors) for excess lifetime cancer risk. Cumulative hazard indices (HI_CUMULATIVE) resulting from these hypothetical exposures have predicted bounds (5th to 95th percentile) between 0.02 and 35 for assumed adult receptors and 0.1 and 146 for assumed child receptors. Predicted health risks are dominated by noncancer endpoints related to ingestion of barium and lithium in impacted groundwater. Hazard indices above unity are largely related to exposure to lithium. Salinity taste thresholds are likely to be exceeded before drinking water exposures result in adverse health effects. The findings provide focus for policy discussions concerning flowback water risk management. They also indicate ways to improve the ability to estimate health risks from drinking water impacted by a flowback water spill (i.e., reducing uncertainty).     

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.     

Uncertain Numbers and Uncertainty in the Selection of Input Distributions—Consequences for a Probabilistic Risk Assessment of Contaminated Land

Risks from exposure to contaminated land are often assessed with the aid of mathematical models. The current probabilistic approach is a considerable improvement on previous deterministic risk assessment practices, in that it attempts to characterize uncertainty and variability. However, some inputs continue to be assigned as precise numbers, while others are characterized as precise probability distributions. Such precision is hard to justify, and we show in this article how rounding errors and distribution assumptions can affect an exposure assessment. The outcome of traditional deterministic point estimates and Monte Carlo simulations were compared to probability bounds analyses. Assigning all scalars as imprecise numbers (intervals prescribed by significant digits) added uncertainty to the deterministic point estimate of about one order of magnitude. Similarly, representing probability distributions as probability boxes added several orders of magnitude to the uncertainty of the probabilistic estimate. This indicates that the size of the uncertainty in such assessments is actually much greater than currently reported. The article suggests that full disclosure of the uncertainty may facilitate decision making in opening up a negotiation window. In the risk analysis process, it is also an ethical obligation to clarify the boundary between the scientific and social domains.     

Comparative Probabilistic Assessment of Occupational Pesticide Exposures Based on Regulatory Assessments

Implementation of probabilistic analyses in exposure assessment can provide valuable insight into the risks of those at the extremes of population distributions, including more vulnerable or sensitive subgroups. Incorporation of these analyses into current regulatory methods for occupational pesticide exposure is enabled by the exposure data sets and associated data currently used in the risk assessment approach of the Environmental Protection Agency (EPA). Monte Carlo simulations were performed on exposure measurements from the Agricultural Handler Exposure Database and the Pesticide Handler Exposure Database along with data from the Exposure Factors Handbook and other sources to calculate exposure rates for three different neurotoxic compounds (azinphos methyl, acetamiprid, emamectin benzoate) across four pesticide‐handling scenarios. Probabilistic estimates of doses were compared with the no observable effect levels used in the EPA occupational risk assessments. Some percentage of workers were predicted to exceed the level of concern for all three compounds: 54% for azinphos methyl, 5% for acetamiprid, and 20% for emamectin benzoate. This finding has implications for pesticide risk assessment and offers an alternative procedure that may be more protective of those at the extremes of exposure than the current approach.     

Mlonte Carlo Techniques for Quantitative Uncertainty Analysis in Public Health Risk Assessments

Most public health risk assessments assume and combine a series of average, conservative, and worst-case values to derive a conservative point estimate of risk. This procedure has major limitations. This paper demonstrates a new methodology for extended uncertainty analyses in public health risk assessments using Monte Carlo techniques. The extended method begins as do some conventional methods--with the preparation of a spreadsheet to estimate exposure and risk. This method, however, continues by modeling key inputs as random variables described by probability density functions (PDFs). Overall, the technique provides a quantitative way to estimate the probability distributions for exposure and health risks within the validity of the model used. As an example, this paper presents a simplified case study for children playing in soils contaminated with benzene and benzo(a)pyrene (BaP).     

Effects of Uncertainties on Exposure Estimates to Methylmercury: A Monte Carlo Analysis of Exposure Biomarkers versus Dietary Recall Estimation

This article presents a general model for estimating population heterogeneity and "lack of knowledge" uncertainty in methylmercury (MeHg) exposure assessments using two-dimensional Monte Carlo analysis. Using data from fish-consuming populations in Bangladesh, Brazil, Sweden, and the United Kingdom, predictive model estimates of dietary MeHg exposures were compared against those derived from biomarkers (i.e., [Hg]hair and [Hg]blood). By disaggregating parameter uncertainty into components (i.e., population heterogeneity, measurement error, recall error, and sampling error) estimates were obtained of the contribution of each component to the overall uncertainty. Steady-state diet:hair and diet:blood MeHg exposure ratios were estimated for each population and were used to develop distributions useful for conducting biomarker-based probabilistic assessments of MeHg exposure. The 5th and 95th percentile modeled MeHg exposure estimates around mean population exposure from each of the four study populations are presented to demonstrate lack of knowledge uncertainty about a best estimate for a true mean. Results from a U.K. study population showed that a predictive dietary model resulted in a 74% lower lack of knowledge uncertainty around a central mean estimate relative to a hair biomarker model, and also in a 31% lower lack of knowledge uncertainty around central mean estimate relative to a blood biomarker model. Similar results were obtained for the Brazil and Bangladesh populations. Such analyses, used here to evaluate alternative models of dietary MeHg exposure, can be used to refine exposure instruments, improve information used in site management and remediation decision making, and identify sources of uncertainty in risk estimates.     

Combining Food Frequency and Survey Data to Quantify Long-Term Dietary Exposure: A Methyl Mercury Case Study

Twenty-four-hour recall data from the Continuing Survey of Food Intake by Individuals (CSFII) are frequently used to estimate dietary exposure for risk assessment. Food frequency questionnaires are traditional instruments of epidemiological research; however, their application in dietary exposure and risk assessment has been limited. This article presents a probabilistic method of bridging the National Health and Nutrition Examination Survey (NHANES) food frequency and the CSFII data to estimate longitudinal (usual) intake, using a case study of seafood mercury exposures for two population subgroups (females 16 to 49 years and children 1 to 5 years). Two hundred forty-nine CSFII food codes were mapped into 28 NHANES fish/shellfish categories. FDA and state/local seafood mercury data were used. A uniform distribution with minimum and maximum blood-diet ratios of 0.66 to 1.07 was assumed. A probabilistic assessment was conducted to estimate distributions of individual 30-day average daily fish/shellfish intakes, methyl mercury exposure, and blood levels. The upper percentile estimates of fish and shellfish intakes based on the 30-day daily averages were lower than those based on two- and three-day daily averages. These results support previous findings that distributions of "usual" intakes based on a small number of consumption days provide overestimates in the upper percentiles. About 10% of the females (16 to 49 years) and children (1 to 5 years) may be exposed to mercury levels above the EPA's RfD. The predicted 75th and 90th percentile blood mercury levels for the females in the 16-to-49-year group were similar to those reported by NHANES. The predicted 90th percentile blood mercury levels for children in the 1-to-5-year subgroup was similar to NHANES and the 75th percentile estimates were slightly above the NHANES.     

Risk Criteria for Nuclear Power Plants: A Pragmatic Proposal1

Criteria are proposed for both an acceptable upper bound of nuclear power plant risk and a lower bound as a design target. Recognizing that the public risk associated with a power plant can be estimated only by probabilistic analysis of the design features, the spread between the lower design target and the upper bound provides a margin for uncertainty in th probabilistic estimate. The combination of a low probabilistic design target and this margin provides a reasonable expectation that the overall performance will be in the domain of an acceptable risk level. Because the exposure to potential risk is chiefly in the locality of the nuclear station, it is also proposed that compensatory benefits should be provided locally and that these be included as a cost of operation. It is suggested that the upper bound be set at a risk level equivalent to those risks of routine living which are normally accepted, i.e., about 10^-4 deaths per year per person (100 deaths/yr/million). The proposed lower design target is 10^-8 (0.1 deaths/yr/million), about one-hundredth of the minimal risk from the natural hazards all people are exposed to.     

Variability in Physical Constant Parameter Values from Standard Data Sources and the Implication of this Variability for Risk Analysis

Physical property values are used in environmental risk assessments to estimate media and risk-based concentrations. Recently, however, considerable variability has been reported with such values. To evaluate potential variability in physical parameter values supporting a variety of regulatory programs, eight data sources were chosen for evaluation, and chemicals appearing in at least four sources were selected. There were 755 chemicals chosen. In addition, chemicals in seven environmentally important subgroups were also identified for evaluation. Nine parameters were selected for analysis--molecular weight (MolWt), melting point (MeltPt), boiling point (BoilPt), vapor pressure (VP), water solubility (AqSOL), Henry's law constant (HLC), octanol-water partition coefficient (Kow), and diffusion coefficients in air (Dair) and water (Dwater). Results show that while 71% of constituents had equal MolWts across data sources, <3% of the constituents had equivalent parameter values across data sources for AqSOL, VP, or HLC. Considerable dissimilarity between certain sources was also observed. Furthermore, measures of dispersion showed considerable variation in data sets for Kow, VP, AqSOL, and HLC compared to measures for MolWt, MeltPt, BoilPt, or Dwater. The magnitude of the observed variability was also noteworthy. For example, the 95th percentile ratio of maximum/minimum parameter values ranged from 1.0 for MolWt to well over 1.0 x 10(6) for VP and HLC. Risk and exposure metrics also varied by similar magnitudes. Results with environmentally important subgroups were similar. These results show that there is considerable variability in physical parameter values from standard sources, and that the observed variability could affect potential risk estimates and perhaps risk management decisions.     

Variability in Physical Constant Parameter Values from Standard Data Sources and the Implication of This Variability for Risk Analysis

Physical property values are used in environmental risk assessments to estimate media and risk-based concentrations. However, considerable variability has recently been reported with such values. To evaluate potential variability in physical parameter values supporting a variety of regulatory programs, eight data sources were chosen for evaluation, and chemicals appearing in at least four sources were selected. There were 755 chemicals chosen. In addition, chemicals in seven environmentally important subgroups were also identified for evaluation. Nine parameters were selected for analysis-molecular weight (MolWt), melting point (MeltPt), boiling point (BoilPt), vapor pressure (VP), water solubility (AqSOL), Henry's law constant (HLC), octanol-water partition coefficient (Kow), and diffusion coefficients in air (Dair) and water (Dwater). Results show that while 71% of constituents had equal MolWts across data sources, <3% of the constituents had equivalent parameter values across data sources for AqSOL, VP, or HLC. Considerable dissimilarity between certain sources was also observed. Furthermore, measures of dispersion showed considerable variation in data sets for Kow, VP, AqSOL, and HLC compared to measures for MolWt, MeltPt, BoilPt, or Dwater. The magnitude of the observed variability was also noteworthy. For example, the 95th percentile ratio of maximum/minimum parameter values ranged from 1.0 for MolWt to well over 1.0E + 06 for VP, and HLC. Risk and exposure metrics also varied by similar magnitudes. Results with environmentally important subgroups were similar. These results show that there is considerable variability in physical parameter values from standard sources, and that the observed variability could affect potential risk estimates and perhaps risk management decisions.     

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.     

Plausible Upper Bounds: Are Their Sums Plausible?1

Quantitative cancer risk assessments are typically expressed as plausible upper bounds rather than estimates of central tendency. In analyses involving several carcinogens, these upper bounds are often summed to estimate overall risk. This begs the question of whether a sum of upper bounds is itself a plausible estimate of overall risk. This question can be asked in two ways: whether the sum yields an improbable estimate of overall risk (that is, is it only remotely possible for the true sum of risks to match the sum of upper bounds), or whether the sum gives a misleading estimate (that is, is the true sum of risks likely to be very different from the sum of upper bounds). Analysis of four case studies shows that as the number of risk estimates increases, their sum becomes increasingly improbable, but not misleading. Though the overall risk depends on the independence, additivity, and number of risk estimates, as well as the shapes of the underlying risk distributions, sums of upper bounds provide useful information about the overall risk and can be adjusted downward to give a more plausible [perhaps probable] upper bound, or even a central estimate of overall risk.     

Does EPA Underestimate Cancer Risks by Ignoring Susceptibility Differences?

A 2009 report of the National Research Council (NRC) recommended that the U.S. Environmental Protection Agency (EPA) increase its estimates of increased cancer risk from exposure to environmental agents by ～7‐fold, due to an approximate ～25‐fold typical ratio between the median and upper 95th percentile persons’ cancer sensitivity assuming approximately lognormally distributed sensitivities. EPA inaction on this issue has raised concerns that cancer risks to environmentally exposed populations remain systematically underestimated. This concern is unwarranted, however, because EPA point estimates of cancer risk have always pertained to the average, not the median, person in each modeled exposure group. Nevertheless, EPA has yet to explain clearly how its risk characterization and risk management policies concerning individual risks from environmental chemical carcinogens do appropriately address broad variability in human cancer susceptibility that has been a focus of two major NRC reports to EPA concerning its risk assessment methods.     

Exposure Assessment: Then, Now, and Quantum Leaps in the Future

Health risk assessments have become so widely accepted in the United States that their conclusions are a major factor in many environmental decisions. Although the risk assessment paradigm is 10 years old, the basic risk assessment process has been used by certain regulatory agencies for nearly 40 years. Each of the four components of the paradigm has undergone significant refinements, particularly during the last 5 years. A recent step in the development of the exposure assessment component can be found in the 1992 EPA Guidelines for Exposure Assessment. Rather than assuming worst-case or hypothetical maximum exposures, these guidelines are designed to lead to an accurate characterization, making use of a number of scientific advances. Many exposure parameters have become better defined, and more sensitive techniques now exist for measuring concentrations of contaminants in the environnment. Statistical procedures for characterizing variability, using Monte Carlo or similar approaches, eliminate the need to select point estimates for all individual exposure parameters. These probabilistic models can more accurately characterize the full range of exposures that may potentially be encountered by a given population at a particular site, reducing the need to select highly conservative values to account for this form of uncertainty in the exposure estimate. Lastly, our awareness of the uncertainties in the exposure assessment as well as our knowledge as to how best to characterize them will almost certainly provide evaluations that will be more credible and, therein, more useful to risk managers. If these refinements are incorporated into future exposure assessments, it is likely that our resources will be devoted to problems that, when resolved, will yield the largest improvement in public health.     

Measures of Compounding Conservatism in Probabilistic Risk Assessment

Concern about the degree of uncertainty and potential conservatism in deterministic point estimates of risk has prompted researchers to turn increasingly to probabilistic methods for risk assessment. With Monte Carlo simulation techniques, distributions of risk reflecting uncertainty and/or variability are generated as an alternative. In this paper the compounding of conservatism⁽¹⁾ between the level associated with point estimate inputs selected from probability distributions and the level associated with the deterministic value of risk calculated using these inputs is explored. Two measures of compounded conservatism are compared and contrasted. The first measure considered, F , is defined as the ratio of the risk value, R _d, calculated deterministically as a function of n inputs each at the j th percentile of its probability distribution, and the risk value, R _j that falls at the j th percentile of the simulated risk distribution (i.e., F=R_d/R_j). The percentile of the simulated risk distribution which corresponds to the deterministic value, R_d , serves as a second measure of compounded conservatism. Analytical results for simple products of lognormal distributions are presented. In addition, a numerical treatment of several complex cases is presented using five simulation analyses from the literature to illustrate. Overall, there are cases in which conservatism compounds dramatically for deterministic point estimates of risk constructed from upper percentiles of input parameters, as well as those for which the effect is less notable. The analytical and numerical techniques discussed are intended to help analysts explore the factors that influence the magnitude of compounding conservatism in specific cases.     

Quantifying Uncertainty in a Risk Assessment Using Human Data

A call for risk assessment approaches that better characterize and quantify uncertainty has been made by the scientific and regulatory community. This paper responds to that call by demonstrating a distributional approach that draws upon human data to derive potency estimates and to identify and quantify important sources of uncertainty. The approach is rooted in the science of decision analysis and employs an influence diagram, a decision tree, probabilistic weights, and a distribution of point estimates of carcinogenic potency. Its results estimate the likelihood of different carcinogenic risks (potencies) for a chemical under a specific scenario. For this exercise, human data on formaldehyde were employed to demonstrate the approach. Sensitivity analyses were performed to determine the relative impact of specific levels and alternatives on the potency distribution. The resulting potency estimates are compared with the results of an exercise using animal data on formaldehyde. The paper demonstrates that distributional risk assessment is readily adapted to situations in which epidemiologic data serve as the basis for potency estimates. Strengths and weaknesses of the distributional approach are discussed. Areas for further application and research are recommended.     

Recommended Distributions for Exposure Factors Frequently Used in Health Risk Assessment

Although there has been nearly complete agreement in the scientific community that Monte Carlo techniques represent a significant improvement in the exposure assessment process, virtually all state and federal risk assessments still rely on the traditional point estimate approach. One of the rate-determining steps to a timely implementation of Monte Carlo techniques to regulatory decision making is the development of "standard" data distributions that are considered applicable to any setting. For many exposure variables, there is no need to wait any longer to adopt Monte Carlo techniques into regulatory policy since there is a wealth of data from which a robust distribution can be developed and ample evidence to indicate that the variable is not significantly influenced by site-specific conditions. In this paper, we propose several distributions that can be considered standard and customary for most settings. Age-specific distributions for soil ingestion rates, inhalation rates, body weights, skin surface area, tapwater and fish consumption, residential occupancy and occupational tenure, and soil-on-skin adherence were developed. For each distribution offered in this paper, we discuss the adequacy of the database, derivation of the distribution, and applicability of the distribution to various settings and conditions.     

Integrated Analysis: Combining Risk and Economic Assessments While Preserving the Separation of Powers

This article presents a process for an integrated policy analysis that combines risk assessment and benefit-cost analysis. This concept, which explicitly combines the two types of related analyses, seems to contradict the long-accepted risk analysis paradigm of separating risk assessment and risk management since benefit-cost analysis is generally considered to be a part of risk management. Yet that separation has become a problem because benefit-cost analysis uses risk assessment results as a starting point and considerable debate over the last several years focused on the incompatibility of the use of upper bounds or "safe" point estimates in many risk assessments with benefit-cost analysis. The problem with these risk assessments is that they ignore probabilistic information. As advanced probabilistic techniques for risk assessment emerge and economic analysts receive distributions of risks instead of point estimates, the artificial separation between risk analysts and the economic/decision analysts complicates the overall analysis. In addition, recent developments in countervailing risk theory suggest that combining the risk and benefit-cost analyses is required to fully understand the complexity of choices and tradeoffs faced by the decisionmaker. This article also argues that the separation of analysis and management is important, but that benefit-cost analysis has been wrongly classified into the risk management category and that the analytical effort associated with understanding the economic impacts of risk reduction actions need to be part of a broader risk assessment process.