Integrated community occupancy models: A framework to assess occurrence and biodiversity dynamics using multiple data sources

The occurrence and distributions of wildlife populations and communities are shifting as a result of global changes. To evaluate whether these shifts are negatively impacting biodiversity processes, it is critical to monitor the status, trends, and effects of environmental variables on entire communities. However, modeling the dynamics of multiple species simultaneously can require large amounts of diverse data, and few modeling approaches exist to simultaneously provide species and community level inferences. We present an integrated community occupancy model (ICOM) that unites principles of data integration and hierarchical community modeling in a single framework to provide inferences on species-specific and community occurrence dynamics using multiple data sources. We use simulations to compare the ICOM to previously developed hierarchical community occupancy models and single species integrated distribution models. We then apply our model to assess the occurrence and biodiversity dynamics of foliage-gleaning birds in the White Mountain National Forest in the northeastern USA from 2010-2018 using three independent data sources. Simulations reveal that integrating multiple data sources in the ICOM increased precision and accuracy of species and community level inferences compared to single data source models, although benefits of integration were dependent on data source quality (e.g., amount of replication). Compared to single species models, the ICOM yielded more precise species-level estimates. Within our case study, the ICOM had the highest out-of-sample predictive performance compared to single species models and models that used only a subset of the three data sources. The ICOM offers an attractive approach to estimate species and biodiversity dynamics, which is additionally valuable to inform management objectives of both individual species and their broader communities.

Integrating automated acoustic vocalization data and point count surveys for estimation of bird abundance

Monitoring wildlife abundance across space and time is an essential task to study their population dynamics and inform effective management. Acoustic recording units are a promising technology for efficiently monitoring bird populations and communities. We present an integrated modeling framework that combines high-quality but temporally sparse bird point count survey data with acoustic recordings. Using simulations, we compare the accuracy and precision of abundance estimates using differing amounts of acoustic vocalizations obtained from a clustering algorithm, point count data, and a subset of manually validated acoustic vocalizations. We also use our modeling framework in a case study to estimate abundance of the Eastern Wood-Pewee (Contopus virens) in Vermont, U.S.A. The simulation study reveals that combining acoustic and point count data via an integrated model improves accuracy and precision of abundance estimates compared with models informed by either acoustic or point count data alone. Combining acoustic data with only a small number of point count surveys yields estimates of abundance without the need for validating any of the identified vocalizations from the acoustic data. Within our case study, the integrated models provided moderate support for a decline of the Eastern Wood-Pewee in this region. Our integrated modeling approach combines dense acoustic data with few point count surveys to deliver reliable estimates of species abundance without the need for manual identification of acoustic vocalizations or a prohibitively expensive large number of repeated point count surveys. Our proposed approach offers an efficient monitoring alternative for large spatio-temporal regions when point count data are difficult to obtain or when monitoring is focused on rare species with low detection probability.

Trends in bird abundance differ among protected forests but not bird guilds

Improved monitoring and associated inferential tools to efficiently identify declining bird populations, particularly of rare or sparsely distributed species, is key to informed conservation and management across large spatio-temporal regions. We assess abundance trends for 106 bird species in a network of eight national park forests located within the northeast USA from 2006-2019 using a novel hierarchical model. We develop a multi-species, multi-region removal sampling model that shares information across species and parks to enable inference on rare species and sparsely sampled parks and to evaluate the effects of local forest structure. Trends in bird abundance over time varied widely across parks, but species showed similar trends within parks. Three parks (Acadia, Marsh-Billings-Rockefeller, and Morristown) decreased in bird abundance across all species, while three parks (Saratoga, Roosevelt-Vanderbilt, and Weir-Farm) increased in abundance. Bird abundance peaked at medium levels of basal area and high levels of percent forest and forest regeneration, with percent forest having the largest effect. Variation in these effects across parks could be a result of differences in forest structural stage and diversity. Our novel hierarchical model enables estimates of abundance at the network, park, guild, and species levels. We found large variation in abundance trends across parks but not across bird guilds, suggesting that local forest condition may have a broad and consistent effect on the entire bird community within a given park. Management should target the three parks with overall decreasing trends in bird abundance to further identify what specific factors are driving observed declines across the bird community. Understanding how bird communities respond to local forest structure and other stressors is crucial for informed and lasting management.

