1 Introduction

Reducing chemical pressure on human and ecological health is an integral part of the global sustainability agenda. This includes the reduction of chemical emissions into the natural environment along product life cycles as well as the reduction of human exposure to chemicals used in consumer goods, as laid out in the United Nations (UN) Sustainable Development Goals (UN 2020) and in the UN Environment Programme’s Strategic Approach to International Chemicals Management (UNEP 2015). Harmonized guidelines are required to consistently quantify life cycle toxicity impacts from chemical emissions as well as from exposure to chemicals in product formulations, materials or articles (hereafter referred to as “products”). Such guidelines should focus especially on providing recommendations for globally applicable and life-cycle–based indicators and underlying methods that are most suitable for the quantitative characterization of human and ecological toxicity impacts associated with chemical emissions and exposure. It is thus important that related recommendations are consistent with the boundary conditions of characterizing toxicity impacts in life cycle impact assessment (LCIA) (Fantke et al. 2018a), while being aligned with chemical-management recommendations from regulatory entities (Saouter et al. 2017a,b).

In response to these needs, the Life Cycle Initiative hosted at the UN Environment Programme, developed and endorsed the scientific consensus model USEtox (Rosenbaum et al. 2008; Westh et al. 2015). Current practice for characterizing human toxicity and freshwater ecotoxicity impacts is implemented in this model. It includes inhalation and ingestion exposure and related health effects from emissions into far-field compartments (air, water, soil) or into a generic indoor compartment. However, despite reflecting—as a scientific consensus model—mature science (Hauschild et al. 2008), the original toxicity characterization framework has limitations, calling for further improvements based on scientific progress. Such improvements are mainly related to increasing the spatiotemporal and population-level resolution of impact estimates and extending the coverage and quality of substance, exposure, and dose-response information (Fantke et al. 2018a,b; Kirchhübel & Fantke 2019; Crenna et al. 2020; Gentil et al. 2020; Holmquist et al. 2020). Furthermore, a series of recent studies has demonstrated that environmentally mediated exposure from chemical emissions is less important for overall exposure than consumer exposure to chemical constituents in products (Shin et al. 2015; Ernstoff et al. 2016; Fantke & Jolliet 2016; Csiszar et al. 2017; Ring et al. 2019; Fantke et al. 2020b; Jolliet et al. 2021). Hence, including pathways related to chemicals in consumer products into toxicity characterization frameworks is crucial for considering all relevant pathways. Another specific need is the improvement of dose-response assessment for non-cancer toxicity that would benefit from the latest progress in stochastic dose-response modeling (WHO 2014, 2017, 2017; Chiu & Slob 2015; Chiu et al. 2018). To address these issues, the Life Cycle Initiative has established taskforces for advancing the original consensus for human toxicity and ecotoxicity characterization, in the frame of a flagship project on global guidance on environmental LCIA indicators (Jolliet et al. 2014, 2018; Frischknecht et al. 2016; Verones et al. 2017). The work conducted by these taskforces has led to a set of recommendations for refining and expanding the original toxicity characterization framework, and for implementing these recommendations into USEtox. In the present paper, we detail the methodological aspects underlying these recommendations and their implementation. We present recommendations for environmental fate, human exposure, and toxicity effects, with focus on including exposure to chemicals in consumer products and non-cancer dose-response modeling, while recommendations for ecotoxicity effects are presented elsewhere. The proposed updated USEtox near-field/far-field exposure and human toxicity modeling framework is evaluated in an illustrative rice production and consumption case study.

2 Consensus-building process

In 2016, the Life Cycle Initiative convened a Human Toxicity taskforce, consisting of leading experts from academia, industry, and government public-health institutions, to expand its guidance for evaluating human toxicity impacts from exposure to chemical substances. The taskforce followed an earlier effort for providing guidance on assessing impacts from exposure to fine particulate matter (Fantke et al. 2015, 2017, 2019b; Hodas et al. 2016). An initial scoping phase identified the need for advancing two major aspects, namely integrating near-field and far-field exposure, and improving human dose-response modeling. To address these aspects, scientific advances were discussed and initial recommendations formulated through three expert workshops held in 2016 and 2017 (Fantke et al. 2018b). Further refinement of these initial recommendations, along with alignment with efforts to update ecotoxicity impact assessment (Fantke et al. 2018a; Owsianiak et al. 2019), led to a set of final recommendations discussed and agreed upon in a scientific expert workshop held in Valencia (Spain) in 2018. While the general recommendations have been summarized in the related workshop report (Frischknecht & Jolliet 2019), the present paper details the underlying methodological framework for human near-field exposure and toxicity characterization.

3 Proposed methodological framework

In USEtox, human toxicity impacts are characterized in terms of population-level lifetime loss associated with cumulative health risk for different health effects (aggregated into cancer and non-cancer effects, each having different average effect severity). Such health effects are associated with exposure to chemicals that distribute among far-field compartments (e.g., outdoor air), finally reaching humans via different exposure pathways (e.g., ingestion of food) (Rosenbaum et al. 2011). Combining these aspects for environmental emissions in a matrix-based framework yields characterization factors expressed in disability-adjusted life years (DALY) lost per unit mass emitted into the environment.

