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.     

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).     

On the Effect of Probability Distributions of Input Variables in Public Health Risk Assessment

A central part of probabilistic public health risk assessment is the selection of probability distributions for the uncertain input variables. In this paper, we apply the first-order reliability method (FORM)^(1–3) as a probabilistic tool to assess the effect of probability distributions of the input random variables on the probability that risk exceeds a threshold level (termed the probability of failure) and on the relevant probabilistic sensitivities. The analysis was applied to a case study given by Thompson et al. ⁽⁴⁾ on cancer risk caused by the ingestion of benzene contaminated soil. Normal, lognormal, and uniform distributions were used in the analysis. The results show that the selection of a probability distribution function for the uncertain variables in this case study had a moderate impact on the probability that values would fall above a given threshold risk when the threshold risk is at the 50th percentile of the original distribution given by Thompson et al. ⁽⁴⁾ The impact was much greater when the threshold risk level was at the 95th percentile. The impact on uncertainty sensitivity, however, showed a reversed trend, where the impact was more appreciable for the 50th percentile of the original distribution of risk given by Thompson et al. ⁴ than for the 95th percentile. Nevertheless, the choice of distribution shape did not alter the order of probabilistic sensitivity of the basic uncertain variables.     

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.     

Residential House and Zone Volumes in the United States: Empirical and Estimated Parametric Distributions

This research reports empirical distributions and estimated univariate parametric probability distributions for house volume and certain zone volumes within households for residential structures in the United States. The author derived the distributions from two separate databases. The volumes were found to be exceptionally well fit by lognormal distributions (adjusted R² >> 0.95) in almost all cases. In addition, data from one database indicates that the correlation between house volume and air changes per hour is very weak.     

First-Order Reliability Analysis of Public Health Risk Assessment

This paper demonstrates a new methodology for probabilistic public health risk assessment using the first-order reliability method. The method provides the probability that incremental lifetime cancer risk exceeds a threshold level, and the probabilistic sensitivity quantifying the relative impact of considering the uncertainty of each random variable on the exceedance probability. The approach is applied to a case study given by Thompson et al. ⁽¹⁾ on cancer risk caused by ingestion of benzene-contaminated soil, and the results are compared to that of the Monte Carlo method. Parametric sensitivity analyses are conducted to assess the sensitivity of the probabilistic event with respect to the distribution parameters of the basic random variables, such as the mean and standard deviation. The technique is a novel approach to probabilistic risk assessment, and can be used in situations when Monte Carlo analysis is computationally expensive, such as when the simulated risk is at the tail of the risk probability distribution.     

Correlated Inputs in Quantitative Risk Assessment: The Effects of Distributional Shape

Application of Monte Carlo simulation methods to quantitative risk assessment are becoming increasingly popular. With this methodology, investigators have become concerned about correlations among input variables which might affect the resulting distribution of risk. We show that the choice of input distributions in these simulations likely has a larger effect on the resultant risk distribution than does the inclusion or exclusion of correlations. Previous investigators have studied the effect of correlated input variables for the addition of variables with any underlying distribution and for the product of lognormally distributed variables. The effects in the main part of the distribution are small unless the correlation and variances are large. We extend this work by considering addition, multiplication and division of two variables with assumed normal, lognormal, uniform and triangular distributions. For all possible pairwise combinations, we find that the effects of correlated input variables are similar to those observed for lognormal distributions, and thus relatively small overall. The effect of using different distributions, however, can be large.     

Parametric Distributions for Soil Ingestion by Children

This note presents parameterized distributions of estimates of the amount of soil ingested by children based on data collected by Binder et al. (1986). Following discussions with Dr. Binder, we modified the Binder study data by using the actual stool weights instead of the 15 g value used in the original study. After testing the data for lognormality, we generated parameterized distributions for use in risk assessment uncertainty analyses such as Monte Carlo simulations.     

Quantifying the Distribution of Inhalation Exposure in Human Populations: 2. Distributions of Time Spent by Adults, Adolescents, and Children at Home, at Work, and at School

