Improving inference for aerial surveys of bears: The importance of assumptions and the cost of unnecessary complexity

Accounting for missing ticks: Use (or lack thereof) of hierarchical models in tick ecology studies

Ixodid (hard) ticks play important ecosystem roles and have significant impacts on animal and human health via tick-borne diseases and physiological stress from parasitism. Tick occurrence, abundance, activity, and key life-history traits are highly influenced by host availability, weather, microclimate, and landscape features. As such, changes in the environment can have profound impacts on ticks, their hosts, and the spread of diseases. Researchers recognize that spatial and temporal factors influence activity and abundance and attempt to account for both by conducting replicate sampling bouts spread over the tick questing period. However, common field methods notoriously underestimate abundance, and it is unclear how (or if) tick studies model the confounding effects of factors influencing activity and abundance. This step is critical as unaccounted variance in detection can lead to biased estimates of occurrence and abundance. We performed a descriptive review to evaluate the extent to which studies account for the detection process while modeling tick data. We also categorized the types of analyses that are commonly used to model tick data. We used hierarchical models (HMs) that account for imperfect detection to analyze simulated and empirical tick data, demonstrating that inference is muddled when detection probability is not accounted for in the modeling process. Our review indicates that only 5 of 412 (1 %) papers explicitly accounted for imperfect detection while modeling ticks. By comparing HMs with the most common approaches used for modeling tick data (e.g., ANOVA), we show that population estimates are biased low for simulated and empirical data when using non-HMs, and that confounding occurs due to not explicitly modeling factors that influenced both detection and abundance. Our review and analysis of simulated and empirical data shows that it is important to account for our ability to detect ticks using field methods with imperfect detection. Not doing so leads to biased estimates of occurrence and abundance which could complicate our understanding of parasite-host relationships and the spread of tick-borne diseases. We highlight the resources available for learning HM approaches and applying them to analyzing tick data.

Distance sampling of duikers in the rainforest: Dealing with transect avoidance

Bushmeat is a major source of protein and income in tropical regions but is often over-harvested. A better monitoring of bushmeat stocks could help achieve sustainability. We used a combination of simulations and transect survey data collected from blue duikers (Philantomba monticola) in the Lomako wildlife reserve, Democratic Republic of the Congo, to investigate the use of transect-based distance sampling to monitor bushmeat stocks. The comparison of dung piles and direct observations of duikers evidenced that animals avoided both the transects in the absence of observers, and the observers themselves. This type of behavioural response appeared common in a literature survey. It causes a negative bias in the estimates of population densities from the standard distance sampling methodology. This negative bias would lead to over-pessimistic predictions of population viability, especially if the behavioural response is more intense in the locations where the animals are hunted. In turn, this would lead to excessively conservative management recommendations. To correct for the effect of the behavioural response of the animals to either the transects or the observers, we recommend recording both the forward and perpendicular distances to the observers (2D distance sampling), not just the perpendicular distance. We also recommend multiple-observer protocols. As a cautionary note, we also demonstrate a scenario where the intensity of the behavioural response is too high to reliably estimate the abundance of the population. As a perspective, we outline the general principles of a local stakeholder-based program combining distance sampling with less intensive types of ecological indicators to monitor wildlife populations.

Steep and deep: Terrain and climate factors explain brown bear (Ursus arctos) alpine den site selection to guide heli-skiing management

Winter recreation and tourism continue to expand worldwide, and where these activities overlap with valuable wildlife habitat, there is greater potential for conservation concerns. Wildlife populations can be particularly vulnerable to disturbance in alpine habitats as helicopters and snowmachines are increasingly used to access remote backcountry terrain. Brown bears (Ursus arctos) have adapted hibernation strategies to survive this period when resources and energy reserves are limited, and disturbance could negatively impact fitness and survival. To help identify areas of potential conflict between helicopter skiing and denning brown bears in Alaska, we developed a model to predict alpine denning habitat and an associated data-based framework for mitigating disturbance activities. Following den emergence in spring, we conducted three annual aerial surveys (2015-2017) and used locations from three GPS-collared bears (2008-2014) to identify 89 brown bear dens above the forest line. We evaluated brown bear den site selection of land cover, terrain, and climate factors using resource selection function (RSF) models. Our top model supported the hypothesis that bears selected dens based on terrain and climate factors that maximized thermal efficiency. Brown bears selected den sites characterized by steep slopes at moderate elevations in smooth, well-drained topographies that promoted vegetation and deep snow. We used the RSF model to map relative probability of den selection and found 85% of dens occurred within terrain predicted as prime denning habitat. Brown bear exposure to helicopter disturbance was evident as moderate to high intensities of helicopter flight tracking data overlapped prime denning habitat, and we quantified where the risk of these impact was greatest. We also documented evidence of late season den abandonment due to disturbance from helicopter skiing. The results from this study provide valuable insights into bear denning habitat requirements in subalpine and alpine landscapes. Our quantitative framework can be used to support conservation planning for winter recreation industries operating in habitats occupied by denning brown bears.

Rejection of Schmidt et al.'s estimators for bear population size

Aerial distance sampling of bears to estimate population size has been used throughout many parts of Alaska. The distance sampling models are complex since they need to account for undetected bears and differences in detection probabilities. This will require covariates and mark-recapture data. The models proposed by Schmidt et al. do not use covariates or mark-recapture data and are inappropriate for these surveys.

Improving inference for aerial surveys of bears: The importance of assumptions and the cost of unnecessary complexity

Obtaining useful estimates of wildlife abundance or density requires thoughtful attention to potential sources of bias and precision, and it is widely understood that addressing incomplete detection is critical to appropriate inference. When the underlying assumptions of sampling approaches are violated, both increased bias and reduced precision of the population estimator may result. Bear (Ursus spp.) populations can be difficult to sample and are often monitored using mark‐recapture distance sampling (MRDS) methods, although obtaining adequate sample sizes can be cost prohibitive. With the goal of improving inference, we examined the underlying methodological assumptions and estimator efficiency of three datasets collected under an MRDS protocol designed specifically for bears. We analyzed these data using MRDS, conventional distance sampling (CDS), and open‐distance sampling approaches to evaluate the apparent bias‐precision tradeoff relative to the assumptions inherent under each approach. We also evaluated the incorporation of informative priors on detection parameters within a Bayesian context. We found that the CDS estimator had low apparent bias and was more efficient than the more complex MRDS estimator. When combined with informative priors on the detection process, precision was increased by >50% compared to the MRDS approach with little apparent bias. In addition, open‐distance sampling models revealed a serious violation of the assumption that all bears were available to be sampled. Inference is directly related to the underlying assumptions of the survey design and the analytical tools employed. We show that for aerial surveys of bears, avoidance of unnecessary model complexity, use of prior information, and the application of open population models can be used to greatly improve estimator performance and simplify field protocols. Although we focused on distance sampling‐based aerial surveys for bears, the general concepts we addressed apply to a variety of wildlife survey contexts.