<<

Behavior and Spatial Ecology of the Harris' Antelope Ground (Ammospermophilus harrisii)

Item Type text; Electronic Thesis

Authors Burnett, Alexandra D.

Publisher The University of .

Rights Copyright © is held by the author. Digital access to this material is made possible by the University Libraries, University of Arizona. Further transmission, reproduction, presentation (such as public display or performance) of protected items is prohibited except with permission of the author.

Download date 28/09/2021 10:13:42

Link to Item http://hdl.handle.net/10150/636507 BEHAVIOR AND SPATIAL ECOLOGY OF THE HARRIS’ ANTELOPE GROUND

SQUIRREL (Ammospermophilus harrisii)

by

Alexandra Burnett

______

Copyright © Alexandra Burnett 2019

A Thesis Submitted to the Faculty of the

SCHOOL OF NATURAL RESOURCES & THE ENVIRONMENT

In Partial Fulfillment of the Requirements

For the Degree of

MASTER OF SCIENCE

NATURAL RESOURCES – WILDLIFE CONSERVATION & MANAGEMENT

In the Graduate College

THE UNIVERSITY OF ARIZONA

2019 2

THE UNIVERSITY OF ARIZONA GRADUATE COLLEGE

As members of the Master’s Committee, we certify that we have read the thesis prepared by: Alexandra Burnett titled:

and recommend that it be accepted as fulfilling the thesis requirement for the Master’s Degree.

______Date: ______Sep 9, 2019 John L Koprowski

______Date: ______Sep 10, 2019 Michael Bogan

______Date: ______Sep 9, 2019 bret pasch

Final approval and acceptance of this thesis is contingent upon the candidate’s submission of the final copies of the thesis to the Graduate College.

I hereby certify that I have read this thesis prepared under my direction and recommend that it be accepted as fulfilling the Master’s requirement.

______Date: ______Sep 9, 2019 John L Koprowski T hesis Committee Chair

School of Natural Resources and the Environment

3

Acknowledgments

I would like to sincerely thank my graduate advisor, Dr. John Koprowski for providing the opportunity for me to complete my master’s degree at the University of Arizona, as well as his support and guidance throughout my thesis project. I would further like to thank my committee members, Dr. Bret Pasch and Dr. Bill Mannan, for their availability and support along the way.

Special thanks to Sarah Hale, Melissa Merrick, and Vicki Greer for training support and field assistants, Alexis Blair, Sandy Slovikosky, Helena Yomantas, Anne-Laure Blanche, Megan

Bethel, Nicole Bokanowski, Brandon Mayer, and Eduardo Gracia for countless hours of field support and data input. Thanks to M. Merrick and M. Mazzamuto for guidance with software and data analysis, and V. Greer for logistical support with permitting, equipment, and vehicles. I am grateful for the support and guidance provided by the Koprowski Conservation Research

Laboratory throughout my time in Arizona. Finally, I thank my partner, Kevin Wagner, and my parents, friends, and family, for their support and patience throughout my thesis project.

Alexandra D. Burnett

4

Table of Contents

List of Figures……………………………………………………………………………………..5

Abstract……………………………………………………………………………………………8

Chapter 1: Introduction…………………………………………………………………………..10

Chapter 2: Present Study…………………………………………………………………………15

Summary of Conclusions…………………………………………………………………………15

Literature Cited……………………………………………………………………………….….17

Appendix A: On using a nonhibernating squirrel to inform models of complexity……..………23

Introduction………………………………………………………………………………………24

Methods…………………………………………………………………………………………..30

Results……………………………………………………………………………………………33

Discussion………………………………………………………………………………………..36

Appendix B: Effects of shrub encroachment on antipredator behavior…………….……………59

Introduction………………………………………………………………………………………60

Methods………………………………………………………………………………………..…62

Results……………………………………………………………………………………………66

Discussion………………………………………………………………………………………..67

5

List of Figures

A.1 Figure 1: Average home range size (ha) of Ammospermophilus harrisii in southeastern

Arizona throughout 2017 and 2018. Data are divided by the mating season (Jan-Apr), juvenile season (May-Aug), and nonbreeding season (Sep-Dec). There are no statistical differences of home ranges sizes among season or sex.………………………………………………………54

A.1 Figure 2: 95% and 50% kernel densities of individual antelope

(Ammospermophilus harrisii) in southeastern Arizona. Home ranges of a) all individuals found within study site 2 and 3 of the Santa Rita Experimental Range in which we collected enough points for a home range. Males are represented by warmer colors (red, orange, yellow) and females are represented by cooler colors (green, blue, purple). Individuals overlapped intra- and intersexually: b) female (green)-male (red) overlap during the months of June-August 2018, b) female-female overlap during the months of June-August 2017, and c) male-male overlap during the months of October and November 2017…………………………………...………………...55

A.1 Figure 3: Average number of antipredator vocalizations given by Harris’ antelope ground squirrels (Ammospermophilus harrisii) detected per hour in each month of 2017 and 2018

(pooled) in southeastern Arizona. Dashed lines represent seasonal boundaries. Both years start with a low hourly rate of antipredator vocalizations until April in the mating season (months 1-

4), remain at a moderate level through the juvenile season (months 5-8), and peak in the nonbreeding season after juvenile dispersal (months 9-12) …………………………………...56

6

A.1 Figure 4: Spectrogram of an adult Harris’ (Ammospermophilus harrisii) call given in southeastern Arizona, looking at frequency (kHz) over time(s). Color represents amplitude; warmer colors represent higher amplitudes. Calls feature a trill with high-frequency overtones that descends in amplitude over time. Squirrels frequently give calls repeatedly, in bouts……………………………………………………………………………………………...57

A.1 Figure 5: Spectrogram of a juvenile Harris’ antelope squirrel (Ammospermophilus harrisii) call in the presence of a Western diamondback rattlesnake (Crotalus atrox) in southeastern

Arizona. Frequency (kHz) is represented on the y-axis and time (s) is represented on the x-axis.

Chirps consist of a short, broadband syllable that is repeated.…………………………………..58

A.2 Figure 1: Selection indices of burrow locations of Ammospermophilus harrisii in southeastern Arizona within each vegetation class based on direct observations compared to average availability across study sites. We included all burrow locations that were used by marked individuals at least once (n=22 , 192 total burrow locations). Negative values indicate selection against, positive values indicate selection for the vegetation class. Burrows beneath Opuntia were used most commonly. and Celtis were the next most frequent genera for burrowing beneath, though at a significantly lower frequency. Individuals also used burrows beneath Senegalia and Cylindropuntia, as well as Ferocactus, Ephedra, and

Heteropogon……………………………………………………………………………………78

A.2 Figure 2: Selection indices of foraging items within each vegetation class based on direct observations of feeding compared to average availability across study sites (n=47 direct

7 observations of feeding). Negative values indicate selection against, positive values indicate selection for the vegetation class. We observed Ammospermophilus harrisii feeding on and fruits of Opuntia the vast majority of the time in southeastern Arizona. The next most frequently observed food source was grass (Poaceae)…………………………………...79

A.2 Figure 3: Selection indices of alarm calling locations within each vegetation class compared to average availability of each vegetation class across study sites. Only observations in which the the individual was calling from was positively identified were included in the analyses

(n=18 observations). Negative values indicate selection against, positive values indicate selection for the vegetation class. Ammospermophilus harrisii called from mesquite (Prosopis) in

38.8% of known alarm call locations in southeastern Arizona. Prickly pear (Opuntia) were the next most common observed alarm call locations (33.3%), and then cholla (Cylindropuntia,

11.1%) and hackberry (Celtis, 11.1%)

…………………………………………………………...80

A.2 Figure 4: Selection indices of active individuals within each vegetation class compared to average availability across study sites (n=22 individuals, 257 observations). Negative values indicate selection against, positive values indicate selection for the vegetation class.

Ammospermophilus harrisii in southeastern Arizona selected for bare ground slightly avoided shrub and cacti………………………………………………………..…………………81

8

Abstract

Ground squirrels provide a valuable model system for studying both theoretical and applied questions in behavioral ecology. High variability among ecology imposes selective forces that result in vast differences in social and communication systems. Resource distribution, body size, and kin selection drive social evolution, leading to hypotheses linking larger body size and shorter active periods to greater sociality, which in turn has implications for strategies that may have interesting implications for models of social and communicative complexity. We conducted a two-year study of a nonhibernating squirrel, Ammospermophilus harrisii, examining social and antipredator behavior to test models of complexity. If nonhibernating squirrels comply with predictions of social hypotheses, A. harrisii should be nonsocial and solitary during the winter due to low body size and long active period. We further tested predictions of communicative complexity relating to call function and variability. Our results suggest that antipredator vocalizations may serve multiple functions and vocalizations contain graded variation that may hold additional layers of information. We found that squirrels were not always solitary during winter months, which could explain evidence of kin selection in calling behavior. Nonhibernating squirrels therefore offer insight into potential mechanisms driving variation in communicative complexity.

We also used A. harrisii to answer conservation questions about rapid environmental changes. Habitat degradation and loss is a significant driver of extinction, limiting species resilience to rapid change and resulting in widespread declines. Behavior is often flexible and context-dependent, allowing researchers to gauge wildlife responses to disturbances and determine whether population decline will result. Monitoring whether novel behaviors arise or survival strategies falter can be the first indicator of a potential threat. Further, understanding

9 how animals respond to environmental stimuli is useful for predicting responses to future disturbance. The Harris’ antelope is representative of many species adapted to open, arid environments in that they rely on vigilance and antipredator signals to evade predators. However, on nearly every continent are shifting from open to mesquite, a phenomenon known as shrub encroachment. Shrub encroachment can alter foraging behavior and communication, and a shift in species composition can ultimately ensue.

We conducted a two-year study of A. harrisii in the Santa Rita Mountains, monitoring individuals in areas of high and low velvet mesquite (Prosopis velutina) cover to assess reliance on mesquite as a resource and quantify differences in antipredatory behavior. We found evidence that the presence of mesquites may alter or reduce antipredator behaviors that may affect adult and juvenile survival. A. harrisii were also found to select against woody vegetation and areas of thick vegetation. Thus, shrub encroachment may have adverse effects on antelope squirrel populations and limit future distribution.

10

Chapter 1: Introduction

Behavioral research has a rich history in wildlife studies that has led to the discovery of a vast array of diversity in survival and reproductive strategies. Decades of behavioral research have provided a framework that allows researchers to use behavior as a tool to understand how environmental selective pressures evolve behavioral adaptations. Many conventional questions regarding human evolution can be answered by turning to systems. For example, the evolution of human society has long been regarded as unique in complexity, however, research into other animal societies has revealed that non-human animals also exhibit a remarkable spectrum in their social relationships and communication. Examining these dynamic characteristics across an array of different ecologies allows evolutionary biologists a means to understand what processes shape these complex behaviors. In an applied setting, behavior is useful as a way of gauging the status of an animal population because it can respond to environmental changes within short time frames but has lasting effects on reproduction and survival. Changes in wildlife behavior that are maladaptive or reduce survival may be indicative of larger issues on the landscape and can serve as a warning system to managers. Behavioral monitoring is becoming increasingly important as climate change brings rapid landscape changes that can decimate populations before multi-year reproductive data can be sufficiently gathered to identify problems. Behavior is therefore an important tool for understanding the underlying principles of animal societies, including that of humans, as well as how anthropogenic actions can put populations at risk.

Of particular interest to many researchers is the evolution of seemingly altruistic or cooperative behaviors, such as sociality and communication. Further study of animal systems has unveiled a large spectrum of social and communication behaviors both within and among taxa

11

(Michener 1984, Blumstein 2007, Schel et al. 2010, Bouchet, Blois-Heulin, and Lemasson

2013). Many of the underlying components of complex communication networks, such as referential signaling (‘words’) and ordering (‘syntax’), are most often seen in species with greater social cohesion (Cheney and Seyfarth 1980, Zuberbühler 2000, Blumstein 2007).

However, hints of network structure are apparent across a variety of social and nonsocial contexts as researchers begin to uncover how interactions between selective pressures can shape communication systems overall. Signal meaning or function can vary based on the presence or absence of additional signals and one signal can have multiple functions directed toward multiple receivers (Hebets et al. 2016). Social and communication behaviors are constantly under selection of physiological, phylogenetic, genetic, and both abiotic and biotic environmental interactions. Thus, a more holistic, systems approach is needed when answering proximate and ultimate questions. The antipredator communication system of ground squirrels has been used to model social (Armitage 1981, Rayor and Armitage 1981, Blumstein and Armitage 1998) and communicative complexity (Blumstein and Armitage 1997b, Pollard and Blumstein 2012) and is frequently cited as evidence for the social complexity hypothesis, in which greater sociality is hypothesized to lead to greater communicative complexity (Freeberg, Dunbar, and Ord 2012,

Sewall 2015). Antipredator vocalizations used by squirrels likely function as both a predator deterrent and a warning to offspring (Blumstein and Armitage 1997a, Digweed and Rendall

2009, Fuong, Maldonado-Chaparro, and Blumstein 2015) and phylogenetic studies show that predator deterrence likely evolved first in Rodentia (Shelley and Blumstein 2005). However, kin selection plays a large role in calling behavior and call structure in a number of species (Sherman

1977, Dunford 1977, Blumstein et al. 1997), resulting in a majority of research focus devoted toward the social implications of antipredator vocalizations in hibernating squirrels that live in

12 large family groups (Ackers and Slobodchikoff 1999, McCowan and Hooper 2002, Blumstein,

Verneyre, and Daniel 2004, Matrosova et al. 2010, Pollard 2011). Nonhibernating squirrel exhibit sociality that may vary throughout the year (Karasov 1983, Williams et al. 2009) and can be readily observed in both social and nonsocial contexts. Antipredator communication systems in nonsocial, nonhibernating could provide insight into how selective pressures other than sociality may shape communication structure due to minimized effects of kin selection.

Behavior not only offers researchers a glimpse into the evolutionary process, but also provides a useful gage for studying the effects of rapid environmental change. Animals rely on communication to find food and mates, minimize predation, and reinforce social bonds.

