Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

A Decision Analytic Approach to Exposure-Based Chemical Prioritization

  • Jade Mitchell,

    Affiliation Biosystems & Agricultural Engineering, Michigan State University, East Lansing, Michigan, United States of America

  • Nicolas Pabon,

    Affiliation Physics Department, Carnegie Mellon University, Pittsburgh, Pennsylvania, United States of America

  • Zachary A. Collier,

    Affiliation Environmental Laboratory, Engineer Research and Development Center, United States Army Corps of Engineers, Concord, Massachusetts, United States of America

  • Peter P. Egeghy,

    Affiliation Office of Research and Development, United States Environmental Protection Agency, Research Triangle Park, North Carolina, United States of America

  • Elaine Cohen-Hubal,

    Affiliation Office of Research and Development, United States Environmental Protection Agency, Research Triangle Park, North Carolina, United States of America

  • Igor Linkov ,

    Igor.Linkov@usace.army.mil

    Affiliation Environmental Laboratory, Engineer Research and Development Center, United States Army Corps of Engineers, Concord, Massachusetts, United States of America

  • Daniel A. Vallero

    Affiliation Office of Research and Development, United States Environmental Protection Agency, Research Triangle Park, North Carolina, United States of America

Abstract

The manufacture of novel synthetic chemicals has increased in volume and variety, but often the environmental and health risks are not fully understood in terms of toxicity and, in particular, exposure. While efforts to assess risks have generally been effective when sufficient data are available, the hazard and exposure data necessary to assess risks adequately are unavailable for the vast majority of chemicals in commerce. The US Environmental Protection Agency has initiated the ExpoCast Program to develop tools for rapid chemical evaluation based on potential for exposure. In this context, a model is presented in which chemicals are evaluated based on inherent chemical properties and behaviorally-based usage characteristics over the chemical’s life cycle. These criteria are assessed and integrated within a decision analytic framework, facilitating rapid assessment and prioritization for future targeted testing and systems modeling. A case study outlines the prioritization process using 51 chemicals. The results show a preliminary relative ranking of chemicals based on exposure potential. The strength of this approach is the ability to integrate relevant statistical and mechanistic data with expert judgment, allowing for an initial tier assessment that can further inform targeted testing and risk management strategies.

Introduction

Manufactured chemicals are widely used in products such as cosmetics, plastics, and electronics, and have applications in almost all industrial processes in sectors including energy, agriculture, and pharmaceuticals [1]. Increasing dependence on manufactured chemicals has not, however, been matched by an adequate increase in our understanding of the risks these may pose to the environment and human health [2]. Many chemicals in U.S. commerce today have unknown environmental fates and poorly understood potential for human exposure, including some of the most ubiquitous commercial chemicals, such as surfactants, fragrances, cleaning agents and pesticides [3], [4]. In this context, exposure is the contact of a stressor (i.e., a chemical agent) with a receptor (i.e., a human or a human population) for a specific duration of time [5]. Because of the lack of resources and sufficient scientific information on toxicity [6] and exposure [3] for the assessment of all chemicals, efforts are typically, and rationally, devoted to assessing those chemicals believed to pose the greatest potential risks based on production volume and chemical properties.

Within the domain of human health risk assessment, toxicity is an indication and measurement of the severity of adverse health effects a chemical causes in relation to an exposure level (dose). We broadly define exposure to be the contact of a stressor with a receptor for a specific duration of time [5]. The stressors of interest are chemical agents that can potentially lead to an adverse impact and the receptors of interest are individuals or population of individuals. Exposure is complex and dynamic in nature due to its spatial and temporal characteristics. For this reason, exposure-based prioritization efforts focus on relative exposure potential as a means to evaluate and rank chemicals. While prioritization is in of itself a risk management strategy, other risk management decisions may follow to include the allocation of scarce resources to complete future risk assessments, collection of additional data or testing, and/or (bio) monitoring. Therefore, the resolution and precision of the data incorporated in these efforts may vary according to the overall objective of the prioritization.

The U.S. EPA Office of Chemical Safety and Pollution Prevention recently performed a chemical prioritization exercise to identify 83 “TSCA Work Plan Chemicals” [7] as candidates for risk assessment during the next few years. Broad stakeholder input was used to identify prioritization and screening criteria and data sources. Chemicals were evaluated based on their combined hazard, exposure potential, and persistence and bioaccumulation characteristics using a two-step process. In the first step, a set of data sources was used to identify 1,235 chemicals meeting one or more criteria suggesting concern, namely: known reproductive or developmental effects; persistent, bioaccumulative, and toxic (PBT) properties; known carcinogenicity; and presence in children’s products. Excluding those chemicals not regulated under TSCA and those with physical and chemical characteristics that do not generally present significant health hazards narrowed the number of chemicals down to 345 candidates. In the second step, a numerical algorithm was used to score each chemical based on three characteristics: hazard, exposure, and potential for persistence or bioaccumulation. Candidate chemicals that ranked highest on the basis of their total score were identified as work plan chemicals; those that could not be scored because of an absence of exposure or hazard data were identified as candidates for information gathering.

Using the methodology described above, EPA has been able to identify a priority set of chemicals for near-term assessment based on criteria widely accepted as warranting concern. The scoring algorithm is transparent and the data sources are well documented. Focusing on chemicals with documented evidence of concern (i.e. “data-rich”) is reasonable in light of limited prototypes for post hoc screening and the paucity of available resources. However, this approach may not adequately address the need to make decisions about the thousands of chemicals in commerce and the hundreds of new chemicals introduced each year for which there is little or no information [1], [3].

To support the development of novel rapid approaches for evaluating potential exposure of both existing and emerging chemicals, the EPA has initiated the ExpoCast research program [8]. This program is keenly interested in characterizing exposures across the chemical life cycle –manufacturing, transportation, product formulation, consumer product usage and finally disposal. EPA seeks to build on current chemical exposure models and knowledge to generate robust new protocols that better support chemical evaluation, risk assessment and risk management. Recent activities under this program have evaluated utility of available approaches for the purpose of rapidly prioritizing large numbers of chemicals on the basis of exposure [9], [10].

