Climate change velocity underestimates climate change exposure in mountainous regions

Abstract

Climate change velocity is a vector depiction of the rate of climate displacement used for assessing climate change impacts. Interpreting velocity requires an assumption that climate trajectory length is proportional to climate change exposure; longer paths suggest greater exposure. However, distance is an imperfect measure of exposure because it does not quantify the extent to which trajectories traverse areas of dissimilar climate. Here we calculate velocity and minimum cumulative exposure (MCE) in degrees Celsius along climate trajectories for North America. We find that velocity is weakly related to MCE; each metric identifies contrasting areas of vulnerability to climate change. Notably, velocity underestimates exposure in mountainous regions where climate trajectories traverse dissimilar climates, resulting in high MCE. In contrast, in flat regions velocity is high where MCE is low, as these areas have negligible climatic resistance to movement. Our results suggest that mountainous regions are more climatically isolated than previously reported.

ARTICLE

Received 16 Feb 2016 | Accepted 23 Jun 2016 | Published 1 Aug 2016

Climate change velocity is a vector depiction of the rate of climate displacement used for assessing climate change impacts. Interpreting velocity requires an assumption that climate trajectory length is proportional to climate change exposure; longer paths suggest greater exposure. However, distance is an imperfect measure of exposure because it does not quantify the extent to which trajectories traverse areas of dissimilar climate. Here we calculate velocity and minimum cumulative exposure (MCE) in degrees Celsius along climate trajectories for North America. We nd that velocity is weakly related to MCE; each metric identies contrasting areas of vulnerability to climate change. Notably, velocity underestimates exposure in mountainous regions where climate trajectories traverse dissimilar climates, resulting in high MCE. In contrast, in at regions velocity is high where MCE is low, as these areas have negligible climatic resistance to movement. Our results suggest that mountainous regions are more climatically isolated than previously reported.

DOI: 10.1038/ncomms12349 OPEN

Climate change velocity underestimates climate change exposure in mountainous regions

Solomon Z. Dobrowski1 & Sean A. Parks2

1 Department of Forest Management, College of Forestry and Conservation, University of Montana. Missoula, Montana 59812, USA. 2 Aldo Leopold Wilderness Research Institute, Rocky Mountain Research Station, US Forest Service. Missoula, Montana 59812, USA. Correspondence and requests for materials should be addressed to S.Z.D. (email: mailto:[email protected]

Web End [email protected] ) or to S.A.P. (email: mailto:[email protected]

Web End [email protected] )

NATURE COMMUNICATIONS | 7:12349 | DOI: 10.1038/ncomms12349 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12349

Climate change impacts on biota depend in part on climate change exposure, the degree to which a system is exposed to climate variations over time or space1,2.

One of the most widely used metrics for estimating exposure is climate change velocity, the direction and rate at which organisms must move to maintain a given climate3. Multiple formulations have been proposed for calculating climate velocity, including approaches based on local climate gradients3 (herein gradient based) and those based on distance to analogue climates (distance based)4. These metrics have been used to describe risk of species extinctions5, climate change responses of marine taxa6, regional patterns in species endemism7, Quaternary range shifts of biota8, the distribution of climate change refugia9 and climate change exposure under observed and projected future conditions3,4,1012. Velocity has also been used to derive climate trajectories, paths that describe the movement of a climate isopleth over a given period of time13. Although velocity estimates are dependent on a number of methodological choices4,10,11, the relative rankings of velocity have been shown to be quite robust to differences in methodology4 and show a general pattern of low values (implying longer climate residence times) in regions with high spatial climate heterogeneity (for example, mountainous regions), and high velocity (shorter residence times) in areas of low topographic relief 3,10,12. Consequently, authors have suggested that at areas are particularly exposed to climate change while areas of complex terrain may act to buffer climate change impacts by allowing organisms to ameliorate climate shifts via short distance dispersal10,1416.

Implicit in the use of velocity is the assumption that climate trajectory length is proportional to exposure; longer paths suggest greater exposure3,4,10,12,17. However, distance is an imperfect measure of exposure in the same way it is an insufcient measure of connectivity18,19. Climatic connectivitythe ability of a landscape to promote or hinder species movement in response to a changing climateis also contingent on the costs of moving through areas of dissimilar climate. These costs may preclude migration in areas with strong climate gradients. For instance, the closest climate analogue for a mountain-top may be found on an adjacent mountain-top. The distance may be short between each

point, but the actual exposure to climate differences along the trajectory through the ensuing valley can be large (Fig. 1). In addition, exposure will depend on assumptions about dispersal preferences of organisms. Current distance-based approaches4,11 to estimating velocity assume that organisms will minimize the distance they travel (Euclidean distanceED) as opposed to minimizing their exposure to dissimilar climate. If organisms minimize exposure and thus deviate from straight line paths, then velocity based on ED (velocityED) will underestimate rates needed to keep pace with changing climate (Fig. 1).

