Success of Grizzly Bear Population Augmentation in Northwest Montana

Abstract: Augmentation of large carnivore populations can be a valuable management and recovery tool, but success of many programs has not been well documented. The Cabinet—Yaak grizzly bear (Ursus arctos) population was located in northwestern Montana and northern Idaho, USA, and was estimated at 30–40 individuals. The Cabinet Mountains portion of this area may be isolated from the remainder of the zone and was the site of a test of grizzly bear population augmentation. Experimental objectives included evaluating site fidelity, reproduction, and long-term survival of the translocated bears. Four subadult females (2–6 yr old) were translocated from southeastern British Columbia, Canada, from 1990 to 1994. Three of 4 transplanted bears remained in the target area for ≥1 year and satisfied the short-term goal for site fidelity. Recent genetic evidence gathered through hair-snagging efforts has determined that at least one of the original transplanted animals has reproduced, thereby providing evidence of success for the long-term goals of survival and reproduction.     

Evaluation of Bear Rub Surveys to Monitor Grizzly Bear Population Trends

Abstract: Wildlife managers need reliable estimates of population size, trend, and distribution to make informed decisions about how to recover at-risk populations, yet obtaining these estimates is costly and often imprecise. The grizzly bear (Ursus arctos) population in northwestern Montana, USA, has been managed for recovery since being listed under the United States Endangered Species Act in 1975, yet no rigorous data were available to evaluate the program's success. We used encounter data from 379 grizzly bears identified through bear rub surveys to parameterize a series of Pradel model simulations in Program MARK to assess the ability of noninvasive genetic sampling to estimate population growth rates. We evaluated model performance in terms of 1) power to detect gender-specific and population-wide declines in population abundance, 2) precision and relative bias of growth rate estimates, and 3) sampling effort required to achieve 80% power to detect a decline within 10 years. Simulations indicated that ecosystem-wide, annual bear rub surveys would exceed 80% power to detect a 3% annual decline within 6 years. Robust-design models with 2 simulated surveys per year provided precise and unbiased annual estimates of trend, abundance, and apparent survival. Designs incorporating one survey per year require less sampling effort but only yield trend and apparent survival estimates. Our results suggest that systematic, annual bear rub surveys may provide a viable complement or alternative to telemetry-based methods for monitoring trends in grizzly bear populations.     

Grizzly bear responses to restrictions of recreation in Yellowstone National Park

Avoiding humans will be more difficult and energetically costly for animals as outdoor recreation increases and people venture farther into wildland areas that provide high-quality habitat for wildlife. Restricting human access can be an attractive management tool to mitigate effects of human recreation activities on wildlife; however, the efficacy of such measures is rarely assessed. In 1982, Yellowstone National Park identified areas important to grizzly bears (Ursus arctos) to help protect critical grizzly bear habitat and reduce the likelihood of human injuries by bears. Referred to as bear management areas (BMAs), human access is restricted in these areas for 2–8 months each year, with timing and type of restrictions varying by area. We examined 2 datasets to evaluate grizzly bear selection of BMAs and differences of bear density in BMAs and non-BMAs. First, we used 17 years of recent global positioning system telemetry data for grizzly bears to assess their selection of BMAs during periods when human access was allowed, and when access was restricted. We used step-selection functions to test the hypothesis that bears spend time in places that allow them to avoid people and select quality food sources. There was support that grizzly bears differentially select for BMAs regardless of whether human access was restricted at the time, compared with areas outside BMAs, and that selection changed with sex and season. Only males during the summer and hyperphagic seasons changed their selection of BMAs based on whether access restrictions were in place, and overall, male bears preferred unrestricted BMAs (BMAs without restrictions in place). Females preferentially selected BMAs regardless of whether the area had access restrictions in place only during the mating season. Individuals varied widely in their preference for BMAs and access restrictions. Bears likely choose to spend time in BMAs based on available food resources rather than restrictions to human access. Supporting this interpretation, our analyses indicated that a greater proportion of BMA in an area was associated with higher densities of grizzly bear. Thus, restrictions to human access likely help reduce the potential for human–bear interactions, accomplishing one of the original objectives for establishing the BMAs.     

