United States
                    Environmental Protection
                    Agency
Environmental Research
Laboratory
Athens GA 30613
-A'
                    Research and Development
EPA-600/S3-83-084  Nov. 1983
ŁEPA          Project  Summary
                    Prediction  of  Ecological  Effects
                    of  Toxic  Chemicals:  Overall
                    Strategy  and  Theoretical   Basis
                    for the  Ecosystem  Model
                    Ray R. Lassiter
                     A strategy was developed for model-
                    ing ecosystems to permit the assess-
                    ment of toxic chemical effects on
                    element cycling and other ecosystem
                    processes. The strategy includes the
                    use of multi-species representations of
                    biotic communities and mathematical
                    descriptions of the processes that are
                    important  in  aquatic ecosystems.
                    Direct effects of toxicants are assigned
                    to the species comprising the biotic
                    community, in a manner suggested by
                    available  toxicological  information.
                    Effects are calculated as the difference
                    between  selected  measures of
                    processes from unaffected systems and
                    systems affected by the presence of a
                    toxic chemical. Ecological effects cal-
                    culated in this manner are considered to
                    be heuristically useful.
                     This Project Summary was developed
                    by  EPA's Environmental  Research
                    Laboratory, Athens. GA, to  announce
                    key findings of the research project that
                    is fully documented in a separate report
                    of  the same  title (see Project Report
                    ordering information at back).

                    Introduction
                     Analysis of environmental risks posed
                    by the release of toxic chemicals neces-
                    sarily would require predictions of the
                    fate of a chemical and its effects. Pre-
                    diction of a chemical's fate presupposes
                    an  understanding   of  transport  and
                    transformation processes, or  "what the
                    environment  does to the  chemical."
                    Prediction of a chemical's effects presup-
                    poses  a  wider knowledge  of  the
                    chemical's fate, and by contrast, "what
the chemical does to the environment." It
is also necessary to know which effects
are of concern, because very many poten-
tial effects exist. Ecological effects are the
concern of this research, and effects that
occur in ecosystems, and whose charac-
teristics  are influenced by ecological
interactions, are considered to be ecolog-
ical effects. A particular class of potential
ecological effects that arguably could be
of major significance are effects on cycles
of elements. Although results  of this
research will not be limited to considera-
tion of the effects of toxicants on element
cycles, they are the class of effects that
have motivated most of the work. A major
reason for  focusing on  element cycle
effects is that there is apparently no on-
going work in this area despite the appar-
ent potential for such effects.

The Nature of Element Cycles
and Importance of Considering
Their Response to Toxicants
 Biotic processes in the proportions and
magnitudes that we know them depend
upon  concentrations of chemicals
remaining approximately at their current
levels. The energy for biotic processes is
derived from catalysis of redox reactions
between reduced and oxidized forms of
elements. The primary source of the dis-
equalibria that drive these reactions is
photosynthesis in which both  highly
reduced organic compounds and
molecular oxygen are produced. The
remaining biotic processes feed off this
disequilibrium, and in doing so contribute
to maintenance of the pool sizes  of geo-

-------
chemicals. These include the atmospher-
ic gases, the composition of the ocean,
dissolved  constituents of  fresh waters,
and the earth's soils.
  These geochemical pools have been
remarkably stable  over geological time
scales. This fact attests to the historical
stability of the contributory processes. No
comparative analysis of the relative roles
of  abiotic  and  biotic  processes  in
maintaining the pools of  geochemicals
has  been made.  The  number  and
importance of the processes known to be
carried out by organisms, however, leaves
little doubt that the characteristics of the
earth depend on the continuation of biotic
processes at  about their current levels.
Whether these processes will  remain
stable  in  the  presence  of  continued
introduction of xenobiotic chemicals is
unknown.  It  is  prudent,  therefore, to
consider potential conditions  that could
lead to serious disturbances of element
cycle transformations.  Conditions that
would  lead to alterations  in ecosystem
structure  and function  are expected to
lead to  corresponding  alterations  in
element cycling.

Strategy of Ecosystem Modeling
for Analyzing Alterations in
Element Cycling
  Element cycling is  a system phenom-
enon. Organisms use pools of chemicals
that  are products  of other organisms'
metabolism. These and other phenomena
that occur as fluxes that change chemical
forms and redox states are processes. For
a given set of elements one can describe
the major processes  that occur to cycle
the elements through their various forms.
They can be grouped into those pertain-
ing to microrganisms or macroorganisms
and to autotrophs  or heterotrophs.  The
level  of  resolution  at  which   these
processes  are  described  is  a  critical
choice  for accomplishing the desired
goals. These  goals are to  represent the
aquatic ecosystem such that  models of
the  processes  respond   to  external
influences, including the introduction of
toxicants approximately the same as do
real systems.  If this goal  is reached, it will
be  possible  to analyze for  ecological
effects of toxic chemicals, including their
effects  on element  cycling processes.
Carrying out analyses of this sort is,  in a
special sense, a predictive (or prognostic)
exercise that is fully comprehensible only
if the usage of the term "effect" is under-
stood,   and  if the  special  sense of
prediction is made clear.
  An effect implies a deviation from some
