skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on August 13, 2026

Title: Quantifying functionally equivalent species and ecological network dissimilarity with optimal transport distances
Abstract Quantifying the structure and dynamics of species interactions in ecological communities is fundamental to studying ecology and evolution. While there are numerous approaches to analysing ecological networks, there is not yet an approach that can (1) quantify dissimilarity in the global structure of ecological networks that range from identical species and interaction composition to zero shared species or interactions and (2) map species between such networks while incorporating additional ecological information, such as species traits or abundances.To address these challenges, we introduce the use of optimal transport distances to quantify ecological network dissimilarity and functionally equivalent species between networks. Specifically, we describe the Gromov–Wasserstein (GW) and Fused Gromov–Wasserstein (FGW) distances. We apply these optimal transport methods to synthetic and empirical data, using mammal food webs throughout sub‐Saharan Africa for illustration. We showcase the application of GW and FGW distances to identify the most functionally similar species between food webs, incorporate additional trait information into network comparisons and quantify food web dissimilarity among geographic regions.Our results demonstrate that GW and FGW distances can effectively differentiate ecological networks based on their topological structure while identifying functionally equivalent species, even when networks have different species. The FGW distance further improves node mapping for basal species by incorporating node‐level traits. We show that these methods allow for a more nuanced understanding of the topological similarities in food web networks among geographic regions compared to an alternative measure of network dissimilarity based on species identities.Optimal transport distances offer a new approach for quantifying functional equivalence between networks and a measure of network dissimilarity suitable for a broader range of uses than existing approaches. OT methods can be harnessed to analyse ecological networks at large spatial scales and compare networks among ecosystems, realms or taxa. Optimal transport‐based distances, therefore, provide a powerful tool for analysing ecological networks with great potential to advance our understanding of ecological community structure and dynamics in a changing world.  more » « less
Award ID(s):
2213568
PAR ID:
10643476
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
British Ecological Society
Date Published:
Journal Name:
Methods in Ecology and Evolution
ISSN:
2041-210X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Earth's biosphere is undergoing drastic reorganization due to the sixth mass extinction brought on by the Anthropocene. Impacts of local and regional extirpation of species have been demonstrated to propagate through the complex interaction networks they are part of, leading to secondary extinctions and exacerbating biodiversity loss. Contemporary ecological theory has developed several measures to analyse the structure and robustness of ecological networks under biodiversity loss. However, a toolbox for directly simulating and quantifying extinction cascades and creating novel interactions (i.e. rewiring) remains absent.Here, we presentNetworkExtinction—a novel R package which we have developed to explore the propagation of species extinction sequences through ecological networks and quantify the effects of rewiring potential in response to primary species extinctions. WithNetworkExtinction, we integrate ecological theory and computational simulations to develop functionality with which users may analyse and visualize the structure and robustness of ecological networks. The core functions introduced withNetworkExtinctionfocus on simulations of sequential primary extinctions and associated secondary extinctions, allowing user‐specified secondary extinction thresholds and realization of rewiring potential.With the packageNetworkExtinction, users can estimate the robustness of ecological networks after performing species extinction routines based on several algorithms. Moreover, users can compare the number of simulated secondary extinctions against a null model of random extinctions. In‐built visualizations enable graphing topological indices calculated by the deletion sequence functions after each simulation step. Finally, the user can estimate the network's degree distribution by fitting different common distributions. Here, we illustrate the use of the package and its outputs by analysing a Chilean coastal marine food web.NetworkExtinctionis a compact and easy‐to‐use R package with which users can quantify changes in ecological network structure in response to different patterns of species loss, thresholds and rewiring potential. Therefore, this package is particularly useful for evaluating ecosystem responses to anthropogenic and environmental perturbations that produce nonrandom and sometimes targeted, species extinctions. 
    more » « less
  2. Abstract Changes in trophic niche—the pathways through which an organism obtains energy and nutrients—are a fundamental way in which organisms respond to environmental conditions. But the capacity for species to alter their trophic niches in response to global change, and the ways they do so when able, remain largely unknown.Here we examine food webs in three long‐term and large‐scale experiments to test how resource availability and nutritional requirements interact to determine an organism's trophic niche in the context of one of the largest global trends in land use—the rise in bioenergy production.We use carbon and nitrogen stable isotope analyses to characterize arthropod food webs across three biofuel crops representing a gradient in plant resource richness (corn monocultures, fields dominated by native switchgrass and restored prairie), and to quantify changes in the trophic niche of a widespread generalist ant species across habitats. In doing so, we measure the effects of basal resource richness on food chain length, niche breadth and trophic position. We frame our results in the context of two hypotheses that explain variation in trophic niche—the niche variation hypothesis which emphasizes the importance of resource diversity and ecological opportunity, and the optimal diet hypothesis which emphasizes dietary constraints and the availability of optimal resources.Increasing plant richness lengthened food chains by 10%–20% compared to monocultures. Niche breadths of generalist ants did not vary with resource richness, suggesting they were limited by optimal diet requirements and constraints rather than by ecological opportunity. The ants instead responded to changes in plant richness by shifting their estimated trophic position. In resource‐poor monocultures, the ants were top predators, sharing a trophic position with predatory spiders. In resource‐rich environments, in contrast, the ants were omnivores, relying on a mix of animal prey and plant‐based resources.In addition to highlighting novel ecosystem impacts of alternate bioenergy landscapes, our results suggest that niche breadth and trophic diversification depend more on the presence of optimal resources than on ecological opportunity alone. 
    more » « less
  3. Abstract The Gromov–Wasserstein distance—a generalization of the usual Wasserstein distance—permits comparing probability measures defined on possibly different metric spaces. Recently, this notion of distance has found several applications in Data Science and in Machine Learning. With the goal of aiding both the interpretability of dissimilarity measures computed through the Gromov–Wasserstein distance and the assessment of the approximation quality of computational techniques designed to estimate the Gromov–Wasserstein distance, we determine the precise value of a certain variant of the Gromov–Wasserstein distance between unit spheres of different dimensions. Indeed, we consider a two-parameter family$$\{d_{{{\text {GW}}}p,q}\}_{p,q=1}^{\infty }$$ { d GW p , q } p , q = 1 of Gromov–Wasserstein distances between metric measure spaces. By exploiting a suitable interaction between specific values of the parameterspandqand the metric of the underlying spaces, we are able to determine the exact value of the distance$$d_{{{\text {GW}}}4,2}$$ d GW 4 , 2 between all pairs of unit spheres of different dimensions endowed with their Euclidean distance and their uniform measure. 
    more » « less
  4. Abstract Global change is impacting biodiversity across all habitats on earth. New selection pressures from changing climatic conditions and other anthropogenic activities are creating heterogeneous ecological and evolutionary responses across many species' geographic ranges. Yet we currently lack standardised and reproducible tools to effectively predict the resulting patterns in species vulnerability to declines or range changes.We developed an informatic toolbox that integrates ecological, environmental and genomic data and analyses (environmental dissimilarity, species distribution models, landscape connectivity, neutral and adaptive genetic diversity, genotype‐environment associations and genomic offset) to estimate population vulnerability. In our toolbox, functions and data structures are coded in a standardised way so that it is applicable to any species or geographic region where appropriate data are available, for example individual or population sampling and genomic datasets (e.g. RAD‐seq, ddRAD‐seq, whole genome sequencing data) representing environmental variation across the species geographic range.To demonstrate multi‐species applicability, we apply our toolbox to three georeferenced genomic datasets for co‐occurring East African spiny reed frogs (Afrixalus fornasini, A. delicatusandA. sylvaticus) to predict their population vulnerability, as well as demonstrating that range loss projections based on adaptive variation can be accurately reproduced from a previous study using data for two European bat species (Myotis escaleraiandM. crypticus).Our framework sets the stage for large scale, multi‐species genomic datasets to be leveraged in a novel climate change vulnerability framework to quantify intraspecific differences in genetic diversity, local adaptation, range shifts and population vulnerability based on exposure, sensitivity and landscape barriers. 
    more » « less
  5. Abstract Optimal transport (OT) methods seek a transformation map (or plan) between two probability measures, such that the transformation has the minimum transportation cost. Such a minimum transport cost, with a certain power transform, is called the Wasserstein distance. Recently, OT methods have drawn great attention in statistics, machine learning, and computer science, especially in deep generative neural networks. Despite its broad applications, the estimation of high‐dimensional Wasserstein distances is a well‐known challenging problem owing to the curse‐of‐dimensionality. There are some cutting‐edge projection‐based techniques that tackle high‐dimensional OT problems. Three major approaches of such techniques are introduced, respectively, the slicing approach, the iterative projection approach, and the projection robust OT approach. Open challenges are discussed at the end of the review. This article is categorized under:Statistical and Graphical Methods of Data Analysis > Dimension ReductionStatistical Learning and Exploratory Methods of the Data Sciences > Manifold Learning 
    more » « less