Skip to main content

Rattlesnake migrations and the implications of thermal landscapes



The importance of thermal resources to terrestrial ectotherms has been well documented but less often considered in larger-scale analyses of habitat use and selection, such as those routinely conducted using standard habitat features such as vegetation and physical structure. Selection of habitat based on thermal attributes may be of particular importance for ectothermic species, especially in colder climates. In Canada, Western Rattlesnakes (Crotalus oreganus) reach their northern limits, with limited time to conduct annual migratory movements between hibernacula and summer habitat. We radio-tracked 35 male snakes departing from 10 different hibernacula. We examined coarse-scale differences in migratory movements across the region, and then compared the route of each snake with thermal landscapes and ruggedness GIS maps generated for different periods of the animals’ active season.


We observed dichotomous habitat use (grasslands versus upland forests) throughout most of the species’ northern range, reflected in different migratory movements of male snakes emanating from different hibernacula. Snakes utilizing higher-elevation forests moved further during the course of their annual migrations, and these snakes were more likely to use warmer areas of the landscape.


In addition to thermal benefits, advantages gained from selective migratory patterns may include prey availability and outbreeding. Testing these alternative hypotheses was beyond the scope of this study, and to collect the data to do so will require overcoming certain challenges. Still, insight into migratory differences between rattlesnake populations and the causal mechanism(s) of migrations will improve our ability to assess the implications of landscape change, management, and efficacy of conservation planning. Our findings suggest that such assessments may need to be tailored to individual dens and the migration strategies of their inhabitants. Additionally, local and landscape-scale migration patterns, as detected in this study, will have repercussions for snakes under climate-induced shifts in ecosystem boundaries and thermal regimes.


The conservation of migratory animals may be particularly challenging, as movements may cross multiple habitats and jurisdictions, including borders between protected and unprotected areas [58]. A starting point is determining habitat use patterns throughout the entire life-history of these species, followed by an understanding of why these patterns exist, and why particular migratory paths are followed. Such information is critical to crafting effective conservation and management strategies.

Migrations occur when animals move explicitly to take advantage of resources distributed through the environment [20, 23, 35, 66]. Migration is differentiated from dispersal movements by way of bidirectional movement, movement of greater duration than ranging, and movement along relatively straight-line paths [20]. Such responses may be apparent when they involve clear, contrasting shifts in habitat use, such as birds migrating between temperate and tropical regions. However, migrations also occur on a finer landscape scale, where habitat factors exert more subtle effects. Food, water or mate availability may drive such movements [2], but for animals inhabiting cooler regions, heat potentially may be another resource influencing migration. Ectothermic animals occupying relatively harsh or variable environments may be particularly responsive to the thermal properties of landscapes, given the lack of a metabolic buffer between the thermal environment and the physiological function of the individual. Although there is consensus that a prominent factor in habitat selection by ectotherms is temperature [7, 17], this relationship has been primarily studied on a fine-scale [11, 17, 33, 64, 71]. Snakes, particularly northern species, provide rare examples of terrestrial ectotherms conducting annual migrations; these migrations can include considerable movements relative to the mobility of the animal and the brevity of the active season (e.g., [36, 47, 73]). However, the potential linkages between these movements and thermal landscape patterns has not been well-studied [72]. By testing for this mechanism, we gain a better understanding of why particular movements (including seasonal migrations) occur. This, in turn, will facilitate more refined habitat mapping and consideration of how other land uses could potentially influence migrations.

Seasonal, relatively lengthy migrations (> 3–5 kms) have been documented for the Western Rattlesnake (Crotalus oreganus) in areas close to the species’ northern limits [10, 29]. Presumably, as in other northern snakes, these movements are partly due to the geographic separation of suitable overwintering hibernacula (‘dens’) and summer habitat conducive to foraging and mating [31]. However, exactly why these animals conduct extensive movements is unclear, particularly when the animals are travelling through areas that contain prey. Further, the extent and directionality of the migratory movements appear highly den specific, with adult individuals taking similar pathways from year to year and, in at least one study, the animals may travel into higher-elevation forests away from the normally associated grassland habitat [29]. While den sites may be afforded protection through land designations (parks, protected areas, etc.) or stewardship by private land owners, migration routes can result in snakes travelling across different jurisdictional boundaries and into areas where protection is not afforded [85]. At present, we lack a clear understanding of why these animals use lengthy and specific migration routes during the summer active season.

Ectotherms, like snakes, employ strategies to maintain their body temperatures at optimal levels. These strategies include physiological mechanisms for temperature control and/or behavioural tactics. Behavioural thermoregulation carries ecological benefits and costs [77]; for example, shifting positions to maintain optimal body temperatures [28, 37, 84] will be worthwhile for an animal only the availability of other resources does not plummet [38] such as food and access to mates [7, 22]. Thermoregulation thus may be considered a form of ‘resource utilization’ [39] and variation of heat as a resource across different habitats likely will have repercussions [76].

Hypothetically, variations in movements and habitat use exhibited by different den populations of rattlesnakes and other northern snakes may be caused by temperature selection at the landscape level. Thermal characteristics on the ground will be influenced by the incident solar radiation that, in turn, is influenced by substrate ruggedness, day length, topographical shadows, and solar azimuth at the given latitude [40]. On an annual timescale, seasonal and daily variation in these dynamics likely do not cause major changes in vegetation communities typically used to delineate habitats; however, they may provide significantly different resources (from a snake’s point of view) at different points in the season and across a heterogeneous landscape. Over the course of the active season, the average thermal properties of the landscape may alter the migratory and large-scale movement patterns of snakes as they leave and return to hibernacula. If snake migration can be linked to thermal attributes of landscapes, it would provide powerful new insight into how the ecology of these animals is influenced on large scales. Being able to predict such movements also would be an important conservation tool.

The goals of this study were (i) to use telemetry to monitor the seasonal migration of Western Rattlesnakes from a relatively large sample of hibernacula across the northern range of the species, (ii) demonstrate more clearly the existence of different migratory phenotypes and associated habitat types, and (iii) examine the role that thermal attributes of the landscape play in large-scale habitat selection. Our working hypothesis was that snakes select warmer areas within the available habitat (‘landscape’) to gain thermoregulatory benefits over the course of the northern summer. We used solar insolation as a proxy for temperature; this metric has been used to predict hibernacula locations [32] but it has not been applied to investigate snake summer habitat use. Thus, we examined whether thermal properties, like other long-term resource distributions, are correlated with habitat use by animals. We conducted our work at the northern limit of rattlesnakes (western Canada), where harsher and more variable temperatures intuitively should cause the thermal attributes of landscapes to exert a relatively strong influence on seasonal migratory movements.

Study species