To facilitate the integration of near-field (consumer and indoor) exposure environments into the existing framework of far-field (emission-related outdoor) exposure environments, we reviewed available approaches that model near-field fate and exposure pathways for chemicals in consumer products (Huang et al. 2017). The integrated framework is needed to address human exposure pathways during and after product use on a consistent mass-balance basis. Following recent recommendations, we build on the original rate constant–based USEtox framework to determine transfer fractions. We extend these recommendations to the near-field environment to obtain by matrix inversion cumulative chemical mass transfer fractions and product–related chemical intake fractions. The latter combine near-field consumer household environments with far-field environments into a metric that incorporates the interactions of humans with both types of environments via dermal, inhalation, and oral exposure, and potential feedback via, e.g., exhalation (Jolliet et al. 2015; Fantke et al. 2016). This allows linking human intake via all exposure routes directly to chemical mass in products. The overall structure of the new recommended integrated near-field/far-field human exposure and toxicity characterization framework is illustrated in Fig. 1.

Fig. 1
figure 1

Recommended human toxicity characterization framework. Emission-based far-field exposure pathways are integrated with product-related near-field consumer-exposure pathways, and linked to human toxicity dose-response and effect severity information. Unit mass emitted along a product life cycle or contained in a product can be expressed as unit inventory mass

In line with the original USEtox matrix framework, we arrange in a single matrix \({\mathbf{TF}}^{\mathbf{cum}} \in {\mathbb{R}}^{c\times c}\) all cumulative chemical mass transfer fractions. We start from products and emissions in the different environments along with intake fractions (relating intake to emissions) and product intake fractions (relating intake to mass in products) (Fantke et al. 2016). We obtain a matrix \(\mathbf{CF} \in {\mathbb{R}}^{e\times c}\) of characterization factors per kg inventory mass [DALY/kginventory mass], relating impacts on humans due to health effects \(e\) per unit mass of a chemical in emission-related or product-related compartment of entry \(c\):

$$\mathbf{CF}=\mathbf{SF} \: \mathbf{DRF} \: {\mathbf{TF}}^{\mathbf{cum}}=\mathbf{EF} \: {\mathbf{TF}}^{\mathbf{cum}}$$
(1)

where the diagonal matrix \(\mathbf{SF} \in {\mathbb{R}}^{e\times e}\) of severity factors [DALY/incidence] multiplies matrix \(\mathbf{DRF} \in {\mathbb{R}}^{e\times r}\) of dose-response factors [incidence risk/kgintake] associated with exposure routes \(r\) (i.e., inhalation, ingestion, and dermal exposure, aggregated over all contributing exposure pathways \(x\in r\), such as ingestion of drinking water and ingestion of food). Matrices \(\mathbf{SF}\) and \(\mathbf{DRF}\) are combined into matrix \(\mathbf{EF}\) of toxicity-related health effect factors [DALY/kgintake], multiplying matrix \({\mathbf{TF}}^{\mathbf{cum}}\). The dimensions of the matrix of dose-response factors are thereby matched to \({\mathbf{TF}}^{\mathbf{cum}}\), i.e., \(\mathbf{DRF} \in {\mathbb{R}}^{e\times c}\), with zeroes for all non-human compartments. From this set of matrices, we derive below-described recommendations for characterizing human toxicity and ecotoxicity impacts for chemical emissions and for chemicals used in products.

Impact scores \({IS}^{\text{H}}\) [DALY/functional unit] representing toxicity-related damages on human health are derived from the chemical inventory mass per functional unit, \({m}_{c}\) [kginventory/functional unit]. Chemicals are either emitted to an environmental or available in a product compartment \(c\), multiplied by the sum of corresponding characterization factors for human toxicity (\(\mathrm{H}\)) impacts, \({CF}_{c,r,e}^{\text{H}}\) [DALY/kginventory], across considered exposure routes \(r\) and effect types \(e\), aggregated over all chemical source compartments \(c\):

$${IS}^{\text{H}}={\sum }_{c}\left({m}_{c}\times {\sum }_{r,e}{CF}_{c,r,e}^{\text{H}}\right)$$
(2)

where the chemical inventory mass corresponds directly to the emitted mass for environmental emissions, and to the mass of chemical in a given product application for chemicals in consumer products. This chemical mass in products is derived from the chemical weight fraction in a product application [kg/kgproduct] and the related product mass \({M}_{c}\) [kgproduct]:

$${m}_{c}={wf}_{c}\times {M}_{c}$$
(3)

Characterization factors \({CF}_{c,r,e}^{\text{H}}\) are derived from human intake–related cumulative chemical mass transfer fractions, \({TF}_{c,u,x\in r}^{\text{cum}}\) [kgintake/kginventory], summed over population groups \(u\) and over exposure pathways \(x\) belonging to the same route \(r\). These cumulative transfer fractions are multiplied by toxicity-related health effect factors, \({EF}_{x,e}^{\text{H}}\) [DALY/kgintake], calculated as the product of the dose-response factor, \({DRF}_{r,e}^{\text{H}}\) [incidence risk/kgintake], for exposure route \(r\) and effect type \(e\), and the effect type–specific severity factor, \({SF}_{e}^{\text{H}}\) [DALY/incidence]:

