Page under development
ECOTRACER: Using Ecosim to predict movement and accumulation of contaminants and tracers in food webs
Ecosim predicts temporal changes flows of biomass among living and detritus pools using nonlinear functional relationships between flow rates and abundances of the interacting species. These flow rates (along with auxiliary information about factors such as isotope decay rate and physical exchange rates) can be used to predict changes in concentrations (and per-biomass burdens) of chemicals like organic contaminants and isotope tracers that ?flow? passively along with the biomass flows. The dynamic equations for such passive flow (and accumulation, e.g. ?bioamplification?) are not the same as the biomass flow rate equations, and in fact are generally linear dynamical equations with time-varying rate coefficients that depend on the biomass flow rates; these linear equations are relatively easy to solve in parallel with the Ecosim biomass dynamics equations.
In Ecosim we allow parallel simulation of one tracer or contaminant type while the biomass dynamics equations are being solved. Tracer molecules are assumed to be either in the ?environment?, or in the biota (in biomass and detritus pools) at any moment. Molecules are assumed to flow between pools at instantaneous rates (i.e. to be sampled along with biomass during biomass flows at rates) equal to the probabilities of being ?sampled? as part of the biomass flow: instantaneous rate=(flow)/(biomass in prey pool). We also allow for direct flows from the environment into pools, representing direct uptake or absorption of the tracer material, and for differential decomposition/decay/export rates by pool and from the environmental pool. Schematically, the flow of tracer molecules through any biomass pool is represented by the components shown below:
In the rate equation for time changes in contaminant concentration in pool i, these components are represented as follows:
Hence Ecotracer users must specify some parameters (using the Ecotracer input form) besides those needed for Ecopath / Ecosim mass balance and biomass dynamics calculations:
The other parameters (GCi, MOi) are already be entered or calculated in Ecopath, and the Qij flows are calculated dynamically in Ecosim. Ecotracer users can also enter inflow, outflow (dispersal/advection exchange to outside systems), and decay rates for the ?environmental? (abiotic) concentration Co.
When the Ecotracer simulation component is enabled (check box on Ecotracer entry form), Ecosim integrates the rate components listed above to generate time concentration patterns Ci (i=0 to number of biomass pools), and time patterns of concentration per biomass Ci/Bi for all Ecopath / Ecosim biomass pools. These results can be compared to data on bioaccumulation or tracer movement. Discrepancies between model and data may help identify weaknesses in the Ecopath/ Ecosim trophic flow rate (Qij) estimates, and/or in the Ecotracer rate parameters (ui,ei,di).