A number of exposure models were recently comparatively evaluated through the EPA Expocast model challenge, where a set of approximately 50 data-rich chemicals of different classes were ranked by several different approaches [10]. The chemicals were chosen to include high interest chemicals with a range of properties. Each modeling approach was capable of analyzing a different number of chemicals from the full set because of varying input requirements. Key findings of the comparative analysis among the prioritization schemes indicated significant differences in chemical ranking as a result of several factors: (1) which processes the model described across the source to effects continuum [11]; (2) the exposure metric or surrogate metric used for prioritization and which statistic (i.e., median, upper bound or lower bound estimate); (3) whether the model inputs included actual, modeled or unit emissions; (4) which exposure pathways were considered (i.e., from aggregated sources or through a dominant pathway); and (5) which type of exposure scenarios were considered (i.e., direct or indirect, diffuse source or concentrated source, etc.) [10]. Only mechanistic models characterizing exposure associated with environmental sources could rapidly evaluate and rank potential exposure for the majority of chemicals. To a great extent, this was due to both the minimum data requirements and the availability of predictive tools (i.e., QSARs) to generate model inputs that could be used to describe fate and transport under steady state and equilibrium conditions. Of the other models evaluated in the EPA Expocast model challenge, those designed for evaluation of chemicals in specific exposure scenarios lacked data for chemical and scenario specific input parameters and were thereby inhibited in their ability to produce ordinal rankings for the 55 chemicals.

Arguably, one of the major limitations of the models evaluated, and perhaps one of the larger knowledge gaps in exposure-based chemical prioritization itself, involves complex social behaviors that determine how humans come in contact with manufactured chemicals, particularly those emanating from near field sources (e.g., residential and consumer products). Thus there is a pressing need for enhancing current approaches with tools and techniques developed for understanding human behaviors, such as human factors engineering and marketing research, to better define scenarios describing how products are used. Accurate use scenarios among population groups of interest are necessary to properly characterize the consumer use component of a chemical’s life cycle.

Decision support tools borne out of the social sciences may also have a place in chemical prioritization. Multi-criteria Decision Analysis (MCDA), a rule-based method of classification for priority setting, is both a set of techniques and an approach for ranking alternatives [12], [13]. MCDA is a promising approach for exposure-based prioritization because it is transparent and understandable, yet complex and rigorous enough to include scenario-based reasoning, stochastic processes and value of information analysis. Moreover, it is amenable to sparse data [14], [15], [16], [17]. These characteristics complement some of the limitations of currently available statistical, mechanistic, or logic models, which provide useful frameworks for gathering relevant data but lack the social and policy context for risk-informed decision making. MCDA can merge a variety of types of exposure metrics from descriptions of physical chemical properties to the socioeconomic measures which characterize human activity, chemical use and contact to ultimately inform screening level risk estimates. Permitting structured integration of different types of information, MCDA methods provide a means for combining quantitative chemical property, production and use data with expert judgments and stakeholder preferences. MCDA assessment criteria can be adaptively weighted and modified in real time to evaluate both data-rich and data-limited chemicals.

Use of MCDA methods to support prioritization decision making under high uncertainty has been demonstrated many times including hazard identification and assessment. Risk management alternatives of industrial hazards or industrial consequences were relatively ranked using an MCDA approach by Paralikas and Lygeros [18]. The method recognizes that a single factor could not be used to define flammability and that different methods, tools, codes and legislation use varying sets of fire hazard properties as an example. Using the MCDA framework, the different decision criteria were successfully integrated using fuzzy logic to deal with linguistic variables and uncertainties allowing broad application for chemical hazard ranking decisions. In another example, life cycle assessment (LCA) was incorporated within a decision framework to prioritize future research and evaluate sensitivities to missing information in an assessment of processes for synthesizing single walled carbon nanotubes [14]. Engineered nanomaterials present uncertainties similar to chemicals in consumer products in terms of unknown environmental and human health across all life stages from formulation to disposal.

This paper demonstrates how analytical tools, such as LCA and MCDA, can offer a versatile and transparent approach to exposure-based prioritization utilizing results from several approaches evaluated in the EPA ExpoCast model challenge. The purpose of prioritization within this context is to focus resources on further evaluation of safety for chemicals with high potential for exposure and risk. A combination of exposure assessment model output with qualitative exposure criteria within such a decision framework has been recommended in the exposure-based waiving protocol within Europe’s REACH Regulation [19] which shares some similar goals for human and environmental health protection.

Materials and Methods

We propose a decision analytic approach for exposure-based chemical prioritization to address the need for novel, rapid exposure potential screening protocols. In this approach, we build on current research and existing models by evaluating relevant chemical exposure criteria within a larger MCDA framework. We employ a two-part prioritization model that incorporates both properties of the chemical itself and properties of the chemical’s life cycle (Figure 1).

thumbnail
Figure 1. MCDA Framework for Exposure-Based Chemical Prioritization.

https://doi.org/10.1371/journal.pone.0070911.g001

The chemical property and life cycle property assessments are structured to analyze exposure-related information associated with specific chemical properties and distinct life cycle phases, respectively. Relevant chemical and life cycle properties are grouped into several criteria based upon the means by which each property contributes to the chemical’s overall exposure potential (e.g., properties associated with a chemical’s ability to bioaccumulate vs. those associated with its ability to be metabolized by the human body). Chemical and life cycle properties in each criterion are then further divided into various sub-criteria. The numerical values associated with these properties for a given chemical serve as inputs to the model. Input data can be obtained from a number of different sources, including existing databases, current literature and expert judgment. The criteria within this decision model were selected by reviewing those used in the models submitted to the ExpoCast model challenge, [10] and then structured into a hierarchical framework based on discussions with exposure science experts.

