Examining the potential effects of species aggregation on the network structure of food webs

One of the key measures that have been used to describe the topological properties of complex networks is the “degree distribution”, which is a measure that describes the frequency distribution of number of links per node. Food webs are complex ecological networks that describe the trophic relationships among species in a community, and the topological properties of empirical food webs, including degree distributions, have been examined previously. Previously, the “niche model” has been shown to accurately predict degree distributions of empirical food webs, however, the niche model-generated food webs were referenced against empirical food webs that had their species grouped together based on their taxonomic and/or trophic relationships (aggregated food webs). Here, we explore the effects of species aggregation on the ability of the niche model to predict the total- (sum of prey and predator links per node), in- (number of predator links per node), and out- (number of prey links per node) degree distributions of empirical food webs by examining two food webs that can be aggregated at different levels of resolution. The results showed that (1) the cumulative total- and out-degree distributions were consistent with the niche model predictions when the species were aggregated, (2) when the species were disaggregated (i.e., higher resolution), there were mixed conclusions with regards to the niche model's ability to predict total- and out-degree distributions, (3) the model's ability to predict the in-degree distributions of the two food webs was generally inadequate. Although it has been argued that universal functional form based on the niche model could describe the degree distribution patterns of empirical food webs, we believe there are some limitations to the model's ability to accurately predict the structural properties of food webs.     

Evidence for the existence of a robust pattern of prey selection in food webs

Food webs aim to provide a thorough representation of the trophic interactions found in an ecosystem. The complexity of empirical food webs, however, is leading many ecologists to focus dynamic ecosystem studies on smaller microcosm or mesocosm studies based upon community modules, which comprise three to five species and the interactions likely to have ecological relevance. We provide here a structural counterpart to community modules. We investigate food-web 'motifs' which are n-species connected subgraphs found within the food web. Remarkably, we find that the over- and under-representation of three-species motifs in empirical food webs can be understood through comparison to a static food-web model, the niche model. Our result conclusively demonstrates that predation upon species with some 'characteristic' niche value is the prey selection mechanism consistent with the structural properties of empirical food webs.     

Current food web models cannot explain the overall topological structure of observed food webs

Topological food webs illustrating “who eats whom” in different systems exhibit similar, non‐random, structures suggesting that general rules govern food web structure. Current food web models correctly predict many measures of food web topology from knowledge of species richness and connectance (fraction of possible predator–prey links that actually occur), together with assumptions about the ecological rules governing “who eats whom”. However, current measures are relatively insensitive to small changes in topology. Improvement of, and discrimination among, current models requires development of new measures of food web structure. Here I examine whether current food web models (cascade, niche, and nested hierarchy models, plus a random null model) can predict a new measure of food web structure, structural stability. Structural stability complements other measures of food web topology because it is sensitive to changes in topology that other measures often miss. The cascade and null models respectively over‐ and underpredict structural stability for a set of 17 high‐quality food webs. While the niche and nested hierarchy models provide unbiased predictions on average, their 95% confidence intervals frequently fail to include the observed data. Observed structural stabilities for all models are overdispersed compared to model predictions, and predicted and observed structural stabilities are uncorrelated, indicating that important sources of variation in structural stability are not captured by the models. Crucially, poor model performance arises because observed variation in structural stability is unrelated to variation in species richness and connectance. In contrast, almost all other measures of food web topology vary with species richness and connectance in natural webs. No model that takes species richness and connectance as the only input parameters can reproduce observed variation in structural stability. Further progress in predicting and explaining food web topology will require fundamentally new models based on different input parameters.     

Latitudinal gradients in biotic niche breadth vary across ecosystem types

