Joint species distribution models quantify drivers of aquatic macroinvertebrate assemblage composition
On this page:
Quantifying drivers of biotic assemblage composition is critical for sound ecological assessment. Potential drivers of composition are often evaluated by aggregating taxon-level responses along environmental gradients. Because this approach neglects biotic interactions, the effects of an environmental variable could be over- or under-estimated. Here, we use joint species distribution models (jSDM) to account for biotic interactions and include taxonomic relatedness and traits information to quantify the effect of several environmental gradients on benthic macroinvertebrate assemblages in streams throughout the contiguous U.S. We fit separate jSDMs to nine major ecoregions using data collected by the U.S. Environmental Protection Agency’s National Rivers and Streams Assessment program. We then evaluated how benthic macroinvertebrate assemblages vary along nutrient, salinity, climate, and physical habitat gradients and whether specific traits facilitate or exclude a taxon from a location. We found that jSDMs sufficiently differentiate between presence and absence locations for most genera and that the relative importance of the environmental gradients varies among ecoregions. However, in many ecoregions a relatively large proportion of genera responded to mean summer temperature and substrate diameter. Further, we found evidence that the expression of some traits by a genus is related to how they respond along an environmental gradient. Generally, genus richness increased with increasing mean substrate diameter but declined with increasing total phosphorous and sulfate concentrations. Our results showcase a new approach to modeling biotic assemblages in aquatic ecosystems and provide insights into how they might change under anthropogenic disturbances.