Here we propose a set of modications for distance-based velocity based on the assumption that organisms will follow paths that minimize their exposure to dissimilar climates. The distance along this trajectory from source to destination pixel (minimum exposure distance; MED) will be greater than or equal to distances calculated using ED. To account for this assumption, we employ least-cost modelling techniques (see Methods section) applied to isotherms within North America (sensu20) for the interval of 19952085. We contrast velocityED with velocity estimates based on MED (velocityMED). In addition, we quantify the minimum cumulative exposure (MCE) in C along each climate trajectory which to our knowledge, is a previously unreported facet of exposure. We nd that velocity is weakly related to MCE; both velocityED and velocityMED underestimate exposure in mountainous regions where even short climate trajectories traverse landscapes with dissimilar climate resulting in high values of MCE. In contrast, velocity is high over at regions of the continent where MCE is low, as these areas have negligible climatic resistance to movement.

ResultsVelocity and Minimum Cumulative Exposure. For North America, the average (geometric mean) climate velocity based on MED (velocityMED) was 3.25 km per year (595th percentile;

0.2830.9 km per year). In comparison, average velocity based on ED (velocityED) was 2.68 km per year (0.2411.8 km per year), which indicates that velocityMED is 21% higher on average than velocityED. The location of the nearest climate analogue varies depending on whether we assume organisms minimize exposure

a b

End

End-2

MED=397 km

Velocity=4.4 km per year

ED=294 km

Velocity=3.3 km per year

MCE: 83.5 C

a b

Path 2

End-1

Path 1

Start-2

MCE: 0 C

Start-1

Mean annual temp. (current period)

Low

High

Start

Figure 1 | Climate trajectories and minimum cumulative exposure. Climate trajectories are dened by a source pixel (start) with a given temperature under current conditions (19812010) and a destination pixel (end) with a similar temperature under future conditions (20712100). An example of the difference between climate trajectories dened by Euclidean distance (ED) and minimum exposure distance (MED) for a site within North America (inset) is presented in a. MED (solid line) minimizes exposure to climates that are dissimilar to that of the source pixel whereas ED (dashed line) minimizes distance travelled. An example of minimum cumulative exposure (MCE) for two trajectories is presented in b. Although the trajectories are similar in length, MCE 0 for path 1 because the trajectory tracks changes in climate over the study period without traversing temperatures that are dissimilar to that

of the source pixel. Path 2 has a high MCE because the trajectory traverses two warm valleys (shading).

2 NATURE COMMUNICATIONS | 7:12349 | DOI: 10.1038/ncomms12349 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12349 ARTICLE

to dissimilar climates or minimize distance travelled; 37% of source pixels had different destination locations when comparing the MED and ED approaches; of these pixels, the locations differed on average (geometric mean) by 236 km. MCE values were strongly right skewed; 71% of the continent exhibited no exposure (MCE 0) despite having widely varying velocityMED. For

cells in which MCE40, the median value was 147 C (3.0 8,005 C). In addition to the results presented herein, which are based on gridded temperature data at 5 km resolution, we also assessed the sensitivity of velocityMED, velocityED and MCE to the resolution of input climate data, the width of temperature bins for dening a climate analogue, and cost penalties associated with the resistance surfaces used in least-cost modelling (see Methods section, Supplementary Fig. 1, Supplementary Fig. 2, and Supplementary Table 1).

There is considerable regional variability in velocityMED and

MCE (summarized in Supplementary Table 2) driven by the rate of temperature change, regional landforms, proximity to water and local physiography. VelocityMED tended to increase with latitude due to greater projected warming, was lower in areas of complex terrain, and was higher in at areas and regions surrounded by large water bodies such as the Great Lakes (Fig. 2, Supplementary Fig. 3). In contrast, MCE values were generally 0 in at regions with high velocities (the Great Plains and Boreal region) and were often 40 in topographically complex areas along the major mountain cordillera of the continent (areas with generally low velocity; Fig. 3). Areas with large MCE values often had climate trajectories that required southward movements (for example, northward oriented peninsulas) or traversed relatively warm valleys (mountainous regions; Figs 1b and 3, Supplementary Fig. 3). The ratio between velocityMED and velocityED was also higher in mountainous regions and on northward oriented peninsulas (Fig. 2). VelocityMED and MCE are only weakly correlated (Spearmans

rank r 0.37); for MCE values 40 this correlation increases

(r 0.81), although there is substantial variation in velocityMED

across all MCE values (Fig. 4).

DiscussionVelocity and MCE describe complementary facets of exposure: the length of potential migration paths and climatic resistance to movement along these paths through time. Our estimates of velocityMED and velocityED mirror existing studies which suggest that mountainous regions have relatively low climate change

exposure (longer climate residence times)3,10,12,14. However, this interpretation assumes climatic resistance to movement does not vary in space. We show that in areas of complex terrain (for example, the cordillera of North America), MCE can be high even for short climate trajectories because they often traverse climatically heterogeneous landscapes (Figs 1b and 3). In this context, spatial climate heterogeneity acts to increase resistance to movement, thus reducing climatic connectivity between a source and destination pixel. In contrast, long climate trajectories over much of the continent pass through climatically homogenous environments with little to no resistance. Under these conditions, velocity is large but MCE is zero (that is, much of the Great Plains and Canadian Boreal region; Fig. 3).

