Grizzly Bear Density in Glacier National Park,Montana

Abstract: We present the first rigorous estimate of grizzly bear (Ursus arctos) population density and distribution in and around Glacier National Park (GNP), Montana, USA. We used genetic analysis to identify individual bears from hair samples collected via 2 concurrent sampling methods: 1) systematically distributed, baited, barbed-wire hair traps and 2) unbaited bear rub trees found along trails. We used Huggins closed mixture models in Program MARK to estimate total population size and developed a method to account for heterogeneity caused by unequal access to rub trees. We corrected our estimate for lack of geographic closure using a new method that utilizes information from radiocollared bears and the distribution of bears captured with DNA sampling. Adjusted for closure, the average number of grizzly bears in our study area was 240.7 (95% CI = 202–303) in 1998 and 240.6 (95% CI = 205–304) in 2000. Average grizzly bear density was 30 bears/1,000 km², with 2.4 times more bears detected per hair trap inside than outside GNP. We provide baseline information important for managing one of the few remaining populations of grizzlies in the contiguous United States.     

Estimates of Abundance and Harvest Rates of Female Black Bears Across a Large Spatial Extent

American black bears (Ursus americanus) are an iconic wildlife species in the southern Appalachian highlands of the eastern United States and have increased in number and range since the early 1980s. Given an increasing number of human-bear conflicts in the region, many management agencies have liberalized harvest regulations to reduce bear populations to socially acceptable levels. Wildlife managers need reliable population data for assessing the effects of management actions for this high-profile species. Our goal was to use DNA extracted from hair collected at barbed-wire enclosures (i.e., hair traps) to identify individual bears and then use spatially explicit capture-recapture methods to estimate female black bear density, abundance, and harvest rate. We established 888 hair traps across 66,678 km² of the southern Appalachian highlands in Georgia, North Carolina, South Carolina, and Tennessee, USA, in 2017 and 2018, arranged in 174 clusters of 2–9 traps/cluster. We collected 9,113 hair samples from those sites over 6 weeks of sampling, of which 1,954 were successfully genotyped to 462 individual female bears. Our spatially explicit estimator included a percent forest covariate to explain inhomogeneous bear density across the region. Densities ranged up to 0.410 female bears/km² and regional abundance was 5,950 (95% CI = 4,988–7,098) female bears. Based on hunter kill data from 2016 to 2018, mean annual harvest rates for females were 12.7% in Georgia, 17.6% in North Carolina, 17.6% in South Carolina, and 22.8% in Tennessee. Our estimated harvest rates for most states approached or exceeded theoretical maximum sustainable levels, and population trend data (i.e., bait-station indices) indicated decreasing growth rates since about 2009. These data suggest that the increased harvest goals and poor hard mast production over a series of prior years reduced bear population abundance in many states. We were able to obtain reasonable population abundance and density estimates because of spatially explicit capture-recapture methods, cluster sampling, and a large spatial extent. Continued monitoring of bear populations (e.g., annual bait-station surveys and periodic population estimation using spatially explicit methods) by state jurisdictions would help to ensure that population trajectories are consistent with management goals. © 2021 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.     

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 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.     

A Comparison of Grizzly Bear Demographic Parameters Estimated from Non-Spatial and Spatial Open Population Capture-Recapture Models

