Translator Disclaimer
19 September 2017 Warming Strengthens the Ecological Role of Intraspecific Variation in a Predator
Author Affiliations +

Recent work shows communities and ecosystems can be shaped by predator intraspecific variation, but it is unclear whether the magnitude and direction of these influences are context-dependent. Temperature is an environmental context of strong ecological influence and widespread relevance given global warming trends. Warming should increase per capita predator effects on prey through increases in predator metabolic rate, potentially exacerbating intraspecific differences in ecological effects. Here, we used two populations of the potent pelagic freshwater predator, Western Mosquitofish (Gambusia affinis), to test how experimental pond temperature mediates the differences between their ecological impacts. Mosquitofish introduction induced a strong pelagic trophic cascade, causing a large reduction of crustacean zooplankton biomass, an increase in phytoplankton biomass, and changes to ecosystem-level response variables. Warming ( 2°C above unwarmed treatments) exacerbated fish-induced reduction of zooplankton biomass, but moderated the cascade to phytoplankton, primary productivity, and nutrient concentrations. Effects of intraspecific variation were apparent only on zooplankton, and only at warmed environmental temperatures. The traits underlying this divergence may be related to the population source thermal environments. Overall, results show that warming may increase the ecological importance of predator intraspecific variation. In general, extrinsic environmental drivers, such as those associated with climate change, may reshape the effects of intraspecific trait variation on ecosystems.

INTRASPECIFIC trait variation in ecologically important species is becoming a widely recognized potential driver of community—and ecosystem—level characteristics and processes (reviewed in Whitham et al., 2003; Hairston et al., 2005; Fussmann et al., 2007; Bailey et al., 2009a; Post and Palkovacs, 2009; Bolnick et al., 2011; Matthews et al., 2011; Schoener, 2011). Intraspecific effects can be strong even contrasted with traditional ecological factors like presence of a dominant species (Bailey et al., 2009b; Palkovacs et al., 2015; Gómez et al., 2016) and habitat size (Farkas et al., 2013). A typical study design is a “common gardening” experiment (sensu Matthews et al., 2011), in which one tests the ecological impacts of intraspecific trait variants in a common environmental context (e.g., Schweitzer et al., 2004; Palkovacs and Post, 2009; Ingram et al., 2011; Lundsgaard-Hansen et al., 2014; Fryxell et al., 2015; Rudman and Schluter, 2016). However, the role of environmental context in determining the strength and direction of intraspecific effects is not well known. Some studies find intraspecific effects depend on biotic context (i.e., presence of another species or another species' particular phenotype; Palkovacs et al., 2009; Ingram et al., 2012; Rudman et al., 2015), but few studies evaluate how intraspecific effects may depend on the abiotic environment, especially in animals (but see El-Sabaawi et al., 2015; Lajoie and Vellend, 2015; Tuckett et al., this volume, 2017). Because abiotic context strongly shapes ecological interactions (Chamberlain et al., 2014), impacts of intraspecific variation likely also depend on abiotic context.

Temperature is an abiotic variable with profound impacts across levels of biological organization. The ecological influence of temperature is fundamental in that it shapes organismal metabolism (Gillooly et al., 2001; Brown et al., 2004), which itself may help explain higher-level ecological patterns such as biodiversity and carbon-flow through ecosystems (Allen et al., 2002; Schramski et al., 2015). Temperature is highly variable across space and through time, and has been increasing rapidly on average across the globe in recent history (IPCC, 2014). Because of its pervasive ecological role and immediate relevance, it is important to understand how temperature may mediate the ecological role of intraspecific variation.

Metabolism increases exponentially with temperature within the range of temperatures typically encountered by an organism (Gillooly et al., 2001), so small increases in body temperature can greatly increase metabolic demand. This increased demand must be met by increased ingestion rates at the individual level (Rall et al., 2012), which, in consumers and predators, could contribute to the widely observed warming-induced strengthening of top-down effects on ecosystems (Sanford, 1999; Barton and Schmitz, 2009; Barton et al., 2009; O'Connor et al., 2009; Hoekman, 2010; Harley, 2011; Kratina et al., 2012; Shurin et al., 2012). If feeding-related trait variation occurs among populations, warming-induced increases in per capita feeding rates could increase the ecological effects differences between populations of ectotherm predators.

Here, we test the prediction that predator intraspecific differences have stronger effects for freshwater communities in a warmed versus an unwarmed environment. We test this prediction using Western Mosquitofish, Gambusia affinis, which prey heavily upon crustacean zooplankton in the pelagic zone of ponds (Hurlbert and Mulla, 1981; Pyke, 2005), causing trophic cascades whereby producer biomass increases, primary productivity increases, and nutrient concentrations decline (Hurlbert et al., 1972; Fryxell et al., 2016). Mosquitofish have been spread globally (Pyke, 2008) and today inhabit a wide diversity of environments to which they have acclimated and adapted (Pyke, 2005). There is considerable trait variation within and among Mosquitofish populations. Body size variation and sex ratio variation are common and can mediate a population's ecological effects. An increasing proportion of females, which are generally larger, can induce stronger trophic cascades (Fryxell et al., 2015). Morphological differences among Mosquitofish populations can emerge via evolutionary responses to predation pressure (Langerhans et al., 2004). Mosquitofish are also known to exhibit rapid evolution of life history traits in response to habitat size variation (Stearns, 1983) and temperature (Stockwell and Weeks, 1999). Such contemporary trait change might also have community and ecosystem effects.

In this experiment, we used wild-caught Mosquitofish from two recently divergent populations of different thermal environment. We crossed three fish treatments (fishless, cool-source, warm-source) with two ecosystem temperature treatments (unwarmed, warmed) to test our predictions that 1) warming exacerbates top-down effects of fish introduction and 2) warming exacerbates the ecological differences between warm- and cool-source fish. Specifically, we expect 1) fish will more strongly suppress zooplankton at warmed versus unwarmed temperatures, which should cascade to affect phytoplankton, productivity, and nutrients, and 2) the ecological differences between populations for these same response variables will be larger at warmed versus unwarmed temperatures. Our use of recently divergent wild-caught fish from populations of different source temperatures may additionally allow us to address how predator trait variation generated along the temperature axis interacts with thermal context to shape ecological conditions.


