Drivers of habitat availability for terrestrial mammals: Unravelling the role of livestock, land conversion and intrinsic traits in the past 50 years.

The global decline of terrestrial species is largely due to the degradation, loss and fragmentation of their habitats. The conversion of natural ecosystems for cropland, rangeland, forest products and human infrastructure are the primary causes of habitat deterioration. Due to the paucity of data on the past distribution of species and the scarcity of fine-scale habitat conversion maps, however, accurate assessment of the recent effects of habitat degradation, loss and fragmentation on the range of mammals has been near impossible. We aim to assess the proportions of available habitat within the lost and retained parts of mammals' distribution ranges, and to identify the drivers of habitat availability. We produced distribution maps for 475 terrestrial mammals for the range they occupied 50 years ago and compared them to current range maps. We then calculated the differences in the percentage of 'area of habitat' (habitat available to a species within its range) between the lost and retained range areas. Finally, we ran generalized linear mixed models to identify which variables were more influential in determining habitat availability in the lost and retained parts of the distribution ranges. We found that 59% of species had a lower proportion of available habitat in the lost range compared to the retained range, thus hypothesizing that habitat loss could have contributed to range declines. The most important factors negatively affecting habitat availability were the conversion of land to rangeland and high density of livestock. Significant intrinsic traits were those related to reproductive timing and output, habitat breadth and medium body size. Our findings emphasize the importance of implementing conservation strategies to mitigate the impacts caused by human activities on the habitats of mammals, and offer evidence indicating which species have the potential to reoccupy portions of their former range if other threats cease to occur.

Evaluating expert-based habitat suitability information of terrestrial mammals with GPS-tracking data.

Aim: Macroecological studies that require habitat suitability data for many species often derive this information from expert opinion. However, expert-based information is inherently subjective and thus prone to errors. The increasing availability of GPS tracking data offers opportunities to evaluate and supplement expert-based information with detailed empirical evidence. Here, we compared expert-based habitat suitability information from the International Union for Conservation of Nature (IUCN) with habitat suitability information derived from GPS-tracking data of 1,498 individuals from 49 mammal species. Location: Worldwide. Time period: 1998-2021. Major taxa studied: Forty-nine terrestrial mammal species. Methods: Using GPS data, we estimated two measures of habitat suitability for each individual animal: proportional habitat use (proportion of GPS locations within a habitat type), and selection ratio (habitat use relative to its availability). For each individual we then evaluated whether the GPS-based habitat suitability measures were in agreement with the IUCN data. To that end, we calculated the probability that the ranking of empirical habitat suitability measures was in agreement with IUCN's classification into suitable, marginal and unsuitable habitat types. Results: IUCN habitat suitability data were in accordance with the GPS data (> 95% probability of agreement) for 33 out of 49 species based on proportional habitat use estimates and for 25 out of 49 species based on selection ratios. In addition, 37 and 34 species had a > 50% probability of agreement based on proportional habitat use and selection ratios, respectively. Main conclusions: We show how GPS-tracking data can be used to evaluate IUCN habitat suitability data. Our findings indicate that for the majority of species included in this study, it is appropriate to use IUCN habitat suitability data in macroecological studies. Furthermore, we show that GPS-tracking data can be used to identify and prioritize species and habitat types for re-evaluation of IUCN habitat suitability data.

Integrating data types to estimate spatial patterns of avian migration across the Western Hemisphere.

For many avian species, spatial migration patterns remain largely undescribed, especially across hemispheric extents. Recent advancements in tracking technologies and high-resolution species distribution models (i.e., eBird Status and Trends products) provide new insights into migratory bird movements and offer a promising opportunity for integrating independent data sources to describe avian migration. Here, we present a three-stage modeling framework for estimating spatial patterns of avian migration. First, we integrate tracking and band re-encounter data to quantify migratory connectivity, defined as the relative proportions of individuals migrating between breeding and nonbreeding regions. Next, we use estimated connectivity proportions along with eBird occurrence probabilities to produce probabilistic least-cost path (LCP) indices. In a final step, we use generalized additive mixed models (GAMMs) both to evaluate the ability of LCP indices to accurately predict (i.e., as a covariate) observed locations derived from tracking and band re-encounter data sets versus pseudo-absence locations during migratory periods and to create a fully integrated (i.e., eBird occurrence, LCP, and tracking/band re-encounter data) spatial prediction index for mapping species-specific seasonal migrations. To illustrate this approach, we apply this framework to describe seasonal migrations of 12 bird species across the Western Hemisphere during pre- and postbreeding migratory periods (i.e., spring and fall, respectively). We found that including LCP indices with eBird occurrence in GAMMs generally improved the ability to accurately predict observed migratory locations compared to models with eBird occurrence alone. Using three performance metrics, the eBird + LCP model demonstrated equivalent or superior fit relative to the eBird-only model for 22 of 24 species-season GAMMs. In particular, the integrated index filled in spatial gaps for species with over-water movements and those that migrated over land where there were few eBird sightings and, thus, low predictive ability of eBird occurrence probabilities (e.g., Amazonian rainforest in South America). This methodology of combining individual-based seasonal movement data with temporally dynamic species distribution models provides a comprehensive approach to integrating multiple data types to describe broad-scale spatial patterns of animal movement. Further development and customization of this approach will continue to advance knowledge about the full annual cycle and conservation of migratory birds.

