56 Introduction and dynamics

Ecotracer is a useful tool within the Ecopath with Ecosim (EwE) modelling approach to track radioisotopes, contaminants, persistent pollutants, or stable isotopes through a food web model.  After achieving a mass‐balanced Ecopath model, Ecotracer can be used with the Ecosim (time dynamic) or Ecospace (spatial‐temporal dynamic) to track the flow of the pollutant through the modelled ecosystem.

EwE uses a mass‐balance approach, and the flows of a contaminant due to predator/prey interactions are tracked within the underlying Ecopath model. However, Ecotracer also needs parameters for groups based on a kinetic toxicology approach to estimate initial conditions. However, similar to Ecopath, Ecotracer can become dynamic through either the use of Ecosim or Ecospace to follow the changes of a contaminant that has different temporal inputs or to variations in temporal spatial concentrations in the water column.

The purpose of this manual is to give a full description of the Ecotracer approach including the dynamic equations that describe the basis for the input parameters, and to familiarize users with the various interfaces for the inputs and outputs. Simulation scenarios are also given to allow users to become more familiar with Ecotracer, and instructions are given on how to navigate through the different interfaces used.

Introduction

Ecotracer is a sub‐routine in the Ecopath with Ecosim (EwE) modeling framework^{[1] [2]} that allows the modeller to follow a contaminant or stable isotope in modelled functional groups and the environment in a balanced Ecopath model.

Many EwE models have been made that focus on fishery‐related questions, but here the focus is on how the Ecotracer routine is used within EwE to trace contaminants such as radionuclides through an aquatic ecosystem. EwE consists of three routines: Ecopath which is a mass balance interpretation of an ecosystem where, in essence, the production of a group in the model is equal to its consumption; Ecosim allows the user to build in time dynamics to the Ecopath model for events such as changes in contaminant loading to an ecosystem; and Ecospace which allows for the spatial‐temporal resolution for such events as the effects of change in loading on marine organisms that result from organisms inhabiting different spatial areas or habitat types that have different environmental concentrations through time.

Typical applications of Ecotracer have been for contaminants such as mercury^[3], ¹⁴C^{[4] [5]}, ¹³⁷Cs^[6] [7], and PCBs^[8] that can have detrimental impacts on human and environmental health. The use of Ecotracer can help to estimate the amount of contaminant or concentration in a group/species of interest, spatial differences in concentration within the same functional groups if Ecospace is used, the fluxes between groups due to trophic interactions, and the importance of diet versus direct environmental uptake. It can also help to estimate functional groups’ concentrations when such data are lacking (i.e., have a starting value of zero) and make forward projections based on changing environmental concentrations. Concentration levels are an important aspect for environmental and human health as, in conjunction with consumption rates, they determine exposure levels that may have detrimental effects. Regulatory limits on the concentration in aquatic products destined for human consumption may also affect trade and fisheries opportunities.

Ecotracer requires a balanced Ecopath model to follow the contaminant, radioisotopes or stable isotopes in the model groups, and environment (e.g., water concentration). Ecotracer when used with Ecosim can provide estimates to important ecotoxicological questions such as,

1. what could be the expected group concentrations if the environmental concentration did not change?
2. what could be the expected results in group concentrations if the environmental concentration changed through time? and
3. is there an effect on concentration levels as a result of changing underlying Ecopath parameters such as fishing mortality?

The first question is useful if many functional groups in the model lack concentration data. The second question can be important to estimate resulting concentrations in biota if the input into the environment changes. The third can help to understand contaminant flows as a result of changes in the dynamics of the underlying structure of an ecosystem.

Ecotracer when used with Ecospace can help to answer whether there are differences in the same species that occur over a large geographic area, and whether different environmental concentrations in different areas impact the resulting concentrations in organisms. In this case, a two‐dimensional representation of the model area is made which has a user defined spatial resolution (i.e., grid cells). Spatial environmental concentrations can be driven by effluents being released as a point source, or from atmospheric deposits that change over space and time. Effluents released from a point source would be affected by currents resulting in different spatial and temporal distributions, and atmospheric releases could be affected by different levels of releases due to industrial activity through time or accident scenarios as well as currents. In the case of large‐scale accidents, such as the Dai‐ichi nuclear accident at Fukushima, spatial differences can result from both point sources and differing atmospheric deposits both of which occurred through time.

Ecotracer dynamics

Ecotracer simulates the contaminant fluxes and resulting amounts and concentrations using a modified transfer contaminant model (e.g.,^{[9] [10]}), and applies it to both the environment and biota. Resulting changes at any time step are dependent upon the gains and losses in functional groups and are described in Walters and Christensen^[11] as,

[latex]\frac{dA_i(t)}{dt} = \alpha_i - \beta_i A_i(t)\tag{1}[/latex]

where α_i represents the gains (Bq∙year^‐1) in each functional group i, β_i represents the rate losses (year^‐1) to each functional group, and A_i represents the amount (e.g., Bq) in each functional group i. This general formulation allows different measurement units of substances (e.g., Bq or µg) to be tracked in the modelled environment, and the resulting concentrations (e.g., Bq∙t^‐1) are computed separately using the biomass output in Ecopath and Ecosim.

