Familiarity confers anti‐predator and foraging advantages on juvenile brown trout

Understanding preferences of animals is of fundamental importance for modelling habitat quality and quantity. Important theoretical developments, for example using ideal free and ideal despotic distributions (IDD), have enabled biologists to build conceptual frameworks for relating habitat preferences of individual animals to distributions and dynamics of populations. At the same time, managers of natural resources have established predictive empirical models as a basis for understanding habitat quality. For example, the Physical Habitat Simulation model (PHABSIM) has been widely applied for managing river flows. The aim of this study was to conduct experiments, using Atlantic salmon parr, to test whether observed population distributions could be predicted using simple behaviour theory and PHABSIM. We show that predictions from PHABSIM depend crucially on population density, discharge and the interaction between density and discharge at the time when the model is parameterized. These findings can, in part, be explained by consideration and application of the IDD. However, the results of the experiments also suggest that models derived from first behavioural principles may need to be unexpectedly complex and species‐specific if they are to capture the population response to variations in water discharge effectively.     

Aggregations from using inadvertent social information: a form of ideal habitat selection

Social information in breeding site selection has received extensive study; however, few attempts have been made to link this process to pre‐existing models. We examine the importance of social information to three pertinent models of habitat selection that describe breeding aggregations and spatial patterns: 1) the ideal despotic distribution (IDD) which considers conspecific competition and habitat availability, 2) the perceptual constraints model which accounts for patch selection when animals experience a threshold of undetectable difference in quality, and 3) the “neighbourhood model” which predicts concordance between resources and settlers can be disrupted by conspecific attraction when resources are patchy. These models all predict initial settlers will select a high quality patch first. However, their predictions of subsequent settlement behaviour in remaining patches differ: the IDD predicts subsequent settlers will be distributed regularly, the perceptual constraints model predicts a random distribution, and the neighbourhood model predicts clustering from conspecific attraction. We examined which model best described settlement patterns of bobolink Dolichonyx oryzivorus and savannah sparrow Passerculus sandwichensis, in the context of social information. We observed settlement timing, quantified available resources, and determined where they occurred in the highest (local population “core”) and lowest densities (local population “periphery”). We then assessed whether individuals in the periphery settled in greater concordance with resources or conspecific presence. Core territories were clustered strongly on relevant resources, and these territory holders were older than in the periphery. Peripheral territories were likewise clustered but did not always co‐occur with the best available resources, matching the neighbourhood model prediction that social information may not always direct them to the best sites available. This suggests older individuals used their own experience to locate ideal habitat, whereas younger individuals attempted to aggregate on seemingly ideal habitat by using conspecific location; such information asymmetry due to age can be viewed as an “ideal aggregative distribution”.     

Simple predictions of instream habitat model outputs for target fish populations

1. Bottom-up approaches based on individual behaviour can help to identify key variables influencing populations at larger scales. Instream habitat models have been developed to predict the consequences, for populations in stream reaches, of fish preferences for particular hydraulic conditions observed at the scale of individuals. Conventional instream habitat models (e.g. PHABSIM) predict habitat values for species or life stages in reaches, and their changes with discharge. Despite their worldwide use, they have been subject to continuing criticism and have been mainly limited to site-specific case studies.
2. We ran conventional instream habitat models in 58 French stream reaches dominated by brown trout. Using non-linear mixed effect models, we demonstrated that the outputs of instream habitat models (habitat values for three trout life stages and five other species) are predictable from average characteristics of reaches (discharge, depth, width and bed particle size).
3. Our models closely reflect variations in habitat values within-reaches (with discharge) and between-reaches. Within-reach changes are linked to the Reynolds number of reaches, while between-reach changes depend mainly on the Froude number at median daily discharge. These two dimensionless variables combine discharge, mean depth and mean width of reaches. Independent model validations showed robust model predictions that are consistent with studies of habitat values for brown trout made in larger streams from western North America.
4. Our results contribute to identifying the main hydraulic variables governing estimates of fish habitat values. They should facilitate habitat studies in multiple streams, at the basin or larger scales, while reducing their cost. They should enhance the biological validation of habitat model predictions, which remains critical.     

Habitat matching: Alternatives and implications to populations and communities

