Evaluation of alternate harvesting strategies using experimental microcosms

Predation yields greater population performance: What are the contributions of density- and trait-mediated effects?

1. Population responses to extrinsic mortality can yield no change in number of survivors (compensation) or an increase in number of survivors (overcompensation) when the population is regulated by negative density-dependence. This intriguing response has been the subject of theoretical studies, but few experiments have explored how the source of extrinsic mortality affects the response. 2. This study tests abilities of three functionally diverse predators, alone and combined, to induce (over)compensation of a prey population. Larval Aedes aegypti (Diptera: Culicidae) were exposed to predation by Mesocyclops longisetus (Crustacea: Copepoda), Anopheles barberi (Diptera: Culicidae), Corethrella appendiculata (Diptera: Corethrellidae), all three in a substitutive design, or no predation. 3. The number of survivors to adulthood, female size and development time, and a composite index of performance (r') were analysed. Predator treatment did not have a significant effect on total number of survivors, nor on number of males, suggesting mortality by predation was compensatory. Predation significantly affected number of female survivors, with a trend of more females produced with predation, though no post hoc tests were significant. Predation significantly increased female development rate and r' relative to no-predator control. 4. A sensitivity analysis indicated that the change in the number of female adults produced was the largest contributing factor to the differences in r' among cohorts. While predation did not significantly increase overall production of adults, it did release survivors from density-dependent effects sufficiently to increase population performance. This study provides an empirical test of mechanisms by which predation may yield positive effects on a population of victims, a phenomenon predicted to occur across many taxa and food webs.

Threshold harvesting as a conservation or exploitation strategy in population management

Threshold harvesting removes the surplus of a population above a set threshold and takes no harvest below the threshold. This harvesting strategy is known to prevent overexploitation while obtaining higher yields than other harvesting strategies. However, the harvest taken can vary over time, including seasons of no harvest at all. While this is undesirable in fisheries or other exploitation activities, it can be an attractive feature of management strategies where removal interventions are costly and desirable only occasionally. In the presence of population fluctuations, the issue of variable harvests and population sizes becomes even more notorious. Here, we investigate the impact of threshold harvesting on the dynamics of both population size and harvests, especially in the presence of population cycles. We take into account semelparous and iteroparous life cycles, Allee effects, observation uncertainty, and demographic as well as environmental stochasticity, using generic mathematical models in discrete time. Our results show that threshold harvesting enhances multiple forms of population stability, namely persistence, constancy, resilience, and dynamic stability. We discuss plausible choices of threshold values, depending on whether the aim is resource exploitation, pest control, or the stabilization of fluctuations.

Success and failure of ecological management is highly variable in an experimental test

When managing natural systems, the importance of recognizing the role of uncertainty has been formalized as the precautionary approach. However, it is difficult to determine the role of stochasticity in the success or failure of management because there is almost always no replication; typically, only a single observation exists for a particular site or management strategy. Yet, assessing the role of stochasticity is important for providing a strong foundation for the precautionary approach, and learning from past outcomes is critical for implementing adaptive management of species or ecosystems. In addition, adaptive management relies on being able to implement a variety of strategies in order to learn-an often difficult task in natural systems. Here, we show that there is large, stochastically driven variability in success for management treatments to control an invasive species, particularly for moderate, and more feasible, management strategies. This is exactly where the precautionary approach should be important. Even when combining management strategies, we show that moderate effort in management either fails or is highly variable in its success. This variability allows some management treatments to, on average, meet their target, even when failure is probable. Our study is an important quantitative replicated experimental test of the precautionary approach and can serve as a way to understand the variability in management outcomes in natural systems which have the potential to be more variable than our tightly controlled system.

Enhancing population stability with combined adaptive limiter control and finding the optimal harvesting-restocking balance

Fluctuations in population size may have negative consequences (e.g., an increased risk of extinction or the occurrence of repeated outbreaks), and many management strategies are aimed at avoiding them by either only restocking or only harvesting the population. Two of these strategies are adaptive limiter control (ALC) and adaptive threshold harvesting (ATH). With ALC the population is controlled by only restocking and with ATH by only harvesting. We propose the strategy of combined adaptive limiter control (CALC) as the combination of ALC and ATH and study the potential advantages of CALC over ALC and ATH. We consider two different population models, namely a stochastic overcompensatory model and a host-pathogen-predator model. For the first model, our results show that the combination of restocking and harvesting under CALC improves the constancy stability of the managed populations when the harvesting and restocking intensities are high enough. Otherwise the effect is marginal or in rare cases negative. For the second model, we show that combining harvesting with restocking reduces the outbreak risk only if the harvesting intensity is low. For medium harvesting intensities the effect is marginal and for high harvesting intensities the risk of outbreaks is increased. In addition, we study the optimal harvesting-restocking balance by considering a proxy of the benefit obtained in terms of the reduction in the outbreak risk and the harvesting and restocking costs.

Finding the sweet spot: What levels of larval mortality lead to compensation or overcompensation in adult production?

Extrinsic mortality impinging on negatively density-dependent populations can result in no change in the number of survivors (compensation) or an increase (overcompensation) by releasing the population from density-dependent effects on survivorship. The relationship between the level of extrinsic mortality (i.e., percentage of mortality) and the level and likelihood of overcompensation is theoretically important, but rarely investigated. We tested the hypothesis that overcompensation occurs below a threshold value of extrinsic mortality that is related to density-dependent mortality rate, and that additive extrinsic mortality occurs above this threshold. This hypothesis predicts that survivorship vs. extrinsic mortality will: 1) be best described by a two-segmented model with a threshold; 2) have a slope >0 below the threshold; and 3) have a slope=-1 above the threshold. We also tested whether mortality imposed by real predators and random harvest have equivalent effects on adult production, and whether magnitude of overcompensation is related to species sensitivity to density-dependence. These hypotheses were tested in the container mosquitoes Aedes aegypti, A. albopictus, A. triseriatus, and Culex pipiens (Diptera: Culicidae). Cohorts of 150 larvae were exposed to random harvest of 0-70% two days after hatch or to predation by 1-3 Mesocyclops longisetus (Crustacea: Copepoda). Overcompensation occurred in A. aegypti in a pattern consistent with predictions. Aedes triseriatus showed strong overcompensation but no evidence of a threshold, whereas A. albopictus and C. pipiens had survival consistent with compensatory mortality but no evidence of a threshold. Compared to random harvest, mortality from predation yielded greater adult production in A. aegypti and A. albopictus, lesser adult production in C. pipiens, and no difference in adult production in A. triseriatus. Our results are largely consistent with our hypothesis about overcompensation, with the caveat that thresholds for additive mortality appear to occur at very high levels of extrinsic mortality. Magnitudes of overcompensation for the three Aedes were inversely related to survival in the 0% mortality treatment, consistent with our hypothesis that overcompensation is related to sensitivity to density-dependence. A broad range of extrinsic mortality levels can yield overcompensation, which may have practical implications for attempts to control pest populations.

