Evidence of age-related improvement in the foraging efficiency of Adélie penguins

Age variation in reproductive performance is well-documented but the mechanisms underlying this variation remain unclear. Foraging efficiency is likely to be a key source of demographic variation as it determines the amount of energy that can be invested in fitness-related activities. Evidence of age-related changes in the foraging efficiency of adult seabirds is scarce and inconsistent. We investigated the effects of age on the foraging efficiency of breeding Adélie penguins, a relatively short-lived seabird species, in order to gain a broader perspective on the processes driving variation in ageing rates. We found support for a positive effect of age, either linear or levelling off at old ages, on both our proxies for daily catch rate and catch per unit effort. Across all age classes, males were more performant foragers than females. We found no strong evidence for differing ageing patterns between sexes or individual quality levels, and no evidence for senescence. We infer that continuous individual improvement could be responsible for a larger amount of the variation in foraging efficiency with age at our study site, compared with selective disappearance of underperforming phenotypes. The different results reported by other studies highlight the need to conduct longitudinal studies across a range of species in different environments.

and the existence of long-term demographic studies. Nonetheless, evidence for age-related changes in foraging efficiency of seabirds remains relatively scarce and inconsistent across species and geographic locations, due in part to the need to couple long-term demographic data and at-sea monitoring of known individuals. Some studies found evidence for lower foraging efficiency of young, first-time breeders compared to older, experienced breeders 17,18,[22][23][24][25]27 but see 26 on the same species 27 . A handful of studies on long-lived procellariiformes showed that old individuals had a lower foraging performance than younger individuals 17,21 . Recently, studies involving large sample sizes detected no 22 (thick-billed murres Uria lomvia) or very little evidence 23 (wandering albatrosses Diomedea exulans) of age-related variation in foraging behaviour, and no sign of senescence. These two studies are also the only ones, to our knowledge, to have quantified the effect of continuous age (rather than age or experience categories) on the foraging performance of breeding male and female seabirds in the wild.
It has previously been suggested that detecting senescence would be difficult in wild animals because high levels of mortality would remove individuals from the population before senescence becomes evident 10 . However, more recent comparative analyses showed that senescence is widespread in birds and mammals and that ageing rates are determined by ranking on the fast-slow life-history continuum 28 . Hence, species with a slower life history should have a lower senescence rate than species with a faster life history, for a given metabolic rate or body size. Interestingly, most studies on the age-related changes in foraging performance (defined as foraging rate or efficiency) of seabirds where senescence has been detected were conducted on species at the slow end of this fast-slow life-history continuum. Here, we investigate the effect of age on foraging efficiency of a seabird that lives a shorter life in harsher environments, the Adélie penguin (Pygoscelis adeliae, mean life expectancy 12-14 years, Dugger et al. unpubl. data, vs. 33 years for the wandering albatross) 29 . Indeed, Adélie penguins exhibit high fecundity and relatively low adult survival 30 for a pelagic seabird, presumably driven by abundant seasonal food resources and high levels of predation on adults 31 . Some populations range widely at sea 32 , though some to a lesser degree 33,34 and spend most of their lives there, yet often raise two chicks, a life-history pattern that is uncommon among pelagic seabirds 31 . Given their relatively fast life-history, we would expect to see faster aging rates in Adélie penguins compared to other seabirds.
In the present cross-sectional study, we investigated the linear and non-linear effects of age on the foraging efficiency of breeding Adélie penguins, using proxies for daily catch rate and catch per unit effort derived from acceleration and dive data. More specifically, we tested the following hypotheses regarding the shape of the age-specific foraging efficiency curve ( Fig. 1): Foraging efficiency (1) linearly increases with age (maturation or progressive selection of more efficient individuals, no senescence); (2) increases in younger age classes then reaches a plateau (pseudo-threshold response with no senescence); (3) increases in younger age classes, reaches a plateau for middle-aged individuals, then increases in older age classes (terminal investment); (4) increases in younger age classes, reaches an optimum for middle-aged individuals, then decreases in older age classes (senescence); or (5) is constant through age classes. As ageing patterns could differ between birds of different sex or different intrinsic quality, we also included interactions between age and sex, and between age and an age-independent index of breeding quality 35 . Once we had determined the most competitive models for foraging efficiency based on intrinsic characteristics of the birds (age, sex, quality), we modelled the additional effect of extrinsic variables (day in the season, year) to account for environmental variation. Schematic representation of the hypotheses regarding the shape of the age-specific foraging efficiency curve in sexually-mature seabirds. Foraging efficiency 1: linearly increases with age; 2: only increases in younger age classes then reaches a plateau; 3: increases in younger age classes, reaches a plateau for middleaged individuals, then increases in older age classes; 4: increases in younger age classes, reaches an optimum for middle-aged individuals, then decreases in older age classes; 5: is constant through age classes. We investigated different intercepts and slopes (ageing rates) for each of these shapes, depending on sex, breeding quality and environmental conditions.