Within each sub-criterion, the constituent chemical or life cycle property is evaluated to determine its contribution to overall exposure potential. Input values for individual properties are compared against established numerical thresholds, which define distinct levels of risk that span the range of possible values for the given sub-criterion. Thresholds are used to score property values based on the indicated level of risk (e.g., a compound with a longer half-life may have higher potential for exposure than a compound with a shorter half-life, all other things being equal).

Following an MCDA approach, sub-criterion scores are then combined according to explicit decision rules to derive scores for their higher-level criterion. Chemical property and life cycle phase criterion scores are then combined to produce a Chemical Properties Exposure Score (CPES) and a Life Cycle Exposure Score (LCES) for each chemical. These scores reflect relative estimates of chemical exposure potential as indicated by available chemical property and life cycle property data, respectively. Exposure scores may then be integrated to derive aggregate measures of exposure potential, which can be used to compare and prioritize chemicals on a relative basis, or can remain separate and be plotted on a risk matrix for a more qualitative assessment.

Chemical property and life cycle phase criteria can be weighted within each assessment to reflect their relevance to the user’s management objectives. Weights may indicate a specific focus of the assessment or reflect expert judgment of a criterion’s predictive reliability or relative importance. Criterion weights can be adjusted to refine the scope of a particular assessment to a particular class of chemicals (e.g., pesticides), a particular exposure scenario (e.g., occupational exposure), or a particular exposure target (e.g., environmental contamination). When eliciting subjective weights, it is important to utilize best practices to avoid potential biases and inconsistencies [20], [21]. Numerous elicitation techniques exist, including rank-based methods and swing-weight methods [13], [21], [22].

Chemical Properties Assessment

As seen in Figure 1, the Chemical Properties Assessment considers four main criteria to estimate potential risk for human exposure: bioaccumulation potential, persistence, ADME (Absorption, Distribution, Metabolism, and Elimination), and physical hazard potential. Each criterion constitutes a unique set of sub-criteria, which define the distinct chemical property data points that serve as inputs to the assessment. Observed chemical properties used to estimate exposure potential are defined by the specific sub-criteria under each of the four main criteria. Using thresholds established for each sub-criterion, individual data points are evaluated and assigned scores representing the potential for exposure indicated by the observed chemical property. Once these initial scores have been calculated, the highest within each set of sub-criteria is assigned as that criterion’s exposure score.

When certain chemical-specific data are unavailable, as is often the case in this context, it may not be possible to assign scores to each sub-criterion. By defining each criterion’s exposure score as the highest of its associated sub-criteria scores, we account for this possibility. By employing this approach, criterion scores can be assigned even in the presence of sparse data.

Each chemical’s bioaccumulation, persistence, ADME, and physical hazard scores are combined with their associated weights. Weighted criteria exposure scores are then summed to produce initial chemical property exposure score for each chemical. Once this has been done for the set of chemicals being assessed, the initial chemical property exposure scores are normalized from 0 to 1 to produce relative rankings.

Bioaccumulation.

Bioaccumulation is a process in which a chemical substance is absorbed by an organism via all routes of exposure in the natural environment, for example through dietary and ambient environmental sources, and increases in concentration over time [23]. Using three bioaccumulation-related sub-criteria, we evaluate surrogate chemical properties in order to predict the compound’s ability to bioaccumulate.

Bioconcentration factor (BCF).

A compound’s BCF is a dimensionless number representing the relative concentration of the compound in organic tissues. In general, chemicals with relatively higher BCFs have greater potential for exposure, and thus are more likely to adversely impact human health and the environment. In this model, four distinct numerical thresholds were used to evaluate chemical BCF data. These thresholds are shown in Table 1, and were used to assign each chemical a BCF sub-criteria score from 1–4 based on the indicated level of bioaccumulation potential. Thresholds are based on previously published values employed by existing exposure assessment models: the EPA Design for the Environment Program [24], and the Clean Production Action’s Green Screen for Safer Chemicals Initiative [25]. To address minor numerical discrepancies, the more conservative thresholds were chosen when values differed between models.

thumbnail
Table 1. Thresholds for Bioaccumulation Potential and Environmental Persistence.

https://doi.org/10.1371/journal.pone.0070911.t001

Log kow.

A compound’s Kow, or octanol-water partition coefficient, describes its ability to transition between water and carbon-based media. Chemical compounds with relatively higher log Kow are capable of greater movement within the environment; they are thus more adaptive and have higher potential for human exposure and absorption. In this model, four distinct numerical thresholds were used to evaluate chemical Kow data. These thresholds are shown in Table 1, and were used to assign each chemical a log Kow sub-criteria score from 1–4 based on the indicated level of bioaccumulation potential. Thresholds are based on previously published values employed by existing exposure assessment models: the EPA Design for the Environment Program [24], and the Clean Production Action’s Green Screen for Safer Chemicals Initiative [25], with the more conservative threshold chosen when values differed between models.

Molecular weight.

Previous studies have identified a significant correlation between a compound’s molecular weight and its ability to bioaccumulate [26], [27]. Results from these studies support the general conclusion that heavy molecules do not easily bioaccumulate, as their size hinders passage through lipid membranes. Lower weight chemicals thus possess a relatively greater potential for human exposure. These and similar findings have been used to inform chemical testing policy and legislation such as the OECD Chemical Substance Control Law (CSCL) in Japan [28] and the EPA Toxic Substances Control Act (TSCA) in the United States [29].

A single cut-off threshold is employed by our model to evaluate molecular weight data. Molecules 1000 amu or greater are given a bioaccumulation criteria score of 1, regardless of their other sub-criteria scores within the bioaccumulation category (BCF & log Kow). The 1000 amu cut-off follows TSCA premanufacture notification policy [29], and is based on current understanding that molecular weights in this range are generally better indicators of chemical bioaccumulation potential than other surrogate properties [26].

