Vol. 8: 269–277, 2010 AQUATIC BIOLOGY Published March 16 doi: 10.3354/ab00213 Aquat Biol

Contribution to the Theme Section ‘Ecophysiology and biomechanics in diving seabirds’ OPENPEN ACCESSCCESS Dive efficiency versus depth in foraging emperor

Ilka Zimmer1,*, Rory P. Wilson2, Michaël Beaulieu3, Yan Ropert-Coudert3, Akiko Kato4, André Ancel3, Joachim Plötz1

1Alfred Wegener Institute for Polar and Marine Research, PO Box 120161, 27515 Bremerhaven, Germany 2Institute of Environmental Sustainability, School of Biological Sciences, Swansea University, Singleton Park, SA2 8PP Swansea, UK 3IPHC—Département Écologie, Physiologie et Éthologie—UMR7178 Centre National de la Recherche Scientifique and Université Louis Pasteur, 23 rue Becquerel, 67087 Strasbourg, France 4National Institute of Polar Research, 1–9–10 Kaga, Itabashi, Tokyo 173–8515, Japan

ABSTRACT: Dive duration generally decreases with dive depth in air-breathing vertebrates. In most species, this occurs due to increasing transit (descent and ascent) durations. However, in U-shaped dives, this is also because the duration of the bottom phase of the dive increases with increasing depth. We considered the time-based efficiency with which can use different diving depths by using depth recorders on 9 foraging emperor penguins Aptenodytes forsteri, studied dur- ing the early and late chick-rearing period in Adélie Land, Antarctica. Dive and post-dive duration positively correlated with diving depth, but particularly long recovery durations were apparent for dives longer than 456 s. This inflection point (i.e. 456 s) corresponded to a theorized maximum rate of gain of prey per dive cycle. By using the number of undulations in the bottom phase of the dive as a proxy for prey capture success, we conclude that the most lucrative dive depths for the birds studied were between 50 and 225 m. Since these depths were also those most often visited, we think that foraging emperor penguins focus on depths where profitability is highest.

KEY WORDS: Emperor penguins · Foraging depth · Bottom dive duration · Dive efficiency

Resale or republication not permitted without written consent of the publisher