The Western Rattlesnake (Crotalus oreganus) is the northernmost viper (F. Viperidae) in the Western Hemisphere, reaching lengths just over 1.2 m in snout-vent length [3]. The species’ range extends across the USA:Canada border into the southern region of British Columbia (BC, Canada), where the snake feeds primarily on rodents and birds [56]. In BC, the species occupies a limited and narrow range, being restricted primarily to the dry semi-arid valleys of the Okanagan and Thompson rivers (see map in [54]). Here, the snakes den communally, displaying high fidelity during their life to specific hibernacula typically located in south-facing rocky crevices, fissures and talus slopes in low-elevation (500–625 m above sea level [MSL]) grasslands [5, 51, 52]. After the snakes emerge from hibernacula in the spring (April, May), males disperse for the summer [53] presumably to feed and locate mates, and then return to the den in autumn (September, October). The movement ecology of females is far less studied, but reproduction is known to be biennial or longer, and there is some evidence that nongravid females also disperse from the den site for summer [52]. Mating takes place away from dens in late summer (July, August) with fertilization delayed until the ensuing spring. As embryos develop, the gravid females remain close to the den site and occupy ‘rookeries’ through gestation and parturition (Larsen pers. obsv.). The lower-valley locations of most hibernacula have led to these northern snakes historically being associated with grassland habitats during the entire active season; however, telemetry work by Gomez et al. [29] revealed that the snakes in one denning population migrated into higher-altitude forests, moving up to 3500 m from the hibernacula over the course of the summer, while snakes from a neighbouring den population remained primarily in the grasslands. Other studies have referenced occasional forest habitat use [10, 13, 14, 29, 50], suggesting the stereotyped habitat association of the species needs further investigation. In particular, the extent of migratory movements by this animal, and why they vary between different hibernacula populations, is a significant gap in our understanding of the species’ ecology in the far north.


Selection of sites

The study was conducted in 2010 and 2011 in the Thompson-Nicola (50.8°N, 120.6°W) and Okanagan-Similkameen (49.3°N, 119.6°W) regions of British Columbia, Canada. This is a semi-arid region containing the northernmost populations of this species, with an elevational range from approximately 300 to 1400 MSL. Details on the ecosystem characteristics are provided in Fig. 1 and in Harvey [34]. Study hibernacula were selected to ensure a diversity of putative thermal landscape conditions were present between hibernacula. This was done using preliminary thermal maps of the area created with GIS (details below). Considerations also were made for logistics, land ownership and access. Further, hibernacula were only considered for the study if the estimated population size was > 12 snakes [9]; this was done to reduce the likelihood that the study would impact the viability of the population. In total, ten different hibernacula were chosen for this study – six in the Thompson-Nicola region (six in 2010, with one repeat in 2011) and four in the Okanagan-Similkameen region (in 2011).

Fig. 1

Typical and atypical habitat associations for the Western Rattlesnake in British Columbia. Open grassland habitats generally occur at elevations of 300 to 800 m in the Bunchgrass and Ponderosa Pine biogeoclimatic zones. Forested habitats generally occur at elevations of 500 to 1200 m in the Interior Douglas-fir biogeoclimatic zone [57]. Hibernacula typically occur at elevations of 500 to 625 m on south-facing, rocky slopes. Photos by lead author

In total 17 snakes from six hibernacula were selected for inclusion in the study in 2010, and 18 snakes from five hibernacula were selected in 2011. We targeted at least three snakes from each hibernaculum in order to achieve replication. However, at two hibernacula only one and two snakes respectively were deemed suitable for the study (based on length and weight). At the remaining nine sites, either three or four snakes were selected as study animals.

Animal capture, processing, selection and surgery

We visited targeted hibernacula repeatedly (i.e., at least twice) during the spring emergence period (April 15 – May 7). Rattlesnakes were captured using tongs and weighed in a canvas bag. Individuals suitable for telemetry (i.e., heavier than 400 g) were ushered into a plexiglass tube for sex determination via hemipenal probing [75]. Because females vary in their migratory behaviour according to their reproductive status (see above - [51]), we selected only male snakes for radio-telemetry. All animal selected for telemetry were transported by vehicle to a veterinary clinic in a towel-lined, aerated rubber container.


Each study animal was surgically implanted with an SB-2 radio-transmitter (Holohil Systems Inc., Ontario, Canada), weighing 3.8–5.2 g. No implant package exceeded 2.7% of the snake’s weight, and transmitter lifespans ranged from 5 to 10 months. Surgeries were carried out by veterinarians following the protocols described by Reinert and Cundall [67] with modifications by Reinert [68]. Implanted snakes were held for approximately 24–48 h post-operation to permit recovery from sedation and allow rehydration. Each animal then was released at its precise point of capture. We tracked and located telemetered snakes every three to 7 days between egress (April/May) and ingress (September/October) using an R-1000 telemetry receiver (Comm-Spec) and RA-159 handheld Yagi directional antenna (Advanced Telemetry Systems). When each snake was located, we recorded date, time, UTM coordinates using a handheld GPS unit (Garmin, GPS 76Cx, accuracy ±3 m), canopy closure using a spherical crown densitometer (Forestry Suppliers, Convex model A), and habitat type (see below). Transmitters were removed from re-captured snakes either when they returned to their hibernacula in autumn of the same year, or as they emerged from hibernation the following spring.

Mapping and analysis

Snake location data were filtered to include only those that constituted independent movements, defined as more than 10 m from the previous location [29]. Each location was assigned to one of two habitat types to enable comparison: those with < 10% canopy closure (bunchgrass and open-canopy Ponderosa Pine) were designated as “Open” habitats, while locations with > 10% canopy closure (Interior Douglas-fir forests) were classified as “Forest” habitat. In the Kamloops region, this threshold for “Forest” was deemed appropriate given the arid nature of the region and the typical canopy closure in the Ponderosa Pine zone. Additionally, each snake was later assigned a category based on the type of habitat reached at the endpoint of its migration, termed the ‘destination habitat’. This classification of habitats allowed comparison between snakes using typical (open grassland) habitats and atypical (forested) habitats (Fig. 1).

The migration path for each individual snake was created by connecting telemetry data points, assuming direct movements between locations. The maximum straight-line distance travelled by each snake from its hibernaculum was calculated using the furthest detected location (apogee). Also, we focused our movement analyses largely on the migration path taken by each snake out to the apogee, terming this the Outward migration path [rattlesnakes tend to simply reverse their Outward migration paths during the ‘homeward’ migration phase [29] when movements appear largely driven by fidelity to the den site]. The Outgoing migration was further divided into an Initial stage and a Late stage containing the movements of each snake in May/June and July/August, respectively.