Using distributions of time spent at various ventilation levels, ranges of inhalation exposure in the population can be established. Distributions of exposure time were determined using results of a study by the California Air Resources Board (CARB) which focused on time spent by humans participating in various activities and the locations where the activities occurred. The daily at-home activities from the CARB study were assigned to one of three ventilation levels, generating aggregate time periods. Distinct age and gender populations were identified, and distributions for aggregate time were established for these populations at each of the ventilation levels. In addition to aggregate time spent at home, distributions for various ages and genders were established for aggregate time spent at school and work. By combining distributions of aggregate time with corresponding ventilation rates, the distribution of inhalation rates can be established for at home, at work, and at school exposures.     

Estimated Distributions for Average Daily Consumption of Total and Self-Caught Fish for Adults in Michigan Angler Households

Fish consumption rates play a critical role in the assessment of human health risks posed by the consumption of fish from chemically contaminated water bodies. Based on data from the 1989 Michigan Sport Anglers Fish Consumption Survey, we examined total fish consumption, consumption of self-caught fish, and consumption of Great Lakes fish for all adults, men, women, and certain higher risk subgroups such as anglers. We present average daily consumption rates as compound probability distributions consisting of a Bernoulli trial (to distinguish those who ate fish from those who did not) combined with a distribution (both empirical and parametric) for those who ate fish. We found that the average daily consumption rates for adults who ate fish are reasonably well fit by lognormal distributions. The compound distributions may be used as input variables for Monte Carlo simulations in public health risk assessments.     

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.     

Lognormal Distributions for Water Intake by Children and Adults

We fit lognormal distributions to data collected in a national survey for both total water intake and tap water intake by children and adults for these age groups in years: 0 less than age less than 1; 1 less than or equal to age less than 11; 11 less than or equal to age less than 20; 20 less than or equal to age less than 65; 65 less than or equal to age; and all people in the survey taken as a single group. These distributions are suitable for use in public health risk assessments.     

Residential Air Exchange Rates in the United States: Empirical and Estimated Parametric Distributions by Season and Climatic Region

The purpose of this paper is to undertake a statistical analysis to specify empirical distributions and to estimate univariate parametric probability distributions for air exchange rates for residential structures in the United States. To achieve this goal, we used data compiled by the Brookhaven National Laboratory using a method known as the perfluorocarbon tracer (PFT) technique. While these data are not fully representative of all areas of the country or all housing types, they are judged to be by far the best available. The analysis is characterized by four key points: the use of data for 2,844 households; a four-region breakdown based on heating degree days, a best available measure of climatic factors affecting air exchange rates; estimation of lognormal distributions as well as provision of empirical (frequency) distributions; and provision of these distributions for all of the data, for the data segmented by the four regions, for the data segmented by the four seasons, and for the data segmented by a 16 region by season breakdown. Except in a few cases, primarily for small sample sizes, air exchange rates were found to be well fit by lognormal distributions (adjusted R² 0.95). The empirical or lognormal distributions may be used in indoor air models or as input variables for probabilistic human health risk assessments.     

A Trichloroethylene Risk Assessment Using a Monte Carlo Analysis of Parameter Uncertainty in Conjunction with Physiologically-Based Pharmacokinetic Modeling

A Monte Carlo simulation is incorporated into a risk assessment for trichloroethylene (TCE) using physiologically-based pharmacokinetic (PBPK) modeling coupled with the linearized multistage model to derive human carcinogenic risk extrapolations. The Monte Carlo technique incorporates physiological parameter variability to produce a statistically derived range of risk estimates which quantifies specific uncertainties associated with PBPK risk assessment approaches. Both inhalation and ingestion exposure routes are addressed. Simulated exposure scenarios were consistent with those used by the Environmental Protection Agency (EPA) in their TCE risk assessment. Mean values of physiological parameters were gathered from the literature for both mice (carcinogenic bioassay subjects) and for humans. Realistic physiological value distributions were assumed using existing data on variability. Mouse cancer bioassay data were correlated to total TCE metabolized and area-under-the-curve (blood concentration) trichloroacetic acid (TCA) as determined by a mouse PBPK model. These internal dose metrics were used in a linearized multistage model analysis to determine dose metric values corresponding to 10^-6 lifetime excess cancer risk. Using a human PBPK model, these metabolized doses were then extrapolated to equivalent human exposures (inhalation and ingestion). The Monte Carlo iterations with varying mouse and human physiological parameters produced a range of human exposure concentrations producing a 10^-6 risk.     