How (co)evolution alters predator responses to increased mortality: extinction thresholds and hydra effects

Population responses to environmental change depend on both the ecological interactions between species and the evolutionary responses of all species. In this study, we explore how evolution in prey, predators, or both species affect the responses of predator populations to a sustained increase in mortality. We use an eco-evolutionary predator-prey model to explore how evolution alters the predator extinction threshold (defined as the minimum mortality rate that prevents population growth at low predator densities) and predator hydra effects (increased predator abundance in response to increased mortality). Our analysis identifies how evolutionary responses of prey and predators individually affect the predator extinction threshold and hydra effects, and how those effects are altered by interactions between the evolutionary responses. Based on our theoretical results, we predict that it is common in natural systems for evolutionary responses in one or both species to allow predators to persist at higher mortality rates than would be possible in the absence of evolution (i.e., evolution increases the predator mortality extinction threshold). We also predict that evolution-driven hydra effects occur in a minority of natural systems, but are not rare. We revisited published eco-evolutionary models and found that evolution causes hydra effects and increases the predator extinction threshold in many studies, but those effects have been overlooked. We discuss the implications of these results for species conservation, predicting population responses to environmental change, and the possibility of evolutionary rescue.

A management case study for a new commercial fishery: brine shrimp harvesting in Great Salt Lake, Utah, USA

A fishery for brine shrimp (Artemia franciscana) cysts to supply the aquaculture industry considerably expanded in the late 1980s in the Great Salt Lake, Utah, USA. With this expansion, concerns emerged in the 1990s about the fishery's sustainability, especially its impact on the abundant western North American waterbirds that use the lake and feed on brine shrimp. We track the development of management strategies using adaptive management by the Utah Division of Wildlife Resources (UDWR), which focused on the biology of the system and development of biology-based harvesting models. The models and their rationale are presented, their success in forecasting is evaluated, and implications for managing the harvest and conserving waterbirds are examined. We view this as an interesting case study because it transpired over a short time in a relatively simple system. This permitted us to clearly track management from the onset of a harvest market, through realization that the harvest had to be managed in the absence of needed biological knowledge, to the adaptive development of management strategies as biological knowledge was accumulated. The outcome illustrates the success that harvest management can attain with careful monitoring of the resource and terminating the harvest when a necessary escapement stock is attained.