Methods
Species and Study Site. Adélie penguins are pelagic seabirds, which spend only 10% of their life on land, where they form breeding colonies distributed around the Antarctic coast and high latitude offshore islands 31 . The study colony, Cape Crozier (77°27′S, 169°12′E), Ross Island, is one of the largest for this species (~275 000 pairs) 36 , and therefore exhibits significant intraspecific competition for food (Ainley et al. 2004;Ballance et al. 2009;Ford et al. 2015). Predation by leopard seals (Hyrurga leptonyx) makes foraging activity hazardous for Adélie penguins year-round 37,38 . Accordingly, Adélie penguins exhibit relatively low adult survival rates for pelagic seabirds (68-88% at Cape Crozier) 30 , and delayed maturity, with age at first breeding ranging 3-7 years for females and 4-8 years for males 31 . An improvement of reproductive success with age among young breeders (from 3-to 7-yr-old) was reported from earlier banding data 31 . Banding studies have shown that some individuals reach 20+ years of age, but few individuals reach the mid-teen years 31 (see also Fig. 2).
Breeders arrive at Cape Crozier in late October/early November, lay (usually two) eggs in mid-November, and feed their chicks between mid-December and early February. During the guard stage, one parent remains with the chick(s) while the other forages at sea. Nest reliefs at Crozier occur every 1-2 days during early chick-rearing and chicks are fed relatively small meals by the attending parent. Afterwards, chick demands are too great for adequate provisioning by one parent, so chicks are left on their own ("crèche" stage) while both parents forage simultaneously. We investigated the effect of age on foraging efficiency during the guard stage of four austral summers (2010-2011 to 2013-2014), from mid-December to early January. Hereafter, we refer to austral summers as seasons, using their initial year (e.g., 2010 refers to the 2010-2011 breeding season).
All penguin survey, capture and handling methods used for data collection were performed following all relevant guidelines and regulations under the approval and oversight of Oregon State University's Institutional Animal Care and Use Committee. Additionally, all work was approved and conducted under Antarctic Conservation Act permits issued by the US National Science Foundation and the U.S. Antarctic Program, and administered by H.T. Harvey & Associates.
Population and Individual Monitoring. From 1996 to 2011, 1000 chicks, and 500 chicks from 2012 onwards, have been banded each year on the left flipper with a numbered stainless steel band (except in 2001 when only 110 chicks were banded due to very low reproduction that year; see 39 for details on band design and effect or lack thereof on foraging behaviour). From 1997 onwards the entire colony was searched annually, on 2-to 7-d intervals (depending on the area) for banded individuals throughout the breeding season, 15 December to 25 January, and 15 November to 25 January beginning in 2002. Bands were read with binoculars from a distance (<10 m). When a banded bird was found breeding (i.e. having a nest with at least one egg or chick), its nest contents was subsequently monitored with binoculars every 2-7 d until chicks crèched (about the first week of January).
Instrumentation. Every season, we randomly selected birds of all available ages among this pool of breeders, paying attention to keep a balanced sex-ratio as much as possible and to time the device deployments appropriately so as to always have birds of different ages and different sexes equipped at the same time. This resulted in attaching time-depth recorders (Mk9, Wildlife Computers, Redmond, Washington, USA; 68 × 17 × 17 mm, mass 30 g; hereafter called TDRs) or acceleration-time-depth recorders (G6A, CEFAS Technology Limited, Lowestoft, UK; 40 × 28 × 15 mm, mass 18 g; hereafter called accelerometers) to a total of 45 birds in 2010, 48 birds in 2011, 46 birds in 2012, and 52 in 2013. Sex was assessed in the field based on extensive observations over the lifetimes of the individuals, comparing relative size of especially head and bill, courtship behaviour (male ecstatic displays and copulation), and also timing of attendance in the colony (which is highly synchronous by sex, particularly during incubation and brood stages) 31 . Post-hoc molecular sexing was also conducted as described below. www.nature.com/scientificreports www.nature.com/scientificreports/ From Dec 19-21 to Jan 2-6 (depending on years), birds were caught by hand on their nests and devices were attached to the lower back using black Tesa ® tape (Hamburg, Germany) 40 , following techniques previously shown to have no detectable impacts on foraging trip duration and reproductive success 41 . Handling lasted <10 min and devices were recovered after one foraging trip at sea, lasting between 0.5-9.2 d. Due to device failure or nest failure before the completion of the foraging trip, we recovered a total of 171 files with usable data (47 TDR data files and 124 accelerometer data files from 85 females and 86 males, age 4-17 yrs, Table 1). TDRs were programmed to record pressure (with a resolution of 0.5 m), temperature and light at 1 Hz, while accelerometers recorded pressure (with a resolution of 0.08 m) and temperature at 1 Hz and 3-axis acceleration at 20 Hz. Of the 124 accelerometer data files, 10 of them failed to record acceleration but contained complete pressure and temperature data.
Diving data were processed using the program divesum (v.7.5.5; G. Ballard, unpublished software, see Appendix A in Lescroël et al. 42 for more details on dive-data processing and analysis). This program corrected the record baseline and computed several individual dive parameters such as maximum dive depth, post-dive interval and number of undulations (number of changes in underwater swimming direction from ascent to descent >1 m, involving at least 3 data points, thus 3 sec).
Measuring Foraging Efficiency. We evaluated foraging efficiency in two ways: (1) catch per unit time (trip duration) or daily catch rate, and (2) catch per unit effort (overall dynamic body acceleration).
As a proxy for prey catch, we used the total number of undulations per foraging trip, which has been shown to be linearly related to the number of prey captures 43,44 for both fish-and krill-eating penguins feeding in the water column (as opposed to benthic feeders), such as Adélie penguins. We acknowledge that this proxy did not allow us to quantify the exact number of prey captures, as a small amount of undulations could occur without prey capture and vice-versa 45,46 . It is nonetheless a useful index to compare foraging success between individuals, in terms of relative number of prey capture events, as long as individuals are targeting the same prey in terms of size and energetic content (see Discussion).
Trip duration (from all data files) was measured from entrance in the water to exit, based on the temperature and dive profiles. Trip duration positively correlates with total energy expenditure during foraging by chick-rearing Adélie penguins (r² = 0.88) 47 .
Vectorial dynamic body acceleration (VeDBA, in units of gravitational force g, from accelerometer data only) is a measure of body motion derived from measurements of acceleration made in all 3 spatial axes 48 . This measure of the mechanical work performed by individuals can be used as a proxy for energy expenditure under certain conditions 49 . Raw acceleration values are composed of a dynamic component and a static component. The static acceleration was estimated by smoothing the data in each channel with a running mean over a period of 2 seconds, following methods by 50 . The dynamic acceleration was derived as follows: where acceleration (A) is the sum of static acceleration (S) resulting from body angle with respect to gravity and dynamic acceleration (D) resulting from body movements along all three spatial axis (x: antero-posterior, y: lateral, z: vertical), i.e. A = S + D. The total mechanical power used within the trip is alluded to by measurements of VeDBA summed over the entire trip duration (cf. 48,51 ).
Based on these measures, we evaluated foraging efficiency as: (1) the total number of undulations per trip divided by trip duration to get undulations per day (Uperday), as an index representing catch per unit time; and (2) the total number of undulations per trip divided by total VeDBA to get undulations per g (Uperg), as the catch per unit effort.
Intrinsic and Extrinsic Covariates. We considered age, sex, and an age-independent index of breeding quality as intrinsic covariates that could influence foraging efficiency or the relationship between foraging efficiency and age. Age was determined by linking an individual's cohort-specific band number to our banding database. Sex of the equipped birds was determined from blood or feather samples collected upon detachment of the devices 52 . Using breeding propensity and success data from our monitoring program, we calculated a Breeding Quality Index (BQI) for each known age individual, based on methods we developed previously 35,42 for birds of known breeding history but unknown age. We first calculated a probability of breeding success for each year and individual using four independent variables (age, previous breeding experience, colony of origin and breeding year). The BQI of each individual was then calculated as the mean per individual of the difference between the actual breeding success and the predicted breeding success for every year during which a given individual had been resighted when at least 3 yr-old, up to the year when we equipped it with a TDR or accelerometer. More details about the BQI calculation can be found in 53 . Negative BQI values indicate lower-than-average long-term breeding performance, while positive values indicate above average long-term breeding performance. BQI for the equipped birds ranged −0.32 to 0.93.
We also considered variation of environmental conditions within and between seasons. At Cape Crozier, foraging range increases with time, requiring longer commutes between the colony and high quality prey patches [47][48][49][50][51][52][53][54][55][56] . Our previous work supports the hypothesis that conditions become more challenging as the breeding  www.nature.com/scientificreports www.nature.com/scientificreports/ season progresses due to prey depletion by foraging Adélie penguins and cetaceans near to the colony 54,55,57 . At the same time, because Adélie Penguins are highly synchronous breeders 31 , offspring needs increase consistently among all parents. Thus, as a proxy for the within-season variation in environmental conditions, we used "study day" (a continuous variable with 0 = 20 December) and evaluated its effects on foraging efficiency. We used "year" as a factor to account for any between-season variation, with 2010 as the reference level.
Statistical Analyses. To test our five main hypotheses (Fig. 1), we evaluated models including each foraging efficiency index (Uperday or Uperg) as the dependent variable and either a linear, pseudo-threshold (−1/Age), quadratic (Age + Age²) or cubic (Age + Age² + Age 3 ) effect of age as independent variables, as well as a null model (intercept only). As ageing patterns could differ between birds of different sex or different intrinsic quality, we also included interactions between age and sex, and between age and BQI. Once we had determined the most competitive models for foraging efficiency based on intrinsic characteristics of the birds (age, sex, quality), we added extrinsic variables (study day, year) to the top intrinsic model(s) including potential interactions with age. A null model was also included in this second model set. Residuals were examined to verify normality, homogeneity of variances, and independence.
To evaluate these models and determine the strength of evidence supporting specific effects, we used an information theoretic approach 58 . Models were ranked using the small-sample-size corrected version of Akaike Information Criterion (AICc), with the best model having the lowest AICc value. We calculated ΔAICc as the difference in AICc between each candidate model and the model with the lowest AICc value, and considered all models within 2 ΔAICc as competitive models (Burnham & Anderson, 2002). We determined the strength of evidence supporting specific effects by examining the unstandardized effect sizes (slope coefficients and differences in means) and the associated 85% confidence intervals (CI), as 85% CI are more compatible with an information theoretic approach using AIC or AIC c than 95% CI 59 . If the 85% CI for a parameter in a competitive model (ΔAIC c < 2.0) included zero, it was considered uninformative 59 . Only models performing better than the null model and within 2 ΔAICc of the best model are reported in Tables 2 and 3.
All statistics were performed using R 3.4.0 60 . Acceleration data were analysed with IGOR Pro 6 (WaveMetrics Inc., USA). Means ± SE are given unless indicated otherwise.