ForestFit : An R package for modeling tree diameter distributions

Modeling the diameter distribution of trees in forest stands is a common forestry task that supports key biologically and economically relevant management decisions. The choice of model used to represent the diameter distribution and how to estimate its parameters has received much attention in the forestry literature; however, accessible software that facilitates comprehensive comparison of the myriad modeling approaches is not available. To this end, we developed an R package called ForestFit that simplifies estimation of common probability distributions used to model tree diameter distributions, including the two- and three-parameter Weibull distributions, Johnsons SB distribution, Birnbaum-Saunders distribution, and finite mixture distributions. Frequentist and Bayesian techniques are provided for individual tree diameter data, as well as grouped data. Additional functionality facilitates fitting growth curves to height-diameter data. The package also provides a set of functions for computing probability distributions and simulating random realizations from common finite mixture models.

Assessing soundscape disturbance through hierarchical models and acoustic indices: a case study on a shelterwood logged northern Michigan forest

Assessing the effects of anthropogenic disturbances on wildlife is a necessary conservation task. The soundscape is a critical habitat component for acoustically communicating organisms, but the use of the soundscape as a tool for assessing disturbance impacts has been relatively unexplored until recently. Here we present a broad modeling framework for assessing disturbance impacts on soundscapes, which we apply to quantify the influence of a shelterwood logging on soundscapes in northern Michigan. Our modeling approach can be broadly applied to assess anthropogenic disturbance impacts on soundscapes. The approach accommodates inherent differences in control and treatment sites to improve inference about treatment effects, while also accounting for extraneous variables (e.g., rain) that influence acoustic indices. Recordings were obtained at 13 sites before and after a shelterwood logging. Four sites were in the logging region and nine sites served as control recordings outside the logging region. We quantify the soundscapes using common acoustic indices (Normalized Difference Soundscape Index (NDSI), Acoustic Entropy (H), Acoustic Complexity Index (ACI), Acoustic Evenness Index (AEI), Welch Power Spectral Density (PSD)) and build two hierarchical Bayesian models to quantify the changes in the soundscape over the study period. Our analysis reveals no long-lasting effects of the shelterwood logging on the soundscape diversity as measured by the NDSI, but analysis of H, AEI, and PSD suggest changes in the evenness of sounds across the frequency spectrum, indicating a potential shift in the avian species communicating in the soundscapes as a result of the logging. Acoustic recordings, in conjunction with this modeling framework, can deliver cost efficient assessment of disturbance impacts on the landscape and underlying biodiversity as represented through the soundscape.

Characterizing functional relationships between anthropogenic and biological sounds: A western New York state soundscape case study

Roads are a widespread feature of landscapes worldwide, and road traffic sound potentially makes nearby habitat unsuitable for acoustically communicating organisms. It is important to understand the influence of roads at the soundscape level to mitigate negative impacts of road sound on individual species as well as subsequent effects on the surrounding landscape. We seek to characterize the relationship between anthropogenic and biological sounds in western New York and assess the extent to which available traffic data explains variability in anthropogenic noise. Recordings were obtained in the spring of 2016 at 18 sites throughout western New York. We used the Welch Power Spectral Density (PSD) at low frequencies (0.5-2 kHz) to represent anthropogenic noise and PSD values at higher frequencies (2-11 kHz) to represent biological sound. Relationships were modeled using a novel two-stage hierarchical Bayesian model utilizing beta regression and basis splines. Model results and map predictions illustrate that anthropogenic noise and biological sound have an inverse relationship, and anthropogenic noise is greatest in close proximity to high traffic volume roads. The predictions have large uncertainty, resulting from the temporal coarseness of public road data used as a proxy for traffic sound. Results suggest that finer temporal resolution traffic sound data, such as crowd-sourced time-indexed traffic data from geographic positioning systems, might better account for observed temporal changes in the soundscape. The use of such data, in combination with the proposed modeling framework, could have important implications for the development of sound management policies.