Biodiversity and ecosystem functioning in food webs: the vertical diversity hypothesis

One challenge in merging community and ecosystem ecology is to integrate the complexity of natural multitrophic communities into concepts of ecosystem functioning. Here, we combine food‐web and allometry theories to demonstrate that primary production, as measured by the total nutrient uptake of the multitrophic community, is determined by vertical diversity (i.e. food web's maximum trophic level) and structure (i.e. distributions of species and their abundances and metabolic rates across trophic levels). In natural ecosystems, the community size distribution determines all these vertical patterns and thus the total nutrient uptake. Our model suggests a vertical diversity hypothesis (VDH) for ecosystem functioning in complex food webs. It predicts that, under a given nutrient supply, the total nutrient uptake increases exponentially with the maximum trophic level in the food web and it increases with its maximum body size according to a power law. The VDH highlights the effect of top–down regulation on plant nutrient uptake, which complements traditional paradigms that emphasised the bottom–up effect of nutrient supply on vertical diversity. We conclude that the VDH contributes to a synthetic framework for understanding the relationship between vertical diversity and ecosystem functioning in food webs and predicting the impacts of global changes on multitrophic ecosystems.     

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.     

Modelling the effects of climate change on the distribution and production of marine fishes: accounting for trophic interactions in a dynamic bioclimate envelope model

Climate change has already altered the distribution of marine fishes. Future predictions of fish distributions and catches based on bioclimate envelope models are available, but to date they have not considered interspecific interactions. We address this by combining the species‐based Dynamic Bioclimate Envelope Model (DBEM) with a size‐based trophic model. The new approach provides spatially and temporally resolved predictions of changes in species' size, abundance and catch potential that account for the effects of ecological interactions. Predicted latitudinal shifts are, on average, reduced by 20% when species interactions are incorporated, compared to DBEM predictions, with pelagic species showing the greatest reductions. Goodness‐of‐fit of biomass data from fish stock assessments in the North Atlantic between 1991 and 2003 is improved slightly by including species interactions. The differences between predictions from the two models may be relatively modest because, at the North Atlantic basin scale, (i) predators and competitors may respond to climate change together; (ii) existing parameterization of the DBEM might implicitly incorporate trophic interactions; and/or (iii) trophic interactions might not be the main driver of responses to climate. Future analyses using ecologically explicit models and data will improve understanding of the effects of inter‐specific interactions on responses to climate change, and better inform managers about plausible ecological and fishery consequences of a changing environment.     

Size-based predictions of food web patterns

We employ size-based theoretical arguments to derive simple analytic predictions of ecological patterns and properties of natural communities: size-spectrum exponent, maximum trophic level, and susceptibility to invasive species. The predictions are brought about by assuming that an infinite number of species are continuously distributed on a size–trait axis. It is, however, an open question whether such predictions are valid for a food web with a finite number of species embedded in a network structure. We address this question by comparing the size-based predictions to results from dynamic food web simulations with varying species richness. To this end, we develop a new size- and trait-based food web model that can be simplified into an analytically solvable size-based model. We confirm existing solutions for the size distribution and derive novel predictions for maximum trophic level and invasion resistance. Our results show that the predicted size-spectrum exponent is borne out in the simulated food webs even with few species, albeit with a systematic bias. The predicted maximum trophic level turns out to be an upper limit since simulated food webs may have a lower number of trophic levels, especially for low species richness, due to structural constraints. The size-based model possesses an evolutionary stable state and is therefore un-invadable. In contrast, the food web simulations show that all communities, irrespective of number of species, are equally open to invasions. We use these results to discuss the validity of size-based predictions in the light of the structural constraints imposed by food webs.     

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.     

Biodiversity in model ecosystems, II: species assembly and food web structure

This is the second of two papers dedicated to the relationship between population models of competition and biodiversity. Here, we consider species assembly models where the population dynamics is kept far from fixed points through the continuous introduction of new species, and generalize to such models the coexistence condition derived for systems at the fixed point. The ecological overlap between species and shared preys, that we define here, provides a quantitative measure of the effective interspecies competition and of the trophic network topology. We obtain distributions of the overlap from simulations of a new model based both on immigration and speciation, and show that they are in good agreement with those measured for three large natural food webs. As discussed in the first paper, rapid environmental fluctuations, interacting with the condition for coexistence of competing species, limit the maximal biodiversity that a trophic level can host. This horizontal limitation to biodiversity is here combined with either dissipation of energy or growth of fluctuations, which in our model limit the length of food webs in the vertical direction. These ingredients yield an effective model of food webs that produce a biodiversity profile with a maximum at an intermediate trophic level, in agreement with field studies.     

Dynamics of springtail and mite populations: the role of density dependence, predation, and weather

Abstract 1. Ecological theory suggests that density‐dependent regulation of organism abundance will vary from exogenous to endogenous factors depending on trophic structure. Changes in abundance of soil arthropods were investigated at three trophic levels, springtails (Collembola), predaceous mites (Acari), and macro‐arthropods (spider, adult and larval beetles, centipedes). Predictions were that springtails are predator regulated and mites are food limited according to the Hairston et al. (1960) model, which predicts alternating regulation by competition and predation from fungi to springtails to mites to macro‐arthropods. The alternate hypothesis was based on the bottom‐up model of trophic dynamics, which predicts that each trophic level is regulated by competition for resources. 2. The relative contributions to springtail and mite population dynamics of endogenous (i.e. density‐dependent population growth related to food availability) and exogenous (i.e. predation and weather) factors were tested using time‐series analysis and experimental manipulation of water conditions. Box patterns were distributed within an aspen forest habitat located in the Canadian prairies and surveyed weekly from May to September 1997–1999. Each box depressed the leaf litter, creating a microhabitat island for soil arthropods that provided counts of invertebrates located immediately beneath the boxes. 3. Strong evidence was found for endogenous control of springtail and mite numbers, indicated by a reduction in population growth related to density in the previous week. Contrary to predictions, no evidence was found for regulation of springtail numbers by mites, or for regulation of mite numbers by macro‐arthropods. Springtail population growth rate was related positively to current springtail density (8 and 23% variation explained) and related negatively to 1‐week lagged density (85 and 58%), and related negatively to temperature (5 and 5%) for time‐series data and for experimental addition of water respectively. Mite population growth rate was related positively to current mite density (54%) and temperature (4%), and negatively to 1‐week lagged mite density (20%) and precipitation (6%) for time‐series analysis. For experimental addition of water, mite growth rate was related positively to current mite density (44%) and temperature (5%), and negatively to 1‐week lagged density (11%). Results differed from the Hairston et al. (1960) model predictions but were consistent with a bottom‐up view that springtail and mite populations were regulated intrinsically by competition for food and secondarily by temperature as a function of reproduction.     

Species loss and the structure and functioning of multitrophic aquatic systems

Experiments and theory in single trophic level systems dominate biodiversity and ecosystem functioning research and recent debates. All natural ecosystems contain communities with multiple trophic levels, however, and this can have important effects on ecosystem structure and functioning. Furthermore, many experiments compare assembled communities, rather than examining loss of species directly. We identify three questions around which to organise an investigation of how species loss affects the structure and functioning of multitrophic systems. 1) What is the distribution of species richness among trophic levels; 2) from which trophic levels are species most often lost; and 3) does loss of species from different trophic levels influence ecosystem functioning differently? Our analyses show that: 1) Relatively few high‐quality data are available concerning the distribution of species richness among trophic levels. A new data‐set provides evidence of a decrease in species richness as trophic height increases. 2) Multiple lines of evidence indicate that species are lost from higher trophic levels more frequently than lower trophic levels. 3) A theoretical model suggests that both the structure of food webs (occurrence of omnivory and the distribution of species richness among trophic levels) and the trophic level from which species are lost determines the impact of species loss on ecosystem functioning, which can even vary in the sign of the effect. These results indicate that, at least for aquatic systems, models of single trophic level ecosystems are insufficient for understanding the functional consequences of extinctions. Knowledge is required of food web structure, which species are likely to be lost, and also whether cascading extinctions will occur.     