Demography and Genetic Structure of a Recovering Grizzly Bear Population

ABSTRACT Grizzly bears (brown bears; Ursus arctos) are imperiled in the southern extent of their range worldwide. The threatened population in northwestern Montana, USA, has been managed for recovery since 1975; yet, no rigorous data were available to monitor program success. We used data from a large noninvasive genetic sampling effort conducted in 2004 and 33 years of physical captures to assess abundance, distribution, and genetic health of this population. We combined data from our 3 sampling methods (hair trap, bear rub, and physical capture) to construct individual bear encounter histories for use in Huggins—Pledger closed mark—recapture models. Our population estimate, Ň = 765 (95% CI = 715–831) was more than double the existing estimate derived from sightings of females with young. Based on our results, the estimated known, human-caused mortality rate in 2004 was 4.6% (95% CI = 4.2–4.9%), slightly above the 4% considered sustainable; however, the high proportion of female mortalities raises concern. We used location data from telemetry, confirmed sightings, and genetic sampling to estimate occupied habitat. We found that grizzly bears occupied 33,480 km² in the Northern Continental Divide Ecosystem (NCDE) during 1994–2007, including 10,340 km² beyond the Recovery Zone. We used factorial correspondence analysis to identify potential barriers to gene flow within this population. Our results suggested that genetic interchange recently increased in areas with low gene flow in the past; however, we also detected evidence of incipient fragmentation across the major transportation corridor in this ecosystem. Our results suggest that the NCDE population is faring better than previously thought, and they highlight the need for a more rigorous monitoring program.     

Components of Grizzly Bear Habitat Selection: Density,Habitats, Roads,and Mortality Risk

Abstract: We used resource selection functions (RSF) to estimate the relative probability of use for grizzly bears (Ursus arctos) adjacent to the Parsnip River, British Columbia, Canada, 1998-2003. We collected data from 30 radiocollared bears on a rolling plateau where a large portion of the landscape had been modified by human activities, primarily forestry. We also monitored 24 radiocollared bears in mountain areas largely inaccessible to humans. Bears that lived on the plateau existed at less than one-quarter the density of bears in the mountains. Plateau bears ate more high-quality food items, such as meat and berries, leading us to conclude that food limitation was not responsible for the differences in densities. We hypothesized that plateau bears were limited by human-caused mortality associated with roads constructed for forestry activities. Independent estimates of bear population size from DNA-based mark-recapture techniques allowed us to link populations to habitats using RSF models to scale habitat use patterns to population density. To evaluate whether differences in land-cover type, roads, or mortality risk could account for the disparity in density we used the mountain RSF model to predict habitat use and number of bears on the plateau and vice versa. We predicted increases ranging from 34 bears to 96 bears on the plateau when switching model coefficients, excluding land-cover types; when exchanging land-cover coefficients, the model predicted that the plateau population would be 9 bears lower than was observed. Large reductions in the numbers of mountain bears were predicted by habitat-selection models of bears using the plateau landscape. Although RSF models estimated in mountain and plateau landscapes could not predict bear use and abundance in the other areas, contrasts in models between areas provided a useful tool for examining the effects of human activities on grizzly bears.     

Grizzly bear population vital rates and trend in the Northern Continental Divide Ecosystem,Montana