The environmental compartment concentrations are also calculated by tracking the gains and losses in the cells representing the environment,

[latex]\frac{dC_o(t)}{dt} = \alpha_o - \beta_o C_o(t)\tag{2}[/latex]

where C_o represents the environmental concentration (e.g., Bq∙km^‐2), α_o, represents the gains and losses in each environmental cell o, and β_o represents the rate losses (year^‐1) in each cell.

Environment

Gains in the environment originate from the release of contaminants into the environment as a base inflow rate, and from the excretion from organisms. Losses originate from the direct uptake from the environment by organisms, physical decay rates, and base volume exchange. In Ecospace, the environment can be represented by multiple grid cells and thus the gains and losses can be considered to be for each environmental compartment o, such that,

[latex]\alpha_o = BI_o + \sum\limits_{i=1}^{n}m_iA_i\tag{3}[/latex]

where BI_o is the base inflow rate (Bq∙km^‐2∙year^‐1) to a grid cell, and m_iA_i are the excretory products for each functional group within each grid cell.

Losses from the environment are due to biological, physical decay processes, environmental volume changes, and direct uptake by organisms, such that,

[latex]\beta_o C_o=(d_i+V_i)C_o+\sum\limits_{i=1}^{n}u_iB_iC_o\tag{4}[/latex]

where d_i represents the physical decay rate (year^‐1), V_i represents the base volume exchange loss (year^‐1), and the second term (u_iB_iC_o) represents the total uptake rate by all functional groups (see below). Temporal changes to the environmental concentration (C_o) can be made by applying a forcing function to the base inflow rate, through a contaminant concentration driver file (Table 1), or by current/advection fields.

Biota

In biota, intake amounts (e.g., Bq∙year‐1) result from direct uptake rates (i.e., respiration) the fraction retained from trophic interactions (i.e., diet), and immigration. i.e.,

[latex]\alpha_i=u_iB_iC_o+AE_i \sum\limits_{i=1}Q_{ij} \frac{A_j}{B_j}+c_iI_i\tag{5}[/latex]

where, C_o represents the environmental concentration (Bq∙km^‐2), Bi is the biomass (t) of group i, u_i represent the intake/biomass/environmental concentration/year (km²∙t^‐1∙year^‐1); AE_i is the assimilation efficiency for each group, Q_ji is the consumption rate (t∙year^‐1) of group j by group i, A_j is the amount of substance in a group (e.g., Bq), B_j is the prey biomass of each prey item j (Bq∙t^‐1); c_i is the group biomass concentration (Bq∙t^‐1) and I_i is the immigrating biomass (t∙year^‐1).

The losses from a group (β_iC_i) are attributed to predation, fisheries, other mortality, excretion and decay, i.e.,

[latex]\beta_iC_i=(\sum\limits_{j=pred}(\frac{Q_{ij}}{B_i})+F_i+MO_i+E_i+m_i+d_i) C_i\tag{6}[/latex]

where Q_ij is the rate of consumption (t∙year^‐1) of group i due to predation by j, F_i is the fishing mortality rate (year_‐1), MO_i (year^‐1) is other mortality rate (i.e., non‐predation mortality), E_i is the emigrating biomass rate (year^‐1), m_i (year^‐1) is the excretion and/or metabolic rate, and d_i (year^‐1) is the physical decay rate. These rates are multiplied by C_i the amount of contaminant (Bq) in each group i. Excretory products that are released from tissues to the environment are added to the environmental concentration.

The solution for finding the equilibrium amount of contaminant in a primary producer with the resulting concentration only being due to direct uptake, losses to due predation, other mortality, metabolism, and decay is given as,

[latex]C_{i,eq}=\frac{u_iB_iC_o}{\sum\limits_{j=pred}(\frac{Q_{ij}}{B_i}) +MO_i+m_i+d_i}\tag{7}[/latex]

whereas for other groups an additional term must be accounted for due to the group’s prey items; in these cases the equilibrium solution can be defined as,

[latex]C_{i,eq}=\frac{u_iB_iC_o+AE_i \frac{Q_{ji}}{B_i}}{\sum\limits_{j=pred}(\frac{Q_{ij}}{B_i}) +MO_i+m_i+d_i}\tag{8}[/latex]

The Ecotracer approach is dynamic and extends the basic concentration ratio (CR) approach, but the CR approach is contained within it as,

[latex]CR_i=\frac{(A_i/B_i)}{C_o}=\frac{u_i+AE_i\frac{Q_{}ji}{B_i}CR_j}{Z_i+m_i}\tag{9}[/latex]

The amount of contaminant in the detritus compartment originates from the unassimilated consumption resulting from predation, as well as non‐predation mortality. Thus, groups feeding on detritus will have exposure levels associated with the contributions from the fraction of unassimilated consumption from all groups. Initial concentrations in the biota and environment are also input parameters that can be used if data is available. For groups lacking contaminant data from field studies or literature data, the model is able to estimate concentration or burdens in the groups, leading to the ability to estimate risk through time and make comparisons to regulatory limits.

.