Ecological associations distribution modelling of marine plankton at a global scale

Marine plankton communities form intricate networks of interacting organisms at the base of the food chain, and play a central role in regulating ocean biogeochemical cycles and climate. However, predicting plankton community shifts in response to climate change remains challenging. While species distribution models are valuable tools for predicting changes in species biogeography under climate change scenarios, they generally overlook the key role of biotic interactions, which can significantly shape ecological processes and ecosystem responses. Here, we introduce a novel statistical framework, association distribution modelling (ADM), designed to model and predict ecological associations distribution in space and time. Applied on a Tara Oceans genome-resolved metagenomics dataset, the present-day biogeography of ADM-inferred marine plankton associations revealed four major biogeographic biomes organized along a latitudinal gradient. We predicted the evolution of these biome-specific communities in response to a climate change scenario, highlighting differential responses to environmental change. Finally, we explored the functional potential of impacted plankton communities, focusing on carbon fixation, outlining the predicted evolution of its geographical distribution and implications for ecosystem function.This article is part of the theme issue 'Connected interactions: enriching food web research by spatial and social interactions'.

A quixotic view of spatial bias in modelling the distribution of species and their diversity

Ecological processes are often spatially and temporally structured, potentially leading to autocorrelation either in environmental variables or species distribution data. Because of that, spatially-biased in-situ samples or predictors might affect the outcomes of ecological models used to infer the geographic distribution of species and diversity. There is a vast heterogeneity of methods and approaches to assess and measure spatial bias; this paper aims at addressing the spatial component of data-driven biases in species distribution modelling, and to propose potential solutions to explicitly test and account for them. Our major goal is not to propose methods to remove spatial bias from the modelling procedure, which would be impossible without proper knowledge of all the processes generating it, but rather to propose alternatives to explore and handle it. In particular, we propose and describe three main strategies that may provide a fair account of spatial bias, namely: (i) how to represent spatial bias; (ii) how to simulate null models based on virtual species for testing biogeographical and species distribution hypotheses; and (iii) how to make use of spatial bias - in particular related to sampling effort - as a leverage instead of a hindrance in species distribution modelling. We link these strategies with good practice in accounting for spatial bias in species distribution modelling.

Sampling from commercial vessel routes can capture marine biodiversity distributions effectively

Collecting fine-scale occurrence data for marine species across large spatial scales is logistically challenging but is important to determine species distributions and for conservation planning. Inaccurate descriptions of species ranges could result in designating protected areas with inappropriate locations or boundaries. Optimizing sampling strategies therefore is a priority for scaling up survey approaches using tools such as environmental DNA (eDNA) to capture species distributions. In a marine context, commercial vessels, such as ferries, could provide sampling platforms allowing access to undersampled areas and repeatable sampling over time to track community changes. However, sample collection from commercial vessels could be biased and may not represent biological and environmental variability. Here, we evaluate whether sampling along Mediterranean ferry routes can yield unbiased biodiversity survey outcomes, based on perfect knowledge from a stacked species distribution model (SSDM) of marine megafauna derived from online data repositories. Simulations to allocate sampling point locations were carried out representing different sampling strategies (random vs regular), frames (ferry routes vs unconstrained), and number of sampling points. SSDMs were remade from different sampling simulations and compared with the "perfect knowledge" SSDM to quantify the bias associated with different sampling strategies. Ferry routes detected more species and were able to recover known patterns in species richness at smaller sample sizes better than unconstrained sampling points. However, to minimize potential bias, ferry routes should be chosen to cover the variability in species composition and its environmental predictors in the SSDMs. The workflow presented here can be used to design effective sampling strategies using commercial vessel routes globally for eDNA and other biodiversity survey techniques. This approach has potential to provide a cost-effective method to access remote oceanic areas on a regular basis and can recover meaningful data on spatiotemporal biodiversity patterns.

Earth and field observations underpin metapopulation dynamics in complex landscapes: Near-term study on carabids

