Joint dynamic species distribution models: a tool for community ordination and spatio‐temporal monitoring

Aim Spatial analysis of the distribution and density of species is of continuing interest within theoretical and applied ecology. Species distribution models (SDMs) are being increasingly used to analyse count, presence–absence and presence-only data sets. There is a growing literature on dynamic SDMs (which incorporate temporal variation in species distribution), joint SDMs (which simultaneously analyse the correlated distribution of multiple species) and geostatistical models (which account for similarity between nearby sites caused by unobserved covariates). However, no previous study has combined all three attributes within a single framework. Innovation We develop spatial dynamic factor analysis for use as a ‘joint, dynamic SDM’ (JDSDM), which uses geostatistical methods to account for spatial similarity when estimating one or more ‘factors’. Each factor evolves over time following a density-dependent (Gompertz) process, and the log-density of each species is approximated as a linear combination of different factors. We demonstrate a JDSDM using two multispecies case studies (an annual survey of bottom-associated species in the Bering Sea and a seasonal survey of butterfly density in the continental USA), and also provide our code publicly as an R package. Main conclusions Case study applications show that that JDSDMs can be used for species ordination, i.e. showing that dynamics for butterfly species within the same genus are significantly more correlated than for species from different genera. We also demonstrate how JDSDMs can rapidly identify dominant patterns in community dynamics, including the decline and recovery of several Bering Sea fishes since 2008, and the ‘flight curves’ typical of early or late-emerging butterflies. We conclude by suggesting future research that could incorporate phylogenetic relatedness or functional similarity, and propose that our approach could be used to monitor community dynamics at large spatial and temporal scales.

[1]  D. Doak,et al.  The Statistical Inevitability of Stability‐Diversity Relationships in Community Ecology , 1998, The American Naturalist.

[2]  Laura J. Pollock,et al.  Understanding co‐occurrence by modelling species simultaneously with a Joint Species Distribution Model (JSDM) , 2014 .

[3]  Brett J Furnas,et al.  Detecting diversity: emerging methods to estimate species diversity. , 2014, Trends in ecology & evolution.

[4]  David B. Roy,et al.  Phenology of British butterflies and climate change , 2000 .

[5]  J. Thorson,et al.  The importance of spatial models for estimating the strength of density dependence. , 2015, Ecology.

[6]  Paul Galpern,et al.  Climate change impacts on bumblebees converge across continents , 2015, Science.

[7]  Kevin A. Thompson,et al.  Evaluating Temporal Consistency in Marine Biodiversity Hotspots , 2015, PloS one.

[8]  Mark H. DeVisser,et al.  A dynamic species distribution model of Glossina subgenus Morsitans: The identification of tsetse reservoirs and refugia , 2010 .

[9]  E. Pollard.  A method for assessing changes in the abundance of butterflies , 1977 .

[10]  Cory Merow,et al.  Developing Dynamic Mechanistic Species Distribution Models: Predicting Bird-Mediated Spread of Invasive Plants across Northeastern North America , 2011, The American Naturalist.

[11]  S. Schneider,et al.  Fingerprints of global warming on wild animals and plants , 2003, Nature.

[12]  J. Andrew Royle,et al.  Estimating Size and Composition of Biological Communities by Modeling the Occurrence of Species , 2005 .

[13]  Brian J. McGill,et al.  Can niche-based distribution models outperform spatial interpolation? , 2007 .

[14]  Antoine Guisan,et al.  Niche dynamics in space and time. , 2008, Trends in ecology & evolution.

[15]  Hugh P Possingham,et al.  Zero tolerance ecology: improving ecological inference by modelling the source of zero observations. , 2005, Ecology letters.

[16]  Otso Ovaskainen,et al.  Making more out of sparse data: hierarchical modeling of species communities. , 2011, Ecology.

[17]  J. Andrew Royle,et al.  Biodiversity of man-made open habitats in an underused country: a class of multispecies abundance models for count data , 2012, Biodiversity and Conservation.

[18]  C. Dormann.  Effects of incorporating spatial autocorrelation into the analysis of species distribution data , 2007 .

[19]  M. StachuraMegan,et al.  Oceanographic influences on patterns in North Pacific salmon abundance , 2014 .