Several properties of food webs—the networks of feeding links between species—are known to vary systematically with the species richness of the underlying community. Under the ‘latitude–niche breadth hypothesis’, which predicts that species in the tropics will tend to evolve narrower niches, one might expect that these scaling relationships could also be affected by latitude. To test this hypothesis, we analysed the scaling relationships between species richness and average generality, vulnerability and links per species across a set of 196 empirical food webs. In estuarine, marine and terrestrial food webs there was no effect of latitude on any scaling relationship, suggesting constant niche breadth in these habitats. In freshwater communities, on the other hand, there were strong effects of latitude on scaling relationships, supporting the latitude–niche breadth hypothesis. These contrasting findings indicate that it may be more important to account for habitat than latitude when exploring gradients in food-web structure.     

Compilation and network analyses of cambrian food webs

A rich body of empirically grounded theory has developed about food webs—the networks of feeding relationships among species within habitats. However, detailed food-web data and analyses are lacking for ancient ecosystems, largely because of the low resolution of taxa coupled with uncertain and incomplete information about feeding interactions. These impediments appear insurmountable for most fossil assemblages; however, a few assemblages with excellent soft-body preservation across trophic levels are candidates for food-web data compilation and topological analysis. Here we present plausible, detailed food webs for the Chengjiang and Burgess Shale assemblages from the Cambrian Period. Analyses of degree distributions and other structural network properties, including sensitivity analyses of the effects of uncertainty associated with Cambrian diet designations, suggest that these early Paleozoic communities share remarkably similar topology with modern food webs. Observed regularities reflect a systematic dependence of structure on the numbers of taxa and links in a web. Most aspects of Cambrian food-web structure are well-characterized by a simple “niche model,” which was developed for modern food webs and takes into account this scale dependence. However, a few aspects of topology differ between the ancient and recent webs: longer path lengths between species and more species in feeding loops in the earlier Chengjiang web, and higher variability in the number of links per species for both Cambrian webs. Our results are relatively insensitive to the exclusion of low-certainty or random links. The many similarities between Cambrian and recent food webs point toward surprisingly strong and enduring constraints on the organization of complex feeding interactions among metazoan species. The few differences could reflect a transition to more strongly integrated and constrained trophic organization within ecosystems following the rapid diversification of species, body plans, and trophic roles during the Cambrian radiation. More research is needed to explore the generality of food-web structure through deep time and across habitats, especially to investigate potential mechanisms that could give rise to similar structure, as well as any differences.     

The Probabilistic Niche Model Reveals the Niche Structure and Role of Body Size in a Complex Food Web

The niche model has been widely used to model the structure of complex food webs, and yet the ecological meaning of the single niche dimension has not been explored. In the niche model, each species has three traits, niche position, diet position and feeding range. Here, a new probabilistic niche model, which allows the maximum likelihood set of trait values to be estimated for each species, is applied to the food web of the Benguela fishery. We also developed the allometric niche model, in which body size is used as the niche dimension. About 80% of the links in the empirical data are predicted by the probabilistic niche model, a significant improvement over recent models. As in the niche model, species are uniformly distributed on the niche axis. Feeding ranges are exponentially distributed, but diet positions are not uniformly distributed below the predator. Species traits are strongly correlated with body size, but the allometric niche model performs significantly worse than the probabilistic niche model. The best-fit parameter set provides a significantly better model of the structure of the Benguela food web than was previously available. The methodology allows the identification of a number of taxa that stand out as outliers either in the model''s poor performance at predicting their predators or prey or in their parameter values. While important, body size alone does not explain the structure of the one-dimensional niche.     

Defining and measuring trophic role similarity in food webs using regular equivalence

We present a graph theoretic model of analysing food web structure called regular equivalence. Regular equivalence is a method for partitioning the species in a food web into "isotrophic classes" that play the same structural roles, even if they are not directly consuming the same prey or if they do not share the same predators. We contrast regular equivalence models, in which two species are members of the same trophic group if they have trophic links to the same set of other trophic groups, with structural equivalence models, in which species are equivalent if they are connected to the exact same other species. Here, the regular equivalence approach is applied to two published food webs: (1) a topological web (Malaysian pitcher plant insect food web) and (2) a carbon-flow web (St. Marks, Florida seagrass ecosystem food web). Regular equivalence produced a more satisfactory set of classes than did the structural approach, grouping basal taxa with other basal taxa and not with top predators. Regular equivalence models provide a way to mathematically formalize trophic position, trophic group and trophic niche. These models are part of a family of models that includes structural models used extensively by ecologists now. Regular equivalence models uncover similarities in trophic roles at a higher level of organization than do the structural models. The approach outlined is useful for measuring the trophic roles of species in food web models, measuring similarity in trophic relations of two or more species, comparing food webs over time and across geographic regions, and aggregating taxa into trophic groups that reduce the complexity of ecosystem feeding relations without obscuring network relationships. In addition, we hope the approach will prove useful in predicting the outcome of predator-prey interactions in experimental studies.     

