Intake Fraction for the Indoor Environment: A Tool for Prioritizing

Aug 25, 2012 - A framework for an alternatives assessment dashboard for evaluating chemical alternatives applied to flame retardants for electronic ap...
0 downloads 15 Views 5MB Size
Article pubs.acs.org/est

Intake Fraction for the Indoor Environment: A Tool for Prioritizing Indoor Chemical Sources Hyeong-Moo Shin,*,† Thomas E. McKone,‡ and Deborah H. Bennett† †

Department of Public Health Sciences, University of CaliforniaDavis, Davis, California, United States School of Public Health, University of CaliforniaBerkeley, Berkeley, California, United States



S Supporting Information *

ABSTRACT: Reliable exposure-based chemical characterization tools are needed to evaluate and prioritize in a rapid and efficient manner the more than tens of thousands of chemicals in current use. This study applies intake fraction (iF), the integrated incremental intake of a chemical per unit of emission, for a suite of indoor released compounds. A fugacity-based indoor mass-balance model was used to simulate the fate and transport of chemicals for three release scenarios: direct emissions to room air and surface applications to carpet and vinyl. Exposure through inhalation, dermal uptake, and nondietary ingestion was estimated. To compute iF, cumulative intake was summed from all exposure pathways for 20 years based on a scenario with two adults and a 1year-old child who ages through the simulation. Overall iFs vary by application modes: air release (3.1 × 10−3 to 6.3 × 10−3), carpet application (3.8 × 10−5 to 6.2 × 10−3), and vinyl application (9.0 × 10−5 to 1.8 × 10−2). These iF values serve as initial estimates that offer important insights on variations among chemicals and the potential relative contribution of each pathway over a suite of compounds. The approach from this study is intended for exposure-based prioritization of chemicals released inside homes. pathways are additional routes of exposure for young children.8 These factors all strengthen the need for quantifying exposures for compounds released to the indoor environment. Intake fraction (iF) is defined as the integrated incremental intake of a pollutant per unit of emission,9,10 and we consider it a useful concept for expressing source-to-intake relationships. Because iF is unitless, it is considered a convenient standardized metric of exposure. A number of researchers have applied iF to particulate matter (PM) and other types of emissions from power plants and mobile sources to determine the relationship among emissions, human intake, and health impacts.11,12 Moreover, iF has been calculated for over 300 organic chemicals released into ambient air and water to examine the relationship between chemical properties and primary exposure routes outdoors.9 Three studies estimated iF through inhalation of PM from indoor sources11,13 and pollutants from episodic indoor emissions.14 However, no other studies to date have included dermal and nondietary exposure pathways to quantify iF values over a wide range for organic compounds found in the indoor environment. Several studies report that exposure through incidental dust ingestion and dermal contact is important in predicting polybrominated diphenyl ether (PBDE) and pesticide body burden,15−19 indicating the importance of

1. INTRODUCTION Over the last several decades, there have been significant increases in the volume of chemicals used in the manufacture of residential consumer products, including plastics, clothing, electronics, and cosmetics.1 Not only the volume of chemicals produced but also the number and type of chemicals have increased. However, information about potential exposure and adverse health effects in humans from residential uses is limited for most chemicals. Thus, the development of both toxicityand exposure-based assessments is needed to evaluate and prioritize chemicals in a rapid and efficient manner.1,2 Approaches for systematically prioritizing chemicals require human exposure data, product use information, and modeled human behavior and must address both near-field (from consumer products used indoors) and far-field (from overall environmental dispersion) exposures.1,3 However, most current prioritization approaches focus mainly on far-field exposure and their output impact metrics are not standardized.1 Because individuals generally spend over 70% of their time indoors at home,4 both the residence times of chemicals in the indoor environment and residential human activity patterns are key inputs to chemical prioritization. Several studies report that the concentrations of volatile organic compounds (VOCs) are consistently higher indoors than outdoors.5,6 Moreover, the levels of pesticides found most frequently in carpet dust tend to be higher than those in outdoor soil, and many pesticides have been found to be much more persistent in the indoor environment relative to outdoors.7 For less volatile and more persistent compounds, dermal uptake or nondietary ingestion © 2012 American Chemical Society

Received: Revised: Accepted: Published: 10063

May 7, August August August

2012 17, 2012 24, 2012 25, 2012

dx.doi.org/10.1021/es3018286 | Environ. Sci. Technol. 2012, 46, 10063−10072

Environmental Science & Technology

Article

Table 1. Physicochemical Properties of 15 Studied Chemicals MW (g/mol) formaldehyde naphthalene diethyl phthalate galaxolide tris(2-chloroethyl) phosphate diazinon chlorpyrifos PBDE-47 bis(2-ethylhexyl) phthalate triphenyl phosphate butyl benzyl phthalate PBDE-99 permethrin benzo[a]pyrene PBDE-209

30.0 128.2 222.2 258.2 284.0 304.0 350.6 485.8 390.5 326.1 312.0 564.7 391.3 252.3 959.2

Kowb 2.2 2.3 5.1 1.1 3.0 5.9 1.0 4.8 3.3 3.9 8.1 2.0 4.4 1.5 2.8

× × × × × × × × × × × × × ×

VPc(Pa) 3

10 102 105 101 103 105 106 108 104 104 107 107 106 109

4.6 2.1 2.2 5.5 1.4 9.4 2.2 5.6 5.3 1.7 9.5 3.3 9.9 2.5 2.2

× × × × × × × × × × × × × × ×

5

10 101 10−1 10−2 10−2 10−3 10−3 10−4 10−4 10−4 10−5 10−5 10−7 10−7 10−10

Hd (Pa·m3/mol) 7.0 1.9 9.3 5.9 5.5 1.3 4.7 1.1 1.8 7.5 3.9 3.0 1.0 3.4 1.5

× × × × × × ×

10 101 10−2 10−1 10−4 10−1 10−1

× × × × × ×

10−3 10−3 10−1 10−2 10−2 10−3

1

Daire (m2/day) 1.4 5.2 2.3 3.0 2.1 1.8 3.3 2.7 1.5 5.0 1.5 1.9 1.7 4.1 1.7

× × × × × × × × × × × × × ×

−1

10 10−1 10−1 10−1 10−1 10−1 10−1 10−1 10−1 10−1 10−1 10−1 10−1 10−1

