Estimating ungulate recruitment and growth rates using age ratios

Woodland caribou calving fidelity: Spatial location, habitat, or both?

Individuals that isolate themselves to give birth can use more than one strategy in choosing birth sites to maximize reproductive success. Previous research has focused on the consistency in the use of the same birth-site across years (i.e., spatial fidelity), but individuals alternatively may use similar habitat conditions across years (i.e., habitat fidelity). Using GPS telemetry, we determined whether woodland caribou expressed spatial or habitat fidelity during calving, and evaluated intrinsic and extrinsic factors associated with expressing either type of fidelity. We identified 56 individuals with ≥2 putative birth events, via a movement-based model, across northern Ontario between 2010 and 2014. Individuals were classified as expressing (1) spatial fidelity by comparing sequential calving locations to a random spatial distribution of available calving locations, (2) habitat fidelity using a logistic use model compared to a null (intercept only) model, (3) no fidelity (neither criterion met), or (4) both spatial and habitat fidelity (both criteria met). Across all individuals, 37% expressed no fidelity (36 of 98), 15% expressed only spatial fidelity (15 of 99), 35% expressed only habitat fidelity (34 of 98), and 14% expressed both spatial and habitat fidelity (14 of 98). Older individuals were more likely to express spatial fidelity, whereas lower availability of upland and lowland conifer forests without linear features increased the probability an individual expressed habitat fidelity. Our results indicate that managing for caribou calving needs to consider protecting both specific, known birthing sites, but also broad-scale areas of preferred habitat for calving. Understanding the mechanisms that influence caribou expressing calving fidelity, and associated fitness costs, is crucial for the conservation of the species.

Temporal and Spatial Influences on Fawn Summer Survival in Pronghorn Populations: Management Implications from Noninvasive Monitoring

Monitoring vital rates allows managers to estimate trends in growth rates of ungulate populations. However, connecting the influence of nutrition on ungulate demography is challenging. Noninvasive sampling offers a low-cost, low-effort alternative for measuring nutritional indices, allowing for an increased understanding of the mechanistic relationships between environmental factors, nutrition, and specific population vital rates. We examined the temporal influence of intrinsic and extrinsic factors on pronghorn (Antilocapra americana) fawn recruitment. We collected fresh fecal samples from adult female pronghorn in five subpopulations spanning three sampling periods associated with critical maternal life-history stages (late gestation, early lactation, breeding season) for 2 years to investigate both intra- and interannual influences. Intrinsic factors were fecal glucocorticoid metabolites (FGMs), nutritional indices (fecal nitrogen (FN) and 2,6-diaminopimelic acid (DAPA)), and dietary composition (protein intake of forbs, graminoids, legumes, other, shrubs), while the extrinsic factor was vegetative greenness (normalized difference vegetation index (NDVI)). We found variations in DAPA, protein intake of forbs, variation in forb protein intake, and protein intake of legumes during late gestation positively influenced fawn recruitment. Fecal nitrogen during early lactation showed the strongest positive influence on the recruitment of any measured parameter. Finally, breeding season NDVI and the variation in DAPA values positively influenced the subsequent year's fawn recruitment. Our longitudinal study enabled us to investigate which parameter was most important to specific periods of fawn development and recruitment. We combined the results across five subpopulations, but interpretation and subsequent management decisions should be made at the subpopulation level such that pronghorn subpopulations with low recruitment can be positively influenced by increasing nitrogen on the landscape available to adult females during the early lactation period. As the use of noninvasive monitoring methods continues to expand, we believe our methodologies and results can be broadly applied to other ungulate monitoring programs.

Evaluating movement-based methods for estimating the frequency and timing of parturition in mule deer

BACKGROUND

Information on reproduction of harvested species such as mule deer (Odocoileus hemionus) is vital for conservation and management. Furthermore, parturition in ungulates may be detected using patterns of movement logged by GPS transmitters. Several movement-based methods have been developed to detect parturition in ungulates including the Peterson method, behavioral change point analysis (BCPA), rolling minimum convex polygons (rMCP), individual-based method (IBM), and population-based method (PBM). Our objectives were to (1) test the accuracy and the precision of each previously described method and (2) develop an improved method optimized for mule deer that incorporated aspects of the other methods.

METHODS

We determined parturition timing and status for female mule deer fitted with GPS collars and implanted with vaginal implant transmitters (VITs). We used movement patterns before and after parturition to set movement thresholds for each movement-based method. Following model training, we used location and birth date data from an external dataset to test the effectiveness of each movement-based method. Additionally, we developed a novel method for detecting parturition called the analysis of parturition indicators (API). We used two regression analyses to determine the accuracy and precision of estimates generated by each method.

RESULTS

The six methods we employed varied in accuracy, with the API, rMCP, and BCPA being most accurate. Precision also varied among methods, with the API, rMCP, and PBM generating the most precise estimates of parturition dates. The API and the rMCP performed similarly and better overall than any of the other existing methods.

CONCLUSIONS

We found that movement-based methods could be used to accurately and precisely detect parturition in mule deer. Further, we determined that the API and rMCP methods had the greatest overall success at detecting parturition in mule deer. The relative success of the API and rMCP may be attributed to the fact that both methods use home range size to detect parturition and are validated using known parturition dates of collared deer. We present the API as an efficient method of estimating birth status and timing of parturition of mule deer fitted with GPS transmitters, as well as affirm the effectiveness of a previously developed method, rMCP.

