skip to main content

Title: A hierarchical N‐mixture model to estimate behavioral variation and a case study of Neotropical birds

Understanding how and why animals use the environments where they occur is both foundational to behavioral ecology and essential to identify critical habitats for species conservation. However, some behaviors are more difficult to observe than others, which can bias analyses of raw observational data. To our knowledge, no method currently exists to model how animals use different environments while accounting for imperfect behavior‐specific detection probability. We developed an extension of a binomial N‐mixture model (hereafter the behavior N‐mixture model) to estimate the probability of a given behavior occurring in a particular environment while accounting for imperfect detection. We then conducted a simulation to validate the model's ability to estimate the effects of environmental covariates on the probabilities of individuals performing different behaviors. We compared our model to a naïve model that does not account for imperfect detection, as well as a traditional N‐mixture model. Finally, we applied the model to a bird observation data set in northwest Costa Rica to quantify how three species behave in forests and farms. Simulations and sensitivity analyses demonstrated that the behavior N‐mixture model produced unbiased estimates of behaviors and their relationships with predictor variables (e.g., forest cover, habitat type). Importantly, the behavior N‐mixture model accurately characterized uncertainty, unlike the naïve model, which often suggested erroneous effects of covariates on behaviors. When applied to field data, the behavior N‐mixture model suggested that Hoffmann's woodpecker (Melanerpes hoffmanii) and Inca dove (Columbina inca) behaved differently in forested versus agricultural habitats, while turquoise‐browed motmot (Eumomota superciliosa) did not. Thus, the behavior N‐mixture model can help identify habitats that are essential to a species' life cycle (e.g., where individuals nest, forage) that nonbehavioral models would miss. Our model can greatly improve the appropriate use of behavioral survey data and conclusions drawn from them. In doing so, it provides a valuable path forward for assessing the conservation value of alternative habitat types.