kOHf cm3/mol·day) 8.1 1.9 3.0 1.9 2.3 1.9 7.9 2.3 1.9 2.3 9.5 2.3 2.8 4.3 2.3

× × × × × × × × × × × × × × ×

10−7 10−6 10−7 10−6 10−6 10−6 10−6 10−6 10−6 10−6 10−7 10−6 10−6 10−6 10−6

a Molecular mass. bOctanol−water partition coefficient (unitless). cVapor pressure. dHenry’s law constant. eDiffusion coefficient in pure air. fOH radical reaction rate constant.

be used to predict time-varying exposures if the emission profile in the indoor residential environment is known. Ideally, a detailed quantitative uncertainty assessment should be applied to the model results to account for uncertainty in chemical property values, variability and uncertainty in housing characteristics, and model uncertainty resulting from various assumptions related to dermal and nondietary exposure pathways. Insufficient information exists to reliably quantify the uncertainty of parameter input values, and there are vast differences between various models available for dermal and nondietary exposure. Therefore, we focused on relatively large modeling uncertainty by determining the sensitivity of results to the selection of two alternative dermal and nondietary models. 2. 2. Indoor Fate and Transport Model. 2.2.1. Model Description and Input Parameters. We used a fugacity-based mass-balance model to simulate the fate and transport of compounds in the indoor environment21 (see Supporting Information for summary of the model). The model includes four compartments that capture the major indoor reservoirs for chemicals: air, carpet, vinyl flooring, and walls. Each compartment in the model is composed of multiple phases, such as gases and particles in the air compartment. We used multiple size fractions of particles, including 0−2.5, 2.5−10, and 10−150 μm, because of the unique physical properties of each particle size fraction. We assumed that each phase in a given compartment was in chemical equilibrium with all other phases in that compartment. However, the compartments are not required to be in chemical equilibrium with each other because there are driving diffusive mass transfers between compartments. Specifically, a compound can partition to carpet after being released to air and then slowly volatilize back out into the air. Since semivolatile organic compounds (SVOCs) tend to partition to dust, advective transport pathways of dust, including resuspension and deposition, were considered. A compound can be removed from the indoor system through ventilation of gas and particle phases or through surface cleaning. A set of differential equations accounts for gains and losses in each compartment as well as transfers between compartments. Input parameters and the properties of particles in different size fractions are listed in Tables S1 and S2 of the Supporting Information. In applying the model, we first defined

these pathways. Quantifying exposure through dermal and nondietary pathways requires a number of assumptions regarding both the primary mechanisms for exposure and parameter values associated with quantifying intake from these mechanisms. The objectives of this study are to (1) estimate iF for a suite of organic compounds from residential exposure via three release scenarios (air emission and carpet and vinyl flooring applications), (2) examine the relationship between intake through different exposure pathways and chemical properties by varying the magnitude of chemical volatility, and (3) address the data and model gaps and their impact for estimating dermal and nondietary ingestion exposure.

2. METHODS 2. 1. Overview. Calculating iF for chemicals released to an indoor residence requires a number of steps. First, the release scenario is defined. Second, the resulting fate and transport of the chemical through multiple indoor compartments is determined by use of an indoor fugacity model. Third, exposure concentrations are calculated from the mass associated with any compartment. Fourth, exposure over time is calculated via multiple exposure models. Finally, the integrated incremental intake of a pollutant per unit of emission, iF, is calculated. Calculating iF, rather than predicting exposure from a broad range of specific source input scenarios, allows simplified release scenarios because the magnitude of iF is not strongly dependent on the details of the release scenario beyond the release compartment. This is critical if the model is to be used to screen chemicals. Compounds were compared on the basis of a pulse−release (i.e., one-time release at time = 0) into three primary household media: air, carpet, and vinyl. A pulse− release may not capture all of the details of different modes of application for the chemicals modeled, but it allows chemical prioritization when calculating iF. For example, if a constant number of individuals is assumed, for a continuous airborne release of a volatile compound, such as naphthalene from a mothball, iF is the same as for a pulse input.20 When chemicals are screened in a prioritization framework, it is unlikely that the emission scenario would be known, particularly if the chemical was assessed prior to its use in commerce. This model can also 10064

dx.doi.org/10.1021/es3018286 | Environ. Sci. Technol. 2012, 46, 10063−10072

Environmental Science & Technology

Article

Table 2. Age-Specific Exposure Parameter Values for Inhalation and Nondietary Ingestion age (years)

inhalation ratea (m3/day)

hand-to-mouth frequencya (events/day)

object-to-mouth frequencya (events/day)

surface area of handsb (m2)

fraction of time spent at residenceb

fraction of time sleeping and nappingb

0.5 1 2 3 4 5 6 7 8 9 10

5.4 8.0 8.9 10.1 10.1 10.1 12.0 12.0 12.0 12.0 12.0

190 200 130 150 150 150 70 70 70 70 70

240 200 100 100 100 100 10 10 10 10 10

0.024 0.030 0.028 0.037 0.037 0.037 0.051 0.051 0.051 0.051 0.051

0.77 0.74 0.68 0.66 0.66 0.66 0.62 0.62 0.62 0.62 0.62

0.53 0.54 0.49 0.47 0.47 0.47 0.43 0.43 0.43 0.43 0.43

a

Child-Specific Exposure Factors Handbook (EPA 2008).23 bExposure Factors Handbook (EPA 2011).26

grams per cubic meter), Rinh is the inhalation uptake rate (cubic meters per day), and f res is the fraction of time spent at the residence (unitless). For inhalation of air particles, we included chemical mass only for particles less than 10 μm in diameter.24 Human activities, such as walking, folding blankets, and vacuuming, can resuspend dust from surface compartments. Ferro et al.25 measured the mean PM2.5 and PM5 concentrations for the three activities mentioned above and reported that personal exposures to PM2.5 and PM5 were respectively 1.4 and 1.6 times as high as direct exposure to average indoor concentrations. Thus, we accounted for both resuspended particles from human activities and the fraction of human activity time at residence in computing exposure through particle inhalation. We estimated the fraction of human activity time at residence by subtracting the fraction of time spent sleeping/napping (fslp) from the fraction of time spent at residence ( f res). The age-specific values for fslp and f res were obtained from the U.S. Environmental Protection Agency (EPA) Exposure Factors Handbook (EFH)26 and are listed in Table 2. The average inhalation uptake rates for male and female adults (Table 3) were applied from the EFH. Those for children under 12 years of age were applied from the ChildSpecific Exposure Factors Handbook (CSEFH).23 The

