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.


Title: NetworkExtinction : An R package to simulate extinction propagation and rewiring potential in ecological networks
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
Award ID(s):
2224915 2129757
PAR ID:
10418656
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley-Blackwell
Date Published:
Journal Name:
Methods in Ecology and Evolution
Volume:
14
Issue:
8
ISSN:
2041-210X
Page Range / eLocation ID:
p. 1952-1966
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract A key challenge in conservation biology is that not all species are equally likely to go extinct when faced with a disturbance, but there are multiple overlapping reasons for such differences in extinction probability. Differences in species extinction risk may represent extinction selectivity, a non‐random process by which species’ risks of extinction are caused by differences in fitness based on traits. Additionally, rare species with low abundances and/or occupancies are more likely to go extinct than common species for reasons of random chance alone, that is, bad luck. Unless ecologists and conservation biologists can disentangle random and selective extinction processes, then the prediction and prevention of future extinctions will continue to be an elusive challenge.We suggest that a modified version of a common null model procedure, rarefaction, can be used to disentangle the influence of stochastic species loss from selective non‐random processes. To this end we applied a rarefaction‐based null model to three published data sets to characterize the influence of species rarity in driving biodiversity loss following three biodiversity loss events: (a) disease‐associated bat declines; (b) disease‐associated amphibian declines; and (c) habitat loss and invasive species‐associated gastropod declines. For each case study, we used rarefaction to generate null expectations of biodiversity loss and species‐specific extinction probabilities.In each of our case studies, we find evidence for both random and non‐random (selective) extinctions. Our findings highlight the importance of explicitly considering that some species extinctions are the result of stochastic processes. In other words, we find significant evidence for bad luck in the extinction process.Policy implications. Our results suggest that rarefaction can be used to disentangle random and non‐random extinctions and guide management decisions. For example, rarefaction can be used retrospectively to identify when declines of at‐risk species are likely to result from selectivity, versus random chance. Rarefaction can also be used prospectively to formulate minimum predictions of species loss in response to hypothetical disturbances. Given its minimal data requirements and familiarity among ecologists, rarefaction may be an efficient and versatile tool for identifying and protecting species that are most vulnerable to global extinction. 
    more » « less
  2. 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
  3. Abstract Habitat loss disrupts species interactions through local extinctions, potentially orphaning species that depend on interacting partners, via mutualisms or commensalisms, and increasing secondary extinction risk. Orphaned species may become functionally or secondarily extinct, increasing the severity of the current biodiversity crisis. While habitat destruction is a major cause of biodiversity loss, the number of secondary extinctions is largely unknown. We investigate the relationship between habitat loss, orphaned species, and bipartite network properties. Using a real seed dispersal network, we simulate habitat loss to estimate the rate at which species are orphaned. To be able to draw general conclusions, we also simulate habitat loss in synthetic networks to quantify how changes in network properties affect orphan rates across broader parameter space. Both real and synthetic network simulations show that even small amounts of habitat loss can cause up to 10% of species to be orphaned. More area loss, less connected networks, and a greater disparity in the species richness of the network's trophic levels generally result in more orphaned species. As habitat is lost to land‐use conversion and climate change, more orphaned species increase the loss of community‐level and ecosystem functions. However, the potential severity of repercussions ranges from minimal (no species orphaned) to catastrophic (up to 60% of species within a network orphaned). Severity of repercussions also depends on how much the interaction richness and intactness of the community affects the degree of redundancy within networks. Orphaned species could add substantially to the loss of ecosystem function and secondary extinction worldwide. 
    more » « less
  4. Abstract Human-driven threats are changing biodiversity, impacting ecosystem services. The loss of one species can trigger secondary extinctions of additional species, because species interact–yet the consequences of these secondary extinctions for services remain underexplored. Herein, we compare robustness of food webs and the ecosystem services (hereafter ‘services’) they provide; and investigate factors determining service responses to secondary extinctions. Simulating twelve extinction scenarios for estuarine food webs with seven services, we find that food web and service robustness are highly correlated, but that robustness varies across services depending on their trophic level and redundancy. Further, we find that species providing services do not play a critical role in stabilizing food webs – whereas species playing supporting roles in services through interactions are critical to the robustness of both food webs and services. Together, our results reveal indirect risks to services through secondary species losses and predictable differences in vulnerability across services. 
    more » « less
  5. Abstract Declines in grassland diversity in response to nutrient addition are a general consequence of global change. This decline in species richness may be driven by multiple underlying processes operating at different time‐scales. Nutrient addition can reduce diversity by enhancing the rate of local extinction via competitive exclusion, or by reducing the rate of colonization by constraining the pool of species able to colonize under new conditions. Partitioning net change into extinction and colonization rates will better delineate the long‐term effect of global change in grasslands.We synthesized changes in richness in response to experimental fertilization with nitrogen, phosphorus and potassium with micronutrients across 30 grasslands. We quantified changes in local richness, colonization, and extinction over 8–10 years of nutrient addition, and compared these rates against control conditions to isolate the effect of nutrient addition from background dynamics.Total richness at steady state in the control plots was the sum of equal, relatively high rates of local colonization and extinction. On aggregate, 30%–35% of initial species were lost and the same proportion of new species were gained at least once over a decade. Absolute turnover increased with site‐level richness but was proportionately greater at lower‐richness sites relative to starting richness. Loss of total richness with nutrient addition, especially N in combination with P or K, was driven by enhanced rates of extinction with a smaller contribution from reduced colonization. Enhanced extinction and reduced colonization were disproportionately among native species, perennials, and forbs. Reduced colonization plateaued after the first few (<5) years after nutrient addition, while enhanced extinction continued throughout the first decade.Synthesis. Our results indicate a high rate of colonizations and extinctions underlying the richness of ambient communities and that nutrient enhancement drives overall declines in diversity primarily by exclusion of previously established species. Moreover, enhanced extinction continues over long time‐scales, suggesting continuous, long‐term community responses and a need for long‐term study to fully realize the extinction impact of increased nutrients on grassland composition. 
    more » « less