Emergence of non-random structure in local food webs generated from randomly structured regional webs

Previous studies have shown that high-resolution, empirical food webs possess a non-random network structure, typically characterized by uniform or exponential degree distributions. However, the empirical food webs that have been investigated for their structural properties represent local communities that are only a subset of a larger pool of regionally coexisting species. Here, we use a simple model to investigate the effects of regional food web structure on local food webs that are assembled by two simple processes: random immigration of species from a source web (regional food web), and random extinction of species within the local web. The model shows that local webs with non-random degree distributions can arise from randomly structured source webs. A comparison of local webs assembled from randomly structured source webs with local webs assembled from source webs generated by the niche model shows that the former have higher species richness at equilibrium, but have a nonlinear response to changing extinction rates. These results imply that the network structure of regional food webs can play a significant role in the assembly and dynamics of local webs in natural ecosystems. With natural landscapes becoming increasingly fragmented, understanding such structure may be a necessary key to understanding the maintenance and stability of local species diversity.     

Exploring the temporal variability of a food web using long‐term biomonitoring data

Ecological communities are constantly being reshaped in the face of environmental change and anthropogenic pressures. Yet, how food webs change over time remains poorly understood. Food web science is characterized by a trade‐off between complexity (in terms of the number of species and feeding links) and dynamics. Topological analysis can use complex, highly resolved empirical food web models to explore the architecture of feeding interactions but is limited to a static view, whereas ecosystem models can be dynamic but use highly aggregated food webs. Here, we explore the temporal dynamics of a highly resolved empirical food web over a time period of 18 years, using the German Bight fish and benthic epifauna community as our case study. We relied on long‐term monitoring ecosystem surveys (from 1998 to 2015) to build a metaweb, i.e. the meta food web containing all species recorded over the time span of our study. We then combined time series of species abundances with topological network analysis to construct annual food web snapshots. We developed a new approach, ‘node‐weighted’ food web metrics by including species abundances to represent the temporal dynamics of food web structure, focusing on generality and vulnerability. Our results suggest that structural food web properties change through time; however, binary food web structural properties may not be as temporally variable as the underlying changes in species composition. Further, the node‐weighted metrics enabled us to detect that food web structure was influenced by changes in species composition during the first half of the time series and more strongly by changes in species dominance during the second half. Our results demonstrate how ecosystem surveys can be used to monitor temporal changes in food web structure, which are important ecosystem indicators for building marine management and conservation plans.     

Plant functional traits reveal the relative contribution of habitat and food preferences to the diet of grasshoppers

Food preferences and food availability are two major determinants of the diet of generalist herbivores and of their spatial distribution. How do these factors interact and eventually lead to diet differentiation in co-occurring herbivores? We quantified the diet of four grasshopper species co-occurring in subalpine grasslands using DNA barcoding of the plants contained in the faeces of individuals sampled in the field. The food preferences of each grasshopper species were assessed by a choice (cafeteria) experiment from among 24 plant species common in five grassland plots, in which the four grasshoppers were collected, while the habitat was described by the relative abundance of plant species in the grassland plots. Plant species were characterised by their leaf economics spectrum (LES), quantifying their nutrient vs. structural tissue content. The grasshoppers’ diet, described by the mean LES of the plants eaten, could be explained by their plant preferences but not by the available plants in their habitat. The diet differed significantly across four grasshopper species pairs out of six, which validates food preferences assessed in standardised conditions as indicators for diet partitioning in nature. In contrast, variation of the functional diversity (FD) for LES in the diet was mostly correlated to the FD of the available plants in the habitat, suggesting that diet mixing depends on the environment and is not an intrinsic property of the grasshopper species. This study sheds light on the mechanisms determining the feeding niche of herbivores, showing that food preferences influence niche position whereas habitat diversity affects niche breadth.     