the fugacity capacity for each compartment, including all phases (Table S3, Supporting Information) and transfer factors between compartments (Table S4, Supporting Information). Then, we solved ordinary differential equations to compute the time-dependent mass in each compartment from the mass balance equations (Table S5, Supporting Information) using MATLAB (Mathworks, Natick, MA). 2.2.2. Chemical Release Scenarios. We selected to study 15 chemicals based on concern for their effects on human health and their large range in vapor pressure values (see Table 1 for chemical properties). The selected chemicals were the following: formaldehyde, naphthalene, galaxolide (a synthetic musk), benzo[a]pyrene, three pesticides commonly found in the indoor environment (diazinon, chlorpyrifos, and permethrin), three common polybrominated diphenyl ether (PBDE) congeners as flame retardants (PBDE-47, PBDE-99, and PBDE-209), two flame retardants newly detected in house dust [tris(2-chloroethyl) phosphate and triphenyl phosphate], and three plasticizers [diethyl phthalate, bis(2-ethylhexyl) phthalate, and butyl benzyl phthalate]. We assumed that each chemical was emitted or applied as a pulse−release to the air, carpet, and vinyl compartments, and we computed timedependent concentration profiles for 20 years. 2.3. Exposure Model. We chose a representative exposure scenario, consisting of two adults and a 1-year-old child who ages through the simulation, based on U.S. Census Bureau data published in 2011, in which the average household was reported to consist of 1.95 adults and 0.63 children.22 In our exposure scenario, the 1-year-old child is tracked for 20 years, during which age-specific exposure factors, including inhalation uptake rates, hand-to-mouth frequencies, and surface area of hands, change over time. A conservative assumption was made that the occupants were present indoors when the chemicals were released and that both adults and the child were exposed through inhalation and dermal uptake. We estimated exposure through nondietary ingestion only for the young child, as the frequency of mouthing activities after surface contact is considerably lower for adults, although we note that this may underestimate adult exposure.23 2.3.1. Inhalation Exposure. For each individual, daily intake (micrograms per day) via inhalation is based on the following equation: Iinh = (Cg + Cp)R inhfres

Table 3. Exposure Parameter Values parameter inhalation uptake rate (m3/day) fraction of time spent at residence fraction of time sleeping/napping transfer efficiency of carpet-to-hand transfer efficiency of vinyl-to-hand transfer efficiency of hand-to-mouth transfer efficiency of object-to-mouth surface area of hands contacted with surface or object (m2/event) surface area of object or hand that is mouthed (m2/event) exposure frequency of hand-to-carpet (events/h) exposure frequency of hand-to-vinyl (events/h) exposure frequency of hand-to-object (events/h) fraction absorbed through skin fraction absorbed through gut

(1)