Capture-recapture studies are frequently used to monitor the status and trends of wildlife populations. Detection histories from individual animals are used to estimate probability of detection and abundance or density. The accuracy of abundance and density estimates depends on the ability to model factors affecting detection probability. Non-spatial capture-recapture models have recently evolved into spatial capture-recapture models that directly include the effect of distances between an animal’s home range centre and trap locations on detection probability. Most studies comparing non-spatial and spatial capture-recapture biases focussed on single year models and no studies have compared the accuracy of demographic parameter estimates from open population models. We applied open population non-spatial and spatial capture-recapture models to three years of grizzly bear DNA-based data from Banff National Park and simulated data sets. The two models produced similar estimates of grizzly bear apparent survival, per capita recruitment, and population growth rates but the spatial capture-recapture models had better fit. Simulations showed that spatial capture-recapture models produced more accurate parameter estimates with better credible interval coverage than non-spatial capture-recapture models. Non-spatial capture-recapture models produced negatively biased estimates of apparent survival and positively biased estimates of per capita recruitment. The spatial capture-recapture grizzly bear population growth rates and 95% highest posterior density averaged across the three years were 0.925 (0.786–1.071) for females, 0.844 (0.703–0.975) for males, and 0.882 (0.779–0.981) for females and males combined. The non-spatial capture-recapture population growth rates were 0.894 (0.758–1.024) for females, 0.825 (0.700–0.948) for males, and 0.863 (0.771–0.957) for both sexes. The combination of low densities, low reproductive rates, and predominantly negative population growth rates suggest that Banff National Park’s population of grizzly bears requires continued conservation-oriented management actions.     

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.     

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.     

Contrasting Activity Patterns of Sympatric and Allopatric Black and Grizzly Bears

ABSTRACT The distribution of grizzly (Ursus arctos) and American black bears (U. americanus) overlaps in western North America. Few studies have detailed activity patterns where the species are sympatric and no studies contrasted patterns where populations are both sympatric and allopatric. We contrasted activity patterns for sympatric black and grizzly bears and for black bears allopatric to grizzly bears, how human influences altered patterns, and rates of grizzly-black bear predation. Activity patterns differed between black bear populations, with those sympatric to grizzly bears more day-active. Activity patterns of black bears allopatric with grizzly bears were similar to those of female grizzly bears; both were crepuscular and day-active. Male grizzly bears were crepuscular and night-active. Both species were more night-active and less day-active when ≤1 km from roads or developments. In our sympatric study area, 2 of 4 black bear mortalities were due to grizzly bear predation. Our results suggested patterns of activity that allowed for intra- and inter-species avoidance. National park management often results in convergence of locally high human densities in quality bear habitat. Our data provide additional understanding into how bears alter their activity patterns in response to other bears and humans and should help park managers minimize undesirable bear-human encounters when considering needs for temporal and spatial management of humans and human developments in bear habitats.     

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.     

Factors associated with black bear density and implications for management

Wildlife density estimates are important to accurately formulate population management objectives and understand the relationship between habitat characteristics and a species’ abundance. Despite advances in density and abundance estimation methods, management of common game species continues to be challenged by a lack of reliable population estimates. In Washington, USA, statewide American black bear (Ursus americanus) abundance estimates are predicated on density estimates derived from research in the 1970s and are hypothesized to be a function of precipitation and vegetation, with higher densities in western Washington. To evaluate current black bear density and landscape relationships in Washington, we conducted a 4-year capture-recapture study in 2 areas of the North Cascade Mountains using 2 detection methods, non-invasive DNA collection and physical capture and deployment of global positioning system (GPS) collars. We integrated GPS telemetry from collared bears with spatial capture-recapture (SCR) data and created a SCR-resource selection model to estimate density as a function of spatial covariates and test the hypothesis that density is higher in areas with greater vegetative food resources. We captured and collared 118 bears 132 times and collected 7,863 hair samples at hair traps where we identified 537 bears from 1,237 detections via DNA. The most-supported model in the western North Cascades depicted a negative relationship between black bear density and an index of human development. We estimated bear density at 20.1 bears/100 km², but density varied from 13.5/100 km² to 27.8 bears/100 km² depending on degree of human development. The model best supported by the data in the eastern North Cascades estimated an average density of 19.2 bears/100 km², which was positively correlated with primary productivity, with resulting density estimates ranging from 7.1/100 km² to 33.6 bears/100 km². The hypothesis that greater precipitation and associated vegetative production in western Washington supports greater bear density compared to eastern Washington was not supported by our data. In western Washington, empirically derived average density estimates (including cubs) were nearly 50% lower than managers expected prior to our research. In eastern Washington average black bear density was predominantly as expected, but localized areas of high primary productivity supported greater than anticipated bear densities. Our findings underscore the importance that black bear density is not likely uniform and management risk may be increased if an average density is applied at too large a scale. Disparities between expected and empirically derived bear density illustrate the need for more rigorous monitoring to understand processes that affect population numbers throughout the jurisdiction, and suggest that management plans may need to be reevaluated to determine if current harvest strategies are achieving population objectives. © 2019 The Wildlife Society.     