Stabilization of Structured Populations via Vector Target-Oriented Control

In contrast with unstructured models, structured discrete population models have been able to fit and predict chaotic experimental data. However, most of the chaos control techniques in the literature have been designed and analyzed in a one-dimensional setting. Here, by introducing target-oriented control for discrete dynamical systems, we prove the possibility to stabilize a chosen state for a wide range of structured population models. The results are illustrated with introducing a control in the celebrated LPA model describing a flour beetle dynamics. Moreover, we show that the new control allows to stabilize periodic solutions for higher-order difference equations, such as the delayed Ricker model, for which previous target-oriented methods were not designed.

Hydra effects in discrete-time models of stable communities

A species exhibits a hydra effect when, counter-intuitively, increased mortality of the species causes an increase in its abundance. Hydra effects have been studied in many continuous time (differential equation) multispecies models, but only rarely have hydra effects been observed in or studied with discrete time (difference equation) multispecies models. In addition most discrete time theory focuses on single-species models. Thus, it is unclear what unifying characteristics determine when hydra effects arise in discrete time models. Here, using discrete time multispecies models (where total abundance is the single variable describing each population), I show that a species exhibits a hydra effect in a stable system only when fixing that species' density at its equilibrium density destabilizes the system. This general characteristic is referred to as subsystem instability. I apply this result to two-species models and identify specific mechanisms that cause hydra effects in stable communities, e.g., in host--parasitoid models, host Allee effects and saturating parasitoid functional responses can cause parasitoid hydra effects. I discuss how the general characteristic can be used to identify mechanisms causing hydra effects in communities with three or more species. I also show that the condition for hydra effects at stable equilibria implies the system is reactive (i.e., density perturbations can grow before ultimately declining). This study extends previous work on conditions for hydra effects in single-species models by identifying necessary conditions for stable systems and sufficient conditions for cyclic systems. In total, these results show that hydra effects can arise in many more communities than previously appreciated and that hydra effects were present, but unrecognized, in previously studied discrete time models.

Hydra effects in stable communities and their implications for system dynamics

A hydra effect occurs when the mean density of a species increases in response to greater mortality. We show that, in a stable multispecies system, a species exhibits a hydra effect only if maintaining that species at its equilibrium density destabilizes the system. The stability of the original system is due to the responses of the hydra-effect species to changes in the other species' densities. If that dynamical feedback is removed by fixing the density of the hydra-effect species, large changes in the community make-up (including the possibility of species extinction) can occur. This general result has several implications: (1) Hydra effects occur in a much wider variety of species and interaction webs than has previously been described, and may occur for multiple species, even in small webs; (2) conditions for hydra effects caused by predators (or diseases) often differ from those caused by other mortality factors; (3) introducing a specialist or a switching predator of a hydra-effect species often causes large changes in the community, which frequently involve extinction of other species; (4) harvest policies that attempt to maintain a constant density of a hydra-effect species may be difficult to implement, and, if successful, are likely to cause large changes in the densities of other species; and (5) trophic cascades and other indirect effects caused by predators of hydra-effect species can exhibit amplification of effects or unexpected directions of change. Although we concentrate on systems that are originally stable and models with no stage-structure or trait variation, the generality of our result suggests that similar responses to mortality will occur in many systems without these simplifying assumptions. In addition, while hydra effects are defined as responses to altered mortality, they also imply counterintuitive responses to changes in immigration and other parameters affecting population growth.

Harvested populations are more variable only in more variable environments

The interaction between environmental variation and population dynamics is of major importance, particularly for managed and economically important species, and especially given contemporary changes in climate variability. Recent analyses of exploited animal populations contested whether exploitation or environmental variation has the greatest influence on the stability of population dynamics, with consequences for variation in yield and extinction risk. Theoretical studies however have shown that harvesting can increase or decrease population variability depending on environmental variation, and requested controlled empirical studies to test predictions. Here, we use an invertebrate model species in experimental microcosms to explore the interaction between selective harvesting and environmental variation in food availability in affecting the variability of stage‐structured animal populations over 20 generations. In a constant food environment, harvesting adults had negligible impact on population variability or population size, but in the variable food environments, harvesting adults increased population variability and reduced its size. The impact of harvesting on population variability differed between proportional and threshold harvesting, between randomly and periodically varying environments, and at different points of the time series. Our study suggests that predicting the responses to selective harvesting is sensitive to the demographic structures and processes that emerge in environments with different patterns of environmental variation.