INTRODUCTION brood survival (cf. Takahashi et al. 2003). The specific behaviour of a diving penguin (e.g. swim speed, rate of Much of the optimal foraging literature deals with change of depth) determines the rate of oxygen use how animals should behave to maximize their rate of underwater (e.g. Kooyman & Ponganis 1994) and sub- energetic gain. The modulation of factors such as the sequently when the is at the surface between time and energy involved in transiting between resting dives (Kooyman & Ponganis 1998). Both physiological and foraging sites plays an important role in this (e.g. processes affect the amount of time that can be spent Shepard et al. 2009 and references therein). Penguins at depth and consequently define the ultimate effi- provisioning chicks are a good example of this: there is ciency of the penguin for hunting prey at depth. considerable selection pressure for them to organise Emperor penguins Aptenodytes forsteri are the the durations and depths underwater and the subse- deepest diving birds, with a maximum recorded dive quent time spent recovering at the surface (Butler & depth of 564 m (Wienecke et al. 2007). This diving abil- Jones 1997), so that they maximize the rate of energy ity allows them to prey on squid and at great acquisition (Weimerskirch 1998, Ropert-Coudert et al. depths in the Southern Ocean where no other avian 2004). This rate of energy acquisition is ultimately species can compete (Ropert-Coudert et al. 2006a). important in chick provisioning rates, and therefore Deeper dives, however, require longer (and energeti-

*Email: [email protected] © Inter-Research 2010 · www.int-res.com 270 Aquat Biol 8: 269–277, 2010

cally more costly) transit phases to depth and may then travelling to or from the colony (when the colony to sea require longer surface recovery time (e.g. Kooyman & distance was about 500 m). The sex of the birds is diffi- Ponganis 1998) and result in poorer diving efficiency cult to determine during this period. Consequently we (Wilson 2003). Such a reduced diving efficiency might inferred sex for one female and one male by their be compensated, however, by appropriately greater vocalisations (Robisson 1992). We attached a satellite- gains in prey acquisition rates at those depths. linked, data-collecting tag (Splash, Wildlife Comput- We examined the foraging behaviour of breeding ers, 7.8 × 5.0 × 2.3 cm, 105 g) with black Tesa-tape emperor penguins. We analysed the depths at which (Wilson et al. 1997b) on the lower back feathers of the these birds operated, assessed how diving efficiency penguins. varied according to foraging depth, and considered The Mk9 archival tag had a relative depth resolution whether our depth-dependent, derived efficiencies of 0.5 m, a memory of 16 MB and was set to record correlated with putative prey ingestion rates and depth every 5 s. The Splash tag had a relative depth whether the birds maximized rate of prey ingestion resolution of 0.1 m, a 14 MB memory, and was set to versus time. record depth every 2 s. Both devices recorded within the range of 0 to 1000 m depth. Data analysis. We analysed dive depth data (cor- MATERIALS AND METHODS rected for surface drift which varied by ± 2 m) using MT-dive, software specially developed for this pur- Study site and data records. We conducted field- pose (Jensen Software). The software analysed dives work at the emperor penguin colony at Pointe Géo- sequentially, writing a number of defining parameters logie, Adélie Land, Antarctica (66° 40’S, 140° 01’E) dur- for each dive into an ASCII output file. These were: the ing the 2005 breeding season, between 20 April and time of the dive initiation; the overall dive duration; the 31 August in winter and 2 and 21 November in spring. maximum depth reached during the dive; the descent-, We attached time-depth recorders (Mk9, Wildlife bottom-, and ascent-phase durations; the vertical velo- Computers, 6.7 × 1.7 × 1.7 cm, 30 g) to 3 female and 2 cities during the descent, bottom and ascent phases; male penguins between 20 April and 11 May, in the the number of rapid succession short ascent/descent pairing period. We glued the instruments to the bird’s phases during the bottom phase; and the post-dive feathers on the midline of the lower back to minimize interval (Fig. 1). drag (Bannasch et al. 1994). We considered dives deeper than 2 m as proper Between 2 and 13 November in the late chick rear- dives for further dive analysis. We defined the bottom ing period, a further 4 adult penguins were captured phase of each dive, during which penguins are most

0 Start of dive End of dive

SPI counted by dive analysis 50

2 to 3 serial points of inflection (SPI) = 1 wiggle ~ 1 prey pursuit 100

150 Dive depth (m)

200 Descent phase Bottom phase Ascent phase

0123 45678 Dive duration (min) Fig. 1. Aptenodytes forsteri. Example of a dive classified into descent, bottom and ascent phase. A wiggle during the bottom phase is defined by 2 to 3 registered serial points of inflection (SPI) counted during dive analysis Zimmer et al.: Dive efficiency versus depth in emperor penguins 271