Summary I evaluate habitat matching rules based on ideal distribution models of density-dependent habitat use. Recent approaches and the ideal free continuous input matching rule on which they depend, are restricted to only those habitats that are jointly occupied across the full range of population sizes. These assumptions may often be inappropriate to field applications of habitat matching. I develop alternatives that can be applied to a wide array of ideal forms of habitat selection, including the ideal free, continuous input example. Input matching can be distinguished from assumptions of consumer-resource models and preemptive habitat use by regressions of density between paired habitats (isodars). Isodars for continuous input models should be linear on a logarithmic scale, while those for consumer-resource models should be linear on an arithmetic scale. Pre-emptive isodars can be distinguished from the others by dramatic non-linearities at both low and high densities. Field data on white-footed mice support the consumer-resource theory. Implications of the rules for population regulation and community organization are highlighted by new models that specify how the fitness of pre-emptive habitat selectors should decline with increasing density. Strong non-linearities produced by comparisons between variable and homogeneous habitats produce reversing source-sink population regulation and a new form of cyclical community dynamics. Variable habitats act as a source of emigrants at low density and a sink for immigrants at high density. Subordinate species may occupy only the variable habitat at both low and high density.     

Integrating individual behavior and population ecology: the potential for habitat-dependent population regulation in a reef fish

We used the predictions of the ideal free and ideal despoticdistributions (IFD and IDD, respectively) as a basis to evaluatethe link between spatial heterogeneity, behavior, and populationdynamics in a Caribbean coral reef fish. Juvenile three-spotdamselfish (Stegastes planifrons) were more closely aggregatedin patch reef habitat than on continuous back reef. Agonisticinteractions were more frequent but feeding rates were lowerin the patch versus the continuous reef habitat. Growth rateswere lower in patch reef habitat than on the continuous reef,but mortality rates did not differ. A separate experiment usingstandard habitat units demonstrated that the patterns observedin natural habitat were the result of the spatial distributionof the habitat patches rather than resource differences between habitats. Our results do not follow the predictions of simpleIFD or IDD models. This deviation from IFD and IDD predictionsmay be the result of a number of factors, including lack ofperfect information about habitat patches, high movement costs,and higher encounter rates of dispersed patches. Our resultsdemonstrate that behavioral interactions are an integral partof population dynamics and that it is necessary to considerthe spatial organization of the habitat in both behavioraland ecological investigations.     

The Dispersal System of a Butterfly: A Test of Source-Sink Theory Suggests the Intermediate-Scale Hypothesis

Theory predicts source-sink dynamics can occur in species with the ideal preemptive distribution but not with the ideal free distribution. Source-sink dynamics can also occur in species with passive dispersal, in which a fixed fraction of the population disperses each generation. However, in nature, dispersal often approximates random diffusion rather than ideal choices or fixed probabilities. Here, I ask which dispersal system occurred in a butterfly (Euphydryas editha) known to have source-sink dynamics. The study used 13 experimental sites, where vacant and occupied habitat patches were juxtaposed. I estimated movement during the flight season and tested hypotheses about the type of dispersal system. Ideal free and ideal preemptive models were rejected because per capita movement rates were density independent. Passive dispersal was rejected because per capita rates were related to patch area and habitat preference. The diffusion model best explained the data because it predicted both the area relationship and an odd feature of the habitat preference: immigration was not higher in preferred habitat; rather, emigration was lower. The diffusion model implied that source-sink dynamics were driven by diffusion from areas of high to low population density. Existing source-sink theory assumes fine-scale patchiness, in which animals have perfect knowledge and ease of mobility. The results from the butterfly suggest that source-sink dynamics arise at coarser spatial scales, where diffusion models apply.     

Evolutionary stability of optimal foraging: Partial preferences in the diet and patch models

In this article the patch and diet choice models of the optimal foraging theory are reanalyzed with respect to evolutionary stability of the optimal foraging strategies. In their original setting these fundamental models consider a single consumer only and the resulting fitness functions are both frequency and density independent. Such fitness functions do not allow us to apply the classical game theoretical methods to study an evolutionary stability of optimal foraging strategies for competing animals. In this article frequency and density dependent fitness functions of optimal foraging are derived by separation of time scales in an underlying population dynamical model and corresponding evolutionarily stable strategies are calculated. Contrary to the classical foraging models the results of the present article predict that partial preferences occur in optimal foraging strategies as a consequence of the ecological feedback of consumer preferences on consumer fitness. In the case of the patch occupation model these partial preferences correspond to the ideal free distribution concept while in the case of the diet choice model they correspond to the partial inclusion of the less profitable prey type in predators diet.     

