Diving Behavior and Fine-Scale Kinematics of Free-Ranging Risso's Dolphins Foraging in Shallow and Deep-Water Habitats

Air-breathing marine predators must balance the conflicting demands of oxygen conservation during breath-hold and the cost of diving and locomotion to capture prey. However, it remains poorly understood how predators modulate foraging performance when feeding at different depths and in response to changes in prey distribution and type. Here, we used high-resolution multi-sensor tags attached to Risso’s dolphins (Grampus griseus) and concurrent prey surveys to quantify their foraging performance over a range of depths and prey types. Dolphins (N=33) foraged in shallow and deep habitats [seabed depths less or more than 560m, respectively] and within the deep habitat, in vertically stratified prey features occurring at several aggregation levels. Generalized linear mixed-effects models indicated that dive kinematics were driven by foraging depth rather than habitat. Bottom-phase duration and number of buzzes (attempts to capture prey) per dive increased with depth. In deep dives, dolphins were gliding for >50% of descent and adopted higher pitch angles both during descent and ascents, which was likely to reduce energetic cost of longer transits. This lower cost of transit was counteracted by the record of highest vertical swim speeds, rolling maneuvers and stroke rates at depth, together with a 4-fold increase in the inter-buzz interval, suggesting higher costs of pursuing and handling prey compared to shallow-water feeding. In spite of the increased capture effort at depth, dolphins managed to keep their estimated overall metabolic rate comparable across dive types. This indicates that adjustments in swimming modes may enable energy balance in deeper dives. If we think of the surface as a central place where divers return to breathe, our data match predictions that central place foragers should increase the number and likely quality of prey items at greater distances. These dolphins forage efficiently from near-shore benthic communities to depth-stratified scattering layers, enabling them to maximize their fitness.

Air-breathing marine predators must balance the conflicting demands of oxygen conservation during breath-hold and the cost of diving and locomotion to capture prey. However, it remains poorly understood how predators modulate foraging performance when feeding at different depths and in response to changes in prey distribution and type. Here, we used high-resolution multi-sensor tags attached to Risso's dolphins (Grampus griseus) and concurrent prey surveys to quantify their foraging performance over a range of depths and prey types. Dolphins (N = 33) foraged in shallow and deep habitats [seabed depths less or more than 560 m, respectively] and within the deep habitat, in vertically stratified prey features occurring at several aggregation levels. Generalized linear mixed-effects models indicated that dive kinematics were driven by foraging depth rather than habitat. Bottom-phase duration and number of buzzes (attempts to capture prey) per dive increased with depth. In deep dives, dolphins were gliding for >50% of descent and adopted higher pitch angles both during descent and ascents, which was likely to reduce energetic cost of longer transits. This lower cost of transit was counteracted by the record of highest vertical swim speeds, rolling maneuvers and stroke rates at depth, together with a 4-fold increase in the inter-buzz interval (IBI), suggesting higher costs of pursuing, and handling prey compared to shallow-water feeding. In spite of the increased capture effort at depth, dolphins managed to keep their estimated overall metabolic rate comparable across dive types. This indicates that adjustments in swimming modes may enable energy balance in deeper dives. If we think of the surface as a central

