Herbivory and pollen limitation at the upper elevational range limit of two forest understory plants of eastern North America

Studies of species' range limits focus most often on abiotic factors, although the strength of biotic interactions might also vary along environmental gradients and have strong demographic effects. For example, pollinator abundance might decrease at range limits due to harsh environmental conditions, and reduced plant density can reduce attractiveness to pollinators and increase or decrease herbivory. We tested for variation in the strength of pollen limitation and herbivory by ungulates along a gradient leading to the upper elevational range limits of Trillium erectum (Melanthiaceae) and Erythronium americanum (Liliaceae) in Mont Mégantic National Park, Québec, Canada. In T. erectum, pollen limitation was higher at the range limit, but seed set decreased only slightly with elevation and only in one of two years. In contrast, herbivory of T. erectum increased from <10% at low elevations to >60% at the upper elevational range limit. In E. americanum, we found no evidence of pollen limitation despite a significant decrease in seed set with elevation, and herbivory was low across the entire gradient. Overall, our results demonstrate the potential for relatively strong negative interactions (herbivory) and weak positive interactions (pollination) at plant range edges, although this was clearly species specific. To the extent that these interactions have important demographic consequences—highly likely for herbivory on Trillium, based on previous studies—such interactions might play a role in determining plant species' range limits along putatively climatic gradients.     

Scale dependence of temperature as an abiotic driver of species' distributions

Aim Scale dependence of patterns and processes remains one of the major unresolved problems in ecology. The responses of ecosystems to environmental stressors are reported to be strongly scale dependent, but projections of the effects of climate change on species' distributions are still restricted to particular scales and knowledge about scale dependence is lacking. Here we propose that the scale dependence of those species' niche dimensions related to climate change is strongly related to the strength of climatic cross‐scale links. More specifically, we hypothesize that the strong cross‐scale links between micro‐ and macroclimatic conditions are related to high cross‐scale similarity (low scale dependence) of species' realized temperature niches and, thus, species' spatial distributions. Location This study covers seven orders of magnitude of spatial scale, ranging from local‐scale (below a metre) and regional‐scale (kilometre) investigations in central European wetland ecosystems to continental‐scale (thousands of kilometres) studies of species' distributions. Methods We combined data on the spatial occurrence of species (vegetation records at local and regional scales, digitized distribution maps at the continental scale) with information about the corresponding temperature regime of vascular plant species occurring in environmentally stable wetland ecosystems characterized by strong cross‐scale links between micro‐ and macroclimatic conditions. Results We observed high cross‐scale similarity of the characteristics of species temperature niches across seven orders of magnitude of spatial scale. However, the importance of temperature as an abiotic driver decreased nonlinearly with decreasing scale, suggesting greater importance of additional (biotic) drivers of species' occurrence at small spatial scales. Main conclusions We report high cross‐scale similarity of realized temperature niches for species inhabiting ecosystems where small‐scale environmental noise is low and cross‐scale links between micro‐ and macroclimatic conditions are strong. By highlighting a strong relationship between abiotic and biotic cross‐scale similarity, our results will help to improve niche‐based species distribution modelling, one of the major assessment tools for determining the ecological effects of climate change.     

Latitudinal variation in herbivory: influences of climatic drivers,herbivore identity and natural enemies

Although a number of investigations have concluded that lower latitudes are associated with increases in herbivore abundance and plant damage, the generality of this pattern is still under debate. Multiple factors may explain the lack of consistency in latitude–herbivory relationships. For instance, latitudinal variation in herbivore pressure may be shaped entirely or not by climatic variables, or vary among herbivore guilds with differing life‐history traits. Additionally, the strength of top–down effects from natural enemies on herbivores might also vary geographically and influence latitude–herbivory patterns. We carried out a field study where we investigated the effects of latitude and climate on herbivory by a seed‐eating caterpillar and leaf chewers, as well as parasitism associated to the former across 30 populations of the perennial herb Ruellia nudiflora (Acanthaceae). These populations were distributed along a 5° latitudinal gradient from northern Yucatan (Mexico) to southern Belize, representing one‐third of the species' latitudinal distribution and the entirety and one‐third of the precipitation and temperature gradient of this species' distribution (respectively). We found opposing latitudinal gradients of seed herbivory and leaf herbivory, and this difference appeared to be mediated by contrasting effects of climate on each guild. Specifically, univariate regressions showed that seed herbivory increased at higher latitudes and with colder temperatures, while leaf herbivory increased toward the equator and with wetter conditions. Multiple regressions including temperature, precipitation and latitude only found significant effects of temperature for seed herbivory and latitude for leaf herbivory. Accordingly, that latitudinal variation in seed herbivory appears to be driven predominantly by variation in temperature whereas latitudinal variation in leaf herbivory was apparently driven by other unexplored correlates of latitude. Parasitism did not exhibit variation with latitude or climatic factors. Overall, these findings underscore that the factors driving latitudinal clines in herbivory might vary even among herbivore species coexisting on the same host plant.     