Niche‐dependent trophic position distributions among primary,secondary and tertiary consumers

Many consumers display flexible feeding strategies that vary among individuals or populations, through life‐history, or spatiotemporally. Despite the recognized influence of flexible feeding on the structure and dynamics of food webs, the consequences of these feeding strategies on the actual shape and characteristics of trophic position distributions have received less attention. We proposed and tested several a priori hypotheses to predict the likely effect of niche‐dependent (e.g. herbivore, secondary consumer) foraging on the shape and statistical properties of consumer trophic position distributions using natural abundance stable isotope data from a diverse dataset of consumers. We found evidence that the structural characteristics of consumer trophic position distributions varied as a function of trophic niche. Herbivores and tertiary consumers tended to be ‘packed’ closely near their mean trophic position, with few individuals realizing trophic positions markedly higher or lower than the mean. Conversely, secondary consumers often displayed broad trophic position distributions with many individuals dispersed away from the center of the distribution. We examined the effect of applying constant versus dynamic isotope trophic fractionation models and found that both models yielded similar although not identical results. Our findings suggest that trophic level omnivory supports a larger fraction of consumer diet at intermediate trophic positions than at either the lowest or the highest positions in aquatic food webs. These results suggest that vertical trophic niche declines among higher order consumers despite general evidence that the range of potential foraging options (i.e. horizontal trophic niche) tends to increase at higher trophic positions. Although further work is needed to test the generality of these patterns in other ecosystems, proactively examining trophic position distributions and reporting appropriate measures of central tendency (e.g. arithmetic versus geometric means) will increase the accuracy of individual trophic studies as well as the applicability of results for meta‐analytical food web models.     

Topological properties of food webs: from real data to community assembly models

We explore patterns of trophic connections between species in the largest and highest-quality empirical food webs to date, introducing a new topological property called the link distribution frequency (i.e. degree distribution), defined as the frequency of species S _L with L links. Non-trivial differences are shown in link distribution frequencies between species-rich and species-poor communities, which might have important consequences for the responses of ecosystems to disturbances. Coarse-grained topological properties observed, as species richness-connectance and number of links-species richness relationships, provide no support for the theory of links-species scaling law or constant connectance across empirical food webs investigated. We further explore these observations by means of simulated food webs resulting from multitrophic assembly models using different functional responses between species. Species richness-connectance and links-species richness relationships of empirical food webs are reproduced by our models, but degree distributions are not properly predicted, suggesting the need of new theoretical approximations to food web assembly. The best agreement between empirical and simulated webs occurs for low values of interaction strength between species, corroborating previous empirical and theoretical findings where weak interactions govern food web dynamics.     

Connectance and parasite diet breadth in flea-mammal webs

The number of links in webs of species interactions, which lies at the heart of the biodiversity-stability debate, has given rise to controversy during the last 20 yr. Studies exploring these web properties have mainly focused on symmetric webs where each species can potentially feed on any other species; asymmetric webs such as host-parasite webs, where one set of species feed on another set of species, have been overlooked. However, food webs are incomplete without parasites and the study of parasite-host sub-web properties deserves attention. Here, using a large database involving 33 regional interaction webs between mammals and their flea parasites, we found a negative relationship between species richness and host-parasite connectance. We suggest that some phylogenetic constraints on flea diet may explain our observed patterns because we found that parasite diet breadth, measured as host taxonomic diversity, was invariant along our host richness gradient. We found that the slope of the logarithmic relationship between the number of realized links and species richness is lower than slope values reported for food webs. We suggest that connectance may not respond to increasing species richness as rapidly in host-parasite webs as in predator-prey food webs due to stronger coevolutionary requirements.     

