Modelling bioaccumulation of oil constituents in aquatic species

Modelling and validation of polycyclic aromatic hydrocarbons emissions from offshore oil production facilities

The development of numerical models for investigating the risks and impact caused by human activities to the marine environment is important. Herein, the recently developed ChemicalDrift Lagrangian dispersion model was coupled to a toxicokinetic model and applied to investigate emissions of polycyclic aromatic hydrocarbons (PAHs) discharged from oil and gas production facilities as produced water. The performance of the model was evaluated with available data from a monitoring survey conducted at two oil fields. The survey provided exposure concentrations by means of passive samplers and bioaccumulation data in caged mussels; multiple depths and locations were assessed. The study included 26 PAHs and alkylated derivatives, showing good agreement between the model and the survey measurements. The compounds dominating the scenario were naphthalenes and phenanthrenes. Model provided contamination gradients were in agreement with the survey results, with levels decreasing with distance away from the main sources and with higher concentrations at 20 m depth. ChemicalDrift and the toxicokinetic model provided detailed time series, showing peaks of C1-naphthalene bioaccumulation significantly higher than values accumulated at the end of the monitored period. The utilised model was able to separate the relative contributions of multiple platforms and to identify the major contamination sources, providing a valuable and versatile tool for assessing the impact of discharges at sea.

Changes in Temperature Alter the Toxicity of Polycyclic Aromatic Compounds to American Lobster (Homarus americanus) Larvae

Polycyclic aromatic compounds (PACs) present in the water column are considered to be one of the primary contaminant groups contributing to the toxicity of a crude oil spill. Because crude oil is a complex mixture composed of thousands of different compounds, oil spill models rely on quantitative structure-activity relationships like the target lipid model to predict the effects of crude oil exposure on aquatic life. These models rely on input provided by single species toxicity studies, which remain insufficient. Although the toxicity of select PACs has been well studied, there is little data available for many, including transformation products such as oxidized hydrocarbons. In addition, the effect of environmental influencing factors such as temperature on PAC toxicity is a wide data gap. In response to these needs, in the present study, Stage I lobster larvae were exposed to six different understudied PACs (naphthalene, fluorenone, methylnaphthalene, phenanthrene, dibenzothiophene, and fluoranthene) at three different relevant temperatures (10, 15, and 20 °C) all within the biological norms for the species during summer when larval releases occur. Lobster larvae were assessed for immobilization as a sublethal effect and mortality following 3, 6, 12, 24, and 48 h of exposure. Higher temperatures increased the rate at which immobilization and mortality were observed for each of the compounds tested and also altered the predicted critical target lipid body burden, incipient median lethal concentration, and elimination rate. Our results demonstrate that temperature has an important influence on PAC toxicity for this species and provides critical data for oil spill modeling. More studies are needed so oil spill models can be appropriately calibrated and to improve their predictive ability. Environ Toxicol Chem 2023;42:2389-2399. © 2023 SETAC.

An impact-based environmental risk assessment model toolbox for offshore produced water discharges

We present a novel approach to environmental risk assessment of produced water discharges based on explicit impact and probability, using a combination of transport, fate and toxicokinetic-toxicodynamic models within a super-individual framework, with a probabilistic element obtained from ensemble simulations. Our approach is motivated by a need for location and species specific tools which also accounts for the dynamic nature of exposure and uptake of produced water components in the sea. Our approach is based on the well-established fate model DREAM, and accounts for time-variable exposure, considers body burden and effects for specific species and stressors, and assesses the probability of impact. Using a produced water discharge in the Barents Sea, with early life stages of spawning haddock, we demonstrate that it is possible to conduct a model-based risk assessment that highlights the effect of natural variations in environmental conditions. The benefits, limitations and potential for further improvements are discussed.

Bridging the lab to field divide: Advancing oil spill biological effects models requires revisiting aquatic toxicity testing