$${CF}_{c,r,e}^\text{H}=\sum\nolimits_{u,x}{TF}_{c,u,x\in r}^\text{cum}\times\overbrace{{DRF}_{r,e}^\text{H}\times{SF}_e^\text{H}}^{{}^{{EF}_{x,e}^\text{H}}}$$
(4)

Human toxicity indicators are most commonly extrapolated from animals to humans (as human data are often unavailable). They express dose-response factors for carcinogenic and non-carcinogenic toxicity effect endpoints that represent a change in human population response per unit change of chemical exposure. To allow for considering the various health endpoints obtained from (human and animal) toxicity studies, results are generally aggregated into cancer and non-cancer effects. Based on the latest work in dose-response modeling (WHO 2014; Chiu & Slob 2015; Chiu et al. 2018), we refine the approaches for the selection of toxicity input data and their extrapolation to derive effect factors for non-cancer effects. Severity factors translate an estimated human response into potential lifetime loss due to death or disease, expressed in DALY. Statistically averaged severity factors for cancer and non-cancer health endpoints were originally based on incidence-weighted DALY from the Global Burden of Disease (GBD) study series for 1990 (Huijbregts et al. 2005), which we update specifically associated with non-cancer responses.

Chemicals emitted into a given environmental compartment or used in a given product application can also reach ecosystems (e.g., via landfilling the product after use or via product use–related emissions into indoor air and subsequent transfer to the outdoor environment). Hence, we also characterize ecotoxicity impacts associated with either the chemical mass emitted to a compartment or the chemical mass in a product. Impact scores \({IS}_{s}^{\text{E}}\) [PDF m3 d/functional unit] for a given ecosystem type \(s\) (e.g., freshwater ecosystems) represent ecotoxicity-related damages on ecosystem quality. They are derived as chemical inventory mass, \({m}_{c}\) [kginventory/functional unit], multiplied by the corresponding characterization factors for ecotoxicity (\(\mathrm{E}\)) impacts, \({CF}_{s,c}^{\text{E}}\) [PDF m3 d/kginventory], aggregated over all chemical source compartments \(c\):

$${IS}_{s}^{\text{E}}={\sum }_{c}\left({m}_{c}\times {CF}_{s,c}^{\text{E}}\right)$$
(5)

Characterization factors \({CF}_{s,c}^{\text{E}}\) are derived from cumulative chemical mass transfer fractions from emission compartment \(c\) to the receiving exposure compartment associated with ecosystem type \(s\), \({TF}_{c,s}^{\text{cum}}\) [kgto compartment/kginventory]. \({TF}_{c,s}^{\text{cum}}\) are multiplied by the first-order residence time in the receiving compartment, which is the inverse of the direct overall removal rate constant of chemical mass from that compartment, \({k}_{s}^{\text{loss}}\) [(kg/d)/kg]. These are further multiplied by the ecosystem exposure factor, \({XF}_{s}\) [kgbioavailable/kgto compartment], representing the bioavailable mass fraction in the compartment of exposed ecosystem type \(s\), and by the ecosystem type–specific, ecotoxicological effect factor, \({EF}_{s}^{\text{E}}\) [PDF m3/kgbioavailable]. \({EF}_{s}^{\text{E}}\) are calculated as the product of the concentration–response factor, \({CRF}_{s}\) [PAF m3/kgbioavailable], and the severity factor, \({SF}_{s}^{\text{E}}\) [PDF/PAF]:

$${CF}_{s,c}^\text{E}={TF}_{c,s}^\text{cum}\times\frac1{k_s^\text{loss}}\times{XF}_s\times\overbrace{{CRF}_s\times{SF}_s^\text{E}}^{{EF}_s^\text{E}}$$
(6)

Since further multimedia transfers are already considered in \({TF}_{c,s}^{\text{cum}}\), we do not use the diagonal of the fate factor matrix as in earlier versions of USEtox, but the first-order chemical residence time (inverse of \({k}_{s}^{\text{loss}}\)), to avoid double counting of feedback into the original emission or product compartment (Margni et al. 2004).

4 Integrated near-field/far-field fate and exposure assessment

4.1 Extended matrix framework

