Post-release survival of translocated fishers: implications for translocation success

As a vital tool for the conservation of species at risk, translocations are also opportunities to identify factors that influence translocation success. We evaluated factors associated with post-release survival of 90 radio-tracked fishers (Pekania pennanti) translocated from central British Columbia, Canada, to the Olympic Peninsula of Washington, USA, from 2008 to 2011. We hypothesized that the survival of translocated fishers would be affected by the same factors that influence the survival of resident, native fishers (i.e., sex, age, season, body condition), and additional factors that were associated with the translocation process (e.g., duration of captivity, release date, yr of release). Fisher survival was most strongly influenced by translocation year (i.e., release-yr cohort), season, sex, and age class of fisher; whereas duration of captivity, standardized body mass, release date, and number of intact canines did not influence survival. Survival was lowest for fishers released in cohort 2 in 2009 and during the breeding season (Mar–Jun), and was greatest for juveniles and males. When combined across release-year cohorts, year 1 survival rates were greatest for juvenile males followed by juvenile females, adult females, and adult males. Sex and age-related differences in survival of translocated fishers were counter to those commonly reported for established fisher populations, where adult females often have the highest survival rates and juveniles the lowest. Predation (40%) and vehicle strikes (20%) were the most common causes of known mortality among the 24 recovered fishers for which cause of death was determined. We speculate that females face higher risks of mortality in translocated populations because their small size makes them more vulnerable to predation and because adult females in resident populations are less likely than males and juveniles to disperse. Our findings support designing translocations that favor releasing a preponderance of female fishers in recognition of their lower survival rates and to ensure adequate breeders are established in the population, and juvenile and young adult fishers to enhance survival of both sexes. Releases conducted over multiple years will minimize the impact of stochastic annual events that may adversely affect survival in any given year. Persistence, widespread distribution, and documented reproduction of fishers within our study area for ≥6 years following the last releases indicate that survival parameters we measured contributed toward successful population establishment over the short term.     

A tale of two populations: vital rates of fishers in British Columbia,Canada

Fishers (Pekania pennanti) are a forest-dependent carnivore of conservation concern in British Columbia, Canada. Ecological, spatial, and genetic evidence suggests that there are 2 distinct populations (Boreal and Columbian) that occur in forests at low to moderate elevations in the boreal and central interior regions of the province. In British Columbia, fishers occur at low densities relative to other parts of their range in North America, are trapped for their fur, and are sensitive to habitat change. Despite these factors, little demographic information exists to assist with management decisions for these populations. We collated and analyzed survival and reproductive data from 100 radio-tagged fishers from 5 independent studies conducted between 1990 and 2012 in British Columbia: 2 in the Boreal population, and 3 in the Columbian population. We also collated litter size data from 1 den box study and a translocation project of fishers from the Columbian population. Annual survival rates were not significantly different between the populations or between males and females; however, adult survival rates were higher than subadults (0.79 and 0.63, respectively). Subadult females had significantly lower survival rates than other sex or age classes. Reproductive rates were significantly different between the 2 populations (denning rate = 0.54 [Columbian], 0.82 [Boreal]; $\overline{x}$ litter size = 1.7 [Columbian], 2.6 [Boreal]). These differences resulted in net reproductive rates in the Columbian population that were less than half of those in the Boreal population (0.92 kits/reproductive season compared to 2.13, respectively). Population growth rates suggest that the Columbian population may have been declining during the studies, whereas the Boreal population may have been increasing (0.96 compared to 1.20). Consequently, we suggest that focused and intensive habitat and population management for fishers are needed in British Columbia to ensure population sustainability, particularly for the Columbian population.     

Occupancy Patterns in a Reintroduced Fisher Population during Reestablishment