more » « less
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Ecological Applications
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Advances in biologging have increased the understanding of how animals interact with their environment, especially for cryptic species. For example, giant armadillos (Priodontes maximus) are the largest extant species of armadillo but are rarely encountered due to their fossorial and nocturnal behavior. Through the analysis of speed, turning angles, and accelerometer activity counts, we estimated behavioral states, characterized activity budgets, and investigated the state-habitat associations exhibited by individuals monitored with GPS telemetry in the Brazilian Pantanal from 2019 to 2020. This methodology is proposed as a useful framework for the identification of priority habitat. Using the non-parametric Bayesian mixture model for movement (M3), we estimated four latent behavioral states that were named ‘vigilance-excavation’, ‘local search’, ‘exploratory’, and ‘transit’. These states appeared to correspond with behavior near burrows or termite mounds, foraging, ranging, and rapid movements, respectively. The first and last hours of activity presented relatively high proportions of the vigilance-excavation state, while most of the activity period was dominated by local search and exploratory states. The vigilance-excavation state occurred more frequently in regions between forest and closed savannas, whereas local search was more likely in high proportions of closed savanna. Exploratory behavior probability increased in areas with high proportions of both forest and closed savanna. Our results establish a baseline for behavioral complexity, activity budgets, and habitat associations in a relatively pristine environment that can be used for future work to investigate anthropogenic impacts on giant armadillo behavior and fitness. The integration of accelerometer and GPS-derived movement data through our mixture model has the potential to become a powerful methodological approach for the conservation of other cryptic species.

    more » « less
  2. Abstract

    Accurate estimates of survival are crucial for many management decisions in translocation programs. Maximizing detection probabilities and reducing sampling biases for released animals can aid in estimates of survival. One important source of sampling bias is an animal’s behavior. For example, individuals that are consistently more exploratory or active may be more likely to be detected visually. Behavioral traits can be related to survival after reintroduction, and because many pre‐release treatments aim to manipulate animal behavior, it is critical to tease apart relationships between behavior and detection probability. Here, we assessed the repeatability (intra‐individual consistency and inter‐individual variation) of behavioral traits for an endangered amphibian, the mountain yellow‐legged frog (Rana muscosa). Because new technological tools offer one potential solution for reducing sampling biases while increasing detection, we also tested whether a long‐range passive integrated transponder (PIT) tag reader could enhance surveys for these individuals after translocation into the wild. After confirming thatex situbredR. muscosaexhibit repeatable behavioral traits (repeatability = 0.25–0.41) and releasing these frogs (N = 196) into the wild, we conducted post‐release surveys visually and with the long‐range PIT tag reader. Integrating the long‐range reader into surveys improved detection probability four‐fold in comparison to visual surveys alone (~0.09 to ~0.36). Moreover, mark–recapture modeling revealed that tag reader detection probability was not biased toward detecting individuals of specific behavioral types, while visual detection was significantly related to behavioral traits. These results will enable a more accurate understanding of individual differences in post‐release success in translocations. This may be particularly important for amphibian species, which can be difficult to detect and are expected to increasingly be involved in human‐managed breeding and translocation programs due to their vulnerable conservation status.

    more » « less
  3. Powell, Roger (Ed.)
    Abstract Island biogeography theory (IBT) explains and estimates large-scale ecological patterns among islands and isolated habitat patches. Specifically, IBT predicts that the number of species per habitat patch differs as a function of area and isolation as a result of local colonization and extinction. Accurate estimates of species richness are essential for testing predictions of IBT, but differences in detectability of species can lead to bias in empirical data. Hierarchical community models correct for imperfect detection by leveraging information from across the community to estimate species-specific occupancy and detection probabilities. Using the fragmented Ozark glades as our model system, we constructed a hierarchical community model to 1) estimate site-level and regional species richness of small mammals while correcting for detection error, and 2) determine environmental covariates driving occupancy. We sampled 16 glades in southwestern Missouri in summer 2016–2017 and quantified mammal community structure within the glade network. The detected species pool included eight species, and the model yielded a regional species estimate of 8.6 species, with a mean of 3.47 species per glade. Species richness increased with patch area but not isolation, and effects of patch shape varied between species in the community. 
    more » « less
  4. Abstract

    Animals that are successful in urban habitats often have reduced antipredator responses toward people (sometimes called “fear” responses). However, few studies test whether sympatric species differ in their responses to humans, which may explain differing sensitivities to urbanization. Here, we quantified the behavioral and physiological responses to humans in two lizard species, side-blotched lizards (Uta stansburiana) and western fence lizards (Sceloporus occidentalis), across three different habitat types that vary in human impact: natural habitats with low levels of human activity, natural habitats with high levels of human activity, and urban habitats. We found that side-blotched lizards had longer flight initiation distances, were found closer to a refuge, and were more likely to hide than fence lizards, behaviors that could indicate greater fearfulness. Both lizard species were found closer to a refuge and were also more likely to hide in the urban habitat than in the natural habitat with low human impact, which could represent adaptive behaviors for increased risks in urban areas (e.g. cats). Western fence lizards exhibited lower body sizes and conditions in the habitats with moderate and high levels of human activity, whereas these traits did not differ among habitats in side-blotched lizards. Baseline and stress-induced corticosterone concentrations did not differ across habitats for both species, suggesting that human-impacted habitats were not stressful or that lizards had undergone habituation-like processes in these habitats. Taken together, our results highlight the importance of standardized measurements across multiple species in the same habitats to understand differential responses to human-induced environmental change.

    more » « less
  5. Abstract

    Vertical movements can expose individuals to rapid changes in physical and trophic environments—for aquatic fauna, dive profiles from biotelemetry data can be used to quantify and categorize vertical movements. Inferences on classes of vertical movement profiles typically rely on subjective summaries of parameters or statistical clustering techniques that utilize Euclidean matching of vertical movement profiles with vertical observation points. These approaches are prone to subjectivity, error, and bias. We used machine learning approaches on a large dataset of vertical time series (N = 28,217 dives) for 31 post‐nesting leatherback turtles (Dermochelys coriacea). We applied dynamic time warp (DTW) clustering to group vertical movement (dive) time series by their metrics (depth and duration) into an optimal number of clusters. We then identified environmental covariates associated with each cluster using a generalized additive mixed‐effects model (GAMM). A convolutional neural network (CNN) model, trained on standard dive shape types from the literature, was used to classify dives within each DTW cluster by their shape. Two clusters were identified with the DTW approach—these varied in their spatial and temporal distributions, with dependence on environmental covariates, sea surface temperature, bathymetry, sea surface height anomaly, and time‐lagged surface chlorophyllaconcentrations. CNN classification accuracy of the five standard dive profiles was 95%. Subsequent analyses revealed that the two clusters differed in their composition of standard dive shapes, with each cluster dominated by shapes indicative of distinct behaviors (pelagic foraging and exploration, respectively). The use of these two machine learning approaches allowed for discrete behaviors to be identified from vertical time series data, first by clustering vertical movements by their movement metrics (DTW) and second by classifying dive profiles within each cluster by their shapes (CNN). Statistical inference for the identified clusters found distinct relationships with environmental covariates, supporting hypotheses of vertical niche switching and vertically structured foraging behavior. This approach could be similarly applied to the time series of other animals utilizing the vertical dimension in their movements, including aerial, arboreal, and other aquatic species, to efficiently identify different movement behaviors and inform habitat models.

    more » « less