Scenario-Based Risk Assessment of Multi-Use Chemicals: Application to Solvents

The risk through chemical exposure is commonly characterized by ratios of exposure concentrations and effect levels (risk quotients). For chemicals with many different applications such as solvents, however, in addition to the risk quotients of different exposure situations it is useful to determine the corresponding numbers of exposed individuals, that is, not only the magnitude but also the extent of the risk. To this end, the Scenario-Based Risk Assessment (SceBRA) method has been developed that makes use of a large set of scenarios, each of which describes a typical situation regarding handling a solvent or solvent-containing product. The scenarios cover the life-cycle steps of production, distribution, and use of solvents. For each scenario, SceBRA provides the risk quotient, r, and the number of exposed individuals, N. This study investigated seven solvents that are used in large amounts in Switzerland. For each solvent, characteristic distributions of r and N values were calculated, making it possible to compare different solvents with respect to their risk profile. Graphical representations of the r, N data provide an informative way for analyzing and communicating the results of SceBRA.     

Sensitivity Analysis, Monte Carlo Risk Analysis, and Bayesian Uncertainty Assessment

Standard statistical methods understate the uncertainty one should attach to effect estimates obtained from observational data. Among the methods used to address this problem are sensitivity analysis, Monte Carlo risk analysis (MCRA), and Bayesian uncertainty assessment. Estimates from MCRAs have been presented as if they were valid frequentist or Bayesian results, but examples show that they need not be either in actual applications. It is concluded that both sensitivity analyses and MCRA should begin with the same type of prior specification effort as Bayesian analysis.     

Quantitative Assessment of the Risk of Lung Cancer Associated with Occupational Exposure to Refractory Ceramic Fibers

We present the results of a quantitative assessment of the lung cancer risk associated with occupational exposure to refractory ceramic fibers (RCF). The primary sources of data for our risk assessment were two long-term oncogenicity studies in male Fischer rats conducted to assess the potential pathogenic effects associated with prolonged inhalation of RCF. An interesting feature of the data was the availability of the temporal profile of fiber burden in the lungs of experimental animals. Because of this information, we were able to conduct both exposure–response and dose–response analyses. Our risk assessment was conducted within the framework of a biologically based model for carcinogenesis, the two-stage clonal expansion model, which allows for the explicit incorporation of the concepts of initiation and promotion in the analyses. We found that a model positing that RCF was an initiator had the highest likelihood. We proposed an approach based on biological considerations for the extrapolation of risk to humans. This approach requires estimation of human lung burdens for specific exposure scenarios, which we did by using an extension of a model due to Yu. Our approach acknowledges that the risk associated with exposure to RCF depends on exposure to other lung carcinogens. We present estimates of risk in two populations: (1) a population of nonsmokers and (2) an occupational cohort of steelworkers not exposed to coke oven emissions, a mixed population that includes both smokers and nonsmokers.     

The Promise of Molecular Epidemiology for Quantitative Risk Assessment

In the long run, molecular epidemiological techniques can provide important insights for understanding a wide variety of important issues in current risk assessment and are applicable across a broad spectrum of adverse effects in addition to carcinogenesis. Unfortunately, current risk assessment practices make very little use of the kind of detailed mechanistic information that molecular epidemiology can provide. Eventually, there is reason to hope that the availability of mechanistic insights provided in part by molecular epidemiology can produce some of the "essential tension" required to reform paradigms for the formulation of quantitative risk assessment models in general.     

Distributions of Total Job Tenure for Men and Women in Selected Industries and Occupations in the United States, February 1996

This article develops and fits probability distributions for the variability in projected (total) job tenure for adult men and women in 31 industries and 22 occupations based on data reported by the U.S. Department of Labor's Bureau of Labor Statistics. It extends previously published results and updates those results from January 1987 to February 1996. The model provides probability distributions for the variability in projected (total) job tenures within the time range of the data, and it extrapolates the distributions beyond the time range of the data, i.e., beyond 25 years.