Balancing Precision and Risk: Should Multiple Detection Methods Be Analyzed Separately in N-Mixture Models?

Using multiple detection methods can increase the number, kind, and distribution of individuals sampled, which may increase accuracy and precision and reduce cost of population abundance estimates. However, when variables influencing abundance are of interest, if individuals detected via different methods are influenced by the landscape differently, separate analysis of multiple detection methods may be more appropriate. We evaluated the effects of combining two detection methods on the identification of variables important to local abundance using detections of grizzly bears with hair traps (systematic) and bear rubs (opportunistic). We used hierarchical abundance models (N-mixture models) with separate model components for each detection method. If both methods sample the same population, the use of either data set alone should (1) lead to the selection of the same variables as important and (2) provide similar estimates of relative local abundance. We hypothesized that the inclusion of 2 detection methods versus either method alone should (3) yield more support for variables identified in single method analyses (i.e. fewer variables and models with greater weight), and (4) improve precision of covariate estimates for variables selected in both separate and combined analyses because sample size is larger. As expected, joint analysis of both methods increased precision as well as certainty in variable and model selection. However, the single-method analyses identified different variables and the resulting predicted abundances had different spatial distributions. We recommend comparing single-method and jointly modeled results to identify the presence of individual heterogeneity between detection methods in N-mixture models, along with consideration of detection probabilities, correlations among variables, and tolerance to risk of failing to identify variables important to a subset of the population. The benefits of increased precision should be weighed against those risks. The analysis framework presented here will be useful for other species exhibiting heterogeneity by detection method.     

Challenges of DNA-Based Mark-Recapture Studies of American Black Bears

Abstract: We explored whether genetic sampling would be feasible to provide a region-wide population estimate for American black bears (Ursus americanus) in the southern Appalachians, USA. Specifically, we determined whether adequate capture probabilities (p > 0.20) and population estimates with a low coefficient of variation (CV < 20%) could be achieved given typical agency budget and personnel constraints. We extracted DNA from hair collected from baited barbed-wire enclosures sampled over a 10-week period on 2 study areas: a high-density black bear population in a portion of Great Smoky Mountains National Park and a lower density population on National Forest lands in North Carolina, South Carolina, and Georgia. We identified individual bears by their unique genotypes obtained from 9 microsatellite loci. We sampled 129 and 60 different bears in the National Park and National Forest study areas, respectively, and applied closed mark-recapture models to estimate population abundance. Capture probabilities and precision of the population estimates were acceptable only for sampling scenarios for which we pooled weekly sampling periods. We detected capture heterogeneity biases, probably because of inadequate spatial coverage by the hair-trapping grid. The logistical challenges of establishing and checking a sufficiently high density of hair traps make DNA-based estimates of black bears impractical for the southern Appalachian region. Alternatives are to estimate population size for smaller areas, estimate population growth rates or survival using mark-recapture methods, or use independent marking and recapturing techniques to reduce capture heterogeneity.     

Using spatially‐explicit capture–recapture models to explain variation in seasonal density patterns of sympatric ursids