Life-history traits and habitat availability shape genomic diversity in birds: implications for conservation.

More than 25% of species assessed by the International Union for Conservation of Nature (IUCN) are threatened with extinction. Understanding how environmental and biological processes have shaped genomic diversity may inform management practices. Using 68 extant avian species, we parsed the effects of habitat availability and life-history traits on genomic diversity over time to provide a baseline for conservation efforts. We used published whole-genome sequence data to estimate overall genomic diversity as indicated by historical long-term effective population sizes (Ne) and current genomic variability (H), then used environmental niche modelling to estimate Pleistocene habitat dynamics for each species. We found that Ne and H were positively correlated with habitat availability and related to key life-history traits (body mass and diet), suggesting the latter contribute to the overall genomic variation. We found that H decreased with increasing species extinction risk, suggesting that H may serve as a leading indicator of demographic trends related to formal IUCN conservation status in birds. Our analyses illustrate that genome-wide summary statistics estimated from sequence data reflect meaningful ecological attributes relevant to species conservation.

A Translational Model for Venous Thromboembolism: MicroRNA Expression in Hibernating Black Bears.

BACKGROUND: Hibernating American black bears have significantly different clotting parameters than their summer active counterparts, affording them protection against venous thromboembolism during prolonged periods of immobility. We sought to evaluate if significant differences exist between the expression of microRNAs in the plasma of hibernating black bears compared with their summer active counterparts, potentially contributing to differences in hemostasis during hibernation. MATERIALS AND METHODS: MicroRNA sequencing was assessed in plasma from 21 American black bears in summer active (n = 11) and hibernating states (n = 10), and microRNA signatures during hibernating and active state were established using both bear and human genome. MicroRNA targets were predicted using messenger RNA (mRNA) transcripts from black bear kidney cells. In vitro studies were performed to confirm the relationship between identified microRNAs and mRNA expression, using artificial microRNA and human liver cells. RESULTS: Using the bear genome, we identified 15 microRNAs differentially expressed in the plasma of hibernating black bears. Of these microRNAs, three were significantly downregulated (miR-141-3p, miR-200a-3p, and miR-200c-3p), were predicted to target SERPINC1, the gene for antithrombin, and demonstrated regulatory control of the gene mRNA expression in cell studies. CONCLUSIONS: Our findings suggest that the hibernating black bears' ability to maintain hemostasis and achieve protection from venous thromboembolism during prolonged periods of immobility may be due to changes in microRNA signatures and possible upregulation of antithrombin expression.

Effects of body size on estimation of mammalian area requirements.

Accurately quantifying species' area requirements is a prerequisite for effective area-based conservation. This typically involves collecting tracking data on species of interest and then conducting home-range analyses. Problematically, autocorrelation in tracking data can result in space needs being severely underestimated. Based on the previous work, we hypothesized the magnitude of underestimation varies with body mass, a relationship that could have serious conservation implications. To evaluate this hypothesis for terrestrial mammals, we estimated home-range areas with global positioning system (GPS) locations from 757 individuals across 61 globally distributed mammalian species with body masses ranging from 0.4 to 4000 kg. We then applied block cross-validation to quantify bias in empirical home-range estimates. Area requirements of mammals <10 kg were underestimated by a mean approximately15%, and species weighing approximately100 kg were underestimated by approximately50% on average. Thus, we found area estimation was subject to autocorrelation-induced bias that was worse for large species. Combined with the fact that extinction risk increases as body mass increases, the allometric scaling of bias we observed suggests the most threatened species are also likely to be those with the least accurate home-range estimates. As a correction, we tested whether data thinning or autocorrelation-informed home-range estimation minimized the scaling effect of autocorrelation on area estimates. Data thinning required an approximately93% data loss to achieve statistical independence with 95% confidence and was, therefore, not a viable solution. In contrast, autocorrelation-informed home-range estimation resulted in consistently accurate estimates irrespective of mass. When relating body mass to home range size, we detected that correcting for autocorrelation resulted in a scaling exponent significantly >1, meaning the scaling of the relationship changed substantially at the upper end of the mass spectrum.

