USGS ScienceSearch

USGS · 70251770

Joint spatial modeling bridges the gap between disparate disease surveillance and population monitoring efforts informing conservation of at-risk bat species

Abstract

White-Nose Syndrome (WNS) is a wildlife disease that has decimated hibernating bats since its introduction in North America in 2006. As the disease spreads westward, assessing the potentially differential impact of the disease on western bat species is an urgent conservation need. The statistical challenge is that the disease surveillance and species response monitoring data are not co-located, available at different spatial resolutions, non-Gaussian, and subject to observation error requiring a novel extension to spatially misaligned regression models for analysis. Previous work motivated by epidemiology applications has proposed two-step approaches that overcome the spatial misalignment while intentionally preventing the human health outcome from informing estimation of exposure. In our application, the impacted animals contribute to spreading the fungus that causes WNS, motivating development of a joint framework that exploits the known biological relationship. We introduce a Bayesian, joint spatial modeling framework that provides inferences about the impact of WNS on measures of relative bat activity and accounts for the uncertainty in estimation of WNS presence at non-surveyed locations. Our simulations demonstrate that the joint model produced more precise estimates of disease occurrence and unbiased estimates of the association between disease presence and the count response relative to competing two-step approaches. Our statistical framework provides a solution that leverages disparate monitoring activities and informs species conservation across large landscapes. Stan code and documentation are provided to facilitate access and adaptation for other wildlife disease applications.

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Christian Stratton, Kathryn Irvine, Katharine M. Banner, Emily S. Almberg, Daniel Bachen, Kristina Smucker. 2024-02-24. Joint spatial modeling bridges the gap between disparate disease surveillance and population monitoring efforts informing conservation of at-risk bat species. https://doi.org/10.1007/s13253-023-00593-8

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

Estimating submersed aquatic vegetation biomass using a hierarchical Bayesian model of ordinal measurements and structural zeros

This study explores the prediction of a continuous target variable that is primarily measured as an ordinal outcome. Interest in such prediction may occur given data from double- or two-phase sampling designs, where an auxiliary variable is measured at all sampling units and a target and predictor variable, X , at a subset of those units. This study focuses on the case where X is nonnegative, is conditional on a zero process, and is sampled using a cluster sampling design. We model the ordinal outcome using cumulative logistic regression and the conditional regressor X , following log transformation, in a Bayesian setting. Estimation of model parameters and prediction of missing X is generally accurate and precise, particularly with sufficiently large Pr(X > 0) and associated variances. This study is motivated by a need to accurately estimate the biomass of submersed aquatic vegetation on the Upper Mississippi River, where two types of measurements for biomass are available: inexpensive, ordinal biomass scores and expensive, continuous diver-harvested biomass data. Here, X represents plant biomass and Pr(X > 0) is the usual species site occupancy measure. When fit using biomass data from the submerged aquatic vegetation species Vallisneria americana Michx, our model estimates parameters within the parameter space region where the model performed acceptably using synthetic data. We expect this model to appeal to investigators with ordered outcomes, cluster designs, and one or more continuous, positive predictors.

Journal of Agricultural, Biological and Environmen

Nonstationary demographic state-space models using unreplicated counts for species undergoing environmental stressors

A fundamental task in ecological statistics is to estimate abundance and growth rate distributions from wildlife monitoring data to inform conservation management. Modeling time series of wildlife populations presents a number of challenges from both statistical and ecological perspectives, including discreteness; lack of replication; nonstationarity; and observation, demographic, and other phenomenological processes. Nonstationary dynamics are often exhibited by populations undergoing environmental stressors. Models must account for these characteristics to produce reliable estimates of abundance and trends, yet estimation can be challenging with unreplicated data. We propose nonstationary demographic state-space models using unreplicated counts for populations undergoing environmental stressors. A reduced growth rate model matches the complexity of the unreplicated count data, and a fecundity bound on growth rate distributions allows the separation of processes affecting growth rates like environmental stressors from those affecting abundance external to growth rates like migration. NDSSMs allow for the embedding of nonstationary model components, and we explore the use of changepoints, volatility clustering, and migration processes. We apply the proposed nonstationary models in case studies of herons affected by predator/competitor reestablishment and three bat species affected by a fungal pathogen causing white-nose syndrome. Nonstationary models outperform stationary models and generalized linear mixed effects models according to model scoring and visual inspection of predictions, and provide estimates more consistent with published values. Incorporating migration improves model fit universally, even with approximate one-way immigration, most likely because populations are extirpated, recolonized, and increase multiple-fold over the upper bound set by species fecundity. In addition, estimates of the timing and severity of the environmental stressor differed for models with migration. Including nonstationary and demographic components in a fecundity-bounded growth rate model improves inference and benefits interpretability of hyperparameters. In turn, this adjusts uncertainties in predictions of abundance and growth rates over time, providing the ingredients needed for informed conservation analysis and for directing future monitoring of at-risk species.

Journal of Agricultural, Biological and Environmen

Bayesian approaches to proxy uncertainty quantification in paleoecology: A mathematical justification and practical integration

Paleoenvironmental data are essential for reconstructing environmental conditions in the distant past, and these reconstructions strongly depend on proxies and age–depth models. Proxies are indirect measurements that substitute for variables that cannot be directly measured, such as past precipitation. Conversely, an age–depth model is a tool that correlates the observed proxy with a specific moment in time. Bayesian age–depth modelling has proved to be a powerful method for estimating sediment ages and their associated uncertainties. However, there remains considerable potential for further integration into proxy analysis. In this paper, we explore a mathematical justification and a computational approach that integrates uncertainty at the age–depth level and propagates it to the proxy scale in the form of a posterior predictive distribution. This method mitigates potential biases and errors by removing the need to assign a single age to a given proxy measurement. It allows for quantifying the likelihood that proxy data values correspond to modelled ages, thus enabling the quantification of uncertainty in both the temporal and proxy value domains. The use of Bayesian statistics in proxy analysis represents a relatively recent advancement. We aim to mathematically justify incorporating the Markov chain Monte Carlo output from age–depth models into proxy analysis and to present a novel methodology for constructing environmental reconstructions using this approach.

Journal of Agricultural, Biological and Environmen