Mutualist‐mediated effects on species' range limits across large geographic scales

Understanding the processes determining species range limits is central to predicting species distributions under climate change. Projected future ranges are extrapolated from distribution models based on climate layers, and few models incorporate the effects of biotic interactions on species' distributions. Here, we show that a positive species interaction ameliorates abiotic stress, and has a profound effect on a species' range limits. Combining field surveys of 92 populations, 10 common garden experiments throughout the range, species distribution models and greenhouse experiments, we show that mutualistic fungal endophytes ameliorate drought stress and broaden the geographic range of their native grass host Bromus laevipes by thousands of square kilometres (~ 20% larger) into drier habitats. Range differentiation between fungal‐associated and fungal‐free grasses was comparable to species‐level range divergence of congeners, indicating large impacts on range limits. Positive biotic interactions may be underappreciated in determining species' ranges and species' responses to future climates across large geographic scales.     

The role of abiotic and biotic factors in determining coexistence of multiple pollinators in the yucca–yucca moth mutualism

The determinants of a species' geographic distribution are a combination of both abiotic and biotic factors. Environmental niche modeling of climatic factors has been instrumental in documenting the role of abiotic factors in a species' niche. Integrating this approach with data from species interactions provides a means to assess the relative roles of abiotic and biotic components. Here, we examine whether the high host specificity typically exhibited in the active pollination mutualism between yuccas and yucca moths is the result of differences in climatic niche requirements that limit yucca moth distributions or the result of competition among mutualistic moths that would co‐occur on the same yucca species. We compared the species distribution models of two Tegeticula pollinator moths that use the geographically widespread plant Yucca filamentosa. Tegeticula yuccasella occurs throughout eastern North America whereas T. cassandra is restricted to the southeastern portion of the range, primarily occurring in Florida. Species distribution models demonstrate that T. cassandra is restricted climatically to the southeastern United States and T. yuccasella is predicted to be able to live across all of eastern North America. Data on moth abundances in Florida demonstrate that both moth species are present on Y. filamentosa; however, T. cassandra is numerically dominant. Taken together, the results suggest that moth geographic distributions are heavily influenced by climate, but competition among pollinating congeners will act to restrict populations of moth species that co‐occur.     

A global distribution of biodiversity inferred from climatic constraints: results from a process‐based modelling study

We investigated the connection between plant species diversity and climate by using a process‐based, generic plant model. Different ‘species' were simulated by different values for certain growth‐related model parameters. Subsequently, a wide range of values were tested in the framework of a ‘Monte Carlo' simulation for success; that is, the capability of each plant with these parameter combinations to reproduce itself during its lifetime. The range of successful parameter combinations approximated species diversity. This method was applied to a global grid, using daily atmospheric forcing from a climate model simulation. The computed distribution of plant ‘species' diversity compares very well with the observed, global‐scale distribution of species diversity, reproducing the majority of ‘hot spot' areas of biodiversity. A sensitivity analysis revealed that the predicted pattern is very robust against changes of fixed model parameters. Analysis of the climatic forcing and of two additional sensitivity simulations demonstrated that the crucial factor leading to this distribution of diversity is the early stage of a plant's life when water availability is highly coupled to the variability in precipitation because in this stage root‐zone storage of water is small. We used cluster analysis in order to extract common sets of species parameters, mean plant properties and biogeographic regions (biomes) from the model output. The successful ‘species' cannot be grouped into typical parameter combinations, which define the plant's functioning. However, the mean simulated plant properties, such as lifetime and growth, can be grouped into a few characteristic plant ‘prototypes', ranging from short‐lived, fast growing plants, similar to grasses, to long‐lived, slow growing plants, similar to trees. The classification of regions with respect to similar combinations of successful ‘species' yields a distribution of biomes similar to the observed distribution. Each biome has typical levels of climatic constraints, expressed for instance by the number of ‘rainy days' and ‘warm days'. The less the number of days favourable for growth, the greater the level of constraints and the less the ‘species' diversity. These results suggest that climate as a fundamental constraint can explain much of the global scale, observed distribution of plant species diversity.     