Oil fate and exposure modeling addresses the complexities of oil composition, weathering, partitioning in the environment, and the distributions and behaviors of aquatic biota to estimate exposure histories, i.e., oil component concentrations and environmental conditions experienced over time. Several approaches with increasing levels of complexity (i.e., aquatic toxicity model tiers, corresponding to varying purposes and applications) have been and continue to be developed to predict adverse effects resulting from these exposures. At Tiers 1 and 2, toxicity-based screening thresholds for assumed representative oil component compositions are used to inform spill response and risk evaluations, requiring limited toxicity data, analytical oil characterizations, and computer resources. Concentration-response relationships are employed in Tier 3 to quantify effects of assumed oil component mixture compositions. Oil spill modeling capabilities presently allow predictions of spatial and temporal compositional changes during exposure, which support mixture-based modeling frameworks. Such approaches rely on summed effects of components using toxic units to enable more realistic analyses (Tier 4). This review provides guidance for toxicological studies to inform the development of, provide input to, and validate Tier 4 aquatic toxicity models for assessing oil spill effects on aquatic biota. Evaluation of organisms' exposure histories using a toxic unit model reflects the current state-of the-science and provides an improved approach for quantifying effects of oil constituents on aquatic organisms. Since the mixture compositions in toxicity tests are not representative of field exposures, modelers rely on studies using single compounds to build toxicity models accounting for the additive effects of dynamic mixture exposures that occur after spills. Single compound toxicity data are needed to quantify the influence of exposure duration and modifying environmental factors (e.g., temperature, light) on observed effects for advancing use of this framework. Well-characterized whole oil bioassay data should be used to validate and refine these models.

North sea produced water PAH exposure and uptake in early life stages of Atlantic Cod

Produced water discharges from offshore oil and gas platforms represent a significant source of petroleum components such as polycyclic aromatic hydrocarbons (PAHs) released to the ocean. High molecular weight PAHs are persistent in the environment and have a potential for bioaccumulation, and the investigation of their fate and uptake pathways in marine life are relevant when assessing environmental risk of produced water discharges. To study the exposure and uptake of 2-5 ring PAHs in early life stages of Atlantic Cod in the North Sea, we run a coupled fate and individual-based numerical model that includes discharges from 26 platforms. We consider 26 different PAH components in produced water which biodegrade with primary degradation rates; intermediate degradation products are not included. Model simulations are run covering multiple years (2009-2012) to study annual exposure variability, while a one-day time slice of spawning products from the peak spawning season are followed. By covering multiple release points and large spatio-temporal scales, we show how individuals can be exposed to produced water from multiple regions in the North Sea. We find that a combination of oceanic fate processes and toxicokinetics lead to markedly different compositions in the predicted internal concentrations of PAHs compared to discharge concentrations; for instance, naphthalene makes up 30% of the total discharged PAHs, but contributes to at most 1% of internal concentrations. In all simulations we find the predicted total internal PAH concentration (26 components) to be below 1.2 nmol/g, a factor of 1000 less than concentrations commonly associated with acute narcotic effects.

Implications of Trophic Variability for Modeling Biomagnification of POPs in Marine Food Webs in the Svalbard Archipelago

The occurrence of persistent organic pollutants (POPs) in the Arctic has been of constant concern, as these chemicals cause reproductive effects and mortality in organisms. The Arctic acts as a chemical sink, which makes this system an interesting case for bioaccumulation studies. However, as conducting empirical studies for all Arctic species and POPs individually is unfeasible, in silico methods have been developed. Existing bioaccumulation models are predominately validated for temperate food chains, and do not account for a large variation in trophic levels. This study applies Monte Carlo simulations to account for variability in trophic ecology on Svalbard when predicting bioaccumulation of POPs using the optimal modeling for ecotoxicological applications (OMEGA) bioaccumulation model. Trophic magnification factors (TMFs) were calculated accordingly. Comparing our model results with monitored POP residues in biota revealed that, on average, all predictions fell within a factor 6 of the monitored POP residues in biota. Trophic variability did not affect model performance tremendously, with up to a 25% variability in performance metrics. To our knowledge, we were the first to include trophic variability in predicting biomagnification in Arctic ecosystems using a mechanistic biomagnification model. However, considerable amounts of data are required to quantify the implications of trophic variability on biomagnification of POPs in Arctic food webs.