Climate trajectories based on ED and MED dene two end-members of a spectrum of potential trajectories that are contingent on a sites physiographic setting and the traits of species that occupy the site. In at regions, for example, there is little difference between these end-members; velocityMED and

velocityED show similar values (Fig. 2). In contrast, velocityMED and velocityED diverge most notably in areas of complex terrain (Fig. 1a and 2). Furthermore, MCE represents the minimum exposure to dissimilar climate along a potential migration path. The actual climate change exposure of individual organisms will also depend on the traits of the species that occupy the site21,22. Vagile organisms may be able to keep pace with large climate velocities, whereas sessile organisms may fall behind changing climate23, resulting in increased exposure as climate change outpaces movements. Moreover, dispersal traits of organisms will likely affect migration routes. For example, volant species such as birds may follow straight line paths that are well approximated by velocityED. Thermal tolerances may also inuence whether an organism tends to follow straight line paths or paths that minimize exposure to dissimilar climates. For instance, thermal specialists (for example, ectotherms from the tropics;) are more likely to minimize their exposure to dissimilar climate (approximated by velocityMED) as compared with thermal generalists (for example, endotherms from high latitudes)22,24,25.

Similarly, the narrow reproductive niche2628 of many plant species will likely require longer migratory paths than those described using velocityED.

The bivariate distribution of velocityMED and MCE provides a logical framework for assessing exposure as a function of species dispersal capacity and thermal tolerance (Fig. 4a). We present a simple classication of the bivariate distribution of velocityMED

and MCE as four groups (Fig. 4b): low velocity, low MCE; low

a b c

km per year Ratio

0.1 1 10 49 1 1.25 1.5+

Figure 2 | Climate velocity for North America for the interval 1995 to 2085. Climate velocity based on minimum exposure distance (velocityMED) is

presented in a. The ratio between velocityMED and velocity based on Euclidean distance (velocityED) is presented in b. An inset (black box) of b resampled to a resolution of 1 km to improve visual interpretation is presented in c. Areas in grey represent pixels whose source or destination pixel is located on an island (see Methods section).

NATURE COMMUNICATIONS | 7:12349 | DOI: 10.1038/ncomms12349 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12349

a b

MCE C

1

100 2,500+

Figure 3 | Minimum cumulative exposure. MCE for North America for the interval 1995 to 2085 is presented in a. An inset (black box) resampled to a resolution of 1 km to improve visual interpretation is presented in b. MCE values40 are shown in colour; MCE 0 is shown in light gray. MCE values for

pixels on islands and pixels whose future climate analogues are on islands are shown in dark grey (see Methods section).

a Increasing exposure for poorly dispersing organisms

b

MCE (C)

5,000

500

50

Increasing exposure for

thermal specialists

5

0

0.05 0.2 1 5 20

Velocity (km per year)

Low velocity low MCE

Low velocity high MCE

High velocity low MCE

High velocity high MCE

Figure 4 | Classication of velocity and minimum cumulative exposure. The relationship between velocityMED and minimum cumulative exposure (MCE) for a random sample of 5000 pixels from North America is presented in a. VelocityMED and MCE values40 are plotted on log-log axis. MCE values 0

(71% of pixels) are plotted in lower panel and are considered low; MCE40 are considered high. The dashed vertical line represents the median velocityMED

which is used to delineate between low and high velocity classes. Map showing the distribution of the four classes based on the bivariate distribution of MCE and velocityMED is presented in b. Dot colour in a corresponds to map color in b. Areas in grey represent pixels whose source or destination pixel is located on an island (see Methods section).

velocity, high MCE; high velocity, low MCE; and high velocity, high MCE. We use the terms low and high to simply qualify the range of our data and do not suggest that the specic thresholds we use (see methods) have biological relevance. Nevertheless, areas with high velocity present increased risk of exposure to poorly dispersing species. Areas of high MCE present increased risk for thermal specialists. Of particular concern are areas with both high velocity and MCE (for example, mountain tops, peninsulas). These areas are climatically isolated and increase the vulnerability of local populations to climate change. Specic examples include the Great Lakes, Yucatan Moist Forests, and the Northern Appalachian/Acadian ecoregions (Supplementary Table 2, Supplementary Fig. 4).

The use of climate velocity and MCE should be done in a manner that is cognizant of the tension between species-specic approaches, which have greater data needs but account for individual species responses to climate change, and coarse-lter approaches that are applied independently of species-specic data,

