Bioconcentration, bioaccumulation, biomagnification and trophic magnification: a modelling perspective.

We present a modelling perspective on quantifying metrics of bio-uptake of organic chemicals in fish. The models can be in concentration, partition ratio, rate constant (CKk) format or fugacity, Z and D value (fZD) format that are shown to be exactly equivalent, each having it merits. For most purposes a simple, parameter-parsimonious one compartment steady-state model containing some 13 parameters is adequate for obtaining an appreciation of the uptake equilibria and kinetics for scientific and regulatory purposes. Such a model is first applied to the bioaccumulation of a series of hypothetical, non-biotransforming chemicals with log KOW (octanol-water partition ratio) values of 4 to 8 in 10 g fish ranging in lipid contents to deduce wet-weight and lipid normalized concentrations, bioaccumulation and biomagnification factors. The sensitivity of biomagnification factors to relative lipid contents is discussed. Second, a hypothetical 5 species linear food chain is simulated to evaluate trophic magnification factors (TMFs) showing the critical roles of KOW and biotransformation rate. It is shown that lipid normalization of concentrations is most insightful for less hydrophobic chemicals (log KOW < 5) when bio-uptake is largely controlled by respiratory intake and equilibrium (equi-fugacity) is approached. For more hydrophobic chemicals when dietary uptake kinetics dominate, wet weight concentrations and BMFs are more insightful. Finally, a preferred strategy is proposed to advance the science of bioaccumulation using a combination of well-designed ecosystem monitoring, laboratory determinations and modelling to confirm that the perceived state of the science contained in the models is consistent with observations.

Benthic invertebrate exposure and chronic toxicity risk analysis for cyclic volatile methylsiloxanes: Comparison of hazard quotient and probabilistic risk assessment approaches.