Generalist predator's niche shifts reveal ecosystem changes in an experimentally fragmented landscape

Habitat fragmentation can alter the trophic structure of communities and environmental conditions, thus driving changes in biodiversity and ecosystem functions. Quantifying niches of generalist predators can reveal how fragmentation alters ecosystems. In a habitat fragmentation experiment, we used stable isotopes of a generalist predator skink to test predictions from spatial theory on trophic structure and to quantify abiotic changes associated with fragmentation among continuous forest, fragments, and matrix habitats. We predicted that in fragments and the matrix, isotopic niches would shift due to decreases in skink trophic positions (δ¹⁵N) from reductions in trophic structure of arthropod food webs and abiotic changes over time (δ¹³C) relative to continuous forest. Contrary to theoretical predictions, we did not find evidence of reductions in trophic structure with fragmentation. In fact, skink δ¹⁵N values were higher in the matrix and fragments than continuous forest, likely due to changes in distributions of a detritivorous prey species. In addition, δ¹³C values in the matrix decreased over years since fragmentation due to abiotic changes associated with matrix tree maturation. We show how isotopic niches are influenced by fragmentation via shifts in biotic and abiotic processes. The potential for either or both spatial and abiotic effects of fragmentation present a challenge for theory to better predict ecological changes in fragmented landscapes.     