where Cg and Cp are the exposure concentrations in the gas and particle phases, respectively, of the air compartment (micro10065

symbol

child

adult

ref

Rih Fresid Fsleep TEc TEv TEh_m TEo_m SAi

0.034 0.033 0.5 0.5 0.098

26 26 26 33 40 45 45 26

SAx

0.005

45

EFc

11

EFv

12

EFo

115

ABs ABg

0.02 0.70

34, 36, 38 34, 36, 38 35, 36, 37 44 41

15.5 0.65 0.36

dx.doi.org/10.1021/es3018286 | Environ. Sci. Technol. 2012, 46, 10063−10072

Environmental Science & Technology

Article

floor, and objects from five published studies34−38 to obtain the values 11, 12, and 115 events/day, respectively, for use in our model (see Supporting Information, Table S7). Although the transfer efficiency from surfaces to human skin is chemicalspecific and influenced by other factors, including condition of skin on the palm, we used the average transfer efficiency of chlorpyrifos from carpet to hand (TEc = 3.35%) from Lu and Fenske33 and Cohen Hubal et al.39 as our default value for carpet and the average of all analytes from vinyl flooring to hand (TEv = 3.15%) from Clothier40 and Cohen Hubal et al.39 as our default value for vinyl. For the absorption fraction through skin (ABSs), we assumed that 2% of the chemical in soil or dust was absorbed through skin.41 In addition to carpet and vinyl flooring, children contact numerous other objects such as toys and furniture, which can contain the chemicals selected for this study. There has been little work to determine the relationship between the concentration on carpet/vinyl surface compartments and the surface concentration on other objects. The existing work is difficult to interpret because toy concentrations were based on an extraction technique rather than a wipe sample.42 In this study, we assumed that the concentrations on the surface of other objects and the transfer efficiency from other objects to hands can be represented by area-weighted values of carpet and vinyl surfaces, acknowledging that this is an uncertain parameter. 2.3.3. Nondietary Ingestion Exposure. As with dermal exposure, there are multiple methods for calculating nondietary ingestion. Cohen-Hubal et al.27 integrated nondietary uptake based on a series of discrete transfers from hand- and object-tomouth activities. Ozkaynak et al.43 estimated 68 mg/day soil and dust ingestion rate via hand and object mouthing, using the residential-SHEDS model. This ingestion rate is similar to the EPA’s recommended soil and dust ingestion rate (= 60 mg/ day) for children.23 Thus, chemical intake through soil and dust ingestion can be estimated by multiplying soil/dust ingestion rate (milligrams per day) by chemical concentrations in soil/ dust (micrograms of chemical/milligram of dust and soil).44 We calculated nondietary exposure using both the Cohen Hubal and Ozkaynak models for comparison, with the Cohen Hubal model as the baseline approach for nondietary exposure. We considered two nondietary ingestion pathways: object-tomouth and hand-to-mouth. In both cases, the equation to estimate the daily intake (micrograms per day) from nondietary ingestion is based on the approach from the model of Cohen Hubal et al.:27

resulting equations for exposure concentrations for all compartments are listed in Table S6 of the Supporting Information. 2.3.2. Dermal Exposure. There are multiple methods available to estimate intake through dermal exposure. Cohen Hubal et al.27 integrated dermal uptake based on a series of discrete transfers from each contact with a contaminated medium. Weschler and Nazaroff28 estimated air-to-skin transdermal uptake of SVOCs by using two permeability coefficients: (1) skin lipid and (2) indoor air transdermal. Zartarian et al.8 developed the residential stochastic human exposure and dose simulation model for pesticides (residentialSHEDS) to estimate dermal exposure and dose by use of microlevel time−location−activity information. We selected to use both the Cohen Hubal and Weschler and Nazaroff models for calculating dermal exposure and comparing the results as a sensitivity analysis. The Cohen Hubal model served as the baseline approach. Because the stochastic approach required to use SHEDS makes it difficult to evaluate pathway contributions over long periods of time, and our goal is to provide a screening-level assessment model, we elected not to include the SHEDS model in our final comparison among alternative modeling approaches. The details for the Cohen Hubal model are presented below, while those for the Weschler and Nazaroff model are presented in the Supporting Information. We used the following equation to estimate the intake (micrograms per day) from dermal exposure using the “microactivity” approach from the model of Cohen Hubal et al.:27 Iderm = Ci(SA)(EF)(TE i i)(ABSs )

(2)

where Ci is the exposure concentration on surface i (carpet or vinyl) (mirograms per square meter), SA is the surface area of the hand contacted (square meters per event), EFi is the exposure frequency of contacts per day for surface i (events/ day), TEi is the efficiency of the contaminant transferring from surface i to the hand (unitless), and ABSs is the fraction of mass absorbed through skin (unitless). For the exposure concentration associated with the carpet compartment (Cc), we assumed that only dust was transferred from carpet and that only a portion of the contaminant residue in that dust was available for transfer to the skin. Dislodgeable fractions of compounds from carpet are often measured by using a polyurethane foam (PUF) roller.7,29 Thus, we multiplied chemical concentrations on carpet dust by the average value of dislodgeable fraction (1.5%) as measured by a PUF roller from Nishioka et al.30 and Fortune et al.31 For the exposure concentration associated with the vinyl compartment (Cv), we assumed that the dust portion of the vinyl compartment and the organic film, which builds up on impermeable surfaces,32 were available for exposure and that 100% of compounds in the organic film and vinyl dust were dislodgeable. This assumption is consistent with measurement techniques, such as surface wipes, used to estimate surface concentrations.33 Based on chemical properties listed in Table 1 and the fugacity approach, the percent of compound mass partitioned to the organic film and vinyl dust relative to vinyl itself ranges from 0.05% to 70%.21 We assumed that mass associated with the walls was not relevant for exposure assessments. We applied the age-specific surface area of hands (SA) from the CSEFH for children under 12 years of age.23 We estimated the average frequency of contacts over one day for carpet, vinyl

Indi = Cx(SA x)(EFx )(TEx)(ABSg )

(3)

where x is the hand or object that is mouthed, Cx is the exposure concentration on hand or object placed in the mouth (micrograms per square meter), SAx is the surface area of hand or object that is mouthed (square meters per event), EFx is the exposure frequency of the mouthing event for the hand or object (events/day), TEx is the transfer efficiency from hand or object to mouth (unitless), and ABSg is the fraction of mass absorbed through gut (unitless). For the exposure concentration on hand (Ch), we multiplied area-weighted exposure concentrations on carpet (Cc) or vinyl (Cv) by transfer efficiencies from carpet to hand (TEc) or vinyl to hand (TEv). For the exposure concentration on objects (Co), we used area-weighted exposure concentrations on carpet or vinyl because objects are inserted directly in the mouth. We assumed that the surface area of an object that is mouthed was 10066

dx.doi.org/10.1021/es3018286 | Environ. Sci. Technol. 2012, 46, 10063−10072

Environmental Science & Technology

Article

Figure 1. Time-dependent relative contribution of each exposure pathway to the total intake and cumulative exposure profiles (solid lines) for chlorpyrifos (left) and PBDE-209 (right).

the same as the surface area of a finger (i.e., 10% of hand) that is mouthed.45 We assumed that the transfer efficiency from hand or object to mouth (TEx) was 0.5, acknowledging that this is an uncertain parameter.45 For the absorption fraction through gut (ABSg), we assumed that 70% of an administered dose was absorbed through the gut.44 For children under 12 years of age, we obtained age-specific exposure parameters, including hand-to-mouth frequency, object-to-mouth frequency, and the surface area of hand (Table 3) from the CSEFH.23

3. RESULTS AND DISCUSSION 3.1. Exposure Profile. We estimated intake for 15 organic compounds commonly found in indoor environmental samples using the fugacity model time-dependent mass distribution among indoor compartments, human activity patterns, our inhalation model, and microactivity approaches of Cohen Hubal et al.27 for dermal and nondietary ingestion exposure. The chemicals were grouped into three categories on the basis of their relative chemical properties. Formaldehyde, naphthalene, and diethyl phthalate were categorized as “volatile” compounds. Galaxolide, tris(2-chloroethyl) phthalate, diazinon, 10067

dx.doi.org/10.1021/es3018286 | Environ. Sci. Technol. 2012, 46, 10063−10072

Environmental Science & Technology

Article

Figure 2. Contribution of intake from each exposure pathway to the total intake (left) and absolute intake fractions (right) for 15 studied chemicals from different chemical application modes. Red triangles and blue diamonds indicate minimum and maximum intake fractions, respectively, depending on the selection of models for dermal and nondietary ingestion.

and chlorpyrifos were categorized as “minimally volatile” compounds. The remaining compounds, benzo[a]pyrene, permethrin, PBDE-47, PBDE-99, PBDE-209, triphenyl phosphate, bis(2-ethylhexyl) phthalate, and butyl benzyl phthalate, were categorized as “nonvolatile” compounds. The primary exposure pathway for volatile compounds is always gas-phase inhalation regardless of application type. Minimally volatile compound exposure is dependent on the release scenario, with gas-phase inhalation dominating for releases to air and nondietary ingestion dominating for applications to surfaces.

For nonvolatile compounds, inhalation of particles is important for releases to air and applications to carpet, while dermal and nondietary ingestion are significant for releases to air and applications to carpet and are the dominant exposure pathways for applications to vinyl. Figure 1 illustrates the relative contribution of each exposure pathway to the total intake and cumulative intake for chlorpyrifos (a, b, and c) and PBDE-209 (e, f, and g) as representative minimally volatile and nonvolatile compounds, respectively. Gas-phase inhalation is a dominant exposure 10068

dx.doi.org/10.1021/es3018286 | Environ. Sci. Technol. 2012, 46, 10063−10072

Environmental Science & Technology

Article

Figure 3. Sensitivity analysis on estimates of intake fractions from different approaches for dermal and nondietary ingestion for diethyl phthalate, chlorpyrifos, permethrin, and benzo[a]pyrene after release to indoor air. “Dermal-Permeation” represents dermal exposure from the permeation model of Weschler and Nazaroff,28 and “Nondietary-Dust Ingest” represents nondietary ingestion from the dust/soil ingestion rate by Ozkaynak et al.43

persistence. In contrast, when these same two compounds are applied to vinyl, it takes about 90 days for them to be removed. 3.2. Overall Intake Fractions for Each Chemical. Figure 2 shows the relative contribution of intake from each exposure pathway to the total intake (left) and the absolute iF (right) for the three chemical application modes. In the case of air release, overall baseline iF ranges from 3.1 × 10−3 to 6.3 × 10−3. For volatile and minimally volatile compounds, the total intake is dominated by inhalation of the gas phase. For nonvolatile compounds, the contribution of particle-phase inhalation, dermal uptake, and nondietary ingestion becomes significant to the total intake. In the case of application to carpet, iF is estimated from 3.8 × 10−5 to 6.2 × 10−3. The overall pattern of contribution of each exposure pathway to the total intake is similar to the results for air release. In the case of application to vinyl, iF is estimated from 9.0 × 10−5 to 1.8 × 10−2. The pattern of exposure pathway contribution for vinyl is clearly different from air and carpet chemical application modes. This arose because of higher dislodgeable fractions of vinyl dust and organic film compared to the carpet compartment and the higher transfer rates from vinyl to air. Because there are no data available to determine the dislodgeable fraction of compounds from vinyl dust and organic carbon to hands, we used an approach that tends toward higher rather than lower exposure estimates. In evaluating the results in Figure 2, we observed that exposure to volatile compounds was primarily driven by gasphase inhalation due to high transfer rates from vinyl to air, while exposure to nonvolatile compounds was driven by dermal and nondietary ingestion due to high exposure concentrations

pathway for chlorpyrifos released into air, whereas particlephase inhalation accounts primarily for PBDE-209 exposure. When applied to carpet, more than 95% of chlorpyrifos exposure is driven by gas-phase inhalation because it moves from carpet to air compartment due to high diffusive and advective transfer rates from carpet to air, while particle-phase inhalation, dermal uptake, and nondietary ingestion account for about 50%, 10%, and 40% of PBDE-209 exposure, respectively. When applied to vinyl, about 50% of chlorpyrifos is taken up from gas-phase inhalation and the rest of the exposure is from dermal and nondietary ingestion. In contrast, dermal uptake and nondietary ingestion account for about 30% and 70% of PBDE-209 exposure, respectively. For this application mode, the contribution of dermal exposure increases for both compounds because the dislodgeable fraction (100%) of the compound associated with the vinyl compartment is much greater than the dislodgeable fraction (1.5%) of the compound associated with the carpet compartment, resulting in higher exposure concentrations on vinyl surfaces. The exposure profile for a compound depends on the actual application compartment. For example, if chlorpyrifos was applied to a vinyl surface in a kitchen or bathroom, a common application room,46 dermal and nondietary exposure would be important, as has been predicted in detailed exposure studies.43,44 In considering the overall persistence of chemicals in the indoor environment, we noted that with respect to the three application modes, both chlorpyrifos and PBDE-209 remained the longest indoors when applied to carpet. This occurs because a significant portion of the compound is in surface reservoirs, such as carpet fibers or foam materials, resulting in long indoor 10069

dx.doi.org/10.1021/es3018286 | Environ. Sci. Technol. 2012, 46, 10063−10072

Environmental Science & Technology

Article

on vinyl surfaces, associated primarily with Kow values. However, tris(2-chloroethyl) phosphate, diazinon, and chlorpyrifos have lower transfer rates from vinyl to air, resulting in relatively lower exposure concentrations in air. These chemicals also have lower Kow values, resulting in lower exposure concentrations on vinyl surfaces, implying that chemical properties determine both total intake and primary exposure pathways. Ventilation, surface cleaning, and degradation due to OH radical reactions are considered potentially important removal pathways for the selected compounds. When released to air, ventilation is a primary removal pathway for all compounds due to the large air exchange rate (Figure S2a, Supporting Information). As the vapor pressure decreases, the contribution of ventilation to the total removal rate decreases while that of surface cleaning increases. When a compound is applied to carpet, ventilation is still a primary removal pathway (>85%) for volatile and minimally volatile compounds (Figure S2b, Supporting Information). However, surface cleaning accounts for more than 70% of total removal for nonvolatile compounds. When a compound is applied to vinyl, ventilation is a dominant removal pathway for formaldehyde and naphthalene, whose vapor pressure values are greater than 10 Pa. On the other hand, diethyl phthalate is removed almost equally by ventilation and surface cleaning, and other chemicals are dominantly removed by surface cleaning (Figure S2c, Supporting Information). 3.3. Sensitivity Analysis: Dermal and Nondietary Ingestion. We calculated iF from the various dermal and nondietary models, which are fundamentally different in underlying assumptions and processes, and found large differences, indicating that model choice is likely one of the greatest contributors to output uncertainties. In order to explore the dependence of iF on chemical properties and alternative modeling approaches, we carried out a sensitivity analysis for four of the compounds (diethyl phthalate, chlorpyrifos, permethrin, and benzo[a]pyrene). Figure 3 illustrates the results of the sensitivity analysis for air release, providing estimates of iF via dermal uptake and nondietary ingestion from different methods. For diethyl phthalate, which has a relatively large inhalation intake, model selection for dermal and nondietary ingestion pathways was not significant on the total iF. For chlorpyrifos and benzo[a]pyrene, although changing a model for dermal uptake (from a touch model to a permeation model) and nondietary ingestion (from a touch model to a dust ingestion model) resulted in increased iF by approximately 0.5−2 orders of magnitude, the total iF did not increased significantly. However, the importance of dermal uptake and nondietary ingestion became comparable to inhalation intake for these compounds. For permethrin, the total iF increased approximately 1.5 orders of magnitude upon changing from a touch model to a dust ingestion model for nondietary ingestion, making nondietary ingestion the driving exposure pathway. Results from carpet and vinyl application are provided in Figures S3 and S4 of the Supporting Information. Using all available models, the minimum (red triangle) and maximum (blue diamond) iF values are shown in Figure 2, indicating that, depending on the chemical application mode, nonvolatile compounds have uncertainties of about 0.5−2 orders of magnitude in the estimates of intake as a result of the uncertainties in dermal and nondietary ingestion. A complete breakdown of iF for both adults and children from all pathways and models, including the percent dermal uptake directly from

air versus uptake from surface contact in the Weschler and Nazaroff model, can be found in Table S8 of the Supporting Information. Use of the alternate model generally increased total iF for emissions to air and applications to carpet. With the baseline models, iF was found to be slightly lower for the lessvolatile compounds than the more volatile compounds, as lessvolatile compounds are more rapidly removed from the air, decreasing inhalation exposure, while dermal and nondietary exposure are not increased accordingly. With the alternate models, increases in dermal and nondietary exposures for the less-volatile compounds far exceed the decrease in exposure through inhalation resulting from partitioning out of air. This results in the less volatile compounds having a higher iF than the more volatile compounds. Thus, the choice of model for dermal and nondietary intake results in different absolute values of total iF. In addition, the model choice dramatically alters the importance of dermal and nondietary exposure, altering the primary exposure pathway for compounds such as chlorpyrifos, permethrin, and benzo[a]pyrene. As a result, the relative ranking of iF shifts. This implies that exposure models should be carefully selected in screening large sets of compounds. Thus, further research is needed to ensure this pathway is properly quantified. 3.4. Implications/Conclusions. It is well established that iF is greater for releases to the indoor environment than for releases to the outdoor environment.9,14 This work is an effort to build more capacity for addressing indoor exposures. However, there has been a lack of understanding about how iF values for indoor releases vary among a broad range of chemical compounds and the relative importance of various exposure pathways for indoor releases. This paper addresses this critical need by providing a first attempt to estimate multipathway residential indoor iF values for a suite of organic compounds. The reported indoor iF values range from 10−3 to 10−1,47 highlighting the critical need to understand releases to the indoor environment relative to the outdoor environment, as iF values for releases to the outdoor environment range from 10−8 to 10−4.47 The 2 orders of magnitude range in indoor iF values also points to the importance of understanding variability between compounds. Estimating dermal and nondietary ingestion via direct observation is a necessary adjunct to this work. However, direct observation is costly and labor-intensive and has many challenges (e.g., different residue-to-skin transfer efficiency with repeated contacts). Thus, it is vital to provide a modeling framework to help guide and interpret these observational studies. In addition, the difficulties or technical issues of collecting urine and blood samples from a young child make it challenging to validate model results associated with exposures driven by dermal and nondietary ingestion. Although the estimates of total intake from these pathways are uncertain, the positive correlations between PBDE biomarkers and indoor dust concentrations suggest that dermal and nondietary ingestion are potentially important pathways contributing to human body burden for PBDEs.15,19,48,49 Moreover, although intake estimates from dermal and nondietary pathways have significant uncertainties, excluding these pathways fails to provide a comprehensive estimation of exposure. For example, Xu et al.50 reported that the total exposure to diethyl phthalate emitted from vinyl flooring was primarily from the oral ingestion of dust. 10070

dx.doi.org/10.1021/es3018286 | Environ. Sci. Technol. 2012, 46, 10063−10072

Environmental Science & Technology

Article

(8) Zartarian, V. G.; Ozkaynak, H.; Burke, J. M.; Zufall, M. J.; Rigas, M. L.; Furtaw, E. J. A modeling framework for estimating children’s residential exposure and dose to chlorpyrifos via dermal residue contact and nondietary ingestion. Environ. Health Perspect. 2000, 108 (6), 505−514. (9) Bennett, D. H.; Margni, M. D.; McKone, T. E.; Jolliet, O. Intake fraction for multimedia pollutants: A tool for life cycle analysis and comparative risk assessment. Risk Anal. 2002, 22 (5), 905−918. (10) Bennett, D. H.; McKone, T. E.; Evans, J. S.; Nazaroff, W. W.; Margni, M. D.; Jolliet, O.; Smith, K. R. Defining intake fraction. Environ. Sci. Technol. 2002, 36 (9), 206A−211A. (11) Humbert, S.; Marshall, J. D.; Shaked, S.; Spadaro, J. V.; Nishioka, Y.; Preiss, P.; McKone, T. E.; Horvath, A.; Jolliet, O. Intake fraction for particulate matter: recommendations for life cycle impact assessment. Environ. Sci. Technol. 2011, 45 (11), 4808−4816. (12) Solomon, P.; Gehr, P.; Bennett, D.; Phalen, R.; Méndez, L.; Rothen-Rutishauser, B.; Clift, M.; Brandenberger, C.; Mühlfeld, C. Macroscopic to microscopic scales of particle dosimetry: from source to fate in the body. Air Qual., Atmos., Health 2012, 1−19. (13) Klepeis, N. E.; Nazaroff, W. W. Modeling residential exposure to secondhand tobacco smoke. Atmos. Environ. 2006, 40 (23), 4393− 4407. (14) Nazaroff, W. W. Inhalation intake fraction of pollutants from episodic indoor emissions. Build. Environ. 2008, 43 (3), 269−277. (15) Stapleton, H. M.; Eagle, S.; Sjödin, A.; Webster, T. F. Serum PBDEs in a North Carolina toddler cohort: associations with hand wipes, house dust and socioeconomic variables. Environ. Health Perspect. 2012, 120 (7), 1049−1054. (16) Harrad, S.; de Wit, C. A.; Abdallah, M. A. E.; Bergh, C.; Bjorklund, J. A.; Covaci, A.; Darnerud, P. O.; de Boer, J.; Diamond, M.; Huber, S.; Leonards, P.; Mandalakis, M.; Oestman, C.; Haug, L. S.; Thomsen, C.; Webster, T. F. Indoor contamination with hexabromocyclododecanes, polybrominated diphenyl ethers, and perfluoroalkyl compounds: an important exposure pathway for people? Environ. Sci. Technol. 2010, 44 (9), 3221−3231. (17) Jones-Otazo, H. A.; Clarke, J. P.; Diamond, M. L.; Archbold, J. A.; Ferguson, G.; Harner, T.; Richardson, G. M.; Ryan, J. J.; Wilford, B. Is house dust the missing exposure pathway for PBDEs? An analysis of the urban fate and human exposure to PBDEs. Environ. Sci. Technol. 2005, 39 (14), 5121−5130. (18) McKone, T. E.; Castorina, R.; Harnly, M. E.; Kuwabara, Y.; Eskenazi, B.; Bradman, A. Merging models and biomonitoring data to characterize sources and pathways of human exposure to organophosphorous pesticides in the Salinas Valley of California. Environ. Sci. Technol. 2007, 41 (9), 3233−3240. (19) Watkins, D. J.; McClean, M. D.; Fraser, A. J.; Weinberg, J.; Stapleton, H. M.; Sjodin, A.; Webster, T. F. Exposure to PBDEs in the office environment: evaluating the relationships between dust, handwipes, and serum. Environ. Health Perspect. 2011, 119 (9), 1247−1252. (20) Nazaroff, W. W. Inhalation intake fraction of pollutants from episodic indoor emissions. Build. Environ. 2008, 43, 269−277. (21) Bennett, D. H.; Furtaw, E. J. Fugacity-based indoor residential pesticide fate model. Environ. Sci. Technol. 2004, 38 (7), 2142−2152. (22) Average number of people per household, by race and Hispanic origin, marital status, age, and education of householder; U.S. Census Bureau, Washington, DC, 2011. (23) Child-Specific Exposure Factors Handbook; U.S. Environmental Protection Agency, Washington, DC, 2008. (24) PM 2.5 objectives and history; U.S. Environmental Protection Agency, Washington, DC, 2008. (25) Ferro, A. R.; Kopperud, R. J.; Hildemann, L. M. Source strengths for indoor human activities that resuspend particulate matter. Environ. Sci. Technol. 2004, 38 (6), 1759−1764. (26) Exposure Factors Handbook; U.S. Environmental Protection Agency, Washington, DC, 2011. (27) Cohen Hubal, E. A.; Sheldon, L. S.; Burke, J. M.; McCurdy, T. R.; Berry, M. R.; Rigas, M. L.; Zartarian, V. G.; Freeman, N. C. G. Children’s exposure assessment: A review of factors influencing

Despite the lack of model evaluation by case studies, the modeled iF from this study offers key insights on the role of chemical properties and chemical application modes in controlling human exposure to chemicals released indoors. We believe this fugacity model can be useful in screening large numbers of chemicals. iF for a pulse−release scenario can be applied to large numbers of chemicals, provided chemical properties are available. Ultimately, for chemical prioritization, this near-field model needs to be integrated with a far-field model, and the portion of the mass of the chemical released to the indoors as compared to the outdoors must be known. Furthermore, models that account for direct consumer product use, such as personal care products applied directly to the skin, also need to be included for complete prioritization. This model fills an important gap in understanding chemical exposures occurring in the indoor environment.



ASSOCIATED CONTENT

S Supporting Information *

Additional text, eight tables, and four figures with input parameters, fugacity capacities, exposure concentration equations, and mass balance equations of an indoor fugacity model; contribution of each removal pathway to total removal amount; and results of this sensitivity analysis for surface applications. This material is available free of charge via the Internet at http://pubs.acs.org.



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]; phone: 949.648.1614; fax: 530.752.5300. Notes