such as the metrics presented here. Distance-based velocity implicitly assumes negligible thermal tolerances or adaptive capacity of local populations and thus represents an upper limit for migration requirements for climate-sensitive biota11. In contrast, sites comprised of organisms with more broadly dened climate tolerances will have lower MCE and velocity values because more of the neighbouring landscape will be suitable for these organisms. This is demonstrated in Supplementary Table 1; mean velocity and MCE decrease as climate bin width increases. Fine-lter approaches that model potential migration paths for individual species have also been developed. For example, circuit theory has been applied to modelling faunal movement routes that track climate changes29. In addition, both distance- and gradient-based approaches for calculating velocity have been adapted to account for individual species climate tolerances; so-called biotic velocity estimates the direction and rate at which a given species must move to track its climatic niche11,30. Biotic velocity suffers from the same limitations we identify here;

4 NATURE COMMUNICATIONS | 7:12349 | DOI: 10.1038/ncomms12349 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12349 ARTICLE

it may not reect the extent to which potential migration routes expose species to dissimilar climate. Nonetheless, the methods we employ can be applied to the calculation of biotic velocity and for estimating species-specic MCE.

VelocityMED, velocityED and MCE are also sensitive to the resolution of input data and to input parameters used in their calculation. VelocityMED, velocityED and MCE increase as data resolution is coarsened whereas they decline as the bin width for dening a climate analogue is broadened (Supplementary Table 1, Supplementary Fig. 1, Supplementary Fig. 2). These results are consistent with previous studies that have examined the scale sensitivity of velocity calculations4,10 and suggest that comparisons between studies and comparisons of species dispersal rates against climate velocities (for example, ref. 23) must be interpreted with explicit consideration to scale of analysis.

Mountains support roughly a quarter of the globes terrestrial biodiversity, contain 32% of protected areas and nearly half of the worlds biodiversity hotspots31. In light of this, our results raise a logical question of how to reconcile the notion of mountains as refugia under changing climate versus mountains as areas of high climatic resistance to movement. Climate change adaptation strategies are often predicated on the former notion which assumes that spatial variability in climate may allow for short distance dispersal to ameliorate climate change impacts1416,32. This paradigm underlies approaches for reserve design33,34 and methods for identifying microrefugia9,16,35,36. Nonetheless, an important and sometimes overlooked consideration is whether sites are accessible to migrating organisms37. Accessibility is partly a function of climatic connectivity to a given site4. For example, velocityMED and MCE calculated for reverse climate trajectories (sensu reverse velocity4,11) provide two measures of climatic accessibility (Fig. 5). Accessibility will also depend on physical barriers to movement, habitat fragmentation, land use and land cover change, among other factors. Temporal scale is another important consideration when trying to reconcile

these two viewpoints. VelocityMED and MCE depend in part on the magnitude of climate shifts (a function of time) relative to the range or extent of the local spatial gradient in climate. Over short time intervals (or in places with sharp climate gradients) local dispersal, within the range of a local monotonic spatial climate gradient, may allow an organism to keep pace with climate change without exposure to dissimilar climates (for example, simple upslope movement; Fig. 1b, path 1). However when the magnitude of climate change exceeds the extent of the local spatial gradient in climate, MCE (and potentially MED) will increase sharply as organisms are required to traverse unsuitable climates en route to a future climate analogue (Fig. 1b, path 2). This implies that at certain spatial and temporal scales, montane sites will act as climate refugia. Beyond these spatio-temporal domains, these sites will provide only temporary holdout habitat (sensu36).

Climate velocity estimates (whether based on MED or ED) for montane regions suggest limited exposure and relatively short migration paths under climate change. However, we show that distance (and velocity) is an imperfect measure of climate connectivity given that it does not account for spatial variability in climatic resistance to movement (that is, MCE). MCE can be high in mountainous regions, which implies that these areas are more climatically isolated than has been previously reported. Conversely, at regions of the continent have high velocities but negligible climatic resistance to movement. These regions may have greater climatic connectivity to their future climate analogues than is currently appreciated. Consequently, a more nuanced assessment of climate exposure is warranted, one which explicitly accounts for climatic resistance to movement. Hence, we advocate for further integration between climate informed connectivity modelling approaches and the development of climate change exposure metrics. This integration will more realistically describe multiple facets of exposure to biota from ongoing climate change.

a b

Figure 5 | Climate trajectories for Yellowstone National Park USA. Forward (a) and reverse (b) climate trajectories for Yellowstone National Park (red polygon) for the period 19952085. Forward trajectories depict routes that minimize the cumulative exposure (MCE) to dissimilar isotherms along paths between a source pixel and its future climate analogue (dot). Forward trajectories (sensu forward velocity4,11) characterize the potential exposure of organisms to climate change. Reverse climate trajectories depict paths with the lowest MCE between a future cell and an analogue under current climate conditions (dot). Reverse trajectories (sensu reverse velocity4,11) characterize how climatically accessible a site is to colonization from other sites. Climate trajectories (forward and reverse) can be characterized by minimum exposure distance (MED), velocityMED (length/time), and MCE.

NATURE COMMUNICATIONS | 7:12349 | DOI: 10.1038/ncomms12349 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12349

Methods