In situ warming strengthens trophic cascades in a coastal food web

Global warming may affect most organisms and their interactions. Theory and simple mesocosm experiments suggest that consumer top–down control over primary producer biomass should strengthen with warming, since consumer respiration increases faster with warming than plant photosynthesis. However, these predictions have so far not been tested on natural communities that have experienced warming over many generations. Natural systems display a higher diversity, heterogeneity and complexity than mesocosms, which could alter predicted effects of warming. Here we used an artificially heated part of the northern Baltic Sea (the Forsmark Biotest basin) to test how warming influences trophic interactions in a shallow coastal food web with four trophic levels: omnivorous fish, invertivorous fish, herbivorous invertebrates, and filamentous macroalgae. Monitoring of fish assemblages over six years showed that small invertivorous fish (gobiids, sticklebacks and minnows) were much less abundant in the heated basin than in unheated references areas. Stomach content analyses of the dominating omnivorous fish – Eurasian perch Perca fluviatilis – revealed a strikingly different diet within and outside the Biotest basin; gammarid crustaceans were the dominating prey at heated sites, whereas invertivorous fish (e.g. gobiids) dominated at unheated sites. A 45‐day cage experiment showed that fish exclusion did not affect the biomass of algal herbivores (gastropods and gammarids), but reduced algal biomass in heated sites (but not unheated). This suggests that warming induced a trophic cascade from fish to algae, and that this effect was mediated by predator‐induced changes in herbivore behavior, rather than number. Overall, our study suggests that warming has effectively compressed the food chain from four to three trophic levels (algae, gammarids and perch), which have benefitted the primary producers by reducing grazing pressure. Consequently, warming appears to have restructured this coastal food web through a combination of direct (physiological) and indirect (species interactions) effects.     

Habitat size influences food web structure in drying streams

Biodiversity in running waters is threatened by an increased severity and incidence of low‐flow extremes resulting from global climate change and a growing human demand for freshwater resources. Although it is unknown how and to what extent riverine communities will change in the face of these threats, considerable insight will be gained from efforts aimed at quantifying habitat size‐related controls on the trophic relationships among taxa in streams experiencing extreme flow loss. Here we report on a detailed space‐for‐time survey of replicate stream food webs sampled along the perennial‐ to‐drying continuum in each of fourteen different intermittent South Island, New Zealand streams. We quantified several structural attributes of food webs at fifty‐eight sites, including two taxonomically‐based metrics (web size, predator:prey ratio) and three stable isotope‐based metrics (food chain length [FCL], trophic area, δ¹³C range); we also quantified habitat size‐, disturbance‐, and resource‐related covariates at each site. Food web structure varied widely across sample sites within and across study streams and much of this variation was explained by habitat size. Consistent with our predictions, we found that food webs became smaller (ca 30 to ca 15 taxa, ca 20‐fold reduction in stable isotope‐based trophic area) and shorter (maximum trophic position [FCL] from 4.1 to 2.0, 25% reduction in predator:prey ratio) as we moved from the largest to smaller habitats. These results, and a comparison of our findings with those from a similar assessment conducted in perennial streams, suggest that there are perturbation thresholds which may trigger food web collapse when exceeded, and further imply that food webs may ultimately be ‘sized’ to minimum flows rather than average flow conditions. Our work provides a basis for making general predictions about how habitat contraction, and flow loss in particular, may affect communities and additionally provides insight on mechanisms warranting further attention.     

Temperature effects on fish production across a natural thermal gradient

Global warming is widely predicted to reduce the biomass production of top predators, or even result in species loss. Several exceptions to this expectation have been identified, however, and it is vital that we understand the underlying mechanisms if we are to improve our ability to predict future trends. Here, we used a natural warming experiment in Iceland and quantitative theoretical predictions to investigate the success of brown trout as top predators across a stream temperature gradient (4–25 °C). Brown trout are at the northern limit of their geographic distribution in this system, with ambient stream temperatures below their optimum for maximal growth, and above it in the warmest streams. A five‐month mark‐recapture study revealed that population abundance, biomass, growth rate, and production of trout all increased with stream temperature. We identified two mechanisms that contributed to these responses: (1) trout became more selective in their diet as stream temperature increased, feeding higher in the food web and increasing in trophic position; and (2) trophic transfer through the food web was more efficient in the warmer streams. We found little evidence to support a third potential mechanism: that external subsidies would play a more important role in the diet of trout with increasing stream temperature. Resource availability was also amplified through the trophic levels with warming, as predicted by metabolic theory in nutrient‐replete systems. These results highlight circumstances in which top predators can thrive in warmer environments and contribute to our knowledge of warming impacts on natural communities and ecosystem functioning.     