Efectos del Tamaño Corporal sobre la Estimación de los Requerimientos de Área de Mamíferos Resumen La cuantificación precisa de los requerimientos de área de una especie es un prerrequisito para que la conservación basada en áreas sea efectiva. Esto comúnmente implica la recolección de datos de rastreo de la especie de interés para después realizar análisis de la distribución local. De manera problemática, la autocorrelación en los datos de rastreo puede resultar en una subestimación grave de las necesidades de espacio. Con base en trabajos previos, formulamos una hipótesis en la que supusimos que la magnitud de la subestimación varía con la masa corporal, una relación que podría tener implicaciones serias para la conservación. Para probar esta hipótesis en mamíferos terrestres, estimamos las áreas de distribución local con las ubicaciones en GPS de 757 individuos de 61 especies de mamíferos distribuidas mundialmente con una masa corporal entre 0.4 y 4,000 kg. Después aplicamos una validación cruzada en bloque para cuantificar el sesgo en estimaciones empíricas de la distribución local. Los requerimientos de área de los mamíferos <10 kg fueron subestimados por una media ∼15% y las especies con una masa ∼100 kg fueron subestimadas en ∼50% en promedio. Por lo tanto, encontramos que la estimación del área estaba sujeta al sesgo inducido por la autocorrelación, el cual era peor para las especies de talla grande. En combinación con el hecho de que el riesgo de extinción incrementa conforme aumenta la masa corporal, el escalamiento alométrico del sesgo que observamos sugiere que la mayoría de las especies amenazadas también tienen la probabilidad de ser aquellas especies con las estimaciones de distribución local menos acertadas. Como corrección, probamos si la reducción de datos o la estimación de la distribución local informada por la autocorrelación minimizan el efecto de escalamiento que tiene la autocorrelación sobre las estimaciones de área. La reducción de datos requirió una pérdida de datos del ∼93% para lograr la independencia estadística con un 95% de confianza y por lo tanto no fue una solución viable. Al contrario, la estimación de la distribución local informada por la autocorrelación resultó en estimaciones constantemente precisas sin importar la masa corporal. Cuando relacionamos la masa corporal con el tamaño de la distribución local, detectamos que la corrección de la autocorrelación resultó en un exponente de escalamiento significativamente >1, lo que significa que el escalamiento de la relación cambió sustancialmente en el extremo superior del espectro de la masa corporal.

Spatial genetic structure in American black bears (Ursus americanus): female philopatry is variable and related to population history.

Previously, American black bears (Ursus americanus) were thought to follow the pattern of female philopatry and male-biased dispersal. However, recent studies have identified deviations from this pattern. Such flexibility in dispersal patterns can allow individuals greater ability to acclimate to changing environments. We explored dispersal and spatial genetic relatedness patterns across ten black bear populations-including long established (historic), with known reproduction >50 years ago, and newly established (recent) populations, with reproduction recorded <50 years ago-in the Interior Highlands and Southern Appalachian Mountains, United States. We used spatially explicit, individual-based genetic simulations to model gene flow under scenarios with varying levels of population density, genetic diversity, and female philopatry. Using measures of genetic distance and spatial autocorrelation, we compared metrics between sexes, between population types (historic and recent), and among simulated scenarios which varied in density, genetic diversity, and sex-biased philopatry. In empirical populations, females in recent populations exhibited stronger patterns of isolation-by-distance (IBD) than females and males in historic populations. In simulated populations, low-density populations had a stronger indication of IBD than medium- to high-density populations; however, this effect varied in empirical populations. Condition-dependent dispersal strategies may permit species to cope with novel conditions and rapidly expand populations. Pattern-process modeling can provide qualitative and quantitative means to explore variable dispersal patterns, and could be employed in other species, particularly to anticipate range shifts in response to changing climate and habitat conditions.