Our first step (detailed in Fig. 2) was to quantify the amount of solar radiation each migratory path would receive in an average year. To do this, we created thermal models of the landscapes surrounding the hibernacula from whence the telemetered snakes originated. These landscape simulations were composed of raster images built on the predicted incident solar radiation for a pixel (cell size) of 25 m2. To quantify the solar radiation each pixel would receive, we used the Area Solar Radiation tool in the Spatial Analyst extension in ArcGIS 9.3 [25]. This algorithm uses slope, aspect, elevation, day length, latitude and solar azimuth to calculate the expected incident solar radiation at a given point on the landscape. The simulation was based on a 25 m digital elevation model (DEM), with the thermal values being expressed in average daily Watt-hours per square meter (Wh/m2) (Fig. 3). It is important to appreciate that these maps did not depict the actual solar radiation a snake would experience during any given year; rather, the modelling process quantified the level of solar radiation (i.e. temperature) each pixel would be expected to receive on average across years, based on general topographic features of the landscape [8, 21, 27].

Fig. 2

Step-by-step process by which assessment of the thermal properties of migration routes was conducted for rattlesnakes in British Columbia, Canada. The entire process outlined was separately carried out on the entire Outward migration pathway, as well as separately on the Initial and Late stages of the same migration pathway. Definitions and additional details are provided in the Methods section of the text

Fig. 3

Examples of Western Rattlesnake migrations from the hibernaculum to destination habitat. In the top illustration, the study hibernaculum is represented by a star (), and different snake movements are represented by different-coloured triangles (Δ) and connecting lines. In the top image, the thermal landscape is represented using modelled average incident solar ration as a proxy for temperature. Telemetered snakes in this example all used warmer and forested areas of the thermal landscape (shown as darker green in the bottom image)

Our second step (Fig. 2) was to produce a ‘thermal score’ for the migratory routes demonstrated by the snakes. To do this, we calculated the total of the solar incidence values assigned to those pixels traversed by each snake’s Outward migration path. We also calculated these values for the Initial and Late portions of the migratory pathways.

To examine the significance of the pathway thermal values, we used a relatively simple analytical approach, given that more elaborate and sophisticated step-selection functions (SSF - [79]) target situations where massive amounts of sequential spatial data are collected on animals, primarily through GPS (e.g. [62]), other remote tracking methods, or constant tracking [80] – all methods we could not employ in this study. These SSF models become more difficult to parameterize and interpret when relocation frequencies are low and vary considerably between animal subjects, such as the situations we faced in our study (see Results). Further, to our knowledge, this type of study has not appeared in the literature, and more complex analyses could make the results more difficult to understand while likely leading to similar conclusions [59].

Our analytical approach started with the creation of random-walk migrations for each of the snakes we tracked (see Step 3, Fig. 2). This allowed comparisons involving our empirical/observed thermal values to those theoretically ‘available’ on the landscape, thus representing a form of resource use [44]. We used Hawth’s tools [6] to generate 100 random walks for each snake, starting at the animal’s hibernaculum. These random walks were parameterized using the actual Outward movements of the snake: i.e. the overall lengths of the random walks were restricted in length to the greatest distance measured from the hibernaculum to the snake (the apogee), and an average turn angle equal to that demonstrated by the animal in the field. For each random walk, we calculated a thermal value for the simulated pathway in the same manner as that previously done for the empirical (observed) movement pathway (see above). Repeating this procedure for the 100 simulated random walks by the snake produced a distribution of thermal values by which to compare the empirical value. From there, we simply determined the percentile score of the observed thermal value within the distributions of values obtained from the random walk simulations (Step 4, Fig. 2). We repeated these calculations on the Initial and Late stages of the migrations. The end result was percentile values reflecting the Outward migrations of the sample of telemetered snakes (n = 30, see Results), with similar datasets derived from the Initial and Late stages of the migrations.

Our statistical analysis was three-fold: We first tested whether the observed thermal values for the snakes were relatively high or low given the landscapes they were migrating through. We simply noted whether they fell above or below the 50th percentile scores in the corresponding distribution of simulated values and applied a χ2 test Ho: 50% of observed observations will fall above the 50th percentile of the 100 simulations (Step 5 – Fig. 2). We then tested if the percentile scores assigned to the snakes differed between animal using forest versus open destination habitats, using a t-test ([41, 88]- Step 6 – Fig. 2). Lastly, Z-tests [87] were used to determine the probability of selecting the observed thermal value of each snake’s migration from the distribution of values drawn from the corresponding simulated walks (Step 7 – Fig. 2). All of these analyses were repeated for each of the three divisions of snake migration (Outgoing, Initial, and Late). Normality of data were confirmed in all cases prior to analyses.

We also conducted an analysis on the role of ‘ruggedness’ in influencing the snake movements. Ruggedness as used in GIS environments in based on differences in elevation between a focal pixel and the eight pixels surrounding it [69]. Using Relative Topographic position and the raster calculator tool, each of these 8 differences is squared (to create all positive values) and a mean value is calculated. The square root of this value then produces the ruggedness index value [42] representing the average elevation change between the focal pixel and the neighbouring area. A ruggedness index value was calculated for each pixel on the landscape using a moving window analysis [42, 69]. The mean value of these pixels for each of the snakes’ outgoing migration paths and for an area with a radius of 4000 m surrounding each hibernaculum was calculated. Ruggedness could theoretically influence snake movements in two manners: (i) ruggedness could impact the mobility and capabilities of a snake to traverse the area, and/or (ii), a high ruggedness value may represent a diversity of microclimates existing on the site, allowing a snake to locate warm areas even on sites where the superficial thermal index is low.

Lastly, to determine the relationship between habitat ruggedness, thermal landscape characteristics, and habitat use, we examined the data in two ways. We first tested for a relationship between the two categories of snakes based on their destination habitat (forest versus open) using a t-test. We then used linear regression to explore the relationship between ruggedness index and the percentile scores calculated for the thermal values of each snake’s migratory pathway.

Statistical considerations

All statistical analyses were performed in the program R version 2.12.1 [65]. Data were tested for normality by examination of histograms and using the Shapiro-Wilk test or the Kolmogorov-Smirnov test [87]. Homogeneity of variances between groups was tested using the Fligner-Killeen test [15, 16]. Percentage data were transformed for normality using an arcsine transformation for analysis. A significance value of α = 0.05 was used to guide the interpretation of the results. Means are reported ±1 standard deviation, unless otherwise stated.


In total, we used 35 snakes for the radio-telemetry study. Thirty of 35 telemetered snakes were tracked through their entire annual migration, resulting in 13 to 27 location data points per individual. Five snakes did not have their entire migration route completely documented either due to mortality or signal loss. These snakes were not considered in the analysis due to incomplete data sets.