Persistence.

Persistence corresponds to the length of time a chemical can exist in the environment before degrading or being transformed by natural processes [23]. Persistent chemicals are more likely to come into contact with humans compared to chemicals that degrade quickly in the environment. We consider the half-life in water, soil, sediment, and air for each chemical as surrogate indicators of persistence for the purpose of evaluating exposure potential.

The numerical thresholds used for evaluating chemical half-life data are shown below in Table 1. Thresholds were used to assign each chemical four distinct half-life sub-criteria scores from 1–4 based on the level of persistence indicated by each of the four half-lives (in water, soil, sediment, and air). Threshold values for water, soil, and sediment are based on previously published values employed by existing exposure assessment models: the EPA Design for the Environment Program [24], and the Clean Production Action’s Green Screen for Safer Chemicals Initiative [25], using the more conservative thresholds. The threshold value for air follows science-based guidance for evaluating chemical long-range transport potential and overall persistence [30]. Chemicals with half-lives in air that are less than two days are assigned an associated sub-criteria score of 1 (“Low”), while those with half-lives in air greater than or equal to two days are assigned an score of 3 (“High”).

ADME.

Properties that describe a chemical’s ability for absorption, distribution, metabolism, and excretion (ADME) are indicators of the potential for biologically relevant human exposure. Chemicals that can be easily absorbed by the body and that are resistive to metabolism or excretion pose a greater threat for extended exposure; therefore it is useful to focus on the entrance and exit of the chemicals within the context of the body. Though recent and current ADME-related research efforts have focused on establishing appropriate surrogate properties and developing predictive models, general consensus has not been reached regarding an accepted approach to ADME assessment for environmental chemicals [10]. Building on current research and existing models, a new ADME assessment protocol intended for screening-level exposure-based chemical prioritization was incorporated into the framework [10]. This method utilizes QikProp software Version 3.0 [31], a QSAR-based model to obtain surrogate chemical property values, which were then integrated to evaluate ADME properties along various sub-criteria briefly discussed below. All QikProp values are based on a 24-hour exposure period. Incidentally, QikProp is a three-dimensionally based structure method, so the SARs depend on the solvent accessible surface area. The properties calculated are dependent on the conformer adopted at the time of calculation and could be sensitive to molecular orientation. In addition, QikProp was designed exclusively to develop organic pharmaceutical compounds, so cannot be used for metals and inorganic compounds. Thus, if the analytics discussed herein are to be applied to metals and inorganic compounds, another QSAR system is needed.

Absorption.

The chemical absorption assessment is based on two QikProp predictors which describe oral availability. The first descriptor represents a qualitative measure of oral absorption potential, and takes values of 1, 2, or 3 for low, medium, or high, respectively. The second descriptor represents a numerical probability of oral absorption on a 0 to 100% scale, with <25% and >80% designating low and high probability, respectively. These values were combined to derive an absorption score (1–3) for each chemical.

Distribution/Excretion: Distribution and excretion-related properties were combined into a single assessment. QikProp predicted octanol/water partition coefficients, serving as surrogates for half-life within the human body, were categorized into bins using subjective thresholds to derive a distribution/excretion score (1–4) for each chemical.

Metabolism.

The assessment of metabolism was derived from the QikProp descriptor representing the number of expected possible metabolites for each chemical over a 24-hour period in the human body. These values were categorized based on the predicted half-life of each chemical in order to represent metabolism via natural degradation in the body. These values were combined to generate average metabolism scores (1–4) for each chemical.

Physical hazard potential.

Highly flammable and reactive chemicals pose human and environmental threats that may not be considered in standard exposure or toxicity-based assessments. Though the properties that determine a given chemical’s flammability and reactivity may be distinct from those that determine its environmental fate and transport, the threat of physical hazard is nonetheless directly related to the likelihood of exposure. The risk of physical hazards (e.g., combustion) is thus an exposure-related risk, and we assess each chemical’s hazard-related properties in order to anticipate threats that may not be considered in other exposure or toxicity-based screenings. In accordance with existing National Fire Protection Association (NFPA) standards and classifications [32], flammability and reactivity were assigned scores of (1–4) using established NFPA thresholds.

Chemical Life Cycle Properties Assessment

Similarly to the assessment of chemical properties, we estimate potential for human exposure by assessing three main life cycle phases of manufactured chemicals: production, consumer use, and disposal. Each phase constitutes a unique subset of exposure-related criteria, which define the distinct life cycle characteristics that serve as inputs to the assessment.

The different criteria associated with each of the three life cycle phases designate the individual life cycle properties that will serve as indicators of a chemical’s exposure potential during the relevant phase. All life cycle criteria are evaluated quantitatively, with higher values indicating higher potential for exposure. Instead of establishing thresholds for each sub-criteria as in the assessment of chemical properties, raw values are used but then normalized across the set of chemicals for each individual sub-criteria. This provides bounds for the range of values and assists in making comparative assessments.

Criteria scores are then calculated by summing the sub-criteria scores. Again, these scores are normalized across the set of chemicals to account for criteria containing more sub-criteria than others, and then multiplied by their weights to produce an initial Life Cycle Properties Exposure Score (LCES). Once initial LCESs have been calculated for all chemicals, we derive final LCESs by normalizing initial scores to the highest and lowest observed scores across all chemicals.

Production.

Number of Potential Exposure Sources: Each chemical is evaluated to determine the possibility for human exposure during processes associated with production of the chemical. We consider one potential source (occupational microenvironments) defined as any workplace environment in which a release might occur during chemical manufacture and/or processing. Each chemical is assigned a score of either 0 or 1 based on whether the compound presents risk of exposure during production.

Projected average annual number of production sites.

A chemical’s exposure risk is increased if it is produced in many locations. Ubiquity classifications for each chemical were used to estimate the amount of chemical production sites [10]. Higher scores indicate increased potential for human exposure during chemical production: very widespread (5), widespread (4), moderate (3), localized (2), low (1).