We estimated grizzly bear (Ursus arctos) population vital rates and trend for the Northern Continental Divide Ecosystem (NCDE), Montana, between 2004 and 2009 by following radio-collared females and observing their fate and reproductive performance. Our estimates of dependent cub and yearling survival were 0.612 (95% CI = 0.300–0.818) and 0.682 (95% CI = 0.258–0.898). Our estimates of subadult and adult female survival were 0.852 (95% CI = 0.628–0.951) and 0.952 (95% CI = 0.892–0.980). From visual observations, we estimated a mean litter size of 2.00 cubs/litter. Accounting for cub mortality prior to the first observations of litters in spring, our adjusted mean litter size was 2.27 cubs/litter. We estimated the probabilities of females transitioning from one reproductive state to another between years. Using the stable state probability of 0.322 (95% CI = 0.262–0.382) for females with cub litters, our adjusted fecundity estimate (m_x) was 0.367 (95% CI = 0.273–0.461). Using our derived rates, we estimated that the population grew at a mean annual rate of approximately 3% (λ = 1.0306, 95% CI = 0.928–1.102), and 71.5% of 10,000 Monte Carlo simulations produced estimates of λ > 1.0. Our results indicate an increasing population trend of grizzly bears in the NCDE. Coupled with concurrent studies of population size, we estimate that over 1,000 grizzly bears reside in and adjacent to this recovery area. We suggest that monitoring of population trend and other vital rates using radioed females be continued. © 2011 The Wildlife Society.     

Dietary adjustability of grizzly bears and American black bears in Yellowstone National Park

Grizzly bears (Ursus arctos) and American black bears (U. americanus) are sympatric in much of Yellowstone National Park. Three primary bear foods, cutthroat trout (Oncorhynchus clarki), whitebark pine (Pinus albicaulis) nuts, and elk (Cervus elaphus), have declined in recent years. Because park managers and the public are concerned about the impact created by reductions in these foods, we quantified bear diets to determine how bears living near Yellowstone Lake are adjusting. We estimated diets using: 1) stable isotope and mercury analyses of hair samples collected from captured bears and from hair collection sites established along cutthroat trout spawning streams and 2) visits to recent locations occupied by bears wearing Global Positioning System collars to identify signs of feeding behavior and to collect scats for macroscopic identification of residues. Approximately 45 ± 22% ( ± SD) of the assimilated nitrogen consumed by male grizzly bears, 38 ± 20% by female grizzly bears, and 23 ± 7% by male and female black bears came from animal matter. These assimilated dietary proportions for female grizzly bears were the same as 10 years earlier in the Lake area and 30 years earlier in the Greater Yellowstone Ecosystem. However, the proportion of meat in the assimilated diet of male grizzly bears decreased over both time frames. The estimated biomass of cutthroat trout consumed by grizzly bears and black bears declined 70% and 95%, respectively, in the decade between 1997–2000 and 2007–2009. Grizzly bears killed an elk calf every 4.3 ± 2.7 days and black bears every 8.0 ± 4.0 days during June. Elk accounted for 84% of all ungulates consumed by both bear species. Whitebark pine nuts continue to be a primary food source for both grizzly bears and black bears when abundant, but are replaced by false-truffles (Rhizopogon spp.) in the diets of female grizzly bears and black bears when nut crops are minimal. Thus, both grizzly bears and black bears continue to adjust to changing resources, with larger grizzly bears continuing to occupy a more carnivorous niche than the smaller, more herbivorous black bear. © 2012 The Wildlife Society.     

Hazards Affecting Grizzly Bear Survival in the Greater Yellowstone Ecosystem