Understanding how environmental factors interact to determine the abundance and distribution of animals is a primary goal of ecology, and fundamental to the conservation of wildlife populations. Studies of these relationships, however, often assume static environmental conditions, and rarely consider effects of competition with ecologically similar species. In many parts of their shared ranges, grizzly bears Ursus arctos and American black bears U. americanus have nearly complete dietary overlap and share similar life history traits. We therefore tested the hypothesis that density patterns of both bear species would reflect seasonal variation in available resources, with areas of higher primary productivity supporting higher densities of both species. We also hypothesized that interspecific competition would influence seasonal density patterns. Specifically, we predicted that grizzly bear density would be locally reduced due to the ability of black bears to more efficiently exploit patchy food resources such as seasonally abundant fruits. To test our hypotheses, we used detections of 309 grizzly and 597 black bears from two independent genetic sampling methods in spatially‐explicit capture–recapture (SECR) models. Our results suggest grizzly bear density was lower in areas of high black bear density during spring and summer, although intraspecific densities were also important, particularly during the breeding season. Black bears had lower densities in areas of high grizzly bear density in spring; however, density of black bears in early and late summer was best explained by primary productivity. Our results are consistent with the hypothesis that smaller‐bodied, more abundant black bears may influence the density patterns of behaviorally‐dominant grizzly bears through exploitative competition. We also suggest that seasonal variation in resource availability be considered in efforts to relate environmental conditions to animal density.     

No long-term effect of black bear removal on elk calf recruitment in the southern Appalachians

In 2001 and 2002, 52 elk (Cervus canadensis; 21 males, 31 females), originally obtained from Elk Island National Park, Alberta, Canada, were transported and released into Cataloochee Valley in the northeastern portion of Great Smoky Mountains National Park (GRSM, Park), North Carolina, USA. The annual population growth rate (λ) was negative (0.996, 95% CI = 0.945–1.047) and predation by black bears (Ursus americanus) on elk calves was identified as an important determinant of population growth. From 2006 to 2008, 49 bears from the primary elk calving area (i.e., Cataloochee Valley) were trapped and translocated about 70 km to the southwestern portion of the Park just prior to elk calving. Per capita recruitment (i.e., the number of calves produced per adult female that survive to 1 year of age) increased from 0.306 prior to bear translocation (2001–2005) to 0.544 during years when bears were translocated (2006–2008) and λ increased to 1.118 (95% CI = 1.096–1.140). Our objective was to determine whether per capita calf recruitment rates after bear removal (2009–2019) at Cataloochee were similar to the higher rates estimated during bear removal (i.e., long-term response) or if they returned to rates before bear removal (i.e., short-term response), and how those rates compared with recruitment from portions of our study area where bears were not relocated. We documented 419 potential elk calving events and monitored 129 yearling and adult elk from 2001 to 2019. Known-fate models based on radio-telemetry and observational data supported calf recruitment returning to pre-2006 levels at Cataloochee (short-term response); recruitment of Cataloochee elk before and after bear relocation was lower (0.184) than during bear relocation (0.492). Recruitment rates of elk outside the removal area during the bear relocation period (0.478) were similar to before and after rates (0.420). In the Cataloochee Valley, cause-specific annual calf mortality rates due to predation by bears were 0.319 before, 0.120 during, and 0.306 after bear relocation. In contrast, the cause-specific annual mortality rate of calves in areas where bears were not relocated was 0.033 after the bear relocation period, with no bear predation on calves before or during bear relocation. The mean annual population growth rate for all monitored elk was 1.062 (95% CI = 0.979–1.140) after bear relocation based on the recruitment and survival data. Even though the effects of bear removal were temporary, the relocations were effective in achieving a short-term increase in elk recruitment, which was important for the reintroduction program given that the elk population was small and vulnerable to extirpation.     

Effects of black bear relocation on elk calf recruitment at Great Smoky Mountains National Park

