Assessing variation in faecal glucocorticoid concentrations in gray whales exposed to anthropogenic stressors

Understanding how individual animals respond to stressors behaviourally and physiologically is a critical step towards quantifying long-term population consequences and informing management efforts. Glucocorticoid (GC) metabolite accumulation in various matrices provides an integrated measure of adrenal activation in baleen whales and could thus be used to investigate physiological changes following exposure to stressors. In this study, we measured GC concentrations in faecal samples of Pacific Coast Feeding Group (PCFG) gray whales (Eschrichtius robustus) collected over seven consecutive years to assess the association between GC content and metrics of exposure to sound levels and vessel traffic at different temporal scales, while controlling for contextual variables such as sex, reproductive status, age, body condition, year, time of year and location. We develop a Bayesian Generalized Additive Modelling approach that accommodates the many complexities of these data, including non-linear variation in hormone concentrations, missing covariate values, repeated samples, sampling variability and some hormone concentrations below the limit of detection. Estimated relationships showed large variability, but emerging patterns indicate a strong context-dependency of physiological variation, depending on sex, body condition and proximity to a port. Our results highlight the need to control for baseline hormone variation related to context, which otherwise can obscure the functional relationship between faecal GCs and stressor exposure. Therefore, extensive data collection to determine sources of baseline variation in well-studied populations, such as PCFG gray whales, could shed light on cetacean stress physiology and be used to extend applicability to less-well-studied taxa. GC analyses may offer greatest utility when employed as part of a suite of markers that, in aggregate, provide a multivariate measure of physiological status, better informing estimates of individuals' health and ultimately the consequences of anthropogenic stressors on populations.

Downsized: gray whales using an alternative foraging ground have smaller morphology

Describing individual morphology and growth is key for identifying ecological niches and monitoring the health and fitness of populations. Eastern North Pacific ((ENP), approximately 16 650 individuals) gray whales primarily feed in the Arctic/sub-Arctic regions, while a small subgroup called the Pacific Coast Feeding Group (PCFG, approximately 212 individuals) instead feeds between northern California, USA and British Columbia, Canada. Evidence suggests PCFG whales have lower body condition than ENP whales. Here we investigate morphological differences (length, skull, and fluke span) and compare length-at-age growth curves between ENP and PCFG whales. We use ENP gray whale length-at-age data comprised of strandings, whaling, and aerial photogrammetry (1926-1997) for comparison to data from PCFG whales collected through non-invasive techniques (2016-2022) to estimate age (photo identification) and length (drone-based photogrammetry). We use Bayesian methods to incorporate uncertainty associated with morphological measurements (manual and photogrammetric) and age estimates. We find that while PCFG and ENP whales have similar growth rates, PCFG whales reach smaller asymptotic lengths. Additionally, PCFG whales have relatively smaller skulls and flukes than ENP whales. These findings represent a striking example of morphological adaptation that may facilitate PCFG whales accessing a foraging niche distinct from the Arctic foraging grounds of the broader ENP population.

Assessment of a non-invasive approach to pregnancy diagnosis in gray whales through drone-based photogrammetry and faecal hormone analysis

Knowledge of baleen whales' reproductive physiology is limited and requires long-term individual-based studies and innovative tools. We used 6 years of individual-level data on the Pacific Coast Feeding Group gray whales to evaluate the utility of faecal progesterone immunoassays and drone-based photogrammetry for pregnancy diagnosis. We explored the variability in faecal progesterone metabolites and body morphology relative to observed reproductive status and estimated the pregnancy probability for mature females of unknown reproductive status using normal mixture models. Individual females had higher faecal progesterone concentrations when pregnant than when presumed non-pregnant. Yet, at the population level, high overlap and variability in progesterone metabolite concentrations occurred between pregnant and non-pregnant groups, limiting this metric for accurate pregnancy diagnosis in gray whales. Alternatively, body width at 50% of the total body length (W50) correctly discriminated pregnant from non-pregnant females at individual and population levels, with high accuracy. Application of the model using W50 metric to mature females of unknown pregnancy status identified eight additional pregnancies with high confidence. Our findings highlight the utility of drone-based photogrammetry to non-invasively diagnose pregnancy in this group of gray whales, and the potential for improved data on reproductive rates for population management of baleen whales generally.

Minke whale feeding rate limitations suggest constraints on the minimum body size for engulfment filtration feeding