For consistency with the original USEtox far-field model, we structure the product compartment as an entry point for chemicals directly emitted to this compartment or for chemicals in products. We link near-field and far-field emission compartments and human receptor compartments in a single matrix of cumulative chemical mass transfer fractions, \({\mathbf{TF}}^{\mathbf{cum}}\) (see Eq. 1), as illustrated in Fig. 2. Matrix \({\mathbf{TF}}^{\mathbf{cum}}\) is derived from inverting the matrix \(\mathbf{TF} \in {\mathbb{R}}^{c\times c}\) of direct transfer fractions, where each element represents a direct chemical mass transfer fraction that is defined between two adjacent or nested source and receiving compartments (Fantke et al. 2016). Matrix columns representing human receptor compartments are only relevant for feedback of chemical mass fractions taken in by humans that are released back to the environment, such as from exhaling a fraction of an inhaled chemical back into indoor air. Blue highlighted cumulative transfer fractions in Fig. 2 divided by feedback factors are used for obtaining the matrix of fate factors as input for modeling ecotoxicity impacts, consistent with the original USEtox approach. Exposure-related matrix elements in Fig. 2 are highlighted in orange. For all rows referring to human receptor compartments \(i\), the first (i.e., product or direct emission source-related) column \(j\) of \({\mathbf{TF}}^{\mathbf{cum}}\) represents the vector \(\overrightarrow{\text{PiF}}={\mathbf{TF}}^{\mathbf{cum}}[{\alpha }_{i},{\beta }_{j=1}]\) of product intake fractions [kgintake/kginventory mass]. \(\overrightarrow{\text{PiF}}\) elements represent cumulative chemical mass transfer from a given product or direct emission compartment to humans, and all other (i.e. emission source-related) columns represent the sub-matrix \(\mathbf{iF}={\mathbf{TF}}^{\mathbf{cum}}[{\alpha }_{i},{\beta }_{j\ne 1}]\) of intake fractions [kgintake/kgemitted], whose elements describe cumulative chemical transfer from different emission compartments to humans. Each \(\overrightarrow{\text{PiF}}\) element is the sum of direct and indirect transfer fractions contributing to an exposure pathway from a compartment of entry to a specific human receptor compartment. Different scenarios can be compared across exposure pathways for different compartments of entry (i.e., for chemicals in different product types), since all other \({\mathbf{TF}}^{\mathbf{cum}}\) columns are product-independent. The compartment of entry can refer to any near- or far-field emission compartment (for direct emissions) or to any product representing a semi-permeable boundary between product and near-field environment, which chemicals enter when installing products or bringing them into the near-field environment.

Fig. 2
figure 2

Representation of the matrix of cumulative chemical mass transfer fractions. It combines a product-related compartment of entry (left-most column) with near-field, far-field, and human receptor compartments. Product intake fractions (PiF) start from the product compartment of entry and intake fractions (iF) start from environmental emission compartments, both linking to pathway- and population-specific human exposure to chemicals

4.2 Direct chemical mass transfer fractions

The first column in the matrix \(\mathbf{T}\mathbf{F}\) of direct chemical mass transfer fractions is specified to be a product type–specific compartment of entry, whose entries vary as function of transfer pathways relevant for different product types or emission scenarios. To calculate transfer fractions for chemicals in an initial emission-related or product-related compartment of entry for a variety of exposure pathways, we incorporated five main underlying models for “direct emission,” and for “article interior,” “skin surface layer,” “object surface,” and “food contact material” products. Table 1 summarizes the direct transfer fractions that are determined by each underlying model.

Table 1 Underlying models for calculating direct chemical mass transfer fractions, \(TF\) [ −], from chemicals in products, including exposure pathways and example products covered

The direct emission model contains a set of fixed direct transfer fractions to environmental compartments. For emissions, direct transfer fractions are derived from original USEtox rate constants describing chemical transfers or losses from compartments representing outdoor environments (Rosenbaum et al. 2008), indoor environments (Rosenbaum et al. 2015), and waste treatment (e.g., Kjeldsen & Christensen 2001). Related direct human exposure fractions include inhalation of indoor air, ingestion of drinking water, and dermal uptake. The article interior model considers direct transfer fractions from articles (e.g., building materials, furniture, toys) to indoor air, to skin via dermal contact, and to the gastrointestinal tract via dust ingestion. Depending on chemical-material properties, different models are used for transfer to indoor air, namely, diffusion-limited models (accounting for chemical diffusion within articles via Fick’s 2nd Law) (Huang & Jolliet 2016) and partition-limited models (assuming that chemicals are evenly distributed inside articles and accounting for indoor sorption) (see Electronic Supplementary Material, ESM, Section S-2a). Direct transfers for dermal contact and dust ingestion assume that chemical concentrations on article surfaces are in equilibrium with those in skin surface lipids and in dust, respectively. The skin surface layer model is applicable to evaluate chemicals in personal care and other products (e.g., hand sanitizers) applied directly on human skin. This model uses a three-compartment mass-balance model that includes skin, indoor air, and product applied, and assumes that volatilization and skin permeation are competing loss processes (Ernstoff et al. 2016). The object surface model is used for chemicals applied on object surfaces (e.g., surface cleaners). This model uses a four-compartment mass-balance approach that includes near-person and far-person surfaces (in the same indoor environment), air, and the applied product (Wang et al. 2016), with transfers between near-person and far-person surfaces to simulate the movement of a person when cleaning surfaces. The food contact material model estimates the fraction of organic chemicals migrating from polymeric packaging materials into food matrices as a function of diffusion coefficients within packaging materials and packaging-food partition coefficients (Ernstoff et al. 2017). Indoor sorption is considered across product-specific models and detailed in ESM, Section S-2d.

5 Human toxicity dose-response framework

We updated the approach for deriving human toxicity dose-response factors (see Eq. 1) for non-cancer endpoints. The updated approach is illustrated in Fig. 3 and consists of two main steps for each chemical: (1) identifying a point of departure (POD) using a decision-tree based on data availability and reliability and (2) deriving human lifetime equivalent dose-response factors from the selected PODs using a probabilistic framework. This framework addresses a number of key limitations of the original approach. First, it significantly broadens the coverage of chemicals for which dose-response factors can be derived by utilizing a wider range of toxicity data sources beyond regulatory human-health-hazard assessments. Second, it incorporates recent guidance from the WHO (2014, 2017) on more explicitly considering both uncertainty and human variability. Although dose-response factors are still fixed slopes based on point estimates, the values themselves are derived from a statistical distribution rather than being the result of deterministic calculations. Finally, the value of the dose-response factor was chosen to better represent the marginal slope of the dose=-response relationship in a dose range that is environmentally relevant.