nominal state or rate. In this model, there-
fore.the effect is referred to as a deviation
from a nominal condition. This deviation
can be obtained by first using the model to
obtain the nominal condition, then using
it again in the same manner, but with the
presence  of  a  toxicant  represented
additionally. To normalize effects so that
measurement units do not influence the
magnitude of the values, the difference
should be divided by the nominal value.
Any quantity of interest can be used to
calculate an effect, from the amount of a
material  present to the rate at which a
process  is  occurring,   regardless   of
whether there is any direct influence of
the toxicant on the quantity of interest.
Because the effect is obtained  by model
calculation in  which  comparable
conditions can be strictly maintained, and
it  is  known that  the only difference
between the nominal and affected condi-
tions is the presence of the chemical, any
differences  can be   attributed to  the
presence of the chemical and referred to,
therefore, as an effect. These calculations
are referred to as predictions of ecological
effect of the toxic chemical.
  The special sense  in which "predic-
tions" is used is that the predicted values
have  heuristic value but do not refer to an
expected effect  in   a  particular  real
ecosystem at some time. The  predicted
value is intended to be heuristically useful
by permitting assessment of the probable
effects in aquatic ecosystems of selected
characteristics without necessarily refer-
encing specific real aquatic systems. The
strategy of modeling is designed to make
this use possible.
  The strategy of modeling makes  use of
representation of multiple species, each
carrying  out a  particular process. It  is
often  observed that  concentrations of
environmental chemicals that are known
to turn  over  rapidly  are fairly  stable,
whereas any particular species that uses
or produces the chemical may fluctuate
widely. This suggests that compensation
occurs among species carrying out  the
same function. That is, reduction in levels
of one species may be compensated by
corresponding increases in the levels of
one or more others, so that the processes
responsible for the level of the chemical
are relatively stable.  The possibility of
stability  resulting  from  species
redundancy  and compensation   is
important when considering the effects
of chemicals  on the cycling  of  major
elements.  If  this   phenomenon   is
important, it is probable that representing
a group of organisms by a single variable
(as opposed to a multi-species represen-
tation) would lead to conclusions that
element  cycling  processes  are  more
fragile than they,  in fact, are.
  Another  important component in this
modeling strategy is the way that element
cycling  processes are  represented by
species populations. In element cycling, it
is the processes that are important. They
are carried out, however, by particular
populations. Each  individual species pop-
ulation is characterized by size, energy
source,   source  of   materials  for
biosynthesis,  motility, and several other
descriptors. The descriptors will define a
particular population as carrying out a
specific  process.   Representations  of
multi-species   populations  will  differ
among  themselves  in   the  particular
values of the descriptors. For example,
one  population may consist of  large
individuals that move slowly and use
biomass  of  microorganisms  for  both
energy and biosynthesis, while another
may differ  in being smaller and  moving
more  rapidly.  These  organisms  are
represented differently from carrying out
different processes by differences in the
energy source and source of materials for
biosynthesis.  For example, another group
of  organisms might use NH4+and 02 for
energy  and  inorganic   chemicals for
biosynthesis.
  Representing the organisms that exist
in  systems to carry out all the important
functions using this approach requires
that a large amount of descriptive data be
available. That amount of data will not be
available for all the organisms necessary
to  represent the processes of any given
ecosystem. Therefore, the use  of this
approach requires an additional strategic
component to make it possible to parame-
terize the model for the  large number of
species representations  necessary to the
approach.  Essentially   this  will  be
accomplished by  first describing a large
number  of species for the initial set of
organisms that occupy the environment.
Selection of the descriptors will be gu ided
by  what is known about the kinds of
organisms carrying  out the   various
ecosystem processes, /. e., they will be
assigned to the  model  representations
over a range that includes what is known.
It  is unlikely  that an assignment can be
made such that the entire initial set will
persist for very long (not approach  zero).
   It is anticipated, however, that some of
the variables will be persistent. These will
represent a  set of organisms  that are
inter-compatible and that are "adapted"
to  the  environment.   It is  extremely
unlikely that  this  set will well represent
real populations  (in the sense of  close
correspondence of their descriptors to

-------
those  of  real  populations).  Yet,  their
descriptors  necessarily will  be in the
neighborhood of those of real populations,
because only such descriptors are to be
used. Repeated selection of compatible
sets of organism representations for a
given environment beginning with differ-
ent initial sets is expected to build a
distribution  of  descriptors that would
include the characteristics of real organ-
isms that would live in the environment
as described.
  This  synthetic approach appears to be
the best  way to  obtain descriptions of
biotic communities  that are appropriate
for  an  environment that is selected  for
assessment in  anticipation of a  toxic
chemical's  possible introduction.  With
this approach  a  nominal  behavior of
element  cycling   processes   can  be
obtained. The same type of uncertainty
exists  for the  direct effects  of a  toxic
chemical  on  the  organisms  of  a
community,  however,  as exists for the
characteristics of types  of organisms that
inhabit the community.  The final element
of the  strategy defines the process by
which toxicological information is used to
obtain  affected  behavior.
  For  each  set of  compatible species