Bulk filter feeding has enabled gigantism throughout evolutionary history. The largest animals, extant rorqual whales, utilize intermittent engulfment filtration feeding (lunge feeding), which increases in efficiency with body size, enabling their gigantism. The smallest extant rorquals (7-10 m minke whales), however, still exhibit short-term foraging efficiencies several times greater than smaller non-filter-feeding cetaceans, raising the question of why smaller animals do not utilize this foraging modality. We collected 437 h of bio-logging data from 23 Antarctic minke whales (Balaenoptera bonaerensis) to test the relationship of feeding rates (λf) to body size. Here, we show that while ultra-high nighttime λf (mean ± s.d.: 165 ± 40 lunges h-1; max: 236 lunges h-1; mean depth: 28 ± 46 m) were indistinguishable from predictions from observations of larger species, daytime λf (mean depth: 72 ± 72 m) were only 25-40% of predicted rates. Both λf were near the maxima allowed by calculated biomechanical, physiological and environmental constraints, but these temporal constraints meant that maximum λf was below the expected λf for animals smaller than ~5 m-the length of weaned minke whales. Our findings suggest that minimum size for specific filter-feeding body plans may relate broadly to temporal restrictions on filtration rate and have implications for the evolution of filter feeding.

Demography of an ice-obligate mysticete in a region of rapid environmental change

Antarctic minke whales (Balaenoptera bonaerensis, AMW) are an abundant, ice-dependent species susceptible to rapid climatic changes occurring in parts of the Antarctic. Here, we used remote biopsy samples and estimates of length derived from unoccupied aircraft system (UAS) to characterize for the first time the sex ratio, maturity, and pregnancy rates of AMWs around the Western Antarctic Peninsula (WAP). DNA profiling of 82 biopsy samples (2013-2020) identified 29 individual males and 40 individual females. Blubber progesterone levels indicated 59% of all sampled females were pregnant, irrespective of maturity. When corrected for sexual maturity, the median pregnancy rate was 92.3%, indicating that most mature females become pregnant each year. We measured 68 individuals by UAS (mean = 8.04 m) and estimated that 66.5% of females were mature. This study provides the first data on the demography of AMWs along the WAP and represents the first use of non-lethal approaches to studying this species. Furthermore, these results provide baselines against which future changes in population status can be assessed in this rapidly changing marine ecosystem.

Fast and Furious: Energetic Tradeoffs and Scaling of High-Speed Foraging in Rorqual Whales

Although gigantic body size and obligate filter feeding mechanisms have evolved in multiple vertebrate lineages (mammals and fishes), intermittent ram (lunge) filter feeding is unique to a specific family of baleen whales: rorquals. Lunge feeding is a high cost, high benefit feeding mechanism that requires the integration of unsteady locomotion (i.e., accelerations and maneuvers); the impact of scale on the biomechanics and energetics of this foraging mode continues to be the subject of intense study. The goal of our investigation was to use a combination of multi-sensor tags paired with UAS footage to determine the impact of morphometrics such as body size on kinematic lunging parameters such as fluking timing, maximum lunging speed, and deceleration during the engulfment period for a range of species from minke to blue whales. Our results show that, in the case of krill-feeding lunges and regardless of size, animals exhibit a skewed gradient between powered and fully unpowered engulfment, with fluking generally ending at the point of both the maximum lunging speed and mouth opening. In all cases, the small amounts of propulsive thrust generated by the tail were unable to overcome the high drag forces experienced during engulfment. Assuming this thrust to be minimal, we predicted the minimum speed of lunging across scale. To minimize the energetic cost of lunge feeding, hydrodynamic theory predicts slower lunge feeding speeds regardless of body size, with a lower boundary set by the ability of the prey to avoid capture. We used empirical data to test this theory and instead found that maximum foraging speeds remain constant and high (∼4 m s^–1) across body size, even as higher speeds result in lower foraging efficiency. Regardless, we found an increasing relationship between body size and this foraging efficiency, estimated as the ratio of energetic gain from prey to energetic cost. This trend held across timescales ranging from a single lunge to a single day and suggests that larger whales are capturing more prey—and more energy—at a lower cost.

Scaling of maneuvering performance in baleen whales: larger whales outperform expectations

Despite their enormous size, whales make their living as voracious predators. To catch their much smaller, more maneuverable prey, they have developed several unique locomotor strategies that require high energetic input, high mechanical power output and a surprising degree of agility. To better understand how body size affects maneuverability at the largest scale, we used bio-logging data, aerial photogrammetry and a high-throughput approach to quantify the maneuvering performance of seven species of free-swimming baleen whale. We found that as body size increases, absolute maneuvering performance decreases: larger whales use lower accelerations and perform slower pitch-changes, rolls and turns than smaller species. We also found that baleen whales exhibit positive allometry of maneuvering performance: relative to their body size, larger whales use higher accelerations, and perform faster pitch-changes, rolls and certain types of turns than smaller species. However, not all maneuvers were impacted by body size in the same way, and we found that larger whales behaviorally adjust for their decreased agility by using turns that they can perform more effectively. The positive allometry of maneuvering performance suggests that large whales have compensated for their increased body size by evolving more effective control surfaces and by preferentially selecting maneuvers that play to their strengths.

Baleen whale prey consumption based on high-resolution foraging measurements