The authors declare no competing financial interest.

■ ■

ACKNOWLEDGMENTS This research is funded by the American Chemistry Council (Grant 3-DBACC01). REFERENCES

(1) Egeghy, P. P.; Vallero, D. A.; Hubal, E. A. C. Exposure-based prioritization of chemicals for risk assessment. Environ. Sci. Policy 2011, 14 (8), 950−964. (2) Dix, D. J.; Houck, K. A.; Martin, M. T.; Richard, A. M.; Setzer, R. W.; Kavlock, R. J. The ToxCast program for prioritizing toxicity testing of environmental chemicals. Toxicol. Sci. 2007, 95 (1), 5−12. (3) Hubal, E. A. C.; Richard, A.; Aylward, L.; Edwards, S.; Gallagher, J.; Goldsmith, M. R.; Isukapalli, S.; Tornero-Velez, R.; Weber, E.; Kavlock, R. Advancing exposure characterization for chemical evaluation and risk assessment. J. Toxicol. Environ. Health B: Crit. Rev. 2010, 13 (2−4), 299−313. (4) Klepeis, N. E.; Nelson, W. C.; Ott, W. R.; Robinson, J. P.; Tsang, A. M.; Switzer, P.; Behar, J. V.; Hern, S. C.; Engelmann, W. H. The National Human Activity Pattern Survey (NHAPS): a resource for assessing exposure to environmental pollutants. J. Exposure Anal. Environ. Epidemiol. 2001, 11 (3), 231−252. (5) Wallace, L. A. Comparison of risks from outdoor and indoor exposure to toxic chemicals. Environ. Health Perspect. 1991, 95, 7−13. (6) Loh, M. M.; Levy, J. I.; Spengler, J. D.; Houseman, E. A.; Bennett, D. H. Ranking cancer risks of organic hazardous air pollutants in the United States. Environ. Health Perspect. 2007, 115 (8), 1160−1168. (7) Lewis, R. G.; Fortmann, R. C.; Camann, D. E. Evaluation of methods for monitoring the potential exposure of small children to pesticides in the residential environment. Arch. Environ. Contam. Toxicol. 1994, 26 (1), 37−46. 10071

