Bulk-filter feeding is an energetically efficient strategy for resource acquisition and assimilation, and facilitates the maintenance of extreme body size as exemplified by baleen whales (Mysticeti) and multiple lineages of bony and cartilaginous fishes. Among mysticetes, rorqual whales (Balaenopteridae) exhibit an intermittent ram filter feeding mode, lunge feeding, which requires the abandonment of body-streamlining in favor of a high-drag, mouth-open configuration aimed at engulfing a very large amount of prey-laden water. Particularly while lunge feeding on krill (the most widespread prey preference among rorquals), the effort required during engulfment involve short bouts of high-intensity muscle activity that demand high metabolic output. We used computational modeling together with morphological and kinematic data on humpback (Megaptera noveaangliae), fin (Balaenoptera physalus), blue (Balaenoptera musculus) and minke (Balaenoptera acutorostrata) whales to estimate engulfment power output in comparison with standard metrics of metabolic rate. The simulations reveal that engulfment metabolism increases across the full body size of the larger rorqual species to nearly 50 times the basal metabolic rate of terrestrial mammals of the same body mass. Moreover, they suggest that the metabolism of the largest body sizes runs with significant oxygen deficits during mouth opening, namely, 20% over maximum at the size of the largest blue whales, thus requiring significant contributions from anaerobic catabolism during a lunge and significant recovery after a lunge. Our analyses show that engulfment metabolism is also significantly lower for smaller adults, typically one-tenth to one-half . These results not only point to a physiological limit on maximum body size in this lineage, but also have major implications for the ontogeny of extant rorquals as well as the evolutionary pathways used by ancestral toothed whales to transition from hunting individual prey items to filter feeding on prey aggregations.
Citation: Potvin J, Goldbogen JA, Shadwick RE (2012) Metabolic Expenditures of Lunge Feeding Rorquals Across Scale: Implications for the Evolution of Filter Feeding and the Limits to Maximum Body Size. PLoS ONE 7(9): e44854. https://doi.org/10.1371/journal.pone.0044854
Editor: Andreas Fahlman, Texas A&M University-Corpus Christi, United States of America
Received: June 15, 2012; Accepted: August 8, 2012; Published: September 14, 2012
Copyright: © Potvin et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: J. Goldbogen was supported by funding from the United States Navy and Office of Naval Research awarded to J. Calambokidis and Cascadia Research Collective. R. Shadwick was funded by the Natural Sciences and Engineering Research Council of Canada (NSERC). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
The extreme body size of baleen whales (Mysticeti) is generally attributed to the overall energetic efficiency of bulk filter feeding . Instead of hunting single prey items which is typified by the raptorial and suction feeding strategies in toothed whales (Odontoceti), baleen whales are obligate suspension filter feeders that engulf large quantities of prey-laden water. By processing vast amounts of small prey in bulk, baleen whales are thought to reap the rewards of an abundant resource using a more efficient feeding mechanism , . However, the energetics of foraging in cetaceans has yet to be fully assessed, especially with respect to microphagy in large mysticetes. Considering that some baleen whale species represent the largest animals to have ever lived, exploring the metabolic expenditures of feeding in baleen whales may shed light on the evolution of gigantism as well as the potential energetic limits to body size.
Baleen whales exhibit a variety of filter feeding mechanisms, with distinct modes occurring in different mysticete lineages , . Bowhead and right whales (Balaenidae) are continuous ram feeders that exploit patches of copepods at slow and steady speeds , using massive tongues to direct incoming prey-laden water along parallel racks of baleen , . Gray whales (Eschrichtiidae), on the other hand, possess relatively smaller, but more mobile tongues that can be depressed to suction-feed on benthic invertebrates . In contrast to both balaenids and eschrichtiids, rorqual whales (Balaenopteridae) exhibit highly expendable ventral pouches along with extensible tongues that invert into capacious sacs, to enable the engulfment of a large and discrete volume of prey and water , . After the jaws have closed around the engulfed water, prey is filtered using plates of baleen as water is purged from the now inflated ventral pouch. Although different rorqual species exhibit a wide diversity in prey preferences (fish, copepods, squid, crabs, etc.), foraging strategy –, and ecological niche , , krill represents one of the most common prey types . In general, engulfment is preceded by prey-approach, and then followed by an obligatory filter phase which may also allow some degrees of metabolic recovery. Depending on the depth and quality of the krill patch, several lunge feeding sequences of pre-approach, engulfment and filtering/recovery may be repeated during a single dive.
Lunge feeding requires a coordinated suite of anatomical and biomechanical adaptations , . In addition to the foldable (and muscularized) ventral pouch –, the rorqual engulfment apparatus is comprised of flexible jaw joints that enable the rotation of the mandibles, which directly increase the area of the mouth ,  and therefore the flux of water into the oropharyngeal cavity . The engulfment capacity of the ventral pouch is ultimately limited by the mechanical properties of the ventral groove blubber (VGB), a specialized blubber layer that is reversibly extensible up to several times its resting length . Both the hard and soft tissue structures of the engulfment apparatus exhibit positive allometry, whereby the dimensions of these elements are relatively longer in larger animals . As a consequence, mass-specific engulfment capacity increases with body size both within and among rorqual species . Because larger rorquals have the ability to engulf relatively larger volumes of prey-laden water, overall feeding efficiency is significantly increased in larger whales .
Hydrodynamic modeling of engulfment parameterized with kinematic and morphological data indicates that lunge feeding comes at high energetic costs which are largely incurred from the engulfment and entrainment of a very large amount of water –. Such high energetic expenditures are significant enough to greatly reduce, in comparison with non lunge-feeding whales of similar size, the diving time in which foraging occurs . This phenomenon has previously been addressed in the context of the Theoretical Aerobic Dive Limit (TADL), which is largely a function of the ratio of the metabolic expenditures during a dive to the oxygen storage capacity in blood, muscle and lungs , , . What TADL doesn’t explain, however, is the fact that even though all large rorquals have similar maximal foraging diving times, namely, anywhere from 9 to 11 min. in humpback (Megaptera noveaangliae), fin (Balaenoptera physalus) and blue (Balaenoptera musculus) whales , their capacity to execute many lunges during a single dive declines significantly with body size, i.e., from 12, to 7, to 4 lunges per dive in those same three species respectively . Although there appears to be enough oxygen storage capacity to meet the energetic demands of foraging at all body sizes and to the same depths, there seems to be a more proximate limit associated with the lunge component of foraging at the largest sizes, in particular with the rate of energy expenditure during the engulfment stage. Here we explore the hypothesis that reduced lunge frequency at the larger sizes is, at least for these three large Rorqual species, a reflection of the metabolic power requirements of engulfment becoming so high to involve significant oxygen deficits and muscle fatigue during each lunge. Such increased power requirements follow from the allometry of the skull which permits the engulfment of ever greater masses of water and prey (relative to body mass ) and at greater energetic costs, but as performed over nearly unchanging engulfment durations constrained by the escape time scales of the prey . Reduced diving capacity and lunge frequency during foraging has major consequences for rorqual ecology and evolution because it begins to decrease the prime benefit of bulk feeding, namely a high energetic efficiency, by limiting access duration to high quality prey resources at depth .
Understanding the relative importance of total energetic expenditures versus power output (the rate of expenditure) requires the consideration of the various time scales characterizing the stages of a lunge. Within the context of a single foraging dive where durations and expenditures for diving, prey-approach, engulfment and filtering/recovery are all accounted for, the (averaged) metabolic rates and corresponding energetic expenditures are found to be only slightly higher than those of steady swimming . When extrapolated over several months of intense foraging during the summer, these power requirements of foraging are still low compared to the rate of energy intake such that it facilitates the deposition of substantial lipid stores needed for long distance migration and reproduction . However, and at the smaller time scale of the engulfment stage (of about 10 seconds or less), there is a requirement for short bouts of intense muscle activity and for a metabolic output during which a large body of water (both external and internal) needs to be quickly set into motion, a process for which muscle fatigue and consequent recovery may become a limitation. Although previous research has focused on estimating the energetic cost of lunge feeding , , , scant attention has been paid to the power output required by this extreme feeding strategy. As a result, both its ecological and evolutionary implications remain poorly understood. Here we address this dearth in our knowledge of rorqual foraging energetics with a new, high time-resolved hydro-mechanical model based on previous work , , but now capable of accurately distinguishing the changing physics of the various sub-stages of engulfment.
The obvious impossibility of studying the energetics of large whales in a laboratory setting makes computer modeling the only tool available for assessing the relevant factors driving the metabolic requirements of engulfment. These include the work done by the locomotor muscles for swimming and by the musculature embedded in the VGB for accelerating the engulfed water mass. An additional factor is the large amount of energy a whale typically loses to hydrodynamic drag. Given the substantial speeds imparted to the engulfed mass, and the need for coordinating the dynamics of a lunge in order to optimize engulfment volume , engulfment modeling must also include a hydrodynamic model that is coupled to the forces generated by VGB musculature. Here we quantify the effects of these factors over the adult sizes of humpback, fin and blue whales, as well as at one size of the much smaller minke whale (Balaenoptera acutorostrata), all investigated at a level of detail that we believe yields the correct physics across all body sizes and over the time scales of prey escape and engulfment . This model now accounts for the water being engulfed anterior to the temporomandibular joint (TMJ) in between the mandibles, in addition to the water engulfed posterior to the TMJ , as well as the effects arising from the wake re-contacting a decelerating whale during mouth closure. The new model thus removes some of the uncertainties that required input sensitivity analyses in previous studies ,  and which inevitably reduced the prediction ability of the model.
The simulated engulfment metabolic rates (EMR) are compared with various standardized measures of metabolic effort, including the Basal Metabolic Rate (BMR) and the Rorqual Average Active Metabolic Rate (RAAMR), a new assessment of active metabolism specific to balaenopterids. New estimates of the metabolic expenditures sustained during the prey-approach (PAMR) and the filter/recovery stages (F/RMR) will be discussed as well. Our analyses show that the metabolic expenditures of engulfment (averaged over mouth-opening time scales) significantly change with body size, from about equal to RAAMR at 10 m (a medium-sized humpback whale or small fin whale), to about 3.7 RAAMR at 27 m (a very large blue whale). Generally, the expenditures sustained prior to, and following engulfment (i.e., PAMR and L/RMR) are similar to RAAMR within 50%. Comparisons of maximum instant EMR (rather than time-averaged EMR) show an even steeper power requirement: namely, and when compared with the BMR of terrestrial mammals of the same mass (BMR(terr)), an increase from about 8 BMR(terr) at 8 m body lengths, to 48 BMR(terr) at 27 m. When further compared in terms of maximum aerobic capacity(as represented by the Maximum Metabolic Rate (MMR)), the magnitude of maximum instantaneous EMR suggests substantial oxygen deficits at the largest body sizes (27 m), with expenditure rates exceeding by 20% during a good portion the mouth opening stage. Characterized as a supramaximal type of effort , we suggest that such expenditure levels can only be met by significant contributions of anaerobic metabolism for which muscular fatigue is more important and metabolic recovery longer (i.e., recovery after each lunge during filtering and after each dive). Interestingly, a simulated 33 m blue whale shows required power outputs as high as 80% above, and maintained over durations that are long enough to explain why it is perhaps non-extant. On the other hand, engulfment metabolism requires significantly lower outputs, namely 0.1 to 0.5, at the small body sizes of all three large species and at the one simulated size of the minke whale, the smallest of the Rorquals (with body lengths of less than 10 m). As often characterized as light submaximal and heavy submaximal , such effort levels can be sustained for longer periods of activity and over several lunges during a single dive. These drastically reduced metabolic demands at smaller body sizes are similar to those of non-feeding swimming, and bring about interesting questions about the ontogeny of extant rorquals, as well as the evolution and physiological limits to different lunge feeding modes.
Materials and Methods
Engulfment metabolic rate calculations are based on a simulation of the forces at play during mouth opening and closure. Determination of the metabolic rates follows after tallying the corresponding (mechanical) energy and power while taking into account assumed metabolic efficiencies. The hydrodynamic model used for engulfment simulations, the Basic Lunge Feeding model (or BLF for short; version 3.0), is a significant upgrade of a model originally devised a few years ago , . Although the BLF dynamically couples both structure (whale) and fluid (flow), its one-dimensional fluid dynamics and simple (bulk) VGB force modeling makes it far less complex than computationally detailed models where changing body shape morphometrics is digitized into structure meshes about which three-dimensional fluid simulations are performed , or for which structural (elastic) strains are computed with finite element methods . Simplicity is allowed here by the physics of lunge feeding, where most of the relevant fluid masses being displaced are moving in only one direction, i.e., forward, along the whale’s trajectory. Moreover, the fluid-pushing forces supplied by the whale’s VGB are active pulls by muscles of (quantitatively) unknown distributions and densities, rather than elastic-based strains of a well-defined rigid structure (the skeleton).
The general principles, upgrades and approximations of the BLF model will be summarized in the following paragraphs and its mathematics discussed in Modeling Details (which follows the Conclusion). Estimates of the metabolic output sustained during prey-approach, as well as that of the active metabolic rate specific to Rorquals, will also be discussed below while the filter/recovery metabolic rates will be described in the Discussion. A list of the symbols and acronyms can be found in Text S1.
Input parameters and scope.
The BLF is informed by a variety of inputs, including six entries based on the known morphology of the VGB and skull (Figure 1 and Tables 1, 2, 3, 4). On the other hand, one table input, as well as the validation of the speeds calculated at each species’ average body size, have relied on the velocity data collected by tags deployed on humpback, fin and blue whales lunging at depth along horizontal or inclined tracks (both “uphill” or “downhill” tracks; see discussions in , , , ). Finally, the overall sequence of body shapes being simulated ,  is similar to what is shown in the film and photographic record of rorquals feeding on krill at or near the surface. Although these data provide many useful insights, they obviously limit the applicability of the model to lunges towards slow-moving prey along horizontal or straight but inclined trajectories. As observed with many rorqual species, lunge-feeding is a plastic behavior where diverse combinations of lunge speeds, body-rolling maneuvers, motion headings and gape angle dynamics are used to exploit different prey types and distributions , , including the bubble-netting technique used by humpback whales to lunge-feed vertically towards fish . However, the types of lunges discussed here represent the most common foraging strategy among rorquals and consequently their simulations should provide a realistic assessment of engulfment energetics.
The model uses the (coupled) Newtonian motion equations of the whale body and engulfed water, as constrained by an engulfment scenario in which, and as suggested by the film record, the filling of the cavity sections posterior and anterior to the TMJ are sequential rather than simultaneous (see Diagram 1 in ). In what is described here as Synchronized Engulfment (SE), a whale is assumed to first fill its cavity post-TMJ until maximum gape. This is followed, during mouth closure, by the filling of the buccal cavity anterior to the TMJ, as well as by the engulfed mass moving at the speed of the whale - a state of motion herein termed as equivelocity. Equivelocity is an important simulation constraint as it implies the absence of flows out of the cavity past the moment of maximum gape, as shown by the film record of surface lunges. On the other hand, the term “synchronized” in SE refers to the coordinated use of swimming muscles and engulfment apparatus (VGB muscles, temporalis, sternomandibularis) being individually modulated to ensure sequential cavity filling , possibly according to the inputs of a recently identified sensory organ located at the mandibular symphysis .
Computer simulations suggest that the expansion of the ventral pouch is actively resisted by the eccentric contraction of the VGB musculature  rather than passively by VGB elasticity . The engulfed water is thus accelerated forward from inside the open mouth and therefore represents an additional source of drag (engulfment drag) to that generated by the flow deflected around the body (shape drag). This concept of active (captured) flow control is supported by photographic evidence showing the lack of sufficient VGB distension at mouth closure , suggesting that the VGB is not stretching enough to enter the high stiffness region of the stress-strain curve obtained during in vitro tests on fin whale VGB . We note, however, that VGB elasticity could play a minor role during and after engulfment as an energy absorber against the sloshing of engulfed water within the ventral pouch.
SE provides important constraints on body and fluid dynamics which simplifies the model without unduly compromising its accuracy. It provides also quantitative relationships among the model’s dynamical variables (such as engulfment duration), in terms of not only time and body dimensions but also of escape modes of the prey. These constraints and relationships are further discussed in  and summarized in Modeling Details.
Forces at Play
The BLF model is a numerical scheme that computes the accelerations and speeds of the whale’s body (ac and Vc) and its engulfed mass (aw and Vw) as both interact with each other and with the surrounding fluid (Figure 2). The forces acting on the body consist in the following: the buccal cavity wall force (FBC), which by virtue of Newton’s third law of motion is equal in magnitude to the engulfment component of hydrodynamic drag (FED); the so-called “shape” component of drag (FSD), arising from the flow deflection around the whale’s body; the tail thrust force (Fthrust), and finally the weight-subtracted buoyancy force (Fext). The forces acting on the engulfed mass include the force FBC generated within the buccal cavity walls and mostly acting on its posterior-most section; and the “ocean-to-engulfed mass” drag Fww, acting on its ocean-facing end (Figure 2). The latter parameterize the effects of fluid pressure buildup under the palate, where the moving engulfed mass is meeting a static ocean. We note that using the force FBC acting in the bulk, rather than being unevenly distributed in magnitude and direction over the surface of the VGB, is a drastic oversimplification of reality, but one allowed by the overall one-dimensional (and incompressible) character of the moving engulfed mass (Figure 2).
As shown in the flow chart of Figure 3, the BLF model simulates the full mouth-opening and closure sequence, beginning with the calculation (and prediction) of the duration of engulfment (tengulf; see equation 8 below) and total engulfed volumetric capacities (equations 9 and 10). This information is then fed into an iteration scheme aimed at calculating the forces, accelerations and speeds applied to the engulfed mass and body via Newton’s 2nd law of motion, over each steps of a temporal sequence encompassing both mouth opening and closure. With the forces and speeds thus known, the mechanical work performed by each relevant force can be calculated at each time step and stored for the subsequent calculation of the metabolic rates.
Modeling Metabolic Power Output
Metabolic power expenditures during engulfment are calculated from the BLF-simulated forces and mechanical energies related to muscle use for swimming thrust (Fthrust) and buccal cavity push (FBC) (or engulfment drag (FED)). Note that in Goldbogen et al. , the expenditures were based on the consideration of total drag and swimming thrust accordingly to a scenario of shape drag possessing an active component (see Appendix S1 of that paper). The expended energy for each force is obtained by summing the increments ΔQTmech = Fthrust ΔXc and ΔQVGBmech = FBC ΔXc over all steps of the simulation, with ΔXc corresponding to the distance travelled during each time step (ΔXc = Vc(t) dt). The translation of mechanical energy into metabolic energy is achieved by dividing the sums ∑ΔQTmech and ∑ΔQVGBmech by efficiency factors, namely by 0.15 and 0.25, respectively . The 0.15 efficiency factor applied to swimming thrust takes into account both losses from muscle activity and propulsive inefficiency, namely 75% and 10%-loss of the total metabolic power generated respectively , ; the 0.25 factor applied to ∑ΔQVGBmech accounts for the 75% metabolic losses incurred during VGB contractions. The rate of energy expenditure is calculated by dividing the work incurred by engulfment time (tengulf), namely, as PTmetab = ∑ΔQTmetab/tengulf and PVGBmetab = ∑ΔQVGBmetab/tengulf. The metabolic power outputs during the mouth opening and closure stages are calculated similarly, but integrated during their respective duration (i.e., over tengulf/2; see Equation 8 below), and the instant EMR (labeled EMR*) integrated over time intervals of 0.1 s.
The calculation of EMR also includes an estimate of the metabolic expenditures incurred by the rest of the body (i.e., besides swimming muscles and VGB musculature). This is achieved by adding to PTmetab + PVGBmetab the metabolic output of the mass fraction X of those body parts that spend metabolic energy at rates obtained from the oft-used Active Metabolic Rate (AMR) , and the output of fraction Y that expend energy at rates similar to the Basal Metabolic Rate:(1)
Previous studies that determined body composition in baleen whales estimates muscle mass at approximately 43% of total mass, in comparison with 25% blubber, 17% bone and 12% viscera . Another study suggested locomotor musculature as representing 15% of total body mass . If VGB musculature comprises about the same proportion as locomotor musculature (i.e., ∼15%), there remains about 43%–30% ∼ 13% of muscle mass which, along with viscera, could represent about 25% of body mass that may also function at higher metabolic rates. Given the presumed intensity of engulfment, and further assuming such tissues operating at rates similar to AMR, Equation 1 would thus use X = 0.25 and Y = 0. Such values are obviously tentative but it turns out that they have minimal effects on the overall value of EMR.
We estimated AMR and BMR from allometric equations relating body mass to power, where AMR = 3 BMR , . Following previous studies of marine mammal metabolism , , , BMR is assumed herein as twice that of terrestrial animals expressed via Kleiber’s scaling formula , . Those two assumptions thus lead to:(2)
AMR-costs related to non-resting metabolic outputs by rorquals have been estimated either via assumed allometric formulations of AMR and BMR , , , or by assuming a specific form of the shape drag force . These AMR-based studies have yielded estimates that differ by as much as 100%, depending on the specific AMR-BMR relationship and BMR-model being used, as discussed by Croll et al. . The drag-based approach has yielded estimates lower than all AMR-based estimates, most likely due to an inadequate drag model that neglects the effects of surface-induced drag (i.e., ventilation and wave drag ), effects which are important in other marine mammals . In the case of blue and fin whales, Equation 2 yields active metabolic rates that are lower than Lockyer  by about 50%, and higher than Croll et al.  by 50%. Further support for Equation 2 is presented below with another measure of active metabolic rate that is more specific to the bio-mechanics and ecology of rorquals.
Rorqual Average Active Metabolic Rate (RAAMR) and Prey-Approach Metabolic Rate (PAMR)
Rorqual Average Active Metabolic Rate (RAAMR) is primarily based on a calculation by Bose and Lien  of the thrust generated by actual flukes using classical 2-dimensional airfoil theory with finite span corrections. Under an assumption of steady travel, power expenditures are calculated from the product of this calculated thrust (Fthrust steady) and average speed of transport (<Vc>), with the latter obtained here from long duration tracking data , . In what follows, <Vc> = 2.16 m/s and 2.4 m/s for fin and blue whales respectively. Being extracted from data spanning days, these speeds reflect motions characteristic of not only long distance transit, but also of shorter events such as resting, mingling (socializing) and lunge-feeding. Although these activities are generally performed at different speeds, the long duration of the monitoring should be dominated mostly by the long-distance travel component of the track.
The calculations by Bose and Lien applied only to the fluke shape and size of a 14.5 m fin whale and, as such, needed to be extrapolated over body size and across species. This was done by first scaling the Bose-Lien result according to body drag, and thus to surface area (or, equivalently to body area Abody (Figure 1 and Tables 1,2,3 4). As non-feeding swimming is often performed in groups of variously-sized individuals moving at the same speed, the transport speed <Vc> is herein assumed as body size-independent. Moreover, non-feeding swimming speed tends to be relatively low across body size in many different taxa , . The Bose-Lien results were also scaled across the body length of blue whales given their morphologic similarity with that of fin whales. This analysis was not extended to humpback whales given their significantly different fluke morphology.
A third extrapolation of the Bose-Lien thrust was performed with respect to travel speed, given the high flow speed these authors used (i.e., 4 to 12 m s−1) in comparison to the average travel speed measured by tracking. Considering the Bose-Lein data below 8 m s−1 and extrapolating down to 2 m s−1 yields the formula Fthrust steady = 1000 (1.216<Vc> +0.0205<Vc>2). The RAAMR that results is computed using a metabolic output model similar to that of Equation 1:(3a)where
(3b)The value of the cross sectional area at 14.5 m (or Abody (14.5 m) ) is 2.81 m2 for both fin and blue whales, and follows from the allometry of Abody shown in Tables 1, 2, 3, 4. The body fraction are set to X = 0.15 and Y = 0.2 but as with Equation 1, their specific choice yield small contributions in comparison to those of Ptransitmetab. Note the efficiency factor of 0.15 discussed previously. Note also that <Vc> = 2.16 m s−1 and 2.4 m s−1 are applied at all body length for fin and blue whales respectively. As shown next, the RAAMR exceeds AMR by factors of only 1.1 to 1.6.
The metabolic rate incurred during prey-approach (PAMR) is calculated by assuming fluking thrust as being much greater than drag and (weight-adjusted) buoyancy (see figure 2 but without the engulfment-specific forces). With thrust being the only force at play, the (metabolic) power output of fluking is computed from the change in (whale) kinetic energy during a stage of duration tpa, namely from the end of the filter/recovery stage of the previous lunge to the beginning of engulfment of the current lunge. Factoring in the metabolic and fluke hydrodynamic efficiency, as well as the “rest-of-the-body” expenditures, one has:(4)
The fractions X and Y are set to X = 0.25 and Y = 0 as with EMR. Values for the whale’s mass (Mc) and speed at the beginning of engulfment (Vc(0)) can be found in Tables 1, 2, 3, 4. The values for tpa are shown Table 5 and extracted from tag studies , , , . These correspond only to each species’ average body size. Finally, the speed at the end of the filter/recovery stage (Vc(tfilter)) is set at 1.0m /s for the three species, again as hinted by digital tag studies , , , .
This BLF upgrade now requires, for each value of Lbody, a total of 22 inputs on body dimensions, body dynamics and hydrodynamics (Tables 1, 2, 3, 4). The former include the relevant body characteristics Lbody, L0, whead, Ljaw, Mc and Abody (Figure 1) which were obtained from reduced major axis regressions  of morphometric studies , – in humpback, fin and blue whales. These regressions yield allometric equations of each body characteristic over body sizes ranging from weaned juveniles to the largest adults. Simulations were performed for minke whales as well, but at only one body size due to limited data.
As the most important body dynamics input, the initial speed of a whale just prior to the mouth opening stage (or Vc(0) in Tables 1, 2, 3, 4) was determined from tag studies of large rorquals lunge-feeding at depth for krill . These data, along with the modeling of possible escape scenarios of the prey , suggest this initial speed to vary as Vc(0) = Vn Lbody, with Vn = 0.159/sec (humpback), 0.148/sec (fin) and 0.146/sec (blue) . Given the unavailability of tag data in the case of minke whales, the initial speed was assumed at Vc(0) = 1.16 m/s (as computed from Vn = 0.15/sec), a value which turned out lower than the reported average minke cruising speeds (∼ 3 m/s ), but one that appears nevertheless within the range of this species’ feeding speeds.
The values of the dynamical and hydrodynamic parameters kopen, φ, χ, CD etc. listed in Table 1 are further discussed in Modeling Details. Note that with minke whales, and again due to lack of tag data, these (dimensionless) parameters were set at values typical of fin and blue whales given their similar morphology. An exception was with the maximum gape angle, which was set to 50° instead of 78°. BLF simulations, along with the film record, seem to suggest the impossibility for minke whales to carry out horizontal engulfments without premature cavity filling or draining at the maximum gape angle of 78° used by the larger Rorqual species (while lunging on krill).
Body Motion, Engulfed Mass, Force Output and Expended Power
The simulation diagram sketched in Figure 3 yields calculations of the engulfed mass, muscular-based forces and whale body speeds, and ultimately of the energy expenditures. Sample outputs are shown below in the case of fin whales (outputs for humpback and blue whales are qualitatively similar). The scaling laws relevant to the graphed results are summarized in Tables 5, 6 and 7 for each species.
Simulations of a 20.2 m fin whale are shown in Figures 4 and 5. These are compared in Figure 4 with averaged speed data collected by digital tags , as 20.2 m fin whales represent the average body size of this species . The fact that the simulated lunges reproduced the measured speeds shouldn’t be too surprising given that several dynamic input parameters were tuned to yield a good match, in particular that of Fthrust (at t < topen). Such parameter tuning was carried out only at the average body size of each species (where the tags data applies), and then scaled with respect to Lbody at all other body sizes according to scaling rules further discussed in Modeling Details.
Simulated speed of a 20.2 m fin whale (continuous line), as compared with data from digital tags  (triangles). The dashed line corresponds to the simulated speed of the engulfed mass. Note the equivelocity constraint operating during mouth closure (i.e., at t >2.85 s).
The forces of engulfment are shown in Figure 5. Comparing with a similar figure produced in a previous version of the BLF shows overall qualitative agreement (see Diagram 10, in Potvin et al.  or Diagram 5 in Goldbogen et al. ). In this new version, however, engulfment drag is dominating over shape drag and fluking thrust during the entire engulfment sequence. Moreover, shape drag becomes very small as a result of the re-contact of the wake onto the whale’s decelerating body. In fact, shape drag can become negative at sufficiently large body size, i.e., a pushing force rather than a resistance (as illustrated further in Modeling Details), in analogy with the wake pushing a power boat upon engine cut-off. But here the effects of wake re-contact on the forward motion are mitigated by fluke thrust which then becomes negative, i.e., when producing a braking rather than a propulsive action. Interestingly, such wake re-contact affects only the larger sized whales given their greater decelerations, a result of having to push forward a larger engulfed mass (relative to body size) as caused by the positive allometry of the skull (more on this below).
Temporal variations of total drag (FSD + FED; continuous line), engulfment drag (FED; dash-dotted), shape drag (FSD; dotted) and combination of fluking thrust and weight-adjusted buoyancy (Fthrust + Fext; dash-doubly-dotted). Mouth closure begins at the 2.85 s mark and is characterized by shape drag cancelling Fthrust + Fext per equation 24.
The variations in intra- and inter- specific scaling among the humpback, fin and blue whales can be summarized by the (mass-specific) peak engulfment drag (i.e., FED at t = tengulf/2 (the time of maximum gape)) and engulfment time tengulf shown in Figures 6 and 7 (see also Table 6). Note that Figure 7 also shows the so-called VGB contraction time scale τ (in contrast to engulfment duration), which is about three times as small, and in the range of 1 to 2 seconds over the entire body dimensions of humpback, fin and blue whales, as further discussed in Modeling Details.
Largest value attained during mouth opening, by humpback (diamonds), fin (squares), blue (triangles) and minke whales (starburst).
Figures 8 and 9 show the calculated metabolic expenditures (mass-specific) during non-feeding transport (RAAMR; Figure 8) and in the course of lunge-feeding (EMR; Figure 9). Unlike the metabolic expenditures of engulfment, RAAMR decreases slightly with body size from about 2.7 to 1.7 W kg−1 (see also Tables 7 and 8). Thus, RAAMR is consistent with the general concept of a lower cost of transport , although we note that RAAMR is somewhat larger than AMR (equation 2) by factors of 1.1 to 1.6. As calculated and applied to averaged sized fin and blue whales, RAAMR is greater than BMR by factors of 4.8 and 3.6 respectively.
In contrast, the estimation of the EMR increases with body size as demonstrated by previous studies . But these simulations show for the first time that metabolic output is highest during the mouth opening stage when the cavity wall forces are at their peak (Figure 5). Overall, such intensity ranges from 3 to 7 W kg−1, where the largest expenditures are incurred at the largest body sizes. This is a reflection of the positive allometry of the skull which enhances buccal cavity volume (relative to body size) and results in larger mass-specific engulfed and deflected water masses .
Figure 9 also compares the metabolic rates involved in each phase of a lunge, namely, prey-approach (PAMR), mouth opening and closure (EMR) and filter/recovery (F/RMR). Although PAMR and EMR are calculated directly in terms of the forces involved (see Eqs. 1, 2, 4), F/RMR is estimated only indirectly, i.e., via PAMR, EMR and RAAMR (Equation 3), and from the average metabolic output (ΔEO2/tmaxdive) that could be obtained from a “maximum duration” foraging dive in which all of the stored oxygen would (theoretically) be used to achieve a maximum number (fmaxlunge) of lunges:(5)
Parameters tpa, tengulf and tfilter correspond to the duration of the prey-approach, engulfment, filter/recovery stages respectively, and tascent, tdecent, tsearch and tmaxdive to the times needed for descending and ascending to and from foraging depth, searching for the krill patch in between lunges (during a same dive), and for performing the longest dive possible with the O2 stores at hand (O2 stored in the lungs, muscle and blood). With the exception of tengulf, most of these parameters are obtained from tag data and listed at the representative average body length of each species in Table 5. The physiological data pertinent to the computation of ΔEO2 are supplied in Table 9. This estimate of the F/RMR is similar to a TADL calculation , ,  except that tmaxdive is not predicted but rather used as an input to yield a metabolic rate. This is a rather speculative estimate given that several input parameters are still poorly known. Nevertheless, the results suggest that the metabolic requirements of engulfment (mouth opening specifically) significantly exceeds those of the other stages of a lunge by 100% or more. Moreover, both PAMR and F/RMR still represent “active” metabolic outputs (in contrast to “basal”), as they turn out to be similar to those of RAAMR. This shouldn’t be too surprising even for the filter/recovery phase, as the significant oxygen deficits incurred during mouth opening stage are likely to require high recovery metabolism.
Ratios calculated for humpback (diamonds), fin (squares) and blue whales (triangles). Engulfment Metabolic Rate (EMR) during mouth opening – black symbols, and during mouth closing – gray symbols. Metabolic rates for prey-approach (PAMR), basal (BMR), filter/recovery (F/RMR) and Rorqual Average Active (RAAMR) – open symbols as indicated.
Body size is one of the most important determinants of energetic efficiency and locomotor performance –. Because locomotion and prey capture are integrated in rorquals, as they are in many vertebrate taxa , both elements are subject to mechanical scaling effects. Here we integrated morphological and kinematic data into a novel hydro-mechanical model, derived from first principles and from a unified theory of predator-prey dynamics , to estimate the forces required for lunge feeding whales. The model output enabled us to quantify the energetics of feeding across an extensive size range of juvenile and adult rorquals. In general, the energy and power output required to lunge feed increases disproportionally with increasing body size, a phenomenon that results from a complex interaction between lunge speed, unsteady hydrodynamics and the allometric scaling of the engulfment apparatus . Because feeding is such a major component of baleen whale life history, these predictions have major consequences for rorqual foraging ecology, ontogeny, and evolution.
The Metabolic Cost of Engulfment
The results of Figure 9 show how significantly more strenuous engulfment is during mouth opening, in comparison with BMR and non-feeding transport. This is emphasized further in figure 10, with a direct comparison with RAAMR at all body sizes in fin and blue whales. Unlike engulfment metabolic expenditures, RAAMR (mass-specific) decreases slightly with body size (Figure 8; see also Tables 7 and 8). This is similar to the Active Metabolic Rate (AMR) of other marine mammals which is approximately three times their Basal Metabolic Rate , –. But the ratio of EMR to RAAMR, much like all other ratios that relate engulfment costs to all other energetic expenditures, increase with body size for all three species up to about 3.7 RAAMR (Figure 10). Most noteworthy is the convergence of EMR and RAAMR at small body sizes, and the conclusion that the metabolic expenditures during engulfment become quite close to those of non-feeding swimming at such sizes. At large body sizes on the other hand, and including the sizes of the non-extant blue whale (33 m), the much greater costs of engulfment become evident.
Ratios calculated for fin (squares) and blue whales (triangles). The 33m blue whale is non-extant.
Engulfment metabolism relative to other standard metrics of power is shown in Figure 11, with a comparison of EMR with the BMR and maximal metabolic rate (MMR) of terrestrial mammals of the same mass (MMR data from ; and BMR(terr) = 4.0 Mc0.75 (in watts)). Marine mammals are different from terrestrial mammals, of course, with their heightened resting metabolic rate on the one hand , and their well-known physiological adaptations for overall BMR reduction during diving on the other (for example, bradycardia and regional vasoconstriction) . But the data available on mitochondrial volume densities  and aerobic scope  on seals and dolphins already point to a similar aerobic metabolic performance for powering locomotor musculature (even during breath-hold), which after all, is the main function of metabolism during strenuous exercise . Thus the physiological demands of lunge-feeding rival those of athletic , or more appropriately “highly active”, terrestrial mammals during maximal effort. What is also interesting is the finding that the relative level of power required by engulfment increases significantly with body size, from that similar to trained human athletes at the smaller sizes (MMR ∼ 20 BMR), to the maximum performance of dogs and horses at the larger sizes (MMR ∼ 30 BMR). Although interesting, this comparison actually underestimates the actual maximum metabolic output that may be required, as EMR represents an averaged metabolic rate, i.e., over the duration of mouth opening. A more informative datum is the “instantaneous EMR” (or EMR*) computed from averaging power output over the significantly shorter time interval of 0.1s, to become time-dependent as shown in Figure 12 (Note: EMR* also neglects the X- and Y-terms using in equation 1). EMR* indeed yields significantly higher maximum values, i.e., EMR*|max/BMR(terr) = 5.8 and 18.2, 12.6 and 39.8, and 24.3 and 48.2, at the smallest and largest sizes in humpback, fin and (extant) blue whales respectively. Instant EMR thus changes the effort picture of Figure 11 somewhat, suggesting the largest extant Rorquals as engulfing prey and water at performance levels rivaling those of the pronghorn antelope. Interestingly, the predicted metabolic performance required from the non-extant 33 m blue whale far exceeds that of the pronghorn antelope at EMR*|max/BMR(terr) = 78.6.
Ratios calculated for humpback (diamonds), fin (squares), blue (triangles) and minke whales (starburst). Note that the 33 m blue whale is non-extant.
Maximum aerobic capacityis assessed with a comparison of instant metabolic rate with the Maximum Metabolic Rate (i.e., MMR = 39.4 Mc0.87 (in watts )), as listed in Table 8 and shown in Figure 12. The figure shows the metabolic demands of engulfment in relation to MMR, which is a direct measure of the maximum aerobic capacity of an air-breathing terrestrial mammals , or in other words, of the limitation of oxidative metabolism of muscle cells to supply energy without recourse to anaerobic glycolysis . Exceeding MMR as shown reveals the magnitude of the oxygen deficit that accumulates during the most demanding phase of engulfment. But oxygen deficits also accumulate from the start, i.e., at small effort level, given the short time scales of engulfment, and also for fueling the anaerobic metabolic reactions that are needed to power the forward push of the engulfed mass by the fast twitch muscle fibers embedded in the VGB (such fibers represent approximately 50% of the VGB muscle fibers [R.E. Shadwick, unpublished data]). Accordingly, overall oxygen deficits may be small and recovered during the latter part of the effort, as during typical light submaximal exercise ; or greater, i.e., as during heavy submaximal exercise, and mostly repaid during longer post-effort recovery (i.e., during filtering or after a dive); or significantly greater (supramaximal exercise), where severe muscle fatigue induced by stored glycogen depletion and inorganic phosphates accumulation can only be recovered during extended rest periods. Again, the comparison suggests that engulfment is disproportionally costly for larger body sizes. At the upper extreme of body mass in the (extant) range of 25 to 27 m, the rates of energy expenditure are most demanding, i.e., reaching MMR-levels over several seconds.
Scaling of Lunge Feeding Power Output: Consequences and Implications
Our analyses show again that the energetic cost of engulfment is positively allometric whereby larger rorquals must expend relatively more energy to lunge feed , , . But the results show for the first time that it is the rate of energy delivery (i.e., the power), rather than the total energy required for lunge feeding, that may be the limiting factor on lunge performance. This is underscored on the one hand by the TADL computed from our simulated energies (and physiology data of ), hinting at maximum dive times in the range of 9 to 12 minutes at all body sizes, including that of the non-extant 33 m blue whale; and on the other, by the power requirements reaching levels of high muscle fatigue and long post-dive recovery. Such a heightened power requirement is a result of the positive allometry of the engulfment apparatus , where larger skull sizes (in relation to body size) yield greater engulfment capacity and increased captured mass kinetic energy costs, as performed over (nearly) unchanging engulfment time scales brought about by the evasion strategies of the prey .
Increased fatigue and metabolic recovery may explain the reduction in lunge frequency seen across species , and presumably, among the largest individuals of each species. However, other factors such as krill patch dimensions may also be an important limiting factor, particularly if they are similar to, or smaller than predator size. But the noted increased post-dive breath numbers in humpback whales after dives encompassing the largest number of lunges  may provide supporting evidence for intensified power delivery rather than for food availability. Interestingly, such ability to deliver the required power may be adversely affected for emaciated whales returning to their feedings grounds after months of fasting, and as a result, experience lower lunge frequencies than at their normal body weight. This is an interesting topic for future simulations of course, but one that will be possible only when the morphology and body mass of these underweight and slimmer animals become known.
Limits on power delivery predictably has a negative impact on filtering capacity over an entire dive , and therefore feeding efficiency , but such power demands of engulfment may also set a hard physiological limit on maximum body size in rorquals. We tested this hypothesis by simulating lunge feeding in a hypothetical 33 m blue whale, and the results indicate that its engulfment power expenditures exceed those of MMR by over 80%, with a peak instant EMR of about 78.6 BMR(terr). In other words, blue whales at this scale that would be able to generate the needed power for engulfment would do so at a cost of requiring extensive rest and recovery even after a single lunge, to obviously limit residence time in high density prey patches at depth. In this way, and as argued for dinosaurs  and blue whales , limits in body size may not only be imposed by a balance between resource availability and energy expenditures, but also by the metabolic power output required for feeding (see also this argument applied to non-feeding transport in ). Nevertheless, either circumstance (power vs. energy limitation) predictably constrains maximum body size in this lineage of baleen whales. However, it is unknown what limits size in other cetacean lineages, as well as other lineages of other aquatic tetrapods that also exhibited gigantism in the past . Clearly more research is needed to explore what factors limit body size in different taxa not only from an energetics perspective, but also with respect to other life history constraints.
Our analyses, including those of the minke whale (a fish- as well as a krill-feeder), also have implications for the lowest body size classes of rorquals, which suggest relatively low power requirements for lunge feeding. This has a significant impact on the ontogeny of rorquals since weaned juveniles reap the advantages of efficiency in a low-cost feeding strategy that will facilitate rapid growth. Such a characteristic, which is exemplified by the large range in intraspecific body size of extant rorqual species, likely played a major role in the evolution of gigantism in balaenopterids. Lower energetic costs associated with engulfment at smaller body sizes may also have implications for how lunge feeding evolved in smaller ancestral baleen whales. Fossil evidence between the late Oligocene and late Miocene indicate that extinct rorquals reached maximum body sizes no larger than extant minke whales (i.e. 10 m long) . These data, together with our minke results (Figure 11), suggest that lunge feeding evolved within a body size range where the cost of engulfment does not appear to be significantly higher than non-feeding swimming (Interestingly, the minke modeling suggests that such costs may be even smaller). We posit that the evolution of baleen in cetaceans  at small body size classes  generated a mode of feeding that exhibited high energetic efficiency where vast amounts of prey could be captured at a relatively low cost . Such an ecological role further led to more specialized types of microphagy, and therefore to several mysticete lineages that each exhibited distinct filter feeding modes, that satisfied particular niches related to differences in prey type. The ability to gulp discrete volumes of prey-laden water likely facilitated the exploitation of more agile zooplankton (i.e. krill, squid, etc.), in contrast to slow moving copepod prey for which continuous ram feeding may be more efficient , and therefore promoted the evolution of extremely large rorqual species fueled by the existence of super-aggregations of krill . Future research should further explore the feeding energetics of the smallest baleen whale species, both extinct and extant, to examine how different feeding modes may have evolved and to what extent it promoted extreme body size.
This paper has shown that the physical requirements of lunge-feeding towards krill, particularly with regards to the rate of energy delivery during engulfment, may present an obstacle to ever increasing body size. This would follow from the allometry of the skull, along with the approach speeds being dictated by the evasion strategies of the prey. But this obstacle may, of course, be altogether avoided if lunges can be performed differently, for example by reducing the maximum gape angle, or by passively engulfing while lunging vertically and/or cooperatively, as already performed by some rorqual species lunging at the surface. Furthermore, limits to body size are also connected to prey type and availability. Large body size obviously limits maneuverability and has most likely forced the blue whale, the largest of the Rorquals, into obligate krill-feeding and into exploiting a resource that has so far been abundant enough to even support these giants in the very large numbers of pre-whaling days. But ultimately, and given that most Rorquals prey on both krill and fish, the linkages between food resources, motion energetics, and body size will not be completely resolved until further modeling is carried out on lunge-feeding towards fish. In general, lunging on fish will demand significantly higher prey approach speeds and accelerations by the predator. In most rorqual species such extra effort will be mitigated by the use of smaller gape angles (as with minke whales), as well as by other maneuvering and dynamic adaptations such as vertical lunging without active VGB muscle action. Although vertical lunging at higher speeds and with smaller gapes can already be simulated with the modeling techniques discussed here, passive engulfment is bound to change the physics and hydrodynamics of the modeled whale-ocean system, and as a result, yield an alternate simulation tool for the study of lunge feeding.
Body and fluid dynamics constraints from synchronized engulfment.
Synchronized Engulfment introduces several constraints on overall engulfed mass motion as well as on specific dynamic variables. The first concerns the eccentric contraction of the VGB musculature, which must impart forward motion to the engulfed mass (as a reflux) without premature draining of the cavity . Here premature cavity draining occurs whenever the engulfed water leaves the mouth aperture prior to complete mouth closure. Moreover, the BLF implements an assumption of premature filling avoidance, or in other words, of preventing complete cavity filling post-TMJ prior to maximum gape, and complete cavity filling ant-TMJ prior to mouth closure. These two constraints are implemented by the use of cavity wall force pushing the engulfed mass with the “right” amount: If the whale exerts a force that is too high, the result is a more rapid reflux, slower whale motions and thus a slower cavity inflation rate or even a negative inflation rate (or draining). Conversely, if the force is too small the reflux is minimized, but the whale speed remains high and the ventral pouch fills too quickly. Note that an extreme of the latter includes “passive” or “compliant” engulfment , where the influx of water is met with little resistance other than the passive mechanical (i.e., elastic) properties of the VGB. In this case, and given the compliant nature of the VGB over most of its allowed strain range , little force is exerted on the engulfed mass during most of the engulfment process.
Because the oropharyngeal cavity (post-TMJ) does not significantly fill during mouth closure, the model invokes a third SE-motivated constraint, namely that of the engulfed mass post-TMJ moving at the instantaneous speed of the whale. This so-called equivelocity constraint is new to the BLF and is further extended to the water captured anterior to the TMJ, as motivated by the water being “scooped” by the rotating mandibles during mouth closure rather than “bagged-in” and pushed forward as during the mouth-opening stage. Most importantly, equivelocity implies the absence of significant surging of engulfed water moments before complete closure, a state confirmed by the film record. Equivelocity, and its consequent equiacceleration, means that the interaction between whale and water represents a perfectly inelastic collision, and leads to a useful derivation of an engulfment hydrodynamic force that would otherwise be very difficult to determine (i.e., Fww in Figure 2). Furthermore, equivelocity provides an additional constraint on drag and thrust and as such helps reduce the uncertainties that have plagued the modeling of previous studies –. The equivelocity of whale and engulfed fluid motions during mouth closure, as well as the acceleration of the reflux during mouth-opening, are clearly visible in Figure 4. (Note that there is no reflux in passive engulfment , ).
Algorithmic flow of the BLF.
The hydrodynamic model is an iterative scheme aimed at computing the forces applied to, and speeds sustained by both whale and engulfed mass. As shown in the flow chart of Figure 3, each iteration involves a calculation of the gape angle, mouth surface area (projected longitudinally) and mass so-far engulfed, which are then used in the computation of all the forces acting on the whale’s body and engulfed mass. The upgrades of this BLF version include the use of engulfed mass rates that are specific to the mouth opening and closing stages – an important ingredient for ensuring physical realism at the needed time scale of metabolic output; other new features implement a formulation of shape drag aimed at approximating the effects of wake re-contact on the whale’s body during mouth closure, and the use of the equivelocity constraint to derive the “ocean-to-engulfed mass” drag Fww, a force that is entirely dynamical in nature.
Mouth opening rates and engulfment duration.
The complete gape angle cycle of engulfment, from opening to closure, lasts several seconds, with the mandibles opening to a maximum gape (θgapemax) of about 78–80° with the humpback, fin and blue whales , , and of about 50° with the minke whale [JAG unpublished data]. Evidence from the film record suggests rates of mouth opening and closure calculated as :(6)(7)with θgape and t corresponding to the gape angle (Figure 1) and time respectively. Here Ljaw, Xjd, hsync, topen and Vc(0) are the length of the mandibles, jaw disarticulation factor, synchronization factor, duration of the mouth-opening stage and whale’s speed just prior to mouth opening, respectively (Figure 1).
The body dimensions used in these formulae (and everywhere else in the BLF) are shown in Tables 1, 2, 3, 4 and characterize the size classes typical of adults in all three Rorqual species . These are the results of reduced major axis regressions  of morphometric studies , –. The initial whale speeds Vc(0) at these average body sizes were obtained from tag data . The duration of mouth opening (topen) follows from the requirement of the (near) maximal extension of the ventral pouch post-TMJ by the time of maximum gape . On the other hand, mouth closure duration is assumed as being the same as mouth opening, as motivated from video footage . When integrated with Equations 6 and 7, this observation yields the means of calculating the durations of engulfment (tengulf), mouth opening (topen) and mouth closure via.(8)
The ratio hsync/Xjd can be shown to scale as hsync/Xjd = (θgapemax/Γ) sin θgapemax Ljaw/(L0 - Ljaw) as a direct result of Synchronized Engulfment , with L0 representing the axial length of the VGB (i.e., from umbilicus to mandibular symphisis) (Tables 1, 2, 3, 4). The value of constant (input) Γ listed in Tables 1, 2, 3, 4 is derived from the values of the ratio hsync/Xjd extracted from the temporal development of the gape angle shown in video footage .
Equation 8 is important as it predicts the duration of engulfment at all body sizes. It also provides an important scale for the expression of the force provided by VGB musculature (see equations 16 and 17 below). Engulfment duration is the first of two fundamental time scales that are relevant to engulfment, as further discussed in the context of cavity wall force (FBC).
Engulfed mass capacity.
The overall capacity of the ventral pouch, as well as the filling rate of each compartment (i.e. ant-TMJ and post-TMJ), are important dynamical components of the model. These volumetric capacities are expressed in terms of quarter-ellipsoidal shapes . When filled during typical horizontal lunges, and excluding local cavity over-extension effects due to sloshing, photogrammetric data suggest the ventral pouch’s extensibility to not expand in width beyond that of the skull and in depth below the length of the mandibles (at maximum gape). From these observations emerge the following “filled” capacity equations:(9)(10)
Parameters ρw and whead correspond to the density of sea water and width of the skull (Figure 1; Tables 1, 2, 3, 4). Like all other body dimensions, the latter is known from morphometric studies. The quantities in parentheses correspond to the three semi-minor radii of the quarter-ellipsoids. Interestingly, the ratio of these capacities amounts to Mwant-TMJ/Mwpost-TMJ ∼ 0.5 for all three sizes of the adult rorqual species represented in Tables 1, 2, 3. Note also that Equations 9 and 10 would not apply to the cavity of bloated dead whales for which the width well-exceeds that of the skull, a result of the decomposition gases stretching the VGB to longitudinal and circumferential strains never reached during actual lunge feeding.
Using the morphological data of Tables 1, 2, 3, 4 in the ratio of engulfed mass to body mass Mwtotal/Mc (with Mwtotal ≡ Mwpost-TMJ + Mwant-TMJ) shows engulfment capacity to exceed body mass by at least 10% at most body sizes, and up to 50% at the largest body sizes (Goldbogen et al ). Such large engulfment volumes are enabled by the mouth opening at wide gape angles, along with the unfolding and stretching of ventral pleats that line the buccal cavity wall. Those pleats and the rest of the VGB are reversibly extensible up to several-times its resting length . The volume of the cavity is also increased by the inversion and distension of a weakly-muscularized and highly elastic tongue (see Diagram 6 in Lambertsen  or Diagram 3 in Goldbogen ), which invades an intermuscular space (the so-called the cavum ventrale) located between the VGB and the rest of the body , , . The result is that at maximum cavity extension, the engulfed water mass runs ventral and posterior to the esophagus .
During mouth opening the filling rate of the cavity post-TMJ, or in other words the amount of fluid mass Mw(t) entering per unit time (or dMw/dt ), is seen as the filling of a quarter-ellipsoidal sac opened at its wide end (see Diagram 2 in ). The opening of the latter is assumed to be shaped as a half-ellipse of surface area a given by:(11)
(see also Diagrams 5a,b,c in Potvin et al ). This instantaneous area is known once the gape angle is calculated from equation 6 or 7. The water flows into the ventral pouch at a speed φ(Vc(t) – Vw(t)), resulting in the flux being given by:(12a)
Here the velocities Vc and Vw are obtained from the equations of motion to be discussed next. Parameter φ is a filling efficiency necessitated by the water entering the cavity (below the TMJ) at speeds that are, relative to the whale, higher than the speed of the ocean ahead, a result of the flow passing through a funnel formed by the palate and buccal cavity ant-TMJ. Specific values were obtained via numerical experimentation constrained to not exceed the volume of equation 9, resulting in φ = 1.60 for t ≤0.66 topen and φ = 0 for 0.66 topen < t ≤ topen (in all three rorqual species). The latter condition, which effectively stops the filling of the cavity 66% into mouth opening duration was necessary for avoiding cavity draining. This constraint leads to the cavity post-TMJ to fill to somewhere between 70% and 80% of maximum engulfment capacity.
The filling rate of the cavity ant-TMJ during mouth closure is based on a rate of (upward) angular sweep corresponding to the motion of the mandibles:(12b)where the θ-dot derivative (or dθ/dt) on the right-hand-side is obtained from Equation 7. Parameter χ represents a capture efficiency quantifying possible fluid loss during mouth closure (χ = 1 here). Note that equation 12b was not implemented in the previous version of the BLF , thus omitting the filling of the buccal cavity ant-TMJ. This omission was partially compensated by letting the cavity post-TMJ fill to 100% capacity throughout engulfment, but still resulted into a 25% underestimate in the computation of the engulfment forces (see Appendix 1 of Goldbogen et al. ).
The instantaneous value of the total mass accumulated so far (Mwtotal) is obtained from a numerical integration of Equation 12a during mouth opening, and of Equation 12b during mouth closure:
Equations of motion.
The fact that most of the engulfed water is moving in the forward direction motivates the use of the simplest fluid dynamics model possible, namely one-dimensional hydrodynamics. This model is then coupled with the trajectory modeling of the whale’s body as follows :(14)(15)
The accelerations ac and aw are those of the whale’s body (Mc) and engulfed mass, respectively. The accelerations are defined from the velocities Vc(t) and Vw(t) measured from a fixed reference frame, via ac = dVc/dt and aw = dVw/dt. The filling rate (M-dot) in equation 15 is given by either Equation 12a or 12b. The forces applied on each body all appear on the right-hand-side of both equations; these are further explained in the following sections.
As shown in Figure 3, equations 14 and 15 are used to compute the accelerations ac and aw (and then the speeds Vc and Vw). This is done once the gape angle and engulfed mass rates have been calculated via equations 6 and 12a respectively (or equations 7 and 12b during mouth closure), and also once the forces shown in Figure 2 are calculated. The numerical scheme used for this step is discussed in further details by Potvin et al. .
These equation of motion also yield a computation of the lead slug position Xc(t) – Xw(t) (i.e., the very first engulfed slug), which is a useful diagnostic for monitoring the filling state of the cavity post-TMJ. The calculation is performed via the numerical integration of Xc(t) – Xw(t) = ∫dt’ (Vc(t’) Vw(t’)). This observable is used to terminate BLF simulations (with an error message) whenever Xc – Xw exceeds the length L0– Ljaw of the VGB posterior to the TMJ, an indication that the lead slug is about to accomplish the biologically impossible feat of traveling past the umbilicus.
VGB push force (FBC) and engulfment drag (FED).
SE involves the active use of musculature within the VGB and tail in order to set (and keep) the engulfed mass into a state of forward motion. Moreover, such push-forward is metered in a manner to avoid premature cavity filling or premature cavity draining. (The former occurs if the push is too weak and the latter if it is too strong). Such control is possible morphologically given the preponderance of both fast- and slow-twitch muscle within VGB tissue, running from the umbilicus to the snout in bundles layered longitudinally and obliquely . This muscle is thought to resist the lengthening of the VGB posterior to the TMJ so to gradually impart with forward speed the mass being engulfed. Moreover, there are two large “paratendinous” cords emerging out of the mandibular symphysis and project posteriorly in a direction parallel to each mandible . Each cord lies within the buccal cavity walls, just in between the blubber and the ventral most muscle layer. This structure, termed the Y-shaped fibrocartilage skeleton , is clearly visible during feeding when the buccal cavity is inflated . It has been hypothesized that this structure adds rigidity to the buccal cavity and acts like a tendon to transmit force from the buccal cavity to the mandibles .
The push-forward of the engulfed mass is modulated so to avoid premature cavity filling or draining. Several mathematical forms of the VGB muscle action FBC have been discussed elsewhere . The most interesting candidate is expressed in terms of the product of an acceleration scale (in brackets) and a mass scale:(16)(17)
Cavity wall force gives rise to engulfment drag FED (FED = FBC) by action-reaction, and adds to the so-called shape drag FSD being produced by the flow deflected around the whale’s body (as discussed in the next section). The function Ac(t) is calculated from equations 11 and 6 (or equation 7 during mouth closure), parameter tengulf from equation 8, and total engulfed mass Mwtotal from equations 13a or 13b.
The reaction constants kopen and kclose (no dimensions) determine how “hard” the engulfed mass is being pushed forward by muscle action: too high of a value and the simulation yields cavity draining as Vc(t) < Vw(t); too small of a value and the cavity post-TMJ fills up too soon, i.e., prior to maximum gape. In this new version of the BLF model, the values for kopen at each body size (Tables 1, 2, 3, 4) are determined via trial simulations so to obtain Vc(t) = Vw(t) at the time of maximum gape, thus ending the filling of the cavity post-TMJ at maximum gape per the equivelocity constraint. Such input kopen are single-valued in that slightly-off values yield either Vc(topen) < Vw(topen) and cavity drainage post-TMJ during the entire mouth closure stage; or Vc(topen) > Vw(topen) and continued filling of the cavity post-TMJ during closure.
Numerical experimentation suggests values of kclose to be set according to ratios independent of body size, but still varies relating to species specific morphology as kopen/kclose = 1.49, 1.26, 1.82 and 1.83 for minke, humpback, fin and blue whales respectively (Tables 1, 2, 3, 4). Different values for kopen and kclose make sense morphologically given that the mouth-opening torque applied to the mandibles (mainly by the longitudinal VGB musculature) is maximum during the mouth-opening stage when buccal cavity wall forces are at their largest . It follows that such action must be much weaker during mouth closure in order to minimize that same (mouth-opening) torque that would then be working against the mouth-closing skull musculature . The required forward-push of the engulfed mass could be provided during mouth-closure by VGB oblique musculature which can generate a force without opening the mouth, as done during the filtering/recovery stage (it is noted that such oblique muscle is likely to also assist the longitudinal muscle in providing cavity wall push during the mouth-opening stage). Note that these kopen/kclose ratios do not guarantee Vc (t) = Vw(t) throughout mouth closure; this requirement is enforced instead through a definition of the force Fww applied on the engulfed mass (Figure 2) as further discussed below. Note also that in previous version of the BLF model – the value of kclose is equal to zero while that of kopen is set to yield maximum capacity Mw/Mc in the cavity post-TMJ (equation 9) by the end of engulfment.
Defining the cavity wall force as a product of an acceleration scale and mass scale leaves a certain degree of arbitrariness in the choice of these scales. Consequently, and although constrained by the equivelocity constraint, the values of kopen and kclose listed in Tables 1, 2, 3, 4 are associated with the particular choice of acceleration and mass scales used in equations 16 and 17. In conjunction with this degree of arbitrariness, it is interesting to realize that the alternative definition FBC = (4 Ac(t))/(π whead τ2)Mw(t) in which τ ≡ tengulf/√ kopen brings up another time scale of interest. Even though the resulting force is obviously identical to that of equations 16 and 17, it suggests another fundamental time scale relevant to engulfment which, as shown in Figure 7, is more insensitive to body size than engulfment time (tengulf) in addition of being three times shorter (∼ 1 to 2 seconds). Whether this time scale τ reflects a property of VGB musculature rather than a constraint on cavity filling (and initial speed Vc(0)) isn’t entirely clear at this time.
Shape drag (FSD).
The drag associated with the flows deflected around an (inflating) body can always be calculated from the general form FSD = ½ ρ C(t)S(t)Vc2(t) in which the values of the time-changing drag coefficient C(t) are chosen to yield the best match with tag data (Potvin et al. ). The time-dependent drag area C(t)S(t) represents the combined effects of the accelerating (and decelerating) external flows created around the whale (i.e., the so-called added or apparent mass), as well as the effects of wake growth and wake turbulence that are generated during the lunge . Although conceptually correct, this approach introduces large uncertainties from the lack of knowledge on the function C(t). Here an alternative but equivalent formulation of shape drag - the so-called apparent mass formulation - is used to reduce such uncertainty by including unsteady flow effects that are known to occur, namely those of the deflected fluid masses co-decelerating with the whale:(18)
The drag area CDnwS(t) is defined from a known reference surface area S(t) and from a drag coefficient CDnw which is typically not equal to the steady-velocity value . This term accounts for the turbulent near-wake found right behind the body. The other two terms in dVc/dt and in Vc describe the effects of the deflected flows moving about the near wake and also co-decelerating with the whale. Parameter Vam(t) is a reference volume for the fluid accelerated around the body. Such volume may change over time as with the cases of opening parachutes . On the other hand, kam parameterizes the fraction of added mass that co-accelerates with the body. Since not all the surrounding fluid is doing so at the rate of dVc/dt, kam measures (roughly) the equivalent amount of fluid that would be doing so in order to achieve the same drag. With decelerating bodies at high Reynolds numbers, kam is a complicated function of time which may be strongly dependent on the motion’s history prior to the deceleration –. This contrasts with the low Reynolds number regime in which kam is a constant and dependent only on the shape of the body . In most cases kam must be determined empirically and verified as time-independent for the system at hand . The added mass formulation has been used previously in comparative bio-mechanics, most notably with the study of shrimp and medusan swimming –.
The BLF makes use of the kam -values associated with ellipsoids, dirigibles and blimps , , , which for the typical rorquals body aspect ratio of 1 to 3 (i.e., aspect ratio defined as buccal cavity length over width (filled)), would yield kam ∼ 0.2. Here the reference volume of the added mass (Vam) is approximated by a half-ellipsoid of revolution defined by semi-minor axes L0/2 and ½(whead/2+ Ljaw).
Note that the added mass formulation embodied by equation 18 is not unique, as both terms in Vc end Vc2 are derived from a deflection rate of the incoming fluid mass . Depending on the specific accounting of what constitutes the near wake (for the Vc2 term) and what constitute the co-decelerated flows off-wake (the kam Vc term), different models, and hence different values of kam, Vam CDnw and S(t), can be used to yield the same force. For example, the kam Vc -term can be shown to behave as a Vc2–term when the cavity is inflating mostly transversely at rates that are proportional to Vc:(19)
(rtransv is the mouth radius and Ξ a proportionality constant).
The added mass drag model embodied in the BLF takes the following forms:
The kam dVc/dt term is neglected in equation 20 given the low body (whale) accelerations and engulfed volumes at play during mouth opening [J. Potvin; unpublished data]. In contrast, this term is quantitatively important during mouth closure and, with dVc/dt being negative (a deceleration), approximately reflects the effects of the near wake re-contact the body by actually reducing shape drag. Parameter Abody is the known maximum cross-sectional area of the body in a closed-mouth, empty-cavity configuration (Figure 1), Ac(t) the instantaneous mouth area (equation 11) and Acmax the maximum mouth area (equation 11 with θgapegape = 50° or 78° ). CDbody corresponds to the drag coefficient for a mouth-closed and empty-body configuration, set to CDbody = 0.05 from a previous hydrodynamic study on fin whale locomotion . Although the drag coefficients CDopen and CDclose were assigned the values shown in Tables 1, 2, 3, 4 (∼ 0.3 to 0.5), the resulting accelerations and speeds turned out to be insensitive to their specific values, in contrast to the simulations discussed in , . This is particularly true during mouth closure when shape drag is exactly cancelled by fluking thrust as discussed next.
Swimming thrust (Fthrust), buoyancy-reduced weight (Fext), and “ocean-to-engulfed mass” drag force (Fww).
The remaining forces acting on the body and engulfed mass are as follows (Figure 2): swimming thrust (Fthrust); buoyancy-reduced weight (Fext; component along the track); and the “ocean-to-engulfed mass” drag force (Fww), which acts on the ocean-facing end of the engulfed mass. Calculating their values during engulfment can only be estimated as they have never been experimentally characterized. The buoyancy-reduced weight is defined as Fext = (FB – Wc) sinΘ, with FB, Wc and Θ corresponding to a whale’s buoyancy, weight and angle of the upward-inclined trajectory respectively (Note: at depth FB < Wc, a case of negative buoyancy). On the other hand, Fww arises from the pressure applied by the oncoming flow under the palate to the ocean-facing end of the engulfed water (Figure 2). In a previous version of the BLF both Fthrust and Fext were estimated from tag data as being quite small relative to engulfment drag (FED) and set to zero values (see Appendix 1 in Goldbogen et al ); similarly, Fww was set to zero. In this paper we use non-zero values which are still small but yield better matches with the (speed) tag data of the three species, in addition to yielding a more realistic picture of the overall fluid motion.
Mouth Opening Stage
The combination Fthrust + Fext is modeled simply by using a constant value throughout that stage, presumably the average of the actual (but unknown) time-varying force. For the size corresponding to the average body length of humpback whales (14.4 m), fin whales (20.2 m), and blue whales (25.2 m), Fthrust + Fext is set to a value that yields the best match with the speed data of tag studies , . Such values (Tables 1, 2, 3, 4) are then scaled all other body dimensions as follows:
These scaling laws are based on assuming Fthrust being much greater than Fext, and by using a simple model in which Fthrust is given by fluke area times flow speed relative to the fluke -squared. As discussed in more details in the very last section of this paper, Fthrust scales mostly with fluke area since the flow speed relative to it is driven mostly by the tangential speed of the fluke along the arc of tail sweep, a variable which turns out to scale weakly with Lbody. Note that the Fthrust-values obtained from equation 21 are generally lower than those obtained from the Bose-Lien model used in the RAAMR calculation (at the same speed and body size). This should follow from the intermittent short bouts of fluking that occurs during mouth opening. On the other hand, the speed of the reflux being generally low during most of the mouth opening stage means that little pressure is felt by the ocean-facing end of the engulfed mass (residing just below the TMJ). Thus the ocean-to-engulfed mass drag force is set to.(25)
Mouth Closure Stage
The equivelocity condition required by SE, as applied to the water in both cavity sections, makes the interaction between whale and engulfed water look like a perfectly inelastic collision. In elementary college physics, such a collision is exemplified by two balls moving towards each other in the absence of external forces, then colliding and moving together while permanently attached. Here, and in the presence of external forces, a perfectly inelastic collision is implemented to yield Vc = Vw and ac = aw, or in other words, equivelocity and equideceleration, by solving equations 14 and 15 under the following constraints:(26)(27)
In this scheme the forces (Fthrust +Fext) and Fww are calculated from the knowledge of the mass so-far engulfed (equation 13), the reflux speed (equation 15), the filling rate (equation 13), and both shape and engulfment drag (equations 16–17 and 20–21).
Although these constraints are not unique in yielding equideceleration, they represent a plausible scenario in which the whale applies swimming thrust in amounts that are just enough to keep up with the deceleration of its engulfed cargo, which itself is caused by the push differential between buccal cavity force (FBC) and the back-end pressure from the ocean (Fww). In this case thrust is modulated to counteract the decelerating effects of shape drag and negative buoyancy. Interestingly, the body size scaling of swimming thrust now reflects that of shape drag as with the case of steady travel. Moreover, even though any shape drag model becomes irrelevant to the actual motion of the whale (being cancelled by Fthrust + Fext), they remain quite relevant to metabolic expenditures as they become a driver of swimming thrust (during mouth closure).
More results; scaling of engulfment time and peak cavity wall force.
Typical BLF results are shown in Figures 5, 6, 7 previously discussed. More are displayed in Figures 13 and 14, which show the temporal changes in the position of the lead slug of water (or first-engulfed slug) and of the mass so-far engulfed. In Figure 13 the lead slug’s motion clearly ends prior to maximum gape (tengulf/2) in accordance with SE. Such motion stops a meter or so away from the aft-end of the cavity located at a distance L0– Ljaw posterior to the TMJ, a feature designed for avoiding cavity drainage during the rest of engulfment. The extent of the partial filling post-TMJ is also shown in Figure 14. Note that the modeling allows the cavity ant-TMJ to fill completely in contrast to the cavity post-TMJ.z
Simulated positions of the first-engulfed slug for a 20.2 m fin whale. Here the lead slug stops traveling beyond the time of maximum gape (i.e. after the 2.85 s mark), per Synchronized Engulfment. Note also that this slug can never travel past the umbilicus, or in other words, cover a distance exceeding L0– Ljaw ( = 7.34 m in this case).
The cavity post-TMJ fills first during the first 66% of the mouth-opening (i.e., until 1.9 s here), to remain filled thereafter. This is followed by the filling of the cavity ant-TMJ during mouth closure. These are compared with the total mass capacity of the cavity post-TMJ (equation 9; dash-dotted line) and total capacity (sum of equations 9 and 10; dash-doubly-dotted). The body mass of this case is 40,705 kg.
The simulations show how engulfment forces and times change with body size (Figures 6, 7 and 15, and Tables 5, 7 and 8), as well as how they are constrained by: 1) the scaling of Synchronized Engulfment, 2) the allometry of the engulfment apparatus , and 3) the initial lunge speed (Vc(0) = Vn Lbody ). Interestingly, these three elements contribute in different ways for each engulfment performance measures.
Temporal development at various body lengths: 24.0 m (black, continuous), 22.7 m (black, dashed), 20.2 m (black, dashed-doubly-dotted), 18.5 m (black, dashed-dotted), 17.7 m (black, dotted), 13.5 m (blue, continuous), 10.0 m (blue, dashed).
Engulfment time reflects the opening and closing (angular) speeds of the mouth as driven by the escape speeds of the prey . Elementary dimensional analysis suggests engulfment time to scale as a ratio of length over speed , , or in the context of engulfment, as tengulf ∼ Ljaw/Vc(0). A recent analysis of tag data for all three species of large rorquals , as well as plausible prey escape scenarios , suggest the initial velocity to scale as Vc(0) ∼ 0.15 Lbody which in turns implies tengulf ∼ Ljaw/Lbody. Using the morphological data of Tables 1, 2, 3, 4 where Ljaw ∼ Lbody1.21, 1.29 and 1.41 for humpback, fin and blue whales, respectively, engulfment time would ultimately scale as tengulf ∼ Lbody0.21, 0.29 and 0.41. Interestingly, this is quite similar to the scaling of SE-constrained engulfment time of equation 8 for which tengulf ∼ Lbody0.20, 0.27 and 0.41 (Table 2; see also Figure 7). This comparison shows that even though SE contributes a non-trivial dynamical scaling factor to engulfment time (namely, the factor Γ(L0– Ljaw)/Ljaw in equation 8), the basic allometry of the skull in which L0 ∼ Lbody1.18, 1.16 and 1.19 (Tables 1, 2, 3, 4) significantly weakens its significance to levels usually expected with isometric body dimensions. This is an important “non-contribution” given the body size range examined here, from weaned juvenile humpback whales (8 m) up to the largest adult blue whales (27 m).
Mass-specific peak cavity push (or peak engulfment drag) at maximum gape is a measure for which SE is somewhat more important. As shown in Figure 6, and per Equation 11, the increase in mass-specific cavity wall force is mostly the result of the positive allometry of the skull which here directly increases mandibular size (Ljaw) as well as engulfment capacity (Mwpost-TMJ). With tengulf scaling as previously discussed, and with Mwpost-TMJ, Mc and Ac/whead scaling as Lbody3.43, 3.66 and 3.80 (equation 9), Lbody3.00, 2.74 and 3.46 (Tables 1, 2, 3, 4) and Lbody1.21, 1.29 and 1.41 (i.e., as ∼Ljaw; see equation 11) respectively, the product (Ac/whead)(Mwpost-TMJ/Mc) thus varies as Lbody1.23, 1.67 and 0.93 and so drives the scaling of the ratio as FBCopen|peak/Mc as ∼ Lbody1.14, 1.09 and 0.75 (Table 5). The difference resides in the scaling of the reaction constant kopen appearing in equation 16, which is tuned at each body size to enforce equivelocity at maximum gape. The relatively weaker scaling of kopen (∼ Lbody−0.36, −0.61, −0.35; Tables 1, 2, 3, 4) shows again the reduced role of SE in comparison to morphology. This is also seen with the mass-specific expended power, where the SE-dependent effects connected to the scaling of kopen are neutralized by those of the (average) speed (also SE-dependent).
Scaling of fluking thrust during engulfment.
The scaling of the fluke thrust (Fthrust) embodied in equations 20 and 21, although seemingly simple, is by no means straightforward. It is usually derived in the context of non-feeding steady travel (at speed Vcruise) during which body drag is equal to fluking thrust (in this context “body drag” is synonymous with “shape drag”). Thus Fthrust ∼ Lbody2 Vcruise2 with an assumption of the drag coefficient being weakly dependent on body size . In cases where the cruise speed is determined by minimal energy expenditure via the use of a resonance related to tail stiffness, one has Vcruise ∼ Lbody1/2 and Fthrust ∼ Lbody3 . On the other hand, where there is a need for the differently-sized members of a group travelling together to move at the same speed, fluking thrust is applied differently by each member so to yield Vcruise ∼ Lbody0 and Fthrust ∼ Lbody2. Using body drag as basis for studying the scaling of fluking thrust no longer applies for engulfment given that the associated motions are fundamentally unsteady, or in other words, given that thrust is equal to body drag in addition to the inertial force (i.e, -Mcac) and engulfment drag (FED).
A much simpler – and fundamental - way to look at thrust is to view the flukes as (flexible) hydrofoils that produce lift and drag forces combining into thrust in the forward direction (Figure 16):(28)
Here the angles Ω and ω measure the directions of the tail and fluke’s tangential speed (Vtan) respectively. CR is the force coefficient resulting from the (vector) addition of both airfoil drag and lift, and Δ the angle made by the resultant with respect to the horizontal. Finally, Sfluke is the reference area defined as fluke span times fluke chord.
The flukes’ tangential speed along the arc of a sweep is estimated as Vtan ∼. Ltail Ωsweep/tsweep, with Ltail ∼ Lbody/2 and Ωsweep = 1.57rads (i.e., 90°). The sweep duration tsweep is informed by the trends seen in tag data collected on the three species , , , , which suggests that one fluking sweep early in the mouth opening stage (either upwards or downwards) lasts about 1/3 of engulfment time (or tsweep ∼ 1/3 tengulf ). In the case of the average body sizes shown in Tables 1, 2, 3, 4, one has tengulf = 4.1s (humpback), 5.69 s (fin) and 6.42s (blue), yielding tsweep = 1.36s (humpback), 1.89s (fin) and 2.33s (blue). With Ltail being equal to 7 m (humpback), 10 m (fin) and 13 m (blue), one has Vtan ∼ 8.1 m/s (humpback), 8.3 m/s (fin) and 8.76 m/s (blue). This tangential velocity turns out to be greater than the whale’s translation speed (Vc) thus leading to Vtan2>> Vc2 and to the neglect of the Vc2–term in equation 26. This result means that the velocity-squared factor of the equation is roughly independent of body size. After again assuming the aerodynamic force coefficient CR as independent of body size, fluke thrust thus scales as Fthrust ∼ Sfluke Vtan2.
Understanding why Vtan scales weakly with body size can be achieved by realizing that Vtan ∼ Lbody/tsweep ∼ Lbody/Lbody ∼ 1. The last step follows from an argument similar to Hill’s estimate of muscle contraction frequency under maximum stress, as discussed by Pennycuick . Thus Fthrust ∼ Sfluke, and although Fthrust ∼ Lbody2 in an isometric world, fluking thrust scales somewhat more weakly in the real world of rorquals . Indeed, scaling the fluke reference area in terms of the fluke span (Kspan) and fluke width at insertion (i.e., KWAI, measured from notch of fluke to medial insertion ) one obtains for both fin and blue whales Sfluke ∼ KspanKWAI ∼ Lbody1.07 Lbody0.84 , and Sfluke ∼ KspanKWAI ∼ (Lbody –0.77) Lbody0.82 for humpback whales, which then lead to the scaling laws of Eqs. 20–22 (for Kspan see ; scaling for KWAI is from J. Goldbogen (unpub. data)).
J.P. thanks C. Scott for useful discussions on the metabolic expenditures during strenuous exercise by humans, and D. Peters for references on the added mass of blimps. The authors thank A. Werth for pointing out the importance of minke whale simulations in the interpretation of low EMR at small body size.
Conceived and designed the experiments: JP JAG RES. Performed the experiments: JP. Analyzed the data: JP JAG RES. Contributed reagents/materials/analysis tools: JP JAG RES. Wrote the paper: JP JAG RES.
- 1. Williams TM (2006) Physiological and ecological consequences of extreme body size in whales. In: Estes JA, DeMaster DP, Doak DF, Williams TM, Brownell RL, editors. Whales, Whaling, and Ocean Ecosystems. Berkeley, CA: University of California Press. 191–201.
- 2. Werth AJ (2000) Feeding in marine mammals. In: Schwenk K, editor. Feeding: Form, Function and Evolution in Tetrapod Vertebrates. New York, NY: Academic Press. 475–514.
- 3. Werth AJ (2001) How do mysticetes remove prey trapped in baleen? Bull Mus Comp Zool 156: 189–203.
- 4. Werth AJ (2007) Adaptations of the cetacean hyolingual apparatus for aquatic feeding and thermoregulation. Anatomical Record-Advances in Integrative Anatomy and Evolutionary Biology 290: 546–568.
- 5. Simon M, Johnson M, Tyack P, Madsen PT (2009) Behaviour and kinematics of continuous ram filtration in bowhead whales (Balaena mysticetus). Proceedings of the Royal Society B-Biological Sciences.
- 6. Werth AJ (2004) Models of hydrodynamic flow in the bowhead whale filter feeding apparatus. Journal of Experimental Biology 207: 3569–3580.
- 7. Lambertsen RH, Rasmussen KJ, Lancaster WC, Hintz RJ (2005) Functional morphology of the mouth of the bowhead whale and its implications for conservation. Journal of Mammalogy 86: 342–352.
- 8. Lambertsen RH (1983) Internal mechanism of rorqual feeding. Journal of Mammalogy 64: 76–88.
- 9. Goldbogen JA (2010) The Ultimate Mouthful: Lunge Feeding in Rorqual Whales. American Scientist 98: 124–131.
- 10. Friedlaender AS, Hazen EL, Nowacek DP, Halpin PN, Ware C, et al. (2009) Diel changes in humpback whale Megaptera novaeangliae feeding behavior in response to sand lance Ammodytes spp. behavior and distribution. Marine Ecology-Progress Series 395: 91–100.
- 11. Ware C, Friedlaender AS, Nowacek DP (2011) Shallow and deep lunge feeding of humpback whales in fjords of the West Antarctic Peninsula. Marine Mammal Science 27: 587–605.
- 12. Wiley D, Ware C, Bocconcelli A, Cholewiak D, Friedlaender A, et al. (2011) Underwater components of humpback whale bubble-net feeding behaviour. Behaviour 148: 575–602.
- 13. Friedlaender AS, Lawson GL, Halpin PN (2009) Evidence of resource partitioning between humpback and minke whales around the western Antarctic Peninsula. Marine Mammal Science 25: 402–415.
- 14. Woodward BL, Winn JP, Fish FE (2006) Morphological specializations of baleen whales associated with hydrodynamic performance and ecological niche. Journal of Morphology 267: 1284–1294.
- 15. Nemoto T (1970) Feeding pattern of baleen whales in the ocean. In: Steele JH, editor. Marine Food Chains. Berkeley: University of California Press. 241–252.
- 16. Pyenson ND, Goldbogen JA, Vogl AW, Szathmary G, Drake RL, et al. (2012) Discovery of a sensory organ that coordinates lunge feeding in rorqual whales. Nature 485: 498–501.
- 17. Brodie PF (1993) Noise Generated by the Jaw Actions of Feeding Fin Whales. Canadian Journal of Zoology-Revue Canadienne De Zoologie 71: 2546–2550.
- 18. Brodie PF (2001) Feeding mechanics of rorquals Balaenoptera sp. In: Mazin JM, de Buffrenil V, editors. Secondary Adaptations of Tetrapods to Life in Water. Munchen, Germany: Verlag. 345–352.
- 19. Orton LS, Brodie PF (1987) Engulfing Mechanics of Fin Whales. Canadian Journal of Zoology 65: 2898–2907.
- 20. Arnold PW, Birtles RA, Sobtzick S, Matthews M, Dunstan A (2005) Gulping behaviour in rorqual whales: underwater observations and functional interpretation. Mem Queensl Mus.
- 21. Lambertsen R, Ulrich N, Straley J (1995) Frontomandibular Stay of Balaenopteridae - a Mechanism for Momentum Recapture During Feeding. Journal of Mammalogy 76: 877–899.
- 22. Goldbogen JA, Pyenson ND, Shadwick RE (2007) Big gulps require high drag for fin whale lunge feeding. Marine Ecology-Progress Series 349: 289–301.
- 23. Goldbogen JA, Potvin J, Shadwick RE (2010) Skull and buccal cavity allometry increase mass-specific engulfment capacity in fin whales. Proceedings of the Royal Society B-Biological Sciences 277: 861–868.
- 24. Goldbogen JA, Calambokidis J, Croll D, McKenna MF, Potvin J, et al. (2012) Scaling of lunge feeding performance in rorqual whales: mass-specific energy expenditure increases with body size and progressively limits diving capacity. Functional Ecology 26: 216–226.
- 25. Goldbogen JA, Calambokidis J, Oleson E, Potvin J, Pyenson ND, et al. (2011) Mechanics, hydrodynamics and energetics of blue whale lunge feeding: efficiency dependence on krill density. Journal of Experimental Biology 214: 131–146.
- 26. Potvin J, Goldbogen JA, Shadwick RE (2009) Passive versus active engulfment: Verdict from trajectory simulations of lunge-feeding fin whales Balaenoptera physalus. Journal of the Royal Society Interface 6: 1005–1025.
- 27. Potvin J, Goldbogen JA, Shadwick RE (2010) Scaling of lunge feeding in rorqual whales: An integrated model of engulfment duration. Journal of Theoretical Biology 267: 437–453.
- 28. Croll DA, Acevedo-Gutierrez A, Tershy BR, Urban-Ramirez J (2001) The diving behavior of blue and fin whales: is dive duration shorter than expected based on oxygen stores? Comparative Biochemistry and Physiology a-Molecular & Integrative Physiology 129: 797–809.
- 29. Shaffer SA, Costa DP, Williams TM, Ridgway SH (1997) Diving and swimming performance of white whales, Delphinapterus leucas: An assessment of plasma lactate and blood gas levels and respiratory rates. Journal of Experimental Biology 200: 3091–3099.
- 30. Wiedenmann J, Cresswell KA, Goldbogen J, Potvin J, Mangel M (2011) Exploring the effects of reductions in krill biomass in the Southern Ocean on blue whales using a state-dependent foraging model. Ecological Modelling 222: 3366–3379.
- 31. Hill RW, Wyse GA, Anderson M (2008) Animal Physiology: Sinauer Associates Inc.
- 32. Van Wassenbergh S, Aerts P (2009) Aquatic suction feeding dynamics: insights from computational modelling. Journal of the Royal Society Interface 6: 149–158.
- 33. Parr WCH, Wroe S, Chamoli U, Richards HS, McCurry MR, et al. (2012) Toward integration of geometric morphometrics and computational biomechanics: New methods for 3D virtual reconstruction and quantitative analysis of Finite Element Models. Journal of Theoretical Biology 301: 1–14.
- 34. Goldbogen JA, Calambokidis J, Croll D, Harvey J, Newton K, et al. (2008) Foraging behavior of humpback whales: kinematic and respiratory patterns suggest a high cost for a lunge. Journal of Experimental Biology 211: 3712–3719.
- 35. Goldbogen JA, Calambokidis J, Shadwick RE, Oleson EM, McDonald MA, et al. (2006) Kinematics of foraging dives and lunge-feeding in fin whales. Journal of Experimental Biology 209: 1231–1244.
- 36. Goldbogen JA, Potvin J, Shadwick RE (2011) Passive elasticity of ventral groove blubber does not limit engulfment volume in lunge-feeding whales. Integrative and Comparative Biology 51: E193–E193.
- 37. Bose N, Lien J (1989) Propulsion of a fin whale (Balaenoptera physalus) - why the fin whale is a fast swimmer. Proceedings of the Royal Society of London Series B-Biological Sciences 237: 175–200.
- 38. Fish FE, Rohr JJ (1999) Review of dolphin hydrodynamics and swimming performance. San Diego, CA.
- 39. Lockyer C (1976) Body Weights of Some Species of Large Whales. ICES Journal of Marine Science 36: 259–273.
- 40. Kermack KA (1948) The propulsive powers of blue and fin whales. Journal of Experimental Biology 25: 237–240.
- 41. Costa DP, Williams TM (1999) Marine Mammal Energetics. In: Reynolds JE, Rommel SA, editors. Biology of Marine Mammals. Washington, D.C.: Smithonsian Institution Press. 176–217.
- 42. Croll D, Kudela R, Tershy BR (2006) Ecosystem impact of the decline of large whales in the North Pacific. In: Estes JAea, editor. Whales, whaling and ocean ecosystems. Berkeley, CA: University of California Press. 202–214.
- 43. Williams TM (1999) The evolution of cost efficient swimming in marine mammals: limits to energetic optimization. Philosophical Transactions of the Royal Society of London Series B-Biological Sciences 354: 193–201.
- 44. Kleiber M (1975) The Fire of Life: an Introduction to Animal Energetics. Huntington, NY: R.E. Kteiger Publishing Co.
- 45. Kleiber M, Rogers TA (1961) Energy metabolism. Annu Rev Physiol 23: 5–36.
- 46. Sigurjonsson J, Vikingsson GA (1997) Seasonal abundance of and Estimated Food Consumption by Cetaceans in Icealandic And Adjacent Waters. J Northw Atl Fish Sci 22: 271–287.
- 47. Lockyer CH (1981) Growth and Energy Budgets of Large Baleen Whales from the Southern Hemisphere. Mammals in the Seas: FAO Fish. 379–487.
- 48. Ahlborn BK (2004) Zoological Physics. Berlin: Springer-Verlag.
- 49. Williams TM (1989) Swimming by sea otters: adaptations for low energetic cost locomotion. Journal of Comparative Physiology A: Sensory Neural and Behavioral Physiology 164: 815–824.
- 50. Watkins WA, Sigurjonsson J, Wartzok D, Maiefski RR, Howey PW, et al. (1996) Fin whale tracked by satellite off Iceland. Marine Mammal Science 12: 564–569.
- 51. Bailey H, Mate B, Irvine L, Palacios DM, Bograd SJ, et al. (2009) Blue whale behavior in the eastern North Pacific inferred from state-space model analysis of satellite tracks. Endangered Species Research 10: 93–106.
- 52. Watanabe YY, Sato K, Watanuki Y, Takahashi A, Mitani Y, et al. (2011) Scaling of swim speed in breath-hold divers. Journal of Animal Ecology 80: 57–68.
- 53. Sato K, Watanuki Y, Takahashi A, Miller PJO, Tanaka H, et al. (2007) Stroke frequency, but not swimming speed, is related to body size in free-ranging seabirds, pinnipeds and cetaceans. Proceedings of the Royal Society B-Biological Sciences 274: 471–477.
- 54. Bohonak AJ, van der Linde K (2004) RMA: Software for Reduced Major Axis regression for Java.
- 55. Víkingsson GA, Sigurjónsson J, Gunnlaugsson T (1988) On the relationship between weight, length and girth dimensions in fin and sei whales caught off Iceland. Report of the International Whaling Commission 38: 323–326.
- 56. Quiring DP (1943) Weight data on five whales. Journal of Mammalogy 24: 39–45.
- 57. Nishiwaki M, Oye T (1951) Biological investigation on blue whales (Balaenoptera musculus) and fin whales (Balaenoptera physalus) caught by the Japanese Antarctic whaling fleets. Sci Rep Whales Res Inst 5: 91–167.
- 58. Mackintosh NA (1942) The southern stocks of whalebone whales. Discovery Reports 22: 197–300.
- 59. Lockyer C, Waters T (1986) Weights and anatomical measurements of Northeastern Atlantic fin (Balaneoptera physalus, Linnaeus) and sei (B. borealis, Lesson) whales. Marine Mammal Science 2: 169–185.
- 60. Ash CE (1953) Weights of Antarctic humpback whales. Norsk Hvalfangsttid 42: 387–391.
- 61. Mackintosh NA, Wheeler JFG (1929) Southern blue and fin whales. Discovery Reports 1: 257–540.
- 62. Matthews LH (1937) The humpback whale, Megaptera nodosa. Discovery Reports 22: 197–300.
- 63. Blix AS, Folkow LP (1995) Daily Energy-Expenditure in Free-Living Minke Whales. Acta Physiologica Scandinavica 153: 61–66.
- 64. Alexander RM (2005) Models and the scaling of energy costs for locomotion. Journal of Experimental Biology 208: 1645–1652.
- 65. Domenici P (2001) The scaling of locomotor performance in predator-prey encounters: from fish to killer whales. Comparative Biochemistry and Physiology a-Molecular and Integrative Physiology 131: 169–182.
- 66. Biewener AA (2005) Biomechanical consequences of scaling. Journal of Experimental Biology 208: 1665–1676.
- 67. Alexander RM (1998) All-time giants: The largest animals and their problems. Palaeontology 41: 1231–1245.
- 68. Alexander RM (2003) Principles of Animal Locomotion. Princeton, NJ: Princeton University Press. 369 p.
- 69. Schmidt-Nielsen K (1984) Scaling: why is animal size so important ? Cambridge: Cambridge University Press.
- 70. Vogel S (2008) Modes and scaling in aquatic locomotion. Integrative and Comparative Biology 48: 702–712.
- 71. Higham TE (2007) The integration of locomotion and prey capture in vertebrates: Morphology, behavior, and performance. Integrative and Comparative Biology 47: 82–95.
- 72. Williams TM, Yeates L (2004) The energetics of foraging in large mammals: a comparison of marine and terrestrial predators. International Congress Series 1275: 351–358.
- 73. Lockyer C (2007) All creatures great and smaller: a study in cetacean life history energetics. Journal of the Marine Biological Association of the United Kingdom 87: 1035–1045.
- 74. Boyd I (2002) Marine mammal energetics. In: Hoelzel A, editor. Marine Mammal Biology: An Evolutionary Approach. Oxford, UK: Blackwell. 247–277.
- 75. Weibel ER, Hoppeler H (2005) Exercise-induced maximal metabolic rate scales with muscle aerobic capacity. Journal of Experimental Biology 208: 1635–1644.
- 76. Kanatous SB, Davis RW, Watson R, Polasek L, Williams TM, et al. (2002) Aerobic capacities in the skeletal muscles of Weddell seals: key to longer dive durations? Journal of Experimental Biology 205: 3601–3608.
- 77. Williams TM, Friedl WA, Haun JE (1993) The physiology of bottlenose dolphins (Tursiops truncatus): heart rate, metabolic rate, and plasma lactate concentration during exercise. Journal of Experimental Biology 179: 31–46.
- 78. McNab BK (2009) Resources and energetics determined dinosaur maximal size. Proceedings of the National Academy of Sciences of the United States of America 106: 12184–12188.
- 79. Pyenson ND, Sponberg SN (2011) Reconstructing Body Size in Extinct Crown Cetacea (Neoceti) Using Allometry, Phylogenetic Methods and Tests from the Fossil Record. Journal of Mammalian Evolution 18: 269–288.
- 80. Demere TA, McGowen MR, Berta A, Gatesy J (2008) Morphological and molecular evidence for a stepwise evolutionary transition from teeth to baleen in mysticete whales. Systematic Biology 57: 15–37.
- 81. Nowacek DP, Friedlaender AS, Halpin PN, Hazen EL, Johnston DW, et al.. (2011) Super-Aggregations of Krill and Humpback Whales in Wilhelmina Bay, Antarctic Peninsula. Plos One 6.
- 82. Pivorunas A (1977) Fibro-Cartilage Skeleton and Related Structures of Ventral Pouch of Balaenopterid Whales. Journal of Morphology 151: 299–313.
- 83. Schulte HVW (1916) Anatomy of a foetus. Balaenoptera borealis. Mem Am Mus Nat Hist 1: 389–502.
- 84. Perrine D (2009) In Hot Pursuit: Bryde's whales rocket through Pacific shallows to gorge on fish. National Geographic 216: 65, 70–71.
- 85. Macha JM (1993) A Simple Approximate Model of Parachute Inflation; 1993; London, United Kingdom.
- 86. Cochran BC, White BR, Macha JM (1991) Experimental Investigation of Added Mass during Parachute Deceleration – Preliminary Results; April 9–11; San Diego, CA. 171–180.
- 87. Sarpkaya T, Isaacson M (1981) Mechanics of Wave Forces on Offshore Structures. New York: Van Nostrand Reinhold Company.
- 88. Strickland JH, Macha M (1990) Preliminary Characterization of Parachute Wake Recontact. Journal of Aircraft 27: 501–506.
- 89. Batchelor GK (2000) An Introduction to Fluid Dynamics: Cambridge University Press.
- 90. Daniel TL (1984) Unsteady aspects of aquatic locomotion. American Zoologist 24: 121–134.
- 91. Daniel TL (1985) Cost of locomotion: unsteady medusan swimming. Journal of Experimental Biology 119: 149–164.
- 92. Daniel TL, Meyhofer E (1989) Size limits in escape locomotion of carridean shrimp Journal of Experimental Biology. 143: 245–265.
- 93. Tuckerman LB (1925) Inertia factors of ellipsoids for use in airship design. National Advisory Committee for Aeronautics (NACA) 210.
- 94. Munk MM (1924) Some tables of the factor of apparent mass. National Advisory Committee for Aeronautics (NACA) 197.
- 95. Wilson JF (2003) Dynamics of offshore structures: John Wiley & Sons Inc.
- 96. Ahlborn BK, Blake RW, Chan KHS (2009) Optimal fineness ratio for minimum drag in large whales. Canadian Journal of Zoology-Revue Canadienne De Zoologie 87: 124–131.
- 97. Pennycuick CJ (1992) Newton Rules Biology. Oxford: Oxford University Press.
- 98. Sousa-Lima RS, Groch KR (2010) Correlation between body length and fluke width in humpback whales, Megaptera novaeangliae. Marine Mammal Science 26: 977–981.