Demography of an ice-obligate mysticete in a region of rapid environmental change

Antarctic minke whales (Balaenoptera bonaerensis, AMW) are an abundant, ice-dependent species susceptible to rapid climatic changes occurring in parts of the Antarctic. Here, we used remote biopsy samples and estimates of length derived from unoccupied aircraft system (UAS) to characterize for the first time the sex ratio, maturity, and pregnancy rates of AMWs around the Western Antarctic Peninsula (WAP). DNA profiling of 82 biopsy samples (2013-2020) identified 29 individual males and 40 individual females. Blubber progesterone levels indicated 59% of all sampled females were pregnant, irrespective of maturity. When corrected for sexual maturity, the median pregnancy rate was 92.3%, indicating that most mature females become pregnant each year. We measured 68 individuals by UAS (mean = 8.04 m) and estimated that 66.5% of females were mature. This study provides the first data on the demography of AMWs along the WAP and represents the first use of non-lethal approaches to studying this species. Furthermore, these results provide baselines against which future changes in population status can be assessed in this rapidly changing marine ecosystem.

Sustainable elk harvests in Alberta with increasing predator populations

Large predators often are believed to cause declines in hunter harvests of ungulates due to direct competition for prey with hunters. In Alberta, predators of elk (Cervus elaphus), including grizzly bear (Ursus arctos), cougar (Puma concolor), and wolf (Canis lupus), have increased in recent years. We used trend analysis replicated by Wildlife Management Unit (WMU) to examine regional trends in elk harvest and hunter success. Over a 26-yr period, average harvest of elk increased by 5.46% per year for unrestricted bull and by 6.64% per year for limited-quota seasons. Also, over the same time frame, average hunter success increased by 0.2% per year for unrestricted bull and by 0.3% per year for limited-quota seasons, but no trend was detected in hunter effort (P>0.05). Our results show that increasing large-predator populations do not necessarily reduce hunter harvest of elk, and we only found evidence for this in Alberta's mountain WMUs where predation on elk calves has reduced recruitment. Furthermore, data indicate that Alberta's elk harvest management has been sustainable, i.e., hunting has continued while populations of elk have increased throughout most of the province. Wildlife agencies can justify commitments to long-term population monitoring because data allow adaptive management and can inform stakeholders on the status of populations.

The abundance and persistence of Caprinae populations

Stable or growing populations may go extinct when their sizes cannot withstand large swings in temporal variation and stochastic forces. Hence, the minimum abundance threshold defining when populations can persist without human intervention forms a key conservation parameter. We identify this threshold for many populations of Caprinae, typically threatened species lacking demographic data. Doing so helps triage conservation and management actions for threatened or harvested populations. Methodologically, we used population projection matrices and simulations, with starting abundance, recruitment, and adult female survival predicting future abundance, growth rate (λ), and population trend. We incorporated mean demographic rates representative of Caprinae populations and corresponding variances from desert bighorn sheep (Ovis canadensis nelsoni), as a proxy for Caprinae sharing similar life histories. We found a population’s minimum abundance resulting in ≤ 0.01 chance of quasi-extinction (QE; population ≤ 5 adult females) in 10 years and ≤ 0.10 QE in 30 years as 50 adult females, or 70 were translocation (removals) pursued. Discovering the threshold required 3 demographic parameters. We show, however, that monitoring populations’ relationships to this threshold requires only abundance and recruitment data. This applied approach avoids the logistical and cost hurdles in measuring female survival, making assays of population persistence more practical.

Multi-objective optimization can balance trade-offs among boreal caribou, biodiversity, and climate change objectives when conservation hotspots do not overlap

The biodiversity and climate change crises have led countries-including Canada-to commit to protect more land and inland waters and to stabilize greenhouse gas concentrations. Canada is also obligated to recover populations of at-risk species, including boreal caribou. Canada has the opportunity to expand its protected areas network to protect hotspots of high value for biodiversity and climate mitigation. However, co-occurrence of hotspots is rare. Here we ask: is it possible to expand the network to simultaneously protect areas important for boreal caribou, other species at risk, climate refugia, and carbon stores? We used linear programming to prioritize areas for protection based on these conservation objectives, and assessed how prioritization for multiple, competing objectives affected the outcome for each individual objective. Our multi-objective approach produced reasonably strong representation of value across objectives. Although trade-offs were required, the multi-objective outcome was almost always better than when we ignored one objective to maximize value for another, highlighting the risk of assuming that a plan based on one objective will also result in strong outcomes for others. Multi-objective optimization approaches could be used to plan for protected areas networks that address biodiversity and climate change objectives, even when hotspots do not co-occur.

Influence of maternity penning on the success and timing of parturition by mountain caribou (Rangifer tarandus caribou)

Invasive conservation actions that require the capture and handling of individual animals are common, but the implications for both survival and reproduction are often not studied. Across North America, most populations of woodland caribou (Rangifer tarandus caribou Gmelin, 1788) are Threatened or Endangered. Maternity penning, where pregnant females are held in an enclosure until the calf is less vulnerable to predation, is one conservation action that is designed to increase population growth. Few studies have reported the influence of maternity penning on the occurrence or timing of parturition and the implications for reproduction. We quantified parturition success and dates of penned and free-ranging caribou within the Klinse-Za population of caribou found across east-central British Columbia, Canada. Parturition dates were identified using daily observations for penned caribou (n=41) and estimated dates for free-ranging caribou (n=27) generated using statistical modelling of GPS collar data. We related parturition outcomes to a range of ecological and environmental variables. We found that the occurrence and date of parturition did not differ between penned and free-ranging caribou. For all monitored animals there was an earlier calving date during years of higher snowfall and warmer winter weather. Our results suggested that maternity penning, a potentially invasive conservation action, did not increase or decrease the probability or date of parturition for this study population.

