Spatial gaussian processes improve multi-species occupancy models when range boundaries are uncertain and non-overlapping

Wright, Wilson; Irvine, Kathryn; Rodhouse, Thomas; Litt, Andrea

1. Species distribution models enable practitioners to analyze large datasets of encounter records and make predictions about species occurrence at unsurveyed locations. In omnibus surveys that record data on multiple species simultaneously, species ranges are often non-overlapping and misaligned with the administrative unit defining the spatial domain of interest (e.g., a state or province). Consequently, some species display differentially restricted extents within a study area and assuming hard boundaries based on expert opinion or published range maps to restrict their occurrence predictions implies a false sense of certainty in model-based inferences. 2. We propose a multi-species occupancy model with a spatial Gaussian process on site specific effects for each species as a model-based solution. Specifying informative Bayesian hyperpriors on the spatial hyperparameters encapsulates broad-scale correlation among site occupancy probabilities for each species. We fit this model to acoustic detection/non-detection data collected with autonomous recording units during summer of 2016 to 2019 throughout Oregon and Washington, USA on 15 bat species. 3. We found vast improvements in spatial predictions of spotted bat (Euderma maculatum), canyon bat (Parastrellus hesperus), and Brazilian free-tailed bat (Tadarida brasiliensis) when the available environmental predictors were insufficient for characterizing their restricted ranges within the region. 4. In contrast, widespread species (Lasionycteris noctivagans, Myotis californicus, Myotis evotis, Myotis volans) were appropriately modeled using only environmental predictors, such as percentage forest cover and cliffs and canyon cover. 5. Utilizing spatial Gaussian processes within a community or multi-species model incorporates uncertainty in range boundaries and allows for simultaneous predictions for the entire faunal assemblage even if species have non-overlapping or restricted ranges within a spatial domain of interest. Such modelling improvements are essential if species distribution models are to accurately inform monitoring, species recovery plans, and other conservation efforts.