This study utilized probabilistic risk assessment techniques to compare field sediment concentrations of the cyclic volatile methylsiloxane (cVMS) materials octamethylcyclotetrasiloxane (D4, CAS # 556-67-2), decamethylcyclopentasiloxane (D5, CAS # 541-02-6), and dodecamethylcyclohexasiloxane (D6, CAS # 540-97-6) to effect levels for these compounds determined in laboratory chronic toxicity tests with benthic organisms. The concentration data for D4/D5/D6 in sediment were individually sorted and the 95th centile concentrations determined in sediment on an organic carbon (OC) fugacity basis. These concentrations were then compared to interpolated 5th centile benthic sediment no-observed effect concentration (NOEC) fugacity levels, calculated from a distribution of chronic D4/D5/D6 toxicologic assays per OECD guidelines using a variety of standard benthic species. The benthic invertebrate fugacity biota NOEC values were then compared to field-measured invertebrate biota fugacity levels to see if risk assessment evaluations were similar on a field sediment and field biota basis. No overlap was noted for D4 and D5 95th centile sediment and biota fugacity levels and their respective 5th centile benthic organism NOEC values. For D6, there was a small level of overlap at the exposure 95th centile sediment fugacity and the 5th centile benthic organism NOEC fugacity value; the sediment fugacities indicate that a negligible risk (1%) exists for benthic species exposed to D6. In contrast, there was no indication of risk when the field invertebrate exposure 95th centile biota fugacity and the 5th centile benthic organism NOEC fugacity values were compared.

Bioaccumulation and trophic transfer of cyclic volatile methylsiloxanes (cVMS) in the aquatic marine food webs of the Oslofjord, Norway.

The trophic transfer of cyclic methylsiloxanes (cVMS) in aquatic ecosystems is an important criterion for assessing bioaccumulation and ecological risk. Bioaccumulation and trophic transfer of cVMS, specifically octamethylcyclotetrasiloxane (D4), decamethylcyclopentasiloxane (D5), and dodecamethylcyclohexasiloxane (D6), were evaluated for the marine food webs of the Inner and Outer Oslofjord, Norway. The sampled food webs included zooplankton, benthic macroinvertebrates, shellfish, and finfish species. Zooplankton, benthic macroinvertebrates, and shellfish occupied the lowest trophic levels (TL ≈2 to 3); northern shrimp (Pandalus borealis) and Atlantic herring (Clupea harengus) occupied the middle trophic levels (TL ≈3 to 4), and Atlantic cod (Gadus morhua) occupied the highest tropic level (TL>4.0). Trophic dynamics in the Oslofjord were best described as a compressed food web defined by demersal and pelagic components that were confounded by a diversity in prey organisms and feeding relationships. Lipid-normalized concentrations of D4, D5, and D6 were greatest in the lowest trophic levels and significantly decreased up the food web, with the lowest concentrations being observed in the highest trophic level species. Trophic magnification factors (TMF) for D4, D5, and D6 were <1.0 (range 0.3 to 0.9) and were consistent between the Inner and Outer Oslofjord, indicating that exposure did not impact TMF across the marine food web. There was no evidence to suggest biomagnification of cVMS in the Oslofjord. Rather, results indicated that trophic dilution of cVMS, not trophic magnification, occurred across the sampled food webs.

Roles of steady-state and dynamic models for regulation of hydrophobic chemicals in aquatic systems: A case study of decamethylcyclopentasiloxane (D5) and PCB-180 in three diverse ecosystems.

We seek to contribute to the improved regulatory use of mass balance models to complement environmental monitoring data by applying the steady-state Quantitative Water Air Sediment Interactive model (QWASI) and a novel unsteady-state QWASI model. A steady-state model can yield not only a useful simulation of chemical fate under near steady-state conditions, but it can provide insights into the likely influences of increasing or decreasing emission rates, temperature changes, and unexpectedly high sensitivities to model parameters that may require additional investigation. We compared the consistency of insights from both types of model, in the expectation that while the dynamic model provides a closer simulation of actual conditions, for many purposes a simple, less computationally demanding, more transparent and less expensive model may be adequate for many regulatory purposes. We investigated the response times of decamethylcyclopentasiloxane (D5) and PCB-180 concentrations in water and sediment under three emission scenarios in three different aquatic systems, namely Lake Ontario, Oslofjord, and Lake Pepin. D5 was predicted to be removed largely by hydrolysis and volatilization in Lake Ontario and Oslofjord whereas it is subject to removal by advective loss in Lake Pepin. The half-times of D5 water concentration to a stepwise reduction in emission were <60 days in all three water bodies. In contrast, the predicted half-times were 0.53, 1.4, and 2.9 years in Lake Pepin, Oslofjord, and Lake Ontario, respectively. We also explored how uncertainties in input parameters propagate into uncertainties of concentrations in water and sediments possibly necessitating more accurate values.

Trophic dilution of cyclic volatile methylsiloxanes (cVMS) in the pelagic marine food web of Tokyo Bay, Japan.

Bioaccumulation and trophic transfer of cyclic volatile methylsiloxanes (cVMS), specifically octamethylcyclotetrasiloxane (D4), decamethylcyclopentasiloxane (D5), and dodecamethylcyclohexasiloxane (D6), were evaluated in the pelagic marine food web of Tokyo Bay, Japan. Polychlorinated biphenyl (PCB) congeners that are "legacy" chemicals known to bioaccumulate in aquatic organisms and biomagnify across aquatic food webs were used as a benchmark chemical (CB-180) to calibrate the sampled food web and as a reference chemical (CB-153) to validate the results. Trophic magnification factors (TMFs) were calculated from slopes of ordinary least-squares (OLS) regression models and slopes of bootstrap regression models, which were used as robust alternatives to the OLS models. Various regression models were developed that incorporated benchmarking to control bias associated with experimental design, food web dynamics, and trophic level structure. There was no evidence from any of the regression models to suggest biomagnification of cVMS in Tokyo Bay. Rather, the regression models indicated that trophic dilution of cVMS, not trophic magnification, occurred across the sampled food web. Comparison of results for Tokyo Bay to results from other studies indicated that bioaccumulation of cVMS was not related to type of food web (pelagic vs demersal), environment (marine vs freshwater), species composition, or location. Rather, results suggested that differences between study areas was likely related to food web dynamics and variable conditions of exposure resulting from non-uniform patterns of organism movement across spatial concentration gradients.

The chemical exposure toxicity space (CETS) model: Displaying exposure time, aqueous and organic concentration, activity, and onset of toxicity.

A 1-compartment toxicokinetic model is used to characterize the chemical exposure toxicity space (CETS), providing a novel graphic tool that can aid in the design of aquatic toxicity tests for fish and for interpreting their results. The graph depicts the solution to the differential equation describing the uptake kinetics of a chemical by a modeled fish under conventional bioassay conditions. The model relates the exposure concentration in the water to a dimensionless time and the onset of toxicity as determined by an estimated or assumed critical body residue or incipient lethal aqueous concentration. These concentration graphs are specific to each chemical and exposure and organism parameters and clearly demonstrate differences in toxicity between chemicals and how factors such as hydrophobicity influence the toxic endpoint. The CETS plots can also be used to assess bioconcentration test conditions to ensure that concentrations are well below toxic levels. Illustrative applications are presented using a recent set of high-quality toxicity data. Conversion of concentrations to chemical activities in the plots enables results for different baseline toxicants to be superimposed. For chemicals that have different modes of toxic action, the increased toxicity then becomes apparent. Implications for design and interpretation of aquatic toxicity tests are discussed. The model, and pictorial visualization of the time-course of aquatic toxicity tests, may contribute to improvements in test design, implementation, and interpretation, and to reduced animal usage. Environ Toxicol Chem 2017;36:1389-1396. © 2016 The Authors. Environmental Toxicology and Chemistry Published by Wiley Periodicals, Inc. on behalf of SETAC.

Processes influencing chemical biomagnification and trophic magnification factors in aquatic ecosystems: Implications for chemical hazard and risk assessment.

Bioconcentration factors (BCFs) and bioaccumulation factors (BAFs) are widely used in scientific and regulatory programs to assess chemical hazards. There is increasing interest in also using biomagnification factors (BMFs) and trophic magnification factors (TMFs) for this purpose, especially for highly hydrophobic substances that may reach high concentrations in predatory species that occupy high trophic level positions in ecosystems. Measurements of TMFs in specific ecosystems can provide invaluable confirmation that biomagnification or biodilution has occurred across food webs, but their use in a regulatory context can be controversial because of uncertainties related to the reliability of measurements and their regulatory interpretation. The objective of this study is to explore some of the recognized uncertainties and dependencies in field BMFs and TMFs. This is accomplished by compiling a set of three simple food web models (pelagic, demersal and combined pelagic-demersal) consisting of up to seven species to simulate field BMFs and TMFs and to explore their dependences on hydrophobicity (expressed as log KOW), rates of biotransformation and growth, sediment-water fugacity ratios, and extent of food web omnivory and issues that arise when chemical concentration gradients exist in aquatic ecosystems. It is shown that empirical TMFs can be highly sensitive to these factors, thus the use of TMFs in a regulatory context must recognize these sensitivities. It is suggested that simple but realistic evaluative food web models could be used to extend BCF and BAF assessments to include BMFs and TMFs, thus providing a tool to address bioaccumulation hazard and the potential risk of exposures to elevated chemical concentrations in organisms at high trophic levels.

Evaluating the roles of biotransformation, spatial concentration differences, organism home range, and field sampling design on trophic magnification factors.

Trophic magnification factors (TMFs) are field-based measurements of the bioaccumulation behavior of chemicals in food-webs. TMFs can provide valuable insights into the bioaccumulation behavior of chemicals. However, bioaccumulation metrics such as TMF may be subject to considerable uncertainty as a consequence of systematic bias and the influence of confounding variables. This study seeks to investigate the role of systematic bias resulting from spatially-variable concentrations in water and sediments and biotransformation rates on the determination of TMF. For this purpose, a multibox food-web bioaccumulation model was developed to account for spatial concentration differences and movement of organisms on chemical concentrations in aquatic biota and TMFs. Model calculated and reported field TMFs showed good agreement for persistent polychlorinated biphenyl (PCB) congeners and biotransformable phthalate esters (PEs) in a marine aquatic food-web. Model testing showed no systematic bias and good precision in the estimation of the TMF for PCB congeners but an apparent underestimation of model calculated TMFs, relative to reported field TMFs, for PEs. A model sensitivity analysis showed that sampling designs that ignore the presence of concentration gradients may cause systematically biased and misleading TMF values. The model demonstrates that field TMFs are most sensitive to concentration gradients and species migration patterns for substances that are subject to a low degree of biomagnification or trophic dilution. The model is useful in anticipating the effect of spatial concentration gradients on the determination of the TMF; guiding species collection strategies in TMF studies; and interpretation of the results of field bioaccumulation studies in study locations where spatial differences in chemical concentration exist.

Bioaccumulation of decamethylpentacyclosiloxane (D5): A review.

Decamethylpentacyclosiloxane (D5) is a widely used, high-production volume personal care product with an octanol-water partition coefficient (log K(OW)) of 8.09. Because of D5's high K(OW) and widespread use, it is subject to bioaccumulation assessments in many countries. The present study provides a compilation and an in-depth, independent review of bioaccumulation studies involving D5. The findings indicate that D5 exhibits depuration rates in fish and mammals that exceed those of extremely hydrophobic, nonbiotransformable substances; that D5 is subject to biotransformation in mammals and fish; that observed bioconcentration factors in fish range between 1040 L/kg and 4920 L/kg wet weight in laboratory studies using non-radiolabeled D5 and between 5900 L/kg and 13 700 L/kg wet weight in an experiment using C(14) radiolabeled D5; and that D5 was not observed to biomagnify in most laboratory experiments and field studies. Review of the available studies shows a high degree of internal consistency among findings from different studies and supports a broad comprehensive approach in bioaccumulation assessments that includes information from studies with a variety of designs and incorporates multiple bioaccumulation measures in addition to the K(OW) and bioconcentration factor.

Bioconcentration and Aquatic Toxicity of Superhydrophobic Chemicals: A Modeling Case Study of Cyclic Volatile Methyl Siloxanes.

Many chemicals in commerce are classified as "superhydrophobic", having log octanol-water partition coefficients (log KOW) approaching or exceeding 7. Examples include long-chain alkanes, halogenated aromatics, and cyclic volatile methylsiloxanes (cVMS). We show that superhydrophobic chemicals present unique assessment challenges because of their sparing solubility in water and difficulties in empirical determinations of bioconcentration factors (BCFs) and aquatic toxicity. Using cVMS as an example, BCFs are considerably lower than expected due to biotransformation. Reviewed aquatic toxicity test data for cVMS in a range of aquatic organisms show little or no toxic effects up to solubility limits in water and sediment. Explanations for this apparent lack of toxicity of cVMS, and by extension to other superhydrophobic chemicals, are explored using a conventional one-compartment uptake model to simulate bioconcentration and toxicity tests using an assumed baseline narcotic critical body residue (CBR) and a range of organism sizes. Because of the low aqueous concentrations, equilibration times are very long and BCFs are sensitive to even very slow rates of biotransformation. Most organisms fail to achieve the assumed CBR during feasible test durations even at the solubility limit. Regulatory evaluation of superhydrophobic substances requires specially designed test protocols addressing biotransformation and dietary uptake.

Decamethylcyclopentasiloxane (D5) environmental sources, fate, transport, and routes of exposure.

The environmental sources, fate, transport, and routes of exposure of decamethylcyclopentasiloxane (D5; CAS no. 541-02-6) are reviewed in the present study, with the objective of contributing to effective risk evaluation and assessment of this and related substances. The present review, which is part of a series of studies discussing aspects of an effective risk evaluation and assessment, was prompted in part by the findings of a Board of Review undertaken to comment on a decision by Environment Canada made in 2008 to subject D5 to regulation as a toxic substance. The present review focuses on the early stages of the assessment process and how information on D5's physical-chemical properties, uses, and fate in the environment can be integrated to give a quantitative description of fate and exposure that is consistent with available monitoring data. Emphasis is placed on long-range atmospheric transport and fate in water bodies receiving effluents from wastewater treatment plants (along with associated sediments) and soils receiving biosolids. The resulting exposure estimates form the basis for assessments of the resulting risk presented in other studies in this series. Recommendations are made for developing an improved process by which D5 and related substances can be evaluated effectively for risk to humans and the environment.

Fugacity and activity analysis of the bioaccumulation and environmental risks of decamethylcyclopentasiloxane (D5).

As part of an initiative to evaluate commercial chemicals for their effects on human and environmental health, Canada recently evaluated decamethylcyclopentasiloxane (D5; CAS no. 541-02-06), a high-volume production chemical used in many personal care products. The evaluation illustrated the challenges encountered in environmental risk assessments and the need for the development of better tools to increase the weight of evidence in environmental risk assessments. The present study presents a new risk analysis method that applies thermodynamic principles of fugacity and activity to express the results of field monitoring and laboratory bioaccumulation and toxicity studies in a comprehensive risk analysis that can support risk assessments. Fugacity and activity ratios of D5 derived from bioaccumulation measures indicate that D5 does not biomagnify in food webs, likely because of biotransformation. The fugacity and activity analysis further demonstrates that reported no-observed-effect concentrations of D5 normally cannot occur in the environment. Observed fugacities and activities in the environment are, without exception, far below those corresponding with no observed effects, in many cases by several orders of magnitude. This analysis supports the conclusion of the Canadian Board of Review and the Minister of the Environment that D5 does not pose a danger to the environment. The present study further illustrates some of the limitations of a persistence-bioaccumulation-toxicity-type criteria-based risk assessment approach and discusses the merits of the fugacity and activity approach to increase the weight of evidence and consistency in environmental risk assessments of commercial chemicals.

Characterization of ecological risks from environmental releases of decamethylcyclopentasiloxane (D5).

Decamethylcyclopentasiloxane (D5) is used in personal care products and industrial applications. The authors summarize the risks to the environment from D5 based on multiple lines of evidence and conclude that it presents negligible risk. Laboratory and field studies show that D5 is not toxic to aquatic organisms or benthic invertebrates up to its solubility limit in water or porewater or its sorptive capacity in sediment. Comparison of lipid-normalized internal concentrations with measured concentrations in benthos indicates that field-collected organisms do not achieve toxic levels of D5 in their tissues, suggesting negligible risk. Exposure to D5 resulted in a slight reduction of root biomass in barley at test concentrations 2 orders of magnitude greater than measured D5 levels in biosolids-amended soils and more than twice as high as the maximum calculated sorptive capacity of the soil. No effects were observed in soil invertebrates exposed to similar concentrations, indicating that D5 poses a de minimis risk to the terrestrial environment. High rates of metabolism and elimination of D5 compared with uptake rates from food results in biodilution in the food web rather than biomagnification, culminating in de minimis risk to higher trophic level organisms via the food chain. A fugacity approach substantiates all conclusions that were made on a concentration basis.

An updated Quantitative Water Air Sediment Interaction (QWASI) model for evaluating chemical fate and input parameter sensitivities in aquatic systems: application to D5 (decamethylcyclopentasiloxane) and PCB-180 in two lakes.

The QWASI fugacity mass balance model has been widely used since 1983 for both scientific and regulatory purposes to estimate the concentrations of organic chemicals in water and sediment, given an assumed rate of chemical emission, advective inflow in water or deposition from the atmosphere. It has become apparent that an updated version is required, especially to incorporate improved methods of obtaining input parameters such as partition coefficients. Accordingly, the model has been revised and it is now available in spreadsheet format. Changes to the model are described and the new version is applied to two chemicals, D5 (decamethylcyclopentasiloxane) and PCB-180, in two lakes, Lake Pepin (MN, USA) and Lake Ontario, showing the model's capability of illustrating both the chemical to chemical differences and lake to lake differences. Since there are now increased regulatory demands for rigorous sensitivity and uncertainty analyses, these aspects are discussed and two approaches are illustrated. It is concluded that the new QWASI water quality model can be of value for both evaluative and simulation purposes, thus providing a tool for obtaining an improved understanding of chemical mass balances in lakes, as a contribution to the assessment of fate and exposure and as a step towards the assessment of risk.

Mathematical relationships between metrics of chemical bioaccumulation in fish.

Five widely used metrics of bioaccumulation in fish are defined and discussed, namely the octanol-water partition coefficient (KOW ), bioconcentration factor (BCF), bioaccumulation factor (BAF), biomagnification factor (BMF), and trophic magnification factor (TMF). Algebraic relationships between these metrics are developed and discussed using conventional expressions for chemical uptake from water and food and first-order losses by respiration, egestion, biotransformation, and growth dilution. Two BCFs may be defined, namely as an equilibrium partition coefficient KFW or as a nonequilibrium BCFK in which egestion losses are included. Bioaccumulation factors are shown to be the product of the BCFK and a novel equilibrium multiplier M containing 2 ratios, namely, the diet-to-water concentration ratio and the ratio of uptake rate constants for respiration and dietary uptake. Biomagnification factors are shown to be proportional to the lipid-normalized ratio of the predator/prey values of BCFK and the ratio of the equilibrium multipliers. Relationships with TMFs are also discussed. The effects of chemical hydrophobicity, biotransformation, and growth are evaluated by applying the relationships to a range of illustrative chemicals of varying KOW in a linear 4-trophic-level food web with typical values for uptake and loss rate constants. The roles of respiratory and dietary intakes are demonstrated, and even slow rates of biotransformation and growth can significantly affect bioaccumulation. The BCFK s and the values of M can be regarded as the fundamental determinants of bioaccumulation and biomagnification in aquatic food webs. Analyzing data from food webs can be enhanced by plotting logarithmic lipid-normalized concentrations or fugacities as a linear function of trophic level to deduce TMFs. Implications for determining bioaccumulation by laboratory tests for regulatory purposes are discussed.

Uncertainty analysis using a fugacity-based multimedia mass-balance model: application of the updated EQC model to decamethylcyclopentasiloxane (D5).

The EQuilibrium Criterion (EQC) model developed and published in 1996 was recently revised to include improved treatment of input partitioning and reactivity data, temperature dependence and an easier sensitivity and uncertainty analysis. This New EQC model was used to evaluate the multimedia, fugacity-based fate of decamethylcyclopentasiloxane (D5; CAS No. 541-02-6) in the environment over a temperature range of 1-25°C. In addition, Monte Carlo uncertainty analysis was used to quantitatively determine the influence of temperature and input partitioning and reactivity data on the behavior of D5 under various emission scenarios. Results indicated that emission mode was the most influential factor determining the fate and distribution of D5 in the model environment. When emitted to air and soil, D5 partitioned to and remained in the air compartment where rates of removal from degradation and advection processes were relatively rapid. In contrast, D5 emitted to water resulted in a substantial mass fraction of D5 being accumulated in the sediment compartment, where rates of removal from degradation and advection processes were slow. The mass distributions and fate of D5 in the model environment were strongly influenced by multiple input parameters, including temperature, the mode of emission (especially emission rate to water), KOC and half-life in air. As temperature decreased from 25°C to 1°C, KOC and half-life in air became increasingly more influential such that the mass distribution of D5 increased in air and decreased in sediment, resulting in decreased overall persistence.