Experimental evidence for alternative stable equilibria in a benthic pond food web

Both local and initial conditions have been implicated in causing variation in the strengths of predator effects among natural communities. In this study, I performed a controlled mesocosm experiment using two common snails, Physella gyrina and Helisoma trivolvis, and their insect predator, Belostoma flumineum, to explicitly quantify and understand the causes of previously observed variation in food web interactions in natural ecosystems. The results from a controlled mesocosm experiment showed that: (1) at low productivity, snails were dominated by small‐vulnerable species (P. gyrina), predator effects were strong and predator effects caused a trophic cascade extending to producers; (2) at high productivity, snails were dominated by large‐invulnerable species (H. trivolvis), predator effects were weak and there was no trophic cascade to producers; (3) at intermediate productivity, alternative stable equilibria were implicated, such that depending on initial conditions, snails were dominated by either small‐vulnerable or large‐invulnerable species and predator effects were either strong or weak. Thus, this study represents one of the first to quantify how variation in food web structure can result from the local environment, initial conditions and their interaction.     

Predictive power of food web models based on body size decreases with trophic complexity

Food web models parameterised using body size show promise to predict trophic interaction strengths (IS) and abundance dynamics. However, this remains to be rigorously tested in food webs beyond simple trophic modules, where indirect and intraguild interactions could be important and driven by traits other than body size. We systematically varied predator body size, guild composition and richness in microcosm insect webs and compared experimental outcomes with predictions of IS from models with allometrically scaled parameters. Body size was a strong predictor of IS in simple modules (r² = 0.92), but with increasing complexity the predictive power decreased, with model IS being consistently overestimated. We quantify the strength of observed trophic interaction modifications, partition this into density‐mediated vs. behaviour‐mediated indirect effects and show that model shortcomings in predicting IS is related to the size of behaviour‐mediated effects. Our findings encourage development of dynamical food web models explicitly including and exploring indirect mechanisms.     

Predicting the consequences of species loss using size‐structured biodiversity approaches

Understanding the consequences of species loss in complex ecological communities is one of the great challenges in current biodiversity research. For a long time, this topic has been addressed by traditional biodiversity experiments. Most of these approaches treat species as trait‐free, taxonomic units characterizing communities only by species number without accounting for species traits. However, extinctions do not occur at random as there is a clear correlation between extinction risk and species traits. In this review, we assume that large species will be most threatened by extinction and use novel allometric and size‐spectrum concepts that include body mass as a primary species trait at the levels of populations and individuals, respectively, to re‐assess three classic debates on the relationships between biodiversity and (i) food‐web structural complexity, (ii) community dynamic stability, and (iii) ecosystem functioning. Contrasting current expectations, size‐structured approaches suggest that the loss of large species, that typically exploit most resource species, may lead to future food webs that are less interwoven and more structured by chains of interactions and compartments. The disruption of natural body‐mass distributions maintaining food‐web stability may trigger avalanches of secondary extinctions and strong trophic cascades with expected knock‐on effects on the functionality of the ecosystems. Therefore, we argue that it is crucial to take into account body size as a species trait when analysing the consequences of biodiversity loss for natural ecosystems. Applying size‐structured approaches provides an integrative ecological concept that enables a better understanding of each species' unique role across communities and the causes and consequences of biodiversity loss.     

Decoupling of cascading trophic interactions in a freshwater,benthic food chain

Food chain theory provides explicit predictions for equilibrium biomasses among trophic levels in food chains of different lengths. Empirical studies on freshwater benthic food chains have typically been performed on chains with up to three levels and in field experiments with limited spatial and temporal scale. Here we use a natural snapshot experiment approach to study equilibrium biomass and abundance among trophic levels in natural ponds differing only with respect to fish assemblage structure. Forty-four ponds were surveyed for their densityand biomass of fish, snails and periphyton. Ponds were divided into three categories based on fish assemblage: ponds with no fish (two trophic levels), ponds with molluscivorous fish (three trophic levels) and ponds that also had piscivorous fish (four trophic levels). Ponds without fish had a high density and biomass of snails and a low biomass of periphyton, whereas snails were scarce and periphyton biomass was high in ponds with molluscivorous fish. In the presence of piscivores, molluscivore populations consisted of low numbers of large individuals. Snail assemblages in piscivore ponds were characterised by relativelyhigh densities of small-bodied detritivorous species and periphyton biomass was not significantlydifferent from ponds with three trophic levels. Thus, predictions from classic food chain theory were upheld in ponds with up to three trophic levels. In ponds with four trophic levels, however, there was a decoupling of the trophic cascade at the piscivore-molluscivore level. Gape-limited piscivory, predation on snails by molluscivores that have reached an absolute size refuge from predation, and changes in food preferences of the dominant snails are suggested to explain the observed patterns.     