Source populations.—We used Mosquitofish from two geothermal sites near Bishop, California, USA. The “warm-source” site was Keough's Hot Ditch (37°15′33.6″N, 118°22′18.5″W), which has a mean temperature of 31.6°C, and the “cool-source” site was an unnamed dammed artesian well (37°21′02.1″N, 118°19′35.7″W), which has a mean temperature 23.7°C, as measured over the same period (Fig. 1). Mosquitofish were introduced to California from a common source population in Texas, USA in 1922 (Lenert, 1923). Fish have occupied these particular sites since at least 1980, though gene flow could have occurred until as recently as 2001, when the translocation of Mosquitofish among sites was discontinued by the Owens Valley Mosquito Abatement Program (Bob Kennedy, pers. comm.). There is little potential for natural gene flow between these sites because they are aquatic islands in a desert landscape and are both dammed pools upstream from separate geothermal tributaries of the Owens River watershed.

Fig. 1. 

The temperature profiles of both fish sources used in the study. “Cool-source” fish were taken from the artesian well and “warm-source” fish were taken from Keough's Hot Ditch, both in Bishop, California, USA. Daily average air temperatures are plotted from the Bishop Airport weather station, Bishop, California (Station ID GHCND:USW00023157), measured over the same period. In these geothermal systems, air temperature has little influence on water temperature.


Mesocosm setup.—We established 36 experimental ponds (1136 L Rubbermaid® stock tanks) in a 6 × 6 array at Long Marine Laboratory in Santa Cruz, California, USA. On 5 August 2013, we filled tanks with city water, evenly spread 19 L of sand each across their bottoms, and placed two cinderblocks (9.3 × 19.0 × 39.3 cm3) adjacently in the center of each. We placed a smaller cinderblock (6.0 × 8.5 × 16.0 cm3) between the two larger cinderblocks in the center of all ponds as cover and habitat structure for Mosquitofish. On 20 August 2013, we established a biological community in each pond by introducing a 4 L aliquot of a homogenized sediment sample and even aliquots of a large zooplankton sample taken from Antonelli Pond (36°57′18.6″N, 122°03′37.8″W), Santa Cruz, California, USA. The zooplankton community was dominated by crustacean zooplankton Daphnia, Ceriodaphnia, Bosmina, calanoid copepods, and cyclopoid copepods.

Design and treatment application.—We assigned treatments in our 3 × 2 factorial design (fishless, cool-source, warm-source x unwarmed, warmed) to the 6 × 6 mesocosm array using latin-squares random assignment. We initiated warming on 12 August 2013 using 300-watt heaters (Fluval Aquatics®) placed on the small central cinderblock in ponds assigned the warmed treatment. Similar warming methods have been used in other experiments, and had maintained a ∼3°C temperature above unwarmed treatments, with similar patterns of temperature variation between treatments (Kratina et al., 2012; Shurin et al., 2012). This temperature increase is within the range of those expected globally over the next century (IPCC, 2014).

We collected experimental fish from both sites using seine and hand nets on 22 August 2013. We transported fish to Santa Cruz, California and immediately introduced a random subsample of ten individuals to each experimental pond following treatment assignments. This density of Mosquitofish was equivalent to densities used in prior mesocosm experiments, which reflect approximate natural summertime densities (Fryxell et al., 2015, 2016).

Mesocosm sampling.—We took repeated samples for water temperature, crustacean zooplankton (hereafter zooplankton) biomass, and phytoplankton biomass. Water temperature was monitored continuously at 15-minute intervals using HOBO Pendant (Onset Corporation®) data loggers. We sampled mesocosm water for zooplankton and phytoplankton two and four weeks following fish introduction. We collected zooplankton from a 1 L water sample using a 40 μm sieve and preserved them in ethanol. Later, each zooplankter was identified to the lowest taxonomic group possible at 100X magnification and measured to get the body length using the “Live measurements” module of Leica Microsystems® Application Suite. Lengths were converted to mass by using published length-weight relationships for each taxonomic group (Bottrell et al., 1976), and masses of all zooplankters within a mesocosm sample were summed to get biomass. We collected phytoplankton from 1 L water samples filtered through Whatman GF/F® filters (pore size 0.7 μm). Filters were placed into a –20°C freezer for later analysis. Chlorophyll a was extracted with acetone for 24 hours at 4°C and its concentration analyzed using the non-acid module on a Trilogy Flourometer (Turner Designs®). Chlorophyll a concentration was used as a proxy for phytoplankton abundance.

We took final time point (week 4) samples for net primary productivity (NPP) and phosphate concentrations. NPP was estimated as the change in dissolved oxygen concentration between dusk and the previous dawn (as in Harmon et al., 2009). Phosphate (soluble reactive phosphorus) was measured from 60 mL water samples filtered (0.7 μm pore size), frozen, and then thawed and analyzed on a LaChat® flow injection analyzer following standard methods (Rice et al., 2012).

After ecological sampling was complete, we removed fish and euthanized them with an overdose of tricaine methanesulfonate. Thereafter, fish were measured for length, counted to get population growth (number of offspring), weighed to get total fish biomass per mesocosm, and sexed to get the adult sex ratio (proportion males).

Analyses.—We used ANOVA to test the significance of warming, fish introduction, and intraspecific variation for each response variable. The “Warming” effect is the difference between the ‘warmed' and ‘unwarmed' treatments. Separate ANOVA analyses were used for the independent contrasts of interest—the “fish introduction effect” (FI) and the “intraspecific variation effect” (IV). The FI effect used all treatments but coded both populations as a single “fish present” factor level. The IV effect used only the treatments with fish to test for differences between the effects of the populations. Specifically, we used MANOVA of repeated measures for zooplankton and phytoplankton responses to account for correlations between sampling points, and used ANOVA for NPP and phosphate. All ecological responses were log10 transformed before analyses to improve adherence to ANOVA assumptions, though zooplankton biomass was log10(1+x) transformed to avoid losing a replicate due to a 0 value. We used Bartlett's tests to ensure approximate equality of variances for the FI effect, since it was an unbalanced comparison (6 replicates without fish, 12 with fish). For zooplankton, we additionally performed post hoc MANOVA tests for warmed and unwarmed treatments separately. We used principal components analysis (PCA) on correlations of the four transformed ecological responses at the final sampling point to understand how treatments arranged themselves in multivariate space. PCAs were run on data from all treatments, and then for warmed and unwarmed treatments separately.