Regional geometric mean production quantity (MQR).

In addition to how widespread production is, estimates are made of the quantity produced. This is estimated using the Regional Geometric Mean Production Quantity (MQR), measured in units of kilotons per year. This is an estimated quantity, but production quantities could also be provided by industry.

Consumer Use

The assessment evaluates several sub-criteria relevant to the consumer use phase in the life cycle of manufactured chemicals. Based on the intended uses of each chemical, primary consumer class is defined as either strictly industrial, or industrial and individual. Chemicals used during industrial processes (e.g., monomers, solvents) and chemicals otherwise noted to have primarily industrial consumers were defined to have a strictly industrial consumer class. Chemicals used in agriculture (e.g., pesticides, insecticides, herbicides) or as food/cosmetic additives (e.g., preservatives, anti-microbials) were defined to have both industrial and individual consumers. Chemicals directly incorporated into consumer products during their production (e.g., plastics, coatings, fabrics, flame retardants) are also defined to have both industrial and individual consumers.

Number of potential exposure sources.

Each chemical was evaluated to determine the possibility for human exposure during processes associated with both industrial and individual consumer uses of the chemical. Ten distinct potential sources associated with consumer exposure were considered (i.e., outdoor air, water, soil, biota, indoor air/dust, in-vehicle air, object contact, tap water, other water, food/beverages) by assigning each chemical a score from 0–10 based on possibility for exposure via each unique source during consumer use of the compound.

Projected average annual number of individual consumers.

Chemicals defined as having industrial and individual consumer classes were assessed to determine their potential for exposure to individual consumers in non-industrial settings. Chemical ubiquity classifications were used to represent the relative size of each chemical’s average, annual, individual consumer base. Chemicals defined as having strictly industrial consumer classes were assigned individual consumer scores of 0. Remaining chemicals were assigned scores from 1–5 based on their ubiquity, with higher scores indicating increased potential for individual consumer exposure during non-industrial use: very widespread (5), widespread (4), moderate (3), localized (2), low (1).

Projected average annual number of industrial consumers.

To assess chemicals’ potential for exposure to industrial consumers, we employ the ubiquity classification to estimate the average, annual size each chemical’s industrial consumer base. As none of the chemicals assessed were defined as having a strictly individual (non-industrial) consumer base, all chemicals were assigned scores from 1–5 based on their ubiquity classification, with higher scores indicating increased potential for industrial consumer exposure during use of the chemical: very widespread (5), widespread (4), moderate (3), localized (2), low (1).

Projected average annual quantity consumed per individual/industrial consumer.

The average annual quantity of each chemical consumed per consumer was predicted using the relative size of the chemical’s total consumer base (including both individual and industrial consumers), and its MQR. Relative measures of consumption quantity per consumer (Q) were calculated by dividing each chemical’s projected mean production volume by their total number of consumers, assuming chemicals with higher consumption quantities to have increased potential for consumer exposure. Projected annual quantities consumed per individual consumer were calculated using the same equation as that for industrial consumers:(1)where (niIndividual+niIndustrial) represents the chemical’s total consumer base, or the number of individual consumers plus the number of industrial consumers.

Susceptible populations.

To determine if there was a heightened exposure risk to susceptible populations (in this case, children), particular processes associated with individual consumer use of the chemical were evaluated. Nine distinct potential sources associated with exposure to children were considered (Outdoor Air, Water, Soil, Indoor Air/Dust, In-Vehicle Air, Object Contact, Tap Water, Other Water, and Food/Beverages), and each chemical was assigned a score from 0–9 based on possibility for exposure via each unique source.

Disposal

Number of potential exposure sources.

Each chemical was evaluated to determine potential for human exposure resulting from disposal events. We consider four distinct disposal-related sources (Outdoor Air, Water, Soil, Biota), assigning each chemical a score from 0–4 based on potential for exposure via each unique source during and after disposal of the compound.

Projected average annual number of disposal events.

Each chemical’s total number of consumers was estimated to determine an annual number of associated chemical disposal events. Assuming that each chemical’s industrial and individual consumers dispose of equal amounts of the compound, we define the projected number of disposal events as each chemical’s total number consumers, and assign scores of 1–10, with higher scores representing greater potential for disposal-related human exposure.

Projected average annual quantity disposed.

To account for assumed variations in the actual quantities disposed during industrial and individual consumer disposal events, we assume that 0.1% of the net production volume of each chemical is disposed of in order to evaluate disposal-related exposure potential. Note that the use of this unit value assumes that no chemical- or product-specific data were available. With larger disposal quantities indicating higher potential for post-disposal chemical exposure, we calculate relative disposal quantities of each chemical (QDISP) as:(2)

Integrating Chemical Properties and Life Cycle Exposure Scores

Once assessments of chemical properties and life cycles have been performed on all chemicals, those chemicals lacking sufficient data to calculate either a chemical properties exposure score or life cycle exposure score are removed from the remainder of the prioritization. Though these chemical’s available scores may indicate significant threat of exposure, they are excluded from the integration process as their scores can skew final exposure potential relationships. The remaining chemicals are renormalized as:(3)where xES denotes the relevant exposure score (either chemical or life cycle). Next, the remaining chemicals’ exposure scores (chemical property and life cycle property) are summed to produce aggregate exposure scores. These scores represent cumulative measures of exposure potential based on each chemical’s distinct properties and characteristics of its projected life cycle. Aggregated exposure scores, which all lie in the range of 0–2, are used to numerically rank chemicals based on their potential for human exposure.

In addition to this quantitative integration, chemical property and life cycle scores can be visualized using a risk-reporting matrix (Figure 2) for a more qualitative assessment of aggregate chemical exposure potential.

thumbnail
Figure 2. Example Chemical Exposure Potential Risk Matrix.

https://doi.org/10.1371/journal.pone.0070911.g002