Developmental effects in fish embryos exposed to oil dispersions - The impact of crude oil micro-droplets

During accidental crude oil spills and permitted discharges of produced water into the marine environment, a large fraction of naturally occurring oil components will be contained in micron-sized oil droplets. Toxicity is assumed to be associated with the dissolved fraction of oil components, however the potential contribution of oil droplets to toxicity is currently not well known. In the present work we wanted to evaluate the contribution of oil droplets to effects on normal development of Atlantic cod (Gadus morhua) through exposing embryos for 96 h to un-filtered (dispersions containing droplets) and filtered (water soluble fractions) dispersions in a flow-through system at dispersion concentrations ranging from 0.14 to 4.34 mg oil/L. After exposure, the embryos were kept in clean seawater until hatch when survival, development and morphology were assessed. The experiment was performed at two different stages of embryonic development to cover two potentially sensitive stages (gastrulation and organogenesis). Exposure of cod embryos to crude oil dispersions caused acute and delayed toxicity, including manifestation of morphological deformations in hatched larvae. Oil droplets appear to contribute to some of the observed effects including mortality, larvae condition (standard length, body surface, and yolk sac size), spinal deformations as well as alterations in craniofacial and jaw development. The timing of exposure may be essential for the development of effects as higher acute mortality was observed when embryos were exposed from the start of gastrulation (Experiment 1) than when exposed during organogenesis (Experiment 2). Even though low mortality was observed when exposed during organogenesis, concentration-dependent mortality was observed during recovery.

Model-based exploration of the variability in lake trout (Salvelinus namaycush) bioaccumulation factors: The influence of physiology and trophic relationships

Because dietary consumption of fish is often a major vector of human exposure to persistent organic pollutants (POPs), much effort is directed toward a quantitative understanding of fish bioaccumulation using mechanistic models. However, many such models fail to explicitly consider how uptake and loss rate constants relate to fish physiology. We calculated the bioaccumulation factors (BAFs) of hypothetical POPs, with octanol-water partition coefficients ranging from 10^4.5 to 10^8.5 , in lake trout (Salvelinus namaycush) with a food-web bioaccumulation model that uses bioenergetics to ensure that physiological parameters applied to a species are internally consistent. We modeled fish in 6 Canadian lakes (Great Slave Lake, Lake Ontario, Source Lake, Happy Isle Lake, Lake Opeongo, and Lake Memphremagog) to identify the factors that cause the BAFs of differently sized lake trout to vary between and within lakes. When comparing differently sized lake trout within a lake, larger fish tend to have the highest BAF because they allocate less energy toward growth than smaller fish and have higher activity levels. When comparing fish from different lakes, the model finds that diet composition and prey energy density become important in determining the BAF, in addition to activity and the amount of total energy allocated to growth. Environ Toxicol Chem 2019;38:831-840. © 2019 SETAC.

Partitioning of PAHs between Crude Oil Microdroplets, Water, and Copepod Biomass in Oil-in-Seawater Dispersions of Different Crude Oils

The impact of oil microdroplets on the partitioning of polycyclic aromatic hydrocarbons (PAHs) between water and marine zooplankton was evaluated. The experimental approach allowed direct comparison of crude oil dispersions (containing both micro-oil droplets and water-soluble fraction; WSF) with the corresponding WSF (without oil droplets). Dispersion concentration and oil type have an impact on the PAH composition of WSFs and therefore affect dispersion bioavailability. Higher T-PAH body residues were observed in copepods treated with dispersions compared to the corresponding WSFs. PAHs with log K_ow 3-4.5 displayed comparable accumulation factors between treatments; however, accumulation factors for less soluble PAHs (log K_ow = 4.5-6) were higher for the WSF than for the dispersions, suggesting low bioavailability for components contained in oil droplets. The higher PAH body residue in dispersion exposures is assumed to result mainly from copepods grazing on oil droplets, which offers an alternative uptake route to passive diffusion. To a large degree this route is controlled by the filtration rates of the copepods, which may be inversely related to droplet concentration.