[20]  R. Hilborn,et al.  Integrating scientific guidance into marine spatial planning , 2014, Proceedings of the Royal Society B: Biological Sciences.

[21]  Francis K. C. Hui,et al.  So Many Variables: Joint Modeling in Community Ecology. , 2015, Trends in ecology & evolution.

[22]  D. Schindler,et al.  Watershed geomorphology and snowmelt control stream thermal sensitivity to air temperature , 2015 .

[23]  David J. Harris.  Generating realistic assemblages with a joint species distribution model , 2015 .

[24]  S. Carpenter,et al.  ESTIMATING COMMUNITY STABILITY AND ECOLOGICAL INTERACTIONS FROM TIME‐SERIES DATA , 2003 .

[25]  A. Zuur,et al.  Dynamic factor analysis to estimate common trends in fisheries time series , 2003 .

[26]  D. Earn,et al.  Coherence and conservation. , 2000, Science.

[27]  Anders Nielsen,et al.  TMB: Automatic Differentiation and Laplace Approximation , 2015, 1509.00660.

[28]  H. Rue,et al.  An explicit link between Gaussian fields and Gaussian Markov random fields: the stochastic partial differential equation approach , 2011 .

[29]  Tavis D. Forrester,et al.  A two-species occupancy model accommodating simultaneous spatial and interspecific dependence. , 2016, Ecology.

[30]  Hans J. Skaug,et al.  Spatial factor analysis: a new tool for estimating joint species distributions and correlations in species range , 2015 .

[31]  K. Gross,et al.  Do growing degree days predict phenology across butterfly species , 2015 .

[32]  Michael J. Fogarty,et al.  Marine Taxa Track Local Climate Velocities , 2013, Science.

[33]  David B. Roy,et al.  Uncovering hidden spatial structure in species communities with spatially explicit joint species distribution models , 2016 .

[34]  Carsten F. Dormann,et al.  Towards novel approaches to modelling biotic interactions in multispecies assemblages at large spatial extents , 2012 .

[35]  David B. Dunson,et al.  Using latent variable models to identify large networks of species‐to‐species associations at different spatial scales , 2016 .

[36]  J. Andrew Royle,et al.  Estimating species richness and accumulation by modeling species occurrence and detectability. , 2006, Ecology.

[37]  J. Elith,et al.  Species Distribution Models: Ecological Explanation and Prediction Across Space and Time , 2009 .

[38]  J. Andrew Royle,et al.  Multi-species occurrence models to evaluate the effects of conservation and management actions , 2010, Biological Conservation.

[39]  J. Thorson,et al.  Spatial semiparametric models improve estimates of species abundance and distribution , 2014 .

[40]  J Andrew Royle,et al.  Models for inference in dynamic metacommunity systems. , 2010, Ecology.

[41]  C. Parmesan,et al.  Poleward shifts in geographical ranges of butterfly species associated with regional warming , 1999, Nature.

[42]  Hans J. Skaug,et al.  Automatic approximation of the marginal likelihood in non-Gaussian hierarchical models , 2006, Comput. Stat. Data Anal..

[43]  M. Hooten,et al.  Statistical Agent-Based Models for Discrete Spatio-Temporal Systems , 2010 .

[44]  Otso Ovaskainen,et al.  Modeling species co-occurrence by multivariate logistic regression generates new hypotheses on fungal interactions. , 2010, Ecology.

[45]  Jay M. Ver Hoef,et al.  Using spatiotemporal statistical models to estimate animal abundance and infer ecological dynamics from survey counts , 2015 .

[46]  Kai Zhu,et al.  More than the sum of the parts: forest climate response from joint species distribution models. , 2014, Ecological applications : a publication of the Ecological Society of America.

[47]  W. Wooster,et al.  Year‐to‐year variations in Bering Sea ice cover and some consequences for fish distributions , 1998 .

[48]  A M Latimer,et al.  Hierarchical models facilitate spatial analysis of large data sets: a case study on invasive plant species in the northeastern United States. , 2009, Ecology letters.

[49]  Jorge Soberón.  Grinnellian and Eltonian niches and geographic distributions of species. , 2007, Ecology letters.

[50]  Peter C. M. Molenaar,et al.  A dynamic factor model for the analysis of multivariate time series , 1985 .