In this method of integration, chemical property and life cycle exposure scores are converted from a scale of 0–1 to a scale of 0–5 by multiplying the initial score by a factor of five to place them within the 5×5 risk matrix, with each chemical’s position representing a qualitative, cumulative measure of exposure potential based on both chemical and life cycle properties. Qualitative exposure potential thresholds (red, yellow, or green) can be defined within the matrix to designate high, moderate, and low risk regions.

Case Study

Data Set

For the case study, a set of 51 chemicals was selected from those presented and evaluated in the model challenge (Table 2), representing a wide variety of chemical classifications (e.g., organics, metals, etc.). Sub-criteria scores for these chemicals were collected from numerous reports and online databases, and the sources for each sub-criterion are listed in Table 3. Case study data can be found in the online File S1.

Prioritization

First, the data for each chemical was compiled. It was found that some chemicals were difficult to assess due to a lack of readily available data. If a chemical did not have any sub-criteria scores for at least one of its criteria, that chemical was removed from the analysis process as having too little data for analysis. Nine of the 51 chemicals (largely metals) were removed for this reason.

Following the MCDA approach outlined above, each of the remaining test chemicals was assessed. Scores for each criterion were weighted by allocating equal weights (i.e., bioaccumulation, persistence, ADME, and physical hazards each weighted 25%; production, consumer use, disposal each weighted 33.33%). The final prioritization under this weighting distribution is shown in Table 4. The risk matrix comparison under this weighting distribution is shown in Figure 3.

thumbnail
Figure 3. Risk Matrix Comparison of Exposure Potential with Even Weighting.

https://doi.org/10.1371/journal.pone.0070911.g003

Discussion

As stated above, one of the major limitations of currently available exposure models involves the inability to fully characterize the influence of complex social behaviors on resulting exposures or contact between humans and manufactured chemical across all life stages of the chemical. This is especially true for chemicals used in residential and consumer products, those arising from near field sources. A multi-criteria decision model was developed to combine typical physiochemical screening level data with measures to characterize human activities. As a proof of concept to show the utility of this approach, a case study was conducted on a small set of chemicals that were also analyzed using higher tiered statistical and mechanistic exposure models in a model challenge [10]. The models used in the model challenge considered different types of exposure scenarios including indirect exposures from diffuse environmental sources and direct, concentrated exposures from micro-environmental sources (i.e. from a personal care product or within a residence), though the latter had significant limitations in terms of necessary data to produce exposure estimates. Ranking results were obtained by three models and the comparative analysis is reported elsewhere [10]. Some agreement between ranking results was observed, but in general these models produced widely incongruous results across a number of different domains of information. Interestingly, some of the results using the MCDA model developed herein coincide with results from these more complex models. The majority of the chemicals (13 of 14) ranked in the top one-third of the list in Table 4 (Rank 1–14), are also ranked in the top one-third of one of the models evaluated in the challenge. In general this agreement is with a “far-field” indirect diffuse source model which does not incorporate human activity at the micro-environmental level. Nonylphenol was the exception as it was ranked low by all other mechanistic models. Similarly, the bottom third of the ranked list in Table 4 (Rank 28–42) shows high agreement with results from a model from the challenge. One model used characterized both far-field and near-field exposures and the other two were far-field models.

Because this analysis was conducted as a proof of concept, an exhaustive search for quality data and subsequent data validation was not conducted independently of the model challenge. However, the absence of the mechanistic relationships involved in the exposure models as well as the equal weighting scheme used in our example would lead to the assumption that the input drivers of the challenge models would be different than the input drivers of MCDA model. To fully explore this assumption and the utility of this methodology for larger scale research prioritization or policy guidance, the results of the case study underscore the need for quality data inputs. Only nine of the chemicals had to be excluded. These chemicals have properties that exclude them from the domain of applicability of the analytics, e.g. models, QSAR type, and other tools. As mentioned, metals and inorganic compounds are not characterized by the ADME models used in this study.

For the majority of compounds that fall within the domain of applicability, the MCDA approach is useful. As shown in Table 4, the majority of the chemicals used in plastics appear in the top half of the ranked list denoting highest exposure potential by highest aggregated exposure score. Plastics are broadly related to exposures that occur in all locations across the life-cycle of the chemicals. The chemicals in the bottom half of the ranked list (lower exposure potential) fit into a number of other of categories, but 11 of 21 are or were used as pesticides/herbicides, agriculturally, in homes or in public and commercial areas. The two pesticides/herbicides, Parathion and Methoxychlor, are ranked relatively low on the list in Table 4. Both chemicals were exclusively used in agriculture only, but have been previously banned or restricted by the EPA and do not have other uses like 1,2,3-trichlorobenzene, ethylene thiourea, and hexachlorobenzene which were also used exclusively in agriculture but are now used as a nonfood commercial additives. The remaining chemical in the agricultural only category is aldicarb. Aldicarb was restricted more recently in 2010 and will not be completely phased out until 2018, so exposure potential may be higher than the others in this category.

It should be noted that the nature of this analysis is to score chemicals in a comparative and relative manner, as opposed to assigning an absolute measure of exposure risk, which would not be practical or appropriate for a screening tool such as this. The relative assessment of chemical exposure potential is therefore dependent upon the set or sub-set of chemicals under consideration, and must be considered when designing the analysis and interpreting the results.

If a risk matrix is used for interpretation or communication of exposure potential results, it is important to note that a chemical with a high chemical property score and low life cycle property score (or vice versa) may be displayed has having a low exposure risk. When the risk matrix is used for score integration, however, these chemicals will appear on the boundaries of the matrix and can easily be identified as outliers that may warrant further assessment. Figure 3 shows the results of the case study on such a risk matrix. The risk matrix approach can be used to graphically visualize qualitative risk categories such as high, medium and low risk. The case study chemicals mostly fall within the same middle risk range of the matrix. Six chemicals fall into the higher exposure risk potential category and seven chemicals fall into the low exposure risk potential category based on the delineations shown in Figure 2. As a high tier screening, this type of representation may be useful for rapid visualization and categorization of large number of chemicals; however risk matrices should be used with caution when guiding risk management decisions [35].