Monitoring population performance in the years following species reintroductions is key to assessing population restoration success and evaluating assumptions made in planning species restoration programs. From 2008–2010 we translocated 90 fishers (Pekania pennanti) from British Columbia, Canada, to Washington's Olympic Peninsula, USA, providing the opportunity to evaluate modeling assumptions used to identify the most suitable reintroduction areas in Washington and enhance understanding of fisher habitat associations in the late-successional forest ecosystems in the coastal Pacific Northwest. From 2013–2016, we deployed 788 motion-sensing cameras and hair (DNA)-snaring devices distributed among 263 24-km² primary sampling units across the Olympic Peninsula. Our objectives were to determine whether occupancy patterns of the reestablishing population supported assumptions of the initial habitat assessment models, whether the population had expanded or shifted in distribution since the initial reintroductions, compare physical habitat attributes among land-management designations, and determine whether the founding fishers had successfully reproduced. We predicted that site occupancy by fishers would be associated with landscapes characterized by high proportional coverage of dense forest canopies and medium-sized and large trees, a diversity of stand structural classes, and area near the administrative boundary separating wilderness from more intensively managed forest lands. We detected fishers across designated wilderness, federal lands outside of wilderness, and other land designations in proportion to land availability on the Peninsula. We found negligible support for predictions that occupancy by fishers was associated with percent forest cover, tree-size class, or structural class diversity. Rather, occupancy was strongly associated with lands near the wilderness boundary on both sides. We speculate that the boundary between wilderness and more intensively managed forest lands provided fishers with the most suitable prey in proximity to contiguous expanses of low- to mid-elevation late-successional forests that provided optimal resting, denning, and security values. Occupancy patterns shifted toward the west and south along a precipitation gradient during the study, indicating that population distribution had not yet stabilized 5–8 years following translocation. Genetic results indicated that ≥2 generations of fishers have been produced on the Peninsula. Annual occupancy rates across the Peninsula (0.08–0.24) were lower than in other previously studied and established fisher populations, indicating that not all habitat was fully occupied or that initial estimates of the extent of habitat was overestimated. The strong selection fishers exhibited for wilderness edge and weak selection against extensive forested wilderness areas suggested that habitat managers should strive for maintaining a suitable interspersion of required forest structures and biotic habitat components, such as prey resource availability. © 2019 The Wildlife Society.     

Estimating occupancy to monitor northern goshawk in the central Rocky Mountains

Designing monitoring programs to evaluate trends in low-density wildlife species at regional scales is challenging given difficulties detecting uncommon organisms distributed in potential habitats over large spatial extents. The northern goshawk (Accipiter gentilis) has been petitioned for listing under the Endangered Species Act and the review of the petition indicated a need for information on population trend. To evaluate trends in goshawk populations, the U.S. Forest Service developed the Northern Goshawk Bioregional Monitoring Design to estimate goshawk occupancy over broad spatial extents. We adapted and implemented this design to approximately 30,600 km² of 88,128 km² of National Forest System lands in the Forest Service Rocky Mountain Region, including portions of Colorado, Wyoming, and South Dakota. We developed a stratified random design to monitor goshawk occupancy in sampling units, defined by primary and secondary habitat quality as well as accessibility. To define habitat quality, we examined a time series for 58 previously located nesting territories. Using logistic regression, we found that the dominant conifer species and status of aspen in postfledging zones best characterized high-quality goshawk nesting habitat. We applied model results to stratify 4,445 sampling units based on habitat quality and further stratified sampling units based on accessibility into easy and difficult access categories. We conducted field sampling during the goshawk breeding season in the summer of 2006 to estimate detection probabilities and occupancy rates. Within our sampling frame, we sampled 51 sampling units and estimated goshawk occupancy of 0.329 (95% CI: 0.213–0.445). Occupancy within primary strata (high quality) sampling units was 0.811 (SE = 0.113), whereas occupancy in secondary strata (lower quality) sampling units was 0.124 (SE = 0.067). Future implementation of this monitoring program can achieve 0.8 power to detect 30–40% declines in with 140 sampling units. Our implementation of a stratified sampling design to monitor occupancy of goshawks at a region-wide scale reduced the number of sampling units in each administrative unit and focused our efforts on those areas most likely to have goshawks. © 2011 The Wildlife Society.     