Abstract: During the past 2 decades, the grizzly bear (Ursus arctos) population in the Greater Yellowstone Ecosystem (GYE) has increased in numbers and expanded its range. Early efforts to model grizzly bear mortality were principally focused within the United States Fish and Wildlife Service Grizzly Bear Recovery Zone, which currently represents only about 61% of known bear distribution in the GYE. A more recent analysis that explored one spatial covariate that encompassed the entire GYE suggested that grizzly bear survival was highest in Yellowstone National Park, followed by areas in the grizzly bear Recovery Zone outside the park, and lowest outside the Recovery Zone. Although management differences within these areas partially explained differences in grizzly bear survival, these simple spatial covariates did not capture site-specific reasons why bears die at higher rates outside the Recovery Zone. Here, we model annual survival of grizzly bears in the GYE to 1) identify landscape features (i.e., foods, land management policies, or human disturbances factors) that best describe spatial heterogeneity among bear mortalities, 2) spatially depict the differences in grizzly bear survival across the GYE, and 3) demonstrate how our spatially explicit model of survival can be linked with demographic parameters to identify source and sink habitats. We used recent data from radiomarked bears to estimate survival (1983–2003) using the known-fate data type in Program MARK. Our top models suggested that survival of independent (age ≥ 2 yr) grizzly bears was best explained by the level of human development of the landscape within the home ranges of bears. Survival improved as secure habitat and elevation increased but declined as road density, number of homes, and site developments increased. Bears living in areas open to fall ungulate hunting suffered higher rates of mortality than bears living in areas closed to hunting. Our top model strongly supported previous research that identified roads and developed sites as hazards to grizzly bear survival. We also demonstrated that rural homes and ungulate hunting negatively affected survival, both new findings. We illustrate how our survival model, when linked with estimates of reproduction and survival of dependent young, can be used to identify demographically the source and sink habitats in the GYE. Finally, we discuss how this demographic model constitutes one component of a habitat-based framework for grizzly bear conservation. Such a framework can spatially depict the areas of risk in otherwise good habitat, providing a focus for resource management in the GYE.     

Brown Bear Density and Estimated Harvest Rates in Northwestern Alaska

Human-caused mortality in general, and unregulated hunting in particular, have been implicated in reductions in brown bear (Ursus arctos) populations throughout much of their range. In northwestern Alaska, USA, bear densities have not been assessed in 20 years while harvest regulations have been liberalized, raising concerns that broad undetected population declines might occur. We used a modified mark-resight approach to estimate brown bear density during 2005–2018 in 4 subareas throughout the region. We also summarized harvest information for each subarea and used our survey results to estimate harvest rates. We estimated densities for independent bears assuming constant or heterogeneous probabilities of detection and occurrence. We present the results of the constant model for more direct comparison with past work and the heterogeneity model results to provide estimates of density that are less likely to be negatively biased. Using the constant model, we estimated the density of independent bears was 17.0, 49.2, 24.9, and 19.4/1,000 km² on portions of the Seward Peninsula, the lower Noatak River, the upper Noatak River, and Gates of the Arctic National Park and Preserve, respectively. These estimates are broadly similar to those from past work in interior and northwestern Alaska, with the exception of the lower Noatak River subarea where our estimates are the highest reported for a bear population in northern Alaska. We estimated that the harvest rate on the Seward Peninsula was approximately 5.2% or 7.7% on average, depending upon the model used. In the remaining areas, we estimated annual harvest rates were <2.5%, well within sustainability guidelines from past work. Overall, our results suggest that brown bear densities are similar or somewhat higher than in the past in much of northwestern Alaska and that current harvest rates are sustainable in most areas, except perhaps the Seward Peninsula. Ongoing survey work will be useful for further evaluating the assumptions of the modified mark-resight survey approach, assessing population trajectory, and determining the effect of harvest on brown bear populations. © 2021 The Wildlife Society.     

Evaluation of Rules to Distinguish Unique Female Grizzly Bears With Cubs in Yellowstone