Both the ranking and risk matrix approaches highlight the potential promise of multi-criteria decision analytic models for exposure-based prioritization, but further development beyond this effort is warranted. Given that the baseline weighting scenario – equal weights distributed among the chemical property and life cycle criteria – is likely an unrealistic one, a sensitivity analysis should be conducted to explore the effects of uncertainty in both the scoring of chemical parameters and the weighting schemes on the final chemical prioritization. This will help identify chemicals which are targets for further exposure assessment and data collection, ideally including better release characterization, proximal exposure assessment, and biomonitoring.

Finally, it is important to recognize that these results are strictly a measure of exposure potential and do not consider toxicological properties. Risk is a function of both hazard and exposure. The means by which organisms are exposed to stressors are complex; with many feedback loops (e.g., an outcome may itself become a stressor or modify other stressors). Risks related to chemical ingredients in products depend not only on the inherent properties of that chemical, but also the manner in which the chemical is formulated and used. Exposure potential therefore might be integrated with computational toxicology to paint a more complete picture of risk and to effectively prioritize the numerous chemicals in commerce.

Conclusions

In this paper, we have presented a decision analytic approach to exposure-based prioritization of manufactured chemicals. The proposed methodology allows for structured and transparent analysis of chemical exposure potential through integration of heterogeneous metrics used to evaluate exposure risk-related information associated with both chemical properties and life cycle phases. The model is scalable to assess as many chemicals as is necessary for the project scope, and the MCDA framework is able to accommodate varied inputs and exposure potential indicators, providing an adaptive and easy-to-use screening tool for rapid prioritization in the face of sparse data. In addition, the use of weighting in the model allows for specific user objectives, expert judgment, and data availability considerations to be explicitly implemented within the assessment.

The proposed approach builds on earlier models and current research relating to rapid evaluation of exposure potential. Specifically, it integrates the results of mechanistic and statistical approaches with semi-quantitative categorical data to describe exposure potential. In this paper, we attempt to address the need for high-level screening tools that (1) are capable of more detailed assessments than those provided by simpler predictive models (i.e., limited to persistence and bioaccumulation as indicators of exposure), and (2) have less intensive data requirements than more complex models, so as to remain efficient at the screening level.

It is important to note that work on this model is ongoing, and that the initial framework presented in this paper is primarily intended to illustrate the application of decision analytic methods to supplement existing exposure potential estimation techniques. Currently, our developmental efforts are focused on: (1) refining ADME assessment criteria and calculations; (2) identifying optimal surrogates for bioaccumulation potential; (3) implementing value of information (VOI) techniques to quantify data gaps and prioritize further research efforts; (4) improving normalization algorithms; and (5) developing a supplemental logic model for more specific exposure scenario evaluation. Additionally, we are working to develop formal means of considering expert judgment and empirical chemical exposure data within our assessments. In the future, we anticipate that the decision analytic approach will be able to provide decision makers with important and reliable information to support efficient, exposure-based prioritization of manufactured chemicals.

Supporting Information

Acknowledgments

This material has been subjected to Environmental Protection Agency review and approved for publication. Mention of trade names or commercial products does not constitute endorsement or recommendation for use. Permission was granted by the US Army Corps of Engineers, Chief of Engineers to publish this material. The views and opinions expressed in this paper are those of the individual authors and not those of the US Army, or other sponsor organizations.

Author Contributions

Conceived and designed the experiments: JM PPE IL. Performed the experiments: JM NP ZAC. Analyzed the data: NP ZAC. Contributed reagents/materials/analysis tools: PPE ECH IL DAV. Wrote the paper: JM NP ZAC PPE ECH IL DAV.