Modelling Size Structured Food Webs Using a Modified Niche Model with Two Predator Traits

The structure of food webs is frequently described using phenomenological stochastic models. A prominent example, the niche model, was found to produce artificial food webs resembling real food webs according to a range of summary statistics. However, the size structure of food webs generated by the niche model and real food webs has not yet been rigorously compared. To fill this void, I use a body mass based version of the niche model and compare prey-predator body mass allometry and predator-prey body mass ratios predicted by the model to empirical data. The results show that the model predicts weaker size structure than observed in many real food webs. I introduce a modified version of the niche model which allows to control the strength of size-dependence of predator-prey links. In this model, optimal prey body mass depends allometrically on predator body mass and on a second trait, such as foraging mode. These empirically motivated extensions of the model allow to represent size structure of real food webs realistically and can be used to generate artificial food webs varying in several aspects of size structure in a controlled way. Hence, by explicitly including the role of species traits, this model provides new opportunities for simulating the consequences of size structure for food web dynamics and stability.     

A modified niche model for generating food webs with stage‐structured consumers: The stabilizing effects of life‐history stages on complex food webs

1. Almost all organisms grow in size during their lifetime and switch diets, trophic positions, and interacting partners as they grow. Such ontogenetic development introduces life‐history stages and flows of biomass between the stages through growth and reproduction. However, current research on complex food webs rarely considers life‐history stages. The few previously proposed methods do not take full advantage of the existing food web structural models that can produce realistic food web topologies.
2. We extended the niche model developed by Williams and Martinez (Nature, 2000, 404, 180–183) to generate food webs that included trophic species with a life‐history stage structure. Our method aggregated trophic species based on niche overlap to form a life‐history structured population; therefore, it largely preserved the topological structure of food webs generated by the niche model. We applied the theory of allometric predator–prey body mass ratio and parameterized an allometric bioenergetic model augmented with biomass flow between stages via growth and reproduction to study the effects of a stage structure on the stability of food webs.
3. When life‐history stages were linked via growth and reproduction, more food webs persisted, and persisting food webs tended to retain more trophic species. Topological differences between persisting linked and unlinked food webs were small to modest. The slopes of biomass spectra were lower, and weak interaction links were more prevalent in the linked food webs than the unlinked ones, suggesting that a life‐history stage structure promotes characteristics that can enhance stability of complex food webs.
4. Our results suggest a positive relationship between the complexity and stability of complex food webs. A life‐history stage structure in food webs may play important roles in dynamics of and diversity in food webs.

     

Carnivory and resource-based niche differentiation in anuran larvae: implications for food web and experimental ecology

1.  Food webs represent the paths of material and energy flow through organisms in an ecosystem. Anuran larvae are important components of pond food webs: they are abundant, consume large quantities of food and serve as prey for many organisms. However, there are very basic uncertainties about the feeding ecology of anuran larvae; for instance, as to which trophic level they belong and whether species differ in resource use. Because anuran larvae have been employed in model systems in experimental ecology for decades, these uncertainties could lead to misinterpretation of published experiments, or inadequate designs of experiments directed at general, conceptual issues in ecology.
2.  Using ¹³C and ¹⁵N stable isotope and gut content analyses of free-ranging and enclosed tadpoles of four ranid species ( Lithobates sylvaticus , L. pipiens , L. clamitans , L. catesbeianus ) in the food webs of six wetlands, we tested the following null hypotheses: (i) that anuran larvae are strict primary consumers; (ii) that they are non-selective feeders and therefore exhibit little feeding niche differentiation; (iii) that they are opportunistic consumers and (iv) that their diet remains unchanged through ontogeny.
3.  All four species consumed and assimilated substantial amounts of animal food; bullfrog larvae, in particular, appear to be predatory. Significant feeding niche differentiation among species occurred with respect to the sources of carbon, consumption of animal matter and nutritional quality of food ingested. We further documented opportunistic feeding habits and ontogenetic shifts in diet.
4.  Collectively, these studies revealed complex trophic relationships that might require a reconsideration of the role of anuran larvae in pond food webs, as well as a reinterpretation of results of previous studies employing anuran larvae in model experimental systems.     