How including ecological realism impacts the assessment of the environmental effect of oil spills at the population level: The application of matrix models for Arctic Calanus species

For oil spill responses, assessment of the potential environmental exposure and impacts of a spill is crucial. Due to a lack of chronic toxicity data, acute data is used together with precautionary assumptions. The effect on the Arctic keystone (copepod) species Calanus hyperboreus and Calanus glacialis populations is compared using two approaches: a precautionary approach where all exposed individuals die above a defined threshold concentration and a refined (full-dose-response) approach. For this purpose a matrix population model parameterised with data from the literature is used. Population effects of continuous exposures with varying durations were modelled on a range of concentrations. Just above the chronic No Observed Effect Concentration (which is field relevant) the estimated population recovery duration of the precautionary approach was more than 300 times that of the refined approach. With increasing exposure concentration and duration, the effect in the refined approach converges to the maximum effect assumed in the precautionary approach.

Numerical evaluation of bioaccumulation and depuration kinetics of PAHs in Mytilus galloprovincialis

Polycyclic aromatic hydrocarbons (PAHs) are important organic pollutants in the aquatic environment due to their persistence and bioaccumulation potential both in organisms and in sediments. Benzo(a)anthracene (BaA) and phenanthrene (PHE), which are in the priority pollutant list of the U.S. EPA (Environmental Protection Agency), are selected as model compounds of the present study. Bioaccumulation and depuration experiments with local Mediterranean mussel species, Mytilus galloprovincialis were used as the basis of the study. Mussels were selected as bioindicator organisms due to their broad geographic distribution, immobility and low enzyme activity. Bioaccumulation and depuration kinetics of selected PAHs in Mytilus galloprovincialis were described using first order kinetic equations in a three compartment model. The compartments were defined as: (1) biota (mussel), (2) surrounding environment (seawater), and (3) algae (Phaeodactylum tricornutum) as food source of the mussels. Experimental study had been performed for three different concentrations. Middle concentration of the experimental data was used as the model input in order to represent other high and low concentrations of selected PAHs. Correlations of the experiment and model data revealed that they are in good agreement. Accumulation and depuration trend of PAHs in mussels regarding also the durations can be estimated effectively with the present study. Thus, this study can be evaluated as a supportive tool for risk assessment in addition to monitoring studies.

Time-varying effects of aromatic oil constituents on the survival of aquatic species: Deviations between model estimates and observations

There is a need to study the time course of toxic chemical effects on organisms because there might be a time lag between the onset of chemical exposure and the corresponding adverse effects. For aquatic organisms, crude oil and oil constituents originating from either natural seeps or human activities can be relevant case studies. In the present study the authors tested a generic toxicokinetic model to quantify the time-varying effects of various oil constituents on the survival of aquatic organisms. The model is based on key parameters applicable to an array of species and compounds with baseline toxicity reflected by a generic, internal toxicity threshold or critical body burden (CBB). They compared model estimates with experimental data on the effects of 8 aromatic oil constituents on the survival of aquatic species including crustaceans and fish. The average model uncertainty, expressed as the root mean square error, was 0.25 (minimum-maximum, 0.04-0.67) on a scale between 0 and 1. The estimated survival was generally lower than the measured survival right after the onset of oil constituent exposure. In contrast, the model underestimated the maximum mortality for crustaceans and fish observed in the laboratory. Thus, the model based on the CBB concept failed to adequately predict the lethal effects of the oil constituents on crustaceans and fish. Possible explanations for the deviations between model estimates and observations may include incorrect assumptions regarding a constant lethal body burden, the absence of biotransformation products, and the steady state of aromatic hydrocarbon concentrations in organisms. Clearly, a more complex model approach than the generic model used in the present study is needed to predict toxicity dynamics of narcotic chemicals. Environ Toxicol Chem 2017;36:128-136. © 2016 SETAC.

Crude oil affecting the biomass of the marine copepod Calanus finmarchicus: Comparing a simple and complex population model