Climate data. We generated estimates of velocity and MCE over the interval (19812010) and (20712100) using mean annual temperature data described by (ref. 38). Projections of future mean annual temperature are based on an ensemble of 15 CMIP5 GCMs under the RCP 8.5 scenario. These data came rounded to the tenth of a degree and had an original spatial resolution of 1 km, which we resampled to a 5 km resolution.

For each 5 km pixel in North America, we identied all pixels with a climate analogue in the future time period. We dene the climate analogue as any pixel that is 0.25 C from the pixel of interest, effectively setting a bin width of 0.5 C. We then used a least-cost algorithm39 to identify the trajectory (potential migration routes) that minimizes exposure to dissimilar climates between each source and destination pixel with a matching future climate analogue. We developed cost surfaces (one for each 0.1 C temperature increment) in a manner that accounts for spatio-temporal changes in climate, allowing isotherms to be traced over landscapes and through time.

Least-cost modelling. The overall workow for our approach to least-cost modelling is presented in Fig. 6. To account for a changing climate, we rst interpolated mean annual temperature grids to represent incremental changes in temperature over the 90 year time period. These interpolated temperature grids (n 31) represent mean annual temperature on about three year increments

between the start (1995) and end date (2085) and assume a linear trend in changes in temperature over this time period (Fig. 6a). Next, we created intermediate cost surfaces (Fig. 6b,c) for each of these interpolated temperature grids such that costs were proportional to the level of dissimilarity between each 0.1 C increment of

interest (in 1995) and all other pixels based on the following equation:

cost i

1 p T T i

1 where costi is the cost assigned to pixel i, P is a dissimilarity penalty (P 2

dimensionless units per oC in this study such that costi increases by 1 for each0.5 C in climate dissimilarity), T is the 0.1 C temperature increment of interest, and Ti is the temperature of pixel i; a value of one is added because least-cost approaches do not allow costs of zero and it ensures that costi distance when

trajectories pass through pixels with an analogous climate. We then computed, among all intermediate cost surfaces, the minimum cost for each pixel (Fig. 6bd); this minimum cost is used to generate the nal cost surfaces (Fig. 6d). Cost surfaces were developed for each 0.1 C temperature increment (n 567 cost surfaces). For

each increment, climate analogues were dened based on bins that were 0.5 C wide so that matching pixels (0.25 C) were assigned a cost of one and cost values increased with the level of dissimilarity as integers (for example,5.3 C0.250.75 C was assigned a cost of two; see Fig. 6, Table 1). We evaluated every 0.1 C increment, yet used a bin width of 0.5 C to dene climate analogues. This was intended to reduce boundary effects; that is, we wanted to ensure that pixels with temperatures of 4.9 and 5.1 C were considered climate analogues as opposed to being treated as separate bins (Table 1). In the nal cost surfaces, lakes and ocean were arbitrarily assigned a cost value of 5,000 to heavily penalize open water, thereby forcing trajectories to avoid water when possible (Supplementary Fig. 3). Excluding water, cost values generally ranged from 1 to B100, reecting the B50 C difference in mean annual temperature across North America.

Accounting for temporal changes in climate (Fig. 6ad) was necessary for two reasons. First, had we based our cost surface on the mean annual temperature in

Mean annual temperature

Intermediate cost surfaces (5.0 C)

Intermediate cost surfaces (5.5 C)

2085

2040

5.0

a c

1

b

2

1

2

6.0

6.5 7.0

5.5

6.0

7.0

8.0

3

5

4

6.0 7.0

6.5 8.0

7.0 7.0

3 4 5

3

2

2 3 4

2

7.5

6

5

4

7.0

5 5 4 5

5

7

6

7.0 7.5

7.5 7.5

5

7

6 6

3

4 4 3 4

6

5 5

8.0 8.5 8.5

8.5 9.0

6 7 8 8

8 9

5 6 7 7

7 8

3.5

3.5

4

4

5

5

4.5

5.0 5.5

4.5

5.5

6.5

3

4.5 5.5

5.0 6.5

5.5 5.5

2

2

1

2 1 2

2

3 2 1

3

6.0

3

2

1

5.5

2 2 1 2

2

2

5.5 6.0

6.0 6.0

4

3 3

4

1

1 1 2 1

3

3

2 2

6.5 7.0 7.0

7.0 7.5

3 4 5 5

5 6

2 3 4 4

4 5

1995

2.0

7

7

8

8

3.0

3.5 4.0

2.0

3.0

4.0

5.0

5

3

5

6

4

3.0 4.0

3.5 5.0

4.0 4.0

5 4 3

6 5 4

6

2

3

3

4

4.0

4.5

3

3 3 4 3

1

4.0 4.5

2 2

4

4 4 5 4

3 3

2

4.5 4.5

5.5 5.5 6.0

2

1

2 2 3

3

2

1 1 2

2

1

5.0

5.5

1

2

Minimum costof three timeperiods = 2 Minimum cost of three time periods = 2

d

Minimum cost of three time periods = 2

Final cost for 5.0 C

Final cost for 5.5 C

