Fine-Scale Microclimatic Variation Can Shape

Integrative and Comparative Biology
Integrative and Comparative Biology, volume 56, number 1, pp. 45–61
doi:10.1093/icb/icw016
Society for Integrative and Comparative Biology
SYMPOSIUM
Fine-Scale Microclimatic Variation Can Shape the Responses of
Organisms to Global Change in Both Natural and Urban
Environments
Sylvain Pincebourde,1,* Courtney C. Murdock,† Mathew Vickers‡ and Michael W. Sears§
*Institut de Recherche sur la Biologie de l’Insecte (IRBI, CNRS UMR 7261), Université François Rabelais, Faculté des
Sciences et Techniques, Tours, 37200, France; †Department of Infectious Diseases, College of Veterinary Medicine, Odum
School of Ecology, University of Georgia, Athens, GA 30602, USA; ‡Station d’Ecologie Théorique Expérimentale, UMR
5321, CNRS et Université Paul Sabatier, 2 route du CNRS, Moulis, 09200, France; §Department of Biological Sciences,
Clemson University, Clemson, South Carolina, 29634
From the symposium ‘‘Beyond the Mean: Biological Impacts of Changing Patterns of Temperature Variation’’ presented
at the annual meeting of the Society for Integrative and Comparative Biology, January 3–7, 2016 at Portland, Oregon.
1
E-mail: [email protected]
Synopsis When predicting the response of organisms to global change, models use measures of climate at a coarse
resolution from general circulation models or from downscaled regional models. Organisms, however, do not experience
climate at such large scales. The climate heterogeneity over a landscape and how much of that landscape an organism can
sample will determine ultimately the microclimates experienced by organisms. This past few decades has seen an important
increase in the number of studies reporting microclimatic patterns at small scales. This synthesis intends to unify studies
reporting microclimatic heterogeneity (mostly temperature) at various spatial scales, to infer any emerging trends, and to
discuss the causes and consequences of such heterogeneity for organismal performance and with respect to changing land
use patterns and climate. First, we identify the environmental drivers of heterogeneity across the various spatial scales that
are pertinent to ectotherms. The thermal heterogeneity at the local and micro-scales is mostly generated by the architecture
or the geometrical features of the microhabitat. Then, the thermal heterogeneity experienced by individuals is modulated
by behavior. Second, we survey the literature to quantify thermal heterogeneity from the micro-scale up to the scale of a
landscape in natural habitats. Despite difficulties in compiling studies that differ much in their design and aims, we found
that there is as much thermal heterogeneity across micro-, local and landscape scales, and that the temperature range is
large in general (49 8C on average, and up to 26 8C). Third, we examine the extent to which urban habitats can be used to
infer the microclimatic patterns of the future. Urban areas generate globally drier and warmer microclimatic patterns and
recent evidence suggest that thermal traits of ectotherms are adapted to them. Fourth, we explore the interplay between
microclimate heterogeneity and the behavioral thermoregulatory abilities of ectotherms in setting their overall performance. We used a random walk framework to show that the thermal heterogeneity allows a more precise behavioral
thermoregulation and a narrower temperature distribution of the ectotherm compared to less heterogeneous microhabitats.
Finally, we discuss the potential impacts of global change on the fine scale mosaics of microclimates. The amplitude of
change may differ between spatial scales. In heterogeneous microhabitats, the amplitude of change at micro-scale, caused
by atmospheric warming, can be substantial while it can be limited at the local and landscape scales. We suggest that the
warming signal will influence species performance and biotic interactions by modulating the mosaic of microclimates.
Introduction
Organisms are not all facing global change equally.
Climatic change can be represented as a signal emitted
by changes in the global atmosphere that are then
received by organisms. The analogy with communication networks can help define how the ‘warming signal’
is transmitted across the environment between global
atmosphere and organisms, and also to understand
Advanced Access publication April 23, 2016
ß The Author 2016. Published by Oxford University Press on behalf of the Society for Integrative and Comparative Biology. All rights reserved.
For permissions please email: [email protected].
46
how this signal, once perceived by the receiver, can be
modified. The signal differs across latitudes and altitudes since the amplitude of warming increases towards the poles (Dillon et al. 2010; IPCC 2014).
Then, the signal can be attenuated, magnified, or
remain unchanged as it is transmitted across different
habitats, depending on their physical properties
(Woods et al. 2015). Finally, the signal is filtered by
organisms (Kingsolver and Watt 1983), biophysically
(e.g., their body coloration, shape, and size; Gates
1980), behaviorally (e.g., thermoregulation; Huey and
Slatkin 1976) and physiologically (e.g., thermal tolerance; Angilletta 2009). Understanding how the warming signal is modified by the habitat and filtered by the
organism can, therefore, be as important as defining
the amplitude of the environmental change itself (Sears
et al. 2011; Sears and Angilletta 2015). In this context,
mechanistic frameworks should be favored over correlative approaches when forecasting the ecological impacts of environmental changes (Buckley et al. 2010;
Kearney et al. 2010; Potter et al. 2013) to truly integrate the physical and biological processes behind habitat and organismal filtering of the climate signal.
When predicting the response of organisms to
global change, models use measures of climate at a
coarse resolution from general circulation models or
from downscaled regional models (Potter et al. 2013;
Kearney et al. 2014). Sophisticated biophysical
models have been developed in the 1960s (Gates
1980) to mechanistically integrate climatic variables
and ecophysiological traits of ectotherms to predict
their body temperature patterns across space and
time, and they have been applied to a great variety
of ectotherms (e.g., Tracy 1976; Gates 1980;
Kingsolver and Moffat 1982; Stevenson 1985; Grant
and Porter 1992; Helmuth 1998; Pincebourde and
Casas 2006; Kearney et al. 2009; Sunday et al.
2014). Nevertheless, most global studies using biophysical models to predict body temperature patterns
also use global or regional climatic inputs (e.g.,
Sunday et al. 2014). Individual ectotherms, however,
do not experience climate at such large scales.
Quoting George Bartholomew (1966), ‘‘. . . climate
in the usual sense of the term is little more than a
crude index to the physical conditions in which most
terrestrial animals live,’’ because these organisms are
living in ‘‘. . . cracks and crevices, holes in logs, dense
underbrush, tunnels, and nests’’ (Bartholomew 1966;
see also Huey and Bennett 2008). The climate heterogeneity over a landscape and how much of that
landscape an organism can sample will determine
ultimately the microclimates that organisms do experience (Grant and Dunham 1988; Sears et al. 2011;
S. Pincebourde et al.
Woods et al. 2015). Further, these microclimates may
magnify or buffer the amplitude of the effects of
climate change on organisms (Woods et al. 2015).
Body temperature integrates all aspects of the microclimate that ectotherms experience at their relevant
scales of space and time (Gates 1980). This past few
decades has seen an important increase in the
number of studies reporting microclimatic patterns
at small scales, focusing mostly on temperature. This
synthesis intends to unify studies reporting microclimatic temperature heterogeneity at various spatial
scales, to infer any emerging trends, and to discuss
the causes and consequences of such thermal heterogeneity for organismal performance and with respect
to changing land use patterns and climate.
In this synthesis, first we identify the environmental
drivers of thermal heterogeneity across the various
spatial scales that are pertinent to ectotherms.
Second, we survey the literature to quantify thermal
heterogeneity from the microscale up to the scale of a
landscape in natural habitats. Third, we examine the
extent to which urban habitats show microscale thermal heterogeneity comparable to natural habitats to
understand how regional warming (i.e., the urban
heat island) causes shifts in spatial heterogeneity of
microclimatic temperatures. Fourth, we explore the
interplay between microclimate temperature heterogeneity and the behavioral thermoregulatory abilities of
ectotherms in setting their overall performance.
Finally, we discuss on the potential impacts of global
change on the fine scale mosaics of microclimatic temperatures. We mainly deal with temperature as it is
one of the most important abiotic factors for ectotherms in general, and also because it is probably
the easiest variable to measure in the field. Other factors including humidity (e.g., Muth 1980) and gas
composition in the microhabitat (e.g., Pincebourde
and Casas 2016) can be especially important. Our approach can be applied to any of these factors, but the
scarcity of data on microclimatic variables other than
temperature may prevent any deep analysis such as
ours. We discuss further this point in our concluding
remarks. In addition, we focus on the spatial scaling of
thermal environments but the temporal variability of
microclimatic temperatures can also influence the response of organisms to environmental changes (e.g.,
Vasseur et al. 2014). These temporal issues are well
covered elsewhere (Dillon et al., this issue).
The drivers of microclimatic temperature heterogeneity across spatial scales
Significant deviations between macro- and microclimates are quite common in various ecosystems
Thermal heterogeneity across space
(Andrewartha 1944; Cloudsley-Thompson 1962;
Helmuth, Kingsolver and Carrington 2005; Scherrer
and Körner 2010; Suggitt et al. 2011). For example,
microhabitat temperatures for intertidal ectotherms
are poorly predicted by macro-environmental conditions (Seabra et al. 2011). Body temperature of mussels and limpets, for example, instead relies on
complex interplay among macro-environment and
local factors such as tidal regime, wave splash, radiation level, and exposure (Helmuth et al. 2006;
Mislan et al. 2009; Denny et al. 2011). Likewise, significant deviations have also been observed between
temperatures at the plant leaf surface (Linacre 1967;
Chelle 2005; Bernard et al. 2013) or ground surface
(Geiger 1959) and in surrounding air. Thus, compiling the drivers of microclimatic heterogeneity requires first that we define the spatial scale(s)
relevant for a given organism.
The environment can be seen as a complex superposition of multiple layers of air across space with
distinct properties. Tiny organisms like small insects
should be influenced by the microclimatic conditions
within the smallest boundary layer of air in their
environment (Kaspari et al. 2015). This layer can
reflect thermal variation near the rocky surface in
the intertidal (Helmuth 1998; Denny and Harley
2006), or along the plant leaf surface (Woods 2010,
2013; Pincebourde and Woods 2012). In contrast,
relatively large organisms experience thermal conditions reflective of larger boundaries of air and at
more local scales. For example, the body temperatures of large caterpillars are shaped by the air temperature within plant canopies (Woods 2013), and
forest canopies produce microclimates that diverge
from atmosphere conditions (Finnigan 2000;
Kaspari et al. 2015). The body temperatures of lizards can mirror variation in air temperature within
the ground boundary layer (Kearney et al. 2014).
Finally, the largest boundary layer reflecting the climate situation over the landscape would matter only
for the largest organisms, such as isolated trees or the
largest ectotherms such as the komodo’s dragon.
These boundary layer considerations help explain
the non-intuitive results of the biophysical modeling
of a spherical lizard, in which a large lizard gets
higher body temperatures than a small size animal
(Fig. 1). Larger sizes of spheres have larger characteristic dimensions which cause the boundary conductance for heat to increase, increasing in turn the
operative temperature. These results are consistent
with data published in Bakken and Gates (1975).
The differences in operative temperature offset the
heating rates of differently sized animals to produce
the patterns in the bottom plot of Fig. 1. Small
47
ectotherms tend to converge toward the surface temperature when this surface is cooler than the potential body temperature when in full sun. This
biophysical modeling approach also highlights the
important role of shading. The body temperature
can potentially be 30 8C higher when a large lizard
is in full sun compared to in the shade under the
vegetation for example (Fig. 1). The shade is not
distributed equally in space and this distribution
can be quite dynamic in space and time (Gates
1980). This spatio-temporal dynamics of shade also
differs between microhabitats, such as on the ground
versus on tree trunks (Huey et al. 1977). Therefore,
the boundary layer thickness should be analyzed relative to the body size of the organism.
Body size can be compared to the boundary layer
thickness to infer the actual conditions that a species
experiences. In general, the thickness of a boundary
layer can be simply calculated using atmospheric
wind speed and the size (characteristic dimension)
of the structure that generates the boundary layer,
assuming laminar air flow over a surface (Oke
1989). Therefore, the smallest boundary layer can
be computed using the smallest structure in the
microchabitat; a rock, a leaf, a branch, a piece of
decaying wood, etc. Typically, the thickness of this
boundary layer ranges from 1 mm to few centimeters
at the most (Oke 1989). This contrasts with the
thickness of more local boundary layers (e.g., plant
canopies) which can be as thick as few meters
(Finnigan 2000), and landscape boundary layers
which can be as thick as tens meters (Maitani 1979).
The drivers of microclimatic heterogeneity and effects of climate change on organism performance will
likely not be the same across different boundary
layers. For example, at the smallest boundary layer,
the drivers of temperature variation for small ectotherms include the thermal properties of the habitat
surface (bare rock, soil, leaf, etc.) as well as its geometrical features. Further, the orientation and the
slope of the surface relative to the position of the
heat source (i.e., the sun) influence greatly the
amount of energy received by the surface (Oke
1989). We propose that the thermal heterogeneity
available to tiny ectotherms can be estimated by measuring the roughness of their microhabitat, that is, the
diversity of the orientation and inclination angles
(Box 1). Even a relatively flat surface is rugose to
some extent, generating potentially diverse thermal
microenvironment for small arthropods (Caillon et
al. 2014). In intertidal organisms, the orientation
and slope of the rocky surfaces generate high levels
of thermal heterogeneity at microscale (Denny et al.
2011; Lathlean et al. 2012). If the habitat surface is
48
S. Pincebourde et al.
Fig. 1 Influence of body size on the body temperature of ectotherms. For a given set of climatic variables, organisms experience the
same conditions differently given features of their morphology, coloration, and thermoregulatory preferences and abilities. For instance,
here, we simulated radiative and thermal conditions for a hypothetical summer day. We used air temperatures, ground temperatures,
and radiation absorbed (in shade or in sun) to simulate the body temperatures of animals of differing sizes (ranging from 1-g to 1000-g,
modeled as spheres). Size alone creates thermal heterogeneity among organisms. Differences of up to 10 8C can be experienced at
midday across a range of body sizes. Thus, caution must be used when applying climatic data from weather stations without transforming the data to organism specific measures such as operative temperatures or potential body temperatures.
associated with another living organism, such as a
leaf, physiology and behavior of both species come
into play. For example, transpiration rates on plant
leaves modulate surface temperatures experienced by
insects (Pincebourde and Woods 2012). At larger
local scales (e.g., forest canopy), the diversity of leaf
angles generates most of the thermal variability found
among leaves exposed to the sun within a tree canopy
(Sinoquet et al. 2001). At even larger scales, the landscape, topography has been shown to induce a high
thermal heterogeneity that can be exploited by large
thermoregulating ectotherms (Porter et al. 2002; Sears
et al. 2011).
Microclimatic temperature heterogeneity across spatial scales
We surveyed the thermal heterogeneity reported by
studies dealing with processes at the various spatial
scales defined above. Our approach did not consist
in analyzing the whole set of published studies available. Our aim was to detail the most important and
representative studies from various biomes to suggest
a general trend and to call for more specific studies.
In this analysis, we excluded extreme thermal environments, such as geyser, for which changes of about
30 8C can be observed within a few millimeters
Thermal heterogeneity across space
49
mountain slope. The probability density function of radiative
energy incident to the surface pixels showed an important heterogeneity of energy distribution over the landscape (B). Globally,
the horizontal landscape remained cooler than the irregular
landscape. We assumed in this example that the whole surface
was exposed to solar radiation. Irregular terrain, however, would
produce shaded areas. Shaded portions would bring more surfaces at low radiative heat flux, possibly generating bimodal surface temperature distributions. For surfaces with variable
microtopography and multiple elements (e.g., bare rock, bark, and
leaves), the spatial heterogeneity in radiative heat absorption can
be estimated from the statistical distribution of slope and absorbance values. The other components of the heat budget, namely
surface evapotranspiration and heat loss via convection, can be
estimated from surface wetness and ambient air temperature. In
general, our approach allows us to estimate microclimate heterogeneity at the scale of organisms even when no spatially explicit information on environmental architecture is known.
Box 1. The major physical drivers of surface temperature are
the amount of solar radiation incident to the substrate (incoming
energy, Sp) and the portion of this energy it stores (absorptivity).
Simple geometrical relationships can be used to compute the
amount of radiation incident to a surface from factors such as
zenith angle (Z), solar azimuth (), and substrate slope (p) and
azimuth (p) (e.g., the cosine law of illumination which is based
on the angle between direct solar beam and the normal to
slope surface) (see Equations (1) and (2) (Oke 1989)). Then, the
amount of radiative energy received by a surface can be calculated for each pixel surface characterized by its angles (p, p)
(A), calculated with the sun at the zenith on June 6 in Tours,
France). The spatial distribution of the angles (p, p) should be
known to estimate the landscape thermal heterogeneity for a
small insect. In the example given here (B), the spatial distribution
of the angles (p, p) was estimated for a flat land, with a level
of rugosity sufficient to create slope with any azimuth angle, by
attributing a normal distribution to angles. We also biased this
normal distribution toward angle values of (458, 1808) for (p,
p) to simulate a similar rugous land but on a south facing
(Dunckel et al. 2009). Overall, it was not possible to
extract explicit temperature ranges across space from
most studies linking thermal tolerance of ectotherms
and environmental temperatures because these studies either focused on extreme temperatures only, or
they did not report metrics of temperature range
across space because that was not the purpose of
their studies. In addition, most studies reported temperature ranges that integrated both space and time,
typically by averaging daily minimal and maximal
temperatures within a given area. This difficulty
raises the importance of agglomerating microclimatic
raw data into a single database for all biomes and
latitudes/altitudes. Finally, we extracted the largest
temperature ranges when different conditions were
studied to remain conservative. Other metrics could
also be used but the temperature range (i.e., maximum minus minimum temperature) is the easiest
metric to compare between studies as it is inherently
standardized for ambient air temperature. Below, we
briefly detail the levels of thermal heterogeneity, and
the factors that generate them, across spatial scales
from the micro-scale to local and landscape scales.
Then, we compiled and analyzed the thermal heterogeneity data we collected.
Thermal heterogeneity within the smallest boundary
layer (microscale)
The expansion of biophysical ecology in the 60s
(Gates 1980), and the continuous developments of
thermal ecology (Angilletta 2009), generated a high
number of studies looking at the thermal processes at
microscale, including the variation in microclimate
and body temperatures and their influence on the
performance of organisms. One of the most striking
examples is the plant leaf microclimate (Jones 1999;
50
Stokes et al. 2006; Pincebourde and Woods 2012;
Pincebourde and Casas 2015). The surface of leaves
is quite heterogeneous with temperature ranges from
4 to 7 8C over single leaf surfaces (Woods 2013;
Caillon et al. 2014; Pincebourde and Suppo, this
issue). The temperature gradient over a leaf surface
is generated by wind effect across the surface (RothNebelsick, 2001), by possible variation in the stomatal behavior (Mott and Buckley 2000), and by the
leaf microtopography (Saudreau and Pincebourde,
personal communication). Therefore, small arthropods can encounter a great variety of surface temperatures by moving over few centimeters (Woods
2013; Caillon et al. 2014).
The rocky intertidal ecosystem provides another
model with high variability at small scale, within a
meter. Rock surface temperature (Lathlean et al.
2012) and mussel body temperature (Denny et al.
2011) vary by as much as 8 8C and 15 8C, respectively,
at a given time of the day and over short distances due to
the microtopography of the rock surface and to variation in orientation of mussel’s body toward the sun
(Dowd, personal communication). Lizards have also
been shown to exploit microscale thermal variation
(Muth 1977). Bakken (1989) showed that simple postural adjustments of a small lizard on an arboreal perch
could result in changes of up to 6.5 8C, and simply
moving from sun to shade could affect the operative
temperature by up to 16 8C. But the record for terrestrial systems probably is the microhabitat selected by a
flat-rock spider in Australia. The underside of rocks
selected for nesting can have a temperature range of
up to 20 8C over a few centimeters (Pike et al. 2012),
which results from the high heat load absorbed by the
rock in full sun and the much lower temperature of the
substrate it is resting on. A temperature range of similar
amplitude (about 17 8C) was also found around rocks
used by snakes as retreat-sites in northeastern
California (Huey et al. 1989). Other studies did
report similar temperature ranges for various types of
systems (Goller et al. 2014; Mota et al. 2015).These
studies demonstrate that ectotherms can access a thermal gradient over very small spatial scales (51 m) that is
as variable as thermal gradients found over larger, geographical scales (Kaspari et al. 2015; e.g.,41600 km in
the intertidal, Denny et al. 2011; 750 km of forest,
Caillon and Pincebourde, personal communication).
Thermal heterogeneity within intermediate boundary layers (local scale)
The local scale, as defined here, can represent spatial
scales similar to a portion of a crop field, a section of
forest with several large trees or the scale of a single
S. Pincebourde et al.
tree canopy. The size of this spatial scale largely surpasses body size, but organisms can move within the
local scale via short-range dispersion. The local scale
can also correspond to the home range of the individual (Briscoe et al. 2014). The recent technological
advances in unmanned aerial vehicles (UAVs)
equipped with infrared cameras has considerably
changed the way we probe the spatial heterogeneity
at local scale. For example, surface temperature
ranges of 20-30 8C were reported over 20 m in
crop fields in Andean agrosystems at the same time
of the day (Faye et al. 2015). Other approaches include thermography in mountains where one can
look at a slope frontally (Scherrer and Körner
2010) or infrared cameras fixed on crane towers
(Leuzinger and Körner 2007). For example, the surface temperature of temperate tree canopies is extremely variable with leaf surface temperatures
varying in a range of 15 8C (Leuzinger and
Körner 2007). This heterogeneity was species specific,
with each single canopy species showing a specific
temperature pattern (Leuzinger and Körner 2007).
In general, authors identified the plant structure as
the cause of this variability. Indeed, the underlying
process is again related to the architecture of the
plant canopies: the temperature distribution should
be related to the distributions in orientation and inclination angles of the leaves or other plant organs
(Sinoquet et al. 2009). The same process applied to
other types of surfaces (Cox and Smith 2011).
Thermal heterogeneity within the largest boundary
layers (landscape scale)
Thermal heterogeneity at the landscape scale integrates
local heterogeneity at the scale of the individual’s home
range to the area encompassing a population of interacting organisms. Largely, this level of heterogeneity is
driven by the attenuation of solar radiation over the
landscape surface as a result of differing slopes and aspects, vegetational cover (shaded vs. open), and underlying surface composition (soil type, rocky exposures,
litter composition, etc.; Sears et al. 2011; Suggitt et al.
2011; Ashcroft et al. 2012). Landscape heterogeneity is
likely to be very important given that many applications, such as species distribution models, assume climate envelopes over large swaths of habitat (30 m to
greater than 1 km pixels). For instance, using a simulation based approach, Sears et al. (2011) showed that
thermal heterogeneity of the substrate could be increased from negligible to greater than 20 8C across differently configured topographies (flat to hilly), despite
using the same climatic drivers. Likewise, using infrared
thermography, Scherrer and Körner (2010) showed an
Thermal heterogeneity across space
51
increase in the variation of surface temperatures from
minimal (up to 5 8C) at night to quite high (420 8C)
during the day in alpine grasslands. These findings carry
implications for predicted consequences of climate
change (but see below).
Comparative analysis
We collected temperature range data from various
ecosystems, including plant surfaces and intertidal
substrates (See Supplementary Table S1). We categorized the temperature range values according to
the spatial scale they represent: micro-scale, localscale and landscape scale. Because one has a priori
expectation of ordering between these three categories, we ran an ordered-heterogeneity test using the
complement of the P-value of a Kruskal–Wallis
analysis of variance (Pc) and the Spearman’s rank
coefficient (rc) (OH tests; Rice and Gaines 1994).
Globally, the mean temperature range at the
micro, local, and landscape scales did not differ
(Fig. 2; OH test: Pc x rc ¼ 0.36550.700 at a probability threshold ¼ 0.05, NS), although mean temperature range tended to increase with area.
Therefore, there is as much thermal heterogeneity
at all spatial scales, at least in terms of composition
(i.e., the ‘quality’ of temperature patches). The temperature range is not informative of the spatial
configuration of the various temperature patches.
Other metrics should be used but the lack of coherence between published studies and the lack of standard methodology impede further analyses on
configuration. Indeed, we lack studies that compared the empirical levels of thermal heterogeneity
across each relevant spatial scale on the same model
system and in a standard way. To our knowledge,
Pincebourde and Suppo (this issue) provided the
only attempt to do so across micro, local, and landscape scales. This study reported a significant increase in thermal variance and temperature range
when scaling up from leaf surface to landscape,
mostly because landscapes are more diverse structurally than single leaf surfaces (Pincebourde and
Suppo, this issue). By contrast, a modeling approach showed that the variance in mussel body
temperatures was larger when looking at local
scale compared to broader scales (Gilman et al.
2006). The tendency for increased heterogeneity
when scaling up is certainly more or less pronounced depending on the type of ecosystem.
Although data from the intertidal ecosystem were
included into the dataset, our analysis focused on
open-air temperature range values. The temperature
heterogeneity at a small spatial scale has been largely
Fig. 2 Mean temperature range (SD) at three different spatial
scales. Sample sizes (number of studies) are N ¼ 3, 11, and 12 at
the scales of landscape, local scale, and micro-scale, respectively
(see Table S1).
disregarded in most aquatic systems (streams, ponds,
etc.), except in the context of the vertical stratification
of deep lakes. The water temperature in lakes varies
substantially with depth, with water temperature differences of more than 10 8C at a depth410 m
(Mortimer 1952; Hondzo and Stefan 1993; Boehrer
and Schultze 2008) as a result of the heat flux balance
at the water surface (Dake and Harleman 1969). The
water temperature range within single ponds was
quantified in rare occasions. Moss (1969) reported
variation in water temperature of up to 5 8C with
depth in a 3 m-deep pond, but this vertical gradient
disappeared during the hottest period of the year,
offering few thermal refuges to ectotherms. Thermal
heterogeneity was also found in rivers and streams,
despite the water flow tends to homogenize water
temperature. The drivers of thermal heterogeneity
are mostly linked to the geo-morphology of the
river basin and to the structural complexity of the
channels (Hawkins et al. 1997; Caissie 2006). Spatial
variations of only a few degrees were found in general
between various biotopes (e.g., rifle vs. run vs. pool)
but even such a small micro-scale heterogeneity was
enough to generate significant deviations in degreedays for instance (Dallas and Rivers-Moore 2011). A
water temperature range of up to nearly 20 8C, however, was recorded in freshwater systems composed of
several types of channels (Wawrzyniak et al. 2013;
Brewitt and Danner 2014). Therefore, aquatic habitats
may show levels of thermal heterogeneity similar to
terrestrial habitats. To our knowledge, only one study
compared thermal metrics between aquatic and
nearby terrestrial habitats, reporting large temperature
52
changes across the aquatic-terrestrial interface
(Tonolla et al. 2010). But this study did not include
metrics such as temperature range. More field-based
research is needed to fully quantify the water temperature heterogeneity at small-scale.
Our analysis reveals that generally there is a high
level of surface temperature heterogeneity within the
micro- and local scales. This result implies that ectotherms can potentially find enough microclimate
diversity locally to thermoregulate properly and
escape lethal temperatures as long as the microclimatic temperature distribution remains mostly
below the temperature threshold of the species.
Indeed, this effect was also found in sessile organisms, with some individuals surviving to extreme
events when they live in the less exposed microhabitats (Dowd et al. 2015). This conclusion contradicts
a recent global analysis which indicated that ectotherms are barely able to buffer the amplitude of
warming by moving over a radius of 50 or 100 km
(Buckley et al. 2013). This study, however, integrated
climatology, topography and behavior at a coarse
grid size (10’ latitude), which may have confounded
the role of local thermal heterogeneity in substrate
and body temperatures (but see below).
Ultimately, the eco-physiology of ectotherms is
driven by their body temperature. Body temperature
can deviate from surface or habitat temperature
somewhat importantly due to the biophysical processes governing the energy budget of relatively
large organisms (Gates 1980). The realized body temperatures achieved by organisms in the field can be
compared to the surface or habitat temperature distributions to estimate the thermoregulatory abilities
or to what extent their biophysical properties magnify or reduce the environmental variance. In general, ectotherms remain at body temperatures close
to their surrounding habitat temperature when in the
shade but body temperature increases when exposed
to solar radiation, causing inevitably an increase in
the body temperature range at local scale compared
to habitat temperature range (Pincebourde et al.
2007; Sunday et al. 2014). Several important studies
showed recently that the body temperature patterns
of ectotherms put them at risk under a warming
scenario, in particular by restricting their activity
levels (Huey et al. 2012; Kearney 2013).
Microclimatic heterogeneity in urban
environments
Human-modified habitats potentially show contrasting microclimatic patterns compared to natural habitats. Variation in socio-economic status and human
S. Pincebourde et al.
behavior shape the environments that organisms experience in urban ecosystems, potentially resulting in
changes in population growth rates, species interactions, disease transmission, and ultimately organism
fitness (Meineke et al. 2013). In general, urban centers tend to have very different temperature (Oke
1982; Arnfield 2003) and precipitation regimes
(Lacke et al. 2009; Shepherd et al. 2010; Niyogi et
al. 2011) relative to surrounding areas. Urban and
suburban landscapes have four main characteristics
that differ from rural areas (Larsen 2015): (1) the use
of dark, dense paving, and building materials, (2)
reduced amounts of vegetative cover, which decreases
natural shading and cooling by evapotranspiration,
(3) three-dimensional buildings that restrict air circulation and absorb solar radiation, and (4) the addition of ‘‘waste heat’’ due to increases in vehicular
traffic, air conditioning units, additional lighting,
and
industrial
equipments
or
processes.
Consequently, in the continental United States, for
example, cities have mean temperatures that are on
average 3 8C higher than non-urban areas, with the
exception of drier ecosystems. Night time temperatures also tend to be higher in urban centers resulting in smaller ranges in diurnal temperature
fluctuation (Kalnay and Cai 2003). These shifts in
temperature associated with urban heat islands can
also translate into changes in other weather variables,
such as precipitation regimes, wind speed, relative
humidity, etc. (USGCRP 2014), all of which will
impact the heat budget of organisms living in urbanized landscapes.
In addition to urban centers being on average
warmer, there is substantial variation in thermal
metrics within and across urban landscapes. The drivers of such variability are diverse and they operate
at different spatial scales. First, the material used for
urban infrastructures greatly impacts surface temperatures. Surface temperature deviations of up to 20 8C
can be observed between different materials under
the same environmental conditions (Doulos et al.
2004). Not only the nature of the material matters,
but also does the texture (roughness) and the surface
color. By contrast, the geometrical feature of the surface (size and thickness) does not generate significant
surface temperature deviations (Doulos et al. 2004).
Second, the urban architecture impacts the net solar
radiation and sensible heat flux, thereby generating
different surface temperatures depending on the orientation of the urban ‘‘canyon’’ (streets) (Todhunter
1990). Therefore, important microclimatic variation
may result from urban systems differing in their
symmetry and orientation. Third, the biophysical
composition of urban landscapes directly causes
Thermal heterogeneity across space
heterogeneity in surface temperatures. In particular,
the relative surface area corresponding to plant cover
and infrastructures impacts surface temperature in a
nonlinear fashion (Stabler et al. 2005; Guo et al.
2015). Interestingly, surface temperatures decrease
with distance from the core of the city on average,
but vegetation patches can modify this relationship
(Stabler et al. 2005). Finally, the amount of ‘‘waste
heat’’ generated by human activities introduces another source of thermal variation in urban areas
(Harlan et al. 2006; Jenerette et al. 2007).
Therefore, these drivers interact with each other to
generate the specific thermal signatures of urban
systems.
The composition of the urban landscape is important, but the spatial configuration of these drivers
matters as well (Li et al. 2013). For example, in
Baltimore (MD, United States), the relative surface
area of buildings increases surface temperatures while
plant cover mitigates the urban heat island (Zhou et
al. 2011). However, keeping the composition of the
urban landscape constant, changing the spatial configuration of land cover features (e.g., vegetation
edge density and complexity of vegetation) led to
different surface temperature patterns (Zhou et al.
2011). The physical connectedness between different
vegetation units (e.g., trees) can also contribute to
decrease the average surface temperature (Rhee et
al. 2014). Overall, edge effects seem to be important
in urban systems and likely contribute to the thermal
heterogeneity in cities.
Urban areas reach a level of complexity that differs
from natural habitats in that the geometrical features
of buildings are much more predictable and constant
than the architecture of vegetation which changes
dynamically throughout the season and from year
to year. The spatial configuration effects as caused
by the relative arrangement of the various infrastructures and vegetation areas would also be more predictable than the configuration effects in natural
systems, simply because the main drivers (architecture and material) are easily identifiable and well
known. Then, the thermal heterogeneity across
space in urban areas would be more predictable
than in natural habitats. The predictability of urban
landscapes makes them ideal systems for modeling
processes in a spatially explicit context (Zhang et
al. 2013).
Are ectotherms adapted to their urban mosaic of
microclimates? Answers to this question may help to
resolve how organisms will adapt to climate change.
Theory predicts that urban populations would exhibit higher and lower tolerance to extreme warm
and cold temperatures, respectively, relative to rural
53
populations because urban environments are characterized by warmer daily minimum, mean, and maximum temperatures (Oke 1982). Yet, there have been
relatively few studies to compare how urban populations differ in thermal tolerance and performance
relative to populations living in natural landscapes.
Angilletta et al. (2007) did find systematic differences
in the ability of urban leaf cutter ants to tolerate a
high temperature (42 8C) relative to rural leaf cutter
ants, suggesting that thermal tolerance of some
urban organisms can respond rapidly to changes in
climate. Whether this difference in tolerance across
urban and rural ant populations was due to phenotypic plasticity or genetic change was not determined. Further, there is some evidence to suggest
that the capacity for evolutionary shifts in upper
thermal limits is more constrained than lower thermal limits (Snyder and Weathers 1975; Araujo et al.
2013), and that a 2-4 8C warming scenario could be
problematic for many ectotherms, especially those
living in mid-latitude continental environments
(Hoffmann et al. 2013). Interestingly, the climatic
gradient from rural areas across urban landscapes is
expected to be larger under temperate compared to
tropical regions (Diamond et al. 2015), suggesting
more difficulties for ectotherms to adapt to urban
thermal environments when coming from surrounding rural habitats at mid-latitudes.
The idea that urban areas can reflect the climates
we will experience in the future is not new. Recently,
a study showed that cities simulate the global warming effect on scale insects (Youngsteadt et al. 2015).
Urban environments may, therefore, be a window to
future microclimatic mosaics for ectotherms. Urban
areas can be used in experiments to simulate the
effects of climate change on organisms (thermal tolerance, water balance, performance, nutrient regulation, etc.). Understanding how microclimate
variation in urban environments affects organism
performance could also provide information on
how we can mitigate the negative effects of climate
warming and extreme events on organisms.
The interplay between spatial heterogeneity and thermoregulatory abilities
The thermal environment is variable in space, and
how variable it is—or the degree of thermal heterogeneity in the environment—is most usefully defined
with reference to a target organism. Size and motility
of the organism are the factors determining the effective spatial heterogeneity of the thermal environment (Porter and Gates 1969; Vickers et al. 2011;
Sears and Angilletta 2015). As organisms become
54
larger (or faster), the grain over which they experience the environment becomes coarser due to their
increased thermal inertia (or to the relatively small
amount of time experiencing any one temperature
over its movement path). Smaller (or slower) organisms, on the other hand, experience the same environment at a much finer scale due to their quickly
changing body temperatures and relatively longer
durations at a particular microclimate. As such,
one could imagine interactions where small/slow organisms experience thermal environments at the
same grain as larger/faster organisms.
Despite the complexity of defining thermal heterogeneity, some properties are general. Spatial heterogeneity in the thermal environment can enable or
hinder behavioral thermoregulation depending on
whether the thermal environment has a temperature
distribution similar to- or different from- the preferred body temperature of the thermoregulator
(Fig. 3). To explore this in detail, consider a
random walk through the habitat that yields the perfect temperature in terms of both accuracy (location,
or mean body temperature), and precision (breadth,
spread, or variance of temperature). We assume our
real organism targets a single temperature value (i.e.,
its preferred body temperature) for optimal performance, and that the more random its walk is, the less
effort it requires. This might be because the organism
can devote energy toward other objectives (e.g., reproduction) while maintaining body temperature for
free. The final assumption is that on account of thermal inertia and the time to sample the habitat, our
random walker takes a certain amount of time to
realize its ultimate or equilibrium temperature
distribution.
If the habitat is as close to optimal as possible, the
random walker can achieve perfect temperature control and, therefore, a real organism in such an environment must invest no further effort in behavioral
thermoregulation (Huey and Slatkin 1976). Consider
that the total range of temperatures and the mean
temperature does not change in the habitat but that
the temperature patches are arranged differently—for
example, they are not aggregated in the heterogeneous habitat but are aggregated in the less heterogeneous. The random walker spends less time at any
one habitat temperature and therefore integrates over
time toward the mean more rapidly in a heterogeneous habitat. In this ergodic system, therefore, the
random walker in habitats of increasing spatial heterogeneity realizes its ultimate temperature distribution sooner, and has a narrower temperature
distribution than in less complex habitats. Thus,
the temperature distribution will have a narrower
S. Pincebourde et al.
temperature range (Fig. 3 top). Our random walker
achieves nearly optimal body temperature more
easily or more efficiently in this case. By corollary,
decreasing thermal heterogeneity increases the time
to realize the distribution and increases the final
temperature distribution breadth (a decrease in precision, Fig. 3 bottom). For our real organism to behaviorally thermoregulate, it must diverge further
from the random walk as thermal heterogeneity decreases: thermoregulation is more difficult. At the
limits of thermal heterogeneity (either perfectly
mixed, or zero heterogeneity), the thermal distribution becomes effectively constant and equal to the
mean, and behavioral thermoregulation becomes logically impossible. Sears and Angilletta (2015) showed
that the relative size, arrangement, and heterogeneity
of thermal patches affect the thermoregulatory performance of individuals similarly as described above.
As the relative abundance, dispersion of patches, and
within patch heterogeneity increased, so did thermoregulatory performance. Further, larger organisms
are less affected by all of these factors. In fact, in
these simulations, heterogeneity and arrangement
seemed to affect thermoregulatory ability as much
as the average temperature of the environment
itself. Taken together, these results demonstrated
that studies of behavioral thermoregulation should
continue to consider the spatial arrangement of temperature patches per se, instead of focusing exclusively on the mean temperature or temperature
range.
Now consider shifting mean temperature with
varying spatial heterogeneity; a random walker in
habitat with suboptimal mean temperature and
high heterogeneity will have a narrow body temperature range that is far from the optimum. A thermoregulator here must invest high effort to correct its
body temperature (Fig. 3, top left). On the other
hand, suboptimal temperature and low heterogeneity
may result in off-center but variable body temperature that still overlaps the thermal optimum (Fig. 3,
bottom left). A thermoregulator might need to invest
relatively less effort here to maintain good control
over body temperature. This interaction of mean
temperature with its spatial heterogeneity can produce complex results.
This thought experiment involves habitats of different thermal heterogeneity, and holds provided the
habitats have similar mean and variance, that is, we
discuss only spatial heterogeneity. The situation will
become more complex as more factors are introduced: increasing the maxima, variance, or skewness
of the environmental temperature distribution will
have further effects on the performance of
Thermal heterogeneity across space
55
heterogeneity should be a priority as shifting thermal
heterogeneity may compound the effect of shifting
location and composition of the distributions.
Which way environments move on the two-way heterogeneity/location axis (Fig. 3) will determine the
severity of climate change to behavioral thermoregulators who, if not plastic enough in their response,
may find themselves suddenly in the wrong patch.
The effects of warming on the microclimate mosaic and ectotherms’
performance
Fig. 3 The interplay between thermal heterogeneity and mean
temperature on temperature experienced by a random walker.
The vertical axis represents some value of spatial heterogeneity
(e.g., agglomerating index, patchiness) while mean temperature
and temperature range remains the same. Therefore, habitat
heterogeneity refers to the spatial configuration of temperature
patches, and the spatial composition remains unchanged. The
horizontal axis represents a variation in mean temperature while
the spatial configuration and the temperature range remain
unchanged. In each quadrant is a realized body temperature
distribution, with optimal temperature indicated (Topt). In a
highly heterogeneous habitat (top), the body temperature experienced will be narrower than an equivalent habitat of lower
heterogeneity (bottom). As the mean temperature of the habitat
moves away from the optimum (left, only colder is shown, but
the same principal applies for hotter), the body temperature
distribution shifts toward ‘‘wrong’’ body temperatures. Then, the
walk in the more heterogeneous habitat suffers more than the
walk in the low heterogeneous environment as a byproduct of
the narrowness of the distribution.
thermoregulation (e.g., Sears and Angilletta 2015). In
particular, increased positive skewness will compound the effect of increased mean temperature
(Vasseur et al. 2014) and increased variance may
alter the variance of realized body temperature
(Huey and Slatkin 1976). Understanding and predicting precisely how global change impacts local
Theory in thermal ecology predicts that the physiological limits of ectotherms track shifts in environmental temperature. The mean tolerance of
ectotherms to elevated temperatures, however,
varies only weakly along latitude and altitude
(Addo-Bediako et al. 2000; Sunday et al. 2011;
Garcı́a-Robledo et al. 2016). These macroecological
studies also reveal a high level of variance between
species or population at a given latitude/altitude.
Physiological limits are certainly more correlated to
very local conditions rather than to global climatic
estimates (e.g., Stillman and Somero 1996; Kaspari et
al. 2015). Indeed, the temperature range we report at
the landscape scale (15 8C) is comparable in amplitude to the global mean temperature difference between the tropics and temperate regions, a metrics
that is often used in global analysis. In addition, the
inter-specific variability of CTmax at temperate latitudes (around 7 8C; Addo-Bediako et al. 2000;
Sunday et al. 2011) compares well to the levels of
thermal heterogeneity at the landscape, local and microscales. For example, the lethal temperature threshold of a temperate leaf mining caterpillar varies in
response to the changes in the microclimatic temperature across larval stages (Pincebourde and Casas
2015). This study shows that, at a given time in
the summer, the CTmax of larvae varies by about
6 8C when the microclimatic temperature range
encompasses about 7 8C within a single apple tree
canopy (Pincebourde et al. 2007). Therefore, we propose that the local microclimatic heterogeneity explains the lack of global coherence in the patterns
of thermal tolerance across space. This suggests that
the warming signal will influence species performance and biotic interactions by modulating the
mosaic of microclimates.
Nevertheless, our understanding of how global
change will modify the spatial distribution and arrangement of the various microclimates in a given
habitat is weak. Most of our knowledge is inferred
from laboratory experiments or simulation studies.
56
For example, warming ambient air reduced thermal
heterogeneity along leaf surfaces (Caillon et al. 2014).
The temperature range was substantially reduced at
the apple leaf surface, and the aggregation index was
increased, when air temperature increased (Caillon et
al. 2014). The thermal homogenization at microscale has consequences for the spider mites that
thermoregulate at the leaf surface: under elevated
temperatures, spider mites tried to remain in the
coldest parts of the leaf surface (Caillon et al.
2014). These organisms probably rely on physiological mechanisms to increase their resistance to high
temperatures. Further, a modeling study showed that
tree architecture influences significantly a leaf miner’s
larval development time and mortality, but that this
effect of manipulating tree structure was weak compared to the impact of climate change by 2100
(Saudreau et al. 2013). This suggests that at some
point the amplitude of warming can surpass the benefits of the thermal heterogeneity at micro and local
scales, especially for ectotherms that already select
microclimates that are warmer than ambient air
(e.g., leaf miners; Pincebourde et al. 2007) and for
species that are already challenged by their current
heterogeneous environment (e.g., lizards that delay
their activity to mornings; Grant and Dunham
1988). This ‘point of no return’ represents a period
at which no thermal refuge exists anymore in space,
causing the extinction of the population. This effect
might occur during an extreme event such as a heat
wave. The generality of this concept to any organism
remains to be assessed.
At the landscape scale, the topography can mitigate the negative impacts of warming on organisms.
For example, a shift in operative temperatures and
activity times for lizards were less extreme as topographic relief increased from flat to hilly (Sears et al.
2011). Similarly, a 2 8C warming in the alpine landscape caused the loss of only 3% of the coolest microhabitats, while 75% of current microhabitats were
reduced in abundance and 22% of microhabitats (the
warmest portion) increased in abundance (Scherrer
and Körner 2011). In this study, Scherrer and Körner
(2011) stated that ‘‘meter-scale thermal contrasts significantly exceed IPCC warming projections for the
next 100 years’’ in alpine ecosystems. While this
implies little risk of temperature-driven extinctions
of species, the competition for cool microclimates
may be exacerbated in a warmer world (Scherrer
and Körner 2011). These suggested that increased
thermal heterogeneity produced through processes
at the landscape scale might mediate the responses
to climate, and that predictions of the negative responses of species to increasing temperatures (e.g.,
S. Pincebourde et al.
multiple extinctions; Sinervo et al. 2010) might be
largely overestimated (Clusella-Trullas and Chown
2011). Yet, global studies largely ignore the influence
of local heterogeneity (but see exceptions; e.g.,
Helmuth et al. 2006).
Concluding remarks
Environmental heterogeneity promotes biodiversity
(Stein et al. 2014). There are multiple ways of quantifying heterogeneity. Our synthesis reveals some
methodological issues in how ecologists quantify
the spatial heterogeneity in surface or body temperatures. Our meta-analysis of spatial temperature
range was considerably limited by the fact that authors routinely report a temperature range that confound both space and time, simply because these
studies did not focus particularly on the spatial
issue alone. We encourage authors to provide the
temperature distributions in space as supplemental
materials, for example, or to upload their raw data
in data repositories. This approach requires that the
spatial resolution of the data is explicitly reported
along with the data. Indeed, the heterogeneity of microclimate and body temperatures at a fine spatial
scale is largely unknown for the vast majority of species. Such knowledge is crucial since the spatial resolution of the data can impact the inference we make
on the effects of climate changes (Potter et al. 2013).
The rise in the use of thermal imaging to explore
thermal landscapes (Faye 2015; Faye et al. 2015,
2016) will also certainly improve our global assessment of spatial heterogeneity in both natural habitats
and urban environments.
The use of environmental architecture or geometry
to assess the potential diversity of surface temperatures at micro and local scales is surprisingly rare.
Advanced technical tools exist, however, to measure
architectural traits in the micro-environment such as
LIDAR (light detection and ranging; Vierling et al.
2008), scanning lasers and electromagnetic digitizers.
Scanning lasers capture the 3D structure of surfaces
at a finer scale than LIDAR, and sometimes authors
coined the term of microtopography to describe the
level of roughness measured with these equipments.
For example, the use of a scanning laser to describe
the microtopography in the rocky intertidal showed
that the local distribution of intertidal ectotherms is
best predicted by the roughness of the rock surfaces
(Peglow 2013). Scanning lasers are robust and can
even be used underwater like in deep-sea ecosystems
(Du Preez and Tunnicliffe 2012). Electromagnetic
digitizing device allows to precisely measure the 3D
architecture of volumic structures such as vegetative
57
Thermal heterogeneity across space
elements within a tree or forest canopy (Sinoquet et
al. 2001). For example, this approach was applied to
predict the distribution of microclimates for insect
pests within apple tree canopies (Pincebourde et al.
2007; Saudreau et al. 2011, 2013). Therefore, the
technology already exists to capture the architectural
properties of any ecosystem, and this is the first step
toward the exhaustive quantification of the thermal
heterogeneity in space.
Our approach developed for temperature can be
applied to any environmental variable. Indeed, the
microclimate of organisms cannot be summarized
by temperature alone. The environment of organisms
is multifactorial and we need to integrate all these
abiotic variables to quantify the net effect of global
change on species. Different factors such as temperature, humidity, rainfall, physical disturbance (e.g.,
wave splash) can vary in their level of variance in
both time and space (Benedetti-Cecchi et al. 2006;
Koussoroplis and Wacker 2016). The spatial configuration of each variable should be considered, although some of them covary strongly (e.g.,
humidity and temperature). An important issue is
that a biological response to temperature can be misinterpreted as a direct effect of temperature while it
results from the effect of a variable that covaries with
temperature. For example, dissolved oxygen varies
with temperature in freshwater systems to the point
that thermal tolerance can be limited by dissolved
oxygen deprivation, a phenomenon that could be
expanded to oxygen limitation in terrestrial invertebrates (Verberk and Calosi 2012; Boardman and
Terblanche 2015). Other factors of anthropogenic
origin such as pollutants can also interact and
lower the performance of ectotherms. For example,
elevated temperatures and cadmium exposure led to
oxygen limitation by impaired performance in
oxygen supply through ventilation and circulation
in oysters, suggesting that the thermal tolerance is
lower in marine ectotherms exposed to polluted
waters compared with pristine environments
(Lannig et al. 2008). In general however, the influence of the covariance structure of the multifactorial
environment on the performance of ectotherms remains to be determined (Koussoroplis and Wacker
2016). Biophysical models of ectotherms’ body temperature have the advantage of integrating environmental variables into a single output, the body
temperature. Therefore, these models seem to be
good candidates to tackle these questions (but see
Denny et al. 2009; Dowd et al. 2015).
The spatial heterogeneity of microclimatic temperatures as described in this review, and combined with
the temporal variability of temperatures (Dillon et al.,
this issue; Sheldon and Dillon, this issue), is ultimately filtered by organisms both behaviorally and
physiologically. The way organisms exploit this
spatio-temporal heterogeneity is influenced by biotic
interactions, including herbivory (Pincebourde and
Casas 2015), competition (e.g., Crespo-Pérez et al.
2015) and predation (e.g., Barton and Schmitz
2009). Climate change can impact species by strengthening these biotic interactions, putting prey, and the
weakest competitors under risk of extinction by modifying the activity of interacting species (Barton and
Schmitz 2009; Huey et al. 2012). We suggest that the
spatial scale at which we shall expect the strongest
shifts in thermal niches and species distributions can
be the scale at which climate change causes the largest
modification in the thermal heterogeneity. Analyzing
the levels of microclimate heterogeneity across spatial
scales promises net progresses in global change
biology.
Acknowledgments
We thank Arthur Woods and Michael Dillon for
organizing the symposium ‘‘Beyond the mean:
Biological impacts of changing patterns of temperature variation’’ at the annual meeting of the SICB in
Portland, OR (January, 2016). We also acknowledge
the other participants to the symposium for fruitful
discussions. We warmly thank Raymond Huey, Wes
Dowd and Olivier Dangles for their feedback on a
previous version of the manuscript. The symposium
was sponsored by the divisions DAB, DCE, DCPB,
DEE, and DIZ.
Funding
This work was partially supported by the SICB,
Agence Nationale pour la Recherche [MicroCliMite,
ANR-2010 BLAN-1706-02 to SP; PROTECTODO,
Région Centre 2016 to SP; LABEX TULIP, ANR10-LABX-41 to MV]; [NDHET to MV]; the
AnaEE-France project [ANR-11-INBS-0001 to MV];
and the ‘‘Investissements d’Avenir’’ program.
Supplementary data
Supplementary Data available at ICB online.
References
Addo-Bediako A, Chown SL, Gaston KJ. 2000. Thermal tolerance, climatic variability and latitude. Proc Roy Soc Lond
B 267:739–45.
Andrewartha HG. 1944. Air temperature records as a guide to
the date of hatching of the nymphs of Austroicetes cruciata
Sauss. (Orthoptera). Bull Entomol Res 35:31–41.
Angilletta MJ. 2009. Thermal adaptation—a theoretical and
empirical synthesis. Oxford, UK: Oxford University Press.
58
Angilletta MJ, Jr, Wilson RS, Niehaus AC, Sears MW, Navas
CA, Ribeiro PL. 2007. Urban physiology: city ants possess
high heat tolerance. PLoS One 2:e258.
Araujo MB, Ferri-Yanez F, Bozinovic F, Marquet PA,
Valladares F, Chown SL. 2013. Heat freezes niche evolution.
Ecol Lett 16:1206–19.
Arnfield AJ. 2003. Two decades of urban climate research: a
review of turbulence, exchanges of energy and water, and
the urban heat island. Inter J Clim 23:1–26.
Ashcroft M, Gollan J, Warton D, Ramp D. 2012. A novel
approach to quantify and locate potential microrefugia
using topoclimate, climate stability, and isolation from
the matrix. Glob Change Biol 18:1866–79.
Bakken GS. 1989. Arboreal perch properties and the operative
temperature experienced by small animals. Ecology 70:922–30.
Bakken GS, Gates DM. 1975. Heat-transfer analysis of animals: some implications for field ecology, physiology, and
evolution. In: Gates DM, Schmerl RB, editors. Perspectives
in biophysical ecology Part IV, vol. 12. Springer Berlin
Heidelberg.
Bartholomew GA. 1966. Interaction of physiology and behavior under natural conditions. In: Bowman RI, editor.
The galapagos.Berkeley, CA: University of California
Press.
Barton BT, Schmitz OJ. 2009. Experimental warming transforms multiple predator effects in a grassland food web.
Ecol Lett 12:1317–25.
Benedetti-Cecchi L, Bertocci I, Vaselli S, Maggi E. 2006.
Temporal variance reverses the impact of high mean intensity of stress in climate change experiments. Ecology
87:2489–99.
Bernard F, Sache I, Suffert F, Chelle M. 2013. The development of a foliar fungal pathogen does react to leaf temperature!. New Phytol 198:232–40.
Boardman L, Terblanche JS. 2015. Oxygen safety margins set
thermal limits in an insect model system. J Exp Biol
218:1677–85.
Boehrer B, Schultze M. 2008. Stratification of lakes. Rev
Geophys 46:RG2005.
Brewitt K, Danner E. 2014. Spatio-temporal temperature variation influences juvenile steelhead (Oncorhynchus mykiss)
use of thermal refuges. Ecosphere 5:92.
Briscoe N, Handasyde K, Griffiths S, Porter WP,
Krockenberger A, Kearney MR. 2014. Tree-hugging koalas
demonstrate a novel thermoregulatory mechanism for arboreal mammals. Biol Lett 10:20140235.
Buckley LB, Tewksbury JJ, Deutsch CA. 2013. Can terrestrial
ectotherms escape the heat of climate change by moving?
Proc Roy Soc Lond B 280:20131149.
Buckley LB, Urban MC, Angilletta MJ, Crozier LG, Rissler LJ,
Sears MW. 2010. Can mechanism inform species’ distribution models? Ecol Lett 13:1041–54.
Caillon R, Suppo C, Casas J, Woods HA, Pincebourde S.
2014. Warming decreases thermal heterogeneity of leaf surfaces: implications for behavioural thermoregulation by arthropods. Func Ecol 28:1449–58.
Caissie D. 2006. The thermal regime of rivers: a review. Fresh
Biol 51:1389–406.
Chelle M. 2005. Phylloclimate or the climate perceived by
individual plant organs: what is it? How to model it?
What for? New Phytol 166:781–90.
S. Pincebourde et al.
Cloudsley-Thompson JL. 1962. Microclimates and the distribution of terrestrial arthropods. Ann Rev Entomol 7:199–222.
Clusella-Trullas S, Chown SL. 2011. Comment on ‘‘Erosion of
lizard diversity by climate change and altered thermal
niches’’. Science 332:537.
Cox TE, Smith CM. 2011. Thermal ecology on an exposed
algal reef: infrared imagery a rapid tool to survey temperature at local spatial scales. Coral Reefs 30:1109–20.
Crespo-Pérez V, Régnière J, Chuine I, Rebaudo F, Dangles O.
2015. Changes in the distribution of multispecies pest assemblages affect levels of crop damage in warming tropical
Andes. Glob Change Biol 21:82–96.
Dake JMK, Harleman DRF. 1969. Thermal stratification in
lakes: analytical and laboratory studies. Water Resour Res
5:484–95.
Dallas H, Rivers-Moore N. 2011. Micro-scale heterogeneity in
water temperature. Water SA 37:505–12.
Denny MW, Harley CDG. 2006. Hot limpets: Predicting body
temperature in a conductance-mediated thermal system.
J Exp Biol 209:2409–19.
Denny MW, Dowd WW, Bilir L, Mach KJ. 2011. Spreading
the risk: Small-scale body temperature variation among intertidal organisms and its implications for species persistence. J Exp Mar Biol Ecol 400:175–90.
Denny MW, Hunt L, Miller L, Harley CDG. 2009. On the
prediction of extreme ecological events. Ecol Monogr
79:397–421.
Diamond S, Dunn RR, Frank SD, Haddad NM, Martin RA.
2015. Shared and unique responses of insects to the interaction of urbanization and background climate. Cur
Opinion Ins Sci 11:71–7.
Dillon ME, Wang G, Huey RB. 2010. Global metabolic impacts of recent climate warming. Nature 467:704–6.
Dillon ME, Woods HA, Wang G, Fey SB, Vasseur DA,
Telemeco RS, Marshall K, Pincebourde S. Life in the frequency domain: the biological impacts of changes in climate variability at multiple time scales. Integr Comp Biol
this issue.
Doulos L, Santamouris M, Livada I. 2004. Passive cooling of
outdoor urban spaces. The role of materials. Solar Energy
77:231–49.
Dowd W, King F, Denny M. 2015. Thermal variation, thermal
extremes and the physiological performance of individuals.
J Exp Biol 218:1956–67.
Du Preez C, Tunnicliffe V. 2012. A new video survey method
of microtopographic laser scanning (MiLS) to measure
small-scale seafloor bottom roughness. Limnol Oceanogr:
Methods 10:899–909.
Dunckel A, Cardenas B, Sawyer A, Bennett P. 2009. Highresolution in-situ thermal imaging of microbial mats at El
Tatio Geyser, Chile shows coupling between community
color and temperature. Geophys Res Lett 36:L23403.
Faye E. 2015. Thermal landscapes and pest dynamic in
Andean tropical agrosystems [PhD dissertation] Paris:
Université Pierre et Marie Curie.
Faye E, Dangles O, Pincebourde S. 2016. Distance makes the
difference in thermography for ecological studies. J Therm
Biol 56:1–9.
Faye E, Rebaudo F, Yanez-Cajo D, Cauvy-Fraunié S, Dangles
O. 2015. A toolbox for studying thermal heterogeneity
across spatial scales: from unmanned aerial vehicle imagery
Thermal heterogeneity across space
to landscape metrics. Meth Ecol Evol doi: 10.1111/2041210X.12488.
Finnigan J. 2000. Turbulence in plant canopies. Ann Rev
Fluid Mech 32:519–71.
Garcı́a-Robledo C, Kuprewicz EK, Staines CL, Erwin TL,
Kress WJ. 2016. Limited tolerance by insects to high temperatures across tropical elevational gradients and the implications of global warming for extinction. PNAS doi:
10.1073/pnas.1507681113.
Gates DM. 1980. Biophysical ecology. New York, USA:
Springer-Verlag.
Geiger R. 1959. The climate near the ground. Cambridge, UK:
Harvard University Press.
Gilman S, Wethey D, Helmuth B. 2006. Variation in the sensitivity of organismal body temperature to climate change
over local and geographic scales. PNAS 103:9560–5.
Goller M, Goller F, French S. 2014. A heterogeneous
thermal environment enables remarkable behavioral
thermoregulation in Uta stansburiana. Ecol Evol
4:3319–29.
Grant BW, Dunham AE. 1988. Thermally imposed time constraints on the activity of the desert lizard Sceloporus merriami. Ecology 69:167–76.
Grant BW, Porter WP. 1992. Modeling global macroclimatic
constraints on ectotherm energy budgets. Am Zool
32:154–78.
Guo G, Wu Z, Xiao R, Chen Y, Liu X, Zhang X. 2015.
Impacts of urban biophysical composition on land surface
temperature in urban heat island clusters. Land Urban Plan
135:1–10.
Harlan SL, Brazel AJ, Prashad L, Stefanov WL, Larsen L. 2006.
Neighborhood microclimates and vulnerability to heat
stress. Soc Sci Med 63:2847–63.
Hawkins CP, Hogue JN, Decker LM, Feminella JW. 1997.
Channel morphology, water temperature, and assemblage
structure of stream insects. J N Am Benthol Soc 16:728–49.
Helmuth B. 1998. Intertidal mussel microclimates: Predicting
the body temperature of a sessile invertebrate. Ecol Monogr
68:51–74.
Helmuth B, Broitman BR, Blanchette C, Gilman SE, Halpin P,
Harley CDG, O’Donnell MJ, Hofmann GE, Menge B,
Strickland D. 2006. Mosaic patterns of thermal stress in
the rocky intertidal zone: implications for climate change.
Ecol Monogr 76:461–79.
Helmuth B, Kingsolver JG, Carrington E. 2005. Biophysics,
physiological ecology, and climate change: Does mechanism
matter? Ann Rev Physiol 67:177–201.
Hoffmann AA, Chown SL, Clusella-Trullas S. 2013. Upper
thermal limits in terrestrial ectotherms: how constrained
they are? Func Ecol 27:934–49.
Hondzo M, Stefan HG. 1993. Regional water temperature
characteristics of lakes subjected to climate change. Clim
Change 24:187–211.
Huey RB, Bennett AF. 2008. Bart’s familiar quotations: The
enduring biological wisdom of George A. Bartholomew.
Physiol Biochem Zool 81:519–25.
Huey RB, Kearney MR, Krockenberger A, Holtum JAM, Jess
M, Williams SE. 2012. Predicting organismal vulnerability
to climate warming: roles of behaviour, physiology and
adaptation. Phil Trans Royal Soc B 367:1665–79.
59
Huey RB, Peterson CR, Arnold SJ, Porter WP. 1989. Hot
rocks and not-so-hot rocks: retreat-site selection by garter
snakes and its thermal consequences. Ecology 70:931–44.
Huey RB, Pianka ER, Hoffmann J. 1977. Seasonal variation in
thermoregulatory behavior and body temperatures of diurnal Kalahari lizards. Ecology 58:1066–75.
Huey RB, Slatkin M. 1976. Cost and benefits of lizard thermoregulation. Quart Rev Biol 51:363–84.
IPCC. 2014. Climate change 2014: impacts, adaptation, and
vulnerability. Cambridge, UK: Cambridge University Press.
Jenerette GD, Harlan SL, Brazel A, Jones N, Larsen L,
Stefanov WL. 2007. Regional relationships between surface
temperature, vegetation, and human settlement in a rapidly
urbanizing ecosystem. Landsc Ecol 22:353–65.
Jones HG. 1999. Use of thermography for quantitative studies
of spatial and temporal variation of stomatal conductance
over leaf surfaces. Plant Cell Env 22:1043–55.
Kalnay E, Cai M. 2003. Impact of urbanization and land-use
change on climate. Nature 423:528–31.
Kaspari M, Clay NA, Lucas J, Yanoviak SP, Kay A. 2015.
Thermal adaptation generates a diversity of thermal limits
in a rainforest ant community. Glob Change Biol
21:1092–102.
Kearney MR. 2013. Activity restriction and the mechanistic
basis for extinctions under climate warming. Eco Lett
16:1470–9.
Kearney MR, Isacc AP, Porter WP. 2014. Microclim: Global
estimates of hourly microclimate based on long-term
monthly climate averages. Sci Data 1:140006.
Kearney MR, Shamakhy A, Tingley R, Karoly DJ, Hoffmann
AA, Briggs PR, Porter WP. 2014. Microclimate modelling at
macro scales: a test of a general microclimate model integrated with gridded continental-scale soil and weather data.
Meth Ecol Evol 5:273–86.
Kearney MR, Shine R, Porter WP. 2009. The potential for
behavioral thermoregulation to buffer ‘‘cold-blooded’’ animals against climate warming. PNAS 106:3835–40.
Kearney MR, Simpson SJ, Raubenheimer D, Helmuth B. 2010.
Modelling the ecological niche from functional traits. Phil
Trans Royal Soc B 365:3469–83.
Kingsolver JG, Moffat RJ. 1982. Thermoregulation and the
determinants of heat transfer in Colias butterflies.
Oecologia 53:27–33.
Kingsolver JG, Watt WB. 1983. Thermoregulatory strategies of
Colias butterflies: thermal stress and the limits to adaptation in temporally varying environments. Am Nat
121:32–55.
Koussoroplis AM, Wacker A. 2016. Covariance modulates the
effect of joint temperature and food variance on ectotherm
life-history traits. Ecol Lett 19:143–52.
Lacke MC, Mote TL, Shepherd JM. 2009. Aerosols and associated precipitation patterns in Atlanta. Atmos Environ
43:4359–73.
Lannig G, Cherkasov AS, Pörtner HO, Bock C, Sokolova IM.
2008. Cadmium-dependent oxygen limitation affects temperature tolerance in eastern oysters (Crassostrea virginica
Gmelin). Am J Physiol 294:1338–46.
Larsen L. 2015. Urban climate and adaptation strategies.
Front Ecol Env 13:486–92.
Lathlean J, Ayre D, Minchinton T. 2012. Using infrared imagery to test for quadrat-level temperature variation and
60
effects on the early life history of a rocky-shore barnacle.
Limnol Oceano 57:1279–91.
Leuzinger S, Körner C. 2007. Tree species diversity affects
canopy leaf temperatures in a mature temperate forest.
Agri For Meteor 146:29–37.
Li X, Zhou W, Ouyang Z. 2013. Relationship between land
surface temperature and spatial pattern of greenspace: what
are the effects of spatial resolution? Land Urban Plan
114:1–8.
Linacre ET. 1967. Further notes on a feature of leaf and air
temperatures. Theor Appl Climatol 15:422–36.
Maitani T. 1979. A comparison of turbulence statistics in the
surface layer over plant canopies with those over several
other surfaces. Bound Layer Meteorol 17:213–22.
Meineke EK, Dunn RR, Sexton JO, Frank SD. 2013. Urban
warming drives insect pest abundance on street trees. PLoS
One 8:1–7.
Mislan KAS, Wethey DS, Helmuth B. 2009. When to worry
about the weather: role of tidal cycle in determining patterns of risk in intertidal ecosystems. Glob Change Biol
15:3056–65.
Mortimer CH. 1952. Water movements in lakes during
summer stratification; evidence from the distribution of
temperature in Windermere. Phil Trans Roy Soc Lond B
236:355–98.
Moss B. 1969. Vertical heterogeneity in the water column of
Abbot’s pond: I. The distribution of temperature and dissolved oxygen. J Ecol 57:381–96.
Mota C, Engelen A, Serrão E, Pearson G. 2015. Some
don’t like it hot: Microhabitat-dependent thermal and
water stresses in a trailing edge population. Func Ecol
29:640–49.
Mott KA, Buckley TN. 2000. Patchy stomatal conductance:
Emergent collective behavior of stomata. Trends Plant Sci
5:258–62.
Muth A. 1977. Body temperatures and associated postures of
the zebra-tailed lizard, Callisaurus draconoices. Copeia
1:122–5.
Muth A. 1980. Physiological ecology of desert iguana
(Dipsosaurus dorsalis) eggs: temperature and water relations. Ecology 61:1335–43.
Niyogi D, Pyle P, Lei M, Arya SP, Kishtawal CM, Shepherd
M, Chen F, Wolfe B. 2011. Urban modification of thunderstorms: An observational storm climatology and model
case study for the Indianapolis urban region. J App
Meteorol Climatol 50:1129–44.
Oke TR. 1982. The energetic basis of the urban heat island.
Quart J Royal Meteorol Soc 108:1–24.
Oke TR. 1989. Boundary layer climates. London: Routledge.
Peglow JR. 2013. Investigating the relationships between community structure and rugosity in the rocky intertidal using
terrestrial laser scanning [Master Thesis]. Faculty of the
Division of Science and Environmental Policy, California
State University, Monterey Bay CA: Monterey Bay.
Pike DA, Webb JK, Shine R. 2012. Hot mothers, cool eggs:
nest-site selection by egg-guarding spiders accommodates
conflicting thermal optima. Funct Eco l26:469–75.
Pincebourde S, Casas J. 2006. Multitrophic biophysical budgets: thermal ecology of an intimate herbivore insect–plant
interaction. Ecol Monogr 76:175–94.
S. Pincebourde et al.
Pincebourde S, Casas J. 2015. Warming tolerance across insect
ontogeny: influence of joint shifts in microclimates and
thermal limits. Ecology 96:986–97.
Pincebourde S, Casas J. 2016. Hypoxia and hypercarbia in
endophagous insects: Larval position in the plant gas exchange network is key. J Insect Physiol 84:137–53.
Pincebourde S, Suppo C. The vulnerability of tropical ectotherms to warming is modulated by the microclimatic heterogeneity. Integr Comp Biol this issue.
Pincebourde S, Sinoquet H, Combes D, Casas J. 2007.
Regional climate modulates the canopy mosaic of favourable and risky microclimates for insects. J Anim Ecol
76:424–38.
Pincebourde S, Woods HA. 2012. Climate uncertainty on leaf
surfaces: the biophysics of leaf microclimates and their consequences for leaf-dwelling organisms. Func Ecol
26:844–53.
Porter WP, Gates DM. 1969. Thermodynamic equilibria of
animals with environment. Ecol Monogr 39:227–44.
Porter WP, Sabo JL, Tracy CR, Reichman OJ, Ramankutty N.
2002. Physiology on a landscape scale: Plant-animal interactions. Integr Comp Biol 43:431–453.
Potter KA, Woods HA, Pincebourde S. 2013. Microclimatic
challenges in global change biology. Glob Chang Biol
19:2932–9.
Rhee J, Park S, Lu Z. 2014. Relationship between land cover
patterns and surface temperature in urban areas. GIScience
Rem Sens 51:521–36.
Rice WR, Gaines SD. 1994. The ordered-heterogeneity family
of tests. Biometrics 50:746–52.
Roth-Nebelsick A. 2001. Computer-based analysis of steady
state and transient heat transfer of small-sized leaves by
free and mixed convection. Plant Cell Env 24:631–40.
Saudreau M, Marquier A, Adam B, Sinoquet H. 2011.
Modelling fruit-temperature dynamics within apple tree
crowns using virtual plants. Ann Bot 108:1111–20.
Saudreau M, Pincebourde S, Dassot M, Adam B, Loxdale HD,
Biron DG. 2013. On the canopy structure manipulation to
buffer climate change effects on insect herbivore development. Trees 27:239–48.
Scherrer D, Körner C. 2010. Infra-red thermometry of alpine
landscapes challenges climatic warming projections. Glob
Change Biol 16:2602–13.
Scherrer D, Körner C. 2011. Topographically controlled thermal-habitat differentiation buffers alpine plant diversity
against climate warming. J Biogeo 38:406–16.
Seabra R, Wethey DS, Santos AM, Lima FP. 2011. Side matters: microhabitat influence on intertidal heat stress over a
large geographical scale. J Exp Mar Biol Ecol 400:200–8.
Sears MW, Angilletta MJ. Jr. 2015. Costs and benefits of thermoregulation revisited: both the heterogeneity and spatial
structure of temperature drive energetic costs. Am Nat
185:E94–E102.
Sears MW, Raskin E, Angilletta MJ. Jr. 2011. The world is not
flat: defining relevant thermal landscapes in the context of
climate change. Integr Comp Biol 51:666–75.
Sheldon KS, Dillon ME. Beyond the mean: biological impacts
of cryptic temperature change. Integr Comp Biol this issue.
Shepherd JM, Carter M, Manyin M, Messen D, Burian S.
2010. The impact of urbanization on current and future
Thermal heterogeneity across space
coastal precipitation: a case study for Houston. Env Plan BPlan Design 37:284–304.
Sinervo B, et al.. 2010. Erosion of lizard diversity by climate
change and altered thermal niches. Science 328:894–99.
Sinoquet H, Le Roux X, Adam B, Amaglio T, Daudet FA.
2001. RATP: a model for simulating the spatial distribution
of radiation absorption, transpiration and photosynthesis
within canopies: application to an isolated tree crown.
Plant Cell Env 24:395–406.
Sinoquet H, Pincebourde S, Adam B, Donès N,
Phattaralerphong J, Combes D, Ploquin S, Sangsing K,
Kasemsap P, Thanisawanyangkura S, et al.. 2009. 3-D
maps of tree canopy geometries at leaf scale. Ecology 90:283.
Snyder GK, Weathers WW. 1975. Temperature adaptations in
amphibians. Am. Nat 109:93–101.
Stabler LB, Martin CA, Brazel AJ. 2005. Microclimates in a
desert city were related to land use and vegetation index.
Urban For Urban Green 3:137–47.
Stein A, Gerstner K, Kreft H. 2014. Environmental heterogeneity as a universal driver of species richness across taxa,
biomes and spatial scales. Ecol Lett 17:866–80.
Stevenson RD. 1985. The relative importance of behavioral
and physiological adjustments controlling body temperature in terrestrial ectotherms. Am Nat 126:362–86.
Stillman J, Somero G. 1996. Adaptation to temperature stress
and aerial exposure in congeneric species of intertidal porcelain crabs (genus Petrolisthes): correlation of physiology,
biochemistry and morphology with vertical distribution.
J Exp Biol 199:1845–55.
Stokes VJ, Morecroft MD, Morison JIL. 2006. Boundary layer
conductance for contrasting leaf shapes in a deciduous
broadleaf forest canopy. Agr For Meteor 139:40–54.
Suggitt A, Gillingham P, Hill J, Huntley B, Kunin W, Roy D,
Thomas C. 2011. Habitat microclimates drive fine-scale
variation in extreme temperatures. Oikos 120:1–8.
Sunday J, Bates AE, Dulvy NK. 2011. Global analysis of thermal tolerance and latitude in ectotherms. Proc Royal Soc
Lond B 278:1823–30.
Sunday JM, Bates AE, Kearney MR, Colwell RK, Dulvy NK,
Longino JT, Huey RB. 2014. Thermal safety margins and
the necessity of thermoregulatory behavior across latitude
and elevation. PNAS 111:5610–5.
Todhunter PE. 1990. Microclimates variations attributable to
urban-canyon asymmetry and orientation. Phys Geogr
11:131–41.
Tonolla D, Acuna V, Uehlinger U, Frank T, Tockner K. 2010.
Thermal heterogeneity in floodplains. Ecosystems 13:727–40.
61
Tracy CR. 1976. A model of the dynamic exchanges of water
and energy between a terrestrial amphibian and its environment. Ecol Monogr 46:293–326.
USGCRP (US Global Change Research Program). 2014.
National climate assessment. Washington, DC: USGCRP.
Vasseur DA, DeLong JP, Gilbert B, Greig HS, Harley CDG,
McCann KS, Savage V, Tunney TD, O’Connor MI. 2014.
Increased temperature variation poses a greater risk to
species than climate warming. Proc Roy Soc B Lond
281:20132612.
Verberk WCEP, Calosi P. 2012. Oxygen limits heat tolerance
and drives heat hardening in the aquatic nymphs of the gill
breathing damselfly Calopteryx virgo (Linnaeus, 1758).
J Therm Biol 37:224–29.
Vickers M, Manicom C, Schwarzkopf L. 2011. Extending the
cost-benefit model of thermoregulation: high-temperature
environments. Am Nat 177:452–61.
Vierling KT, Vierling LA, Gould WA, Martinuzzi S,
Clawges RM. 2008. Lidar: shedding new light on habitat characterization and modeling. Front Ecol Environ
6:90–98.
Wawrzyniak V, Piégay H, Allemand P, Vaudor L, Grandjean
P. 2013. Prediction of water temperature heterogeneity of
braider rivers using very high resolution thermal infrared
(TIR) images. Intern J Remote Sens 34:4812–31.
Woods HA. 2010. Water loss and gas exchange by eggs of
Manduca sexta: trading off costs and benefits. J Insect
Physiol 56:480–7.
Woods HA. 2013. Ontogenetic changes in the body temperature of an insect herbivore. Func Ecol 27:1322–31.
Woods HA, Dillon EM, Pincebourde S. 2015. The roles of
microclimatic diversity and of behaviour in mediating the
responses of ectotherms to climate change. J Therm Biol
54:86–97.
Youngsteadt E, Dale AG, Terando AJ, Dunn RR, Frank SD.
2015. Do cities simulate climate change? A comparison of
herbivore response to urban and global warming. Glob
Change Biol 21:97–105.
Zhang C, Wu J, Grimm NB, McHale M, Buyantuyev A. 2013.
A hierarchical patch mosaic ecosystem model for urban
landscapes: model development and evaluation. Ecol Mod
250:81–100.
Zhou W, Huang G, Cadenasso ML. 2011. Does spatial configuration matters? Understanding the effects of land cover
pattern on land surface temperature in urban landscapes.
Land Urban Plan 102:54–63.