Previous research from 2001 to 2006 on an experimentally released elk (Cervus elaphus) population at Great Smoky Mountains National Park (GSMNP or Park) indicated that calf recruitment (i.e., calves reaching 1 yr of age per adult female elk) was low (0.306, total SE = 0.090) resulting in low or negative population growth (λ = 0.996, 95% CI = 0.945–1.047). Black bear (Ursus americanus) predation was the primary calf mortality factor. From 2006 to 2008, we trapped and relocated 49 bears (30 of which were radiocollared) from the primary calving areas in the Park and radiomonitored 67 (28 M:39 F) adult elk and 42 calves to compare vital rates and population growth with the earlier study. A model with annual calf recruitment rate correlating with the number of bears relocated each year was supported (ΔAIC_c = 0.000; β = 0.070, 95% CI = 0.028–0.112) and a model with annual calf recruitment differing from before to during bear relocation revealed an increase to 0.544 (total SE = 0.098; β = −1.092, 95% CI = −1.180 to −0.375). Using vital rates and estimates of process standard errors observed during our study, 25-yr simulations maintained a mean positive growth rate in 100% of the stochastic trials with λ averaging 1.118 (95% CI = 1.096–1.140), an increase compared with rates before bear relocation. A life table response experiment revealed that increases in population growth were mostly (67.1%) due to changes in calf recruitment. We speculate that behavioral adaptation of the elk since release also contributed to the observed increases in recruitment and population growth. Our results suggest that managers interested in elk reintroduction within bear range should consider bear relocation as a temporary means of increasing calf recruitment. © 2011 The Wildlife Society.     

Prevalence of Trichinella spp. in black bears, grizzly bears, and wolves in the Dehcho Region, Northwest Territories, Canada, including the first report of T. nativa in a grizzly bear from Canada

Samples of muscle from 120 black bears (Ursus americanus), 11 grizzly bears (Ursus arctos), and 27 wolves (Canis lupus) collected in the Dehcho Region of the Northwest Territories from 2001 to 2010 were examined for the presence of Trichinella spp. larvae using a pepsin-HCl digestion assay. Trichinella spp. larvae were found in eight of 11 (73%) grizzly bears, 14 of 27 (52%) wolves, and seven of 120 (5.8%) black bears. The average age of positive grizzly bears, black bears, and wolves was 13.5, 9.9, and approximately 4 yr, respectively. Larvae from 11 wolves, six black bears, and seven grizzly bears were genotyped. Six wolves were infected with T. nativa and five with Trichinella T6, four black bears were infected with T. nativa and two with Trichinella T6, and all seven grizzly bears were infected with Trichinella T6 and one of them had a coinfection with T. nativa. This is the first report of T. nativa in a grizzly bear from Canada. Bears have been linked to trichinellosis outbreaks in humans in Canada, and black bears are a subsistence food source for residents of the Dehcho region. In order to assess food safety risk it is important to monitor the prevalence of Trichinella spp. in both species of bear and their cohabiting mammalian food sources.     

Assessing Population Viability of Black Bears using Spatial Capture-Recapture Models