Results
Effects of Age on Foraging Efficiency. We found no evidence for a quadratic (i.e., #4 in Fig. 1) or cubic (#3 in Fig. 1) effect of age on either of the foraging efficiency proxies (Tables 1 and 2). Instead, all the top models included either a linear (#1 in Fig. 1) or a pseudo-threshold (#2 in Fig. 1) effect of age on foraging efficiency. For both Uperday (Table 2) and Uperg (Table 3), models including a linear (models 2 and 15) and a pseudo-threshold (models 1 and 16) effect of age were equally competitive, therefore supporting the hypothesis that foraging efficiency increases in younger age classes (maturation or progressive selection of more efficient individuals), then either continues to increase in older age classes or reaches a plateau.
For Uperday, the best model based on intrinsic characteristics only (

Effects of sex and Individual Quality on Foraging Efficiency. Higher foraging efficiency of male
Adélie penguins relative to females was supported in terms of Uperday (Fig. 3) and male Adélie penguins performed 157 more undulations per day (Table 2, model 1, βˆ = 157.2, SE = 61.5, 85% CI: 68.2 to 246.2) than their female counterparts. Sex was also included in 6 of the 9 top models based on intrinsic variables for Uperg ( Table 3)  Higher foraging efficiency of higher quality individuals was not strongly supported. BQI was included in 2 of the 6 top models based on intrinsic variables for Uperday (Table 2) and in 2 of the 9 top models for Uperg (Table 3). However, in all models, the associated 85% CI included zero (e.g.   www.nature.com/scientificreports www.nature.com/scientificreports/ 2012 (βˆ = −1914.9 compared to 2010, SE = 943.6, 85% CI: −3283.0 to −546.9). Interactions between year and age were not supported.
A negative effect of study day on foraging efficiency was included in 4 of the 8 top models including extrinsic variables for Uperday (Table 2, models 7-10) and 2 of the 4 top models for Uperg (

Discussion
We found support for a positive effect of age, either linear or levelling off at old ages, on both indices of foraging efficiency in breeding Adélie penguins. Younger birds exhibited lower foraging efficiency, both in terms of prey capture per unit of time (Uperday) and per unit of dynamic acceleration (Uperg). Other variables had different effects on these two indices. Across all age classes, males showed higher foraging efficiency than females in terms of Uperday but this effect was weaker for Uperg. Uperg was marked by strong inter-annual variations, while Uperday was more driven by an increase in trip duration throughout the breeding season (study day effect). Overall, extrinsic variables, year and study day, tended to explain more of the variation in foraging efficiency than age.
Individual Improvement and Selective Disappearance of Phenotypes. Our results are consistent with both hypotheses (Hyp. 1 and 2) involving a continuous improvement of foraging efficiency and/or the progressive selection of more efficient individuals, possibly levelling off in older age classes. These two mechanisms will be difficult to tease apart and to quantify without a long-term, longitudinal study of foraging efficiency. However, it is very likely that both individual improvement and selective disappearance of individuals from the population play a role in the observed pattern. Adélie penguins live in an extremely harsh and selective environment and experience high predation rates 31 . In a previous study 35 , we showed the existence of demographic heterogeneity (i.e. some types of individuals having greater or lesser chances of surviving or reproducing) 61 in the Cape Crozier population. We also showed that the first reproductive event was likely acting as a selective process leading to a more homogeneous pool of experienced breeders, highlighting the fact that some selective disappearance is indeed happening, at least during the early stages of reproductive life. On the other hand, recent studies on reproductive performance accounting for unobserved heterogeneity (see 62 and references therein) showed that average individual improvement may be a widespread phenomenon in birds and mammals. In their study on common terns, Rebke et al. 62 suggested that the observed age-specific patterns in average improvement in reproductive performance among survivors between age 3 and 14 yrs may be due to a gain in experience, possibly through improved foraging ability. Clearly, in brown pelicans, older birds have fewer accidents, i.e. taking water depth into consideration before plunging into the ocean, than do naïve ones 27 . Our results on Adélie penguins confirm the possibility of continuous improvement of foraging ability through old age classes.
We did not find strong support for a direct relationship between breeding quality and foraging efficiency. This is in accordance with our previous results, showing that this relationship is only apparent under harsh environmental conditions 42 , with this current study not including years of conditions as harsh as those encountered in 2001-2005, when giant icebergs were present in the foraging area. Our focus on breeding individuals, i.e. only a portion of the diversity of phenotypes at fledging, also likely restrained the range of inter-individual variation. Additionally, we found no evidence for different age-related variations according to breeding quality, i.e. low and high-quality breeders do not age differently, at least as far as foraging efficiency is concerned. Therefore, it seems that under "normal" environmental conditions, individuals of higher intrinsic quality do not exhibit significantly higher foraging efficiency than lower quality individuals or maintain better their foraging efficiency into old ages (no effect of the interaction between age and BQI). Altogether, these results suggest that within-individual improvement, more so than selective disappearance of underperforming phenotypes and at least in older ages, may be responsible for a large amount of the observed increase in foraging efficiency with age. Interestingly, a recent longitudinal demographic analysis using long-term data gathered at all 3 colonies of Ross Island from 1996 to 2013 (Kappes et al. In prep) does show that within-individual improvement is the main mechanism underlying patterns of age-related changes in reproductive performance at Cape Crozier. Given that the average age at first reproduction (shown to act as a selective process) is just under 6 yr old at Cape Crozier (Kappes et al. In prep), and that most of the observed variation in Uperday (but not in Uperg) occurs up to this age, it is possible that selective disappearance could be at least partially responsible of the observed improvement of foraging efficiency in early ages.
Behavioral senescence. As in the most recent studies on age-related variation in foraging behaviour of seabirds 22,23 , we did not find any sign of senescence on global indices of foraging efficiency. On thick-billed murres, Elliott et al. 22 were able to find signs of physiological senescence in parameters linked to diving capacity (haematocrit, resting metabolic rate and thyroid hormone levels), while Cunningham et al. 63 detected signs of deterioration of the diving performance at the scale of the dive. The fact that it was possible to detect early signs of senescence at a very fine scale but not on the resulting foraging efficiency means that either old birds are able to adjust their behaviour to compensate for decreased fine-scale performance and/or that they die shortly after the onset of senescence. In common terns, only 20% of all mature birds reached the age when reproductive senescence may become important 62 . In Adélie penguins, high mortality levels due to predation 37 and harsh environmental conditions may explain why we could not detect signs of senescence in foraging efficiency.
Sex Differences in Foraging Efficiency. Sexual differences in foraging behaviour and provisioning rates have been reported for many seabird species, including species with no or little sexual dimorphism such as Adélie www.nature.com/scientificreports www.nature.com/scientificreports/      www.nature.com/scientificreports www.nature.com/scientificreports/ penguins. Such differences have been explained by a sex-specific allocation of foraging effort between parents and offspring 64 , differences in time spent defending the nest [65][66][67] , or spatial segregation due to intra-specific competition or habitat selection 68 . Here and in previous studies 42 , e.g. 69 , we showed that female Adélie penguins from Cape Crozier stayed longer at sea for the same amount of food caught (using the number of dive undulations) or brought back to the colony (using an automated weighing system). This pattern was also found at Béchervaise Island, where female Adélie penguins made longer foraging trips than males, ranged greater distances more frequently and consumed larger quantities of krill 70 . In Adélie Land, however, Angelier et al. 71 did not find any significant sex difference in foraging success, trip duration, or maximal foraging range, although females tended to range farther than males (94.72 ± 6.44 vs. 65.20 ± 16.46 km).
In our study colony, males bring back the same amount of food as females after shorter foraging trips but do not make more frequent trips or feed their chicks more frequently 42 . It seems therefore unlikely that the observed sex differences would stem from a sex-specific allocation of foraging effort between parents and offspring. Rather, it appears that male and female foraging areas are spatially segregated, with females going further away from the colony on longer trips and diving shallower 42 , while males exploit waters closer to the colony by diving deeper, thus allowing them to spend more time defending their territory. Interestingly, we found a much weaker effect of sex on Uperg than Uperday, suggesting that if the male strategy was indeed more efficient in terms of number of prey captures per unit of time spent at sea (Uperday), it was about as costly in terms of mechanical work (deep dives for males vs. longer trips for females). In an optimal foraging theory framework 72 , it seems that male Adélie penguins tend to maximize net energy gain per unit time, while females tend to maximize net energy gain per energy spent. However, more work is needed to examine fine-scale foraging habitat and diet preferences of males and females so as to determine whether this difference in foraging habitat use is based on based on niche specialization 68 , competitive exclusion 73 , or differences in nutritional requirements 74 .
We also acknowledge that our index of foraging success, i.e. the number of undulations in the dive profile, while linearly related to the number of prey captures in penguins, is an imperfect measure and could represent vastly different energy intake depending on whether it relates to crystal krill capture (Euphausia crystallorophias, average individual wet weights for subadults to adults:0.02 to 0.2 g, average individual energy content: 0.05 to 0.46 kJ) [75][76][77] or silverfish (Pleuragramma antarcticum, average individual wet weight for older subadults: 3.42 g, average individual energy content: 8.92 kJ) 75,78 , the two main prey of Adélie penguins in the southern Ross Sea 79 . An individual foraging primarily on silverfish might undulate significantly less (in another study, Watanabe & Takahashi 46 showed that Adélie penguins catch about 7 krill individuals for 1 fish Pagothenia borchgrevinki) and still gain the same amount of energy. However, krill is by far the dominant item of the diet of breeding Adélie penguins during the early chick-rearing season 57,75 , when we equipped them in this study. Additionally, the average number of undulations per trip is 2403 ± 1151 (this study) and average food load (i.e. food brought back to the chicks by a parent) is 592 ± 285 g 50 . This translates into 0.25 g per undulation. If we account for 15% of false positives (i.e. undulations without actual prey capture) 78 , it gives 2043 successful undulations and 0.29 g per undulation. Both estimates are consistent with a krill-dominated diet.
Drivers of Foraging Efficiency. In our study, intrinsic characteristics of the birds, such as age, sex, or breeding quality explained a maximum of 6% (for Uperg) to 10% (for Uperday) of the variability in foraging efficiency (3 to 5% for the age only), compared to 9% (for Uperday) to 18% (for Uperg) for extrinsic variables (year, study day) only (Tables 2 and 3). We confirmed that foraging efficiency of birds of both sexes decreased as the season progressed, likely related to resources becoming depleted or inaccessible close to the very large Cape Crozier colony 54,55,57 . We identified large inter-annual variations in average foraging efficiency as proxied by Uperg, which may at least partially be driven by inter-annual variations in sea ice concentration (SIC) within the summer foraging area of Adélie penguins. At our study colony, maximum foraging efficiency during summer is known to occur at relatively low SIC, peaking at 6-12% and decreasing with higher SIC 80 . In 2011, the average SIC over monitored penguin foraging days was especially high (18.4%) and coincided with low foraging efficiency as well as an attenuation of the age effect. Under high SIC (13-38%), foraging efficiency is constrained to low levels 80 as penguins have to walk over extended, continuous stretches of ice, or hop from floe to floe to avoid being crushed or attacked by predators hunting in the leads. In these conditions, foraging efficiency becomes essentially uncoupled from SIC variations as phenotypic plasticity is suppressed 80 , older age does not seem to provide a significant advantage anymore, but inter-individual differences linked to breeding quality are exacerbated 50 . Interestingly, the fact that Uperg was marked by stronger inter-annual variations than Uperday suggests that Adélie penguins first respond to more difficult environmental conditions by intensifying their foraging effort before increasing their trip duration, which can only be stretched to a given extent before chick survival is jeopardized.
As the amount of variation explained has rarely been reported in previous studies about the effect of age on foraging efficiency, it is difficult to put our results in perspective. As an indication, Elliott et al. 22 report R² of 0.19 to 0.29 for the effect of age on physiological parameters in thick-billed murres, and Cunningham et al. 63 report R² of 0.14 to 0.46 for the effect of age on energy expenditure (reported as Partial Dynamic Body Acceleration) during different parts of the dive in the same species. Compared to these, the 3-5% variation attributed to age in this study appears to be small but is in accord with behavioural traits generally being more labile than physiological or morphological ones 81 . Seabird behaviour in particular is sensitive to the variations in environmental conditions (hence their use as ecological indicators 82 and also very plastic so that shifts in phenotypic plasticity can affect or blur the functional relationships between behavioural traits and environmental variables 80,83 . It is therefore not too surprising that age would account for only a small part of the variation in foraging efficiency, especially measured over 4 different years. Assuming that one undulation represents the capture of 0.29 g of food (at the minimum, not taking into account prey digestion during the early part of the trip), then age alone explains differences in prey catch up to 86.71 g (or 14.7% of the total food load) per day between the youngest and oldest birds in our www.nature.com/scientificreports www.nature.com/scientificreports/ study. Under competitive conditions at large colonies such as Cape Crozier 47,54,57 , even a small advantage due to age or intrinsic "quality" may make a difference in terms of reproductive success and/or survival.

Data Availability
Data collected are available at California Avian Data Center (CADC) hosted by Point Blue Conservation Science and metadata registered with the "Antarctic Master Directory" (http://gcmd.nasa.gov/KeywordSearch/Home. do?Portal=amd&MetadataType=0). Data are and will be available at CADC: http://data.prbo.org/apps/penguinscience/.