Spatial synchrony of recruitment in mountain-dwelling woodland caribou

Spatial synchrony in population dynamics is a ubiquitous feature across a range of taxa. Understanding factors influencing this synchrony may shed light on important drivers of population dynamics. Three mechanisms influence the degree of spatial synchrony between populations: dispersal, shared predators, and spatial environmental covariance (the Moran effect). We assessed demographic spatial synchrony in recruitment (calf:cow ratio) of 10 northern mountain caribou herds in the Yukon Territory, Canada (1982–2008). Shared predators and dispersal were ruled out as causal mechanisms of spatial recruitment synchrony in these herds and therefore any spatial synchrony should be due to the Moran effect. We also assessed the degree of spatial synchrony in April snow depth to represent environmental variability. The regional average spatial synchrony in detrended residuals of April snow depth was 0.46 (95% CI 0.37 to 0.55). Spatial synchrony in caribou recruitment was weak at 0.13 (95% CI −0.06 to 0.32). The spatial scale of synchrony in April snow depth and caribou recruitment was 330.2 km (95% CI 236.3 to 370.0 km) and 170.0 km (95% CI 69.5 to 282.8 km), respectively. We also investigated how the similarity in terrain features between herds influenced the degree of spatial synchrony using exponential decay models. Only the difference in elevation variability between herds during calving was supported by the data. Herds with more similar elevation variability may track snowmelt ablation patterns in a more similar fashion, which would subsequently result in more synchronized predation rates on calves and/or nutritional effects impacting juvenile survival. Interspecific interactions with predators and alternate prey may also influence spatial synchrony of recruitment in these herds.     

Woodland Caribou Extirpation and Anthropogenic Landscape Disturbance in Ontario

Abstract: The decline of woodland caribou (Rangifer tarandus caribou) has been attributed to anthropogenic landscape disturbances, but critical distance thresholds and time lags between disturbance and extirpation are unknown. Using a database of caribou presence and extirpation for northern Ontario, Canada, geo-coded to 10 times 10-km cells, we constructed logistic regression models to predict caribou extirpation based on distance to the nearest of each of 9 disturbance types: forest cutovers, fires, roads, utility corridors, mines, pits and quarries, lakes, trails, and rail lines. We used Akaike's Information Criterion to select parsimonious models and Receiver-Operating Characteristic curves to derive optimal thresholds. To deal with the effects of spatial autocorrelation on estimates of model significance, we used subsampling and restricted randomizations. Forest cutovers were the best predictor of caribou occupancy, with a tolerance threshold of 13 km to nearest cutover and a time lag of 2 decades between disturbance by cutting and caribou extirpation. Management of woodland caribou should incorporate buffers around habitat and requires long-term monitoring of range occupancy.     

Impacts of spruce budworm defoliation on the habitat of woodland caribou,moose, and their main predators

Forest logging has contributed to the decline of several woodland caribou populations by causing the fragmentation of mature coniferous stands. Such habitat alterations could be worsened by spruce budworm (SBW) outbreaks. Using 6201 vegetation plots from provincial inventories conducted after the last SBW outbreak (1968–1992) in boreal forests of Québec (Canada), we investigated the influence of SBW‐caused tree defoliation and mortality on understory vegetation layers relevant to woodland caribou and its main predators. We found a positive association between severe outbreaks and the cover of most groups of understory plant species, especially in stands that were dominated by balsam fir before the outbreak, where a high canopy openness particularly benefited relatively fast‐growing deciduous plants. Such increases in early successional vegetation could provide high‐quality forage for moose, which is likely to promote higher wolf densities and increase predation pressure on caribou. SBW outbreaks may thus negatively affect woodland caribou by increasing predation risk, the main factor limiting caribou populations in managed forests. For the near future, we recommend updating the criteria used to define critical caribou habitat to consider the potential impacts of spruce budworm defoliation.     

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.     

Boreal Caribou Can Coexist with Natural but Not Industrial Disturbances

