Below are the first 10 and last 10 pages of uncorrected machine-read text (when available) of this chapter, followed by the top 30 algorithmically extracted key phrases from the chapter as a whole.
Intended to provide our own search engines and external engines with highly rich, chapter-representative searchable text on the opening pages of each chapter.
Because it is UNCORRECTED material, please consider the following text as a useful but insufficient proxy for the authoritative book pages.
Do not use for reproduction, copying, pasting, or reading; exclusively for search engines.
OCR for page 72
4
Relationship between Propagule Pressure
and Establishment Risk
Chapter 3 highlighted the many factors that influence invasion risk, namely
(in broad terms) propagule pressure, species traits, abiotic environmental charac-
teristics, and biotic interactions. Managing invasion risk by setting a discharge
standard assumes that, despite these powerful modifying factors, organism den-
sity alone is a reasonable predictor of establishment probability. Consequently,
this chapter examines the relationship between organism density and invasion
risk, and considers how this relationship might help inform an organism-based
discharge standard. Subsequent chapters examine other, non-modeling-based
approaches to setting a discharge standard (Chapter 5) and evaluate the data re-
quirements and limitations in estimating the relationship between invasion risk
and propagule pressure, including uncertainty, variability over space and time,
the relative merits of historical, survey, and experimental data, and the use of
proxy variables (Chapter 6).
THE RISK–RELEASE RELATIONSHIP
Concepts and Terms
There are many definitions of invasion risk (e.g., Drake and Jerde, 2009).
This report uses the term invasion risk interchangeably with establishment prob-
ability to refer to the chance that an introduced group of individuals establishes a
self-maintaining population. In formal risk assessment, risk is defined a func-
tion of both exposure (the probability of a harmful event) and hazard (the effect
of the harmful event). In this framework, invasion risk could be defined as a
function of the probability of a species establishing (exposure) and its expected
impact (hazard); however, such terminology is not used here. Rather, the term
invasion risk is defined simply as the probability of establishment. For repeated
introductions of invasive species, it is important to consider the time scale for
establishment probability because, over infinite time, any invasion with a non-
zero probability will eventually occur. The term invasion rate refers to the
72
OCR for page 72
Relationship Between Propagule Pressure and Establishment Risk 73
number of nonindigenous species that establish in a given region per unit time.
It is straightforward to convert invasion rate to establishment probability.
A propagule is any biological material (such as particles, cells, spores, eggs,
larvae, and mature organisms) that is or may become a mature organism. Prop-
agule pressure is a general term expressing the quantity, quality, and frequency
with which propagules are introduced to a given location. As discussed in
Chapter 3, propagule pressure is a function of a suite of variables reflecting the
nature of the species and the transport vector. The remainder of this chapter,
however, focuses on the quantity of propagules alone.
In the context of ballast water, it is useful to distinguish two measures of
propagule quantity. Following Minton et al. (2005), inoculum density is defined
as the density of organisms in released ballast water. Inoculum density (denoted
as DI in equations) is given simply as the total number of organisms in the in-
oculum (NI) divided by the inoculum, or ballast water, volume (VI):
DI = NI / VI.
The initial population size is the initial number of organisms released into
the environment in a given location at a given time, i.e., the inoculum number
(NI). As these organisms will tend to spread out into their new habitat, their
density in the environment (DE) is, in the simplest formulation, given as the
number of organisms released (NI) divided by the volume of water in the envi-
ronment (VE):
DE = NI / VE.
It is important to recognize that it is the inoculum density (DI) that is subject
to a ballast water discharge standard. However, both the initial number of or-
ganisms NI (conventionally denoted N0 in population modeling) and their densi-
ty in the new environment DE are expected to affect establishment probability.
In other words, the variable to be managed to reduce invasion risk is clearly dis-
tinct from the variable that is typically used in predicting invasion risk. This
disconnect is the central conceptual challenge in converting empirical and theo-
retical results in population establishment to an operational discharge standard
(see section below).
The relationship between invasion risk and propagule pressure is the risk–
release relationship. Understanding the risk–release relationship is essential to
predicting and comparing the invasion risk associated with different discharge
standards. However, understanding this relationship is not a straightforward
proposition. It is easiest to define, model, and estimate this relationship for sin-
gle species, focusing on the relationship between the number of individuals re-
leased at a given time in a given location, and the probability of that population
establishing. For larger-scale analyses of invasions by multiple species, the rela-
tionship becomes less easy to define, model, and measure. In general, however,
it can be thought of as the number of established species as a function of the
OCR for page 72
74 Propagule Pressure and Invasion Risk in Ballast Water
number of released species, organisms, or combination thereof, over a given
time period. In the case of ballast water, the central but somewhat indirect risk–
release relationship is the number of species that establish as a function of the
large-scale release of a varying number of varying species at varying densities.
The Hypothetical Risk–Release Relationship
In general, the relationship between invasion risk and propagule pressure is
expected to be positive, although its shape is unknown. A priori, it might take
any of a number of standard shapes including linear, exponential, hyperbolic,
and sigmoid (Ruiz and Carlton, 2003) as shown in Figure 4-1. The shape of the
risk–release relationship has important implications for managing invasion risk.
If the relationship were linear, then a given reduction in release density
would always lead to a proportional reduction in invasion risk (Figure 4-2A). If
the relationship were nonlinear, with one or more inflection or slope-balance
points, then interesting management thresholds would emerge (Figure 4-2B). If
the relationship were exponential, a reduction from high to moderate release
density leads to a much greater reduction in invasion risk than a similar-sized
reduction from moderate to low release density. The opposite would hold for
FIGURE 4‐1 Common shapes for relationships between two variables: (a) hyperbolic, (b)
sigmoid, (c) linear, (d) exponential (both axes linear). A priori, any of these could
represent the relationship between invasion risk (probability of a species establishing)
and propagule pressure (e.g., number of individuals released). SOURCE: Adapted, with
permission, from Ruiz and Carlton (2003). ©2003 by Island Press.
OCR for page 72
Relationship Between Propagule Pressure and Establishment Risk 75
A B
FIGURE 4‐2 Conceptual application of a risk–release relationship to inform ballast water
organism discharge standards. S is the observed risk of ballast‐mediated species inva‐
sion; N is the observed number of organisms released. (Axis units depend on whether
the model represents a single species or multiple species.) S* is the target invasion risk;
N* is the corresponding target release value. Assuming a robust risk–release relation‐
ship, reducing the ballast water release by the proportion RN is predicted to reduce the
invasion risk by the proportion RS. (A) Under the assumption of a linear risk‐release
relationship, a given reduction in the release rate is predicted to give the same propor‐
tional reduction in invasion rate (i.e., RN=RS). (B) Under the assumption of a sigmoid
relationship, the same reduction in release (RN1=1‐[N*1/N1] = RN2=1‐[N*2/N2]) is pre‐
dicted to give a much lesser (RS1 =1‐[S*1/S1]) or a much greater (RS2 =1‐[S*2/S2]) reduc‐
tion in invasion risk, depending on the range over which RN occurs. Solid dot indicates
the inflection point around which the greatest reduction in risk is obtained for the least
reduction in release. Open dots indicate the points at which the slope passes through a
o
45 angle: outside these bounds, increasingly less risk reduction is obtained for the same
release reduction. Panel (A) provides a graphical illustration of the multi‐species linear
model proposed by Cohen (2005, 2010) and Reusser (2010), where the y‐axis is the in‐
vasion risk (characterized as invasion rate, the number of species established per unit
time) and the x‐axis is organism release (characterized as release rate, the total number
organisms per unit time). Reusser (2010) defined as the per capita invasion risk the
number of introduced species that establish per organism released, i.e., the slope of the
straight‐line relationship. Cohen (2005, 2010) defined as the reduction factor the equal
proportional reductions in release and risk from current to target levels.
OCR for page 72
76 Propagule Pressure and Invasion Risk in Ballast Water
a hyperbolic relationship. For an S-shaped curve, a reduction over the middle
range of release densities would to the greatest reduction in invasion risk (Figure
4-2B). Because of qualitative change in the risk–release relationship at the in-
flection point and at the slope-balance points (where the tangent to the curve
passes through 45o), non-linear relationships present influential management
thresholds.
Theory tells us that for a single population, only two of these four shapes—
the hyperbolic and sigmoid curves—can represent the overall risk–release rela-
tionship (see Box 4-1). The combined curves of multiple populations would
also be expected to be nonlinear. Thus, it is expected that there should be at
least one threshold in the risk–release relationship that could in principle prove
useful in informing discharge standards.
Despite understanding that the overall theoretical shape is hyperbolic or
sigmoid for a single species, it is possible for a given set of risk–release data to
be better characterized by a linear or even an altogether different model. This
apparent discrepancy could emerge for two main reasons. First, there may be
insufficient data points to support a curved line over a straight line. This diffi-
culty will be exacerbated as the true slope decreases (for example, the lower-left
or upper-right ends of the hyperbolic or sigmoid curves). Second, any underly-
ing theoretical risk–release relationship may be swamped out by other more im-
portant sources of variation that affect establishment probability (see Chapter 3),
such that it cannot be recovered from the data.
To quantitatively predict the effects of a discharge standard on invasion
risk, and to compare the risk associated with different discharge standards, it is
essential to understand the shape and strength of the risk–release relationship.
The following section reviews a range of approaches that have been taken to
fitting risk–release curves to empirical data.
Modeling the Risk–Release Relationship
An ideal analysis of the risk–release relationship would involve developing
and testing a suite of candidate theoretical models, collecting multiple rigorous
empirical datasets, and comparing the fit of the models to the data to determine
(1) which model best captures the risk–release relationship and (2) how strong
this relationship is relative to other potential explanatory variables.
Several approaches have been taken to modeling the risk–release relation-
ship, categorized in Table 4-1 along two axes. First, models can range from
descriptive models that simply represent the shape of the relationship to mecha-
nistic models that define the processes generating the relationship (e.g., Drake
and Jerde, 2009). These are not mutually exclusive categories—a given model
may include both mechanistic and descriptive components—but at their ex-
tremes they represent very different modeling philosophies and goals, and they
define a useful spectrum for organizing modeling approaches. Second, Table 4-
1 distinguishes models that focus on the establishment of a single species vs.
OCR for page 72
Relationship Between Propagule Pressure and Establishment Risk 77
TABLE 4‐1 Approaches Taken to Modeling the Relationship Between Invasion Risk and
Propagule Pressure
Scale
Sample structure1
Type Single species Multiple species
examples examples
Descriptive Statistical Memmott et al. Lonsdale (1999)
Logistic regression: (2005) Levine and D’Antonio
y Jongejans et al. (2003)
b0 b1 x (2007)
ln
1 y
Ricciardi (2006)
Bertolino (2009) Cohen (2005, 2010)*
Functional form Drake and Jerde Reusser (2010)*
Species‐area curve: (2009)
z
y=cx
Mechanistic Probabilistic Shea and Possing‐ Drake et al. (2005)*
N
pE 1 e ham (2000) Costello et al. (2007)*
Leung et al. (2004) USCG (2008)*
Dynamic demographic Memmott et al. See text
dN (2005)
rN
dt Drake et al. (2006)
USCG (2008)*
Drake and Jerde
(2009)
Jerde et al. (2009)
Bailey et al. (2009)
Kramer and Drake
(2010)
Notes: Models are categorized as being on the descriptive or the mechanistic end of a spectrum,
and as representing single or multiple invading species. Sample structures show simple, generic
forms of these model types. Most of these approaches have been widely developed and imple‐
mented throughout the biological literature; only a few recent examples, further discussed in the
text, are listed here. Some studies illustrate more than one modeling approach.
1
Parameters: y, dependent variable (invasion risk); x, independent variable (organism discharge or
proxy variable); b0, b1, c, z, shape parameters; pE, population establishment probability; α,
ln(individual establishment probability); N, number of individuals; r, per capita population growth
rate.
* Proposed application of modeling approach to ballast water management.
OCR for page 72
78 Propagule Pressure and Invasion Risk in Ballast Water
that of multiple species. The following sections highlight illustrative examples
of each approach and outline their key advantages and disadvantages. Since
population establishment theory applies across all species, habitats, and vectors,
many of the given examples fall outside the immediate realm of ballast-
mediated invasions. Nevertheless, the approaches illustrated are readily appli-
cable to the risk–release relationship for ballast water.
SINGLE-SPECIES MODELS
It is informative to examine the risk–release relationship at the scale of a
single species, for two main reasons. First, this approach allows examination of
invasion scenarios for certain model species, such as fast growing, high impact,
or commonly released invaders, which could be used to obtain upper bounds for
discharge standards under best-case (for invasion) scenarios. Second, it allows
for clarification of model structures and assumptions before scaling up to the
more realistic scenario of multi-species releases. The primary disadvantage of
the single-species approach, in the context of managing ballast water, is that it
does not represent the reality of the simultaneous and continuous release of
many species from ballast water. The greatest challenge in this approach is in
converting experimental and theoretical results premised on N0 to a discharge
standard applicable to DI.
A ballast-mediated invasion may be expected to begin from the introduction
of relatively few individuals. Three factors are particularly relevant to small
population dynamics—demographic stochasticity, positive density dependence,
and the spatial environment—and must be considered in developing an effective
risk–release model. Their net effect can be captured by a descriptive model;
their individual effects can be tested by incorporating them explicitly into a me-
chanistic model.
Demographic stochasticity is the natural variability in individual survival
and reproduction that occurs in populations of any size, but that in small popula-
tions can lead to large fluctuations in population growth rate. On average, de-
mographic stochasticity makes extinction more likely than the equivalent deter-
ministic model would predict; on the other hand, in a given realization, stochas-
ticity can lead to establishment when a deterministic model would otherwise
predict extinction (Morris and Doak, 2002; Drake, 2004; Andersen, 2005; Drake
et al., 2006; Jerde et al., 2009). Demographic stochasticity is expected to lower
the risk–release curve at low density.
Positive density dependence, or Allee effects, is the intuitively logical notion
that as organism density decreases, individuals may suffer increasing difficulty
finding mates or foraging. Population growth rate would thus be expected de-
cline at low density, rather than to increase as per an exponential growth model
(Allee, 1931; Courchamp et al., 2009). Thus, across the range of low initial
densities, it takes a higher density to achieve the same invasion risk when Allee
effects are at work than when they are not. Allee effects lower the risk–release
OCR for page 72
Relationship Between Propagule Pressure and Establishment Risk 79
curve at low density, leading to a characteristic sigmoid curve (see Curve B in
Figure 4-1). Allee effects are logically appealing and have been applied to mod-
els of both sexually reproducing and parthenogenetic organisms (Drake, 2004).
However, they have proved elusive to document empirically (see discussion and
examples in Dennis, 2002; Morris and Doak, 2002; Leung et al., 2004; Drake et
al., 2006; Courchamp et al., 2009; Jerde et al., 2009; Kramer and Drake, 2010).
The third feature is the effect of the spatial environment on population den-
sity. One of the greatest challenges in population modeling is that organisms
released into an environment will tend to spread out, through both passive and
active dispersal. Depending on the relative scales of dispersal and population
growth, the effective initial population size may be very much lower than the
original number of organisms released. This reduction in density will presuma-
bly exacerbate the impacts of demographic stochasticity and Allee effects.
Thus, in general it is expected that both individual and population establishment
probabilities to be much lower in the wild than in contained laboratory experi-
ments. (It is possible, of course, that hydrodynamic features or aggregative be-
havior will have the opposite effect, tending to concentrate organisms in a locale
and enhancing their chances of establishment; see Chapter 3).
Since establishment probability seems generally likely to be dominated by
the dynamics of small populations, the most rigorous modeling approach is to
construct models that incorporate (or can phenomenologically reflect the effects
of) demographic stochasticity, Allee effects, and their modification by dispersal,
and to allow the empirical data to indicate on a case-by-case basis the impor-
tance of these additional features. The following examples consider none, some,
or all of these features.
Descriptive Models
Descriptive models, such as regression and similar statistical techniques, of-
fer a phenomenological characterization of the risk–release relationship. That is,
they can be formulated and parameterized without having to understand or spe-
cify the underlying mechanisms by which the independent variables explain the
dependent variable (Drake and Jerde, 2009). As a result, one’s confidence in
their predictive ability is limited. Familiar descriptive models include statistical
models such as regression and functional forms such as species-area curves and
behavioral responses.
These models have been applied to the results of both experimental and his-
torical survey data. In a simple example, Drake and Jerde (2009) fit a spline, or
a series of local regressions, to establishment probability as a function of propa-
gule pressure in the scentless chamomile (Matricaria perforata) (Figure 4-3A).
In this case, establishment was defined as survival simply from seed to flower-
ing, but the same method could be applied to a longer-term study of population
establishment. The same data were also fit with a probabilistic model (see later
section). In a slightly more complex field experiment, Jongejans et al. (2007)
OCR for page 72
80 Propagule Pressure and Invasion Risk in Ballast Water
estimated the establishment probability (defined as persistence over six years) of
the European thistle (Carduus acanthoides) as a function of propagule pressure
and native plant biomass. Using a generalized linear model, they found that
establishment probability increased significantly with higher initial seed num-
ber, and tended to increase with reduced native biomass; together these variables
accounted for 37 percent of the variation (Figure 4-3B). Bertolino (2009) mod-
eled the success of global historical squirrel introductions (defined as persistence
to the present day of populations introduced over a >100-year period) as a func-
tion of propagule pressure, environmental matching, native diversity, and the
invaders’ biogeographical origin. For the genus Sciurus (squirrels), a logistic
regression fit to the initial number of individuals alone explained 55 percent of
the variation in establishment probability (Figure 4-3C).
Statistical models can be made increasingly complex by adding ever more
independent variables, and have been used to describe invasion risk over a spa-
tial domain using species distribution (environmental niche) modeling (Peterson
and Vieglais, 2001; Herborg et al., 2007, 2009; Dullinger et al., 2009).
Mechanistic Models
In contrast to descriptive statistical models, mechanistic models represent
invasion establishment as a function of parameters that have a readily defined
biological meaning. Whereas statistical models describe a relationship only
over the range of data to which they are fit, mechanistic models are presumed to
extrapolate well over the entire biologically realistic parameter space. Further-
more, descriptive models allow one to investigate the shape of a relationship,
while mechanistic models force the user to specify the processes driving the
relationship and to link causative variables explicitly.
It is useful to distinguish two general classes of mechanistic models: proba-
bilistic statements and dynamic, demographic models. In their simpler forms,
these two model classes possess different mathematical structures and require
different data to parameterize and validate. In more complex models of popula-
tion establishment, this distinction blurs and a given model may incorporate
elements of both classes (e.g., Jerde and Lewis, 2007; Leung and Mandrak,
2007; Jerde et al., 2009).
Probabilistic Models
In the context of the risk–release relationship, the probabilistic models con-
sidered here are composed of probability statements beginning with the proba-
bility of an individual’s establishment probability and scaling up to a population
level. A probabilistic model is written immediately in terms of its solution,
namely, in terms of population establishment probability. When the model’s
constituent probabilities are represented as frequency distributions rather than as
OCR for page 72
Relationship Between Propagule Pressure and Establishment Risk 81
A
B
C
FIGURE 4‐3 Single‐species risk–release relationships obtained from descriptive models.
(A) Spline fit to short‐term establishment probability of scentless chamomile (Matricaria
perforata) (Drake and Jerde, 2009). (B) Generalized linear model fit to establishment
outcome (success or failure) for European thistle (Carduus acanthoides) across three
levels of native plant biomass reduced by clipping (Jongejans et al., 2007). (C) Logistic
regression fit to establishment probability of squirrel Sciurus. Replotted from Bertolino
(2009) with data generously provide by S. Bertolino. SOURCES: (A) Reprinted, with per‐
mission from, Drake and Jerde (2009). © 2009 by Oxford University Press. (B) Re‐
printed, with permission, from Jongejans et al. (2009). © by Springer. (C) Reprinted,
with permission, from Bertolino (2009). © 2009 by John Wiley and Sons.
OCR for page 72
82 Propagule Pressure and Invasion Risk in Ballast Water
fixed points, it is known as a hierarchical probability model (HPM). The HPM
approach to representing stochastic events has been extensively applied to medi-
cal, engineering, and ecological problems, and it allows the explicit incorpora-
tion and analysis of uncertainty (Dennis and Lele, 2009; Ponciano et al., 2009).
Parameterizing and validating probabilistic models requires comparatively sim-
ple data: the outcome, as success or failure, of a series of introductions inocu-
lated over a range of initial organism numbers. As will be discussed later, prob-
abilistic models can readily be expanded to represent multiple species and envi-
ronmental conditions that have different associated probabilities of establish-
ment.
Probabilistic models of population establishment have been developed to
serve as the basis for a metapopulation model of biocontrol release (Shea and
Possingham, 2000) and for a gravity model of zebra mussel spread (Leung et al.,
2004). In their simplest, non-spatial form, they contain a sole parameter—the
probability of a single individual producing an established population. This
value is then scaled up to obtain the probability of a group of individuals leading
to an established population (Shea and Possingham, 2000; Leung et al., 2004,
Leung and Mandrak, 2007; Jerde and Lewis, 2007; Jerde et al., 2009). The dif-
fusion approximation to exponential growth shares the same core probability
structure, and has likewise been used to model population establishment (Drake
et al., 2006; Bailey et al., 2009).
The basic construction of a simple probabilistic model is outlined in Box 4-
1, and its application is illustrated in several examples below.
To implement these probabilistic models, studies have examined population
establishment over a range of initial population sizes either from directed expe-
riments or from descriptive population data. Memmott et al. (2005) fit both a
logistic regression and a probabilistic model to the success after six years of
biocontrol insect releases (Arytainilla spartiophila; Figure 4-5A). Drake and
Jerde (2009) used short-term data for the success of the prairie weed scentless
chamomile from seed to flowering (Matricaria perforata; Jerde and Lewis,
2007) to scale up to a population-level establishment model. This model was
subsequently extended by Jerde et al. (2009) to incorporate mate-finding limita-
tions that led to a biologically driven Allee effect, and it was used to predict
invasion risk for Chinese mitten crabs (Eriocheir sinensis) and apple snails
(Pomacea canaliculata) (Jerde et al., 2009; Figure 4-5B). Drake et al. (2006)
used descriptive population growth data of the spiny waterflea Bythotrephes
longimanus in three lakes over four years, and Bailey et al. (2009) conducted
100-day mesocosm studies of a variety of cladocerans, to parameterize diffusion
approximations to exponential growth models (Figure 4-5C,D).
Probabilistic models have also been constructed to investigate the accumu-
lation of nonindigenous species over time (Solow and Costello, 2004; Wonham
and Pachepsky, 2005); this type of model has been extended to examine the rela-
tionship between propagule pressure and invasion risk (Costello et al., 2007; see
later section). Basic probabilistic models can represent the invasion outcome
alone, without necessarily representing the population dynamics leading to
OCR for page 72
Relationship Between Propagule Pressure and Establishment Risk 103
BOX 4‐4
Outline of a Simple Multispecies Probabilistic Model
The single‐species probabilistic model developed in Box 4‐1 is readily ex‐
tended to a multispecies probabilistic model. Equation 4‐3 in Box 4‐1 defines
the establishment probability for a single species as pE = 1-e-N. This equation
can be modified to accommodate S species, each with its own establishment
probabilty ps. Following the same general approaches described in Shea and
Possingham (2000) and USCG (2008), the model then describes SE, the ex‐
pected number of species that establish, as
S
SE 1 e i Ni
ci
s 1 (4-8)
and ps, the probability that at least one species establishes, as
s
pS 1 ei Ni (4-9)
ci
i 1
The same principles can be used to extend the model to consider variation
in establishment probability across multiple locations (L), multiple vessels arriv‐
ing in those locations (VL), and so on. In this case, the expected number of spe‐
cies that establish can be written as
VL
S L
SE 1 e s,l ,v Ns,l ,v
cs , l , v
l 1 v 1
s 1 (4-10)
and the probability that at least one species establishes as
VL
S L
pS 1 e s,l ,v Ns,l ,v
cs ,l , v
s 1 l 1 v 1 (4-11)
box continues next page
OCR for page 72
104 Propagule Pressure and Invasion Risk in Ballast Water
BOX 4‐4 Continued
These models follow the same logic as the basic one (4‐3), where the
shape parameters and c describe the species establishment probability as a
function of the initial number of individuals N; the subscripts s, l, and v allow
variation among species, locations, and vessels; and the complement of all the
propagules failing to establish gives the final probability of establishment. This
model does not account for variation in parameter values over time or for po‐
tential interactions among species. Nevertheless, parameterizing such a mod‐
el, particularly in a hierarchical structure where each parameter is characte‐
rized as a frequency distribution, would require a tremendous amount of data.
Qualitatively, the overall establishment probability ps obtained from a mul‐
ti‐species model can only be the same as, or greater than, the largest estab‐
lishment probability pE of the constituent species.
To the Committee’s knowledge, this approach has not yet been applied to
predicting the success of multiple nonindigenous species. However, a related
approach using a reaction-diffusion model, which is a standard spatial extension
of a demographic model, offers interesting insights. Reaction-Diffusion (R-D)
models represent a class of models that were originally developed to model the
spread of organisms in continuous time and space (Skellam, 1951). These mod-
els were later developed to model the spread of invading organisms across a
one- or two-dimensional landscape with the goal of defining the rate of spread
as a travelling wave and so provide a description of the rate of spread and the
area occupied by the invasion (Okubo et al., 1989). The classical version of
these models typically involves several restrictive assumptions including spatial
homogeneity and random movement (at least at the population level). However,
they are comparatively easy to parameterize, requiring only estimates of per
capita rate of population increase and the mean squared displacement per unit
time of individuals in the population. Among the advantages of this approach
include the ability to approximate the spread of the invading population as a
travelling wave (Okubo et al., 1989). This permits the estimation of the rate of
spread as a linear function of time, so that the arrival of an invader at a new site
could be reasonably estimated under the model assumptions.
These models typically do not provide an estimate of the rate of establish-
ment, although recent applications have attempted this for multiple species using
some simplifying assumptions (Drake et al., 2005). In Drake et al. (2005), the
authors combined a generic exponential-growth reaction-diffusion model with
an allometric relationship between body size and population growth rate to ex-
OCR for page 72
Relationship Between Propagule Pressure and Establishment Risk 105
amine the establishment probability of a variety of aquatic species. Their goal
was not to predict establishment probabilities of any particular species, but to
predict invasion rates over a range of species. Under a series of assumptions,
particularly concerning Allee dynamics, this model estimates the risk–release
relationship in terms of the proportion of species of a given body size that estab-
lish vs. the volume of water released. From this output, the chance of a single
invasion by a size class of organism can be predicted as a function of the num-
ber and volume of releases, independent of the number of individual organisms.
This proposed approach is specific to a size class, and has not yet been vali-
dated, but seems reasonable.
Obtaining a Discharge Standard from Multiple-Species Models
Multi-species models represent an attempt to capture the complexity of
wholesale ballast water release. Descriptive and mechanistic models can readily
be formulated in conceptual and mathematical terms at this scale. In the current
absence of data for parameterizing and validating mechanistic models, descrip-
tive statistical models can be developed. However, these must be interpreted
with caution given the uncertainty in the estimates, and disconnect between the
scales of the independent and dependent variables. The use of proxy variables
introduces a further challenge: although well-fitting models may be obtained,
proxy variables (such as ballast volume, shipping tonnage, vessel abundance)
must be evaluated for their relationship to the direct variable of interest (dis-
charge density) before ascribing any mechanistic meaning to their relationships.
To summarize, a multi-species approach focuses on the assemblage of spe-
cies released from ballast water. Because of the associated data requirements,
descriptive models are more likely than mechanistic models to yield estimates of
the risk–release relationship at this scale. Even so, given the uncertainty and
mismatch in both the independent and dependent variables, the applicability of
any apparent relationships is questionable. Relative to the single-species ap-
proach, a multi-species approach has the advantage of being conceptually more
realistic in the context of ballast water release, and the disadvantage of being
more complex and more difficult to ground in the relevant empirical data.
CONCLUSIONS AND RECOMMENDATIONS
Models are generally useful in environmental management because they
provide a transparent framework, force an explicit statement of assumptions,
allow us to predict and compare future projections under different management
scenarios, and can be updated in their structure and parameter estimates as new
information emerges. In principle, a well-supported model of the relationship
between invasion risk and organism release could be used to inform a ballast
water discharge standard. For a given discharge standard, the corresponding
OCR for page 72
106 Propagule Pressure and Invasion Risk in Ballast Water
invasion risk could be predicted, or, for a given target invasion risk, the corres-
ponding target release level could be obtained.
Candidate risk–release models developed to date include single- and mul-
tiple-species scales, and extend along the spectrum from descriptive to mecha-
nistic in their construction. Mechanistic single-species models require fewer
data to parameterize than do mechanistic multi-species models, but do not
represent the more realistic scenario of ballast discharge of an assemblage of
species. Descriptive single-species models are simpler, but offer none of the
predictive advantages of mechanistic ones. Descriptive multi-species models
are an appealing tool for investigating large correlative datasets, but are ham-
pered by a current lack of appropriate data.
The rigorous use of models requires that multiple candidate models be for-
mulated and compared in their ability to represent the data. This approach is
well established in the population dynamic literature at the single-species scale.
However, currently there are insufficient data to distinguish among risk–release
relationship models at the multi-species scale. The following conclusions and
recommendations identify how models might be put to use at present, and in the
future, to help inform a discharge standard.
Ballast water discharge standards should be based on models, and be
explicitly expressed in an adaptive framework to allow the models to be
updated in the future with new information. Before being applied, it is es-
sential that candidate models be tested and compared, and their compounded
uncertainty be explicitly analyzed. Only a handful of quantitative analyses of
invasion risk–release relationships thus far have tested multiple models and
quantified uncertainty.
The predicted shape of the risk–release relationship is non-linear. In-
flection points and slope-balance points could provide natural breakpoints for
informing a discharge standard. However, the apparent shape of the relation-
ship for a given system will depend on the quantity, error, and parameter range
of the empirical data, as well as the biology of the species and the nature of the
environment.
In the short term, mechanistic single-species models are recommended
to examine risk–release relationships for best-case (for invasion) scenario
species. This approach makes sense biologically because in general concerns
are only about the small subset of released species that establish as high-impact
invaders. Such an approach to setting a standard is conservative and would pro-
vide maximum safety against invader establishment.
Candidate best-case-scenario species should be those with life histories that
would favor establishment with the smallest inoculum density. Species with the
highest probability of establishment relative to inoculum density will have the
greatest influence in determining the shape of the risk–release curve. Life histo-
ry traits promoting such sensitivity to small inoculum density possibly include
OCR for page 72
Relationship Between Propagule Pressure and Establishment Risk 107
fast-growth, parthenogenetic or other asexual reproductive abilities, lecitho-
trophic larvae, etc. Other considerations of best-case species might include
those that have a high ecological or economic impact, or are frequently intro-
duced. The greatest challenge in this approach will be converting the re-
sults of small-scale studies to an operational discharge standard.
Developing a mechanistic multi-species model of risk and release, parame-
terized for an assemblage of best-case scenario species, would only be recom-
mended over the longer term. This model would allow a detailed theoretical
investigation of the relationship between total organism number and invasion
risk, by permitting the analysis of the risk associated with different species rich-
ness and frequency distributions that sum to the same total organism number.
The challenges in this approach include the massive time and effort needed to
gather the necessary data as well as converting model results to a fleet-wide dis-
charge standard. The implications of these models would therefore be highly
specific, no more (and possibly less) informative than those of single-species
models, and will require more data and computational effort to construct, para-
meterize, and validate.
Developing a robust statistical model of the risk–release relationship is
recommended. It is unclear whether the current lack of a clear pattern across
ports reflects a true absence of pattern, or the absence of appropriate data to test
this model. Nevertheless, given spatial variation in shipping patterns and envi-
ronmental variables, it is anticipated that this approach will be more fruitful at a
local scale than a nation-wide scale.
Within a region, this relationship should be estimated across multiple time
intervals, rather than from a single point. The effect of temporal bin sizes on the
shape of the relationship must be examined. The choice of independent variable
must be carefully considered. Since long-term historical data on ballast-
organism density do not exist, the committee recommends an extremely careful
analysis and validation of any proxy variables. The greatest challenge in this
approach is the currently insufficient scope and scale of the data. There is no
evidence that any proxy variable used thus far is a reliable stand-in for or-
ganism density.
Finally, models of any kind are only as informative as their input data.
In the case of ballast water, both invasion risk and organism density dis-
charged from ballast water are characterized by considerable and largely
unquantified, uncertainty. At the multi-species scale in particular, the existing
data (historical invasion records vs. recent ballast surveys) are substantially
mismatched in time, and patchy in time, space, and taxonomy; current statistical
relationships with these or proxy variables are of dubious value. The judicious
use of an appropriate model combined with robust data may help inform stan-
dard setting in the future.
OCR for page 72
108 Propagule Pressure and Invasion Risk in Ballast Water
REFERENCES
Allee, W. 1931. Animal Aggregations. A Study in General Sociology. Chicago: Uni-
versity of Chicago Press.
Andersen, M. 2005. Potential applications of population viability analysis to risk as-
sessments for invasive species. Human and Ecological Risk Assessment 11:1083–
1095.
Bailey, S., L. Velez-Espino, O. Johannsson, M. Koops, and C. Wiley. 2009. Estimating
establishment probabilities of Cladocera introduced at low density: an evaluation of
the proposed ballast water discharge standards. Canadian Journal of Fisheries and
Aquatic Sciences 66:261–276.
Barry, K. L., and C. D. Levings. 2002. Feasibility of using the RAMAS metapopulation
model to assess the risk of a non-indigenous copepod (Pseudodiaptomus marinus)
establishing in Vancouver Harbour from ballast water. Canadian Technical Reports
Fisheries Aquatic Science 2401.
Bartell, S., and S. Nair. 2004. The establishment of invasive species: an interface be-
tween risk analysis and theoretical population ecology. Risk Analysis 20:833–846.
Bertolino, S. 2009. Animal trade and non-indigenous species introduction: the world-
wide spread of squirrels. Diversity and Distributions 15:701–708.
Cameron, E., and E. Bayne. 2009. Road age and its importance in earthworm invasion
of northern boreal forests. Journal of Applied Ecology 46:28–36.
Carlton, J., and J. Geller. 1993. Ecological roulette: The global transport of nonindigen-
ous marine organisms. Science 261:78–82.
Case, T. 1990. Invasion resistance arises in strongly interacting species-rich model
competition communities. Proceedings of the National Academy of Sciences of the
United States 87:9610–9614.
Case, T. 1995. Surprising behavior from a familiar model and implications for competi-
tion theory. The American Naturalist 146:961–966.
Castro, S. A., and F. M. Jaksic. 2008. How general are global trends in biotic homoge-
nization? Floristic tracking in Chile, South America. Global Ecology and Biogeo-
graphy 17:524–531.
Caswell, H. 1989. Matrix Population Models: Construction, Analysis and Interpretation.
Sunderland, MA: Sinauer Assoc.
Chapman, J., and J. Carlton. 1991. A test of criteria for introduced species: The global
invasion by the isopod Synidotea laevidorsalis (Miers, 1881). Journal of Crustacean
Biology 11:386–400.
Chiron, F., S. Shirley, and S. Kark. 2009. Human-related processes drive the richness of
exotic birds in Europe. Proceedings of the Royal Society of London Series B
276:47–53.
Cohen, A. 2005. Memo on a natural invasion rate standard. California State Lands
Commission.
Cohen, A. N. 2010. The natural invasion rate model and California’s ballast water dis-
charge standards. Presentation to the National Academy of Science/National Rese-
arch Council Committee on Assessing Numeric Limits for Living Organisms in Ba-
llast Water Washington, DC. June 2, 2010.
Cohen, A. N., and J. T. Carlton. 1998. Accelerating invasion rate in a highly invaded
estuary. Science 279(5350):555–558.
Copp, G. H., M. Templeton, and R. E. Gozlan. 2007. Propagule pressure and the inva-
sion risks of non-native freshwater fishes: a case study in England. Journal of Fish
Biology 71:148–159, Suppl. D.
OCR for page 72
Relationship Between Propagule Pressure and Establishment Risk 109
Copp, G., L. Vilizzi, and R. Gozlan. 2010. The demography of introduction pathways,
propagule pressure and occurrences of non-native freshwater fish in England. Aqua-
tic Conservation: Marine and Freshwater Ecosystems 20:595–601.
Cordell, J., D. Lawrence, N. Ferm, L. Tear, S. Smith, and R. Herwig. 2009. Factors
influencing densities of non-indigenous species in the ballast water of ships arriving
at ports in Puget Sound, Washington, United States. Aquatic Conservation: Marine
and Freshwater Ecosystems 19:322–343.
Costello, C., and A. Solow. 2003. On the pattern of discovery of introduced species.
Proceedings of the National Academy of Sciences 100:3321–3323.
Courchamp, F., L. Berec, and J. Gascoigne. 2009. Allee Effects. In: Ecology and Con-
servation. Oxford: Oxford University Press.
Dawson, W., D. Burslem, and P. Hulme. 2009. Factors explaining alien plant invasion
success in a tropical ecosystem differ at each stage of invasion. Journal of Ecology
97:657–665.
Dehnen-Schmutz, K., J. Touza, C. Perrings, and M. Williamson. 2007. The horticultural
trade and ornamental plant invasions in Britain. Conservation Biology 21:224–231.
Dennis, B. 2002. Allee effects in stochastic populations. Oikos 96:389–401.
Dennis, B., and S. Lele. 2009. Bayesian methods for hierarchical models: are ecologists
making a Faustian bargain? Ecological Applications 19:581–584.
Drake, J. 2004. Allee effects and the risk of biological invasion. Risk Analysis 24:795–
802.
Drake, J., K. Drury, D. Lodge, A. Blukacz, and N. Yan. 2006. Demographic stochastici-
ty, environmental variability, and windows of invasion risk for Bythotrephes longi-
manus in North America. Biological Invasions 8:843–861.
Drake, J., and C. Jerde. 2009. Stochastic models of propagule pressure and establish-
ment. Pp. 83–102 In: Bioeconomics of Invasive Species. R. Keller, D. Lodge, M.
Lewis, and J. Shogren, editors. Oxford: Oxford University Press.
Drake, J. M., and D. M. Lodge. 2004. Global hot spots of biological invasions: evaluat-
ing options for ballast-water management. Proceedings of the Royal Society of
London Series B-Biological Sciences 271(1539):575–580.
Drake, J., D. Lodge, and M. Lewis. 2005. Theory and preliminary analysis of species
invasions from ballast water: controlling discharge volume and location. American
Midland Naturalist 154:459–470.
Duggan, I., C. van Overdijk, S. Bailey, P. Jenkins, H. Limén, and H. MacIsaac. 2005.
Invertebrates associated with residual ballast water and sediments of cargo-carrying
ships entering the Great Lakes. Canadian Journal of Fisheries and Aquatic Sciences
62:2463–2474.
Dullinger, S., I. Kleinbauer, J. Peterseil, M. Smolik, and F. Essl. 2009. Niche based
distribution modelling of an invasive alien plant: effects of population status, propa-
gule pressure and invasion history. Biological Invasions 11:2401–2414.
Dunstan, P., and N. Bax. 2007. How far can marine species go? Influence of population
biology and larval movement on future range limits. Marine Ecology Progress Se-
ries 344:15–28.
Fitzpatrick, M., E. Preisser, A. Ellison, and J. Elkinton. 2009. Observer bias and the
detection of low-density populations. Ecological Applications 19:1673–1679.
Fofonoff, P. W., G. M. Ruiz, A. H. Hines, B. D. Steves, and J. T. Carlton. 2009. Four
Centuries of Estuarine Biological Invasions in the Chesapeake Bay Region. Pp.
479–506 In: Biological Invasions in Marine Ecosystems: Ecological, Management,
and Geographic Perspectives. Rilov, G. and J. Crooks (editors). Berlin: Springer-
Verlag.
OCR for page 72
110 Propagule Pressure and Invasion Risk in Ballast Water
Gotelli, N. J. 1991. Metapopulation Models: The Rescue Effect, the Propagule Rain, and
the Core-Satellite Hypothesis. American Naturalist 138(3):768–776.
Grigorovich, I. A., R. I. Colautti, E. L. Mills, K. Holeck, and H. J. MacIsaac. 2003. Bal-
last-mediated animal introductions in the Laurentian Great Lakes: retrospective and
prospective analysis. Canadian Journal of Fisheries and Aquatic Sciences 60:740–
756.
Haccou, P., and Y. Iwasa. 1996. Establishment probability in fluctuating environments:
a branching process model. Theoretical Population Biology 50:254–280.
Hastings, A., K. Cuddington, K. Davies, C. Dugaw, S. Elmendorf, A. Freestone, S. Harri-
son, M. Holland, J. Lambrinos, U. Malvadkar, B. Melbourne, K. Moore, C. Taylor,
and D. Thomson. 2005. The spatial spread of invasions: new developments in
theory and evidence. Ecology Letters 8:91–101.
Herborg, L.-M., C. Jerde, D. Lodge, G. Ruiz, and H. MacIsaac. 2007. Predicting inva-
sion risk using measures of introduction effort and environmental niche models.
Ecological Applications 17:663–674.
Herborg, L. M., P. O'Hara, and T. W. Therriault. 2009. Forecasting the potential distri-
bution of the invasive tunicate Didemnum vexillum. Journal of Applied Ecology
46(1):64–72.
Hines, A. H., and G. M. Ruiz. 2000. Biological invasions of cold-water coastal
ecosystems: Ballast-mediated introductions in Port Valdez / Prince William Sound,
Alaska. Final Report, Regional Citizens Advisory Council of Prince William
Sound.
Holyoak, M., M. Leibold, and R. Holt. 2005. Metacommunities: Spatial Dynamics and
Ecological Communities. Chicago: University of Chicago Press.
Hulme, P. 2009. Trade, transport and trouble: managing invasive species pathways in an
era of globalization. Journal of Applied Ecology 46:10–18.
Jerde, C., C. Bampfylde, and M. Lewis. 2009. Chance establishment for sexual, semel-
parous species: Overcoming the Allee effect. American Naturalist 173:734–746.
Jerde, C., and J. Bossenbroek. 2009. Uncertain invasions: a biological perspective. In:
Bioeconomics of Invasive Species. R. Keller, D. Lodge, M. Lewis, and J. Shogren,
editors. Oxford: Oxford University Press.
Jerde, C., and M. Lewis. 2007. Waiting for invasions: A framework for the arrival of
nonindigenous species. American Naturalist 170:1–9.
Jeschke, J., and D. Strayer. 2006. Determinants of vertebrate invasion success in Europe
and North America. Global Change Biology 12:1608–1619.
Jongejans, E., O. Skarpaas, P. Tipping, and K. Shea. 2007. Establishment and spread of
founding populations of an invasive thistle: The role of competition and seed limita-
tion. Biological Invasions 9:317–325.
Klein, G., K. MacIntosh, I. Kaczmarska, and J. Ehrman. 2010. Diatom survivorship in
ballast water during trans-Pacific crossings. Biological Invasions 12:1031–1044.
Kramer, A., and J. Drake. 2010. Experimental demonstration of population extinction
due to a predator-driven Allee effect. Journal of Animal Ecology 79:633–639.
Lavoie, D., L. Smith, and G. Ruiz. 1999. The potential for intracoastal transfer of non-
indigenous species in the ballast water of ships. Estuarine Coastal and Shelf Science
48:551–564.
Lawrence, D., and J. Cordell. 2010. Relative contributions of domestic and foreign
sourced ballast water to propagule pressure in Puget Sound, Washington, USA.
Biological Conservation 143:700–709.
Lee II, H., D. Reusser, M. Frazier, and G. Ruiz. 2010. Density Matters: Review of Ap-
proaches to Setting Organism-Based Ballast Water Discharge Standards.
OCR for page 72
Relationship Between Propagule Pressure and Establishment Risk 111
EPA/600/R-10/031. EPA Office of Research and Development, National Health
and Environmental Effects Research Laboratory, Western Ecology Division.
Lee, S., and J. Bruno. 2009. Propagule supply controls grazer community structure and
primary production in a benthic marine ecosystem. Proceedings of the National
Academy of Sciences of the United States 106:7052–7057.
Leung, B., J. Drake, and D. Lodge. 2004. Predicting invasions: propagule pressure and
the gravity of Allee effects. Ecology 85:1651–1660.
Leung, B., and N. Mandrak. 2007. The risk of establishment of aquatic invasive species:
joining invasibility and propagule pressure. Proceedings of the Royal Society of
London Series B 274:2603–2609.
Levine, J. M., and C. M. D’Antonio. 2003. Forecasting biological invasions with in-
creasing international trade. Conservation Biology 17(1):322–326.
Levings, C., J. Cordell, S. Ong, and G. Piercey. 2004. The origin and identity of inver-
tebrate organisms being transported to Canada’s Pacific coast by ballast water. Ca-
nadian Journal of Fisheries and Aquatic Sciences 61:1–11.
Lewis, M., M. Neubert, H. Caswell, J. Clark, and K. Shea. 2005. A guide to calculating
discrete time invasion rates from data. In: Conceptual Ecology and Invasions Biol-
ogy: Reciprocal Approaches to Nature. M. Cadotte, S. McMahon, and T. Fukami,
editors. New York: Springer.
Locke, A., D. Reid, H. Van Leeuwen, W. Sprules, and J. Carlton. 1993. Ballast water
exchange as a means of controlling dispersal of fresh-water organisms by ships.
Canadian Journal of Fisheries and Aquatic Sciences 50:2086–2093.
Lockwood, J., P. Cassey, and T. Blackburn. 2009. The more you introduce the more you
get: the role of colonization pressure and propagule pressure in invasion ecology.
Diversity and Distributions 15:904–910.
Lonsdale, W. 1999. Global patterns of plant invasions and the concept of invasibility.
Ecology 80:1522–1536.
Lubina, J., and S. Levin. 1988. The spread of reinvading species: range expansion in the
California sea otter. American Naturalist 131:526–543.
Lutscher, F., E. McCauley, and M. Lewis. 2007. Spatial patterns and coexistence me-
chanisms in systems with unidirectional flow. Theoretical Population Biology
71:267–277.
Lutscher, F., R. Nisbet, and E. Pachepsky. 2010. Population persistence in the face of
advection. Theoretical Ecology 3:271–284.
MacArthur, R., and E. Wilson. 1967. The Theory of Island Biogeography. Monographs
in Population Biology. Princeton, NJ: Princeton University Press.
Memmott, J., P. G. Craze, H. M. Harman, P. Syrett, and S. V. Fowler. 2005. The effect
of propagule size on the invasion of an alien insect. Journal of Animal Ecology
74:50–62.
Minton, M., E. Verling, A. Miller, and G. Ruiz. 2005. Reducing propagule supply and
coastal invasions via ships: effects of emerging strategies. Frontiers in Ecology and
the Environment 3:304–308.
Morris, W., and D. Doak. 2002. Quantitative Conservation Biology: Theory and Prac-
tice of Population Viability Analysis. Sunderland, MA: Sinauer Associates.
Neubert, M., and I. Parker. 2004. Projecting rates of spread for invasive species. Risk
Analysis 24:817–831.
Okubo, A., P. Maini, M. Williamson, and J. Murray. 1989. On the spatial spread of the
gray squirrel in Britain. Proceedings of the Royal Society of London Series B
238:113–125.
Pachepsky, E., F. Lutscher, R. Nisbet, and M. Lewis. 2005. Persistence, spread and the
OCR for page 72
112 Propagule Pressure and Invasion Risk in Ballast Water
drift paradox. Theoretical Population Biology 67:61–73.
Peterson, A., and D. Vieglais. 2001. Predicting species invasions using ecological niche
modeling: new approaches from bioinformatics attack a pressing problem. BioS-
cience 51:363–371.
Parker, I. 2000. Invasion dynamics of Cytisus scoparius: A matrix model approach.
Ecological Applications 10(3):726–743.
Ponciano, J. M., M. L. Taper, B. Dennis, and S. R. Lele. 2009. Hierarchical models in
ecology: confidence intervals, hypothesis testing, and model selection using data
cloning. Ecology 90(2):356–362.
Reusser, D. 2010. Biogeography of Nonindigenous Species: From Description to Pre-
diction. Corvallis, OR: Oregon State University.
Ricciardi A. 2001. Facilitative interactions among aquatic invaders: is an “invasional
meltdown” occurring in the Great Lakes? Canadian Journal of Fisheries and Aqua-
tic Sciences 58(12):2513–2525.
Ricciardi, A. 2006. Patterns of invasion in the Laurentian Great Lakes in relation to
changes in vector activity. Diversity and Distributions 12:425–433.
Ruiz, G., and J. Carlton. 2003. Invasion vectors: A conceptual framework for manage-
ment. In: Invasive species: Vectors and management strategies. G. Ruiz and J.
Carlton (editors). Washington, DC: Island Press.
Ruiz, G., P. Fofonoff, J. Carlton, M. Wonham, and A. Hines. 2000. Invasion of coastal
marine communities in North America: Apparent patterns, processes, and biases.
Annual Reviews in Ecology and Systematics 31:481–531.
Shea, K., and H. Possingham. 2000. Optimal release strategies for biological control
agents: An application of stochastic dynamic programming to population manage-
ment. Journal of Applied Ecology 37(1):77–86.
Shigesada, N., and K. Kawasaki. 1997. Biological Invasions: Theory and Practice. Ox-
ford: Oxford University Press.
Shurin, J. 2000. Dispersal limitation, invasion resistance, and the structure of pond
zooplankton communities. Ecology 81:3074–3066.
Shurin, J., J. Havel, M. Leibold, and B. Pinel-Alloul. 2000. Local and regional zoop-
lankton species richness: a scale-independent test for saturation. Ecology 81:3062–
3073.
Skellam, J. 1951. Random dispersal in theoretical populations. Biometrika 38:196–218.
Smith, L., M. Wonham, L. McCann, G. Ruiz, A. Hines, and J. Carlton. 1999. Invasion
pressure to a ballast-flooded estuary and an assessment of inoculant survival. Bio-
logical Invasions 1:67–87.
Solow, A., and C. Costello. 2004. Estimating the rate of species introductions from the
discovery record. Ecology 85:1822–1825.
Suarez, A., D. Holway, and P. Ward. 2005. The role of opportunity in the unintentional
introduction of nonnative ants. Proceedings of the National Academy of Sciences of
the United States 102:17032–17035.
Taylor, B., and R. Irwin. 2004. Linking economic activities to the distribution of exotic
plants. Proceedings of the National Academy of Sciences of the United States
101:17725–17730.
Tilman, D. 1997. Community invasibility, recruitment limitation, and grassland biodi-
versity. Ecology 78:81–92.
USCG. 2008. United States Coast Guard Ballast Water Discharge Standards Draft Pro-
grammatic Environmental Impact Statement.
OCR for page 72
Relationship Between Propagule Pressure and Establishment Risk 113
Verling, E., G. M. Ruiz, L. D. Smith, et al. 2005. Supply-side invasion ecology: charac-
terizing propagule pressure in coastal ecosystems. Proceedings of the Royal Society
B-Biological Sciences 272(1569):1249–1256.
Wonham, M., and J. Carlton. 2005. Trends in marine biological invasions at local and
regional scales: the Northeast Pacific Ocean as a model system. Biological Inva-
sions 7:369–392.
Wonham, M., and M. Lewis. 2009. Modeling marine invasions: current and future ap-
proaches. Pp. 71–105 In: Marine Bioinvasions: Ecology, Conservation, and Man-
agement Perspectives. G. Rilov and J. Crooks, editors. Berlin: Springer-Verlag.
Wonham, M., and E. Pachepsky. 2005. A null model of temporal trends in biological
invasion records. Ecology Letters 9:663–672.
Wonham, M. J., W. C. Walton, G. M. Ruiz, A. M. Frese, and B. S. Galil. 2001. Going to
the source: role of the invasion pathway in determining potential invaders. Marine
Ecology-Progress Series 215:1–12.
Wood, D., J. Koprowski, and P. Lurz. 2007. Tree squirrel introduction: A theoretical
approach with population viability analysis. Journal of Mammalogy 88(5):1271–
1279.