Impacts of biomass production at civil airports on grassland bird conservation and aviation strike risk.

Growing concerns about climate change, foreign oil dependency, and environmental quality have fostered interest in perennial native grasses (e.g., switchgrass [Panicum virgatum]) for bioenergy production while also maintaining biodiversity and ecosystem function. However, biomass cultivation in marginal landscapes such as airport grasslands may have detrimental effects on aviation safety as well as conservation efforts for grassland birds. In 2011-2013, we investigated effects of vegetation composition and harvest frequency on seasonal species richness and habitat use of grassland birds and modeled relative abundance, aviation risk, and conservation value of birds associated with biomass crops. Avian relative abundance was greater in switchgrass monoculture plots during the winter months, whereas Native Warm-Season Grass (NWSG) mixed species plantings were favored by species during the breeding season. Conversely, treatment differences in aviation risk and conservation value were not biologically significant. Only 2.6% of observations included avian species of high hazard to aircraft, providing support for semi-natural grasslands as a feasible landcover option at civil airports. Additionally, varied harvest frequencies across a mosaic of switchgrass monocultures and NWSG plots allows for biomass production with multiple vegetation structure options for grassland birds to increase seasonal avian biodiversity and habitat use.

Influence of drift and admixture on population structure of American black bears (Ursus americanus) in the Central Interior Highlands, USA, 50 years after translocation.

Bottlenecks, founder events, and genetic drift often result in decreased genetic diversity and increased population differentiation. These events may follow abundance declines due to natural or anthropogenic perturbations, where translocations may be an effective conservation strategy to increase population size. American black bears (Ursus americanus) were nearly extirpated from the Central Interior Highlands, USA by 1920. In an effort to restore bears, 254 individuals were translocated from Minnesota, USA, and Manitoba, Canada, into the Ouachita and Ozark Mountains from 1958 to 1968. Using 15 microsatellites and mitochondrial haplotypes, we observed contemporary genetic diversity and differentiation between the source and supplemented populations. We inferred four genetic clusters: Source, Ouachitas, Ozarks, and a cluster in Missouri where no individuals were translocated. Coalescent models using approximate Bayesian computation identified an admixture model as having the highest posterior probability (0.942) over models where the translocation was unsuccessful or acted as a founder event. Nuclear genetic diversity was highest in the source (AR = 9.11) and significantly lower in the translocated populations (AR = 7.07-7.34; P = 0.004). The Missouri cluster had the lowest genetic diversity (AR = 5.48) and served as a natural experiment showing the utility of translocations to increase genetic diversity following demographic bottlenecks. Differentiation was greater between the two admixed populations than either compared to the source, suggesting that genetic drift acted strongly over the eight generations since the translocation. The Ouachitas and Missouri were previously hypothesized to be remnant lineages. We observed a pretranslocation remnant signature in Missouri but not in the Ouachitas.

Behavioral traits and airport type affect mammal incidents with U.S. civil aircraft.

Wildlife incidents with aircraft cost the United States (U.S.) civil aviation industry >US$1.4 billion in estimated damages and loss of revenue from 1990 to 2009. Although terrestrial mammals represented only 2.3 % of wildlife incidents, damage to aircraft occurred in 59 % of mammal incidents. We examined mammal incidents (excluding bats) at all airports in the Federal Aviation Administration (FAA) National Wildlife Strike Database from 1990 to 2010 to characterize these incidents by airport type: Part-139 certified (certificated) and general aviation (GA). We also calculated relative hazard scores for species most frequently involved in incidents. We found certificated airports had more than twice as many incidents as GA airports. Incidents were most frequent in October (n = 215 of 1,764 total) at certificated airports and November (n = 111 of 741 total) at GA airports. Most (63.2 %) incidents at all airports (n = 1,523) occurred at night but the greatest incident rate occurred at dusk (177.3 incidents/hr). More incidents with damage (n = 1,594) occurred at GA airports (38.6 %) than certificated airports (19.0 %). Artiodactyla (even-toed ungulates) incidents incurred greatest (92.4 %) damage costs (n = 326; US$51.8 million) overall and mule deer (Odocoileus hemionus) was the most hazardous species. Overall, relative hazard score increased with increasing log body mass. Frequency of incidents was influenced by species relative seasonal abundance and behavior. We recommend airport wildlife officials evaluate the risks mammal species pose to aircraft based on the hazard information we provide and consider prioritizing management strategies that emphasize reducing their occurrence on airport property.