Disruption in communication systems can therefore lead to a breakdown in reproduction and survival. Successful transmission of signals between signalers and recipients depends on the surrounding environment, however, vast landscape changes over very short time intervals affect the transmission of visual, acoustic, and olfactory signals and alter the context in which individuals may communicate. Antipredator signals typically depend on specific environmental qualities in order for both parties to benefit from the signal (Wiley and Richards 1978, Forrest

1994, Briefer et al. 2009). Stotting is an antipredator behavior in Artiodactyla that signals to predators when they’ve been spotted and communicates their fitness to the predator, deterring further pursuit (Caro 2004). Perceived predation risk can influence habitat use and may prompt herbivores to use open environments preferentially (Valeix et al. 2009). Rodents living in open habitats that experienced environmental changes resulting in less visibility increased perceived predation risk and altered foraging behavior (Arenz and Leger 1997, Wheeler and Hik 2014).

Thus, prey may rely on open visual environments to detect and display to the predator as well as for the predator to detect and assess such signals (Lima 1988). Rodents are unlikely to signal

13 unless they are able to track the predator before betraying their location with attention-grabbing behavior such as tail flicking or alarm calling (Blumstein and Armitage 1997a, Shelley and

Blumstein 2005, Owings 2010). Rodents in grasslands or desert scrub that have evolved antipredator cues for open environments may therefore experience detrimental effects when environmental change prevents normal antipredator behavior from occurring (Arenz and Leger

1997, Wheeler and Hik 2014, Akunke and O’Connell 2017). Shrub encroachment is an ongoing problem in arid grasslands worldwide (Archer 1995, Van Auken et al. 2010), and many areas, including Southwest North America, are already experiencing biodiversity losses due to shrub encroached habitat (Ratajczak et al. 2012, Smit and Prins 2015). Examining the effects of shrub encroachment on antipredator behavior and habitat selection may shed light onto the mechanism of population decline and possible management strategies that may help mitigate or minimize species decline before local extirpation.

Ammospermophilus harrisii offers an excellent model to answer the theoretical and conservation questions outlined above. Ammospermophilus harrisii is a desert-dwelling sciurid that emits variable alarm calls but is thought to be solitary, providing a good system to fill gaps in current models of sociality and communication. Individuals are active year-round, providing an opportunity to study antipredator vocalizations in both social contexts, in which juveniles are present, as well as nonsocial contexts. By observing calling behavior with and without juvenile presence, we can identify signal receivers and call function, offering insight into selective forces other than sociality could shape communication systems. A. harrisii is also endemic to the

Sonoran Desert but is generally associated with more upland desert scrub, which is currently facing shrub encroachment of velvet mesquite (Prosopis velutina) (Archer 1995, Bahre and

Shelton 1993). In the Santa Rita Experimental Range, A. harrisii occur across an elevational

14 gradient that exhibits variable levels of shrub encroachment. This provides a natural experiment allowing us to compare calling behavior between shrub encroached and areas that remain open.

If antipredator behavior is reduced in shrub encroached areas, predation may increase and population numbers may decline. These findings would have broader implications for other desert and grassland species that have evolved antipredator strategies that rely on open environments, many of which threatened or endangered.

15

Chapter 2: Present Study

This thesis is a compilation of two manuscripts prepared for publication, located in

Appendices A and B. Appendix A: On using a nonhibernating squirrel to inform models of complexity, formatted for Animal Behaviour, discusses the sociality and communication of

Harris’ antelope squirrel (Ammospermophilus harrisii) in relation to the study of complexity. We focus on the function of antipredator vocalizations and subsequent implications for the current structure of the antipredator communication system. Appendix B: Effects of shrub encroachment on antipredator behavior, formatted for Conservation Biology, describes avoidance of shrub cover by squirrels and shows how shrub encroachment can reduce predator-deterring behavior, which could indicate higher foraging costs and altered predator-prey interactions in dense areas.

Below is a summary of the conclusions of these studies.

Summary of Conclusions

We used the Harris’ antelope ground squirrel (Ammospermophilus harrisii) to inform models of social and communicative complexity as well as management decisions relating to shrub encroachment. In our first study, we found that A. harrisii vocalizes year-round but that callers are female-biased. We also found that A. harrisii emits at least two different vocalizations in antipredator contexts and that trill vocalizations contained a large amount of variation within and between call bouts. Our findings coincide with other research that suggests that antipredator vocalizations are likely used for both a predator deterrent and a warning to offspring, and vocalizations may therefore be shaped by kin selection. Further, species that are typically classified as solitary may still benefit from neighbors and maintain social bonds throughout their lifetime. Thus, social or kin selection may explain some of the graded variation observed. In a separate study, we found that antipredator behavior in a low-density stud site differed from our

16 shrub encroached study site. We also found that squirrels preferred desert vegetation and bare ground, avoiding woody patches. Shrub encroachment may therefore reduce antipredator behavior adapted to open environments, potentially decreasing survival. This could further imperil grassland species that are already threatened or endangered. Additionally, avoidance of woody patches could indicate that shrub encroachment may further alter population distribution if squirrels select against such areas when dispersing. Other open-habitat species may also exhibit similar avoidance, potentially due to increased costs of vigilance (Arenz and Leger 1997) or decreased likelihood of survival (Akunke and O’Connell 2017) due to reduced visibility.

Managers should employ a variety of tactics to maintain open habitat for grassland and upland desert species (Thomas and Goodson 1992, Twidwell et al. 2013). Further, travel corridors between areas of high biodiversity should be created and maintained to encourage gene flow and increase resilience to climate change.

17

Literature Cited

Ackers, S., & Slobodchikoff, C. (1999). Communication of stimulus size and shape in alarm

calls of Gunnison’s prairie dogs, Cynomys gunnisoni. Ethology, 105, 149–162.

Akunke, F., & O’Connell, T. J. (2017). The landscape of fear as an emergent property of

heterogeneity: Contrasting patterns of predation risk in grassland ecosystems. Ecology and

Evolution, 7, 4782–4793.

Archer, S., Schimel, D. S., & Holland, E. A. (1995). Mechanisms of shrubland expansion: Land

use, climate or CO2? Climatic Change, 29, 91–99.

Arenz, C., & Leger, D. W. (1997). Artificial visual obstruction, antipredator vigilance, and

predator detection in the thirteen0lined ground squirrel ( tridecemlineatus).

Behaviour, 134, 1101–1114.

Armitage, K. B. (1981). Sociality as a life-history tactic of ground squirrels. Oecologia, 48, 36–

49.

Bahre, C. J., & Shelton, M. L. (1993). Historic vegetation change, mesquite increases, and

climate in southeastern Arizona. Journal of Biogeography, 20, 489–504.

Blumstein, D. T., & Armitage, K. B. (1997). Alarm calling in yellow-bellied : I. The

meaning of situationally variable alarm calls. Animal Behaviour, 53, 143–171.

Blumstein, D. T., & Armitage, K. B. (1997). Does sociality drive the evolution of

communicative complexity? A comparative test with ground-dwelling sciurid alarm calls.

The American Naturalist, 150, 179–200.

Blumstein, D. T., & Armitage, K. B. (1998). Life history consequences of social complexity: a

comparative study of ground dwelling sciurids. Behavioral Ecology, 9, 1–7.

18

Blumstein, D. T., Steinmetz, J., Armitage, K. B., & Daniel, J. C. (1997). Alarm calling in

yellow-bellied marmots: II. The importance of direct fitness. Animal Behaviour, 53, 173–

184.

Blumstein, D. T., Verneyre, L., & Daniel, J. C. (2004). Reliability and the adaptive utility of

discrimination among alarm callers. Proceedings of the Royal Society of London B:

Biological Sciences, 271, 1851–1857.

Blumstein, D. T. (2007). The evolution of alarm communication in rodents: structure, function,

and the puzzle of apparently altruistic calling. In J. O. Wolf & P. W. Sherman (Eds.),

Rodent Societies (pp. 317–327). University of Chicago Press, 1427 E 60th St, Chicago, IL

60637-2954 USA.

Bouchet, H., Blois-Heulin, C., & Lemasson, A. (2013). Social complexity parallels vocal

complexity: a comparison of three non-human primate species. Frontiers in Psychology, 4,

1–15.

Briefer, E., Osiejuk, T., Rybak, F., & Aubin, T. (2009). Are song complexity and song

sharing shaped by habitat structure? An information theory and statistical approach. Journal

of Theoretical Biology, 262, 151.

Caro, T. M., Graham, C. M., Stoner, C. J., & Vargas, J. K. (2004). Adaptive significance of

antipredator behaviour in artiodactyls. Animal Behaviour, 67, 205–228.

Cheney, D. L., & Seyfarth, R. M. (1980). Vocal recognition in free-ranging vervet monkeys.

Animal Behaviour, 28, 362–367.

Digweed, S. M., & Rendall, D. (2009). Predator-associated vocalizations in North American red

squirrels (Tamiasciurus hudsonicus): To whom are alarm calls addressed and how do they

function? Ethology, 115, 1190–1199.

19

Dunford, C. (1977). Kin selection for ground squirrel alarm calls. The American Naturalist, 111,

782–785.

Forrest, T. G. (1994). From sender to receiver: Propagation and environmental effects on

acoustic signals. American Zoologist, 34, 644–654.

Freeberg, T. M., Dunbar, R. I. M., & Ord, T. J. (2012). Social complexity as a proximate and

ultimate factor in communicative complexity. Philosophical Transactions of the Royal

Society B, 367, 1785–1801.

Fuong, H., Maldonado-Chaparro, A., & Blumstein, D. T. (2015). Are social attributes associated

with alarm calling propensity? Behavioral Ecology, 26, 587–592.

Hebets, E. A., Barron, A. B., Balakrishnan, C. N., Hauber, M. E., Mason, P. H., & Hoke, K. L.

(2016). A systems approach to animal communication. Proceedings of the Royal Society B:

Biological Sciences, 283, 1–10.

Karasov, W. H. (1983). Wintertime energy conservation by huddling in antelope ground

squirrels (Ammospermophilus leucurus). Journal of Mammalogy, 64, 341–345.

Lima, S. L. (1988). Initiation and termination of daily feeding in dark-eyed juncos: influences of

predation risk and energy reserves. Oikos, 53, 3–11.

Matrosova, V. A., Volodin, I. A., Volodina, E. V, & Vasilieva, N. A. (2010). Stability of

acoustic individuality in the alarm calls of wild yellow ground squirrels Spermophilus fulvus

and contrasting calls from trapped and free-ranging callers. Naturwissenshaften, 97, 707–

715.

McCowan, B., & Hooper, S. L. (2002). Individual acoustic variation in Belding’s ground squirrel

alarm chirps in the High Sierra Nevada. The Journal of the Acoustical Society of America,

111, 1157.

20

Michener, G. R. (1984). Age, sex, and species differences in the annual cycles of ground-

dwelling sciurids: implications for sociality. In J. Murie & G. R. Michener (Eds.), The

biology of ground-dwelling squirrels: Annual cycles, behavioral ecology, and sociality. (pp.

81–107). University of Nebraska Press. Lincoln, NE

Owings, D. H. (2010). Tonic communication in the antipredator behavior of ground squirrels. In

J. H. Mitani, J. Brockmann, T. Roper, M. Naguib, & K. Wynne-Edwards (Eds.), Advances

in the Study of Behavior (Vol. 41, pp. 119–149). Elsevier Inc.

Pollard, K. A., & Blumstein, D. T. (2011). Report social group size predicts the evolution of

individuality. Current Biology, 21, 413–417.

Pollard, K. A., & Blumstein, D. T. (2012). Evolving communicative complexity: insights from

rodents and beyond. Philosophical Transactions of the Royal Society B: Biological

Sciences, 367, 1869–1878.

Ratajczak, Z., Nippert, J. B., & Collins, S. L. (2012). Woody encroachment decreases diversity

across North American grasslands and savannas. Ecology, 93, 697–703.

Rayor, L. S., & Armitage, K. B. (1991). Social behavior and space-use of young of ground-

dwelling squirrel species with different levels of sociality. Ethology Ecology & Evolution,

3, 185–205.

Schel, A. M., Candiotti, A., & Zuberbühler, K. (2010). Predator-deterring alarm call sequences

in Guereza colobus monkeys are meaningful to conspecifics. Animal Behaviour, 80, 799–

808.

Sewall, K. B. (2015). Social complexity as a driver of communication and cognition. Integrative

and Comparative Biology, 55, 384–395.

21

Shelley, E. L., & Blumstein, D. T. (2005). The evolution of vocal alarm communication in

rodents. Behavioral Ecology, 16, 169–177.

Sherman, P. W. (1977). Nepotism and the evolution of alarm calls. Science, 197, 1246–1253.

Smit, I. P. J., & Prins, H. H. T. (2015). Predicting the effects of woody encroachment on

communities, grazing biomass and fire frequency in African savannas. PLoS ONE,

10, 1–16.

Thomas, P. A., & Goodson, P. (1992). Conservation of succulents in desert grasslands managed

by fire. Biological Conservation, 60, 91–100.

Twidwell, D., Fuhlendorf, S. D., Taylor, C. A., & Rogers, W. E. (2013). Refining thresholds in

coupled fire-vegetation models to improve management of encroaching woody in

grasslands. Journal of Applied Ecology, 50, 603–613

Valeix, M., Loveridge, A. J., Chamaille-Jammes, S., Davidson, Z., Murindagomo, F., Fritz, H.,

& MacDonald, D. W. (2009). Behavioral adjustments of African herbivores to predation

risk by lions: spatiotemporal variations influence habitat use. Ecology, 90, 23–30.

Van Auken, O. W. (2009). Causes and consequences of woody plant encroachment into western

North American grasslands. Journal of Environmental Management, 90, 2931–2942.

Wheeler, H. C., & Hik, D. S. (2014). Giving-up densities and foraging behaviour indicate

possible effects of shrub encroachment on arctic ground squirrels. Animal Behaviour, 95, 1–

8.

Wiley, R. H., & Richards, D. G. (1978). Physical constraints on acoustical communication in the

atmosphere: implications for the evolution of animal vocalizations. Behavioral Ecology and

Sociobiology, 3, 69–94.

22

Williams, C. T., Gorrell, J. C., Lane, J. E., McAdam, A. G., Humphries, M. M., & Boutin, S.

(2013). Communal nesting in an “asocial” mammal: social thermoregulation among

spatially dispersed kin. Behavioral Ecology and Sociobiology, 67, 757–763.

Zuberbühler, K. (2000). Referential labelling in Diana monkeys. Animal Behaviour, 59, 917–

927.

23

Appendix A: Use of a nonhibernating squirrel (Ammospermophilus harrisii) to inform models of complexity

Burnett, Alexandra D.*, Koprowski, John L.*