For species at risk, it is important that demographic models be consistent with our most recent knowledge because alternate model versions can have differing predictions for wildlife and natural resource management. To establish and maintain this consistency, we can compare predicted model values to current or past observations and demographic knowledge. When novel predictor information becomes available, testing for consistency between modeled and observed values ensures the best models are used for robust, evidence-based, wildlife management. We combine novel information on the extent of historical disturbance regimes (industrial and fire) to an existing demographic model and predict historical and projected demographics of woodland caribou (Rangifer tarandus caribou). Exploring 6 simulation experiments across 5 populations in Alberta, Canada, we identify the relative importance of industrial disturbance, fire, and population density to observed population size and growth rate. We confirm the onset of significant declines across all 5 populations began approximately 30 years ago, demonstrate these declines have been consistent, and conclude they are more likely due to industrial disturbance from the oil and gas sector within contemporary population ranges than historical fire regimes. These findings reinforce recent research on the cause of woodland caribou declines. Testing for consistency between observations and models prescribed for species recovery is paramount for assessing the cause of declines, projecting population trends, and refining recovery strategies for effective wildlife management. We provide a novel simulation method for conducting these tests. © 2020 The Authors. The Journal of Wildlife Management published by Wiley Periodicals LLC on behalf of The Wildlife Society.     

Validation of a Demographic Model for Woodland Caribou

Abstract: Wildlife population models are potentially valuable for conservation planning. Validation is necessary to ensure that models are sufficiently robust for predicting management outcomes consistent with conservation objectives. Sorensen et al. (2008) produced a model of woodland caribou (Rangifer tarandus) population growth rate that was recently modified and used as a predictive tool at several scales. We computed confidence intervals and evaluated the performance of this model using novel data. Confidence intervals were wide, and results suggested that the model may have a positive bias, resulting in over-estimation of population growth rates, as well as low predictive power. Wide confidence intervals mean that current understanding of factors governing woodland caribou herd dynamics is not sufficient for wildlife managers to make reliable projections of responses to management.     

Genetic variation in domestic reindeer and wild caribou in Alaska

Reindeer ( Rangifer tarandus tarandus ) were introduced into Alaska 100 years ago and have been maintained as semidomestic livestock. They have had contact with wild caribou ( R. t. granti ) herds, including deliberate crossbreeding and mixing in the wild. Reindeer have considerable potential as a domestic animal for meat or velvet antler production, and wild caribou are important to subsistence and sport hunters. Our objective was to quantify the genetic relationships of reindeer and caribou in Alaska. We identified allelic variation among five herds of wild caribou and three herds of reindeer with DNA sequencing and restriction enzymes for three loci: a DQA locus of the major histocompatibility complex ( Rata-DQA1 ), K-casein and the D-loop of mitochondrial DNA. These loci are of interest because of their potential influence on domestic animal performance and the fitness of wild populations. There is considerable genetic variation in reindeer and caribou for all three loci, including five, three and six alleles for DQA , K-casein and D-loop respectively. Most alleles occur in both reindeer and caribou, which may be the result of recent common ancestry or genetic introgression in either direction. However, allele frequencies differ considerably between reindeer and caribou, which suggests that gene flow has been limited.     

Global declines of caribou and reindeer

Caribou and reindeer herds are declining across their circumpolar range, coincident with increasing arctic temperatures and precipitation, and anthropogenic landscape change. Here, we examine the mechanisms by which climate warming and anthropogenic landscape change influence caribou and reindeer population dynamics, namely changes in phenology, spatiotemporal changes in species overlap, and increased frequency of extreme weather events, and demonstrate that many caribou and reindeer herds show demographic signals consistent with these changes. While many caribou and reindeer populations historically fluctuated, the current, synchronous population declines emphasize the species' vulnerability to global change. Loss of caribou and reindeer will have significant, negative socioeconomic consequences for northern indigenous cultures.     

Impacts of climate change on the seasonal distribution of migratory caribou

