skip to main content


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
NSF-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

    Mutualistic relationships, such as those between plants and pollinators, may be vulnerable to the local extinctions predicted under global environmental change. However, network theory predicts that plant–pollinator networks can withstand species loss if pollinators switch to alternative floral resources (rewiring). Whether rewiring occurs following species loss in natural communities is poorly known because replicated species exclusions are difficult to implement at appropriate spatial scales.

    We experimentally removed a hummingbird‐pollinated plant,Heliconia tortuosa, from within tropical forest fragments to investigate how hummingbirds respond to temporary loss of an abundant resource. Under therewiring hypothesis, we expected that behavioural flexibility would allow hummingbirds to use alternative resources, leading to decreased ecological specialization and reorganization of the network structure (i.e. pairwise interactions). Alternatively, morphological or behavioural constraints—such as trait‐matching or interspecific competition—might limit the extent to which hummingbirds alter their foraging behaviour.

    We employed a replicated Before‐After‐Control‐Impact experimental design and quantified plant–hummingbird interactions using two parallel sampling methods: pollen collected from individual hummingbirds (‘pollen networks’, created from >300 pollen samples) and observations of hummingbirds visiting focal plants (‘camera networks’, created from >19,000 observation hours). To assess the extent of rewiring, we quantified ecological specialization at the individual, species and network levels and examined interaction turnover (i.e. gain/loss of pairwise interactions).

    H. tortuosaremoval caused some reorganization of pairwise interactions but did not prompt large changes in specialization, despite the large magnitude of our manipulation (on average, >100 inflorescences removed in exclusion areas of >1 ha). Although some individual hummingbirds sampled through time showed modest increases in niche breadth followingHeliconiaremoval (relative to birds that did not experience resource loss), these changes were not reflected in species‐ and network‐level specialization metrics.

    Our results suggest that, at least over short time‐scales, animals may not necessarily shift to alternative resources after losing an abundant food resource—even in species thought to be highly opportunistic foragers, such as hummingbirds. Given that rewiring contributes to theoretical predictions of network stability, future studies should investigate why pollinators might not expand their diets after a local resource extinction.

     
    more » « less
  2. Abstract

    Biodiversity studies rely heavily on estimates of species' distributions often obtained through ecological niche modelling. Numerous software packages exist that allow users to model ecological niches using machine learning and statistical methods. However, no existing package with a graphical user interface allows users to perform model calibration and selection based on convex forms such as ellipsoids, which may match fundamental ecological niche shapes better, incorporating tools for exploring, modelling, and evaluating niches and distributions that are intuitive for both novice and proficient users.

    Here we describe anrpackage, NicheToolBox(ntbox), that allows users to conduct all processing steps involved in ecological niche modelling: downloading and curating occurrence data, obtaining and transforming environmental data layers, selecting environmental variables, exploring relationships between geographic and environmental spaces, calibrating and selecting ellipsoid models, evaluating models using binomial and partial ROC tests, assessing extrapolation risk, and performing geographic information system operations via a graphical user interface. A summary of the entire workflow is produced for use as a stand‐alone algorithm or as part of research reports.

    The method is explained in detail and tested via modelling the threatened feline speciesLeopardus wiedii. Georeferenced occurrence data for this species are queried to display both point occurrences and the IUCN extent of occurrence polygon (IUCN, 2007). This information is used to illustrate tools available for accessing, processing and exploring biodiversity data (e.g. number of occurrences and chronology of collecting) and transforming environmental data (e.g. a summary PCA for 19 bioclimatic layers). Visualizations of three‐dimensional ecological niches modelled as minimum volume ellipsoids are developed with ancillary statistics. This niche model is then projected to geographic space, to represent a corresponding potential suitability map.

    Usingntboxallows a fast and straightforward means by which to retrieve and manipulate occurrence and environmental data, which can then be implemented in model calibration, projection and evaluation for assessing distributions of species in geographic space and their corresponding environmental combinations.

     
    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

    Quantitative evaluations to optimize complexity have become standard for avoiding overfitting of ecological niche models (ENMs) that estimate species’ potential geographic distributions.ENMevalwas the first R package to make such evaluations (often termed model tuning) widely accessible for the Maxent algorithm. It also provided multiple methods for partitioning occurrence data and reported various performance metrics.

    Requests by users, recent developments in the field, and needs for software compatibility led to a major redesign and expansion. We additionally conducted a literature review to investigate trends inENMevaluse (2015–2019).

    ENMeval2.0 has a new object‐oriented structure for adding other algorithms, enables customizing algorithmic settings and performance metrics, generates extensive metadata, implements a null‐model approach to quantify significance and effect sizes, and includes features to increase the breadth of analyses and visualizations. In our literature review, we found insufficient reporting of model performance and parameterization, heavy reliance on model selection with AICc and low utilization of spatial cross‐validation; we explain howENMeval2.0 can help address these issues.

    This redesigned and expanded version can promote progress in the field and improve the information available for decision‐making.

     
    more » « less
  5. 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