dx.doi.org/10.1021/es3018286 | Environ. Sci. Technol. 2012, 46, 10063−10072

Environmental Science & Technology

Article

dose for the Minnesota and Arizona NHEXAS populations. J. Exposure Anal. Environ. Epidemiol. 2001, 11 (3), 253−268. (45) Tulve, N. S.; Suggs, J. C.; McCurdy, T.; Hubal, E. A. C.; Moya, J. Frequency of mouthing behavior in young children. J. Exposure Anal. Environ. Epidemiol. 2002, 12 (4), 259−264. (46) Wu, X. M.; Bennett, D. H.; Ritz, B.; Frost, J.; Cassady, D.; Lee, K.; Hertz-Picciotto, I. Residential insecticide usage in northern California homes with young children. J. Exposure Anal. Environ. Epidemiol 2011, 21 (4), 427−436. (47) Hellweg, S.; Demou, E.; Bruzzi, R.; Meijer, A.; Rosenbaum, R. K.; Huijbregts, M. A. J.; McKone, T. E. Integrating human indoor air pollutant exposure within life cycle impact assessment. Environ. Sci. Technol. 2009, 43 (6), 1670−1679. (48) Wu, N.; Herrmann, T.; Paepke, O.; Tickner, J.; Hale, R.; Harvey, E.; La Guardia, M.; McClean, M. D.; Webster, T. F. Human exposure to PBDEs: Associations of PBDE body burdens with food consumption and house dust concentrations. Environ. Sci. Technol. 2007, 41 (5), 1584−1589. (49) Watkins, D. J.; McClean, M. D.; Fraser, A. J.; Weinberg, J.; Stapleton, H. M.; Sjodin, A.; Webster, T. F. Impact of dust from multiple microenvironments and diet on PentaBDE body burden. Environ. Sci. Technol. 2012, 46 (2), 1192−1200. (50) Xu, Y.; Cohen Hubal, E. A.; Little, J. C. Predicting residential exposure to phthalate plasticizer emitted from vinyl flooring: sensitivity, uncertainty, and implications for biomonitoring. Environ. Health Perspect. 2010, 118 (2), 253−258.