likely to hunt (e.g. Wilson et al. 1995) and appear to with post-dive durations of less than 25 min (as delin- capture most of their prey (e.g. Simeone & Wilson eated by a point of inflection in the frequency of post- 2003, Ropert-Coudert et al. 2006b), by 3 conditions: it dive duration graph; cf. bout-ending criterion Gentry could only occur (1) at depths >85% of the maximum & Kooyman 1986). We processed data using Origin depth of the dive (dive analysis package DA, Wildlife 7.5, TableCurve 2D and JMP 6.0, SAS 2005. Computers); (2) if it was bounded by 2 points of inflec- tion in the rate of change of depth (indicating the start and end of a clear bottom period); and (3) if the overall RESULTS rate of change of depth for the whole of the putative bottom period did not exceed 0.2 m s–1 (Rodary et al. We analysed dive records of 9 emperor penguins, 2000a). Short ascent and descent phases that exceeded consisting of one complete foraging trip for each. We 2 m during the bottom phase of a dive were quantified also analysed dives of a second foraging trip for one of according to the number of points of inflection during the females in winter. Those trips provided 277 d of the ascents and descents (Fig. 1). We described 2 or foraging time at sea and 24 403 dives. The maximum 3 serial points of inflection (SPI) as a ‘wiggle’. Single dive depth was 438 m (mean: 68.83 m ± 0.52, median: wiggles of this type reported for Magellanic penguins 32.76 m, range: 2.00–438.37 m). The maximum dive Spheniscus magellanicus have been shown, overall, to duration was 16 min (mean: 3.00 min ± 0.02, median: result in the capture of a single prey item during the 2.67 min, range: 0.10–16.08 min). bottom phases of their dives (Simeone & Wilson 2003). We examined dive duration and post-dive duration Moreover, wiggles are generally considered to be indi- for all dives that had post-dive intervals <25 min (n = cative of prey pursuit in penguins (e.g. Kirkwood & 24 068 dives). Both parameters were positively corre- Robertson 1997, Rodary et al. 2000b). Consequently, lated with maximum dive depth (log dive duration: 2 we considered the number of wiggles occurring in the r = 0.88, F1,24063 = 165011.4, p < 0.001 [Fig. 2a] and log 2 bottom phase of emperor penguin dives to be linearly post-dive duration: r = 0.40, F1,24063 = 13551, p < 0.001 related to prey abundance in that they represented a [Fig. 2b]) and there were no significant differences be- pursuit with, or without, prey capture. We note that the tween seasons (log dive duration: F2,6 = 1.59, p = 0.28, validity of this assumption is critical to our assessment and log post-dive duration: F2,6 = 2.48, p = 0.16). of ‘catch per unit effort’. Our estimate of prey abun- The rate of increase of the post-dive duration was dance via ‘catch per unit effort’ (CPUE, SPI min–1) sta- higher for deep dives between 225 to 450 m, being tistics was obtained by dividing the number of serial 1.7 s m–1, than for shallower dives between 0 to 225 m, points of inflection during the bottom phase by the bot- being 0.35 s m–1 (Fig. 2b). Log-transformed post dive tom dive duration. Likewise, the diving efficiency (DE) duration and dive duration were positively correlated 2 was calculated following the equation of Ydenberg & (r = 0.37, F1,24062 = 11759.6, p < 0.001 [Fig. 3a]) and not Clark (1989): different between seasons (F2,6 = 3.31, p = 0.11).

DE = bottom duration/ (1) (dive duration + post-dive duration) being expressed as a non-dimensional fraction. As an 600 Total dive duration Post-dive duration overall measure for DE and CPUE, the product of these parameters was calculated as the CPUE×DE index, 500 being expressed as SPI min–1. We analysed dive data for each data set. The 2 dif- 400 ferent study periods of winter and spring in the 300 breeding cycle of emperor penguins resulted in dif- ferent foraging trip characteristics (e.g. trip direction Duration (s) 200 and range, trip duration and, thus, the overall num- ber of dives conducted; see Zimmer et al. 2007). We 100 log-transformed results of the dive parameters (e.g. as dive durations, dive efficiency and the CPUE, 0 being variables that change with depth) in order to 0 50100 150 200 250 300 350 400 450 satisfy normality rules before being investigated Dive depth (m) using General Linear Models with bird as the random Fig. 2. Aptenodytes forsteri. Dive and post-dive duration as factor. For clarity, the data are presented as grand function of maximum dive depth for dives with post-dive means ± 1 SE and averaged per 25 m maximum dive durations <25 min; n = 24 068 dives; n = 9 penguins. Values depth interval. We calculated all parameters for dives are grand means ± 1 SE 272 Aquat Biol 8: 269–277, 2010

y = 28.24 + 3.46 x 10–5x2.5 a 250 y = 0.69 + 0.01 (max. depth) 600 r2 = 0.98 r2 = 0.97 p < 0.001 p < 0.0001 F 200 F = 567.95 500 12 = 455.46 17 y = 175.21 + 1.33 x 10–72x26.9 r2 = 0.99 150 400 p < 0.001 F 5 = 241.61 300 100

200 50 Post-dive duration (s)

100 Descent 0