Loss of adaptive variation during evolutionary responses to climate change

The changes in species' geographical distribution demanded by climate change are often critically limited by the availability of key interacting species. In such cases, species' persistence will depend on the rapid evolution of biotic interactions. Understanding evolutionary limits to such adaptation is therefore crucial for predicting biological responses to environmental change. The recent poleward range expansion of the UK brown argus butterfly has been associated with a shift in female preference from its main host plant, rockrose (Cistaceae), onto Geraniaceae host plants throughout its new distribution. Using reciprocal transplants onto natural host plants across the UK range, we demonstrate reduced fitness of females from recently colonised Geraniaceae‐dominated habitat when moved to ancestral rockrose habitats. By contrast, individuals from ancestral rockrose habitats show no reduction in fitness on Geraniaceae. Climate‐driven range expansion in this species is therefore associated with the rapid evolution of biotic interactions and a significant loss of adaptive variation.     

Altitude‐mediated soil properties,not geography or climatic distance,explain the distribution of a tropical endemic herb

Understanding the ecological processes that govern species'' range margins is a fundamental question in ecology with practical implications in conservation biology. The center‐periphery hypothesis predicts that organisms have higher abundance at the center of their geographic range. However, most tests of this hypothesis often used raster data, assuming that climatic conditions are consistent across one square km. This assumption is not always justified, particularly for mountainous species for which climatic conditions can vary widely across a small spatial scale. Previous studies rarely evenly sample occurrence data across the species'' distribution. In this study, we sampled an endemic perennial herb, Thunbergia atacorensis (Acanthanceae), throughout its range in West Africa using 54 plots and collected data on (a)biotic variables, the species density, leaf mass per area, and basal diameter. We built a structural equation model to test the direct and indirect effects of distance from geographic and climatic niche centers, and altitude on Thunbergia density as mediated by abiotic and biotic factors, population demographic structure, and individual size. Contrary to the prediction of the center‐periphery hypothesis, we found no significant effect of distance from geographic or climatic niche centers on plant density. This indicates that even the climatic center does not necessarily have optimal ecological conditions. In contrast, plant density varied with altitudinal gradient, but this was mediated by the effect of soil nitrogen and potassium which had positive effect on plant size. Surprisingly, we found no direct or mediating effect of interspecific competition on plant density. Altogether, our results highlight the role of geography, climatic, and ecological mismatch in predicting species distribution. Our study highlights that where altitudinal gradient is strong local‐scale heterogeneity in abiotic factors can play important role in shaping species range limits.     

Potential distribution of the invasive tree Triadica sebifera (Euphorbiaceae) in the United States: evaluating climex predictions with field trials

Models that project introduced species distributions based on the climates in native and potential introduced ranges can provide valuable insights on the extent of a species' future spread. Yet, the lack of direct field evaluation of these range projections remains a major limitation. We evaluated results from the climex model in conjunction with results from seed and plant field trials in assessing environmental constraints to spread of the invasive tree Triadica sebifera (Chinese tallow tree) in the southeastern USA. climex incorporates key climatic parameters to generate large‐scale projections of potential distributions based on the climate across the species' current distribution. By employing field trials within microhabitats within and beyond the tree's current range, we were able to determine seed and young plants' response to the heterogeneity of the environment at regional scales. Based on projections of the climex model, T. sebifera has the potential to spread 500 km northward beyond its current distribution in the southeastern USA; minimum temperature and limited precipitation are the key climatic constraints in the eastern and western USA, respectively. climex results correlate strongly with seed germination across sites in the southeastern USA. These results do not however correlate with plant growth rates, which were often higher in sites with low projected climatic suitability. Competition and herbivory were not constraints on the growth of T. sebifera in our field trials and were therefore not responsible for the lack of correlation between model results and plant growth rates. If the minimum and maximum temperatures were to rise by 2 °C, the range of T. sebifera could extend northward 700 km beyond its current distribution. While both climex and the field trials indicate that T. sebifera is capable of extensive northward spread in the eastern USA, results of field trials indicate that the patterns of invasion within the region are likely to vary substantially with local site conditions.     

