- Open Access
The power of oral and nasal calls to discriminate individual mothers and offspring in red deer, Cervus elaphus
Frontiers in Zoologyvolume 12, Article number: 2 (2015)
In most species, acoustical cues are crucial for mother-offspring recognition. Studies of a few species of ungulates showed that potential for individual recognition may differ between nasal and oral contact calls.
Vocalizations of 28 hinds and 31 calves of farmed Iberian red deer (Cervus elaphus hispanicus) were examined with discriminant function analyses (DFA) to determine whether acoustic structure of their oral and nasal contact calls encodes information about the caller’s identity. Contact calls were elicited by brief separation of individually identified animals by a distance over 10 m or by a bar fence. Both oral and nasal calls of both hinds and calves showed high potential to discriminate individuals. In hinds, individuality was significantly higher in the oral than in the nasal calls, whereas in calves, individuality was equally well expressed in both oral and nasal calls. For calves, the maximum fundamental frequency was higher and the duration was longer in oral calls than in nasal calls. For hinds, the maximum fundamental frequency and the duration were indistinguishable between oral and nasal calls. Compared to the pooled sample of oral and nasal calls, separate oral or nasal call samples provided better classifying accuracy to individual in either hinds or calves. Nevertheless, in both hinds and calves, even in the pooled sample of oral and nasal calls, the degree of individual identity was 2–3 times greater than expected by chance. For hinds that provided calls in both years, cross-validation of calls collected in 2012 with discriminant functions created with calls from 2011 showed a strong decrease of classifying accuracy to individual.
These results suggest different potentials of nasal and oral calls to allow the discrimination of individuals among hinds, but not among red deer calves. The high potential of individual recognition even with the pooled sample of oral and nasal calls allows mother and young to remember only one set of acoustic variables for mutual vocal recognition. Poor between-year stability of individual characteristics of hind oral and nasal calls would require updating keys to individual recognition each calving season.
Mother-offspring recognition is critically important for the survival of the young for many taxa: ungulates [1,2], pinnipedes [3-5], bats [6,7], penguins [8,9], alcids [10,11], larids [12,13] and cranes . Along with the visual and olfactory cues vocalizations play an important role in this process [15-18]. Visual and olfactory cues are generally perceived across a limited area relative to acoustic cues , which can also be perceived at night, and through dense vegetation . Vocal recognition is important for selective feeding one’s own offspring [1,3,8,19], protection against predators that are dangerous for the young but not for the mother [21-23], and for maintaining of spatial proximity between mother and young in herds and flocks .
Vocal recognition between mother and young is based on individual features of calls. In mammals, the calls are generated by vibrations of vocal folds in the larynx. After that, the sound passes through the vocal tract, which filters the signal, accenting the resonant frequencies of the vocal tract, called formant frequencies [25,26]. So, the vocal output is the joint product of the vocal fold vibrations, determining the call fundamental frequency (f0), and of the work of the vocal tract, determining the values of formant frequencies. According to the source-filter theory, source and filter variables are independent of each other [26,27]. Vocal individual identity can be encoded by call source variables and filter variables. In ungulates, the fundamental frequency has primarily contributed to vocal identity, as exemplified by mother and young fallow deer Dama dama , red deer Cervus elaphus , domestic sheep Ovis aries ; goats Capra hircus [17,18], in the young of reindeer Rangifer tarandus , white-tailed deer Odocoileus virginianus , mule deer O. hemionus , and goitred gazelles Gazella subgutturosa [31,32]. Formant frequencies play a role in encoding individual identity in the young of goitred gazelles  and reindeer , in red deer hinds , and in mother and young domestic goats . Call duration proved significant in encoding individual identity in mother and young red deer  and domestic sheep .
Encoding individuality with acoustic variables may be complicated due to variable modes of vocal production. Mother and young ungulates produce two types of contact calls: oral calls that are made through a widely opened mouth, and nasal calls that are made through the nose with a closed mouth. The oral and nasal modes of vocal production were previously reported for the young of white-tailed deer  and goitred gazelles [31,34], for mother domestic sheep , for mother and young saiga anthelopes Saiga tatarica , and for mother and young red deer . The acoustic structure of the oral and nasal calls is substantially different [31,35-37]. Oral and nasal calls can be produced by callers in the same series [36,37]. Although it is commonly considered that the oral calls are emitted in situations of higher tension compared to the nasal calls, the communicative potential of these two different call types remains unclear [31,35,36]. In particular, it is unclear whether oral and nasal calls share common cues identifying individuals across types, or hinds and calves must remember two unique sets of cues for encoding individuality specifically in the oral and nasal calls.
Individual vocal traits can be used for recognition only if individual differences remain stable over time. In some species, vocalizations are stable over time [38-40], and mothers can recognize their offspring calls [3,41]. In contrast in others, individual vocal traits change over a short time [42-45]. Although stability of acoustic individuality of offspring calls has been studied , to date, no information is available on the stability of mother vocal signature in mammals.
In red deer, individual features were reported for rutting roars of farmed and free-ranging stags [46-48] and for contact calls of mother and offspring . Playbacks showed that calves respond stronger to contact calls of their own mothers compared to alien mothers, whereas hinds respond stronger to contact calls of their offspring compared to calls of alien calves . In a previous study, we described the structure of oral and nasal calls of hinds and calves of Iberian red deer C.e. hispanicus during the rut period; however, callers were not individually identified .
This study investigates the structure of the oral and nasal calls of individually identified mother and young Iberian red deer in calving period. We compared the potential of variables of the oral and nasal calls in both mother and young (Figure 1) to encode individual identity. We also tested whether body mass and sex affected the acoustics of calls in calves. We estimate which acoustic variables primarily contribute to vocal identity in the oral and in the nasal calls and across call types. In addition, for hinds, we estimate the between-year stability of individual identity in their oral and nasal contact calls.
Comparison of oral and nasal calls
For hinds, a repeated measures ANOVA controlling for individuality revealed no significant differences in the duration of oral and nasal calls, but revealed a significantly shorter dur-to-max in the oral compared to the nasal calls (Table 1). None of the variables associated with f0 differed between oral and nasal calls with the exception of f0beg, which was significantly higher in oral calls. The fpeak, q25, power-f0 and peak-harm were also significantly higher in the oral than in the nasal calls (Table 1).
For calves, a repeated measures ANOVA controlling for individuality revealed significantly higher values in all measured variables in the oral calls compared to the nasal calls with the exception of dur-to-max (Table 1). Unlike hinds, in calves the oral calls were longer and had higher values of f0 variables compared to the nasal calls. Overall, oral and nasal calls differed by structure more strongly in calves than in hinds.
Individual discrimination with DFA
For the same sample of 22 hinds, we estimated the values of correct classification to individual for oral calls, for nasal calls, and for the pooled sample of oral and nasal calls (Figure 2). In all three DFAs, the average value of correct assignment to individual (77.0% for oral calls, 61.8% for nasal calls, 53.9% for the pooled sample of oral and nasal calls) significantly exceeded our random expectation (24.1 ± 2.5%, 23.5 ± 2.6% and 15.5 ± 1.6% respectively, all p < 0.001). The average value of correct assignment to individual was higher in oral than in nasal calls (χ 2 1 = 10.80, p = 0.001) and in oral calls than in calls of the pooled sample of oral and nasal calls (χ 2 1 = 34.57, p < 0.001), although it did not differ significantly between nasal calls and the pooled sample of oral and nasal calls (χ 2 1 = 3.24, p = 0.07) (Figure 2).
Because the same variables were included in a DFA for both oral and nasal calls in the same sample of 22 hinds, we could compare percentages of correct assignment for nasal and oral calls to particular individuals with repeated measures ANOVA. The ANOVA revealed significantly higher percentages of correct assignment to individual for oral calls compared to nasal calls (F 1,21 = 10.76; p = 0.004) and for oral calls compared to the pooled sample of oral and nasal calls (F 1,21 = 28.63; p < 0.001), although it did not reveal statistical differences in classifying accuracy between nasal calls and the pooled sample of oral and nasal calls (F 1,21 = 4.03; p = 0.06).
For hind oral calls, the f0max, f0end and f0mean (in order of decreasing importance) were mainly responsible for discrimination of individuals (Table 2). For hind nasal calls, the duration, f0beg and f0mean (in order of decreasing importance) were mainly responsible for discrimination of individuals (Table 2). For hind pooled sample of oral and nasal calls, the f0max, f0mean and duration (in order of decreasing importance) were mainly responsible for discrimination of individuals (Table 2). Therefore, in all three DFAs for hinds, similar sets of cue discriminating variables were found.
For the same sample of 17 calves, we estimated the values of correct classification to individual for oral calls, for nasal calls, and for the pooled sample of oral and nasal calls (Figure 3). In all the three DFAs, the average value of correct assignment to individual (61.1% for oral calls, 64.2% for nasal calls, 49.5% for the pooled sample of oral and nasal calls) exceeded significantly the random value (respectively 30.8 ± 3.5%, 32.2 ± 3.6%, 20.3 ± 2.0%), calculated with randomization procedure (all differences are significant, p < 0.001). Unlike hinds, in calves the average value of correct assignment to individual did not differ between oral and nasal calls (χ 2 1 = 0.17, p = 0.68) and was significantly higher either in oral or in nasal calls compared to the pooled sample of oral and nasal calls (χ 2 1 = 4.83, p = 0.03 and χ 2 1 = 7.34, p = 0.007 respectively) (Figure 3).
Because the same variables were included in a DFA for both oral and nasal calls in the same sample of 17 calves, we could compare percentages of correct assignment for nasal and oral calls to particular individuals with repeated measures ANOVA. The ANOVA did not reveal significant differences in percentages of correct assignment to individual between oral and nasal calls (F 1,16 = 0.08; p = 0.78). For the pooled sample of oral and nasal calls, the percentages of correct assignment to individual were significantly lower than with unpooled samples of either oral (F 1,16 = 5.09; p = 0.04) or nasal calls (F 1,16 = 5.82; p = 0.03).
For calf oral calls, the f0beg, duration and power-f0 (in order of decreasing importance) were primarily responsible for discrimination of individuals (Table 3). For calf nasal calls and for the pooled sample of oral and nasal calls, the f0beg, duration and q75 (in order of decreasing importance) were primarily responsible for discrimination of individuals (Table 3). Thus, in all three DFAs for calves, the sets of variables primarily responsible for discriminating individuals were similar.
Body mass of 22 hinds whose calls were included in DFA was mean ± SD = 104.8 ± 12.2 kg and ranged from 85.3 to 121.5 kg, and body mass of 17 calves whose calls were included in DFA was 16.0 ± 4.4 kg and ranged from 10.0 to 25.8 kg. To approximate probable effects of body mass on the acoustic variables, for hinds we calculated correlations between log3 body mass and five variables (durat, f0beg, f0end, f0max, f0mean, see Table 2) primarily contributing to discrimination of individuals by calls. We did not find effects of body mass on variables of either oral calls (durat: r = 0.12, p =0.61; f0beg: r = −0.20, p =0.38; f0end: r = −0.27, p =0.24; f0max: r = −0.13, p =0.55; f0mean: r = −0.23, p = 0.31) or nasal calls (durat: r = 0.08, p = 0.74; f0beg: r = −0.26, p = 0.24; f0end: r = −0.20, p = 0.38; f0max: r = −0.14, p = 0.55; f0mean: r = −0.21, p = 0.34).
For calves, we estimated the effects of sex, body mass and the conjoint effect of sex*body mass on four acoustic variables (durat, f0beg, q75, power-f0, see Table 3) primarily contributing to discrimination of individuals by calls. ANCOVA with sex as a fixed categorical factor and body mass as a continuous factor revealed significant effects of body mass on the duration, f0beg and q75 of the oral calls and on the duration of the nasal calls (Table 4). For all these acoustic variables the values decrease with increase of body mass. Only for the nasal calls ANCOVA revealed significant effects of sex and sex*body mass on the duration and q75 (Table 4), which values were lower in male than in female calves (duration: males 0.227 ± 0.019 s, females 0.232 ± 0.031 s; q75: males 3588 ± 258 Hz, females 3642 ± 243 Hz). At the same time, body mass did not differ significantly (one-way ANOVA, F 1,16 = 1.15, p = 0.30) between male and female calves (males 17.0 ± 4.9 kg, females 14.7 ± 3.7 kg).
Between-year stability of hind calls
For hinds that provided sufficient number of calls in both study years (2011 and 2012) we compared the stability of vocal individual identity between years (Figure 4). Within years, individuality was very high in either oral or nasal calls. Among oral calls, the average values of correct classification to individual (93.8% in 2011 and 86.0% in 2012) significantly exceeded the random expectation (49.9 ± 5.5% in 2011 and 52.0 ± 6.2% in 2012, p < 0.001 in both cases) and did not differ between years (χ 2 1 = 1.33, p = 0.25). Among nasal calls, the average values of correct classification (69.6% in 2011 and 68.6% in 2012) significantly exceeded the random expectation (33.6 ± 3.6% in 2011 and 34.7 ± 3.8% in 2012, p < 0.001 in both cases) and did not differ between years (χ 2 1 = 0, p = 0.97).
However, cross-validation of calls recorded in 2012 using discriminant functions created for calls recorded in 2011 revealed a strong decrease in the correct classification of individuals with either oral or nasal calls (Figure 4). For the oral calls, the average value of correct classification dropped to the level expected by chance alone (56.1%), and became significantly lower compared to call samples from either 2011 (χ 2 1 = 21.79, p < 0.001) or from 2012 (χ 2 1 = 10.92, p = 0.001). Similarly, for the nasal calls, the value of correct classification dropped to the level expected by chance alone (36.4%) and became significantly lower compared to call samples from either 2011 (χ 2 1 = 25.96, p < 0.001) or from 2012 (χ 2 1 = 23.93, p < 0.001). Thus, in hinds, individual identity of oral and nasal calls was unstable between years.
The acoustics of oral and nasal calls
Our study revealed that among calves, the fundamental frequency was higher in oral than in nasal calls. The same relationship was previously found in calls of juvenile goitred gazelles , in rumbles of adolescent African elephants Loxodonta africana , in calls of adult female domestic sheep , and in calls of newborn and adult female saiga antelopes . The higher fundamental frequency of oral than nasal calls may represent a common rule, in that oral calls are made at higher arousal. The anatomical and productional basis of this difference in fundamental frequency between oral and nasal calls in mammals has been discussed in detail in previous studies [31,36]. Nevertheless, in Iberian red deer hinds in this study, the fundamental frequency did not differ between oral and nasal calls. These results contradict those obtained for other ungulates [31,35,36] and confirm our previous data for unmarked Iberian red deer hinds during the rut .
The values of fundamental frequency variables of Iberian red deer calf oral calls in our study (f0max = 875 Hz, Table 1) were similar to those reported for oral distress calls of neonate Middle-European red deer C.e. hippelaphus (f0max = 844 Hz ), and were slightly higher in frequency than distress calls of 2–4 day-old Middle-European red deer calves (f0max = 737 Hz ) and contacts calls of 4 month-old Corsican red deer calves C.e. corsicanus (f0max = 710 Hz ). At the same time, power quartiles of calf oral calls in our study (q25 = 952 Hz, q50 = 1731 Hz, q75 = 2471 Hz, Table 1) were substantially lower than on  (q25 = 2385 Hz, q50 = 4669 Hz, q75 = 7398 Hz). These differences could be attributable to degree of emotional arousal of a caller, as we collected contact calls from calves that were separated from their mothers, whereas Teichroeb et al.  collected distress calls made during captures, ear tagging and human handling. The energy shifts towards higher frequencies (and respective increase of power quartile values) represent a common rule of acoustics changes attending increased emotional arousal of mammalian callers [52-54].
Compared to calls of 3–4 month-old calves recorded on the same farm during the rut , calls during the calving period were shorter and higher-pitched, which is a product of decreasing of fundamental frequency with age. Hind calls were higher in fundamental frequency during the rut relative to the calving period, which may result from the lower disturbance of animals in the current study because of the longer period of data collection. In both the rut  and calving period (this study) we found a higher fundamental frequency in oral than in nasal calf calls and did not observe significant differences between oral and nasal calls of hinds. Our study also confirmed the findings of previous studies of Iberian red deer, that revealed that the fundamental frequency of hind contact calls is lower than that in rutting roars of stags [37,55].
Individuality of oral and nasal calls
In hinds, individual identity was more pronounced in oral than in nasal calls. The same relationship was previously reported for calls of juvenile goitred gazelles . One possible explanation for the higher potential of oral calls to allow deer to discriminate different individuals is that such calls are produced in situations of higher tension than nasal calls [31,35]. Among Iberian red deer calves in our study, however, the expression of individual identity detected with DFA, did not differ between oral and nasal calls. No comparative data is available to explain this discrepancy.
In hinds and in calves, classification accuracy of callers was two to three times greater than the random expected value. This was observed even for the pooled sample of oral and nasal calls, probably because the same acoustic variables were responsible for individuality of either oral or nasal calls. This should allow mothers and their young to remember only a single set of acoustic variables for mutual vocal recognition instead of remembering two different sets for the oral and for the nasal calls. At the same time, variables mostly contributing to discrimination differed between mother and young. In mothers, discriminating power of oral and nasal calls was mainly based on variables related to fundamental frequency, whereas in calves, it included f0beg along to the duration and power-f0, that is, was based on variables of all modalities, temporal, frequency and power. To reveal the communicative functions of nasal and oral calls, further playback study would be necessary. Consistent with our demonstration of individually-distinct contact calls among hinds and calves, orally produced rutting roars of red deer stags also have a potential to encode individual identity. In two studies, the degree of correct classification of individual stags by their rutting roars was about 70% for a sample of 5–6 stags, in spite of different sets of measured acoustic variables [47,48].
Estimating vocal individuality was complicated by differences in age ranges among calves during the data collection period. Call parameters, in particular the fundamental frequency, depend on the caller’s body size among ungulates [34,56], although in adult red deer this dependence is lacking [37,57]. The fundamental frequency parameters contributed substantially to classifying individuals with DFA in both calves and hinds. Therefore, to address this issue, we repeatedly recorded individuals across the period of data collection. Nevertheless, we found significant effects of body mass on the start fundamental frequency and duration in the oral calls and on duration in the nasal calls of calves. Thus, we could not totally exclude the effects of the body mass on individual distinctiveness of calf contact calls in this study. Nevertheless, the vocal individuality of calves was comparable with that of hinds (in which effects of body mass on call variables were not found) and with values in other ungulates, e.g. 20 goitred gazelle calves, with correctly classified to individual 75% oral calls and 65% nasal calls .
Between-year stability of hind calls
Individual identity of hind oral and nasal calls was unstable between years. Thus, if the ability to discriminate among individuals based on calls alone is to persist across years, cues to individual identity in hind oral and nasal calls must develop anew each year. Degradation of vocal individuality with time has been demonstrated for common roars of red deer stags  and for groans of fallow deer bucks . These results support the hypothesis that vocal individuality does not result exclusively from individual vocal anatomy and that there is some vocal tuning and vocal flexibility in vocalizations of ungulates [32,58,59].
Materials and methods
Study site, subjects and dates of recordings
Calls of hinds and calves were recorded in June 2011 and 2012 at the farm of the University of Castilla-La Mancha (Albacete, Spain, 38°57’10”N, 1°47’00”W, 690 m a.s.l.). The population originated in 1994 from 15 male and 50 female pure Iberian deer from a nearby Las Dehesas public game reserve in Alpera (Albacete) and from Cabaneros National Park (Toledo). The animals used in this study were born and kept in four 10,000 m2 enclosures on an irrigated pasture. They were fed ad libitum with a diet of barley straw and meal from barley, alfalfa, oats and sugar beets .
All mothers and calves were captive-born and kept together in permanent groups (4 groups in 2011 and 3 groups in 2012) separately from adult stags and yearlings. The groups ranged in size from 6 hinds and 2 calves to 30 hinds and 24 calves (mean ± SD = 17.4 ± 7.6 hinds and 15.0 ± 7.7 calves per group). The entire population of animals from which we collected calls counted 61 hinds and 52 calves in 2011 and 61 hinds and 53 calves in 2012 (45 hinds were the same in both years). The age of mothers was 2–19 years (mean ± SD = 10.0 ± 4.6 years). The calves were born from 6 May to 23 June in 2011 and from 13 May to 14 June in 2012. During data collection, the age of calves varied of 1–52 days. All study calves were singletons, excluding two sibling pairs, one in 2011 and one in 2012.
All hinds and calves were individually labeled with Allflex (Palmerston North, New Zealand) plastic ear tags and all hinds were additionally marked with the Allflex color collars with numbers. Three times a day, the farm staff inspected each enclosure with hinds and calves for newborns and labeled them individually with an ear tag, which allowed us to establish a mother for each calf.
Call and body mass collection
For acoustic recordings of hind and calf (48 kHz, 16 bit), we used a solid state recorder Marantz PMD-660 (D&M Professional, Kanagawa, Japan) with a Sennheiser K6-ME66 cardioid electret condenser microphone (Sennheiser electronic, Wedemark, Germany). The distance from the hand-held microphone to the animals was 5–35 m, the level of recording was adjusted during the recordings accordingly to the intensity of the produced calls.
We recorded calls daily, from 6:00–7:00 to 12:00–13:00, often with synchronous video, using a digital camcorder Panasonic HDC-HS100 (Panasonic Corp., Kadoma, Japan). During recordings, individual identities of callers producing calls through the mouth and through the nose were labeled by voice. Recordings have been conducted both inside and outside the outdoor enclosures during different contexts: everyday routine activity, when mothers searched for their offspring which were hidden in the enclosures; at translocations to small paddocks, from where the animals were taken for weightings; during temporal separations of hinds and calves after the weighing and at short separations by presence of researchers in enclosure on the way between a mother and her calf. In total, in 2011 and 2012, we collected 30 hours of audio recordings (16 hours in 2011 and 14 hours in 2012) from 28 individual hinds (9 hinds were recorded in both 2011 and 2012) and from 31 individual calves. All animals were weighted one time with Mettler-Toledo ID1 scales (Mettler-Toledo S.A.E., Barcelona, Spain) as the part of routine farm management  during the periods of acoustic recordings. The age of calves at weightings ranged from 1 to 40 days.
For acoustic analyses, we took only calls of good quality with high signal-to-noise ratios sufficient for analysis of all acoustic variables, measured in this study that were not disrupted by wind or the calls of other animals or overloaded during the recordings. We analysed only individually identified calls of known call type (nasal or oral). Calls were classified to nasal and oral call types based on voice comments of researchers made during recording, by video clips, made synchronously with the recordings, by the obvious nasal quality of sound within a recording (Additional files 1 and 2) and by the difference in call energy distribution, shifted towards higher frequencies in oral calls due to the shortening the vocal tract at opening of the mouth. These methods of classification to nasal and oral call types were previously applied for the Iberian red deer , for goitred gazelles [31,32] and for saiga antelopes . Two researchers (OS and IV) independently classified all calls, and we took for analysis only calls where both researchers were concordant in their judgments concerning their type. To reduce pseudoreplication, we took calls from different recording sessions per animal and from different parts within session, because calls from the same sequence are commonly more similar in their acoustic structure than calls from different sequences . The mean ± SD number of sessions per animal was 5.3 ± 4.4, and we took up to 15 high-quality calls per individual per call type (nasal and oral) from 28 mothers and 31 young (16 males, 15 females). We analysed 801 calls of mothers (354 oral and 447 nasal) and 469 calls of calves (281 oral and 188 nasal), for 1270 calls in total.
Acoustic analyses were conducted in the same way for hinds and calves and for the oral and nasal calls. For each nasal and oral call, we measured the same set of 14 acoustic variables: 2 temporal, 6 variables of fundamental frequency (f0) and 6 power variables, as they proved their use for estimating vocal individual identity in red deer [28,47,48]. Before analysis, the calls were downsampled to 11025 Hz and high-pass filtered at 50 Hz, to increase frequency resolution and to reduce the low-frequency background noise. We measured the duration of each call and the duration from call onset to the point of maximum f0 (dur-to-max) manually on the screen with the reticule cursor in the spectrogram window (Hamming window, FFT 1024 points, frame 50% and overlap 96.87%) by using Avisoft SASLab Pro software (Avisoft Bioacoustics, Berlin, Germany). Then we performed manual measurements on the screen with the standard marker cursor of the start (f0beg), maximum (f0max) and end (f0end) fundamental frequencies of each call (Figure 1). Measurements were exported automatically to Microsoft Excel (Microsoft Corp., Redmond, WA, USA).
In a 50 ms call fragment symmetrical about f0 maximum, we created the power spectrum, from which we automatically measured fpeak, representing the value of the frequency of maximum amplitude and the q25, q50 and q75, representing the lower, medium and upper quartiles, covering 25%, 50% and 75% of the energy of the call spectrum respectively (Figure 1). On the same spectrum, we estimated (in dB) the power-f0, representing the relative power of the f0 band compared to the peak harmonic, on the screen using two harmonic cursors (Figure 1). The power-f0 was equal to 0 when the f0 band coincided with the fpeak band. In addition, we recorded the peak-harm, representing the order number of the harmonic with the maximum energy.
In addition, we measured the f0 variables following Reby and McComb  by using the Praat DSP package . The f0 contour was extracted by using a cross-correlation algorithm (to Pitch (cc) command in Praat). The time steps in the analysis were 0.01 s for hinds and 0.005 s for calves; the lower and upper limits of the f0 range were 50–400 Hz for hinds and 100–1200 Hz for calves. A preliminary visual analysis of the spectrograms in Avisoft showed that the lower limit was lower than the minimum f0 for calls of either hinds or calves. Spurious values and octave jumps in the f0 contour were corrected manually on the basis of the spectrograms. Values of f0min, f0max, the depth of frequency modulation f0 (∆f0 = f0max - f0min) and average f0 of a call (f0mean) were taken automatically by using by using the Pitch info command in the Pitch edit window.
Two different methods of measuring f0max (one using Avisoft and another using Praat) applied to the same calls, resulted in very similar values. Coefficients of correlation, calculated separately for the oral and for the nasal calls of hinds and calves, ranged between 0.994 and 0.996 (0.988 < R2 < 0.992). Thus, for subsequent acoustic analyses we could select between these methods and we used the f0 values measured with Praat.
We did not measure formants, as they can only be measured either in low-frequency calls with closely spaced harmonics or in noisy calls (e.g. stag harsh roars), where the sound energy is dispersed over the call spectrum [26,37,55]. In most calls of hinds and in all calls of calves, formant frequencies could not be measured because they fell either between the harmonics and were thus invisible, or they were indistinguishable from the harmonics because they coincided with them in frequency.
Call samples and statistical analyses
Statistical analyses were conducted using STATISTICA v. 6.0 (StatSoft, Tulsa, OK, USA) and R v.3.0.1 . Means are given as mean ± SD, all tests were two-tailed, and differences were considered significant whenever p < 0.05. Distributions of most measured parameter values (excepting fpeak and peak-harm) did not depart from normality and distributions of all mean parameter values did not depart from normality (Kolmogorov-Smirnov test, p > 0.05). As parametric ANOVA and discriminant function analysis (DFA) are relatively robust to departures from normality , this was not an obstacle to the application of these tests.
We used a repeated measures ANOVA controlled for individuality, to compare the parameter values between oral and nasal calls. We used 1–15 (9.3 ± 5.2) calls per animal per call type, from 28 hinds (300 oral and 325 nasal calls) and from 31 calves (16 males, 15 females; 281 oral and 188 nasal calls). From hinds which provided calls in both years, we took calls recorded only in one of the years (2011 or 2012). For oral and nasal calls of these 28 hinds and 31 calves, we calculated average values of acoustic variables.
We used DFA to calculate the probability of the assignment of calls to the correct individual for three call samples in both hinds and calves (of nasal calls, oral calls and the pooled sample of oral and nasal calls). We took 5–10 (9.0 ± 1.5) calls per animal per call type from 22 hinds (209 oral and 212 nasal calls in total) and from 17 calves (10 males, 7 females; 149 oral and 134 nasal calls in total). DFA requires balanced sample sizes per group, and thus we excluded from the analyses all the animals with less than 5 oral/nasal calls. We took all the calls from all the animals with from 5 to 10 oral/nasal calls and randomly selected 10 calls per type from animals with more than 10 measured calls of each type. We included 11 of the 14 measured call variables in the DFA, excluding fpeak and peak-harm (for not meeting the criterion of normality), and f0min (because it was used for calculating another variable).
Then we investigated the stability of acoustic individuality of hind oral and nasal calls between years for hinds that provided calls in both years. We classified hind calls from 2012 year with DFA functions derived from 2011, considering the value of the correct cross-validation as a measure of the retention of individuality over time [39,43,44,67]. We used 7–15 (12.2 ± 3.4) oral calls per animal per year from 5 hinds, which provided a sufficient number of oral calls in both years (in total 65 oral calls in 2011 and 57 oral calls in 2012). Also, we used 10–15 (13.7 ± 1.8) nasal calls per animal from 9 hinds which provided sufficient numbers of nasal calls in both years (in total 125 nasal calls in 2011 and 121 nasal calls in 2012).
We used Wilks’ Lambda values to estimate how strongly acoustic variables of calls contribute to the discrimination of individuals. With a 2×2 Yates’ chi-squared test, we compared the values of correct assignment of nasal and oral calls to the correct individual. We used a repeated measures ANOVA to compare the percentages of correct assignment of oral and nasal calls to particular individuals. To validate our DFA results, we calculated the random values of correct assignment of calls to individual by applying randomization procedure with macros, created in R. The random values were averaged from DFAs performed on 1000 randomized permutations on the data sets as described by . For example, to calculate the random value of classifying oral calls to individual hinds, each permutation procedure included the random permutation of 209 calls among 22 randomization groups, respectively to 22 individual hinds which were examined, and followed by DFA standard procedure built-in in STATISTICA. All other permutation procedures were made similarly. Using a distribution obtained by the permutations, we noted whether the observed value exceeded 95%, 99% or 99.9% of the values within the distribution . If the observed value exceeded 95%, 99% or 99.9% of values within this distribution, we established that the observed value did differ significantly from the random one with a probability p < 0.05, p < 0.01 or p < 0.001 respectively [44,68,69].
Because body mass should theoretically be proportional to the cube of a linear dimension like body size, for hinds we used log3 body mass to calculate Pearson’s correlation between body mass and acoustic variables of oral and nasal calls. For calves we used ANCOVA (with sex as a fixed categorical factor and body mass as continuous factor) for estimating effects of sex, body mass and the conjoint effect of sex*body mass on acoustic variables of oral and nasal calls. We used one-way ANOVA to compare body mass values between male and female calves.
Differences in fundamental frequency between oral and nasal calls are age-specific in Iberian red deer. This finding is new for mammals and probably may be found in other Cervidae species. The higher fundamental frequency of oral than nasal calls of red deer young is consistent with available data on three Bovidae and one Elephantidae species [31,35,36,50], whereas the oral and nasal calls of hinds do not fit to this proposal common rule. The role of factors responsible for this discrepancy, e.g. differences in emotional arousal between mother and offspring red deer during brief separations, should be tested by further studies with controlled experimental conditions [e.g. 24,35].
Differences in power to discriminate individuals between oral and nasal calls are age-specific in Iberian red deer. Whereas in both hinds and calves, both oral and nasal call types allow accurate discrimination of individuals with DFA, the more pronounced individuality of oral calls compared to nasal calls was found only in hinds. Whereas the higher individuality of oral calls was found also in one Bovidae species , the same individuality of oral and nasal calls has no analogies in mammals.
Hind calls were unstable between seasons. These results indicate that the vocal individuality is not supported automatically even in mature animals which stopped their growth. This may point also on some costs of supporting vocal individuality over long time. Thus, vocal individuality is not permanently supported in hinds, but re-stores on the basis of new combinations of acoustic traits .
Discriminant function analysis
The fundamental frequency
The duration from call onset to the point of the maximum fundamental frequency
The fundamental frequency at the onset of a call
The fundamental frequency at the end of a call
The maximum fundamental frequency of a call
The minimum fundamental frequency of a call
The average fundamental frequency of a call
The depth of frequency modulation, calculated as the difference between f0max and f0 min
The frequency of maximum amplitude within a call
The lower quartile of a call, covering 25% of call energy
The medium quartile of a call covering 50% of call energy
The upper quartiles of a call covering 75% of call energy
The relative power of the f0 band compared to the peak frequency band
The order number of the harmonic with the maximum energy
Nowak R, Porter RH, Levy F, Orgeur P, Schaal B. Role of mother-young interactions in the survival of offspring in domestic mammals. Rev Reprod. 2000;5:153–63.
Lingle S, Wyman MT, Kotrba R, Teichroeb LJ, Romanov CA. What makes a cry a cry? a review of infant distress vocalizations. Curr Zool. 2012;58:698–726.
Charrier I, Mathevon N, Joventin P. Fur seal mother memorises growing pup’s voice steps: adaptation to long-term recognition or evolutionary by-product? Biol J Linn Soc. 2003;80:305–12.
Insley SJ, Phillips AV, Charrier I. A review of social recognition in pinnipeds. Aquatic Mammals. 2003;29:181–201.
Pitcher BJ, Harcourt RG, Charrier I. Rapid onset of maternal vocal recognition in a colonially breeding mammal, the Australian sea lion. PLoS One. 2010;5(8):e12195.
Bohn KM, Wilkinson GS, Moss CF. Discrimination of infant isolation calls by female greater spear-nosed bats, Phyllostomus hastatus. Anim Behav. 2007;73:423–32.
Knörnschild M, Feifel M, Kalko EKV. Mother-offspring recognition in the bat Carollia perspicillata. Anim Behav. 2013;86:941–8.
Jouventin P, Aubin T, Lengagne T. Finding a parent in a king penguin colony: the acoustic system of individual recognition. Anim Behav. 1999;57:1175–83.
Aubin T, Jouventin P, Hildebrand C. Penguins use the two-voice system to recognise each other. Proc Roy Soc Lond B. 2000;267:1081–7.
Jones IL, Falls JB, Gaston AJ. Vocal recognition between parents and young of Ancient Murrelets, Synthliboramphus antiquus (Aves: Alcidae). Anim Behav. 1987;35:1405–15.
Insley SJ, Paredes R, Jones IL. Sex differences in razorbill Alca torda parent-offspring vocal recognition. J Exp Biol. 2003;206:25–31.
Charrier I, Mathevon N, Jouventin P, Aubin T. Acoustic communication in a black-headed gull colony: how do chicks identify their parents? Ethology. 2001;107:961–74.
Mathevon N, Charrier I, Jouventin P. Potential of individual recognition in acoustic signals: a comparative study of two gulls with different nesting patterns. C R Biologies. 2003;326:329–37.
Klenova AV, Volodin IA, Volodina EV. The variation in reliability of individual vocal signature throughout ontogenesis in the red-crowned crane Grus japonensis. Acta Ethol. 2009;12:29–36.
Beecher MD. Successes and failures of parent-offspring recognition in animals. In: Hepper PG, editor. Kin recognition. Cambridge: Cambridge University Press; 1991. p. 94–124.
Nowak R. Senses involved in the discrimination of Merino ewes at close contact and from a distance by their newborn lambs. Anim Behav. 1991;42:357–66.
Terrazas A, Serafin N, Hernandez H, Nowak R, Poindron P. Early recognition of newborn goat kids by their mother: II. auditory recognition and evidence of an individual acoustic signature in the neonate. Dev Psychobiol. 2003;43:311–20.
Briefer EF, McElligott AG. Mutual mother-offspring recognition in an ungulate hider species (Capra hircus). Anim Cogn. 2011;14:585–98.
Keller M, Meurisse M, Poindron P, Nowak R, Ferreira G, Shayit M, et al. Maternal experience influences the establishment of visual/auditory, but not olfactory recognition of the newborn lamb by ewes at parturition. Devel Psychobiol. 2003;43:167–76.
Torriani MVG, Vannoni E, McElligott AG. Mother-young recognition in an ungulate hider species: a unidirectional process. Amer Nat. 2006;168:412–20.
Lingle S, Rendall D, Pellis SM. Altruism and recognition in the antipredator defence of deer: 1. species and individual variation in fawn distress calls. Anim Behav. 2007;73:897–905.
Lingle S, Rendall D, Wilson WF, Deyoung RW, Pellis SM. Altruism and recognition in the antipredator defence of deer: 2. why mule deer help nonoffspring fawns. Anim Behav. 2007;73:907–16.
Teichroeb LJ, Riede T, Kotrba R, Lingle S. Fundamental frequency is key to response of female deer to juvenile distress calls. Behav Process. 2013;92:15–23.
Sebe F, Nowak R, Poindron P, Aubin T. Establishment of vocal communication and discrimination between ewes and their lamb in the first two days after parturition. Dev Psychobiol. 2007;49:375–86.
Fitch WT, Hauser MD. Unpacking “honesty”: vertebrate vocal production and the evolution of acoustic signals. In: Simmons A, Fay RR, Popper AN, editors. Acoustic communication, springer handbook of auditory research. New York: Springer; 2002. p. 65–137.
Taylor AM, Reby D. The contribution of source–filter theory to mammal vocal communication research. J Zool. 2010;280:221–36.
Fant G. Acoustic theory of speech production. Mouton & Co: The Hague, Netherlands; 1960.
Vankova D, Malek J. Characteristics of the vocalizations of red deer Cervus elaphus hinds and calves. Bioacoustics. 1997;7:281–9.
Searby A, Jouventin P. Mother-lamb acoustic recognition in sheep: a frequency coding. Proc Roy Soc Lond B. 2003;270:1765–71.
Espmark Y. Individual characteristics in the calls of reindeer calves. Behaviour. 1975;54:50–9.
Volodin IA, Lapshina EN, Volodina EV, Frey R, Soldatova NV. Nasal and oral calls in juvenile goitred gazelles (Gazella subgutturosa) and their potential to encode sex and identity. Ethology. 2011;117:294–308.
Lapshina EN, Volodin IA, Volodina EV, Frey R, Efremova KO, Soldatova NV. The ontogeny of acoustic individuality in the nasal calls of captive goitred gazelles, Gazella subgutturosa. Behav Process. 2012;90:323–30.
Richardson LW, Jacobson HA, Muncy RJ, Perkins CJ. Acoustics of white-tailed deer (Odocoileus virginianus). J Mammal. 1983;64:245–52.
Efremova KO, Volodin IA, Volodina EV, Frey R, Lapshina EN, Soldatova NV. Developmental changes of nasal and oral calls in the goitred gazelle Gazella subgutturosa, a nonhuman mammal with a sexually dimorphic and descended larynx. Naturwissenschaften. 2011;98:919–31.
Sebe F, Duboscq J, Aubin T, Ligout S, Poindron P. Early vocal recognition of mother by lambs: contribution of low- and high-frequency vocalizations. Anim Behav. 2010;79:1055–66.
Volodin IA, Sibiryakova OV, Kokshunova LE, Frey R, Volodina EV. Nasal and oral calls in mother and young trunk-nosed saiga antelopes, Saiga tatarica. Bioacoustics. 2014;23:79–98.
Volodin I, Matrosova V, Volodina E, Garcia AJ, Gallego L, Márquez R, et al. Sex and age-class differences in calls of Iberian red deer during rut: reversed sex dimorphism of pitch and contrasting roars from farmed and wild stags. Acta Ethol. 2014. doi:10.1007/s10211-013-0179-8.
Lengagne T. Temporal stability in the individual features in the calls of eagle owls (Bubo bubo). Behaviour. 2001;138:1407–19.
Klenova AV, Volodin IA, Volodina EV. Examination of pair-duet stability to promote long-term monitoring of the endangered red-crowned crane (Grus japonensis). J Ethol. 2009;27:401–6.
Klenova AV, Zubakin VA, Zubakina EV. Inter- and intra-season stability of vocal individual signatures in a social seabird, the crested auklet. Acta Ethol. 2012;15:41–52.
Briefer EF, de la Torre MP, McElligott AG. Mother goats do not forget their kids’ calls. Proc Roy Soc Lond B. 2012;279:3749–55.
Puglisi L, Adamo C. Discrimination of individual voices in male great bitterns (Botaurus stellaris) in Italy. Auk. 2004;121:541–7.
Matrosova VA, Volodin IA, Volodina EV. The short-term and long-term individuality in speckled ground squirrel alarm calls. J Mammal. 2009;90:158–66.
Matrosova VA, Volodin IA, Volodina EV, Vasilieva NA, Kochetkova AA. Between-year stability of individual alarm calls in the yellow ground squirrel Spermophilus fulvus. J Mammal. 2010;91:620–7.
Briefer E, Vannoni E, McElligott AG. Quality prevails over identity in the sexually selected vocalisations of an ageing mammal. BMC Biol. 2010;8:35.
Reby D, Andre-Obrecht R, Galinier A, Farinas G, Cargnelutti B. Cepstral coefficients and hidden markov models reveal idiosyncratic voice characteristics in red deer (Cervus elaphus) stags. J Acoust Soc Am. 2006;120:4080–9.
Bocci A, Telford M, Laiolo P. Determinants of the acoustic behaviour of red deer during breeding in a wild alpine population, and implications for species survey. Ethol Ecol Evol. 2013;25:52–69.
Passilongo D, Reby D, Carranza J, Apollonio M. Roaring high and low: composition and possible functions of the iberian stag’s vocal repertoire. PLoS One. 2013;8:e63841.
Vankova D, Bartoš L, Malek J. The role of vocalisations in the communication between red deer hinds and calves. Ethology. 1997;103:795–808.
Stoeger AS, Heilmann G, Zeppelzauer M, Ganswindt A, Hensman S, Charlton BD. Visualizing sound emission of elephant vocalizations: evidence for two rumble production types. PLoS One. 2012;7:e48907.
Kidjo N, Cargnelutti B, Charlton BD, Wilson C, Reby D. Vocal behaviour in the endangered corsican deer: description and phylogenetic implications. Bioacoustics. 2008;18:159–81.
Gogoleva SS, Volodina EV, Volodin IA, Kharlamova AV, Trut LN. The gradual vocal responses to human-provoked discomfort in farmed silver foxes. Acta Ethol. 2010;13:75–85.
Gogoleva SS, Volodin IA, Volodina EV, Kharlamova AV, Trut LN. Sign and strength of emotional arousal: vocal correlates of positive and negative attitudes to humans in silver foxes (Vulpes vulpes). Behaviour. 2010;147:1713–36.
Briefer EF. Vocal expression of emotions in mammals: mechanisms of production and evidence. J Zool. 2012;288:1–20.
Frey R, Volodin I, Volodina E, Carranza J, Torres-Porras J. Vocal anatomy, tongue protrusion behaviour and the acoustics of rutting roars in free-ranging iberian red deer stags (Cervus elaphus hispanicus). J Anat. 2012;220:271–92.
Briefer E, McElligott AG. Indicators of age, body size and sex in goat kid calls revealed using the source–filter theory. Appl Anim Behav Sci. 2011;133:175–85.
Volodin I, Volodina E, Frey R, Carranza J, Torres-Porras J. Spectrographic analysis points to source-filter coupling in rutting roars of Iberian red deer. Acta Ethol. 2013;16:57–63.
Briefer EF, McElligott AG. Social effects on vocal ontogeny in an ungulate, the goat, Capra hircus. Anim Behav. 2012;83:991–1000.
Volodin IA, Volodina EV, Lapshina EN, Efremova KO, Soldatova NV. Vocal group signatures in the goitred gazelle Gazella subgutturosa. Anim Cogn. 2014;17:349–57.
Landete-Castillejos T, García A, Gómez JA, Molina A, Gallego L. Subspecies and body size allometry affect milk production and composition, and calf growth in red deer: comparison of Cervus elaphus hispanicus and C. e. scoticus. Physiol Biochem Zool. 2003;76:594–602.
Landete-Castillejos T, Garcia A, Gallego L. Calf growth in captive iberian red deer (Cervus elaphus hispanicus): effect of birth date and hind milk production and composition. J Anim Sci. 2001;79:1085–92.
Durbin LS. Individuality in the whistle call of the asiatic wild dog Cuon alpinus. Bioacoustics. 1998;9:197–206.
Reby D, McComb K. Anatomical constraints generate honesty: acoustic cues to age and weight in the roars of red deer stags. Anim Behav. 2003;65:519–30.
Boersma P, Weenink D. Praat: doing phonetics by computer. [http://www.praat.org/]
R foundation for statistical computing. [http://www.R-project.org]
Dillon WR, Goldstein M. Multivariate analysis: methods and applications. New York: Wiley; 1984.
Tripp TM, Otter KA. Vocal individuality as a potential long-term monitoring tool for western screech-owls, Megascops kennicottii. Can J Zool. 2006;84:744–53.
Solow AR. A randomization test for misclassification probability in discriminant analysis. Ecology. 1990;71:2379–82.
Matrosova VA, Volodin IA, Volodina EV, Vasilieva NA. Stability of acoustic individuality in the alarm calls of wild yellow ground squirrels Spermophilus fulvus and contrasting calls from trapped and free-ranging callers. Naturwissenschaften. 2010;97:707–15.
We thank S.S. Gogoleva and K.O. Efremova for help with data collection. We thank the staff of Albacete red deer farm for their help and support. We thank Maya Pollard for her courteous correction of English in the early draft of the manuscript. We thank James Hare and the anonymous reviewer for their most helpful and constructive comments. We adhered to the ‘Guidelines for the treatment of animals in behavioural research and teaching’ (Anim. Behav., 2006, 71, 245–253) and to the laws on animal welfare for scientific research of Spain and the Russian Federation, where the study was conducted. Data collection protocol # 2011–36 was approved by the Committee of Bio-ethics of Lomonosov Moscow State University. This study was supported by the Russian Scientific Foundation, grant No 14-14-00237 (for IV, EV, OS and VM) and by Ministerio de Economía y Competitividad, Spain, project AGL2012-38898 (for TL, AG and LG).
The authors declare that they have no competing interests.
OS, IV, VM, EV and AG collected the data. OS, IV and VM analysed the data. IV, EV, TL, AG and LG designed the study. IV, EV, OS, TL, AG and LG wrote the manuscript. All authors read and approved the final manuscript.
About this article
- Acoustic communication
- Mother-offspring recognition
- Iberian red deer
- Cervus elaphus hispanicus
- Separation calls