A multistate dynamic site occupancy model for spatially aggregated sessile communities

TL;DR

This paper introduces a multistate dynamic site occupancy model that accounts for local community structure and resampling errors, improving estimates of transition probabilities in sessile communities.

Contribution

It develops a novel nonparametric kernel smoothing approach to incorporate local community structure into occupancy models, addressing biases caused by resampling errors.

Findings

Resampling errors can bias transition probability estimates.

Accounting for local community structure alters inferred community dynamics.

The model can handle anisotropic spatial correlations.

Abstract

Markov community models have been applied to sessile organisms because such models facilitate estimation of transition probabilities by tracking species occupancy at many fixed observation points over multiple periods of time. Estimation of transition probabilities of sessile communities seems easy in principle but may still be difficult in practice because resampling error (i.e., a failure to resample exactly the same location at fixed points) may cause significant estimation bias. Previous studies have developed novel analytical methods to correct for this estimation bias. However, they did not consider the local structure of community composition induced by the aggregated distribution of organisms that is typically observed in sessile assemblages and is very likely to affect observations. In this study, we developed a multistate dynamic site occupancy model to estimate transition probabilities that accounts for resampling errors associated with local community structure. The model applies a nonparametric multivariate kernel smoothing methodology to the latent occupancy component to estimate the local state composition near each observation point, which is assumed to determine the probability distribution of data conditional on the occurrence of resampling error. By using computer simulations, we confirmed that an observation process that depends on local community structure may bias inferences about transition probabilities. By applying the proposed model to a real dataset of intertidal sessile communities, we also showed that estimates of transition probabilities and of the properties of community dynamics may differ considerably when spatial dependence is taken into account. Our approach can even accommodate an anisotropic spatial correlation of species composition, and may serve as a basis for inferring complex nonlinear ecological dynamics.