skip to main content


Title: Landscape analyses using eDNA metabarcoding and Earth observation predict community biodiversity in California
Abstract

Ecosystems globally are under threat from ongoing anthropogenic environmental change. Effective conservation management requires more thorough biodiversity surveys that can reveal system‐level patterns and that can be applied rapidly across space and time. Using modern ecological models and community science, we integrate environmental DNA and Earth observations to produce a time snapshot of regional biodiversity patterns and provide multi‐scalar community‐level characterization. We collected 278 samples in spring 2017 from coastal, shrub, and lowland forest sites in California, a complex ecosystem and biodiversity hotspot. We recovered 16,118 taxonomic entries from eDNA analyses and compiled associated traditional observations and environmental data to assess how well they predicted alpha, beta, and zeta diversity. We found that local habitat classification was diagnostic of community composition and distinct communities and organisms in different kingdoms are predicted by different environmental variables. Nonetheless, gradient forest models of 915 families recovered by eDNA analysis and using BIOCLIM variables, Sentinel‐2 satellite data, human impact, and topographical features as predictors, explained 35% of the variance in community turnover. Elevation, sand percentage, and photosynthetic activities (NDVI32) were the top predictors. In addition to this signal of environmental filtering, we found a positive relationship between environmentally predicted families and their numbers of biotic interactions, suggesting environmental change could have a disproportionate effect on community networks. Together, these analyses show that coupling eDNA with environmental predictors including remote sensing data has capacity to test proposed Essential Biodiversity Variables and create new landscape biodiversity baselines that span the tree of life.

 
more » « less
Award ID(s):
1759756
NSF-PAR ID:
10375741
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  more » ;  ;  ;  ;  ;  ;   « less
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Ecological Applications
Volume:
31
Issue:
6
ISSN:
1051-0761
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    1. The occurrence and distributions of wildlife populations and communities are shifting as a result of global changes. To evaluate whether these shifts are negatively impacting biodiversity processes, it is critical to monitor the status, trends and effects of environmental variables on entire communities. However, modelling the dynamics of multiple species simultaneously can require large amounts of diverse data, and few modelling approaches exist to simultaneously provide species and community‐level inferences.

    2. We present an ‘integrated community occupancy model’ (ICOM) that unites principles of data integration and hierarchical community modelling in a single framework to provide inferences on species‐specific and community occurrence dynamics using multiple data sources. The ICOM combines replicated and nonreplicated detection–nondetection data sources using a hierarchical framework that explicitly accounts for different detection and sampling processes across data sources. We use simulations to compare the ICOM to previously developed hierarchical community occupancy models and single species integrated distribution models. We then apply our model to assess the occurrence and biodiversity dynamics of foliage‐gleaning birds in the White Mountain National Forest in the northeastern USA from 2010 to 2018 using three independent data sources.

    3. Simulations reveal that integrating multiple data sources in the ICOM increased precision and accuracy of species and community‐level inferences compared to single data source models, although benefits of integration were dependent on the information content of individual data sources (e.g. amount of replication). Compared to single species models, the ICOM yielded more precise species‐level estimates. Within our case study, the ICOM had the highest out‐of‐sample predictive performance compared to single species models and models that used only a subset of the three data sources.

    4. The ICOM provides more precise estimates of occurrence dynamics compared to multi‐species models using single data sources or integrated single‐species models. We further found that the ICOM had improved predictive performance across a broad region of interest with an empirical case study of forest birds. The ICOM offers an attractive approach to estimate species and biodiversity dynamics, which is additionally valuable to inform management objectives of both individual species and their broader communities.

     
    more » « less
  2. Abstract

    Correlative species distribution models are widely used to quantify past shifts in ranges or communities, and to predict future outcomes under ongoing global change. Practitioners confront a wide range of potentially plausible models for ecological dynamics, but most specific applications only consider a narrow set. Here, we clarify that certain model structures can embed restrictive assumptions about key sources of forecast uncertainty into an analysis. To evaluate forecast uncertainties and our ability to explain community change, we fit and compared 39 candidate multi‐ or joint species occupancy models to avian incidence data collected at 320 sites across California during the early 20th century and resurveyed a century later. We found massive (>20,000 LOOIC) differences in within‐time information criterion across models. Poorer fitting models omitting multivariate random effects predicted less variation in species richness changes and smaller contemporary communities, with considerable variation in predicted spatial patterns in richness changes across models. The top models suggested avian environmental associations changed across time, contemporary avian occupancy was influenced by previous site‐specific occupancy states, and that both latent site variables and species associations with these variables also varied over time. Collectively, our results recapitulate that simplified model assumptions not only impact predictive fit but may mask important sources of forecast uncertainty and mischaracterize the current state of system understanding when seeking to describe or project community responses to global change. We recommend that researchers seeking to make long‐term forecasts prioritize characterizing forecast uncertainty over seeking to present a single best guess. To do so reliably, we urge practitioners to employ models capable of characterizing the key sources of forecast uncertainty, where predictors, parameters and random effects may vary over time or further interact with previous occurrence states.

     
    more » « less
  3. Abstract Aim

    We may be able to buffer biodiversity against the effects of ongoing climate change by prioritizing the protection of habitat with diverse physical features (high geodiversity) associated with ecological and evolutionary mechanisms that maintain high biodiversity. Nonetheless, the relationships between biodiversity and habitat vary with spatial and biological context. In this study, we compare how well habitat geodiversity (spatial variation in abiotic processes and features) and climate explain biodiversity patterns of birds and trees. We also evaluate the consistency of biodiversity–geodiversity relationships across ecoregions.

    Location

    Contiguous USA.

    Time period

    2007–2016.

    Taxa studied

    Birds and trees.

    Methods

    We quantified geodiversity with remotely sensed data and generated biodiversity maps from the Forest Inventory and Analysis and Breeding Bird Survey datasets. We fitted multivariate regressions to alpha, beta and gamma diversity, accounting for spatial autocorrelation among Nature Conservancy ecoregions and relationships among taxonomic, phylogenetic and functional biodiversity. We fitted models including climate alone (temperature and precipitation), geodiversity alone (topography, soil and geology) and climate plus geodiversity.

    Results

    A combination of geodiversity and climate predictor variables fitted most forms of bird and tree biodiversity with < 10% relative error. Models using geodiversity and climate performed better for local (alpha) and regional (gamma) diversity than for turnover‐based (beta) diversity. Among geodiversity predictors, variability of elevation fitted biodiversity best; interestingly, topographically diverse places tended to have higher tree diversity but lower bird diversity.

    Main conclusions

    Although climatic predictors tended to have larger individual effects than geodiversity, adding geodiversity improved climate‐only models of biodiversity. Geodiversity was correlated with biodiversity more consistently than with climate across ecoregions, but models tended to have a poor fit in ecoregions held out of the training dataset. Patterns of geodiversity could help to prioritize conservation efforts within ecoregions. However, we need to understand the underlying mechanisms more fully before we can build models transferable across ecoregions.

     
    more » « less
  4. ABSTRACT

    Tropical forests harbour the highest levels of terrestrial biodiversity and represent some of the most complex ecosystems on Earth, with a significant portion of this diversity above ground. Although the vertical dimension is a central aspect of the ecology of forest communities, there is little consensus as to prominence, evenness, and consistency of community‐level stratification from ground to canopy. Here, we gather the results of 62 studies across the tropics to synthesise and assess broad patterns of vertical stratification of abundance and richness in vertebrates, the best studied taxonomic group for which results have not been collated previously. Our review of the literature yielded sufficient data for bats, small mammals, birds and amphibians. We show that variation in the stratification of abundance and richness exists within and among all taxa considered. Bat richness stratification was variable among studies, although bat abundance was weighted towards the canopy. Both bird richness and abundance stratification were variable, with no overriding pattern. On the contrary, both amphibians and small mammals showed consistent patterns of decline in abundance and richness towards the canopy. We descriptively characterise research trends in drivers of stratification cited or investigated within studies, finding local habitat structure and food distribution/foraging to be the most commonly attributed drivers. Further, we analyse the influence of macroecological variables on stratification patterns, finding latitude and elevation to be key predictors of bird stratification in particular. Prominent differences among taxa are likely due to taxon‐specific interactions with local drivers such as vertical habitat structure, food distribution, and vertical climate gradients, which may vary considerably across macroecological gradients such as elevation and biogeographic realm. Our study showcases the complexity with which animal communities organise within tropical forest ecosystems, while demonstrating the canopy as a critical niche space for tropical vertebrates, thereby highlighting the inherent vulnerability of tropical vertebrate communities to forest loss and canopy disturbance. We recognise that analyses were constrained due to variation in study designs and methods which produced a variety of abundance and richness metrics recorded across different arrangements of vertical strata. We therefore suggest the application of best practices for data reporting and highlight the significant effort required to fill research gaps in terms of under‐sampled regions, taxa, and environments.

     
    more » « less
  5. Abstract

    Data deficiencies among rare or cryptic species preclude assessment of community‐level processes using many existing approaches, limiting our understanding of the trends and stressors for large numbers of species. Yet evaluating the dynamics of whole communities, not just common or charismatic species, is critical to understanding and the responses of biodiversity to ongoing environmental pressures.

    A recent surge in both public science and government‐funded data collection efforts has led to a wealth of biodiversity data. However, these data collection programmes use a wide range of sampling protocols (from unstructured, opportunistic observations of wildlife to well‐structured, design‐based programmes) and record information at a variety of spatiotemporal scales. As a result, available biodiversity data vary substantially in quantity and information content, which must be carefully reconciled for meaningful ecological analysis.

    Hierarchical modelling, including single‐species integrated models and hierarchical community models, has improved our ability to assess and predict biodiversity trends and processes. Here, we highlight the emerging ‘integrated community modelling’ framework that combines both data integration and community modelling to improve inferences on species‐ and community‐level dynamics.

    We illustrate the framework with a series of worked examples. Our three case studies demonstrate how integrated community models can be used to extend the geographic scope when evaluating species distributions and community‐level richness patterns; discern population and community trends over time; and estimate demographic rates and population growth for communities of sympatric species. We implemented these worked examples using multiple software methods through the R platform via packages with formula‐based interfaces and through development of custom code in JAGS, NIMBLE and Stan.

    Integrated community models provide an exciting approach to model biological and observational processes for multiple species using multiple data types and sources simultaneously, thus accounting for uncertainty and sampling error within a unified framework. By leveraging the combined benefits of both data integration and community modelling, integrated community models can produce valuable information about both common and rare species as well as community‐level dynamics, allowing for holistic evaluation of the effects of global change on biodiversity.

     
    more » « less