*School of Natural Resources and the Environment, University of Arizona, Tucson, AZ, USA

Ground squirrels offer a fascinating model for studying the evolution of sociality and communication. High variability among species ecology imposes selective forces that result in vast differences in social and communication systems. Resource distribution, body size, and kin selection as main drivers of social evolution, which in turn has implications for the evolution of communication. Conspecifics in many ground squirrel species respond to antipredator vocalizations and sociality functions as a strong selective force favoring more informative antipredator vocalizations. However, comparatively little research exists on identifying other selective forces acting on call structure and function, particularly on communication systems with multiple potential receivers. If predators are the intended receiver of antipredator vocalizations, sociality may not account for all observed variation in antipredator communication systems. We conducted a two-year study to determine the potential function of antipredator vocalizations in a non-hibernating squirrel, Ammospermophilus harrisii. If non-hibernating squirrels comply with predictions of sociality models, A. harrisii should remain solitary throughout the year, exhibit non-overlapping home ranges, and display neutral or avoidant behavior toward other individuals. We further hypothesized that if calls function as a predator- deterrent and are only given when predators are present, callers should be of equal sex ratio and vocalize year-round. If calls function primarily as a warning to offspring, callers should be predominantly female and vocalize only when juveniles are above ground. Additionally,

24 vocalizations should exhibit little variation due to low sociality. We found that callers were predominantly female but vocalized throughout the year. We also found that call bouts varied considerably in duration, amplitude, rate, and call type, which could hold additional layers of information. Our results suggest that antipredator vocalizations may function as a predator deterrent, however, calling behavior may still be shaped by kin selection in solitary species. If antipredator vocalizations have multiple functions and multiple receivers, a number of drivers are likely affecting the evolution of complexity in antipredator vocalizations in ground squirrels.

Organisms exhibit a rich diversity of strategies and behaviors to maximize fitness and ensure survival. In many animal taxa, developing effective, evolutionarily stable strategies depends on coping with tradeoffs that result in the evolution of complex behaviors (Alexander

1974, Arnold 1992). A recent surge of interest in the evolution of complexity and the resulting emergent interactions between system components (Fischer, Wadewitz, and Hammerschmidt

2017, Kappeler et al. 2019, Patricelli and Hebets 2016) has resulted in a systems approach to the study of complex behaviors (Bradbury and Vehrencamp 2014, Fischer et al. 2017, Hebets et al.

2016). Complexity is composed of nonlinear, dynamic systems with interacting emergent properties (Bradbury and Vehrencamp 2014). Comparative analyses within and across phylogenetic clades have uncovered common patterns underlying complexity, focusing on behaviors related to cognition (Seyfarth and Cheney 2015), sociality (Armitage 1981, Benson-

Amram et al. 2016, Blumstein and Armitage 1998, Rayor and Armitage 1991), as well as communication (Bouchet, Blois-Heulin, and Lemasson 2013, Blumstein and Armitage 1997b,

Freeberg, Dunbar, and Ord 2012, Pollard and Blumstein 2012). Though some debate exists over the definition (Kappeler et al. 2019), quantification (Bergman and Beehner 2015, Fischer et al.

25

2017, Kappeler 2019, Peckre, Kappeler, and Fichtel 2019) and comparison of complexity across clades (Kappeler et al. 2019), communication systems that exhibit larger, more graded repertoires or vocalizations that serve multiple functions are generally considered more complex

(Blumstein 2007, Freeberg et al. 2012, Fischer et al. 2017, Hebets et al. 2016).

Complex behaviors, including social group formation, may form due to phylogenetics

(Shelley and Blumstein 2005), physiological and environmental constraints (Armitage 1981,

Blumstein and Armitage 1998, Jarman 1974, Michener 1984), and genetic influences (Alexander

1974, Hamilton 1964, Maynard Smith 1965, Sherman 1977). Because steep trade-offs must be overcome for sociality to be evolutionarily stable (Alexander 1974), predicting the circumstances under which complex social systems evolve is challenging. Sociality can evolve as a strategy to secure resources and ensure survival through times of scarcity, and resource distribution and availability have significant influences on reproductive ecology and social structure in ground squirrels. Mammalian species exhibit larger and more complex groups in geographic locations with patchy resources (Clutton-Brock 1989, Travis and Slobodchikoff 1993, Verdolin and

Slobodchikoff 2009) or short seasonal availability of resources (Armitage 1981, Michener 1984).

Armitage (1981) and Michener (1984) looked comparatively across species, finding that larger bodied species or species with shorter active periods delay dispersal, developing multi-year relationships with parents and siblings. These models focused predominantly on northern latitude species that use hibernation as a strategy to endure resource constraints. Species in arid environments with more variable resource distribution and availability, in turn, exhibit more variable behavioral strategies among species that may result in different social structures

(Dunford 1975, Karasov 1983, Munroe and Koprowski 2014, Waterman 1996). Sympatric species in arid environments can display alternative strategies to cope with seasonal resource

26 depletion, exhibiting both hibernating (i.e. tereticaudus; Neal 1965) or nonhibernating behaviors (i.e. Ammospermophilus harrisii; Neal 1965). Nonhibernating squirrels often still experience seasonal energetic constraints but remain active throughout the year, evolving unique social behaviors. squirrels, which do not hibernate, face resource depletion and thermoregulatory constraints during winter months, that prompt solitary species to nest with neighboring squirrels during winter months (Edelman and Koprowski, 2007, Halloran and

Bekoff, 1994; Koprowski, 1996; Ramos-Lara and Koprowski 2012, Williams et al. 2013).

Ammospermophilus leucurus, a nonhibernating ground squirrel found as far north as Idaho, likewise share burrows with neighbors during winter months, likely to defray the cost of thermoregulation (Karasov 1983). Nonhibernating ground squirrels living in arid regions of

Africa and North America may live in small groups with isolated burrows (Xerus rutilus, O’Shea

1976; Ammospermophilus harrisii, Chew and Chew 1930, Neal 1965; Xerospermophilus tereticaudus, Neal 1965, Bradley 1967) or large family groups (Cynomys ludovicianus,

Hoogland 1983; Xerus inauris, Waterman 1995). Though fewer studies are available, nonhibernating species offer an opportunity to inform models of sociality using species that lack obligatory constraints due to hibernation but must still cope with energetic constraints during wintertime.

The antipredator vocalizations of ground squirrels can similarly serve as a model system for the evolution of communication (Blumstein and Armitage 1997b, Krams et al. 2012,

Pollard and Blumstein 2012). Comparative analyses show support for the social complexity hypothesis, which predicts vocal complexity to increase with greater sociality (Bouchet, Blois-

Heulin, and Lemasson 2013, Blumstein and Armitage 1997b, Krams et al. 2012, Pollard and

Blumstein 2012). Conspecifics glean information such as the type, direction, or urgency of threat

27 from variability in antipredator vocalizations (Rodentia: Leger and Owings 1978, Ackers and

Slobodchikoff 1999, Blumstein and Armitage 1997a, Blumstein and Daniel 2004, Weary and

Kramer 1995; Primate: Schel, Candiotti, and Zuberbühler 2010, Zuberbühler 2001). Many squirrel species further exhibit a female-biased ratio of callers, and reproductive females or females with relatives nearby are more likely to call (Dunford 1977, Hoogland 1983, Leger and

Owings 1978, Schwagmeyer 1980, Sherman 1977). These findings provide evidence that calls function as a warning to kin and kin selection is the principal mechanism driving call structure.

Though there is strong support for the social complexity hypothesis, communication structure and variability are constrained by physiological characteristics of both the caller and the receiver(s) (Guilford and Dawkins 1991), as well as the surrounding environment (Marten and

Marler 1977, Wiley and Richard 1978), which must be considered when evaluating proximate and ultimate drivers of vocalization structure (Patricelli and Hebets 2016, Peckre et al. 2019).

Signals must be detectable by receivers (Guilford and Dawkins 1991, Rowe 2013) and depending on the function or audience of signals, signal structure may be subject to opposing or coinciding evolutionary forces (Arnold 1992, Hunt et al. 2009, Krebs and Dawkins 1984).

Animal communication systems often include signals that have multiple functions (pluripotency:

Guindre-Parker et al. 2012, Manno et al, 2007), increasing the number of selective pressures shaping signal structure (Hebets et al. 2016, Hunt et al. 2009). In particular, signals directed toward both conspecific and heterospecific receivers would impose complex selection pressures that may constrain or select for variation. Thus, studying the signal of interest in a variety of contexts is imperative for understanding the structure and evolution of communication systems and evaluating complexity.

28

Research on antipredator vocalizations has largely focused on conspecific receivers, however, visual and acoustic antipredator signals across taxa are capable of altering predator behavior, prompting movement of the predator to a new location (Artiodactyla: Caro 2004;

Anseriformes: Woodland 1980; Columbiformes: Amorim and Dias 2019; Primate: Adams and

Kitchen 2018, Isbell and Binder 2016, Zuberbühler, Jenny, and Bshary 1999; Rodentia: Barbour and Clark 2012, Clark 2005). An additional hypothesis explaining alarm calling behavior is that prey can decrease the chances of a costly chase by signaling when the predator has been detected

(Hasson 1991, Sherman 1977, Woodland 1980). Phylogenetic studies that map antipredator vocalization behavior in Rodentia found that alarm calling was largely explained by the evolution of diurnality rather than sociality, further supporting the predator-deterrent hypothesis

(Shelley and Blumstein 2005). Antipredator vocalizations may accomplish both functions

(Blumstein and Armitage 1997a, Caro 2004, Digweed and Rendall 2009, Fuong, Maldonado-

Chaparro, and Blumstein 2015, Shelley and Blumstein 2005, Zuberbühler, Jenny, and Bshary

1999). Unfortunately, antipredator behavior in ground squirrels is rarely evaluated to assess how vocalizations affect predator behavior, and studies are largely conducted on social, hibernating species in which juveniles and conspecifics are always present. Teasing apart whether antipredator vocalizations are given regardless of conspecific presence and identifying other potential receivers is therefore extremely difficult. To evaluate antipredator vocalizations in

Sciuridae, development of a systems approach to study antipredator vocalizations in different contexts is necessary. Examining antipredator behavior in a nonsocial context may provide insight into how intraspecific receivers may shape signal structure (Rowe 2013) as well as other possible drivers of complexity (Patricelli and Hebets 2016).

29

Nonhibernating species in desert environments exhibit more variable reproduction and sociality (Waterman 1996), offering an opportunity to examine whether antipredator vocalizations exhibit pluripotency and explore other selective forces imposed by multiple receivers or the environment. Ammospermophilus harrisii occupy the , a semidesert in southwest North America (Audubon and Bachman 1854). A. harrisii is a non- hibernating species of small body size (~120g) that lives at low densities (Chew and Chew 1930) in simple family groups composed of mother and offspring of the year (Neal 1965). Juveniles disperse within their first year and relatedness of neighboring individuals is not known. Sociality models (Armitage 1981, Michener 1984) classify squirrels that obtain sexual maturity and disperse within one year as nonsocial if they live solitarily in the breeding season. Few studies have characterized the social relationships and antipredator vocalizations of Ammospermophilus

(but see Bolles 1988 and Neal 1965) and the was not considered in previous studies modeling complexity, offering an opportunity to test hypotheses in sociality and communication.

If A. harrisii are nonsocial, individuals should remain solitary throughout the year, maintain non- overlapping home ranges, and individuals should exhibit neutral or avoidant behavior toward neighbors.

A. harrisii are active year-round but occur at low densities during the winter season after offspring dispersal (Chew and Chew 1930), allowing researchers to gain insight into the structure and function of antipredator vocalizations in this species. Nonsocial squirrels are predicted by models of communicative complexity to exhibit low variability in antipredator vocalization systems (Pollard and Blumstein 2012). Under the assumption that adults are largely solitary during the nonbreeding season, individuals should call only during the breeding season and caller ratio should be female biased if the primary function of antipredator vocalizations is a

30 warning to offspring. We restrict our hypotheses to offspring and do not include relatives as potential beneficiaries because dispersal patterns and relatedness between neighbors is unknown.

Further, direct fitness likely acts as a stronger driver than indirect fitness (Blumstein et al. 1997), particularly in solitary species. If antipredator vocalizations function to communicate with predators, animals should continue to give antipredator vocalizations throughout the year and caller ratio should be equal to the sex ratio of males and females in the population. Antipredator vocalizations may therefore be shaped by multiple receivers, shedding light on how complexity may arise.

Methods

Study Area

We conducted a two-year study from January 2017 to January 2019 in the Santa Rita

Experimental Range (SRER), a 21,000-ha area in Pima County, Arizona about 32 km south of

Tucson. Located at the western edge of the Sonoran Desert, the SRER boundary lies at ~900m elevation and receives ~300mm of rainfall yearly in a bimodal pattern, driving seasonal vegetation blooms and spikes in wildlife activity (McClaran 2003). This ecotone is largely dominated by cacti and woody shrubs, including prickly pear (Opuntia spp.), mesquite (Prosopis spp.), and desert hackberry (Celtis spp.).

Data Collection

We used Tomahawk live traps (model No. 201) baited with peanut butter to capture individuals, following the American Society of Mammalogists guidelines for humane trapping

(Sikes et al. 2016). We placed traps at localities where we identified squirrel presence. We used a cloth handling cone (Koprowski 2002) when handling squirrels to ensure efficient handling and minimal stress to the squirrels. We recorded weight, sex, life stage, and reproductive condition of

31 captured animals. We fit a VHF radio collar (Wildlife Materials; 3-4g, <5% body weight) to each captured adult and monitored 20 adult A. harrisii of both sexes. We used the homing technique to locate radio collared individuals at least twice weekly throughout the year, and up to eight times weekly during the summer field season until the radio collar died or the individual disappeared. A. harrisii move frequently and widely between locations, such that allowing one hour to pass between tracking an individual was enough time to allow movement to any location within the home range. Because squirrels were small and well camouflaged, individuals were often not seen before running from their original location, indicated by a large change in radio signal strength. In these circumstances, a point was recorded with a GPS unit where the signal was initially strongest to avoid chasing the squirrel and affecting location. We collected behavior and location data for all marked and unmarked squirrels detected visually or audibly. Detections of antipredator vocalizations and general location were collected opportunistically throughout the day. All antipredator vocalizations detected were documented with notes on the general location, number of calls included in the bout, and whether any variation in amplitude, trill duration, or call rate was detected by ear. If possible, we attempted to localize the caller by slowly approaching to within 15m of the caller while recording with a ME66 Sennheiser omnidirectional microphone and a Marantz PMD661 MKIII Recorder. We also used this microphone-recorder system to record any calls given from individuals while confined in a trap or during handling.