Adaptive threshold harvesting and the suppression of transients

Fluctuations in population size are in many cases undesirable, as they can induce outbreaks and extinctions or impede the optimal management of populations. We propose the strategy of adaptive threshold harvesting (ATH) to control fluctuations in population size. In this strategy, the population is harvested whenever population size has grown beyond a certain proportion in comparison to the previous generation. Taking such population increases into account, ATH intervenes also at smaller population sizes than the strategy of threshold harvesting. Moreover, ATH is the harvesting version of adaptive limiter control (ALC) that has recently been shown to stabilize population oscillations in both experiments and theoretical studies. We find that ATH has similar stabilization properties as ALC and thus offers itself as a harvesting alternative for the control of pests, exploitation of biological resources, or when restocking interventions required from ALC are unfeasible. We present numerical simulations of ATH to illustrate its performance in the presence of noise, lattice effect, and Allee effect. In addition, we propose an adjustment to both ATH and ALC that restricts interventions when control seems unnecessary, i.e. when population size is too small or too large, respectively. This adjustment cancels prolonged transients.

Stochastic life history modeling for managing regional-scale freshwater fisheries: an experimental study of brook trout

Environmental heterogeneity can combine with evolutionary responses to create very dynamic and often locally independent populations across a landscape. Such complexity creates difficulties for managers trying to conserve populations across large areas. This study develops, applies, and tests the use of stochastic life history modeling and Monte Carlo simulation to assess management scenarios related to the realities of regional fisheries management and conservation. We apply this approach to the management of recreational brook trout (Salvelinus fontinalis) fishing; an activity that can severely impact species balance, abundance, and the size structure of fish communities. Specifically, the model incorporates population-specific life-history information (e.g., growth rate, reproductive effort, and survival) to allow forecasts of the impact of various management strategies and/or changes to environmental conditions on a population's ecological characteristics (e.g., size structure, abundance, and probability of persistence). Sampling was carried out in 16 water bodies spread across four sites in Atlantic Canada. Each water body was sampled in 2005 and reassessed in 2008. This sampling had two primary objectives: (1) define a significant proportion of life-history variation of brook trout in Atlantic Canada, and (2) to test the precision and accuracy of model predictions of population responses to experimental exploitation and management changes. The model successfully predicted population responses to changes in adult survival in 12 of 13 populations having sufficient data for validation testing, while also proving to be a useful tool when engaging stakeholders regarding management options and their associated risk. We suggest that such models are cost-effective and have great potential for informing proactive management of jurisdictions with numerous and diverse populations.

A comparison of six methods for stabilizing population dynamics

Over the last two decades, several methods have been proposed for stabilizing the dynamics of biological populations. However, these methods have typically been evaluated using different population dynamics models and in the context of very different concepts of stability, which makes it difficult to compare their relative efficiencies. Moreover, since the dynamics of populations are dependent on the life-history of the species and its environment, it is conceivable that the stabilizing effects of control methods would also be affected by such factors, a complication that has typically not been investigated. In this study, we compare six different control methods with respect to their efficiency at inducing a common level of enhancement (defined as 50% increase) for two kinds of stability (constancy and persistence) under four different life-history/environment combinations. Since these methods have been analytically studied elsewhere, we concentrate on an intuitive understanding of realistic simulations incorporating noise, extinction probability and lattice effect. We show that for these six methods, even when the magnitude of stabilization attained is the same, other aspects of the dynamics like population size distribution can be very different. Consequently, correlated aspects of stability, like the amount of persistence for a given degree of constancy stability (and vice versa) or the corresponding effective population size (a measure of resistance to genetic drift) vary widely among the methods. Moreover, the number of organisms needed to be added or removed to attain similar levels of stabilization also varies for these methods, a fact that has economic implications. Finally, we compare the relative efficiencies of these methods through a composite index of various stability related measures. Our results suggest that Lower Limiter Control (LLC) seems to be the optimal method under most conditions, with the recently proposed Adaptive Limiter Control (ALC) being a close second.

Stabilizing spatially-structured populations through adaptive Limiter Control