Losing your edge: climate change and the conservation value of range‐edge populations

Populations occurring at species' range edges can be locally adapted to unique environmental conditions. From a species' perspective, range‐edge environments generally have higher severity and frequency of extreme climatic events relative to the range core. Under future climates, extreme climatic events are predicted to become increasingly important in defining species' distributions. Therefore, range‐edge genotypes that are better adapted to extreme climates relative to core populations may be essential to species' persistence during periods of rapid climate change. We use relatively simple conceptual models to highlight the importance of locally adapted range‐edge populations (leading and trailing edges) for determining the ability of species to persist under future climates. Using trees as an example, we show how locally adapted populations at species' range edges may expand under future climate change and become more common relative to range‐core populations. We also highlight how large‐scale habitat destruction occurring in some geographic areas where many species range edge converge, such as biome boundaries and ecotones (e.g., the arc of deforestation along the rainforest‐cerrado ecotone in the southern Amazonia), can have major implications for global biodiversity. As climate changes, range‐edge populations will play key roles in helping species to maintain or expand their geographic distributions. The loss of these locally adapted range‐edge populations through anthropogenic disturbance is therefore hypothesized to reduce the ability of species to persist in the face of rapid future climate change.     

Macro‐ and microclimatic interactions can drive variation in species' habitat associations

Many species are more restricted in their habitat associations at the leading edges of their range margins, but some species have broadened their habitat associations in these regions during recent climate change. We examine the effects of multiple, interacting climatic variables on spatial and temporal patterns of species' habitat associations, using the speckled wood butterfly, Pararge aegeria, in Britain, as our model taxon. Our analyses reveal that this species, traditionally regarded as a woodland‐dependent insect, is less restricted to woodland in regions with warmer winters and warmer and wetter summers. In addition, over the past 40 years of climate change, the species has become less restricted to woodland in locations where temperature and summer rainfall have increased most. We show that these patterns arise mechanistically because larval growth rates are slower in open (i.e. nonwoodland) habitats associated with colder microclimates in winter and greater host plant desiccation in summer. We conclude that macro‐ and microclimatic interactions drive variation in species' habitat associations, which for our study species resulted predominantly in a widening of habitat associations under climate change. However, species vary in their climatic and nonclimatic requirements, and so complex spatial and temporal patterns of changes in habitat associations are likely to be observed in future as the climate changes.     

Concurrent biotic interactions influence plant performance at their altitudinal distribution margins

Recent studies have shown that biotic interactions can shape species' distributions, but empirical data on multiple biotic interactions are scarce. Therefore, we examined effects of plant–plant and plant–herbivore interactions on plant survival, growth and reproduction at different altitudes. For these purposes we conducted a factorial neighbor removal and large herbivore exclusion experiment with six transplant species (three tall forbs with their main distribution at low altitudes and three small forbs with their main distribution at high altitudes) on Låkta?ohkka Mountain, northern Sweden, replicated at two altitudes (ca 600 and 900 m a.s.l.) and consequently a 2.1°C difference in summer air temperatures. Overall transplant survival was 93%. Two out of three tall forbs grew better at low than at high altitudes, while no significant differences in growth between altitudes were found for any of the three small forbs. Since the main difference in abiotic conditions between the altitudes was most likely in temperature (as the sites were topographically and edaphically matched as closely as possible), this result indicates that climatic warming could induce upward migration of tall low‐altitude forbs. Negative plant–plant interactions prevailed at both altitudes, and we found indications that competition may set the lower altitudinal limits of some small tundra forbs. Thus, increased competition in response to climate warming may potentially shift the lower margins of high‐altitude forbs' distributions upward. Large mammalian grazers reduced the growth of tall forbs and enhanced the flowering of small forbs, and grazers could thus at least partly counteract the anticipated warming‐induced distribution shifts.     

The effect of competition on species' distributions depends on coexistence,rather than scale alone