In the current study differences were evaluated between a complex 3D multistage population model (SINMOD) and a simpler consumer-resource population model for estimating the effects of crude oil on the marine copepod Calanus finmarchicus. The SINTEF OSCAR model was used to simulate hypothetical oil spills in the Lofoten area in 1995, 1997, and 2001. Both population models simulated a negligible effect of crude oil on the Calanus' biomass when assuming low species sensitivity. The simple model estimated a larger effect on the biomass (up to a 100% decline) compared to the complex model (maximum decline of 60-80%) at high species sensitivity to crude oil. These differences may be related to the inclusion of copepod advection in the complex model. Our study showed that if little data is available to parameterize a model, or if computational resources are scarce, the simple model could be used for risk screening. Nevertheless, the possibility of including a dilution factor for time-varying biomass should be examined to improve the estimations of the simple model. The complex model should be used for a more in depth risk analysis, as it includes physical processes such as the drift of organisms and differentiation between developmental stages.

A fugacity-based toxicokinetic model for narcotic organic chemicals in fish

A novel dynamic fugacity-based model is described, developed, and tested that simulates the uptake of narcotic organic chemicals in fish from water as occurs in aquatic bioconcentration and toxicity tests. The physiologically based toxicokinetic model treats the time course of chemical distribution in 4 compartments (tissue groups) in the fish, including the liver, in which biotransformation may occur. In addition to calculating bioconcentration and toxicokinetics, 5 possible toxic endpoints are defined corresponding to chemical concentration, fugacity, or activity reaching a critical value that causes 50% mortality. The mathematical description of multicompartment uptake is simplified by expressing the equations in the fugacity format. The model is parameterized and tested against reported empirical data for the bioconcentration of pentachloroethane in rainbow trout and for uptake and mortality from aquatic exposures to naphthalene and 1,2,4-trichlorobenzene in fathead minnows. Model performance is evaluated, and it is concluded that with suitable parameterization it has potential for application for assessment of both bioconcentration and toxicity expressed as median lethal concentrations, critical body residues, and chemical activity as a function of time to death.

Release, deposition and elimination of radiocesium ((137)Cs) in the terrestrial environment

Radionuclide contamination in terrestrial ecosystems has reached a dangerous level. The major artificial radionuclide present in the environment is (137)Cs, which is released as a result of weapon production related to atomic projects, accidental explosions of nuclear power plants and other sources, such as reactors, evaporation ponds, liquid storage tanks, and burial grounds. The release of potentially hazardous radionuclides (radiocesium) in recent years has provided the opportunity to conduct multidisciplinary studies on their fate and transport. Radiocesium's high fission yield and ease of detection made it a prime candidate for early radio-ecological investigations. The facility setting provides a diverse background for the improved understanding of various factors that contribute toward the fate and transfer of radionuclides in the terrestrial ecosystem. In this review, we summarize the significant environmental radiocesium transfer factors to determine the damaging effects of radiocesium on terrestrial ecosystem. It has been found that (137)Cs can trace the transport of other radionuclides that have a high affinity for binding to soil particles (silts and clays). Possible remedial methods are also discussed for contaminated terrestrial systems. This review will serve as a guideline for future studies of the fate and transport of (137)Cs in terrestrial environments in the wake of the Fukushima Nuclear Power Plant disaster in 2011.

Toxicity data for modeling impacts of oil components in an Arctic ecosystem

Ecological impact assessment modeling systems are valuable support tools for managing impacts from commercial activities on marine habitats and species. The inclusion of toxic effects modeling in these systems is predicated on the availability and quality of ecotoxicology data. Here we report on a data gathering exercise to obtain toxic effects data on oil compounds for a selection of cold-water marine species of fish and plankton associated with the Barents Sea ecosystem. Effects data were collated from historical and contemporary literature resources for the endpoints mortality, development, growth, bioaccumulation and reproduction. Evaluating the utility and applicability of these data for modeling, we find that data coverage is limited to a sub-set of the required endpoints. There is a need for new experimental studies for zooplankton focused on the endpoints development and bioaccumulation and for larvae and juvenile fish focused on growth and development.