Stabilizing the dynamics of complex, non-linear systems is a major concern across several scientific disciplines including ecology and conservation biology. Unfortunately, most methods proposed to reduce the fluctuations in chaotic systems are not applicable to real, biological populations. This is because such methods typically require detailed knowledge of system specific parameters and the ability to manipulate them in real time; conditions often not met by most real populations. Moreover, real populations are often noisy and extinction-prone, which can sometimes render such methods ineffective. Here, we investigate a control strategy, which works by perturbing the population size, and is robust to reasonable amounts of noise and extinction probability. This strategy, called the Adaptive Limiter Control (ALC), has been previously shown to increase constancy and persistence of laboratory populations and metapopulations of Drosophila melanogaster. Here, we present a detailed numerical investigation of the effects of ALC on the fluctuations and persistence of metapopulations. We show that at high migration rates, application of ALC does not require a priori information about the population growth rates. We also show that ALC can stabilize metapopulations even when applied to as low as one-tenth of the total number of subpopulations. Moreover, ALC is effective even when the subpopulations have high extinction rates: conditions under which another control algorithm had previously failed to attain stability. Importantly, ALC not only reduces the fluctuation in metapopulation sizes, but also the global extinction probability. Finally, the method is robust to moderate levels of noise in the dynamics and the carrying capacity of the environment. These results, coupled with our earlier empirical findings, establish ALC to be a strong candidate for stabilizing real biological metapopulations.

Spatial climate patterns explain negligible variation in strength of compensatory density feedbacks in birds and mammals

The use of long-term population data to separate the demographic role of climate from density-modified demographic processes has become a major topic of ecological investigation over the last two decades. Although the ecological and evolutionary mechanisms that determine the strength of density feedbacks are now well understood, the degree to which climate gradients shape those processes across taxa and broad spatial scales remains unclear. Intuitively, harsh or highly variable environmental conditions should weaken compensatory density feedbacks because populations are hypothetically unable to achieve or maintain densities at which social and trophic interactions (e.g., competition, parasitism, predation, disease) might systematically reduce population growth. Here we investigate variation in the strength of compensatory density feedback, from long-term time series of abundance over 146 species of birds and mammals, in response to spatial gradients of broad-scale temperature precipitation variables covering 97 localities in 28 countries. We use information-theoretic metrics to rank phylogenetic generalized least-squares regression models that control for sample size (time-series length) and phylogenetic non-independence. Climatic factors explained < 1% of the remaining variation in density-feedback strength across species, with the highest non-control, model-averaged effect sizes related to extreme precipitation variables. We could not link our results directly to other published studies, because ecologists use contrasting responses, predictors and statistical approaches to correlate density feedback and climate--at the expense of comparability in a macroecological context. Censuses of multiple populations within a given species, and a priori knowledge of the spatial scales at which density feedbacks interact with climate, seem to be necessary to determine cross-taxa variation in this phenomenon. Despite the availability of robust modelling tools, the appropriate data have not yet been gathered for most species, meaning that we cannot yet make any robust generalisations about how demographic feedbacks interact with climate.

Trend and fluctuations: analysis and design of population dynamics measurements in replicate ecosystems

The dynamical evolution of complex systems is often intrinsically stochastic and subject to external random forces. In order to study the resulting variability in dynamics, it is essential to make measurements on replicate systems and to separate arbitrary variation of the average dynamics of these replicates from fluctuations around the average dynamics. Here we do so for population time-series data from replicate ecosystems sharing a common average dynamics or common trend. We explain how model parameters, including the effective interactions between species and dynamical noise, can be estimated from the data and how replication reduces errors in these estimates. For this, it is essential that the model can fit a variety of average dynamics. We then show how one can judge the quality of a model, compare alternate models, and determine which combinations of parameters are poorly determined by the data. In addition we show how replicate population dynamics experiments could be designed to optimize the acquired information of interest about the systems. Our approach is illustrated on a set of time series gathered from replicate microbial closed ecosystems.

Stabilizing biological populations and metapopulations through Adaptive Limiter Control

Despite great interest in techniques for stabilizing the dynamics of biological populations and metapopulations, very few practicable methods have been developed or empirically tested. We propose an easily implementable method, Adaptive Limiter Control (ALC), for reducing the magnitude of fluctuation in population sizes and extinction frequencies and demonstrate its efficacy in stabilizing laboratory populations and metapopulations of Drosophila melanogaster. Metapopulation stability was attained through a combination of reduced size fluctuations however, and synchrony at the subpopulation level. Simulations indicated that ALC was effective over a range of maximal population growth rates, migration rates and population dynamics models. Since simulations using broadly applicable, non-species-specific models of population dynamics were able to capture most features of the experimental data, we expect our results to be applicable to a wide range of species.

Decoupling of component and ensemble density feedbacks in birds and mammals