INTRODUCTION
When animals are foraging, their efficiency can be defined as the difference in energy gained from ingesting prey relative to the energy expenditures associated with searching for and capturing prey (Parker et al., 1996). Therefore, to increase efficiency, animals must minimize the cost of prey capture (Williams et al., 2000) and/or increase the energetic benefits from prey (Watanabe and Takahashi, 2013;Watanabe et al., 2014). To minimize locomotor costs, aquatic animals adopt gliding gaits, and slow swimming speeds to decrease drag and the energetic requirements of searching for food (Williams, 2001;Fahlman et al., 2008Fahlman et al., , 2013Watanabe et al., 2011). To increase energy intake, they can increase feeding rates by extending the time spent in a prey patch (Doniol-Valcroze et al., 2011;Wilson et al., 2011), or selecting prey patches of higher quality (Orians and Pearson, 1979). Higher quality prey can take the form of larger bodied prey for animals that take prey individually (i.e., particle feeders, hunting predators), greater densities of the same type for multiple-prey loaders (i.e., filter feeders or grazers), or other prey changes that decrease handling time or increase catch per unit effort. Moreover, the cumulative benefit from a foraging bout may increase due to a combination of these processes and mechanisms (Daniel et al., 2008;Doniol-Valcroze et al., 2011;Gibb et al., 2016). However, what remains poorly understood is how predators modulate foraging performance when faced with heterogeneous, dynamic prey-scapes .
Many animals across taxonomic groups and environments are central place foragers where feeding bouts occur far from a return location (Houston and McNamara, 1985). Air-breathing divers can be described as central place foragers, where the sea surface is the return location that provides a critical source of oxygen replenishment (Kramer, 1988). In this context, divers face the conflicting demands of oxygen conservation during apnea and the energetic costs associated with diving, swimming, and feeding (Davis et al., 2004). Therefore, diving foragers must balance managing oxygen stores while maximizing net energy intake during feeding at depth (Hazen et al., 2015). If the high energetic costs of feeding deplete body oxygen stores, this in turn would reduce aerobic diving capacity and thus impose a potential constraint on foraging performance (Croll et al., 2001;Acevedo-Gutiérrez et al., 2002;Goldbogen et al., 2012).
Central place foraging theory predicts that predators should only swim farther from the central place to encounter increasing numbers or quality of prey (Carbone and Houston, 1996). When the same prey resource is targeted, diminishing energetic returns are expected at increasing distances from the central place, because of increased travel time. Thus, at greater distances animals should feed on higher value prey to optimize foraging efficiency (Schoener, 1969(Schoener, , 1979. It is known that toothed whales (Odontoceti) hunt and capture individual prey across a wide size range (Clarke, 1996;Andersen et al., 2016); however it remains unclear whether odontocetes selectively vary feeding rates when feeding at different depths and/or modulate their locomotor and diving performance in response to changes in prey distribution and availability. To test the above hypothesis, we explore the finescale foraging behavior and kinematics of the cephalopod-eating Risso's dolphin using motion-sensing and acoustic recording tag data to quantify the foraging performance of this species in the wild.

Dolphin Data Collection
Archival tags (Johnson and Tyack, 2003) were attached dorsally to free-ranging Risso's dolphins off Santa Catalina Island, California during field efforts from 2010 through 2016. The tags had a suction-cup system that released at local sunset, or unintentionally due to movement of the dolphins. They were positively buoyant and transmitted a VHF signal when at the surface, allowing tracking, and recovery. Stereo acoustic data were sampled with 16-bit resolution at a sampling rate of 240 kHz, except for tag gg11_216a where 120 kHz was used. Data from the other sensors were sampled at 200 Hz per channel and converted into depth, pitch, roll, and heading of the tagged animal, following methods described in Johnson and Tyack (2003). Pressure and accelerometer data were decimated to 20 and 5 Hz, respectively, before analysis. Analyses were carried out in Matlab 9.1.0 (http://www.mathworks.es/) using the DTAG toolbox (https://www.soundtags.org/dtags/dtag-toolbox/) and custom-made scripts. A subset of 18 tagged dolphins were subject to controlled acoustic playback exposure experiments (CEE) ( Table 1). Detailed exposure protocols are described in Southall et al. (2012). Dives were classified as "non-CEE, " "pre-exposure, " "CEE, " or "post-exposure, " depending on the time overlap with the onset and end of the acoustic stimuli. Generalized additive mixed models (GAMMs) were used to assess the effect of CEE mode and type on two representative metrics of foraging performance. We fitted GAMMs using maximum dive depth and buzz rate as response variables, including CEE categorical data as a factor and individual dolphin as random effect, to identify potential behavioral responses as a function of CEE mode (non-CEE, pre-CEE, during-CEE, and post-CEE) and type (non-CEE, Control, PRN, Simul-MFA, Real MFA). As we were testing a specific hypothesis no model selection was needed. This statistical approach allowed us to assess whether there was a behavioral response as a consequence of CEE, which was important to determine if the full Grampus tag dataset could be pooled across CEE modes and types for further analysis. Foraging models were fitted as a functions of dive metrics and of CEE modes and types, using the mgcv 1.8-15 package (Wood, 2011) in R 3.3.2 "Sincere Pumpkin Patch" (http://www.R-project.org/).
Focal follows of tagged animals were conducted whenever possible, avoiding approaches closer than 25 m. Observations were made from the tag boat with the aid of binoculars and VHF radio tracking equipment. Individuals were identified via photos of their dorsal fin and scar pattern, which facilitated tagging different individuals every time. Location of the tagged animals was recorded with a GPS on board and seabed depth at dolphin sighting locations was extracted using NOAA's NGDC bathymetric charts in ArcMap 9.2 (Environmental Systems Resource Institute, Redlands, CA, USA) with 3arc second resolution (https://www.ngdc.noaa.gov/mgg/coastal/ grddas06/grddas06.htm). Average seabed depth was estimated from all positions gathered during each tag deployment and the maximum depth of the deepest dive recorded among the 33 dolphins. This was 560 m, considering that in the deepwater habitat the seabed must be deeper than the deepest dive recorded. In this way dolphins were broadly classified as foraging in a shallow or deep-water habitat ( Table 1). Only two of the dolphins (data sets gg14_223a and gg14_253a) were close to the boundary defined by this depth threshold, as they remained above average seabed depths of 450 and 500 m, respectively. Those depths may have allowed deep diving but if so, foraging would likely have occurred close to the seabed, in the benthic boundary layer (Angel and Boxshall, 1990). Given the differences between benthopelagic and pelagic organisms, these dolphins were classified as in a shallow-water habitat. Moreover, data for each dive were analyzed with respect to which habitat the dolphin was in, based on the seabed depth of the closest focal follow position gathered within 15 min of the start or end of a given dive. We acknowledge the low accuracy of these estimates but the seabed depth at the nearest focal follow location was considered a representative measure of the type of foraging habitat used by the animal. Dives without focal follow information within the 15 min time window was not considered for kinematic analyses.

Prey Data Collection
Data on the distribution of prey in the Catalina Basin, off the eastern coast of Santa Catalina Island, California, comprising the deep-water habitat of the dolphins, were obtained from ship and underwater autonomous vehicle (AUV) based hydroacoustic surveys. The ship transects covered seabed depths ranging from 300 to 900 m and the AUV sampled pre-defined depths between 50 and 500 m, depending on the location of prey layers. Both platforms were equipped with Simrad EK60 echo sounders at 38 and 120 kHz. The AUV survey provides a 15 × 10 cm sampling resolution in horizontal and vertical planes, respectively. This allowed individual animals to be observed within scattering features, whereas the ship-based sensors provided a view of entire features. Comprehensive explanation of the sensors and platforms used is provided by Moline et al. (2015) while Benoit-Bird et al. (2017) provide a detailed description of biological sampling and active acoustics methodology. Briefly, acoustic scattering data were processed using Echoview and individual prey detected from the AUV within scattering features were identified as single targets (Sawada et al., 1993), providing measurements of target strength at two frequencies. Frequency differencing was used to facilitate coarse taxonomic classification (i.e., fish, squid, or crustacean) while target strength was used as a proxy of length within each taxonomic class. Calculation of inter-individual spacing of prey in layers was measured using the nearest neighbor distance for each individual target sampled by the AUV, both in the beam and along the track. The dominant composition of layers identified from the active acoustics was determined using a net towed at relatively high speed (1-2.2 m s −1 ) that captured mobile organisms between 1 and 35 cm body length. This catch matched the composition identified acoustically (Benoit-Bird et al., 2017). Sampling of the deep-water habitat of the dolphins was coincident in space and time with tagging of two of the 33 dolphins (tag id gg13_266b and gg13_267a, Table 1). The spatial coverage of the prey mapping and its overlap with the tag data is presented in Figure 2 in Arranz et al. (2018). Another 15 dolphins were tagged in the same general habitat, 6 of them in the same year. Prey fields in the shallow-water habitat were not sampled in this study. Dolphins in shallow-water habitat were mostly tagged in 2013, when increased upwelling productivity conditions may have influenced prey availability, particularly of market squid (Doryteuthis opalescens) in coastal waters off the island. Hence the dolphin's foraging opportunities were higher close to shore (Vanderzee et al. pers. comm.). This was supported by visual observations during tagging in 2013, when market squid were observed at the surface and in the mouths of Risso's dolphins. Dives occurring within the shallow-water habitat were classified based on their depth distribution as Shallow (<100 m depth) or Intermediate (>100 m), for comparison of dive parameters.

Kinematic Analyses
Dives were defined as vertical excursions exceeding 20 m depth. Dives recorded within 15 min of tag-on were excluded from the analysis, to remove data potentially affected by the tagging procedure (the 15 min duration is roughly equivalent to two deep dives or 5 shallow dives by this species). Incomplete dives at the start or end of the record (caused by a delay in triggering the salt water switch or release at depth) were also excluded. Dives performed in the deep-water habitat were classified according to their maximum depth and the prey feature targeted in the bottom phase (sensu Arranz et al., 2018). Echolocation clicks and buzzes from the tagged dolphin, respectively, indicative of prey search and capture attempts within dives, were isolated using a supervised click detector together with spectrogram visualizations of acoustic recordings, following the methods described in Arranz et al. (2016). The interpretation that buzzes are associated with prey capture attempts has been confirmed for several echolocating marine mammals, including sperm whales (Miller et al., 2004), beaked whales , pilot whales (Aguilar de Soto et al., 2008), porpoises (Deruiter et al.,  2009), and belugas (Ridgway et al., 2014). Foraging dives were defined as including one or more buzzes. For each dive, the mean descent and ascent depth rates and the proportion of descent time spent gliding were computed. Fluke strokes were identified in the recordings as cyclic variations in the pitch with a magnitude >3 degrees and with a period between 0.4 and 5 s, estimated from a nominal fluking period of 1.5 s for Risso's dolphin. The stroke count was verified visually by checking the pitch angle in random sections of the record. Glides were defined as 10 s time intervals with no stroking. Descent, bottom and ascent phases in dives were defined following a 70% rule, as in Arranz et al. (2016). The descent phase of each dive was considered the period from when the dolphin left the surface to the first time the depth exceeded 70% of the maximum dive depth. The ascent phase started at the last time the depth exceeded 70% of maximum dive depth and ended when the dolphin reached the surface. The bottom phase was defined as the period from the first to the last time the depth exceeded 70% of the maximum dive depth. The inter-buzz interval (IBI) was estimated as the time elapsed since the end of a buzz and the start of the next one. The vertical swim speed was computed as the rate of change in depth (as in Aguilar de Soto et al., 2008).

Energetic Analyses
As a proxy for the activity level of the dolphins we computed the overall dynamic body acceleration (ODBA, Wilson et al., 2006) as the norm of the high-pass filtered acceleration resulting from the movement of the animal and recorded by the tag sensors (ODBA tool from the soundtags toolbox). Cut-on frequency required for ODBA estimation was set to 0.35 Hz (about half of the normal stroking rate for these dolphins). Previous studies have shown that ODBA is a reasonable proxy for metabolic rate in marine mammals (Fahlman et al., 2013) although other studies report difficulty in observing such a relationship (Halsey et al., 2011;Halsey, 2017). An exploratory analysis of this relationship in Risso's dolphins is presented here, assuming that for each individual dolphin, the lowest ODBA recorded on each dive (ODBA min ) corresponded to theVO 2 min , or estimated resting metabolic rate (RMR) and the highest ODBA for each dive (ODBA max ) corresponded to itsVO 2 max or estimated maximum metabolic rate.
For comparative purposes, we modeledVO 2 min in three ways ( Table S1 in Data Sheet S1): (1) using the measuredVO 2 min from the smaller bottlenose dolphin (Tursiops truncatus) (Equation 3.2 in Fahlman et al., 2018) and scaled allometrically for the larger Risso's dolphin, (2) assuming an average mass-specificVO 2 min of 3.9 ml O 2 kg min −1 ( Figure 3B in Fahlman et al., 2018) and multiplying by the estimated body mass (M b ), or (3) using Kleiber's equation (Equation 1) for terrestrial mammals and the estimated M b of the tagged dolphin (Kleiber, 1947). The M b of tagged Risso's dolphins was approximated from individual age class (Table 1), implying a mass of 500 kg for adults and 300 kg for sub-adults and juveniles (Jefferson et al., 2008).
VO 2 max was modeled in two different ways ( Table S1 in Data Sheet S1): scaling upVO 2 min (1) by 5 times or Metabolic Equivalents of Task (MET's, Savage et al., 2007), a level of exercise that one would expect animals to undertake routinely; or (2) by the Hoppeler equation (Equation 2) (Weibel and Hoppeler, 2005).
A regression equation was built using estimatedVO 2 min , aVO 2 max for each dolphin, in one of the ways described above, estimated ODBA min and ODBA max for a given dive. We acknowledge reasonable concerns about generating a regression equation from two data points, but it allowed us to estimate theVO 2 (or field metabolic rate) of the dive based on the overall activity. Regression lines were re-calculated for each dive to account for different tag placements within the same individual (due to potential tag moves) and likely variation of minimum and maximum ODBA values amongst dolphins ( Figure S1 in Data Sheet S1). The resulting intercept and slope of the regression line obtained for a given dive were multiplied either by the average ODBA of the dive+ surface interval, to estimate the averageVO 2 per dive, or by the average ODBA of the dive phase (descent, bottom and ascent) to estimate the averageVO 2 of each phase, as follows: Recorded ODBA min and ODBA max per dive were regressed as a function of duration of the deployment, to rule out correlation between these variables ( Figure S2 in Data Sheet S1).
The ratio between the estimated metabolic costs of each dive phase in different dive types was used comparatively, to assess relative differences in dolphin metabolic costs when diving to different depths.

Statistical Modeling
To assess differences in dive parameters by dive type, we used generalized linear mixed-effects models fitted in R, using the package glmmTMB (Brooks et al., 2017). We modeled each dive parameter as a function of habitat and dive type, controlling for effects of age class (via fixed effects) and individual differences and temporal autocorrelation (via a random intercept of tag ID). We used Gaussian models with identity link, except for the following parameters. Inter-buzz interval and buzz duration were modeled using gamma family with log link, to account for rightskew and overdispersion of residuals, and for number of buzzes a negative binomial type I with log link and an offset for dive duration, since data are count data, and to account for right skew and overdispersion of residuals. We used ANOVA to compare models with and without habitat and dive type as predictors, and in the case of a significant result, we used Tukey's method to assess pairwise differences between dive types [using R package emmeans (Lenth, 2018)].

RESULTS
A total of 124 h of acoustic and movement data from 33 Risso's dolphins (20 adults, 5 sub-adults, 1 juvenile, and 7 individuals of unknown age) resulted in 331 dives during 22 deployments as analyzed here. Most tags recorded data during daytime only but two documented 2 h 19 min and 3 h 30 min after local sunset, respectively, resulting in an overall time coverage from 7:25 a.m. to 9:57 p.m. local (UTC−7) (Table 1, Figure 1). No clear timeof-day dependent pattern of feeding behavior or inshore-offshore movements were detected from plots of the seabed depth at the focal follow positions of the dolphins in relation to time of day. These results are based on a small sample size without full day/night periods from tag records. GAMMs revealed no significant effect of CEEs on the buzz rate per dive nor the maximum dive depth in any of the 18 out of 33 tagged dolphins subject to the experiments ( Table 2).

Habitat
Hydroacoustic sampling of the biotic structure in the water column throughout the Santa Catalina Basin >560 m seabed depth, i.e., deep-water habitat, revealed prey fields vertically stratified into four prey features comprising three soundscattering layers. Each prey layer was horizontally identified as follows: Shallow

Foraging Behavior
Dolphins spent 83 (SD 48)% of the time near (<20 m) or at the surface and 17 (SD 35)% of the time diving. The pattern of dive cycles differed among individuals, with variable-duration bouts of foraging dives interspersed with periods at or near the surface (entailing behaviors like resting, traveling, socializing) (Figure 1). Tagged dolphins foraged throughout the day, with search clicks recorded in 76% of dives and buzzes recorded in 62% of dives. Foraging dives lasted 5.0 [0.7-10] min on average [range], covering a broad range of depths up to 588 m (Figure 1) . Deep dives occurred mostly in the afternoon, since before noon most data were recorded within the shallow habitat. Diurnal patterns of inshore-offshore movement or vice versa were absent when visualizing plots of seabed depth in relation to time of day (average Pearson's correlation across tags R= −0.14). Seabed depths from focal follows revealed that dolphins were foraging near the seabed when in shallowwaters, although benthopelagic feeding may also take place in the deep-water habitat (Figure 2).

Dive Kinematics
Dolphins' descent and ascent rates were up to twice as fast when diving deeper (Spearman's correlation ρ = 0.59, p-value < 0.005, and ρ = 0.65, p-value < 0.005, respectively). They also tended to adopt higher absolute body pitch angles (Spearman's correlation ρ = 0.51, p-value < 0.005) with active stroking on ascent from deeper dives ( Table 3). The proportion of time spent gliding relative to the duration of the descent phase of each dive increased up to 5 times with maximum depth of the dive (Spearman's correlation ρ = 0.53, p-value < 0.005), with an associated dramatic drop in fluke stroking rate ( Figure 4C). Stroke frequency decreased over the first 30 s of the descent during deep dives, from 15 to < 5 strokes every 20 s, after which they glided for about a minute. After that, they increased stroke frequency again, resulting in high vertical velocities (up to 4.5 m s −1 ) during the final phase of most descents of deep dives (Figures 1, 4E). Stroking patterns during descents of shallow, intermediate and midwater dives were more variable, with fewer episodes of high vertical swim speed. Deep dives featured on average twice the number of prey capture attempts in the bottom phase (7 ± 3 buzzes) as shallow dives (3 ± 6 buzzes), whereas in scattered patches and midwater dives dolphins featured 5 ± 4 prey capture attempts, respectively. There was a significant correlation between maximum dive depth and number of buzzes performed at the bottom phase of the dive FIGURE 2 | Foraging activity of Risso's dolphins relative to the seabed. Seabed depth from the closest focal follow position gathered within 6 min of the given buzz (x-axis) and difference between buzz depth and seabed depth (y-axis), colored according to individual. Positive differences are the result of the offset between focal follow position from where seabed depth is extrapolated and depth where dolphins are feeding. Feeding during descent and/or ascent of dives increased the number of prey items targeted in shallow dives by 17% (4 ± 7 buzzes), 40% in scattered patches (6 ± 5 buzzes), 60% in midwater dives (8 ± 3 buzzes), and 28% in deep dives (9 ± 3 buzzes). Overall foraging rates, i.e., the cumulative number of buzzes emitted per dive divided by duration of dives, revealed between 28 and 55 prey attempts per hour, depending on the foraging depth. Dolphins pursued on average 32 ± 52 prey items h −1 in shallow dives, 36 ± 26 items h −1 in intermediate dives, 55 ± 40 items h −1 in midwater dives and 28 ± 20 items h −1 in deep dives.
GLMMs revealed significant differences in buzz rate for different dive types (chisq = 9.83, p-value = 0.02), although pairwise comparisons indicated significant differences only for shallow vs. intermediate dives (Data Sheet S2). Differences in the duration of the bottom-phase of the dive were significant across dive types (chisq = 110.01, p-value = 1.08 −23 ) with differences between intermediate and deep dives, as well as shallow dives and all other dive types. Similarly, mean stroke rate at the bottom phase of the dive differ significantly for all dive types  (chisq = 3, p-value = 3.62 −8 ). According to the full model, IBI did not differ by dive type (chisq = 2.69, p-value = 0.44), however pairwise comparison for individual dive types reported significant differences between shallow and deep dives. There was also evidence of a difference in roll variance by dive type (chisq = 8.56, p-value = 0.03), with shallow dives being different from midwater and deep dives. Buzz duration did not vary among dive types (chisq = 4.28, p-value = 0.23). None of the models indicated significant differences in dive response variables when fitted with habitat type as predictor. Figure 3 shows kernel density of the above modeled parameters separated by age class, habitat and dive type.

DISCUSSION
The Risso's dolphins in this study appear to employ several mechanisms to maximize foraging efficiency, depending on their feeding depth. The dolphins varied their diving kinematics FIGURE 3 | Summary statistics of Risso's dolphins foraging parameters by dive type. Multi-panel violin plot of kernel probability density colored by individual dolphin, different shapes correspond to age class (dots represent adults and triangles non-adults) and separate panels correspond to each habitat. Results of the GLMM (Generalized linear mixed model) fitted for the above parameters with and withouth habitat and dive type as predictors are reported in the text and Data Sheet S2.
and prey quantity when foraging in prey layers with different composition. They foraged throughout the day in shallow and deep-water habitats, attempting to capture prey at depths up to 488 m. Foraging dives had maximum dive depths up to 560 m and occurred in bouts representing about 20% of a dolphin's daytime activity, although tag-on durations averaged only a few hours and at most 13.5 h. These foraging dives targeted epipelagic to mesopelagic prey resources and may have also included near-benthic and coastal water foraging. Foraging bouts were interspersed by traveling/resting/socializing periods of variable length at the surface, which may be needed for recovery. Diving patterns were stereotyped when foraging at particular depths irrespective of the habitat type. They altered their activity for different types of dives in response to metabolic demands, indicating that their foraging tactics were influenced by foraging costs and benefits. When feeding in the shallow-water habitat, buzz depth and seabed depth at focal follow positions point to dolphins targeting benthic organisms during some of these dives. Opportunistic observations report Risso's dolphins feeding on non-gregarious benthic prey, like octopods in very shallow-waters (Ruiz et al., 2011). Stomach contents from stranded specimens from another location reveal that the benthic octopus Eledone cirrosa can represent up to 55% of the cephalopod prey consumed by this species (Clarke and Pascoe, 1985;Blanco et al., 2006), which would support the potential benthic feeding strategy inferred from data presented here. Nevertheless, foraging patterns were comparable when diving at similar depths but in different habitats.
Furthermore, dolphins feeding in deep-water exhibited several adaptations to reduce the energy cost of locomotion when searching for food at different depths, such as the use of gliding gaits and higher pitch angles during ascent and descent (Table 3). Therefore, despite the longer transit distance to forage for deep-water resources, the dolphins managed to keep their average field metabolic rate similar across dive types (i.e.,

FIGURE 4 |
Representative examples of fine-scale kinematics of a Risso's dolphin feeding in a deep-water habitat during deep, midwater, and shallow dives (dolphin with tag id gg13_267a). Dive profiles with buzzes in red (A) and associated measurements of body roll (B), fluke stroke rate (C), ODBA rms (D) and vertical swim speed (E) during the first dive of each type recorded on this tag (15 min after the tag on). Red dots indicate prey capture attempts (buzzes). ODBA rms corresponds to the root-mean-square (rms) of the 2-norm overall dynamic body acceleration (ODBA). Shallow and midwater dives started 5 min before and 12 min after the start of the deep dive, respectively, by the same dolphin and within the same general habitat. Note the gliding periods during the second part of the descent during deep and midwater dives.
shallow, intermediate, midwater, and deep). While our approach to estimate field metabolic rate relies on the assumption that the lowest and highest ODBA for each dive represent resting and maximal metabolic rate. This assumption may not be true for each individual dive, but this analysis allows us to explore the relative metabolic costs within dives, and accounts for possible changes in tag placement on overall ODBA. In addition, the heterogeneity in the regression estimate between ODBA and metabolic rate within an individual between dives was acceptable (C.V. of 32% for slope and of 0.7% for intercept), suggesting that there were not large variation in estimated costs between dives. It has been suggested that glide and stroke behavior helps conserve energy (Williams et al., 2017). Gliding during dive descent may help to reduce the metabolic costs of deep dives up to 40% to balance their need to perform long ascents while fluking actively. Similarly, Steller sea-lions diving to artificial prey patches from 5 to 50 m in experimental studies have been observed to decrease their overall (average) diving metabolic rate and activity level with dive depth (Fahlman et al., 2009(Fahlman et al., , 2013. These adaptations may allow them to extend their time at depth in deeper dives and to make a greater number of prey attempts at the bottom of the dive; although bottom buzz rates remained comparable for shallow vs. deep dives. Foraging during transits from (and to a lesser extent to) foraging depths slightly increases the overall foraging rate during midwater dives compared to intermediate dives, but not in deep dives. This represented about a quarter (25%) of prey caught in deep dives, however if dolphins were still diving deep to forage, it suggests intermediate and midwater prey is of lower quality, with respect to prey energy content and/or prey density per unit volume, compared to deeper prey.
Hydroacoustic sampling within the same general area did not find greater prey density at depth (Benoit-Bird et al., 2017). Active acoustics data revealed that inter-individual distances in schools targeted by Risso's dolphins increased with total body length of prey items (Benoit-Bird et al., 2017). Moreover, the scattering layers were stratified vertically in terms of composition and size. The average length of individuals constituting mono-specific prey patches in the Catalina basin generally increased with depth, while the average density of prey in patches decreased. It therefore seems that the dolphins benefit from deeper foraging trips through accessing prey of varied quality rather than greater prey density. The longer IBI associated with deeper buzzes is consonant with dolphins targeting more widely spaced prey items (i.e., larger distance between individuals) during deep dives or that it requires longer handling times, an interpretation compatible with prey data reported here.
The records of higherVO 2 in the bottom phase of deep dives is consistent with the registered increase in body roll and stroke rate and suggest a greater effort involved in the pursuit of deep-water prey items. These results support the notion of shift in prey type with increasing foraging depth and could be related to the observed escalation of prey length as  (Fernández et al., 2017). Such squid, including species from genera Todarodes, have strong locomotor abilities (O'dor and Webber, 1991) and present high energy density per unit of wet weight (3-4.5 kJ g −1 wet weight , constituting a potential high-quality target for Risso's dolphins during deep dives. Data on swimming speeds of the dolphins while foraging at the bottom could not be estimated reliably in this study (but see Cade et al., 2017), however the relatively high vertical swimming speeds recorded in the bottom phase of deep dives (up to 4.5 m s −1 ) (Figure 4), are coherent with dolphins chasing prey at depth. Deep-water sprints associated with capturing large squid have been described in short-finned pilot whales (Globicephala macrorhynchus) (Aguilar de . Similarly, Aoki et al. (2012) reported bursts of speed in sperm whales (Physeter macrocephalus), suggesting that they performed these bursts to catch powerful and nutritious deep-water prey (i.e., large and/or muscular). In deep dives, we found no strong evidence of longer buzzes that have been associated with the capture of larger prey items in beaked whales , although it may depend on the clicking rate and closing speed at which targets are approached. Nevertheless, the apparently high levels of exercise undertaken by Risso's dolphins during deep-water feeding recorded here are highly significant, together with the measured increase of total body length of prey items in schools targeted by these dolphins as a function of depth. This evidence points to a potential common strategy of some deep-diving toothed whales, in which the capture of larger, more nutritious prey in deep-waters may compensate for the higher metabolic costs of performing longer foraging dives at depth. Some deep diving species, such as sperm whales and pilot whales, adapt to circadian migrations of prey distribution (Watwood et al., 2006;Aoki et al., 2007;Aguilar de Soto et al., 2008) and exploit deep-water resources when they are available at shallower depths at night. Other species, like Blainville's beaked whales (Mesoplodon densirostris) seem to adjust little (Arranz et al., 2011). All prey layers in the deep-water habitat of Risso's dolphins, other than the deepest, have been described to perform diel vertical migrations (Benoit-Bird et al., 2017). It is therefore possible that these dolphins benefit from the availability of deep-water prey in shallower waters at night, so as to reduce the apparently higher energetic costs involved in foraging at depth. Soldevilla et al. (2010) reported higher click detection rates of Risso's dolphins at night within the Catalina basin, suggesting either a higher feeding rate or shallower foraging depths. A few data about night-time foraging behavior of Risso's dolphins presented here are inconclusive, limiting interpretation of whether, and to what extent, this species adapts to circadian changes in prey distribution. Collection of night-time tag data on foraging dolphins would be needed to better understand the feeding patterns of this species.
A number of methods have been developed to estimate metabolic rate in free-ranging marine mammals. In pinnipeds and smaller cetaceans, indirect measurements have been made using doubly labeled water, or species-specific calibrations of proxies of energy expenditure, such as heart rate and activity (Butler et al., 2004;Maresh et al., 2014). In larger cetaceans, heart rate has seldom been measured in situ, nor has the relationship between activity and metabolic cost been validated. However, activity correlates reasonably well with diving metabolic rate in Steller sea-lions (Fahlman et al., 2008(Fahlman et al., , 2013. Recent work has also determined the RMR during restrained near-shore shallow diving or off-shore deep-diving in bottlenose dolphins undergoing capture-release health assessments (Fahlman et al., 2018). Indeed, there is a relative similarity in the deep diving capacity, anatomy, morphology, and body structure of bottlenose dolphins and Risso's dolphins. Despite the fact bottlenose dolphins are within Delphininae and Risso's dolphins within Globicephalinae, groups separated by about 10 million years of evolution, we hypothesize that a comparable relationship exists in Risso's dolphin. While we have made the assumption that the maximal and minimal ODBA represent maximal and restingVO 2 , respectively, we cannot confirm that these are true estimates of these two measures. This is especially true for the resting ODBA, where the variability was higher for tag durations that were < 4 h (Figures S1, S2). However, for maximal ODBA and minimal ODBA for tag durations > 4 h and up to 13.5 h there were no systematic trends with tag duration. We therefore argue that if these estimated minimal or maximal values were grossly over-or underestimated, we would have seen at least one outlier for 13 dolphins with an overall tag duration of 92 h. We also used two independent methods, based on Kleiber's equation (Kleiber, 1947) and Fahlman et al. (2018), to estimate the RMR in free-ranging Risso's dolphins. This was to address the caveats associated with the application of respirometry data from restrained bottlenose dolphins in relatively warm water, since temperature differences have a potential impact. Consequently, the estimated absolute metabolic rate is higher in the larger Risso's dolphin, whereas the mass-specific metabolic rate is lower (Lavigne et al., 1986). Although the results of the three methods varied between 1 and 3 l O 2 min −1 , they all reveal a consistent overall consumption rate across dive types (Figure 5), indicating that Risso's dolphins compensate for energetic costs associated with diving, swimming and feeding at different depths (Davis et al., 2004).
The calculated aerobic dive limit (cADL), the total O 2 stores divided by the divingVO 2 , provides a useful index to indicate whether the animal mainly forages aerobically or possibly uses anaerobic pathways (Butler and Jones, 1997). The majority of diving animals are thought to perform most dives well within their cADL, since this maximizes long-term foraging efficiency (Costa, 2001). Nevertheless, there are some notable examples like fur seals that appear to regularly dive outside their cADL ). In the current study, we used the estimated mass-specific O 2 stores from the bottlenose dolphin [36 ml kg −1 , see Table 5 in Ponganis (2011)] and the average divingVO 2 for each dive type (via Equation 3.2 in Fahlman et al., 2018) to estimate cADL. This resulted in a cADL ranging between 14.8 and 16.2 min and 8.9-9.7 min for adult and non-adult Risso's dolphins, respectively. While our estimated divingVO 2 entail a number of uncertainties, they provide an initial estimate of energy use and indicate that all dives in the current study are apparently within the estimated cADL.
The Risso's dolphins in this study fed in shallow and deepwater habitats, foraging in different depth layers in deep-waters (occasionally within the same dive) and possibly benthically near the shore when in shallow-waters. In order to balance the energetic needs and costs of foraging at increased distances from the surface, these dolphins altered their dive kinematics to minimize energy consumption during longer transits and extended feeding time with increased depth. They thus allowed for a higher number of prey encounters per dive. Moreover, in spite of the increased overall cost of deep dives, they left longer intervals between prey capture attempts during deep dives. Even faced with the increased need to manage oxygen during deeper and longer dives, prey capture appeared to be more energetic during deep foraging. To maintain energetic returns at increasing distances from the surface, Risso's dolphins appear to choose to feed on larger prey at increasing depth, which points to a potential increased prey size-distance relationship for these predators. This strategy is however difficult to measure in free-ranging diving animals. This adaptive strategy appears to increase their foraging efficiency both within and between foraging bouts, thus enhancing their fitness throughout different habitats. Together with this, these predators show an apparent flexibility in using a variety of resources while maintaining their average metabolic rates per dive similar across foraging depths.

DATA AVAILABILITY
Datasets are available on request. The raw data supporting the conclusions of this manuscript will be made available by the authors, without undue reservation, to any qualified researcher.