ABSTRACT The United States Fish and Wildlife Service uses counts of unduplicated female grizzly bears (Ursus arctos) with cubs-of-the-year to establish limits of sustainable mortality in the Greater Yellowstone Ecosystem, USA. Sightings are clustered into observations of unique bears based on an empirically derived rule set. The method has never been tested or verified. To evaluate the rule set, we used data from radiocollared females obtained during 1975–2004 to simulate populations under varying densities, distributions, and sighting frequencies. We tested individual rules and rule-set performance, using custom software to apply the rule-set and cluster sightings. Results indicated most rules were violated to some degree, and rule-based clustering consistently underestimated the minimum number of females and total population size derived from a nonparametric estimator (Chao2). We conclude that the current rule set returns conservative estimates, but with minor improvements, counts of unduplicated females-with-cubs can serve as a reasonable index of population size useful for establishing annual mortality limits. For the Yellowstone population, the index is more practical and cost-effective than capture-mark-recapture using either DNA hair snagging or aerial surveys with radiomarked bears. The method has useful application in other ecosystems, but we recommend rules used to distinguish unique females be adapted to local conditions and tested.     

Noninvasive Estimation of Black Bear Abundance Incorporating Genotyping Errors and Harvested Bear

ABSTRACT Estimating black bear (Ursus americanus) population size is a difficult but important requirement when justifying harvest quotas and managing populations. Advancements in genetic techniques provide a means to identify individual bears using DNA contained in tissue and hair samples, thereby permitting estimates of population abundance based on established mark-capture-recapture methodology. We expand on previous noninvasive population-estimation work by geographically extending sampling areas (36,848 km²) to include the entire Northern Lower Peninsula (NLP) of Michigan, USA. We selected sampling locations randomly within biologically relevant bear habitat and used barbed wire hair snares to collect hair samples. Unlike previous noninvasive studies, we used tissue samples from harvested bears as an additional sampling occasion to increase recapture probabilities. We developed subsampling protocols to account for both spatial and temporal variance in sample distribution and variation in sample quality using recently published quality control protocols using 5 microsatellite loci. We quantified genotyping errors using samples from harvested bears and estimated abundance using statistical models that accounted for genotyping error. We estimated the population of yearling and adult black bears in the NLP to be 1,882 bears (95% CI = 1,389-2,551 bears). The derived population estimate with a 15% coefficient of variation was used by wildlife managers to examine the sustainability of harvest over a large geographic area.     

Spatial variation in density of American black bears in northern Yellowstone National Park

The quality and availability of resources are known to influence spatial patterns of animal density. In Yellowstone National Park, relationships between the availability of resources and the distribution of grizzly bears (Ursus arctos) have been explored but have yet to be examined in American black bears (Ursus americanus). We conducted non-invasive genetic sampling during 2017–2018 (mid-May to mid-July) and applied spatially explicit capture-recapture models to estimate density of black bears and examine associations with landscape features. In both years, density estimates were higher in forested vegetation communities, which provide food resources and thermal and security cover preferred by black bears, compared with non-forested areas. In 2017, density also varied by sex, with female densities being higher than males. Based on our estimates, the northern range of Yellowstone National Park supports one of the highest densities of black bears (20 black bears/100 km²) in the northern Rocky Mountains (6–12 black bears/100 km² in other regions). Given these high densities, black bears could influence other wildlife populations more than previously thought, such as through displacement of sympatric predators from kills. Our study provides the first spatially explicit estimates of density for black bears within an ecosystem that contains the majority of North America's large mammal species. Our density estimates provide a baseline that can be used for future research and management decisions of black bears, including efforts to reduce human–bear conflicts.     

DNA-Based Population Demographics of Black Bears in Coastal North Carolina and Virginia

ABSTRACT Noninvasive genetic sampling has become a popular method for obtaining population parameter estimates for black (Ursus americanus) and brown (U. arctos) bears. These estimates allow wildlife managers to develop appropriate management strategies for populations of concern. Black bear populations at Great Dismal Swamp (GDSNWR), Pocosin Lakes (PLNWR), and Alligator River (ARNWR) National Wildlife Refuges in coastal Virginia and North Carolina, USA, were perceived by refuge biologists to be at or above cultural and perhaps biological carrying capacity, but managers had no reliable abundance estimates upon which to base population management. We derived density estimates from 3,150 hair samples collected noninvasively at each of the 3 refuges, using 6–7 microsatellite markers to obtain multilocus genotypes for individual bears. We used Program MARK to calculate population estimates from capture histories at each refuge. We estimated densities using both traditional buffer strip methods and Program DENSITY. Estimated densities were some of the highest reported in the literature and ranged from 0.46 bears/km² at GDSNWR to 1.30 bears/km² at PLNWR. Sex ratios were male-biased at all refuges. Our estimates can be directly utilized by biologists to develop effective strategies for managing and maintaining bears at these refuges, and noninvasive methods may also be effective for monitoring bear populations over the long term.     