A component density feedback represents the effect of change in population size on single demographic rates, whereas an ensemble density feedback captures that effect on the overall growth rate of a population. Given that a population's growth rate is a synthesis of the interplay of all demographic rates operating in a population, we test the hypothesis that the strength of ensemble density feedback must augment with increasing strength of component density feedback, using long-term censuses of population size, fertility, and survival rates of 109 bird and mammal populations (97 species). We found that compensatory and depensatory component feedbacks were common (each detected in approximately 50% of the demographic rates). However, component feedback strength only explained <10% of the variation in ensemble feedback strength. To explain why, we illustrate the different sources of decoupling between component and ensemble feedbacks. We argue that the management of anthropogenic impacts on populations using component feedbacks alone is ill-advised, just as managing on the basis of ensemble feedbacks without a mechanistic understanding of the contributions made by its components and environmental variability can lead to suboptimal decisions.

Strength of density feedback in census data increases from slow to fast life histories

Life-history theory predicts an increasing rate of population growth among species arranged along a continuum from slow to fast life histories. We examine the effects of this continuum on density-feedback strength estimated using long-term census data from >700 vertebrates, invertebrates, and plants. Four life-history traits (Age at first reproduction, Body size, Fertility, Longevity) were related statistically to Gompertz strength of density feedback using generalized linear mixed-effects models and multi-model inference. Life-history traits alone explained 10 to 30% of the variation in strength across species (after controlling for time-series length and phylogenetic nonindependence). Effect sizes were largest for body size in mammals and longevity in birds, and density feedback was consistently stronger for smaller-bodied and shorter-lived species. Overcompensatory density feedback (strength <-1) occurred in 20% of species, predominantly at the fast end of the life-history continuum, implying relatively high population variability. These results support the idea that life history leaves an evolutionary signal in long-term population trends as inferred from census data. Where there is a lack of detailed demographic data, broad life-history information can inform management and conservation decisions about rebound capacity from low numbers, and propensity to fluctuate, of arrays of species in areas planned for development, harvesting, protection, and population recovery.

Is hunting mortality additive or compensatory to natural mortality? Effects of experimental harvest on the survival and cause-specific mortality of willow ptarmigan

1. The effects of harvest on the annual and seasonal survival of willow ptarmigan Lagopus lagopus L. were tested in a large-scale harvest experiment. Management units were randomly assigned to one of three experimental treatments: 0%, 15% or 30% harvest. Seasonal quotas were based on the experimental treatment and estimates of bird density before the hunting season. Survival rates and hazard functions for radio-marked ptarmigan were then estimated under the competing risks of harvest and natural mortality. 2. The partially compensatory mortality hypothesis was supported: annual survival of ptarmigan was 0·54 ± 0·08 SE under 0% harvest, 0·47 ± 0·06 under 15% harvest, and was reduced to 0·30 ± 0·05 under 30% harvest. Harvest mortality increased linearly from 0·08 ± 0·05, 0·27 ± 0·05 and 0·42 ± 0·06 from 0% to 30% harvest, whereas natural mortality was 0·38 ± 0·08, 0·25 ± 0·05 and 0·28 ± 0·06 under the same treatments. 3. Realized risk of harvest mortality was 0·08-0·12 points higher than our set harvest treatments of 0-30% because birds were exposed to risk if they moved out of protected areas. The superadditive hypothesis was supported because birds in the 30% harvest treatment had higher natural mortality during winter after the hunting season. 4. Natural mortality was mainly because of raptor predation, with two seasonal peaks in fall and spring. Natural and harvest mortality coincided during early autumn with little potential for compensation during winter months. Peak risk of harvest mortality was 5× higher than natural mortality. Low natural mortality during winter suggests that most late season harvest would be additive mortality. 5. Environmental correlates of natural mortality of ptarmigan included seasonal changes in snow cover, onset of juvenile dispersal, and periods of territorial activity. Natural mortality of ptarmigan was highest during autumn movements and nesting by gyrfalcons Falco rusticolus L. Mortality was low when gyrfalcons had departed for coastal wintering sites, and during summer when ptarmigan were attending nests and broods. 6. Our experimental results have important implications for harvest management of upland gamebirds. Seasonal quotas based on proportional harvest were effective and should be set at ≤ 15% of August populations for regional management plans. Under threshold harvest of a reproductive surplus, 15% harvest would be sustainable at productivity rates ≥ 2·5 young per pair. Impacts of winter harvest could be minimized by closing the hunting season in early November or by reducing late season quotas.

Determining the functional form of density dependence: deductive approaches for consumer-resource systems having a single resource