Statistical Analysis

We used R version 3.5.0 (R Development Core Team) software to complete analyses.

We installed the “move” package (Kranstauber, Smolla, and Scharf 2018) to determine minimum locations per individual needed to compute a home range. A bootstrap analysis shows home

32 range size plateaus at 20 points, which informed our threshold for minimum number of locations required per individual to estimate a reliable home range. We used the “adehabitathr” package

(Calenge 2006) to calculate 95% and 50% kernel densities for all individuals with enough points

We calculated the smoothing factor calculated with least square cross validation, LSCV, in all cases except to calculate nonbreeding seasonal home ranges, in which an ad hoc method was used because cross-validation criterion could not be minimized. To avoid including points in which squirrel location was influenced by observer presence, all points in which the squirrel moved just before homing, or in which a point had been taken within the prior hour, were eliminated from home range analyses. We used locations and home range boundaries to analyze home range size and overlap.

To estimate proximity to another known individual and therefore get a sense of interindividual overlap both temporally and spatially, we used Ranges software (v 9.1.2) to run an analysis of dynamic interactions. This analysis provides overlap percentages between individual home ranges as well as a Jacob’s index: ( DO - DE / DE + DO) comparing observed distances between 2 individuals of interest (DO) with the mean distance to all possible locations of individual 2 (DE) (Jacobs 1974, Kenward et al. 1993b). The resulting value ranges from -1 to

1, in which a positive value indicates that the animals are found within close proximity relative to home range boundaries and a negative value indicates avoidance between individuals.

We used Welch’s nonparametric t-test for unequal sample sizes to test null hypotheses regarding our predictions of home range size, home range overlap and calling rate between seasons. To compare home range sizes and calling rate among seasons, we divided the year into three seasons: mating (Jan-Apr), juvenile (May-Aug), and nonbreeding (Sep-Dec). We used the earliest and latest dates of juvenile (<90g) detections to inform seasonal boundaries. Seasonal

33 divisions allowed us to examine the effect of juvenile presence on antipredator vocalizations and/or whether propensity to call could be associated with reproductive activity. The three seasons represent times when squirrels are reproductively active (as evidenced by visual assessment of reproductive condition) but no juveniles are aboveground, females are reproductively active and juveniles are above ground, and no squirrels are reproductively active and juveniles have dispersed, respectively. Juveniles were considered dispersed from the natal burrow the week following the last day seen at the burrow or with an adult female. Juveniles and subadults are distinguished from adults by relative size and pelage color, which dulls with age. In some cases, new unmarked individuals appeared in nearby areas and/or within the home range of the reproductive females after juveniles had dispersed from the natal burrow, and so these juveniles were thought to have dispersed to an adjacent or overlapping home range.

Results

Social & Spatial Organization

Adult squirrels interacted or were within close proximity of neighboring adults throughout the year, albeit at low frequency. Another adult squirrel was detected within 20m in

5.5% of observed squirrel locations (n=1654 total observations). Juveniles emerged above ground in mid to late May and dispersed away from the natal burrow in late July or August.

During the juvenile season, unknown female adults were seen within 5m of a reproductively active collared females’ burrow on eight different occasions. On one occasion, we observed an adult female in October interacting with two subadults believed to be her offspring from May that year. Two females shared a burrow together during the winter months of 2017. We observed an adult female emerge from another marked adult female’s burrow on two occasions within one week in December 2017, during which the marked female was either in the burrow or interacting

34 with the unknown adult female. We captured and collared an adult female, believed to be the unknown above, outside of the same marked female’s burrow approximately one month later in

January 2018. Within 7 days of collar deployment, her collar was chewed so we were unable to obtain further movement data. A week later, we observed the unknown female with a scrotal male and one other squirrel of unidentified sex. On a separate occasion, an adult female was also captured outside the burrow of a marked adult female in the morning in October 2018. Collars chewed at the base of the antenna or on the radio transmitter, out of the owner’s reach, provides further evidence that adults interact. Throughout the study period, we observed chewing on six collars deployed and collected outside of the juvenile season (after the month of August and before the month of May).

Home ranges did not differ between sex, season, or reproductive stage (all p values

>0.5). Therefore, subgroups were pooled for remaining analyses. Home range average (±SE) of

95% kernel densities averaged 2.76±0.47 ha and 50% kernel densities of home range averaged

0.63±0.15 ha (n=20 individuals) (Fig. 1), however, home range size varied considerably from 0.5 to 7.5 ha. In the field, we observed individuals shift areas of core activity from week to week and did not always use the entire home range every day. Individuals that were tracked for several months with >50 points often showed a secondary plateau at ~40 points, indicative of a gradual shift in home range over time, consistent with field observations.

Home ranges overlapped intra and intersexually in time and space (Fig. 2). Overall overlap between individuals averaged 44.75 ± 19.8% overlap among 95% kernel density areas and 27.47 ± 5.01% among 50% kernel density areas. Overlap did not differ between male- female, male-male, or female-female subgroups (all p>0.06), except for female-female 50% kernel density overlap, which was smaller than overlap between 50% kernel densities for pooled

35 sexes (Welch’s t-test, t=-2.37 df= 27.095, p=0.025). The mean Jacob’s index (Jacobs 1974,

Kenward et al 1993b) was 0.028 ± 0.019, with no significant differences between male-male, male-female, or female-female indices (all p>0.3).

Antipredator Vocalizations System

A. harrisii gave antipredator vocalizations throughout the year, though the average number of antipredator vocalizations detected per field visit varied between season (Fig. 3). We pooled the breeding and nonbreeding season and compared the resulting call rate with that of the juvenile season to test whether offspring presence affected calling behavior. There was no statistical difference between the juvenile season and the remainder of the year (t= -0.169, df=206.85, p=0.865). Average hourly rates of call detections were initially low in the mating season (mean=0.06 ± 0.02 bouts/hour), increasing in April to a mean of 0.17 ± 0.03 bouts/hour throughout the juvenile season, and peaking in the nonbreeding season (mean=0.27 ±0.06 bouts/hour). Calls differed between mating and juvenile season (Welch’s t-test: t=2.92, df=148.89, p=0.004) as well as the mating and nonbreeding season (Welch’s t-test: t= -3.45, df=76.04, p<0.001), but not between the juvenile and the nonbreeding season (Welch’s t-test: t=

-1.66, df=87.14, p= 0.099).

Detections of spontaneous antipredator vocalizations in which the sex of the caller was positively identified (n=18) were female-biased (88.8% of callers). Although squirrels of both sexes emitted calls from within traps (14.1% of 73 trapped squirrels), callers were female-biased

(73% female). In comparison, the sex ratio of the population based on trapping data was roughly

1:1 female to males, or the proportion of squirrels caught that were female was 0.49 (n=73 total squirrels caught). Juveniles and subadults also produced antipredator vocalizations

36 spontaneously (n=3) and when confined (n=4). All juvenile and subadult callers with positively identified sex were female (n=5).

Harris’ antelope squirrels emitted antipredator vocalizations both singly (n=43 detections) and in bouts (n=111 detections). Antipredator vocalizations (Fig. 4) that included more than one trill call from a consistent distance and direction were defined as an alarm bout.

The average (±SE) number of calls for a bout was 5.7±0.52 vocalizations. Average vocalizations per bout differed between the strategy used to collect vocalizations. When we only listened and tallied the calls (n=67 call bouts), the average (±SE) was 4.3±0.39 vocalizations per bout. When we attempted to localize the caller (n=42 call bouts), however, the average was 8.3±1.12 vocalizations per bout (Welch’s t-test: t=3.41, df=50.08, p=0.001). Calls within bouts varied considerably in amplitude, duration, and call rate. We detected variation in at least one of these elements in 54.4% of call bouts. Alarm call bouts were followed by one or more bouts given by other individuals in 33% of call detections, resulting in up to four individuals calling at the same time. Finally, we documented two additional call types that were previously undescribed. We have observed juveniles and an adult female give a repeated chirp-like call toward threats close to the burrow (<20m) on five occasions (Fig. 5). We documented an additional instance in which an adult female emitted a trill, followed by several chirps, and repeating this sequence four times.

She emitted this call combination in October after her juveniles dispersed.

Discussion

Social & Spatial Organization

Adult social and spatial organization resembles that of other nonhibernating ground squirrels living in scrub or grassland environments (i.e. White-tailed antelope squirrel (A. leucurus, Karasov 1983), African unstriped ground squirrels (Xerus rutilus; O’Shea 1976)), in

37 that adults live in large, overlapping home ranges composed of several isolated burrows and adults may share burrows and interact occasionally. Warm temperatures throughout the year permit extended breeding seasons, such that there are only a few months throughout the year in which neither males or females are reproductively active, however, females produce only one litter a year. These conditions allow for continued parental care or late dispersal among subadults. Breeding and rearing of young were consistent with patterns described by Neal

(1965), who sampled A. harrisii throughout southeastern Arizona. Juveniles emerged above ground in May and dispersed between August and October. Similarly, Nelson’s antelope squirrels (Ammospermophilus nelson) in the San Joaquin Desert emerge aboveground in April, but may not disperse until as late as November (Hawbecker 1958). Whorley and Kenagy (2007) found that breeding and lactation periods were longer in southern populations of white-tailed antelope squirrels (A. leucurus) than northern populations.

We observed evidence of burrow sharing in several marked adults during winter months. Although winter burrow sharing was previously unknown for this species, it has been observed in their sister species, A. leucurus (Karasov 1983), and may be a common behavior for

‘asocial’ species. Tree squirrels, which are nonhibernating and solitary during the breeding season, also exhibit communal nesting during winter months (Edelman and Koprowski, 2007,

Halloran and Bekoff, 1994; Koprowski, 1996; Ramos-Lara and Koprowski 2012, Williams et al.

2013). Burrow sharing could function as a means of thermoregulation during cooler winter months (Karasov 1983, Williams et al 2013) and implies that squirrels may benefit from maintaining social relationships. We also observed two adult females within 5m of a single active burrow on several occasions throughout the study period. Burrow sharing and adult persistence at a reproductively active female’s burrow could also be evidence of extended

38 parental care (Williams et al. 2013). Adults with more individuals at a burrow may also benefit from greater collective vigilance (van der Marel, Lopez-Darias, and Waterman 2019) or decreased foraging trade-offs (Ortiz et al. 2019). A genetic study analyzing the relatedness between neighboring individuals is needed to determine dispersal patterns, informing social organization and whether kin selection is likely to have a strong influence.

Compared to home ranges of other North American ground squirrels, which typically fall below 1 hectare (Evans and Holdenreid 1943, Evans 1951, Smith and Johnson 1985, Drabek

1973), A. harrisii home ranges are quite large. Other antelope squirrels in neighboring also maintain expansive home ranges of similar size, however, ranging from 3 to 6 hectares

(Hawbecker 1958; Bradley 1967). African ground squirrels living in arid environments exhibit large home ranges (O’Shea 1976) and home ranges in arid environments may be larger than that of northern ground squirrels due to lower resource density (Bradley 1967, Fridell and Litvaitis

1991, Fisher and Fisher 2000, Doherty et al. 2019). Home ranges overlapped intra and intersexually in A. harrisii and did not differ between sexes, characteristic of polygamous mating systems (Clutton Brock 1989) in arid environments (Fisher and Fisher 2000).

Antipredator Vocalization System

We found that although both sexes gave antipredator vocalizations, alarm callers were disproportionately female when calls were given freely as well as when confined. Calls were heard year-round, but peaked during the nonbreeding season, after juvenile dispersal.

Antipredator vocalizations also persisted at approach, suggesting that antipredator vocalizations are intended for the predatory threat (Hasson 1991). However, because our initial assumption of winter solitary behavior may not hold, it is possible that antipredator vocalizations do function as a warning to conspecifics and continue to provide social benefit in the winter time, in which case

39 alarm calling behavior would not cease entirely during nonbreeding months. Antipredator vocalization rate peaks during nonbreeding seasons and presence of calling across age and sex indicate that endocrine cycles are not a likely mechanism driving calling behavior.

Communication systems that include multiple vocalizations and/or vocalizations with multiple functions and graded variation represent higher complexity (Freeberg et al. 2012,

Fischer et al. 2017). Our findings suggest that antipredator vocalizations function as both a predator deterrent and a warning to offspring, indicating that selective forces other than sociality shape antipredator vocalizations. A. harrisii communication system exhibited high variation within and among call bouts. Variation may represent a graded response to the predator’s approach or transmit information to the receiver. Graded variation in trill length or amplitude may be perceived as different signals by the receiver (Fischer et al. 2017, Peckre et al. 2019,

Marler 1974) or prevent the intended receiver from habituating to antipredator vocalizations