Species–area relationships across four trophic levels – decreasing island size truncates food chains

That larger areas will typically host more diverse ecological assemblages than small ones has been regarded as one of the few fundamental ‘laws’ in ecology. Yet, area may affect not only species diversity, but also the trophic structure of the local ecological assemblage. In this context, recent theory on trophic island biogeography offers two clear‐cut predictions: that the slope of the species–area relationship should increase with trophic rank, and that food chain length (i.e. the number of trophic levels) should increase with area. These predictions have rarely been verified in terrestrial systems. To offer a stringent test of key theory, we focused on local food chains consisting of trophic specialists: plants, lepidopteran herbivores, and their primary and secondary parasitoids. For each of these four trophic levels, we surveyed species richness across a set of 20 off‐shore continental islands spanning a hundred‐fold range in size. We then tested three specific hypotheses: that species richness is affected by island size, that the slope of the species–area curve is related to trophic rank, and that such differences in slope translate into variation in food chain length with island size. Consistent with these predictions, estimates of the species–area slope steepened from plants through herbivores and primary parasitoids to secondary parasitoids. As a result of the elevated sensitivity of top consumers to island size, food chain length decreased from large to small islands. Since island size did not detectably affect the ratio between generalists and specialists among either herbivores (polyphages vs oligophages) or parasitoids (idiobionts vs koinobionts), the patterns observed seemed more reflective of changes in the overall number of nodes and levels in local food webs than of changes in their linking structure. Overall, our results support the trophic‐level hypothesis of island biogeography. Per extension, they suggest that landscape modification may imperil food web integrity and vital biotic interactions.     

Predators,fire or resources: What drives the distribution of herbivores in fragmented mesic forests?

Trophic interactions and disturbance events can shape the structure and function of ecosystems. However, the effects of drivers such as predation, fire and climatic variables on species distributions are rarely considered concurrently. We used a replicated landscape‐scale predator management experiment to compare the effects of red fox Vulpes vulpes control, time‐since‐fire, vegetation type and other environmental variables on native herbivore distributions. Occurrence data for four native herbivores and an invasive predator – the red fox – were collected from 240 sites across three baited (for lethal fox control) and three unbaited forest blocks (4659–9750 ha) in south‐western Victoria, Australia, and used to build species distribution models. The herbivore taxa were as follows: red‐necked wallaby Macropus rufogriseus, black wallaby Wallabia bicolour, grey kangaroo Macropus fuligenosus and Macropus giganteus and common brushtail possum Trichosurus vulpecula. Fox control and fire had little effect on herbivore occurrence, despite the literature suggesting it can influence abundance, while climate, proximity to farmland and topography were more influential. This may be because the region’s high productivity and agricultural pastures subsidise food resources for both predators and prey within the forest blocks and so dampen trophic interactions. Alternatively, these drivers may affect herbivore abundance, but not herbivore occurrence. Understanding the drivers of herbivore distributions is an important step in predicting the effects of herbivory on other species, particularly after management interventions such as predator control and prescribed burns.     

Introduced invertebrates are important prey for a generalist predator

Concern over biological invasions has drawn increased attention to the impacts of introduced predators or competitors, but not to the importance of introduced prey. North American forests are rich in introduced invertebrates, including species that represent relatively novel taxonomic or trophic guilds and show biased distributions among forest types. We analysed the diets of red‐backed salamanders, Plethodon cinereus, from three upland and three lowland forests to determine whether introduced prey are important contributors to geographical or temporal variation in salamander food resources. We found several introduced species were volumetrically important salamander prey, and were responsible for resource differences between forest types and much of the seasonal fluctuation in food resources in both forest types. In lowland forests, rain had a stronger effect on salamander predation on non‐native earthworms than native taxa, creating more dynamic resource fluctuations in resource levels than was observed in upland forests where earthworms were absent. With one exception, predation on non‐native species was positively associated with predation on native species, suggesting non‐native prey have added to salamander resources rather than replaced salamander predation on native taxa. We hypothesize that the novel resource gradients created by non‐native prey introductions are contributing to patterns of geographical and temporal phenotypic variation among salamander populations.