Consumer-resource models are used to deduce the functional form of density dependence in the consumer population. A general approach to determining the form of consumer density dependence is proposed; this involves determining the equilibrium (or average) population size for a series of different harvest rates. The relationship between a consumer's mortality and its equilibrium population size is explored for several one-consumer/one-resource models. The shape of density dependence in the resource and the shape of the numerical and functional responses all tend to be "inherited" by the consumer's density dependence. Consumer-resource models suggest that density dependence will very often have both concave and convex segments, something that is impossible under the commonly used theta-logistic model. A range of consumer-resource models predicts that consumer population size often declines at a decelerating rate with mortality at low mortality rates, is insensitive to or increases with mortality over a wide range of intermediate mortalities, and declines at a rapidly accelerating rate with increased mortality when mortality is high. This has important implications for management and conservation of natural populations.

Adaptive changes in prey vulnerability shape the response of predator populations to mortality

Simple models are used to explore how adaptive changes in prey vulnerability alter the population response of their predator to increased mortality. If the mortality is an imposed harvest, the change in prey vulnerability also influences the relationship between harvest effort and yield of the predator. The models assume that different prey phenotypes share a single resource, but have different vulnerabilities to the predator. Decreased vulnerability is assumed to decrease resource consumption rate. Adaptive change may occur by phenotypic changes in the traits of a single species or by shifts in the abundances of a pair of coexisting species or morphs. The response of the predator population is influenced by the shape of the predator's functional response, the shape of resource density dependence, and the shape of the tradeoff between vulnerability and food intake in the prey. Given a linear predator functional response, adaptive prey defense tends to produce a decelerating decline in predator population size with increased mortality. Prey defense may also greatly increase the range of mortality rates that allow predator persistence. If the predator has a type-2 response with a significant handling time, adaptive prey defense may have a greater variety of effects on the predator's response to mortality, sometimes producing alternative attractors, population cycles, or increased mean predator density. Situations in which there is disruptive selection on prey defense often imply a bimodal change in yield as a function of harvesting effort, with a minimum at intermediate effort. These results argue against using single-species models of density dependent growth to manage predatory species, and illustrate the importance of incorporating anti-predator behavior into models in applied population ecology.

Community extinction patterns in coloured environments

Understanding community responses to environmental variation is a fundamental aspect of ecological research, with direct ecological, conservation and economic implications. Here, we examined the role of the magnitude, correlation and autocorrelation structures of environmental variation on species' extinction risk (ER), and the probability of actual extinction events in model competitive communities. Both ER and probability increased with increasing positive autocorrelation when species responded independently to the environment, yet both decreased with a strong correlation between species-specific responses. These results are framed in terms of the synchrony between--and magnitude of variation within--species population sizes and are explained in terms of differences in noise amplification under different conditions. The simulation results are robust to changes in the strength of interspecific density dependence, and whether noise affects density-independent or density-dependent population processes. Similar patterns arose under different ranges of noise severity when these different model assumptions were examined. We compared our results with those from an analytically derived solution, which failed to capture many features of the simulation results.

Testing the mechanisms by which source-sink dynamics alter competitive outcomes in a model system

Dispersal among sites can affect within-site competitive outcomes via source-sink dynamics. Source-sink dynamics are thought to affect competitive outcomes primarily via spatial subsidies: by redistributing individuals from sources to sinks, source-sink dynamics can alter competitive outcomes in both sources and sinks. However, dispersal also can affect competitive outcomes via demography modification, which occurs when dispersal alters the parameters governing species' per capita demographic rates. For instance, dispersal of exploitative competitors might cause extinction of some of the resources for which competition occurs, thereby altering the competition coefficients. I used protist microcosms as a model system to test whether spatial subsidies alone could explain the effects of source-sink dynamics on competitive outcomes. I examined the long-term outcome of exploitative competition among three bacterivorous ciliate protists in microcosms of high enrichment (sources) and low enrichment (sinks) in both the presence and the absence of dispersal. Dispersal altered competitive outcomes. Fitting mathematical models to the population dynamics revealed that spatial subsidies were insufficient to account for the effects of dispersal. Fitting alternative models strongly suggested that demography modification was an important determinant of competitive outcomes. These results provide the first evidence that dispersal does not simply redistribute competitors but can alter their per capita demographic rates.

Experimental simulations about the effects of overexploitation and habitat fragmentation on populations facing environmental warming