One of the key problems in ecology is our need to anticipate the set of locations in which a species will be found (hereafter species' distributions). A major source of uncertainty in these models is the role of interactions among species (hereafter biotic interactions). Unfortunately, it is difficult to directly study this problem at large spatial scales and we lack a clear understanding of when biotic interactions shape species' distributions. We show a simple, direct link between the ease of species' coexistence and the importance of competition for shaping species' distributions. We show that increasing the ease of species' coexistence reduces the influence of biotic interactions. Changing the spatial scale of the analysis can reduce the influence of species interactions, but only when it promotes regional coexistence. Using these ideas, we analyze the conditions under which biotic interactions alter species' distributions in a Lotka–Volterra model of competition along an environmental gradient and argue that coexistence theory, rather than scale alone, provides a guide to the influence of species interactions. Our results provide a guide to the facets of biotic interactions that are necessary to anticipate their effects on species distributions. As such, we expect our work will help the development of more realistic distribution models.     

Competition and facilitation may lead to asymmetric range shift dynamics with climate change

Forecasts of widespread range shifts with climate change stem from assumptions that climate drives species' distributions. However, local adaptation and biotic interactions also influence range limits and thus may impact range shifts. Despite the potential importance of these factors, few studies have directly tested their effects on performance at range limits. We address how population‐level variation and biotic interactions may affect range shifts by transplanting seeds and seedlings of western North American conifers of different origin populations into different competitive neighborhoods within and beyond their elevational ranges and monitoring their performance. We find evidence that competition with neighboring trees limits performance within current ranges, but that interactions between adults and juveniles switch from competitive to facilitative at upper range limits. Local adaptation had weaker effects on performance that did not predictably vary with range position or seed origin. Our findings suggest that competitive interactions may slow species turnover within forests at lower range limits, whereas facilitative interactions may accelerate the pace of tree expansions upward near timberline.     

Lags in the response of mountain plant communities to climate change

Rapid climatic changes and increasing human influence at high elevations around the world will have profound impacts on mountain biodiversity. However, forecasts from statistical models (e.g. species distribution models) rarely consider that plant community changes could substantially lag behind climatic changes, hindering our ability to make temporally realistic projections for the coming century. Indeed, the magnitudes of lags, and the relative importance of the different factors giving rise to them, remain poorly understood. We review evidence for three types of lag: “dispersal lags” affecting plant species’ spread along elevational gradients, “establishment lags” following their arrival in recipient communities, and “extinction lags” of resident species. Variation in lags is explained by variation among species in physiological and demographic responses, by effects of altered biotic interactions, and by aspects of the physical environment. Of these, altered biotic interactions could contribute substantially to establishment and extinction lags, yet impacts of biotic interactions on range dynamics are poorly understood. We develop a mechanistic community model to illustrate how species turnover in future communities might lag behind simple expectations based on species’ range shifts with unlimited dispersal. The model shows a combined contribution of altered biotic interactions and dispersal lags to plant community turnover along an elevational gradient following climate warming. Our review and simulation support the view that accounting for disequilibrium range dynamics will be essential for realistic forecasts of patterns of biodiversity under climate change, with implications for the conservation of mountain species and the ecosystem functions they provide.     

Population dynamics can be more important than physiological limits for determining range shifts under climate change

Evidence is accumulating that species' responses to climate changes are best predicted by modelling the interaction of physiological limits, biotic processes and the effects of dispersal‐limitation. Using commercially harvested blacklip (Haliotis rubra) and greenlip abalone (Haliotis laevigata) as case studies, we determine the relative importance of accounting for interactions among physiology, metapopulation dynamics and exploitation in predictions of range (geographical occupancy) and abundance (spatially explicit density) under various climate change scenarios. Traditional correlative ecological niche models (ENM) predict that climate change will benefit the commercial exploitation of abalone by promoting increased abundances without any reduction in range size. However, models that account simultaneously for demographic processes and physiological responses to climate‐related factors result in future (and present) estimates of area of occupancy (AOO) and abundance that differ from those generated by ENMs alone. Range expansion and population growth are unlikely for blacklip abalone because of important interactions between climate‐dependent mortality and metapopulation processes; in contrast, greenlip abalone should increase in abundance despite a contraction in AOO. The strongly non‐linear relationship between abalone population size and AOO has important ramifications for the use of ENM predictions that rely on metrics describing change in habitat area as proxies for extinction risk. These results show that predicting species' responses to climate change often require physiological information to understand climatic range determinants, and a metapopulation model that can make full use of this data to more realistically account for processes such as local extirpation, demographic rescue, source‐sink dynamics and dispersal‐limitation.     

Historical and event‐based bioclimatic suitability predicts regional forest vulnerability to compound effects of severe drought and bark beetle infestation

Vulnerability to climate change, and particularly to climate extreme events, is expected to vary across species ranges. Thus, we need tools to standardize the variability in regional climatic legacy and extreme climate across populations and species. Extreme climate events (e.g., droughts) can erode populations close to the limits of species' climatic tolerance. Populations in climatic‐core locations may also become vulnerable because they have developed a greater demand for resources (i.e., water) that cannot be enough satisfied during the periods of scarcity. These mechanisms can become exacerbated in tree populations when combined with antagonistic biotic interactions, such as insect infestation. We used climatic suitability indices derived from Species Distribution Models (SDMs) to standardize the climatic conditions experienced across Pinus edulis populations in southwestern North America, during a historical period (1972–2000) and during an extreme event (2001–2007), when the compound effect of hot drought and bark beetle infestation caused widespread die‐off and mortality. Pinus edulis climatic suitability diminished dramatically during the die‐off period, with remarkable variation between years. P. edulis die‐off occurred mainly not just in sites that experienced lower climatic suitability during the drought but also where climatic suitability was higher during the historical period. The combined effect of historically high climatic suitability and a marked decrease in the climatic suitability during the drought best explained the range‐wide mortality. Lagged effects of climatic suitability loss in previous years and co‐occurrence of Juniperus monosperma also explained P. edulis die‐off in particular years. Overall, the study shows that past climatic legacy, likely determining acclimation, together with competitive interactions plays a major role in responses to extreme drought. It also provides a new approach to standardize the magnitude of climatic variability across populations using SDMs, improving our capacity to predict population's or species' vulnerability to climatic change.     

Differential effects of climate and species interactions on range limits at a hybrid zone: potential direct and indirect impacts of climate change

The relative contributions of climate versus interspecific interactions in shaping species distributions have important implications for closely related species at contact zones. When hybridization occurs within a contact zone, these factors regulate hybrid zone location and movement. While a hybrid zone''s position may depend on both climate and interactions between the hybridizing species, little is known about how these factors interact to affect hybrid zone dynamics. Here, we utilize SDM (species distribution modeling) both to characterize the factors affecting the current location of a moving North American avian hybrid zone and to predict potential direct and indirect effects of climate change on future distributions. We focus on two passerine species that hybridize where their ranges meet, the Black‐capped (Poecile atricapillus) and Carolina (P. carolinensis) chickadee. Our contemporary climate models predict the occurrence of climatically suitable habitat extending beyond the hybrid zone for P. atricapillus only, suggesting that interspecific interactions primarily regulate this range boundary in P. atricapillus, while climatic factors regulate P. carolinensis. Year 2050 climate models predict a drastic northward shift in suitable habitat for P. carolinensis. Because of the greater importance of interspecific interactions for regulating the southern range limit of P. atricapillus, these climate‐mediated shifts in the distribution of P. carolinensis may indirectly lead to a range retraction in P. atricapillus. Together, our results highlight the ways climate change can both directly and indirectly affect species distributions and hybrid zone location. In addition, our study lends support to the longstanding hypothesis that abiotic factors regulate species'' poleward range limits, while biotic factors shape equatorial range limits.     

Effects of experimental watering but not warming on herbivory vary across a gradient of precipitation

Climate change can affect biotic interactions, and the impacts of climate on biotic interactions may vary across climate gradients. Climate affects biotic interactions through multiple drivers, although few studies have investigated multiple climate drivers in experiments. We examined the effects of experimental watering, warming, and predator access on leaf water content and herbivory rates of woolly bear caterpillars (Arctia virginalis) on a native perennial plant, pacific silverweed (Argentina anserina ssp. pacifica), at two sites across a gradient of precipitation in coastal California. Based on theory, we predicted that watering should increase herbivory at the drier end of the gradient, predation should decrease herbivory, and watering and warming should have positive interacting effects on herbivory. Consistent with our predictions, we found that watering only increased herbivory under drier conditions. However, watering increased leaf water content at both wetter and drier sites. Warming increased herbivory irrespective of local climate and did not interact with watering. Predation did not affect herbivory rates. Given predictions that the study locales will become warmer and drier with climate change, our results suggest that the effects of future warming and drying on herbivory may counteract each other in drier regions of the range of Argentina anserina. Our findings suggest a useful role for range‐limit theory and the stress‐gradient hypothesis in predicting climate change effects on herbivory across stress gradients. Specifically, if climate change decreases stress, herbivory may increase, and vice versa for increasing stress. In addition, our work supports previous suggestions that multiple climate drivers are likely to have dampening effects on biotic interactions due to effects in different directions, though this is context‐dependent.