Functional links and robustness in food webs

The robustness of ecosystems to species losses is a central question in ecology, given the current pace of extinctions and the many species threatened by human impacts, including habitat destruction and climate change. Robustness from the perspective of secondary extinctions has been addressed in the context of food webs to consider the complex network of species interactions that underlie responses to perturbations. In-silico removal experiments have examined the structural properties of food webs that enhance or hamper the robustness of ecosystems to species losses, with a focus on the role of hubs, the most connected species. Here we take a different approach and focus on the role of the connections themselves. We show that trophic links can be divided into functional and redundant based on their contribution to robustness. The analysis of empirical webs shows that hubs are not necessarily the most important species as they may hold many redundant links. Furthermore, the fraction of functional connections is high and constant across systems regardless of size and interconnectedness. The main consequence of this scaling pattern is that ecosystem robustness can be considerably reduced by species extinctions even when these do not result in any secondary extinctions. This introduces the possibility of tipping points in the collapse of ecosystems.     

Parasites in food webs: the ultimate missing links

Parasitism is the most common consumer strategy among organisms, yet only recently has there been a call for the inclusion of infectious disease agents in food webs. The value of this effort hinges on whether parasites affect food‐web properties. Increasing evidence suggests that parasites have the potential to uniquely alter food‐web topology in terms of chain length, connectance and robustness. In addition, parasites might affect food‐web stability, interaction strength and energy flow. Food‐web structure also affects infectious disease dynamics because parasites depend on the ecological networks in which they live. Empirically, incorporating parasites into food webs is straightforward. We may start with existing food webs and add parasites as nodes, or we may try to build food webs around systems for which we already have a good understanding of infectious processes. In the future, perhaps researchers will add parasites while they construct food webs. Less clear is how food‐web theory can accommodate parasites. This is a deep and central problem in theoretical biology and applied mathematics. For instance, is representing parasites with complex life cycles as a single node equivalent to representing other species with ontogenetic niche shifts as a single node? Can parasitism fit into fundamental frameworks such as the niche model? Can we integrate infectious disease models into the emerging field of dynamic food‐web modelling? Future progress will benefit from interdisciplinary collaborations between ecologists and infectious disease biologists.     

Structural Degradation in Mediterranean Sea Food Webs: Testing Ecological Hypotheses Using Stochastic and Mass-Balance Modelling

Human-mediated disturbances such as fishing, habitat modification, and pollution have resulted in significant shifts in species composition and abundance in marine ecosystems which translate into degradation of food-web structure. Here, we used a comparative ecological modelling approach and data from two food webs (North-Central Adriatic and South Catalan Sea) and two time periods (mid-late 1970s and 1990s) in the Mediterranean Sea to evaluate how changes in species composition and biomass have affected food-web properties and the extent of ecosystem degradation. We assembled species lists and ecological information for both regions and time periods into stochastic structural and mass-balance food-web models, and compared the outcomes of 22 food-web properties. Our results show strong similarities in structural food-web properties between the North-Central Adriatic and South Catalan Seas indicating similar ecosystem structure and levels of ecological degradation between regions and time periods. In contrast, a comparison with other published marine food webs (Caribbean, Benguela, and US continental shelf) suggested that Mediterranean webs are in an advanced state of ecological degradation. This was reflected by lower trophic height, linkage density, connectance, omnivory, species involved in looping, trophic chain length and fraction of biomass at higher trophic levels, as well as higher generality and fraction of biomass at lower trophic levels. An analysis of robustness to simulated species extinction revealed lower robustness to species removals in Mediterranean webs and corroborated their advanced state of degradation. Importantly, the two modelling approaches used delivered comparable results suggesting that they both capture fundamental information about how food webs are structured. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.