Fig. 3
figure 3

Recommended approach for deriving human toxicity dose-response information based on a hierarchy of different toxicity input data to identify a suitable individual point of departure (POD), and a probabilistic extrapolation toward human population lifetime toxicity dose-response factors

5.1 Identifying PODs

As shown in Fig. 3 (top panel), our updated framework substantially expands the acceptable data sources for points of departure (PODs). Animal or human toxicity values (e.g., benchmark dose lower confidence limit (BMDL), no-observed adverse effect level (NOAEL), or lowest observed adverse effect level (LOAEL)) used in regulatory human health hazard assessments are the preferred source. If such data are not available, other experimental animal data, new approach methods (NAMs), such as in vitro and in silico models, or Thresholds of Toxicological Concern (TTC) are used to derive an individual animal or human POD, as detailed in ESM, Section S-3a. However, even with regulatory assessments, it is often the case that different countries or organizations come to different conclusions and select different PODs, rendering it necessary to have an explicit hierarchy of sources. As an example, for USA-based regulatory values, we apply the hierarchy developed by Wignall et al. (2014), (2018) (see ESM, Section S-3b).

5.2 Deriving human lifetime equivalent dose-response factors

Once PODs are identified along with accompanying meta-data, we use the automated workflow developed by Chiu et al. (2018) to implement the WHO (2014, 2017) probabilistic framework and derive population-level dose-response functions. This workflow involves the following steps for each POD. First, we assign a “conceptual model” based on the type of effect (continuous or dichotomous, reflecting a stochastic process or not) to each animal or human POD, using criteria listed in ESM, Table S6. Next, uncertainty distributions are generated for each POD using the default distributions derived by WHO (2014, 2017) and informed by the POD meta-data, as outlined in ESM, Table S7. With that, a median individual benchmark dose (BMD) is obtained from POD-specific extrapolations:

$$\mathrm{BMD}= \frac{\mathrm{POD}}{{ef}_{\text{BMD}}}$$
(7)

where \({ef}_{\text{BMD}}\) is a probabilistic factor for NOAEL, LOAEL, or BMDL to BMD extrapolation (for POD = BMD, the central estimate of \({ef}_{\text{BMD}}=1\)). From the BMD, we derive a human-equivalent individual dose-response level for a given effect magnitude (e.g., % change in blood pressure) expressed in mg/(kgBW/d) for oral exposure and mg/m3 for inhalation exposure. We define this dose-response level as \({\text{ED5}}{0}\), which is the effect dose for the median (i.e., 50th %-ile) individual of the human population, thus inducing a 50% human population response:

$${\text{ED50}}=\frac{\mathrm{BMD}}{{ef}_{\text{BW}}\times {ef}_{\text{TKTD}}\times {ef}_{t}}$$
(8)

where \({ef}_{\text{BW}}\), \({ef}_{{\text{TKT}}{\text{D}}}\), and \({ef}_{t}\) are probabilistic factors for, respectively, interspecies body weight (BW) scaling, interspecies toxicokinetic (TK)/toxicodynamic (TD) differences, and exposure duration extrapolation (ESM, Table S7). In case of human-based POD, central estimates for \({ef}_{\text{BW}}={ef}_{\text{TKTD}}=1\), but their uncertainty remains. These distributions are applied using the approach illustrated in Fig. 3 (bottom panel), based on the WHO model. In this model, the human-population incidence-response level \(I\) (e.g., 10%) can be expressed as a cumulative lognormal distribution reflecting intra-species variability around the exposure dose level \(X\), namely, \(I(X)=\Phi\lbrack z(X)\rbrack,z(X)=(\text{ln}X-\text{ln}\text{ED5}0)/\sigma_\text{H}\). In this notation, \({\text{ED5}}{0}\) is defined as in Eq. 8, \({\sigma }_{\text{H}}\) is the log-scale standard deviation that reflects human variability in sensitivity, and the dose level \(X={\text{ED}}{10}\) corresponds to the population incidence response level \(I=10\%\), which is then used to derive a human-equivalent population lifetime dose by following the steps described below. This model is implemented, using the approach proposed by Chiu and Slob (2015), dividing the \({\text{ED5}}{0}\) by a probabilistic factor \({ef}_{{\text{H}},I}\), calculated for an incidence level of \(I=10\%\):

$${\text{ED}}{10}= \frac{{\text{ED}}50}{{ef}_{{\text{H}},I=0.1}}$$
(9)

Specifically, \({ef}_{{\text{H}},I=0.1}=\mathrm{exp}(1.282\times {\sigma }_{\text{H}})\), where 1.282 corresponds to the absolute value of the \(z\)-score for a 10% human-population response level. Not only the POD but also all extrapolation factors in Eqs. 78 are assigned distributions, from which \({\text{ED}}{10}\) is derived probabilistically.

Next, to derive an effect metric compatible with USEtox, the central tendency (median estimate) of the \({\text{ED}}{10}\) is converted into a human-equivalent lifetime dose, \({\text{DL}}{\text{T}}{10}\) (kg/lifetime):