Demographic responses of a threatened, low-density ungulate to annual variation in meteorological and phenological conditions

As global climate change progresses, wildlife management will benefit from knowledge of demographic responses to climatic variation, particularly for species already endangered by other stressors. In Canada, climate change is expected to increasingly impact populations of threatened woodland caribou (Rangifer tarandus caribou) and much focus has been placed on how a warming climate has potentially facilitated the northward expansion of apparent competitors and novel predators. Climate change, however, may also exert more direct effects on caribou populations that are not mediated by predation. These effects include meteorological changes that influence resource availability and energy expenditure. Research on other ungulates suggests that climatic variation may have minimal impact on low-density populations such as woodland caribou because per-capita resources may remain sufficient even in “bad” years. We evaluated this prediction using demographic data from 21 populations in western Canada that were monitored for various intervals between 1994 and 2015. We specifically assessed whether juvenile recruitment and adult female survival were correlated with annual variation in meteorological metrics and plant phenology. Against expectations, we found that both vital rates appeared to be influenced by annual climatic variation. Juvenile recruitment was primarily correlated with variation in phenological conditions in the year prior to birth. Adult female survival was more strongly correlated with meteorological conditions and declined during colder, more variable winters. These responses may be influenced by the life history of woodland caribou, which reside in low-productivity refugia where small climatic changes may result in changes to resources that are sufficient to elicit strong demographic effects. Across all models, explained variation in vital rates was low, suggesting that other factors had greater influence on caribou demography. Nonetheless, given the declining trajectories of many woodland caribou populations, our results highlight the increased relevance of recovery actions when adverse climatic conditions are likely to negatively affect caribou demography.

Variations, validations, degradations, and noninvasive determination of pregnancy using fecal steroid metabolites in free-ranging pronghorn

Pregnancy status is a key parameter used to assess reproductive performance of a species as it represents a starting point for measuring vital rates. Vital rates allow managers to determine trends in populations such as neonate survival and recruitment; two important factors in ungulate population growth rates. Techniques to determine pregnancy have generally involved capture and restraint of the animal to obtain blood samples for determining serum hormone levels. Non-invasive pregnancy assessment, via feces, eliminates any hazards between handler and animal, as well as removes handling-induced physiological biases. Using noninvasive fecal sampling, we conducted hormone validations, investigated pregnancy rates, and determined hormone degradation rates across five subpopulations of pronghorn (Antilocapra americana) in Idaho. Samples were collected during April-May of 2018 and 2019 from adult pronghorn of known sex and age class. Metabolites of testosterone, cortisol, 17β-estradiol, and progesterone were measured in fecal samples, and concentrations of estradiol and progesterone were examined for pregnancy determination. Average fecal progesterone metabolite (FPM) levels of pregnant females were more than double compared to levels of nonpregnant females. Fecal estrogen metabolite (FEM) levels did not differ during concurrent sampling. The largest difference in FPM levels between pregnant and nonpregnant females began on 28 April. Pregnancy determination sampling showed average FPM levels for all five subpopulations were significantly different than the nonpregnant female validation group. Nonetheless, pregnancy rates for some subpopulations lacked conclusive estimates due to early fecal sampling. Fecal glucocorticoid metabolites (FGM) levels significantly differed between pregnant females and male pronghorn, but did not differ from nonpregnant females. Degradation rates of FPM and FGM differed across days, with values for FPM from Day 1 being significantly different from all subsequent days, and after Day 9 for FGM, demonstrating the requirement of fresh samples to accurately measure hormone concentrations. We concluded that a noninvasive method to diagnosis pregnancy is possible in pronghorn via progesterone metabolites if fresh samples are collected during late gestation.

Determining Optimal Stock Density of Punjab Urial (<i>Ovis vignei punjabiensis</i>) in Captivity for Breeding, Population Growth and Reintroduction Potential

BACKGROUND AND OBJECTIVE

The Punjab urial (Ovis vignei punjabiensis) is an endangered wild sheep of Pakistan, raised in captivity with the aim of re-introduction. To date, no information is available about population trends of this species in captivity. The current study was conducted with the aim to evaluate the population trend to better guide captive breeding for improved productivity and conservation value.

MATERIALS AND METHODS

Annual population data recorded and maintained by the Wildlife Department Khyber Pakhtunkhwa, Pakistan, were used. The data were compiled and analyzed in Microsoft Excel 2010 for determining growth rates and package Growthcurver in R-version 3.5.1 was used to produce a graphical representation of the population trend.

RESULTS

The overall average annual population growth rate was rN = 0.22. Results revealed a fast initial growth rate with an average value of rN = 0.4 per year. Birth rates of bN = 0.45 for the first nine years were considerably higher than the death rates dN = 0.22 and the population increased with exponential growth. In the subsequent year, very high mortalities rates (dN = 1.2), likely attributed to the clumping of the population, resulted in the collapse of the population, leaving it in a state of unstable equilibrium.

CONCLUSION

Results support the evaluation of management data to reveal carrying capacity in captive populations, to guide and inform appropriate release of surplus animals into natural habitats.