Populations of many species are dramatically declining worldwide, but the causal mechanism remains debated among different human-related threats. Coping with this uncertainty is critical to several issues about the conservation and future of biodiversity, but remains challenging due to difficulties associated with the experimental manipulation and/or isolation of the effects of such threats under field conditions. Using controlled microcosm populations, we quantified the individual and combined effects of environmental warming, overexploitation and habitat fragmentation on population persistence. Individually, each of these threats produced similar and significant population declines, which were accelerated to different degrees depending upon particular interactions. The interaction between habitat fragmentation and harvesting generated an additive decline in population size. However, both of these threats reduced population resistance causing synergistic declines in populations also facing environmental warming. Declines in population size were up to 50 times faster when all threats acted together. These results indicate that species may be facing risks of extinction higher than those anticipated from single threat analyses and suggest that all threats should be mitigated simultaneously, if current biodiversity declines are to be reversed.

Harvesting induced fluctuations: Insights from a threshold management policy

In this work, it is shown that in a deterministic context, a threshold policy can induce cyclic behavior in an otherwise exploited stable population. These dynamics ensue as a result of the combination of the degree of harvesting pressure and more protective threshold densities. Virtual equilibrium in variable structure systems plays a determinant role in this dynamical outcome.

Harvest reserves reduce extinction risk in an experimental microcosm

Overharvesting by humans threatens a substantial fraction of endangered species. Reserves have recently received enormous attention as a means of better conserving harvested resources, despite limited empirical evidence of their efficacy. We used manipulated microcosms to test whether reserves reduce extinction risk in mobile populations of harvested Tetrahymena thermophila, a ciliate. Here we show that patterns of population distribution inside and outside reserves can be accurately predicted on the basis of simple models of diffusion coupled with logistic controls on local population growth. No extinctions occurred in eight experimental trials with reserves, whereas extinction occurred in seven of eight trials without reserves, as predicted by population viability models based on stochastic population processes. These results suggest that marine reserves may be an effective means of improving long-term viability in heavily harvested fish species.

Incorporating known sources of uncertainty to determine precautionary harvests of saltwater crocodiles

It has been demonstrated repeatedly that the degree to which regulation operates and the magnitude of environmental variation in an exploited population will together dictate the type of sustainable harvest achievable. Yet typically, harvest models fail to incorporate uncertainty in the underlying dynamics of the target population by assuming a particular (unknown) form of endogenous control. We use a novel approach to estimate the sustainable yield of saltwater crocodile (Crocodylus porosus) populations from major river systems in the Northern Territory, Australia, as an example of a system with high uncertainty. We used multimodel inference to incorporate three levels of uncertainty in yield estimation: (1) uncertainty in the choice of the underlying model(s) used to describe population dynamics, (2) the error associated with the precision and bias of model parameter estimation, and (3) environmental fluctuation (process error). We demonstrate varying strength of evidence for density regulation (1.3-96.7%) for crocodiles among 19 river systems by applying a continuum of five dynamical models (density-independent with and without drift and three alternative density-dependent models) to time series of density estimates. Evidence for density dependence increased with the number of yearly transitions over which each river system was monitored. Deterministic proportional maximum sustainable yield (PMSY) models varied widely among river systems (0.042-0.611), and there was strong evidence for an increasing PMSY as support for density dependence rose. However, there was also a large discrepancy between PMSY values and those produced by the full stochastic simulation projection incorporating all forms of uncertainty, which can be explained by the contribution of process error to estimates of sustainable harvest. We also determined that a fixed-quota harvest strategy (up to 0.2K, where K is the carrying capacity) reduces population size much more rapidly than proportional harvest (the latter strategy requiring temporal monitoring of population size to adjust harvest quotas) and greatly inflates the risk of resource depletion. Using an iconic species recovering from recent extreme overexploitation to examine the potential for renewed sustainable harvest, we have demonstrated that incorporating major forms of uncertainty into a single quantitative framework provides a robust approach to modeling the dynamics of exploited populations.

Minimum viable population sizes and global extinction risk are unrelated

Theoretical and empirical work has shown that once reduced in size and geographical range, species face a considerably elevated risk of extinction. We predict minimum viable population sizes (MVP) for 1198 species based on long-term time-series data and model-averaged population dynamics simulations. The median MVP estimate was 1377 individuals (90% probability of persistence over 100 years) but the overall distribution was wide and strongly positively skewed. Factors commonly cited as correlating with extinction risk failed to predict MVP but were able to predict successfully the probability of World Conservation Union Listing. MVPs were most strongly related to local environmental variation rather than a species' intrinsic ecological and life history attributes. Further, the large variation in MVP across species is unrelated to (or at least dwarfed by) the anthropogenic threats that drive the global biodiversity crisis by causing once-abundant species to decline.