In total, the migration data from 30 snakes were included in the analysis. As expected, the snakes travelled away from the hibernaculum to summer habitats (Outgoing migration), reaching their apogee on an average date of August 8. The earliest apogee occurred on June 23 and the latest on September 21. Most snakes returned to the hibernaculum (Homeward migration) along a very similar path as their Outgoing migration. The mean maximum straight-line distance measured for Outgoing migration (from the hibernaculum to the apogee) of the telemetered snakes was 1847.8 m ± 930.0 m (n = 30, range = 373.0 m to 3985.7 m). The migration metrics are presented in Table 1. There were no significant differences observed between straight-line migration distances between the Thompson-Nicola and Okanagan-Similkameen regions [migration (t26 = 1.42, P = 0.170).

Table 1 Distances and trends of movements by Western Rattlesnakes radio-tracked from ten hibernacula in British Colombia (2010 and 2011). The movements of snakes originating at the same den were considered ‘Directional’ if all snakes tracked from the hibernaculum travelled trajectories within a 40° range of each other and ‘Random’ if the trajectories ranged over more the 40°. The ruggedness index was calculated for an area surrounding each hibernaculum using Relative Topographic position. Means are reported ±1 standard deviation, unless otherwise stated

All study hibernacula (and thus the starting points of all monitored migrations) were located in open habitats. From these sites, the mean distance to forest habitat was 656 m (± 958 m). Fifteen of the 30 snakes, henceforth termed ‘Forest snakes’, travelled to and used forests as a Destination habitat (i.e., in the latter part of their outgoing migrations during July and August), while the use of open habitats through the entire migration was observed in the other 15 study snakes (henceforth termed ‘Open-Habitat snakes’ - see Table 1). Maximum straight-line distances reached from the hibernaculum were significantly longer for Forest snakes (2359 m ± 837.0 m) than Open-Habitat snakes (1337 m ± 729.0 m; t27 = 3.57, P = 0.001) for the entire season.

Due to the limited sample size of snakes at each study site (<four from the majority of the study hibernacula), statistical analysis of migration directionality by the snakes from each hibernaculum was not possible. We, therefore, used a less-rigorous approach by classifying migrations from a particular hibernaculum as directional when the telemetered snakes leaving that hibernaculum displayed mean migration bearings within 40° of one another (Table 1). Within-hibernaculum groups of snakes whose migration bearings were more than 40° from one another were considered to have a random distribution. This distinction was made based on a natural break in the data and qualitative judgement of snakes’ travel directions.

The observed thermal metrics of the Forest snakes tended to occur in the upper half of the distribution of simulated migration paths (i.e., relatively warmer pathways) more often than Open-Habitat snakes, in all three categories of the migration (Outgoing, Initial and Late stages; see Fig. 3, Table 2A).

Table 2 Comparisons of observed thermal metrics to the thermal metrics derived from 100 random-walk migration path simulations for Western Rattlesnakes in British Columbia. Using thermal landscape maps, both observed thermal metrics and those of the simulated migrations were derived from the average incident solar radiation along the migration path for each migration category (Outgoing, Initial and Late stages). Forest snakes utilized forested habitats as the destination for their migration, while Open-Habitat snakes remained in sparsely-treed or open grasslands throughout the active season. Means are reported ±1 standard deviation, unless otherwise stated

In all three categories of the migration (Outgoing, Initial, Late) the observed thermal values derived from the migratory pathways of the Forest snakes were significantly higher than the same values for the Open-Habitat snakes, as compared to their respective simulated movements (Table 2 B, Fig. 4). These differences were most noticeable during the Late stage of migration, when the thermal values of the observed (observed) pathways for Forest snakes averaged near the 80th percentile, compared to a 50th percentile average for the snakes that remained in open habitats. In fact, thermal values for the Open-Habitat snake migrations averaged close to 50th percentile scores in all three migration categories (Table 2 B). The average probability values (from z-test scores) as determined for the Outgoing, Initial and Late stages of the migration were significantly lower for Forest snakes than Open-Habitat snake (Table 2 C).

Fig. 4

Average percentile scores of migration path values within values derived from 100 random-walk simulations, for each category of migration, compared to migration path ruggedness for telemetered Western Rattlesnakes in British Columbia. Trend lines (shown) were fit for the pooled group of snakes during the entire Outgoing migration (R2 = 0.13), the Initial migration (R2 = 0.09) and the Late stage of migration (R2 = 0.13). Forest snakes are represented by solid markers; Open-Habitat snakes are represented by open markers

The ruggedness index values for migration paths were significantly higher for Forest snakes ( =77.7 ± 23.0) than for Open-Habitat snakes (x̄ =56.6 ± 13.2; t22 = 3.09, P < 0.01). A significant effect of ruggedness was found for the percentile score of the observed migration path values for the Late migration category (F1,28 = 4.31, P = 0.047; R2 = 0.13), but not for either the Outgoing migration (F1,28 = 4.05, P = 0.054; R2 = 0.13) or Initial migration stages (F1,28 = 2.97, P = 0.095; R2 = 0.09; Fig. 3). Rattlesnake hibernacula with higher average ruggedness values were more likely to produce snakes that migrated to forest habitat (t25 = 3.99, P = < 0.001).


Overall, the results of our study indicate that the annual migrations of these northern snakes are dictated at least in part by thermal attributes of landscapes, at a relatively large-scale, particularly for snakes that utilize forested habitat. This is intuitive for terrestrial ectotherms occurring at high latitudes, but what is more interesting is the fact the animals appear to travel relatively longer distances to access this habitat, and their movements take them out of the lower arid grassland valleys that are considered typical summer habitat. This pattern also was not universal: exactly half of the animals we followed undertook the longer migrations into the higher elevation forests, and slightly greater than half of the animals travelled into relatively warmer areas on the landscape. Thus, the thermal parameters we examined in this study do not fully explain the dichotomy of movements exhibited by snakes in this and a previous study [29], but they do shed light on the factors influencing migration patterns in northern herpetofauna.

Use of forest habitat by rattlesnakes is well-documented but in locations further south and east than our study location [63, 81, 83]. Similarly, long-distance movements also have been reported elsewhere [4, 24, 45] and within the same region of our study [29]. Long-distance migratory movements for northern snakes, in general, have often been hypothesized to reflect widely separated resources, such as hibernacula and summer foraging sites. The mean migration distances in our study were similar to those previously reported in this region; however, several of the maximum distances recorded in this study were longer than those previously reported [5, 13, 14, 29, 49, 52]. The prevalence of forest habitat use [first detected by Gomez et al. [29] and now well-demonstrated by our study] is somewhat unexpected, given the presumption that the animals at their northern limits should be strongly tied to the arid, open grassland habitat of the valley bottoms [52, 55].

Our analysis partially explains the patterns of summer habitat use by these animals. The variable terrain of British Columbia coincides with the northern limit of the species, providing energetic and thermoregulatory challenges for the animals and raising the benefits of using habitat with optimal thermal attributes. In particular, the thermal landscape properties of the outgoing migration as a whole and the late stage of migration differed the most from the simulated random walks. During this time, most of the snakes travelled along warmer pathways, and snakes heading to or occupying forest habitats tended towards warmer pathways. Snakes appeared to move through a variety of conditions to reach destination habitats with ideal thermal properties, as evidenced by fewer snakes travelling warmer pathways during the initial migration stage. The consequences of using cooler paths are not known; however, it is likely that microhabitat selection may compensate for temperature changes [11, 77, 86].

Rattlesnakes migrating through landscapes with higher ruggedness (more elevation variation in the terrain) were more likely to use forested habitats and had higher migration path percentiles during the late stage of migration. Although this relationship was significant, the amount of variation actually attributable to ruggedness was low. Still, ruggedness has been included as an important attribute in habitat selection for a variety of wildlife including caribou [60], bighorn sheep [74], badger [1] and grouse [12], and our data support the assertion that this metric should be considered an important influence in snake habitat use [26, 30], at least in northern areas with noticeable variation in topography.

Thermal patterns of the landscape (as we measured them) appeared to be influencing the migrations observed in our study; however, there may be additional factors dictating migration routes that work in combination or separate from thermal attributes of the landscape. Animals may migrate in search of resources such as prey, mates or other suitable habitat conditions. Several studies link altered spatial behaviours to prey availability [23, 82], whereas others indicate limited support for this effect [61, 78]. In Wyoming, the movement of male rattlesnakes has been attributed to mate-searching [24]. There is, however, insufficient knowledge to extrapolate these effects to other locations, such as our northern study site. While our results indicate that there may be landscape thermal influences on snake movements, investigation of other putative factors is warranted.

As with all ecological models, thermal landscape simulations are simplistic representations of complex systems. The thermal models developed in this study provided insight into the role that thermal attributes of the landscape play in rattlesnake habitat use during the active season; however, they may be constrained by a spatial database resolution of only 25 × 25 m pixels. As the thermal landscape used in the analysis is based on this resolution, any variation occurring below this scale is not captured. A higher resolution digital elevation model, perhaps 3 × 3 m, in concert with a ground-cover mapping layer such as LIDAR, could be used to examine landscape dynamics through a finer lens.

The results of this study demonstrate that the relationship between hibernaculum location, migration distance and direction, and summer habitat utilization for these animals is far more complex than initially suspected. This is particularly important given our understanding of northern rattlesnake ecology historically has been based largely on one study that collected detailed information on one population of snakes [52]. Clearly, widespread migratory differences exist both within and between hibernating populations of these animals. To further our understanding of the ecology of migratory animals, we need to realize di- and polytomous habitat use may require a reassessment of “typical” habitat [18, 70].

Climate-induced shifts in ecosystem boundaries and thermal regimes, such as encroachment of forests into grassland habitats and warming climate, may have significant effects on the habitat use by reptiles, including the subjects of this study [43, 48]. The effects on northern snakes are likely two-fold: First, changes in the vegetation will modify the structure and heterogeneity of the landscape, requiring behavioral modifications to small and large-scale thermoregulation [46]. On a landscape scale, expansion of one habitat type (e.g. the bunchgrass zone in this study) could mean longer migrations to reach other habitat (forested habitat as detected in this study). This will lead to increased energy expenditure on travel, perhaps reducing overall fitness. This potentially will have significant impacts on reptile species that use these habitats in temperate regions [19].


Although the thermal landscape seems to explain some variation with respect to forest migrants, grassland migrants selected pathways that seemed random with respect to the thermal landscape, at least at the resolution of our analysis. Still, this work provides new insight into the factors dictating seasonal snake movements in locations where extensive migrations occur. On a landscape scale, snakes use habitats that provide a thermal advantage through the short, northern summer. The role of thermal landscape attributes in colder environments, and how they affect migratory pathways of snakes (and other species) warrants consideration along with other resource values in assessing habitat. For northern snakes, developing a tool that precisely predicts summer movement patterns of snakes emanating from specific dens may present a lofty challenge, but any insight or mapping tool that can delineate probable movement pathways will greatly improve our ability to manage ecosystems. This may be particularly important where economic drivers of landscape use (e.g. ranching, forestry) coincide with high-value species habitat. Additionally, understanding local and landscape-scale movement patterns has obvious management applications for species likely to be affected by climate-induced shifts in ecosystem boundaries and thermal regimes.

Availability of data and materials

Digital elevation model and biogeoclimatic zone spatial data are available upon request from the BC Ministry of Forests, Land, Natural Resource Operations and Rural Development. Telemetry data are archived with the BC Conservation Data Centre as part of the Wildlife Species Inventory (WHI) managed by the BC Ministry of Environment.



Geographic Information System


Standard Deviation


Universal Transverse Mercator


  1. 1.

    Apps CD, Newhouse NJ, Kinley TA. Habitat associations of American badgers in southeastern British Columbia. Can J Zool. 2002;80:1228–39.

    Article  Google Scholar 

  2. 2.

    Ashton KG. Movements and mating behavior of adult male midget faded rattlesnakes, Crotalus oreganus concolor, in Wyoming. Copeia. 2003;2003:190–4.

    Article  Google Scholar 

  3. 3.

    Ashton KG, de Queiroz A. Molecular systematics of the Western rattlesnake, Crotalus viridis (Viperidae), with comments on the utility of the D-loop in phylogenetic studies of snakes. Mol Phylogenet Evol. 2001;21:176–89.

    CAS  PubMed  Article  Google Scholar 

  4. 4.

    Bauder JM, Akenson H, Peterson CR. Movement patterns of prairie rattlesnakes (Crotalus v. viridis) across a mountainous landscape in a designated wilderness area. J Herpetol. 2015;49:377–87.

    Article  Google Scholar 

  5. 5.

    Bertram N, Larsen KW, Surgenor J. Identification of critical habitats and conservation issues for the western rattlesnake and great basin gopher snake within the Thompson-Nicola region of British Columbia. Report prepared for the British Columbia Ministry of Water, Land and Air Protection and the Habitat Conservation Trust Fund of British Columbia, Kamloops, BC, Canada. 2001. Accessed 24 Jul 2018.

  6. 6.

    Beyer HL. 2004. Hawth's analysis tools for ArcGIS. 2004 http://wwwspatialecologycom/htools Accessed 5 May 2010.

    Google Scholar 

  7. 7.

    Blouin-Demers G, Weatherhead PJ. An experimental test of the link between foraging, habitat selection and thermoregulation in black rat snakes Elaphe obsoleta obsoleta. J Anim Ecol. 2001;70:1006–13.

    Article  Google Scholar 

  8. 8.

    Bode CA, Limm MP, Power ME, Finlay JC. Subcanopy Solar Radiation model: Predicting solar radiation across a heavily vegetated landscape using LiDAR and GIS solar radiation models. Remote Sens Environ. 2014;154:387–97.

    Article  Google Scholar 

  9. 9.

    British Columbia Conservation Data Centre. Endangered Species and Ecosystems – Confidential Sensitive Occurrences (digital file). Victoria; 2009.

  10. 10.

    Brown JR, Bishop CA, Brooks RJ. Effectiveness of short-distance translocation and its effects on western rattlesnake. J Wildl Manag. 2009;73:419–25.

    Article  Google Scholar 

  11. 11.

    Brown WS, Pyle DW, Greene KR, Friedlaender JB. Movements and temperature relationships of timber rattlesnakes (Crotalus horridus) in northeastern New York. J Herpetol. 1982;16:151–61.

    Article  Google Scholar 

  12. 12.

    Carpenter J, Aldridge C, Boyce MS. Sage-grouse habitat selection during winter in Alberta. J Wildl Manag. 2010;74:1806–14.

    Article  Google Scholar 

  13. 13.

    Charland MB, Nelson KJ, Gregory PT. Status of northern pacific rattlesnake in British Columbia. Wildlife Working Report No. WR-54. Victoria: Ministry of Environment, Lands and Parks, Wildlife Branch; 1993a. Accessed 24 Jul 2018.

    Google Scholar 

  14. 14.

    Charland MB, Nelson KJ, Gregory PT. Status of Northern Pacific Rattlesnake in British Columbia. Wildlife Working Report No. WR-54. Victoria: Ministry of Environment, Lands and Parks, Wildlife Branch; 1993b.

    Google Scholar 

  15. 15.

    Conover WJ, Johnson ME, Johnson MM. A comparative study of tests for homogeneity of variances, with applications to the outer continental shelf bidding data. Technometrics. 1981;23:351–61.

    Article  Google Scholar 

  16. 16.

    Crawley MJ. The R Book. West Sussex, London: Wiley; 2007.

    Google Scholar 

  17. 17.

    Diaz JA. Ecological correlates of the thermal quality of an ectotherm's habitat: a comparison between two temperate lizard species. Funct Ecol. 1997;11:79–89.

    Article  Google Scholar 

  18. 18.

    Diggins CA, Kelly CA, Ford WM. Atypical den use of Carolina northern flying squirrels (Glaucomys sabrinus coloratus) in the southern Appalachian Mountains. Southeast Nat. 2015;14:N44–9.

    Article  Google Scholar 

  19. 19.

    Dingle H. Evolutionary genetics of animal migration. Am Zool. 1991;31:253–64.

    Article  Google Scholar 

  20. 20.

    Dingle H. Migration: the biology of life on the move. Oxford: Oxford University Press; 2014.

    Google Scholar 

  21. 21.

    Dubayah R, Rich PM. Topographic solar radiation models for GIS. Internat J Geograph Inform Systems. 1995;9(4):405–19.

    Google Scholar 

  22. 22.

    Dunham AE, Grant BW, Overall KL. Interfaces between biophysical and physiological ecology and the population ecology of terrestrial vertebrate ectotherms. Physiol Zool. 1989;62:335–55.

    Article  Google Scholar 

  23. 23.

    Duvall D, Goode MJ, Hayes WK, Leonhardt JK, Brown DG. Prairie rattlesnake vernal migration: field experimental analysis and survival value. Natl Geogr Res. 1990;6:457–69.

    Google Scholar 

  24. 24.

    Duvall D, Schuett GW. Straight-line movement and competitive mate searching in prairie rattlesnakes, Crotalus viridis viridis. Anim Behav. 1997;54:329–34.

    CAS  PubMed  Article  Google Scholar 

  25. 25.

    Environmental Systems Research Institute. ArcGIS, version 9.3. Redlands: Environmental Systems Research Institute; 2009.

    Google Scholar 

  26. 26.

    Fitzgerald M, Shine R, Lemckert F, Towerton A. Habitat requirements of the threatened snake species Hoplocephalus stephensii (Elapidae) in eastern Australia. Austral Ecol. 2005;30:465–74.

    Article  Google Scholar 

  27. 27.

    Fu P, Rich PM. A geometric solar radiation model with applications in agriculture and forestry. Comput Electron Agric. 2002;37:25–35.

    Article  Google Scholar 

  28. 28.

    Gannon VPJ, Secoy DM. Seasonal and daily activity patterns in a Canadian population of the prairie rattlesnake, Crotalus viridis viridis. Can J Zool. 1985;63:86–91.

    Article  Google Scholar 

  29. 29.

    Gomez L, Larsen KW, Gregory PT. Contrasting patterns of migration and habitat use in neighboring rattlesnake populations. J Herpetol. 2015;49:1–6.

    Article  Google Scholar 

  30. 30.

    Greenberg DB, McClintock WJ. Remember the third dimension: terrain modeling improves estimates of snake home range size. Copeia. 2008;2008:801–6.

    Article  Google Scholar 

  31. 31.

    Gregory PT. Northern lights and seasonal sex: the reproductive ecology of cool climate snakes. Herpetologica. 2009;65:1–13.

    Article  Google Scholar 

  32. 32.

    Hamilton BT, Nowak EM. Relationships between insolation and rattlesnake hibernacula. West N Am Naturalist. 2009;69:319–28.

    Article  Google Scholar 

  33. 33.

    Harvey DS, Weatherhead PJ. Hibernation site selection by eastern massasauga rattlesnakes (Sistrurus catenatus catenatus) near their northern range limit. J Herpetol. 2006;40:66–73.

    Article  Google Scholar 

  34. 34.

    Harvey JA. Thermal influences on summer habitat use by western rattlesnakes (Crotalus oreganus) in British Columbia. M.Sc. thesis, Thompson Rivers University. Kamloops, BC, Canada. 2015.

    Google Scholar 

  35. 35.

    Hopcraft JGC, Morales JM, Beyer HL, Borner M, Mwangomo E, Sinclair ARE, Olff H, Haydon DT. Competition, predation, and migration: individual choice patterns of Serengeti migrants captured by hierarchical models. Ecol Monogr. 2014;84:355–72.

    Article  Google Scholar 

  36. 36.

    Hsiung AC, Boyle WA, Cooper RJ, Chandler RB. Altitudinal migration: ecological drivers, knowledge gaps, and conservation implications. Biol Rev. 2018;93:2049–70.

    PubMed  Article  Google Scholar 

  37. 37.

    Huey RB, Peterson CR, Arnold SJ, Porter WP. Hot rocks and not-so-hot rocks: retreat-site selection by garter snakes and its thermal consequences. Ecology. 1989;70:931–44.

    Article  Google Scholar 

  38. 38.

    Huey RB, Slatkin M. Cost and benefit of lizard thermoregulation. Q Rev Biol. 1976;51:363–84.

    CAS  PubMed  Article  Google Scholar 

  39. 39.

    Huey RB. Physiological consequences of habitat selection. Am Nat. 1991;137:S91–S115.

    Article  Google Scholar 

  40. 40.

    Iqbal M. An introduction to solar radiation. Don Mills: Academic Press Canada; 1983.

    Google Scholar 

  41. 41.

    James ARC, Stuart-Smith AK. Distribution of caribou and wolves in relation to linear corridors. J Wildl Manag. 2000;64(1):154–9.

    Article  Google Scholar 

  42. 42.

    Jenness JS. Calculating landscape surface area from digital elevation models. Wildl Soc Bull. 2004;32:829–39.

    Article  Google Scholar 

  43. 43.

    Jezkova T, Jaeger JR, Oláh-Hemmings V, Jones KB, Lara-Resendiz RA, Mulcahy DG, Riddle BR. Range and niche shifts in response to past climate change in the desert horned lizard Phrynosoma platyrhinos. Ecography. 2015;38:001–12.

    Article  Google Scholar 

  44. 44.

    Johnson DH. The comparison of usage and availability measurements for evaluating resource preference. Ecology. 1980;61(1):65–71.

    Article  Google Scholar 

  45. 45.

    Jorgensen D. Annual migrations of female prairie rattlesnakes, Crotalus v. viridis, in Alberta. M.Sc. Calgary: Thesis, University of Calgary; 2009.

  46. 46.

    Kearney M, Shine R, Porter WP. The potential for behavioral thermoregulation to buffer “cold-blooded” animals against climate warming. PNAS. 2009;106(10):3835–40.

    CAS  PubMed  Article  Google Scholar 

  47. 47.

    King MB, Duvall D. Prairie rattlesnake seasonal migrations: episodes of movement, vernal foraging and sex differences. Anim Behav. 1990;39(5):924–35.

    Article  Google Scholar 

  48. 48.

    La Sorte FA, Jetz W. Tracking climate niche boundaries under recent climate change. J Anim Ecol. 2012;8(4):914–25.

    Article  Google Scholar 

  49. 49.

    Lomas E, Maida JR, Bishop CA, Larsen KW. Movement ecology of northern Pacific rattlesnakes (Crotalus o. oreganus) in response to disturbance. Herpetologica. 2019;75(2):153–61.

    Article  Google Scholar 

  50. 50.

    Lomas E, Larsen KW, Bishop CA. Persistence of northern Pacific rattlesnakes masks the impact of human disturbance on weight and body condition. Anim Conserv. 2015;18:548–56.

    Article  Google Scholar 

  51. 51.

    Macartney JM, Gregory PT. Reproductive biology of female rattlesnakes (Crotalus viridis) in British Columbia. Copeia. 1988;1988:47–57.

    Article  Google Scholar 

  52. 52.

    Macartney JM. The ecology of the northern Pacific rattlesnake, Crotalus viridis oreganus, in British Columbia. MSc Thesis. Victoria: University of Victoria; 1985.

  53. 53.

    Maida JR, Bishop CA, Larsen KW. Migration and disturbance: impact of fencing and development on rattlesnake spring movements in British Columbia. Can J Zool. 2020;98:1–12.

    Article  Google Scholar 

  54. 54.

    Maida JR, Kirk DA, McKibbin O, Row JR, Larsen KW, Stringham C, Bishop CA. Population estimate, survivorship, and generation time of northern Pacific rattlesnakes (Crotalus o. oreganus) at its northern-most range limits. Herpetol Conserv Biol. 2018;13:662–72.

    Google Scholar 

  55. 55.

    Matsuda BM, Green DM, Gregory PT. Amphibians and reptiles of British Columbia. Victoria: Royal British Columbia Museum; 2006.

    Google Scholar 

  56. 56.

    McAllister JM, Maida JR, Dyer O, Larsen KW. Diet of Roadkilled Western rattlesnakes (Crotalus oreganus) and Gophersnakes (Pituophis catenifer) in southern British Columbia. Northwest Nat. 2016;97(3):181–9.

    Article  Google Scholar 

  57. 57.

    Meidinger DV, Pojar J. Ecosystems of British Columbia. Victoria: Research Branch BC Ministry of Forests; 2011.

  58. 58.

    Milner-Gulland EJ, Fryxell JM, Sinclair ARE. Animal migration: a synthesis. Oxford: Oxford University Press; 2011.

  59. 59.

    Murtaugh PA. Simplicity and complexity in ecological data analysis. Ecology. 2007;88(1):56–62.

    PubMed  Article  Google Scholar 

  60. 60.

    Nellemann C, Fry G. Quantitative analysis of terrain ruggedness in reindeer winter grounds. Arctic. 1995;48:172–6.

    Article  Google Scholar 

  61. 61.

    Nowak EM, Schuett GW, Theimer TC, Sisk TD, Nishikawa K. Does short-term provisioning of resources to prey result in behavioral shifts by rattlesnakes? J Wildlife Manage. 2015;79:357–72.

    Article  Google Scholar 

  62. 62.

    Osipova L, Okjello MM, Njumbi SJ, Ngene S, Western D, Hayward MW, Blakenhol N. Using step-selection functions to model landscape connectivity for African elephants: accounts for variability across individuals and seasons. Anim Conserv. 2018;22:35–48.

    Article  Google Scholar 

  63. 63.

    Parker JM, Anderson SH. Ecology and behavior of the midget faded rattlesnake (Crotalus oreganus concolor) in Wyoming. J Herpetol. 2007;41:41–51.

    Article  Google Scholar 

  64. 64.

    Price-Rees SJ, Brown GP, Shine R. Habitat selection by bluetongue lizards (Tiliqua, Scincidae) in tropical Australia: a study using GPS telemetry. Anim Biotelem. 2013;1:7.

    Article  Google Scholar 

  65. 65.

    R Development Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2011. Accessed 24 Jul 2018.

    Google Scholar 

  66. 66.

    Ramenofsky M, Wingfield JC. Regulation of migration. BioScience. 2007;57:135–43.

    Article  Google Scholar 

  67. 67.

    Reinert HK, Cundall D. An improved surgical implantation method for radiotracking snakes. Copeia. 1982;1982:702–5.

    Article  Google Scholar 

  68. 68.

    Reinert HK. Radiotelemetric field studies of pitvipers: data acquisition and analysis. In: Campbell JA, Brodie Jr ED, editors. Biology of the pitvipers. Tyler: Selva Publishing; 1992. p. 185–97.

    Google Scholar 

  69. 69.

    Riley SJ, DeGloria SD, Elliot R. A terrain ruggedness index that quantifies topographic heterogeneity. Intermt J Sci. 1999;5:23–7.

    Google Scholar 

  70. 70.

    Robson KM. American pika population genetic structure, demographic history, and behavior in an atypical environment. MSc. Thesis. University of British Columbia (Okanagan), Kelowna, BC, Canada. 2013.

    Google Scholar 

  71. 71.

    Row JR, Blouin-Demers G. Thermal quality influences effectiveness of thermoregulation, habitat use, and behaviour in milk snakes. Oecologia. 2006;148:1–11.

    PubMed  Article  Google Scholar 

  72. 72.

    Row JR, Blouin-Demers G. Thermal quality influences habitat selection at multiple spatial scales in milksnakes. Écoscience. 2009;13(4):443–50.

    Article  Google Scholar 

  73. 73.

    Russell AP, Bauer AM, Johnson MK. Migration in amphibians and reptiles: an overview of patterns and orientation mechanisms in relation to life history strategies. In: Elewa AMT, editor. Migration of organisms. Berlin, Heidelberg: Springer; 2005.

    Google Scholar 

  74. 74.

    Sappington JM, Longshore KM, Thompson DB. Quantifying landscape ruggedness for animal habitat analysis: a case study using bighorn sheep in the Mojave Desert. J Wildlife Manage. 2007;71:1419–26.

    Article  Google Scholar 

  75. 75.

    Schaefer WH. Diagnosis of sex in snakes. Copeia. 1934;1934:181.

    Article  Google Scholar 

  76. 76.

    Sears MW, Angilletta MJ Jr, Schuler MS, Borchert J, Dilliplane KF, Stegman M, Rusch TW, Mitchell WA. Configuration of the thermal landscape determines thermoregulatory performance of ectotherms. PNAS. 2016;113(38):10595–600.

    CAS  PubMed  Article  Google Scholar 

  77. 77.

    Shoemaker KT, Gibbs JP. Evaluating basking-habitat deficiency in the threatened eastern Massasauga rattlesnake. J Wildl Manag. 2010;74(3):504–13.

    Article  Google Scholar 

  78. 78.

    Taylor EN, Malaway MA, Browning DM, Lemar SV, DeNardo DF. Effects of food supplementation on the physiological ecology of female western diamond-backed rattlesnakes (Crotalus atrox). Oecologia. 2005;144:206–13.

    PubMed  Article  Google Scholar 

  79. 79.

    Thurjfell H, Ciuti S, Boyce MS. Applications of step-selection functions in ecology and conservation. Mov Ecol. 2014;2:4.

    Article  Google Scholar 

  80. 80.

    Volpe NL, Hadley AS, Robinson WD, Betts MG. Functional connectivity experiments reflect routine movement behavior of a tropical hummingbird species. Ecol Appl. 2014;24:2122–31.

    PubMed  Article  Google Scholar 

  81. 81.

    Waldron JL, Bennett SH, Welch SM, Dorcas ME, Lanham JD, Kalinowsky W. Habitat specificity and home-range size as attributes of species vulnerability to extinction: a case study using sympatric rattlesnakes. Anim Conserv. 2006;9:414–20.

    Article  Google Scholar 

  82. 82.

    Wasko DK, Sasa M. Food resources influence spatial ecology, habitat selection, and foraging behavior in an ambush-hunting snake (Viperidae: Bothrops asper): an experimental study. Zoology. 2012;115:179–87.

    PubMed  Article  Google Scholar 

  83. 83.

    Weatherhead PJ, Prior KA. Preliminary observations of habitat use and movements of the eastern Massasauga rattlesnake (Sistrurus c. catenatus). J Herpetol. 1992;26:447–52.

    Article  Google Scholar 

  84. 84.

    Webb JK, Shine R. Using thermal ecology to predict retreat-site selection by an endangered snake species. Biol Conserv. 1998;86:233–42.

    Article  Google Scholar 

  85. 85.

    Williams KE, Hodges KE, Bishop CA. Small reserves around hibernacula sites many not adequately protect mobile snakes: the example of Great Basin Gophersnakes (Pituophis catenifer deserticola) in British Columbia. Can J Zool. 2012;90:304–12.

    Article  Google Scholar 

  86. 86.

    Wills CA, Beaupre SJ. An application of randomization for detecting evidence of thermoregulation in timber rattlesnakes (Crotalus horridus) from northern Arkansas. Physiol Biochem Zoology. 2000;73:325–34.

    CAS  Article  Google Scholar 

  87. 87.

    Zar JH. Biostatistical Analysis. 3th ed. Prentice Hall: Upper Saddle River; 1999.

  88. 88.

    Zimmerman DW, Zumbo BD. Can percentiles replace raw scores in the statistical analysis of test data? Educ Psychol Meas. 2005;65:616–38.

    Article  Google Scholar 

Download references


We thank the BC Wildlife Park as well as veterinarians David Sedgmann, Bruce Maricle, and Malcolm McAdie for their expertise and facilities. We also thank field assistants Andy Greschner, Melany Leontowich, and Cara Haywood-Farmer. The anonymous reviewers and editor (R. Kays) were extremely helpful and diligent in providing useful comments and suggestions seeking clarity of our methods.


The Grasslands Conservation Council of British Columbia, the Habitat Conservation Trust Fund, and the British Columbia Government’s Future Forest Ecosystems Science Committee funded this work.

Author information




This work was part of the MSc thesis research of JAH under the supervision of KWL. The work was conceived by KWL with JAH making modifications and selecting sites. JAH analyzed the data. Both authors were involved in the writing of this manuscript and the final approval for publication.

Corresponding author

Correspondence to Karl W. Larsen.

Ethics declarations

Ethics approval and consent to participate

Surgeries and field procedures were performed under Protocols AUP 2010–05 and AUP 2011-04R and Standard Operating Procedure SOP 2010–02 as approved by the Animal Ethics Committee at Thompson Rivers University. General Wildlife Permits KA10–61565 and VI10–61565 were issued by the British Columbia Ministry of Environment. No consent to participate was required.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Harvey, J.A., Larsen, K.W. Rattlesnake migrations and the implications of thermal landscapes. Mov Ecol 8, 21 (2020).

Download citation


  • Thermal ecology
  • Crotalus oreganus
  • Ectotherm
  • Incident solar radiation
  • Thermoregulation
  • Habitat selection
  • Random walk
  • Migration strategies