Rattlesnake migrations and the implications of thermal landscapes
Movement Ecology volume 8, Article number: 21 (2020)
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 . 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 . 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 , 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 . 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 . 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 . 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 . 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 ; 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  such as food and access to mates [7, 22]. Thermoregulation thus may be considered a form of ‘resource utilization’  and variation of heat as a resource across different habitats likely will have repercussions .
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 . 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  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.
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 . 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 . 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 ). 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  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 . 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.  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 . 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 ; 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).
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 . Because females vary in their migratory behaviour according to their reproductive status (see above - ), 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  with modifications by Reinert . 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 . 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  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 . 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].
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 - ) target situations where massive amounts of sequential spatial data are collected on animals, primarily through GPS (e.g. ), other remote tracking methods, or constant tracking  – 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 .
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 . We used Hawth’s tools  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  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 . 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  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.
All statistical analyses were performed in the program R version 2.12.1 . Data were tested for normality by examination of histograms and using the Shapiro-Wilk test or the Kolmogorov-Smirnov test . 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).
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).
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).
The ruggedness index values for migration paths were significantly higher for Forest snakes (x̄ =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 , 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 . 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.  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 , bighorn sheep , badger  and grouse , 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 . 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 . 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 . 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 .
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
Universal Transverse Mercator
Apps CD, Newhouse NJ, Kinley TA. Habitat associations of American badgers in southeastern British Columbia. Can J Zool. 2002;80:1228–39.
Ashton KG. Movements and mating behavior of adult male midget faded rattlesnakes, Crotalus oreganus concolor, in Wyoming. Copeia. 2003;2003:190–4.
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.
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.
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. http://www.bcreptiles.ca/docs/bertram%20et%20al%202001.pdf Accessed 24 Jul 2018.
Beyer HL. 2004. Hawth's analysis tools for ArcGIS. 2004 http://wwwspatialecologycom/htools Accessed 5 May 2010.
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.
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.
British Columbia Conservation Data Centre. Endangered Species and Ecosystems – Confidential Sensitive Occurrences (digital file). Victoria; 2009. https://www2.gov.bc.ca/gov/content/environment/plants-animals-ecosystems/conservation-data-centre.
Brown JR, Bishop CA, Brooks RJ. Effectiveness of short-distance translocation and its effects on western rattlesnake. J Wildl Manag. 2009;73:419–25.
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.
Carpenter J, Aldridge C, Boyce MS. Sage-grouse habitat selection during winter in Alberta. J Wildl Manag. 2010;74:1806–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; 1993a. http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.71.7160&rep=rep1&type=pdf Accessed 24 Jul 2018.
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.
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.
Crawley MJ. The R Book. West Sussex, London: Wiley; 2007.
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.
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.
Dingle H. Evolutionary genetics of animal migration. Am Zool. 1991;31:253–64.
Dingle H. Migration: the biology of life on the move. Oxford: Oxford University Press; 2014.
Dubayah R, Rich PM. Topographic solar radiation models for GIS. Internat J Geograph Inform Systems. 1995;9(4):405–19.
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.
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.
Duvall D, Schuett GW. Straight-line movement and competitive mate searching in prairie rattlesnakes, Crotalus viridis viridis. Anim Behav. 1997;54:329–34.
Environmental Systems Research Institute. ArcGIS, version 9.3. Redlands: Environmental Systems Research Institute; 2009.
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.
Fu P, Rich PM. A geometric solar radiation model with applications in agriculture and forestry. Comput Electron Agric. 2002;37:25–35.
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.
Gomez L, Larsen KW, Gregory PT. Contrasting patterns of migration and habitat use in neighboring rattlesnake populations. J Herpetol. 2015;49:1–6.
Greenberg DB, McClintock WJ. Remember the third dimension: terrain modeling improves estimates of snake home range size. Copeia. 2008;2008:801–6.
Gregory PT. Northern lights and seasonal sex: the reproductive ecology of cool climate snakes. Herpetologica. 2009;65:1–13.
Hamilton BT, Nowak EM. Relationships between insolation and rattlesnake hibernacula. West N Am Naturalist. 2009;69:319–28.
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.
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.
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.
Hsiung AC, Boyle WA, Cooper RJ, Chandler RB. Altitudinal migration: ecological drivers, knowledge gaps, and conservation implications. Biol Rev. 2018;93:2049–70.
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.
Huey RB, Slatkin M. Cost and benefit of lizard thermoregulation. Q Rev Biol. 1976;51:363–84.
Huey RB. Physiological consequences of habitat selection. Am Nat. 1991;137:S91–S115.
Iqbal M. An introduction to solar radiation. Don Mills: Academic Press Canada; 1983.
James ARC, Stuart-Smith AK. Distribution of caribou and wolves in relation to linear corridors. J Wildl Manag. 2000;64(1):154–9.
Jenness JS. Calculating landscape surface area from digital elevation models. Wildl Soc Bull. 2004;32:829–39.
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.
Johnson DH. The comparison of usage and availability measurements for evaluating resource preference. Ecology. 1980;61(1):65–71.
Jorgensen D. Annual migrations of female prairie rattlesnakes, Crotalus v. viridis, in Alberta. M.Sc. Calgary: Thesis, University of Calgary; 2009.
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.
King MB, Duvall D. Prairie rattlesnake seasonal migrations: episodes of movement, vernal foraging and sex differences. Anim Behav. 1990;39(5):924–35.
La Sorte FA, Jetz W. Tracking climate niche boundaries under recent climate change. J Anim Ecol. 2012;8(4):914–25.
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.
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.
Macartney JM, Gregory PT. Reproductive biology of female rattlesnakes (Crotalus viridis) in British Columbia. Copeia. 1988;1988:47–57.
Macartney JM. The ecology of the northern Pacific rattlesnake, Crotalus viridis oreganus, in British Columbia. MSc Thesis. Victoria: University of Victoria; 1985.
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.
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.
Matsuda BM, Green DM, Gregory PT. Amphibians and reptiles of British Columbia. Victoria: Royal British Columbia Museum; 2006.
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.
Meidinger DV, Pojar J. Ecosystems of British Columbia. Victoria: Research Branch BC Ministry of Forests; 2011.
Milner-Gulland EJ, Fryxell JM, Sinclair ARE. Animal migration: a synthesis. Oxford: Oxford University Press; 2011.
Murtaugh PA. Simplicity and complexity in ecological data analysis. Ecology. 2007;88(1):56–62.
Nellemann C, Fry G. Quantitative analysis of terrain ruggedness in reindeer winter grounds. Arctic. 1995;48:172–6.
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.
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.
Parker JM, Anderson SH. Ecology and behavior of the midget faded rattlesnake (Crotalus oreganus concolor) in Wyoming. J Herpetol. 2007;41:41–51.
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.
R Development Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2011. http://www.R-project.org/. Accessed 24 Jul 2018.
Ramenofsky M, Wingfield JC. Regulation of migration. BioScience. 2007;57:135–43.
Reinert HK, Cundall D. An improved surgical implantation method for radiotracking snakes. Copeia. 1982;1982:702–5.
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.
Riley SJ, DeGloria SD, Elliot R. A terrain ruggedness index that quantifies topographic heterogeneity. Intermt J Sci. 1999;5:23–7.
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.
Row JR, Blouin-Demers G. Thermal quality influences effectiveness of thermoregulation, habitat use, and behaviour in milk snakes. Oecologia. 2006;148:1–11.
Row JR, Blouin-Demers G. Thermal quality influences habitat selection at multiple spatial scales in milksnakes. Écoscience. 2009;13(4):443–50.
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.
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.
Schaefer WH. Diagnosis of sex in snakes. Copeia. 1934;1934:181.
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.
Shoemaker KT, Gibbs JP. Evaluating basking-habitat deficiency in the threatened eastern Massasauga rattlesnake. J Wildl Manag. 2010;74(3):504–13.
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.
Thurjfell H, Ciuti S, Boyce MS. Applications of step-selection functions in ecology and conservation. Mov Ecol. 2014;2:4.
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.
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.
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.
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.
Webb JK, Shine R. Using thermal ecology to predict retreat-site selection by an endangered snake species. Biol Conserv. 1998;86:233–42.
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.
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.
Zar JH. Biostatistical Analysis. 3th ed. Prentice Hall: Upper Saddle River; 1999.
Zimmerman DW, Zumbo BD. Can percentiles replace raw scores in the statistical analysis of test data? Educ Psychol Meas. 2005;65:616–38.
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.
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
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Harvey, J.A., Larsen, K.W. Rattlesnake migrations and the implications of thermal landscapes. Mov Ecol 8, 21 (2020). https://doi.org/10.1186/s40462-020-00202-0
- Thermal ecology
- Crotalus oreganus
- Incident solar radiation
- Habitat selection
- Random walk
- Migration strategies