(Loughry and McDonough 1988, Owings 2010. The number of trills included in a call bout also varied considerably, and call repetition may be used to incite increased vigilance in conspecifics

(i.e. tonic calling) (Loughry and McDonough 1988, Owings and Hennessey 1984, Owings 2010).

Alternatively, variation may represent levels of information that are used to inform conspecifics of threat details (Ackers and Slobodchikoff 1999, Blumstein and Armitage 1997a, van der Marel,

López-Darias, and Waterman 2019). The different call types and call combinations we observed provide further evidence of greater complexity (Freeberg et al. 2012, Hebets et al. 2016, Peckre et al. 2019) and could be designed to alert different species of predators (Guilford and Dawkins

1991) or communicate threat urgency to conspecifics (Blumstein and Armitage 1997a). Finally, we observed multiple individual callers, which could further communicate risk (Weary and

Kramer 1995) or location of the predator (Thompson and 2010).

40

We observed non-aggressive adult interactions, possible burrow sharing, and home range overlap throughout the year, such that kin selection may play a larger role in shaping this species’ sociality and communication than previously thought. Thus, species living in low densities with the ability to communicate long distances could exhibit a higher degree of social selective pressures than has been appreciated traditionally. Research in tree squirrels and ground squirrels has revealed a number of alternative social tactics to offset energetic constraints. A number of tree squirrel species gain thermoregulatory benefits from engaging in winter social behavior and territorial North American red squirrels (Tamiasciurus hudsonicus) benefit from having familiar neighbors (Siracusa et al. 2019), even when unrelated (Siracusa et al. 2017).

Thus, communication with conspecifics may be selected for in species that are primarily regarded as solitary and may result in the higher variation observed than expected given low levels of sociality. Further study is needed to determine whether conspecifics alter behavior in response to antipredator vocalizations and if so, whether call variation is meaningful to conspecifics.

Conclusions

Species considered to be nonsocial by classic definitions may exhibit more complex sociality and communication than models currently predict. Traditional assumptions about seasonal energetic constraints based on Nearctic fauna are not met throughout much of the planet’s geographic area, which can lead to vastly different mammalian social systems (Jarman

1974, Rayor and Armitage 1991, Waterman 1995). Therefore, a broader spectrum of species ecology that includes those in tropical or arid locations, should be included in future comparative analyses. Though A. harrisii is typically seen solitarily and at low densities, as predicted in an arid biome, individuals maintain large, overlapping home ranges, are seen interacting non-

41 aggressively with neighboring squirrels on occasion, and may benefit from long-distance vocalizations. Other non-hibernating squirrels also benefit from neighboring relationships (Eason

2010, Karasov 1983, Waterman 1995), including territorial, ‘asocial’ species (Siracusa et al.

2017, 2019, Williams et al. 2013), suggesting that traditional social systems are not required for social relationships to evolve.

We also found that A. harrisii may direct antipredator vocalizations toward predators, and alarm calling may serve a dual function of deterring predators and warning nearby relatives.

Alarm signals alter predator behavior across taxa (Clark 2005, Isbell and Bidner 2016, Woodland

1980), however other prey animals, both conspecific (Ackers and Slobodchikoff 1999, Blumstein and Daniels 2004, Leger and Owings 1978) and heterospecific (Aschemeier and Maher 2011,

Igic et al. 2019, Magrath et al. 2009, Magrath et al. 2015, Zuberbühler 2000), monitor such signals. Antipredator signals therefore likely serve multiple functions, introducing a number of potential drivers constraining or selecting for variation. If alarm signals do serve both functions, a more expansive view of alarm communication systems that accounts for nonlinear communication involving three or more parties (Bradbury and Vehrencamp 2014) should be evaluated when considering alarm behavior and call structure. By implementing a systems approach, future studies of antipredator vocalizations may find that communicative complexity is not restricted to more social species as initially concluded. Supporting this view, we found evidence of a graded repertoire and multiple call types with combination. Prey species exposed to a variety of predator strategies may result in variation that provides a framework for information transfer to conspecifics to flourish, thereby increasing the benefits of sociality and imposing selection for social and communicative complexity to evolve.

42

Acknowledgements

We would like to thank our reviewers, B. Mannan and B. Pasch, for their comments that

improved the manuscript. We thank our field assistants for their countless hours of tracking and

observing. Finally, we thank T&E Inc., and the American Society of Mammalogists for

providing funding for equipment and field assistance.

Literature Cited

Ackers, S., & Slobodchikoff, C. (1999). Communication of stimulus size and shape in alarm calls

of Gunnison’s prairie dogs, Cynomys gunnisoni. Ethology, 105, 149–162.

Adams, D. B., & Kitchen, D. M. (2018). Experimental evidence that titi and saki monkey alarm

calls deter an ambush predator. Animal Behaviour, 145, 141–147.

Alexander, R. D. (1974). The evolution of social behavior. Annual Review of Ecological Systems,

5, 325–383.

Amorim, P. S., & Dias, R. I. (2019). Non-vocal communication as an anti-predator strategy in

scaled doves (Columbina squammata). Journal of Ethology, 37, 157–165.

Armitage, K. B. (1981). Sociality as a life-history tactic of ground squirrels. Oecologia, 48, 36–

49.

Arnold, S. J. (1992). Constraints on phenotypic evolution. The American Naturalist, 140, 85–107.

Aschemeier, L. M., & Maher, C. R. (2011). Eavesdropping of woodchucks (Marmota monax) and

eastern ( striatus) on heterospecific alarm calls. Journal of Mammalogy,

92, 493–499.

Audubon, J.J. and Bachman, J. (1854). Spermophilus harrisii. In V.G. Audubon (Ed), The

quadrupeds of North America (pp 267-269). New York, NY

43

Barbour, M. A., & Clark, R. W. (2012). Ground squirrel tail-flag displays alter both predatory

strike and ambush site selection behaviours of rattlesnakes. Proceedings of the Royal

Society B: Biological Sciences, 279, 3827–3833.

Benson-Amram, S., Dantzer, B., Stricker, G., Swanson, E. M., & Holekamp, K. E. (2016). Brain

size predicts problem-solving ability in mammalian carnivores. Proceedings of the National

Academy of Sciences, 113, 2532–2537.

Bergman, T. J., & Beehner, J. C. (2015). Measuring social complexity. Animal Behaviour, 103,

203–209.

Blumstein, D. T. (2007). The evolution of alarm communication in rodents: Structure, function,

and the puzzle of apparently altruistic calling. In J. O. Wolff & P. W. Sherman (Eds.),

Rodent Societies (pp. 317–327). University of Chicago Press, 1427 E 60th St, Chicago, IL

60637-2954 USA.

Blumstein, D. T., & Armitage, K. B. (1997). Alarm calling in yellow-bellied marmots: I. The

meaning of situationally variable alarm calls. Animal Behaviour, 53, 143–171.

Blumstein, D. T., & Armitage, K. B. (1997). Does sociality drive the evolution of communicative

complexity? A comparative test with ground-dwelling sciurid alarm calls. The American

Naturalist, 150, 179–200.

Blumstein, D. T., & Armitage, K. B. (1998). Life history consequences of social complexity: a

comparative study of ground dwelling sciurids. Behavioral Ecology, 9, 1–7.

Blumstein, D. T., & Daniel, J. C. (2004). Yellow-bellied marmots discriminate between the alarm

calls of individuals and are more responsive to calls from juveniles. Animal Behaviour, 68,

1257–1265.

44

Blumstein, D. T., Steinmetz, J., Armitage, K. B., & Daniel, J. C. (1997). Alarm calling in yellow-

bellied marmots: II. The importance of direct fitness. Animal Behaviour, 53, 173–184.

Bouchet, H., Blois-Heulin, C., & Lemasson, A. (2013). Social complexity parallels vocal

complexity: a comparison of three non-human primate species. Frontiers in Psychology, 4,

1–15.

Bolles, K. (1988). Evolution and variation of anti-predator vocalization of antelope squirrels,

Ammospermophilus (Rodentia: Sciuridae). International Journal of Mammalian Biology,

53, 129–147.

Bradbury, J. W., & Vehrencamp, S. L. (2014). Complexity and behavioral ecology. Behavioral

Ecology, 25, 435–442.

Bradley, W. G. (1967). Home range, activity patterns, and ecology of the antelope ground squirrel

in southern Nevada. The Southwestern Naturalist, 12, 231–251.

Calenge, C. (2006). The package adehabitat for the R software: a tool for the analysis of space

and habitat use by animals. Ecological Modelling, 197, 516-519

Caro, T. M. (1995). Pursuit deterrence revisited. Trends in Ecology and Evolution, 10, 500–503.

Caro, T. M., Graham, C. M., Stoner, C. J., & Vargas, J. K. (2004). Adaptive significance of

antipredator behaviour in Artiodactyls. Animal Behaviour, 67, 205–228.

Chew, R. M., & Chew, A. E. (1970). Energy relationships of the of a desert shrub

(Larrea tridentata) community. Ecological Monographs, 40, 1–21.

Clark, R. W. (2005). Pursuit-deterrent communication between prey animals and timber

rattlesnakes (Crotalus horridus): the response of snakes to harassment displays. Behavioral

Ecology and Sociobiology, 59, 258–261.

45

Clutton-Brock, T. H. (1989). Mammalian mating systems. Proceedings of the Royal Society of

London B: Biological Sciences, 236, 339–372.

Digweed, S. M., & Rendall, D. (2009). Predator-associated vocalizations in North American red

squirrels (Tamiasciurus hudsonicus): To whom are alarm calls addressed and how do they

function? Ethology, 115, 1190–1199.

Doherty, T. S., Fist, C. N., & Driscoll, D. A. (2019). Animal movement varies with resource

availability, landscape configuration and body size: a conceptual model and empirical

example. Landscape Ecology, 34, 603–614.

Dunford, C. (1977). Behavioral limitation of round-tailed ground squirrel density. Ecology,

58, 1254-1268

Dunford, C. (1977). Kin selection for ground squirrel alarm calls. The American Naturalist, 111,

782–785.

Eason, P. (2010). Alarm signaling in a facultatively social mammal, the southern Amazon red

squirrel Sciurus spadiceus. Mammalia, 74, 343–345.

Edelman, A. J., & Koprowski, J. L. (2007). Communal nesting in asocial Abert’s squirrels: The

role of social thermoregulation and breeding strategy. Ethology, 113, 147–154.

Evans, F. C., & Holdenried, R. (1943). A population study of the Beechey Ground Squirrel in

central . Journal of Mammalogy, 24, 231–260.

Fischer, J., Wadewitz, P., & Hammerschmidt, K. (2017). Structural variability and

communicative complexity in acoustic communication. Animal Behaviour, 134, 229–237.

Fisher, D. O., Owens, I. P. F., & Fisher, D. (2000). Female home range size and the evolution of

social organization in macropod marsupials. Journal of Animal Ecology, 69, 1083–1098.

46

Freeberg, T. M., Dunbar, R. I. M., & Ord, T. J. (2012). Social complexity as a proximate and

ultimate factor in communicative complexity. Philosophical Transactions of the Royal

Society B, 367, 1785–1801.

Fridell, R., & Litvaitis, J. A. (1991). Influence of resource distribution and abundance on home-

range characteristics of southern flying squirrels. Canadian Journal of Zoology, 69, 2589–

2593.

Guilford, T., & Dawkins, M. S. (1991). Receiver psychology and the evolution of animal signals.

Animal Behaviour, 42, 1–14.

Guindre-Parker, S., Gilchrist, H. G., Baldo, S., Doucet, S. M., & Love, O. P. (2012). Multiple

achromatic plumage ornaments signal to multiple receivers. Behavioral Ecology, 24, 672–

682.

Halloran, M. E., & Bekoff, M. (1994). Nesting behaviour of Abert squirrels (Sciurus aberti).

Ethology, 97, 236–248.

Hamilton, W. D. (1964). The genetical evolution of social behaviour. I. Journal of Theoretical

Biology, 7, 1–16.

Hasson, O. (1991). Pursuit-deterrent signals: communication between prey and predator. Trends

in Ecology & Evolution, 6, 325–329.

Hawbecker, A. C. (1958). Survival and home range in the Nelson antelope ground squirrel.

Journal of Mammalogy. 39, 207–215

Hebets, E. A., Barron, A. B., Balakrishnan, C. N., Hauber, M. E., Mason, P. H., & Hoke, K. L.

(2016). A systems approach to animal communication. Proceedings of the Royal Society B:

Biological Sciences, 283, 1–10.

Hoogland, J. L. (1986). Nepotism in prairie dogs (Cynomys ludovicianus) varies with

47

competition but not with kinship. Animal Behaviour, 34, 263–270.

Hunt, J., Breuker, C. J., Sadowskià, J. A., & Moore, A. J. (2009). Male-male competition, female

mate choice and their interaction: determining total sexual selection. Journal of

Evolutionary Biology, 22, 13–26.

Igic, B., Ratnayake, C. P., Radford, A. N., & Magrath, R. D. (2019). Eavesdropping magpies

respond to the number of heterospecifics giving alarm calls but not the number of species

calling. Animal Behaviour, 148, 133–143.

Isbell, L. A., & Bidner, L. R. (2016). Vervet monkey (Chlorocebus pygerythrus) alarm calls to

leopards (Panthera pardus) function as a predator deterrent. Behaviour, 153, 591–606.

Jacobs, J. (1974). Quantitative measurement of food selection a modification of the forage ratio

and Ivlev’s electivity index. Oecologia, 14, 413–417.

Jarman, P. J. (1974). The social organisation of antelope in relation to their ecology. Behaviour,

48, 215–267.

Kappeler, P. M., Clutton-Brock, T., Shultz, S., & Lukas, D. (2019). Social complexity: patterns,

processes, and evolution. Behavioral Ecology and Sociobiology, 73, 1–6.

Karasov, W. H. (1983). Wintertime energy conservation by huddling in antelope ground

squirrels (Ammospermophilus leucurus). Journal of Mammalogy, 64, 341–345.

Kenward, R. E., Marstrom, V., & Karlbom, M. (1993). Post-nestling behaviour in goshawks,

Accipiter gentilis: II. Sex differences in sociality and nest-switching. Animal Behaviour, 46,

371–378.

Koprowski, J. L. (1996). Natal philopatry, communal nesting, and kinship in fox squirrels and

gray squirrels. Journal of Mammalogy, 77, 1006–1016.

48

Koprowski, J. L. (2002). Handling tree squirrels with a safe and efficient restraint. The Wildlife

Society Bulletin, 30, 101–103.

Krams, I., Krama, T., Freeberg, T. M., Kullberg, C., & Lucas, J. R. (2012). Linking social

complexity and vocal complexity: a parid perspective. Philosophical Transactions of the

Royal Society B, 367, 1879–1891.

Kranstauber, B., Smolla, M., and Scharf, A. (2018). move: Visualizing and Analyzing Animal

Track Data. R package version 3.1.0.

Krebs, J., & Dawkins, R. (1984). Animal signals: Mind-reading and manipulation. In J. R. Krebs

& N. B. Davies (Eds.), Behavioural ecology: An evolutionary approach (2nd ed., pp. 380–

402). Oxford University Press, Oxford, UK

Leger, D. W., & Owings, D. H. (1978). Responses to alarm calls by California ground squirrels:

Effects of call structure and maternal status. Behavioral Ecology and Sociobiology, 3, 177–

186.

Loughry, W. J., & McDonough, C. M. (1988). Calling and vigilance in California ground

squirrels: a test of the tonic communication hypothesis. Animal Behaviour, 36, 1533–1540.

Magrath, R. D., Haff, T. M., Fallow, P. M., & Radford, A. N. (2015). Eavesdropping on

heterospecific alarm calls: from mechanisms to consequences. Biological Reviews, 90, 560–

586.

Magrath, R. D., Pitcher, B. J., & Gardner, J. L. (2009). Recognition of other species’ aerial alarm

calls: speaking the same language or learning another? Proceedings of the Royal Society B:

Biological Sciences, 276, 769–774.

Marler P. (1974) Animal Communication. In Krames L., Pliner P., Alloway T. (Eds.),

49

Nonverbal Communication. Advances in the Study of Communication and Affect, vol 62

(pp. 25–50). Plenum Press, 227 West 17th Street, New York, N.Y. 10011

Marten, K., & Marler, P. (1977). Sound transmission for animal vocalization I. Temperate

Habitats. Behavioral Ecology and Sociobiology, 2, 271–290

Maynard Smith, J. (1965). The evolution of alarm calls. The American Naturalist, 99, 59–63.

McClaran, M. P. (2003). A Century of Vegetation Change on the Santa Rita Experimental Range.

USDA Forest Service Proceedings RMRS-P-30, (520), 16–33.

Michener, G. R. (1984). Age, sex, and species differences in the annual cycles of ground-dwelling

sciurids: implications for sociality. In J. Murie & G. R. Michener (Eds.), The biology of

ground-dwelling squirrels: Annual cycles, behavioral ecology, and sociality. (pp. 81–107).

University of Nebraska Press. 111 Lincoln Mall #400, Lincoln, NE. 68508

Munroe, K. E., & Koprowski, J. L. (2014). Levels of social behaviors and genetic structure in a

population of round-tailed ground squirrels (Xerospermophilus tereticaudus). Behavioral

Ecology and Sociobiology, 68, 629–638.

Neal, B. J. (1965). Reproductive habits of round tailed and Harris antelope ground squirrels.

Journal of Mammalogy, 46, 200–206.

O’Shea, T. J. (1976). Home range, social behavior, and dominance relationships in the African

Unstriped Ground Squirrel, Xerus rutilus. Journal of Mammalogy, 57, 450–460.

Ortiz, C. A., Pendleton, E. L., Newcomb, K. L., & Smith, J. E. (2019). Conspecific presence and

microhabitat features influence foraging decisions across ontogeny in a facultatively social

mammal. Behavioral Ecology and Sociobiology, 73, 42

50

Owings, D. H. (2010). Tonic communication in the antipredator behavior of ground squirrels. In

J. H. Mitani, J. Brockmann, T. Roper, M. Naguib, & K. Wynne-Edwards (Eds.), Advances

in the Study of Behavior (Vol. 41, pp. 119–149). Elsevier Inc.

Owings, D. H., & Hennessy, D. F. (1984). The importance of variation in sciurid visual and vocal

communication. In J. O. Murie, and G. R. Michener (Eds.), The biology of ground-

dwelling squirrels (pp. 169–200). University of Nebraska Press, 1111 Lincoln Mall,

Lincoln, NE 68508

Patricelli, G., Hebets, E. A., & Patricelli, G. L. (2016). New dimensions in animal

communication: the case for complexity. Current Opinion in Behavioral Sciences, 12, 80–

89.

Peckre, L., Kappeler, P. M., Fichtel, C., Kappeler, P., Shultz, S., Clutton-Brock, T., & Lukas, D.

(2019). Clarifying and expanding the social complexity hypothesis for communicative

complexity. Behavioral Ecology and Sociobiology, 73, 1–19.

Pollard, K. A., & Blumstein, D. T. (2012). Evolving communicative complexity: insights from

rodents and beyond. Philosophical Transactions of the Royal Society B: Biological

Sciences, 367, 1869–1878.

R Core Team (2017). R: A language and environment for statistical computing. R Foundation for

Statistical Computing, Vienna, Austria.

Ramos-Lara, N., & Koprowski, J. L. (2012). Communal Nesting Behavior in Mearns’s Squirrels

(Tamiasciurus mearnsi). The Southwestern Naturalist, 57, 195–225.

Rayor, L. S., & Armitage, K. B. (1991). Social behavior and space-use of young of ground-

dwelling squirrel species with different levels of sociality. Ethology Ecology & Evolution,

3, 185–205.

51

Rowe, C. (2013). Receiver psychology: a receiver’s perspective. Animal Behaviour, 85, 517–523.

Rowe, C. (1999). Receiver psychology and the evolution of multicomponent signals. Animal

Behaviour, 58, 921–931.

Seyfarth, R. M., & Cheney, D. L. (2015). Social cognition. Animal Behaviour, 103, 191–202.

Shelley, E. L., & Blumstein, D. T. (2005). The evolution of vocal alarm communication in

rodents. Behavioral Ecology, 16, 169–177.

Sherman, P. W. (1977). Nepotism and the evolution of alarm calls. Science, 197, 1246–1253.

Sikes, R. S., Gannon, W. L., & Animal Care and Use Committee. (2011). Guidelines of the

American Society of Mammalogists for the use of wild mammals in research. Journal of

Mammalogy, 92, 235–253.

Siracusa, E., Boutin, S., Humphries, M. M., Gorrell, J. C., Coltman, D. W., Dantzer, B., Lane,

J.E., McAdam, A. G. (2017). Familiarity with neighbours affects intrusion risk in territorial

red squirrels. Animal Behaviour, 133, 11–20.

Siracusa, E. R., Wilson, D. R., Studd, E. K., Boutin, S., Humphries, M. M., Dantzer, B., Lane,

J.E., and McAdam, A. G. (2019). North American red squirrels mitigate costs of territory

defence through social plasticity. Animal Behaviour, 151, 29–42.

Thompson, A. B., & Hare, J. F. (2010). Neighbourhood watch: multiple alarm callers

communicate directional predator movement in Richardson’s ground squirrels,

Spermophilus richardsonii. Animal Behaviour, 80, 269–275.

Travis, S. E., & Slobodchikoff, C. N. (1993). Effects of food resource distribution on the social

system of Gunnison (Cynomys gunnisoni). Canadian Journal of Zoology, 71,

1186–1192.

52

Valeix, M., Loveridge, A. J., Chamaille-Jammès, S., Davidson, Z., Murindagomo, F., Fritz, H., &

MacDonald, D. W. (2009). Behavioral adjustments of African herbivores to predation risk

by lions: spatiotemporal variations influence habitat use. Ecology, 90, 23–30. van der Marel, A., López-Darias, M., & Waterman, J. M. (2019). Group-enhanced predator

detection and quality of vigilance in a social ground squirrel. Animal Behaviour, 151, 43–

52.

Verdolin, J. L., & Slobodchikoff, C. N. (2009). Resources, not kinship, determine social

patterning in the territorial Gunnison’s prairie dog (Cynomys gunnisoni). Ethology, 115, 59–

69.

Waterman, J. (1995). The social organization of the Cape ground squirrel. Ethology, 101, 130–

147.

Waterman, J. (1996). Reproductive biology of a tropical, non-hibernating ground squirrel.

Journal of Mammalogy, 77, 134–146.

Weary, D. M., & Kramer, D. L. (1995). Response of eastern chipmunks to conspecific alarm

calls. Animal Behaviour, 49, 81–93.

Whorley, J. R., & Kenagy, G. J. (2007). Variation in reproductive patterns of antelope ground

squirrels, Ammospermophilus leucurus, from Oregon to Baja California. Journal of

Mammalogy, 88, 1404–1411.

Wiley, R. H., & Richards, D. G. (1978). Physical constraints on acoustical communication in the

atmosphere: Implications for the evolution of animal vocalizations. Behavioral Ecology and

Sociobiology, 3, 69–94.

53

Williams, C. T., Gorrell, J. C., Lane, J. E., McAdam, A. G., Humphries, M. M., & Boutin, S.

(2013). Communal nesting in an “asocial” mammal: social thermoregulation among

spatially dispersed kin. Behavioral Ecology and Sociobiology, 67, 757–763.

Woodland, J., Jaafar, Z., & Knight, M. L. (1980). The “pursuit deterrent” function of alarm

signals. The American Naturalist, 115, 748–753.

Zuberbühler, K. (2000). Referential labelling in Diana monkeys. Animal Behaviour, 59, 917–927.

Zuberbühler, K. (2001). Predator-specific alarm calls in Campbell’s monkeys. Behavioral

Ecology and Sociobiology, 50, 414–422.

Zuberbühler, K., Jenny, D., & Bshary, R. (1999). The predator deterrence function of primate

alarm calls. Ethology, 105, 477–490.

54

Figure 1: Average home range size (ha) of Ammospermophilus harrisii in southeastern Arizona throughout 2017 and 2018. Data are divided by the mating season (Jan-Apr), juvenile season

(May-Aug), and nonbreeding season (Sep-Dec). There are no statistical differences of home ranges sizes among season or sex.

55

a) b)