0 0 50100 150 200 250 300 350 400 450 0 100200 300 400 500 600 Dive duration (s) b Fig. 3. Aptenodytes forsteri. Post-dive duration as function of 100 dive duration for dives with post-dive durations <25 min; n = 24 068 dives, n = 9 penguins. This relation is described by two different correlations due to a point of inflection at 455 ± 9 s, 80 causing a change of the slope. Values are grand means ± 1 SE 60 We visually determined a point of inflection at 456 ±

9 s from plots of dive duration in relation to post-dive Duration (s) 40 duration (Fig. 3a) and tested the relationship sepa- rately before and after this same point of inflection. For 20 the shorter log dive durations (dives <456 s) the log 2 Bottom post-dive duration was positively correlated (r = 0.36, 0 F = 11189, p < 0.001) and did not differ between 1,23389 0 50100 150 200 250 300 350 400 450 seasons (F2,6 = 2.59, p = 0.15 [Fig. 3b]). For the longer log dive durations (of dives >456 s) the log post-dive 300 duration was also positively correlated but with a c y = 0.75 + 0.01 (max. depth) 2 2 higher rate of change with depth (r = 0.16, F1,663 = r = 0.96 250 p < 0.0001 18.7, p < 0.001) while, similarly, not differing between F 17 = 406.29 seasons (F2,5 = 3.08, p = 0.13 [Fig. 3c]). Duration of descent and ascent was positively corre- 200 2 lated to maximum depth (r = 0.85, F1,24063 = 130625, 2 150 p < 0.001 [Fig. 4a] and r = 0.84, F1,24063 = 116114, p < 0.001 [Fig. 4c]) and not different between seasons 100 (F2,6 = 0.035, p = 0.97 and F2,6 = 1.24, p = 0.35). Eighty percent of 24 068 dives (with post-dive inter- 50 vals <25 min; n = 19 248 dives) had a bottom phase. Ascent Bottom duration was shorter than descent and ascent 0 during dives and accounted for 30 ± 0.4% of the dive duration on average overall. Bottom duration in- 0 50100 150 200 250 300 350 400 450 creased with increasing maximum dive depth for dives Maximum dive depth (m) 2 <225 m (r = 0.61, F1,17621 = 26134, p < 0.001 [Fig. 4b]) Fig. 4. Aptenodytes forsteri. Durations of (a) descent, (b) bot- and was not different between seasons (F2,6 = 2.12, p = tom, and (c) ascent phases of dives as function of maximum 0.20), but then decreased with dives >225 m (r2 = 0.14, dive depth for dives with post-dive durations <25 min; n = 24 068 dives, n = 9 penguins. Values are grand means ± 1 SE F1,1576 = 49.5, p < 0.001 [Fig. 4b]), again without sea- sonal difference (F2,2 = 1.20, p = 0.45). Dives with serial points of inflection (SPI dives, n = apparent peaks, one for dives with maximum depths 13 915), which suggest prey pursuit movements, oc- between 25 and 75 m (23% of dives) and another for curred in 72% of dives with a bottom phase. The fre- dives with maximum depths between 125 and 225 m quency distribution of all dives with SPIs showed 2 (41% of dives, Fig. 5a). Zimmer et al.: Dive efficiency versus depth in emperor penguins 273

a 0.30 12 0.25 10 0.20 8

0.15

6 DE

4 0.10 SPI frequency (%) 2 0.05

0 0.00 0 50100 150 200 250 300 350 400 450 0.0 0.51.0 1.5 2.0 2.5 3.0 3.5 4.0 CPUE b 4 Fig. 6. Aptenodytes forsteri. Dive efficiency (DE) as function of catch per unit effort (CPUE, SPI min–1) for dives with post- dive durations <25 min and bottom phase; n = 19 248 dives,

) 3 (n = 9 penguins). Values are grand means ± 1 SE –1

tionship was different between seasons (F = 9.17, p = 2 2,6 0.016). The dive efficiency (DE) only increased with maxi- mum depth for dives up to 50 m before decreasing CPUE (SPI min 1 substantially with increasing dive depth (r2 = 0.011,

F1,19239 = 55.0, p < 0.001) to a minimum mean value of 0.0