References

  1. 1. Wilson MP, Schwarzman MR (2009) Toward a new U.S. chemicals policy: rebuilding the foundation to advance new science, green chemistry, and environmental health. Environ Health Perspect 117 (8): 1202–1209.
  2. 2. Dellarco V, Henry T, Sayre P, Seed J, Bradbury S (2010) Meeting the common needs of a more effective and efficient testing and assessment paradigm for chemical risk management. J Toxicol Environ Health B Crit Rev 13: 347–60.
  3. 3. Egeghy PP, Judson R, Gangwal S, Mosher S, Smith D, et al. (2012) The exposure data landscape for manufactured chemicals. Sci Total Environ 414: 159–166.
  4. 4. Muir DCG, Howard PH (2006) Are there other persistent organic pollutants? A challenge for environmental chemists. Environ Sci Technol 40: 7157–7166.
  5. 5. Zartarian V, Bahadori T, McKone T (2005) Adoption of an official ISEA glossary. J Expo Anal Environ Epidemiol 15: 1–5.
  6. 6. Judson R, Riachrd A, Dix DJ, Houk K, Martin M, et al. (2009) The toxicity landscape for environmental chemicals. Environ Health Perspect 117: 685–695.
  7. 7. USEPA (2012) TSCA Work Plan Chemicals: Methods Document. Office of Pollution Prevention and Toxics. Washington, DC, USA. http://www.epa.gov/oppt/existingchemicals/pubs/wpmethods.pdf.
  8. 8. Cohen Hubal EA, Richard A, Aylward L, Edwards S, Gallagher J, et al. (2010) Advancing exposure characterization for chemical evaluation and risk assessment. J Toxicol Environ Health B Crit Rev 13: 299–313.
  9. 9. Egeghy PP, Vallero DA, Cohen Hubal (2011) EA (2011) Exposure-based prioritization of chemicals for risk assessment. Environ Sci Pol 14 (8): 950–64.
  10. 10. Mitchell J, Arnot J, Jolliet O, Georgopolous P, Isukapalli S, et al. (2013) Comparison of modeling approaches to prioritize chemicals based on estimates of exposure and exposure potential. Sci Total Environ 458–460: 555–567.
  11. 11. USEPA (2009) A Conceptual Framework for U.S. EPA’s National Exposure Research Laboratory. Office of Research and Development. Washington, DC, USA.
  12. 12. Shatkin J, Qian S (2005) Classification schemes for priority setting and decision making. In: Linkov I, Ramadan AB. Comparative risk assessment and environmental decision making. Dordrecht: Kluwer Academic Publishers. 213–243.
  13. 13. Linkov I, Moberg E (2011) Multi-criteria decision analysis: environmental applications and case studies. Boca Raton: CRC Press.
  14. 14. Canis L, Linkov I, Seager TP (2010) Application of stochastic multiattribute analysis to assessment of single walled carbon nanotube synthesis processes. Environ Sci Technol 44: 8704–8711.
  15. 15. Linkov I, Bates ME, Canis LJ, Seager TP, Keisler JM (2011) A decision-directed approach to prioritizing research into the impact of nanomaterials on the environment and human health. Nat Nanotechnol 6: 784–787.
  16. 16. Yokota F, Gray G, Hammitt JK, Thompson KM (2004) Tiered chemical testing: s value of information approach. Risk Anal 24 (6): 1625–1639.
  17. 17. Gabbert S, Weikard H-P (2010) A theory of chemicals regulation and testing. Nat Resour Forum 34: 155–164.
  18. 18. Paralikas AN, Lygeros AI (2005) A multi-criteria and fuzzy logic based methodology for the relative ranking of the fire hazard of chemical substances and installations. Process Saf Environ Prot 83 (2): 122–134.
  19. 19. Vermeire T, van de Bovenkamp M, Bruinen de Bruin Y, Delmaar C, van Engelen J, et al. (2010) Exposure-based waiving under REACH. Regul Toxicol Pharmacol 58: 408–420.
  20. 20. von Winterfeldt, D, Edwards W (1986) Decision analysis and behavioral research. Cambridge: Cambridge University Press.
  21. 21. Belton V, Stewart T (2002) Multiple criteria decision analysis: an integrated approach. Boston: Kluwer Academic Publishers.
  22. 22. Stillwell WG, Seaver DA, Edwards W (1981) A comparison of weight approximation techniques in multi-attribute utility decision making. Organ Behav Hum Perform 28 (1): 62–77.
  23. 23. USEPA (2011a) Estimation Programs Interface Suite™ for Microsoft® Windows, v 4.10. United States Environmental Protection Agency. Washington, DC, USA. http://www.epa.gov/opptintr/exposure/pubs/episuite.htm. Accessed 2012 Feb 14.
  24. 24. USEPA (2011b) Design for the Environment Program Alternative Assessment Criteria for Hazard Evaluation, v 2.0. Office of Pollution Prevention & Toxics, United States Environmental Protection Agency. Washington DC, USA. http://www.epa.gov/dfe/alternatives_assessment_criteria_for_hazard_eval.pdf.
  25. 25. Rossi M, Heine L (2009) The Green Screen for Safer Chemicals, v 1.0. North American Clean Production Action. Spring Brook, NY, USA. http://www.cleanproduction.org/library/cpa-fact%20grscreen_Jan09_final.pdf.
  26. 26. Dow Corning Corporation (1999) Polydimethylsiloxanes Do Not Bioaccumulate. Health Environment & Regulatory Affairs (HERA) Ref. no. 01-1128-01. Midland, MI, USA. http://www.dowcorning.com/content/publishedlit/01-1128-01.pdf.
  27. 27. Muller M, Nendza M (2007) Literature Study: Effects of Molecular Size and Lipid Solubility on Bioaccumulation Potential. Fraunhofer Institut Molekularbiologie and Angewandte Oekologie, FKZ 360 01 043. Dessau, Germany.http://www.reach-info.de/dokumente/bcf.pdf.
  28. 28. OECD. (1992) Test No. 301: Ready Biodegradability, OECD Guidelines for the Testing of Chemicals, Section 3, OECD Publishing. https://doi.org/10.1787/9789264070349-enhttp://search.oecd.org/officialdocuments/displaydocumentpdf/?cote=ocde/gd(95)43&doclanguage=en
  29. 29. USEPA (1995) Premanufacture Notification Exemption; Revision of Exemption for Chemical Substances Manufactured in Small Quantities; Low Release and Exposure Exemption; Final Rule. (60 FR 16336–16351).U.S. Environmental Protection Agency, Office of Pollution Prevention and Toxics. Washington, DC, USA. http://www.smallbiz-enviroweb.org/Resources/sbopubs/edocs/e04.pdf.
  30. 30. Scheringer M, Jones KC, Matthies M, Simonich S, van de Meent D (2009) Multimedia partitioning, overall persistence, and long-range transport potential in the context of POPs and PBT chemical sssessments. Integr Environ Assess Manag 5 (4): 557–76.
  31. 31. Schrödinger LLC (2007) QikProp 3.0 User Manual. http://www.schrodinger.com/supportdocs/18/17.
  32. 32. National Fire Protection Agency (NFPA) (2007) NFPA 704: Standard System for the Identification of the Hazards of Materials for Emergency Response, 2007 Edition.
  33. 33. Environmental Science Center (2007) PBT Profiler. U.S. Environmental Protection Agency, Environmental Science Center, Office of Chemical Safety and Pollution Prevention. Washington, DC, USA. <http://www.pbtprofiler.net/> Accessed 2013 Mar 18.
  34. 34. USEPA (2011c) Exposure-Based Chemical Prioritization Workshop II: Extending Capabilities for High Throughput Assessment. http://www.epa.gov/ncct/expocast/exposure_based_2.html. Accessed 2013 Mar 18.
  35. 35. Cox Jr (2008) LA (2008) What’s Wrong with Risk Matrices? Risk Anal 28 (2): 497–512.