Think Before They Squeak: Vocalizations of the Squirrel Family

Squirrels (Sciuridae) are a diverse group in behavior, morphology, and ecology. This variation is typified by the wide range of vocalizations spanning ground squirrels (Marmotini and Xerini), tree squirrels (Callosciurinae and Sciurini), and flying squirrels (Pteromyini). Squirrels produce calls that range in frequency, modulation, and function, with a complex set of social calls occurring across the family. We review the history of recording methods used in the development of squirrel vocalization repertoires, with emphasis on how the ecology and methodology impact the frequency values reported. The fundamental (F0 – the mean frequency of the fundamental harmonic), dominant (FDom – the frequency of maximum energy or amplitude), minimum (FMin – the minimum frequency of the fundamental harmonic), maximum (FMax – the maximum frequency of the dominant harmonic), and highest harmonic (FHarm – the mean frequency of the highest visible harmonic) frequencies were considered against popular hypotheses that have attempted to explain the evolution of vocal frequency characteristics in terrestrial mammals. These hypotheses include body size, predator avoidance, habitat type, and diel activity pattern. Phylogenetic generalized least squared modeling revealed that body mass and the frequency limits of the methods were the strongest drivers of high-frequency communication. Consistent with popular hypotheses, social squirrels exhibited a broader range of F0, FDom, and FMax than solitary squirrels while habitat openness promoted higher FDom and FHarm. Additionally, nocturnality was significantly associated with higher F0, FDom, and FMax, suggesting that flying squirrels, the only nocturnal squirrels, commonly use high-frequency acoustic signals, a finding that merits further investigation. In conclusion, our review provides a unique insight into the role of behavioral ecology on vocal repertoires and the importance of accurate equipment selection for sampling across a diverse taxon.


