Open Access
Translator Disclaimer
21 July 2022 Neighbors matter: Vocal variation in Gentoo Penguins depends on the species composition of their colony
Helen Rößler, Maureen Lynch, Sara Torres Ortiz, Ole Næsbye Larsen, Michaël Beaulieu
Author Affiliations +
Abstract

When suitable nesting habitat is rare, birds may have to share it with heterospecific individuals with similar nesting requirements. The resulting species mosaic may in turn affect how breeding birds communicate vocally with each other. For instance, they may modify their vocalizations if the ambient noise produced by heterospecifics interferes with their own vocalizations. The colonies of Gentoo Penguins (Pygoscelis papua ellsworthi) are interesting to examine whether such a vocal variation occurs in colonial seabirds, as these birds frequently breed in the same colony as Adélie (P. adeliae) or Chinstrap (P. antarcticus) penguins. We investigated whether the presence of congeneric penguins breeding in the same colony as Gentoo Penguins affected the acoustic characteristics of their most common vocalization, the ecstatic display call. Based on vocalizations recorded in 23 breeding colonies along the Antarctic Peninsula during 3 field seasons, we found that the frequency of ecstatic display calls of Gentoo Penguins was consistently lower (average exhale phrase ∼460 Hz lower, average inhale phrase ∼370 Hz lower) as well as decreased in energy distribution when they bred in mixed colonies with Adélie Penguins than when they bred only in proximity to conspecifics. In contrast, the frequency of both inhale and exhale phrases was unaffected by the presence of Chinstrap Penguins, potentially due to the already greater frequency difference between these 2 species. The apparent vocal plasticity, in mixed colonies with Adélie Penguins, likely allows Gentoo Penguins to adjust the quality of their vocalizations depending on the acoustic space available in their colonies, and hence possibly enhance signal transmission between conspecifics. However, whether this acoustic adjustment is sufficient to increase the detectability of vocalizations and localization of individuals has yet to be determined.

LAY SUMMARY

  • We examined whether Gentoo Penguin calls differ between breeding colonies that included multiple species of penguins and those where only Gentoo Penguins bred.

  • We compared call parameters including frequency, structure, and duration between individual Gentoo Penguins breeding in colonies with and without the presence of Adélie or Chinstrap penguins along the Antarctic Peninsula during three austral summers.

  • We found that the calls of Gentoo Penguins showed a lower frequency and narrower energy distribution when they bred in proximity to Adélie Penguins but did not change when they bred in proximity to Chinstrap Penguins with higher main frequency vocalizations..

  • This suggests that penguins can adjust the quality of their vocalizations according to the social (and presumably acoustic) environment of their colonies.

Cuando el hábitat de anidación adecuado es raro, las aves pueden tener que compartirlo con individuos heteroespecíficos con requisitos de anidación similares. El mosaico de especies resultante puede, a su vez, afectar la forma en que las aves reproductoras se comunican vocalmente entre sí. Por ejemplo, pueden modificar sus vocalizaciones si el ruido ambiental producido por los individuos heteroespecíficos interfiere con sus propias vocalizaciones. Las colonias de Pygoscelis papua ellsworthi son interesantes para examinar si tal variación vocal ocurre en las aves marinas coloniales, ya que estas aves frecuentemente se reproducen en la misma colonia que P. adeliae o P. antarcticus. Investigamos si la presencia de pingüinos congéneres que se reproducen en la misma colonia que P. p. ellsworthi afecta las características acústicas de su vocalización más común, la llamada de despliegue extático. Con base en las vocalizaciones registradas en 23 colonias de reproducción a lo largo de la Península Antártica durante tres temporadas de campo, encontramos que la frecuencia de las llamadas de despliegue extático de P. p. ellsworthi fue consistentemente más baja (frase de exhalación promedio ∼ 460 Hz más baja, frase de inhalación promedio ∼370 Hz más baja) y disminuyó en la distribución de energía cuando se reproduce en colonias mixtas con P. adeliae que cuando se reproduce solo en la proximidad de sus congéneres. En contraste, la frecuencia de las frases de inhalación y exhalación no se vio afectada por la presencia de P. antarcticus, posiblemente debido a la diferencia de frecuencia ya mayor entre estas dos especies. La plasticidad vocal aparente, en colonias mixtas con P. adeliae, probablemente le permite a P. p. ellsworthi ajustar la calidad de sus vocalizaciones según el espacio acústico disponible en sus colonias y, por lo tanto, posiblemente mejorar la transmisión de señales entre congéneres. Sin embargo, aún no se ha determinado si este ajuste acústico es suficiente para aumentar la detectabilidad de las vocalizaciones y la localización de los individuos.

INTRODUCTION

Birds may nest in association with other bird species that simultaneously share similar nesting requirements because of heterospecific attraction or the paucity of available breeding sites (Mönkkönen and Forsman 2002, Morinay et al. 2020). The costs and benefits of such breeding associations between birds appear to be context dependent (Brussee et al. 2016, Swift et al. 2018, 2020) but for such associations to persist, the related benefits (e.g., better access to resources, predator detection, nest defense, information transfer) have to overcome the related costs (resource competition, disease transmission, nest parasitism; Farine et al. 2014).

These preponderant advantages likely underlie the existence of large bird colonies (i.e., aggregates of breeding birds nesting at a given site at a higher probability than predicted by a free distribution) that can be not only monospecific but also heterospecific (Evans et al. 2015). Specifically in seabirds, coloniality has been shown to be the most favored breeding structure, with 98% of seabird species being colonial species (Evans et al. 2015).

Seabirds, including penguins, exchange information with conspecifics mostly through vocalizations with specific acoustic characteristics. Within breeding grounds, these vocalizations are essential for partner recognition and bonding, feeding, and parent–chick recognition (Boucaud et al. 2016, Tyson et al. 2017). However, the active space of vocalizations (i.e., the area around the signal source over which the signal remains detectable and recognizable) may extend beyond the nest, which may help conspecifics nesting in the surroundings to synchronize their phenology and ultimately optimize their breeding performance. Such a situation can be observed in seabird colonies where nests are close to each other (Waas et al. 2005). In such dense colonies, acoustic space is limited and has to be shared between individuals occupying the same active space if the characteristics of their vocalizations overlap. Within breeding grounds shared by multiple species, competition for acoustic space may occur between birds of different species. To reduce such vocal competition, coexisting species may show distinct vocalizations, reduce distances between conspecifics (through the creation of meeting points) and a “courtesy rule” (i.e., abstention from calling when another individual is calling) may apply (Aubin and Jouventin 2002). If acoustic overlap persists, one of the species competing within the acoustic space may have to modify the quality of its vocalizations to optimize vocal transmission between conspecifics (acoustic avoidance hypothesis; Brumm and Slabbekoorn 2005). Just as birds may increase the frequency of their vocalizations in the presence of anthropogenic noise (Katti and Warren 2004, Patricelli and Blickley 2006, Nemeth et al. 2013), birds may also have to alter the characteristics of their vocalizations when exposed to heterospecific vocalizations. Thus, within breeding grounds, the available acoustic environment may be partitioned between bird species.

Suitable rocky habitats for breeding are relatively rare and fragmented along the Antarctic coastline, and as such, many Antarctic seabirds breed in close proximity and sometimes in multispecies colonies, as observed for Pygoscelis spp. penguins—Gentoo (P. papua), Adélie (P. adeliae) and Chinstrap (P. antarctica) penguins, though all 3 species are rarely found breeding in the same area. This association is made possible because of the use of different ecological feeding niches by each species, with Gentoo Penguins feeding mostly inshore on fish whereas Adélie and Chinstrap penguins forage more extensively on krill in more oceanic areas (Trivelpiece et al. 1987). Moreover, Gentoo Penguins show a higher foraging plasticity than Adélie and Chinstrap penguins (Miller et al. 2009, Casanovas et al. 2015) as well as a greater plasticity in breeding phenology (Lynch et al. 2012a), which may make them more resilient to environmental variation. This higher plasticity may explain why Gentoo Penguins represent the most abundant penguin species along the western Antarctic peninsula (Herman et al. 2020).

Similar to many seabirds, penguins use vocalizations mainly for the recognition and localization of their partner and offspring within colonies, and to maintain social relationships and group cohesion (Jouventin 1982, Favaro et al. 2014). They may show some vocal plasticity, as demonstrated in King Penguins (Aptenodytes patagonicus) that increase their number of calls and syllables per call to be better audible to their partner in response to noisy windy conditions (Lengagne et al. 1999). In Pygoscelis penguins, the ecstatic display call is the most frequent call used during the breeding period. It is hypothesized that it serves both to attract and identify partners and as an advertisement display of nest occupancy (Jouventin 1982, Favaro et al. 2014).

To our knowledge, how sympatric Pygoscelis penguins share their acoustic space in mixed colonies has not been examined so far. Despite very different display vocalizations between species and even though Pygoscelis penguins use frequency bands for individual recognition (Aubin and Jouventin 2002), vocalizations overlap considerably in frequency among Gentoo, Adélie, and Chinstrap penguins (mean main frequency: ∼1,285, 1,430, 2,000 Hz, respectively; maximum frequency: ∼ 6,535, 4,325, 4,000 Hz, respectively; Jouventin and Dobson 2017). Because of this acoustic overlap, one or more of these species might be limited in the use of the available acoustic space when breeding in contact with other species and have to alter the characteristics of its vocalizations to occupy a new vocal niche. It was recently found that colonies of Gentoo Penguins breeding along the Antarctic Peninsula differed in the frequency of their ecstatic display calls, with much variations coming from differences among colonies and individuals within colonies (Lynch and Lynch 2017). Because no clear correlation was found between the characteristics of these vocalizations and the geographical distribution of colonies, the reasons for this vocal variation between colonies remained unclear (Lynch and Lynch 2017). As these differences may lie in the acoustic properties of each colony (Nottebohm 1975) as well as in the interrelationship of phylogenetic affinities, ecological requirements, and acoustic similarities (Littlejohn 1977, Duellman and Pyles 1983), we hypothesize here that the presence of other sympatrically breeding Pygoscelis penguins may contribute to vocal differences between colonies of Gentoo Penguins. Such vocal adjustments may be highly adaptive, as Gentoo Penguins occupying mixed colonies are constantly exposed to the ambient noise caused by heterospecific vocalizations. Specifically, we expect Gentoo Penguins, on which we focused in our study, to be limited to using the frequency bands that are mostly inhabited by other species. Because the main frequency of Gentoo Penguins' vocalizations is closer to that of Adélie Penguins' vocalizations than to that of Chinstrap Penguins' vocalizations (Jouventin 1982, Jouventin and Dobson 2017), we predict that this frequency shift should be more dramatic when Gentoo Penguins breed in proximity to Adélie Penguins.

METHODS

Acoustic Recordings

Passive audio recordings were conducted during three austral summers (2014/2015, 2018/2019, 2019/2020) in a total of 23 Gentoo Penguin colonies (P. papua ellsworthi; Pertierra et al. 2020, Tyler et al. 2020) along the Western Antarctic Peninsula and the South Shetland Islands. The colonies were classified as exclusive species colonies with only Gentoo Penguins (G), or multispecies colonies of Gentoo and Adélie penguins (GA) or Gentoo and Chinstrap penguins (GC, Figure 1). Recordings were conducted at the end of the incubation stage and at the beginning of the chick-rearing period. The spatial distribution of the multispecies colonies that we sampled in our study reflects the natural distribution of the different penguin species along the Antarctic Peninsula (Black 2016).

During the breeding season of 2014/2015, passive soundscape recorders (Song Meter SM2+, 24,000 Hz sampling rate, stereo recordings) were placed 3–5 m from one or more small subgroups of nesting penguins (for more details see Lynch and Lynch 2017 from which part of the data was used). In the breeding seasons of 2018/2019 and 2019/2020, recordings were carried out by a person positioned at the edge of a subgroup (i.e., the individuals whose vocalizations we could record above the noise during a recording session, ∼10 nests; Mustafa et al. 2017) of nesting Gentoo Penguins and using a unidirectional microphone (Sennheiser ME66 and ME67, Wedemark, Germany, a frequency response rate of 40–20,000 Hz, the sensitivity of 50 mV/Pa ± 2.5 dB at 1kHz) connected to a recorder (Olympus LS 100, Tokyo, Japan; internal amplification level 2–5 (5 = +4.4 dB), 48–96 kHz sampling rate of 16–24-bit resolution, respectively) and directed towards a given penguin showing ecstatic displays in the absence of its mate. Each location was only recorded at one resolution rate in the same year. To additionally make sure that bit resolution did not have an effect on the recordings or measurements, we ran a playback experiment, recording exactly the same call, one time with 16-bit resolution, the other time with 24-bit resolution. We could not find any difference. The microphone was mounted on a boom pole to better approach and target specific penguins at 3–5 m distance and was protected with a windjammer. The sensitivity of the microphones was confirmed before fieldwork by the manufacturer and by using a Brüel & Kjaer 4230 sound calibrator (SPL = 94 dB–1kHz) both before and after fieldwork. In addition to the recording of Gentoo Penguins, we also recorded some ecstatic display calls of Adélie (n = 20) and Chinstrap (n = 37) penguins, to estimate how the calls emitted by Gentoo Penguins acoustically compared with those of other species.

FIGURE 1.

Map showing the sampling locations along the Antarctic peninsula and the South Shetland Islands. Navy colored dots indicate exclusive Gentoo colonies (G), red dots indicate Gentoo colonies mixed with Chinstrap Penguins (GC) and yellow dots indicate Gentoo colonies mixed with Adélie Penguins (GA) (map created in R with marmap and getNOAA.bathy; Pante and Simon-Bouhet 2013).

img-z4-4_01.jpg

Analysis of Acoustic Sequences

Recordings were screened using Adobe Audition* 3.0 (Adobe Systems Inc., San Jose, CA, USA) and ecstatic display calls were manually selected from the recordings. An ecstatic display call was defined as any display call performed by a single individual with a repeated series of long, low-frequency exhale phrases and short, higher-frequency inhale phrases (Jouventin 1982, Lynch and Lynch 2017; Figure 2). Signal to noise ratio (SNR) was obtained in Raven sound analysis software (Raven Pro 1.5, Center for Conservation Bioacoustics, Cornell Lab of Ornithology, Ithaca, NY) by taking the difference between the measured energy of the selected area with the signal and the selected area without any signal at the same recording period. For further analysis, we only selected high-quality calls (SNR > 11 dB, without any background interference, e.g., of strong wind noise, calling chicks, and other individuals calling at the same time) resulting in a total of 706 individual ecstatic display calls. Among these calls, 419 came from 13 exclusive Gentoo Penguin colonies, 104 came from three colonies mixed with Adélie Penguins, and 183 from 7 colonies mixed with Chinstrap Penguins (Table 1).

To exclude low-frequency background wind noise, each selected call was filtered, with Raven sound analysis software (at 100 or 150 Hz) for the 2014/2015 recordings and for better control of filter settings with a custom-written script (MATLAB R2019a, 7th order Butterworth high pass filter at 100 Hz) for the 2018/2019 and 2019/2020 recordings. All adjusted ecstatic display calls were analyzed with Raven sound analysis software (window size = 625 samples, overlap = 65%, discrete Fourier transform (DFT) size = 2,048 samples). We measured energy [dB re 20 µPa], 90% duration [s], center frequency [Hz], 5% frequency [Hz], 95% frequency [Hz], peak frequency contour (PFC), PFC slope and peak frequency inflection points (for definitions, see Table 2; Charif et al. 2010). After filtering, 95% frequency of background noise was consistently below the 5% frequency of any ecstatic call measured, minimizing the possibility that background noise or variation in filtering interfered with our analysis. Due to the highly repetitive and variable number of phrases, we included in the analysis only measurements for the whole call and the first exhale and first inhale phrase. The first phrases are the clearest and are sufficient for individual identification (Jouventin and Aubin 2002). Received levels were only determined to rule out major transmission loss from off-axis or far distant calls but were not further analyzed because they could not be standardized between recordings due to variable distance between penguins and microphone, and variable orientation of the penguin relative to the microphone. Energy distribution within the ecstatic display calls was measured as aggregate entropy from the spectrogram (Raven Pro 1.5, Center for Conservation Bioacoustics, Cornell Lab of Ornithology, Ithaca, NY) of Gentoo Penguins in exclusive colonies (n = 364), Gentoo Penguins in mixed colonies (n = 275), as well as Adélie Penguins (n = 22) and Chinstrap Penguins (n = 33) from opportunistically taken mixed colonies, as well as Chinstrap Penguins (n = 9) in an exclusive colony. Aggregate entropy measures disorder in a sound by analyzing the energy distribution within a selection. Higher entropy values correspond to greater disorder in the sound whereas a pure tone with energy in only one frequency bin would have zero entropy. The size of a frequency bin is determined by the spectrogram parameters (Charif et al. 2010). As the DFT size was kept the same throughout the measurements independent of sampling rate, we compared 6 randomly selected calls of 48 kHz and 96 kHz analyzed with a DFT size of 2,048 samples to the adjusted DFT sizes of 4,096 samples or 8,192 samples, respectively and statistically tested for differences. We could not find any significant differences for our here analyzed parameters (paired t-test, df = 11, P > 0.05), except for peak frequency inflection points but this was also accounted for in the statistical models. The additional aggregate entropy measurements data were adjusted for differences in sampling rate by a statistically tested consistent factor (determined by the mean difference between randomly chosen original sampled calls and the down sampled version, respectively, 48 kHz + 0.9 and 96 kHz + 1.8; paired t-test, t = 0.41, df = 9, p = 0.7).

FIGURE 2.

Shown are spectrograms of a selected typical complete ecstatic display call by a Gentoo Penguin with the repeating series of the starting longer exhale (broadband pulses) and shorter following inhale phrases (harmonic stacks) and a zoomed in display of only the first exhale and inhale phrases in an exclusive Gentoo colony (G) and in a mixed Gentoo and Adélie penguin colony (GA). Darker hues indicate more power at that frequency. Also visible in the zoomed in spectrograms are the 5% frequency (bottom green line), center frequency (middle yellow line), 95% frequency (top red line), and the PFC (light blue squiggly line). Blue and red boxes show the selection windows. Spectrogram details: Window type: Hann, window size = 625 samples, overlap = 65%, DFT size = 8,192 and 4,096 samples with a 96 and 48 kHz sampling rate, respectively.

img-z5-4_01.jpg

TABLE 1.

Summary of the recording sessions of 3 field seasons along the Western Antarctic Peninsula and at the South Shetland Islands. Location, number of Gentoo Penguin ecstatic display calls recorded, and the species breeding are noted for each colony; G = Gentoo Penguins, A = Adélie Penguins, C = Chinstrap Penguins.

img-z5-7_01.gif

TABLE 2.

Descriptions of spectrogram measurements used in our analysis.

img-z6-5_01.gif

Statistical Analysis

In 2018–2019 and 2019–2020, we could assign a given individual to a given call, as we actively selected given individuals during recording sessions. In contrast, it was not possible to assign a given individual to a given call for vocalizations recorded in 2014–2015, as individuals were passively recorded (Lynch and Lynch 2017). However, in this previous study, a significant variation within colonies in a manner typical of inter-individual variation was found. Combined with a known low rate of recorded repeated calls of identified individuals within one colony, we therefore considered each of the calls analyzed from 2014/2015 in our study was also produced by a given individual.

Due to the highly correlated number of inhale and exhale phrases per call (rs = 0.84, P < 0.001), both variables were combined into one (total number of phrases per call) for further analyses. To reduce the number of vocal parameters to analyze, we then conducted principal component analyses (PCA). PCA loadings were limited to components with an eigenvalue ≥ 1 and varimax rotated for easier visual differentiation. The PCA model that was used for further analysis contained components explaining the highest variance (83%, Table 3) and variables showing no cross-loadings on different axes. We then used PC1 to PC5 scores in generalized linear mixed-effect models (GLMMs; Gaussian distribution) with colony composition as an explanatory variable (fixed factor) and, year and location nested within the year as random variables (to also account for differences in recording methods between years). We also repeated these analyses within each year to examine the repeatability of our results despite the use of different approaches in the field. The normality of residuals was checked through the visual inspection of QQ-plots and by conducting Shapiro–Wilk tests. Data were log-transformed (and mirrored in case of left skewed data) when their distribution was not normal. When this transformation was not sufficient for residuals to be normally distributed, a Gamma error distribution (log link function) was used in the GLMM (lme4 package; Bates et al. 2015). Using an analysis of variance (ANOVA), we tested the significance of a “full model” comprising both the fixed and random factors against a “null model” comprising the random factor only. We additionally re-conducted statistical tests for relevant acoustic parameters (95% frequency and aggregate entropy; see Results) independently to examine more precisely how colony composition affected individual acoustic parameters (when those reflected PC scores or parameters that were affected by colony composition). A nonparametric comparison (Kruskal–Wallis test) was used to test for significant differences between means and a Wilcoxon test for variation across years. To ensure that geographic distance alone did not explain the acoustic variation described here, we performed a Spearman's rank correlation test to assess whether a correlation existed between divergence in Gentoo acoustic features (using the difference in mean PC1 scores) and latitude or longitude across all colonies.

All statistical tests were conducted in R for Windows (The R Foundation 2022) with R Studio* 1.1.456. Data are provided in the result section as median ± SD.

TABLE 3.

Listed are loadings of variables of the best-fitting principal component analysis (PCA) where the first 5 components (PC1–PC5) explain 83% of variance.

img-z7-5_01.gif

RESULTS

Colony composition significantly affected PC1 scores (ANOVAPC1, χ2 = 11.034, P = 0.004; Figure 3), and these results were consistent across years when calls were recorded in the 3 colony compositions, even though the scores were overall lower in 2014/2015 than in 2019/2020 (2014/2015: loge (WWilcoxon) = 9.26, P = < 0.001, r̂ = 0.30, 95% CI: 0.20–0.40, nobs = 280, 2019/2020: loge (WWilcoxon) = 7.94, P = < 0.001, r̂ = 0.42, 95% CI: 0.27–0.58, nobs = 145; Figure 4). Specifically, lower PC1 scores were found when Gentoo Penguins bred in mixed colonies with Adélie Penguins (GLMM, estimate = –0.106, P < 0.001; Figure 3). Such a pattern was not found when Gentoo penguins bred in mixed colonies with Chinstrap Penguins (GLMM, estimate = 0.009, p = 0.75; Figure 3). Colony composition did not affect any other components of the ecstatic display calls of Gentoo Penguins (ANOVAPC2, χ2 = 0.52, P = 0.77, ANOVAPC3, χ2 = 0.55, P = 0.76, ANOVAPC4, χ2 = 1.50, P = 0.47, ANOVAPC5, χ2 = 3.23, P = 0.20).

As PC1 scores mostly reflected the 95% frequency of ecstatic display calls, we re-conducted our analyses for this given parameter. The 95% frequency of the exhale and the inhale phrases was on average 458 Hz and 367 Hz lower, respectively, when Gentoo Penguins bred in mixed colonies with Adélie Penguins than when they bred by themselves or in mixed colonies with Chinstrap Penguins (exhale: χ2 = 29.32, P = < 0.001, ε2 = 0.04, 95% CI: 0.02–0.09, nobs = 706; inhale: χ2 = 11.83, P = 0.003, ε2 = 0.02, 95% CI: 0.003–0.06, nobs = 706; Figure 5). Specifically, the complete ecstatic display calls of Gentoo Penguins showed a median 95% frequency of 2,730 ± 691 Hz in the exclusive Gentoo Penguin colonies (n = 419) and of 2625 ± 584 Hz in the colonies with Chinstrap Penguins (n = 183), whereas it showed a median 95% frequency of 2,344 ± 697 Hz in the colonies shared with Adélie Penguins (n = 104). In comparison, random selected ecstatic display calls of Adélie (n = 22) and Chinstrap (n = 33) penguins showed a median 95% frequency of 2,625 ± 482 and 2,953 ± 666 Hz, respectively (Figures 2 and 7).

FIGURE 3.

PC1 scores of Gentoo calls compared by breeding colony composition (G = exclusive Gentoo Penguins breeding colony, GA = Gentoo and Adélie penguins breeding colony, GC = Gentoo and Chinstrap penguins breeding colony).

img-z7-9_01.jpg

Energy distribution in the complete ecstatic display calls of Gentoo Penguins narrowed down to fewer frequencies used when Gentoo Penguins bred in mixed colonies with Adélie Penguins than when they bred by themselves or in mixed colonies with Chinstrap Penguins (χ2 = 52.86, df = 2, P < 0.0001, ε2 = 0.083, 95% CI: 0.042–0.132, nobs = 639; Figures 6A and 7). Energy distribution in the complete ecstatic display calls of Gentoo Penguins in mixed colonies with either Adélie or Chinstrap Penguins is lower compared to calls of Adélie Penguins (χ2 = 21.98, P < 0.0001, df = 1, nobs = 127) and Chinstrap Penguins (χ2 = 4.22, P = 0.04, df = 1, nobs = 203; Figures 6B–C and 7). Chinstrap Penguin ecstatic display calls in an exclusive Chinstrap penguin colony showed no change in energy distribution compared to when they bred in mixed colonies with Gentoo Penguins (χ2 = 0.12, P = 0.72, df = 1, nobs = 42; Figure 6D).

FIGURE 4.

PC1 scores plotted as grouped by year and by colony composition (G = exclusive Gentoo Penguins breeding colony [blue dots], GA = Gentoo and Adélie penguins breeding colony [yellow dots], GC = Gentoo and Chinstrap penguins breeding colony [red dots]).

img-z8-3_01.jpg

There was no significant correlation between geographic locations of all colonies by latitude or longitude, and acoustic parameters (mean PC1 scores) of the Gentoo Penguins ecstatic display calls (latitude: P = 0.67, rs [1,836] = 0.09; longitude: P = 0.54, rs [2,294] = –0.13).

DISCUSSION

We found that the presence of other penguin species breeding in colonies of Gentoo Penguins was associated with a qualitative change in their vocalizations (frequency, duration, contour). Interestingly, this difference between colonies, mostly reflecting a change in the frequency and energy distribution of vocalizations, was only observed when Gentoo Penguins bred in proximity to Adélie Penguins but not when they bred in proximity to Chinstrap Penguins, presumably because of a stronger acoustical overlap between the vocalizations of Gentoo and Adélie penguins. The geographic locations of colonies along the Antarctic Peninsula did not explain the vocal variation described here. Based on these results, our original hypothesis stipulating that birds sympatrically breeding with other birds need to adjust the use of their acoustic space by modulating their vocalizations therefore cannot be rejected.

FIGURE 5.

Box-violin plots showing the median and quartiles (A) 95% frequency of the first exhale and (B) first inhale phrase of ecstatic display calls of Gentoo Penguins plotted as grouped by colony composition (G = exclusive Gentoo Penguins breeding colony, GA = Gentoo and Adélie penguins breeding colony, GC = Gentoo and Chinstrap penguins breeding colony).

img-z8-11_01.jpg

In mixed colonies with Adélie and Gentoo penguins, Adélie Penguins occupy similar acoustic space as Gentoo Penguins use in exclusive colonies. In mixed colonies, the average frequency difference of ∼300 Hz that we detected between both penguin species might be biologically meaningful, as it corresponds to the critical ratio bandwidth difference for each species to distinguish conspecific vocalizations (Langemann et al. 1995). Average Chinstrap vocalizations already exceed that ∼300 Hz frequency difference to Gentoo Penguin ecstatic display calls and would render relative frequency change unnecessary. Frequency values of harmonics and timbre have been shown to be important factors for parent chick recognition in penguins, and the associated critical bandwidth can be as low as 25–50 Hz in the main frequency of Gentoo and Adélie penguins (Jouventin and Aubin 2002). Due to characteristics of cochlear processing, the detection of lower frequencies has a smaller critical bandwidth than higher frequencies, which might explain the lowering in Gentoo Penguin call frequency as well as the narrowing of energy distribution we observed in mixed colonies (instead of increased frequencies; Langemann et al. 1995, Dooling et al. 2000, Shera et al. 2002, Henry and Lucas 2010). Especially at fundamental frequencies of 200–400 Hz, masking is less pronounced in songbirds, which was suggested to reflect general avian anatomical and physiological mechanisms in the early peripheral auditory system (Dooling et al. 2001). Alternatively, it may be hypothesized that Gentoo Penguins avoid using high frequencies when breeding with Adélie Penguins, as such vocalizations may be associated with a higher cost/benefit ratio under such conditions. Indeed, vocalizations with higher frequencies may be more costly to produce (as a result of a stronger contraction of the syringeal muscles; Suthers et al. 1999) whereas their detectability may be limited when overlapping with other vocalizations with similarly high frequencies. The presumed higher cost of vocalizations with higher frequencies may explain why Gentoo Penguins do not raise the frequency of their vocalizations when breeding with Adélie Penguins, but rather use lower and narrower frequency bands. This vocal plasticity would fit quite well with the general behavioral plasticity of Gentoo Penguins (Miller et al. 2009, Lynch et al. 2012a, Casanovas et al. 2015).

FIGURE 6.

Box-violin plots showing the median and quartiles aggregate entropy (A) of ecstatic display calls of Gentoo Penguins grouped by colony composition (G = exclusive Gentoo Penguins breeding colony, GA = Gentoo and Adélie penguins breeding colony, GC = Gentoo and Chinstrap penguins breeding colony); (B) of ecstatic display calls grouped by penguin species in mixed GA colonies (A = Adélie Penguins, G = Gentoo Penguins); (C) of ecstatic display calls grouped by penguin species in mixed GC colonies (C = Chinstrap Penguins and G = Gentoo Penguins); (D) of ecstatic display calls of Chinstrap Penguins plotted as grouped by colony composition (C = exclusive Chinstrap Penguins breeding colony). Statistical categories shown are: ns (not significant) is p > 0.05, *p ≤ 0.05, **p ≤ 0.01, ***P ≤ 0.001, ****P ≤ 0.0001.

img-z9-4_01.jpg

Even though penguins are not known to experience vocal production learning (Tyack 2020), our results suggest that penguins might be able to learn vocal production to some limited extent, by adjusting the acoustic properties of their inherited vocal motor patterns based on the social structure and the acoustic environment of the colony where they were raised (Searcy et al. 2021). Discrimination of neighbors, strangers, and mates in Adélie Penguins supports the learning of calls in high nest fidelity birds (Speirs and Davis 1991, Beaulieu et al. 2009). It would therefore be interesting to examine the ontogeny of calls across the development of chicks in the different types of colonies we considered.

It might also be hypothesized that Gentoo Penguins able to breed in proximity to Adélie Penguins differ in some traits affecting vocalizations from Gentoo Penguins breeding in exclusive colonies. For instance, body size can affect frequency parameters in birds, with larger individuals typically producing lower frequencies (Ratcliffe and Grant 1985, Podos 1997, 2001). It is therefore possible that only larger Gentoo Penguins, using lower frequencies, are able to breed in colonies shared with Adélie Penguins, thereby explaining the pattern we observed. Thus, the selection for stronger and larger animals nesting in mixed colonies would lead to a lower frequency purely through larger vocal organs and vocal tracts which are better capable to produce and resonate low frequencies, and coupling them to the outside medium (Friis et al. 2021). Such a selection for larger individuals might be possible precisely because the deeper vocalizations of larger individuals provide them with an adaptive advantage by making their vocalizations more distinguishable, thereby facilitating intraspecific communication (Aubin and Jouventin 1998, Jouventin and Aubin 2002). In that case, smaller Gentoo Penguins with higher-frequency vocalizations might be counter-selected from colonies where Adélie Penguins breed because such vocalizations do not allow a proper intraspecific communication under such acoustic conditions. This scenario suggests that Adélie Penguins may limit smaller Gentoo penguins to breed where they could otherwise. Such a potential limitation may explain why colonies of Gentoo Penguins are rarer in association with Adélie Penguins than exclusive Gentoo colonies or colonies of Gentoo Penguins in association with Chinstrap Penguins.

FIGURE 7.

Power spectra of examples of typical Gentoo Penguin calls in all different kinds of examined colony types (G exclusive, GA mixed, GC mixed) as well as of Adélie Penguin calls and Chinstrap Penguin calls in the mixed colony type. Green vertical line = 5% frequency measurement, yellow vertical line = center frequency measurement, red vertical line = 95% frequency measurement. Spectrum details: Window type: Hann, window size = 625 samples, overlap = 65%, DFT size = 2,048, 4,096 or 8,192 samples at a 24, 48 or 96 kHz sampling rate, respectively.

img-z10-4_01.jpg

Our correlative results did not enable us to disentangle the possible mechanisms underlying the vocal differences we found between colony compositions (selection for given individuals depending on their size vs. vocal plasticity) and may be affected by some confounding factors (e.g., variable distance to heterospecific nests, colony topography, and structure). Moreover, the number of mixed colonies included in our study was limited (because their number is naturally limited). Finally, because of the few exclusive colonies of Adélie and Chinstrap penguins present along the Antarctic Peninsula, our correlative study focused solely on more frequent colonies of Gentoo Penguins. It is for instance possible that Gentoo Penguins did not modify their vocalizations in the presence of Chinstrap Penguins because Chinstrap Penguins show themselves some vocal plasticity, leading to less acoustic overlap in those mixed colonies. All these limitations call for the use of experimental approaches, such as approaches based on the use of heterospecific playbacks and their effects on the vocalizations of focal individuals. The use of such approaches would allow us to examine whether penguins can directly alter the quality of their vocalizations in response to their exposure to other vocalizations. Moreover, this experimental approach could also be used to examine whether the vocal plasticity of penguins differs between the three Pygoscelis species, or even further to assess the applicability of frequency changes as limited vocal learning mechanisms in other colonial bird nesting associations.

The coexistence of the 3 Pygoscelis species in colonies along the coast of the Antarctic Peninsula is made possible through the use of different foraging niches (Croxall and Prince 1980, Trivelpiece et al. 1987, Wilson 2010) and breeding phenology niches (Trivelpiece et al. 1987). However, climate change is strongly affecting Pygoscelis penguin populations and distributions along the Antarctic Peninsula. As Gentoo Penguins increase in abundance and expand their range southward (Hinke et al. 2007, Lynch et al. 2012b, Herman et al. 2020), interspecies competition for physical space, and subsequently acoustic space may also increase or shift. How climate change and competition for acoustic space interact in bird nesting associations therefore remains an important aspect to explore.

ACKNOWLEDGEMENTS

We deeply thank L. Smith and F. Guerrero of Quixote Expeditions to invite us (M. Beaulieu and H. Rößler) as guest scientists on board Ocean Tramp, which crews and passengers actively supported us in the field. For additional administrative support, we thank Michael Dähne and Magnus Wahlberg. We also thank Oceanwide Expeditions and Troels Jakobsen for taking us on board (S. T. Ortiz and O. N. Larsen). For additional active support in the field, we thank Morgan Martin and Lee Miller. Finally, we are grateful to Patrick Schwarzbach for his programming support.

Funding statement: The expenses of this work were shared by the Umweltbundesamt (UBA; German Environment Agency) and the Deutsches Meeresmuseum (DMM, German Oceanographic Museum). No funders had input into the content of the manuscript or required approval of the manuscript before submission or publication.

Ethics statement: This research did not require the handling of any animals. Data collection was conducted in 2014–2015 under Stony Brook University IACUC no. 237420, Antarctic Conservation Act no. ACA2014-024. In 2018–2019 and 2019–2020, data collection in the field was permitted based on the Act on the Implementation of the Environmental Protection Protocol to the Antarctic Treaty (AUG) Scientific research according to §6 AUG no. II 2.8-94033/141 (2018-2019) and II 2.8-94033/148 (2019-2020) granted by the German Environment Agency. In addition, in 2019–2020 recordings were granted by the authorities of the Netherlands Permit no. RWS-2019/10270.

Conflict of interest statement: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Author contributions: HR and MB conceived the idea for the study; ML, MB, HR, STO, and ONL collected data; HR compiled and analyzed the data, made tables and figures HR wrote a paper under the guidance of MB; ML, STO, ONL further contributed to the study through additional writing and editing of the manuscript.

Data availability: Analyses reported in this article can be reproduced using the data provided by Rößler et al. (2022).

LITERATURES CITED

1.

Aubin, T., and P. Jouventin (1998). Cocktail–party effect in king penguin colonies. Proceedings of the Royal Society of London, Series B: Biological Sciences 265:1665–1673. Google Scholar

2.

Aubin, T., and P. Jouventin (2002). How to vocally identify kin in a crowd: The penguin model. Advances in the Study of Behavior 31:243–277. Google Scholar

3.

Bates, D., M. Mächler, B. Bolker, and S. Walker (2015). Fitting linear mixed-effects models using lme4. Journal of Statistical Software 67:1–48. Google Scholar

4.

Beaulieu, M., A. -M. Thierry, Y. Le Maho, Y. Ropert-Coudert, and A. Ancel (2009). Alloparental feeding in Adélie Penguins: Why is it uncommon? Journal of Ornithology 150:637–643. Google Scholar

5.

Black, C. E. (2016). A comprehensive review of the phenology of Pygoscelis penguins. Polar Biology 39:405–432. Google Scholar

6.

Boucaud, I. C. A., M. M. Mariette, A. S. Villain, and C. Vignal (2016). Vocal negotiation over parental care? Acoustic communication at the nest predicts partners' incubation share. Biological Journal of the Linnean Society 117:322–336. Google Scholar

7.

Brumm, H., and H. Slabbekoorn (2005). Acoustic communication in noise. Advances in the Study of Behavior 35:151–209. Google Scholar

8.

Brussee, B. E., P. S. Coates, R. L. Hothem, K. B. Howe, M. L. Casazza, and J. M. Eadie (2016). Nest survival is influenced by parental behaviour and heterospecifics in a mixed-species colony. Ibis 158:315–326. Google Scholar

9.

Casanovas, P., R. Naveen, S. Forrest, J. Poncet, and H. J. Lynch (2015). A comprehensive coastal seabird survey maps out the front lines of ecological change on the western Antarctic Peninsula. Polar Biology 38:927–940. Google Scholar

10.

Charif, R., A. Waack, and L. Strickman (2010). Raven Pro 1.4 User's Manual.  https://ravensoundsoftware.com/wp-content/uploads/2017/11/Raven14UsersManual.pdfGoogle Scholar

11.

Croxall, J. P., and P. A. Prince (1980). Food, feeding ecology and ecological segregation of seabirds at South Georgia. Biological Journal of the Linnean Society 14:103–131. Google Scholar

12.

Dooling, R. J., M. L. Dent, M. R. Leek, and O. Gleich (2001). Masking by harmonic complexes in birds: Behavioral thresholds and cochlear responses. Hearing Research 152:159–172. Google Scholar

13.

Dooling, R. J., B. Lohr, and M. L. Dent (2000). Hearing in birds and reptiles. InComparative Hearing: Birds and Reptiles ( R. J. Dooling, R. R. Fay and A. N. Popper, Editors). Springer New York, New York, NY, USA. pp. 308–359. Google Scholar

14.

Duellman, W. E., and R. A. Pyles (1983). Acoustic resource partitioning in Anuran communities. Copeia 1983:639. Google Scholar

15.

Evans, J., S. Votier, and S. Dall (2015). Information use in colonial living. Biological Reviews of the Cambridge Philosophical Society 91:658–672. Google Scholar

16.

Farine, D. R., C. P. Downing, and P. A. Downing (2014). Mixed-species associations can arise without heterospecific attraction. Behavioral Ecology 25:574–581. Google Scholar

17.

Favaro, L., L. Ozella, and D. Pessani (2014). The vocal repertoire of the African Penguin (Spheniscus demersus): Structure and function of calls. PLoS One 9:e103460. Google Scholar

18.

Friis, J. I., T. Dabelsteen, and G. C. Cardoso (2021). Contingency and determinism in the evolution of bird song sound frequency. Scientific Reports 11:11600. Google Scholar

19.

Henry, K. S., and J. R. Lucas (2010). Habitat-related differences in the frequency selectivity of auditory filters in songbirds: Songbird auditory filters. Functional Ecology 24:614–624. Google Scholar

20.

Herman, R., A. Borowicz, M. Lynch, P. Trathan, T. Hart, and H. Lynch (2020). Update on the global abundance and distribution of breeding Gentoo Penguins (Pygoscelis papua). Polar Biology 43:1947–1956. Google Scholar

21.

Hinke, J. T., K. Salwicka, S. G. Trivelpiece, G. M. Watters, and W. Z. Trivelpiece (2007). Divergent responses of Pygoscelis penguins reveal a common environmental driver. Oecologia 153:845. Google Scholar

22.

Jouventin, P. (1982). Visual and vocal signals in penguins, their evolution and adaptive characters. Fortschritte der Verhaltensforschung 24:148–148. Google Scholar

23.

Jouventin, P., and T. Aubin (2002). Acoustic systems are adapted to breeding ecologies: Individual recognition in nesting penguins. Animal Behaviour 64:747–757. Google Scholar

24.

Jouventin, P., and F. S. Dobson (2017). Why Penguins Communicate: The Evolution of Visual and Vocal Signals. Academic Press, New York, NY, USA. Google Scholar

25.

Katti, M., and P. S. Warren (2004). Tits, noise and urban bioacoustics. Trends in Ecology & Evolution 19:19–20. Google Scholar

26.

Langemann, U., G. M. Klump, and R. J. Dooling (1995). Critical bands and critical-ratio bandwidth in the European starling. Hearing Research 84:167–176. Google Scholar

27.

Lengagne, T., T. Aubin, J. Lauga, and P. Jouventin (1999). How do King Penguins (Aptenodytes patagonicus) apply the mathematical theory of information to communicate in windy conditions? Proceedings of the Royal Society B 266:1623–1628. Google Scholar

28.

Littlejohn, M. J. (1977). Long-range acoustic communication in Anurans: An integrated and evolutionary approach. InThe Reproductive Biology of Amphibians ( D. H. Taylor and S. I. Guttman, Editors). Springer US, Boston, MA, USA. pp. 263–294. Google Scholar

29.

Lynch, H. J., W. F. Fagan, R. Naveen, S. G. Trivelpiece, and W. Z. Trivelpiece (2012a). Differential advancement of breeding phenology in response to climate may alter staggered breeding among sympatric pygoscelid penguins. Marine Ecology Progress Series 454:135–145. Google Scholar

30.

Lynch, H. J., R. Naveen, P. N. Trathan, and W. F. Fagan (2012b). Spatially integrated assessment reveals widespread changes in penguin populations on the Antarctic Peninsula. Ecology 93:1367–1377. Google Scholar

31.

Lynch, M. A., and H. J. Lynch (2017). Variation in the ecstatic display call of the Gentoo Penguin (Pygoscelis papua) across regional geographic scales. The Auk: Ornithological Advances 134:894–903. Google Scholar

32.

Miller, A. K., N. J. Karnovsky, and W. Z. Trivelpiece (2009). Flexible foraging strategies of Gentoo Penguins Pygoscelis papua over 5 years in the South Shetland Islands, Antarctica. Marine Biology 156:2527–2537. Google Scholar

33.

Mönkkönen, M., and J. T. Forsman (2002). Heterospecific attraction among forest birds: A review. Ornithological Science 1:41–51. Google Scholar

34.

Morinay, J., J. T. Forsman, and B. Doligez (2020). Heterospecific song quality as social information for settlement decisions: An experimental approach in a wild bird. Animal Behaviour 161:103–113. Google Scholar

35.