2

1

Minimum cost of three time periods = 3

2

Cost surfaces for two temperatures

2

2

2 1 2

2

1

2

1

2

2 2 1

2

2

2 1 1 2

1

1

2

2

1

2 2

1

1 1 2 1

2

2

3 1

2 1 2 2

2 3

2 2 1 1

1 2

e

Target pixel (climate analogue)(5.0 C in 2085)

Target pixel (climate analogue)(5.5 C in 2085)

Climate trajectory (least-cost path)

Source pixel (5.0 C in 1995)

Source pixel (5.5 C in 1995)

Figure 6 | Workow for least-cost modelling. Illustration of how nal cost surfaces were generated for two temperature increments (5.0 and 5.5 C) and two examples of climate trajectories for the interval between 1995 and 2085. To account for changing temperature elds through time, we produced intermediate cost surfaces between the start and end years (a). Here we show a single intermediate temperature eld for simplicity; however our analysis utilized 31 intermediate temperature surfaces which were linearly interpolated between 1995 and 2085. Cost surfaces were generated based on equation 1 for the 5.0 C isotherm (b) and 5.5 C isotherm (c) (see Methods section). For each pixel, the minimum cost over all time periods is retained for the nal cost surface (d), which is then used to generate climate trajectories using a least-cost algorithm (e). Note that the example climate trajectories have the lowest accumulated cost compared to any other potential trajectories. This approach accounts for a changing climate by allowing trajectories to follow the climate as it warms. Cost surfaces were generated in this fashion for each 0.1 C increment (Table 1).

6 NATURE COMMUNICATIONS | 7:12349 | DOI: 10.1038/ncomms12349 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12349 ARTICLE

Table 1 | Examples of how climate analogues and costs were dened in the study.

Temperature of source pixel Cost 1 (climate analogue) Cost 2 Cost 3 Cost 4

5.0 C 5.0 C0.25 5.0 C(0.250.75) 5.0 C(0.751.25) 5.0 C(1.251.75)

4.755.25 C 4.254.75 C 3.754.25 C 3.253.75 C

5.255.75 C 5.756.25 C 6.256.75 C5.1 C 5.1 C0.25 5.1 C(0.250.75) 5.1 C(0.751.25) 5.1 C(1.251.75)

4.855.35 C 4.354.85 C 3.854.35 C 3.353.85 C

5.355.85 C 5.856.35 C 6.356.85 C5.2 C 5.2 C0.25 5.2 C(0.250.75) 5.2 C(0.751.25) 5.2 C(1.251.75)

4.955.45 C 4.454.95 C 3.954.45 C 3.453.95 C

5.455.95 C 5.956.45 C 6.456.95 C

Examples show increasing temperature dissimilarities for source pixels with 5.0, 5.1 and 5.2 C.

the resolution of the climate data from 2 to 20 km, (2) the bin width for dening a climate analogue and (3) the cost penalty associated with dissimilar climates (equation 1) for a subset of N. America. Specically, we used the state of Montana, USA as a subset as it has both mountainous and at terrain. For each permutation, we report the geometric mean, and the 5th, 50th and 95th percentiles for velocityMED, velocityED, the ratio between velocityMED and velocityED, and MCE (Supplementary Table 1, Supplementary Fig. 1, Supplementary Fig. 2).

Code availability. R code for implementing the methods is available at: http://adaptwest.databasin.org/pages/adaptwest-velocitymed

Web End =http://adaptwest.databasin.org/pages/adaptwest-velocitymed

Data availability. Outputs from the analysis are available at: http://adaptwest.databasin.org/pages/adaptwest-velocitymed

Web End =http://adaptwest.- http://adaptwest.databasin.org/pages/adaptwest-velocitymed

Web End =databasin.org/pages/adaptwest-velocitymed

References

1. IPCC. Climate Change 2001:The scientic basis. Contribution of working group I to the third assessment report of the intergovermental panel on climate change. (eds Houghton et al.) 881 (Cambridge Univeristy Press, 2001).

2. Garcia, R. A., Cabeza, M., Rahbek, C. & Araujo, M. B. Multiple dimensions of climate change and their implications for biodiversity. Science 344, 12475791247579 (2014).

3. Loarie, S. R. et al. The velocity of climate change. Nature 462, 10521055 (2009).

4. Hamann, A., Roberts, D. R., Barber, Q. E., Carroll, C. & Nielsen, S. E. Velocity of climate change algorithms for guiding conservation and management. Global Change Biol. 21, 9971004 (2015).

5. Saltr, F. et al. Climate change not to blame for late Quaternary megafauna extinctions in Australia. Nat. Commun. 7, 10511 (2016).

6. Pinsky, M. L., Worm, B., Fogarty, M. J., Sarmiento, J. L. & Levin, S. A. Marine taxa track local climate velocities. Science 341, 12391242 (2013).

7. Sandel, B. et al. The inuence of late quaternary climate-change velocity on species endemism. Science 334, 660664 (2011).