Baleen whales influence their ecosystems through immense prey consumption and nutrient recycling^1-3. It is difficult to accurately gauge the magnitude of their current or historic ecosystem role without measuring feeding rates and prey consumed. To date, prey consumption of the largest species has been estimated using metabolic models^3-9 based on extrapolations that lack empirical validation. Here, we used tags deployed on seven baleen whale (Mysticeti) species (n = 321 tag deployments) in conjunction with acoustic measurements of prey density to calculate prey consumption at daily to annual scales from the Atlantic, Pacific, and Southern Oceans. Our results suggest that previous studies^3-9 have underestimated baleen whale prey consumption by threefold or more in some ecosystems. In the Southern Ocean alone, we calculate that pre-whaling populations of mysticetes annually consumed 430 million tonnes of Antarctic krill (Euphausia superba), twice the current estimated total biomass of E. superba¹⁰, and more than twice the global catch of marine fisheries today¹¹. Larger whale populations may have supported higher productivity in large marine regions through enhanced nutrient recycling: our findings suggest mysticetes recycled 1.2 × 10⁴ tonnes iron yr^-1 in the Southern Ocean before whaling compared to 1.2 × 10³ tonnes iron yr^-1 recycled by whales today. The recovery of baleen whales and their nutrient recycling services^2,3,7 could augment productivity and restore ecosystem function lost during 20th century whaling^12,13.

Scaling of oscillatory kinematics and Froude efficiency in baleen whales

High efficiency lunate-tail swimming with high-aspect-ratio lifting surfaces has evolved in many vertebrate lineages, from fish to cetaceans. Baleen whales (Mysticeti) are the largest swimming animals that exhibit this locomotor strategy, and present an ideal study system to examine how morphology and the kinematics of swimming scale to the largest body sizes. We used data from whale-borne inertial sensors coupled with morphometric measurements from aerial drones to calculate the hydrodynamic performance of oscillatory swimming in six baleen whale species ranging in body length from 5 to 25 m (fin whale, Balaenoptera physalus; Bryde's whale, Balaenoptera edeni; sei whale, Balaenoptera borealis; Antarctic minke whale, Balaenoptera bonaerensis; humpback whale, Megaptera novaeangliae; and blue whale, Balaenoptera musculus). We found that mass-specific thrust increased with both swimming speed and body size. Froude efficiency, defined as the ratio of useful power output to the rate of energy input ( Sloop, 1978), generally increased with swimming speed but decreased on average with increasing body size. This finding is contrary to previous results in smaller animals, where Froude efficiency increased with body size. Although our empirically parameterized estimates for swimming baleen whale drag were higher than those of a simple gliding model, oscillatory locomotion at this scale exhibits generally high Froude efficiency as in other adept swimmers. Our results quantify the fine-scale kinematics and estimate the hydrodynamics of routine and energetically expensive swimming modes at the largest scale.

Lunge filter feeding biomechanics constrain rorqual foraging ecology across scale

Fundamental scaling relationships influence the physiology of vital rates, which in turn shape the ecology and evolution of organisms. For diving mammals, benefits conferred by large body size include reduced transport costs and enhanced breath-holding capacity, thereby increasing overall foraging efficiency. Rorqual whales feed by engulfing a large mass of prey-laden water at high speed and filtering it through baleen plates. However, as engulfment capacity increases with body length (engulfment volume∝body length3.57), the surface area of the baleen filter does not increase proportionally (baleen area∝body length1.82), and thus the filtration time of larger rorquals predictably increases as the baleen surface area must filter a disproportionally large amount of water. We predicted that filtration time should scale with body length to the power of 1.75 (filter time∝body length1.75). We tested this hypothesis on four rorqual species using multi-sensor tags with corresponding unoccupied aircraft systems-based body length estimates. We found that filter time scales with body length to the power of 1.79 (95% CI: 1.61-1.97). This result highlights a scale-dependent trade-off between engulfment capacity and baleen area that creates a biomechanical constraint to foraging through increased filtration time. Consequently, larger whales must target high-density prey patches commensurate to the gulp size to meet their increased energetic demands. If these optimal patches are absent, larger rorquals may experience reduced foraging efficiency compared with smaller whales if they do not match their engulfment capacity to the size of targeted prey aggregations.

Energetic and physical limitations on the breaching performance of large whales

The considerable power needed for large whales to leap out of the water may represent the single most expensive burst maneuver found in nature. However, the mechanics and energetic costs associated with the breaching behaviors of large whales remain poorly understood. In this study we deployed whale-borne tags to measure the kinematics of breaching to test the hypothesis that these spectacular aerial displays are metabolically expensive. We found that breaching whales use variable underwater trajectories, and that high-emergence breaches are faster and require more energy than predatory lunges. The most expensive breaches approach the upper limits of vertebrate muscle performance, and the energetic cost of breaching is high enough that repeated breaching events may serve as honest signaling of body condition. Furthermore, the confluence of muscle contractile properties, hydrodynamics, and the high speeds required likely impose an upper limit to the body size and effectiveness of breaching whales.

Scaling of swimming performance in baleen whales

The scale-dependence of locomotor factors have long been studied in comparative biomechanics, but remain 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.