Technical Comments

Response to Comment on "Declining Wild Salmon Populations in Relation to Parasites from Farm Salmon"

See allHide authors and affiliations

Science  19 Dec 2008:
Vol. 322, Issue 5909, pp. 1790
DOI: 10.1126/science.1156578

Abstract

We evaluated the effect of sea lice (Lepeophtheirus salmonis) infestations on wild pink salmon (Oncorhynchus gorbuscha) populations in the Broughton Archipelago, British Columbia. Riddell et al. suggest that we ignored factors and selectively used data. Here, we clarify misunderstandings and provide analysis to test the strength of our conclusions.

We agree with Riddell et al. (1) that many factors affect salmon survival. In our study of the effects of sea lice, we controlled for nonlouse factors, using the stochastic Ricker model in a comparative analysis (2). The model accounts for density-dependent mortality and environmental variation in survival (3), two mortality factors not associated with sea lice. We applied the model comparatively to exposed pre-infestation, exposed infested, and unexposed populations. The exposed pre-infestation and unexposed populations share many factors affecting pink salmon population dynamics, evidenced by similar population growth rates (2) and synchronous population dynamics (4). The main difference between unexposed populations and exposed populations is that the former did not experience infestations (5, 6), whereas the latter experienced a series of infestations (5, 7, 8). Because the exposed infested populations have a depressed population growth rate (r), this indicates that sea lice infestations have driven the difference in r between exposed and unexposed populations. Further inclusion of sea lice abundance estimates improved the model fit and showed that increased louse abundance is associated with lower survival (2).

The spatial time series data we analyzed have temporal and spatial correlation, raising the issue of pseudoreplication and its impact on statistical inference. We controlled for temporal correlation, using parametric bootstrapping of the stochastic Ricker model to generate confidence intervals for r (3). Spatial correlation is known for salmon populations but occurs at a large scale that encompasses both exposed and unexposed populations, indicating synchronous regional environmental variation (4). An appropriate statistical analysis in this situation follows a matched case-control design in epidemiology (9), where subjects experiencing the same environmental variation are divided into a treatment group (exposed populations) and control group (unexposed populations). After conditioning on the shared common environmental variation, the spatial replicates can be considered independent observations on differences in population growth rate between exposed and unexposed populations. The crucial result is that growth rate r is similar for control and treatment groups before infestations and then declines significantly for the exposed group during infestations but not the unexposed group (Table 1). A potential flaw in this analysis is that the assignment to control and treatment groups was not random but rather was based on proximity to salmon farms and sea lice infestation. Although the assumption of common environmental variation is well supported, it is not impossible that a nonlouse factor changed systematically in a negative way for exposed populations but not for unexposed populations during infestation years. However, neither Riddell et al. (1) nor we have been able to identify such a confounding factor. Furthermore, the estimated effect of sea lice is consistent with predictions from other independent analyses that established and quantified the underlying mechanisms of transmission and mortality (8, 10).

Table 1.

Point estimates and 95% confidence intervals for the population growth rate, r, of pink salmon for exposed and unexposed populations before and during the sea lice infestations and the fallow. Small deviation in estimates from those in the corrected results (16, 17) are due to dividing the unexposed populations into Before Infestations, During Infestations, and Fallow groups in the model fits to get the point estimates of the parameters and also due to stochasticity in the bootstrap simulations to get the confidence intervals.

View this table:

Riddell et al. suggest that it was inappropriate to exclude the 2004 data from our analysis. These data were not excluded but rather were assigned to a separate “fallow” category in accordance with the management intervention. In spring 2003, the provincial and federal governments implemented the Pink Salmon Action Plan, which involved fallowing the primary migration route we identified in (2). This constitutes a major management intervention that had not been conducted before or replicated since. The fallowing management action was associated with a significant reduction in L. salmonis abundance in 2003 (11) and high marine survival for that cohort (12) (Table 1). Contrary to Riddell et al.'s interpretation, our analysis estimated the survival for pink salmon over their entire life cycle, not just the early out-migration period. Further fallowing interventions may help restore pink salmon, but the fallow treatment needs replication before its effects can be robustly assessed.