Research priorities to fill knowledge gaps in wild boar management measures that could improve the control of African swine fever in wild boar populations

The European Commission asked EFSA to provide study designs for the investigation of four research domains (RDs) according to major gaps in knowledge identified by EFSA in a report published in 2019: (RD 1) African swine fever (ASF) epidemiology in wild boar; (RD 2) ASF transmission by vectors; (RD 3) African swine fever virus (ASFV) survival in the environment, and (RD 4) the patterns of seasonality of ASF in wild boar and domestic pigs in the EU. In this Scientific Opinion, the second RD on ASF epidemiology in wild boar is addressed. Twenty-nine research objectives were proposed by the working group and broader ASF expert networks and 23 of these research objectives met a prespecified inclusion criterion. Fourteen of these 23 research objectives met the predefined threshold for selection and so were prioritised based on the following set of criteria: (1) the impact on ASF management; (2) the feasibility or practicality to carry out the study; (3) the potential implementation of study results in practice; (4) a possible short time-frame study (< 1 year); (5) the novelty of the study; and (6) if it was a priority for risk managers. Finally, after further elimination of three of the proposed research objectives due to overlapping scope of studies published during the development of this opinion, 11 research priorities were elaborated into short research proposals, considering the potential impact on ASF management and the period of one year for the research activities.

Fall supplemental feeding increases population growth rate of an endangered caribou herd

Most woodland caribou (Rangifer tarandus caribou) populations are declining primarily because of unsustainable predation resulting from habitat-mediated apparent competition. Wolf (Canis lupus) reduction is an effective recovery option because it addresses the direct effect of predation. We considered the possibility that the indirect effects of predation might also affect caribou population dynamics by adversely affecting summer foraging behaviour. If spring and/or summer nutrition was inadequate, then supplemental feeding in fall might compensate for that limitation and contribute to population growth. Improved nutrition and therefore body condition going into winter could increase adult survival and lead to improved reproductive success the next spring. To test that hypothesis, we fed high-quality food pellets to free-ranging caribou in the Kennedy Siding caribou herd each fall for six years, starting in 2014, to see if population growth rate increased. Beginning in winter 2015–16, the Province of British Columbia began a concurrent annual program to promote caribou population increase by attempting to remove most wolves within the Kennedy Siding and the adjacent caribou herds’ ranges. To evaluate the impact of feeding, we compared lambdas before and after feeding began, and to the population trend in the adjacent Quintette herd over the subsequent four years. Supplemental feeding appeared to have an incremental effect on population growth. Population growth of the Kennedy Siding herd was higher in the year after feeding began (λ = 1.06) compared to previous years (λ = 0.91) and to the untreated Quintette herd (λ = 0.95). Average annual growth rate of the Kennedy Siding herd over the subsequent four years, where both feeding and wolf reduction occurred concurrently, was higher than in the Quintette herd where the only management action in those years was wolf reduction (λ = 1.16 vs. λ = 1.08). The higher growth rate of the Kennedy Siding herd was due to higher female survival (96.2%/yr vs. 88.9%/yr). Many caribou were in relatively poor condition in the fall. Consumption of supplemental food probably improved their nutritional status which ultimately led to population growth. Further feeding experiments on other caribou herds using an adaptive management approach would verify the effect of feeding as a population recovery tool. Our results support the recommendation that multiple management actions should be implemented to improve recovery prospects for caribou.

Trophic consequences of terrestrial eutrophication for a threatened ungulate

Changes in primary productivity have the potential to substantially alter food webs, with positive outcomes for some species and negative outcomes for others. Understanding the environmental context and species traits that give rise to these divergent outcomes is a major challenge to the generality of both theoretical and applied ecology. In aquatic systems, nutrient-mediated eutrophication has led to major declines in species diversity, motivating us to seek terrestrial analogues using a large-mammal system across 598 000 km2 of the Canadian boreal forest. These forests are undergoing some of the most rapid rates of land-use change on Earth and are home to declining caribou (Rangifer tarandus caribou) populations. Using satellite-derived estimates of primary productivity, coupled with estimates of moose (Alces alces) and wolf (Canis lupus) abundance, we used path analyses to discriminate among hypotheses explaining how habitat alteration can affect caribou population growth. Hypotheses included food limitation, resource dominance by moose over caribou, and apparent competition with predators shared between moose and caribou. Results support apparent competition and yield estimates of wolf densities (1.8 individuals 1000 km-2) above which caribou populations decline. Our multi-trophic analysis provides insight into the cascading effects of habitat alteration from forest cutting that destabilize terrestrial predator-prey dynamics. Finally, the path analysis highlights why conservation actions directed at the proximate cause of caribou decline have been more successful in the near term than those directed further along the trophic chain.

Line-transect data may not produce reliable estimates of interannual sex-ratio and age structure variation in West African savannah ungulates

Adult sex ratios and age structures are important wildlife population parameters, but they have been poorly investigated in ungulate species in West African savannahs. We used line transects to investigate these parameters in 11 ungulates from a protected area in south-western Burkina Faso during the period 2010-2018. We created an empirical model of “detectability” for each species based on its main ecological characteristics (habitat and group size) and body size, and then compared the observed interannual inconsistency in sex ratios and age structures with the a priori detectability score. Six out of 11 species showed low interannual inconsistency in sex ratio and age structure. In 82% of the study species, however, the predicted detectability score matched the observed score, with two exceptions being Tragelaphus scriptus and Sincerus caffer.