Wolverine Occupancy,Spatial Distribution,and Monitoring Design

In the western United States, wolverines (Gulo gulo) typically occupy high-elevation habitats. Because wolverine populations occur in vast, remote areas across multiple states, biologists have an imperfect understanding of this species' current distribution and population status. The historical extirpation of the wolverine, a subsequent period of recovery, and the lack of a coordinated monitoring program in the western United States to determine their current distribution further complicate understanding of their population status. We sought to define the limits to the current distribution, identify potential gaps in distribution, and provide a baseline dataset for future monitoring and analysis of factors contributing to changes in distribution of wolverines across 4 western states. We used remotely triggered camera stations and hair snares to detect wolverines across randomly selected 15-km × 15-km cells in Idaho, Montana, Washington, and Wyoming, USA, during winters 2016 and 2017. We used spatial occupancy models to examine patterns in wolverine distribution. We also examined the influence of proportion of the cell containing predicted wolverine habitat, human-modified land, and green vegetation, and area of the cluster of contiguous sampling cells. We sampled 183 (28.9%) of 633 cells that comprised a suspected wolverine range in these 4 states and we detected wolverines in 59 (32.2%) of these 183 sampled cells. We estimated that 268 cells (42.3%; 95% CI = 182–347) of the 633 cells were used by wolverines. Proportion of the cell containing modeled wolverine habitat was weakly positively correlated with wolverine occupancy, but no other covariates examined were correlated with wolverine occupancy. Occupancy rates (ψ) were highest in the Northern Continental Divide Ecosystem (ψ range = 0.8–1), intermediate in the Cascades and Central Mountains of Idaho (ψ range = 0.4–0.6), and lower in the Greater Yellowstone Ecosystem (ψ range = 0.1–0.3). We provide baseline data for future surveys of wolverine along with a design and protocol to conduct those surveys. © 2020 The Authors. The Journal of Wildlife Management published by Wiley Periodicals, Inc. on behalf of The Wildlife Society.     

Seasonal Resource Selection of Canada Lynx in Managed Forests of the Northern Rocky Mountains

ABSTRACT We investigated seasonal patterns in resource selection of Canada lynx (Lynx canadensis) in the northern Rockies (western MT, USA) from 1998 to 2002 based on backtracking in winter (577 km; 10 M, 7 F) and radiotelemetry (630 locations; 16 M, 11 F) in summer. During winter, lynx preferentially foraged in mature, multilayer forests with Engelmann spruce (Picea engelmannii) and subalpine fir (Abies lasiocarpa) in the overstory and midstory. Forests used during winter were composed of larger diameter trees with higher horizontal cover, more abundant snowshoe hares (Lepus americanus), and deeper snow compared to random availability; multilayer, spruce-fir forests provided high horizontal cover with tree branching that touched the snow surface. During winter, lynx killed prey at sites with higher horizontal cover than that along foraging paths. Lynx were insensitive to snow depth or penetrability in determining where they killed prey. During summer, lynx broadened their resource use to select younger forests with high horizontal cover, abundant total shrubs, abundant small-diameter trees, and dense saplings, especially spruce-fir saplings. Based on multivariate logistic-regression models, resource selection occurred primarily at a fine spatial scale as was consistent with a sight-hunting predator in dense forests. However, univariate comparisons of patch-level metrics indicated that lynx selected homogenous spruce-fir patches, and avoided recent clear-cuts or other open patches. Given that lynx in Montana exhibit seasonal differences in resource selection, we encourage managers to maintain habitat mosaics. Because winter habitat may be most limiting for lynx, these mosaics should include abundant multistory, mature spruce-fir forests with high horizontal cover that are spatially well-distributed.     

