ABSTRACT
The scale dependence of locomotor factors has long been studied in comparative biomechanics, but remains poorly understood for animals at the upper extremes of body size. Rorqual baleen whales include the largest animals, but we lack basic kinematic data about their movements and behavior below the ocean surface. Here, we combined morphometrics from aerial drone photogrammetry, whale-borne inertial sensing tag data and hydrodynamic modeling to study the locomotion of five rorqual species. We quantified changes in tail oscillatory frequency and cruising speed for individual whales spanning a threefold variation in body length, corresponding to an order of magnitude variation in estimated body mass. Our results showed that oscillatory frequency decreases with body length (∝length−0.53) while cruising speed remains roughly invariant (∝length0.08) at 2 m s−1. We compared these measured results for oscillatory frequency against simplified models of an oscillating cantilever beam (∝length−1) and an optimized oscillating Strouhal vortex generator (∝length−1). The difference between our length-scaling exponent and the simplified models suggests that animals are often swimming non-optimally in order to feed or perform other routine behaviors. Cruising speed aligned more closely with an estimate of the optimal speed required to minimize the energetic cost of swimming (∝length0.07). Our results are among the first to elucidate the relationships between both oscillatory frequency and cruising speed and body size for free-swimming animals at the largest scale.
INTRODUCTION
For most animals, locomotion and maneuverability factor into critical life functions such as prey capture, predator avoidance and migratory success, and therefore can influence the margin between life and death (Dakin et al., 2018; Fish et al., 2003; Goldbogen et al., 2012; Hein et al., 2012; Williams et al., 2014). These locomotor functions manifest as complex maneuvers or behaviors that can be deconstructed into a series of simple elemental changes in body conformation, rotation or displacement (Full et al., 2002; Pennycuick, 1975; Segre et al., 2016, 2018; Webb, 1997). Many of these elements have evolved similarly among different taxa with distinct body plans as a result of common physical constraints or functional requirements (Donley et al., 2004). As a result, animals have evolved a wide range of convergent locomotor and maneuvering strategies across varying body sizes (Gleiss et al., 2011, 2017). Nowhere are the effects of scale more interesting than in the ocean, where release from gravitational constraints allows species with similar body plans and locomotory styles to range across several orders of magnitude in body mass, from the smallest fishes up to the blue whale (Balaenoptera musculus) (Motani, 2002; Pennycuick, 1992). Both toothed whales (Odontoceti) and baleen whales (Mysticeti) evolved gigantism after the evolutionary appearance of specialized foraging strategies (i.e. echolocation and filter feeding) that required the coupling of locomotion and feeding (Goldbogen and Madsen, 2018; Higham et al., 2016; Pyenson, 2017; Slater et al., 2017). However, what remains poorly understood is the influence of gigantism on the locomotor performance of whales and the ways in which it has enabled their foraging strategies and overall success as predators in a wide range of ocean ecosystems.
The body plans for cetaceans (whales and dolphins) typically consist of a fusiform body shape with appendages modified into streamlined flippers, fins and tails that oscillate or undulate to generate thrust, adjust attitude and help overcome drag (Fish and Lauder, 2006; Webb and De Buffrénil, 1990). At the largest scales, many large aquatic vertebrates (sharks, tuna, cetaceans) have evolved high aspect ratio lunate tails that are specialized for efficient locomotion and maneuvering (Bose and Lien, 1989; Dewar and Graham, 1994; Donley et al., 2004; Fish, 1998; Fish et al., 2014; Graham et al., 1990; Webb and De Buffrénil, 1990; Woodward et al., 2006). Paired, anteriorly placed flippers or fins are generally used as control surfaces to stabilize body trim or generate asymmetrical forces during turns, rolls or pitching motions (Fish and Lauder, 2017; Fish and Shannahan, 2000; Fish et al., 2006; Segre et al., 2016, 2018). Cetaceans (whales, dolphins and porpoises) include the largest animals that have ever lived (McClain et al., 2015), making them a unique and interesting organismal model to study the extremes of swimming performance.
Previous research on cetaceans has focused primarily on small, highly maneuverable toothed whales that can be studied in captivity (Fish, 1993, 1998; Fish and Rohr, 1999; Fish et al., 2014, 2018; Rohr and Fish, 2004). In contrast, the logistical difficulties of working with larger baleen whales has made it more difficult to study the biomechanics and locomotion of these animals in the wild. As a result, much of the progress in this area has come about recently with the advent of archival, motion-sensing tags, which can be attached to wild animals for hours to days (Gleiss et al., 2011; Goldbogen et al., 2017; Martín López et al., 2016; Sato et al., 2007).
In considering the effects of scale on locomotion, previous studies have predicted that the oscillatory frequency of propulsive structures decreases in proportion to increasing body length (Hill, 1950; Pennycuick, 1992; Sato et al., 2007). Using simple muscle models, Hill (1950) predicted that the contraction speed of geometrically similar but differently sized muscles should scale inversely with their physical dimensions in order to balance the production of useful work and wasteful kinetic energy at a given movement speed. Accordingly, steady-state cruising by most vertebrates has been shown to occur in the 1–2 m s−1 range regardless of body size (Fish and Rohr, 1999; Sato et al., 2007; Watanabe et al., 2011). In contrast, dependence on body size has been shown to occur in other swimming performance measures such as maximum speed (Meyer-Vernet and Rospars, 2016; Shadwick and Gemballa, 2005; Wardle, 1975). Thus, the choice of movement speed for routine swimming may not reflect maximally achievable performance but instead relate to economy of transport (or the minimum cost of transport, COTmin) (Williams and Noren, 2009; Williams et al., 1993). However, recent research on free-ranging largemouth bass showed that swimmers exhibited slower speeds than their predicted COTmin, perhaps in order to balance competing functional demands such as prey capture (Han et al., 2017), suggesting that swimming performance is highly context dependent.
The hypothesis that stroke frequency is inversely correlated with body length is also supported by analyses of the dimensionless Strouhal number St, a ratio of the oscillatory frequency f and tailbeat amplitude A over swimming speed U [St=(f·A)/U]. High propulsive efficiency occurs at St≈0.25–0.35 (Anderson et al., 1998; Rohr and Fish, 2004; Taylor et al., 2003; Triantafyllou et al., 1991). If swimming speed is constant and tailbeat amplitude scales isometrically with body length, then the oscillatory frequency must scale inversely with body length (∼L−1.0) in order to achieve optimal vortex generation (Fish, 1998; Pennycuick, 1992).
These size-dependent factors have played a major role in shaping the functional ecological theory of gigantic marine predators, but, until now, we have lacked the data for a comprehensive study on the kinematics of their oscillatory swimming. Specifically, it remains to be seen how swimming speed and oscillatory frequency scale with body length, particularly at the largest scale. To answer this question, we used tags with inertial sensors paired with unoccupied aerial system (UAS, or drones) imagery to quantify how key locomotor kinematics scale across a body length range of an order of magnitude, from morphologically similar Antarctic minke whales (∼9 m; Balaenoptera bonaerensis, Burmeister 1867) to blue whales (∼30 m) (Fig. 1). We used these data to test the hypothesis that stroke frequency decreases in inverse proportion with body length but instead found it scales to the approximate power of −0.5, while cruising speed remains near 2 m s−1 in order to minimize energy expenditure and thus is expected to be largely invariant with size (Sato et al., 2007; Watanabe et al., 2011, 2015).
Representative overhead unoccupied aerial system (UAS) photographs and scale bars for each study species. (A) Minke, (B) humpback, (C) Bryde's, (D) fin and (E) blue whale.
MATERIALS AND METHODS
All procedures in the USA were conducted under approval of the National Marine Fisheries Service (Permits 781-1824, 16163, 14809, 16111, 19116, 15271, 20430), Canada DFO SARA/MML 2010-01/SARA-106B, National Marine Sanctuaries (MULTI-2017-007) and institutional IACUC committees. South Africa operations were conducted under approval from the Department of Environmental Affairs (Permit RES 2018/63) and Nelson Mandela University (Animal Ethics Approval A18-SCI-ICMR_001).
Tag data collection
We used multi-sensor suction-cup tags (Customized Animal Tracking Solutions, CATS; www.cats.is) to collect high sample rate kinematic and behavioral data from humpback whales (Megaptera novaeangliae, Borowski 1781), blue whales [Balaenoptera musculus (Linnaeus 1758)], fin whales [Balaenoptera physalus (Linnaeus 1758)], Bryde's whales (Balaenoptera brydei, Olsen 1913) and Antarctic minke whales (Balaenoptera bonaerensis, Burmeister 1867). Each tag included tri-axial accelerometers, magnetometers and gyroscopes sampling at 40–400 Hz, and pressure sensors sampling at 10 Hz. Pressure and high sample rate accelerometer data recorded by the tag allowed estimation of swimming speed using the root-mean-square (RMS) amplitude of tag vibrations (Cade et al., 2018). Tag sensors were rotated into whale-body frame orientation, and animal pitch, roll and heading were derived (Cade et al., 2016; Johnson and Tyack, 2003). Cameras collected video at 30 frames s−1 and were integrated into the tags in different orientations: (1) forward along the mid-line axis of the tag, (2) offset at ∼45 deg from the mid-line axis, and (3) forward and backward along the mid-line axis. Additional information about tag design and deployment can be found in Cade et al. (2016). Tag deployments occurred on whales from July 2014 to October 2018 in a wide range of geographic locations including the Antarctic Peninsula for B. bonaerensis, the Azores, Portugal, for B. physalus and B. musculus, Eastern Greenland for B. physalus and M. novaeangliae, Monterey Bay, USA, for B. physalus, B. musculus and M. novaeangliae, Plettenburg Bay, South Africa, for B. brydei, Santa Barbara Channel, USA, for B. musculus and M. novaeangliae, and Stellwagen Bank, USA, for M. novaeangliae (Table S1).
Whale photogrammetry and morphometrics
We collected high-resolution aerial images via UAS between 2017 and 2018 (Fig. 1). Nadir images of whales were collected during a subset of tag deployments for B. bonaerensis (Antarctic Peninsula), M. novaeangliae (Antarctic Peninsula, Eastern Greenland, Monterey Bay and Santa Barbara Channel), B. brydei (Plettenburg Bay), B. physalus (Azores, Eastern Greenland, Monterey Bay) and B. musculus (Monterey Bay, Santa Barbara Channel, Azores). Specifically, we used a DJI Phantom 3 Professional quadcopter, a DJI Phantom 4 Pro quadcopter, and two types of hexacopters, the FreeFly Alta 6 and a custom Mikrokopter-based LemHex-44. The Phantom 3 Professional quadcopter was fitted with a Sony EXMOR 1/2.3 inch camera, 4000×3000 pixel resolution and a 94 deg field of view. The Phantom 4 Pro quadcopter was fitted with a 1 inch CMOS camera sensor, 5472×3078 pixel resolution and an 84 deg field of view. Both hexacopters were fitted with a Lightware SF11/C laser altimeter and a Sony Alpha A5100 camera with an APS-C (23.5×15.6 mm) sensor, 6000×4000 pixel resolution and either a Sony SEL50 50 mm or SEL35 35 mm focal length low distortion lens. The laser altimeter and cameras were co-located on a 2-axis gimbal with pitch angle controlled via remote control to aid in positioning and ensure image collection at the nadir. Measurement error for both quadcopters was calculated as <1.10% by Putch (2017). Measurement error for both hexacopters was calculated by measuring an object of known size placed in the water (i.e. a paddle) at various altitudes. The Alta 6 had a measurement error of <0.30%, while the LemHex-44 measurement error was <1.50%. We used similar methods for hand launch and recovery from small boats to those described in Durban et al. (2015), with the addition of a first-person view (FPV) screen attached to each flight controller, giving the pilot a live feed from the photogrammetry camera. The LemHex-44 required a single operator where the pilot manually controlled the gimbal and camera's shutter, whereas the Alta 6 required two operators, a pilot and camera operator. We collected images in bursts (LemHex-44 and Alta 6) or on a 2 s timer (Phantom 3, Phantom 4) with the whale in full frame lengthwise with a high shutter rate of 6 frames s−1 as the animal surfaced or was just below the surface (Fig. 1).
We manually selected UAS images for measuring the total length and maximum diameter of individual whales if the lower jaw, fluke notch and sides were clearly visible and if the animal appeared straight with minimum curvature at the surface or just below. We measured total length from the tip of the lower jaw to the fluke notch and manually selected maximum diameter as the widest width of the whale posterior of the axilla. We performed all measurements using ImageJ 1.5i (Schindelin et al., 2012). We used the segmented or straight-line tool to draw a line from the tip of the lower jaw to the fluke notch, or from each side at the widest width, to measure the distance in number of pixels. We calculated total length and maximum diameter using an approach similar to Fearnbach et al. (2012), where the number of pixels was multiplied by the ground sampling distance (GSD):
(1)
(2)where L is the total length of the animal (m), npixels is the number of pixels, a is altitude (m), lfocal is focal length (mm), wS is the width of the sensor (mm) and wp is the width of the image resolution (pixels). The width was used for the sensor size and image resolution because the whales were captured full frame widthwise.
Tailbeat measurement
We chose the transverse axis of the gyroscope signal to quantify rotational body movements along the pitch axis of the animal. This improved upon previous studies (Martin López et al., 2015, 2016) which relied on the longitudinal axis of the accelerometer because it allowed us to disregard the orientation of the body and it was highly effective during high pitch angle segments that would otherwise be prone to gimbal lock. For 5 out of 143 deployments, we relied on the longitudinal axis accelerometer method with a high-pass filter because of malfunctioning gyroscope signals. We compared gyroscope and accelerometer signals from the same deployment to quantify uncertainties and ensure that slight differences would not affect our results. For each deployment, we calculated the average value for the oscillatory signal, then used a customized Matlab script to identify full tailbeats, here defined as one consecutive upstroke and downstroke, using a zero-crossing method. We set a consistent and conservative series of thresholds aimed at ignoring incomplete or asymmetrical tailbeats in the data sets. To be included in our analyses, a tailbeat had to have a period measured at less than 10 s, feature exactly two peaks (one upstroke and one downstroke), with a magnitude towards either the upstroke or downstroke being less than 15 times the magnitude in the other direction. We visually inspected the resultant tailbeats to confirm that our thresholds were performing correctly. We removed deployments from subsequent analyses if they identified fewer than 200 full tailbeats (Table 1). We calculated multiple time-synchronized kinematic and behavioral parameters for the duration of each full tailbeat (Fig. 2) including the oscillatory frequency, f, as 1/period, the average swimming speed (Cade et al., 2018), the average depth and the modal behavioral state (described below).
Metadata related to each included species
Section of CATS tag data showing time-synced data streams. The video stills are faced caudally from a common reference point and illustrate a single tailbeat cycle. These stills are correlated in time with the vertical red bar intersecting each data stream. The green horizontal lines at the zero line of the first data stream (Gyroscope) correspond to symmetrical tailbeats that fit our detection thresholds and were included in our analyses. The three lines in the third data stream (Orientation) correspond to pitch (right y-axis; orange), roll (left y-axis; pink) and heading (left y-axis; blue) changes.
Measurement of speed
Measurement of the swim speed relied on turbulent flow to vibrate the tag and, as a result, could not resolve speeds below ∼1 m s−1 (Cade et al., 2018). For each deployment, we regressed the oscillatory frequency of each individual tailbeat period against the average speed during that stroke. To ensure that our regressions were accurately capturing the relationship between oscillatory frequency and swimming speed, we only included tailbeat periods with an average speed higher than one standard deviation above the lowest recorded speed. So as not to bias results with artificially high minimum speeds as a result of the method's low-speed threshold, we used median swimming speed values for comparative analyses among species.
Modeling optimal speed
We combined measurements of body total length and maximum diameter with a simple model of non-feeding transport to predict the swimming speed associated with the minimum energetic cost of transport (Weihs, 1973). The modeling starts with estimating the (metabolic) energy expended (ε; J) to travel a given distance (d) as follows:
(3)The second term corresponds to the energy expended by the locomotor musculature to compensate for the drag produced (as propulsive drag equals drag in this instance). Here, d is the distance traveled (set to 1 body length), U is the swimming speed (m s−1), ρ is the density of seawater (=1025 kg m−3), S is the wetted surface area (m2) of the animal calculated using coefficients from computational fluid dynamics modeling or literature sources for each included species (Bose and Lien, 1989; Kermack, 1948; Parry, 1949), Cd is the dimensionless drag coefficient of the body in a rigid configuration (∼0.002–0.005) (Kooyman, 1989), λ is the drag correction due to the heaving body and caudal tail (∼2–3) (Fish, 1993, 1998), ηm is the (dimensionless) metabolic efficiency and ηp is the (dimensionless) propulsive efficiency. Both of the efficiency terms are considered to be constant across body size ranges (Fish, 1996, 2001; Fish and Rohr, 1999). The first term in Eqn 3 accounts for the energy expenditure by the organ and tissue in the rest of the body, which here is estimated by the Kleiber (1975) basal metabolic rate (BMR) correlation, corrected by a factor of 1.6, that is based upon previous estimates of swimming or dive recovery metabolic rates (Costa and Williams, 1999; Lockyer, 1981; Scholander, 1940) and accounts for a higher level of metabolic activity during active swimming over long durations:
(4)where m is body mass (kg). The model yields an optimal swimming speed (Uopt) for which a minimum of energy stored in the body is used to cover distance d (Weihs, 1973):
(5)
Tail-attached tag measurements
As for fish swimming in the wild (Han et al., 2017; Watanabe et al., 2019), stroke amplitude is difficult to measure in tagged cetaceans (van der Hoop et al., 2017). However, during a single humpback whale deployment in 2017, the tag slipped from its original placement on the back of the animal (between the dorsal fin and blowhole) and attached to the dorsal aspect of the tail near the right fluke-blade midline. As shown in Fig. 3, the tag was aimed anteriorly towards the body and remained so for 68 min. Such a fortuitous event allowed us to measure the tailbeat amplitude as a local, short-duration change in depth from the tag's pressure transducer. We included sequences with at least three tailbeats that had an overall depth change near zero (n=14). For each sequence, we calculated the midline – from the midpoint of the first tailbeat to the midpoint of the final tailbeat – and then averaged the distance from that midpoint to the top of each upstroke and to the bottom of each downstroke (Fig. 3 inset). For these same tailbeat sequences, we measured the oscillatory frequency and swimming speed, and estimated the Strouhal number for each sequence.
Tailbeat amplitude measurement. A video still from the deployment of interest. (A) Rostrum, (B) dorsal fin, (C) left flipper, (D) right flipper, (E) peduncle, (F) left fluke blade and (G) right fluke blade. The inset is a schematic diagram of our peak-to-peak amplitude measurement method.
Comparison with two simple oscillatory models
We used a derivation of the Strouhal number equation to estimate the oscillatory frequency of an optimally oscillating vortex generator:
(6)with U set at 2 (m s−1) based upon the average speeds in our observed dataset, St set at 0.3 to simulate optimality and A set at one-fifth of body length (Bainbridge, 1958; Fish and Rohr, 1999; Rohr and Fish, 2004). We applied this model over a 10–20 m range of total body lengths to facilitate a direct comparison with our observed dataset.
An alternative to Strouhal number scaling is the idea that, to sustain a form of low-cost swimming, the oscillatory frequency is set to the resonant frequency of the caudal tail, here seen as a driven and damped oscillator (Ahlborn et al., 2009; Hamilton et al., 2004; Pabst, 1996, 2000). In this model, drag provides the damping, the tail musculature provides the driving force, and the sub-dermal sheath and blubber surrounding the tail stock acts as the spring-like oscillatory restoring force. Swimming near the ‘optimal’ oscillatory frequency means that the driving force will be set at the oscillatory frequency that generates the greatest tailbeat amplitude. In low-damping conditions, this oscillatory frequency will coincide with the so-called ‘natural’ frequency of the tail tissue and bone under the sole influence of the spring-like restoring force provided by the sub-dermal sheath (Kreyszig, 2016). With the sub-dermal sheath forces acting like a cantilever beam oscillating in seawater, the oscillatory frequency would thus be modeled as:
(7)where E parameterizes the effective Young's modulus of cetacean tissue (approximated as 2 MPa from Long et al., 2002), I is the second moment of area (m4) and me is the effective mass (kg) of the beam (Denny, 1988); we doubled that effective mass to account for the added mass of the body accelerating in a fluid (Lighthill, 1970; Webb, 1975). Here, the model assumed that the cranial 40% of body length was stationary during oscillation, with the caudal 60% of body length oscillating (Los). We applied this model over the same 10–20 m range of total body lengths as our Strouhal model.
Statistical analyses
The datasets supporting this article are available from the corresponding author on request. We performed all statistical analyses using R (version 3.4.1, package nlme) with a significance level of 0.05. For each analysis, continuous variables (total length, oscillatory frequency, median cruising speed, coefficient of variation) were log10 transformed before inclusion as predictors or response variables to normalize our data and conform to the basic model of scaling as a power function. For the subset of deployments with morphometric measurements (n=55), we created a linear mixed-effects model with body length as the predictor, the mean oscillatory frequency of a deployment as the response and species as a random effect. The results of this model were used to determine the regression line and standard error of the mean for our comparison of body length and oscillatory frequency. Next, we created a linear mixed-effects model with body length as the predictor, the coefficient of variation for the oscillatory frequency of a deployment as the response and species as a random effect. Finally, we created a linear mixed-effects model with body length as the predictor, median swimming speed as the response and species as a random effect.
We used our complete dataset to create a linear mixed model comparing the oscillatory frequency between dive descent and dive ascent stages across all species. The model included species and individuals within species as random effects. We also performed identical analyses on a species-specific basis for B. bonaerensis, M. novaeangliae and B. musculus. These models included individuals within species as a random effect. Our final models included swimming phase as the predictor, coefficient of variation for the oscillatory frequency of a deployment as the response and species as a random effect.
RESULTS
Variation in oscillatory kinematics
Metadata and summary statistics for tagged whales (n=143 total, n=60 with UAS-measured morphometrics) are shown in Table 1. The number of tailbeats in a deployment ranged from 260 to 9458 with a mean (±s.d.) of 1927±1784. Swim speeds ranged up to 8.28 m s−1 (lower bound was unresolved) with a median of 1.89±0.038 (median±standard error) and oscillatory frequencies ranged from 0.10 to 1.00 Hz with a mean (±s.e.m.) of 0.23±0.004 Hz (Fig. 4). Swim speed increased with an increase in oscillatory frequency in all species (Fig. 5). For all species combined, we calculated a single regression equation of swim speed on oscillatory frequency (
=0.54x−0.80; R2=0.24). For the lone humpback whale tail-attached deployment, we found the mean (±s.d.) tailbeat amplitude (2.63±0.79 m), oscillatory frequency (0.25±0.06 m), swimming speed (2.61±0.25 m) and Strouhal number (0.24±0.04).
Density distributions of the oscillatory frequency calculated for each included tailbeat for a single whale. Distributions have been grouped by species. The box plots show the mean, 25% and 75% percentiles, and whiskers moving out to the farthest outliers.
Regression between the oscillatory frequency (Hz) and average speed (m s−1) of each whale. Individuals of each species are represented by the colored regression lines. The black regression lines are the averages for each species. In the ‘All’ plot (bottom right), each of the species-averaged regression lines is matched to the initial line color chosen for that species. The number of individuals (N) and the mean (±s.d.) number of tailbeats is shown in each panel. Additional information (R2 and slope) is given in the form of density distributions grouped by species.
Effect of body length on oscillatory frequency and swimming speed
Fig. 6 compares the tag-based measurements of the heaving frequency with predictions from optimal Strouhal scaling (Eqn 6) and oscillating beam models (Eqn 7). The same figure also shows a comparison between measurements of steady-state swimming speed and the predicted optimal speed as per Eqn 5. There was a clear decrease in oscillatory frequency (fobs) with increasing total length in both an intraspecific and interspecific context, but no relationship between cruising speed and total length (Table 1). The regression for observed speed (Uobs) corresponded very closely with the regression for optimal speed (Uopt) (Table 2). Overall, the measured swimming speed followed the insensitivity to body length exhibited by the model. However, the scaling of our oscillatory frequency data f≈L−0.5 does not agree with the f≈L−1.0 trend suggested by both the St and oscillating cantilever models.
Relationships between body length and both oscillatory frequency and cruising speed. (A) Interspecific log-transformed mean oscillatory frequency (Hz) for all species. (B) Intraspecific humpback mean oscillatory frequency. (C) Interspecific log-transformed median cruising speed (m s−1) for all species. (D) Intraspecific humpback median cruising speed. Solid lines denote the observed regressions and gray shading represents the standard error. We calculated the regression for A using our linear mixed model with species as a random effect. The dot-dash and dotted lines represent the regressions for our cantilever beam and Strouhal models, respectively. The open circles shown in C represent the modeled Uopt values for each whale. For each point in C and D, a single vertical line extends from the 25th percentile (bottom of the line) to the 75th percentile (top of the line), with the point in the middle denoting the median.
Equations and R2 values calculated for each regression shown in Fig. 6
More specifically, we found a significant dependence on total body length in the mean oscillatory frequency for an individual whale of any species (d.f.=54; t=−6.69; P<0.001). The effect of total length on the coefficient of variation of the oscillatory frequency for an individual whale was not significant (d.f.=54; t=−0.62; P=0.54) and, similarly, the effect of total length on the median cruising speed for an individual whale was not significant (d.f.=54; t=0.90; P=0.37).
DISCUSSION
The scale dependence of locomotor performance
Our measurements from baleen whales of varying body lengths indicate that while oscillatory frequency decreases with body length to the power of −0.53, cruising speed is largely invariant with body length (to the power of 0.08). These results are broadly consistent with previous studies for other swimming animals (Sato et al., 2007), particularly endotherms (Watanabe et al., 2015). The empirically derived oscillatory frequencies and cruising speeds found in our investigation can inform studies attempting to model these kinematic parameters across wide body length ranges (Alexander, 2005; Vogel, 2008). Many studies have included inexact measurements for large swimming animals or excluded them completely from the modeling (Bejan and Marden, 2006; Gazzola et al., 2014). The model put forth by Gazzola et al. (2014) includes an average blue whale cruising speed of 6 m s−1 and a resulting oscillatory frequency of 0.36 Hz. These values are roughly triple what we found in the present study using empirical data for the same species performing natural behaviors.
Our observed scaling exponent of −0.53, derived by least-squares regression of f onto L, differs from that predicted by the two simple scaling models. The cantilevered beam model (Eqn 7) – taking as inputs constant values of E and I, and variable values of L and M – predicts a scaling exponent of −1, nearly double what we observed. We obtained the same exponent, −1, in the Strouhal model, which takes as inputs a constant speed (2 m s−1), constant Strouhal number (0.3) and the variable L. The fact that the observed exponent differs from the predicted suggests that these whales are not operating as simple cantilevered beams or as optimally oscillating vortex generators. The cantilevered beam model lacks damping forces that would reduce the predicted exponent. The Strouhal model assumes optimal fluking based on efficient vortex shedding and interaction. The apparent failure of these models points us toward appropriate ways to complicate subsequent models through added damping or relaxation of optimality constraints.
Our results corroborate previous controlled studies and lend support for a set of hydrodynamic principles governing the locomotion of free-swimming animals down to the level of an individual stroke (Bainbridge, 1958; Fish, 1998; Kohannim and Iwasaki, 2014). The difficulty of obtaining stroke amplitude measures for our study animals is a limitation that precludes the type of precise hydrodynamic analysis found in previous studies, but the stroke amplitude determined from the tail-attached deployment corresponds well with previous estimates of stroke amplitude as one-fifth of body length (Bainbridge, 1958; Fish and Rohr, 1999; Rohr and Fish, 2004). Although we lacked UAS-based measures of the tail-attached animal, a calculation for our average amplitude measure of 2.37 m suggests a total body length of 11.85 m. For comparison, the mean (±s.d.) body length for humpback whales within our dataset was 11.09±1.81 m. Previous studies have shown amplitude to be invariant with changes in forward speed, so the lack of amplitude measures for each stroke does not considerably lessen our understanding of the animal's overall hydrodynamic scope (Bainbridge, 1958; Fish, 1998; Kohannim and Iwasaki, 2014). Using the tail-attached deployment, we also obtained the first empirical measurements of Strouhal number for a large cetacean. These values fall near the predicted range of St≈0.25–0.35 that has been shown to correspond to high hydrodynamic efficiency (Anderson et al., 1998; Rohr and Fish, 2004; Taylor et al., 2003; Triantafyllou et al., 1991). Throughout the duration of this deployment, the animal transited (moved steadily from one location to another) and performed a single feeding lunge. Of all the behaviors performed by animals in our dataset, transiting is most similar to the steady-state swimming performed by animals in controlled, laboratory and aquarium studies (Bainbridge, 1958; Fish and Rohr, 1999; Fish et al., 2014). This could explain why the Strouhal number for this individual fell so close to the theoretical optimum.
One of the common methods used to quantify the efficiency of locomotor modes is calculating the cost of transport over a set distance, either as a mass-specific or as an objective estimation (Alexander, 1999; Schmidt-Nielsen, 1972; Watanabe et al., 2011; Williams, 1999; Williams et al., 2014). Calculation of the cost of transport typically requires an estimate of metabolic rate using oxygen consumption or some other metric of metabolism (Schmidt-Nielsen, 1972; Sumich, 1983; Williams, 1999; Williams et al., 1993). Here, from first principles we predicted the forward speed at which the energetic cost of transport is minimized. This is a common optimization estimate and can be useful as a comparison against measured speed (Alexander, 1999; Watanabe et al., 2011). The results of our modeled optimum speed agree with previous research showing that the preferred speed of many endothermic swimmers is 1–2 m s−1 (Sato et al., 2007; Watanabe et al., 2011, 2015). The slight increase in observed swimming speed above our modeled values could be due to increased behavioral variation throughout the course of the deployment. Our optimum swimming speed equation assumes that an animal is in a steady-state, non-feeding locomotor mode. Conversely, our observed values include maneuvers and non-steady swimming related to feeding or other ecologically relevant functions, many of which are being performed above typical cruising speeds (Cade et al., 2016; Goldbogen et al., 2006, 2011, 2012; Simon et al., 2012). Han et al. (2017) performed a similar comparison between optimal and observed swimming speeds in largemouth bass (Micropterus salmoides) and found the opposite result. Bass routinely swam slower than their predicted optimum for cost of transport. It was suggested by Han et al. (2017) that other considerations may be prioritized in the wild for largemouth bass, such as the need for slower speeds in order to efficiently detect and capture food.
Conclusion
Our study is the first to quantify the kinematics of free-swimming baleen whales over an order of magnitude in body length. Our results for oscillatory frequency and cruising speed are more accurate than previous models and show that, while oscillatory frequency decreases with increasing body length, cruising speed remains largely invariant. A model predicting the optimal swimming speed aligned closely with the cruising speeds observed in our study animals. Models predicated on Strouhal number and an oscillating cantilever beam could not predict the full variation in oscillatory frequency within our dataset, but a small number of empirical measurements of the Strouhal number fell within the range for high hydrodynamic efficiency. Major deviations from theory may reflect the competing demands for different locomotor functions during common behaviors such as foraging.
Acknowledgements
Our fieldwork efforts would not have been successful without our boat crew and logistics personnel. Dr Nikolai Liebsch and Dr Peter Kraft from Customized Animal Tracking Solutions have been instrumental in the design and maintenance of our inertial sensing tags. Xitlaali Castellanos performed the Strouhal measurements for the humpback tail-attached deployment. Dr Robin Elahi assisted with the design and implementation of our statistical analyses.
FOOTNOTES
Competing interests
The authors declare no competing or financial interests.
Author contributions
Conceptualization: W.T.G., P.S.S., D.E.C., J.P., F.E.F., A.S.F., J.H.L., M.S.S., J.A.G.; Methodology: W.T.G., P.S.S., K.C.B., D.E.C., J.P., F.E.F., A.S.F., J.K., J.H.L., M.S.S., J.A.G.; Software: W.T.G., D.E.C., J.K.; Formal analysis: W.T.G., K.C.B., J.P., J.K.; Investigation: W.T.G., P.S.S., K.C.B., D.E.C., J.P., F.E.F., J.D., J.d.C., A.S.F., D.W.J., S.R.K.-R., M.O., G.P., M.S.S., M.S., S.K.V., F.V., D.N.W., J.A.G.; Resources: K.C.B., D.E.C., J.P., J.D., J.d.C., A.S.F., D.W.J., J.K., M.O., G.P., M.S., S.K.A.V., F.V., D.N.W., J.A.G.; Data curation: W.T.G., P.S.S., K.C.B., D.E.C., J.D., A.S.F., D.W.J., S.R.K.-R., J.K., M.O., G.P., M.S.S., M.S., S.K.A.V., F.V., D.N.W.; Writing - original draft: W.T.G.; Writing - review & editing: W.T.G., P.S.S., K.C.B., D.E.C., J.P., F.E.F., J.D., J.d.C., A.S.F., D.W.J., S.R.K.-R., J.K., J.H.L., M.O., G.P., M.S.S., M.S., S.K.A.V., F.V., D.N.W., J.A.G.; Visualization: W.T.G., M.S.S., J.A.G.; Supervision: J.P., F.E.F., A.S.F., J.H.L., M.S.S., J.A.G.; Project administration: A.S.F., J.A.G.; Funding acquisition: J.P., F.E.F., A.S.F., J.A.G.
Funding
This research was funded in part by grants from the National Science Foundation (IOS-1656691, IOS-1656676, IOS-1656656; OPP-1644209), the Office of Naval Research (N000141612477), a Terman Fellowship from Stanford University, and a National Science Foundation INSPIRE Special Projects grant (1344227). Additional Monterey Bay CATS tag deployments were funded by grants from the American Cetacean Society Monterey and San Francisco Bay chapters, and by the Meyers Trust. Bryde's whale tagging operations in South Africa were funded by the Percy Sladen Memorial Fund, PADI Foundation, Torben and Alice Frimodts Fund, and Society for Marine Mammology.
Supplementary information
Supplementary information available online at http://jeb.biologists.org/lookup/doi/10.1242/jeb.204172.supplemental
- Received March 29, 2019.
- Accepted September 24, 2019.
- © 2019. Published by The Company of Biologists Ltd
References
Article navigation
Related articles
Cited by...
More in this TOC section
Similar articles
Subject collections
Other journals from The Company of Biologists
Meet the Editors at SICB Virtual 2021
Reserve your place to join some of the journal editors, including Editor-in-Chief Craig Franklin, at our Meet the Editor session on 17 February at 2pm (EST). Don’t forget to view our SICB Subject Collection, featuring relevant JEB papers relating to some of the symposia sessions.
2020 at The Company of Biologists
Despite 2020's challenges, we were able to bring a number of long-term projects and new ventures to fruition. As we enter a new year, join us as we reflect on the triumphs of the last 12 months.
Critical temperature window sends migratory black-headed buntings on their travels
The spring rise in temperature at black-headed bunting overwintering sites is essential for triggering the physical changes that they undergo before embarking on their spring migration – read more.
Developmental and reproductive physiology of small mammals at high altitude
Cayleih Robertson and Kathryn Wilsterman focus on high-altitude populations of the North American deer mouse in their review of the challenges and evolutionary innovations of pregnant and nursing small mammals at high altitude.
Read & Publish participation extends worldwide
“Being able to publish Open Access articles free of charge means that my article gets maximum exposure and has maximum impact, and that all my peers can read it regardless of the agreements that their universities have with publishers.”
Professor Roi Holzman (Tel Aviv University) shares his experience of publishing Open Access as part of our growing Read & Publish initiative. We now have over 60 institutions in 12 countries taking part – find out more and view our full list of participating institutions.