Arctic ecosystems are especially vulnerable to global climate change as temperature and precipitation regimes are altered. An ecologically and socially highly important northern terrestrial species that may be impacted by climate change is the caribou, Rangifer tarandus . We predicted the current and potential future occurrence of two migratory herds of caribou [Rivière George herd (RG) and Rivière-aux-Feuilles (RAF) herd] under a Canadian General Circulation Model climate change scenario, across all seasons in the Québec–Labrador peninsula, using climatic and habitat predictor variables. Argos satellite-tracking collars have been deployed on 213 caribou between 1988 and 2003 with locations recorded every 4–5 days. In addition, we assembled a database of climate (temperature, precipitation, snowfall, timing and length of growing season) and habitat data obtained from the SPOT VEGETATION satellite sensor. Logistic regression models indicated that both climatic and physical habitat variables were significant predictors of current migratory caribou occurrence. Migratory caribou appeared to prefer regions with higher snowfall and lichen availability in the fall and winter. In the summer, caribou preferred cooler areas likely corresponding to a lower prevalence of insects, and they avoided disturbed and recently burnt areas. Climate change projections using climate data predicted an increased range for the RAF herd and decreased range for the RG herd during 2040–2069, limiting the herds to northeastern regions of the Québec–Labrador peninsula. Direct and indirect consequences of climate change on these migratory caribou herds possibly include alteration in habitat use, migration patterns, foraging behaviour, and demography, in addition to social and economic stress to arctic and subarctic native human populations.     

Fire,grazing history,lichen abundance,and winter distribution of caribou in Alaska's taiga

In the early 1990s the Nelchina Caribou (Rangifer tarandus) Herd (NCH) began a dramatic shift to its current winter range, migrating at least an additional 100 km beyond its historic range. We evaluated the impacts of fire and grazing history on lichen abundance and subsequent use and distribution by the NCH. Historic (prior to 1990) and current (2002) winter ranges of the NCH had similar vascular vegetation, lichen cover (P = 0.491), and fire histories (P = 0.535), but the former range had significantly less forage lichen biomass as a result of grazing by caribou. Biomass of forage lichens was twice as great overall (P = 0.031) and 4 times greater in caribou selected sites on the current range than in the historic range, greatly increasing availability to caribou. Caribou on the current range selected for stands with >20% lichen cover (P < 0.001), greater than 1,250 kg/ha (P < 0.001) forage lichen biomass and stands older than 80 yr postfire (P < 0.001). After fires, forage lichen cover and biomass seldom recovered sufficiently to attract caribou grazing until after ≥60 yr, and, as a group, primary forage lichen species did not reach maximum abundance until 180 yr postfire. Recovery following overgrazing can occur much more quickly because lichen cover, albeit mostly fragments, and organic substrates remain present. Our results provide benchmarks for wildlife managers assessing condition of caribou winter range and predicting effects of fires on lichen abundance and caribou distribution. Of our measurements of cover and biomass by species, densities and heights of trees, elevation, slope and aspect, only percentage cover by Cladonia amaurocraea, Cladina rangiferina, Flavocetraria cuculata, and lowbush cranberry (Vaccinium vitis-idaea) were necessary for predicting caribou use of winter range. © 2011 The Wildlife Society     

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.     

Spatial differences in genetic diversity and northward migration suggest genetic erosion along the boreal caribou southern range limit and continued range retraction

With increasing human activities and associated landscape changes, distributions of terrestrial mammals become fragmented. These changes in distribution are often associated with reduced population sizes and loss of genetic connectivity and diversity (i.e., genetic erosion) which may further diminish a species' ability to respond to changing environmental conditions and lead to local population extinctions. We studied threatened boreal caribou (Rangifer tarandus caribou) populations across their distribution in Ontario/Manitoba (Canada) to assess changes in genetic diversity and connectivity in areas of high and low anthropogenic activity. Using data from >1,000 caribou and nine microsatellite loci, we assessed population genetic structure, genetic diversity, and recent migration rates using a combination of network and population genetic analyses. We used Bayesian clustering analyses to identify population genetic structure and explored spatial and temporal variation in those patterns by assembling networks based on R_ST and F_ST as historical and contemporary genetic edge distances, respectively. The Bayesian clustering analyses identified broad‐scale patterns of genetic structure and closely aligned with the R_ST network. The F_ST network revealed substantial contemporary genetic differentiation, particularly in areas presenting contemporary anthropogenic disturbances and habitat fragmentation. In general, relatively lower genetic diversity and greater genetic differentiation were detected along the southern range limit, differing from areas in the northern parts of the distribution. Moreover, estimation of migration rates suggested a northward movement of animals away from the southern range limit. The patterns of genetic erosion revealed in our study suggest ongoing range retraction of boreal caribou in central Canada.     