ednaoccupancy: An r package for multiscale occupancy modelling of environmental DNA data

In this article, we describe ednaoccupancy , an r package for fitting Bayesian, multiscale occupancy models. These models are appropriate for occupancy surveys that include three nested levels of sampling: primary sample units within a study area, secondary sample units collected from each primary unit and replicates of each secondary sample unit. This design is commonly used in occupancy surveys of environmental DNA (eDNA). ednaoccupancy allows users to specify and fit multiscale occupancy models with or without covariates, to estimate posterior summaries of occurrence and detection probabilities, and to compare different models using Bayesian model‐selection criteria. We illustrate these features by analysing two published data sets: eDNA surveys of a fungal pathogen of amphibians and eDNA surveys of an endangered fish species.     

Tree recruitment at the treeline across the Continental Divide in the Northern Rocky Mountains,USA: the role of spring snow and autumn climate

ABSTRACT

Background: Topoclimate can influence tree establishment within treeline ecotones. Yet much less is known about how regional topography, such as the Continental Divide, Rocky Mountains, mediates the role of climate in governing treeline dynamics.

Aims: To utilise the Continental Divide to test whether contrasts in growing-season moisture regimes to the west (summer-dry) and east (summer-wet) impact the spatio-temporal patterns of tree establishment and rates of treeline advance in the Northern Rocky Mountains.

Methods: We sampled trees at sites on north- and south-facing slopes, west and east of the Continental Divide. We used dendroecological techniques to reconstruct patterns of tree establishment. Age-structure data were quantitatively compared with climate to evaluate possible mechanistic linkages.

Results: Across all sites, 96% of trees established after 1950. There was a treeline advance (range = 39–140 m) accompanied by increases in tree density. Significantly more trees established during wet springs on both sides of the Divide.

Conclusions: Overall, snow duration in spring and autumn temperatures appear to influence patterns of tree recruitment at the treeline. Continued warming will likely amplify the role of autumn climate in regulating tree establishment throughout treeline ecotones in the Northern Rocky Mountains, particularly west of the Divide where summer-dry conditions persist.     

Wolf population dynamics in the U.S. Northern Rocky Mountains are affected by recruitment and human-caused mortality

Reliable analyses can help wildlife managers make good decisions, which are particularly critical for controversial decisions such as wolf (Canis lupus) harvest. Creel and Rotella (2010) recently predicted substantial population declines in Montana wolf populations due to harvest, in contrast to predictions made by Montana Fish, Wildlife and Parks (MFWP). We replicated their analyses considering only those years in which field monitoring was consistent, and we considered the effect of annual variation in recruitment on wolf population growth. Rather than assuming constant rates, we used model selection methods to evaluate and incorporate models of factors driving recruitment and human-caused mortality rates in wolf populations in the Northern Rocky Mountains. Using data from 27 area-years of intensive wolf monitoring, we show that variation in both recruitment and human-caused mortality affect annual wolf population growth rates and that human-caused mortality rates have increased with the sizes of wolf populations. We document that recruitment rates have decreased over time, and we speculate that rates have decreased with increasing population sizes and/or that the ability of current field resources to document recruitment rates has recently become less successful as the number of wolves in the region has increased. Estimates of positive wolf population growth in Montana from our top models are consistent with field observations and estimates previously made by MFWP for 2008–2010, whereas the predictions for declining wolf populations of Creel and Rotella (2010) are not. Familiarity with limitations of raw data, obtained first-hand or through consultation with scientists who collected the data, helps generate more reliable inferences and conclusions in analyses of publicly available datasets. Additionally, development of efficient monitoring methods for wolves is a pressing need, so that analyses such as ours will be possible in future years when fewer resources will be available for monitoring. © 2011 The Wildlife Society.     

Optimal sampling design and the accuracy of occupancy models