d) c)

Figure 2: 95% and 50% kernel densities of individual antelope squirrels (Ammospermophilus harrisii) in southeastern Arizona. Home ranges of a) all individuals found within study site 2 and

3 of the Santa Rita Experimental Range in which we collected enough points for a home range.

Males are represented by warmer colors (red, orange, yellow) and females are represented by cooler colors (green, blue, purple). Individuals overlapped intra- and intersexually: b) female

(green)-male (red) overlap during the months of June-August 2018, b) female-female overlap during the months of June-August 2017, and c) male-male overlap during the months of October and November 2017.

56

Figure 3: Average number of antipredator vocalizations given by Harris’ antelope ground squirrels (Ammospermophilus harrisii) detected per hour in each month of 2017 and 2018

(pooled) in southeastern Arizona. Dashed lines represent seasonal boundaries. Both years start with a low hourly rate of antipredator vocalizations until April in the mating season (months 1-

4), remain at a moderate level through the juvenile season (months 5-8), and peak in the nonbreeding season after juvenile dispersal (months 9-12).

57

Figure 4: Spectrogram of an adult Harris’ antelope squirrel (Ammospermophilus harrisii) call given in southeastern Arizona, looking at frequency (kHz) over time(s). Color represents amplitude; warmer colors represent higher amplitudes. Calls feature a trill with high-frequency overtones that descends in amplitude over time. Squirrels frequently give calls repeatedly, in bouts.

58

Figure 5: Spectrogram of a juvenile Harris’ antelope squirrel (Ammospermophilus harrisii) call in the presence of a Western diamondback rattlesnake (Crotalus atrox) in southeastern Arizona.

Frequency (kHz) is represented on the y-axis and time (s) is represented on the x-axis. Chirps consist of a short, broadband syllable that is repeated.

59

Appendix B: Effects of shrub encroachment on antipredator behavior

Burnett, Alexandra D.*, Koprowski, John L.*

*School of Natural Resources and the Environment, University of Arizona, Tucson, AZ, USA

Habitat degradation and loss is a significant driver of extinction, limiting species resilience to rapid change and resulting in widespread declines. Monitoring whether behaviors arise or survival strategies falter can be the first indicator of a potential threat. Management agencies often struggle against time and lack of resources when saving a species already in critical decline, and so recognizing and mitigating threats in a timely fashion is paramount. Further, understanding how animals respond to environmental stimuli is useful for predicting responses to future disturbance. The Harris’ antelope ground squirrel (Ammospermophilus harrisii) is representative of many species adapted to open, arid environments in that they rely on vigilance and alarm signals to evade predators. However, many of these areas are shifting from open desert grassland to mesquite, a phenomenon known as shrub encroachment. We conducted a two-year study of A. harrisii in the Santa Rita Mountains, monitoring individuals in areas of high- and low-density velvet mesquite (Prosopis velutina) cover to assess reliance on mesquite as a resource and quantify differences in antipredatory behavior. We found evidence that the presence of mesquites may alter or reduce anti-predatory behaviors that could affect adult and juvenile survival. A. harrisii were also found to select against woody vegetation and areas of thick vegetation. Thus, shrub encroachment may have adverse effects on antelope squirrel populations and limit future distribution.

60

Introduction

Habitat degradation and loss are significant drivers of extinction (Heinrichs, Bender, and Schumaker 2016), limiting species resilience to rapid change and resulting in widespread declines (Thomas 2004; Mantyka-Pringle, Martin, and Rhodes 2011). Widespread landscape change over short periods of time wildlife isolated from dwindling refuges and unable to adapt before population numbers are decimated beyond recovery (Thomas 2004). Environments with high temperatures that are prone to , including arid grasslands and savannas, are particularly threatened by habitat degradation, fragmentation, and climate change (Mantyka-

Pringle et al. 2011). Shrub encroachment, generally defined as the increase of woody vegetation in areas previously dominated by herbaceous vegetation, is a global phenomenon threatening grassland communities (Archer 1995). The main drivers of shrub encroachment vary by region, climate, and community composition (Archer 1995; Barger et al. 2011; Stevens et al. 2016,

Archer et al. 2017). General trends show that anthropogenic activities such as overgrazing and fire suppression alter plant- interactions that are further complicated by increasing concentrations of atmospheric carbon dioxide and changing hydrological regimes, triggering ecosystem shifts (Archer et al. 1995; Van Auken 2009). Shrub encroachment significantly alters hydrological trends (Pierini et al. 2014), soil composition (Barger et al. 2011), and population distributions, leading to declines in local biodiversity (Bakker 2003; Dutra et al. 2011,

Ratajczack 2012; Smit and Prins 2015). One such mechanism of wildlife decline may be the breakdown of predator-prey dynamics due to shrub encroachment. Animals reliant on open environments for anti-predatory strategies and communication may alter or increase vigilance behavior in response to lateral obstruction (Arenz and Leger 1997; Wheeler and Hik 2014), thereby increasing the costs of foraging (Metcalfe 1984; Edwards and Waterman 2011;

61

Jayadevan, Mukherjee, and Vanak 2018; Ortiz et al. 2019). Many taxa in open, arid environments evolve alarm signals that deter predator pursuit (Primates: Isbell and Binder 2016;

Artiodactyla: Caro et al. 2004; Lagomorpha: Kamlar and Ballard 2006; Rodentia: Clark 2005;

Barbour and Clark 2012; Gruiformes: Woodland et al. 1980; Passeriformes: Jones and

Whittingham 2008), indicating that these signals may function as a means for prey to communicate that the predator is detected. Once the prey is alert, further pursuit is likely a waste of the predator’s energy and the individual gives up the chase. However, this defense is contingent on the ability of the prey to identify and track the predator through the environment before revealing its location (Shelley & Blumstein 2005) as well as the ability of the predator to detect the prey’s signal (Guilford and Dawkins 1991; Rowe 2013; Patricelli and Hebets 2016).

