Environmental Health Indicators: Jake Rice
Environmental Health Indicators: Jake Rice
Abstract
Ecologists have proposed hundreds of quantitative indicators of the status of ecosystems for
evaluation of and reporting on the status of marine ecosystems. The talk applies a common
approach to classifying indicators, summarises the main properties of each class of indicator,
and provides some illustrations.
Indicators of ecosystem status have roles in both communication and decision support. For
both roles, strengths and weaknesses of indicators are usually only partially known. Few have
been tested systematically for sensitivity and robustness across a range of contextsy
However, to determine best practices for selecting indicators for specific uses, one must have a
fairly complete understanding of the information content of the various indicators. The paper
explores some alternative approaches to documenting the information content of various
indicators of ecosystem status.
As the Precautionary Approach becomes broadly used as the basis for management
decision-making, the role of indicators of ecosystem status becomes central. The PA is made
operational in decision-making through use of indicators and reference points. This means
that it is necessary to identify values of an ecosystem indicator associated with harm to the
environment that is serious or difficult to reverse. Methods for identifying and justifying
reference points for ecosystem indicators are being developed and tested, but the task is
turning out to be complex. Alternative strategies for identifying reference points are reviewed.
When choices must be made from a suite of candidate indicators, it is desirable make the
selection on objective grounds. This requires explicit a priori criteria, on which there is not yet
scientific consensus. Recent developments in this area are reviewed, and again a way forward
proposed.
r 2003 Elsevier Science Ltd. All rights reserved.
0964-5691/03/$ - see front matter r 2003 Elsevier Science Ltd. All rights reserved.
doi:10.1016/S0964-5691(03)00006-1
236 J. Rice / Ocean & Coastal Management 46 (2003) 235–259
1. Introduction
Even with those things known, choice of the correct indicators is not
straightforward. It requires explicit criteria for performance of the candidate
indicators, and objective tools for evaluating performance against the criteria.
The talk will develop these considerations and criteria more fully. It will not
provide the definitive answer to the question of which indicator(s) to use in
Integrated Coastal Zone Management. Rather, it will help users of Indicators of
Ecosystem Status ask better questions when selecting the indicator(s) to use. If the
question ‘‘Which Indicator(s) to use?’’ is asked in a better way, better answers may
be found.
Ecosystem indicators have been classified in many ways. Four classes have proved
effective at organising discussion on ecosystem effects of fishing [2]. Indicator species
are also included because the class is used frequently in some environmental quality
studies, although in resource management applications single-species indicators,
although common, are rarely thought of as indicators of larger ecosystem properties.
J. Rice / Ocean & Coastal Management 46 (2003) 235–259 237
Diversity is a joint construct of both how many species are present in a collection
(richness), and how similar their abundances are (evenness). Of the dozens of
diversity indices (and their close relatives, similarity indices) that have been proposed
since Preston [11], many vary only in the relative weight given to richness vs.
evenness [12]. Some indices additionally try to emphasise dominance; the role of the
most important species in a community [13]. Recently, Warwick and Clarke [14]
have proposed diversity indices that also include taxonomic similarity as well as
information about the number and relative abundance of species in a collection.
In trying to collapse all information on species richness and relative abundance
into a single number, any index of diversity can be misleading in (at least) two ways.
Two communities with different numbers of species can have similar estimated
‘‘diversities’’, if the distribution of abundances across species differs in a reciprocal
manner (i.e. the richer community is dominated by a few very common species;
whereas species in the species-poor community are similar in abundance).
Alternatively, two communities which are quite similar in all their common species
can still have different estimated ‘‘diversities’’, if many rare species are recorded in
one of the communities, whereas few are recorded in the other. Another possible
source of concern is that diversity indices intrinsically treat all species as equally
informative with regard to community structure and impacts.
None of these problems are insurmountable, because the relative sensitivity of the
various indices to richness, evenness, and taxonomic relatedness are known [15]. This
knowledge can be used to select the appropriate index for specific applications,
although the selection offers a number of meta-problems. To choose wisely which
diversity indices to use, one must know which of the major forcers perturbs the
community most strongly. To illustrate, fishing rarely extirpates marine species [16],
so if fishing is an important pressure, diversity indices which weight richness more
strongly than evenness will be insensitive. However, pollution often causes loss of
species intolerant to the chemical while allowing the few most tolerant species to
thrive, so indices emphasising richness are often recommended [17].
In summary, if a researcher or manager is interested in one specific type of
perturbation, and knows a substantial amount about the ecosystem of interest, it
may be possible to choose a diversity index with the proper weighting of richness,
evenness, and dominance. However, such a strategy allows one’s expectations to
influence strongly the analytical results to be obtained, which is rarely desirable. The
strategy also requires that the perturbation of interest be the dominant one affecting
the system, whereas the point of integrated coastal zone management is to address
multiple interacting forcing factors.
small number of dimensions whose axes are mathematically translatable into to the
original community matrix. In this reduced space sites with similar species
compositions are close together and sites with very different species compositions
are far apart. Statistically independent gradients usually are represented as
orthogonal ordination axes. This gives the scores for individual species or sites on
the ordination axes some desirable statistical properties [18]. The computational
steps of identifying the underlying gradients and allocating cases on the gradients are
almost always automated. It is the ecological interpretation of the ordination axes
that requires skill and judgements.
The various classes of ordination methods differ in the assumptions made about
the underlying statistical distributions of the abundances of species and the
functional relationship between the species’ abundances and the underlying
gradients. For example, the frequently encountered principal components analysis
(PCA), assumes that the abundance data have normal error structure, zero
abundances are rare, and gradient(s) are well sampled across their full range(s).
As a consequence, PCA and its close relatives are only applicable when the
ecological range sampled is narrow, with similar species present over the full range of
sites. Correspondence analyses can be used when the gradient is longer and the
abundance matrix is sparser [19]. However, correspondence analysis performs poorly
when all species do not share a common, monotonic functional form for their
distribution; a problem when dealing with a mixture of both widely distributed
species and specialists, or which species that aggregate.
Non-metric scaling (MDS) provides an alternative to dealing with problematic
underlying distributions of abundances. With MDS the similarities or dissimilarities
among sites are estimated using the weaker assumption that the rank order of
abundance of a species across sites is informative, but the actual quantitative
estimates of abundance may not be [20]. MDS then represents the species by site
matrix in a space defined by many fewer dimensions, but preserving as well as
possible the (dis)similarities of cases in the smaller dimensional space. This weaker
assumption of ordinal rather than interval information increases robustness in the
face of irregular distributions of abundance and high sampling variance. As a result
MDS has become a preferred technique for ecological ordinations of marine
communities [21]. However, it is not without its problems. MDS methods require
specifying a priori the number of axes which exist and a starting configuration for
patterns in the data. These a priori requirements create substantial opportunity for
the analyst’s preconceptions to influence, if not dominate, the analytical results.
MDS also requires selecting a measure of (dis)similarity, which brings with it all the
complexities discussed under diversity indices.
Choosing among ordination methods can be done objectively, based on statistical
properties of the site-abundance matrix. Interpreting results in applied contexts may
not be so easy. The simplest interpretations are by narrative; the positions of cases
on the extracted axes are examined, and ecological inferences are drawn from
knowledge of the cases. Opportunity for expectation bias is high. More rigorously,
the values of cases on the ordination axis can be regressed on some environmental
variable, reflecting intensity of perturbation due the forcing factor. Unfortunately,
240 J. Rice / Ocean & Coastal Management 46 (2003) 235–259
unless the forcing factor is such a dominant signal that it describes a community
gradient of its own, its effect may be masked altogether when extracting orthogonal
axes.
If there were no risk of committing errors when using indicators, then the choice of
indicators would be a matter of taste and cost-effectiveness. On the other hand, if
errors can be made, then their nature and the associated consequences should be
242 J. Rice / Ocean & Coastal Management 46 (2003) 235–259
only part of the story, even if to the technical expert that part of the story is
explained very well.
These types of communications errors are not unique to indicators of ecosystem
status, nor to indicators of any type. Neither is the protection against them unique to
this context. Clarity and care in language is paramount, as is awareness of the target
audience. Avoidance of value-laden terms is a virtue. Testing messages with focal
groups is effective when there is a specific message that needs to be conveyed to a
large audience. All these things are good common sense. Although communications
problems will occur, vigilance and co-operation with appropriate experts in
communications arts are the keys to keeping them infrequent and minimally
damaging.
How these probabilities are used then becomes a function of the governance
system. In areas with long management histories and strong regulatory capacities,
there can be quite rigid control rules. In advance of crisis situations managers and
stakeholders agree on specific management actions, tied to specific probabilities of
violating limit reference points and/or achieving targets. The decision support
consists of estimating the probability that the current state of the indicator violates
the limit (or achieves the target), and taking the predetermined action corresponding
to the estimated risk. This approach occurs in a few places in resource management
[38], and with many human health issues, where, for example beaches may be closed
to recreation or shellfish harvesting when coliform bacteria levels exceed pre-
identified reference points. In resource management contexts, though, the manage-
ment system often uses softer rules. The scientific advice still reports the probability
of violating limits and achieving targets, and may even recommend actions
consistent with a precautionary approach [39]. However, the decision-making uses
this information as one of many factors, although often an influential one, in the
final decision.
Regardless of the governance use of the information, estimation of two properties
associated with indicators is fundamental to effective decision support. One is
reference points that are reliably indicative of conservation or human health
concerns, and of provision of desired benefits. The other is the probability that the
current state of the environment complies with the reference points.
signals about anything else. Managers need to be warned when environmental events
have occurred in the real world so that measures to minimise harm or maximise
opportunity can be implemented promptly. Managers do not want frequent
erroneous warnings, however, prompting them to act unnecessary, and possibly
creating problems where none existed.
Signal detection theory in built on the 2 2 matrix of real-world events and
indicator status, using the four possible combinations, labelled:
* hits—event occurred and the signal says something happened,
* misses—event occurred but the signal is indistinguishable from background noise,
* false alarms—nothing actually happened but the signal says something did,
* true negatives—nothing happened, and there was only background noise in the
signal.
A perfect signal produces only hits and true negatives, but in a noisy environment,
a perfect signal is impossible. The theoretical framework was developed to optimise
performance of imperfect signal detection systems. In this optimisation signal
detection theorists focused on several factors that can guide in the selection of
indicators.
(1) The costs of the two types of errors may not be symmetrical, and the asymmetry
usually can be specified in advance. By ‘‘cost’’ in our context, we mean the
consequences of managers not acting when needed vs. acting when the actions
are unnecessary.
(2) Using knowledge of the relative cost of the two types of errors, indicators can be
selected whose error patterns best match the cost ratio.
(3) Knowing the absolute cost of errors guides investments in various monitoring
and signal detection systems. Do the benefits of fewer errors from a more
reliable indicator justify the cost of improved monitoring?
The parallel between developing reliable signals and selecting informative
environmental indicators is direct. Informative environmental indicators tell
managers about changes in ecosystem status. There will always be background
variation (noise) in the indicator, from which signals requiring action must be
differentiated. The differentiation can be done probabilistically in a risk management
framework, or as a binary decision of a reference point either being exceeded or not.
Science advisors, coastal managers, and stakeholders can explore in advance what
the possible consequences of misses and false alarms are, both to the ecosystem and
to those using the coastal zone. Misses are likely to have very different costs than
false alarms, with the costs of false alarms often borne by society, and costs of misses
borne by the ecosystem (and, of course. later by society). The interactions between
error tolerance rates and costs to improve indicator performance can be discussed
among all parties with a role in governance.
Thus, we have a framework for objective screening of candidate indicators of
environmental status. Select case histories where, with all the benefits of hindsight,
we know when some events occurred of magnitudes that coastal zone managers
should have reacted. Calculate the values of the suite of candidate indicators for the
246 J. Rice / Ocean & Coastal Management 46 (2003) 235–259
areas, and see when the risk quantified for each indicator relative to its reference
point would be large enough that technical advisors would have recommended
action. Which ones have error rates within our tolerances? Which ones have miss and
false alarm ratios that match our cost tolerances? Are the best error rates achievable
good enough to provide the desired likelihood of achieving coastal zone manage-
ment objectives?
Using a signal detection framework to select environmental indicators does not
guarantee that we will find an effective basis for science advice on coastal zone
management. Informative metrics with acceptable error rates and distributions may
not exist, or be prohibitively expensive. It may be impossible to get consensus among
diverse, legitimate stakeholders in what are acceptable error rates and distributions.
These situations would still be successes for the approach, because it would have
informed managers and stakeholders what to expect from the indicators being used
in making decisions about coastal zone management, that they assume are science-
based and risk averse.
There is cause for cautious optimism here. When an indicator species was chosen
specifically because of its relationship to a particular forcer, one would expect a
priori that fluctuations in the indicator species relative to its reference point would be
responding to changes in the forcer, and be informative to managers about
appropriate actions. This should be true when the indicator species itself is the
ecosystem property of interest, which is the case with species on protected or
endangered lists. It is less likely to be met when the species was chosen primarily for
convenience or to be a focus for stakeholder engagement. When the presence of a
species is diagnostic of particular conditions, such as a pollution-tolerant species,
false alarms may be rare, but miss rates might be high if the species has additional
ecological restrictions or slow dispersal.
The experience with fisheries is informative here. Fisheries are a particularly strong
forcer, affecting the target species directly (in addition to all their other effects [41]).
J. Rice / Ocean & Coastal Management 46 (2003) 235–259 247
Ordination methods are vulnerable to both misses and false alarms. Management
applications ordinate the multispecies samples from a block of years, using the
248 J. Rice / Ocean & Coastal Management 46 (2003) 235–259
resultant species and site scores on major axes as a reference framework. Scores of
monitoring sites are then tracked through ordination space over time. When sites
make big excursions, big ecosystem changes at those sites are inferred. Because
ordination methods are variance-structuring tools, the initial ordination is
dominated by the species that were most variable in the reference years. If these
are species whose abundances are inherently highly variable in space and time or
there is a lot of sampling error, ordinations will give many false alarms.
Alternatively, species which were initially either uncommon everywhere, or
abundant and widespread but without clearly defined optima will get little weight
in the initial ordination axes. Hence major increases in rare species, or declines in
common but eruptive species, will not show up as big excursions of sites in
ordination space, and be misses.
It is, of course, possible to ordinate and interpret a full time series. If there were
long-term trends, either driven by environmental forcing or incremental effects of
anthropogenic perturbations, they could be apparent in the systematic ordering
of samples over time along latent axes. Misses would be rare if the forcing factor of
interest contributed a meaningful and systematic portion of the total variance, so it
was an independent axis itself. Misses would be common if different parts of the
system responded to the forcing factor at different rates. Not all extracted axes
would be informative about the effects of the forcing factor. False alarms could be
common or rare, depending on the ability to interpret the information content of
the axes in non-circular ways (that is, in ways that used information independent of
the ordination itself to evaluate which axes were reflected effects of forces to which
managers should respond).
Compared to diversity indices calculated from the same data, size spectrum
methods should provide few misses or false alarms of changes in overall productivity,
transfer efficiency, or mortality rate. False alarms should be uncommon
because ln(number) spectra reduce distortions due to occasional sampling anomalies
or strong year classes of one or a few species. If strong recruitments occur for
several species or years, this is an increase in system productivity, and will affect the
metric’s value. The same is true, although possibly to a lesser extent, for dominance
curves. They may become more markedly humped with sampling or recruitment
anomalies, but unless productivity and survivorship of the rest of the species being
sampled is affected at the same time, the dominance curve is unlikely to change
substantially.
Both summary indicators will miss some types of ecosystem change, such as high
turnover rates of uncommon species. However, few indicators will have a high hit
rate for such types of community changes, where meta-populations of many
uncommon species are being lost and re-colonised over time. Moreover, system
changes due to several uncommon species becoming even rarer will be captured well
by dominance curves, and if they are larger individuals, by size spectra as well. These
types of changes might be of particular concern to coastal zone managers.
J. Rice / Ocean & Coastal Management 46 (2003) 235–259 249
There is a large literature and extensive experience with setting reference points for
indicator species. For exploited species, many science advisory and management
agencies use limit and precautionary or buffer reference points as core tools (ICES
NAFO, NASCO). Different jurisdictions may use somewhat different criteria for
positioning the reference points on abundance and mortality rate indicators, some
using more aggressive strategies than the Fmsy and Bmsy suggested by UNFA. What
matters though, is that once the criteria are made explicit, estimating both the
reference point and the probability that a current value violates the reference point,
is just one more task in the assessment of stock status. If the data for the stock make
a traditional assessment straightforward, it will be straightforward to estimate
the reference point and risk as well. If the assessment is difficult and unstable, so will
be the task of estimating reference points and risks.
The threatened or endangered species community also has extensive experience
with using quantitative reference points. The IUCN Red List of 1996 [47] is a
benchmark in using quantitative reference points for determining the appropriate
category for species being evaluated. Although the IUCN SSG work refers to
guidelines and criteria, in practice they are using reference points in combination
with soft decision rules. Great effort when into determining what the criteria should
be, and the values have generated substantial controversy [48]. The experience
illustrates again, though, that estimation of reference points for indicator species is
practical, as long as historic data are available, and are informative about what
states comprise harm that is serious or difficult to reverse.
There is another message for coastal zone management in the experience of
reference points used with fisheries and endangered species. These are the simplest
cases, but they have been neither easy nor free of controversy. Stock assessments
have not inspired uniform belief outside (or even within) the fisheries community.
Application of IUCN criteria or their derivatives has been hotly contests for some
types of species, and even the experts can argue at length about the proper category
for a species, given the quantitative criteria, reference points, and the same data
available to all. We should expect the task to be only harder for more complex
classes of indicators.
because the components of diversity, particularly richness, vary with many factors;
latitude, productivity, substrate, exposure, depth, etc. This makes it hard to justify
taking a value that is well justified as a reference point for one area, and using it as a
de facto reference point for other areas.
Ordination scores give slightly more opportunity for objectively setting reference
points than do diversity indices. When the only information available is the data
used in the ordination then the situation is identical for the two types of indicators.
Historic data can be used to provide estimates of ordination scores characteristic of
previous states of the ecosystem. However, only if some of those states were
unacceptable are they informative about where in ordination space reference limits
(points on one axis, boundaries in spaces of more than one axis) should be located.
Resampling methods can provide estimates of uncertainty of specific values of
ordination scores, and when all the assumptions of the ordination are met,
uncertainties can be estimated parametrically as well, so risk quantification should
usually be possible.
The potential gain with ordination scores comes when there are data available for
even a subset of sites, regarding independent variables that may be directly indicative
of the forcer of interest; for example contaminant levels, nutrient loading, etc. In
such cases, methods such as direct gradient analysis or canonical correspondence
analysis [49] can be used to relate trends in the forcer to trends in the ordination
scores. As long as one can justify extrapolation beyond the parameterisation data,
the results can then be used to identify the location(s) in ordination space associated
with unacceptable levels of the forcer and/or state of the ecosystem, even if the
conditions have not yet been observed. These then can be used as objectively
determined reference points.
Aggregate indicators like size spectra and dominance curves have not been used in
formal decision support contexts, although k-dominance curves have been
considered in evaluating pollution effects. The strong theoretical basis for size
spectra, and the corresponding empirical and modelling support for systematic
behaviour of the slope of the spectrum as a function of fishing mortality both suggest
that in principle it should be possible to establish a slope for a spectrum that would
correspond to an unsustainable mortality rate for the community as a whole. This
would require some initialisation data, from sampling done when the community
was considered healthy. The mortality rate associated with the observed slope of that
community could be the basis for modelling to establish the mortality rate that
would be unsustainable. The slope associated with that mortality would then be
the reference point. Typical regression methods can estimate the uncertainty in
the estimate of the slope from annual samples, so the risk quantification relative to
the reference slope would be possible. Reference points for the intercept of the size
252 J. Rice / Ocean & Coastal Management 46 (2003) 235–259
In Sections 4–6 I have laid out a rigorous framework for objectively selecting
indicators whose performance standards would be known. All the steps in applying
the framework are practical, but it reflects a noteworthy workload. Groups like
SCOR Working Group 119 (Ecosystem Indicators of Fisheries Effects), and a large
J. Rice / Ocean & Coastal Management 46 (2003) 235–259 253
research team in Australia, are pursuing initiatives this ambitious. In the medium
term it is realistic to expect the tools to be developed that would make operating with
the necessary rigour widely possible. However, in the short term, it is likely that
indicators will continue to be chosen on the basis of expert opinion and group
consensus, rather than objectively determined performance standards. Towards that
end a Task Group within SCOR WG 119 (Marie-Joelle Rochet, Philippe Cury, and
I) have been trying to codify practical screening criteria for other Task Groups to use
in selecting among indicators of ecosystem status for various system properties. This
work extended work reported by the ICES Ecosystem Effects of Fishing Working
Group, in their recommendations for scientific advice on Ecosystem Qualities and
Ecosystem Quality Objectives to OSPAR and the North Sea Council of Ministers.
Those developing selection criteria would be relevant for coastal zone managers,
faced with the need to choose among a large suite of possible indicators.
theories that are hotly contested among professionals may pose a high risk of mis-
informing users, if parts of the contested theory must be revised in future. In
evaluating this property for an indicator, one must evaluate not only the
consistency of an indicator with ecological theory, but also the degree to which
the diversity of professional views all accept the theoretical arguments.
7. Conclusions
References
[1] Belfiore S, Balgos M, Galofre J, McLean B, Blaydes M, Tesch D. Background Paper; International
Workshop on the Role of Indicators in Integrated Coastal Management. IOC document, 2002.
[2] Rice JC. Evaluating fishery impacts using metrics of community structure. ICES Journal of Marine
Science 2000;57:682–8.
[3] Jennings S, Kaiser MJ. The effects of fishing on marine ecosystems. Advances in Marine Biology
1998;34:201–352.
[4] ICES. Report of the ICES advisory committee on fisheries management, 2001; ICES Cooperative
Research Report 246(3 volumes);
Philippert CJM. Long-term impact of bottom fisheries on several bycatch species of demersal fish
and benthic invertebrates in the south-eastern North Sea. ICES Journal of Marine Science
1998;55:342–52.
[5] Blacker RW. Benthic animals as indicators of hydrographic conditions and climatic change in
Svalbard waters. Fishery Investigations (London) Ser. II 1957;20:1–49.
[6] ICES. Report of the ICES advisory committee on the marine environment, 1999. ICES Cooperative
Research Report No. 239, 1999. 277pp.
[7] Spellerberg IF. Monitoring ecological change. Cambridge, UK: Cambridge University Press, 1991.
334pp.
[8] Menge BA. Indirect effects in marine rocky intertidal interaction webs: patterns and importance.
Ecological Monographs 1995;65:21–74;
Zacharias MA, Roff JC. Use of focal species in marine conservation and management: a review and
critique. Aquatic Conservation: Marine and Freshwater Ecosystems 2001;11:59–76.
[9] ICES. Report of the working group on ecosystem effects of fishing. ICES CM 1996/Assess/Env:1,
1996.
[10] OSPAR BDC. Testing of the Faial criteria for the selection of species and habitats which need to be
protected with the deep-water and habitat-forming coral reef species Lophelia pertusa. Paper 00/6/5
presented by Norway to OSPAR Biodiversity Committee, 2000.
[11] Preston FW. On the commonness and rarity of species. Ecology 1948;29:254–83.
[12] Hill MO. Diversity and evenness: a unifying notation and its consequences. Ecology 1973;54:427–32;
Magurran AE. Ecological diversity and its measurement. London: Chapman & Hall, 1988.
[13] McNaughton J. Relationship among functional properties of California grassland. Nature
1967;216:168–9.
[14] Warwick RM, Clarke KR. New ‘‘biodiversity’’ measures reveal a decrease in taxonomic distinctness
with increasing stress. Marine Ecology Progress Series 1995;129:301–5.
[15] Ludwig JA, Reynolds JF. Statistical ecology. New York: Wiley, 1988. 337pp;
Magurran AE. Ecological diversity and its measurement. London: Chapman & Hall, 1988.
[16] Powles H, Bradford MJ, Bradford RG, Doubleday WG, Innes S, Levings CD. Assessing and
protecting endangered marine species. ICES Journal of Marine Science 2000;57:669–76;
FAO. An appraisal of the suitability of the CITES criteria for listing commercially exploited aquatic
species. FAO Fisheries Circular, vol. 954, 2000. 66pp.
[17] Spellerberg IF. Monitoring ecological change. Cambridge, UK: Cambridge University Press, 1991.
334pp;
Diez I, Secilla A, Santolaria A, Gorostiagia JM. Phytobenthic intertidal community structure along
an environmental pollution gradient. Marine Pollution Bulletin 1999;38:463–72.
[18] Gauch HG. Multivariate analysis in community ecology. Cambridge UK: Cambridge University
Press, 1982. 298pp;
Jongman RHG, ter Braak CJF, van Tongeren OFR. Data analysis in community and landscape
ecology. Pudoc-Wageningen, Den Haag, 1987. 299pp.
[19] Hill MO. Correspondence analysis: a neglected multivariate method. Applied Statistics 1974;23:
340–54.
[20] Kruskal JB. Multi-dimensional scaling by optimizing goodness of fit to a non-parametric hypothesis.
Psychometrica 1964;29:1–27;
J. Rice / Ocean & Coastal Management 46 (2003) 235–259 257
McRae G, Camp DK, Lyons WG, Dix TL. Relating benthic infauna community structure to
environmental variables in estuaries using non-metric multi-dimensional scaling and similarity
indices. Environmental Monitoring and Assessment 1998;51:233–46.
[21] Clarke KR, Ainsworth M. A method of linking multivariate community structure to environmental
variables. Marine Ecology Progress Series 1993;92:205–19;
Cao Y, Bark AW, Williams WP. Measuring the response of macroinvertebrate communities to water
pollution: a comparison of multivariate approaches, biotic and diversity indices. Hydrobiologia
1996;341:1–19.
[22] Rice JC, Gislason H. Patterns of change in the size spectra of numbers and diversity of the
North Sea fish assemblage, as reflected in surveys and models. ICES Journal of Marine Science
1996;53:1214–25;
Bianchi G, Gislason H, Graham K, Hill L, Jin X, Koranteng K, Manckchand-Heileman S, Paya I,
Sainsbury K, Sanchez F, Zwanenburd K. Impact of fishing on sinze composition and diversity of
demersal fish communities. ICES Journal of Marine Science 2000;57:558–71.
[23] Clarke KR. Comparisons of dominance curves. Journal of Experimental Marine Biology and
Ecology 1990;138:143–57;
Spellerberg IF. Monitoring ecological change. Cambridge, UK: Cambridge University Press, 1991.
334pp.
[24] Boudreau PA, Dickie LM. Biomass spectra of aquatic ecosystems in relation to fisheries yield.
Canadian Journal of Fisheries and Aquatic Science 1992;49:1528–38;
Thibaux ML, Dickey LM. Structure of the body-size spectrum of the biomass in aquatic ecosystems:
a consequence of allometry in predator–prey interactions. Canadian Journal of Fisheries and Aquatic
Science. 1993;50:1308–17.
[25] Gislason H, Rice JC. Modelling the response of size and diversity spectra of fish assemblages to
changes in exploitation. ICES Journal of Marine Science 1998;55:362–70;
Gislason H, Lassen H. On the linear relationship between fishing effort and the slope of the size
spectrum. ICES CM 1997/DD:05, 1997;
Bianchi G, Gislason H, Graham K, Hill L, Jin X, Koranteng K, Manckchand-Heileman S, Paya I,
Sainsbury K, Sanchez F, Zwanenburd K. Impact of fishing on sinze composition and diversity of
demersal fish communities. ICES Journal of Marine Science 2000;57:558–71.
[26] Beukema JJ. An evaluation of the ABC method (abundance biomass comparison) as applied to
macrozoobenthic communities living on tidal flats in the Dutch Wadden Sea. Marine Biology
1988;99:425–33;
Clarke KR. Comparisons of dominance curves. Journal of Experimental Marine Biology and
Ecology 1990;138:143–57.
[27] Lambshead PJD, Platt HM, Shaw KM. The detection of differences among assemblages of marine
benthic species based on an assessment of dominance and diversity. Journal of Natural History
1983;17:859–74;
Clarke KR. Comparisons of dominance curves. Journal of Experimental Marine Biology and
Ecology 1990;138:143–57.
[28] ICES. Report of the working group on ecosystem effects of fishing. ICES CM 1998/ACFM/ACME:1;
ICES 2001. Report of the working group on ecosystem effects of fishing. ICES CM ACME:01, 1998.
[29] Pimm SL. Food webs. New York: Chapman & Hall, 1982;
Cohen JE, Briand F, Newman CM. Community food webs: data and theory. Berlin: Springer, 1990.
308pp.
[30] Yodzis P. Introduction to theoretical ecology. New York: Harper & Row, 1989;
Yodzis P. Predator-prey theory and management of multispecies fisheries. Ecological Applications
1994;4:51–8.
[31] Pimm SL, Rice JC. The dynamics of multi-species, multi-life stage models of aquatic food webs.
Theoretical Population Biology 1987;32:303–25.
[32] Breitburg DL, Sanders JG, Gilmore CC, Hatfield CA, Osmun RW, Riedel GF, Seitzinger SP, Sellner
KG. Variability in responses to nutrients and trace elements, and transmission of stressor effects
through an estuarine food web. Limnology and Oceanography 1999;44:837–63.
258 J. Rice / Ocean & Coastal Management 46 (2003) 235–259
[33] Pauly D, Christensen V, Walters C. Ecopath, ecosim, and ecospace as tools for evaluating ecosystem
impact of fishing. ICES Journal of Marine Science 2000;57:697–706;
Walters C, Christensen V, Pauly D. Structuring dynamic models of exploited ecosystems from trophic
mass-balance assessments. Reviews in Fisheries and Fish Biology 1997;7:139–72.
[34] Jarre-Teichmann A. Seasonal mass-balance models of carbon flow in the central Baltic Sea with
emphasis on upper trophic levels. ICES CM 1995/T:6, 1995. 26pp;
Jarre-Teichmann A. The potential use of mass-balance models for the management of upwelling
ecosystems. Ecolological Applications 1998;8: 93–103;
Pauly D, Christensen V, Walters C. Ecopath, Ecosim, and Ecospace as tools for evaluating ecosystem
impact of fishing. ICES Journal of Marine Science 2000;57:697–706.
[35] Dalsgaard J, Jarre-Teichmann A, Walters C, Pauly D. An approach to the modelling of persistent
pollutants in marine ecosystems. ICES CM/1998, vol. 10, 1998;
Dalsgaard J, Wallace SS, Salas S, Preikshot D. Mass-balance model reconstruction of the Strait of
Georgia: the present, one hundred, and five hundred years ago. UBC Fisheries Centre Research
Reports 1998;6:72–91.
[36] Townsley P. Social issues in fisheries. FAO Fisheries Technical Report 375, 1998. 93pp.
[37] Garcia SM. The precautionary approach to fisheries, its implications for fishery research, technology,
and management: an updated review. FAO Fisheries Technical Paper No 350, Part 2, 1996. pp. 1–76.
[38] Bell ED, Stefansson G. Performance of some harvest control rules. NAFO Scientific Council
Research Document 98/7, 1998. 19pp.
[39] ICES. Report of the ICES advisory committee on fisheries management, 2001; ICES Cooperative
Research Report 246(3 volumes), 2001.
[40] Helstrom CW. Statistical theory of signal detection. Oxford: Pergamon Press, 1968. 470pp.
[41] Hall SJ. The effects of fishing on marine ecosystems and communities. Oxford: Blackwell Science,
1999.
[42] Rice JC. Changes to the large marine ecosystem of the Newfoundland—Labrador shelf. In: Sherman
KS, Skoldal HR, editors. Large marine ecosystems of the North Altantic. Oxford: Blackwell Scientific
Press, in press.
[43] Rose GA, Kulka DW. Hyperaggregation of fish and fisheries: how catch-per-unit-effort increased as
the northern cod (Gadus morhua) declined. Canadian Journal Fisheries and Aquatic Sciences.
1999;56(Su):118–27.
[44] Anderson PJ, Piatt JF. Community reorganization in the Gulf of Alaska following ocean climate
regime shift. Marine Ecology Progress Series 1999;189:117–23.
[45] Daan N. A spatial, temporal diversity index taking into account species rarity, with an application to
the North Sea fish community. ICES CM 2001/T:04, 2001. 6pp.
[46] Rice JC, Gislason H. Patterns of change in the size spectra of numbers and diversity of the North Sea
fish assemblage, as reflected in surveys and models. ICES Journal of Marine Science 1996;53:1214–25;
ICES. Report of the working group on ecosystem effects of fishing. ICES CM 1996/Assess/Env:1,
1996.
[47] Baillie J, Groomsbridge B. 1996 IUCN red list of threatened animals. Cambridge UK: IUCN
Publishing, 1996.
[48] Mace GM, Collar N, Cook J, Gaston K, Ginsberg J, Leader-Williams N, Maunder M, Miller-
Gulland EJ. The development of the new criteria for listing species on the IUCN Redlist. Species
1992;19:16–22;
Musick JA. Endangered marine fishes: criteria and identification of North American stocks at risk.
‘‘Fisheries’’ (Bethesda) 1998;23(2):28–30;
FAO. An appraisal of the suitability of the CITES criteria for listing commercially exploited aquatic
species. FAO Fisheries Circular, vol. 954, 2000. 66pp.
[49] Whitaker RH. Gradient analysis of vegetation. Biological Reviews 1967;49:207–64;
Gauch HG. Multivariate analysis in community ecology. Cambridge, UK: Cambridge University
Press, 1982. 298pp;
Clarke KR, Ainsworth M. A method of linking multivariate community structure to environmental
variables. Marine Ecology Progress Series 1993;92:205–19.
J. Rice / Ocean & Coastal Management 46 (2003) 235–259 259
[50] Dickie LM, Kerr SR, Boudreau PR. Size-dependent processes underlying regularities in ecosystem
structure. Ecological Monographs 1987;57:233–50.
[51] Daan N, Richardson K. Changes in the North Sea Ecosystem and their causes: Aarhus 1875 revisited.
ICES Journal of Marine Science 1996;53:875–1225.
[52] Drinkwater KF, Myers RA. Testing predictions of marine fish and shellfish landings from
environmental variables. Canadian Journal of Fisheries and Aquatic Science. 1987;44:1568–73.