Comparing survey and multiple recruitment-mortality models to assess growth rates and population projections

Estimation of population trends and demographic parameters is important to our understanding of fundamental ecology and species management, yet these data are often difficult to obtain without the use of data from population surveys or marking animals. The northeastern Minnesota moose (Alces alces Linnaeus, 1758) population declined 58% during 2006-2017, yet aerial surveys indicated stability during 2012-2017. In response to the decline, the Minnesota Department of Natural Resources (MNDNR) initiated studies of adult and calf survival to better understand cause-specific mortality, calf recruitment, and factors influencing the population trajectory. We estimated population growth rate (λ) using adult survival and calf recruitment data from demographic studies and the recruitment-mortality (R-M) Equation and compared these estimates to those calculated using data from aerial surveys. We then projected population dynamics 50 years using each resulting λ and used a stochastic model to project population dynamics 30 years using data from the MNDNR's studies. Calculations of λ derived from 2012 to 2017 survey data, and the R-M Equation indicated growth (1.02 ± 0.16 [SE] and 1.01 ± 0.04, respectively). However, the stochastic model indicated a decline in the population over 30 years (λ = 0.91 ± 0.004; 2014-2044). The R-M Equation has utility for estimating λ, and the supporting information from demographic collaring studies also helps to better address management questions. Furthermore, estimates of λ calculated using collaring data were more certain and reflective of current conditions. Long-term monitoring using collars would better inform population performance predictions and demographic responses to environmental variability.

Web-based application for threatened woodland caribou population modeling

Woodland caribou (Rangifer tarandus caribou) are threatened in Canada, with population and distribution declines evident in most regions of the country. Causes of declines are linked to landscape change from forest fires and human development, notably forestry oil and gas activities, which result in caribou habitat loss and affect ecosystem food webs. The Federal Species at Risk Act requires effective protection and restoration of caribou habitat, with actions to increase caribou survival. These requirements call for effective monitoring of caribou population trends to gauge success. Many woodland caribou populations are nearly impossible to count using traditional aerial survey methods, but demographic‐based monitoring approaches can be used to estimate population trends based on population modeling of vital rates from marked animals. Monitoring programs have used a well‐known simple population model (the Recruitment‐Mortality [R/M] equation) to estimate demographic rates for woodland caribou, but have faced challenges in managing large data streams and providing transparency in the demographic estimation process. We present a stand‐alone statistical software application using open‐source software to permit efficient, transparent, and replicable demographic estimation for woodland caribou populations. We developed an easy‐to‐use, interactive web‐based application for the R/M population model that uses a Bayesian estimation approach and provides the user flexibility in choice of prior distributions and other output features. We illustrate the web‐application to the A la Pêche Southern Mountain (Central Group) woodland caribou population in west‐central Alberta, Canada, during 1998–2017. Our estimates of population demographics are consistent with previous research on this population and highlight the utility of the application in assessing caribou population responses to species recovery actions. We provide example data, computer code, the web‐based application package, and a user manual to guide installation and use. We also review underlying assumptions and challenges of population monitoring in this case study. We expect our software will contribute to efficient monitoring of woodland caribou and help in the assessment of recovery actions for this species. © 2019 The Authors. Wildlife Society Bulletin Published by Wiley Periodicals, Inc.

We developed a stand‐alone Web‐Application to support population trend estimation for endangered and threatened woodland caribou populations.

Saving endangered species using adaptive management

Adaptive management is a powerful means of learning about complex ecosystems, but is rarely used for recovering endangered species. Here, we demonstrate how it can benefit woodland caribou, which became the first large mammal extirpated from the contiguous United States in recent history. The continental scale of forest alteration and extended time needed for forest recovery means that relying only on habitat protection and restoration will likely fail. Therefore, population management is also needed as an emergency measure to avoid further extirpation. Reductions of predators and overabundant prey, translocations, and creating safe havens have been applied in a design covering >90,000 km² Combinations of treatments that increased multiple vital rates produced the highest population growth. Moreover, the degree of ecosystem alteration did not influence this pattern. By coordinating recovery involving scientists, governments, and First Nations, treatments were applied across vast scales to benefit this iconic species.

The Importance of Reproduction for the Conservation of Slow-Growing Animal Populations

Both survival and reproduction are important fitness components, and thus critical to the viability of wildlife populations. Preventing one death (survival) or contributing one newborn (reproduction), has arguably the same effect on population dynamics-in each instance the population grows or is maintained by one additional member. However, for the conservation of slow-growing animal populations, the importance of reproduction is sometimes overlooked when evaluating wildlife management options. This has to do with the use of demographic sensitivity analyses, which quantify the relative contribution of vital rates to population growth. For slow-growing populations, the results of such analyses typically show that growth rates are more sensitive to changes in survival than to equal proportional changes in reproduction. Consequently, for slow-growing taxa, survival has been labelled a better fitness surrogate than reproduction. However, such a generalization, derived from conventional sensitivity analyses, is based on flawed approaches, such as omitting appropriate scaling of vital rates, and sometimes misinterpretations. In this chapter, I make the case that for the conservation of slow-growing species the role of reproduction is considerably greater than conventional sensitivity analyses would suggest. This is illustrated by case studies on wildlife populations that underscore the importance of reproduction for the conservation of slow-growing birds, ungulates, carnivores, and cetaceans.

Population Dynamics Model to Inform Harvest Management of a Small Elk Herd in Central New Mexico

