| | 1 | == (To be inserted in TOC) Contaminant Tracing == |
| | 2 | |
| | 3 | '''''Page under development''''' |
| | 4 | |
| | 5 | ECOTRACER: Using Ecosim to predict movement and accumulation of contaminants and tracers in food webs |
| | 6 | |
| | 7 | 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. |
| | 8 | |
| | 9 | 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: |
| | 10 | |
| | 11 | Your browser may not support display of this image. In the rate equation for time changes in contaminant concentration in pool i, these components are represented as follows: |
| | 12 | |
| | 13 | 1. Uptake from food: C,,''j'',,GC,,''i'',,Q,,''ji'',,/B,,''j'',, where C,,''j'',,=conc in food ,,''j'',,, GC,,''i'',,=proportion of food assimilated by type ,,''i'',, organisms; Q,,''ji'',,=biomass flow rate from ,,''j'',, to ,,''i'',,, B,,''j'',,=food ,,''j'',, biomass; |
| | 14 | 2. Direct uptake from environment: u,,''i'',,B,,''i'',,C,,''o'',,, where u,,''i'',,=parameter representing uptake per biomass per time, per unit environmental concentration, B,,''i'',,=biomass, C,,''o'',,=environmental concentration; |
| | 15 | 3. Concentration in immigrating organisms: c,,''i'',,I,,''i'',,, where c,,''i'',,=parameter (tracer per unit biomass in immigrating biomass), I,,''i'',,=biomass of pool I immigrants per time; |
| | 16 | 4. Predation: C,,''i'',,Q,,''ij'',,/B,,''i'',,, where C,,''i'',,=concentration in pool ,,''i'',,, Q,,''ij'',,=consumption rate of type ,,''i'',, organisms by predator type ,,''j'',,, B,,''i'',,=biomass in pool ,,''i'',,. |
| | 17 | 5. Detritus: C,,''i'',,MO,,''i'',,+(1-GC,,''i'',,)S,,''j'',,C,,''j'',,Q,,''ji'',,/B,,''j'',,), where MO,,''i'',,=non-predation death rate of type ,,''i'',, (yr-1), GC,,''i'',,=fraction of food intake assimilated, Q,,''ji'',,=intake rate if type ,,''j'',, biomass by type ,,''i'',,; |
| | 18 | 6. Emigration: e,,''i'',,C,,''i'',,, where e,,''i'',,=emigration rate (yr-1); |
| | 19 | 7. Metabolism: d,,''i'',,C,,''i'',,, where d,,''i'',,=metabolism+decay rate for the material while in pool ,,''i'',,. |
| | 20 | |
| | 21 | Hence Ecotracer users must specify some parameters (using the Ecotracer input form) besides those needed for Ecopath/ecosim mass balance and biomass dynamics calculations: |
| | 22 | |
| | 23 | 1. initial pool concentrations C,,''i'',,, including environmental concentration C,,''o'',,; |
| | 24 | 2. direct uptake rate parameters u,,''i'',, as rates per time per biomass per unit C,,''o'',,; |
| | 25 | 3. conentrations per biomass c,,''i'',, in immigrating organisms; |
| | 26 | 4. metabolism/decay rates d,,''i'',,. |
| | 27 | |
| | 28 | The other parameters (GC,,''i'',,, MO,,''i'',,) are already be entered or calculated in Ecopath, and the Q,,''ij'',, 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 C,,''o'',,. |
| | 29 | |
| | 30 | 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 C,,''i'',, (,,''i'',,=0 to number of biomass pools), and time patterns of concentration per biomass C,,''i'',,/B,,''i'',, 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 (Q,,''ij'',,) estimates, and/or in the Ecotracer rate parameters (u,,''i'',,,e,,''i'',,,d,,''i'',,). |
| | 31 | |
| | 32 | |