Riddell et al. (1) argue that the beginning of the infestation period (2001 out-migration and 2002 return) is confounded by high spawner abundance in 2000. We agree that density-dependent mortality likely contributed to the 2002 collapse, and we controlled for this in our analysis by using a density-dependent population growth model. The Ricker model shows that pink salmon are commonly in the overcompensation range where high spawner abundance leads to low returns (Fig. 1). Even in the absence of sea louse infestation, the Ricker curve predicts a drop from an average of 2.6 times historical abundance in 2000 to a mean of 0.78 times historical abundance in 2002, a decline by a factor of 3.3. The observed mean abundance in 2002 was 0.085 times historical abundance, which is a decline by a factor of 31, indicating that a density-independent factor was the primary cause of the collapse. Riddell et al. (1) are incorrect that our results depend on setting 2000 as the start of the infestation period. Annual survival during all the infestation years was mostly negative (Table 2), and excluding the 2000 to 2002 collapses from the data yields an estimated five generations to reach 99% loss (r = –0.85; 95% CI: –1.48 to –0.25 for exposed populations during infestation years), as opposed to four generations, as estimated in our original report (2).

Fig. 1.

Plot of the Ricker map n(t) = n(t – 2)exp [rbn(t – 2)] parameterized by pre-infestation (1970 to 2001) escapement data for Broughton Archipelago pink salmon [n(t) is an escapement estimate in year t normalized by its time-series mean (2)]. The dotted line is the 1:1 line showing the carrying capacity at the intersection with the Ricker curve. Overcompensation occurs where the curve has negative slope, indicating that high spawner abundance leads to low returns. Comparison with figure 3 in (2) shows that pink salmon were commonly in the overcompensatory regime, including the exposed infested group.

Table 2.

Annual survival estimates, ln[n(t)/n(t – 2)] per pink salmon population in the Broughton Archipelago during 2002 to 2006.

View this table:

Riddell et al. (1) are concerned that excluding the Glendale but including the Kakweiken in the analysis may affect the results. This criticism has already been addressed and shown to have no effect on the results (13). Following the method in (2), with Glendale included, the population growth rate for Broughton pink salmon populations during the sea lice infestations is –1.00 (95% CI: –1.52 to –0.52), and without Glendale and Kakweiken, the population growth rate is –1.23 (95% CI: –1.80 to –0.62). The population growth rate for Broughton pink salmon populations during the infestations as reported in our paper is r = –1.17 (95% CI: –1.71 to –0.59).

Riddell et al. (1) also suggest that competition with fish from the Glendale spawning channel contributed to declines in other Broughton pink salmon populations. Unfortunately, they provide no quantitative support for this hypothesis and fail to explain why it is apparent in only half the data (odd but not even years in their figure 1) and only in pre-infestation years. Figure 2 in (2) shows that the marked decline in pink salmon during the infestations is unprecedented. Figure 1 in Riddell et al. (1) indicates the Glendale populations fluctuated synchronously with the other populations in response to the infestations. If their hypothesis were correct, the Glendale populations would have fluctuated opposite to—not synchronously with—the other populations. The common factors that explain recent differences in the fluctuations of pink salmon abundance in the Broughton Archipelago relative to unexposed populations are the sea lice infestations and the fallow treatment.

Riddell et al. (1) argue from a genetic basis that extirpation of exposed pink salmon populations is acceptable because recolonization could come from other populations in the Broughton Archipelago. Our statistical sample for the Broughton Archipelago was 14 pink salmon populations. The statistical population—and hence the scope of inference—in our analysis is all the Broughton Archipelago pink salmon populations subjected to the infestations. This includes populations that are similarly exposed but excluded from analysis because of enhancement (Glendale) or many missing data points (several rivers). The conclusion that sea lice infestations have depressed wild pink salmon in the Broughton Archipelago applies generally to populations in the Broughton Archipelago (the statistical population), not only to the populations in the analysis (the statistical sample). Should populations be lost, colonization must come from outside the Broughton rather than populations within it; recovery would be more greatly compromised by genetic loss than Riddell et al. have suggested.

Referring to extirpation rather than local extinction implies that recovery is possible. This is not the case under a regime of infestations and negative population growth rates; extant populations will be lost and recovery efforts will fail. Riddell et al. supply encouraging news that 2007 escapement estimates indicate improved survival. This is only inconsistent with our analysis if 2006 infestations remained unchanged from recent non-fallow years, but Riddell et al. provide no supporting data. The impacts of salmon aquaculture on sympatric wild salmon stocks are now known to be widespread (14), and solutions are possible given sufficient political will (15). Our analysis indicates a critical sea lice threshold of r*/α = 1.3 motile lice per juvenile pink salmon, below which population declines can be reversed (2). We agree with Riddell et al. that management needs to consider many factors affecting wild salmon. Prevention of sea lice infestations of wild juvenile salmon is a management and policy option that may help restore wild salmon.

References and Notes

View Abstract

Navigate This Article