Predicting current and future habitat of Indian pangolin (Manis crassicaudata) under climate change.

Climate change is among the greatest drivers of biodiversity loss, threatening up to 15-30% of described species by the end of the twenty-first century. We estimated the current suitable habitat and forecasted future distribution ranges of Indian pangolin (Manis crassicaudata) under climate change scenarios. We collected occurrence records of Indian pangolin using burrow counts, remote camera records and previously published literature in Pakistan during 2021-2023. We downloaded bioclimatic data for current (1970-2000) and future (2041-2060, 2061-2080, 2081-2100) climate scenarios from the WorldClim database using the Hadley Global Environment Model (HadGEM3-GC31-LL). We used MaxEnt software to predict current and future distributions of Indian pangolin, then computed the amount of habitat lost, gained, and unchanged across periods. We obtained 560 Indian pangolin occurrences overall, 175 during the study, and 385 from our literature search. Model accuracy was very good (AUC = 0.885, TSS = 0.695), and jackknife tests of variable importance showed that the contribution of annual mean temperature (bio1) was greatest (33.4%), followed by the mean temperature of the coldest quarter (bio-12, 29.3%), temperature seasonality (bio 4, 25.9%), and precipitation seasonality (bio 15, 11.5%). The maxent model predicted that during the current time period (1970-2000) highly suitable habitat for Indian pangolin was (7270 km2, 2.2%), followed by moderately suitable (12,418 km2, 3.7%), less suitable (49,846 km2, 14.8%), and unsuitable habitat (268,355 km2, 79.4%). Highly suitable habitat decreased in the western part of the study area under most SSPs and in the central parts it declined under all SSPs and in future time periods. The predicted loss in the suitable habitat of the Indian pangolin was greatest (26.97%) under SSP 585 followed by SSP 126 (23.67%) during the time 2061-2080. The gain in suitable habitat of Indian pangolin was less than that of losses on average which ranged between 1.91 and 13.11% under all SSPs during all time periods. While the stable habitat of the Indian pangolin ranged between 64.60 and 83.85% under all SSPs during all time periods. Our study provides the current and future habitat ranges of Indian pangolin in the face of a changing climate. The findings of our study could be helpful for policymakers to set up conservation strategies for Indian pangolin in Pakistan.

Spatial overlap of gray wolves and ungulate prey changes seasonally corresponding to prey migration.

BACKGROUND: Prey are more vulnerable during migration due to decreased familiarity with their surroundings and spatially concentrated movements. Predators may respond to increased prey vulnerability by shifting their ranges to match prey. Moose (Alces alces) and white-tailed deer (Odocoileus virginianus) are primary gray wolf (Canis lupus) prey and important subsistence species for Indigenous communities. We hypothesized wolves would increase use of ungulate migration corridors during migrations and predicted wolf distributions would overlap primary available prey. METHODS: We examined seasonal gray wolf, moose, and white-tailed deer movements on and near the Grand Portage Indian Reservation, Minnesota, USA. We analyzed GPS collar data during 2012-2021 using Brownian bridge movement models (BBMM) in Migration Mapper and mechanistic range shift analysis (MRSA) to estimate individual- and population-level occurrence distributions and determine the status and timing of range shifts. We estimated proportional overlap of wolf distributions with moose and deer distributions and tested for differences among seasons, prey populations, and wolf sex and pack affiliations. RESULTS: We identified a single migration corridor through which white-tailed deer synchronously departed in April and returned in October-November. Gray wolf distributions overlapped the deer migration corridor similarly year-round, but wolves altered within-range distributions seasonally corresponding to prey distributions. Seasonal wolf distributions had the greatest overlap with deer during fall migration (10 October-28 November) and greatest overlap with moose during summer (3 May-9 October). CONCLUSIONS: Gray wolves did not increase their use of the white-tailed deer migration corridor but altered distributions within their territories in response to seasonal prey distributions. Greater overlap of wolves and white-tailed deer in fall may be due to greater predation success facilitated by asynchronous deer migration movements. Greater summer overlap between wolves and moose may be linked to moose calf vulnerability, American beaver (Castor canadensis) co-occurrence, and reduced deer abundance associated with migration. Our results suggest increases in predation pressure on deer in fall and moose in summer, which can inform Indigenous conservation efforts. We observed seasonal plasticity of wolf distributions suggestive of prey switching; that wolves did not exhibit migratory coupling was likely due to spatial constraints resulting from territoriality.

