Northeast Fisheries Science Center Reference Document 04-02
SalmonPVA:
A Population Viability Analysis Model
for Atlantic Salmon in the Maine Distinct Population Segment
by Christopher M. Legault
National Marine Fisheries Serv., Woods Hole Lab., 166 Water St., Woods
Hole, MA 02543
Print
publication date January 2004;
web version posted January 16, 2004
Citation: Legault, C.M. 2004. Salmon PVA: a population viability analysis model for Atlantic salmon in the Maine
Distinct Population Segment. U.S. Dep. Commer. Northeast Fish. Sci. Cent. Ref. Doc. 04-02; 88 p.
Download complete PDF/print version
EXECUTIVE SUMMARY
A population viability analysis (PVA) model has been developed for
Atlantic salmon (Salmo salar) in Maine. This model incorporates
uncertainty in juvenile and adult survival rates, direct and indirect
linkages among populations in different rivers, and accounts for a
number of sources of potential human removals or stocking in a flexible,
modular Fortran program named SalmonPVA. The structure of the model
is based on a state-space approach with a detailed life history cycle.
Multiple cohorts in multiple rivers progress through their life history
based on stage-specific survival rates and fecundity, with limits imposed
by riverine habitat capacity. The model projects the populations forward
in time (usually 100 years) numerous times, with stochastic variables
selected based on a Monte Carlo approach to calculate the probability
of extinction. This model is being developed with input from scientists
and policy makers from NOAA Fisheries, US Fish and Wildlife Service,
Atlantic Salmon Commission, and the University of Maine. Results from
this model will form the basis for informing the selection of delisting
criteria in the Recovery Plan for the Maine Distinct Population Segment,
which was listed as endangered in 1999.
The SalmonPVA model is structured to represent Atlantic salmon life history
characteristics in the US. For example, most fish spend three years in
the river and two years at sea before returning to the river to spawn.
However, the model allows for the possibility to return from sea after
one or three years, and for some spawners to survive (termed kelts) and
spawn again two years later. Inputs to the model allow for a wide range
of simulations incorporating various combinations of mechanisms and parameter
variability. The number of rivers is a dynamic variable limited only by
the capacity of the computer running the program. Linkages among rivers
are determined as input and allow for various straying hypotheses as well
as for linkages among juvenile survival rates due to year effects. Habitat
capacity limits are input and combined with the approach used for fecundity
cause a Beverton and Holt type spawner-recruitment relationship. This underestimates
the probability of extinction when populations are large relative to a
Ricker type spawner-recruitment relationship. However, the salmon populations
in Maine are currently so low that this concern is not of immediate consequence.
A number of sources of human-induced removals from the populations are
allowed (but not required) by the model including interception fisheries
at sea, in river fishing, and broodstock removals of either returning adults
or parr. Stocking of any life stage during any year of the simulation is
also allowed within the model. Stocked fish are followed in a separate
matrix in the program from the natural fish to allow for different survival
rates or removals. The offspring from the hatchery fish are added to the
natural population so that hatchery populations disapper if stocking is
discontinued. The model allows direct examination of specific simulations
as well as aggregate summary results from all of the simulations. The probability
of extinction is the most important output, but trends in adult returns
can also be examined. Population trend information is useful in combination
with risk of extinction because a five percent chance of extinction in
one hundred years has different implications if the overall trend for the
population is increasing or decreasing over the projected time series.
1. Introduction
In November 2000, the Gulf of Maine Distinct Population
Segment (DPS) of Atlantic salmon was listed as endangered under the Endangered
Species Act (ESA) by the Fish and Wildlife Service and the National Marine
Fisheries Service (Fed. Regist. 65(223):69459-69483). Eight coastal rivers
in Maine contained extant populations of Atlantic salmon within the DPS
covered by the ESA listing (Figure 1). This listing
brought additional protections for the species as well as requirements
for creating objective and measurable delisting criteria. One approach
to developing delisting criteria is by using a population viability analysis
(PVA), a stochastic modeling technique for predicting changes in population
abundance given uncertain biological parameters (Beissinger 2002). This
approach was selected for the Atlantic salmon in the Gulf of Maine DPS
to project future abundance and to provide information for establishing
recovery criteria. PVA can also be used to help direct research by quantifying
the relative improvement in population viability due to different management
decisions or the reduction of uncertainty in key biological parameters.
This document describes the PVA developed for the Atlantic salmon DPS in
terms of the modeling decisions as well as the available input data.
Decisions made in constructing the PVA and in evaluating which data to
use were made by a Working Group comprised of experts of salmon biology.
This group consisted of representatives from the National Marine Fisheries
Service (NMFS), US Fish and Wildlife Service (FWS), and Maine Atlantic
Salmon Commission (ASC) (see Appendix for participant
list). The group met as needed to discuss model development and input data
selection and reached agreement through consensus. This activity benefitted
greatly from the many years of experience of the Working Group members
and the collegial atmosphere at these meetings.
2. Theory
Population viability analysis is a technique to estimate the probability
of a stock attaining given sizes, usually zero or very low, sometime in
the future (Gilpin and Soulé, 1986). A wide range of modeling approaches
are used in PVA, from simple extrapolation of current trends to complex
individual based models (Beissinger and McCullough, 2002). Software to
conduct PVAs is widely available, e.g. RAMAS and Vortex, but models built
specifically for a given species have also been utilized (Nickelson and
Lawson, 1998; Reed et al., 2002). Whatever approach is taken, the purpose
is the same, to predict the probability of the population persisting into
the future. Because predictions are difficult, PVA results should not be
taken as absolute truth. Rather, the forecasts should be used to explore
potential consequences of management actions in the light of an uncertain
future and variation in input assumption and data.
A life history modeling approach was selected for the Atlantic salmon
Gulf of Maine DPS due to the large amount of data available for the species.
This approach has the benefit of higher biological realism but requires
many more input parameters and distributions relative to simpler PVA models.
Complex features of Atlantic salmon biology (such as anadromy, precocious
parr, kelting, and hatchery supplementation) are captured in the model,
but at the cost of having to make many assumptions (such as how survival
in riverine life stages is related among the DPS rivers). These features
and assumptions are described below. Verification of the model and input
assumptions was conducted using historical data for the Narraguagus river.
The status of the DPS was projected under the most likely input assumptions,
and sensitivity analyses were conducted for a number of input parameters.
Finally, recovery criteria were explored with the model.
3. Program SalmonPVA
SalmonPVA, written in Fortran90 using IMSL numerical
routines wherever possible, uses a modular approach to track the fate of
populations from multiple rivers separately through time. The core structure
is based on Atlantic salmon life history characteristics in Maine, but
is flexible in terms of number of rivers and years as well as survival
rates, stocking, fishing rates, straying, and other components. The program
consists of input and output sections with two main loops: a simulation
loop that builds the resulting distributions of extinction probability,
number of spawners, habitat limitation, etc., and a year loop that projects
the multiple populations into the future (Figure 2).
SalmonPVA is a state-space model in which the number
of fish alive in a cohort is simulated based on the number in the previous
life history stage and time step and the rates of removals between steps
(Figure 3). The structure of the model is fixed,
but vital rates are input as distributions and Monte Carlo sampling is
utilized to simulate many realizations of the populations forward in time.
Summary statistics from these realizations are used to generate the probability
of events such as extinction in a river or group of rivers, or the distribution
of response variables such as the replacement rate. Populations are considered
unique to rivers, with fish that stray from one river to another taking
on all the vital rates of the fish found in the latter river.
Life history stages in this model were set such that
year time steps could advance cohorts through most stages (Table
1). Eggs are laid in October of year 1 and become fry 8 months later
in June of year 2. Two months later, the exception to the year advancement
per life stage, the fish become parr0+. In August of the following year
the fish advance to parr1+ and then become smolts in May of year 4. Adults
return to the rivers in June of years 5, 6, and 7 as 1 sea winter, 2 sea
winter and 3 sea winter adults, respectively. Eggs are laid in the same
year as the adults return to the rivers, completing the life cycle. Atlantic
salmon are iteroparous, surviving adults become kelts which move to the
sea and return to spawn two years later. Only integer values of fish and
eggs are allowed in the simulation.
For bookkeeping purposes, SalmonPVA uses six separate
matrices (arranged in two sets of three matrices) to store the number of
fish in each life stage for each river and year. The two sets are for the
natural and hatchery populations while the three matrices are for the different
life history stages. In the model, natural fish are hatched from eggs in
the rivers and spend their entire life either in the river or at sea. In
contrast, the hatchery population consists of fish stocked at some life
stage. Note that the progeny of hatchery fish are considered natural fish
in this model, and the same survival rates are applied to both the natural
and hatchery populations. The juvenile matrices contain the egg through
smolt life stages. The adult matrices keep track of the number of fish
of each sea winter age group (1 through 3) in four separate locations or
states within a given year: at sea, return to river prior to straying,
return to river after straying, and spawners. The kelt matrices keep track
of which fish have just become kelts and which are returning to spawn.
The kelt matrices are the only ones which keep track of the sex of the
fish. The progression of fish through life stages, along with linkages
among rivers, are described below and shown schematically in Figure
4.
3.1 Juvenile Survival
Juvenile survival is modeled for each of the four life cycle transitions:
egg to fry, fry to parr0+, parr0+ to parr1+, and parr1+ to smolt. For each
transition, a minimum and maximum main effect is input along with a noise
level. The minimum and maximum values for the main effect define a uniform
distribution for survival at that stage. The noise level is assumed to
be a uniform distribution centered on zero with width two times the noise
level (i.e., from negative NL to positive NL, where NL is the input noise
level). To create a survival rate for a specific juvenile transition, river,
and year, a random value is selected from both the main effect and noise
level distributions and these two values added together.
The number of fish that pass from one stage to the next is randomly selected
from a binomial distribution, with initial population abundance as the
number of events and the survival rate as the probability of success. Thus,
the number of fish at the next life stage in the next time step is bounded
by zero and the initial population abundance, with the expected value equal
to the product of the initial population abundance and the survival rate.
This assumes that each fish has an equal and independent probability of
surviving through the time step. When the number of fish in the starting
life stage is large (thousands), then the number of fish in the next life
stage at the next time step approximately follows a normal distribution.
However, when the population abundance in the starting population is low
(tens), then the variability due to the binomial process critically affects
the persistence of the stocks.
The random survival values for each life stage within a realization are
generated prior to the year loop from one call to the IMSL routine DRNUN
which returns a vector of uniform random numbers of length equal to the
number of years being simulated. These random numbers are then scaled by
the range of survival for that life stage. This is just a time saving feature
relative to calling the random number generator within the year loop. The
survival rate is independent of population abundance for all juvenile life
stages, and so can be generated at the start of the simulation. However,
the application of the survival rate through the binomial distribution
must occur within the year loop as it depends upon the population abundance
at that life stage. The application of the binomial distribution both here
and elsewhere is accomplished using the IMSL routine RNBIN. For example,
in Figure 3 the uniform distribution for survival of eggs to fry is (0,25,
0.35). A value for the year effect would be chosen randomly from this distribution,
say 0.31, and then the random noise for the river effect, say -0.01, would
produce the value of 0.30 shown. Thus, the survival rate for the particular
year and river is formed from the sum of the year effect and the river
effect, both of which are uniformly distributed random values. Application
of this 30% survival rate to the 15,000 eggs would be expected to produce
4,500 fish (15,000 * 0.30 = 4,500), but the use of the binomial distribution
instead produced 4,417 fry as survivors in year 2 in this example. For
large sample sizes, the binomial distribution is similar to the normal
distribution, but for small sample sizes skewed distributions can be generated
because only integer values are possible. The added noise at the river
effect level can cause the selected survival value to fall outside the
original range, as in this example for year 3 egg to fry survival. The
program checks to ensure that the survival rate is bounded by zero and
one and exits with an error message during input if a survival rate outside
of zero and one is possible.
Survival rate linkages are allowed both among juvenile stages and among
rivers, but these linkages are not required. The linkages are defined by
the main effect. In the case of the juvenile stage linkages, the same proportional
value from within the two main effect distributions is used. For example,
if the egg survival rate main effect is the interval (0.15,0.35) and the
fry survival rate main effect is the interval (0.43,0.60) then linking
these stages could produce an egg survival rate of 0.20 and a fry survival
rate of 0.4725, as both are 25% from the start of the interval. The program
allows none, all, or some of the juvenile stages to be linked. The none
and all options require only a single input flag, while the some option
requires definition of the linkages. This is accomplished by entering the
number of linkages followed by the pairs of linked life stages. For example,
if only the two parr stage survivals were desired to be linked, then the
user enters the flag for some linkages followed by the number 1 for one
linkage followed by the set of numbers 3, 4 for the third and fourth juvenile
survivals. Any linkages among life stages are applied to all rivers. In
a similar manner, linkages are allowed among the rivers; none, all, or
some. The survival rate main effects for the linked rivers are all the
same. If river 1 and river 2 are linked and river 2 and river 3 are linked,
then rivers 1 and 3 will be linked whether or not this linkage is explicitly
input.
The relative level of noise to main effect can overcome
any linkages among rivers. The linkages allow "good" and "bad" years to
occur for linked life stages and rivers, but the linkages must be combined
with the noise level to determine the specific survival rate for a stage,
river, and year. By varying the relative levels of main effects and noise,
various life stages and rivers can vary between completely linked to completely
independent (Figure 5).
3.2 Marine Survival
Marine survival was modeled assuming a cyclic pattern (sine wave) that
varies in mean, period, amplitude, and phase over time. The first year
at sea is assumed to have a lower survival rate than subsequent years at
sea, with the difference between the two annual rates established through
an input variable. The lower marine survival during the first year at sea
is thought to be due to additional mortality caused during the transition
from freshwater to the marine environment, but the actual mechanism is
unknown. A component of random noise is added to the expected survival
rate from the sine waves to more closely match the non-smooth variability
observed in marine survival rates throughout the world. These cycles may
be related to environmental cycles, such as the North Atlantic Oscillation,
but no causal mechanism has been demonstrated for this relationship and
there are still not enough data to adequately judge long term correlations.
Cycles of marine survival are modeled in SalmonPVA by assuming an empirical
distribution for the cycle period, and uniform distributions for the cycle
mean, amplitude, and phase. The empirical distribution for the cycle period
is used to generate a generic cumulative distribution function from which
values are sampled. A sufficient number of periods are collected during
each realization to cover the range of years simulated. For each cycle,
random realizations from the uniform distributions of the cycle mean, amplitude
and phase are selected and applied to that cycle for the appropriate period
of years. In this way, marine survival cycles between high and low values
but for random durations, at random levels, with random ranges, and in
random sequences that match the variability observed in the environment.
By varying both the relative ranges assumed for the cycles and the additional
noise, marine survival can be modeled as noise about a fixed mean, a constant
sine wave, a near chaotic series, or any intermediate pattern.
As in the juvenile population, the marine survival rate in any year and
sea winter age combination is applied using a binomial distribution, with
the probability of success determined by the survival rate. The survival
rates for all years and for both the first and later sea winters are computed
for each realization prior to the year loop, but the application of the
survival rates through the binomial distribution occurs for each year.
Note that the initial number of 2 and 3 sea winter fish at sea is the number
of these cohorts in the sea during the previous year minus the number of
these fish that returned to rivers the previous year. The marine survival
rate is the same for all stocks as they all inhabit the same environment.
However, the use of the binomial distribution to determine the actual number
of survivors means that two rivers with the same number of emigrating smolts
in a given year can produce different numbers of one sea winter adults
at sea.
3.3 Habitat Limitations
The capacity of rivers to support Atlantic salmon is limited by the amount
of appropriate habitat. The abundance of all juvenile life stages, except
smolts, may be limited by habitat. Smolts are not limited because they
leave the rivers. Within the simulation, the habitat limitation occurs
after natural mortality and human induced removals have occurred. All four
juvenile stages can have a habitat cap, but a cap is not required for each
stage. Data availability and mechanistic hypotheses should determine which
stage(s) are limited. Elliot (2001) presented evidence of density-dependence
in Atlantic salmon, but the stage at which the limitation occurred was
not defined. Elson (1975) found a linear relationship between egg density
and fry but found a density-dependent relationship between egg density
and smolts. The habitat capacity approach used in SalmonPVA is equivalent
to a Beverton and Holt type stock-recruitment relationship. The use of
a Ricker type stock-recruitment relationship would be more pessimistic
in terms of recovery potential, but there are no observations for recruitment
levels at high stock sizes for Atlantic salmon in Maine upon which to base
a Ricker type relationship.
The habitat capacity for a given juvenile stage in a particular river
is determined from a mean level for each river, a random year effect applied
to all rivers, and a random river specific effect. Both the year and river
effects are multiplicative and defined by uniform distributions such that
if set to zero cause no change from the mean. For example, if the year
effect is desired to change the habitat capacity by plus or minus 10%,
the entered range for the uniform distribution would be -0.1, 0.1. In this
way, all rivers can experience "good" and "bad" years or fluctuate independently,
depending upon the relative variability in the year and river effects.
For each realization, the habitat limitations for each year and river are
derived prior to the year loop as they are independent of the salmon population
dynamics. These "good" and "bad" years of habitat limitation are independent
of the "good" and "bad" years of juvenile survival. The habitat caps for
each stage are chosen independently such that "good" and "bad" years can
occur simultaneously for different stages within a river.
Habitat limitation for any juvenile stage can occur in one of three ways:
hatchery fish are removed first, natural fish are removed first, or the
two sets of fish are removed in proportion to their relative abundance.
As such, hypotheses related to ecological effects can be examined. For
example, it has been hypothesized that hatchery fish are more susceptible
to predation than are naturally spawned fish. In all cases, the total number
of fish that will become smolts in a given year and river are summed from
the hatchery and natural populations. If this sum is less that the corresponding
habitat limitation, then nothing happens. However, if this sum is greater
than the habitat limitation then the number of fish that are removed is
the difference between the sum and the cap. If one set of fish, hatchery
or natural, is selectively removed then the number of fish in this set
is compared to the number to be removed. If the number to be removed exceeds
the number of fish in this set, then the rest of the removals are taken
from the other set. The proportional option reduces each set by the ratio
of the number to be removed and the number present. For example, if 60
hatchery fish and 500 natural fish are present in a river that has a habitat
limitation of 400 fish that year then under the three options there could
be 0 hatchery and 400 natural fish (hatchery first), 60 hatchery and 340
natural fish (natural first), or 42 hatchery and 358 natural fish (proportional)
that remain after application of the habitat limitation.
3.4 Population Initializations
Two options exist for initializing the populations in each river: entering
a number of years of data for a specific life stage, or entering all life
stages in a single year. Both options utilize uniform distributions to
reflect the uncertainty in the abundance estimates of most life stages.
Different options can be used for different rivers. All life stages are
projected forward to the last year with initialization data using all the
processes described in this section, with the exception of straying which
is not used during the initialization. Input data should be sufficient
to ensure fish are present at all life stages in the first true year for
projections. However, if there are uninitialized life stages, these values
will be set to zero by the program. If enough data are entered such that
a cohort is projected into an initialized life stage, then the initialized
value will be used for the natural population. In contrast, previously
stocked life stages are additive when projected forward into an initialized
life stage because hatchery input is the number stocked not a population
estimate. Note that since the fry to parr0+ transition occurs within a
year time step, and this time step is done during the projection part of
the time loop, if annual fry values are entered for initialization, then
one more year than expected will be required to fill the cohorts. If adult
stages are input, then ranges for each sea winter age should be input.
Also note that when the second option for initialization is chosen, all
life stages in a single year, the program will not allow the number of
returns to be greater than the number of fish at sea for any of the sea
winter ages. This is to ensure that the calculation of fish at sea in the
next year does not start at a negative value.
3.5 Human Induced Removals
3.5.1 Interception Fisheries
Fish at sea can be caught in directed or non-directed fisheries. The interception
fisheries remove fish from the river-specific populations randomly in proportion
to their relative abundance. The input for interception fishery removals
is a uniform distribution for each year, which reflects the uncertainty
in estimating the number of Maine salmon caught in the mixed fisheries
of Greenland or St. Pierre et Miquelon. The number of fish caught in a
year is selected randomly from this range and if it exceeds the total number
of fish at sea from all rivers, then all these fish are removed. If the
total removal is less than the sum of the population currently existing
in the run, then random fish are selected from the populations of natural
and hatchery fish. The removal occurs temporally prior to the return of
fish to the rivers in a given year.
3.5.2 Broodstock
Some adults that return to rivers to spawn are collected as broodstock
for restoration efforts. These collections are based on two rules applied
in the model: the number desired and the maximum percent of the run that
can be removed as broodstock. The number of returning fish to a river is
first compared to the desired total collection. The smaller of these two
values becomes the first check. The second check is calculated as the maximum
percentage of the returning fish allowed for broodstock collection. The
lower value of these two checks is taken as the broodstock for that river,
year and sea winter age. Note that since each sea winter age is treated
separately, the total broodstock collection may be less than if all sea
winter fish were aggregated. However, broodstock collection focuses on
certain sea winter ages normally, and so the model treats them separately.
Only natural fish are collected for broodstock.
In a similar manner, natural parr1+ are collected from some rivers through
electrofishing for raising in hatcheries to serve as broodstock for future
restoration efforts. In this case, only a maximum number of fish are entered
and it is assumed that electrofishing operations occur only in stream reaches
where natural fish are expected to occur. If the value entered is larger
than the natural population size in a particular river for a given year,
then all the natural fish are removed. Note that for projections, this
option should only be used when a restoration stocking program is employed.
3.5.3 River Adults
When salmon return to a river to spawn, there is the potential for fishery
removals to occur, through a directed fishery, as bycatch in a non-directed
fishery, or through poaching. In the model, these in river removals are
assumed to occur before broodstock collection. The river fishery removals
are entered for each sea winter age separately to allow for directed fishing
on different age classes or poaching of different age classes at different
rates. The input for river adult removals is from a uniform distribution
of the probability of a fish being caught. This probability is used with
the binomial distribution to randomly remove a fraction of fish due to
fishing. This approach mimics the expectation that larger returns will
have larger removals, but at a given rate determined by management (for
a directed fishery) or bycatch/poaching (for non-directed or illegal fisheries).
The catch of natural and hatchery fish are treated separately to allow
for a directed fishery on stocked fish.
3.6 Straying
Although Atlantic salmon have a remarkable ability to return to their
natal river, there are occasions when fish return to a different river,
termed straying. The model uses two input values to determine the proportion
of natural and hatchery fish that stray (these can be set to different
values). Two separate input values determine the fraction of natural and
hatchery strays that will die, respectively. Finally, a single input matrix
determines to which other rivers the fish will stray and in what expected
proportions. In the model, the stray rate is defined as the proportion
of fish from the total that is expected to return to a specific river that
instead returns to a different river. For example, if 10 of 100 Narraguagus
fish return to the Pleasant river, the stray rate for the Narraguagus river
is 10% independent of the number of fish returning to the Pleasant river.
The stray rate is applied to all fish within the population, natural or
hatchery, but applied to each sea winter age separately using a binomial
distribution to determine the number of strays from each river. This initial
number of strays is then decreased according to the stray mortality fraction
through the use of the binomial distribution. Once the number of strays
that will survive is determined for each river, the river of return for
these fish is determined by the straying matrix. The straying matrix has
river of origin as the rows and river of return as the columns. The diagonal
of the matrix should always be zero because the stray rate assumes that
an actual stray occurs and the diagonal would put the fish back into the
river of origin. However, the program will allow positive values along
the diagonal such that stray rates can be made different among the rivers.
The program checks that the sum of each row in the stray matrix is equal
to one, and rescales to sum to one if not. In this way, physical distance
or some other measure of straying can be used during input, and the program
will calculate the straying matrix. The allocation of strays from a given
river to all other rivers is done randomly but is based on the probability
associated with the transfer in the straying matrix. Straying occurs after
fish return to the river, but before any broodstock collection or fishing
mortality occur, and before the counting of returns to the river.
3.7 Maturity and Fecundity
Atlantic salmon first return from the sea to spawn after spending from
one to three years at sea. The maturity rate (proportion of fish at returning
at each sea winter age) is stock dependent, with DPS fish predominantly
returning as 2SW fish. A maturity vector is input that is held constant
throughout the simulation, but applied through a binomial distribution
to incorporate variability in the maturity rate.
Spawning fish are determined to be either male or female based on an input
probability by sea winter age. For spawning to occur, both males and females
must be present in the river. Males can be either adults or precocious
parr (parr 1+) but females must be adults. The females spawn a quantity
of eggs that is chosen randomly from an input normal distribution, with
mean and variance set by sea winter age. Each female is treated independently,
such that the resulting average number of eggs per female will be close
to the mean when the number of females is large, due to the central value
theorem, but will be more variable when the number of females is low. Hatchery
fecundity for each female spawner is multiplied by an input constant that
allows for either decreased or increased fecundity relative to natural
spawners. The eggs produced by both natural and hatchery adults are placed
into the natural population. Thus, after one life cycle of the ending of
stocking, the hatchery populations will consist only of surviving kelts.
3.8 Stocking
Restoration efforts for Atlantic salmon have used stocking of hatchery
reared fish to supplement natural populations. Different life stages of
fish have been stocked in different rivers and at different times. The
model allows user input of fish into the hatchery population at any stage
and in any year. The numbers of fish stocked can also be varied over time,
although the input format is designed for blocks of years with the same
number of fish stocked at a given life stage. As noted above in the fecundity
section, the hatchery population is not self-sustaining, fish only remain
in this population until they spawn at which point their progeny are considered
natural. There is not a random component to the stocking. This is because
the levels of stocking are considered to be so high and discrete that any
randomness would be artificial. Basically the questions are whether or
not to stock, or which life stage to stock, but not variation in how many
fish are stocked.
3.9 Counters
Given the large number of inputs, with imprecise estimation, to the Atlantic
salmon population viability analysis, the model has to be run many times
to produce distributions of response variables under different parameter
combinations. The variables of interest collected from each realization
are primarily related to population size and extinction, although habitat
limitation effects and replacement rates are also evaluated. The number
of realizations in which a given river went extinct or was limited by the
habitat capacity is counted and reported in the program output. The average
number of spawners over an input defined number of years is also computed
and reported for each realization. This averaging allows the effect of
cyclic marine survival to be mitigated in terms of expected numbers of
spawners in the future, as using just a single year would induce more variability
because the output would then depend upon where in the marine cycle the
observation of survival was taken. Similarly, the replacement rate for
a realization is calculated as the median of the last twenty years in the
simulation. The replacement rate for year t is defined as the sum of adult
offspring for that cohort in years t+4 to t+7 divided by year t adult returns
(Rago, 2001). The rivers can be grouped into any number of geographical
or size combinations to examine patterns in extinction, habitat capacity,
and average number of spawners. A river can also be included in multiple
groups.
The complete output from individual realizations can also be reported
by the model including: natural and hatchery population abundance for all
rivers, years, and life stages; removals due to intercept fishing, river
fishing, and broodstock; and habitat limitations for each river and year
and whether or not the cap was utilized for each river and year. Since
a given random number seed will always produce the same output, multiple
realizations from a single scenario can be collected by using the same
random number seed and changing which realization is output.
3.10 Components Currently Not Included in Model
Catastrophic mortalities are currently not included in the model. This
phenomenon could be added, but would require many more model realizations
due to being, by definition, a low probability event. It is easier to consider
the impacts of catastrophic events a posteriori, by taking the
results and increasing the probability of extinction by some small amount
corresponding to the hypothesized probability of catastrophic events.
Genetics are not considered at all in the model and would be difficult
to incorporate unless a qualitative aspect is used. For example, decreasing
fecundity if the population size goes below a certain limit.
Mortality associated with dam passage is not included directly in the
model but can be aliased using the river fishing removals. The user would
have to ensure that these removals were attributed to dams and not fishing
when describing results.
Disease is not included directly in the model. It would be difficult to
incorporate this as a transmissible vector, but disease might be considered
as more of a catastrophic effect.
The impact of aquaculture escapees is also not included directly in the
model but could be approximated using river fishing removals. Again, when
describing results the user would have to ensure that these removals were
attributed to aquaculture interactions and not fishing.
Predation is not modeled directly, but rather as a consequence of the
survival rates. Although many predators of Atlantic salmon have been identified,
the ability to predict future predator abundance, or that of alternative
prey, makes direct incorporation into the model difficult.
4. Input for Base Case
4.1 Juvenile Survival
Estimates of survival rates for the juvenile life stages of Atlantic salmon
were obtained from the literature and combined using an objective process
that accounted for the uncertainty in each study. Some subjectiveness entered
the process through decisions to include or exclude specific results in
the process. The combination process first standardized survival values
for all studies for a particular life stage to the same time interval assuming
that the reported survival rate would be constant in the new time period.
For example, if the reported survival was 20% over 10 months, then a standardized
survival for 12 months would be 14.5%, calculated as 0.2(12/10).
This conversion was done to the mean as well as the minimum and maximum
survival rates from each study. The standardized triplets of minimum, mean,
and maximum survival for each study were then combined assuming triangular
distributions with zero probability at the minimum and maximum values and
the probability associated with the mean fixed to give an unit area under
the curve, calculated as 2/(maximum-minimum). The triangles from studies
selected for inclusion in the final estimate were then simply added together
and rescaled to form a new distribution of the probability of each survival
value for that life stage. This new probability distribution function was
then converted to a cumulative distribution function and the 10th and
90th percentiles selected as the limits of a uniform distribution
to describe the uncertainty associated with survival of that life stage.
The studies included in the generation of the uniform distribution for
survival were based on group discussions involving representatives from
NMFS, FWS, and ASC. The choices made for each life stage are detailed below.
The survival rates for egg to fry life stages came directly
from a study conducted on Maine DPS Atlantic salmon (Jordan and Beland,
1981). The resulting uniform distribution for survival is 15-35% which
was taken directly from the study instead of applying the objective process
described above. This range in survival covers most other estimates available
in the literature (Table 2) and is thought to
best represent survival of Atlantic salmon in Maine.
The survival rates for fry to parr0+
life stages were derived using the objective process described above, with
the standard time period of two months. There were 13 studies found in
the literature for this life stage transition which were used to generate
a total of 24 possible survival triangles (Table 3).
Of these 24 possibilities, five were selected for use in the objective
procedure to create the uniform distribution of 43-60% survival. The 19
other possibilities were not considered representative of survival of Atlantic
salmon in Maine for a variety of reasons. The duration of a number of studies
could not be determined, which was deemed too important to overlook for
this life stage (Stewart, 1963; Mills, 1969; Greenwood, 1981; Knight et
al., 1982; Kennedy, 1984). One study was conducted specifically as part
of a low competition experiment (Heggenes and Borgstroem, 1991). Some studies
were conducted in low productivity streams in Vermont (Whalen and LaBar,
1994; Whalen and LaBar, 1998) that are not considered representative of
Maine DPS rivers. The two parts of the McMenemy (1995) study were averaged
after adjusting for the different time periods to prevent this one study
from having too much influence on the overall calculation of survival for
this stage. The five selected studies had mean standardized survival rates
ranging from 48.6-59.2% (Egglishaw and Shackley, 1973; Egglishaw and Shackley,
1980; Gardiner and Shackley, 1991; Orciari et al., 1994; McMenemy, 1995).
The range of survival rates in these studies were quite similar and resulted
in a relatively narrow distribution of survival rates for this life stage
(Figure 6).
The survival rates for parr0+ to parr1+
life stages were derived using the objective process described above with
a standard time period of twelve months. Ten studies were found in the
literature for this life stage transition, and these were used to generate
12 possible survival triangles (Table 4). Of these
12 possibilities, seven were selected for use in the objective procedure
to create a uniform distribution of 12-58% survival. One study was conducted
specifically as part of a low competition experiment (Heggenes and Borgstroem,
1991). The unknown duration of two studies precluded their selection (Symons,
1979; Knight et al., 1982). The survival rates from the two time periods
in the Orciari et al. (1994) study were averaged after standardization.
The seven selected studies had mean standardized survival rates ranging
from 11.3-50.2% (Meister, 1962; Egglishaw and Shackley, 1980; Kennedy and
Strange, 1980; Kennedy and Strange, 1986; Gardiner and Shackley, 1991;
Orciari et al., 1994; Cunjak et al., 1998). The range of survival rates
in these studies was quite large, resulting in a wide distribution of survival
for this life stage (Figure 7).
The survival rates for parr1+ to smolt
life stages were derived using the objective process described above with
a standard time period of nine months. Eight studies were found in the
literature for this life stage transition and together with one set of
data from the Narraguagus river (J. Kocik, NOAA Fisheries, pers. comm.)
were used to generate 16 possible survival triangles (Table
5). Of these 16 possibilities, five were used to create a uniform distribution
of 17-50% survival. The unknown duration of two studies obviated their
selection (Elson, 1957; Symons, 1979). The both used Average standardized
survival rates were calculated from the different sets of fish in the Myers
(1984) and Orciari et al. (1994) studies to prevent one study from having
too much affect in the overall calculations. Two studies were conducted
in Vermont rivers, not considered to be representative of Maine DPS rivers
(Whalen, 1998; Whalen et al., 2000). The five selected studies had mean
standardized survival rates ranging from 16.8-45.8% (Meister, 1962; Myers,
1984; Orciari et al., 1994; Cunjak et al., 1998; Kocik pers. comm.). The
range of survival rates in these studies was quite large and resulted in
a relatively wide distribution of survival for this life stage (Figure
8).
Combining these survival rates produced
a possible range of egg to smolt survival of 0.13-6.09% (Table
6). However, the distribution of egg to smolt survival is not uniform,
but approximates a lognormal distribution (Figure
9). This is because the sum of random values from any distribution
is approximately normal for large sample sizes and egg to smolt survival
can be expressed as the sum of the natural logs of each survival rate.
Survival values between 0.5 and 3.5% occur within the 90% confidence interval
of this distribution, which corresponds to the general impression expressed
by the working group and found in the literature (Bley and Moring 1988)
that egg to smolt survival should be around 1-2%.
Linkages in survival rates were assumed to occur among all rivers and
between the two parr life stages. The noise level for each life stage was
set at 0.05 to generate relatively tight linkages, relative to the variability
in the main effects. Relatively tight linkages among rivers are justified
based on hydrological studies where, for example, data for the Machias
River are used to supplement and extend the record for the Dennys River.
The relatively tight linkage between the two parr survival rates is based
on the similarity in habitat requirements of these two stages.
4.2 Marine Survival
Correlations exist between sea surface temperature and marine survival
of Atlantic salmon (Scarnecchia 1984; Martin and Mitchell 1985; Scarnecchia
et al. 1989; Friedland et al. 1993; Friedland et al. 1996; Friedland 1998;
Friedland et al. 2003). One possible causal mechanism for these relationships
for Maine DPS Atlantic salmon is the North Atlantic Oscillation (NAO).
Although no direct link has yet been established between the NAO and Atlantic
salmon marine survival, a cyclic pattern in marine survival is common to
many species of salmon throughout the world. For example, marine survival
of various species of Pacific salmon has been shown to vary with El Niño
Southern Oscillation (ENSO) events (Johnson 1988; Beamish and Bouillon
1993; Francis and Hare 1994).
The NAO is defined as the normalized pressure difference
between the Azores and Iceland. Monthly averages are computed and the deviations
from December to March are averaged to create an annual winter value. NAO
values from 1824 to 2001 have been estimated to contain thirteen cycles
of length 7 to 23 years (Figure 10).
The NAO cycle lengths were used to generate a cumulative distribution
function from which random cycle lengths were resampled, as described in
section 3.2. The random cycle lengths were scaled to cover the range of
marine survival through the use of the mean, amplitude and added noise.
Marine survival from smolt to 2 sea winter adults of Maine Atlantic salmon
has been estimated to be on the order of 0.5 - 4% based on tagging studies
(Baum, 1983) and returns of stocked hatchery smolts (USASAC, 2002). Survival
is hypothesized to be lower during the first year at sea due to the stress
of moving from freshwater to a fully marine environment. In the SalmonPVA
program, survival during the first sea winter is randomly selected from
a uniform distribution ranging from 0.08 to 0.10 and the survival during
that year for two sea winter and older fish is set to the random value
plus 0.10. This produces an expected survival of 1.71% from smolt to 2
SW adult. Combining this expected survival with the average egg to smolt
survival of 1.51% and an expected fecundity of a 2 SW female (described
below) of 7,560 eggs yields an expectation of 1.95 adults surviving from
the eggs of one female to their second sea winter. Given that Atlantic
salmon can spawn in multiple years due to kelting, approximately 2.02 spawning
adults would be expected to survive from the eggs of one 2SW female and
thus, the populations are sustainable on average if there are no human
induced removals. However, variability in survival rates at different life
stages could still cause populations to decline and small populations to
go extinct.
Variability in marine survival is produced by both the
cycles and added random noise. Amplitudes of the marine survival cycles
are randomly selected from a uniform distribution ranging from 0.02 to
0.05. The shift parameter (timing of the cycle) for the marine cycles is
randomly selected from a uniform distribution ranging from -0.3 to -0.2
to reflect the decreasing trend in marine survival observed in recent years.
The mean, amplitude, and shift values produce survival rates for smolt
to 1 SW adults of 3% to 15%, and for 2 SW adults and older of 13% to 25%.
The added noise is randomly selected from a uniform distribution ranging
from -0.04 to 0.04. The survival rate of a smolt through to the 2 SW stage
covers the range 0.05% to 4%, suggested by the working group as an appropriate
level, with a few values outside of this range (Figure
11). Note that periods of good and poor marine survival are present
in this approach and exhibit some of the same characteristics as the NAO,
such as variable period length, large sudden spikes, and changes in the
amount of variability at different times (compare Figure 10 with Figure
11).
4.3 Habitat Limitations
The base case habitat limitation in SalmonPVA operates
on the survivors of the parr 1+ stage, those fish that will be smolts in
that year. Insufficient data are available to reliably estimate habitat
limits for other juvenile life stages. The amount of juvenile rearing habitat
for parr1+ salmon in each river is calculated based upon habitat surveys
(Table 7). These surveys determined the amount
of suitable habitat based upon river characteristics such as flow, stream
width, bottom type and aquatic vegetation. Following standard practice
it was assumed that each unit of suitable juvenile rearing habitat was
capable of producing seven large parr. Thus, multiplying the number of
units of juvenile rearing habitat from the surveys by seven generates the
average habitat limitation for each river (Table 7). Variability is introduced
into the simulations through a year effect of plus or minus 20% and a river
noise effect of plus or minus 5%. These variability levels are considered
appropriate given the geographic proximity of the Maine river systems (Figure
1) and the annual variations in water levels and temperatures experienced
in recent years. One possible exception is the Dennys River in which flows
can be controlled to some extent and therefore may not be as tightly linked
to good and bad years as the other rivers.
4.4 Population Initializations
The numbers of returning Atlantic salmon for years 1995 to 2002 were derived
from either trap counts or a regression between trap counts and redd counts
(Kocik and Trial, pers. comm.). Fish counted at weirs were identified for
sea winter age, and the average sea winter age distribution from all available
returns in a year were applied to those rivers which only had redd counts
(USASAC, 2002). The 90% confidence interval from the regression between
trap counts and redds was used to derive the variability in these initializations.
The trap counts were assumed to be minimum estimates and a 10% increase,
rounded to the nearest integer, was used as the upper bound for the actual
run sizes of each sea winter age. For Cove Brook, neither weir nor redd
counts exist for the years 1995 to 1999, so an assumption was made that
the data from 2000 were applicable to these years as well. The initial
populations produced a full population structure, including kelts, in 2003,
the beginning of the projection period.
As returning fish produce natural eggs
in the year of their return, both hatchery and natural origin fish were
combined in the population initialization. Aquaculture escapees are not
included in these totals because these fish would be stopped at the weirs
and cannot be estimated at those rivers without weirs. Thus, the initial
populations are entered as returning adults, after straying, by sea winter
age with a range of possible abundance values from which a random selection
is made (Table 8). The overall trend during this
period has been declining (Figure 12).
4.5 Human Induced Removals
4.5.1 Interception Fisheries
Maine DPS Atlantic salmon while at sea can potentially be caught in the
Greenland mixed-stock fishery or in the St. Pierre et Miquelon fishery.
The Greenland fishery is currently small and captures salmon primarily
from European and Canadian rivers. Genetic analysis of salmon caught off
Greenland suggests that the number of DPS fish caught in 2002 was no greater
than one to two fish (Sheehan et al. 2003). Assuming an equally small catch
in the St. Pierre et Miquelon fishery produces a range of zero to four
fish caught. The Working Group agreed that this removal range of DPS fish
by the interception fisheries would be used as the base case for the years
2003 through 2012. The river and sea age composition of these fish depends
upon the relative size of the runs from each river and sea winter age,
implying that mainly 1SW fish will be caught.
4.5.2 Broodstock
Currently, there are no broodstock removals of returning adults from any
of the Maine DPS rivers. Instead, electrofishing for parr is used to gather
juvenile fish which are then raised in a hatchery to become broodstock.
Returning adults are removed from the Penobscot river for broodstock. The
annual number of parr1+ removed by electrofishing for broodstock has been
300 fish from both the Machias and Narraguagus River and 200 fish apiece
from the East Machias, Sheepscot, Dennys, and Pleasant Rivers. In the model,
these collections occur for years 2003 through 2014 to allow for the stocking
of the rivers (see Section 4.8).
4.5.3 River Adults
No directed river fisheries occur on any of the Maine DPS rivers. Although
illegal fishing may occur, it is totally unquantified. For the base case,
removal of river adults due to legal or illegal fishing is set to zero
for all rivers.
4.6 Straying
Straying in Atlantic salmon is thought to occur at a low rate, but is
higher for hatchery fish than wild fish (Quinn 1993). In the base case
simulations, stray rates of 1% for natural fish and 2% for hatchery fish
are assumed. No mortality is assumed to occur because of straying for either
natural or hatchery fish. To create the straying matrix, three separate
components were combined: distance between river mouths, relative river
size, and order of encounter. Each of the three components has an associated
straying matrix. These matrices are combined in a simple spreadsheet application
using weights to scale the relative importance of each factor.
The distance between river mouths straying
matrix was computed using cost weighted distances between rivers estimated
from a bathymetric scale in GIS (Marty Anderson, NOAA Fisheries, pers.
comm.). A bathymetric grid for the Gulf of Maine was obtained from the
Maine Office of GIS and stratified into three depth intervals: intertidal,
0-30 feet, and greater than 30 feet. These depth intervals were assigned
a cost, related to the perceived ability of Atlantic salmon to navigate
through these depths given the amount of obstructions expected. Overlaying
the DPS rivers on the bathymetry grid allowed the calculation of the minimum
cost path to travel from one river mouth to another using ARC VIEW GIS
Spatial Analyst. For example, in Figure 13 Cove
Brook is the starting river and the shades of color denote increasing cost
as a fish travels to the other rivers. The inverse of each cost weighted
distance was computed and the sum of the inverse distances was used to
scale each river of origin such that the sum of each row was one (Table
9a).
The straying matrix for river size is computed by taking the inverse of
the difference between river sizes and rescaling to sum to one (Table
9b). Since all the rivers in the DPS are small coastal river, the relative
differences in size are magnified. If a much larger river, such as the
Penobscot, was included in the analysis, the stray rates among the DPS
rivers would be much more similar. For this reason, river size was downweighted
relative to the other two components when constructing the overall straying
matrix.
Returning Atlantic salmon stage south of their natal river and are more
likely to stray into rivers located to the south than to the north (Collette
and Klein-MacPhee 2002). The straying matrix for river order is thus created
by assuming a linear relationship between probability of straying and order
of encounter, and then rescaling to sum to one (Table
9c).
The final base case straying matrix for SalmonPVA was created by giving
equal weight to the distance between river mouths and river order components
and only half weight to the river size component. Each matrix was multiplied
by its weighting factor and the matrices summed, element by element, and
then rescaled to sum to one (Table 9d).
4.7 Maturity and Fecundity
The maturity rate for DPS Atlantic salmon was calculated
based on the observed sea age composition of returns (USASAC 2002) and
reflects the predominance of 2SW fish (Table 10).
The proportion of male and female spawners varies by sea winter age, with
1 SW spawners mainly males and older fish more evenly split between the
sexes (Baum 1997; Table 10). The expected number of eggs produced per female
increases with sea winter age (Baum 1997; Table 10). The mean eggs per
female for kelts is an approximate average from multi-sea winter females.
The variability in numbers of eggs per female was assumed to follow a constant
coefficient of variation of 15%. The hatchery fecundity multiplier is set
to one, meaning no difference between natural and hatchery fecundity.
4.8 Stocking
In the base case simulations, the rate of stocking in
2001 was assumed to continue during the years 2000 to 2015 (USASAC 2002; Table
11). Stocked fish remain in the appropriate hatchery matrix (juvenile,
adult, or kelt) throughout their life, but their eggs get transferred to
the natural juvenile matrix (Figure 4).
5. Testing and Verification
A number of types of testing and verification were conducted for SalmonPVA,
not all of which are reported herein. The first tests compared the use
of the point estimates (mid-points of the uniform distributions described
above for the base case) for a single river, the Narraguagus, using the
full range of uncertainty in the parameters. The point estimates case should
show a self-sustaining population if the initial population size is not
too small, based on the average survival calculations described above.
The second set of tests expanded the comparison to include all eight DPS
rivers. The impact of including multiple rivers was also assessed by comparing
the results of the Narraguagus River runs to the runs with the eight DPS
rivers. The third set of tests were run multiple times with different seed
values for the random number generator to determine how many realizations
are sufficient to get a consistent result. The fourth test was a verification
approach that used historical data for the Narraguagus River to initialize
the population. This single population was followed to determine if the
model could produce a confidence interval for returning adults in years
1975 through 2002 that contained the observed point estimates. In this
verification test, the appropriate number of realizations determined from
the previous test were run and the full range of uncertainty in the base
case input parameters was used.
5.1 Single Population Parameter Point Estimates vs Range
The midpoint of each range of uncertain
parameter inputs for the Narraguagus River was simulated 1,000 times. In
none of the 1,000 simulations did the population go extinct. The average
number of natural spawners during the years 2082 to 2102 ranged from 112
to 575, with a median of 287 (Figure 14). In
all 1,000 simulations the habitat limitation constrained the number of
fish that became smolts. The habitat limitation was reached due to the
stocking of 353,000 fry in the years 2000 to 2015, but could also be reached
in later years when stocking was no longer conducted. The population was
replacing itself, with a median replacement rate of 1.057 and range of
0.7 to 1.5 (Figure 15). When the ranges of uncertainty
from the base case, described above, were input for the Narraguagus River
and simulated 1,000 times, in 89 of the 1,000 simulations the population
went extinct. The average number of natural spawners during the years 2082
to 2102 ranged from zero to 572 with a median of 36 (Figure 14). A large
shift in the frequency distribution of average spawners was induced by
the addition of uncertainty to the input parameters. Although all the uncertainties
were evenly distributed about the midpoint, the response of the number
of spawners is skewed and shifted. Decreasing the range of uncertainty
in survival rates would therefore be expected to decrease the probability
of extinction even if the average survival rate remained constant. In all
1,000 simulations, the habitat limitation was used to reduce the number
of fish that became smolts, mostly due to stocking in the early years of
the projections. In this case the population was not replacing itself,
with a median replacement rate of 0.653 and range from zero to 1.6 (Figure
15). However, some of the zero replacement rates were due to the stock
going extinct. Discounting the zero replacement rates, the median is between
0.7 and 0.8, indicative of a stock still not replacing itself.
5.2 Multiple Populations Parameter Point Estimates vs Range
The base case parameter inputs for all eight DPS rivers
were simulated 10,000 times under two scenarios; (1) the parameter point
estimates and (2) the full range of uncertainty. The output was generated
for each river separately, all eight rivers combined, and as two special
groups: Downeast rivers (DE, EM, MC, PL, and NG) and Southwest rivers (CB,
DT, and SHP). Additional groupings are possible, but were not done for
these tests. As in the single population example, inclusion of uncertainty
in the parameters increased the probability of extinction, decreased the
average numbers of spawners in years 2082 to 2102, and decreased the replacement
rate (Table 12). The estimates from the single
population examples were quite similar to those for the same river in the
multiple population examples. For example, the median replacement rates
from the single population and multiple population examples for the Narraguagus
River were within 0.02 (3%) of each other.
5.3 Determination of the Appropriate Number of Realizations
When performing Monte Carlo simulations, the appropriate number of simulations
to conduct must be found by trial and error. One method to accomplish this
is to conduct multiple trials using a different number of realizations.
For example, one might conduct ten trials of 1,000, 10,000, and 100,000
realizations and compare the variability in results among the three levels.
Conducting more simulations would produce more consistent results but take
more computation time. The tradeoff between consistency and time must be
made for each study, but once determined can usually be used for all sensitivity
analyses and changes in input parameters.
Ten trials were made using three levels of number of
realizations: 1,000, 10,000 and 100,000. Each trial used a different random
number generator seed value. Using 1,000 realizations, the ten trials had:
(a) relatively large river specific ranges (one to five percent) for the
probability of extinction; (b) one to three percent probability of habitat
limitation for those rivers not always limited; and (c) zero to seven mean
spawners for years 2082 to 2102 (Table 13a).
Given the base case settings, 1,000 simulations were insufficient to adequately
capture the amount of variability in the model. Increasing the number of
simulations to 10,000 decreased the range of the ten results to less than
two for all rivers and all three outputs described for the 1,000 simulations
(Table 13b). Further increasing the number
of simulations to 100,000 decreased the range of the ten results to less
than one in almost all cases (Table 13c). Both
the 10,000 and 100,000 simulation cases sufficiently captured the uncertainty
in results. The additional time needed to run the 100,000 simulations compared
to the 10,000 simulations does not seem to be justified by the relatively
minor decrease in variability. Thus, 10,000 simulations was used as the
standard for the base case and all further sensitivity runs.
5.4 Narraguagus River Verification
Historical numbers of adult Atlantic salmon returning
annually to the Narraguagus River, measured either through rod catches
or at adult traps, were estimated for the years 1967 to 2002 and compared
to model predictions (Table 14). The Cherryfield
adult trap was used to count returning Atlantic salmon to the Narraguagus
River from 1962 to 1974, and from 1991 to the present (Baum 1997). In the
early period fish could bypass the trap by jumping over the dam. A modification
to the spillway was made in 1991 and video monitors added to eliminate
this undercount. Approximately 25% of the returning fish in the first period
were thought to bypass the trap, so the trap records were divided by 0.75
to produce expanded trap catch estimates. Rod catch was collected during
the period 1967 to 1995 by state agencies (Baum 1997). The average ratio
of rod to expanded trap catch in the early period was applied to years
without trap catch data to estimate the number of fish expected at the
trap for years 1975 to 1990. These values represent escapement as spawners
for each year and were considered to be of the correct order of magnitude,
based on occasional redd counts. There were a few years when complete redd
counts were made, and these values were used instead of the rod to trap
ratio estimates. Ratios of hatchery and natural fish in the catch or trap
were calculated from the US Atlantic Salmon Assessment Committee database
for each year and applied to the predicted number of fish at the trap to
allow comparison with model estimates of hatchery and natural fish passing
the trap. The simulated populations were initialized based on the expanded
catch using a range of 15% to 35% bypass of the trap and then applying
the proportion at each sea winter age from the US Atlantic Salmon Assessment
Committee database.
Three sources of removals were modeled in the verification
test: (1) parr1+ for broodstock; (2) adult river fishing; and (3) the adult
interception fisheries at sea. Each year 300 parr1+ salmon were collected
by electrofishing for use as broodstock, a minor removal relative to the
parr population abundance. In contrast, the adult removals have been large
relative to abundance. In river fishing for adults was assumed to remove
between 15% and 25% of the annual returns during 1967 to 1992. Based on
catch at age information collected from the rod fishery, approximately
90% of the catch was 2 sea winter fish. This targeting of 2 sea winter
fish was modeled by assuming only 2 sea winter fish were caught by the
in river fishery. The interception fisheries were modeled by assuming that
catches occurred in proportion to relative abundance. The number of fish
caught in the interception fisheries was assumed to be equal to the run
size for years 1967 to 1975 based on tag recapture data (range of 45% to
55% of total run caught) and then was linearly reduced over time to the
current low levels (range of 0 to 10% of total run caught) (Table
15).
The Narraguagus River has been stocked with all juvenile
life stages of Atlantic salmon at some point in its history (Table
16). In the years 1967 through 1982, only smolts were stocked. During
1983 through 1991, all life stages were stocked, although smolt stocking
predominated. No stocking occurred in 1992, 1993, or 1994. Since 1995,
fry have been the predominant life stage stocked, with relatively minor
stocking of other juvenile life stages. As fry stocked fish cannot be distinguished
from natural fish based on scale readings, simulated returns of hatchery
fish for the years 1998 through 2002 were added to the natural returns
for comparison to the predicted trap data. The returns from fry stocking
in years 1985 though 1989 could not be distinguished from the larger expected
returns of parr and smolt stocked fish for these cohorts. Thus, the plots
of model natural + hatchery fry stocked returns for years 1988 through
1995 are slight underestimates and the associated plots of model hatchery
(non fry stocked) will be overestimated for this period.
The model reports the median and 80% confidence intervals of the number
of returns in each year for the natural population, hatchery population,
and total. These confidence intervals were compared to the observations
of predicted returns at the Cherryfield trap. For the verification test,
only the Narraguagus River was simulated, so no straying was included.
The annual medians
and 80% confidence intervals from 10,000 simulations are plotted with the
observed returns in Figure 16. These SalmonPVA
results did not mirror the observations of total returns; more than half
of the observed points are outside the 80% confidence intervals and the
trends are not similar. The observed data show a clear decline in total
returns from the mid 1980s to 2002 while the model predictions show an
overall increase during that period. The returns in the early part of the
series are underestimated by the model while the returns from the later
part of the series are overestimated. The early underestimation could be
due to lower interception fishery catches than assumed or higher marine
survival during this time period. Reducing the interception fishery catch
by 90% did not substantially improve the predictions for the early part
of the time series, and still produced large overestimates of more recent
returns (results not shown). A possible explanation for the overestimation
of recent returns is that stocked fish may have significantly lower survival
than natural fish. Overestimation of hatchery returns impacts not only
the hatchery population estimates in any year, but also the natural population
return estimates in subsequent years, because the offspring of hatchery
returns are classified as natural. Reducing the number of fish stocked
by 50%, but using the original interception fishery catches, reduced the
estimates of recent returns, but still produced a larger underestimation
of returns early in the series and the increasing trend in run size for
recent years remained (Figure 17). Reducing the
number of fish stocked by 90% caused the predicted returns to be quite
similar to those observed in the recent years, but still produced an increasing
trend in recent years and greatly underestimated returns in the early part
of the time series (Figure 18).
These results imply that something changed during the 1967 - 2002 period
to cause the survival rate to decrease. Early in the time series large
catches, both at sea and in the river, were supported and returns to the
trap were high. In contrast, the fisheries in recent years have been almost
completely eliminated, but the number of returns continues to decline.
Marine survival of hatchery marked smolts has declined in many rivers during
this time period. Thus, the assumption of stationarity in marine survival
is probably not applicable over the 35 year period. It is likely that there
has been a protracted downward trend in marine survival. Using SalmonPVA
with the base case inputs does not allow for non-stationarity in survival
rates. Although marine survival is cyclic, it will average the same level
at the start and end of the projection over many realizations.
While the model can predict the correct order of magnitude of returns,
the trend is not matched under the scenarios examined. The questions remain,
what has changed and what are the appropriate levels for use in forecasts?
The verification test is limited in that a model used for long-term projections
is compared to a limited set of actual observations. Even if all the model
processes and data inputs were correct, there is a chance that a single
realization will fall outside the confidence intervals. Furthermore, there
are many ways that model processes and data input can be changed to produce
similar confidence intervals. No attempts were made to tune the model to
Narraguagus River data.
6. Output for Base Case
The sections of the output file are described here with
some examples (see Figure 19). The first part
of the SalmonPVA output file is a record of the input file formatted so
that the output file can be used directly as an input file to confirm results.
This also allows confirmation of input selections to identify any input
errors.
Following the input is an optional section containing
the results from any one of the realizations, determined by user input.
This section can be quite large, approximately 5,000 lines for the base
case, but can be skipped by entering an integer less than 1 or greater
than the number of simulations. This output serves two purposes: 1) it
allows for quick troubleshooting of input when results are obviously different
from those expected and 2) it allows for a detailed evaluation of one of
the simulation trajectories. Using the same random number generator seed
value allows the user to examine multiple simulation trajectories from
a single set of results. In this section of output, the number of natural
fish at each life stage and year are reported by river and as a total over
all rivers (Table 17). The population matrix
shows how the number of fish in a cohort decreases as it ages and passes
through the different life stages. The hatchery population numbers are
presented next in the same format as the natural populations. A matrix
of habitat capacity follows with years as rows and rivers as columns. These
values are the randomly chosen habitat limitations imposed for each river
in each year. Whether or not these habitat capacity values were limiting
follows as a similar matrix with entries of 1 when habitat was limiting
and 0 when it was not. If even one year is limited for a given river in
a simulation, then that simulation is recorded as having been limited by
habitat. The next set of matrices output are the interception fishery removals
(the number of fish caught in the Greenland and St. Pierre et Miquelon
fisheries) presented by river with years as rows and classified by natural,
hatchery, and total, each of which is further stratified by sea winter
age. The next set of matrices output are the river adult fishery removals,
presented in the same format as the interception fishery removals. These
are followed by a set of matrices of natural fish removed for broodstock,
one for adult and another for parr1+. The adult matrices have years as
rows and sea winter age as columns, while the parr1+ matrices have years
as rows and stock as columns. Using this output, the populations of fish
can be tracked in conjunction with human induced removals to determine
the impact of human activity on the populations. The impact of straying
can also be derived by comparing the number of fish at each sea winter
age before and after straying, for example 1SWret and 1SWstray, to assess
the net gain or loss from each river.
Following the specific results from one realization is an optional section
that lists the occurrence of extinction or habitat limitation for each
stock and group for every realization, together with the average number
of spawners for each stock and group in every realization. This section
can quickly become huge as it grows at the rate of three times the number
of realizations. For the base case, this section would contain over 30,000
lines. This section can be used to create histograms of outputs or to determine
in which realization a certain stock went extinct. Combined with the optional
realization specific results, these outputs can show why the salmon population
in a river went extinct.
The remaining output
is always provided. The river names and components of river groups are
defined first, followed by the probability of extinction and habitat limitation
for each stock and group, along with the minimum, median, and maximum number
of years of habitat limitation for each stock and group (Table
18). These tables are followed by the average number of natural spawners
in the specified years for each stock and group and the spawner exceedence
tables (Table 19). The time series of returns
are next given as medians and 80% confidence intervals for the natural,
hatchery, and total populations (not shown here). Finally, the median replacement
rate and table of replacement rate distributions are reported (Table
20).
The output in this section contains the bottom line values of interest
from the SalmonPVA model. Here it can be seen that under the base case
assumptions, the probability of salmon becoming extinct in the next 100
years ranges from less than 10% (Narraguagus and Machias rivers) to over
50% (Cove Brook and Ducktrap). As a group, salmon in the Downeast rivers
have a 3% chance of all going extinct while the Southwest rivers have approximately
a 15% chance of going extinct. The probability of the entire DPS going
extinct in the next 100 years is less than 3%. Conversely, habitat limitations
occur frequently for all rivers, with only three rivers having less than
98% occurrence of habitat capacity limitations (Cove Brook, Pleasant, and
Ducktrap). The salmon populations in these rivers have the highest probability
of extinction, the smallest initial population sizes, and no stocking,
but are also the smallest rivers with the lowest habitat capacity. The
effect of stocking is seen in the number of years when habitat is limiting,
with the medians for stocked rivers corresponding to the approximate length
of stocking. The average number of spawners in years 2082 to 2102 reflect
this with Cove Brook, Pleasant, and Ducktrap having the lowest medians
and maximums.
Although the entire
DPS has less than 5% probability of going extinct in the next 100 years,
the time trend and replacement rate results counter this suggestion of
viability. Using only 1,000 simulations (due to computer memory limitations
with 10,000 simulations) the trajectories of returns were examined as in
the Narraguagus verification test (Figure 20, Figure 21, and Figure 22).
All the rivers and groups of rivers showed the same pattern of an increase
in returns due to stocked fish, followed by an overall decrease during
2020 to 2102. Even rivers that were not stocked, such as Cove Brook, were
impacted by the stocking of other rivers, due to straying. More importantly,
the arbitrary threshold of 5% probability of extinction in 100 years may
not be appropriate due to the continued downward trend in run sizes during
2020 to 2102. The trajectories also show the overly optimistic impact of
stocking in the base case scenario; projected returns under stocking are
much higher than currently observed at these stocking levels.
The median replacement rates for all rivers and river
groups are all below 0.75 and there is at least an 85% probability that
the replacement rate for each river and river group is less than one. This
indicates that the populations are not viable, based on the definition
provided in the Viable Salmonid Populations document (McElhany et al. 2000).
The replacement rates are mainly determined by the survival rates assumed
in the base case, but the low initial population sizes are also a decisive
factor because in a number of rivers salmon go extinct early in the projection
of many simulations, causing a replacement rate of zero (as manifest by
unusually high frequencies of zero in the replacement rate distributions)
(Figure 23). Even discounting the simulations
with zero replacement rates, the distributions reveal that base case survival
rates will not produce viable populations.
7. Sensitivity Analyses
7.1 Increased Survival
A series of five
sensitivity analyses was conducted in which survival rates of certain life
stages were increased by 20%. For each of the four juvenile life stages,
the midpoint of the survival range was increased by 20% but the range of
uncertainty held constant. This shifts the survival rates without a change
in variability. In a similar fashion, the midpoint of the average marine
survival was increased by 20% with the range remaining constant. The amplitude,
frequency, and phase of the marine cycle were not changed in the sensitivity
analyses. As expected, probability of extinction decreased, the median
value of the average spawners for years 2082 to 2102 increased, and replacement
rate increased (Table 21). The more unexpected
result was that all four juvenile life stage sensitivity results were virtually
identical. Furthermore, the magnitude of response was much greater than
20% (the amount of change in the survival rate) demonstrating the non-linear
response of the model to differences in survival rates. Changing the marine
survival rate had a larger effect than changes in any of the juvenile survival
rates for all three response variables. This is not unexpected because
the marine survival rate affects multiple years of an individual cohort,
as opposed to a single year in applying the juvenile life stage survival
rates. The increase in marine survival has a larger impact on the total
returns to the rivers than the increase in juvenile survival of any one
stage (Figure 24). However, all five sensitivity
analyses produced declining populations (as seen in the figure and in the
table of replacement rates). The consistent fluctuations in the median
returns, especially for the increased marine survival scenario, reflect
the interaction of the marine cycles with the cyclic nature of Atlantic
salmon returns; a good return year frequently produces good returns six
years later. The fluctuations are not caused by the simulated trajectories
all following the same path, as evidenced by comparing three trajectories
selected at random with the median and 80% confidence interval generated
from 10,000 simulations (Figure 25).
7.2 Stocking
When stocking did not occur in any rivers, the probability
of extinction increased to high levels for all rivers and groups of rivers
(Table 22). This demonstrates the importance
of stocking in perpetuating the currently small stocks, the more so given
that the stocking projections are optimistic, as noted above. When stocking
did not occur, the probability of habitat limitation dropped below 50%
for all eight rivers. Both the average number of spawners in the years
2080 to 2100 and the replacement rate dropped to zero for many of the rivers
due to the high probabilities of extinction. These results suggest that
if the survival rates input to the model are correct, there is a low probability
of Atlantic salmon persisting in the DPS rivers without hatchery supplementation.
7.3 Straying
The impact of straying on probability of extinction
and replacement rates was examined under two scenarios. The first sensitivity
run eliminated straying. The second sensitivity run maintained the stray
rate matrix as in the base case, but increased the straying rates for both
natural and hatchery populations five fold, from 1% and 2% to 5% and 10%,
respectively. As expected, straying caused the rivers with smaller initial
populations to have lower probabilities of extinction at the expense of
the larger rivers (Table 23). The smaller rivers
also had a higher probability of being limited by habitat when straying
occurred. The overall effect of straying is a reallocation of fish from
large populations to small populations, as seen in the DPS totals. This
component of the model is important for the rivers with small populations,
but less so for the rivers with large populations. However, this dependence
upon strays in the small river populations is not seen in genetic analyses,
where Cove Brook (the smallest of the DPS river) has consistently been
found to be the most genetically distinct river in the DPS (King et al.
2000). The recent low abundance in all the DPS rivers makes drawing a conclusion
regarding straying rates difficult.
7.4 Habitat Limitation
The sensitivity of results to the habitat limitation
was examined by changing the number of large parr per unit of habitat to
reflect potential changes in riverine productivity. The productivity of
all rivers was either decreased from 7 to 3 large parr per unit, or increased
to 11 large parr per unit, as approximate 50% changes in habitat limitation.
As expected, increasing the number of large parr that the river can hold
decreases the probability of extinction, decreases the probability of habitat
limitation, increases the average number of spawners in years 2082 to 2102,
and increases the replacement rate (Table 24).
The changes are not linear; the change from 7 parr/unit to 3 parr/unit
has a larger impact than the change from 7 parr/unit to 11 parr/unit. This
is because the habitat cap prevents Malthusian growth, and the smaller
the cap the more likely the population is to decrease back to zero. The
changes in replacement rate are due to the higher probability of extinction
because the juvenile and marine survival rates are the same among these
runs. These results show the potential for habitat restoration to decrease
the probability of extinction by allowing larger populations during periods
of good survival.
7.5 Survival Linkages Among Rivers
As noted previously, SalmonPVA allows the user to determine
the strength of the juvenile survival linkages among rivers. The full range
of possibilities include no linkages (where each river has the survival
value for that year drawn independently from all the others) through complete
linkages (where only one survival rate is drawn for a year and applied
to all rivers). These two extremes were simulated and compared to the base
case assumption of relatively strong linkages among the rivers and a noise
level of 0.05. Each scenario had the same maximum range of survival values
for a given juvenile life stage, from the minimum minus noise to maximum
plus noise. Surprisingly, the base case analysis did not, in general, produce
results intermediate to the two extreme cases (Table
25). The probability of extinction was lower for all the individual
rivers in the base case compared with both the complete linkage and complete
independence among rivers. The entire DPS probability of extinction in
the base case was intermediate to the two extremes. This apparent contradiction
between individual rivers and all rivers combined arises because all rivers
must be extinct for the DPS to become extinct. When juvenile survival rates
are completely linked among rivers, there is no difference among the rivers
in terms of good years and bad years. A period of bad years will cause
all the rivers to simultaneously go extinct because there is no population
unaffected by the poor survival rates that could otherwise provide strays.
At the other extreme, when juvenile survival is completely independent
among rivers, there are almost always some rivers with poor survival preventing
the populations in all the rivers from rebuilding together. So although
the probability of extinction over all rivers is lowered, the probability
of extinction for any given river is increased because of the loss of strong
cohorts spread over many rivers. These effects can also be seen in both
the average numbers of spawners and replacement rates.
7.6 Initial Population Abundance
The initial population abundance for each river from
1995 through 2002 was doubled to examine the sensitivity of results. Both
the lower and upper bounds of the input ranges were doubled, which also
increased the variance of the input ranges. This approach was chosen to
avoid the need to round to whole fish, as would be required by doubling
the average of the range while maintaining the spread of the input range.
Doubling the initial populations had almost no effect on: (a) the probabilities
of extinction or habitat limitation; (b) the number of average spawners
in years 2082 to 2102; or (c) the replacement rate (Table
26). The differences between the two scenarios were within the range
of noise expected using 10,000 simulations (as shown above) which explains
the counter-intuitive increases in some probabilities of extinction under
larger starting population abundances.
7.7 Sensitivity Analyses Summary
Considering both the sensitivity analyses conducted in section 7 and the
tests conducted in sections 5.1 and 5.2, the survival rates are the most
influential parameters determining probability of population persistence.
Both the level and amount of uncertainty in the survival rates determine
the equilibrium condition to which the population proceeds when projected
many years into the future. This is not surprising as the exploitation
rates on these stocks are currently thought to be close to zero and are
expected to remain low. The importance of stocking appears to be exaggerated
by the model, with more returns per stocking event than have been observed
historically. Some discounting of the amount of fish stocked therefore
seems appropriate. Straying and survival linkages among rivers are not
important for the DPS as a unit, but can be quite important for individual
rivers with small populations. Habitat capacity does cause a limiting effect
and increases the probability of extinction, but does not appear to have
a major impact in the base case conditions. The potential for habitat restoration
to improve the probability of population persistence, even under the same
survival rates, demonstrates the utility of this restoration work. Initial
population abundance is the least influential of the input parameters,
but is obviously important for rivers with very small populations.
8. Viability Analyses
8.1 Relative Population Sizes
River size is a major factor in determining the number
of Atlantic salmon that can be supported. Population viability analysis
requires that relative sizes of populations reflect differences in river
size. Given that adult returns are the most easily measured life stage,
it is desirable to conduct viability analyses in terms of numbers of returning
adults. One way to integrate river size and population size is to compute
the number of adults that would fully seed the available juvenile habitat.
These calculations have been performed for all the DPS rivers and are reported
by the US Atlantic Salmon Assessment Committee where they are denoted Conservation
Spawning Escapement (CSE) levels (USASAC 2002; Table
27). The viability analyses presented here used these values (or some
multiple), expressed as the number of 2 sea winter returns during 1995
through 2002 to initialize population sizes used in the model.
8.2 Replacement Rate
One of the conditions necessary for a viable population
is that it at least replace itself (McElhany et al. 2000). Using the CSE
levels as initial population sizes, but with no stocking, the base case
scenario results in replacement rates well below one for all rivers and
combinations of rivers (Table 28). In a highly
parameterized model, such as SalmonPVA, numerous inputs can be changed
that will affect the replacement rate. Based on the sensitivity analyses,
increases in survival rates are most likely to be required to achieve a
replacement rate of one. Since increases in each of the four juvenile life
stage survival rates produced identical percent results in the sensitivity
analyses, only one of the survival rates (the parr1+ to smolt transition)
was modified, under the assumption that a similar percentage change in
any of the other juvenile life stage survival rates would produce a similar
change in replacement rate. Average marine survival rate was also examined
to determine the amount of change in this parameter needed to generate
a replacement rate of one. The replacement rate for the entire DPS was
chosen as the DPS is the entity protected under the Endangered Species
Act. An alternative approach would be to determine the survival rates that
produced a minimum replacement rate of one for all the rivers.
The PVA results showed that the median
replacement rate for the DPS responded more quickly to relative changes
in marine survival than juvenile survival (Figure
26). A 21% increase in marine survival or a 35% increase in juvenile
survival each produced a median replacement rate of one for the DPS. The
median replacement rates for all the individual rivers were slightly below
one but the distributions were normal and centered near one (Table 28; Figure
27). Due to the effect of population abundance on the replacement rate
calculations and the effect of grouping rivers (a group of rivers has a
lower probability of extinction than any individual river in the group),
it is not possible to have all the rivers produce the same replacement
rate.
8.3 Minimum Viable Population Abundances
How small can the initial populations be and still have
a low probability of going extinct? SalmonPVA results using multiples of
CSE levels of initial population abundance (10%CSE to 50%CSE) and the marine
and juvenile survival scenarios that produced a replacement rate of 1.0
for the DPS (see above) reveal that small populations (10-20% CSE) can
have low probabilities of extinction but that the smallest rivers have
the highest probability of extinction (Table 29).
These extinction risks are true extinction, not quasi-extinction of dropping
below a certain number of returns over a given number of generations. Managers
need to decide how much risk is acceptable when setting the recovery targets
for these endangered populations.
8.4 Recovery Criteria
When setting recovery criteria for endangered species, other features
besides population abundance and replacement rates must be considered.
For Atlantic salmon, the distribution of fish among rivers is an important
attribute which mitigates against catastrophic risks and provides diversity
within the DPS. Genetic bottlenecks should also be considered in establishing
recovery criteria, although these factors cannot be addressed directly
with the current version of SalmonPVA. However, the model can be used to
guide management decisions when setting recovery criteria by depicting
the consequences of different management actions under various assumptions
about future environmental conditions.
9. Acknowledgments
Many of the decisions made while constructing the model or deciding which
data to use were made by an informal working group. This group consisted
on representatives from the National Marine Fisheries Service (NMFS; M.
Colligan, J. Hawkes, J. Kocik, M. Minton, R. Saunders, and T. Sheehan),
US Fish and Wildlife Service (FWS; C. Burger, A. Hecht, D. Kimball, J.
Marancik, M. Parkin, and D. Perkins), Maine Atlantic Salmon Commission
(ASC; K. Beland, G. Mackey, and J. Trial) and the University of Maine (M.
Kinnison). The group met as needed to discuss model development and input
data selection and reached agreements through consensus. This exercise
benefitted greatly from the many years of experience and collegial atmosphere
at these meetings. I also thank members of the US Atlantic Salmon Assessment
Committee and the ICES Working Group on North Atlantic Salmon for feedback
on this work as well as Fred Serchuk and Steve Murawski for many helpful
comments which improved this manuscript.
10. Literature Cited
Baum, E.T. 1983. The Penobscot River: an Atlantic salmon river management
report. Maine Atlantic Sea Run Commission, Bangor, Maine. 67 p.
Baum, E.T. 1997. Maine Atlantic Salmon: A National Treasure. Atlantic
Salmon Unlimited. Hermon, Maine. 224 p.
Beamish, R.J. and D.R. Bouillon. 1993. Pacific salmon production trends
in relation to climate. Canadian Journal of Fisheries and Aquatic Sciences
50: 1002-1016.
Beissinger, S.R. 2002. Population viability analysis: past, present and
future. in Beissinger, S.R. and D.R. McCullough (eds.) Population Viability
Analysis. The University of Chicago Press. Chicago.
Beissinger, S.R. and D.R. McCullough. 2002. Population Viability Analysis.
The University of Chicago Press. Chicago.
Bley, P.W. and J.R. Moring. 1988. Freshwater and ocean survival of Atlantic
salmon and steelhead: A synopsis. U.S. Fish and Wildlife Service, Biological
Report 88(9). 23 p.
Collette, B.B. and G. Kleine-MacPhee. 2002. Bigelow and Schroeder's Fishes
of the Gulf of Maine. 3rd Edition. Smithsonian Institution Press,
Washington.
Cunjak, R.A., T.D. Prowse, and D.L. Parrish. 1998. Atlantic salmon (Salmo
salar) in winter: 'The season of parr discontent? Canadian Journal
of Fisheries and Aquatic Sciences 55(Suppl.1): 161-180.
Egglishaw, H.J. and P.E. Shackley. 1973. An experiment on faster growth
of salmon (Salmo salar L.) in a Scottish stream. Journal of Fish
Biology 5: 197-204.
Egglishaw, H. and P. Shackley. 1980. Survival and growth of salmon, Salmo
salar (L.), planted in a Scottish stream. Journal of Fish Biology
16: 565-584.
Elliott, J.M. 2001. The relative role of density in the stock-recruitment
relationship of salmonids. In Stock, Recruitment and Reference Points:
Assessment and Management of Atlantic Salmon. E. Prevost and G. Chaput
(Eds). INRA, Paris. pp. 25-66.
Elson, P.F. 1957. The Role of Hatcheries in Securing Maritime Stocks of
Atlantic Salmon. Canadian Fish Culture 21.
Elson, P.F. 1975. Atlantic salmon rivers, smolt production and optimal
spawnings: an overview of natural productions. In J.R. Bohhe and L. Sochasky
(Eds.) New England Atlantic salmon restoration conference 1975. International
Atlantic Salmon Foundation Special Publication Series 6: 96-119.
Francis, R.C. and S.R. Hare. 1994. Decadal-scale regime shifts in the
large marine ecosystems of the North-east Pacific: a case for historical
science. Fisheries Oceanography 3: 279-291.
Friedland, K.D. 1998. Ocean climate influences on critical Atlantic salmon
(Salmo salar) life history events . Canadian Journal of Fisheries
and Aquatic Sciences 55 (Suppl 1): 119-130 .
Friedland, K.D., R.E. Haas, and T.F. Sheehan. 1996. Post-smolt growth,
maturation, and survival of two stocks of Atlantic salmon. Fishery Bulletin.
94: 654-663.
Friedland, K.D., D.G. Reddin, and J.F. Kocik. 1993. Marine survival of
North American and European Atlantic salmon: effects of growth and environment.
ICES Journal of Marine Science 50: 481-492.
Friedland, K.D., D.G. Reddin, J.R. McMenemy, and K.F. Drinkwater. 2003.
Multidecadal trends in North American Atlantic salmon (Salmo salar)
stocks and climate trends relevant to juvenile survival. Canadian Journal
of Fisheries and Aquatic Sciences 60: 563-583.
Gardiner, R. and P. Shackley. 1991. Stock and recruitment and inversely
density-dependent growth of salmon, Salmo salar L., in a Scottish
stream. Journal of Fish Biology 38: 691-696.
Gilpin, M.E. and M.E. Soulé. 1986. Minimum viable populations,
processes of extinction. In Soulé, M.E. (Ed.), Conservation Biology:
The Science of Scarcity and Diversity, Sinauer Associates. Sunderland,
MA, pp. 19-34.
Greenwood, J. C. 1981. Baker River: potential production of nursery habitat
for young Atlantic salmon. New Hampshire Fish and Game Dept., Proj. AFS-10-R.
Concord, New Hampshire. 28 p.
Heggenes, J. and R. Borgstroem. 1991. Effect of habitat types on survival,
spatial distribution and production of an allopatric cohort of Atlantic
salmon, Salmo salar L., under conditions of low competition. Journal
of Fish Biology 38: 267-280.
Johnson, S.L. 1988. The effects of the 1983 El Niño on Oregon's
coho (Oncorhynchus kisutch ) and chinook (O. tshawytscha )
salmon. Fisheries Research 6: 105-123.
Jordan, R. M. and K.F. Beland. 1981. Atlantic salmon spawning and evaluation
of natural spawning success. Maine Atlantic Sea-Run Salmon Commission.
Augusta, Maine. 25 p.
Kennedy, G. J. A. 1984. Factors affecting the survival and distribution
of salmon (Salmo salar L.) stocked in upland trout (Salmo
trutta L.) streams in Northern Ireland. EIFAC Technical Paper 42 suppl.
1: 227-242.
Kennedy, G. and C. Strange. 1980. Population changes after two years of
salmon (Salmo salar L.) stocking in upland trout (Salmo trutta L.)
streams. Journal of Fish Biology 17: 577-586.
Kennedy, G. and C. Strange. 1986. The effects of intra- and inter-specific
competition on the distribution of stocked juvenile Atlantic salmon, Salmo
salar L., in relation to depth and gradient in an upland trout, Salmo
trutta L., stream. Journal of Fish Biology 29: 199-214.
King, T.L., A.P. Spidle, M.S. Eackles, B.A. Lubinski, and W.B. Schill.
2000. Mitochondrial DNA diversity in North American and European Atlantic
salmon with emphasis on the Downeast rivers of Maine. Journal of Fish Biology
3: 614-630.
Knight, A. E., G. Marancik, and J.C. Greenwood. 1982. Atlantic salmon
production potential of the Mad River, New Hampshire, 1975-1980. United
States Fish and Wildlife Service. Laconia, New Hampshire. 14 p.
Martin, J.H.A. and K.A. Mitchell. 1985. Influence of sea temperature upon
the numbers of grilse and multi-sea-winter Atlantic salmon (Salmo salar)
caught in the vicinity of the River Dee (Aberdeenshire). Canadian Journal
of Fisheries and Aquatic Sciences 42: 1513-1521.
McElhany, P., M.H. Rucklelshaus, M.J. Ford, T.C. Wainwright, and E.P.
Bjorkstedt. 2000. Viable salmonid populations and the recovery of evolutionarily
significant units. NOAA Technical Memorandum NMFS-NWFSC-42. 156 p.
McMenemy, J.R. 1995. Survival of Atlantic salmon fry stocked at low density
in the West River, Vermont. North American Journal of Fisheries Management
15: 366-374.
Meister, A.L. 1962. Atlantic salmon production in Cove Brook, Maine. Transactions
of the American Fisheries Society 91: 208-212.
Mills, D.H. 1969. The survival of juvenile Atlantic salmon and brown trout
in some Scottish streams. Freshwater Salmon Fish. Res. 39. 12 p.
Myers, R.A. 1984. Demographic consequences of precocious maturation of
Atlantic salmon (Salmo salar). Canadian Journal of Fisheries and
Aquatic Sciences 41: 1349-1353.
Nickelson, T.E. and P.W. Lawson. 1998. Population viability of coho salmon, Oncorhynchus
kisutch, in Oregon coastal basins: Application of a habitat-based
life cycle model. Canadian Journal of Fisheries and Aquatic Science 55:
2383-2392.
Orciari, R.D., G.H. Leonard, D.J. Mysling, and E.C. Schluntz. 1994. Survival,
growth, and smolt production of Atlantic salmon stocked as fry in a southern
New England stream. North American Journal of Fisheries Management 14:
588-606.
Quinn, T.P. 1993. A review of homing and straying of wild and hatchery-produced
salmon. Fisheries Research 18: 29-44.
Rago, P.J. 2001. Index measures and stock assessment in Atlantic salmon.
In Stock, Recruitment and Reference Points: Assessment and Management of
Atlantic Salmon. E. Prevost and G. Chaput (Eds). INRA, Paris. pp. 137-176.
Reed, J.M, L.S. Mills, J.B. Dunning Jr., E.S. Menges, K.S. McKelvey, R.
Frye, S.R. Bessinger, M.-C. Anstett, and P. Miller. 2002. Emerging issues
in population viability analysis. Conservation Biology 16: 7-19.
Scarnecchia, D.L. 1984. Climatic and oceanic variations affecting yield
of Icelandic stocks of Atlantic salmon (Salmo salar). Canadian
Journal of Fisheries and Aquatic Sciences 41: 917-935.
Scarnecchia, D.L., A. Isaksson, and S.E. White. 1989. Effects of oceanic
variations and the West Greenland fishery on age at maturity of Icelandic
west coast stocks of Atlantic salmon (Salmo salar). Canadian Journal
of Fisheries and Aquatic Sciences 46: 16-27.
Sheehan T., C. Legault, R. Brown, A. Spidle, and T. King. 2003. Identifying
the Gulf of Maine distinct population segment component of the West Greenland
Atlantic salmon catch. ICES North Atlantic Working Group Working Paper
7.
Stewart, L. 1963. Investigations into migratory fish propagation in the
area of the Lancashire River Board. Lancaster, Barber. 80 p.
Symons, P. 1979. Estimated escapement of Atlantic salmon (Salmo salar )
for maximum smolt production in rivers of different productivity. Journal
of the Fisheries Research Board of Canada 36: 132-140.
USASAC (US Atlantic Salmon Assessment Committee). 2002. Annual Report
of the U.S. Atlantic Salmon Assessment Committee Report No. 14 - 2001 Activities.
Concord, New Hampshire. 69 p.
Whalen, K.G. 1988. Smolt production and overwinter mortality of Atlantic
salmon (Salmo salar) stocked as fry. Thesis, University of Massachusetts,
Amherst, Massachusetts, USA.
Whalen, K.G. and G.W. LaBar. 1994. Survival and growth of Atlantic salmon
(Salmo salar) fry stocked at varying densities in the White River,
Vermont. Canadian Journal of Fisheries and Aquatic Sciences 51: 2164-2169.
Whalen, K.G. and G.W. LaBar. 1998. Survival and Growth of Unfed and Fed
Atlantic Salmon Fry Stocked in a Vermont Tributary of the Connecticut River.
North American Journal of Fisheries Management 18(4): 931-935.
Whalen, K.G., D.L. Parrish, M.E. Mather, and J.R. McMenemy. 2000. Cross-tributary
analysis of parr to smolt recruitment of Atlantic salmon (Salmo salar).
Canadian Journal of Fisheries and Aquatic Sciences 57: 1607-1616.
Appendix
|
List of Working Group Participants and Affiliations |
Anderson, Marty |
NMFS |
Beland, Ken |
ASC |
Bjorkstedt, Eric |
NMFS |
Burger, Carl |
FWS |
Colligan, Mary |
NMFS |
Cooney, Tom |
NMFS |
Hawkes, Jim |
NMFS |
Hecht, Anne |
FWS |
ICES WGNAS |
various |
Kimball, Dan |
FWS |
Kinnison, Mike |
UMaine |
Kocik, John |
NMFS |
Legault, Chris |
NMFS |
Lindley, Steve |
NMFS |
Mackey, Greg |
AFS |
Marancik, Jerry |
FWS |
McElhany, Paul |
NMFS |
Minton, Mark |
NMFS |
Nichols, Henry |
ASC |
Nickerson, Paul |
FWS |
Parkin, Mary |
FWS |
Perkins, Dave |
FWS |
Saunders, Rory |
NMFS |
Scida, Pat |
NMFS |
Sheehan, Tim |
NMFS |
Trial, Joan |
ASC |
Wainwright, Tom |
NMFS |
Waples, Robin |
NMFS |
|
|
ASC = Maine Atlantic Salmon Commission
FWS = US Fish and Wildlife Service
NMFS = National Marine Fisheries Service
UMaine = University of Maine
WGNAS = Working Group on North Atlantic Salmon |