Crop depredation by wildlife is a frequent concern for natural resource managers and mitigation of this issue is often an important task for wildlife agencies. Elk Cervus elaphus and other ungulate species have depredated corn Zea mays at Bosque del Apache National Wildlife Refuge, New Mexico, USA, interfering with the ability of the Refuge to provide sufficient supplemental nutrition to overwintering sandhill cranes Antigone canadensis and geese (Anatidae). We estimated annual adult survival and calf recruitment rates of elk from 2011 to 2013 at Bosque del Apache National Wildlife Refuge. Natural adult survival (excludes human-related mortalities) was high (mean = 98.3%; 95% CI = 95.0–100.0%). Calf recruitment was lower than in some populations, and ranged from 13.0 to 36.7 calves : 100 cows at time of recruitment (March and April) with a mean of 21.9 (SD =12.9). Using this information, we constructed a harvest management model to determine annual harvest quotas required to stabilize the growth of the elk herd on the Refuge. The female segment of the herd is growing at an annual rate of 9.0% (95% CI = −1.1–24.1%). To stabilize the growth rate of the female elk population, 8.0% (95% CI = −1.1–19.4%) of the cows would need to be harvested annually. We estimated an adult elk abundance of 40.0 (SE = 4.57; 95% CI = 33.8–52.6) in 2012 and 61.1 (SE = 7.21; 95% CI = 49.9–78.8) in 2013. Our harvest management model provides Refuge staff, who ultimately intend to improve corn yield, with valuable information needed to stabilize the elk herd. Further, our approach outlines a simple, easily implemented modeling technique that can be used for the management of other ungulate herds.

Nowhere to hide: Effects of linear features on predator-prey dynamics in a large mammal system

Rapid landscape alteration associated with human activity is currently challenging the evolved dynamical stability of many predator-prey systems by forcing species to behaviourally respond to novel environmental stimuli. In many forested systems, linear features (LFs) such as roads, pipelines and resource exploration lines (i.e. seismic lines) are a ubiquitous form of landscape alteration that have been implicated in altering predator-prey dynamics. One hypothesized effect is that LFs facilitate predator movement into and within prey refugia, thereby increasing predator-prey spatial overlap. We evaluated this hypothesis in a large mammal system, focusing on the interactions between boreal woodland caribou (Rangifer tarandus caribou) and their two main predators, wolves (Canis lupus) and black bears (Ursus americanus), during the calving season of caribou. In this system, LFs extend into and occur within peatlands (i.e. bogs and nutrient-poor fens), a habitat type highly used by caribou due to its refugia effects. Using resource selection analyses, we found that LFs increased predator selection of peatlands. Female caribou appeared to respond by avoiding LFs and areas with high LF density. However, in our study area, most caribou cannot completely avoid exposure to LFs and variation in female response had demographic effects. In particular, increasing proportional use of LFs by females negatively impacted survival of their neonate calves. Collectively, these results demonstrate how LFs can reduce the efficacy of prey refugia. Mitigating such effects will require limiting or restoring LFs within prey refugia, although the effectiveness of mitigation efforts will depend upon spatial scale, which in turn will be influenced by the life-history traits of predator and prey.

Experimental moose reduction lowers wolf density and stops decline of endangered caribou

The expansion of moose into southern British Columbia caused the decline and extirpation of woodland caribou due to their shared predators, a process commonly referred to as apparent competition. Using an adaptive management experiment, we tested the hypothesis that reducing moose to historic levels would reduce apparent competition and therefor recover caribou populations. Nested within this broad hypothesis were three specific hypotheses: (1) sport hunting could be used to substantially reduce moose numbers to an ecological target; (2) wolves in this ecosystem were primarily limited by moose abundance; and (3) caribou were limited by wolf predation. These hypotheses were evaluated with a before-after control-impact (BACI) design that included response metrics such as population trends and vital rates of caribou, moose, and wolves. Three caribou subpopulations were subject to the moose reduction treatment and two were in a reference area where moose were not reduced. When the moose harvest was increased, the moose population declined substantially in the treatment area (by 70%) but not the reference area, suggesting that the policy had the desired effect and was not caused by a broader climatic process. Wolf numbers subsequently declined in the treatment area, with wolf dispersal rates 2.5× greater, meaning that dispersal was the likely mechanism behind the wolf numerical response, though reduced recruitment and starvation was also documented in the treatment area. Caribou adult survival increased from 0.78 to 0.88 in the treatment area, but declined in the reference. Caribou recruitment was unaffected by the treatment. The largest caribou subpopulation stabilized in the treatment area, but declined in the reference area. The observed population stability is comparable to other studies that used intensive wolf control, but is insufficient to achieve recovery, suggesting that multiple limiting factors and corresponding management tools must be addressed simultaneously to achieve population growth.

Assessing the importance of demographic parameters for population dynamics using Bayesian integrated population modeling