Effects of lethal management on gray wolf pack persistence and reproduction in Wisconsin, USA.

Direct human-caused mortality accounts for about half of all large mammal mortality in North America. For social species like gray wolves (Canis lupus), the death of pack members can disrupt pack structure and cause pack dissolution, and mortality of breeding adults or wolves during reproduction and pup-rearing can decrease pup recruitment. We estimated minimum and maximum probability of wolf pack persistence in Wisconsin, USA, during biological years (15 April-14 April) 2011-2019 and evaluated the influence of pack size and legal harvest mortality on pack persistence during 2012-2014. Harvests comprised 75-161 mortalities within 194 monitored packs during 2012-2014, with 56-74% of packs having no wolves harvested each year. As an index of reproduction during 2013-2019, we also estimated the proportion of packs where pups responded to howl surveys. We evaluated the influence of pack size, legal harvest, and agency removal on reproduction during 2013-2015. Annual maximum pack persistence probability was uniformly high (0.95-1.00), and annual minimum pack persistence probability ranged from 0.86-0.98 with a possible decline during years of harvest. Reproduction was similar in years following harvest and agency removal (2013-2015, pup response = 0.27-0.40), and years without harvest or agency removal the year prior (2016-2019, pup response = 0.28-0.66). Pack size had a positive effect on pack persistence and reproduction. Total number of wolf mortalities and number of adult male and females removed did not influence pack persistence or reproduction. We suggest that low per-pack mortality, timing of harvest and agency removal, and harvest characteristics during 2012-2014 supported stable pack persistence and reproduction.

Occurrence and temporal overlap of sympatric jungle cats and leopard cats in Parsa‒Koshi Complex, Nepal.

Co-occurrence and spatial and temporal overlap of sympatric jungle and leopard cats are influenced by habitat preferences, and interspecific competition. Understanding these factors influence is crucial for developing effective conservation strategies. We conducted a camera survey in Parsa‒Koshi Complex (PKC), Nepal during December 2022-March 2023 to investigate factors influencing occupancy and spatial and temporal overlap between jungle cats (Felis chaus) and leopard cats (Prionailurus bengalensis). The mean detection probability (t = 0.664, p = 0.507) did not differ between jungle cats (p = 0.500 ± 0.289) and leopard cats (p = 0.501 ± 0.288); however, occupancy (t = 31.008, p < 0.001) was greater for jungle cats (ψ = 0.247 ± 0.020) than leopard cats (ψ = 0.178 ± 0.019). Jungle cats and leopard cats were positively associated with large predators, and jungle cats were positively associated with human presence and negatively associated with canopy cover. We observed high diel overlap between leopard cats and jungle cats (Dhat1 = 0.802, norm0CI: 0.720-0.884), with both species largely nocturnal. Co-existence of jungle cats and leopard cats in PKC appears to be facilitated by spatial segregation. These findings provide valuable insights into the complex ecological dynamics and interactions between sympatric jungle and leopard cats.

Uncertainty and precaution in hunting wolves twice in a year: Reanalysis of Treves and Louchouarn.

Management of wolves is controversial in many jurisdictions where wolves live, which underscores the importance of rigor, transparency, and reproducibility when evaluating outcomes of management actions. Treves and Louchouarn 2022 (hereafter TL) predicted outcomes for various fall 2021 hunting scenarios following Wisconsin's judicially mandated hunting and trapping season in spring 2021, and concluded that even a zero harvest scenario could result in the wolf population declining below the population goal of 350 wolves specified in the 1999 Wisconsin wolf management plan. TL further concluded that with a fall harvest of > 16 wolves there was a "better than average possibility" that the wolf population size would decline below that 350-wolf threshold. We show that these conclusions are incorrect and that they resulted from mathematical errors and selected parameterizations that were consistently biased in the direction that maximized mortality and minimized reproduction (i.e., positively biased adult mortality, negatively biased pup survival, further halving pup survival to November, negatively biased number of breeding packs, and counting harvested wolves twice among the dead). These errors systematically exaggerated declines in predicted population size and resulted in erroneous conclusions that were not based on the best available or unbiased science. Corrected mathematical calculations and more rigorous parameterization resulted in predicted outcomes for the zero harvest scenario that more closely coincided with the empirical population estimates in 2022 following a judicially prevented fall hunt in 2021. Only in scenarios with simulated harvest of 300 or more wolves did probability of crossing the 350-wolf population threshold exceed zero. TL suggested that proponents of some policy positions bear a greater burden of proof than proponents of other positions to show that "their estimates are accurate, precise, and reproducible". In their analysis, TL failed to meet this standard that they demanded of others.