Without adequate visibility, the conditions required for effective communication between the predator and prey are degraded and anti-predatory defenses may be ineffective, resulting in wasted energy and injury potential for both parties.

Species inhabiting North American deserts and grasslands evolved in an open environment in which wildlife employ predator-deterrent strategies that capitalize on high visibility and allow animals to forage efficiently in these risky environments (Berger et al. 1983;

Kamler and Ballard 2006). Geographical barriers such as mountain ranges or expanses of desert also result in hundreds of endemic species in southeastern Arizona that, after epochs of isolation have evolved unique genetic diversity and adaptations to extreme environments (Koprowski et al. 2013). Despite their impressive resilience to arid conditions, southwestern species face many new challenges that threaten their populations, including urban expansion, grazing pressure, and shrub encroachment. While the southwestern United States has been a focus of research on shrub encroachment for decades, behavioral studies examining Sonoran wildlife response to such

62 landscape changes are lacking. If shrub encroachment does affect anti-predatory behavior, there could be negative consequences for offspring recruitment that result in local extirpation. Further, understanding whether open-adapted species select for or against woody vegetation may provide insight into how population distributions may shift with future shrub encroachment.

The Harris’ Antelope Ground Squirrel (Ammospermophilus harrisii) is a small-bodied

(~120g), diurnal animal that exhibits anti-predatory communication suited to open, desert scrubland (Bolles 1988; Best 1990). While A. harrisii are considered common, their range is limited to upland Sonoran Desert scrub and like many endemic Sonoran species, the land they occupy has faced significant disturbance in the past century due to grazing pressure, urban expansion, and shrub encroachment. Much of their core distribution in Southeastern Arizona is scattered across mid-elevation desertscrub and semidesert grasslands currently threatened by encroachment of native velvet mesquites (Prosopis velutina) (Bahre and Shelton 1993). Ground squirrels are particularly sensitive to lateral occlusion, and must compensate for poor visibility by increasing vigilance and scanning events (Arenz and Leger 1997; Wheeler and Hik 2014).

Thus, A. harrisii behavior and space use may differ between areas of high and low woody cover.

To better understand A. harrisii behavioral response to shrub encroachment and predict how encroachment might affect future distributions, we gathered spatial and behavioral data to inform three questions: 1) Are shrubs used as a resource by squirrels? 2) How is A. harrisii antipredator behavior affected by shrub density? and 3) Do squirrels select for or against shrub patches during daily activity?

Methods

Study Area

63

The Santa Rita Experimental Range (SRER) is a 21,000-ha area in Pima County,

Arizona, about 32 km south of Tucson. The SRER boundary lies at the base of the Santa Rita

Mountains at ~900 m elevation and gradually inclines to ~2500m elevation. This area receives

~300 mm of rainfall yearly in a bimodal pattern, driving seasonal vegetation blooms and peaks in wildlife activity (Krausman and Morrison 2003, McClaran 2003). Generally, vegetation shifts from open Sonoran desertscrub to woodland as elevation increases, however, the SRER has a well-documented history of shrub encroachment over the past century and velvet mesquite (P. velutina) is the dominant shrub between 1000m and 1200m (McClaran 2003). Knowledge regarding wildlife response to such disturbance is warranted but severely lacking (Krausman and

Morrison 2003). We focused on a specific study area of the Santa Rita experimental range covering ~5km2 area between a 950 and 1150m elevation gradient where squirrels were most abundant. The vegetation between these elevations is largely dominated by desert scrub and cacti, including prickly pear (Opuntia spp.), mesquite (Prosopis spp.), and desert hackberry

(Celtis spp). We further broke down the area into five main study sites (~0.1km2 each) that generally encompassed the home ranges of collared squirrels to compare sites of high and low mesquite cover.

Data Collection

We used Tomahawk live traps (model No. 201) baited with peanut butter to capture individuals, following the American Society of Mammalogists guidelines for humane trapping

(Sikes et al. 2016). Squirrels were trapped opportunistically due to low squirrel density and trapping success. We used a cloth handling cone (Koprowski 2002) to handle squirrels and record weight, sex, life stage, and reproductive condition of captured animals. We fit a VHF radio collar (Wildlife Materials; 3-4g, <5% body weight) to 20 adult A. harrisii of both sexes.

64

We located individuals at least twice weekly throughout the year, and up to eight times weekly during the summer field season until the radio collar died or the individual disappeared. A. harrisii move frequently and widely between locations, such that allowing one hour to pass between tracking an individual was enough time to allow movement to any location within the home range. Data were not collected when temperatures rose above 40℃, due to safety concerns, resulting in more locations during morning hours. To avoid influencing results, a point was taken where the signal was initially strongest in circumstances where the squirrel ran from its initial location. We collected behavior and location data for all unmarked squirrels detected visually or audibly. Locations with visual detections used in vegetation analyses may be biased toward open areas, where the observer was more likely to detect the squirrel.

Statistical Analysis

We used R version 3.5.0 (R Development Core Team) software to complete home range analyses. We installed the “move” package (Kranstauber, Smolla, and Scharf 2018) to determine minimum locations per individual needed to compute a home range. We then used the

“adehabitathr” package (Calenge 2006) to calculate 95% and 50% kernel densities for all individuals with enough points, (smoothing factor calculated using least square cross validation,

LSCV). To avoid including points in which squirrel location was influenced by observer presence, all points in which the squirrel moved just before homing, or in which a point had been taken within the prior hour, were eliminated from home range analyses. We used locations and home range boundaries to analyze vegetation associations and selection behaviors.

For resource selection analyses, we quantified squirrel use of burrows, food items, and alarm locations within each vegetation class as well as selection behaviors during general activity. We loaded high resolution (10cm) imagery of the study area provided by NEON (2019)

65 into ArcMap (v.10.6.1) and ran a supervised classification analysis to divide the region into four classes: bare ground, herbaceous cover, cacti, and woody plants. Classification output was used to quantify the proportion of each vegetation class over the entire study area as well as within each of our five study sites. Based on maximum likelihood estimates of supervised classification, the study area is composed of ~17% woody shrubs, 21% cacti, 48% bare ground, and 14% herbaceous cover. We used this output to distinguish between high-density mesquite (greater than 25% woody cover) and low-density mesquite study sites (less than 25% proportion woody cover). We calculated Jacobs index to compare the proportions of observations within each vegetation class to the average proportion throughout the entire study area for burrow, forage, and antipredator vocalization locations (O-E)/(O+E) (Jacobs 1974). For these analyses, observed proportions of vegetation class were informed by direct observations. Burrow locations were informed by radio telemetry data on collared individuals and classified into vegetation class based on direct observations. We included observations of both marked and unmarked individuals in selection analyses of forage and alarm calling locations. To compare behavioral differences between high- and low-density sites, we included observations from study site 3, a low-density mesquite site, and study site 5, a high-density mesquite site of equal area. We used

Welch’s t-test for unequal sample sizes to test for differences between visual and auditory detections and trapping success between the two study sites. To quantify second order habitat selection (Johnson 1980), we compared the proportion of each vegetation class available in the study site (~0.1km2) against the proportion of the marked squirrel’s active locations in each vegetation class using Jacob’s index. For this analysis, we included only GPS locations of marked individuals in which we had a visual confirmation of the squirrel. We estimated

66 significance of observed patterns with Pearson’s chi squared test for analyses of all four selection datasets.

Results

Are shrubs used as a resource by squirrels?

Squirrels dug burrows most frequently beneath mature prickly pear plants (63.9% of observations) but also occasionally used woody shrubs such as mesquite (10.4% of observations)

(n= 22 animals, 192 observations). Squirrels selected for cacti and shrubs and against bare ground and herbaceous cover for burrows compared to expected values that assume neutral selection toward vegetation classes (Chi-squared test: χ2= 316.09 under the null hypothesis, p<<0.0001) (Fig. 1). We observed A. harrisii foraging on a variety of food sources (n= 47 observations total), including buds, grass seeds, berries, mesquite beans, and cacti fruits.

Observations of foraging behavior most frequently involved prickly pear fruits and flowers

(69.6% of observations). In contrast, individuals were observed foraging on mesquite beans in

5.3% of observations. Individuals consumed cacti and grasses greater than measured availability and consumed shrubs less than expected given availability (Chi-squared test: χ 2=26.71, p<<0.0001; Fig. 2). When alarm calling, squirrels selected for both cacti and woody vegetation but avoided bare ground and herbaceous cover (Chi-squared test: χ 2=32.33, p<<0.0001, n=18 observations) (Fig 3). Individuals were seen alarm calling in mesquite trees and prickly pear in

38.8% and 33.3% of known locations, respectively. We also observed squirrels alarm calling from cholla (11.1%) and hackberry (11.1%).

How is squirrel behavior affected by shrub density?

67

We compared two years of visual and auditory detections between a low-density (<25% woody cover) study site (study site 3: 251.05 total hours of observation) to a high-density (>25% woody cover) study site (study site 5: 116.7 total hours of observation), both measuring approximately 0.1km2. We observed fewer squirrel detections at our high-density study site compared to our low-density study site (Welch’s t-test: t=2.16, df=267.19, p=0.03). We documented squirrels at a rate of 1.19±0.12 (SE) detections/hour in areas of low-density mesquite and 0.8 ±0.09 detections/hour in areas of high density. Call detection rate also varied between study sites of high and low mesquite density (n=68 call detections, t= 3.94, df=273.62, p=0.0001). Squirrels gave 0.24 ±0.04 antipredator vocalizations per hour at open study sites with

<25% woody cover. At the study site with >25% woody cover, squirrels gave only 0.03±0.02 antipredator vocalizations per hour, totaling to only three detections of antipredator vocalizations over the course of the two-year study. Trapping success did not vary between site (Welch’s t- test: t=0.21, df=56.67, p=0.83), averaging (±SE) 0.09±0.005 successful captures per trap hour.

Do squirrels select for or against shrub patches?

We observed active individuals <1 m from cover in 58.6% of total observations (n=22 animals, 257 observations), 41.4% of active observations were recorded ≧1 m from cover.

Squirrels selected for slightly more bare ground and herbaceous cover during daily activity than given availability but selected slightly against woody vegetation and cacti (Chi-squared test: χ

2=9.92, p=0.019) (see Fig. 4).

Discussion

Shrubs as a resource for squirrels

We observed A. harrisii using woody vegetation, including mesquite trees, as a resource for burrows, forage, and antipredator strategies, however, cacti were favored over

68 shrubs for burrow locations and forage items. Cacti and shrubs act as food and water sources, cover from predators, a watch post to scan for threats, and vegetation for burrowing beneath. A. harrisii are adept at climbing woody vegetation and will spend several minutes foraging, vigilant, or alarm calling in tall vegetation. As the tallest vegetation on the landscape, mesquites offer a valuable vantage point traditionally provided by large cacti, which may explain why shrubs were preferred over cacti for alarm calling locations.

Shrub density affects squirrel antipredator behavior

Despite the resources that mesquite trees provide, the antipredatory behavior of A. harrisii was altered in areas with high-density mesquite. We recorded fewer visual and auditory detections of squirrels. Due to low trapping reliability, we cannot confirm whether fewer detections of the squirrel correlates with a lower density or whether squirrels instead employed more cryptic behavioral strategy to avoid detection. Regardless, antipredator vocalizations were detected on only three occasions in high-density mesquite areas but on hundreds of occasions in low-density mesquite sites, a difference that time spent in each area cannot account for. This result suggests that individuals occupying dense vegetation may alter antipredatory strategies and inhibit predator-deterring behavior. These changes could be in part due to lack of visibility or inability to track the predator through thick vegetation (Arenz and Leger 2000; Wheeler and Hik

2014; Ortiz et al. 2019). Squirrels may also fail to detect predators until they are within a proximity that most individuals choose not to call (Hasson 1991). Alternatively, dense areas could represent a different density or composition of predators that affects perceived predation risk or antipredator strategy. Another possible reason that squirrels fail to call is the lack of reproductive females present in the area, however, our trapping records indicate that at least four reproductive females inhabited the high-density site during the study period.

69

Squirrels select against shrub cover

The Santa Rita population of A. harrisii is reliant on cover for burrows and predator avoidance, but also uses areas of bare ground for foraging while maintaining visibility. Squirrels exhibited positive selection for cacti and woody vegetation during alarm calling, however, for day-to-day activity, squirrels selected against woody vegetation but selected for areas of bare ground. Thus, while mesquite trees could serve as an important resource going forward, squirrels may prefer open areas to maintain visibility while foraging, as seen in S. tridecemlineatus (Arenz and Leger 1997) and other herbivores (Valeix et al. 2009). High woody cover presents challenges for maintaining visibility for predator detection and avoidance, and may require a higher state of vigilance while foraging, decreasing foraging efficiency (Metcalfe 1984; Edwards and Waterman 2011; Jayadevan, Mukherjee, and Vanak 2018; Ortiz et al. 2019). A. harrisii therefore may face steep trade-offs as mesquite density increases. Poor quality habitat could act as a sink if squirrels experience reduced survival (Pulliam 1988) or deter recolonization by new residents, which may ultimately exclude A. harrisii from shrub encroached areas and limit their distribution.

Management Implications

Small mammals are known to play an integral role in community composition by affecting plant distribution and cover (Walsberg 2000; Beck & Wall 2010; Hale, Koprowski, and

Hicks 2013; Ewacha 2016), soil characteristics (Davidson, Detling, and Brown 2012; Ewacha

2016), and insect abundances (Ewacha 2016). These effects may be exacerbated in desert environments with low resource density (Walsberg 2000; Ewacha et al. 2016). A. harrisii consume and transport an abundance of fruits and seeds. We frequently observed A. harrisii carrying fruits and seeds several meters and occasionally brought food items into the burrow.

70

This foraging behavior was found to be extremely important for plant distribution in A. harrisii’s sister species, Ammospermophilus leucurus (Beck & Wall 2010). A. harrisii likely benefits the

Sonoran Desert vegetation in a manner similar to that of A. leucurus and other ground squirrels.

If shrub encroachment does result in a lower density of desert rodents such as A. harrisii, it may compound on shrub competition with cacti by limiting dispersal.