To successfully respond to changing habitat, climate or harvest, managers need to identify the most effective strategies to reverse population trends of declining species and/or manage harvest of game species. A classic approach in conservation biology for the last two decades has been the use of matrix population models to determine the most important vital rates affecting population growth rate (λ), that is, sensitivity. Ecologists quickly realized the critical role of environmental variability in vital rates affecting λ by developing approaches such as life-stage simulation analysis (LSA) that account for both sensitivity and variability of a vital rate. These LSA methods used matrix-population modeling and Monte Carlo simulation methods, but faced challenges in integrating data from different sources, disentangling process and sampling variation, and in their flexibility. Here, we developed a Bayesian integrated population model (IPM) for two populations of a large herbivore, elk (Cervus canadensis) in Montana, USA. We then extended the IPM to evaluate sensitivity in a Bayesian framework. We integrated known-fate survival data from radio-marked adults and juveniles, fecundity data, and population counts in a hierarchical population model that explicitly accounted for process and sampling variance. Next, we tested the prevailing paradigm in large herbivore population ecology that juvenile survival of neonates <90 d old drives λ using our Bayesian LSA approach. In contrast to the prevailing paradigm in large herbivore ecology, we found that adult female survival explained more of the variation in λ than elk calf survival, and that summer and winter elk calf survival periods were nearly equivalent in importance for λ. Our Bayesian IPM improved precision of our vital rate estimates and highlighted discrepancies between count and vital rate data that could refine population monitoring, demonstrating that combining sensitivity analysis with population modeling in a Bayesian framework can provide multiple advantages. Our Bayesian LSA framework will provide a useful approach to addressing conservation challenges across a variety of species and data types.

Seeking sanctuary: the neonatal calving period among central mountain woodland caribou (Rangifer tarandus caribou)

Woodland caribou (Rangifer tarandus caribou (Gmelin, 1788)) populations have declined throughout their range. With the goal of better understanding habitat selection and fidelity during the neonatal calving period (0–4 weeks), we applied a noninvasive method that estimates calving events and subsequent survival based on changes in movement rates among GPS-collared female caribou. We examined a long-term GPS-collar data set (1998–2014) collected from 81 adult female caribou in two central mountain herds in Alberta and British Columbia, Canada. Although we were unable to validate our results with aerial surveys and pregnancy tests, our estimates of parturition rates, survival rates, calving dates, and habitat selection were consistent with previous studies. We identified 83 calving sites. Female caribou selected calving sites and postparturition habitat on high-elevation ridgetops with gradual slopes and avoided anthropogenic linear features. Female caribou displayed low fidelity to interannual calving ranges with a mean distance of 8.7 km between calving ranges. Fidelity was lower in areas with high seismic-line density. Conservation of high-elevation habitat with limited anthropogenic disturbance is likely to provide the greatest benefit to central mountain caribou during the neonatal calving period, and represents a potential management strategy for population recovery efforts.

Testing predator-prey theory using broad-scale manipulations and independent validation

A robust test of ecological theory is to gauge the predictive accuracy of general relationships parameterized from multiple systems but applied to a new area. To address this goal, we used an ecosystem-level experiment to test predator-prey theory by manipulating prey abundance to determine whether predation was density dependent, density independent, compensatory or depensatory (inversely density dependent) on prey populations. Understanding the nature of predation is of primary importance in community ecology because it establishes whether predation has little effect on prey abundance (compensatory), whether it promotes coexistence (density dependent) and reduces the equilibrium of prey (density independent) or whether it can be destabilizing (depensatory). We used theoretical predictions consisting of functional and numerical equations parameterized independently from meta-analyses on wolves (Canis lupus) and moose (Alces alces), but applied to our specific wolf-moose system. Predictions were tested by experimentally reducing moose abundance across 6500 km(2) as a novel way of evaluating the nature of predation. Depensatory predation of wolves on moose was the best explanation of the population dynamic - a mechanism that has been hypothesized to occur but has rarely been evaluated. Adding locally obtained kill rates and numerical estimates to the independent data provided no benefit to model predictions, suggesting that the theory was robust to local variation. These findings have critical implications for any organism that is preyed upon but that also has, or will be, subject to increased human exploitation or perturbations from environmental change. If depensatory predation is not accounted for in harvest models, predicted yields will be excessive and lead to further population decline.

Reply to the comment by Harron on “Widespread declines in woodland caribou (Rangifer tarandus caribou) continue in Alberta”

Estimation of demographic trends from vital rates provides a powerful means to estimate population trends in cryptic or difficult to study species such as woodland caribou (Rangifer tarandus caribou (Gmelin, 1788)). Using such methods, Hervieux et al. (2013; Can. J. Zool. 91(12): 872–882) recently showed 11 of 14 woodland caribou populations in Alberta were declining at ∼8%/year following up to 18 years of monitoring. Harron (2015; Can. J. Zool. 93(2): 149–150) critiques our original study, claiming that negative biases in our demographic monitoring exaggerate our conclusions of widespread caribou declines. Here, we systematically review each of Harron’s claims of bias, rejecting each of his claims upon careful review of the mechanisms by which his purported claims would manifest in our population trend estimation. Therefore, we conclude that Harron’s scientific critique was superficial and misleading. Delays in conservation actions raised by Harron’s critique risk diminishing opportunities to conserve and recover this federally and provincially protected species.

Managing wolves (Canis lupus) to recover threatened woodland caribou (Rangifer tarandus caribou) in Alberta

Across Canada, woodland caribou (Rangifer tarandus caribou (Gmelin, 1788)) populations are declining because of human-induced changes to food webs that are resulting in apparent competition-induced increases in predator-caused caribou mortality. We tested the hypothesis that wolf (Canis lupus L., 1758) population reduction could reverse declines in a woodland caribou population following a BACI (before-after-control-impact) design conducted over a 12-year period in west-central Alberta, Canada. We monitored annual survival for 172 adult female caribou and calf recruitment from 2000 through 2012 and conducted a provincial government delivered wolf population reduction program annually during the winters of 2005–2006 to 2012 (inclusive) in an area centered on the Little Smoky range. Wolf removal translated to a 4.6% increase in mean population growth rate of the Little Smoky population mostly through improvements in calf recruitment. In contrast, the Red Rock Prairie Creek control population exhibited a 4.7% decline. Although the wolf population reduction program appeared to stabilize the Little Smoky population, it did not lead to population increase, however, with λ remaining approximately equal to 1. Therefore, we recommend, if required, predation management be combined with effective habitat conservation and long-term planning to effect the recovery of species, such as woodland caribou, which are declining as a result of habitat-mediated apparent competition.

Widespread declines in woodland caribou (Rangifer tarandus caribou) continue in Alberta

Nowhere across Canada is the continued persistence of woodland caribou (Rangifer tarandus caribou (Gmelin, 1788)) more uncertain than in Alberta due to widespread industrial development. A recent Government of Canada critical habitat review determined that habitat conditions within all Alberta boreal ecotype caribou ranges are unlikely or very unlikely to allow for self-sustaining caribou populations. This habitat-based assessment was based only indirectly on empirical population trends. Here, we estimated empirical population trend and growth rate (λ) for 13 of Alberta’s 16 remaining woodland caribou populations (plus one adjacent population from Saskatchewan) from 1994 to 2012 using demographic monitoring of adult female survival and calf recruitment. We captured and radio-collared a total of 1337 adult female caribou in 14 populations and estimated the mean annual adult female survival across all populations as 0.851. We conducted 158 late-winter calf recruitment surveys across the 14 populations classifying 20 872 caribou and estimated mean recruitment of 0.154 calves/cow (i.e., 0.077 female calves/cow). We then combined annual estimates of adult female survival and female calf recruitment within each population in a simple age-structured population model to estimate population trend. Annual population growth rate across caribou populations averaged 0.918, and was significantly declining, λ < 1.0, for 10 of the 14 caribou populations. Our results confirm that woodland caribou are declining rapidly (with a realized decline of approximately 50% every 8 years) across Alberta, and support conclusions from previous habitat-based assessments that the population viability of caribou is currently compromised, supporting development and implementation of recovery actions to reverse trends.

Linking habitat selection and predation risk to spatial variation in survival

1. A central assumption underlying the study of habitat selection is that selected habitats confer enhanced fitness. Unfortunately, this assumption is rarely tested, and in some systems, gradients of predation risk may more accurately characterize spatial variation in vital rates than gradients described by habitat selection studies.

2. Here, we separately measured spatial patterns of both resource selection and predation risk and tested their relationships with a key demographic trait, adult female survival, for a threatened ungulate, woodland caribou (Rangifer tarandus caribou Gmelin). We also evaluated whether exposure to gradients in both predation risk and resource selection value was manifested temporally through instantaneous or seasonal effects on survival outcomes.

3. We used Cox proportional hazards spatial survival modelling to assess the relative support for 5 selection- and risk-based definitions of habitat quality, as quantified by woodland caribou adult female survival. These hypotheses included scenarios in which selection ideally mirrored survival, risk entirely drove survival, non-ideal selection correlated with survival but with additive risk effects, an ecological trap with maladaptive selection and a non-spatial effect of annual variation in weather.

4. Indeed, we found positive relationships between the predicted values of a resource selection function (RSF) and survival, yet subsequently incorporating an additional negative effect of predation risk greatly improved models further. This revealed a positive, but non-ideal relationship between selection and survival. Gradients in these covariates were also shown to affect individual survival probability at multiple temporal scales. Exposure to increased predation risk had a relatively instantaneous effect on survival outcomes, whereas variation in habitat suitability predicted by an RSF had both instantaneous and longer-term seasonal effects on survival.

5. Predation risk was an additive source of hazard beyond that detected through selection alone, and woodland caribou selection thus was shown to be non-ideal. Furthermore, by combining spatial adult female survival models with herd-specific estimates of recruitment in matrix population models, we estimated a spatially explicit landscape of population growth predictions for this endangered species.

Inferring parturition and neonate survival from movement patterns of female ungulates: a case study using woodland caribou

Analyses of animal movement data have primarily focused on understanding patterns of space use and the behavioural processes driving them. Here, we analyzed animal movement data to infer components of individual fitness, specifically parturition and neonate survival. We predicted that parturition and neonate loss events could be identified by sudden and marked changes in female movement patterns. Using GPS radio-telemetry data from female woodland caribou (Rangifer tarandus caribou), we developed and tested two novel movement-based methods for inferring parturition and neonate survival. The first method estimated movement thresholds indicative of parturition and neonate loss from population-level data then applied these thresholds in a moving-window analysis on individual time-series data. The second method used an individual-based approach that discriminated among three a priori models representing the movement patterns of non-parturient females, females with surviving offspring, and females losing offspring. The models assumed that step lengths (the distance between successive GPS locations) were exponentially distributed and that abrupt changes in the scale parameter of the exponential distribution were indicative of parturition and offspring loss. Both methods predicted parturition with near certainty (>97% accuracy) and produced appropriate predictions of parturition dates. Prediction of neonate survival was affected by data quality for both methods; however, when using high quality data (i.e., with few missing GPS locations), the individual-based method performed better, predicting neonate survival status with an accuracy rate of 87%. Understanding ungulate population dynamics often requires estimates of parturition and neonate survival rates. With GPS radio-collars increasingly being used in research and management of ungulates, our movement-based methods represent a viable approach for estimating rates of both parameters.