8. Ordonez, A. & Williams, J. W. Climatic and biotic velocities for woody taxa distributions over the last 16000 years in eastern North America. Ecol. Lett. 16, 773781 (2013).

9. Roberts, D. R. & Hamann, A. Climate refugia and migration requirements in complex landscapes. Ecography. doi: 10.1111/ecog.01998 (26th February 2016).

10. Dobrowski, S. Z. et al. The climate velocity of the contiguous United States during the 20th century. Global Change Biol. 19, 241251 (2013).

11. Carroll, C., Lawler, J. J., Roberts, D. R. & Hamann, A. Biotic and climatic velocity identify contrasting areas of vulnerability to climate change. PLoS ONE 10, e0140486 (2015).

12. Burrows, M. T. et al. The pace of shifting climate in marine and terrestrial ecosystems. Science 334, 652655 (2011).

13. Burrows, M. T. et al. Geographical limits to species-range shifts are suggested by climate velocity. Nature 507, 492495 (2014).

14. Ackerly, D. D. et al. The geography of climate change: implications for conservation biogeography. Divers. Distrib. 16, 476487 (2010).

15. Lenoir, J. et al. Local temperatures inferred from plant communities suggest strong spatial buffering of climate warming across Northern Europe. Global Change Biol. 19, 14701481 (2013).

16. Dobrowski, S. Z. A climatic basis for microrefugia: the inuence of terrain on climate. Global Change Biol. 17, 10221035 (2011).

17. Ohlemller, R., Gritti, E. S., Sykes, M. T. & Thomas, C. D. Towards European climate risk surfaces: the extent and distribution of analogous and non-analogous climates 19312100. Glob. Ecol. Biogeogr. 15, 395405 (2006).

18. Etherington, T. R. & Penelope Holland, E. Least-cost path length versus accumulated-cost as connectivity measures. Landsc. Ecol. 28, 12231229 (2013).

one time period (for example, 1995), trajectories would reect the minimum exposure in that time period only and would not reect the optimum path based on a changing climate. The second reason is one of accounting: MCE would always be 40 if we did not base the nal cost surfaces on the interpolated temperature grids. For example, had we based cost on the mean annual temperature in 1995, the cost associated with the destination location, which by denition is a climate analogue and should have a cost equal to one (equation 1), would be high because the cost would reect 1995 climate and not 2085 climate (compare top and bottom panels in Fig. 6b,c). Consequently, our approach accounts for a changing climate by allowing trajectories to follow the climate as it warms and acknowledges that paths between source and destination locations with no intervening topographic variation should not traverse dissimilar climates and will have an MCE 0

(for example, simple northward migration in at terrain).

For each source pixel in each temperature increment, we used the costDistance function of the gdistance package40 in the R statistical environment to calculate the accumulated cost to all pixels with the corresponding future climate analogue. We then identied the destination pixel with the least-accumulated cost and used the shortestPath function of the gdistance package to delineate the trajectory from source to destination pixel (Fig. 6e). By denition, trajectories delineated with these methods minimize exposure to dissimilar climates (as dened by our cost surfaces) and the length (km) of these trajectories is termed the MED. We subsequently calculated velocityMED (km year 1) as:

velocityMED

MEDtime 2 where time is the elapsed time between current and future time periods (time 90

years in this study). MCE was calculated as:

MCE

Pds costi li

MED

p 3

The term Pd

s costi li

is the least-accumulated cost, a standard output of least-

cost methods where costi is the cost assigned by equation 1 and is summed from s, the source pixel to d, the destination pixel, li is the length (km) of the trajectory though pixel i (this acknowledges that a diagonal trajectory through pixel i is longer than the horizontal or vertical equivalent; Fig. 6d), MED is the minimum exposure distance (which we subtract from total accumulated cost because it implicitly includes path distance; this is a result of assigning a value of 1 to all pixels with the same climate analogue; equation 1), and p is a climate dissimilarity penalty (as described in equation 1). In summary, MCE is a cumulative tally of climate dissimilarities that were encountered along each path. For example, MCE 2.5 C

indicates that the path, regardless of its total length, traversed 5 km (that is, one pixel; the resolution of our climate data) that was 0.5 C from the pixel of interest; MCE 5.0 C indicates that the path, regardless of its length, traversed

10 km (that is, two pixels) that was 0.5 C from the pixel of interest OR traversed 5 km (one pixel) that was 1.0 C from the pixel of interest. For comparative purposes, we also identied the ED nearest neighbour and ED-based velocity (cf., ref. 4) using the same climate analogue and binning approach described above.

We calculated velocityMED, velocityED, the ratio of velocityMED to velocityED and MCE for North America and produced corresponding maps. We excluded pixels from these maps (and all statistical analysis) where the source or destination pixel was located on an island; MCE values are not interpretable for these areas because we assigned large costs to water to ensure that trajectories circumvented water when possible (Supplementary Fig. 3).

To help visualize spatial patterns in exposure, we classied pixels within North America using a simple four-category scheme based on the bivariate distribution of velocityMED and MCE. Given that 71% of the MCE values were 0, we chose