Understanding risks to biodiversity requires predictions of the spatial distribution of species adapting to changing ecosystems and, to that end, Earth observations integrating field surveys prove essential as they provide key numbers for assessing landscape-wide biodiversity scenarios. Here, we develop, and apply to a relevant case study, a method suited to merge Earth/field observations with spatially explicit stochastic metapopulation models to study the near-term ecological dynamics of target species in complex terrains. Our framework incorporates the use of species distribution models for a reasoned estimation of the initial presence of the target species and accounts for imperfect and incomplete detection of the species presence in the study area. It also uses a metapopulation fitness function derived from Earth observation data subsuming the ecological niche of the target species. This framework is applied to contrast occupancy of two species of carabids (Pterostichus flavofemoratus, Carabus depressus) observed in the context of a large ecological monitoring program carried out within the Gran Paradiso National Park (GPNP, Italy). Results suggest that the proposed framework may indeed exploit the hallmarks of spatially explicit ecological approaches and of remote Earth observations. The model reproduces well the observed in situ data. Moreover, it projects in the near term the two species' presence both in space and in time, highlighting the features of the metapopulation dynamics of colonization and extinction, and their expected trends within verifiable timeframes.

Competition with insectivorous ants as a contributor to low songbird diversity at low elevations in the eastern Himalaya

Competitive interactions between distantly related clades could cause complementary diversity patterns of these clades over large spatial scales. One such example might be ants and birds in the eastern Himalaya; ants are very common at low elevations but almost absent at mid-elevations where the abundance of other arthropods and insectivorous bird diversity peaks. Here, we ask if ants at low elevations could compete with birds for arthropod prey. Specifically, we studied the impact of the Asian weaver ant (Oecophylla smaragdina), a common aggressive ant at low elevations. Diet analysis using molecular methods demonstrate extensive diet overlap between weaver ants and songbirds at both low and mid-elevations. Trees without weaver ants have greater non-ant arthropod abundance and leaf damage. Experimental removal of weaver ants results in an increase in the abundance of non-ant arthropods. Notably, numbers of Coleoptera and Lepidoptera were most affected by removal experiments and were prominent components of both bird and weaver ant diets. Our results suggest that songbirds and weaver ants might potentially compete with each other for arthropod prey at low elevations, thereby contributing to lower insectivorous bird diversity at low elevations in eastern Himalaya. Competition with ants may shape vertebrate diversity patterns across broad biodiversity gradients.

A pathway for multivariate analysis of ecological communities using copulas

We describe a new pathway for multivariate analysis of data consisting of counts of species abundances that includes two key components: copulas, to provide a flexible joint model of individual species, and dissimilarity-based methods, to integrate information across species and provide a holistic view of the community. Individual species are characterized using suitable (marginal) statistical distributions, with the mean, the degree of over-dispersion, and/or zero-inflation being allowed to vary among a priori groups of sampling units. Associations among species are then modeled using copulas, which allow any pair of disparate types of variables to be coupled through their cumulative distribution function, while maintaining entirely the separate individual marginal distributions appropriate for each species. A Gaussian copula smoothly captures changes in an index of association that excludes joint absences in the space of the original species variables. A permutation-based filter with exact family-wise error can optionally be used a priori to reduce the dimensionality of the copula estimation problem. We describe in detail a Monte Carlo expectation maximization algorithm for efficient estimation of the copula correlation matrix with discrete marginal distributions (counts). The resulting fully parameterized copula models can be used to simulate realistic ecological community data under fully specified null or alternative hypotheses. Distributions of community centroids derived from simulated data can then be visualized in ordinations of ecologically meaningful dissimilarity spaces. Multinomial mixtures of data drawn from copula models also yield smooth power curves in dissimilarity-based settings. Our proposed analysis pathway provides new opportunities to combine model-based approaches with dissimilarity-based methods to enhance understanding of ecological systems. We demonstrate implementation of the pathway through an ecological example, where associations among fish species were found to increase after the establishment of a marine reserve.

Biodiversity Models: What If Unsaturation Is the Rule?

Improving biodiversity predictions is essential if we are to meet the challenges posed by global change. As knowledge is key to feed models, we need to evaluate how debated theory can affect models. An important ongoing debate is whether environmental constraints limit the number of species that can coexist in a community (saturation), with recent findings suggesting that species richness in many communities might be unsaturated. Here, we propose that biodiversity models could address this issue by accounting for a duality: considering communities as unsaturated but where species composition is constrained by different scale-dependent biodiversity drivers. We identify a variety of promising advances for incorporating this duality into commonly applied biodiversity modelling approaches and improving their spatial predictions.

The majority of biodiversity modelling approaches do not explicitly address the question of saturation.

Theoretical and methodological implications of saturation or unsaturation in biodiversity modelling.

Addressing saturation or unsaturation is vital to produce more reliable conservation strategies.

Integrative community modelling frameworks may be the way forward.