IMPACTS OF FENCING AND DEVELOPMENT ON WESTERN RATTLESNAKE (CROTALUS OREGANUS) SPRING MOVEMENTS IN BRITISH COLUMBIA by JARED RYAN MAIDA BNRSc, Thompson Rivers University, 2014 A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF MASTER OF SCIENCE IN ENVIRONMENTAL SCIENCE In the Department of Natural Resource Science Thompson Rivers University Kamloops, British Columbia, Canada December, 2018 Thesis examining committee: Karl Larsen (PhD), Professor and Thesis Supervisor, Department of Natural Resource Sciences, Thompson Rivers University Christine Bishop (PhD), Adjunct Faculty and Thesis Co-Supervisor, Research Scientist, Environment and Climate Change Canada, Wildlife Research Division, Government of Canada Mark Rakobowchuck (PhD), Assistant Professor and Committee Member, Biological Sciences, Thompson Rivers University Mike Cardwell (MSc), External Examiner, Adjunct Researcher, Department of Biology, San Diego State University i Thesis Supervisor(s): Professor Karl Larsen (PhD) and Adjunct Faculty Christine Bishop (PhD) ABSTRACT Due to increasing anthropogenic pressures, the process of natural animal migration is undergoing alterations across many different taxa. To investigate these impacts is difficult on large-scale migrations, yet small-scale migrants will provide useful and tractable systems to understand the effects of disturbance and landscape barriers on natural movement patterns and migrations. The Western Rattlesnake (Crotalus oreganus) in British Columbia is a small, migratory predator that undertakes annual spring movements from communal hibernaculum to summer hunting grounds. In this thesis, I examined changes to these migration movements in male rattlesnakes that encountered both snake exclusion fencing barriers and disturbed habitats from 2011-2016. My thesis field work included 2015 and 2016, and data from 20112014 was drawn from a long-term project database. Individuals moving through disturbed habitats or intercepted by snake exclusion fencing demonstrated shorter migration distances and reduce spring path sinuosity (more crooked routes) compared to individuals migrating in undisturbed habitats. Specifically, individuals encountering a snake fence during spring movements had shorter total spring migration path lengths and smaller home ranges. Regardless, migration distance was strongly associated with the distance individuals travelled until they first encountered human disturbance. Despite this difference, duration of the migration did not differ between rattlesnakes that encountered a snake exclusion fence, experienced other forms of disturbance, or migrated through undisturbed habitats. Overall, my work reveals the importance of looking closely at seasonal movement patterns, including migration, of smaller-scale migrants in response to abrupt and long-term landscape changes and barriers. Small-scale migrants can be used as an important indicator of landscape and community connectivity and can be used as an important reference for managing ecosystem function and health on a changing, anthropocentric landscape. Keywords: British Columbia, Conservation, Crotalus oreganus, Migration, Radio-telemetry, Species-at-Risk ii TABLE OF CONTENTS ABSTRACT .................................................................................................................................. I TABLE OF CONTENTS ................................................................................................................. II LIST OF FIGURES ....................................................................................................................... IV LIST OF TABLES ....................................................................................................................... VII ACKNOWLEDGEMENTS ............................................................................................................ VIII CHAPTER ONE: ANIMALS, BEHAVIOUR AND DISTURBANCE....................................................1 DISTURBANCE AND WILDLIFE ............................................................................................................ 1 DISTURBANCE AND SNAKES ............................................................................................................... 3 SNAKES AND FENCING ........................................................................................................................ 5 THE WESTERN RATTLESNAKE IN BRITISH COLUMBIA –DEALING WITH DISTURBANCE ................... 7 OSOYOOS INDIAN RESERVE (OIR) AND THE NK’MIP SNAKE PROGRAM ......................................... 10 THESIS OBJECTIVES .......................................................................................................................... 18 LITERATURE CITED........................................................................................................................... 21 CHAPTER TWO: IMPACTS OF FENCING AND DISTURBANCE ON WESTERN RATTLESNAKE (CROTALUS OREGANUS) SPRING MOVEMENTS ........................................................................32 INTRODUCTION ................................................................................................................................. 32 METHODOLOGY ................................................................................................................................ 36 Study Site...................................................................................................................................... 36 Radio Telemetry ........................................................................................................................... 36 Spring Migration Parameters ........................................................................................................ 37 Home Range Parameters .............................................................................................................. 38 Body Condition and Growth ........................................................................................................ 40 Quantifying Migration Groups ..................................................................................................... 40 Prey Population Monitoring ......................................................................................................... 41 Statistical Analysis ....................................................................................................................... 43 RESULTS ........................................................................................................................................... 43 Migration Parameters ................................................................................................................... 44 Home Range Parameters .............................................................................................................. 48 Body Characteristics ..................................................................................................................... 48 Prey Population Monitoring ......................................................................................................... 52 iii DISCUSSION ...................................................................................................................................... 56 LITERATURE CITED........................................................................................................................... 61 CHAPTER THREE: CONCLUSION AND MANAGEMENT IMPLICATIONS ...................................70 SUMMARY ......................................................................................................................................... 70 CONSERVATION AND MANAGEMENT IMPLICATIONS ....................................................................... 71 KNOWLEDGE GAPS AND FUTURE RESEARCH CONSIDERATIONS ..................................................... 73 CONCLUSION .................................................................................................................................... 76 LITERATURE CITED........................................................................................................................... 77 iv LIST OF FIGURES Figure 1.1. Western Rattlesnake (Crotalus oreganus) range in the Southern Interior of British Columbia (B.C.) with inset showing the species’ northern range limits in North America (Produced by Robyn Reudink). Osoyoos, B.C. (indicated by red star) was the location of my study site (see below). B.C. map produced by Rene McKibbin……………………………………………………………………………...9 Figure 1.2. Above: Photograph of the Osoyoos valley looking north over Osoyoos Lake in 1935 (photograph from the Okanagan Archive Society; www.oldphotos.ca); Below: same area in 2012 (photograph from Emily Lomas)………………………..……….11 Figure 1.3. Map of the Osoyoos Indian Reserve (OIR) study site (indicated in orange) near Osoyoos, B.C., Canada. The northern and eastern portions remained undeveloped with developed areas concentrated to the southern and western portions of the study site. The condominium community and associated snake fence is located at the southeast corner of the study site. Inset shows the northern extent of the Western Rattlesnake (Crotalus oreganus) range with study site location (indicated by star)…………………………………………………………………………………..12 Figure 1.4. Examples of the dominant ecotypes on the Osoyoos Indian Reserve (OIR): (A) open Ponderosa Pine (Pinus ponderosa) and Douglas-fir (Pseudotsuga menziezii) forests, (B) Shrub-steppe habitat dominated by Big Sagebrush (Artemesia tridentata) and Antelope Brush (Purshia tridentata), (C) and mid-elevation rocky slopes……..13 Figure 1.5. Contrast in habitat quality and landscape disturbance between the southern developed (right) and northern undisturbed (left) portions of the Osoyoos Indian Reserve (OIR) study site. The northern portions containing no habitat barriers and limited human presence. Conversely, disturbance within the southern portions were geared towards tourism and the landscape fragmented by roads, condominiums, vineyards, fencing, golf course, hiking trials and a campground…...…………….....14 Figure 1.6. Mean monthly temperatures in Osoyoos, B.C. during the years of this study (2011– 2016) compared to the historical 30-year (1981–2010) mean monthly temperatures (Environmental Canada 2018). The Osoyoos Indian Reserve (OIR) study site was located approximately 2 km north-west of the location where the weather data was recorded…………………………………………………………...15 Figure 1.7. Timeline showing progression of development and habitat loss in the southern section of the Osoyoos Indian Reserve (OIR) study site from 2004–2016. The campground, golf course and vineyards were established in 2004, and 2004–2016 saw the construction of the Spirit Ridge Resort, Nk’Mip Dersert Cultural Centre, hydro substation, expansion of roads, 4 km long snake exclusion fence, expansion of campground and construction of the condominium community………..…………...17 v Figure 1.8. Condominium community and construction surrounded by a wooden, permanent snake exclusion fence………………………………………………………..………20 Figure 2.1. Map of the southern, developed portion of the Osoyoos Indian Reserve (OIR) study site near Osoyoos, B.C. Development included a golf course, vineyards, resort, cultural centre, hiking trails, campground, condominium community and construction, and associated roads and parking lots. Snake fencing represented in red. Small mammal trapping gird locations used from 2012–2016 are as follows: (1) golf course, (2) hillside, (3) vineyard, and (4) undisturbed shrub-steppe..……………….42 Figure 2.2. Mean ± SE of log-transformed migration distance (MD) of the three migration categories (top) and relationship between log-transformed migration distance and distance to disturbance (bottom) of male Western Rattlesnakes (Crotalus oreganus) on the Osoyoos Indian Reserve (OIR) near Osoyoos, B.C., Canada. Sample sizes for the categories are FENCE (n= 7), DSTB (n= 10), and UN (n= 10). Means with different letters indicate a significant difference at α= 0.05. Individuals in the FENCE category encountered the condominium snake fence during their spring migration, while DSTB snakes encountered other forms of disturbance. Thirdly, UN snakes were individuals in the northern portion of the study site that did not encounter any type of disturbance during spring migration or throughout entire active season….................46 Figure 2.3. Mean ± SE of migration path sinuosity (MPS) of the three migration categories (top) and the relationship between migration sinuosity and distance to disturbance (bottom) of male Westner Rattlesnakes (Crotalus oreganus) on the Osoyoos Indian Reserve (OIR) near Osoyoos, B.C., Canada. Sample sizes for each group are FENCE (n= 7), DSTB (n=10), and UN (n=10). Means with different letters indicate a significant difference at α= 0.05. See Figure 2.2 caption for description of migration categories………………………………………………………………..…………...47 Figure 2.4. Mean ± SE of log-transformed minimum convex polygon (MCP) home range estimates of the three migratory categories (top) and the relationship between logtransformed MCP and distance to disturbance (bottom) of male Western Rattlesnakes (Crotalus oreganus) on the Osoyoos Indian Reserve (OIR) near Osoyoos, B.C., Canada. Sample sizes for the migration categories are FENCE (n= 7), DSTB (n= 10), and UN (n= 10). Means with different letters indicate a significant difference at α= 0.05. See Figure 2.2 caption for description of migration categories..………………49 Figure 2.5. Size class distribution of male Western Rattlesnake (Crotalus oreganus) captures (n= 489) on the Osoyoos Indian Reserve (OIR) near Osoyoos, B.C., Canada from 2011-2016……………………………………………………………………………50 vi Figure 2.6. Mean ± SE of spring snout-vent-length (SVL) between the three migrations categories (top) and the relationship between spring SVL and distance to disturbance (bottom) of male Western Rattlesnakes (Crotalus oreganus) on the Osoyoos Indian Reserve (OIR) near Osoyoos, B.C., Canada. Samples sizes of each category are FENCE (n= 7), DSTB (n= 10), and UN (n= 10). Means with different letters indicate a significant difference at α= 0.05. See Figure 2.2 caption for description of migration categories…………………………………………………………………………….51 Figure 2.7. Mean ± SE instantaneous growth rate between the three migration categories (top) and the relationship between instantaneous growth rate and distance to disturbance (bottom) of male Western Rattlesnakes (Crotalus oreganus) on the Osoyoos Indian Reserve (OIR) near Osoyoos, B.C., Canada. Sample sizes for each category are FENCE (n= 6), DSTB (n=10), and UN (n=9). Means with different letters indicate a significant difference at α= 0.05. See Figure 2.2 caption for description of migration categories…………………………………………………..53 vii LIST OF TABLES Table 2.1. Spring migration and home range parameter descriptions and overview of methods used to calculate each parameter for Western Rattlesnakes from 2011-2016 on the Osoyoos Indian Reserve (OIR) near Osoyoos, B.C……………………..…………...39 Table 2.2. P-values for Pearson’s correlation tests between male Western Rattlesnake (Crotalus oreganus) end of year body condition and instantaneous growth and each migration parameter used (MD, MPL, DTD, DAYS, and MPS)……………………54 Table 2.3. Small mammal population density estimates (animals/ha) of each trapping grid on the Osoyoos Indian Reserve (OIR) near Osoyoos, B.C., Canada from 2012-2016. Small mammals represented within the population were Great Basin Pocket Mice (Perognathus parvus) and Deer Mice (Peromyscus maniculatus)…...……………...55 viii ACKNOWLEDGEMENTS I am eternally grateful for the support I received throughout the process of completing this thesis; how individuals stepped up during times of personal crisis was overwhelming. Firstly, I’d like to thank my supervisor, Dr. Karl Larsen, for his unwavering support, guidance, friendship and his relentless passion and excitement of rattlesnakes and wildlife research. I would also like to thank Dr. Christine Bishop for her logistical support and expertise as well as my third committee member Dr. Mark Rakobowchuk for his assistance. Thank you to my external examiner, Mike Cardwell, for his time and insight through the editorial and defense process. Hiking in the desert all day, every day is not easy, but my field assistants made those long days enjoyable. Thank you Cole Hooper, Kirsten Hales, Kirsten Eyolfson, Jim Song and Ryan Richter for your hard work, passion and putting up with me constantly singing to myself throughout the day. Without your hard work and dedication this research would not have been possible. I will forever appreciate my time at the Nk’Mip Desert Cultural Centre. I will look back on countless memories and lifelong friendships made. I always felt welcome and (more importantly) also discovered my love for bannock. I met and worked with so many amazing people. I would especially like to thank John Herbert and Barb Sabyan who always seemed to have time to offer assistance or cater to needs or requests I may have had. Thank you to Charlotte Stringam and Margaret Holm for their logistical support and communicating with OIB council. Thank you to the Osoyoos Indian Band for allowing me access to an amazing field site and their continued leadership and dedication to snake conservation and education. I would also like to the Dr. Henry Kleinhofmeyer and the rest of the awesome staff at the South Okanagan Vetrinary Clinic for their availability and welcoming energy (and at times delicious baked goods) during long hours of snake surgeries. To my fellow graduate students and band of misfits of the TRU Research Centre: you are awesome. To Michelle Phillips, Roman Chapman, Hilary Hamilton, Dan Doutaz and Doug Terspsma, I will always remember Friday afternoon ramen, patio beers, dog birthday parties and making hours of struggling with R somehow bearable! I (“mod”) look forward to a lifelong of friendship with you all. To my amazing friend, Stephanie Winton, I could not ix have picked a better partner to go on this rollercoaster of a ride with, thanks for the countless memories. Thank you to my incredible family and friends for their support throughout my graduate school journey. Finally, to Valerie Blow, you are my inspiration. I cannot thank you enough for your unwavering support and patience throughout this long process (not to mention having to listen to me obsessively talk about snakes and stats all the time). This work was supported by Environment and Climate Change Canada, Bellstar Hotels and Resorts, Aboriginal Fund for Species at Risk, Interdepartmental Recovery Fund, the Natural Sciences and Engineering Research Council of Canada, the Osoyoos Indian Band, and Nk’Mip Desert Cultral Centre. All field-work, handling procedures and surgeries for my thesis work were performed under the following permits: PE15-164975, PE15171661, SARA-PYR-2015-0316 and Thompson Rivers University Animal Care Protocol Numbers 100343 and 100344. THANKS FOR THE MEMORIES. 1 CHAPTER ONE: ANIMALS, BEHAVIOUR AND DISTURBANCE DISTURBANCE AND WILDLIFE With the human population continuing to increase, habitat loss, fragmentation and alteration have emerged as the dominant threats to global biodiversity and wildlife species persistence (Wilcove et al. 1998; Gibbons et al. 2000; Saha et al. 2018). Habitat loss and fragmentation due to anthropogenic activities are responsible for the majority of the rapidly increasing extinction events globally, and remain a primary cause for the displacement and elimination of local wildlife species (Czech et al. 2000; Kjoss and Litvaitis 2001; Marzluff 2001; Mckinney 2002). In Canada, federal legislations exist to protect species, specifically those at risk, from human disturbance. However, urban sprawl, resource extraction, tourism and outdoor recreation continue to present novel stressors and selection pressures on wildlife populations. Individuals or species unable to cope with these novel pressures may suffer both short and long-term consequences including direct mortality and/or indirect consequences (Nattrass and Lusseau 2016). The direct effects (mortality) of disturbance and habitat loss often are readily apparent, and at times, simple to quantify. Road mortality (Souza et al. 2015), illegal and/or over-harvesting (Wittemyer et al. 2014) and human persecution (Venter et al. 2006) are among the obvious sources of mortality for wildlife populations globally, and such impacts may result in the immediate extirpation of a species. Conversely, indirect effects of disturbance on wildlife often are challenging to quantify and understand, but may represent severe, long-term consequences on individual health and species persistence (Ohashi et al. 2013). For example, Cougars (Puma concolor) avoid roadways and trails during high traffic seasons, likely altering resource selection (Morrison et al. 2014). Bird population densities decrease along busy roadways, and males in these areas may have difficulties attracting mates (Reijnen et al. 1997; Kociolek et al. 2011). In the presence of ecotourism activity, Woodland Caribou (Rangider tarandus caribou) are more vigilant at the expense of time spent resting and foraging (Duchesne et al. 2000; Fortin et al. 2016). Repeated alterations in normal behaviour due to disturbance, such as presented here, may have negative impacts on 2 individuals - especially when the individuals affected have difficulties compensating for potential lost opportunities (Nattrass and Lusseau 2016). Behavioural alterations in response to disturbance or human presence can result in long-term biologically-significant impacts on individuals, and therefore populations (Gill et al. 2001; Christiansen and Lusseau 2015). These impacts may include elevated stress levels (e.g., Brown Bears, Ursus arctus - Støen et al. 2015), increased nestling mortality (e.g., European Storm Petrels, Hydrobates pelagicus - Watson et al. 2014), reduced juvenile growth and body mass (e.g., Hoatzins, Opisthocomus hoazin - Müllner et al. 2004), and negative effects of genetic drift through a reduction in dispersal (e.g., Spanish Imperial Eagle, Aquila adalberti - Martínez-Cruz et al. 2004). Therefore, the persistence per se of animals on disturbed landscapes may not always be indicative of long-term population health and stability (Lomas et al. 2015). Despite recent research, there are very few clear patterns in the behavioural (and physiological) responses animals demonstrate to disturbed habitats and direct human interactions (Tablado and Jenni 2017). This lack of clarity may presumably reflect each species’ specific evolutionary adaptations and furthermore, not all species appear to be negatively impacted by disturbance (Devictor et al. 2008). How animals respond and behave on disturbed landscapes is predicated by disturbance type, species-specific tolerance and life history, as well as the local conditions (Madsen 1988; Béchet et al. 2004; Marini et al. 2017). Specifically, niche breadth has been hypothesized as the main contributor in an animal’s ability to respond to disturbed and changing landscapes (Kitahara and Fujii 1994; Vásquez and Simberloff 2001). Species that follow generalist strategies are typically better suited to adapt to novel environments created by disturbance due to their broad habitat preferences, behavioural flexibility, and/or malleable life history strategies (Kark et al. 2007; Lowry et al. 2013; Šálek et al. 2014). Conversely, specialists rely on a narrow range of habitat requirements and resources (Rousseau et al. 2012; Matthews et al. 2014), and/or dispersal and migration to locate specific resources and avoid genetic isolation (Martin and Fahrig 2018). These specialists typically are more prone to negative consequences when facing novel environments and anthropogenic pressures. The persistence (and at times influx) of species 3 may mask the decline of specialists in fragmented and disturbed areas (Matthews et al. 2014), providing more evidence that species persistence should not be a sole indicator of population or community health. Rather, tracking behavioural responses of individuals as an index of disturbance is extremely important for wildlife managers. However, assessing behavioural responses alone may have a number of potential limitations and may not identify the underlying issues (Beale and Monaghan 2004). Overall, understanding the adaptive limits and behavioural niches of a species is incredibly important to consider in wildlife conservation and management. DISTURBANCE AND SNAKES Snakes are relatively small-bodied, meso-predators that typically specialize in annual site fidelity to specific habitats throughout the course of a year, or over the course of their lifetime (Waldron et al. 2013; Bauder et al. 2015; Gomez et al. 2015). It is likely advantageous for individuals to exhibit fidelity to areas where resource availability and habitats are stable through time (Waldron et al. 2013), especially if the benefits of reusing these areas exceed the cost of searching for resources in other areas (Switzer 1993). Further, snakes as ectotherms rely directly on the thermal properties of their environment to ensure natural physiological processes (i.e., hibernation, digestion, and reproduction; Putmann and Clark 2017). To re-locate and use areas of value (i.e., hibernacula, migration corridor, summer foraging habitat), snakes have well-developed orientation abilities (Landreth 1973; Larsen 1987; Lawson 1991), and their movement patterns may change substantially if reorientation is needed after natural movement patterns are blocked and/or individuals are displaced (Butler et al. 2005; Brown et al. 2009). Species that exhibit high spatial fidelity are associated with stable habitats and are likely not suited to effectively respond to abrupt changes on the landscape. For example, adult snakes have typically demonstrated the inability to re-orientate themselves on the landscape in response to habitat removal and/or barriers limiting movement and access to desired areas (Waldron et al. 2013). Northern snakes likely are particularly vulnerable to these abrupt anthropogenic-influenced changes on the landscape due to the restricted and shortened active season (Gregory 2007). Snake migrations and seasonal movements appear associated with accessing and/or searching for prey (King and Duvall 1990; Bauder et al. 4 2015). Therefore, abrupt changes in the environment may have serious repercussions: short active seasons likely do not leave enough time for individuals to adjust to abrupt landscape changes and still successfully find prey (and mates, and effective hibernacula) compared to their southern counterparts. Furthermore, most northern snake species exhibit life history attributes that include delayed maturation, high adult survival, and infrequent parturition (Macartney and Gregory 1988; Gregory 2007; Maida et al. 2018); novel and abrupt changes in anthropogenic-influenced landscapes can abruptly change selection pressures, and these animals cannot rapidly respond to sudden, dynamic changes (Waldron et al. 2013). In northern regions, snakes den communally and undertake seasonal migrations primarily due to the long, cold winters they face (for example: Western Rattlesnakes in British Columbia, see below). Communal denning and migration between overwintering habitat (hibernacula) and summer foraging sites are closely linked behaviours associated with cold climates and scarce suitable hibernacula habitat (Gregory 2007; Bauder et al. 2015). In contrast, snakes at more southern latitudes typically show reduced movement distances and more flexible overwintering behaviours (Bauder et al. 2015). Within these northern regions, the lack of sufficient summer resources (i.e., prey availability) in the proximity of hibernacula obligates individuals to undertake seasonal migrations between overwintering sites and summer activity areas. Therefore, novel barriers on the landscape, or habitat loss, that interferes with these important movements may impact the survival (much less fitness) of animals. The impacts of human disturbance on snake populations are most apparent through road mortality (Hartmann et al. 2011; Garrah et al. 2015; Marsh and Jaeger 2015; Meek 2015) and accidental or intentional killing from humans (Bonnet et al. 1999; Gomez et al. 2004; Pitts et al. 2017). However, disturbance may have other indirect behavioural and physiological implications for snakes. For example, snake behavioural patterns change in the presence of roads and golf course fairways (Shine et al. 2004; Robson and Blouin-Demers 2013; Goode 2010). Further, rattlesnake home range sizes and migration path straightness tend to become altered in fragmented landscapes (Breininger et al. 2011; Martin et al. 2017). Roads limit dispersal and genetic diversity (Clark et al. 2010) and may cause alterations to physiological stress responses (Owen et al. 2014). Increased human interactions and outdoor 5 recreation activities can increase concealment behaviour by snakes (Beale et al. 2016) and reduce movement frequencies (Parent and Weatherhead 2000). Urban snakes also are less likely to have prey present in their stomachs and are relatively lighter and smaller than nonurban snakes (Wolfe et al. 2017). Lastly, landscape disturbances result in reduced growth rates, fecundity and lower offspring body conditions in rattlesnakes (Jenkins et al. 2009). Snakes are an important component of the natural communities they reside in and have strong effects on community dynamics (Lind et al. 2005), yet they are often viewed as a nuisance animal in urban and other human-dominated environments (Pitts et al. 2017); not surprisingly, persecution has been identified as a leading threat to a number of at-risk snake species. In Canada, and the rest of North America, human social issues play a large role in the conservation of snakes (Gregory 2007). Fear and misconceptions for snakes compounds the challenges of conservation initiatives for these animals as they typically receive limited sympathy and support compared to more charismatic animals (Gomez et al. 2004; Ballouard et al. 2012). Monitoring snake populations is hampered by their small size, cryptic and nocturnal behaviour, and often low-population densities (Durso et al. 2011; Ward et al. 2017). Despite this, snakes on the landscape are critically important for ecosystem health and biodiversity. Specifically, as meso-predators, snakes are a crucial component of energy flow in the complex food web, and can be used as an important indicator species to monitor ecosystem health and function (Lind et al. 2005; Beaupre and Douglas 2009). SNAKES AND FENCING Historically, wildlife fencing has predominantly been used to reduce large mammalvehicle collisions on major highways (Sawyer et al. 2012; Huijser et al. 2016). These types of fencing structures are typically very large (e.g. 2.4 m high for large ungulates), but do not inhibit the passage of small animals like snakes, and can be used over a large scale (Evans and Wood 1980). However, more recently, the use of snake exclusion fencing has been used to mitigate negative interactions between humans and snakes and other reptiles (BaxterGilbert et al. 2015; Jackson et al. 2015; Colley et al. 2017; Markle et al 2017), and have often proved to be effective at reducing mortality rates. For small-bodied animals such as snakes that are highly susceptible to mortality on roads and other human-influenced areas, fencing is valuable for the persistence and viability 6 of local populations. Yet these structures still pose an abrupt, physical barrier for individuals. Encountering abrupt barriers such as snake exclusion fences likely has a large impact on snake behaviour, especially in areas where individuals must annually migrate between overwintering hibernacula and summer foraging areas, as these structures will effectively exclude animals from accessing areas they previously utilized. For example, the Mojave Desert Tortoise (Gopherus agassizzii), who have similar annual site fidelity traits as many snake species, demonstrate a decrease in home range size and an alteration of normal behaviours in response to impenetrable roadside fencing (Peaden et al. 2017). Based on adaptive traits and behaviour, snakes likely are unable to avoid barrier structures, and therefore impermeable snake exclusion fences likely exert negative influences on individual behaviour and home range attributes. There are currently two main management approaches to reducing negative humansnake interactions on the landscape: relocating animals and erecting exclusion/barrier fencing. Efforts to manage and limit snake-human interactions historically have included the relocation of individuals (Brown et al. 2009; Holding et al. 2014). However this strategy often is only a short-term solution (Brown et al. 2009), and comes with potential negative consequences on translocated snakes such as longer movements, increased home range sizes, and a reduction in survival (Nowak et al. 2002; Brown et al. 2009; Wolfe et al. 2018). Fencing structures are intended to redirect movement of animals away from certain areas and/or eliminate access to specific areas (i.e., roads, urban areas), and can be useful for research and inventory purposes (Willson and Gibbons 2010; Hanson and McElroy 2015; Markle et al. 2017). However, although fencing can be very effective at limiting negative impacts, particularly road mortality (Colley et al. 2017), it may come with unintended consequences for individual snakes (Eye et al. 2018) and other herpetofauna species (Ferronato et al. 2014). Globally, and for a suite of different target species, fencing is becoming a very common conservation tool. Yet considering the prevalence of these structures, there has been an obvious lack of research-based initiatives focusing on the effects of these structures on individuals and populations (Jakes et al. 2018). In order to better understand the long-term implications (both positive and negative) of wildlife fences, efforts need to also be directed at 7 understanding the potential shift in the ecology of individuals and populations after fences have been installed. THE WESTERN RATTLESNAKE IN BRITISH COLUMBIA –DEALING WITH DISTURBANCE In British Columbia (B.C.), the Western Rattlesnake (Crotalus oreganus) inhabits areas undergoing particularly rapid development and widespread habitat disturbance. Coincidentally, the rattlesnakes in this province represent the northernmost extent of the clade that extends from northern California through Oregon, Washington and central Idaho. In B.C., suitable rattlesnake habitat is represented by the dry, arid valley bottoms of the province’s southern interior, occurring in the Fraser, Thompson, Nicola, Okanagan, Similkameen and Kettle watersheds within four distinct areas (Southern Interior Reptile and Amphibian Working Group 2016; Figure 1.1). Rattlesnakes are found in the Bunchgrass (BG), Ponderosa Pine (PP), and Interior Douglas-fir (IDF) biogeoclimatic zones of the province (BEC; Biogeoclimatic Zones of British Columbia 2018), and are strongly associated with shrub-steppe, open forest, riparian and rocky ecotypes. Overwintering occurs in communal hibernacula between October–April with known hibernacula ranging between ~400–800m elevation within B.C. (Southern Interior Reptile and Amphibian Working Group 2016). During the early spring and into the summer (April– May), rattlesnakes embark on annual migrations moving away from their hibernacula to their summer foraging habitat (Macartney and Gregory 1988; Gomez et al. 2015; Lomas et al. 2015), with maximum distances from hibernacula recorded up to 4.0 km (Harvey 2015). At the same time, rattlesnake home range attributes, sizes, and movement patterns can vary dramatically across populations, individuals within a hibernacula and available habitat (Lomas 2013; Gomez et al. 2015; Harvey 2015). Rattlesnakes in B.C. show life-history attributes that likely render them particularly susceptible to the impacts of disturbance. The Okanagan Valley represents a large part of the rattlesnake range in B.C., and contains one of the fastest rates of human population growth and development in Canada (Statistics Canada 2014). In some areas of the Okanagan Valley, up to 90 percent of natural vegetation communities have been lost (Lea 2008). Specifically, rapid increases in urbanization, outdoor recreation activities and agriculture (primarily vineyards) are placing increased pressures on these animals. As fallout from the 8 transformation of natural landscapes, recent population estimates suggest a dramatic decrease in rattlesnake densities within the south Okanagan Valley over the last 30 years (Maida et al. 2018). Consequently, Western Rattlesnakes have become federally listed as ‘threatened’ (Species At Risk Public Registry 2018) and provincially the species is considered of ‘special concern’ (Southern Interior Reptile and Amphibian Working Group 2016). Various life history strategies and constraints of this species limits potential recovery from habitat loss and renders them vulnerable to declines including delayed sexual maturity of females (5–8 years; Maida et al. 2018), infrequent parturition (typically every 3 years or longer for females; Macartney and Gregory 1988), and reduced juvenile growth and survival (Charland 1989). Although rattlesnake life histories often include high adult survival rates, novel threats and the dynamics of anthropogenic landscapes (i.e. road mortality or persecution) can quickly change the selection process, potentially resulting in normal strategies becoming dysfunctional and negatively influencing the viability of populations. The high annual fidelity these animals show towards their hibernacula, movement corridors and summer foraging grounds likely further compounds the potential impacts disturbance may have on these rattlesnakes (similar to Eastern Diamondback Rattlesnakes, see: Waldron et al. 2013). 9 Figure 1.1. Western Rattlesnake (Crotalus oreganus) range in the southern interior of British Columbia (B.C.) with inset showing the species’ northern range limits in North America (Produced by Robyn Reudink). Osoyoos, B.C. (indicated by red star) was the location of my study site (see below). B.C. map produced by Rene McKibbin. 10 OSOYOOS INDIAN RESERVE (OIR) AND THE NK’MIP SNAKE PROGRAM The Osoyoos Indian Reserve (OIR) is located in the south Okanagan Valley near Osoyoos, B.C. approximately 4.5 km north of the Canada-USA border (Figure 1.1). Like most of the Okanagan Valley, these lands have been exposed to drastic changes to the natural vegetation communities and landscape (see Figure 1.2). The site encompasses the extreme southern portion of rattlesnake range in B.C (Figure 1.3), and the habitat is primarily comprised of low elevation (300–400 m) arid shrub-steppe habitat characterized by Big Sagebrush (Artemesia tridentata), Antelope Brush (Prushia tridentata) and native grasses. The area is bordered by Osoyoos Lake to the west and by open Ponderosa Pine (Pinus ponderosa) and Douglas-fir (Pseudotsuga menziezii) forests and rocky mountain slopes to the east, which contained rattlesnake communal hibernacula for the study area (~450–600m hibernaculum elevation; Maida et al. 2017; Figure 1.4). To the south, the OIR is bordered by a major highway (Highway 3), extensive agriculture, primarily orchards and vineyards, and housing. At the time of this study, the OIR contained drastic habitat quality contrasts, from heavily developed and fragmented portions, to those in near-pristine condition with minimal human involvement or habitat fragmentation (Figure 1.5). Intensive human activity and development has been primarily situated in the southern areas of the OIR, targeting tourism by including a golf course, winery and vineyard, cultural centre, walking trails, a large campground and condominium resort with associated parking lots and roads. A ~ 4 km snake exclusion fence surrounded most of the resort, parking lots and roadways, as well as the campground in efforts to reduce negative human-snake conflict. The north and eastern portions of the study area had limited human influence and lacked any landscape alterations or disturbance. Mean daily temperatures during this study were similar to historical means over the previous 30 years (1981–2010; Figure 1.6). 11 Figure 1.2. Above: Photograph of the Osoyoos Valley looking north over Osoyoos Lake in 1935 (photograph from the Okanagan Archive Society; www.oldphotos.ca); Below: same area in 2012 (photograph from Emily Lomas) 12 Figure 1.3. Map of the Osoyoos Indian Reserve (OIR) study site (indicated in orange) near Osoyoos, B.C., Canada. The northern and eastern portions remained undeveloped with developed areas concentrated to the southern and western portions of the study site. The condominium community and associated snake exclusion fence is located at the southeast corner of the study site. Inset shows the northern extent of the Western Rattlesnake (Crotalus oreganus) range with study site location (indicated by star). 13 Figure 1.4. Examples of the dominant ecotypes on the Osoyoos Indian Reserve (OIR): (A) open Ponderosa Pine (Pinus ponderosa) and Douglas-Fir (Pseudotsuga menziezii) forests, (B) Shrub-steppe habitat dominated by Big Sagebrush (Artemesia tridentata) and Antelope Brush (Purshia tridentata), (C) and mid-elevation rocky slopes. 14 Figure 1.5. Contrast in habitat quality and landscape disturbance between the southern developed (right) and northern undisturbed (left) portions of the Osoyoos Indian Reserve (OIR) study site. The northern portions containing no habitat barriers and limited human presence. Conversely, disturbance within the southern portions were geared towards tourism and the landscape fragmented by roads, condominiums, vineyards, snake exclusion fencing, golf course, hiking trials and a campground. 15 Figure 1.6. Mean monthly temperatures in Osoyoos, B.C. during the years of this study (2011– 2016) compared to the historical 30-year (1981–2010) mean monthly temperatures (Environmental and Climate Change Canada 2018). The Osoyoos Indian Reserve (OIR) study site was located approximately 2 km north-west of the location where the weather data was recorded. 16 The Nk’Mip Snake Program is a long-term rattlesnake population-monitoring program initiated in 2002 on the OIR, through a partnership with the Osoyoos Indian Band, Nk’Mip Desert Cultural Centre, Environment and Climate Change Canada, and later Thompson Rivers University. As such, the Nk’Mip Snake Program is one of the longest running mark-recapture (since 2002), radio-telemetry (since 2004), and public outreach and education snake programs in western Canada. The overarching objective of the program is to monitor the local snake community as development and human intrusions continue to increase, particularly in the southern areas of the OIR study site (Figure 1.7). The ongoing research has provided information on critical habitat (e.g., hibernaculum, shedding and gestation sites) within the OIR as well as long-term estimates of population size, density and survivorship (Maida et al. 2018). Research also identified that rattlesnakes frequenting the disturbed areas of the OIR demonstrate lower body conditions compared to individuals in the undisturbed portions of the study area (Lomas et al. 2015), and individuals residing in undisturbed habitats have larger home ranges and longer range lengths than animals in minimal to highly disturbed habitats (Lomas et al. in press). Furthermore, ongoing assessment of the effectiveness and implications of management strategies and techniques such as exclusion fencing (Eye et al. 2018) and short-distance translocation (Brown et al. 2009) continue to be explored. In addition to rattlesnakes, the Nk’Mip Snake program has provided inventory information for the conservation of other local snake species in B.C., including the Great Basin Gophersnake (Pituophis catenifer desertocola) and Western Yellow-bellied Racer (Coluber constrictor mormon), thus serving as a reference for recovery and management strategies (Southern Interior Reptile and Amphibian Working Group 2016). The analyses presented in this thesis includes data from my own field seasons on the OIR (2015–2016) as well as radio telemetry and capture data from 2011–2014. 17 Figure 1.7. Timeline showing progression of development and habitat loss in the southern section of the Osoyoos Indian Reserve (OIR) study site from 2004–2016. The campground, golf course and vineyards were established in 2004, and 2004–2016 saw the construction of the Spirit Ridge Resort, Nk’Mip Dersert Cultural Centre, hydro substation, expansion of roads, 4 km long snake fence, expansion of campground and construction of the condominium community. 18 THESIS OBJECTIVES In 2012, construction began on a condominium community located between the east edge of the golf course and the rocky slopes bordering the study site to the east (labeled as “24-unit condominium community” in Figure 1.4; Figure 1.8). To minimize human-snake conflict, an associated permanent, wooden drift fence was built surrounding the ~4 ha community to limit snake access into the community and to re-direct snakes into neighbouring areas. This fence (~0.9 m tall) was approximately 800 m (~400 m circumference) in total length and was composed of galvanized mesh hardware cloth with 0.6 cm openings that was buried approximately 15 cm under the surface (see Eye et al. 2018 for more detail on the entire snake fence structure on the OIR study site), and wooden posts, and top and bottom rails supported the fence structure. Based on previous research and the location of the new development and fence (well within 400 m of some hibernacula), rattlesnakes embarking on spring migration into the valley bottom were expected to encounter the fence early in their annual spring migration. In the following chapter (Chapter 2), I investigate how these barriers and disturbances influenced the spring migration patterns of the snakes, with a specific focus on interactions with the snake fencing structure. I used radio-telemetry, spatial mapping and statistical analysis to quantify rattlesnake spring migration parameters, home range attributes and body characteristics. More specifically, I address the following questions: 1. How does disturbance and varying degrees of land use impact Western Rattlesnake spring migration? 2. How does mitigative snake fencing spatially and behaviourally impact Western Rattlesnake spring migration? 3. How does spring migration influence home range attributes of the snakes during the course of the active season? 4. Is there a relationship between the initial spring migration experience and rattlesnake body characteristics at the end of the active season? The novelty and conservation relevance of this thesis is two-fold. Firstly, spring migration is a notable and important life-history attribute of rattlesnakes in B.C. (see Chapter 2). However, spring migration specifically has never been investigated on this species 19 intensively in B.C. Secondly, snake exclusion fencing is becoming an important and effective strategy to limit negative human-snake encounters, primarily road mortality (Colley et al. 2017; Markle et al. 2017). Although effective, these fencing barriers have recently been documented to have negative impacts on various herpetofauna species (Peaden et al. 2017; Wilson and Topham 2009; Eye et al. 2018). It thus has become of paramount importance to quantify the spatial and behavioural impacts these physical barriers have on individual snakes, and in doing so, assess the long-term consequences of this management technique. Quantifying behavioural shifts amongst individuals intercepted by fencing may lead to uncovering the potential trade-offs of reducing mortality rates while shifting and limiting animal behaviour. In the last chapter (Chapter 3), I briefly summarize the key points of my study, identify important knowledge gaps and make recommendations for the OIR study site and northern snake management and research in general. 20 Figure 1.8. Condominium community and construction surrounded by a wooden, permanent snake exclusion fence. 21 LITERATURE CITED Ballouard, J-M., G. Provost, D. Barré, and Xavier Bonnet. 2012. Influence of a field trip on the attitude of schoolchildren towards unpopular organisms: An experience with snakes. Journal of Herpetology 46: 423–228. Bauder, J.M., H. Akenson, and C.R. Peterson. 2015. Movement patterns of Prairie Rattlesnakes (Crotalus v. viridis) across a mountainous landscape in a designated wilderness area. Journal of Herpetology 49: 377–387. Baxter-Gilbert, J.H., J. Riley, D. Lesbarrères, and J.D. Litzgus. 2015. Mitigating reptile road mortality: Fence failures compromises ecopassage effectiveness. PLoS ONE 10: 15 pp. Beale, C.M., and P. Monaghan. 2004. Behavioural responses to human disturbance: A matter of choice? Animal Behaviour 68: 1056–1069. Beale, M., S. Poulin, C. Ivanyi, and G. Blouin-Demers. 2016. Anthropogenic disturbance affects movement and increases concealment in Western Diamondback Rattlesnakes (Crotalus atrox). Journal of Herpetology 50: 216–221. Beaupre, S.J., and L.E. Douglas. 2009. Snakes as indicators and monitors of ecosystem properties. Pp. 334–368 In Mullin, S.J., and R.A. Siegel (Eds.). Snakes: Ecology and Conservation. Cornell University Press, Ithaca, New York, USA. Béchet, A. J-F Giroux, and G. Gauthier. 2004. The effects of disturbance on beahviour, habitat use and energy of spring staging Snow Geese. Journal of Applied Ecology 41: 689–700. Biogeoclimatic Zones of British Columbia. 2018. Accessed 17 July 2018. https://www.sfu.ca/geog/geog351fall07/Group06/index.html Bonnet, X., G. Naulleau, and R. Shine. 1999. The dangers of leaving home: Dispersal and mortality in snakes. Biological Conservation 89: 39–50. Breininger, D.R., M.R. Bolt, M.L. Legare, J.H. Drese, and E.D. Stolen. 2011. Factors influencing home-range sizes of Eastern Indigo Snakes in central Florida. Journal of Herpetology 45: 484–490. 22 Brown, J.R., C.A. Bishop, and R.J. Brooks. 2009. Effectiveness of short-distance translocations and its effects on Western Rattlesnakes. The Journal of Wildlife Management 73: 419–425. Butler, H., B. Malone, and N. Clemann. 2005. The effects of translocation on the spatial ecology of Tiger Snakes (Notechis scutatus) in a suburban landscape. Wildlife Research 32: 165–171. Charland, M B. 1989. Size and winter survivirship in neonatal Western Rattlesnakes (Crotalus Viridis). Canadian Journal of Zoology 67: 1620–1625. Clark, R.W., W.S. Brown, R. Stechert, and K.R. Zamudio. 2010. Roads, interrupted dispersal, and genetic diversity in Timber Rattlesnakes. Conservation Biology 24: 1059–1069. Colley, M., S.C. Lougheed, K. Otterbein, and J.D. Litzgus. 2017. Mitigation reduces road mortality of a threatened rattlesnake. Wildlife Research 44: 48–59. Christiansen, F., and D. Lusseau. 2015. Linking behavior to vital rates to measure the effects of non-lethal disturbance on wildlife. Conservation Letters 8: 424-431. Czech, K., P.R. Krausman, and P.K. Devers. 2000. Economic associations among causes of species endangerment in the United States. Bio Science 50: 593–601. Devictor, V., R. Julliard, and F. Jiguet. 2008. Distribution of specialist and generalist species along spatial gradients of habitat disturbance and fragmentation. Oikos 117: 507–514. Duchesne, M., S.D. Côte, and C. Barrette. 2000. Responses of Woodland Caribou to winter ecotourism in the Charlevoix Biosphere Reserve, Canada. Biological Conservation 96: 311–317. Durso, A.M., J.D. Wilson, and C.T. Winne. 2011. Needles in haystacks: estimating detection probability and occupancy of rare and cryptic snakes. Biological Conservation 144: 1508–1515. Environment and Climate Change Canada. 2018. National climate data and information archive. Accessed 11 May 2018. http://climate.weatheroffice.gc.ca/climateData/. 23 Evans, G.W., and K.W. Wood. 1980. Assessment of environmental aesthetics in scenic highway corridors. Environment and Behavior 12: 255–273 Eye, D.M., J.R. Maida, O.M. McKibbin, K.W. Larsen, and C.A. Bishop. 2018. Snake mortality and cover board effectiveness along exclusion fencing in British Columbia, Canada. Canadian Field-Naturalist 132:30–35. Ferronato, B.O., J.H. Roe, and A. Georges. 2014. Reptile bycatch in a pest-exclusion fence established for wildlife reintroductions. Journal for Nature Conservation 22: 577–585. Fortin, J.K., K.D. Rode, G.V. Hilderbrand, J. Wilder, S. Farley, C. Jorgensen, and B.G. Marcot. 2016. Impacts of human recreation on Brown Bears (Ursus arctos): A review and management tool. PLoS ONE 11: 26pp. Garrah, E., R.K. Danby, E. Eberhardt, G.M. Cunnington, and S. Mitchell. 2015. Hot spots and hot times: Wildlife road mortality in a regional conservation corridor. Environmental Management 54: 874–889. Gibbons, J.W., D.E. Scott, T.J. Ryan, K.A. Buhlmann, T.D. Tuberville, B.S. Metts, J.L. Greene, T. Mills, S. Poppy, and C.T. Winne. 2000. The global decline of reptiles, déjà vu amphibians. BioScience 50: 653–666. Gill, J.A., K.S. Norris, and W.J. Sutherland. 2001. Why behavioural responses may not reflect the population consequences of human disturbance. Biological Conservation 97: 265–268. Gomez, L.M., K.W. Larsen, and P. Walton. 2004. 'Snake talks' in the classroom: Do they influence children's attitudes? Herpetological Review 34: 338–341. Gomez, L., K.W. Larsen, and P.T. Gregory. 2015. Contrasting patterns of migration and habitat use in neighboring rattlesnake populations. Journal of Herpetology 49: 371– 376. Goode, M. 2010. Population and community responses of reptiles to golf courses. USGA Green Section Record 48: 28–31. 24 Gregory, P.T. 2007. Biology and conservation of a cold-climate snake fauna. Pp, 41–56 In Sebrun, C., and C. Bishop (Eds.). Ecology, conservation and status of reptiles in Canada. Society for the Study of Amphibian and Reptiles, Salt Lake City, Utah, USA. Hanson, K.M., and E.J. McElroy. 2015. Anthropogenic impacts and long-term changes in herpetofauna diversity and community composition on a barrier island in the southeastern United States. Herpetological Conservation and Biology 10: 765–780. Harvey, J.A. 2015. Thermal influences on summer habitat use by Western Rattlesnakes (Crotalus oreganus) in British Columbia. MSc. Thesis. Thompson Rivers University, Kamloops, BC. 94pp. Hartmann, P.A., M. T. Hartmann, and M. Martins. 2011. Snake road mortality in a protected area in the Atlantic Forest of southeastern Brazil. South American Journal of Herpetology 6: 25–42. Holding, M.L., J.A. Frazier, S.W. Dorr, S.N. Henningsen, I.T. Moore and E.N. Taylor. 2014. Physiological and behavioral effects of repeated handling and short-distance translocations on free-ranging Northern Pacific Rattlesnakes (Crotalus oreganus oreganus). Journal of Herpetology 48: 233–239. Huijser, M.P., E.R. Fairbank, W. Camel-Means, J. Graham, V. Watson, P. Basting, and D. Becker. 2016. Effectiveness of short sections of wildlife fencing and crossing structures along highways in reducing wildlife-vehicle collisions and providing safe crossing opportunities for large mammals. Biological Conservation 2016 197: 61–68. Jackson, S.D., D.J. Smith, and K.E. Gunson. 2015. Mitigating road effects on small animals. Pp. 177–207 In Andrews, K.M., P. Nanjappa, and S.P.D. Riley (Eds.). Roads and ecological infrastructure. Johns Hopkins University Press, Baltimore, Maryland, USA. Jakes, A.F., P.F. Jones, L.C. Paige, R.G. Seidler, and M.P. Huijser. 2018. A fence runs through it: A call for greater attention to the influence of fences on wildlife ecosystems. Biological Conservation 227: 310–318. 25 Jenkins, C.L., C.R. Peterson, S.C. Doering, and V.A. Cobb. 2009. Microgeographic variation in reproductive characteristics among Western Rattlesnake (Crotalus oreganus) populations. Copeia 2009: 774–780. Kark, S., A. Iwaniuk, A. Schalimtzek, and E. Banker. 2007. Living in the city: Can anyone become an 'urban exploiter'? Journal of biogeography 34: 638–651. King, M.B, and D. Duvall. 1990. Prairie Rattlesnake seasonal migrations: Episodes of movement, vernal foraging and sex differences. Animal Behaviour 39: 924–935. Kitahara, M., and K. Fujii. 1994. Biodiversity and community structure of temperate butterfly species within a gradient of human disturbance: An analysis on the concept of generalist vs. specialist strategies. Research on Population Ecology 36: 187–199. Kjoss, V.A., and J.A. Litvaitis. 2001. Community structure of snakes in a human-dominated landscape. Biological Conservation 98: 285–292. Kociolek, A.V., A.P. Clevenger, C.C. St. Clair, and D.S. Proppe. 2011. Effects of road networks on bird populations. Conservation Biology 25: 241–249. Landreth, H.F. 1973. Orientation and behavior of the rattlesnake, Crotalus atrox. Copeia 1973: 26–31. Larsen, K.W. 1987. Movements and behavior of migratory garter snakes, Thamnophis sirtalis. Canadian Journal of Zoology 65: 2241–2247. Lawson, P.A., and D.M. Secoy. 1991. The use of solar cues as migratory orientation guides by the Plains Garter Snake, Thamnophis radix. Canadian Journal of Zoology 69: 2700–2702. Lea, T. 2008. Historical (pre-settlement) ecosystems of the Okanagan Valley and lower Similkameen valley of British Columbia - pre-Europe contact to the present. Davidsonia 19: 33p. Lind, A.J., H.H. Welsh, and D.A. Tallmon. 2005. Garter snake population dynamics from a 16-year study: Considerations for ecological monitoring. Ecological Applications 15: 294–303. 26 Lomas, E.V. 2013. Effects of disturbance on the Northern Pacific Rattlesnake (Crotalus oreganus oreganus) in British Columbia. MSc. Thesis. Thompson Rivers University, Kamloops, BC. 117pp. Lomas, E., K.W. Larsen, and C.A. Bishop. 2015. Persistence of Northern Pacific Rattlesnakes masks the impact of human Disturbance on weight and body condition. Animal Conservation 18: 548–556. Lomas, E., C.A. Bishop, J.R. Maida, and K.W. Larsen. In press. Movement ecology of the Northern Pacific Rattlesnake (Crotalus oreganus oreganus) in response to disturbance. Herpetelogica. Lowry, H., A. Lill, and B.B.M. Wong. 2013. Behavioural responses of wildlife to urban environments. Biological Review 88: 537–549. Macartney, M.J., and P.T Gregory. 1988. Reproductive biology of female rattlesnakes (Crotalus viridis) in British Columbia. Copeia 1988: 47–57. Madsen J. 1988. Experimental refuges for migratory waterfowl in Danish wetlands. II. Tests of hunting disturbance effects. Journal of Applied Ecology 35: 398–417 Maida, J.R., D.A. Kirk, O. McKibbin, J.R. Row, K.W. Larsen, C. Stringam, and C.A. Bishop. 2018. Population estimate, survivorship and generation time of the Northern Pacific Rattlesnake (Crotalus o. oreganus) at its northern-most range limit. Herpetological Conservation and Biology 13: 662–672. Maida, J., C. Hooper, K. Larsen, and C. Bishop. 2017. CROTALUS OREGANUS (Western Rattlesnake) Mortality. Herpetological Review 48: 447. Marini, K.L.D., K.A. Otter, S.E. LaZerte, and M.W. Reudink. 2017. Urban environments are associated with earlier clutches and faster nestling feather growth compared to natural habitats. Urban Ecosystems 20: 1291–1300. Markle, C.E., S.D. Gillingwater, R. Levick, and P. Chow-Fraser. 2017. The true cost of partial fencing: Evaluating strategies to reduce reptile road mortality. Wildlife Society Bulletin 41: 342–350. 27 Marsh, D.M., and J.A.G. Jaeger. 2015. Direct effects of roads on small animal populations. Pp. 42–56 In Andrews, K.M., P. Nanjappa, and S.P.D. Riley (Eds.). Roads and ecological infrastructure. Johns Hopkins University Press, Baltimore, Maryland, USA. Martin, A.E., D. Jorgensen, and C.C. Gates. 2017. Costs and benefits of straight versus tortuous migration paths for Prairie Rattlesnakes (Crotalus viridis viridis) in seminatural and human-dominated landscapes. Canadian Journal of Zoology 928: 921– 928. Martin, A.E., and L. Fahrig. 2018. Habitat specialist birds disperse farther and are more migratory than habitat generalist birds. Ecology (online): 9pp. Martínez-Cruz, B., J.A. Goody, and J.J. Negro. 2004. Population genetics after fragmentation: The case of the endangered Spanish imperial eagle (Aquila adalberti). Molecular Ecology 13: 2243–2255. Marzluff, J.M. 2001. Worldwide urbanization and its effects on birds. Pp. 19-47 In Marzluff, J.M., R. Bowman, and R. Donnelly (Eds.). Avian ecology in an urbanizing world. Kluwer Academic Publisher, Norwell, Massachusetts, USA. Matthews, T.J., H.E. Cottee-Jones, and R.J. Whittaker. 2014. Habitat fragmentation and the species-area relationship: A focus on total species richness obscures the impacts oh habitat loss on habitat specialists. Diversity and Distributions 20: 1136–1146. McKinney, M.L. 2002. Urbanization, biodiversity and conservation: The impacts of urbanization on native species are poorly studied, but educating a highly urbanized human population about these impacts can greatly improve species conservation in all ecosystems. BioScience 52: 883–890. Meek, R. 2015. Where do snakes cross roads? Habitat associated road crossings and mortalities in a fragmented landscape in western France. Herpetological Journal 25: 15–19. 28 Morrison, C.D., M.S. Boyce, S.E. Neilsen, and M.M. Bacon. 2014. Habitat selection of a recolonized cougar population in response to seasonal fluctuations of human activitiy. The Journal of Wildlife Management 78:1394–1403. Müllner, A., E.K. Linsenmair, and M. Wikelski. 2004. Exposure to ecotourism reduces survival and affects stress response in hoatzin chicks (Opisthocomus hoazin). Biological Conservation 118: 549–558. Nattrass, S., and D. Lusseau. 2016. Using resilience to predict the effects of disturbance. Scientific Reports 6: 9pp. Nowak, E.M., T. Hare, and J. McNally. 2002. Management of "nuisance" vipers: Effects of translocation on Western Diamond-backed Rattlesnakes (Crotalus atrox). Pp. 533– 560 In Schuett, G.W., M. Höggren, M. Douglas, and H. Greene (Eds.). Biology of Vipers. Eagle Mountain Publishing, Eagle Mountain, Utah, USA. Ohashi, H., M. Saito, R. Horie, H. Tsunnoda, H. Noba, H. Ishii, T. Kuwabara, Y. Hiroshige, S. Koike, Y. Hoshino, H. Toda, and K. Kaji. 2013. Differences in the activity pattern of the Wild Boar, Sus scrofa, related to human disturbance. European Journal of Wildlife Research 59: 167–177. Owen, D.A.S, E.T. Carter, M.L. Holding, K. Islam, and I.T. Moore. 2014. Roads are associated with a blunted stress response in a North American pit viper. General and Comparative Endocrinology 202: 87–92. Parent, C., and P.J. Weatherhead. 2000. Behavioral and life history responses of Eastern Massasauga Rattlesnakes (Sistrurus catenatus catenatus) to human disturbance. Oecologia 125 (2): 170–178. Peaden, J.M., A.J. Nowkowski, T.D. Tuberville, and K.A. Buhlmann. 2017. Effects of roadside fencing on the movements, space, use and carapace length of a threatened tortoise. Biological Conservation 214: 12–22. Pitts, S.L., B.D. Hughes, and I. Mali. 2017. Rattlesnake nuisance removals and urban expansion in Phoenix, Arizona. Western North American Naturalist 77: 309–16. 29 Putman, B.J., and R.W. Clark. 2017. Behavioral thermal tolerances of free-ranging rattlesnakes (Crotalus oreganus) during the summer foraging season. Journal of Thermal Biology 65: 8–15. Reijnen, R., R. Foppen, and G. Veenbaas. 1997. Disturbance by traffic of breeding birds: Evaluating the effect and considerations in planning and managing road corridors. Biodiversity and Conservation 6: 567–581. Robson, L.E., and G. Blouin-Demers. 2013. Eastern Hognose Snakes (Heterodon platirhinos) avoid crossing paved roads, but not unpaved roads. Copeia 2013: 507– 511. Rousseau, P., A. Desrochers, and A.S. Hadley. 2012. Habitat selection and fidelity by Whitethroated Sparrows (Zonotrichia albicollis): Generalist species, specialist individuals? Canada Journal of Zoology 90: 595–601. Saha, A., L. McRae, C.K. Dodd JR., H. Gadsen, K.M. Hare, V. Lukoschek, and M. Böhm. 2018. Tracking global population trends: Population time-series data and a Living Planet Index for reptiles. Journal of Herpetology 52: 259–268. Šálek, M., L. Drahníková, and E. Tkadlec. 2014. Changes in home range sizes and population densities of carnivore species along the natural to urban habitat gradient. Mammal Review 45: 14pp. Sawyer, H., C. LeBeau, and T. Hart. 2012. Mitigating roadway impacts to migratory Mule Deer- a case study with underpasses and continuous fencing. Wildlife Society Bulletin 36: 492–498. Shine, R., M. Lemaster, M.Wall, T. Langkilde, and R. Mason. 2004. Why did the snake cross the road? Effects of roads on movement and location of mates by garter snakes (Thamnophis sirtalis parietalis). Ecology and Society 9: 9pp. Southern Interior Reptile and Amphibian Working Group. 2016. Recovery plan for the Western Rattlesnake (Crotalus oreganus) in British Columbia. Prepared for the B.C. Ministry of Environment, Victoria, B.C. 37pp. 30 Souza, A.M., R.C. Pires, V. Borges, and P.C. Eterovick. 2015. Road mortality of the herpetofauna in a Cerrado ecosystem, central Brazil. Herpetological Journal 25: 141– 148. Species at Risk Public Registry. 2018. Species profile: Western Rattlesnake, Accessed 22 June 2018. http://www.registrelepsararegistry.gc.ca/species/speciesDetails_e.cfm?sid=808 Statistics Canada. 2014. Canada’s population estimates: Subprovincial areas, July 1, 2014. Støen, O-G., A. Ordiz, A.L. Evans, T.G. Laske, J. Kindberg, O. Fröbert, J.E. Swenson, and J.M. Arnemo. 2015. Physiological evidence for human-induced landscape of fear in brown bears (Ursus arctos). Physiology and Behavior 152: 244–248. Switzer, P.V. 1993. Site fidelity in predictable and unpredicatble habitats. Evolutionary Ecology 7: 533–555. Tablado, Z., and L. Jenni. 2017. Determinants of uncertainty in wildlife responses to human disturbance. Biological Reviews 92: 216–233. Vázquez, D.P., and D. Simberloff. 2001. Ecological specialtion and suceptibility to disturbance: Conjectures and refutations. The American Naturalist 159: 606–623. Venter, O., N.N. Brodeur, L. Nemiroff, B. Belland, I.J. Dolinsek, and J.W. Grant. 2006. Threats to endangered species in Canada. Bioscience 56: 903–910. Waldron, J.L., S.M. Welch, S.H. Bennett, W.G. Kalinowsky, and T.A. Mousseau. 2013. Life history constraints contribute to the vulnerability of a declining North American rattlesnake. Biological Conservation 159: 530–538. Ward, R.J., R.A. Giffiths, J.K. Wilkinson, and N. Cornish. 2017. Optimizing monitoring efforts for secretive snakes: A comparison of occupancy and N-mixture models for assessment of population status. Scientific Reports 7: 12pp. Watson, H., M. Bolton, and P. Monaghan. 2014. Out of sight but not out of harm’s way: Human disturbance reduces reproductive success of a cavity-nesting seabird. Biological Conservation 174: 127–133. 31 Wilcove, D.S., D. Rothstein, J. Dubow, A. Phillips, and E. Losos. 1998. Quantifying threats to imperiled species in the United States: Assessing the relative importance of habitat destruction, alien species, pollution, overexploitation, and disease. BiosScience 48: 607–615. Willson, J.D., and J.W. Gibbons. 2009. Drift fences, cover boards and other traps. Amphibian Ecology and Conservation: A handbook of techniques. Pp. 229–245 In Dodd, K, Jr. (Ed.). Amphibian Ecology and Conservation: A handbook of techniques. Oxford University Press, New York, New York, USA. Wilson, J.S., and S. Topham. 2009. The negative effects of barrier fencing on the Desert Tortoise (Gopherus agassizii) and non-target species: is there room for improvement? Contemporary Herpetology 3: 1–4. Wittemyer. G., J.M. Northup, J. Blanc, I. Douglas-Hamilton, R. Omondi, and K.P. Burnham. 2014. Illegal killing for ivory drives global decline in African elephants. Proceedings of the National Academy of Science of the United States of America. 111: 13117– 13121. Wolfe, A.K., P.W. Bateman, and P.A. Fleming. 2017. Does urbanization influence the diet of a large snake? Current Zoology 64: 311–318. Wolfe, A.K., P.A. Flemming, and P.W. Bateman. 2018. Impacts of translocation on a large urban-adapted venomous snake. Wildlife Research 45: 316–324. 32 CHAPTER TWO: IMPACTS OF FENCING AND DISTURBANCE ON WESTERN RATTLESNAKE (CROTALUS OREGANUS) SPRING MOVEMENTS INTRODUCTION How animals distribute themselves on the landscape is regulated by various factors including habitat quality and availability, access to food, predator avoidance, and mate acquisition (DeGregorio et al. 2011; Christiansen et al. 2017). One notable and high-profile example of this process is migration. Migratory behaviour is the seasonal movement by individuals within a population (Wilcove et al. 1998; Naidoo et al. 2016; Martin et al. 2017) and is categorized by having a distinguishable objective with a defined route. However, migration as a biological concept tends to be fluid and at times elusive, complicated by the variation seen in distances travelled, timing (annual to daily), and frequency within individual lifetimes (Hoare 2009). Migration is critical for species survival (Hoare 2009), allowing individuals to respond to seasonal changes in resource availability and exploit habitats beneficial or critical for different life history processes such as breeding, overwintering and grazing or hunting prey (Alerstam et al. 2003). Therefore, interruptions to this phenomenon can restrict individuals from acquiring resources, in turn impacting fitness and ultimately population persistence (Gubili et al. 2016). Understanding the patterns and implications of migration within wildlife populations may be a critical component of maintaining functional and connected landscapes (Epps et al. 2011; Naidoo et al. 2016) crucial for both species conservation and community/landscape preservation. Due to increasing anthropogenic pressures and continued land-use transformations globally, the process of animal migration has been undergoing alterations in many different taxa (Wilcove et al. 1998; Dobson et al. 2010; Naidoo et al. 2016). The specific drivers of these changes vary across species, locale and temporal scales, but the leading identifiable threats are the destruction of habitat, construction of physical barriers (e.g., dams and fences), animal exploitation and climate change (Harris et al. 2009). Disturbed and fragmented areas typically exhibit ecological traits that differ from natural habitat, specifically altering resource composition, abundance and access for animals (Shine et al. 2004). Furthermore, 33 features and situations that present a barrier or restrict animal movement run the risk of fragmenting populations and presenting obstacles for individuals (Gubili et al. 2016). For example, weirs that fragment streams and tributaries prevent migratory trout (Salmo trutta) from reaching their spawning destinations (Gosset et al. 2006), and Reindeer (Rangifer tarandus tarandus) avoid crossing parallel power lines and roads, reducing migration area and access to grazing habitat (Vistnes et al. 2004). Geese (Anser brachyrhynchus) that stage their migration within disturbed, human-influenced farmland habitat display relatively lower growth rates and reduced reproductive success (Madsen 1995). Migrations involving long distances and the mass movement of animal aggregations tend to attract relatively more attention (e.g. Wildebeest- Hopcraft et al. 2014; Pacific salmon- Kovach et al. 2015) compared to smaller scale, individual migrations such as performed by many temperate amphibians and reptiles (Duvall et al. 1990; Chiszar et al. 2014; Yermokhin et al. 2015). However, human disturbance and landscape barriers may have a large and less visible impact on the seasonal movement of these animals despite smaller annual travel distances and the variation of spatial use on the landscape. Smaller migratory species with small-scale home ranges and movement patterns therefore may provide useful systems for investigating the implications of disturbance and barriers on natural migration behaviour. Many temperate snakes, particularly those at northern latitudes, undertake seasonal migrations between overwintering hibernacula to summer foraging and mating grounds (Landreth 1973; Duvall et al. 1990; Jørgensen et al. 2008; Chiszar et al. 2014; Gomez et al. 2015). Such migrations are dangerous and energetically expensive but necessary in cold climates where dependable suitable hibernacula are scarce. These scarce winter shelters are often used by individuals from a wide area and the lack of sufficient summer resources near communal hibernacula obligates these animals to migrate seasonally. These snakes typically exhibit three types of movement and behavioural patterns throughout the course of an active season, namely (1) spring migration away from hibernaculum, (2) mid-summer movements within established hunting and mating grounds, and (3) fall migration back to hibernaculum for overwintering. A number of factors have been posited to drive snake movements away from hibernacula, namely the spatial separation between hibernacula and foraging/mating 34 grounds, limited suitable habitat for hibernacula, reduced intra-competition within a communal hibernacula community (Larsen 1987; Bauder et al. 2015), and the unequal dispersion of appropriate thermal habitat across the landscape (Huey 1991; Shine et al. 2004; Harvey 2015). Spring migration can result in snakes travelling considerable distances from their hibernaculum while typically demonstrating linear pathways (Martino et al. 2012; Chiszar et al. 2014). These long and straight movements likely optimize encounter rates with irregularly distributed resources on the landscape such as prey and mates (Duvall and Schuett 1997; Duvall et al. 1997). Individual Prairie Rattlesnakes (Crotalus viridis viridis) presented with food, during a food supplementation study, had reduced movement frequencies and distance during spring migration (Duvall et al. 1990), suggesting a primary function of spring migration for this species was to locate and successfully hunt prey. However, detailed assessments of spring migration patterns in temperate climate snakes are rare. Migratory effects aside, the impacts of human disturbance, habitat loss and fragmentation on snake behaviour and movements are widely recognized (Parent and Weatherhead 2000; Shine et al. 2004; Breininger et al. 2011; Lomas et al. in press). For example, Eastern Indigo Snakes (Drymarchon couperi) have smaller home ranges in fragmented landscapes in Florida (Breininger et al. 2011), while Western Diamondback Rattlesnakes (Crotalus atrox) and Eastern Massasauga Rattlesnakes (Sistrurus c. catenatus) both move less frequently in disturbed, human-influenced habitats (Parent and Weatherhead 2000; Beale et al. 2016). Furthermore, garter snakes (Thamnophis sirtalis parietalis- Shine et al. 2004), and Eastern Hognose Snakes (Heterodon platirhinos; Robson and Blouin-Demers 2013) appear to exhibit roadway avoidance. Prairie Rattlesnakes (Crotalus v. viridis) show more tortuous (crooked) movements and reduced body conditions in human-dominated landscapes (Martin et al. 2017). Fencing infrastructure on the landscape is a common cause of habitat fragmentation. But, it also is used to reduce direct impacts of disturbance to herpetofauna, including snakes (Colley et al. 2017; Markle et al. 2017). Fencing is used to restrict snake access or movement into a specific area (Gregory 2007), or it can be used to deflect snake movement away from certain areas and/or towards favourable habitat (Willson and Gibbons 2009). These 35 mitigative barriers can be very effective at reducing human-snake encounters in residential areas or roadways (Colley et al. 2017). However, fencing still constitutes a physical barrier obstructing natural movement patterns and behaviour that may have consequences for individuals (i.e., mortality- Ferronato et al. 2014; Eye et al. 2018). The effects on spring migration patterns and movements of snakes by fencing or other anthropogenic barriers have not been adequately assessed. The Western Rattlesnake (Crotalus oreganus) exists at the northern limits of its range in British Columbia (B.C.) (Gomez et al. 2015; Lomas et al. 2015; Maida et al. 2018). Listed as “threatened” in Canada (COSEWIC 2015), the primary threats to Western Rattlesnakes are habitat loss, fragmentation, road mortality and human persecution. Rattlesnakes in B.C. occupy hibernacula between October and April (Macartney and Gregory 1988; Brown et al. 2009; Maida et al. 2017) and can move up to 4 km away from their overwintering habitat during the active season (Harvey 2015). Within B.C., this snake resides in southern, semiarid valley bottoms that are subjected to increasing human development and some of the fastest rates of human population growth in Canada (Statistics Canada 2014; Lomas et al. 2015). Due to strong fidelity to hibernacula, summer foraging grounds and migration corridors (Brown et al. 2009; Gomez et al. 2015) these animals appear to lack the behavioural plasticity to adjust movement patterns and spatial distributions to avoid new development or other land-use changes, similarly to that documented for Eastern Diamondback Rattlesnakes (Crotalus adamanteus- Waldron et al. 2013). Lomas et al (2015; in press) found rattlesnakes frequenting disturbed habitats demonstrated reduced body conditions, smaller home ranges and shorter range lengths at the extreme southern portion of their range in Canada. At that same study site, I investigated how human disturbance and fragmented landscapes altered rattlesnake spring migration patterns. I was specifically interested in shifts of spring migration behaviour by those snakes encountering mitigative fencing. I predicted that migrating rattlesnakes encountering barriers (landscape and fencing) would have shorter spring migration distances, shorter migration duration, and more tortuous movements compared to individuals moving through undisturbed landscapes. Further, I also predicted that alterations to spring migration would 36 impact rattlesnake home range parameters, and have a negative influence on individual body characteristics (growth and body condition) over the course of the active season. METHODOLOGY Study Site This study took place on the Osoyoos Indian Reserve (OIR) near Osoyoos, B.C., Canada (119.4° W, 49.28° N). The 450 ha area mainly was comprised of low elevation (300– 400 m) arid shrub-steppe habitat characterized by Big Sagebrush (Artemesia tridentata), Antelope Brush (Purhsia trdentata) and native grasses (see Brown et al. 2009; Lomas et al. 2015). The site was bordered by Osoyoos Lake to the west and by open Ponderosa Pine (Pinus ponderosa) forests and rocky slopes to the east that contained rattlesnake communal hibernacula at ~500–650 m elevation (Maida et al. 2017). Hibernacula on the OIR study site ranged between ~280–1880 m to human disturbance. On the study site, rattlesnakes demonstrate two main spring movement patterns: (1) move west and downslope into the valley bottom or (2) move east into the higher elevation forest (Lomas 2013). In this study, I focused on the animals moving from their hibernaculum into the valley bottoms in both the southern and northern portions of the study site. See Chapter 1 for a more detailed study site description. Radio Telemetry From 2011–2016 only adult male rattlesnakes were radio-tracked due to the extreme variations in movement patterns exhibited by females of different reproductive status (Macartney and Gregory 1988). Snakes were captured each year during egress (April-May) at or near hibernacula for transmitter implantation. Animals were transported to a nearby veterinary clinic for surgical implantation of radio-transmitters (SB-2; Holohil Sytems Ltd., Carp, Ontario, Canada). Surgical procedures followed Reinert and Cundall (1982) and pharmaceutical procedures followed Brown et al. (2009). Transmitters weighing on average 2.5% of total body weight (range: 1.6–3.7%) were implanted into the coelomic cavity. Rattlesnakes were held in captivity 24–48 hours following surgery to allow recovery and rehydration before being released at their original capture location. Each telemetered snake was located approximately every 2–3 days throughout the 37 entire active season (April–October). In doing so, snakes were tracked from hibernacula in the spring, throughout the summer season and then back to hibernacula in the fall (ingress), unless mortality and/or transmitter failure occurred. A 3-element yagi antenna and a portable radio-telemetry receiver (TRX-1000S; Wildlife Materials Inc., Murphysboro, Illinois) was used to track individuals and the Universal Transvers Mercator (UTM) coordinates of each location using a MobileMapper 6 device (Magellan Professional Inc, Santa Clara, California, USA) was recorded. Each snake was briefly recaptured approximately once a month and inspected to monitor health and ensure the surgery incision was healing appropriately; during these ‘check-ups’ individual snakes were weighed (g) and measured by directing their heads inside plastic tubing to safely record snout-vent length (SVL - cm). Apart from that, care was taken to collect location data of animals from 2–5 m away in order to limit disruption and changes in snake behaviour. Spring Migration Parameters I calculated a series of parameters to quantify rattlesnake spring migration on the study site (Table 2.1). Outbound spring migration was identified as ending when movements away from hibernaculum became < 50 m between consecutive locations and/or four days spanned between movements (Lomas et al. in press). I assessed spring migration distance (MD) as the total straight-line distance between the start of migration (hibernaculum) to the end (start of summer foraging). Spring migration path length (MPL) was the total sum distance of each yearly movement event until the start of the summer foraging season (i.e., shorter, more infrequent movements), calculated from all locational fixes. I calculated spring migration path sinuosity (MPS) as MD/MPL, creating a ratio ranging between 0 and 1, with values approaching 1 indicating straighter movements and values approaching 0 equating to more crooked movements (Saumure et al. 2010; Martin et al. 2017). Distance to nearest source of disturbance (DTD) was the straight-line distance from each individual’s starting point in the spring (hibernacula) to the first confirmed contact with human disturbance. Rattlesnakes from the same hibernacula used dissimilar migration routes and corridors and may have come into first contact with disturbance at different times. Therefore, snakes from the same hibernacula often had different DTD values. Spring migration duration (DAYS) was calculated from the first move away from hibernaculum in the spring until the beginning of summer foraging. I calculated MD, MPL, and DTD using the measuring tool in Garmin 38 Basecamp (version 4.6.2). Home Range Parameters For the entire active season, I estimated rattlesnake home ranges using the 100% minimum convex polygon (MCP). The MCP method creates a polygon around the outermost points plotted on a map, including 100% of the individual’s telemetry locations. I chose the MCP method to maximize comparisons with existing literature on snake activity studies, including many studies at the northern extent of the Western Rattlesnake’s range in B.C. (Brown et al. 2009; Holding et al. 2014; Harvey 2015; Zappalorti et al. 2015; Lomas et al. in press). This method also has been suggested to best reflect herpetofauna home range size (Row and Blouin-Demers 2006; Shipley et al. 2013; MacGowan et al. 2017), compared to kernel estimators which are more appropriate for investigating habitat use (Row and BlouinDemers 2006; Mata-Silva et al. 2018). I calculated the core area of use within each snake’s home range using 50% fixed-kernel (KD) isopeths (Tiebout and Cary 1987; Mata-Silva et al. 2018; Lomas et al. in press). For this calculation, I used an ad hoc method to select the appropriate smoothing factor (hadhoc) to prevent over or under-smoothing (Berger and Gese 2007). I decreased the reference bandwidth (href) by 0.1 until I found the home range estimate that included all telemetry locations and represented the smallest continuous polygon with no lacuna (Berger and Gese 2007; Kie 2013; Bauder et al. 2015), and I then used the 50% isopeths outputs from those estimates. I calculated both 100% MCP and 50% fixed kernel isopeths using the HRT extension (Home Range Tools; Rodgers et al. 2007) in ArcGis 10.2.2. 39 Table 2.1. Spring migration and home range parameter descriptions and overview of methods used to calculate each parameter for Western Rattlesnakes from 2011–2016 on the Osoyoos Indian Reserve (OIR) near Osoyoos, B.C. Parameter Abbr. Description Spring Migration Parameters Spring Migration Distance (m) Spring Migration Path Length (m) Spring Migration Path Sinuosity Distance to Disturbance (m) Spring Migration Duration MD MPL MPS DTD DAYS Straight line distance (m) between hibernacula and beginning of summer foraging behaviour. Total sum distance of each movement between start of migration and beginning of summer foraging behaviour. (MD/MPL) Estimating straightness of path. Values closer to 1 indicated straighter paths, values closer to 0 indicate more crooked paths The total straight line distance from hibernacula to the first confirmed contact with human disturbance Duration of spring migration, calculated from the first move from hibernacula until start of summer foraging behaviour. Home Range Parameters Home Range (100% Maximum Convex Polygon) (ha) MCP Core Area of Use (ha) CA Home range created by a polygon around all the outermost points of all telemetry locations throughout summer active season Calculated using the 50% kernel isopeths of a rattlesnake's home range 40 Body Condition and Growth I calculated body condition index scores for each telemetered rattlesnake by using the residuals from a log-transformed regression between weight and SVL (cf. Parent and Weatherhead 2000; Taylor et al. 2005; Shipley et al. 2013; Lomas et al. 2015). A positive residual value indicated the mass for an individual snake was higher than predicted by its SVL, and vice versa. I considered weight, SVL and body condition data collected near the end of each active season to reflect summer foraging success and health. To this end, I used August measurements because not all snakes were successfully captured later in the active season (September/October). To determine individual growth, I used a modification of Brody’s formula to standardize instantaneous growth over the active season for each telemetered rattlesnake (Brody 1945; Maida et al. 2018): ∆GR= (logeSVL2 – logeSVL1)/((t2-t1)/167) In this equation, loge is the natural base of the logarithm, SVL1 is the SVL at the beginning of the active season (April/May), SVL2 is the SVL at the end of the active season (August), t1 is the date of first capture and t2 represents the date at second (last) capture. Based on previous radio telemetry studies at this study site, the estimated active season for rattlesnakes is approximately 1 April – 15 September (Brown et al. 2009; Lomas et al. in press). Therefore, the constant (167) in the equation represents in days the estimated length of the active season or ‘growth season’ (cf. Maida at el. in press). We assumed that no growth occurred during hibernation (King et al. 2016). Quantifying Migration Groups I assigned individual telemetered rattlesnakes into one of three distinct categories based on their spring migration experience: Individuals in the FENCE category encountered the condominium snake fence during their spring migration, while DSTB snakes encountered other forms of disturbance (i.e. a hiking trail, golf course fairway or road) during migration but did not contact the snake fence. Lastly, UN snakes were individuals in the northern portion of the OIR study site that did not encounter any type of disturbance during migration 41 or throughout the entire active season. For these snakes, I simply used the straight-line distance from the individual hibernaculum to the nearest source of human disturbance to quantify DTD. Prey Population Monitoring Small mammal live-trapping was performed on the study site to monitor the prey base population for rattlesnakes from 2012–2016. Sampling was performed once a year (end of June/early July) over four trapping grids. Fixed grid locations were used each year to encompass the main land use types on the landscape (Figure 2.1). To this end, a single grid was established in the golf course (in the shrub-steppe between fairways), a vineyard, undisturbed shrub-steppe, and the toe slope of the eastern mountains (‘Hillside’- below hibernacula locations). Each grid was sampled using Longworth-style traps (Little Critter Live Traps, Rogers Manufacturing, West Kelowna, BC) set up in a 6 x 6 grid (36 traps in total) spaced 15 m apart, and each trap was covered with a board (~ 15 cm x 30 cm) to protect them from sun and rain (Larsen et al. 2007). Three baiting nights were followed by three trapping nights. Traps were baited with a piece of apple (source of water) along with rolled oats and sunflower seeds and were provided with synthetic bedding for warmth. During trapping nights, traps were baited and armed immediately prior to sunset and then checked the following morning starting at sunrise (~ 0500h). Captured small mammals were identified to species, sexed, ear-tagged (Monel #1; Kent Scientific Corporation, Torrington CT) and released immediately at the point of capture. Due to low capture rates, mark-recapture models could not be fitted so I used counts (i.e. number of unique individuals captured within three days) to represent the small mammal population density within each grid (Torre et al. 2016). I assumed that the uncaptured proportion of the small mammal population is constant, and that individual counts within the grids would yield a representative density estimate (Slade and Blair 2000; Torre et al. 2016). 42 Figure 2.1. Map of the southern, developed portion of the Osoyoos Indian Reserve (OIR) study site near Osoyoos, B.C. Development included a golf course, vineyards, resort, cultural centre, hiking trails, campground, condominium community and construction, and associated roads and parking lots. Snake fencing represented in red. Small mammal trapping gird locations used from 2012–2016 are as follows: (1) golf course, (2) hillside, (3) vineyard, and (4) undisturbed shrub-steppe. 43 Statistical Analysis I used R version 3.4.3 (R Development Team Core 2017) to perform all statistical analysis. I log-transformed migration distance (MD), 100% MCP and 50% KD to meet assumptions of normality. I used univariate one-way ANOVA to compare each migration and home range parameter as well as end-of-year body condition and instantaneous growth between the three migration categories. I used Tukey’s honesty significant difference (HSD) test for post hoc comparisons. I also used linear regressions to investigate the relationship between migration and home range attributes and the distance to nearest source of disturbance. Data on telemetered snakes were used only if the animals were tracked for the entire spring migration (started at or within ~50 m of hibernacula), and for more than 75% of the active season, or they had provided at least 20 locations (Lomas et al. in press). However, due to mortality and/or transmitter failure, not all snakes were tracked for the entire active season. Given strong site fidelity to hibernacula (Maida et al. 2017), for animals not successfully tracked back to their hibernacula I used the location of their spring hibernacula as the predicted endpoint in home range calculations (note that every snake tracked during this study for an entire active season returned to the same hibernacula from which it originated). Considering GPS error was generally ± 5m, I assigned consecutive locations that were <10 m apart as identical, i.e. representing no movement. All statistical interpretation was guided by =0.05 and all means reported include ± 1 standard error (SE). RESULTS A total of 78 adult male rattlesnakes were equipped with radio-transmitters from 2011–2016. Of those, 27 individual rattlesnakes (FENCE n = 7; DSTB n = 10; UN n = 10) were tracked for their entire spring migration (starting at or near hibernaculum) and for more than 75% of the active season, each furnishing more than 20 location points. No snakes from 2013 met this criteria and were thus excluded in my analysis (telemetry started in June that year). Four snakes tracked for their entire spring migration did not provide data through the entire active season due to road mortality (n = 1), predation (n = 1) and transmitter failure/unknown (n = 2). Both spring migration (F2,24= 0.056, P= 0.95) and entire active season (F2,24= 0.30, P= 0.74) monitoring effort (number of telemetry locations) did not differ across the three spring migration categories. During their spring migration, rattlesnakes 44 within the DSTB or FENCE categories encountered either a snake exclusion fence, hiking trail or a golf course fairway. However, snakes in disturbed habitats also encountered roads, vineyards, fencing, fairways and hiking trails throughout the rest of the active summer period. Migration Parameters Monitoring effort (i.e., the number of times telemetry locations were recorded for each snake) did not influence spring migration distance (MD; F1,25= 0.10, P= 0.75), spring migration path length (MPL; F1,25= 2.53, P= 0.12), or spring migration path sinuosity (MPS; F1,25= 2.96, P= 0.09) estimates. Mean MD was 892.6 m (± 92.1 m; FENCE: 𝑥̅ = 421m ±30.8 m; DSTB: 𝑥̅ = 780 m ±80.8 m; UN: 𝑥̅ = 1335 m ± 137.0 m), and there was a significant relationship with distance to disturbance (DTD; F1,25= 58.7, P< 0.01). Log-transformed MD differed significantly between categories (F2,24= 34.78, P< 0.01; Figure 2.2), with snakes in the UN category displaying significantly greater MD lengths than both FENCE and DSTB categories (FENCE: P< 0.01; DSTB: P< 0.01). In addition, snakes in the FENCE category had shorter MD than snakes in the DSTB category (P< 0.01). The mean migration path length (MPL) for all snakes in my study was 1443.8 m (± 87.4 m), and there was a significant relationship between MPL and DTD (F1,25= 14.1, P< 0.01). Rattlesnakes within the UN category (𝑥̅ = 1714.1 m ±148.6 m) has significantly longer MPL than individuals in the FENCE category (𝑥̅ = 1056 m ±142.9 m; P= 0.01). The mean migration path sinuosity index (MPS) for all snakes in my study was 0.60 (± 0.04), and there was a significant correlation between this measurement and DTD (F1,25= 23.1, P< 0.01; Figure 2.3). Furthermore, MPS differed significantly between the three migration categories (F2,24= 13.62, P< 0.01) with the UN category (𝑥̅ = 0.78 ±0.03) MPS being significantly straighter than both the FENCE and DSTB categories (FENCE: 𝑥̅ = 0.43 ±0.04; P< 0.01; DSTB: 𝑥̅ = 0.55 ±0.06; P< 0.01). Mean values indicate the FENCE category migration path was more tortuous than straight, whereas DSTB category sinuosity index indicated slightly more straight than crooked and UN snakes had almost completely straight migration movements. Furthermore, the lowest sinuosity index value for snakes in the UN category was 0.54 (range 0.54–0.89), compared to the highest sinuosity index values registered in the FENCE category (0.53, range 0.28 –0.53). 45 The mean spring migration duration (DAYS) across the study site was 41 days (± 2.4 days), and there was no relationship between DAYS and DTD (F1,25= 0.1, P= 0.48). Moreover, there was no difference in DAYS between the three groups (F2,24= 0.09, P= 0.91), nor was there a relationship in MD (the length of the individual’s migration) and DAYS (F1,25= 0.03, P= 0.87). 46 Figure 2.2. Mean ± SE of log-transformed migration distance (MD) of the three migration categories (top) and relationship between log-transformed migration distance and distance to disturbance (bottom) of male Western Rattlesnakes (Crotalus oreganus) on the Osoyoos Indian Reserve (OIR) near Osoyoos, B.C., Canada. Sample sizes for the categories are FENCE (n= 7), DSTB (n= 10), and UN (n= 10). Means with different letters indicate a significant difference at α= 0.05. Individuals in the FENCE category encountered the condominium snake fence during their spring migration, while DSTB snakes encountered other forms of disturbance. Thirdly, UN snakes were individuals in the northern portion of the study site that did not encounter any type of disturbance during spring migration or throughout entire active season. 47 Figure 2.3. Mean ± SE of migration path sinuosity (MPS) of the three migration categories (top) and the relationship between migration sinuosity and distance to disturbance (bottom) of male Westner Rattlesnakes (Crotalus oreganus) on the Osoyoos Indian Reserve (OIR) near Osoyoos, B.C., Canada. Sample sizes for each group are FENCE (n= 7), DSTB (n=10), and UN (n=10). Means with different letters indicate a significant difference at α= 0.05. See Figure 2.2 caption for description of migration categories. 48 Home Range Parameters Rattlesnake MCP home range sizes varied from 4.4–103 ha (FENCE: 𝑥̅ = 10.1 ±1.4 ha; DSTB: 𝑥̅ = 21.1 ±2.9 ha; UN: 𝑥̅ = 36.8 ±8.2 ha) on the OIR study site, and there was a significant relationship between rattlesnake log-transformed (MCP) home range size and both DTD and MD (DTD: F1,25= 23.3, P< 0.01; MD: F1,25= 51.66, P< 0.01). Logtransformed home range size differed significantly between the three migration categories (F2,24= 9.50, P< 0.01; Figure 2.4) with the FENCE category home range being significantly smaller than both DSTB and UN categories (DSTB: P= 0.02; UN: P< 0.01). Furthermore, rattlesnake home ranges in the FENCE category only extended an average of 106.9 m (±39.4 m; range 0–205 m) beyond the fencing structure (distance away from hibernacula). Mean core area of use was 8.9 ha (± 1.63 ha; FENCE: 𝑥̅ = 2.7 ±0.4 ha; DSTB: 𝑥̅ = 7.4 ±0.8 ha; UN: 𝑥̅ = 13.9 ±10.2 ha) and log-transformed core area was significantly related to DTD (F1,25=23.9, P< 0.01). Log-transformed core area of use differed significantly between the three groups (F2,24=18.3, P< 0.01), where the FENCE category core area was significantly lower than DSTB and UN categories (DSTB: P< 0.01; UN: P< 0.01). Lastly, snakes with larger home ranges (MCP) contained larger core areas (F1,25=137.5, P< 0.01). Body Characteristics The mean initial spring body weight of telemetered snakes was 209.5 g (± 8.7 g), and mean SVL was 66.8 cm (± 1.3 cm); mass was strongly associated with SVL (F1,25= 43.49, P< 0.01). The range of SVL in telemetered snakes was 54.5–78.5 cm, reflecting the bulk of the larger male rattlesnake sizes at the OIR study site (see Figure 2.5). From captures in 2011– 2016 (n= 489), male rattlesnakes larger than 78.5 cm SVL represented 6.7% of total male captures. Furthermore, Maida et al (2018) estimated asymptotic SVL at this site as 73.4 cm SVL for the entire population. Undisturbed snakes (UN: 𝑥̅ = 71.7 ±1.5 cm) were longer than both FENCE (𝑥̅ = 59.8 ±1.8 cm; P< 0.01) and DSTB snakes (𝑥̅ = 65.6 ±1.9 cm; P= 0.04; Figure 2.6), and there was a positive, significant relationship between rattlesnake size (SVL) and DTD (F1,25= 9.7, P< 0.01; i.e. smaller snakes travelled shorter distances to first encounter disturbance than bigger snakes). Due to snake mortality and transmitter failure, size and weight data were not recorded on two snakes within the sample in the month of August resulting in a body characteristics (body condition and growth) sample size of 25 individuals (FENCE n= 6; DSTB n= 10, UN n= 9). 49 Figure 2.4. Mean ± SE of log-transformed minimum convex polygon (MCP) home range estimates of the three migratory categories (top) and the relationship between logtransformed MCP and distance to disturbance (bottom) of male Western Rattlesnakes (Crotalus oreganus) on the Osoyoos Indian Reserve (OIR) near Osoyoos, B.C., Canada. Sample sizes for the migration categories are FENCE (n= 7), DSTB (n= 10), and UN (n= 10). Means with different letters indicate a significant difference at α= 0.05. See Figure 2.2 caption for a description of migration categories. 50 Figure 2.5. Size class distribution of male Western Rattlesnake (Crotalus oreganus) captures (n= 489) on the Osoyoos Indian Reserve (OIR) near Osoyoos, B.C., Canada from 2011– 2016. 51 Figure 2.6. Mean ± SE of spring snout-vent-length (SVL) between the three migrations categories (top) and the relationship between spring SVL and distance to disturbance (bottom) of male Western Rattlesnakes (Crotalus oreganus) on the Osoyoos Indian Reserve (OIR) near Osoyoos, B.C., Canada. Samples sizes of each category are FENCE (n= 7), DSTB (n= 10), and UN (n= 10). Means with different letters indicate a significant difference at α= 0.05. See Figure 2.2 for a description of the migration categories. 52 There was no significant relationship between DTD and end-of-year body condition (F1,23= 0.58, P= 0.46), and body condition did not differ between the three migration categories (F2,22= 0.19, P= 0.83). A post hoc analysis showed the observed power for this test was 11% (f= 0.18, Power (1-β)= 0.11). Given the variances in the samples, to achieve reasonable power (>0.80) in the analysis, the total sample size between the three migration groups would require 303 individuals. Rattlesnakes with higher end-of-year body condition had both larger home ranges and core areas of use (Home range: F1,23= 7.5, P= 0.01; Core area: F1,23= 5.0, P= 0.04). Not surprisingly, smaller snakes grew at a faster rate than larger snakes (F1,23= 5.5, P= 0.02) but instantaneous growth did not appear to differ between the three migration categories (F2,22= 1.17, P= 0.33; Figure 2.7). Post hoc power for this analysis was 35% (f= 0.39, Power (1-β)= 0.35), and to obtain a power estimate > 0.80 would require 69 individuals in the analysis, compared to my own sample size of 25 individuals. Furthermore, there was no significant relationship between instantaneous growth and DTD (F1,23=0.86, P= 0.36). Overall, a Pearson’s correlation matrix showed no significant correlations between each spring migration parameter (MD, MPL, DTD, DAYS and MPS) with either rattlesnake end of year body condition or instantaneous growth (Table 2.2). Prey Population Monitoring Over the 2,160 trap nights conducted during 2012–2016 only two species of small mammals were captured: Great Basin Pocket Mice (Perognathus parvus) were the prominent species, with 179 animals caught (97.2% of total), compared to five individual Deer Mice (Peromyscus maniculatus). This equated to a mean annual population density of 16.4 animals/ha (range: 10.7–26.3 animals/ha; Table 2.3) across the study area. Furthermore, small mammal density appeared homogenous across the landscape (Table 2.3), with no significant differences in capture numbers occurring between the grids (F3,16=1.03, P= 0.41). 53 Figure 2.7. Mean ± SE instantaneous growth rate between the three migration categories (top) and the relationship between instantaneous growth rate and distance to disturbance (bottom) of male Western Rattlesnakes (Crotalus oreganus) on the Osoyoos Indian Reserve (OIR) near Osoyoos, B.C., Canada. Sample sizes for each category are FENCE (n= 6), DSTB (n=10), and UN (n=9). Means with different letters indicate a significant difference at α= 0.05. See Figure 2.2 for a description of the migration categories. 54 Table 2.2. P-values for Pearson’s correlation tests between male Western Rattlesnake (Crotalus oreganus) end of year body condition and instantaneous growth and each migration parameter used (MD, MPL, DTD, DAYS, and MPS). Migration Parameter End of Year Body Condition Instantaneous Growth Spring Migration Distance (MD) P= 0.20 P= 0.58 Total Spring Migration Path (MPL) P= 0.14 P= 0.85 Distance to Disturbance (DTD) P= 0.46 P= 0.77 Spring Migration Duration (DAYS) P= 0.23 P= 0.59 Spring Migration Sinuosity (MPS) P= 0.49 P= 0.44 55 Table 2.3. Small mammal population density estimates (animals/ha) of each trapping grid on the Osoyoos Indian Reserve (OIR) near Osoyoos, B.C., Canada from 2012–2016. Small mammals represented within the population were Great Basin Pocket Mice (Perognathus parvus) and Deer Mice (Peromyscus maniculatus). Year Hillside Golf Course Vineyard Shrub-Steppe Mean 2012 10.7 3.6 25.0 3.6 10.7 2013 28.6 0.0 5.4 8.9 10.7 2014 23.2 35.7 19.6 26.8 26.3 2015 30.4 8.9 21.4 12.5 18.3 2016 Mean 14.3 21.4 7.1 11.1 23.2 18.9 19.6 14.3 16.1 16.4 56 DISCUSSION This is the first detailed study investigating rattlesnake spring migration in relation to human disturbance. More specifically, this is the first study addressing how fence barriers affect rattlesnake spring migratory behaviour and the consequences for summer home ranges and body characteristics. I observed a strong correlation between distance travelled to disturbance with all of the spring migration metrics calculated, home range estimates and animal body size (SVL). I also found that disturbed and fragmented habitats had a strong, consistent influence on rattlesnake spring migration distance and migration sinuosity compared to snakes migrating through undisturbed landscapes during the spring. Furthermore, encountering the fence structure during spring migration had a significantly negative influence on rattlesnake migration distance as well as home range and core area size, compared to snakes in disturbed areas that did not encounter the fencing structure. Physical construction was taking place on site during all/part of the study, and that may have had potential proximate impacts on snakes moving through adjacent habitat or along the fencing structure. Tremors created from the use of heavy machinery may have negative impacts on animals, such as snakes, that use vibrations as a source of communication and sensory input (Hartline 1971; Hill 2001; Lovich and Ennen 2017). However, I did not have the equipment to measure and quantify vibration intensity from the construction zone. In contrast, variations in spring migration distance and path length did not correlate with migration durations within the three categories; in other words, individuals migrating shorter distances in disturbed areas took the same amount of time to complete their spring movements as individuals migrating further distances in undisturbed habitats. In addition, there did not appear to be a clear, consistent relationship between the animals’ body characteristics and the spring migration metrics considered in this study. On the study site, there were two general types of barriers on the landscape interrupting rattlesnake spring migration: physical (fence) and permeable ‘landscape’ barriers (golf course and hiking trails), and both appeared to influence rattlesnake spring migration compared to snakes travelling through undisturbed landscapes. Snakes and other reptiles have been previously observed avoiding or altering their behaviour along open areas including hiking trails and golf course fairways (Parent and Weatherhead 2000; Andrews and 57 Gibbons 2005; Goode 2010;). In Arizona, Gila Monsters (Heloderma suspectum) tend to avoid golf course fairways, while Tiger Rattlesnake (Crotalus tigris) home range sizes and shapes alter in golf course habitats (Goode 2010). Furthermore, Massassauga Rattlesnakes (Sisturus c. catenatus) and Western Diamondback Rattlesnakes (Crotalus atrox) exhibit a reduction in movement frequency within areas of high human use (Parent and Weatherhead 2002; Beale et al. 2016) including hiking trails. The scale of my analysis was insufficient to determine if rattlesnakes avoided golf course fairways and hiking trails on the OIR study site. However, reductions in spring migration distances, total migration path lengths and lower path sinuosity (more crooked movements) in these ‘habitats’ suggest the animals altered their spring migration behaviour compared to individuals in areas lacking human involvement and landscape disturbance. Shorter, more tortuous movements may be a behavioural response to perceived threats or predator avoidance caused by human influence and landscape alterations, and may help decrease the probability of encountering predators and/or humans (Martin et al. 2017). Encountering the snake fence appeared to have a large influence on rattlesnake spring migration, specifically migration distance, migration path length and migration sinuosity compared to individuals migrating in undisturbed habitats. Furthermore, individuals encountering the fence had greatly reduced spring migration distances compared to all other snakes in our study (including both disturbed and undisturbed). Previous studies (see above) show avoidance of open areas associated with landscape disturbance; however, individuals still have the ability to cross these features. Animals encountering barriers such as a golf course fairway or hiking trail face trade-offs between the costs of crossing these impedances versus moving around them (Beyer et al. 2016). On the other hand, individuals encountering a fence do not have this option and must either halt movement or circumnavigate the barrier. Having alternative movement options when encountering an impediment (golf course and/or hiking trails) or a non-permeable barrier (fence) appears to strongly influence the total distance a rattlesnake travels during its spring migration on the OIR study site. Furthermore, rattlesnakes exhibit a high level of fidelity to various areas/habitats throughout their life history (hibernacula, corridors, summer foraging areas). Specifically, adult rattlesnakes seem to contain high spatial and site fidelity on the landscape and are not likely to redistribute themselves on the landscape (Waldron et al. 2013). A novel obstruction, such as a fencing 58 structure, likely impedes access to areas used in previous years. Again, in response to this, individuals may or may not attempt to skirt the barrier (altering path sinuosity and likely distance travelled). Overall, home range patterns and sizes I observed appear to be similar to that previously reported from rattlesnakes at the OIR study site and in other areas at the northern extent of the species’ range (B.C.). Harvey (2015) found MCP home ranges in undisturbed habitats throughout B.C. ranged from 1.5–184.7 ha (𝑥̅ = 52 ha). Furthermore, at the OIR study site, Brown et al. (2009) calculated average rattlesnake MCP to be 25.1 ha and Lomas et al. (in press) calculated mean MCP home ranges at 20.0 ha. Previously, Macartney (1985) determined that den population home ranges (estimated size of area used by individuals in a given den) averaged 26.0–122.3 ha in a study area approximately 140 km north of the OIR study site. Although home range estimates among previous studies, as well as results presented herein, are comparable, it is also indicative of the range of home range sizes across individuals and populations. This is likely determined by a variety of factors including habitat, level of disturbance on the landscape and potentially latitudinal position. Using a disturbance rating (DR) system (Parent and Weatherhead 2002), Lomas et al. (in press) found that rattlesnakes in the undisturbed areas of the OIR study site had larger home ranges than animals within minimal to highly-disturbed areas. I found similar trends, however, interactions with the fencing structure during spring migration had a drastic influence on rattlesnake summer home range size and core areas of use compared to both undisturbed and disturbed snakes. Regardless of disturbance, my results further indicate that home range size was strongly associated with the length of an individual’s spring migration distance, and with respect to disturbance, rattlesnake home range size decreased when in closer proximity to disturbance (shorter DTD = shorter MD). Similarly, Mojave Desert Tortoise (Gopherus agassizii) home ranges decrease with proximity to roads and mitigative road-side fencing (Peaden et al. 2017). The authors inferred this may be a result of road avoidance or increased resources along road sides (water run-off or forbs). Based on anecdotal observations of vegetation and resources adjacent to the construction site and snake fence, this does not appear to be the case on the OIR study site, and changes to home range size likely result from the interaction and behavioural shifts in migration of individuals in 59 response to landscape disturbance and fencing. Growth and body condition are the two most common metrics used to quantify rattlesnake physiology, health and individual/population success (Moore et al. 2000; Jenkins et al. 2009; Lomas et al. 2015; King et al. 2016). One of the primary functions of spring migration in communally-denning rattlesnakes appears to be finding habitats with food (Duvall et al. 1990). Prairie Rattlesnakes (Crotalus v. viridis) with straighter, longer range lengths had higher end of year body condition (Martin et al. 2017), while increased body condition in male Timber Rattlesnakes (Crotalus horridus) allowed individuals to allocate more time to finding mates than hunting (Lind and Beaupre 2015). Furthermore, in B.C., Harvey (2015) observed rattlesnakes in forested habitats moved farther distances and contained higher body conditions than rattlesnakes within the valley bottoms. In support of this, I found individuals with higher body conditions at the end of the active season had maintained larger home ranges and larger core areas. Being a sit-and-wait ambush predator, larger home ranges and core areas, regardless of habitat quality, likely reflect more time for rattlesnakes to search for mates during the summer months, rather than hunting. However, I did not notice any significant trends or differences in end-of-year body condition among the three migration categories. The lack of significant difference and the large sample size required (303 individuals) to achieve statistical power suggests body condition index scores may not be a feasible physiological metric for gauging the effects of different migratory tactics, at least at this specific study area. My observations of rattlesnakes being smaller (SVL) when travelling shorter distances to disturbance mirrors previous observations by Lomas et al. (2015) at the OIR study site; similar trends also were noticed in a Western Rattlesnake population south of the OIR study site in Idaho (rattlesnakes living in disturbed areas were smaller: Jenkins et al. 2009). Furthermore, this trend has been observed in other snake species, such as urban-living Dugites (Pseudonaja affinis) being smaller than non-urban individuals in Australia (Wolfe et al. 2017). Typical snake growth rates involve a highly correlated, linear, negative relationship between size and growth (i.e. smaller snakes grow faster than larger snakes; cf. Macartney et al. 1990; King et al. 2016; Dreslik et al. 2017; Maida et al. 2018). With mean SVL of snakes in disturbed areas (both fence and disturbed categories) being significantly smaller than 60 undisturbed snakes and previously estimated population asymptotic SVL (Maida et al. 2018), we would expect to find an increased rate of growth in snakes in the two disturbed migration categories compared to snakes in the undisturbed category, but this appeared not to be the case in this study. However, due to the lack of power in the analysis, further investigation the impacts on growth (and body condition) based on rattlesnake spring movements and behaviour may be required. Previously, an analysis of road-killed rattlesnakes from the broader region of the south Okanagan showed a relatively large suite of mammalian prey species (McAllister et al. 2016). In contrast, the OIR study site appears to support only two microtine mammal species, Deer Mice and Great Basin Pocket Mice, and in relatively low densities. The homogenous array of low densities and a low species diversity of mammalian prey across the OIR study site suggests the ability for rattlesnakes to successfully find prey may be of particular significance. However, the similarity in prey density estimates between grids in natural landscapes and disturbed areas (i.e., golf course, vineyard) suggests that prey abundances does not appear to be a likely contributor to changes in home ranges and movements that I documented above. Overall, spring migration and behaviour of these animals appear to be highly influenced by human disturbance including interactions with associated exclusion fencing. Similar to the effects of disturbance on rattlesnake body condition outlined by Lomas et al (2015), these effects are not readily apparent when rattlesnakes remain persistent on the landscape. Furthermore, this study does not measure the long-term implications of these impacts. Snake exclusion fencing is becoming a very important conservation tool to reduce the direct impacts of road mortality and other negative consequences of human interactions with snakes (and a variety of species and taxa). However, my findings show that fencing also needs to be considered as a disturbance in itself. Fencing can be effective at reducing direct impacts (i.e., road mortality) of human interactions on snakes. Yet, these structures also pose influences on animals that have been overlooked to date. 61 LITERATURE CITED Alerstam, T., A. Hedenström, and S. Åkesson. 2003. Long-distance migration: Evolution and determinants. Oikos 103: 247–260. Andrews, K.M., and J.W. Gibbons. 2005. Ho do highways influence snake movement? Behavioral responses to roads and vehicles. Copeia 2005: 772 –782. Bauder, J.M., H. Akenson, and C.R. Peterson. 2015. Movement patterns of Prairie Rattlesnakes (Crotalus v. viridis) across a mountainous landscape in a designated wilderness area. Journal of Herpetology 49: 377–387. Beale, M., S. Poulin, C. Ivanyi, and G. Blouin-Demers. 2016. Antorpogenic disturbance affects movement and increases concealment in western diamondback rattlesnakes (Crotalus atrox). Journal of Herpetology 50: 216–221. Berger, K.M., and E.M. Gese. 2007. Does interference competition with wolves limit the distribution and abundance of coyotes? Journal of Animmal Ecology 67: 1075–1085. Beyer, H.L., E. Gurarie, L. Börger, M. Panzacchi, M. Basille, I. Herfindal, B.V. Moorter, S.R. Lele, and J. Matthiopoulos. 2016. 'You shall not pass!': Quantifying barrier permeability and proximity avoidance by animals. Journal of Animal Ecology 85: 43–53. Breininger, D.R., M.R. Bolt, M.L. Legare, J.H. Drese, and E.D. Stolen. 2011. Factors influencing home-range sizes of Eastern Indigo Snakes in central Florida. Journal of Herpetology 45: 484–490. Brody, S. 1945. Bioenergetics and growth with specific reference to the energetic efficiency complex in domestic animals. New York: Reinhold Publishing Corporation. Brown, J.R., C.A. Bishop, and R.J. Brooks. 2009. Effectiveness of short-distance translocations and its effects on Western Rattlesnakes. Journal of Wildlife Management 73: 419–425. Chiszar, D., B.K. Shipley, H.M. Smith, K. Fitgerald, and A.J. Saviola. 2014. Straightness-ofpath during and after vernal migration in Prairie Rattlesnakes, Crotalus viridis, in eastern Colorado. Herpetology Notes 7: 425–436. Christiansen, F., N. Esteban, J.A. Mortimer, A.M. Dujon, and G.C. Hays. 2017. Diel and 62 seasonal patterns in activity and home range size of Green Turtles on their foraging ground revealed by extended Fastloc-GPS tracking. Marine Biology 164: 10p. Colley, M., S.C. Lougheed, K. Otterbein, and J.D. Litzgus. 2017. Mitigation reduces road mortality of a threatened rattlesnake. Wildlife Research 44: 48–59. COSEWIC. 2015. COSEWIC assessment and status report on the western rattlesnake Crotalus oreganus in Canada. Committee on the Status of Endangered Wildlife in Canada. Ottawa. xi + 44pp. (www.registrelep-sararegistry.gc.ca/default_e.cfm). DeGregorio, B.A, J.V. Manning, N. Bieser, and B.A. Kingsbury. 2011. The spatial ecology of the Eastern Massasauga (Sistrurus c. catenatus) in northern Michigan. Herpetologica 67: 71–79. Dobson, A.P., M. Borner, and A.R.E. Sinclair 2010. Road will ruin Serengeti. Nature 467: 272– 273. Dreslik, M.J., D.B. Shepard, S.J. Baker, B.C. Jellen, C.A. and Phillips. 2017. Body size, growth and sexual dimorphism in the Eastern Massasauga (Sistrusus catenatus). Pp 65– 77 In Dreslik, M.J., W.K. Hayes, S.J. Beaupre, and S.P. Mackessy (Eds.). The Biology of Rattlesnakes II. ECO Herpetological Publishing and Distribution. Rodeo, New Mexico, USA. Duvall, D., M.J. Goode, W.K. Hayes, J.K. Leonhardt, and D.G. Brown. 1990. Prairie Rattlesnake vernal migration: Field experimental analyses and survival value. National Geogrgraphic Research 64: 457–469. Duvall, D., D. Chiszar, R.A. Mintzer, and N.O. Roster. 1997. Experimental simulation in behavioral ecology: a multimedia approach with the spatial searching simulation RattleSnake ©. Experimental Biology Online 2: 11pp. Duvall, D., and G.W. Schuett. 1997. Straight-line movement and competitve mate seraching in Prarie Rattlesnakes, Crotalus viridis viridis. Animal Behaviour 54: 329–334. Epps, C.W., B.M. Mutayoba, L. Gwin, and J.S. Brashares. 2011. An empeirical evaluation of the African Elephant as a focal species for connectivity planning in east Africa. Diversisty and Distribrutions 17: 603–612. 63 Eye, D.M., J.R. Maida, O.M. McKibbin, K.W. Larsen, and C.A. Bishop. 2018. Snake mortality and cover board effectiveness along exclusion fencing in British Columbia, Canada. Canadian Field Naturalist 132: 30–35 Ferronato, B.O., J.H. Roe, and A. Georges. 2014. Reptile bycatch in a pest-exclusion fence established for wildlife reintroductions. Journal for Nature Conservation 22: 577–585. Gomez, L., K.W. Larsen, and P.T. Gregory. 2015. Contrasting patterns of migration and habitat use in neighboring rattlesnake populations. Journal of Herpetology 49: 371–376. Goode, M. 2010. Population and community responses of reptiles to golf courses. USGA Green Section Record 48: 28–31. Gosset, C., J. Rives, and J. Labonne. 2006. Effect of habitat fragmentation on spawning migration of brown trout (Salmo trutta L.). Ecology of Freshwater Fish 15: 247–254. Gregory, P.T. 2007. Biology and conservation of a cold-climate snake fauna. Pp. 41–56 In Sebrun, C., and C. Bishop (Eds.). Ecology, conservation and status of reptiles in Canada. Society for the Study of Amphibian and Reptiles, Salt Lake City, Utah, USA. Gubili, C., S. Mriani, B.V. Weckworth, P. Galpern, A.D. McDevitt, M. Hebblewhite, B. Nickel, and M. Musiani. 2016. Environmental and anthropogenic drivers of connectivity patterns: A basis for prioritizing conservation efforts for threatened populations. Evolutionary Applications 10: 199–211. Harris, G., S. Thirgood, J.G.C. Hopcraft, J.P.G.M. Cromsigt, and J. Berger. 2009. Global declines in aggregated migrations of large mammals. Endangered Species Research. 7: 55–76. Hartline, P.H. 1971. Physiological basis for detection of sound and vibration in snakes. Journal of Experimental Biolgy 54: 349–371. Harvey, J.A. 2015. Thermal influences on summer habitat use by Western Rattlesnakes (Crotalus oreganus) in British Columbia. MSc Thesis, Thompson Rivers University, Kamloops, British Columbia. 94 pp. 64 Hill, P.S.M. 2001. Viration and animal communication: a review. Integrative and Comparative Biology 41: 1135-1142. Hoare, B. 2009. Animal migration: Remarkable journeys in the wild. University of California Press. Holding, M.L., J.A. Frazier, S.W. Dorr, S.N. Sloanne, I.T. Moore, and E.N. Taylor. 2014. Physiological and behavioral effects of repeated handling and short-distance translocations on free-ranging Northern Pacific Rattlesnakes (Crotalus oreganus oreganus). Journal of Herpetology 48: 233–239. Hopcracft, J.G.C., J.M. Morales, H.L. Beyer, M. Borner, E. Mwangomo, A.R.E. Sinclair, H. Olff, and D.T. Haydon. 2014. Competition, predation, and migration: Individual choice patterns of Serengeti migrants captured by heirarchical models. Ecological Monographs 84: 355–372. Huey, R.B. 1991. Physiological consequences of habitat selection. The American Naturalist 137: 91–115. Jenkins, C.L., C.R. Peterson, S.C. Doering, and V.A. Cobb. 2009. Microgeographic variation in reproductive characterisitics among Western Rattlesnake (Crotalus oreganus) populations. Copeia 2009: 774–780. Jørgensen, D., Guyer, C., Smith, L.L., and Schuett, G.W., 2008. Movements, migrations, and mechanisms: a review of radiotelemetry studies of Prairie (Crotalus v. viridis) and Western (Crotalus oreganus) Rattlesnakes. Pp 303–316 In Hayes, W.K, Beaman, K.R., Cardwell, M.D. and Bush, S.P. (Eds.). pp 303 –316. The Biology of Rattlesnakes. Loma Linda University Press. Linda Loma, California, USA. Kie, J.G. 2013. A rule-based ad hoc method for selecting a bandwidth in kernel home-range analysis. Animal Biotelemetry 1: 13p King, R.B., K.M. Stanford, P.C. Jones, and K. Bekker. 2016. Size matters: Individual variation in ectotherm growth and asymptotic size. PLoS One 11: e0146299. 65 Kovach, R.P., S.C. Ellison, S. Pyare, and S.A. Tallmon. 2015. Temporal patterns in adult salmon migration timing across southeast Alaska. Global Change Biology 21: 1821– 1833. Landreth, H.F. 1973. Orientation and behavior of the rattlesnakes, Crotalus atrox. Copeia 1973: 26–31. Larsen, K W. 1987. Movements and behavior of migratory garter snakes, Thamnophis sirtalis. Canadian Journal of Zoolology 65: 2241–2247. Larsen, K.W., I.T. Adams, and D.L. Haughland. 2007. Small mammal communities in a pyrogenic habitat mosaic. International Journal of Wildland Fire 16: 728–740. Lind, C.M., and S.J. Beaupre. 2015. Male snakes allocate time and energy according to individual energetic status: Body condition, steroid hormones, and reproductive beahvior in Timber Rattlesnakes, Crotalus horridus. Physiological and Biochemical Zoology 88: 624–633. Lomas, E.V. 2013. Effects of disturbance on the Northern Pacific Rattlesnake (Crotalus oreganus oreganus) in British Columbia. MSc Thesis, Thompson Rivers University, Kamloops, British Columbia. 117 pp. Lomas, E., K.W. Larsen, and C.A. Bishop. 2015. Persistence of Northern Pacific Rattlesnakes masks the impact of human disturbance on weight and body condition. Animal Conservation 18: 548–556. Lomas, E., C.A. Bishop, J.R. Maida, and K.W. Larsen. In press. Movement ecology of the Northern Pacific Rattlesnake (Crotalus oreganus oreganus) in response to disturbance. Herpetologica. Lovich, J.E., and J.R. Ennen. 2017. Reptiles and amphibians. Pp. 97–118 In Perrow, R.M. (Ed). Wildlife and Wind Farms, Conflicts and Solutions. Pelagic Publishing. Exeter, UK. Macartney, J.M. 1985. The ecology of the Northern Pacific Rattlesnake, Crotalus viridis oreganus, in British Columbia. MSc Thesis, Thompson Rivers University, Victoria, British Columbia. 289 pp. 66 Macartney, J.M., and P.T. Gregory. 1988. Reproductive biology of female rattlesnakes (Crotalus viridis) in British Columbia. Copeia 1988: 47–57. Macartney, J.M., P.T. Gregory, and M.B. Charland. 1990. Growth and sexual maturity of the Western Rattlesnake, Crotalus viridis, in British Columbia. Copeia 1990: 528–542. MacGowan, B.J., A.F.T. Currylow, and J.E. MacNeil. 2017. Short-term responses of Timber Rattlesnakes (Crotalus horridus) to even-aged timber harvests in Indiana. Forest Ecology and Management 387: 30–36. Madsen, J. 1995. Impacts of disturbance on migratory waterfowl. International Journal of Avian Science 137: S67–S74. Maida, J.R., D.A. Kirk, O. McKibbin, J.R. Row, K.W. Larsen, C. Stringam, and C.A. Bishop. 2018. Population estimate, survivorship and generation time of the Northern Pacific Rattlesnake (Crotalus o. oreganus) at its northern-most range. Herpetological Conservation and Biology 13: 662–672. Maida, J., K. Larsen, C. Hooper, and C.A. Bishop. 2017. CROTALUS OREGANUS (Western Rattlesnake) Mortality. Herpetological Review 48: 487. Markle, C.E., S.D. Gillingwater, R. Levick, and P. Chow-Fraser. 2017. The true cost of partial fencing: Evaluating strategies to reduce reptile road mortality. Wildlife Society Bulletin 41: 342–350. Martin, A.E., D. Jørgensen, and C.C. Gates. 2017. Costs and benefits of straight versus tortuous migration paths for Prairie Rattlesnakes (Crotalus viridis viridis) in seminatural and human-dominated landscapes. Canadian Journal of Zoology 95: 921–928. Martino, J.A., R.G. Poulin, D.L. Parker, and C.M. Somers. 2012. Habitat selection by grassland snakes at northern range limits: Implications for conservation. Journal of Wildlife Management 76: 759–767. Mata-Silva, V., D.L. DeSantis, A.E. Wagler, and J.D. Johnson. 2018. Spatial ecology of Rock Rattlesnakes (Crotalus lepidus) in far west Texas. Herpetologica 74: 245–254. McAllister, J.M., J.R. Maida, O. Dyer, and K.W. Larsen. 2016. Diet of roadkilled Western 67 Rattlesnakes (Crotalus oreganus) and Gophersnakes (Pituophis catenifer) in southern British Columbia. Northwestern Naturalist 97: 181–189. Moore, I.T., J.P. Lerner, D.T. Lerner, and R.T. Mason. 2000. Relationships between annual cycles of testosterone, corticosterone, and body condition in male Red-Spotted Garter Snakes, Thamnophis sirtalis concinnus. Physiological and Biochemical Zoology 73: 307–312. Naidoo, R., M. J. Chase, P. Beytell, P. Du Preez, K. Landen, G. Stuart-Hill, and R. Taylor. 2016. A newly discovered wildlife migration in Namibia and Botswana is the longest in Africa. Oryx 50: 138–146. Peaden, J.M., A.J. Nowakowski, T.D. Tuberville, K.A. Buhlmann, and B.D. Todd. 2017. Effects of roads and roadside fencing on movements, space use, and carapace temperatures of a threatened tortoise. Biological Conservation 214: 13–22. Parent, C., and P.J. Weatherhead. 2000. Behavioral and life history responses of Eastern Massasauga Rattlesnakes (Sistrurus catenatus catenatus) to human disturbance. Oecologia 125: 170–178. R Development Tream Core. 2017. R: A language and environment for statistical computing. R Foundation for Statisitical Computing, Vienna, Austria. https://www.r-project.org/. Reinert, H.K., and D. Cundall. 1982. An imporved surgical implantation method for radiotracking snakes. Copeia 1982: 702–705. Robson, L.E., and G. Blouin-Demers. 2013. Eastern Hognose Snakes (Heterodon platirhinos) avoid crossing paved roads, but not unpaved roads. Copeia 2013: 507–511. Rodgers, A.R., A.P. Carr, H.L. Beyer, L. Smith, and J.G. Kie. 2007. HRT: Home range tools for ArcGIS. Ontario Ministry of Natrual Resources, Centre for Northern Forest Ecosystems Research. Row, J.R, and G. Blouin-Demers. 2006. Kernels are not accurate estimators of home-range size for herpetofouna. Copeia 2006: 797–802. 68 Saumure, R.A., T.B. Herman, and R.D. Titman. 2010. Effects of patch size and habitat structure on the movements of adult male Wood Turtles, Glyptemys insulpta. Herpetological Conservation Biology 5: 403–413. Shine, R., M. Lemaster, M. Wall, T. Langkilde, and R. Mason. 2004. Why did the snake cross the road? Effects of roads on movement and location of mates by garter snakes (Thamnophis sirtalis parietalis). Ecology and Society 9: 9p. Shipley, B.K., D. Chiszar, K.T. Fitzgerald, and A.J. Saviola. 2013. Spatial ecology of Prairie Rattlesnake (Crotalus viridis) assocaited with Black-tailed Prairie Dog (Cynomys ludovicianus) colonies in Colorado. Herpetological Conservation and Biology 8: 240– 250. Slade, N.A., and S.M. Blair. 2000. An empirical test of using counts of individuals captured as indices of population size. Journal of Mammalogy 81: 1035–1045. Statistics Canada. 2014. Canada's population estimates: Subprovincial areas, July 1, 2014. Taylor, E.N., M.A. Malawy, D.M. Browning, S.V. Lemar, and D.F. DeNardo. 2005. Effects of food supplemenation on the physiological ecology of female Western Diamond-backed Rattlesnakes (Crotalus atrox). Oecologia 144: 206–213. Tiebout, H.M., and J.R. Cary. 1987. Dynamic spatial ecology of the water snake, Nerodia sipedon. Copeia. 1987: 1–18. Torre, I., L. Freixas, A. Arrizabalga, and M. Diaz. 2016. The efficiency of two widely used commercial live-traps to develop monitoring protocols for small mammal biodiversity. Ecological Indicators 66: 481–487. Vistnes, I., C. Nellemann, P. Jordhøy, and O. Strand. 2004. Effects of infrastructure on migration and range use of wild reindeer. Journal of Wildlife Management 68: 101–108. Waldron, J.L., S.M. Welch, S.H. Bennett, W.G. Kalinowsky, and T.A. Mousseau. 2013. Life history constraints contribute to the vulnerability of a declining North America rattlesnake. Biological Conservation 159: 530–538. 69 Wilcove, D.S., D. Rothstein, J. Dubow, A. Phillips, and E. Losos. 1998. Quantifying threats to imperiled species in the United States. BioScience 48: 607–615. Willson, J.D., and J.W. Gibbons. 2009. Drift fences, cover boards and other traps. Pp. 229–245 In Dodd, K, Jr. (Ed.). Amphibian Ecology and Conservation: A handbook of techniques. Oxford University Press, New York, New York, USA. Wolfe, A.K., P.W. Bateman, and P.A. Fleming. 2017. Does urbanization influence the diet of a large snake? Current Zoology 64: 311–318. Yermokhin, M.V., V.G. and C.A. Ivanov. 2015. Spawning migration phenology of the Spadefoot Toad Pelobates fuscus (Pelobatidae, Amphibia) in the valley of the Nedveditsa River (Saratov oblast). Biology Bulletin 42: 931–936. Zappalorti, R., J. Burger, and F. Peterson. 2015. Home range size and distance traveled from hibernacula in Northern Pinesnakes in the New Jersey pine barrens. Herpetologica 71: 26–36. 70 CHAPTER THREE: CONCLUSION AND MANAGEMENT IMPLICATIONS SUMMARY My thesis contributes to the broad understanding of the impacts of disturbance, habitat loss and fragmentation on animal behaviour. Specifically, I used a small-bodied, migrant predator to investigate shifts in migration patterns due to disturbance and humaninfluenced changes on the landscape. Using radio-telemetry and GIS analysis I investigated alterations to spring migration patterns of Western Rattlesnakes in the southern interior of British Columbia, and the potential implications these movement alterations had on individual body characteristics and home ranges throughout the active season. To recap, the general objectives of my thesis were to: (1) Attain baseline information of rattlesnake spring migration movements in shrubsteppe ecosystems in B.C., (2) Outline how disturbance and barrier fencing influence rattlesnake spring migration movements, and (3) Describe the influence of alterations in spring migration has on rattlesnake body characteristics and home range attributes throughout the active season. The major results of my thesis are:  The distance individual rattlesnakes traveled until they first encountered a disturbance had a large influence on spring migration patterns and home range sizes.  Encountering disturbance (both landscape and barrier fencing) reduced both the total distance moved in the spring, and path sinuosity of spring migration in rattlesnakes.  Individuals who encountered the snake fence barrier had reduced spring migration distances, home range and core area sizes compared to all other snakes in the study. 71  Rattlesnakes migrating shorter distances to first encounter disturbance were smaller than snakes moving further distances to first encounter disturbance in the spring. Contrary to my predictions, I did not observe a significant influence on end of year body condition for rattlesnakes based on their migration experience. My results also showed no significant difference in growth rates between rattlesnakes encountering a fence, other forms of disturbance or moving through undisturbed habitats during their spring movements. However, modest sample sizes indicated relatively low power (in both body condition and growth rate) in these analysis suggesting further work may be warranted. CONSERVATION AND MANAGEMENT IMPLICATIONS The province-wide threat assessment for Western Rattlesnakes in B.C. has been classified as high (Southern Interior Reptile and Amphibian Recovery Team 2016), with the primary threats identified include direct population loss from road mortality and human persecution. Furthermore, lower-ranked threats include habitat loss and fragmentation (specifically from agriculture, housing, and recreation). These threats are magnified by the south Okanagan Valley representing one of the fastest growing regions in Canada (Statistics Canada 2014) resulting in a considerable reduction of natural grassland and shrub steppe ecotypes associated with many species-at-risk, including rattlesnakes (Lea 2008). Moreover, less than 15% of the remaining land area within the rattlesnake’s range in B.C. is represented by provincial parks, protected areas, and ecological reserves (Grassland Conservation Council of British Columbia 2004). Animals do not abide by jurisdictional boundaries, so within protected areas, animals with large home ranges and/or migration patterns, like rattlesnakes, may leave these protected areas (Thirgood et al. 2006; Harvey 2015; Allen and Singh 2016). This results in areas of reduced habitat quality (i.e., golf courses, urban areas, vineyard and orchards, etc.) becoming a common aspect of rattlesnake home ranges and habitats used throughout the active season and lifetime of individuals. Therefore, understanding shifts in rattlesnake ecology within disturbed landscapes will lead to more thoughtful and effective management strategies for these animals within less desirable habitat and high human influence. 72 My thesis work provides important baseline ecological information to help address several key objectives outlined for the recovery of this species, namely: (1) “Continue to clarify the impact of human disturbance (both direct and through habitat alteration) on snakes and develop effective mitigative strategies”, (2) “Address biological knowledge gaps regarding distribution, movements, population structure, metapopulation dynamics, prey relationships, genetic and landscape connectivity, and health impacts from threats”, and (3) “Identify and strategically reduce movement barriers in terrestrial habitat where loss of habitat and connectivity is seriously affecting population viability”. The objectives listed above are all identified as essential priorities for rattlesnake conservation in B.C. (Southern Interior Reptile and Amphibian Recovery Team 2016). Based on the results of my thesis, I make the following specific management recommendations: 1) Ensure habitat complexity along fencing structures Barrier fencing has been proven to be effective at reducing direct negative snake-human conflict (Colley et al. 2017) as well as protecting other herpetofauna species primarily from road mortality (Aresco 2005; Baxter-Gilbert et al. 2015; Markle et al. 2017). However, undoubtedly these fences also likely prevent animals from accessing previously used areas prior to construction (Peaden et al. 2017). Fences can also have direct consequences to animals as previously observed along the exclusion fencing near the campground (on the OIR study site), where dead snakes were observed likely dying due to overheating while trying to navigate this feature (Eye et. al. 2018). Ensuring a complex habitat mosaic along these structures should be considered during the planning and construction phases of establishing new exclusion fencing. The addition of complex vegetation and rock cover provides an important mosaic of different thermal opportunities along these structures (Croak et al. 2010). Secondly, it ensures snakes would have access to important habitats for basking, shedding and mating. Lastly, the addition of cover habitat would increase habitat opportunities for small mammals and other species, increasing the hunting opportunities for snakes while navigating exclusion fencing. 73 2) Continue long-term population monitoring and research on OIR The Nk’Mip Snake Program is one of the longest running research programs of its kind in western Canada. But, their relatively small sizes, cryptic nature, often irregular distribution and predominantly nocturnal behaviour makes studying snakes challenging and obtaining long-term population data difficult (Maida et al. in 2018). Long-term data accumulation and research is a crucial component to develop insightful conservation plans for long-lived species-at-risk (Blouin-Demers et al. 2002; Schneider et al. 2018). The OIR study site also contains a unique situation of a drastic landscape gradient from highly disturbed and fragmented with high human use to near-pristine with limited human influence on the landscape. Furthermore, the OIR study site provides the opportunity to monitor the rattlesnake population as development and human involvement (likely) continues to increase in many areas. Considering habitat loss and fragmentation accounts for one of the main threats to rattlesnakes in B.C., this makes the Nk’Mip program very valuable for rattlesnake conservation moving forward. Furthermore, the OIR study site contains two other at-risk snake species: the Great Basin Gopher Snake (Pituophis catenifer deserticola) and the Western Yellow-bellied Racer (Coluber contrictor) providing a great opportunity to determine important baseline demographic assessments for these species as well as species’ niche partitioning and comparing influences of disturbance on these three meso-predator snakes. KNOWLEDGE GAPS AND FUTURE RESEARCH CONSIDERATIONS 1) Fine-scale behavioural and physiological assessment of fencing The scale of my data collection and analysis did not allow me to assess the fine-scale movements of rattlesnakes directly associated with responding to the fencing barriers. Therefore, an important next step is investigating how rattlesnakes respond and behave when they are directly trying to navigate a fencing barrier. Comparing and contrasting fine-scale behaviour as well as physiological impacts such as stress hormones and body temperatures on individuals navigating fencing will help lead to more effective implementations of these structures. Considering the results of my thesis, it is clear that fencing alters rattlesnake movement and behaviour. 74 2) Long term spatial and behavioural assessment of individuals along fencing Adult rattlesnakes do not appear to have the behavioural plasticity to redistribute themselves on the landscape (Waldron et al. 2013). My results outline a short term behavioural response to encountering fencing during rattlesnake spring movements. Investigating the long term (annual) spatial and behavioural movements of individual snakes coming into contact with this fencing would be insightful. Understanding how migratory and spatial behaviour potentially changes (or remains static) over time in response to fencing and other influences or on the landscape would be beneficial to conservation and management techniques of this species. 3) Rattlesnake ecology in different land-use types Urban development and agriculture are two primary land-use disturbances on the landscape in the Okanagan Valley. Specifically, wineries and vineyards account for approximately 4,000 ha of the valley bottoms of the Okanagan Valley (not including orchards), and there are approximately 30 golf courses between Osoyoos and Vernon (approximately 140 km). Vineyards, orchards and golf courses have become three common land-use types that rattlesnakes either move through or remain in throughout their entire active season, however previous information on this animal in these land use-types are lacking. Detailed and targeted ecological information on rattlesnakes in these areas will help promote improved conservation strategies for these snakes and will likely enhance both snake and human safety. 4) Hibernaculum proximity to disturbance My results indicated that individual snakes travelling shorter distances to disturbance were smaller (SVL) than individual snakes that travel further to disturbance or didn’t encounter disturbance at all on the OIR. Previously, Lomas (2013) observed that snakes frequenting undisturbed areas appeared larger than rattlesnakes in various disturbance intensities on the same OIR study site. Both of these examples suggest that disturbance appears to consistently affect rattlesnake size, and the proximity of hibernacula in regards to disturbance may have a large impact on individuals within that community. Assessing the influence that proximity to disturbance has on individual hibernacula communities may assist with targeted management approaches and priorities. 75 5) Ecological drivers of movement/migration in rattlesnakes in B.C. As described previously in this thesis, a number of factors have been hypothesized to drive rattlesnake movement in temperate regions (Bauder et al. 2015; Harvey 2015; Duvall et al. 1997). Previous studies suggest Prairie Rattlesnake (Crotalus viridis viridis) vernal (spring) movement is dedicated to finding food (Duvall et al. 1990) based on a food supplementation experiment. Whereas, water supplementation does not seem to influence Western Rattlesnake spatial ecology in the southern limits of the species’ range (California Capehart et al. 2016). Individuals in some hibernacula populations in B.C. show different migration and home range characteristics than individuals on the OIR study site. Within the valley bottoms, individuals in the OIR study site typically establish identifiable summer home ranges (Lomas et al. in press). Conversely, in other areas of B.C., rattlesnakes embark on a directional outward migration from their hibernaculum, reach an identifiable “turnaround” point, and then return back to the hibernaculum (Gomez et al. 2015; Harvey 2015). Identifying the motivations for rattlesnake movement in B.C. can lead to the discovery of this varying phenomena and improve managing for those potential ecological drivers on disturbed and fragmented landscapes. 6) “Mountain snakes” Rattlesnakes are associated with both the dry, arid valley bottoms and also higher elevation forests throughout their range in southern B.C. (Lomas 2013; Gomez et al. 2015; Harvey 2015). Individuals moving from their hibernaculum into higher elevations typically move further distances on average than individuals moving into the valley bottoms (Lomas 2013; Gomez et al. 2015; Harvey 2015), and “mountain” snakes on the OIR study site grow larger and exhibit higher body conditions than valley individuals (Lomas 2013). However, intensive radio-telemetry and mark-recapture has been limited on these individuals mainly due to logistical and accessibility hurdles. Rattlesnake moving into the higher elevations appear to consist of a large portion of the OIR population (For example six of 15 telemetry rattlesnakes in 2015 moved into the mountains). Based on differences in landscape and vegetation structures, prey composition (Larsen personal communications) and thermal properties, these individuals likely have drastically different experiences throughout the summer active season than snakes in the arid valley bottoms. However, the “mountain” 76 segment of the OIR rattlesnake population has not been investigated other than to define their migration direction in research-based initiatives. CONCLUSION Spring migration is an important aspect of snake ecology in temperate regions, and my results are the first to provide baseline information on spring migration patterns and behaviour of Western Rattlesnakes. However, like many other migratory species (large or small), landscape disturbance and barriers appear to be restricting natural spring migration behaviour in rattlesnakes on the OIR study site. For a species with a slow life history, elevated adult survival is crucial for the viability of this species. Adult Western Rattlesnakes likely do not have the behavioural plasticity to respond to disturbance and anthropogenic threats by redistributing themselves on the landscape (Waldron et al. 2013). Therefore, barrier fencing is crucial for reducing adult mortality due to these anthropogenic threats (i.e., roads, urban communities, etc.). My thesis can act as an important point of reference for managers to consider when implementing mitigative barrier fencing for snakes or other reptile species moving forward, as it is clear that these fencing structures are not without impact to individuals. Moving forward, special attention must be given to understanding the trade-offs fencing has on individuals to ensure long-term efficacy of these important structures. Overall, the results, conclusion and recommendations moving forward outlined in this thesis can help contribute to more thoughtful and meaningful conservation and management strategies for this species and help promote a healthy coexistence between rattlesnakes and humans in B.C. 77 LITERATURE CITED Allen, A.M., and N.J. Singh. 2016. Linking movement ecology with wildlife management and conservation. Frontiers in Ecology and Evolution 3: 13pp. Aresco, M.J. 2005. Mitigation measures to reduce highway mortality of turtles and other herpetofauna at a north Florida lake. Journal of Wildlife Management 69: 549–560. Bauder, J.M., H. Akenson, and C.R. Peterson. 2015. Movement patterns of Prairie Rattlesnakes (Crotalus v. viridis) across a mountainous landscape in a designated wilderness area. Journal of Herpetology 49: 377–387. Baxter-Gilbert, J.H., J.L. Riley, D. Lesbarrères, and J.D. Litzgus. 2015. Mitigating reptile road mortality: Fence failures compromise ecopassage effectiveness. Plos One 10: 15pp. Blouin-Demers, G. K.A Prior, and P.J. Weatherhead. 2002. Comparative demography of Black Rat Snakes (Elaphe obsoleta) in Ontario and Maryland. Journal of Zoology 256:1–10. Capehart, G.D., C. Escallon, B.J. Vernasco, I.T. Moore, and E.N. Taylor. 2016. No drought about it: Effects of supplemental hydration on the ecology, behavior and physiology of free-ranging rattlesnakes. Journal of Arid Environments 134: 79–86. Colley, M., S.C. Lougheed, K. Otterbein, and J.D. Litzgus. 2017. Mitigation reduces road mortality of a threatened rattlesnake. Wildlife Research 44: 48–59. Croak, B.M., D.A. Pike, J.K. Webb, and R. Shine. 2010. Using artificial rocks to restore nonrenewable shelter sites in human-degraded systems: Colonizing by fauna. Restoration Ecology 18: 428-438. Duvall, D., M.J. Goode, W.K. Hayes, J.K. Leonhardt, and D.G. Brown. 1990. Prairie Rattlesnake vernal migration: Field experimental analysis and survival value. National Geographic Research 64: 457–469. Duvall, D., D. Chiszar, R.A. Mintzer, and N.O. Roster. 1997. Experimental simulation in behavioral ecology: a multimedia approach with the spatial searching simulation RattleSnake ©. Experimental Biology Online 2: 11pp. 78 Eye, D.M., J.R. Maida, O.M. McKibbin, K.W. Larsen, and C.A. Bishop. 2018. Snake mortality and cover board effectiveness along exclusion fencing in British Columbia, Canada. Canadian Field Naturalist 132: 30–35. Grassland Conservation Council of British Columbia. 2004. Grasslands mapping project: A conservation risk assessment. 116 pp. Gomez, L., K.W. Larsen, and P.T. Gregory. 2015. Contrasting patterns of migration and habitat use in neighboring rattlesnake populations. Journal of Herpetology 49: 371– 376. Harvey, J.A. 2015. Thermal influences on summer habitat use by Western Rattlesnakes (Crotalus oreganus) in British Columbia. MSc. Thesis, Thompson Rivers University, Kamloops, BC. 94pp. Lea, T. 2008. Historical (pre-settlement) ecosystems of the Okanagan Valley and lower Similkameen valley of British Columbia- pre-Europe contact to the present. Davidsonia 19: 33p. Lomas, E.V. 2013. Effects of disturbance on the Northern Pacific Rattlesnake (Crotalus oreganus oreganus) in British Columbia. MSc. Thesis, Thompson Rivers University, Kamloops, BC. 117pp. Lomas, E., C.A. Bishop, J.R. Maida, and K.W. Larsen. In press. Movement ecology of the Northern Pacific Rattlesnake (Crotalus oreganus oreganus) in response to disturbance. Herpetelogica. Maida, J.R., D.A. Kirk, O. McKibbin, J.R. Row, K.W. Larsen, C. Stringam, and C.A. Bishop. 2018. Population estimate, survivorship and generation time of the Northern Pacific Rattlesnake (Crotalus o. oreganus) at its northern-most range limits. Herpetological Conservation and Biology 13: 662–672. Markle, C.E., S.D. Gillingwater, R. Levick, and P. Chow-Fraser. 2017. The true cost of partial fencing: Evaluating strategies to reduce reptile road mortality. Wildlife Society Bulletin 41: 342–350. 79 Peaden, M.J., J.A. Nowakowski, T.D. Tuberville, K.A. Buhlmann, and B.D. Todd. 2017. Effects of roads and roadside fencing on movements, space use, and carapace temperatures of a threatened tortoise. Biological Conservation 214: 13–22. Schneider, A.C., T.W. Arnold, P.W. Huber, and T.L. Lewis. 2018. An 18-year markrecapture study of Wood Turtles (Glyptemys insculpta) in Michigan. Journal of Herpetology 52: 193–200. Statistics Canada. 2014. Canada’s population estimates: Subprovincial areas, July 1, 2014. Southern Interior Reptile and Amphibian Recovery Team. 2016. Recovery strategy for the Western Rattlesnake (Crotalus oreganus) in British Columbia. 37 pp. Thirgood, S., A. Mosser, S. Tham, G. Hopcraft, E. Mwangomo, T. Mlengeya, M. Kilewo, J. Fryxell, A.R.E. Sinclair, and M. Borner. 2006. Can parks protect migratory ungulates? The case of the Seregenti wildebeest. Animal Conservation 7: 113–120. Waldron, J.L., S.M. Welch, S.H. Bennett, W.G. Kalinowsky, and T.A. Mousseau. 2013. Life history constraints contribute to the vulnerability of a declining North American rattlesnake. Biological Conservation 159: 530–538.