We calculated standardized effect sizes (Hedge's G) for the effects of FI and IV separately for unwarmed and warmed treatments to compare the magnitude of effects between community-level responses variables. Effect sizes were calculated across all data (i.e., including data from both time-points). Effect sizes and confidence intervals were calculated using the “cohen.d” function with arguments “hedges.correction=T” and “pooled=T” in the package “effsize” (Torchiano, 2016) on the R platform (R Core Team, 2015).

To associate population traits with effects, we performed a number of tests. First, we used a t-test to test for population differences in sex ratio at the start of the experiment. Next, we tested how temperature and source population affected final fish biomass and population growth using ANOVA. Lastly, we tested how initial sex ratio, final fish biomass, and population growth may have affected zooplankton biomass using standard least-squares linear regressions, separately for unwarmed and warmed ponds. All ANOVA- and regression-related analyses, and PCAs were performed in JMP Pro 12 (SAS Institute®).


Warmed treatments (overall mean 25.0°C) were an average 2.21°C warmer than unwarmed treatments (overall mean 22.8°C) and tracked unwarmed diel temperature patterns closely (Fig. 2A). The treatment difference did vary with time, from as little as 1.30°C to as much as 3.60°C, over the course of the experiment (Fig. 2B). Diel patterns in the magnitude of warming also emerged, where treatment warming was strongest in morning as temperatures naturally rose, and treatment warming was weakest in the evenings as temperatures naturally cooled (Fig. 2).

Fig. 2. 

(A) The mean of each logged temperature for both the warmed (black) and unwarmed (gray) treatments. (B) The difference between treatment means at each time-point, with the hatched line showing the overall mean temperature increase in the warming treatment: 2.21°C.


Fish introduction (FI) effects were significant at community and ecosystem levels, though there was little evidence for a FI x Warming interaction in most cases. Bartlett's test identified no violations (P > 0.01) of the heteroscedasticity assumption of ANOVA (Supplemental Appendix 1; see Data Accessibility), allowing us to proceed with FI tests. FI and Warming reduced zooplankton biomass (FI: F1,32 = 2.70, P < 0.0001; Warming: F1,32 = 7.03, P = 0.0124; Fig. 3A, B), without a FI x Warming interaction (F1,32 = 1.06, P = 0.3104; Fig. 3A, B). Univariate tests revealed a reduction in the FI effect through time (Time x FI: F1,32 = 8.42, P = 0.0067), which was dominated by the unwarmed treatment (post hoc unwarmed only Time x FI: F1,16 = 9.18, P = 0.008; post hoc warmed only Time x FI: F1,16 = 1.08, P = 0.3143; Fig. 3A, B), though the three-way interaction was nonsignificant (Time x FI x Warming: F1,32 = 2.13, P = 0.1540). We found evidence of a trophic cascade associated with FI, as phytoplankton biomass greatly increased with fish introduction (F1,32 = 38.04, P < 0.0001; Fig. 3C, D). Warming increased phytoplankton (F1,32 = 6.06, P = 0.0194; Fig. 3C, D), without an interaction with the FI effect (F1,32 = 0.9645, P = 0.3334; Fig. 3C, D), despite the result that the FI effect was weaker on average at warmer temperatures (Fig. 3C, D). NPP increased as a result of both FI and Warming, though Warming weakened FI effects (FI x Warming: F1,32 = 7.56, P = 0.0097; Fig. 4A). Lastly, FI reduced phosphate concentrations (F1,32 = 10.2826, P = 0.0030; Fig. 4B), and Warming may have moderated this impact (Warming: F1,32 = 3.4948, P = 0.0707; FI x Warming: F1,32 = 2.7594, P = 0.1065; Fig. 4B).

Fig. 3. 

Responses of zooplankton (A, B) and phytoplankton (C, D) to warming and fish treatments for both sampling points. Y-axes are plotted on the log10 scale.


Fig. 4. 

Response of net primary production (A) and phosphate concentration (B) to warming and fish treatments at the final sampling point of the experiment—four weeks after fish introduction. Y-axes are plotted on the log10 scale.


Zooplankton biomass was affected by intraspecific variation (IV), though this effect only emerged at warmed temperatures (IV: F1,20 = 2.92, P = 0.1031; post hoc unwarmed only: F1,10 = 0.2634, P = 0.6189; post hoc warmed only: F1,10 = 6.04, P = 0.0338; Fig. 3A, B). We did not detect cascading impacts of this effect, as IV effects were nonsignificant for phytoplankton biomass (IV x Warming: F1,20 = 0.03, P = 0.8571; IV: F1,20 = 0.001, P = 0.9717; Fig. 3C, D), NPP (IV x Warming: F1,20 = 0.01, P = 0.9067; IV: F1,20 = 0.04, P = 0.8531; Fig. 4A, B), and phosphate concentration (IV x Warming: F1,20 = 0.22, P = 0.6406; IV: F1,20 = 0.42, P = 0.5244; Fig. 4C, D).

The first principal axes of the three PCAs (across all mesocosms, or for warmed and unwarmed mesocosms separately) each explained ∼60% of the variation in the data, while the second principal axes explained ∼20% of the variation. Phytoplankton and NPP had positive loadings on the first axis, while phosphate and zooplankton had negative loadings on the first axis. The second principal axes showed variable loadings across PCAs, suggesting different relationships among responses due to warming. PCAs showed spatial separation of fish from fishless treatments, but significant overlap of the two treatments with fish (Supplemental Appendix 2; see Data Accessibility).

Effect size calculations showed very strong (>>1) effects of FI and little mediation of FI and IV effects by warming (Fig. 5A, B). For the IV effect, 90 percent confidence intervals overlap with zero except for effects on zooplankton in warmed environments (Fig. 5B). Effect sizes generally declined down trophic levels from zooplankton to phytoplankton.

Fig. 5. 

The absolute value of effect sizes for all responses calculated over all sampling points, separately for (A) the fish introduction effect and (B) the intraspecific variation effect.


Fish survival was high. Only 2 of the 240 individuals that were stocked at the start of the experiment died over the course of the experiment. Because of this high survival, and because few fish offspring reached maturity by the end of the experiment, we estimated initial sex ratio as the final adult sex ratio in experimental ponds. The cool-source population (M = 0.37, SD = 0.16) had a significantly higher proportion males than the warm-source population (M = 0.21, SD = 0.16; t22 = 2.34, P = 0.0285), though a range of 0.0 to 0.6 proportion males was present across mesocosms of both populations. Despite this, the cool-source population tended to have higher population growth (F1,18 = 9.90, P = 0.0056, outliers removed; Fig. 6A), though this was apparent only after removal of two large, influential outlier values, or a log transformation (Supplemental Appendix 3; see Data Accessibility). Both populations had similar final fish biomasses (F1,20 = 0.78, P = 0.388; Fig. 6B). Warming increased fish biomass (F1,20 = 4.53, P = 0.0458; Fig. 6B) and population growth (F1,18 = 8.06, P = 0.0109; Fig. 6A), with no interaction with source population (F1,20 = 0.20, P = 0.6599; F1,18 = 0.06, P = 0.9638; Supplemental Appendix 3; see Data Accessibility). Lastly, our linear regressions of zooplankton biomass on sex ratio, fish biomass, and population growth found no significant effects on zooplankton biomass (P > 0.10; Supplemental Appendix 1; see Data Accessibility). However, the variance in final zooplankton biomass explained by sex ratio (unwarmed: R2 = 0.01; warmed: R2 = 0.03), fish biomass (unwarmed: R2 = 0.04; warmed: R2 = 0.16) and population growth (unwarmed: R2 = 0.02; warmed: R2 = 0.14) was greater in warmed than unwarmed treatments. A summary of all regression and ANOVA-related statistical tests is available in Supplemental Appendix 1 (see Data Accessibility).

Fig. 6. 

Trait differences between populations at the end of the experiment. Fish population growth (A) was the number of offspring present at the end of the experiment, and total fish biomass (B) was the dry mass of all fish within a mesocosm at the end of the experiment. Means and standard errors were calculated without removing outliers (Supplemental Appendix 3; see Data Accessibility). Raw data points each reflect a mesocosm (triangles are warm-source, and circles are cool-source).



Warming mediates fish effects.—Mosquitofish have strong top-down impacts in pond ecosystems and show considerable trait variation within and among populations across their global range, making them good candidates for tests of the ecological effects of intraspecific trait variation (Pyke, 2005). As predicted, Mosquitofish introduction reduced crustacean zooplankton biomass, which induced a trophic cascade whereby phytoplankton bloomed, productivity spiked, and nutrient concentrations declined (Figs. 3, 4). These cascading impacts are consistent with prior studies of Mosquitofish-driven trophic cascades (Hurlbert et al., 1972; Hurlbert and Mulla, 1981; Fryxell et al., 2015, 2016) and set the stage for effects of intraspecific variation and their interaction with warming.

We predicted warming would increase the strength of the top-down effects of fish introduction. For zooplankton, this prediction was supported by the final time-point, despite the overall effect of fish introduction having no significant interaction with warming (Fig. 3A, B). This time-dependence was driven by changes in unwarmed treatments only, where fishless and fish present treatments converged through time. In warmed treatments, fish effects were stable. This difference may reflect accelerated biological and ecological processes at warmed temperatures (Brown et al., 2004). Beyond four weeks, it is unclear whether zooplankton biomasses would have further converged in unwarmed treatments.

Despite warming-induced increases in fish introduction effects on zooplankton biomass, we did not find that warming increased pelagic trophic cascades from fish to phytoplankton, because warming directly increased phytoplankton biomass (Fig. 3C, D). This result is contrary to another warming study in pelagic freshwater systems, which showed strengthened trophic cascades from fish to phytoplankton under warming, and hence, enhanced top-down control (Kratina et al., 2012). However, in that experiment, warming did not enhance zooplankton biomass reduction by fish (Shurin et al., 2012). Perhaps fish excretion or zooplankton body size were primary drivers of the observed warming-induced trophic cascade in that study. Warming generally increases excretion rates (Vanni and McIntyre, 2016), which can increase nutrient supply for phytoplankton growth (Vanni and Layne, 1997). Zooplankton community size structure mediates trophic cascades through relationships between body size and grazing rates (Carpenter et al., 1985). In our study, increased fish excretion would not explain warming-induced phytoplankton blooms, as these effects were primarily seen in fishless treatments (Fig. 3C, D). Simultaneous warming-induced increases in phytoplankton and decreases in zooplankton in fishless ecosystems could occur if warming favored defended or toxic phytoplankton phenotypes (e.g., small phytoplankton, gelatinous-sheathed phytoplankton, cyanobacteria). These phenotypes may be favored via increased per capita predation by zooplankton (Peter and Sommer, 2012) and via some combination of nutrients (Vanni, 1987) and temperature (Rigosi et al., 2014; Yvon-Durocher et al., 2015). Without phytoplankton trait and/or community composition data, we are unable to directly address this hypothesis.

Increased productivity is a long-hypothesized (Carpenter et al., 1992) and common response to warming (Yvon-Durocher et al., 2015) and fish introduction (Carpenter et al., 1985). Ecosystem-level primary productivity is a process driven by metabolic rates of producers and producer standing stock, and is ultimately fueled by light and nutrient availability. As in similar studies (Shurin et al., 2012), we found warming and fish introduction increased NPP and decreased phosphate (Fig. 4). However, warming slightly reduced the effect of fish on these responses, perhaps because trophic cascades from fish to phytoplankton biomass were weaker with warming. Overall, we found support for our prediction that warming increases top-down effects of fish on zooplankton, but without cascading impacts for phytoplankton biomass and ecosystem characteristics. Bottom-up factors appear to have moderated this impact across trophic levels and to higher levels of organization.

Warming mediates intraspecific effects.—We predicted that ecological effects of intraspecific variation would be stronger with warming if top-down effects were stronger with warming. We found intraspecific variation had stronger effects at warmed versus unwarmed temperatures for crustacean zooplankton, but not other responses. Bottom-up factors at warmed temperatures likely prevented cascading effects of intraspecific variation, as they did for fish introduction. However, the effect size of intraspecific variation was considerably smaller than fish introduction (Fig. 5B), so our power to detect its cascading influences was also smaller (Cohen, 1992). Even at warmed temperatures, the effect of intraspecific variation was less than half that of fish introduction (Fig. 5B), which is a notably weaker effect than found in many other studies (Palkovacs et al., 2015; Hendry, 2017).

The intraspecific effect from our study could be relatively weak for many reasons. Our focal populations are recently divergent (∼90 years maximum, likely less), with human-assisted gene flow possible until very recently. Many other studies of intraspecific effects in fishes compare traits variants with a much longer or larger divergence—e.g., Alewives at >300 years (Palkovacs et al., 2008), stickleback, which are post-Pleistocene incipient species (Rundle et al., 2000), and guppies, which are probably longer-derived, showing considerable within-drainage genetic divergence and incomplete phenotypic transformation after a 50-year-long transplant experiment (Willing et al., 2010). Second, populations of many fish species exhibit countergradient variation with respect to temperature—where natural selection favors phenotypic similarity across populations in different environments (Conover and Schultz, 1995; Angilletta, 2009). Countergradient variation may have minimized trait differences between our focal populations and minimized ecological effects differences despite potential genetic differences (Kinnison et al., 2015). Third, other experiments are designed with a more thorough understanding of the trait differences between populations (e.g., Alewives, guppies, stickleback), allowing implementation of methods aimed specifically at detecting or maximizing the ecological effects of the functional trait differences. Investigators might choose maximally divergent focal populations for ecological effects tests, might focus on contexts (i.e., experimental arenas) hypothesized to maximize effects, and/or might control for many aspects of ecologically relevant natural variation like sex ratio, body size, and density, so as not to “swamp out” intraspecific effects. Lastly, studies with incomplete knowledge of ecological effects of focal population differences tend to measure a large number of ecological responses, and may emphasize responses showing relatively large intraspecific effects.

The effect sizes we measure in this study may more accurately reflect general intraspecific effect sizes among populations in nature. With little understanding of specific trait differences between our focal populations, we were unable to cater our experimental arena for maximizing intraspecific effects. We chose to control only for initial fish density, so other sources of natural intraspecific variation were present. Lastly, we focused on measuring a few responses with well-known linkages and broad ecological relevance (i.e., trophic cascades). Using such a less biased strategy in picking focal populations and focusing on one or a few ecological mechanisms of general relevance may better inform as to how strong intraspecific effects are more generally in nature.

Traits underlying intraspecific effects.—After testing whether the populations diverged in ecological effects, we asked whether we could predict those effects from traits and attributes of known ecological relevance across experimental populations (i.e., mesocosms)—sex ratio, population growth (a reflection of life history and phenological differences), and final fish biomass. Temperature is known to have strong phenotypic effects on Mosquitofish traits like body size and life history (e.g., Vondracek et al., 1988; Meffe, 1991; Meffe et al., 1995; Stockwell and Weeks, 1999; Stockwell and Vinyard, 2000), and may affect sex ratio via sex-specific responses to temperature (e.g., see Seebacher et al., 2015). In our experiment, cool-source fish ended with the same biomass of fish as warm-source populations, but had more male-biased initial sex ratios and had greater population growth overall (Fig. 6). Therefore, at a given temperature treatment, cool-source populations had the same biomass of fish but had more and smaller individuals. Nevertheless, initial sex ratio, population growth, and biomass were each independently not associated with final zooplankton biomass at either warmed or unwarmed temperatures. Despite these factors being nonsignificant predictors of zooplankton biomass, we did find stronger associations (i.e., explained variance) between them at warmed versus unwarmed temperatures, in line with our overall prediction.

We hypothesized that intraspecific effects would be greater at warmed versus unwarmed temperatures because of increased metabolic demand. This general prediction is an over-simplification. Metabolic demand does not increase exponentially over the range of all possible temperatures, and does not do so with the same rate across evolutionary lineages (Gillooly et al., 2001). Further, feedings rates and other ecological rates do not perfectly track this metabolic demand (Rall et al., 2012), usually increasing approximately linearly up until a performance maximum, beyond which performance declines (Angilletta, 2009). In this study, experimental temperatures were approximately at or below population source temperatures. Theory therefore predicts focal populations were tested within the range that performance should increase with increasing temperature. Cool-source fish were nearer the source temperature of their natural pond (23.7°C) in experimental treatments (unwarmed: 22.8°C; warmed: 25.0°C) than warm-source fish (31.6°C), and greater reduced zooplankton (Fig. 3A, B). Cool-source fish also greater suppressed zooplankton at warmed versus unwarmed temperatures, perhaps because organisms are commonly found to perform optimally at a temperature slightly above that to which they are acclimated or adapted (Angilletta, 2009). Intraspecific variation in thermal performance curves likely contributes to ecological effects differences at different temperatures.

Eco-evo dynamics and global change.—“Eco-evolutionary dynamics” is an emerging worldview that posits evolutionary trait change and ecological change can interact on contemporary timescales (Hendry, 2017). Trait responses to anthropogenic drivers can be particularly fast (Hendry et al., 2008), so eco-evolutionary dynamics may be important in the face of anthropogenic change (Lavergne et al., 2010; Norberg et al., 2012; Urban et al., 2016). A question recently posed by researchers is whether global change drives eco-evolutionary dynamics such that population trait responses buffer or exacerbate further ecological or environmental change (Palkovacs et al., 2012). Tuckett et al. (this volume, 2017) found that lake eutrophication causes trait changes in white perch that exacerbate nutrient loading. Urban et al. (this volume, 2017) found that adaptive responses of anuran prey to climate-induced predator range expansion mitigated the predator's impact on prey survival. Trait responses may thus exacerbate or mitigate further change in the face of global change drivers. Our study focused on another potential global change driver of eco-evolutionary dynamics—temperature. We found the population from temperatures more similar to experimental conditions more strongly affected prey community biomass. If trait differences reflected recent adaptation, these results may suggest that predator thermal adaptation can mitigate ecological change by maintaining top-down effects on prey across different temperatures. However, our study only shows how a presumably relatively adapted versus a relatively non-adapted population might impact prey in one temperature range (i.e., both experimental temperatures were near the “home” temperature of one of the populations). Future studies might profitably test whether predator populations recently adapted to different temperatures have ecological effects more similar in their specific respective “home” temperatures than in their “away” temperatures.


Supplemental material is available at


We thank A. Hendry and D. Green for organizing the JMIH session on eco-evolutionary dynamics and organizing and editing these symposium proceedings. We thank T. Apgar, S. Parmenter, C. Wickham (and the Owens Valley Mosquito Abatement Program), K. Roessler, E. Argo, D. Hasselman, M. Sabal, and C. Ray for help in the lab and the field. P. Thompson and P. Kratina assisted with implementation of warming methods. D.C.F. was supported by a NSF GRF. E.P.P. was supported by NSF-DEB #1457333. Fish collections were completed under California Department of Fish and Wildlife Permit 12726. Animals were handled ethically following UC Santa Cruz IACUC protocol Palke1306-2.



Allen, A. P., J. H. Brown, and J. F. Gillooly. 2002. Global biodiversity, biochemical kinetics, and the energetic-equivalence rule. Science 297:1545–1548. Google Scholar


Angilletta, M. J., Jr. 2009. Thermal Adaptation: A Theoretical and Empirical Synthesis. Oxford University Press, New York. Google Scholar


Bailey, J. K., A. P. Hendry, M. T. Kinnison, D. M. Post, E. P. Palkovacs, F. Pelletier, L. J. Harmon, and J. A. Schweitzer. 2009a. From genes to ecosystems: an emerging synthesis of eco-evolutionary dynamics. New Phytologist 184:746–749. Google Scholar


Bailey, J. K., J. A. Schweitzer, F. Úbeda, J. Koricheva, C. J. LeRoy, M. D. Madritch, B. J. Rehill, R. K. Bangert, D. G. Fischer, G. J. Allan, and T. G. Whitham. 2009b. From genes to ecosystems: a synthesis of the effects of plant genetic factors across levels of organization. Philosophical Transactions of the Royal Society B: Biological Sciences 364:1607–1616. Google Scholar


Barton, B. T., A. P. Beckerman, and O. J. Schmitz. 2009. Climate warming strengthens indirect interactions in an old-field food web. Ecology 90:2346–2351. Google Scholar


Barton, B. T., and O. J. Schmitz. 2009. Experimental warming transforms multiple predator effects in a grassland food web. Ecology Letters 12:1317–1325. Google Scholar


Bolnick, D. I., P. Amarasekare, M. S. Araújo, R. Bürger, J. M. Levine, M. Novak, V. H. W. Rudolf, S. J. Schreiber, M. C. Urban, and D. A. Vasseur. 2011. Why intraspecific trait variation matters in community ecology. Trends in Ecology & Evolution 26:183–192. Google Scholar


Bottrell, H. H., A. Duncan, Z. M. Gliwicz, E. Grygierek, A. Herzig, A. Hillbricht-Ilkowska, H. Kurasawa, P. Larsson, and T. Weglenska. 1976. Review of some problems in zooplankton production studies. Norwegian Journal of Zoology 24:419–456. Google Scholar


Brown, J. H., J. F. Gillooly, A. P. Allen, V. M. Savage, and G. B. West. 2004. Toward a metabolic theory of ecology. Ecology 85:1771–1789. Google Scholar


Carpenter, S. R., S. G. Fisher, N. B. Grimm, and J. F. Kitchell. 1992. Global change and freshwater ecosystems. Annual Review of Ecology and Systematics 23:119–139. Google Scholar


Carpenter, S. R., J. F. Kitchell, and J. R. Hodgson. 1985. Cascading trophic interactions and lake productivity. BioScience 35:634–639. Google Scholar


Chamberlain, S. A., J. L. Bronstein, and J. A. Rudgers. 2014. How context dependent are species interactions? Ecology Letters 17:881–890. Google Scholar


Cohen, J. 1992. A power primer. Psychological Bulletin 112:155–159. Google Scholar


Conover, D. O., and E. T. Schultz. 1995. Phenotypic similarity and the evolutionary significance of countergradient variation. Trends in Ecology & Evolution 10:248–252. Google Scholar


El-Sabaawi, R. W., R. D. Bassar, C. Rakowski, M. C. Marshall, B. L. Bryan, S. N. Thomas, C. Pringle, D. N. Reznick, and A. S. Flecker. 2015. Intraspecific phenotypic differences in fish affect ecosystem processes as much as bottom-up factors. Oikos 124:1181–1191. Google Scholar


Farkas, T. E., T. Mononen, A. A. Comeault, I. Hanski, and P. Nosil. 2013. Evolution of camouflage drives rapid ecological change in an insect community. Current Biology 23:1835–1843. Google Scholar


Fryxell, D. C., H. A. Arnett, T. M. Apgar, M. T. Kinnison, and E. P. Palkovacs. 2015. Sex ratio variation shapes the ecological effects of a globally introduced freshwater fish. Proceedings of the Royal Society B: Biological Sciences 282:20151970. Google Scholar


Fryxell, D. C., A. R. Diluzio, M. A. Friedman, N. A. Menge, and E. P. Palkovacs. 2016. Cross-habitat effects shape the ecosystem consequences of co-invasion by a pelagic and a benthic consumer. Oecologia 182:519–528. Google Scholar


Fussmann, G. F., M. Loreau, and P. A. Abrams. 2007. Eco-evolutionary dynamics of communities and ecosystems. Functional Ecology 21:465–477. Google Scholar


Gillooly, J. F., J. H. Brown, G. B. West, V. M. Savage, and E. L. Charnov. 2001. Effects of size and temperature on metabolic rate. Science 293:2248–2251. Google Scholar


Gómez, P., S. Paterson, L. De Meester, X. Liu, L. Lenzi, M. D. Sharma, K. McElroy, and A. Buckling. 2016. Local adaptation of a bacterium is as important as its presence in structuring a natural microbial community. Nature Communications 7:12453. Google Scholar


Hairston, N. G., S. P. Ellner, M. A. Geber, T. Yoshida, and J. A. Fox. 2005. Rapid evolution and the convergence of ecological and evolutionary time. Ecology Letters 8:1114–1127. Google Scholar


Harley, C. D. G. 2011. Climate change, keystone predation, and biodiversity loss. Science 334:1124–1127. Google Scholar


Harmon, L. J., B. Matthews, S. Des Roches, J. M. Chase, J. B. Shurin, and D. Schluter. 2009. Evolutionary diversification in stickleback affects ecosystem functioning. Nature 458:1167–1170. Google Scholar


Hendry, A. 2017. Eco-evolutionary Dynamics. Princeton University Press, Princeton, New Jersey. Google Scholar


Hendry, A. P., T. J. Farrugia, and M. T. Kinnison. 2008. Human influences on rates of phenotypic change in wild animal populations. Molecular Ecology 17:20–29. Google Scholar


Hoekman, D. 2010. Turning up the heat: temperature influences the relative importance of top-down and bottom-up effects. Ecology 91:2819–2825. Google Scholar


Hurlbert, S. H., and M. S. Mulla. 1981. Impacts of mosquitofish (Gambusia affinis) predation on plankton communities. Hydrobiologia 83:125–151. Google Scholar


Hurlbert, S. H., J. Zedler, and D. Fairbanks. 1972. Ecosystem alteration by mosquitofish (Gambusia affinis) predation. Science 175:639–641. Google Scholar


Ingram, T., W. E. Stutz, and D. I. Bolnick. 2011. Does intraspecific size variation in a predator affect its diet diversity and top-down control of prey? PLoS ONE 6:e20782. Google Scholar


Ingram, T., R. Svanbäck, N. J. B. Kraft, P. Kratina, L. Southcott, D. Schluter, and R. Svanb. 2012. Intraguild predation drives evolutionary nich shift in threespine stickleback. Evolution 66:1819–1832. Google Scholar


IPCC. 2014. Climate Change 2014: Synthesis Report. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change [Core Writing Team, R. K. Pachauriand L. A. Meyer (eds.)]. Geneva, Switzerland. Google Scholar


Kinnison, M. T., N. G. Hairston, and A. P. Hendry. 2015. Cryptic eco-evolutionary dynamics. Annals of the New York Academy of Sciences 1360:120–144. Google Scholar


Kratina, P., H. S. Greig, P. L. Thompson, T. S. A. Carvalho-Pereira, and J. B. Shurin. 2012. Warming modifies trophic cascades and eutrophication in experimental freshwater communities. Ecology 93:1421–1430. Google Scholar


Lajoie, G., and M. Vellend. 2015. Understanding context dependence in the contribution of intraspecific variation to community trait-environment matching. Ecology 96:2912–2922. Google Scholar


Langerhans, R. B., C. A. Layman, A. M. Shokrollahi, and T. J. DeWitt. 2004. Predator-driven phenotypic diversification in Gambusia affinis. Evolution 58:2305–2318. Google Scholar


Lavergne, S., N. Mouquet, W. Thuiller, and O. Ronce. 2010. Biodiversity and climate change: integrating evolutionary and ecological responses of species and communities. Annual Review of Ecology, Evolution, and Systematics 41:321–350. Google Scholar


Lenert, L. G. 1923. Gambusia affinis: its use in mosquito control. California State Board of Health Weekly Bulletin 2. Google Scholar


Lundsgaard-Hansen, B., B. Matthews, and O. Seehausen. 2014. Ecological speciation and phenotypic plasticity affect ecosystems. Ecology 95:2723–2735. Google Scholar


Matthews, B., A. Narwani, S. Hausch, E. Nonaka, H. Peter, M. Yamamichi, K. E. Sullam, K. C. Bird, M. K. Thomas, T. C. Hanley, and C. B. Turner. 2011. Toward an integration of evolutionary biology and ecosystem science. Ecology Letters 14:690–701. Google Scholar


Meffe, G. K. 1991. Life history changes in eastern mosquitofish (Gambusia holbrooki) induced by thermal elevation. Canadian Journal of Fisheries and Aquatic Sciences 48:60–66. Google Scholar


Meffe, G. K., S. C. Weeks, M. Mulvey, and K. L. Kandl. 1995. Genetic differences in thermal tolerance of eastern mosquitofish (Gambusia holbrooki; Poeciliidae) from ambient and thermal ponds. Canadian Journal of Fisheries and Aquatic Sciences 52:2704–2711. Google Scholar


Norberg, J., M. C. Urban, M. Vellend, C. A. Klausmeier, and N. Loeuille. 2012. Eco-evolutionary responses of biodiversity to climate change. Nature Climate Change 2:747–751. Google Scholar


O'Connor, M. I., M. F. Piehler, D. M. Leech, A. Anton, and J. F. Bruno. 2009. Warming and resource availability shift food web structure and metabolism. PLoS Biology 7:e1000178. Google Scholar


Palkovacs, E. P., K. B. Dion, D. M. Post, and A. Caccone. 2008. Independent evolutionary origins of landlocked alewife populations and rapid parallel evolution of phenotypic traits. Molecular Ecology 17:582–597. Google Scholar


Palkovacs, E. P., D. C. Fryxell, N. E. Turley, and D. M. Post. 2015. Ecological effects of intraspecific consumer biodiversity for aquatic communities and ecosystems, p. 37–51. In:Aquatic Functional Biodiversity: An Ecological and Evolutionary Perspective. A. Belgrano, G. Woodward, and U. Jacob (eds.). Academic Press, Elsevier, London. Google Scholar


Palkovacs, E. P., M. T. Kinnison, C. Correa, C. M. Dalton, and A. P. Hendry. 2012. Fates beyond traits: ecological consequences of human-induced trait change. Evolutionary Applications 5:183–191. Google Scholar


Palkovacs, E. P., M. C. Marshall, B. A. Lamphere, B. R. Lynch, D. J. Weese, D. F. Fraser, D. N. Reznick, C. M. Pringle, and M. T. Kinnison. 2009. Experimental evaluation of evolution and coevolution as agents of ecosystem change in Trinidadian streams. Philosophical Transactions of the Royal Society B: Biological Sciences 364:1617–1628. Google Scholar


Palkovacs, E. P., and D. M. Post. 2009. Experimental evidence that phenotypic divergence in predators drives community divergence in prey. Ecology 90:300–305. Google Scholar


Peter, K. H., and U. Sommer. 2012. Phytoplankton cell size: intra- and interspecific effects of warming and grazing. PLoS ONE 7:e49632. Google Scholar


Post, D. M., and E. P. Palkovacs. 2009. Eco-evolutionary feedbacks in community and ecosystem ecology: interactions between the ecological theatre and the evolutionary play. Philosophical Transactions of the Royal Society B: Biological Sciences 364:1629–1640. Google Scholar


Pyke, G. H. 2005. A review of the biology of Gambusia affinis and G. holbrooki. Reviews in Fish Biology and Fisheries 15:339–365. Google Scholar


Pyke, G. H. 2008. Plague minnow or mosquito fish? A review of the biology and impacts of introduced Gambusia species. Annual Review of Ecology, Evolution, and Systematics 39:171–191. Google Scholar


R Core Team. 2015. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. Google Scholar


Rall, B. C., U. Brose, M. Hartvig, G. Kalinkat, F. Schwarzmuller, O. Vucic-Pestic, and O. L. Petchey. 2012. Universal temperature and body-mass scaling of feeding rates. Philosophical Transactions of the Royal Society B: Biological Sciences 367:2923–2934. Google Scholar


Rice, E. W., R. B. Baird, A. D. Eaton, and L. S. Clesceri (Eds.). 2012. Standard Methods for the Examination of Water and Wastewater. 22nd edition. American Public Health Association, Washington, D.C. Google Scholar


Rigosi, A., C. C. Carey, B. W. Ibelings, and J. D. Brookes. 2014. The interaction between climate warming and eutrophication to promote cyanobacteria is dependent on trophic state and varies among taxa. Limnology and Oceanography 59:99–114. Google Scholar


Rudman, S. M., M. A. Rodriguez-Cabal, A. Stier, T. Sato, J. Heavyside, R. W. El-Sabaawi, and G. M. Crutsinger. 2015. Adaptive genetic variation mediates bottom-up and top-down control in an aquatic ecosystem. Proceedings of the Royal Society B: Biological Sciences 282:20151234. Google Scholar


Rudman, S. M., and D. Schluter. 2016. Ecological impacts of reverse speciation in threespine stickleback. Current Biology 26:490–495. Google Scholar


Rundle, H. D., L. Nagel, J. W. Boughman, and D. Schluter. 2000. Natural selection and parallel speciation in sympatric sticklebacks. Science 287:306–308. Google Scholar


Sanford, E. 1999. Regulation of keystone predation by small changes in ocean temperature. Science 283:2095–2097. Google Scholar


Schoener, T. W. 2011. The newest synthesis: understanding the interplay of evolutionary and ecological dynamics. Science 331:426–429. Google Scholar


Schramski, J. R., A. I. Dell, J. M. Grady, R. M. Sibly, and J. H. Brown. 2015. Metabolic theory predicts whole-ecosystem properties. Proceedings of the National Academy of Sciences of the United States of America 112:2617–2622. Google Scholar


Schweitzer, J. A., J. K. Bailey, B. J. Rehill, G. D. Martinsen, S. C. Hart, R. L. Lindroth, P. Keim, and T. G. Whitham. 2004. Genetically based trait in a dominant tree affects ecosystem processes. Ecology Letters 7:127–134. Google Scholar


Seebacher, F., V. Ducret, A. G. Little, and B. Adriaenssens. 2015. Generalist–specialist trade-off during thermal acclimation. Royal Society Open Science 2:140251. Google Scholar


Shurin, J. B., J. L. Clasen, H. S. Greig, P. Kratina, and P. L. Thompson. 2012. Warming shifts top-down and bottom-up control of pond food web structure and function. Philosophical Transactions of the Royal Society B: Biological Sciences 367:3008–3017. Google Scholar


Stearns, S. 1983. The genetic basis of differences in life-history traits among six populations of mosquitofish (Gambusia affinis) that shared ancestors in 1905. Evolution 37:618–627. Google Scholar


Stockwell, C. A., and G. L. Vinyard. 2000. Life history variation in recently established populations of western mosquitofish (Gambusia affinis). Western North American Naturalist 60:273–280. Google Scholar


Stockwell, C. A., and S. C. Weeks. 1999. Translocations and rapid evolutionary responses in recently established populations of western mosquitofish (Gambusia affinis). Animal Conservation 2:103–110. Google Scholar


Torchiano, M. 2016. effsize: efficient effect size computation(R package). Google Scholar


Tuckett, Q. M., K. S. Simon, and M. T. Kinnison. 2017. Cultural eutrophication mediates context-dependent eco-evolutionary feedbacks of a fish invader. Copeia 105. Google Scholar


Urban, M. C., G. Bocedi, A. P. Hendry, J.-B. Mihoub, G. Peer, A. Singer, J. R. Bridle, L. G. Crozier, L. De Meester, W. Godsoe, A. Gonzalez, J. J. Hellmann, R. D. Holt, A. Huth, K. Johst, C. B. Krug, P. W. Leadley, S. C. F. Palmer, J. H. Pantel, A. Schmitz, P. A. Zollner, and J. M. J. Travis. 2016. Improving the forecast for biodiversity under climate change. Science 353:aad8466. Google Scholar


Urban, M. C., J. L. Richardson, N. A. Freidenfelds, D. L. Drake, J. F. Fischer, and P. P. Saunders. 2017. Microgeographic adaptation of Wood Frog tadpoles to an apex predator. Copeia 105. Google Scholar


Vanni, M. J. 1987. Effects of nutrients and zooplankton size on the structure of a phytoplankton community. Ecology 68:624–635. Google Scholar


Vanni, M. J., and C. D. Layne. 1997. Nutrient recycling and herbivory as mechanisms in the “top-down” effect of fish on algae in lakes. Ecology 78:21–40. Google Scholar


Vanni, M. J., and P. B. McIntyre. 2016. Predicting nutrient excretion of aquatic animals with metabolic ecology and ecological stoichiometry: a global synthesis. Ecology 97:3460–3471. Google Scholar


Vondracek, B., W. A. Wurtsbaugh, and J. J. Cech. 1988. Growth and reproduction of the mosquitofish, Gambusia affinis, in relation to temperature and ration level: consequences for life history. Environmental Biology of Fishes 21:45–57. Google Scholar


Whitham, T. G., W. P. Young, G. D. Martinsen, C. A. Gehring, J. A. Schweitzer, S. M. Shuster, G. M. Wimp, D. G. Fischer, J. K. Bailey, R. L. Lindroth, S. Woolbright, and C. R. Kuske. 2003. Community and ecosystem genetics: a consequence of the extended phenotype. Ecology 84:559–573. Google Scholar


Willing, E., P. Bentzen, C. van Oosterhout, M. Hoffmann, J. Cable, F. Breden, D. Weigel, and C. Dreyer. 2010. Genome-wide single nucleotide polymorphisms reveal population history and adaptive divergence in wild guppies. Molecular Ecology 19:968–984. Google Scholar


Yvon-Durocher, G., A. P. Allen, M. Cellamare, M. Dossena, K. J. Gaston, M. Leitao, J. M. Montoya, D. C. Reuman, G. Woodward, and M. Trimmer. 2015. Five years of experimental warming increases the biodiversity and productivity of phytoplankton. PLoS Biology 13:e1002324. Google Scholar
© 2017 by the American Society of Ichthyologists and Herpetologists
David C. Fryxell and Eric P. Palkovacs "Warming Strengthens the Ecological Role of Intraspecific Variation in a Predator," Copeia 105(3), 523-532, (19 September 2017).
Received: 13 October 2016; Accepted: 1 May 2017; Published: 19 September 2017

Back to Top