$${\text{DLT}}{10}= \frac{{\text{ED}}{10}\times {F}_{\left\{BW,BR\right\}}\times LT\times {cf}_{\text{d/yr}}}{{cf}_{\text{mg/kg}}}$$
(10)

where for oral exposure, \({F}_{\left\{BW,BR\right\}}=70\) kg is the average human body weight and for inhalation exposure, \({F}_{\left\{BW,BR\right\}}=16\) m3/d is the average human individual breathing rate; \(LT=70\) years is the average human lifetime, and \({cf}_{\text{d/yr}}=365.25\) d/yr and \({cf}_{\text{mg/kg}}={10}^{6}\) mg/kg are conversion factors for, respectively, days per year and mg per kg. Finally, the human-equivalent lifetime dose-response factor, \(DRF\) (incidence risk/kg), is determined as linear slope to \({\text{DLT}}{10}\):

$$DRF= \frac{0.1}{{\text{DLT}}{10}}$$
(11)

This central-tendency linearly extrapolated slope from \({\text{ED}}{10}\) is also approximately equal to that of the marginal slope at 1% human-population response level as illustrated in Fig. 3, bottom panel (Fantke et al. 2019a), a response level that corresponds to a dose range that is environmentally relevant, especially when considering near-field exposure to chemicals in products.

The derived \(DRF\) is meant to represent the incremental (or marginal) increase in cumulative population risk for each unit of additional exposure, which is the slope of the dose-response relationship. The original LCIA approach uses the line between the point at 50% response and zero to represent this slope—using the above notation, this slope is simply \(0.5/{\text{ED}}{50}\). A key limitation of the previous approach is that it does not incorporate any information about population variability, which can be substantial, and can vary widely across chemicals (Zeise et al. 2013; Chiu et al. 2018). Indeed, the true slope of the dose-response relationship is non-linear, with a value that varies continuously depending on dose.

Ideally, each study would establish a “working point” exposure at which this slope is evaluated and the distribution to calculate the \(DRF\) in each case. However, there are several difficulties with implementing such a case-specific \(DRF\). First, due to lack of exposure data, it is very difficult to establish the current “background” exposure level for any given compound and spatially varying emission, even at the population level, let alone its population distribution. Second, as pointed out by Zeise et al. (2013), chemicals are additive to background disease, rather than background exposure to a specific compound. Indeed, Huijbregts et al. (2005) and Fantke et al. (2018b) have suggested that the “background exposure” could be set as an “effective” exposure at the equivalent background incidence rate. In sum, there is substantial uncertainty in the appropriate exposure dose at which to evaluate the slope of the dose-response relationship. Given these uncertainties, the linear slope from the 10% response level was found to be adequately representative of a reasonable range of possible marginal slopes at 1% response as further detailed in ESM, Section S-3d.

6 Human health effects severity

To allow comparisons of human toxicity impacts with other impacts contributing to damages on human health, levels of severity are assigned to the predicted population response levels, yielding DALY estimates associated with incremental exposure. Huijbregts et al. (2005) provided estimates of the incidence-weighted-average DALY associated with a range of both cancer and non-cancer health endpoints. They proposed that impacts could be assessed using these average DALY values (albeit with high uncertainty, particularly for non-cancer endpoints), even though responses estimated from animal toxicity data can rarely be mapped to specific human diseases.

We propose that the previous approach be refined to address at least one additional subset of non-cancer responses separately from other non-cancer responses. Substances that cause birth defects may be of special interest, because of the clear dichotomous nature of the response, the presence of directly analogous disease states in humans, and the severity and duration of the outcome (US-EPA 1991). This includes a group of disease categories designated in the Global Burden of Disease (GBD) databases as congenital anomalies (birth defects). However, inclusion of these endpoints in the broader non-cancer category can severely underestimate such outcomes. Separation of this category of effects is justifiable from both a mechanistic and a statistical point of view, given the heterogeneity in DALY between this category and other non-cancer effect categories (Hay et al. 2017). The term congenital anomalies as used in the GBD refers, in the context of toxicity data, to “reproductive/developmental toxicity” effects. Developmental outcomes are effects that manifest in the offspring, while reproductive effects are those that affect the fertility or function of a parent for reproduction. This entire category of reproductive/developmental toxicity draws special attention in the regulatory world (along with cancer and mutagenesis). An example is the CMR (carcinogenicity, mutagenicity, reproductive/developmental toxicity) designation used in Europe (EC 2008). Developmental outcomes thereby range in their severity from mild to extreme. However, the important point for all reproductive or developmental effects is that they have the potential to adversely affect humans for their entire lifetime, either because they were never born (effects on reproduction) or because they were born with either functional or morphological deficits (effects on development).

We recommend using incidence-weighted DALY values for all cancer endpoints from Huijbregts et al. (2005), as well as separate values for the reproductive or developmental effects and other non-cancer diseases as presented in ESM, Table S8. Separation of the reproductive/developmental category from other non-cancer responses results in a substantial decrease in the uncertainty associated with the average non-cancer (other than reproductive or developmental effects) DALY values. The reproductive/developmental DALY should be applied for any substance for which the effect factor is derived from reproductive or developmental outcomes. In addition to this separation, we recommend that DALY values per incidence for both non-cancer categories be updated with the most recent GBD statistics (Salomon et al. 2015; Stanaway et al. 2018).

