The predictive relationship between heart rate (fH) and oxygen consumption () has been derived for several species of marine mammals swimming horizontally or diving in tanks to shallow depths. However, it is unclear how dive activity affects the relationship and whether the existing equations apply to animals diving to deeper depths. We investigated these questions by simultaneously measuring the fH and of Steller sea lions (Eumetopias jubatus) under different activity states (surface resting or diving), types of dives (single dives or dive bouts), and depths (10 or 40 m). We examined the relationship over dives only and also over dive cycles (dive + surface interval). We found that fH could only predict over a complete single dive cycle or dive bout cycle (i.e. surface intervals had to be included). The predictive equation derived for sea lions resting on the surface did not differ from that for single dive cycles. However, the equation derived over dive bout cycles (multiple dives + surface intervals) differed from those for single dive cycles or surface resting, with similar fH for multiple dive bout equations yielding higher predicted than that for single dive bout cycles (or resting). The relationships were not significantly affected by dive duration, dive depth, water temperature or cumulative food consumed under the conditions tested. Ultimately, our results demonstrate that fH can be used to predict activity-specific metabolic rates of diving Steller sea lions, but only over complete dive cycles that include a post-dive surface recovery period.
While there are numerous methods for measuring the energy expenditure of marine mammals in the field, not all are equal in reliability, precision or ease of use. The three most common current methods employ measures of doubly labelled water (DLW) turnover, body acceleration metrics and heart rate (fH). The DLW method provides a mean estimate of field metabolic rate over a finite period that is not activity specific, with individual error estimates in pinnipeds ranging from –39 to +44% (Sparling et al., 2008). Body acceleration metrics, such as overall dynamic body acceleration (ODBA) and flipper stroking, are limited to predicting metabolic rate during active behaviours, and cannot account for changes in physiological state such as digestion (Bevan et al., 1997; Hays et al., 2004; Williams et al., 2004; Fahlman et al., 2008b; Green et al., 2009; Halsey et al., 2009). The fH method estimates oxygen consumption rate (, an indirect measure of energy expenditure) from measured fH (Fick, 1870). fH can therefore provide activity-specific estimates of energy expenditure on a much finer time scale and for longer periods of time than the DLW method (Boyd et al., 2004; Butler et al., 2004; Ponganis, 2007), with comparable error estimates to ODBA under steady-state conditions (Green et al., 2009; Halsey et al., 2009).
The fH method requires deriving species-specific predictive equations between fH and in a controlled environment before the technique can be used to predict in the field (for a review, see Green, 2010). Previous studies in marine mammals and birds have made comparisons while animals were submerged in a shallow swim mill (Woakes and Butler, 1983; Williams et al., 1991; Butler et al., 1992; Boyd et al., 1995; Ponganis et al., 1997; McPhee et al., 2003), walking on a treadmill (Froget et al., 2001; Green, 2001; Froget et al., 2002; Green et al., 2005), swimming horizontally in open water (Williams et al., 1993) or diving in a shallow tank (Bevan et al., 1992; Webb et al., 1998a). However, the artificial modes of locomotion and environments employed in these studies raise questions of applicability to animals that spend a considerable time diving to depth.
Despite the current use of fH to predict of diving marine mammals in the wild (i.e. Hindell and Lea, 1998; Boyd et al., 1999), it is not clear whether the fH method is accurate for pinnipeds foraging at natural depths and for realistic dive durations. Accurate estimates of metabolic rate in the wild require relationships that are derived under controlled conditions which encompass dive durations and dive depths that are representative of free-ranging animals. Previous studies that have investigated fH and in diving marine mammals have been limited by maximum tank depth (Webb et al., 1998a; Sparling and Fedak, 2004). Our study derived relationships between fH and in trained Steller sea lions freely diving in the open ocean to depths up to 40 m and for durations of 1–6 min, which reflect dive characteristics comparable to those of free-ranging animals (Merrick and Loughlin, 1997; Rehberg et al., 2009). Our specific objectives were therefore to (1) simultaneously measure and determine the relationship between fH and while Steller sea lions were foraging and diving to depths of up to 40 m and (2) determine whether fH could be used to predict average metabolic rate (AMR) or diving metabolic rate (DMR) over either a single dive or a series of continuous dives (dive bout).
MATERIALS AND METHODS
Our study was conducted from August to September, 2008, using three trained female Steller sea lions (Eumetopias jubatus, Schreber 1776) housed at the UBC Open Water Research Laboratory (Port Moody, British Columbia, Canada). Animals were housed in a specially designed floating pen that provided access to seawater and haulout space [for a full description, see Hastie et al. (Hastie et al., 2007)]. Two of the sea lions were 11 years old (body mass of F97SI, 218±4.4 kg; F97HA, 172±0.6 kg) and one was 8 years old (F00BO, 145±5.1 kg; total mass range 130–224 kg). All sea lions had been raised in captivity and had previously been trained to use all experimental apparatus. On non-trial days, animals were fed a diet of herring (Clupea pallassi) supplemented with vitamin tablets. Sea lions were fasted overnight, and then weighed on a platform scale prior to trials (±0.5 kg). All procedures and protocols were conducted in accordance with the guidelines of the Canadian Council on Animal Care (University of British Columbia no. A07-0413), and under permits from the Department of Fisheries and Oceans Canada (no. MML 2007-0001) and the Vancouver Aquarium. All animal work was conducted voluntarily under trainer control.
Measurement of fH
fH was measured using subcutaneous electrodes connected to a wireless transmitter, coupled to a datalogger. fH monitoring equipment composition and placement are detailed elsewhere (Young, 2010; Young et al., 2010). Briefly, animals were outfitted with fH monitoring equipment while anaesthetized (0–5% isoflurane in O2), and all equipment was removed after each trial while the animals were under trainer control. The fH monitoring system consisted of (1) a fH datalogger (HTR, Wildlife Computers, Redmond, WA, USA) which recorded the R–R interval and (2) a fH transmitter (HRX, Wildlife Computers) that had two ∼80 cm leads. To reduce infection risk, we spliced 30 gauge, 99.9% pure silver Teflon-coated wire (Grass Technologies, Longueuil, QC, Canada) to the terminal ends of the electrode leads, and these ends were sterilized and inserted subcutaneously prior to each trial. The fH transmitter and fH datalogger were carried by the sea lions on a custom-fitted harness. The harness also held a time–depth recorder (TDR) that measured dive time, surface time and dive depth (sampling frequency 1 Hz, SU-05272, ReefNet Inc., Mississauga, ON, Canada). Animals were allowed to fully recover from anaesthesia in a dry area before the start of the trial. fH was continuously measured immediately after both fH electrodes were inserted.
We measured using open-circuit gas respirometry as previously described for Steller sea lions (Hastie et al., 2007). Fractional oxygen and carbon dioxide concentrations within a desiccated (via CaSO4) subsample of the excurrent airstream were measured using Sable System FC-1B and CA-1B analysers, coupled to a 500H mass flow generator and controller (Sable Systems Inc., Las Vegas, NV, USA).
Trials were conducted either at a dive site next to the holding pen or after the animals were transported to a nearby dive site in a 7 m research boat. A second boat towing a floating barge carried the respirometry equipment to the dive site. The floating barge had an opening in the middle containing a cage (152×152×250 cm) used to contain the animal at the surface before and after dives, and a floating transparent Plexiglas® respirometry dome. Air was drawn through the respirometry dome at a constant rate of 475 l min–1. The excurrent airstream was continuously subsampled, and averaged every 0.5 s (Sable Data Acquisition System, Sable Systems Inc.). Fractional gas concentration readings were corrected for electronic drift against ambient air before and after each trial. Barometric pressure and excurrent air temperature were also recorded (Airguide Instruments, Chicago, IL, USA) to correct readings to standard temperature and pressure, dry (STPD). Water temperature ranged between 14.4 and 17.5°C (mean 15.8±1.0°C, N=15 trials).
The sea lions were trained to remain stationary within the floating respirometry dome, and then dive between two feeding tubes positioned at the same set depth (10 or 40 m) ∼6 m apart (mean animal length ∼2 m). The feeding tubes were used to maintain swimming activity at depth. These trial depths reflect the lower and upper range of dive depths commonly noted for free-ranging females (Merrick and Loughlin, 1997; Rehberg et al., 2009). Each complete dive trial was repeated three times at each depth in random order for each animal (N=6 trials per animal), with only one trial per animal per day. Sea lions were not tested on sequential days. In the end it was necessary to discard one entire trial from each animal because of noise in the dataset (N=5 trials analysed per animal, N=15 bouts, N=25 single dives).
A trial consisted of (1) a 6–10 min resting period where sea lions were stationed in the respirometry dome that provided the baseline measurements of surface metabolic rate (MRs), (2) a single dive with recovery period (see below), (3) a 4-series dive bout with recovery period at the end of the bout, and (4) a second single dive with recovery period, of longer dive duration than the first dive (40 m only, Fig. 1). Animals remained at the surface in the respirometry dome during post-dive recovery surface intervals (recovery periods), that ended when measured returned to within 2% of MRs (Fig. 1).
Animals were fed small pieces of herring (20 g) during trials to facilitate cooperation. A tube system allowed fish delivery to set depths and into the respiratory dome without altering the integrity of the flow-through measurements. Sea lions were encouraged (via altered underwater fish delivery rates) to execute dives of 0.5–2 min duration; field data demonstrate that 45% of free-ranging dives are 0–1 min and 33% are 1–2 min for female Steller sea lions (Merrick and Loughlin, 1997). For the 40 m trials only, an additional, longer single dive was encouraged at the end of the trial by providing food underwater at higher rates, which enabled consideration of the influence of dive duration independent of depth on the relationship.
The cumulative amount of food consumed from the time the animals were outfitted with the experimental gear until the end of each specific event (i.e. end of dive or end of dive + surface interval) ranged from 0.4 to 2.8 kg. The amount of food sea lions were fed during a dive trial was intentionally minimized as previous work has indicated that the consumption of a 4 or 6 kg meal changes the relationship in Steller sea lions resting in water (Young et al., 2010). However, the previous study did not account for the potential effects of diving or feeding on smaller amounts of prey at depth on the relationship. Therefore, in this study we explicitly tested the effect of cumulative food consumption on the relationship.
Calculating metabolic rate
was calculated from fractional O2 and CO2 concentrations [see eqn 3b in Withers (Withers, 1977)] using Datacan Data Analysis software (V 1.0.24; Sable Systems Inc.) and then exported into Microsoft Excel. For resting MRs, was measured for 6–10 min after had reached a steady state and averaged over the 2 min preceding the first dive of a trial. We analysed underwater O2 consumption in two ways to maximize the comparability of our data to previously published work. First, we calculated AMR across a dive cycle (dive + surface interval). Second, we calculated DMR over solely the underwater time. Both analyses were performed on both single dives and dive bouts. AMR over a single dive cycle or dive bout cycle was calculated by dividing the total integrated volume (l) of O2 consumed during the dive and surface interval by the entire dive cycle or dive bout cycle duration. DMR was calculated as the integrated post-dive above the MRs baseline, divided by time underwater. For single dives, time underwater was the actual dive duration. For dive bouts, time underwater was the cumulative dive duration.
Unfortunately, there is not a consensus on whether should be scaled with body mass and, if so, which exponent to use (Packard and Boardman, 1999; Brown and West, 2005; White and Seymour, 2005; Savage et al., 2007). Overall, model results with non-mass-corrected (l O2 min–1) and mass-corrected metabolic (ml O2 min–1 kg–0.66 or ml O2 min–1 kg–0.75) did not differ from each other (i.e. different equations, but the same factors and models were significant). In addition, mass could not be incorporated into the equations as a covariate because of either convergence errors or a lack of significance as a factor.
Ultimately, we chose the most parsimonious approach and decided to use units that did not assume a scaling exponent (l O2 min–1). Therefore, all equations presented should only be applied to wild Steller sea lions within our experimental mass range (130–224 kg). To facilitate comparison with previous studies that did scale metabolic rate with body mass [including Young et al. (Young et al., 2010)], all model equations and figures using (ml O2 min–1 kg–0.75) are available elsewhere [see pp. 50–71 of Young (Young, 2010)].
fH measurements stored in the fH datalogger were downloaded and analysed with Microsoft Excel and R 2.9.2 (R Core Development Team, 2009). First, R–R intervals were converted to instantaneous fH (fH,inst in beats min–1). A series of algorithms were applied to systematically remove any fH,inst data that were artifacts of muscle or wire movement [detailed in Young et al. (Young et al., 2010)]. Finally, fH,inst values were averaged across appropriate intervals to match the integrations of the same time periods, yielding mean fH (beats min–1; Fig. 1). Specifically, resting fH was calculated over the 2 min directly preceding the first dive. Mean fH was calculated over the entire dive and surface period, from the start of the dive (or first dive in a bout) until the end of the surface recovery interval. Diving fH was calculated over either the underwater portion of a single dive or the cumulative underwater duration for dive bouts (Fig. 1).
It has previously been suggested that otariid fH may not scale linearly with body mass (Castellini and Zenteno-Savin, 1997). Furthermore, preliminary analysis of our data showed that there were no significant linear relationships between body mass and mean fH over a single dive cycle (ANOVA F-test, F1,21=0.44, P=0.52 for slope) or over a dive bout cycle (ANOVA F-test, F1,11=0.90, P=0.36 for slope). These results, combined with our goal of consistency between fH and units, contributed to our decision not to scale fH to body mass.
Data from each animal within a trial and data from each animal across trial types were treated as a repeated measures dataset using linear mixed-effects (LME) models in R 2.12.1 [nlme library from Pinheiro and Bates (Pinheiro and Bates, 2000)]. LME models consider the correlation between repeated measurements within and among animals, while also characterizing individual animal variation relative to the mean of the population (Pinheiro and Bates, 2000). As animals were allowed to recover to pre-dive levels between single dives and dive bouts (i.e. returned to within 2% of MRs), single dives and dive bouts within the same trial were considered to be statistically independent. All models were run using the maximum likelihood method, and the slope and intercept were allowed to vary for each animal during model optimization. Animal ID was treated as a random effect for all models (which permitted inference from the sample captive population to the free-ranging population) and fixed effects explored included: dive duration (1.0–1.99 min, 2.0–3.99 min, >4.0 min), dive depth (10 m, 40 m), activity (resting at the surface or diving) and cumulative amount of food consumed (0–1.0 kg or >1.0 kg, max. 2.8 kg).
For each analysis, the best model in terms of fixed effect factors was determined using an ANOVA. Execution on a single model generated a conditional F-test to determine the significance of model slope, intercept and fixed effects. As all intercepts proved to be highly significant (P≤0.001), the F-values reported here are for model slopes only. An ANOVA performed on two nested models (the fixed effect model hierarchically nested within the model without fixed effects) produced likelihood ratio tests (LRT) that compared the two models (i.e. whether addition of fixed effect significantly improved null model). All models presented only had one fixed effect applied at a time; therefore, for all LRT tests d.f.=1. Within LME analysis, both models must be linear (i.e. significant slope and intercept) as determined with an F-test before being compared with an LRT test. Therefore, unless noted otherwise, all models compared via LRT tests were linear and statistics reported are for LRT tests. All reported data are presented as means ± s.e.m., and statistical significance was set at α=0.05.
Confidence intervals (95% CI) for final predictive equations were calculated using standard bootstrapping methods with 1000 replicates per model (Venables et al., 1999; R Core Development Team, 2009; Whitlock and Schluter, 2009). The ordered 24th and 976th bootstrapped replicates were then plotted to represent the 95% CI for a model.
The error associated with the predictive equations is not constant and increases with the distance from the mean fH value, but many studies evaluating other techniques to estimate energy expenditure report a single, average error (e.g. Boyd et al., 1995). To facilitate comparison among studies, we estimated model error using both standard error of the estimate (s.e.e.) as well as ‘average percent error’. The s.e.e. for was calculated using the minimum fH for each model (Zar, 1999) [see eqn 10 in Green (Green, 2001)]. It is important to note that the s.e.e. values are maximum estimates because the error at the tails of each model is greater than the error at the mean fH value (see Appendix Table A1 for s.e.e. mean and maximum fH values). All s.e.e. equations included the error associated with variation between sea lions in slope and intercept during model optimization. We also derived the average percent error associated with each model, which was a single representative percent of the error based on the standard residuals. We calculated the average percent error by dividing the mean absolute fixed residual by the median predicted value. We then multiplied this by 100 to calculate a percent error. Standard error of the mean (s.e.m.) for the slope and intercept is also presented in Table 2.
Dive duration for all single dives (N=25 dives) recorded from the TDR ranged from 1.0 to 6.3 min, and cumulative dive time for bouts ranged from 2.6 to 8.0 min (N=13 bouts, Table 1). Maximum dive depth ranged from 11 to 14 m (10 m trials) or from 41 to 58 m (40 m trials) as animals occasionally swam below the feeding stations temporarily before feeding. Our dive characteristics were comparable to those observed in free-ranging female Steller sea lions (Merrick and Loughlin, 1997; Rehberg et al., 2009).
Relationship between fH and resting metabolic rate
There was a significant relationship between fH and for animals resting at the surface in open water (MRs). This relationship does not significantly differ from the relationship derived for four other Steller sea lions resting in a swim mill (∼16 h postprandial) (Young et al., 2010). Further, combining data from these two sets of trials creates a stronger composite predictive equation. This composite equation that described all animals resting in water (fasted in swim mill + MRs in open water) was designated resting metabolic rate (RMR), and was used to make all comparisons to the diving relationships. Mean RMR ranged from 0.81 to 2.30 l O2 min–1 and the corresponding fH ranged from 57 to 108 beats min–1.
Relationship between fH and DMR
Overall, fH and DMR data were normally distributed, had normally distributed errors, and had homogeneous variance across single dives and dive bouts. Diving fH during a single dive ranged from 27 to 68 beats min–1, and DMR ranged from 1.21 to 3.21 l O2 min–1. Over dive bouts, diving fH ranged from 34 to 75 beats min–1 (subsurface time only), and DMR ranged from 1.40 to 3.00 l O2 min–1. Bradycardia was notable for both dive types (on average, dive fH was reduced from resting by 64% in single dives and 44% in bouts). However, no significant differences in DMR (2.06±0.56 l O2 min–1 for single dives, 2.22±0.47 l O2 min–1 for bouts) relative to RMR (1.41±0.35 l O2 min–1) were noted.
None of the predictive relationships for DMR were significantly linear under any analysis conditions (either without fixed effects or including appropriate fixed effects such as depth, dive duration, cumulative food consumption and dive type, or after log transformation, Fig. 2). For example, diving relationships were not linear for dive bout data alone (F1,11=0.89, P=0.365), when combined with data from animals resting in water (F1,77=0, P=1.00) or when combined with single dive data (F1,32=1.18, P=0.283). Given that no predictive equations for DMR were linear, further comparisons (such as DMR vs AMR or RMR) were not undertaken.
Relationship between fH and AMR
There was a significant relationship between fH and AMR, but this was not significantly affected by dive duration, dive trial depth, maximum dive depth or water temperature across single dive or dive bout cycles (dive + surface recovery interval). Cumulative meal size did not change the fH:AMRsingle relationship over dive cycles (LRT=1.14, P=0.89). Animals consumed less food during the first single dive cycle (mean 0.36±0.19 kg, 0–0.62 kg, N=17) compared with the last single dive cycle (mean 1.38±0.38 kg; 0.88–2.08 kg, N=8). Despite this, adding dive cycle order as a fixed factor (i.e. first or last dive bout) did not significantly change the fH:AMRsingle relationship (LRT=6.16, P=0.19). We therefore combined data from all single dive cycles for both trial depths for analysis.
Over a single dive cycle, average fH ranged from 46 to 103 beats min–1 (mean 75±18.1 beats min–1), and AMRsingle ranged from 1.23 to 2.44 l O2 min–1 (mean 1.82±0.35 l O2 min–1 water surface (Table 2: Eqn 1) were statistically similar (LRT=10.73, P=0.030). Therefore, a single equation could be used to describe both data sets (Table 2: Eqn 3, Fig. 3A).
Over a dive bout cycle, average fH ranged from 46 to 97 beats min–1 (mean 71±18.3 beats min–1) and AMRbout ranged from 1.59 to 2.70 l O2 min–1 (mean 2.03±0.33 l O2 min–1). During dive bout cycles, sea lions were fed a maximum of 1.8 kg herring to facilitate cooperation (cumulative mean 1.0±0.33 kg). However, feeding did not affect the significantly linear relationship over a dive bout cycle (LRT=7.93, P=0.09; Table 2: Eqn 4). Combining data sets to test different factors yields slightly different predictive equations with mixed-effects models because of variability in degrees of freedom and replicate datasets from the same individuals across more than one trial type (i.e. Eqns 4, 5a and 6a in Table 2 for dive bouts are similar but not exactly the same). The relationship over a dive bout cycle differed from that for both RMR (LRT=34.66, P<0.0001; Fig. 3B; Table 2: Eqns 5a and 5b) and a single dive cycle (LRT=23.75, P=0.0001, Fig. 3C; Table 2: Eqns 6a and 6b). One aspect of these differences was that predicted AMR was often greater for bout cycles than for either single dive cycles or RMR at a given fH (Fig. 3). Model convergence errors prevented direct comparison of the relationship between dive bout cycles (Table 2: Eqn 4) and the equation for pooled RMR and AMRsingle data (Table 2: Eqn 3, Fig. 3A). However, given that dive bout equations were distinct from both single dive and RMR equations independently (Fig. 3A,B), we suggest that logically dive bouts are also different from the composite RMR + AMRsingle line (Table 2: Eqn 3, Fig. 3A).
The predictive relationship between fH and has been explored in several marine mammal species submerged in a shallow swim mill (Woakes and Butler, 1983; Williams et al., 1991; Butler et al., 1992; Boyd et al., 1995; Ponganis et al., 1997; McPhee et al., 2003), swimming horizontally in open water (Williams et al., 1993) or diving in a shallow tank (Webb et al., 1998b; Sparling and Fedak, 2004). The implicit or underlying assumption in applying these equations to animals in the wild is that these relationships remain constant independent of the animals' specific aquatic behaviour. This is an ambitious assumption given the range of cardiovascular and metabolic adjustments known to be employed by diving pinnipeds.
Our results demonstrate that fH can be used reliably to predict in diving Steller sea lions, but only when averaged over single dive cycles or dive bout cycles where animals are allowed to recover fully from any potential O2 debt incurred during diving. Further, while the predictive equation for single dive cycles did not differ from that for animals resting in water, the relationship for dive bout cycles was significantly different from that for both single dives and resting. In contrast to previous studies that have simply pooled the data for several animals together, our conclusions are strengthened by using mixed-effects models within a repeated measures framework, with each animal treated as a random sample from the population. The statistical strength of our predictive models – despite our sample size – strongly indicates that detectable differences in the relationships represent a true biological change rather than a statistical artifact.
The results of our study also imply that predictive relationships must be derived from controlled studies that incorporate representative dive characteristics of free-ranging animals if fH is to accurately estimate the metabolic rates of diving endotherms in the field. While many previous studies investigating the relationship in diving marine mammals have been limited by tank depth (Webb et al., 1998a; Sparling and Fedak, 2004), our open ocean site allowed sea lions to dive freely with dive characteristics (depths and dive durations) similar to those of sea lions in the wild (Merrick and Loughlin, 1997). The duration of such dives is probably shorter than the estimated aerobic dive limit of adult female Steller sea lions [7.5 min (Richmond et al., 2006); mean 7.1–7.6 min (Hastie et al., 2007)]. It is also important to collect calibration data from foraging behaviour approximating field conditions, given that locomotion can potentially increase the total energetic costs of diving (Williams et al., 2004; Hastie et al., 2007). Our use of two separated feeding tubes created active feeding and underwater exercise, which appeared to elevate DMR (cf. Hastie et al., 2007) to levels that were comparable to pre-dive values.
Relationship between fH and DMR
We determined that fH could be used to predict AMR over a dive cycle, for either single dives or dive bouts, but could not predict DMR. The lack of predictive power for DMR could be a function of the temporal disconnect between measurements of (measured as post-dive excess O2 consumption above resting levels) and fH (measured directly during submergence) (Fig. 1). This comparison assumes that the post-dive increase in is independent of any post-dive tachycardia that can facilitate increased gas exchange upon surfacing (see below). By contrast, the fH:AMR relationship was formulated from data with a greater temporal overlap (Fig. 1). Specifically, the fH component encompassed the dive itself and any surface recovery periods, and thus incorporated both diving bradycardia and post-dive tachycardia.
The lack of a linear relationship between fH and DMR supports the concept that the dive cycle is a discrete and fundamental unit for diving that cannot be partitioned physiologically into underwater and surface elements (Fedak et al., 1988; Butler and Jones, 1997). Several similar studies on diving vertebrates concluded that fH and should be measured over complete dive cycles because their relationships over dives alone were inconsistent (Fedak, 1986; Fedak et al., 1988; Bevan et al., 1992; Butler, 1993). Dives and subsequent surface intervals are physiologically linked because gas exchange and removal of metabolic byproducts (i.e. CO2) only occur at the surface (Kramer, 1988). Furthermore, repetitive dives in many species are associated with a progressive O2 debt, which is not repaid until the bout is completed (Kooyman et al., 1973; Ponganis et al., 1993; Fahlman et al., 2008a). It is therefore not surprising that mean diving fH was not useful in predicting DMR for Steller sea lions.
Relationship between fH and AMR
Our study found significant relationships from which AMR can be predicted using fH over periods of rest, single dive or multiple dive bouts. The relationship was similar when the animals were resting at the surface or performing single dive bouts. The similar relationship between single dive cycles and surface resting suggests that the O2 debt incurred over a single dive cycle was not sufficient to alter the relationship.
However, a unique equation was needed to describe the relationship during multiple dive bouts (Fig. 3). While the range and pattern of fH and across single dive cycles and dive bout cycles overlapped, dive bouts yielded slightly higher for a given fH because of a greater intercept value for the dive bouts (Fig. 3C). More detailed analyses of the dive bradycardia data collected here further demonstrated that mean diving fH was similar across single dives and dive bouts (Hindle et al., 2010), as has been demonstrated in California sea lions (Zalophus californianus, dive durations 1–3 min) (Ponganis et al., 1997).
Although single dives and dive bout cycles are subject to similar physiological and environmental effects, dive bouts incur the cumulative effects of multiple dives, most notably a progressive O2 debt (Kooyman et al., 1973; Ponganis et al., 1993; Fahlman et al., 2008a). Unfortunately, the relationship between fH and is not constant across levels of oxygen debt. This is partly because rates of reoxygenation are affected by the PO2 gradient (Kramer, 1988), such that oxygen is taken up faster into the bloodstream with greater levels of depletion. Therefore, a greater oxygen debt (as from dive bouts) will result in a greater effective upon resurfacing without similar changes in fH.
Whereas single dive cycles probably incorporate minimal exercise (and limited O2 debt), dive bouts often encompass longer cumulative dive time, and more frequent changes in depth and pressure relative to single dives. Although dive depth and dive duration were not significant factors (when averaged over a complete dive cycle or dive bout cycle), it is possible that instantaneous changes in these parameters during the dive bout influenced the overall relationship. Ultimately, the oxygen debt appears sufficiently pronounced in multiple dive bouts to result in a relationship that is distinguishable from both AMR during surface resting and metabolic rate over a single dive cycle.
Errors associated with using the fH technique
The model error for our ‘best’ models over single dive cycles and periods of rest was 17% average percent error (12% s.e.e., Table 2, Appendix Table A1), and the model error over all of the dive bout cycles was 9% average percent error (4% s.e.e., Table 2, Appendix Table A1). The model error for our recommended predictive equations was greater than the error recorded for equations predicting from fH of submerged swimming California sea lions [mean ∼3%, range –28 to +23% (Boyd et al., 1995)], and slightly greater than the estimates derived from ODBA on Steller sea lions [range –4 to +7% (Fahlman et al., 2008b)]. Both of our types of error estimates were also greater than those of the DLW method for Antarctic fur seals [mean ∼3%, range –11 to +20% (Costa and Trillmich, 1988)], but considerably less than error estimates for California sea lions using the same technique [mean ∼36%, range –10 to +86% (Boyd et al., 1995)]. It is relevant to note that caution should be used when comparing error estimates from different studies as calculation methods vary widely among studies.
Our percent error estimates assume that the appropriate predictive equation will be applied to fH data collected in the field. Dataloggers can be used to measure fH, dive depth and dive duration and to determine when an animal is resting at the surface of the water, executing a single dive, or executing multiple dives in series (Ponganis, 2007). These data are essential for applying the most appropriate equation from Table 2 to estimate metabolic rate in the field.
Given the absence of statistical differences between predictive equations for animals resting in water or completing a single dive cycle, we recommend using the composite equation encompassing all trials for animals resting at the surface or when performing single dives (Table 2: Eqn 3). Given that statistical differences exist between the resting and dive bout equations, it follows that the impact of using the ‘wrong’ equation to predict AMR of dive bouts should be considered. We therefore calculated the percent error for using the ‘wrong’ equation by comparing the predicted AMRs for a mean resting fH of 100 beats min–1 among appropriate equations in Table 2. Incorrectly using the predictive equation for animals resting in water (Table 2: Eqn 1) to estimate of a free-ranging Steller sea lion performing multiple dives in a series (Table 2: Eqn 4) would overestimate by ∼25% (i.e. 0.90 vs 1.20 l O2 min–1).
Given the large potential error associated with using an inappropriate equation, and knowing that it is reasonably easy to determine dive behaviour via small TDRs in situations were fH is already being telemetered, Eqn 3 in Table 2 should be used when animals are resting at the surface or executing single dives, and Eqn 4 should be used when animals are executing a series of multiple dives (Table 2). The errors incurred by applying the wrong equation will only alter calculation of energy expenditure during the time when that specific activity is occurring. Therefore, the cumulative error for estimates of daily energy expenditure will depend directly on the activity budget. That is, if energy expenditure is miscalculated for an infrequent activity, the resulting error on daily expenditure calculations will be small. For example, of the 48% of the time female Steller sea lions spend foraging at sea, only ∼22% is spent diving (Rehberg et al., 2009).
Logistically, it is not always possible to distinguish single recovery dive cycles from dive bout cycles in free-ranging animals (i.e. Hindell and Lea, 1998; Boyd et al., 1999). Wild female Steller sea lions dive quite frequently and dive behaviour may differ between the declining western Alaska sea lion population (mean of 17.5 dives h–1 for approximately 5 h day–1 (Merrick and Loughlin, 1997) or the slowly increasing eastern Alaska population (mean of 13.5 dives h–1 for approximately 3 h day–1 (Rehberg et al., 2009). Consequently, dive bout cycles are more likely to accompany successful foraging compared with single dive cycles; single isolated dives with a full surface recovery of O2 stores (AMRsingle) are less likely to occur than dive bouts when sea lions are exploiting a prey patch. We therefore recommend using the predictive equation for AMR over a dive bout cycle (AMRbout, Table 2: Eqn 4) for any ambiguous diving behaviour when using fH to predict field metabolic rate of diving Steller sea lions.
Prudent application of the appropriate predictive equations will facilitate obtaining reasonable estimates of activity-specific energy expenditure in wild Steller sea lions. This information can then be used to quantify the energetic cost of observed or predicted changes in sea lion behaviour related to changes in their biotic or abiotic environment.
LIST OF SYMBOLS AND ABBREVIATIONS
average metabolic rate () averaged over the dive cycle (dive + TMRs)
doubly labelled water
diving metabolic rate () averaged over dive
heart rate (beats min–1)
instantaneous heart rate (beats min–1)
resting at the surface in open water (pre-dive)
overall dynamic body acceleration
resting metabolic rate (MRs + resting swim mill trials) (Young et al., 2010)
post-dive surface interval (min)
SI until returned to ±2% of pre-dive level (MRs)
oxygen consumption rate (l O2 min–1)
Financial support was provided by a grant from the North Pacific Marine Science Foundation to the North Pacific Universities Marine Mammal Research Consortium (NA05NMF4391068), with additional financial support from the US National Oceanic and Atmospheric Administration.
We thank the technicians and training staff at the UBC Open Water Research Laboratory and the Vancouver Aquarium for assisting with data collection and training the sea lions. We are also grateful for the assistance of Ruth Joy and Tiphaine Jeanniard du Dot with statistical analyses.