Our findings provide evidence of potential consequences for the survival and distribution of A. harrisii resulting from shrub encroachment. Other Sonoran wildlife reliant on open environments for antipredator behavior, such as Sonoran pronghorn (Antilocapra americana sonoriensis, Berger et al. 1983), antelope jackrabbit (Lepus alleni, Best and Henry 1993), or black-tailed jackrabbits (Lepus californicus, Kamler and Ballard 2006), could be affected in a similar manner. Shrub encroachment creates a more complicated landscape that can instill a higher depredation of small mammals by aerial predators, or provide more cover for predators to hide (Akunke and O’Connell 2017). Animals that have evolved in open environments increase vigilance behavior in closed areas or increase giving up densities in closed areas (Wheeler and

Hik 2014; Jayadevan et al. 2018) which can result in steep foraging tradeoffs that may alter spatial and temporal use of encroached areas or distribution across the landscape (Valeix et al.

2009). To mitigate the negative consequences outlined above, managers should consider employing controlled burns at regular frequencies to maintain grasslands (Thomas and Goodson

1992; Twidwell et al. 2013). Where shrub encroachment persists, managers may consider thinning and pruning to maintain visibility needed for antipredator strategies. Travel corridors between grasslands should be created and maintained to encourage gene flow between populations and increase resilience to climate change.

71

Acknowledgements

We thank our reviewers, B. Mannan and B. Pasch, for their input in improving the manuscript. We also thank the Koprowski Conservation lab for their guidance and support in addition to our numerous field assistants for their contributions. We thank our funding sources

T&E, ASM, and UA for their financial and logistical support. The National Ecological

Observatory Network is a program sponsored by the National Science Foundation and operated under cooperative agreement by Battelle Memorial Institute. This material is based in part upon work supported by the National Science Foundation through the NEON Program.

Literature Cited

Archer, S., Schimel, D. S., & Holland, E. A. (1995). Mechanisms of Shrubland Expansion: Land

Use, Climate or CO2? Climatic Change, 29, 91–99.

Arenz, C. L., & Leger, D. (1997). The Antipredator Vigilance of Adult and Juvenile Thirteen-

lined Ground Squirrels (Sciuridae: Spermophilus tridecemlineatus): Visual Obstruction

and Simulated Hawk Attacks. Ethology, 103, 945–953.

Bakker, K. (2003). A synthesis of the effect of woody vegetation on grassland nesting .

Proceedings of the South Dakota Academy of Science, 82, 119–141

Barbour, M. A., & Clark, R. W. (2012). Ground squirrel tail-flag displays alter both predatory

strike and ambush site selection behaviours of rattlesnakes. Proceedings of the Royal

Society B: Biological Sciences, 279, 3827–3833.

Barger, N. N., Archer, S. R., Campbell, J. L., Huang, C., Morton, J. A., & Knapp, A. K. (2011).

Woody plant proliferation in North American drylands: A synthesis of impacts on

ecosystem carbon balance. Journal of Geophysical Resources, 116, 1–7.

72

Beck, M. J., & Vander Wall, S. B. (2010). by scatter-hoarding rodents in arid

environments. Journal of Ecology, 98, 1300–1309.

Berger, J., Daneke, D., Johnson, J., & Berwick, S. H. (1983). Pronghorn foraging economy and

predator avoidance in a desert ecosystem: implications for the conservation of large

mammalian herbivores. Biological Conservation, 25, 193–208.

Best, T. L., & Henry, T. H. (1993). Lepus alleni. Mammalian Species, 5, 1–8.

Best, T. L., Titus, A. S., Caesar, K., & Lewis, C. L. (1990). Ammospermophilus harrisii.

Mammalian Species, 366, 1–7.

Blumstein, D. T., & Armitage, K. B. (1997). Does sociality drive the evolution of

communicative complexity? A comparative test with ground-dwelling sciurid alarm calls.

The American Naturalist, 150, 179–200.

Bolles, K. (1988). Evolution and variation of antipredator vocalization of antelope squirrels,

Ammospermophilus (Rodentia: Sciuridae). International Journal of Mammalian Biology,

53, 129–147.

Bradley, W. G. (1967). Home range, activity patterns, and ecology of the antelope ground

squirrel in southern Nevada. The Southwestern Naturalist, 12, 231–251.

Caro, T. M., Graham, C. M., Stoner, C. J., & Vargas, J. K. (2004). Adaptive significance of

antipredator behaviour in artiodactyls. Animal Behaviour, 67, 205–228.

Clark, R. W. (2005). Pursuit-deterrent communication between prey animals and timber

rattlesnakes (Crotalus horridus): the response of snakes to harassment displays.

Behavioral Ecology and Sociobiology, 59, 258–261.

73

Davidson, A., Detling, J., & Brown, J. (2012). Ecological roles and conservation challenges of

social, burrowing herbivorous mammals in the world’s grasslands. Frontiers in Ecology

and the Environment, 10, 477–486.

Dutra, H. P., Barnett, K., Reinhardt, J. R., Marquis, R. J., & Orrock, J. L. (2011). Invasive plant

species alters consumer behavior by providing refuge from predation. Oecologia, 166,

649–657.

Ewacha, M. V. A., Kaapehi, C., Waterman, J. M., & Roth, J. D. (2016). Cape ground squirrels as

ecosystem engineers: modifying habitat for plants, small mammals and beetles in Namib

Desert grasslands. African Journal of Ecology. 54, 68–75

Guilford, T., & Dawkins, M. S. (1991). Receiver psychology and the evolution of animal signals.

Animal Behaviour, 42, 1–14.

Hale, S. L., J. L. Koprowski, and H. Hicks. 2013. Review of black-tailed prairie dog

reintroduction strategies and site selection: Arizona reintroduction. In: Gottfried, Gerald

J.; Ffolliott, Peter F.; Gebow, Brooke S.; Eskew, Lane G.; Collins, Loa C. Merging

science and management in a rapidly changing world: Biodiversity and management of

the Madrean Archipelago III and 7th Conference on Research and Resource Management

in the Southwestern Deserts; 2012 May 1-5; Tucson, AZ. Proceedings. RMRS-P-67. Fort

Collins, CO: U.S. Department of Agriculture, Forest Service, Rocky Mountain Research

Station. p. 310-315.

Hawbecker, A. C. (1958). Survival and home range in the Nelson antelope ground squirrel.

Journal of Mammalogy, 39, 207–215.

Heinrichs, J. A., Bender, D. J., & Schumaker, N. H. (2016). Habitat degradation and loss as key

drivers of regional population extinction. Ecological Modelling, 335, 64–73.

74

Jayadevan, A., Mukherjee, S., & Tamim Vanak, A. (2018). Bush encroachment influences

nocturnal rodent community and behaviour in a semi-arid grassland in Gujarat, India.

Journal of Arid Environments, 153, 32–38.

Johnson, D. H. (1980). The comparison of usage and availability measurements for evaluating

resource preference. Ecology, 61, 65–71.

Jones, K. A., & Whittingham, M. J. (2008). Anti-predator signals in the chaffinch Fringilla

coelebs in response to habitat structure and different predator types. Ethology, 114,

1033–1043

Kamler, J. F., & Ballard, W. B. (2006). Ear flashing behavior of black-tailed jackrabbits (Lepus

californicus). American Midland Naturalist, 155, 402–403.

Kenward, R.E., Casey, N.M., Walls, S.S. & South A.B. (2014) Ranges9: For the analysis of

tracking and location data. Online manual. Anatrack Ltd. Wareham, UK.

Koprowski, J. L. (2002). Handling tree squirrels with a safe and efficient restraint. Wildlife

Society Bulletin, 30, 101–103.

Koprowski, J. L., S. L. Doumas, M. J. Merrick, B. Oleson, E. E. Posthumus, T. G. Jessen, and R.

N. Gwinn. 2013. It's lonely at the top: Biodiversity at risk to loss from climate change. In

Gottfried, Gerald J.; Ffolliott, Peter F.; Gebow, Brooke S.; Eskew, Lane G.; Collins, Loa

C., comps. 2013. Merging science and management in a rapidly changing world:

Biodiversity and management of the Madrean Archipelago III and 7th Conference on

Research and Resource Management in the Southwestern Deserts; 2012 May 1-5;

Tucson, AZ. Proceedings. RMRS-P-67. Fort Collins, CO: U.S. Department of

Agriculture, Forest Service, Rocky Mountain Research Station. 593 p. 53–59.

75

Kranstauber, B., Smolla, M., and Scharf, A. (2018). move: Visualizing and analyzing animal

track data. R package version 3.1.0.

Krausman, P and M. Morrison. 2003. Wildlife ecology and management, Santa Rita

Experimental Range (1903 to 2002). In: McClaran, Mitchel P.; Ffolliott, Peter F.;

Edminster, Carleton B., tech. coords. Santa Rita Experimental Range: 100 years (1903 to

2003) of accomplishments and contributions; conference proceedings; 2003 October 30–

November 1; Tucson, AZ. Proc. RMRS-P-30. Ogden, UT: U.S. Department of

Agriculture, Forest Service, Rocky Mountain Research Station.

Kupfer, J. A., Balmat, J., & Smith, J. L. (2005). Shifts in the potential distribution of Sky Island

plant communities in response to climate change. Connecting Mountain Islands and Desert

Seas: Biodiversity and Management of the Madrean Archipelago II. Proceedings RMRS-P-

36, 485–490.

Mantyka-Pringle, C., Martin, T., & Rhodes, J. (2011). Interactions between climate and habitat

loss effects on biodiversity: a systematic review and meta-analysis. Global Change Biology,

18, 1239–1252.

McClaran, M. P. (2003). A Century of Vegetation Change on the Santa Rita Experimental

Range. USDA Forest Service Proceedings RMRS-P-30, 520, 16–33.

Metcalfe, N. B. (1984). The effects of habitat on the vigilance of shorebirds: is visibility

important? Animal Behaviour. 32, 981–985

National Ecological Observatory Network. 2019. Provisional data downloaded from

http://data.neonscience.org on [01-23-2019]. Battelle, Boulder, CO, USA

76

Pierini, N., Vivoni, E., Robles-Morua, A., Scott, R., & Nearing, M. (2014). Using observations

and a distributed hydrolic model to explore runoff thresholds linked with mesquite

encroachment in the Sonoran Desert. Water Resources Research, 50, 8191–8215.

Pulliam, H. (1988). Sources, sinks, and population regulation. The American Naturalist, 132,

652–661.

Ratajczak, Z., Nippert, J. B., & Collins, S. L. (2012). Woody Encroachment Decreases Diversity

Across North American Grasslands and Savannas. Ecology, 93, 697–703.

Rowe, C. (2013). Receiver psychology: a receiver’s perspective. Animal Behaviour, 85, 517–

523.

Shelley, E. L., & Blumstein, D. T. (2005). The evolution of vocal alarm communication in

rodents. Behavioral Ecology, 16, 169–177.

Smit, I. P. J., & Prins, H. H. T. (2015). Predicting the effects of woody encroachment on

mammal communities, grazing biomass and fire frequency in African savannas. PLoS

ONE, 10, 1–16.

Stevens, N., Erasmus, B. F. N., Archibald, S., & Bond, W. J. (2016). Woody encroachment over

70 years in South African savannahs: overgrazing, global change or extinction

aftershock? Philosophical Transactions of the Royal Society B: Biological Sciences, 371,

20150437.

Thomas, P. A., & Goodson, P. (1992). Conservation of succulents in desert grasslands managed

by fire. Biological Conservation, 60, 91–100.

Thomas, C. D., Cameron, A., Green, R. E., Bakkenes, M., Beaumont, L. J., Collingham, Y. C.,

… Williams, S. E. (2004). Extinction risk from climate change. Nature, 427, 145–148.

77

Twidwell, D., Fuhlendorf, S. D., Taylor, C. A., & Rogers, W. E. (2013). Refining thresholds in

coupled fire-vegetation models to improve management of encroaching woody plants in

grasslands. Journal of Applied Ecology, 50, 603–613.

Valeix, M., Loveridge, A. J., Chamaille-Jammès, S., Davidson, Z., Murindagomo, F., Fritz, H.,

& MacDonald, D. W. (2009). Behavioral adjustments of African herbivores to predation

risk by lions: spatiotemporal variations influence habitat use. Ecology, 90, 23–30.

Van Auken, O. W. (2009). Causes and consequences of woody plant encroachment into western

North American grasslands. Journal of Environmental Management, 90, 2931–2942.

Walsberg, G. E. (2000). Small mammals in hot deserts: Some generalizations revisited.

Bioscience, 50, 109–120.

Wheeler, H. C., & Hik, D. S. (2014). Giving-up densities and foraging behaviour indicate

possible effects of shrub encroachment on arctic ground squirrels. Animal Behaviour, 95,

1–8.

Woodland, J., Jaafar, Z., & Knight, M.-L. (1980). The “pursuit deterrent” function of alarm

signals. The American Naturalist, 115, 748–753.

78

Figure 1: Selection indices of burrow locations of Ammospermophilus harrisii in southeastern

Arizona within each vegetation class based on direct observations compared to average availability across study sites. We included all burrow locations that were used by marked individuals at least once (n=22 animals, 192 total burrow locations). Negative values indicate selection against, positive values indicate selection for the vegetation class. Burrows beneath

Opuntia were used most commonly. Prosopis and Celtis were the next most frequent genera for burrowing beneath, though at a significantly lower frequency. Individuals also used burrows beneath Senegalia and Cylindropuntia, as well as Ferocactus, Ephedra, and Heteropogon.

79

Figure 2: Selection indices of foraging items within each vegetation class based on direct observations of feeding compared to average availability across study sites (n=47 direct observations of feeding). Negative values indicate selection against, positive values indicate selection for the vegetation class. We observed Ammospermophilus harrisii feeding on flowers and fruits of Opuntia the vast majority of the time in southeastern Arizona. The next most frequently observed food source was grass seeds (Poaceae).

80

Figure 3: Selection indices of alarm calling locations within each vegetation class compared to average availability of each vegetation class across study sites. Only observations in which the plant the individual was calling from was positively identified were included in the analyses

(n=18 observations). Negative values indicate selection against, positive values indicate selection for the vegetation class. Ammospermophilus harrisii called from mesquite trees (Prosopis) in

38.8% of known alarm call locations in southeastern Arizona. Prickly pear (Opuntia) were the next most common observed alarm call locations (33.3%), and then cholla (Cylindropuntia,

11.1%) and hackberry (Celtis, 11.1%).

81

Figure 4: Selection indices of active individuals within each vegetation class compared to average availability across study sites (n=22 individuals, 257 observations). Negative values indicate selection against, positive values indicate selection for the vegetation class.

Ammospermophilus harrisii in southeastern Arizona selected for bare ground slightly avoided shrub and cacti.