A piecewise linear modeling approach for testing competing theories of habitat selection: an example with mule deer in northern winter ranges

Habitat selection fundamentally drives the distribution of organisms across landscapes; density-dependent habitat selection (DDHS) is considered a central component of ecological theories explaining habitat use and population regulation. A preponderance of DDHS theories is based on ideal distributions, such that organisms select habitat according to either the ideal free, despotic, or pre-emptive distributions. Models that can be used to simultaneously test competing DDHS theories are desirable to help improve our understanding of habitat selection. We developed hierarchical, piecewise linear models that allow for simultaneous testing of DDHS theories and accommodate densities from multiple habitats and regional populations, environmental covariates, and random effects. We demonstrate the use of these models with data on mule deer (Odocoileus hemionus) abundance and net energy costs in different snow depths within winter ranges of five regional populations in western Idaho, USA. Regional population density explained 40 % of the variation in population growth, and we found that deer were ideal free in winter ranges. Deer occupied habitats with lowest net energy costs at higher densities and at a higher rate than compared to habitats with intermediate and high energy costs. The proportion of a regional population in low energy cost habitat the previous year accounted for a significant amount of variation in population growth (17 %), demonstrating the importance of winter habitat selection in regulating deer populations. These linear models are most appropriate for empirical data collected from centralized habitat patches within the local range of a species where individuals are either year-round residents or migratory (but have already arrived from migration).     

Ideal free distributions when individuals differ in competitive ability: phenotype-limited ideal free models

A series of prospective models is developed to investigate ideal free distributions in populations where individuals differ in competitive ability. The models are of three types. In the continuous-input models, there is continuous arrival of food or mates into each habitat patch, and competitors scramble to obtain as large a share as possible. In the interference models, the prey density in a particular patch stays constant but the presence of competitors slows down the rate at which prey are captured. In the kleptoparasitism model, individuals have food or females stolen from them by competitors higher in the dominance hierarchy, and in turn steal items from subordinates. A general result of the continuous-input and interference models is that the population of competitors can be truncated between patches so that the individuals with the highest competitive ability occur in the best patches, or in the patches where competitive differences are greatest. Individuals of lowest competitive ability occur in the poorest patches or where competitive differences are least, and intermediate phenotypes are ranked between these two extremes. Thus the ideal free prediction that all individuals will achieve equal fitness will not apply. However, in continuous-input cases where competitive differences between phenotypes remain constant across patches, this solution is only neutrally stable, and forms only one element of a set of equilibrium distributions. The fact that many empirical studies of continuous-input have found approximately equal mean fitness across patches may relate to this finding. Most interference studies contradict the simple ideal free solution by having different mean intake rates across patches; this may relate to the predicted positive correlation of competitive ability with patch quality. The kleptoparasitism model usually generated continuous cycling of individuals between habitat patches, though some correlation could be found between competitive ability and patch quality.     

Density-dependent habitat selection: evaluating isoleg theory with a Lotka-Volterra model

Contemporary models of density-dependent habitat selection generally focus on long-term evolutionary consequences of intraspecific or interspecific competition and/or patterns of resource use in patchy environments. A primary goal of such studies often is to elucidate evolutionary stable strategies based on steady-state dynamics of population growth. In contrast, we developed a simulation model to explore short-term movements of interspecific competitors among fine-grained habitats of differing attributes, as might result from field manipulations of habitat quality or population densities. In this model, habitat quality is expressed in terms of mean individual fitness, represented by average per capita growth rate calculated according to the Lotka-Volterra equations describing interspecific competition. This model provides a mechanism for quantifying the effects of habitat quality, patterns of resource use and competition on distributions of individuals. Results demonstrate the heuristic value of this model in corroborating predictions derived from the ideal free distribution and isodar theory, and in generating isolegs to test the predictions of isoleg theory. Results indicate that small changes in model parameters have substantial impacts on patterns of habitat use and co-occurrence between species. The model identifies a variety of conditions under which isolegs for a given type of community organization deviate from predictions of contemporary isoleg theory, potentially expanding the universe of possible interspecific behaviors underlying the development of evolutionary stable strategies.     

The ideal free pike: 50 years of fitness-maximizing dispersal in Windermere

