Daily Soil Ingestion Estimates for Children at a Superfund Site

Ingestion of contaminated soil by children may result in significant exposure to toxic substances at contaminated sites. Estimates of such exposure are based on extrapolation of short-term-exposure estimates to longer time periods. This article provides daily estimates of soil ingestion on 64 children between the ages of 1 and 4 residing at a Superfund site; these values are employed to estimate the distribution of 7-day average soil ingestion exposures (mean, 31 mg/day; median, 17 mg/day) at a contaminated site over different time periods. Best linear unbiased predictors of the 95th-percentile of soil ingestion over 7 days, 30 days, 90 days, and 365 days are 133 mg/day, 112 mg/day, 108 mg/day and 106 mg/day, respectively. Variance components estimates (excluding titanium and outliers, based on Tukey's far-out criteria) are given for soil ingestion between subjects (59 mg/day)2, between days on a subject (95 mg/day)2, and for uncertainty on a subject-day (132 mg/day)2. These results expand knowledge of potential exposure to contaminants among young children from soil ingestion at contaminated sites. They also provide basic distributions that serve as a starting point for use in Monte Carlo risk assessments.     

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.     

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.     

A Meta-Analysis of Children's Hand-to-Mouth Frequency Data for Estimating Nondietary Ingestion Exposure

Because of their mouthing behaviors, children have a higher potential for exposure to available chemicals through the nondietary ingestion route; thus, frequency of hand-to-mouth activity is an important variable for exposure assessments. Such data are limited and difficult to collect. Few published studies report such information, and the studies that have been conducted used different data collection approaches (e.g., videography versus real-time observation), data analysis and reporting methods, ages of children, locations, and even definitions of "mouthing." For this article, hand-to-mouth frequency data were gathered from 9 available studies representing 429 subjects and more than 2,000 hours of behavior observation. A meta-analysis was conducted to study differences in hand-to-mouth frequency based on study, age group, gender, and location (indoor vs. outdoor), to fit variability and uncertainty distributions that can be used in probabilistic exposure assessments, and to identify any data gaps. Results of this analysis indicate that age and location are important for hand-to-mouth frequency, but study and gender are not. As age increases, both indoor and outdoor hand-to-mouth frequencies decrease. Hand-to-mouth behavior is significantly greater indoors than outdoors. For both indoor and outdoor hand-to-mouth frequencies, interpersonal, and intra-personal variability are approximately 60% and approximately 30%, respectively. The variance difference among different studies is much bigger than its mean, indicating that different studies with different methodologies have similar central values. Weibull distributions best fit the observed data for the different variables considered and are presented in this article by study, age group, and location. Average indoor hand-to-mouth behavior ranged from 6.7 to 28.0 contacts/hour, with the lowest value corresponding to the 6 to <11 year olds and the highest value corresponding to the 3 to <6 month olds. Average outdoor hand-to-mouth frequency ranged from 2.9 to 14.5 contacts/hour, with the lowest value corresponding to the 6 to <11 year olds and the highest value corresponding to the 6 to <12 month olds. The analysis highlights the need for additional hand-to-mouth data for the <3 months, 3 to <6 months, and 3 to <6 year age groups using standardized collection and analysis because of lack of data or high uncertainty in available data. This is the first publication to report Weibull distributions as the best fitting distribution for hand-to-mouth frequency; using the best fitting exposure factor distribution will help improve estimates of exposure. The analyses also represent a first comprehensive effort to fit hand-to-mouth frequency variability and uncertainty distributions by indoor/outdoor location and by age groups, using the new standard set of age groups recommended by the U.S. Environmental Protection Agency for assessing childhood exposures. Thus, the data presented in this article can be used to update the U.S. EPA's Child-Specific Exposure Factors Handbook and to improve estimates of nondietary ingestion in probabilistic exposure modeling.     

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.     

Field Measurement of Dermal Soil Loading Attributable to Various Activities: Implications for Exposure Assessment

Estimates of soil adherence to skin are required for assessment of dermal exposures to contaminants in soils. Previously available estimates depend heavily on indirect measurements and/or artificial activities and reflect sampling of hands only. Results are presented here from direct measurement of soil loading on skin surfaces of volunteers before and after normal occupational and recreational activities that might reasonably be expected to lead to soil contact. Skin surfaces assayed included hands, forearms, lower legs, faces and/or feet. Observed hand loadings vary over five orders of magnitude (roughly from 10^–3 to 10² mg/cm²) and are dependent upon type of activity. Hand loadings within the current default range of 0.2 to 1.0 mg/cm² were produced by activities providing opportunity for relatively vigorous soil contact (rugby, farming). Loadings less than 0.2 mg/cm² were found on hands following activities presenting less opportunity for direct soil contact (soccer, professional grounds maintenance) and on other body parts under many conditions. The default range does not, however, represent a worst case. Children playing in mud on the shore of a lake generated geometric mean loadings well in excess of 1 mg/cm² on hands, arms, legs, and feet. Post-activity average loadings on hands were typically higher than average loadings on other body parts resulting from the same activity. Hand data from limited activities cannot, however, be used to conservatively predict loadings that might occur on other body surfaces without regard to activity since non-hand loadings attributable to higher contact activities exceeded hand loadings resulting from lower contact activities. Differences between pre- and post-activity loadings also demonstrate that dermal contact with soil is episodic. Typical background (pre-activity) geometric mean loadings appear to be on the order of 10^-2 mg/cm² or less. Because exposures are activity dependent, quantification of dermal exposure to soil will remain inadequate until data describing relevant human behavior (type of activity, frequency, duration including interval before bathing, clothing worn, etc.) are generated.     

Effect of Particle Size on Risk Assessment of Direct Soil Ingestion and Metals Adhered to Children's Hands at Playgrounds

The quantity of heavy metals in soil is measured after 2‐mm sieving in Japan for risk assessment of direct soil ingestion. A study was conducted on the relationship between the size of soil particles and quantity of heavy metals in soil, and the particle‐size distribution of soil adhered to children's hands, and the risks of direct intake of soil considering the particle sizes ingested were evaluated. The results showed that smaller particles had a tendency to contain more heavy metals than bigger ones, that the particle size of approximately 90% of the soil particles from playgrounds adhered to children's hands was less than 100 μm, and that 2‐mm sieving in preparation for measuring heavy metal content caused underestimation of the risk of direct soil intake. The amount of heavy metals on children's hands after playing outside was investigated. Various metals and soil were adhered to their hands, and the amount of soil adhered could be estimated from the concentration of metals. To develop accurate risk assessment, the particle‐size distribution of ingested soil and more detailed scenarios of soil intake are necessary.     

Derivation of a Target Level of Lead in Soil at Residential Sites Corresponding to a De Minimis Contribution to Blood Lead Concentration1

Inability to define either a clear toxicologic threshold or a stochastic all-or-nothing (cancer-type) response model for the noncarcinogenic effects of lead (Pb) in young children has posed difficulties for derivation of risk-based target levels of Pb in residential soil. Approaches based on empirical relationships between Pb levels in blood (PbB) and Pb in soil suffer from inability to specify the numerous variables which mediate between these two quantities. Approaches based on achieving a toxicologically de minimis target PbB level (e.g., 10 μg/dl) are subject to large uncertainty in estimating the distribution of existing PbB levels in a specific exposed population and in estimating the relative contribution from nonsoil sources of Pb. The multisource contribution to the distribution of PbB makes this approach unsuited for determination of a target Pb level in a single medium. An alternative approach is presented based on achieving a de minimis contribution to PbB (ΔPbB) from soil. Contributions to Pb exposure from outdoor soil and indoor soil-derived dust (ISDD) are modeled and appropriate values are suggested for input parameters. This analysis predicts that chronic exposure of young children to 200 μg Pb/g (ppm) in residential soil will result in a ΔPbB of 2 μg Pb/dl blood. This concentration of Pb in soil may provide an appropriate target level for residential soil when other significant sources of Pb exposure are present. In other cases, this approach can be used to predict a soil concentration of Pb corresponding to an appropriate non-de minimisΔPbB.     

Assessment of the Effect of Population and Diary Sampling Methods on Estimation of School‐Age Children Exposure to Fine Particles

Population and diary sampling methods are employed in exposure models to sample simulated individuals and their daily activity on each simulation day. Different sampling methods may lead to variations in estimated human exposure. In this study, two population sampling methods (stratified‐random and random‐random) and three diary sampling methods (random resampling, diversity and autocorrelation, and Markov‐chain cluster [MCC]) are evaluated. Their impacts on estimated children's exposure to ambient fine particulate matter (PM_2.5) are quantified via case studies for children in Wake County, NC for July 2002. The estimated mean daily average exposure is 12.9 μg/m³ for simulated children using the stratified population sampling method, and 12.2 μg/m³ using the random sampling method. These minor differences are caused by the random sampling among ages within census tracts. Among the three diary sampling methods, there are differences in the estimated number of individuals with multiple days of exposures exceeding a benchmark of concern of 25 μg/m³ due to differences in how multiday longitudinal diaries are estimated. The MCC method is relatively more conservative. In case studies evaluated here, the MCC method led to 10% higher estimation of the number of individuals with repeated exposures exceeding the benchmark. The comparisons help to identify and contrast the capabilities of each method and to offer insight regarding implications of method choice. Exposure simulation results are robust to the two population sampling methods evaluated, and are sensitive to the choice of method for simulating longitudinal diaries, particularly when analyzing results for specific microenvironments or for exposures exceeding a benchmark of concern.     

Statistical Distributions of Daily Breathing Rates for Narrow Age Groups of Infants and Children

Children may be more susceptible to toxicity from some environmental chemicals than adults. This susceptibility may occur during narrow age periods (windows), which can last from days to years depending on the toxicant. Breathing rates specific to narrow age periods are useful to assess inhalation dose during suspected windows of susceptibility. Because existing breathing rates used in risk assessment are typically for broad age ranges or are based on data not representative of the population, we derived daily breathing rates for narrow age ranges of children designed to be more representative of the current U.S. children's population. These rates were derived using the metabolic conversion method of Layton (1993) and energy intake data adjusted to represent the U.S. population from a relatively recent dietary survey (CSFII 1994–1996, 1998). We calculated conversion factors more specific to children than those previously used. Both nonnormalized (L/day) and normalized (L/kg-day) breathing rates were derived and found comparable to rates derived using energy estimates that are accurate for the individuals sampled but not representative of the population. Estimates of breathing rate variability within a population can be used with stochastic techniques to characterize the range of risk in the population from inhalation exposures. For each age and age-gender group, we present the mean, standard error of the mean, percentiles (50th, 90th, and 95th), geometric mean, standard deviation, 95th percentile, and best-fit parametric models of the breathing rate distributions. The standard errors characterize uncertainty in the parameter estimate, while the percentiles describe the combined interindividual and intra-individual variability of the sampled population. These breathing rates can be used for risk assessment of subchronic and chronic inhalation exposures of narrow age groups of children.     

Small Parts Aspiration, Ingestion, and Choking in Small Children: Findings of the Small Parts Research Project

Obtaining information on current child injury trends and their associated issues is an important factor in developing products that meet or surpass acceptable toy safety boundaries. Understanding these boundaries helps determine safe product design characteristics that reduce the risk of product-related injury. Inchcape Testing Services developed a Small Parts Aspiration, Ingestion, and Choking Hazards Research Project, independent of an ongoing consumer Product Safety Commission (CPSC) project. The project's purpose was to determine, independent of CPSC, whether a more stringent small parts regulatory standard was necessary and, if so, to ensure that the standard was determined objectively. This article reports on the project's findings relating to critical characteristics (size, shape, and consistency) relative to the victim's age, of objects responsible for child choking injuries and fatalities.     

Exposures from Chrysotile‐Containing Joint Compound: Evaluation of New Model Relating Respirable Dust to Fiber Concentrations

The potential for fiber exposure during historical use of chrysotile‐containing joint compounds (JCC) has been documented, but the published data are of limited use for reconstructing exposures and assessing worker risk. Consequently, fiber concentration distributions for workers sanding JCC were independently derived by applying a recently developed model based on published dust measurements from sanding modern‐day (asbestos‐free) joint compound and compared to fiber concentration distributions based on limited historical measurements. This new procedure relies on factors that account for (i) differences in emission rates between modern‐day and JCC and (ii) the number of fibers (quantified by phase contrast microscopy [PCM]) per mass of dust generated by sanding JCC, as determined in a bench‐scale chamber study using a recreated JCC, that convert respirable dust concentrations to fiber concentrations. Airborne respirable PCM‐fiber concentration medians (and 95% confidence intervals) derived for output variables using the new procedure were 0.26 (0.039, 1.7) f/cm³ and 0.078 (0.013, 0.47) f/cm³, and corresponding total fiber concentrations were 1.2 (0.17, 9.2) f/cm³ and 0.37 (0.056, 2.5) f/cm³, in enclosed and nonenclosed environments, respectively. Corresponding estimates of respirable and total PCM fiber concentrations measured historically during sanding of asbestos‐containing joint compound—adjusted for differences between peak and time‐weighted average (TWA) concentrations and documented analytical preparation and sampling artifacts—were 0.15 (0.019, 0.95) f/cm³ and 0.86 (0.11, 5.4) f/cm³, respectively. The PCM‐fiber concentration distributions estimated using the new procedure bound the distribution estimated from adjusted TWA historical fiber measurements, suggesting reasonable consistency of these estimates taking into account uncertainties addressed in this study.     

Time- and Loading-Dependence in the McKone Model for Dermal Uptake of Organic Chemicals from a Soil Matrix

McKone has recently proposed an innovative two-layer model for dermal uptake of organic chemicals from a soil matrix that explicitly includes variables for properties of the chemical, the soil, the skin, and the exposure. In this note, we investigate the joint time- and loading-dependencies implicit in the model by using MATHEMATICA to find and plot a closed-form function for the uptake fraction for six aromatic hydrocarbons.     

Risk Analysis of Dust Explosion Scenarios Using Bayesian Networks

In this study, a methodology has been proposed for risk analysis of dust explosion scenarios based on Bayesian network. Our methodology also benefits from a bow‐tie diagram to better represent the logical relationships existing among contributing factors and consequences of dust explosions. In this study, the risks of dust explosion scenarios are evaluated, taking into account common cause failures and dependencies among root events and possible consequences. Using a diagnostic analysis, dust particle properties, oxygen concentration, and safety training of staff are identified as the most critical root events leading to dust explosions. The probability adaptation concept is also used for sequential updating and thus learning from past dust explosion accidents, which is of great importance in dynamic risk assessment and management. We also apply the proposed methodology to a case study to model dust explosion scenarios, to estimate the envisaged risks, and to identify the vulnerable parts of the system that need additional safety measures.     

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.     

Assessment of Health Risk from Exposure to Contaminated Soil

The risk to human health posed by contaminated soil in a residential area depends on the potential extent of exposure to soil and on the toxic properties of the contaminants. A detailed soil exposure analysis is presented for young children, older children, and adults living in a house surrounded by contaminated soil. From this analysis, a lifetime exposure model is derived and used to assess chronic health risks.     

Dermal Uptake of Organic Chemicals from a Soil Matrix

Uptake of chemicals from soil on human skin is considered. Based on a review of literature on the structure of human skin, the processes by which chemicals pass through this boundary, and experiments that reveal the rate and magnitude of this transport process; a two-layer model is presented for estimating how chemical uptake through the stratum corneum depends on chemical properties, skin properties, soil properties and exposure conditions. The model is applied to two limiting scenarios--(1) continuous deposition and removal of soil on the skin surface and (2) a one-time deposition of soil onto the skin surface. The fraction of soil-bound chemical that passes through the stratum corneum is dependent on the skin-soil layer thickness; the dimensionless Henry's law constant, Kh and the octanol-water partition coefficient, Kow of the soil-bound chemical. The nature of this dependence is discussed.     

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.     

A Probabilistic Arsenic Exposure Assessment for Children Who Contact Chromated Copper Arsenate (CCA)-Treated Playsets and Decks, Part 2: Sensitivity and Uncertainty Analyses

A probabilistic model (SHEDS-Wood) was developed to examine children's exposure and dose to chromated copper arsenate (CCA)-treated wood, as described in Part 1 of this two-part article. This Part 2 article discusses sensitivity and uncertainty analyses conducted to assess the key model inputs and areas of needed research for children's exposure to CCA-treated playsets and decks. The following types of analyses were conducted: (1) sensitivity analyses using a percentile scaling approach and multiple stepwise regression; and (2) uncertainty analyses using the bootstrap and two-stage Monte Carlo techniques. The five most important variables, based on both sensitivity and uncertainty analyses, were: wood surface residue-to-skin transfer efficiency; wood surface residue levels; fraction of hand surface area mouthed per mouthing event; average fraction of nonresidential outdoor time a child plays on/around CCA-treated public playsets; and frequency of hand washing. In general, there was a factor of 8 for the 5th and 95th percentiles and a factor of 4 for the 50th percentile in the uncertainty of predicted population dose estimates due to parameter uncertainty. Data were available for most of the key model inputs identified with sensitivity and uncertainty analyses; however, there were few or no data for some key inputs. To evaluate and improve the accuracy of model results, future measurement studies should obtain longitudinal time-activity diary information on children, spatial and temporal measurements of residue and soil concentrations on or near CCA-treated playsets and decks, and key exposure factors. Future studies should also address other sources of uncertainty in addition to parameter uncertainty, such as scenario and model uncertainty.