The Central Georgia Bear Population (CGP) is the least abundant and most isolated of Georgia's 3 American black bear (Ursus americanus) populations. Beginning in 2011, changes to regulations governing harvest of the CGP resulted in an increase in female bear harvest, creating concern that future harvest could be an important influence on population viability. Hence, our objective was to assess viability of the CGP under various levels of female mortality. During 2012–2016, we used barbed-wire hair snares to collect bear hair samples from within the range of the CGP in Georgia, USA. We used microsatellite genotyping to identify individual bears and created robust-design, spatial detection histories for all female bears detected. We fit open population spatial capture-recapture (SCR) models to the detection histories in a Bayesian framework. We used the Widely Applicable Information Criterion (WAIC) to rank models that varied with respect to sources of variation in detection probability, survival, and per capita recruitment, and used the model with the lowest WAIC to forecast dynamics of the CGP 50 years into the future under various levels of female mortality. We assessed the 50-year extinction probability under a continuation of mortality levels documented during 2012–2016, and under incremental increases in female mortality above this baseline. The top model included density-dependent per capita recruitment, annual variation in detection probability, and a trap-level behavioral response. Abundance increased from 106 (95% CI = 86–132) females in 2012 to 136 (95% CI = 113–161) females in 2013 and remained relatively stable thereafter. Annual female survival was 0.75 (95% CI = 0.69–0.82) and did not vary among years. The per capita recruitment rate decreased over time as density increased, and was 0.49 (95% CI = 0.33–0.66) during the first time interval and 0.29 (95% CI = 0.20–0.38) during the final time interval. Annual growth rate () was 1.28 (95% CI = 1.07–1.52) between 2012 and 2013 but decreased throughout the study, ending at 1.04 (95% CI = 0.93–1.17). Forecasts indicated continuation of the female mortality levels experienced from 2012–2016 were sustainable over 50 years, with the estimated extinction risk being <0.001%. Increasing annual harvest by 5 females introduced a negligible increase in the 50-year probability of extinction, but harvesting an additional 10 females/year caused extinction risk to rise to 1.15%. We recommend that harvest regulations are structured such that mortality rates remain at current levels or do not increase by more than an annual average of 5 females above levels observed during our study. Furthermore, we recommend that managers continue to monitor the population so that harvest regulations and population models can be refined over time. © 2020 The Wildlife Society.     

A study of predominant aerobic microflora of black bears (Ursus americanus) and grizzly bears (Ursus arctos) in northwestern Alberta

Swab specimens were obtained from nasal, rectal, and preputial or vaginal areas of 37 grizzly and 17 black bears, captured during May to June of 1981 to 1983, to determine the types and frequency of predominant aerobic microflora. Bacterial genera most frequently isolated from bears were Escherichia, Citrobacter, Hafnia, Proteus, Staphylococcus, and Streptococcus species, comprising about 65% of the isolates. Erwinia, Xanthomonas, Agrobacterium, Rhizobium, and Gluconobacter/Acetobacter were also isolated but at lower frequencies (less than 5%). Comparison of bacterial generic composition using similarity quotient values showed no appreciable differences between grizzly and black bear flora. Also, no outstanding differences in bacterial generic composition were observed among grizzly bear samples; however, differences were noted among black bear samples. Fungal genera most commonly encountered included Cryptococcus, Rhodotorula, Cladosporium, Penicillium, Sporobolomyces, and Candida. In general, the microflora of both bear types were marked by generic diversity and random distribution. The majority of microorganisms isolated from the plant samples in the study area were also found in bear samples. This observation and the presence of certain water and soil bacteria in samples from bears suggest that the predominant microflora of both grizzly and black bears were transient and probably influenced by their foraging habits and surrounding environments.     

The Impact of Roads on the Demography of Grizzly Bears in Alberta

One of the principal factors that have reduced grizzly bear populations has been the creation of human access into grizzly bear habitat by roads built for resource extraction. Past studies have documented mortality and distributional changes of bears relative to roads but none have attempted to estimate the direct demographic impact of roads in terms of both survival rates, reproductive rates, and the interaction of reproductive state of female bears with survival rate. We applied a combination of survival and reproductive models to estimate demographic parameters for threatened grizzly bear populations in Alberta. Instead of attempting to estimate mean trend we explored factors which caused biological and spatial variation in population trend. We found that sex and age class survival was related to road density with subadult bears being most vulnerable to road-based mortality. A multi-state reproduction model found that females accompanied by cubs of the year and/or yearling cubs had lower survival rates compared to females with two year olds or no cubs. A demographic model found strong spatial gradients in population trend based upon road density. Threshold road densities needed to ensure population stability were estimated to further refine targets for population recovery of grizzly bears in Alberta. Models that considered lowered survival of females with dependant offspring resulted in lower road density thresholds to ensure stable bear populations. Our results demonstrate likely spatial variation in population trend and provide an example how demographic analysis can be used to refine and direct conservation measures for threatened species.