Estimating ungulate recruitment and growth rates using age ratios

Trends in population growth can be monitored with data for key vital rates without knowledge of abundance. Although adult female survival has the highest elasticity for ungulate population dynamics, the more variable recruitment rates are commonly monitored to track local variation in growth rates. Specifically, recruitment is often measured using late winter young:adult age ratios, though these age ratios are difficult to reliably interpret given the contribution of multiple vital rates to annual ratios. We show that the supplementation of age ratio data with concurrent radio-telemetry monitoring of adult female survival allows both retrospective estimation of empirical population growth rates and the decomposition of recruitment-specific vital rates. We demonstrate the estimation of recruitment and population growth rates for 1 woodland caribou population using these methods, including elasticity and life-stage simulation analysis of the relative contribution of adult female survival and recruitment rates to variation in population growth. We show, for this woodland caribou population, that adult female survival and recruitment rates were nearly equivalent drivers of population growth. We recommend the concurrent monitoring of adult female survival to reliably interpret age ratios when managing caribou and other ungulates. © 2011 The Wildlife Society.     

Population size and major valleys explain microsatellite variation better than taxonomic units for caribou in western Canada

Identifying conservation units below the species level is becoming increasingly important, particularly when limited resources necessitate prioritization for conservation among such units. This problem is exemplified with caribou, a mammal with a circum-Arctic distribution that is exposed to a broad spectrum of ecological conditions, but is also declining in many parts of its range. We used microsatellite markers to evaluate the suitability of existing intra-specific taxonomic designations to act as population units for conservation and contrasted this with landscape features that were independent of taxonomy. We also quantified the relationship between genetic differentiation and subpopulation size, a factor that has been under-represented in landscape genetic research. Our data set included three subspecies and three ecotypes of caribou that varied in population size by five orders of magnitude. Our results indicated that genetic structure did not correspond to existing taxonomic designation, particularly at the level of ecotype. Instead, we found that major valleys and population size were the strongest factors associated with substructure. There was a negative exponential relationship between population size and F(ST) between pairs of adjacent subpopulations, suggesting that genetic drift was the mechanism causing the structure among the smallest subpopulations. A genetic assignment test revealed that movement among subpopulations was a fraction of the level needed to stabilize smaller subpopulations, indicating little chance for demographic rescue. Such results may be broadly applicable to landscape genetic studies, because population size and corresponding rates of drift have the potential to confound interpretations of landscape effects on population structure.     

A data-driven demographic model to explore the decline of the Bathurst caribou herd

The Bathurst herd of barren-ground caribou (Rangifer tarandus groenlandicus) in the Canadian central arctic declined from an estimated 203,800 to 16,400 breeding females from 1986 to 2009, with the most rapid decline from 2006 to 2009. A key research and management question was whether the decline was mainly due to decreases in productivity alone or also due to reduced adult female survival. Investigating causes of the decline was hampered by a lack of direct estimates of caribou demographic parameters. We developed a demographic model that could be objectively fitted to field data to explore the mechanisms for the Bathurst decline, with a focus on the recent accelerated decline from 2006 to 2009. Our modeling indicated that the decline was driven by increasing negative trends in adult female and calf survival rates and possibly reduced fecundity The effect of a constant hunter harvest on the declining herd was one potential cause for the recent accelerated decline in adult survival. The demographic model detected negative trends in adult female survival that were not detected using standalone analyses of collar-based survival data. The model allowed rigorous interpretation of trends in productivity by controlling for the simultaneous influence of trends in adult, calf, and yearling survival and adult fecundity on field-based calf–cow ratios. Stochastic simulations suggested that large increases in adult survival and productivity would be needed for the herd to recover. Our methods enable objective modeling of caribou demography that can assist in caribou management based upon all sources of available data. © 2011 The Wildlife Society.