Geometry of the ideal free distribution: individual behavioural variation and annual reproductive success in aggregations of a social ungulate.

Variation in social environment can mitigate risks and rewards associated with occupying a particular patch. We aim to integrate Ideal Free Distribution (IFD) and Geometry of the Selfish Herd (GSH) to address an apparent conflict in their predictions of equal mean fitness between patches (IFD) and declining fitness benefits within a patch (GSH). We tested these hypotheses in a socio-spatial context using individual caribou that were aggregated or disaggregated during calving and varied in their annual reproductive success (ARS). We then tested individual consistency of these spatial tactics. We reveal that two socio-spatial tactics accorded similar mean ARS (IFD); however, ARS for aggregated individuals declined near the periphery (GSH). Individuals near the aggregation periphery exhibited flexibility, whereas others were consistent. The integration of classical theories through a contemporary lens of consistent individual differences provides evidence for an integrated GSH and IFD strategy that may represent an evolutionary stable state.

Spatiotemporal heterogeneity in prey abundance and vulnerability shapes the foraging tactics of an omnivore.

Prey abundance and prey vulnerability vary across space and time, but we know little about how they mediate predator-prey interactions and predator foraging tactics. To evaluate the interplay between prey abundance, prey vulnerability and predator space use, we examined patterns of black bear (Ursus americanus) predation of caribou (Rangifer tarandus) neonates in Newfoundland, Canada using data from 317 collared individuals (9 bears, 34 adult female caribou, 274 caribou calves). During the caribou calving season, we predicted that landscape features would influence calf vulnerability to bear predation, and that bears would actively hunt calves by selecting areas associated with increased calf vulnerability. Further, we hypothesized that bears would dynamically adjust their foraging tactics in response to spatiotemporal changes in calf abundance and vulnerability (collectively, calf availability). Accordingly, we expected bears to actively hunt calves when they were most abundant and vulnerable, but switch to foraging on other resources as calf availability declined. As predicted, landscape heterogeneity influenced risk of mortality, and bears displayed the strongest selection for areas where they were most likely to kill calves, which suggested they were actively hunting caribou. Initially, the per-capita rate at which bears killed calves followed a type-I functional response, but as the calving season progressed and calf vulnerability declined, kill rates dissociated from calf abundance. In support of our hypothesis, bears adjusted their foraging tactics when they were less efficient at catching calves, highlighting the influence that predation phenology may have on predator space use. Contrary to our expectations, however, bears appeared to continue to hunt caribou as calf availability declined, but switched from a tactic of selecting areas of increased calf vulnerability to a tactic that maximized encounter rates with calves. Our results reveal that generalist predators can dynamically adjust their foraging tactics over short time-scales in response to changing prey abundance and vulnerability. Further, they demonstrate the utility of integrating temporal dynamics of prey availability into investigations of predator-prey interactions, and move towards a mechanistic understanding of the dynamic foraging tactics of a large omnivore.

Phase-dependent climate-predator interactions explain three decades of variation in neonatal caribou survival.

Climate can have direct and indirect effects on population dynamics via changes in resource competition or predation risk, but this influence may be modulated by density- or phase-dependent processes. We hypothesized that for ungulates, climatic conditions close to parturition have a greater influence on the predation risk of neonates during population declines, when females are already under nutritional stress triggered by food limitation. We examined the presence of phase-dependent climate-predator (PDCP) interactions on neonatal ungulate survival by comparing spatial and temporal fluctuations in climatic conditions, cause-specific mortality and per capita resource limitation. We determined cause-specific fates of 1384 caribou (Rangifer tarandus) from 10 herds in Newfoundland, spanning more than 30 years during periods of numerical increase and decline, while exposed to predation from black bears (Ursus americanus) and coyotes (Canis latrans). We conducted Cox proportional hazards analysis for competing risks, fit as a function of weather metrics, to assess pre- and post-partum climatic influences on survival on herds in population increase and decline phases. We used cumulative incidence functions to compare temporal changes in risk from predators. Our results support our main hypothesis; when caribou populations increased, weather conditions preceding calving were the main determinants of cause-specific mortality, but when populations declined, weather conditions during calving also influenced predator-driven mortality. Cause-specific analysis showed that weather conditions can differentially affect predation risk between black bears and coyotes with specific variables increasing the risk from one species and decreasing the risk from the other. For caribou, nutritional stress appears to increase predation risk on neonates, an interaction which is exacerbated by susceptibility to climatic events. These findings support the PDCP interactions framework, where maternal body condition influences susceptibility to climate-related events and, subsequently, risk from predation.

Divergent estimates of herd-wide caribou calf survival: Ecological factors and methodological biases.