MCE40 to delineate high values and MCE 0 to delineate low values. For

velocityMED, we used the median as the boundary between high and low values.

Sensitivity analysis. To evaluate the sensitivity of our results to input parameters and resolution of climate data, we conducted a sensitivity analyses by altering (1)

NATURE COMMUNICATIONS | 7:12349 | DOI: 10.1038/ncomms12349 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12349

19. Parks, S. A., Mckelvey, K. S. & Schwartz, M. K. Effects of weighting schemes on the identication of wildlife corridors generated with least cost methods. Conserv. Biol. 27, 145154 (2013).

20. Nunez, T. A. et al. Connectivity planning to address climate change. Conserv. Biol. 27, 407416 (2013).

21. Corlett, R. T. & Westcott, D. A. Will plant movements keep up with climate change? Trends Ecol. Evol. 28, 482488 (2013).

22. Huey, R. B. et al. Predicting organismal vulnerability to climate warming: roles of behaviour, physiology and adaptation. Philos. Trans. R. Soc. Lond. B Biol. Sci. 367, 16651679 (2012).

23. Schloss, C. A., Nunez, T. A. & Lawler, J. J. Dispersal will limit ability of mammals to track climate change in the Western Hemisphere. Proc. Natl Acad. Sci. USA 109, 86068611 (2012).

24. Tewksbury, J. J., Huey, R. B. & Deutsch, C. A. Putting the heat on tropical animals. Science 320, 12961297 (2008).

25. Deutsch, C. A. et al. Impacts of climate warming on terrestrial ectotherms across latitude. Proc. Natl Acad. Sci. USA 105, 66686672 (2008).

26. Grubb, P. J. The maintenance of species-richness in plant communities: the importance of the regeneration niche. Biol. Rev. 52, 107145 (1977).

27. Bell, D. M., Bradford, J. B. & Lauenroth, W. K. Early indicators of change: divergent climate envelopes between tree life stages imply range shifts in the western United States. Glob. Ecol. Biogeogr. 23, 168180 (2014).

28. Dobrowski, S. Z. et al. Forest structure and species traits mediate projected recruitment declines in western US tree species. Glob. Ecol. Biogeogr. 24, 917927 (2015).

29. Lawler, J. J., Ruesch, A. S., Olden, J. D. & Mcrae, B. H. Projected climate-driven faunal movement routes. Ecol. Lett. 16, 10141022 (2013).

30. Serra-Diaz, J. M. et al. Bioclimatic velocity: the pace of species exposure to climate change. Divers. Distrib. 20, 169180 (2014).

31. Korner, C. & Ohsawa, M. Ecosystem and Human Well-Being I (Island Press, 2005).

32. Scherrer, D. & Korner, C. Infra-red thermometry of alpine landscapes challenges climatic warming projections. Global Change Biol. 16, 26022613 (2010).

33. Anderson, M. G. & Ferree, C. E. Conserving the stage: climate change and the geophysical underpinnings of species diversity. PLoS ONE 5, e11554 (2010).

34. Lawler, J. J. et al. The theory behind, and the challenges of, conserving natures stage in a time of rapid change. Conserv. Biol. 29, 618629 (2015).

35. Keppel, G. et al. Refugia: Identifying and understanding safe havens for biodiversity under climate change. Glob. Ecol. Biogeogr. 21, 393404 (2012).

36. Hannah, L. et al. Fine-grain modeling of species response to climate change: Holdouts, stepping-stones, and microrefugia. Trends Ecol. Evol. 29, 390397 (2014).

37. Keppel, G. et al. The capacity of refugia for conservation planning under climate change. Front. Ecol. Environ. 13, 106112 (2015).

38. Hamann, A., Wang, T., Spittlehouse, D. L. & Murdock, T. Q. A comprehensive, high-resolution database of historical and projected climate surfaces for western North America. Bull. Am. Meteorol. Soc. 94, 13071309 (2013).

39. Dijkstra, E. W. A note on two problems in connexion with graphs. Numer. Math. 1, 269271 (1959).

40. van Etten, J. gdistance: Distances and Routes on Geographical Grids. R package version 1.1-9. (2015).

Acknowledgements

S.Z.D. was supported by the National Science Foundation (DEB 1145985; BCS 1461576) and USFS Rocky Mountain Research Station (Agreement number: 15-JV-11221639-119). We thank Brady Allred for providing computing resources, and John Abatzoglou and three anonymous reviewers for constructive comments on previous drafts of this manuscript.

Author contributions

S.Z.D. and S.A.P. contributed equally to the design, analysis, and writing of this paper.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Dobrowski, S. Z. et al. Climate change velocity underestimates climate change exposure in mountainous regions. Nat. Commun. 7:12349 doi: 10.1038/ncomms12349 (2016).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

r The Author(s) 2016

8 NATURE COMMUNICATIONS | 7:12349 | DOI: 10.1038/ncomms12349 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications