Introduction
Spatial distribution of seabirds is often determined by biological and physical attributes of the environment that determine appropriate foraging grounds. These areas can change through time according to modifications of the environment and prey distribution (
Weimerskirch 2007), as well as the stage of the seabirds’ annual cycle (
Herman et al. 2017). Knowing the at-sea distribution of seabirds is important to comprehend their ecological requirements (
Thiebot et al. 2012), their dispersion patterns (
Carneiro et al. 2020), and possible threats in their foraging areas (
Catry et al. 2013). Thus, it provides information that lead us to identify accurately areas to target conservation efforts and for marine spatial planning (
Sansom et al. 2018).
Direct tracking of seabirds provides location information on individuals through time and allows us to develop distribution maps that reflect individuals’ use of space at a specific moment (
Sansom et al. 2018). Moreover, some devices simultaneously record location and dive depth data, combining both the horizontal and the vertical movements of the birds, which is relevant for diving seabirds (e.g.,
Rosciano et al. 2016). Complementary, stable isotope analysis is a reliable tool that provides information on the trophic ecology and assesses patterns of habitat use by organisms (
Hobson 2005;
Newsome et al. 2012). In marine systems, stable isotope values of carbon (δ
13C) reflect primary carbon sources within a food web and can be used to trace trends in marine habitat use by consumers (inshore or benthic vs. offshore or pelagic;
Cherel and Hobson 2007). Stable isotope values of nitrogen (δ
15N) reflect the trophic position of consumers due to a stepwise enrichment of
15N between trophic levels (
Minagawa and Wada 1984).
The largest colony of Southern Rockhopper Penguins (
Eudyptes chrysocome (J.R. Forster, 1781); henceforth Rockhopper Penguins) is located at Isla de los Estados (Franklin Bay), Tierra del Fuego, Argentina (
Schiavini 2000). It is home to approximately 130 000 breeding pairs stretching about 5 km along the western coast (
Raya Rey et al. 2014). Rockhopper Penguins in this colony were tracked during different stages of their annual cycle (inter-breeding and incubation periods), and showed a wide range of use of the South Atlantic Ocean waters, even overlapping foraging ranges with conspecifics from the closest colonies at the Falkland/Malvinas Islands (∼500 km;
Pütz et al. 2002,
2006,
2009;
Raya Rey et al. 2007). During the breeding season, after the chicks have hatched, female Rockhopper Penguins forage at sea and feed their offspring while males remain at the nest to guard them (
Warham 1975).
Here, we aimed to investigate the at-sea distribution of breeding female Rockhopper Penguins during the guard stage. We evaluated this by exploring a multivariate clustering approach for statistically identifying different groups according to the foraging trip and diet (δ
15N) characteristics of each female. Across three breeding seasons, we used a combination of tracking and diving data from GPS and temperature and depth (GPS–TD) data loggers to obtain the track and diving profile of the penguins (e.g.,
Rosciano et al. 2018), and δ
13C and δ
15N stable isotope analysis on whole blood samples to reveal the assimilated diet of the penguins for an extended period of time (
Barquete et al. 2013).
Results
We equipped a total of 36 female Rockhopper Penguins with GPS–TD data loggers over the 3-year study period (2011: n = 13; 2012: n = 10; 2013: n = 13). We obtained a complete data set covering one foraging trip (i.e., diving, GPS, and isotope data) from 24 female Rockhopper Penguins (9 in 2011, 4 in 2012, and 11 in 2013), and incomplete data sets (only GPS and isotope data) for 5 female Rockhopper Penguins in 2012. Differences in data recovered were due to logger failure (3 in 2011, 1 in 2012, and 1 in 2013), and loggers that were lost because penguins did not return to the colony (1 in 2011 and 1 in 2013).
The k-means cluster analysis using k = 3 fitted well to the data (within-group square sum per group: group 1 = 7.31 (n = 4), group 2 = 24.82 (n = 9), and group 3 = 59.18 (n = 11); within-groups square sum/total square sum = 50.4%). When comparing with the result of the grouping factor k = 4, we did not observe a substantial difference in the amount of variability explained (within-groups square sum/total square sum = 57.8%); thus, we chose to use the simplest cluster classification k = 3.
Based on the
k-means cluster analysis, we mapped all individuals by group and identified three different recurrent foraging locations supporting the specified grouping number
k = 3. Cluster 1 included the female Rockhopper Penguins that headed north towards the coast of Tierra del Fuego Island; cluster 2 consisted of individuals that foraged in areas around Isla de los Estados; and cluster 3 consisted of females that headed south towards the continental slope (
Fig. 1). The three groups were well separated by the bathymetry of the respective foraging areas: while individuals in cluster 1 “neritic” remained in waters on the continental shelf (<200 m), individuals in cluster 2 “neritic–oceanic” used both shallow waters around the continental shelf and deeper waters around Isla de los Estados and individuals in cluster 3 “oceanic” foraged in more oceanic waters deeper than 1000 m (
Fig. 1,
Table 1). Females in cluster 1 “neritic” and cluster 3 “oceanic” traveled longer distances from the colony to feed (36.3 ± 21.3 km (cluster 1), 40.3 ± 14.0 km (cluster 3)), though in opposite directions, whereas females in cluster 2 “neritic–oceanic” fed near the colony (16.8 ± 7.8 km) (
Table 1).
Although the incomplete data sets obtained in 2012 had to be excluded (
n = 5) from the
k-means analysis, the GPS data recorded allowed the calculation of the foraging trip duration and thus the association of those trips to the different foraging strategies obtained from the
k-means cluster analysis. Overall, we observed that none of the females used the cluster 2 “neritic–oceanic” strategy in 2011, whereas 8 out of 11 females tracked were found in this cluster in 2013. In 2012, females used mostly the cluster 2 “neritic–oceanic” strategy (4 out of 9) and cluster 3 “oceanic” strategy (3 out of 9) (
Table 2).
Foraging trip duration significantly differed between clusters (
F[2,19] = 8.99,
p = 0.002), with cluster 2 “neritic–oceanic” presenting the shortest foraging trips and cluster 1 “neritic” and cluster 3 “oceanic” revealing comparable foraging trip durations (
Table 1). Females in cluster 3 “oceanic” dived deeper (
F[2,19] = 4.05,
p = 0.0002;
Table 1) and spent more time in the bottom phase per trip hour (
F[2,19] = 7.45,
p = 0.004) compared with cluster 1 “neritic” and cluster 2 “neritic–oceanic” (
Table 1). We did not obtain significant inter-annual differences for any of the variables mentioned above.
The δ
13C values were similar among clusters (
F[2,19] = 2.31,
p = 0.13) and all δ
13C values were lower than −20.0‰ (
Fig. 2,
Table 1). The δ
15N values were significantly higher in cluster 1 compared with cluster 2 and cluster 3 (
F[2,19] = 11.92,
p = 0.0004;
Fig. 2,
Table 1), and also differed significantly between years (
F[2,19] = 4.82,
p = 0.02) with penguins in year 2013 (δ
15N = 8.6‰ ± 0.6‰) having lower values compared with penguins in year 2011 (δ
15N = 9.5‰ ± 0.5‰) and year 2012 (δ
15N = 9.4‰ ± 0.4‰). Chlorophyll
a concentration (Chl
a; mg/m
3) differed significantly between cluster 2 and cluster 3 (
F[2,19] = 8.07,
p = 0.003;
Table 3). Similarly, sea surface temperature (SST; °C) was significantly different between cluster 2 and cluster 3 (
F[2,19] = 9.37,
p = 0.001;
Table 3). Cluster 3 had the lowest Chl
a concentrations and SST. No differences were observed between years either for Chl
a (
F[2,19] = 0.89,
p = 0.43;
Table 3) nor for SST (
F[2,19] = 0.58,
p = 0.58;
Table 3) values. Finally, we observed that chick survival was similar between years (
F[2,19] = 2.45,
p = 0.11;
Table 4) and clusters (
F[2,19] = 0.88,
p = 0.43;
Table 4).
All the cluster × year interactions in the linear models were not significant and thus were discarded to simplify the models.
Discussion
Female Rockhopper Penguins breeding on Isla de los Estados, Argentina, exhibited differential at-sea distribution during their foraging trips at the early chick-rearing period with regards to trip duration, foraging location and associated bathymetry, and dive depth. Spatial segregation in seabirds can be beneficial within a large population as a strategy to relax intra-specific competition (
Wakefield et al. 2011;
Pütz et al. 2018), but may also be linked to environmental characteristics or prey distribution (
Weimerskirch 2007). Furthermore, it could vary in relation to age, sex, reproductive condition and experience, and to some level of individual preferences (
Masello et al. 2013;
Pelletier et al. 2014;
Ceia and Ramos 2015), although this is beyond the scope of this study.
Rockhopper Penguins in this study differed in the areas used to forage (horizontal axis) and in the depth ranges explored (vertical axis). The observed differences in the diving depths may be attributed to a different distribution of the prey and (or) prey sizes within the water column in the distinct foraging areas. Female Rockhopper Penguins that explored greater water depths spent more time in the bottom phase per hour of trip (cluster 3 “oceanic”). Although
Bost et al. (2008) suggested that those parameters combined could be indicating a higher foraging success for those penguins,
Scioscia et al. (2016) showed an association between a greater foraging effort of Magellanic Penguins (
Spheniscus magellanicus (J.R. Forster, 1781)) during their foraging trips (e.g., greater % bottom time) for years in which the consumption of Fuegian sprat (
Sprattus fuegensis (Jenyns, 1842)) (their main prey) decreased. Therefore, further exploration of the effects on the differences presented in the diving profile of female Rockhopper Penguins at the different locations used to forage are desirable (e.g., differences in chick growth and (or) food delivered to chicks). There is evidence for benthic foraging in Rockhopper Penguins (
Ludynia et al. 2013;
Pütz et al. 2018) and for female
Eudyptes chrysocome filholi Hutton, 1879 (
Tremblay and Cherel 2000) during incubation, suggesting that female Rockhopper Penguins may be able to explore the sea floor. However, and in agreement with previous studies at this location, the dives performed by female Rockhopper Penguins in this study were always pelagic (
Schiavini and Raya Rey 2004), since they were too shallow to reach the sea floor, regardless of the foraging area.
The differential at-sea distribution observed (clusters) for female Rockhopper Penguins may also reflect different environmental conditions in the foraging areas and subsequently on the availability of different prey type. For instance, cluster 1 “neritic” and cluster 2 “neritic–oceanic” had similar values for both Chl a concentration and SST, which were higher than those in cluster 3 “oceanic”. Differences in Chl a concentration and SST between years were not statistically significant, and cluster selection by penguins did not depend on year (interaction was not significant), which may indicate that the observed variation in the at-sea distribution of female Rockhopper Penguins during their foraging trips across years may not have been related to annual changes in those environmental characteristics. Moreover, though environmental conditions at sea can rapidly change, female Rockhopper Penguins that were equipped with loggers at the same time (e.g., same day) chose different areas to forage. For example, two females equipped with devices on 9 December 2011 grouped in cluster 1 and cluster 3. In another example, of three females equipped with loggers on 7 December 2012, two females grouped in cluster 3 and one female grouped in cluster 2. Both examples illustrate that the females departed from the colony the same day, but they chose different at-sea areas to forage. This revealed that the areas used to forage did not follow a daily pattern.
Small differences among clusters in δ
15N values suggest that the different foraging areas selected by female Rockhopper Penguin presented different availability of prey types and (or) sizes. Consumption of fish can be associated with higher δ
15N values compared with crustaceans (e.g., euphasiids;
Dehnhard et al. 2016). Thus, the higher δ
15N values observed in cluster 1 “neritic” probably indicate higher consumption of fish (e.g.,
S. fuegensis), whereas the lower δ
15N isotope value in cluster 3 “oceanic” more likely indicate crustaceans as the main prey source (
S. fuegensis usually presents higher δ
15N values compared with crustaceans, e.g., copepods and euphausiids;
Ciancio et al. 2008;
Riccialdelli et al. 2017,
2020). Accordingly, blood samples of Rockhopper Penguin chicks at the Falkland/Malvinas Islands showed that the higher δ
15N values were associated with greater consumption of fish and squid, and lower δ
15N values were associated with a higher prevalence of euphausiids (
Dehnhard et al. 2016).
The available information on prey distribution, abundance, and availability in the study area confirm that it is very productive in terms of biodiversity (
Sánchez et al. 1995;
Hansen 1999;
Ivanovic 2010;
Padovani et al. 2012;
Diez et al. 2016;
Riccialdelli et al. 2020), mainly as a result of the frontal areas that occur in the region (
Acha et al. 2004). For instance, swarms of swarming squat lobsters (
Munida gregaria (Fabricius, 1793)) were recently reported occurring on the continental shelf possibly following a recent expansion process occurring throughout their distribution range (
Diez et al. 2016).
Sprattus fuegensis was also reported in the coastal region off southern Patagonia (
Sánchez et al. 1995;
Hansen 1999); moreover, both species (
S. fuegensis and
M. gregaria) usually overlap spatially and share an ecological niche (
Diez et al. 2018). Some studies also showed high concentration of zooplankton south of 45° (
Sabatini and Álvarez Colombo 2001;
Sabatini et al. 2004;
Romero et al. 2006) coupled with the higher productivity (Chl
a concentration) that occurs during spring and summer (October–March). High heterogeneous densities of amphipods (e.g.,
Themisto gaudichaudii Guérin-Méneville, 1825) associated with the large abundance of zooplankton in the Southern Ocean were also reported for the study area (
Padovani et al. 2012,
2015). The distribution pattern of these prey types is usually patchy and presents peaks according to productivity of the waters. Further studies that focus on prey distribution and availability specifically in the areas that female Rockhopper Penguin use to forage will contribute to our understanding of the differences in the isotope values observed in this study and if those areas actually present differential prey options for this predators.
Although the most parsimonious explanation to the variation in δ
15N values between clusters is related to different prey selected and (or) available in the different foraging areas, we do not fully discard the alternative hypotheses. For instance, the spatial, temporal, and (or) ontogenetic variation in stable isotope values at the base of the food web also has the potential to propagate up the food chain and influence the stable isotope values of marine predators (
Polito et al. 2019). Thus, it may be possible that differences in δ
15N values in each cluster are related to differences in foraging habitats and their baseline δ
15N values. Recent information regarding the variation in the baseline of the food web in the area and the isotopic composition of the prey species for Rockhopper Penguins indicate some degree of spatial variation in δ
13C and δ
15N values at the base of the food web in the study area (
Riccialdelli et al. 2020), thereby emphasizing the need for further studies to fully understand how it may propagate across the food web in the study area.
The δ
13C values in this study were low and similar among clusters, indicating that all individuals fed in pelagic and oceanic environments. This is not totally unexpected, because use of oceanic areas was observed, mainly in female Rockhopper Penguins grouped in cluster 3 “oceanic”. However, we would have expected higher δ
13C values for female Rockhopper Penguins foraging in cluster 1 “neritic” using near-shore foraging habitats, because, for instance, Magellanic Penguins foraging in similar areas presented higher δ
13C values in comparison (
Rosciano et al. 2018). Yet, the stable isotope data in this study integrate primarily the early chick-rearing period, so they may also include a portion from the late incubation period, during which female Rockhopper Penguins breeding at Isla de los Estados were located foraging to the south and east of the Island, reaching oceanic areas near the Polar Front and the Burwood Bank (Area Protegida Marina Namuncura), respectively (
Pütz et al. 2006). These oceanic areas are most likely depleted in
13C (
Cherel and Hobson 2007;
Lara et al. 2010).
Traveling farther to find food for the chicks can be a trade-off for parents: they will need more time to commute between the colony and the foraging areas, which may affect their reproductive success and hence chick survival (
Boersma and Rebstock 2009). This constraint may be compensated with the possibility to prey on higher quality items and (or) more available prey in areas farther away from the colony (
Burke and Montevecchi 2009). To have an overview of the breeding success of the female equipped with GPS loggers, we accounted for chick survival through the early chick-rearing period, and numbers were similar among clusters and years. Therefore, during the guard stage, the length of the foraging trips, in terms of time and distance to the colony, had no effect on chick survival of the birds studied. However, the number of nests followed in our study was quite small; thus, we suggest increasing the sample size in further studies. For example, nesting Rockhopper Penguins at the Falkland/Malvinas Islands colonies showed a growth rate and survival of chicks during the early and late chick-rearing periods that were not affected by the distance from the colony traveled during their foraging trips (
Dehnhard et al. 2016). Also, the quality of the adult breeders and their capacity to efficiently forage under harsh environmental conditions (e.g., prey depletion) can determine the efficiency in energy gain that can be translated into successful reproduction, as was demonstrated for Adélie Penguins (
Pygoscelis adeliae (Hombron and Jacquinot, 1841)) (
Lescroël et al. 2010). Future studies are warranted to further explore the observed pattern in this context.
Different foraging areas and diving patterns within a population can mitigate intra-specific competition (
Grémillet et al. 2004). The high number of penguins nesting in the area (
Schiavini 2000;
Raya Rey et al. 2014) may suggest that choosing different foraging areas could act as a mechanism to reduce intra-specific competition as has been reported between neighboring colonies of the species at other locations (
Masello et al. 2010) as well as other seabird species (
Lewis et al. 2001;
Grémillet et al. 2004). This study provides initial knowledge on the at-sea distribution of Rockhopper Penguins during the early chick-rearing period and allows to identify and to plan how the sample size needs to be increased to answer this question (
Sequeira et al. 2019). Also, individual specialization and (or) consistency among Rockhopper Penguins may reveal other traits associated with the foraging behavior of the species. Individual consistency in seabirds was reported to reduce intra-specific competition (
Ceia and Ramos 2015), revealing the possibility of individual-level differences in foraging strategies (
Dingemanse and Dochtermann 2013). Further studies at the Rockhopper Penguin colony located on Isla de los Estados should focus in understanding this behavior, e.g., deploying loggers to record two or more foraging trips per penguin (e.g.,
Traisnel and Pichegru 2019).
Conclusions
This study provides important data on the spatial distribution of foraging female Rockhopper Penguins during one central phase of their annual cycle, the early chick-rearing period. It enhances our understanding of the distribution of Rockhopper Penguins in the Southwest Atlantic Ocean. It outlines the different foraging areas that female Rockhopper Penguins use during the early chick-rearing period, associated with dive depths explored and most likely different type of prey selected, which may reflect different strategies selected by female Rockhopper Penguins. Although the underlying mechanisms of this differential use of the water masses around Isla de los Estados remain unclear, the prospect of flexibility in the foraging behavior and the availability of several foraging areas within reach of the colony are important features in a scenario of constant environmental changes for a penguin population such as the Rockhopper Penguin (e.g.,
Pütz et al. 2018). Further studies are needed to fully explore individual specialization and (or) consistency among Rockhopper Penguins that may reveal other traits associated with the foraging behavior of the species. This study fills gaps in our knowledge of the spatial distribution of this penguin species during its complete annual cycle (
Pütz et al. 2006,
2018;
Raya Rey et al. 2007), strengthening that female Rockhopper Penguins foraged inside the northern limits of the recently created Marine Protected Area “Yaganes” (particularly, female Rockhopper Penguins in cluster 3 “oceanic”;
Fig. 1) to protect the vast wildlife occurring in the area, preventing fisheries and (or) other extractive activities in the area (
DNC/APN 2017). Thus, this study complements available information on the use of the Southern Ocean by marine mesopredators and top predators and contributes to marine spatial planning in the area.