GRIZZLY BEAR DEMOGRAPHICS IN AND AROUND BANFF NATIONAL PARK AND KANANASKIS COUNTRY,ALBERTA

Abstract: The area in and around Banff National Park (BNP) in southwestern Alberta, Canada, is 1 of the most heavily used and developed areas where grizzly bears (Ursus arctos) still exist. During 1994–2002, we radiomarked and monitored 37 female and 34 male bears in this area to estimate rates of survival, reproduction, and population growth. Annual survival rates of bears other than dependent young averaged 95% for females and 81–85% for males. Although this area was largely unhunted, humans caused 75% of female mortality and 86% of male mortality. Females produced their first surviving litter at 6–12 years of age (x̄ = 8.4 years). Litters averaged 1.84 cubs spaced at 4.4-year intervals. Adult (≥6-years-old) females produced 0.24 female cubs per year and were expected to produce an average of 1.7 female cubs in their lifetime, based on rates of reproduction and survival. Cub survival was 79%, yearling survival was 91%, and survival through independence at 2.5–5.5 years of age was 72%, as no dependent young older than yearlings died. Although this is the slowest-reproducing grizzly bear population yet studied, high rates of survival seem to have enabled positive population growth (Λ = 1.04, 95% CI = 0.99–1.09), based on analyses using Leslie matrices. Current management practices, instituted in the late 1980s, focus on alleviating human-caused bear mortality. If the 1970–1980s style of management had continued, we estimated that an average of 1 more radiomarked female would have been killed each year, reducing female survival to the point that the population would have declined.     

Evaluation of noninvasive genetic sampling methods for cougars in Yellowstone National Park

Conventional methods for monitoring cougar, Puma concolor, populations involve capture, tagging, and radio-collaring, but these methods are time-consuming, expensive, and logistically challenging. For difficult-to-study species such as cougars, noninvasive genetic sampling (NGS) may be a useful alternative. The ability to identify individuals from samples collected through NGS methods provides many opportunities for developing population-monitoring tools, but the utility of these survey methods is dependent upon collection of samples and accurate genotyping of those samples. In January 2003, we initiated a 3-yr evaluation of NGS methods for cougars using a radio-collared population in Yellowstone National Park (YNP), USA. Our goals were to: 1) determine which DNA collection method, hair snares or snow tracking, provided a better method for obtaining samples for genetic analysis, 2) evaluate reliability of the genetic data derived from hair samples collected in the field, and 3) evaluate the potential of NGS for demographic monitoring of cougar populations. Snow tracking yielded more hair samples and was more cost effective than snagging hair with rub pads. Samples collected from bed sites and natural hair snags (e.g., branch tips, thorn bushes) while snow tracking accurately identified and sexed 22 individuals (9 F, 13 M). The ratio of the count from snow tracking to the count from radio-telemetry was 15:24 in winter 2004, 13:12 in 2005, and 22:29 for both years combined. Annual capture probabilities for obtaining DNA from snow tracking varied considerably between years for females (0.42 in 2004 and 0.88 in 2005) but were more consistent for males (0.77 in 2004 and 0.88 in 2005). Our results indicate that snow tracking can be an efficient, reliable NGS method for cougars in YNP and has potential for estimating demographic and genetic parameters of other carnivore populations in similar climates. © 2011 The Wildlife Society.