The ideal free distribution (IFD) theory is one of the most influential theories in evolutionary ecology. It predicts how animals ought to distribute themselves within a heterogeneous habitat in order to maximize lifetime fitness. We test the population level consequence of the IFD theory using 40-year worth data on pike (Esox lucius) living in a natural lake divided into two basins. We do so by employing empirically derived density-dependent survival, dispersal and fecundity functions in the estimation of basin-specific density-dependent fitness surfaces. The intersection of the fitness surfaces for the two basins is used for deriving expected spatial distributions of pike. Comparing the derived expected spatial distributions with 50 years data of the actual spatial distribution demonstrated that pike is ideal free distributed within the lake. In general, there was a net migration from the less productive north basin to the more productive south basin. However, a pike density-manipulation experiment imposing shifting pike density gradients between the two basins managed to switch the net migration direction and hence clearly demonstrated that the Windermere pike choose their habitat in an ideal free manner. Demonstration of ideal free habitat selection on an operational field scale like this has never been undertaken before.     

Simple predictions of instream habitat model outputs for fish habitat guilds in large streams

1. In the context of a generalised modification of hydraulic conditions in medium to large streams, modelling the impacts of stream regulation on fish communities in multiple streams is an important challenge for basic and applied freshwater ecology. Conventional instream habitat models such as PHABSIM link a hydraulic model with preference curves for various species to estimate habitat value changes with discharge in stream reaches. Despite world‐wide applications, they have been scarcely used in multiple sites with multiple species. 2. We assigned 21 size classes of European fish species to four habitat guilds (cluster analysis grouping size classes with comparable microhabitat preference curves). Then, we ran a conventional instream habitat model on 28 French stream reaches belonging to the `barbel zone', to estimate habitat values versus discharge curves for the 21 size classes. We summarised the outputs as mean habitat values for guilds, and tested if they were predictable from average characteristics of reaches (discharge, depth, width, particle size). 3. As was obtained elsewhere for populations, habitat values for guilds were strongly related to average, dimensionless characteristics of reaches. The Reynolds number of reaches, equivalent to a discharge per width unit, reflected most of the discharge‐dependent changes in habitat values (within reaches). In particular, habitat values of species preferring bank (respectively midstream) microhabitats decreased (respectively increased) with increasing Reynolds number. The Froude number at median discharge was the major predictor of reach‐dependent but discharge‐independent variations in habitat values. Habitat values of species preferring riffle versus pool or bank microhabitats were higher in reaches with high Froude numbers. These relationships were consistent with existing knowledge on the different species. 4. Such results suggest that the input variables required to estimate habitat values for fish communities can be greatly simplified, as illustrated by a general estimation of the sensitivity of species preferring midstream habitats to discharge changes in any reach. Cost‐efficient alternatives to conventional instream habitat models should facilitate their validation in multiple sites, a point that remains critical in instream habitat modelling of fish communities.     

Investigating a major assumption of predictive instream habitat models: is water velocity preference of juvenile Atlantic salmon independent of discharge?

The mean column velocity preference of juvenile Atlantic salmon Salmo salar (LF 30–55 mm) was investigated by observing their spatial pattern of habitat use in a laboratory flume while varying discharge (Q) over a 18‐fold range (Q=2·6–46·8l s^‐1). Based on 341 fish observations at three discharges (Q=2·6, 15·0 and 46·8l s^‐1), three separate velocity preference curves were developed using standard procedures. The mean column velocities measured at 0·6 depth for the fish positions at the set low, medium and high discharges had medians of 7, 10 and 24 cm s^‐1, respectively, and varied significantly between the discharges. Across the range of flows, the fish utilized mean column velocities between 0 and 56 cm s^‐1, but the three velocity preference curves differed. Differences between juvenile Atlantic salmon use of habitat, defined according to mean column velocities at different discharges, were greatest at the lower end of the available range of velocities (<20 cm s^‐1). Weighted usable area (WUA), the output of the instream flow model PHABSIM that is used to describe the available habitat at a given discharge, was calculated for the flume using the preference curves built at the three set discharges. The model was highly sensitive to differences between the three preference curves and WUA varied by up to a two‐fold difference. Furthermore, habitat‐discharge relationships derived from the three preference curves were very different. Predicted habitat losses across the modelled range of discharges varied by up to 150% depending upon which velocity preference curve was used in the model. Thus, the assumption that a single preference curve can be applied across a range of discharges is not valid and is likely to result in large errors when employing PHABSIM and other models that use similar principles.     

Dynamical facilitation of the ideal free distribution in nonideal populations

The ideal free distribution (IFD) requires that individuals can accurately perceive density‐dependent habitat quality, while failure to discern quality differences below a given perception threshold results in distributions approaching spatial uniformity. Here, we investigate the role of population growth in restoring a nonideal population to the IFD. We place a simple model of discrete patch choice under limits to the resolution by which patch quality is perceived and include population growth driven by that underlying quality. Our model follows the population's distribution through both breeding and dispersal seasons when perception limits differ in their likely influence. We demonstrate that populations of perception limited movers can approximate an IFD provided sufficient population growth; however, the emergent IFD would be temporally inconstant and correspond to reproductive events. The time to emergence of the IFD during breeding is shorter under exponential growth than under logistic growth. The IFD during early colonization of a community persists longer when more patches are available to individuals. As the population matures and dispersal becomes increasingly random, there is an oscillation in the observance of IFD, with peaks most closely approximating the IFD occurring immediately after reproductive events, and higher reproductive rates producing distributions closer to the IFD.     

Modeling ecological traps for the control of feral pigs

Ecological traps are habitat sinks that are preferred by dispersing animals but have higher mortality or reduced fecundity compared to source habitats. Theory suggests that if mortality rates are sufficiently high, then ecological traps can result in extinction. An ecological trap may be created when pest animals are controlled in one area, but not in another area of equal habitat quality, and when there is density‐dependent immigration from the high‐density uncontrolled area to the low‐density controlled area. We used a logistic population model to explore how varying the proportion of habitat controlled, control mortality rate, and strength of density‐dependent immigration for feral pigs could affect the long‐term population abundance and time to extinction. Increasing control mortality, the proportion of habitat controlled and the strength of density‐dependent immigration decreased abundance both within and outside the area controlled. At higher levels of these parameters, extinction was achieved for feral pigs. We extended the analysis with a more complex stochastic, interactive model of feral pig dynamics in the Australian rangelands to examine how the same variables as the logistic model affected long‐term abundance in the controlled and uncontrolled area and time to extinction. Compared to the logistic model of feral pig dynamics, the stochastic interactive model predicted lower abundances and extinction at lower control mortalities and proportions of habitat controlled. To improve the realism of the stochastic interactive model, we substituted fixed mortality rates with a density‐dependent control mortality function, empirically derived from helicopter shooting exercises in Australia. Compared to the stochastic interactive model with fixed mortality rates, the model with the density‐dependent control mortality function did not predict as substantial decline in abundance in controlled or uncontrolled areas or extinction for any combination of variables. These models demonstrate that pest eradication is theoretically possible without the pest being controlled throughout its range because of density‐dependent immigration into the area controlled. The stronger the density‐dependent immigration, the better the overall control in controlled and uncontrolled habitat combined. However, the stronger the density‐dependent immigration, the poorer the control in the area controlled. For feral pigs, incorporating environmental stochasticity improves the prospects for eradication, but adding a realistic density‐dependent control function eliminates these prospects.     

Process-based approach reveals directional effects of environmental factors on movement between habitats

1.?Understanding the effects of environmental factors on animal distributions is a central issue in ecology. However, movement rules inferred from distribution patterns do not reveal the processes through which animal distribution is realized. 2.?We investigated individual movement rules using a process-based approach. In experiments, coastal fish larvae (red drum, Sciaenops ocellatus) were matched with an intraspecific competitor of different sizes, and time series of habitat transition of individuals were fitted with a continuous-time Markov chain model to evaluate the effects of the presence of a competitor, behavioural interactions and habitat quality on the likelihoods of habitat transition. 3.?The process-based approach revealed that these factors did not simply act as a 'slope' between habitats that makes it easier to go in one direction and more difficult to return. Rather, these factors modify the movement rules differently depending on the directions of the movement. 4.?Individuals were less likely to enter a better habitat in the presence of a larger conspecific, more likely to shift to a poorer habitat when they received aggressive behaviour and more likely to stay in a better habitat in the presence of food. However, no effect was found on the transition intensity for moving in the opposite direction. 5.?The process-based approach to evaluating movement rules of animals allowed us to see the contrasting directional effects of different factors on the underlying movement rules used by animals, as opposed to pattern-based fitting of observed distributions. Consideration of these rules would improve the existing habitat-choice models.     

Bayesian joint models with INLA exploring marine mobile predator–prey and competitor species habitat overlap

Understanding spatial physical habitat selection driven by competition and/or predator–prey interactions of mobile marine species is a fundamental goal of spatial ecology. However, spatial counts or density data for highly mobile animals often (1) include excess zeros, (2) have spatial correlation, and (3) have highly nonlinear relationships with physical habitat variables, which results in the need for complex joint spatial models. In this paper, we test the use of Bayesian hierarchical hurdle and zero‐inflated joint models with integrated nested Laplace approximation (INLA), to fit complex joint models to spatial patterns of eight mobile marine species (grey seal, harbor seal, harbor porpoise, common guillemot, black‐legged kittiwake, northern gannet, herring, and sandeels). For each joint model, we specified nonlinear smoothed effect of physical habitat covariates and selected either competing species or predator–prey interactions. Out of a range of six ecologically important physical and biologic variables that are predicted to change with climate change and large‐scale energy extraction, we identified the most important habitat variables for each species and present the relationships between these bio/physical variables and species distributions. In particular, we found that net primary production played a significant role in determining habitat preferences of all the selected mobile marine species. We have shown that the INLA method is well‐suited for modeling spatially correlated data with excessive zeros and is an efficient approach to fit complex joint spatial models with nonlinear effects of covariates. Our approach has demonstrated its ability to define joint habitat selection for both competing and prey–predator species that can be relevant to numerous issues in the management and conservation of mobile marine species.     

Is conservation based on best available science creating an ecological trap for an imperiled lagomorph?

Habitat quality regulates fitness and population density, making it a key driver of population size. Hence, increasing habitat quality is often a primary goal of species conservation. Yet, assessments of fitness and density are difficult and costly to obtain. Therefore, species conservation often uses “best available science,” extending inferences across taxa, space, or time, and inferring habitat quality from studies of habitat selection. However, there are scenarios where habitat selection is not reflective of habitat quality, and this can lead to maladaptive management strategies. The New England cottontail (Sylvilagus transitionalis) is an imperiled shrubland obligate lagomorph whose successful recovery hinges on creation of suitable habitat. Recovery of this species is also negatively impacted by the non‐native eastern cottontail (Sylvilagus floridanus), which can competitively exclude New England cottontails from preferred habitat. Herein, we evaluate habitat quality for adult and juvenile New England and eastern cottontails using survival and density as indicators. Our findings did not support selection following an ideal free distribution by New England cottontails. Instead, selected resources, which are a target of habitat management, were associated with low survival and density and pointed to a complex trade‐off between density, survival, habitat, and the presence of eastern cottontails. Further, movement distance was inversely correlated with survival in both species, suggesting that habitat fragmentation limits the ability of cottontails to freely distribute based on habitat quality. While habitat did not directly regulate survival of juvenile cottontails, tick burden had a strong negative impact on juvenile cottontails in poor body condition. Given the complex interactions among New England cottontails, eastern cottontails, and habitat, directly assessing and accounting for factors that limit New England cottontail habitat quality in management plans is vital to their recovery. Our study demonstrates an example of management for possible ecological trap conditions via the application of incomplete knowledge.     

Sampling Considerations for Disease Surveillance in Wildlife Populations

Abstract Disease surveillance in wildlife populations involves detecting the presence of a disease, characterizing its prevalence and spread, and subsequent monitoring. A probability sample of animals selected from the population and corresponding estimators of disease prevalence and detection provide estimates with quantifiable statistical properties, but this approach is rarely used. Although wildlife scientists often assume probability sampling and random disease distributions to calculate sample sizes, convenience samples (i.e., samples of readily available animals) are typically used, and disease distributions are rarely random. We demonstrate how landscape-based simulation can be used to explore properties of estimators from convenience samples in relation to probability samples. We used simulation methods to model what is known about the habitat preferences of the wildlife population, the disease distribution, and the potential biases of the convenience-sample approach. Using chronic wasting disease in free-ranging deer (Odocoileus virginianus) as a simple illustration, we show that using probability sample designs with appropriate estimators provides unbiased surveillance parameter estimates but that the selection bias and coverage errors associated with convenience samples can lead to biased and misleading results. We also suggest practical alternatives to convenience samples that mix probability and convenience sampling. For example, a sample of land areas can be selected using a probability design that oversamples areas with larger animal populations, followed by harvesting of individual animals within sampled areas using a convenience sampling method.