7 Case study on rice production and consumption

The proposed updated USEtox near-field/far-field exposure and human toxicity modeling framework was demonstrated in an illustrative rice production and consumption case study based on three scenarios, using 1 kg of cooked, white rice as functional unit (FU). In the first scenario, rice is produced/processed in rural areas and distributed/cooked in urban China (CN); in the second scenario, rice production-to-cooking takes place in rural India (IN); and in the third scenario, rice is produced/processed in rural USA and distributed/cooked in urban Switzerland (US/CH). Further case study details are provided in Frischknecht et al. (2016). For human toxicity impacts, we considered 115 organic chemicals for which we could determine non-zero emissions along the rice production life cycle for at least one scenario. To specifically evaluate our near-field exposure component, we additionally considered consumer exposure to six chemicals found in rice packaging (see ESM, Table S9), namely a single recycled cardboard package (CN), a single polyethylene rice package (IN), and multiple polyethylene rice cooking bags stored in a single recycled cardboard package (US/CH).

We applied the article interior model to calculate chemical mass released from rice packaging to indoor air and the food contact material model to obtain the chemical mass migrated from rice packaging into rice during storage and boiling (see Table 1 for model descriptions). On the toxicity side, we derived non-cancer \(DRF\) values according to the approach presented in Fig. 3 for 71 chemicals based on regulatory toxicity values, 47 based on in vivo experimental animal toxicity studies, and 3 utilizing QSAR (see ESM, Fig. S5 for ingestion exposure and Fig. S6 for inhalation exposure). Furthermore, we obtained information on carcinogenicity for 58 chemicals. Among those, we derived cancer \(DRF\) values for 37 chemicals, while the data for 21 chemicals showed no carcinogenic effect indication in the underlying experimental studies.

All individual chemicals across inventory, intake, effect, and impact-score levels were compared and ranked by decreasing impact scores for the US/CH scenario (ESM, Fig. S7). Packaging-related manufacturing and disposal emissions are for most chemicals one to three orders of magnitude lower than total cradle-to-gate emissions. Overall, inventory results span over 10 orders of magnitude, with varying contributions of emission compartments across chemicals and scenarios. Intake fractions, aggregated over all contributing exposure pathways, span over 8 orders of magnitude, ranging from \(5\times {10}^{-9}\) to almost 1. For individual pathways, emission-based intake fractions range from \(<{10}^{-15}\), i.e., 1 pg inhaled per kg halosulfuron-methyl, up to \(3\times {10}^{-3}\), i.e., 3 g ingested per kg 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), both being emitted to freshwater, spanning 16 orders of magnitude. Product intake fractions for chemicals in rice packaging range from 2 ppt of laurolactam inhaled by household members after volatilization during storage to > 90% ingested by users after migration to rice during boiling, spanning 12 orders of magnitude across exposed sub-populations and pathways for the same chemical. Human-toxicity effect factors combining dose-response and effect severity span over more than 10 orders of magnitude, from \({10}^{-4}\) DALY per kg propene inhaled to 7 million DALY per kg TCDD ingested. Inventory data, (product) intake fractions, and effect factors for all case study chemicals are provided in ESM, Table S10.

Combining emission-based intake fractions and product intake fractions with effect factors yields characterization factors across rice case study chemicals, aggregated over all pathways per exposure route (Fig. 4a). Substantial characterization factors are found for all six chemicals in rice packaging (ranging from 0.0003 to 0.07 DALY/kg in packaging) and for some of the 115 emitted chemicals (highest being TCDD with 21,000 DALY/kg emitted to freshwater). Characterization factors span 18 orders of magnitude, with varying contributions of effect types (cancer, reproductive/developmental, and other non-cancer effects) and exposure routes (inhalation, ingestion, and dermal exposure) across chemicals.

Fig. 4
figure 4

Human toxicity characterization factors (a) and impact scores (bd) for the rice case study. Characterization factors (right-side y-axis) are plotted as function of human (product) intake fractions combining fate and exposure (x-axis) and corresponding human toxicity-related effect factors combining dose-response and effect severity (left-side y-axis). Impact scores (right-side y-axis) are plotted as function of life cycle inventory mass (x-axis) and characterization factors (left-side y-axis). Industrial chemicals are emitted to continental rural air, freshwater, and natural soil, and pesticides are emitted to agricultural soil. TCDD 2,3,7,8-tetrachlorodibenzo-p-dioxin