Mustafa, O., J. Esefeld, H. Grämer, J. Maercker, M.-C. Rümmler, M. Senf, H.-U. Peter, and C. Pfeifer (2017). Monitoring penguin colonies in the Antarctic using remote sensing data. Umweltbundesamt/German Environment Agency, Germany.  https://www.umweltbundesamt.de/en/publikationen/monitoring-penguin-colonies-in-the-antarctic-using  Google Scholar

36.

Nemeth, E., N. Pieretti, S. A. Zollinger, N. Geberzahn, J. Partecke, A. C. Miranda, and H. Brumm (2013). Bird song and anthropogenic noise: vocal constraints may explain why birds sing higher-frequency songs in cities. Proceedings of the Royal Society B: Biological Sciences 280:20122798. Google Scholar

37.

Nottebohm, F. (1975). Vocal behavior in birds. InAvian Biology ( J. R. King and D. S. Famer, Editors). Academic Press, New York, NY, USA. pp. 287–332. Google Scholar

38.

Pante, E., and B. Simon-Bouhet (2013). marmap: A package for importing, plotting and analyzing bathymetric and topographic data in R. PLoS One 8:e73051. Google Scholar

39.

Patricelli, G. L., and J. L. Blickley (2006). Avian communication in urban noise: Causes and consequences of vocal adjustments. The Auk 123:639–649. Google Scholar

40.

Pertierra, L. R., N. I. Segovia, D. Noll, P. A. Martinez, P. Pliscoff, A. Barbosa, P. Aragón, A. R. Rey, P. Pistorius, P. Trathan, et al. (2020). Cryptic speciation in Gentoo Penguins is driven by geographic isolation and regional marine conditions: Unforeseen vulnerabilities to global change. Diversity and Distributions 26:958–975. Google Scholar

41.

Podos, J. (1997). A performance constraint on the evolution of trilled vocalizations in a songbird family (Passeriformes: Emberizidae). Evolution 51:537–551. Google Scholar

42.

Podos, J. (2001). Correlated evolution of morphology and vocal signal structure in Darwin's finches. Nature 409:185–188. Google Scholar

43.

Ratcliffe, L. M., and P. R. Grant (1985). Species recognition in Darwin's finches (Geospiza, Gould). III. Male responses to playback of different song types, dialects and heterospecific songs. Animal Behaviour 33:290–307. Google Scholar

44.

Rößler, H., M. Lynch, S. T. Ortiz, O. N. Larsen, and M. Beaulieu (2022). Data from: Neighbors matter: Vocal variation in Gentoo Penguins depends on the species composition of their colony. Ornithology 139:ukac031.  https://doi.org/10.5061/dryad.wwpzgmsn1  Google Scholar

45.

Searcy, W. A., J. Soha, S. Peters, and S. Nowicki (2021). Variation in vocal production learning across songbirds. Philosophical Transactions of the Royal Society B 376:0257. Google Scholar

46.

Shera, C. A., J. J. Guinan, and A. J. Oxenham (2002). Revised estimates of human cochlear tuning from otoacoustic and behavioral measurements. Proceedings of the National Academy of Sciences USA 99:3318–3323. Google Scholar

47.

Speirs, E. A. H., and L. S. Davis (1991). Discrimination by Adélie Penguins, Pygoscelis adeliae, between the loud mutual calls of mates, neighbours and strangers. Animal Behaviour 41:937–944. Google Scholar

48.

Suthers, R. A., F. Goller, and C. Pytte (1999). The neuromuscular control of birdsong. Philosophical Transactions of the Royal Society of London, Series B: Biological Sciences 354:927–939. Google Scholar

49.

Swift, R. J., M. J. Anteau, E. A. Roche, M. H. Sherfy, D. L. Toy, and M. M. Ring (2020). Asymmetric benefits of a heterospecific breeding association vary with habitat, conspecific abundance and breeding stage. Oikos 129:1504–1520. Google Scholar

50.

Swift, R. J., A. D. Rodewald, and N. R. Senner (2018). Context-dependent costs and benefits of a heterospecific nesting association. Behavioral Ecology 29:974–983. Google Scholar

51.

The R Foundation (2022). R: The R Project for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria.  https://www.r-project.org/  Google Scholar

52.

Trivelpiece, W. Z., S. G. Trivelpiece, and N. J. Volkman (1987). Ecological segregation of Adelie, Gentoo, and Chinstrap penguins at King George Island, Antarctica. Ecology 68:351–361. Google Scholar

53.

Tyack, P. L. (2020). A taxonomy for vocal learning. Philosophical Transactions of the Royal Society B: Biological Sciences 375:20180406. Google Scholar

54.

Tyler, J., M. T. Bonfitto, G. V. Clucas, S. Reddy, and J. L. Younger (2020). Morphometric and genetic evidence for four species of Gentoo Penguin. Ecology and Evolution 10:13836–13846. Google Scholar

55.

Tyson, C., H. Kirk, A. Fayet, E. E. Van Loon, A. Shoji, B. Dean, C. Perrins, R. Freeman, and T. Guilford (2017). Coordinated provisioning in a dual-foraging pelagic seabird. Animal Behaviour 132:73–79. Google Scholar

56.

Waas, J. R., P. W. Colgan, and P. T. Boag (2005). Playback of colony sound alters the breeding schedule and clutch size in Zebra Finch (Taeniopygia guttata) colonies. Proceedings of the Royal Society B: Biological Sciences 272:383–388. Google Scholar

57.

Wilson, R. P. (2010). Resource partitioning and niche hyper-volume overlap in free-living Pygoscelid penguins. Functional Ecology 24:646–657. Google Scholar
© American Ornithological Society 2022. Published by Oxford University Press for the American Ornithological Society.
Helen Rößler, Maureen Lynch, Sara Torres Ortiz, Ole Næsbye Larsen, and Michaël Beaulieu "Neighbors matter: Vocal variation in Gentoo Penguins depends on the species composition of their colony," Ornithology 139(4), 1-13, (21 July 2022). https://doi.org/10.1093/ornithology/ukac031
Received: 10 February 2022; Accepted: 24 June 2022; Published: 21 July 2022
JOURNAL ARTICLE
13 PAGES


Share
SHARE
KEYWORDS
Antarctic Peninsula
colonias mixtas
ecstatic display call
Gentoo penguins
llamada de despliegue extático
mixed colonies
Península Antártica
RIGHTS & PERMISSIONS
Get copyright permission
Back to Top