representations, a set of direct effects will
be  assigned from  the range of direct
effects  suggested   by  the  available
toxicological  information.  Recalculation
of the  behavior in the presence of the
toxicant,  then,  will   build  another
distribution  of behaviors.  Normalized
differences between the two behaviors,
as already discussed, are the ecological
effects of the chemical.
  It is  readily apparent that ecological
effects calculated in this manner are not
expected to be  observed exactly in any
particular  system.  They are  expected,
however,  to include  values  for   real
systems  similar  to  the one  analyzed.
Furthermore, if the nature of the distribu-
tion of  descriptors of the compatible sets
of species and of the direct effects of toxi-
cants on them  is  realistic, a measure of
the  uncertainty  associated  with  the
predicted effects  can be obtained. The
value of the calculated effects will not be
primarily  quantitative, however,   but
rather, will more appropriately be used as
a learning exercise for identification of
possible ecological effects that are too far
removed  from  the  direct  effects to  be
apparent without  an ecosystem model.
Descriptions of the Processes
  Much  of the  research  effort  was
devoted to assimilation of existing and
development of new mathematical de-
scriptions of the processes. Uptake of
materials by microorganisms is described
by a rectangle hyperbola that is derived in
a way that rationalizes its use for all
uptake processes. Depletion of materials
being taken up from a zone in the vicinity
of cells is considered. Photoautotrophy is
described in terms of both bioenergetics
and uptake of materials, so that limitation
by light and materials is possible. Chemo-
autotrophs and photoautotrophs share
the same biosynthetic pathways. Chemo-
autotrophs,  however, use a variety of
terminal electron acceptors. The bioener-
getics of their energy reactions is used to
calculate the rate of energy production,
and as with photoautotrophs and for that
matter all organisms, limitation of growth
is considered to be the minimum of the
potential  growth rates for each of the
requirements. Photoautotrophs produce
organic  matter  using  solar   energy,
whereas,   chemoautotrophs   produce
organic  matter  using energy of redox
reactions between  reduced  forms of
inorganic compounds and oxygen. In the
process, oxidized forms of elements are
formed.  The  reduced  forms  of  the
inorganic compounds are produced by
heterotrophs using  oxidized species of
elements in anaerobic zones.
  Heterotrophic  processes  are  repre-
sented  for  both  microorganisms  and
macroorganisms. All use organic com-
pounds for both energy and for synthesis.
Microorganisms  are  found  in  both
aerobic  and anaerobic zones of  aquatic
systems.  In  aerobic zones  they  use
molecular  oxygen   as  their  terminal
electron  acceptor  (oxidizing  agent) for
oxidizing organic compounds. In anaero-
bic  zones oxidized forms of  the major
elements are used as terminal electron
acceptors. These oxidized forms are those
produced  by chemoautotrophy.  In  this
manner,  the   element  cycles  are
completed. Growth  of  heterotrophs is
described as the difference between the
available energy and materials and their
need for maintenance and activity. Macro-
heterotrophs differ  from  microhetero-
trophs in being limited  to aerobic envi-
ronments and  m their use of specific
structures and behaviors to obtain food.
Two major kinds of feeding are described
for  the model: filter feeding and  pursuit
feeding. General model descriptions for
both types are developed for the general
case of any number of consumers feeding
on  any  number of prey. Rates of con-
sumption  are calculated  directly from
characteristics  of  predator  and prey,
rather than assigned as coupling constants.
Outlook
  The  strategy  of modeling  and  the
requisite theory to carry it  out exist.
Additional work is needed to bring the
expression of toxicological theory to the
state of development  (in the  context of
this model) as the  theory for ecosystem
modeling. Further work is  needed to
represent fermentation as an ecological
process.  Coding is underway  for  the
ecosystem model.  Whether the  sets of
compatible  representations  of  biotic
communities and the effects as computed
will be usefully close to real systems
needs to be tested. The fact that no speci-
fic real system is represented presents a
difficulty in that regard. Laboratory eco-
systems, however,  can be run repeatedly
under the same conditions and with or
without the  presence  of a toxicant. One
potential testing scheme would make use
of this capability and  compare distribu-
tions and effects from  laboratory ecosys-
tems with distributions and effects from
the model as described.

-------
     The Project Summary was authored by Ray R. Lassiter who is also the EPA
       Project Officer (see below).
     The complete report, entitled "Prediction of Ecological Effects of Toxic Chemicals:
       Overall Strategy and Theoretical Basis for the Ecosystem Model," (Order No. PB
       83-261 685; Cost: $10.00, subject to change) will be available only from:
             National Technical Information Service
             5285 Port Royal Road
             Springfield, VA 22161
              Telephone: 703-487-4650
     The EPA Project Officer can be contacted at:
             Environmental Research Laboratory
             U.S. Environmental Protection Agency
             Athens, GA 30613
                                                   4U.S. GOVERNMENT PRINTING OFFICE: 1M3-459-017/7227
United States
Environmental Protection
Agency
Center for Environmental Research
Information
Cincinnati OH 45268
  BULK RATE
 U.S. POSTAG
     PAID
Cincinnati, Oh
 Permit No. G3
Official Business
Penalty for Private Use $300

-------