Combining characterization factors with chemical mass in packaging and cradle-to-gate emissions yields impact scores expressed as DALY per functional unit (Fig. 4b–d). In these cumulative scores, results from near-field exposure to chemicals in consumer products (in our case rice packaging) are directly comparable with results from far-field population exposure to chemical emissions. Chemicals in packaging all show high impact scores, mostly via ingestion following migration from rice packaging into rice, and to a lesser extent via inhalation after volatilization from packaging into indoor air. This emphasizes the importance of including exposure to chemicals in consumer products and incorporating an adequate indoor environment in toxicity characterization methods. While this primarily illustrates how chemicals in food packaging and near-field exposure are relevant, it would also be interesting in future studies to compare the impacts of packaging, accounting for the entire supply chain of the packaging and a larger number of contained chemicals. Impact scores range from \({10}^{-23}\) DALY for trichlorofluoromethane to \(2\times {10}^{-5}\) DALY for dibutyl phthalate, spanning 18 orders of magnitude. Aggregating impact scores yields between \(5\times {10}^{-8}\) DALY for US/CH and \(2\times {10}^{-7}\) DALY for IN, dominated by pesticide emissions to soil across scenarios. Rice packaging manufacturing and disposal contributes < 1% to overall impact scores across scenarios. Impacts from direct exposure to packaging chemicals exceed cradle-to-gate impacts by up to > 2 orders of magnitude, ranging from \(2\times {10}^{-6}\) DALY for IN to \(3\times {10}^{-5}\) DALY for CN, dominated by exposure of consumers via ingestion of packaging chemicals migrating into rice. When translating the functional unit into daily use, we would arrive at impacts per individual usage per person. With that, our framework can also be applied in chemical exposure assessment and risk prioritization frameworks (Fantke and Illner 2019; Fantke et al. 2020a; Jolliet et al. 2021).

Based on the application of our framework in the rice case study, we developed additional recommendations. First, we recommend expressing human toxicity characterization factors as DALY per unit mass emitted (for environmentally mediated exposure) or as DALY per unit mass in a product application (for product use-related exposure), which can be interpreted as relative capacity to cause harm. We further recommend presenting impact scores separately for organic substances and metals/metalloids (the latter have not been discussed in the present document), to acknowledge important differences in characterizing fate, exposure, and effects for these substance groups. Finally, we recommend presenting human toxicity impact results on log10-scale and allow for no more than two significant digits (e.g., \(3.4\times {10}^{-5}\)) due to the substantial variability across and uncertainty in characterization factors.

8 Conclusions and next steps

Following up on recommendations from an expert task force for refining and expanding the original LCIA toxicity characterization approach, we present the methodological framework underlying these recommendations and describe their implementation in USEtox. We considered recommendations for environmental fate, human exposure, and toxicity effects, with a focus on expanding the exposure assessment to consider consumer exposure to chemicals in products and refining the non-cancer dose-response assessment by adapting a probabilistic framework. These adaptations enable future studies to (a) go beyond life cycle emissions by also considering exposure to chemicals related to various product use settings; (b) obtain inhalation, ingestion, and dermal exposure estimates for various chemical-product combinations, and delineating between product users, other relevant sub-populations, and the rest of the human population; (c) derive dose-response information for a much broader range of chemicals; (d) improve the reliability of human effect factors by applying a non-linear extrapolation from systematically selected points of departure; and (e) discriminate severity for three indicators, namely cancer, developmental/reproductive non-cancer, and other non-cancer effects. In all, the presented methodological advancements of the USEtox toxicity characterization model address several key limitations in earlier approaches. This includes considering both emission-based and product-based exposure along with their respective health effects in a consistent matrix framework. This is supported by our case study results, illustrating that product use–related exposure dominates overall life cycle exposure. Sub-models for product-related exposure already allow for considering dynamics in fate and exposure settings. To consider such dynamics also for environmental far-field emissions, we would need to define an exposure duration and introduce a dynamic solution for the entire multi-compartment chemical mass transfer fractions system.

Applying a single overarching framework for considering various exposure and toxicity effects for a wide range of chemicals found in multiple products and their life cycles in a comparable way is very beneficial. However, it requires some necessary simplifications and assumptions (especially with respect to low-dose exposure and the consideration of differences in exposure and toxicity), which all contribute to a substantial level of uncertainty. In a comparative context with results ranging over many orders of magnitude and with a limited number of chemicals often dominating overall impact profiles, such uncertainties are acceptable while guiding further research efforts for refining our framework.

For impact of food contact substances, useful insights on maximal levels and potential material matrices for intended conditions of use can be found in regulatory databases (e.g., the US-FDA “Inventory of Environmental Impact Decisions for Food Contact Substance Notifications” with more than 2000 substance-material combinations). To study specific types of packaging, foreground data from manufacturing companies might be necessary to determine the specific food contact substances used in packaging, also considering the life cycle of the packaging. In the longer term, LCI databases might add data on typical chemical masses in products for a range of household products, such as food packaging, cosmetics, and paints. However, these are likely to be customized to each specific product and usage type, thus the importance of making the USEtox tool and its near-field component available to practitioners to determine exposure and impacts for the specific product studied.

Table 2 Primary application areas of the recommended near-field/far-field USEtox modeling framework. The presented framework is suitable for comparative evaluations of chemicals emitted along product life cycles and chemicals in various types of product applications

Next steps are to test the recommended characterization framework in various case studies and to close the remaining research gaps related to, for example, missing chemical content in product information, introducing additional exposure pathways and product types, as well as including the latest effect test data and effect severity information. The new toxicity data that have recently been made available also open the possibility to substantially extend the training set for in silico extrapolation and improve the quality of route-of-exposure–specific dose-response extrapolations. Our framework is applicable not only for evaluating chemical emissions and product-related exposure in LCIA but also in chemical alternatives assessment (CAA) and chemical substitution, consumer exposure and risk screening, and high-throughput chemical prioritization (see Textbox 2). The formal USEtox near-field/far-field model is freely available at https://usetox.org.