children’s exposure, and the data available to characterize and assess that exposure. Environ. Health Perspect. 2000, 108 (6), 475−486. (28) Weschler, C. J.; Nazaroff, W. W. SVOC exposure indoors: fresh look at dermal pathways. Indoor Air 2012, in press. (29) Nishioka, M. G.; Burkholder, H. M.; Brinkman, M. C.; Gordon, S. M.; Lewis, R. G. Measuring transport of lawn-applied herbicide acids from turf to home: Correlation of dislodgeable 2,4-D tuff residues with carpet dust and carpet surface residues. Environ. Sci. Technol. 1996, 30 (11), 3313−3320. (30) Nishioka, M. G.; Burkholder, H. M.; Brinkman, M. C.; Lewis, R. G. Distribution of 2,4-dichlorophenoxyacetic acid in floor dust throughout homes following homeowner and commercial lawn applications: Quantitative effects of children, pets, and shoes. Environ. Sci. Technol. 1999, 33 (9), 1359−1365. (31) Fortune, C. R.; Blanchard, F. T.; Ellenson, W. D. Analysis of aged in-home carpeting to determine the distribution of pesticide residues between dust, carpet, and pad compartments; U.S. Environmental Protection Agency, Research Triangle Park, NC, 2000. (32) Diamond, M. L.; Gingrich, S. E.; Fertuck, K.; McCarry, B. E.; Stern, G. A.; Billeck, B.; Grift, B.; Brooker, D.; Yager, T. D. Evidence for organic film on an impervious urban surface: Characterization and potential teratogenic effects. Environ. Sci. Technol. 2000, 34 (14), 2900−2908. (33) Lu, C. S.; Fenske, R. A. Dermal transfer of chlorpyrifos residues from residential surfaces: Comparison of hand press, hand drag, wipe, and polyurethane foam roller measurements after broadcast and aerosol pesticide applications. Environ. Health Perspect. 1999, 107 (6), 463−467. (34) Zartarian, V. G.; Ferguson, A. C.; Leckie, J. O. Quantified dermal activity data from a four-child pilot field study. J. Exposure Anal. Environ. Epidemiol. 1997, 7 (4), 543−552. (35) Reed, K. J.; Jimenez, M.; Freeman, N. C. G.; Lioy, P. J. Quantification of children’s hand and mouthing activities through a videotaping methodology. J. Exposure Anal. Environ. Epidemiol. 1999, 9 (5), 513−520. (36) Freeman, N. C. G.; Hore, P.; Black, K.; Jimenez, M.; Sheldon, L.; Tulve, N.; Lioy, P. J. Contributions of children’s activities to pesticide hand loadings following residential pesticide application. J. Exposure Anal. Environ. Epidemiol. 2005, 15 (1), 81−88. (37) Freeman, N. C. G.; Jimenez, M.; Reed, K. J.; Gurunathan, S.; Edwards, R. D.; Roy, A.; Adgate, J. L.; Pellizzari, E. D.; Quackenboss, J.; Sexton, K.; Lioy, P. J. Quantitative analysis of children’s microactivity patterns: The Minnesota Children’s Pesticide Exposure Study. J. Exposure Anal. Environ. Epidemiol. 2001, 11 (6), 501−509. (38) Auyeung, W.; Canales, R. A.; Beamer, P.; Ferguson, A. C.; Leckie, J. O. Young children’s hand contact activities: An observational study via videotaping in primarily outdoor residential settings. J. Exposure Sci. Environ. Epidemiol. 2006, 16 (5), 434−446. (39) Cohen Hubal, E. A.; Nishioka, M. G.; Ivancic, W. A.; Morara, M.; Egeghy, P. P. Comparing surface residue transfer efficiencies to hands using polar and nonpolar fluorescent tracers. Environ. Sci. Technol. 2008, 42 (3), 934−939. (40) Clothier, J. M. Dermal transfer efficiency of pesticides from new, vinyl sheet flooring to dry and wetted palms; U.S. Environmental Protection Agency, Research Triangle Park, NC, 2000. (41) Nolan, R. J.; Rick, D. L.; Freshour, N. L.; Saunders, J. H. Chlorpyrifos: pharmacokinetics in human volunteers. Toxicol. Appl. Pharmacol. 1984, 73, 8−15. (42) Hore, P.; Zartarian, V.; Xue, J. P.; Ozkaynak, H.; Wang, S. W.; Yang, Y. C.; Chu, P. L.; Sheldon, L.; Robson, M.; Needham, L.; Barr, D.; Freeman, N.; Georgopoulos, P.; Lioy, P. J. Children’s residential exposure to chlorpyrifos: Application of CPPAES field measurements of chlorpyrifos and TCPy within MENTOR/SHEDS-pesticides model. Sci. Total Environ. 2006, 366 (2−3), 525−537. (43) Ozkaynak, H.; Xue, J. P.; Zartarian, V. G.; Glen, G.; Smith, L. Modeled estimates of soil and dust ingestion rates for children. Risk Anal. 2011, 31 (4), 592−608. (44) Buck, R. J.; Ozkaynak, H.; Xue, J. P.; Zartarian, V. G.; Hammerstrom, K. Modeled estimates of chlorpyrifos exposure and 10072

dx.doi.org/10.1021/es3018286 | Environ. Sci. Technol. 2012, 46, 10063−10072