Population monitoring is a critical part of effective wildlife management, but methods are prone to biases that can hinder our ability to accurately track changes in populations through time. Calf survival plays an important role in ungulate population dynamics and can be monitored using telemetry and herd composition surveys. These methods, however, are susceptible to unrepresentative sampling and violations of the assumption of equal detectability, respectively. Here, we capitalized on 55 herd-wide estimates of woodland caribou (Rangifer tarandus caribou) calf survival in Newfoundland, Canada, using telemetry (n = 1,175 calves) and 249 herd-wide estimates of calf:cow ratios (C:C) using herd composition surveys to investigate these potential biases. These data included 17 herd-wide estimates replicated from both methods concurrently (n = 448 calves and n = 17 surveys) which we used to understand which processes and sampling biases contributed to disagreement between estimates of herd-wide calf survival. We used Cox proportional hazards models to determine whether estimates of calf mortality risk were biased by the date a calf was collared. We also used linear mixed-effects models to determine whether estimates of C:C ratios were biased by survey date and herd size. We found that calves collared later in the calving season had a higher mortality risk and that C:C tended to be higher for surveys conducted later in the autumn. When we used these relationships to modify estimates of herd-wide calf survival derived from telemetry and herd composition surveys concurrently, we found that formerly disparate estimates of woodland caribou calf survival now overlapped (within a 95% confidence interval) in a majority of cases. Our case study highlights the potential of under-appreciated biases to impact our understanding of population dynamics and suggests ways that managers can limit the influence of these biases in the two widely applied methods for estimating herd-wide survival.

SimSurvey: An R package for comparing the design and analysis of surveys by simulating spatially-correlated populations.

Populations often show complex spatial and temporal dynamics, creating challenges in designing and implementing effective surveys. Inappropriate sampling designs can potentially lead to both under-sampling (reducing precision) and over-sampling (through the extensive and potentially expensive sampling of correlated metrics). These issues can be difficult to identify and avoid in sample surveys of fish populations as they tend to be costly and comprised of multiple levels of sampling. Population estimates are therefore affected by each level of sampling as well as the pathway taken to analyze such data. Though simulations are a useful tool for exploring the efficacy of specific sampling strategies and statistical methods, there are a limited number of tools that facilitate the simulation testing of a range of sampling and analytical pathways for multi-stage survey data. Here we introduce the R package SimSurvey, which has been designed to simplify the process of simulating surveys of age-structured and spatially-distributed populations. The package allows the user to simulate age-structured populations that vary in space and time and explore the efficacy of a range of built-in or user-defined sampling protocols to reproduce the population parameters of the known population. SimSurvey also includes a function for estimating the stratified mean and variance of the population from the simulated survey data. We demonstrate the use of this package using a case study and show that it can reveal unexpected sources of bias and be used to explore design-based solutions to such problems. In summary, SimSurvey can serve as a convenient, accessible and flexible platform for simulating a wide range of sampling strategies for fish stocks and other populations that show complex structuring. Various statistical approaches can then be applied to the results to test the efficacy of different analytical approaches.

Implementing a novel movement-based approach to inferring parturition and neonate caribou calf survival.

In ungulates, parturition is correlated with a reduction in movement rate. With advances in movement-based technologies comes an opportunity to develop new techniques to assess reproduction in wild ungulates that are less invasive and reduce biases. DeMars et al. (2013, Ecology and Evolution 3:4149-4160) proposed two promising new methods (individual- and population-based; the DeMars model) that use GPS inter-fix step length of adult female caribou (Rangifer tarandus caribou) to infer parturition and neonate survival. Our objective was to apply the DeMars model to caribou populations that may violate model assumptions for retrospective analysis of parturition and calf survival. We extended the use of the DeMars model after assigning parturition and calf mortality status by examining herd-wide distributions of parturition date, calf mortality date, and survival. We used the DeMars model to estimate parturition and calf mortality events and compared them with the known parturition and calf mortality events from collared adult females (n = 19). We also used the DeMars model to estimate parturition and calf mortality events for collared female caribou with unknown parturition and calf mortality events (n = 43) and instead derived herd-wide estimates of calf survival as well as distributions of parturition and calf mortality dates and compared them to herd-wide estimates generated from calves fitted with VHF collars (n = 134). For our data, the individual-based method was effective at predicting calf mortality, but was not effective at predicting parturition. The population-based method was more effective at predicting parturition but was not effective at predicting calf mortality. At the herd-level, the predicted distributions of parturition date from both methods differed from each other and from the distribution derived from the parturition dates of VHF-collared calves (log-ranked test: χ2 = 40.5, df = 2, p < 0.01). The predicted distributions of calf mortality dates from both methods were similar to the observed distribution derived from VHF-collared calves. Both methods underestimated herd-wide calf survival based on VHF-collared calves, however, a combination of the individual- and population-based methods produced herd-wide survival estimates similar to estimates generated from collared calves. The limitations we experienced when applying the DeMars model could result from the shortcomings in our data violating model assumptions. However despite the differences in our caribou systems, with proper validation techniques the framework in the DeMars model is sufficient to make inferences on parturition and calf mortality.