INTRODUCTION
Mammalian bioacoustics is an expansive field in which biology, physics, and psychology are used to study the taxonomic diversity of signalers and interpret the complexity of vocal communication. Over the last century or so, the bioacoustic community has seen the development, implementation, and specialization of their research technology. Ever-improving recording devices have allowed advancements from phonetic descriptions of non-human animal calls to the first spectrographic visualizations to the development of software allowing quantitative analysis of calls in minute detail (Terhune, 2011). While phonetic descriptions of the calls served an important purpose at one point in time, the recording of calls provided the first opportunity for a detailed comparison of the differences between species and the evolution of vocal production in mammals through spectrographic analyses. The first acoustic spectrograph was developed in 1946 to study human speech, though the authors suggested that this new device could be used to study environmental sounds, musical instruments and orchestral composition, and potentially animal sounds under the right conditions (Koenig et al., 1946).
Spectrographs were, and continue to be, important because they allowed researchers to study the shape and structure of the sounds being produced, but they also allowed the visualization of calls that cannot be heard by the human ear. Humans are limited to a hearing frequency range of about 20 Hz to a maximum of 20 kHz, though 15-18 kHz is the cut-off for the average human (Masterton et al., 1969). Given that it is difficult for even the best human ears to hear sounds near the upper limit (Masterton et al., 1969), some have argued that 15 kHz should be the cutoff for differentiating between the audible range (what can be heard by humans) and the ultrasonic range (calls that exceed the upper limit of human hearing). However, most publications describe 20 kHz, the absolute maximum of human hearing, as the frequency at which calls can be described as ultrasonic.
Ultrasonic vocalizations (USVs), calls that completely lie in the ultrasonic range (>20 kHz herein), were first studied in cetaceans and bats that had been observed navigating and hunting without visual cues. USVs were first recorded with the aid of frequency-altering heterodynes in bats (Griffin and Galambos, 1941) and were first visualized in porpoises (Wood, 1952). Bats and cetaceans largely remained the focus of ultrasonic studies for decades, due in particular to interest in echolocation. The intensive echolocation focus has likely overshadowed ultrasound-production in many other wild species such as shrews (Blarina and Sorex: Gould et al., 1964;Sorex: Tomasi, 1979), singing mice (Baiomys and Scotinomys: Miller and Engstrom, 2007;Scotinomys: Fernández-Vargas et al., 2011), and even several species of squirrel (Urocitellus richardsonii: Wilson and Hare, 2004;Glaucomys: Murrant et al., 2013). Additionally, the widespread availability of lab mice and rats has led to many studies investigating the function and biomechanical processes involved with USV production. The study of USV production in these rodents has been quite active since the early 1970s (Oswalt and Meier, 1975;Shair, 1978, 1980;Corrigan and Flannelly, 1979;Geyer, 1979;Smith, 1979), leading to a wealth of publications, but also, along with bats, a taxonomic bias in publications about ultrasound use. While still crucial in understanding the physiological and cognitive processes associated with USVs, captive rats and mice have been isolated from natural processes for generations and are not known to produce well-studied wild calls such as predator alarm calls (Blanchard et al., 1991;Hahn and Lavooy, 2005;Sirotin et al., 2014;Kalcounis-Rueppell et al., 2018).
Ultrasound can be a useful communication tool because the short sound waves are highly directional and are not readily heard by most predatory birds, reptiles, and amphibians (Knudsen, 1981;Yamazaki et al., 2004). It has been proposed that the main drivers of ultrasound use for all terrestrial vertebrates are body size, predator avoidance, and increasing the signal-to-noise ratio (Arch and Narins, 2008). Several studies support the hypothesis that body size has a strong relationship to the frequency of many vocalizing animals (avian mass: Ryan and Brenowitz, 1985; mammalian mass: Charlton and Reby, 2016;Martin et al., 2016; mass of mole rats: [Cryptomys, Heterocephalus, Spalax], Credner et al., 1997;length of ground squirrels: Spermophilus, Nikol'skii, 2017). The mass-call frequency relationship is influenced directly by the physical mechanisms of vocal production (i.e., lung capacity, larynx size) and through indirect effects like species ecology (e.g., larger animals using different environments than smaller animals; Ryan and Brenowitz, 1985). In contrast, predator-avoidance is more difficult to quantify and this hypothesis seems to be more supported by anecdotal evidence of vocal shifts toward ultrasonic frequencies when animals are in the presence of predatory birds (Kotler, 1984;Tamura and Yong, 1993;Digweed and Rendall, 2009). Finally, the signal-to-noise ratio hypothesis has been supported by only a few examples, such as the Kihansi spray toad (Nectophyrynoides asperginis) which has shifted all vocal production to USVs to reduce interference from nearby roaring, low-frequency waterfalls (Arch et al., 2011). Similarly, the acoustic adaptation hypothesis, a version of the signal-to-noise idea that focuses on sound propagation relative to environmental characteristics (e.g., open versus closed) instead of ambient sound, postulates that higher frequencies will be employed in open habitats where the shorter sound waves are less likely to be attenuated by the environment before reaching the receiver (Ey and Fischer, 2009). In contrast to the acoustic adaptation hypothesis, it has also been shown that some mammals inhabiting forested habitats have improved hearing sensitivities to continue to use high-frequency calls in closed habitats (Charlton et al., 2019).
Our review will focus on literature covering the recording and reporting of the highly vocal family, Sciuridae. With spectrographic reports of squirrel vocalizations beginning in the 1960s, squirrel genera have been well sampled, and many vocal libraries have been developed and updated with improving technologies. We apply popular and novel hypotheses related to vocal characteristics ( Table 1) to observe which behavioral and ecological traits drive frequency use in squirrels. We also challenge the currently established repertoires in light of new reports of USVs in flying squirrels and call for future research to employ equipment that is intended for non-human animals, and not subject to the limitations of equipment designed for use with humans.

SQUIRREL VOCALIZATIONS -A SHORT HISTORY
The first published spectrographic analyses of squirrels appeared in 1966; the authors used these visualizations to describe a variety of calls of the Uinta ground squirrel (Urocitellus armatus; Balph and Balph, 1966) and the yellow-bellied (Marmota flaviventris) and hoary marmots (M. caligata; Waring, 1966). 1 | Summary of parameters used in the phylogenetic least square (PGLS) analyses and how they relate to popular or novel vocal range hypotheses.

Parameter Hypothesis Citation
Body mass (-) Body size controls the size of vocal producing structures; body size estimated through body mass. Charlton and Reby, 2016;Martin et al., 2016;Arch and Narins, 2008 Sociality: social (+) Social species communicate more often and rely on conspecifics for predator avoidance; higher frequencies are less likely to detected by predatory birds, reptiles, and amphibians. Ramsier et al., 2012;Arch and Narins, 2008 Habitat: openness (+) Acoustic adaptation hypothesis: species will alter call structure to maximize transmission in different environments; closed habitats attenuate high-frequency sound waves rapidly while open habitats do not inhibit transmission. Ey and Fischer, 2009 Diel activity pattern: nocturnality (+) Gliding and nocturnality are intertwined for Sciuridae; nocturnality facilitates the use of echolocation as visual cues are greatly reduced.

Novel
Method limits (+) The frequency limits of the microphone and spectrograph restrict the detection range of the equipment.

Novel
The predicted direction of the relationship to frequency is indicated as positive (+) or negative (−).
While the authors published separately, they used the same equipment to record the squirrels: a unidirectional microphone with a relatively narrow frequency bandwidth (Electrovoice-644 Sound Spot, frequency response of ∼0.05-12 kHz; Balph and Balph, 1966;Waring, 1966) attached to a relatively large portable tape recorder. Waring (1966) used a W-610-EV batteryoperated tape-recorder (weight ∼7.7 kg; Amplifier Corporation of America, 1959) whereas Balph and Balph (1966) used an 11.8 kg Wollensak T 1700 tape recorder (Revere-Wollsenak Division, 2010). The fundamental frequency (F0 -the vibration frequency of the vocal folds which is determined by the shape of the larynx as well as airflow; Fitch et al., 2002;Charlton and Reby, 2016;Riede et al., 2017) and the dominant frequency (F Dom -the frequency with the maximum amplitude or energy; Figure 1) of the calls produced by the ground squirrels and marmots were captured within this narrow frequency range as demonstrated on the spectrographs. However, Waring noted at least one call produced by the yellow-bellied marmot had strong harmonics (resonant frequencies produced at intervals above F 0 ) exceeding 40 kHz, well beyond the dynamic range of the microphone and only detectable with a heterodyne. The dynamic range of the recording microphone is limited by the physical shape of the transducer; this element vibrates when sound is produced, translating the physical sound wave into an electronic signal (Obrist et al., 2010). Similarly, heterodynes (handheld ultrasonic detectors that transpose USVs down to the human hearing range so that researchers can listen to USVs in real-time) can detect a broad range of frequencies with a specialized transducer. However, the transposer (which translates the vibration into an audible sound) is restricted to a narrow bandwidth and the researcher must tune the transducer to listen to different frequencies (Ahlén, 2004).
As Waring (1966) demonstrated, harmonics can lead to vocalizations that partially exceed human hearing abilities (i.e., ultrasonic harmonics that exceed 20 kHz). F 0 is produced as air exits the lungs (called the glottal flow) and without any interruptions, this tone would only consist of F 0 (Titze, 2009). However, if the glottal flow collides with the vocal tract (vocal folds, upper larynx, tongue, mouth, etc.), infinite harmonics will be produced, decreasing in strength as they increase in frequency (Monsen and Engebretson, 1977;Titze, 2009). Harmonics can be strengthened or dampened, depending on the shape of the vocal tract (like the addition of harmonic resonance chambers; Riede et al., 2017) and the presence of strong high-frequency harmonics can be indicative of active manipulation of the vocal tract (Titze, 2009;Fenton et al., 2011). Harmonics can allow animals to extort a larger bandwidth (as in bats; Kingston and Rossiter, 2004;Fenton et al., 2011) or to communicate at a higher frequency (as in anurans and birds; Narins et al., 2004;Arch et al., 2011), though the function of harmonics in most squirrels remains unclear. The yellow-bellied marmot is not the only squirrel for which ultrasonic harmonics have been described. Other genera include Callospermophilus (Eiler and Banack, 2004), Glaucomys (Gilley, 2013;Murrant et al., 2013;Gilley et al., 2019), Sciurus (Zelley, 1971;Lishak, 1982), Spermophilus (Koshev and Pandourski, 2008;Schneiderová, 2008Schneiderová, , 2012Matrosova et al., 2011Matrosova et al., , 2012Matrosova et al., , 2016Schneiderová and Policht, 2012;Schneiderová et al., 2020), Tamias (Blake, 1992), and Urocitellus (Melchior, 1971;Manno et al., 2007). The function of harmonics has been discussed in the context of a few species of ground squirrel (Spermophilus). California ground squirrels (Spermophilus beecheyi) will dampen the strength of higher harmonics to increase the clarity of F 0 and improve the ability for receivers to locate the sender . As well, California ground squirrels habituating to noisy environments have developed stronger harmonics than their counterparts inhabiting quieter environments, improving the signal-to-noise ratio of their calls (Rabin et al., 2003). In contrast, some ground squirrels with larger geographic ranges or with multiple subspecies may alter the strength or number of harmonics to accommodate different habitat or social constraints, resulting in distinct dialects (Eiler and Banack, 2004). So, while all species would have infinite ultrasonic harmonics if recorded under ideal conditions (i.e., facing the recorder directly, maintaining close proximity to reduce attenuation of higher frequencies, remove all background noises, microphone with infinite sensitivity), we consider it of some interest if strong ultrasonic harmonics have only been detected in a subset of free-roaming squirrels.
Since the first spectrographic publications, the calls of over 70 extant squirrel species have been reported. There has been a strong emphasis on alarm calls (reported in 74/89 publications), likely because these are often the highest-amplitude calls, they FIGURE 1 | Examples highlighting the variation of how Sciuridae calls are presented as well as how frequency characteristics were measured on different spectrographs and call types. Freqeuncy characteristics measured for all call types: F 0 (fundamental frequency: the mean frequency of the primary vibrational frequency of the vocal membrane; when the dominant frequency occurs on the fundamental harmonic, F 0 = (F Max + F Min )/2), F Dom (dominant frequency: the frequency with the greatest energy, power or amplitude), F Min (minimum frequency: the minimum frequency of the fundamental frequency), F Max [maximum frequency: the maximum frequency of the fundamental frequency (or of harmonic on which F Dom is measured)] and F Harm (frequency of highest visible harmonic: mean frequency of the highest complete harmonic visible on the spectrograph). (A) Broadband noise produced by a female southern flying squirrel (Glaucomys volans; Murrant et al., 2013). (B) Alarm call ("shrill chirp") recorded in a colony of Columbian ground squirrels (Urocitellus columbianus; Betts, 1976). (C) Alarm call (two call elements identified by the initial researchers and therefore recorded separately, red = element 1, blue = element 2) recorded in a population of European ground squirrels in the Ponor Mountains (Spermophilus citellus; Koshev and Pandourski, 2008). (D) Alarm call (also with two distinct elements identified by red and blue) emitted by a female little ground squirrel (Spermophilus pygmaeus; Nikol'skii, 2019). Attribution details for previously published figures available in Supplementary Table 3. are produced in response to human observers, and they are most obvious to humans given our own sensory bias (Masterton et al., 1969). Therefore alarm calls are likely the easiest to record. Alarm calls function to notify conspecifics of the presence of a predator and usually prompt a behavioral response from the receivers. Blumstein (1999b) spent over a decade looking at the role of alarm calls in marmots (Marmota spp.), noting that several species can communicate the threat level of a predator by modifying the speed or repetitive pattern of the alarm calls (Blumstein and Armitage, 1997a;Blumstein, 1999a). Additionally, many authors have reported that squirrels will produce different alarm calls (distinct in structure and frequency) to identify unique predator types and to elicit a specific behavioral response in the receivers . For example, Richardson's ground squirrels (Urocitellus richardsonii) produce a short chirp (lasting only 75 ms and rapidly decreasing in frequency from 8 to 3.5 kHz) to warn of aerial predators and evoke an immediate retreat response in conspecifics. In contract, a terrestrial predator is identified by a long whistle (lasting around 400 ms with a constant frequency of 10 kHz) which provokes conspecifics to become for alert stand erect (Davis, 1984). Additionally, these squirrels can add a chuck syllable (only 22.5 ms long and decreasing in frequency from about 6 to 2 kHz) to the end of either the chirp or whistle when predators are nearby, eliciting increased vigilance from receivers after their initial behavioral response (Sloan et al., 2005) and like marmots, Richardson's ground squirrels will increase the repetition of their calls as the threat level increases (Warkentin et al., 2001). When predator-specific alarm calls are used in a species, the caller most often only distinguishes between terrestrial and aerial predators [Callosciurus (Tamura and Yong, 1993;Tamura, 1995), Cynomys (Ackers and Slobodchikoff, 1999;Loughry et al., 2019), Marmota (Taulman, 1977;Davis, 1991;Blumstein and Arnold, 1995;Blumstein, 1999a), Otospermophilus (Owings and Virginia, 1978;Owings and Leger, 1980), and Urocitellus (Balph and Balph, 1966;Melchior, 1971)], though some Asian tree squirrels have a third call which elicits snake mobbing behavior in conspecifics (Tamura and Yong, 1993).
While alarm calls have been the focus of most studies, squirrels have many distinct call types ( Table 2) varying in function from mating (male courtship, pre-and postcopulatory calls, and female mating calls), territoriality, facilitating aggressive (threat, combat) or amicable (maternal, group-cohesion) interactions, and communicating internal states (distress). Several species of flying squirrel have also been observed producing calls during their glides and while exploring their environment (Glaucomys sabrinus and G. volans: Murrant et al., 2013;Petaurista leucogenys: Ando and Kuramochi, 2008), suggesting this nocturnal tribe (Pteromyini; Brandt, 1855) may be using echo-based orientation (Siemers et al., 2009). While this hypothesis has been explored once (Chattin, 1969), the recent discovery of ultrasonic vocalizations in North American flying squirrels (Gilley, 2013;Murrant et al., 2013;Eisinger et al., 2016;Gilley et al., 2019) and the vast improvement in ultrasonic recording since the 1960s warrants a re-examination of the hypothesis.
It was not until 2004 that a USV (i.e., a call where F 0 entirely exceeds 20 kHz) was first reported in a squirrel species: the Richardson's ground squirrel (Urocitellus richardsonii; Wilson and Hare, 2004). Richardson's ground squirrels are known to both produce predator-specific calls that differentiate between terrestrial and aerial predators (Davis, 1984) and vary the repetition and call structure to communicate perceived risk (Warkentin et al., 2001;Sloan et al., 2005). USVs were found to be another way for this species to convey the perceived predation risk by varying their calls between weaker, ultrasonic vocalizations and louder, audible vocalizations as the threat imposed by the presumptive predator increases (Wilson and Hare, 2006). When a predator is farther away, it is less likely to detect highly directional USVs and therefore prey individuals can call altruistically while remaining inconspicuous. As the predator becomes more of a threat, the colony benefits from louder, more omnidirectional alarm calls to ensure that the maximum number of receivers are notified. In addition to varying USV use between perceived threats, it was found that there is a significant increase in the frequency of ultrasonic calls when the emerging juveniles reach their highest density, implying that ultrasonic communication may be most beneficial for social cohesion (Wilson and Hare, 2006). It is key to note that the authors did not use a secondary microphone in the audible range (<20 kHz) to confirm that the calls were not biphonic (two unrelated fundamental frequencies are produced at the same time; Matrosova et al., 2012) with one fundamental in the audible range and a second in the ultrasonic range. However, the authors published a follow-up report in which they attempted to quantity biophonic versus ultrasonic calls in the same population by noting ultrasonic call behaviors (thoracic cavity expanding and mouth held in communicative position, sometimes a very faint whistle can be heard by the observer) while recording with an audible microphone and they were able to differentiate between ultrasonic and biphonic calls (termed "mixed calls" by the authors; Wilson and Hare, 2006). Wilson and Hare's (2004) publication challenges the squirrel vocal literature because most of the published repertoires have used microphones that are not sensitive to the ultrasonic range (85% of studies reviewed in our report did not use equipment that encompassed ultrasonic frequencies). By constraining the microphones to ranges that maximize human voice recording, the published repertoires have been unintentionally biased by methodology. While the detection of ultrasonic calls can be extremely rare even when the correct equipment is being used (as in some Spermophilus; Matrosova et al., 2012), since 2004, three new Sciuridae species have been observed to produce USVs. USVs have now been reported in northern (Glaucomys sabrinus: Murrant et al., 2013) and southern flying squirrels (G. volans: Murrant et al., 2013;Eisinger et al., 2016) and noted, though rarely, in the speckled ground squirrel (Spermophilus suslicus: Matrosova et al., 2012). Calls where F 0 partially extend into the ultrasonic range have been also been noted in both the southern (G. volans: Gilley, 2013;Murrant et al., 2013;Eisinger et al., 2016;Gilley et al., 2019) and northern flying squirrels (G. sabrinus: Gilley, 2013;Gilley et al., 2019) as well as the Carolina flying squirrel (G. sabrinus coloratus: Gilley, 2013).
Despite the underrepresentation of USVs in published squirrel vocal repertoires owing to technological limitations and difficulties associated with recording free-ranging individuals, squirrels are a useful taxonomic group to preliminarily evaluate factors controlling the frequency limits of mammalian vocal signals. Squirrels have been well sampled across the taxa and represent a diverse range of social structure, habitat use, and behavior. We used phylogenetic methods to investigate the following five popular and novel hypotheses relating to different frequency characteristics of the vocalizing members of Sciuridae (Table 1). (1) The role of body size, which has been shown to drive frequencies higher as body size decreases, will be explored (Arch and Narins, 2008;Martin et al., 2016). However, it has already been noted in Sciuridae that body size differences attributed to age class (pup versus adult) do not result in within-species differences in vocal characteristics (Matrosova et al., 2007;Swan and Hare, 2008), therefore this relationship is only expected to be significant when making cross-species comparisons. (2) To investigate the potential role of predator avoidance (prey species will favor acoustic ranges that exceed the detection of common predators; Arch and Narins, 2008), we used sociality as a proxy because species that exhibit social complexity are more likely to exhibit increased vocal complexity (Blumstein and Armitage, 1997b;Blumstein, 2007). Primates that exhibit more complex social systems are sensitive to higher frequencies (Ramsier et al., 2012) and mammals with complex social systems are more likely to produce higher frequency vocalizations associated with submissiveness and appeasement (Charlton and Reby, 2016). Additionally, Wilson and Hare (2006) suggest that squirrel social systems in which kin are more clustered, either spatially or temporally from juvenile emergence, are more likely to favor USV production which is less likely to be detected at a distance by an eavesdropping predator. Therefore, we predicted that more socially complex squirrels will vocalize at higher frequencies.
(3) We explored the acoustic-adaptation hypothesis (terrestrial animals will modify call structure depending on their environment; Blumstein, 2007;Ey and Fischer, 2009) by comparing the frequency limits of squirrels in open versus closed habitats, with open habitats expected to promote the use of higher frequencies.
(4) The role of diel activity patterns is interesting as the only nocturnal squirrels are the flying squirrels which might be using vocalizations for echolocation as echo-based orientation (or "echonavigation") is associated with reduced visual cues in other rodents (Ancillotto et al., 2014;Panyutina et al., 2017). Therefore, we predicted that nocturnal species would have higher frequencies than diurnal species. (5) Finally, the relevance of the equipment used will be analyzed by incorporating the method detection limits (defined as the limits of either the microphone or spectrographic analysis, hereinafter "method limits") of the publication methods into the models defined below. We predicted that method limits, particularly the upper-frequency limit, will constrain the perceived vocal ranges of recorded species.

REVIEWING FREQUENCY CHARACTERISTICS
We developed a database beginning with a list of publications that described the vocalizations of squirrels (summarized in Supplementary Table 1). The minimum requirement for each publication was the description of at least one call with either a spectrographic analysis or numerical data, though the majority of publications described multiple calls per species or described multiple species per publication (493 calls from 72 species represented in 89 publications; summarized in Supplementary  Table 1). The databases used to search for these publications were Google Scholar, JSTOR, Web of Science, and Wiley Online Library. We used the keywords acoustics, acoustic repertoire, calls, frequency, Hz, vocalizations, and ultrasound paired with Sciuridae, squirrel, or an exhaustive list of currently valid and invalid genera (the most updated nomenclature was taken from the Integrated Taxonomic Information System 1 ). For each call described in the selected publication, the following characteristics were taken: the fundamental frequency (F 0 : the mean frequency of the primary vibrational frequency of the vocal membrane; kHz), dominant frequency (F Dom : the frequency with the greatest energy, power or amplitude; kHz), minimum frequency (F Min : the minimum frequency of the fundamental frequency; kHz), maximum frequency (F Max : the maximum frequency of the fundamental frequency or of harmonic on which F Dom is measured; kHz), and the highest visible harmonic (F Harm : mean frequency of the highest complete harmonic visible on the spectrograph; kHz) (Figure 1). Once our review of vocalization publications was complete, we searched for the body mass (g), diel activity pattern (diurnal or nocturnal), social complexity, and habitat openness of the dominant habitat (open or closed) of each species from the relevant vocalization papers. If not provided, other resources including Mammalian Species accounts, PanTHERIA (Jones et al., 2009), and the Animal Diversity Web (Myers et al., 2020) were reviewed. Both male and female body masses were initially recorded, but male body size could not be found for Spermophilus taurensis. Male and female body mass were strongly correlated (r = 0.98, p < 0.001), therefore female body mass was chosen to represent squirrel body size. Because we could only assign an adult female body mass to all species, calls that are exclusively produced by males or pups were removed from the dataset before analysis. We pooled all other calls (calls produced by both sexes or females only as well as calls produced by juveniles and adults) as there is little evidence to suggest that juveniles and adults produce acoustically distinct calls across the family (Matrosova et al., 2007(Matrosova et al., , 2011Volodina et al., 2010;Schneiderová, 2012;but see: Nikol'skii, 2007). While the initial database included a fivetiered social classification ranging from solitary to colonial [based on the social grades of ground squirrels described by Matějů et al. (2016)], social classes were reduced to social or solitary living to reduce model parameters. Species that exhibit dynamic social structures, such as flying squirrels that engage in social nesting to a greater extent during one portion of the year , were treated as socially living. Two subspecies (Marmota baibacina centralis and Tamias dorsalis dorsalis) could not be used in the subsequent analyses because ecological data and body mass-specific to each subspecies could not be found; similarly, the species Spermophilus pallidicauda could not be included as body mass for either sex could not be found.

Phylogeny
VertLife, an online resource that allows the user to extract pruned trees from vertebrate supertrees, was used to produce 100 pruned trees from the Mammalian supertree (Figure 2; Upham et al., 2019). Three subspecies had to be incorporated under their parent species, so branch tips were broken in two and subspecies were treated as equivalent to parent species, with branch lengths identical between the parent and subspecies (the addition of a subspecies did not create any polytomies in the tree). Three species are represented by subspecies only: Sciurus aberti kaibensis, Sciurus niger rufiventer, and Callosciurus erythraeus thaiwanensis.

Statistical Analysis
Phylogenetic generalized least square (PGLS) modeling was used to account for the variation in acoustic repertoire that may be explained by phylogenetic relatedness. PGLS models produce a lambda parameter, λ, that represents the degree to which the variance of traits is explained by the phylogenetic relationships in the model. The λ parameter varies between 0 and 1, with 0 representing no phylogenetic trace and 1 representing absolute Brownian motion (Freckleton et al., 2002;Martin et al., 2016).
PGLS modeling restricts each species to a single observation (i.e., no subsampling of species permitted). Therefore, the numerous data entries per species had to be reduced. For the fundamental, dominant, maximum, and highest harmonic frequencies, the absolute maximum value for each characteristic reported among all publications was chosen. Likewise, for minimum frequency, the absolute minimum reported frequency was chosen. We use maximum and minimum values rather than the median for a more rigorous test of our hypothesis about method limits.
Body mass and all frequency characteristics were logtransformed to achieve normal distributions. Additive models were built for each frequency type [β 0 + body mass (β Mass) + diel activity pattern (β Diel ) + sociality (β Sociality ) + habitat openness (β Open ) + method limits (β Lim )] using the caper package in R (ver 3.6.2). We reported the test statistics of the regression to evaluate significance and effect size (F-statistic, p-value, and adjusted R 2 ).

EVALUATING CORRELATES OF FREQUENCY CHARACTERISTICS
The effect of body mass and method limits were significant for all five frequency characteristics (F 0 , F Dom , F Min , F Max , and F Harm ) whereas habitat openness only had a significant effect for F Dom and F Harm ( Table 3). The PGLS models fit the data well for most frequency characteristics, with effect sizes (adjusted R 2 ) of 0.64 (p < 0.001), 0.57, (p < 0.001), 0.61 (p < 0.001), and 0.41 (p < 0.001) for F 0 , F Dom , F Max , and F Harm accordingly. The minimum frequency PGLS model did not fit the data well, yielding an effect size of only 0.12 (p < 0.05), however, it was the only model that yielded a significant phylogenetic signal (Pagel's λ with 95% confidence intervals) of 0.65 [0, 0.88]. In contrast, the phylogenetic signal was estimated to be 0 for  Table 2). Notably, harmonic frequency values could only be found for 43 of the 73 species, with species reported in only one publication being significantly less likely to have reported harmonics (Chi-square = 17.1, p < 0.001) than those appearing in multiple publications.

Body Mass
Body mass (Figure 3) had a significant effect on F 0 (slope = −0.22 (±SE = 0.05; Table 4), F 1,58 = 23.14, p < 0.001), F Dom (slope = −0.20 (±0.05), F 1,60 = 12.37, p < 0.001), F Min (slope = −0.19 (±0.09), F 1,58 = 6.24, p < 0.05), F Max (slope = −0.25 (±0.05), F 1,59 = 37.61, p < 0.001), and F Harm (slope = −0.29 (±0.08), F 1,37 = 9.18, p < 0.01). These results are consistent with our prediction that increased body mass results in the production of lower frequencies as already shown in mammals (Charlton and Reby, 2016) and birds (Ryan and Brenowitz, 1985). The avian mass-frequency relationship focused on the "emphasized" frequency, similar to F Dom reported here (Ryan and Brenowitz, 1985). Our results are somewhat consistent with Martin et al. (2016) who also found that F Min was constrained by mammalian body mass, but the mass-F Min relationship reported was twice as strong as reported here (slope = 0.41). This stronger relationship likely arose because their dataset ranged from microbats (weighing less than 0.1 kg) to whales (larger than 100,000 kg), representing a much larger range of body sizes and therefore showing a stronger effect of body size on frequency types. In comparison, this dataset covered only a range of about 36 g (Tamias alpinus) to 8.0 kg (Marmota sibirica). In contrast to our results, Martin et al. (2016) found that F max was not constrained by body size, though the lack of relationship was largely attributed to the interaction between environment and body size. They suggest that because sound waves travel differently in the water, aquatic species are freed from the constraints of body size on the maximum frequency compared to their terrestrial relatives. Like Martin et al. (2016), Ryan and Brenowitz (1985) also found that there was an interaction between environment and body size in avian song frequency, but FIGURE 2 | Phylogeny of Sciuridae reduced to all species used in analysis (n = 70). Species present represent the subfamily Callosciurinae (gray: southern Asian tree squirrels), Sciurinae (orange: Pteromyini, flying squirrels; yellow: Sciurini, Holarctic tree squirrels tree squirrels) and Xerinae (pink: tribe Xerini, African ground squirrels; Marmotini, Holarctic ground squirrels). Marmotini is subdivided into the genera Tamias (green: chipmunks), Marmota (light blue: marmots) and the remaining Holarctic ground squirrels (dark blue). Edge lengths to scale; scale bar represents 5 million years from present. Tree downloaded from Vertlife.org. Attribution details for the animal images available in Supplementary Table 3.  Notably, sociality had a significant effect on these models while yielding a non-significant slope; investigating the models further, there was no evidence of collinearity between sociality and other variables. However, there was a larger variation in frequencies used by social squirrels compared to solitary squirrels. These results are consistent with Blumstein and Armitage (1997b) who found that social complexity was related to increased use of alarm calls in the repertoire in marmots. In our results, we found that alarm calls are widespread across the taxa, appearing in 64 of the 73 vocalizing Sciuridae species. While social squirrels might not be more likely to have alarm calls in their vocal repertoire than solitary squirrels (Chi-Square = 3.51, p = 0.06), they likely maintain shorter distances between the sender and receiver and increase their production of alarm calls (e.g., Urocitellus beldingi; Sherman, 1985). This reduced distance could facilitate the incorporation of higher frequencies into the repertoire as these frequencies are highly directional but are more likely to be attenuated over longer distances (Smith, 1979). As well, the increased proximity not only results in more altruistic alarm calling (Sherman, 1985) but also the production of more vocalizations to facilitate group cohesion and maintain social hierarchies, particularly between emerging juveniles and established adults. For example, calls in adult female yellowbellied marmots (Marmota flaviventris) are directly correlated to the emergence of juveniles (Blumstein et al., 1997). Primates, another group of highly communicative and socially variable mammals, also trend toward increased vocal complexity and overall improved hearing sensitivities in more social species and larger groups (Ramsier et al., 2012). To compensate for increased exposure via sound, social animals benefit from communicating in ranges less readily heard by predatory species (Arch and Narins, 2008). In contrast, solitary squirrels would benefit from louder calls if they are warning surrounding conspecifics of potential danger, especially if the species holds relatively large home ranges. However, the maintenance of these territories would likely lead to vocal production that maximizes the effectiveness of territorial displays (which we only found reports of in solitary and socially alternating species), meaning that lower frequencies would be more beneficial to maximize amplitude over longer distances and to increase omnidirectionality. In addition Significant estimates are given in bold. *p < 0.05; **p < 0.01; ***p < 0.0001.
to social structure, mating systems may be complicating the results observed. We grouped monogamous and polygynous species. However, Charlton and Reby (2016) suggests that F 0 is a sexually selected trait and is higher in polygynous species and lower in monogamous species. While we excluded calls produced exclusively by males, we kept calls produced by both sexes, therefore not eliminating calls produced by both sexes that are also under sexual selection. We suggest further research into the role of sociality and mating systems in vocal complexity across Sciuridae.

Habitat
The  (Smith, 1979) resulting in a shift of vocal range (minimum to maximum frequency). Instead, we found that squirrels in open habitats spend their maximum acoustic energy (F Dom ) on higher frequencies than squirrels in closed habitats, with little effect on range, as F Min and F Max were not significantly influenced by habitat openness (Tables 3, 4). Charlton et al. (2019) have shown that forested animals may not be hindered by their ability to produce or detect high frequencies as many have been shown to co-evolve improved acoustic structures to not lose their access to these frequencies in a closed environment. They suggest that forested animals may boost the amplitude and increase the frequency band of their calls to help counteract the increased attenuation of their environment and improve the ability of receivers to locate the sender (Charlton et al., 2019). This hypothesis would explain why squirrels in closed habitats continue to use high frequencies, despite the limitations. Based on this hypothesis, we would also predict F Harm would not be affected by habitat openness as increased harmonic strength would improve the bandwidth of calls in closed environments and high harmonics would be unhindered in open habitats, making them favorable in both conditions. We found instead that higher harmonics are detected in open habitats, supporting the more traditional acoustic adaptation hypothesis (Blumstein, 2007;Ey and Fischer, 2009). However, we think it is important to note that the detection of the highest harmonic is largely dependent on the recording techniques of the researcher and how the calls have been displayed on spectrographs. It may likely be that high harmonics are easier to record in open environments, whether because the acoustics are better for recording (reduced attenuation of harmonics before reaching the microphone as in bats; Fenton et al., 2011) or the individual squirrel ecology lends to better study subjects (i.e., congregated ground squirrels allow microphones to be placed close to predictable burrow holes compared to free-ranging tree squirrels which have many temporary nests and retreats across a forest canopy). While not explored in this study, environments with greater humidity transmit higher frequencies better than arid environments where calls are attenuated faster (Nikol'skii, 1984). The importance of habitat ariditity was demonstrated in ground squirrels (Spermophilus; Nikol'skii, 2017) when the strength of the relationship between the F Dom and body mass was significantly improved when habitat aridity was incorporated as a factor into the models. Overall, our findings suggest that higher dominant frequencies are more easily incorporated and recorded in squirrels in open habitats, but high frequencies are not necessarily unused or completely inhibited by closed habitats. all frequency characteristics trend toward nocturnal species using higher frequencies than diurnal species (Figure 4). Unfortunately, the lack of behavioral observations associated with the published flying squirrel repertoires, the only nocturnal squirrels, means that interpretation of these results is mostly hypothetical (as nocturnal patterning is only represented by 5 different (sub)species). Both species of Asian giant flying squirrel (Petaurista lena and P. leucogenys) stand out as the only squirrels with movement-associated vocalizations that were recorded while observing free-ranging individuals (Ando and Kuramochi, 2008;Shen, 2013). Additionally, one call associated with the exploration of a novel environment was reported in the southern flying squirrel (Glaucomys volans; Murrant et al., 2013). While there has yet to be any direct evidence of echobased navigation (Siemers et al., 2009) or echolocation in flying squirrels (Chattin, 1969), we propose that the strong effect of nocturnality on multiple frequency characteristics, correlated to calls in higher frequencies, as well as the recent discovery of USVs in North American flying squirrels (Gilley, 2013;Murrant et al., 2013), warrants a re-examination of the function of USVs in these squirrels with specific echo-based navigation experiments included in the research.

Method Limits
As we predicted, the limits of the methods (microphone or analysis ranges) did significantly influence all the frequency characteristics as shown in the PGLS models. . We conclude that the type of equipment being used is crucial for detecting all frequency characteristics and that using equipment with large bandwidths (the entire range of frequencies that could be detected) is important for all squirrel acoustics research as we continue to discover novel USV production in Sciuridae.

Summary
We predicted that body mass, sociality, habitat openness, diel activity pattern, and method limits would influence the  fundamental, dominant, minimum, maximum, and highest harmonic frequencies (kHz ; Table 1). We found at least some support for all of these hypotheses. We found that across the squirrel family, body mass significantly affected all frequency characteristics, a finding that is mostly consistent with other literature (Arch and Narins, 2008;Charlton and Reby, 2016;Martin et al., 2016). Social squirrels used higher fundamental, dominant, and maximum frequencies which aligns with the predator-avoidance hypothesis (Arch and Narins, 2008). Social squirrels also used a larger range of frequencies which can be partially attributed to increased communication amongst kin (Blumstein and Armitage, 1997b) or sexual selection acting on vocal traits (Charlton and Reby, 2016). As predicted by the acoustic adaptation hypothesis (Blumstein, 2007;Ey and Fischer, 2009), we found that species inhabiting open habitats used higher dominant and harmonic frequencies than those in closed habitats where high frequencies are more rapidly attenuated by the environment (Smith, 1979); however, the maximum frequency was not affected by habitat type which suggests that squirrels in closed habitats may have co-evolved structures to still maintain access to USVs despite the environmental complications (Charlton et al., 2019). Nocturnality was associated with higher frequencies for the fundamental, dominant, and maximum frequencies, suggesting that nocturnal, gliding squirrels are under unique selection pressures for using USVs. An explanation for this general pattern of higher frequencies in nocturnal squirrels remains elusive. Finally, the significant effect of method limits on all frequency characteristics and the presence of method limits in all candidate models (Supplementary Table 2) suggest that method limits are important for framing the ranges detected during recording sessions.

THINK BEFORE THEY SQUEAK -CONCLUSION
Squirrels are vocally complex taxa, displaying a wide range of call frequencies and functions. While there has been a representative sampling across the genera, there are large gaps in the geographic sampling of squirrels ( Figure 5). The United States alone accounts for over 50% of the published vocal repertoires. In contrast, only a single study has occurred in Africa and no repertoires have been published from Central or South America. In the face of global extinction crises and habitat loss, it is more important than ever to archive the traits of rare and cryptic species so that we can maximize the utility of all available monitoring and conservation tools. Bioacoustic monitoring is an effective tool for cryptic species (Walker, 1964;Gilbert et al., 1994;Tripp and Otter, 2006) and has been shown to improve the monitoring of flying squirrels which are notoriously difficult to detect (Diggins et al., 2016). Considering these analyses, we encourage researches to record squirrel vocalizations in under-sampled locations. We further suggest the use of recording equipment that can resolve frequencies that extend into the ultrasonic range so that repertoires accurately reflect the ranges used by the species being studied. Additionally, we caution that researchers studying the vocalizations of previously reported squirrels should consider the accuracy of previous literature, especially when there is possible evidence of partial USVs or ultrasonic harmonics present. Squirrels that are socially complex or inhabiting open areas are of interest for studies focused on high-frequency vocalizations; we also strongly encourage the examination of the role of USVs in nocturnal, flying squirrels. We caution that while USV research requires specialized technology, this work also requires ideal recording conditions and patience as USVs can be rare and difficult to detect in some Sciuridae species (Matrosova et al., 2012). In conclusion, while rodent bioacoustic studies are typified by controlled environments and lab animals, studying the vocalizations of free-ranging rodents can reveal new information about the natural world, providing informative findings and new insights that apply to real-time conservation efforts.

DATA AVAILABILITY STATEMENT
The data and R scripts that support the findings of this study are openly available on Dryad at https://doi.org/10.5061/dryad. vt4b8gtpm.

AUTHOR CONTRIBUTIONS
JB and SN conceptualized the manuscript. SN collected data, conducted analyses, and prepared the manuscript with guidance and critical feedback from JB. All authors contributed to the article and approved the submitted version.

FUNDING
The authors' institutions provided support for this publication.

ACKNOWLEDGMENTS
We thank Trent University for access to online databases, P. O'Brien for assistance in article collection and valuable feedback during conceptualization, and PO'B, K. Solmundson, and R. Dillion for their insightful comments during the final editing stages.