Factors influencing mammalian community occupancy in Dhorpatan Hunting Reserve, Nepal.

The decline in mammalian species diversity is increasing worldwide, with areas characterized by high human activities experiencing more prominent effects. Knowledge of spatial distributions of species and factors acting on them is necessary for effective management. We evaluated community-level occupancy of mammal species in Dhorpatan Hunting Reserve (DHR), Nepal using remote cameras during 15 March-15 June 2022. We used mammal species detections from remote cameras and multispecies hierarchical occupancy modeling to assess the effects of environmental and anthropogenic variables on community-level occupancy of detected mammal species. We identified a highly heterogeneous mammal species community at DHR with greatest detection probability (0.21) for red fox (Vulpes vulpes) and lowest (0.08) for blue sheep (Pseudois nayaur). Naïve occupancy ranged from 0.31 for giant-flying squirrel (Petaurista magnificus) to 0.84 for red fox. Mammal community occupancy increased with increasing canopy cover and number of livestock detections, but overall occupancy declined close to human settlements. The findings of this study can be used for developing policy at DHR for the management of mammal species through reducing the potential increase of human settlements or livestock grazing.

Factors influencing frequency and severity of human-American black bear conflicts in New York, USA.

Free-ranging large carnivores are involved in human-wildlife conflicts which can result in economic costs. Understanding factors that lead to human-wildlife conflicts is important to mitigate these negative effects and facilitate human-carnivore coexistence. We used a human-American black bear (Ursus americanus) conflict database maintained by the New York State Department of Environmental Conservation to determine whether drought, conflicts within the Adirondack and Catskill Parks as compared to outside of these parks, mild severity (Class 3) conflicts early in the year (April-June), and bear harvest in the previous year (as an index of bear abundance), were associated with greater frequency of high or moderate severity (Class 1-2) conflicts later in the year (July-September) across New York, USA. During 2006-2019, we obtained 3,782 mild severity conflict records early in the year, and 1,042 high or moderate severity records later in the year. We found that a one standard deviation increase in the cumulative precipitation difference from mean early in the year (about 7.59 cm) coincided with a 20% decrease in conflicts, and that Wildlife Management Units (WMUs) within the parks were predicted to have 5.61 times as many high or moderate severity conflicts as WMUs outside the parks. We also found that a one standard deviation increase in the frequency of mild severity conflicts (equivalent to 5.68 conflicts) early in the year coincided with an increase in the frequency of high or moderate severity conflicts in a WMU later in the year by 49%, while a one standard deviation increase in the bear abundance index in the previous year (0.14 bears/10 km2) coincided with a 23% increase in high or moderate severity conflicts. To reduce the frequency and severity of conflicts to facilitate human-black bear coexistence, we recommend the following measures to be taken in place consistently and build over time in local communities: (i) further reducing black bear access to anthropogenic foods and other attractants, (ii) non-lethal measures including bear-resistant waste management, (iii) electric fencing, and (iv) modifying placement or configuration of field crops.

A review of caracal and jungle cat diets across their geographical ranges during 1842-2021.

Co-occurring carnivore species that are phylogenetically related or of similar size, morphology, and ecological needs often reduce competition by partitioning shared resources through temporal, spatial, and dietary niche segregation via behavioral adaptations. Caracals (Caracal caracal) and jungle cats (Felis chaus) co-occur in portions of their geographical ranges and are expected to display resource segregation in these ranges. We compiled scat, stomach content, and prey remains found data from published and unpublished sources to summarize information on the diets of caracals and jungle cats across their geographical ranges during 1842-2021. We obtained 63 sources from 26 countries in Europe, Asia, and Africa, in which caracal diet included 151 species while jungle cat diet included 61 species. We found that caracals and jungle cats did not exhibit dietary niche partitioning and had greater dietary similarities in areas of range overlap. We also found that caracals consumed more diverse prey species including prey with greater average body mass compared to jungle cats. Our results suggest that greater prey diversity in areas of range overlap, caracal predation on wide range of prey, and opportunistic feeding behavior that facilitates consumption of more diverse prey species compared to jungle cats, may facilitate co-occurrence between these two felid species.