We present general theoretical limits on the possible accuracy (mean squared error or MSE) of occupancy estimates for a large range of occupancy study designs with imperfect detection and confirm our theoretical results via a simulation study. In particular, we show that for a given total survey effort, the best possible MSE is driven by two design-related factors: the fraction of visits made at occupied sites (regardless of whether that occupancy status is known or not) and the number of visits made to each site with unknown occupancy status (ie, sites with no detections). The limits reveal that there is very little room for improvement over optimal implementations of the three existing occupancy design paradigms: standard design (visit S sites K times each), removal design (visit S sites up to K times each, halting visits to each site following a positive detection), and conditional design (visit S sites once, then resurvey sites with a positive detection an additional $K-1$ times). For the small portion of the occupancy-detection parameter space where improvement can be achieved, we introduce a new hybrid survey design with accuracy closer to the theoretical limit, which we illustrate by reanalyzing an existing coyote (Canis latrans) camera trap dataset. Our results provide new clarity and intuition regarding key factors of occupancy study design.     

Marten and fisher responses to fluctuations in prey populations and mast crops in the northern hardwood forest

Many forest tree species produce seed (mast) crops that are consumed by a variety of wildlife species and these pulsed resources may mediate interactions among predator and prey populations. In the northern hardwood forests of New York, we investigated interactions among mast production, prey abundance, and harvests of American martens (Martes americana) and fishers (Martes pennanti) during 1988–2009. Mast production for beech (Fagus grandifolia), sugar maple (Acer saccharum), and mountain ash (Sorbus americana) was synchronous and an alternate-year pattern in production was evident for most of the time series. We documented considerable temporal variation in summer small mammal relative abundance and our numerical response models received substantial support for 5 of the 8 species, indicating lagged responses to autumn mast crops. Trap response of martens to the autumn production of beech mast and mountain ash berries was immediate and numerical responses to the relative abundance of small mammal prey occurred during the preceding summer. The age structure of the marten harvest differed based on the dominant alternate-year pattern of summer prey relative abundance and autumn mast production (χ²₄ = 33.06, P < 0.001). The proportion of juvenile marten in the autumn harvest was 52% and 34% following summers when small mammal relative abundance was high and low, respectively and these differences resulted in a persistent cohort effect that was apparent until age 3.5. Trap response of fishers to the autumn production of beech mast was immediate and numerical responses to the relative abundance of Sciurid prey occurred during the preceding summer. Marten and fisher harvests fluctuated similarly among New York, Maine, and New Brunswick, which may indicate regional synchronization of mast crops and responses of martens and fishers to similar prey dynamics. A better understanding of how food availability influences demographic responses and trapping vulnerability of martens and fishers would aid our ability to manage harvests of these species on a sustained yield basis. © 2011 The Wildlife Society.     

Optimizing occupancy surveys by maximizing detection probability: application to amphibian monitoring in the Mediterranean region

• Setting up effective conservation strategies requires the precise determination of the targeted species’ distribution area and, if possible, its local abundance. However, detection issues make these objectives complex for most vertebrates. The detection probability is usually <1 and is highly dependent on species phenology and other environmental variables. The aim of this study was to define an optimized survey protocol for the Mediterranean amphibian community, that is, to determine the most favorable periods and the most effective sampling techniques for detecting all species present on a site in a minimum number of field sessions and a minimum amount of prospecting effort. We visited 49 ponds located in the Languedoc region of southern France on four occasions between February and June 2011. Amphibians were detected using three methods: nighttime call count, nighttime visual encounter, and daytime netting. The detection nondetection data obtained was then modeled using site‐occupancy models. The detection probability of amphibians sharply differed between species, the survey method used and the date of the survey. These three covariates also interacted. Thus, a minimum of three visits spread over the breeding season, using a combination of all three survey methods, is needed to reach a 95% detection level for all species in the Mediterranean region. Synthesis and applications: detection nondetection surveys combined to site occupancy modeling approach are powerful methods that can be used to estimate the detection probability and to determine the prospecting effort necessary to assert that a species is absent from a site.