Baleen whale prey consumption based on high-resolution foraging measurements.

Baleen whales influence their ecosystems through immense prey consumption and nutrient recycling1-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 models3-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 studies3-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. superba10, and more than twice the global catch of marine fisheries today11. Larger whale populations may have supported higher productivity in large marine regions through enhanced nutrient recycling: our findings suggest mysticetes recycled 1.2 × 104 tonnes iron yr-1 in the Southern Ocean before whaling compared to 1.2 × 103 tonnes iron yr-1 recycled by whales today. The recovery of baleen whales and their nutrient recycling services2,3,7 could augment productivity and restore ecosystem function lost during 20th century whaling12,13.

Modeling individual growth reveals decreasing gray whale body length and correlations with ocean climate indices at multiple scales.

Changes in body size have been documented across taxa in response to human activities and climate change. Body size influences many aspects of an individual's physiology, behavior, and ecology, ultimately affecting life history performance and resilience to stressors. In this study, we developed an analytical approach to model individual growth patterns using aerial imagery collected via drones, which can be used to investigate shifts in body size in a population and the associated drivers. We applied the method to a large morphological dataset of gray whales (Eschrichtius robustus) using a distinct foraging ground along the NE Pacific coast, and found that the asymptotic length of these whales has declined since around the year 2000 at an average rate of 0.05-0.12 m/y. The decline has been stronger in females, which are estimated to be now comparable in size to males, minimizing sexual dimorphism. We show that the decline in asymptotic length is correlated with two oceanographic metrics acting as proxies of habitat quality at different scales: the mean Pacific Decadal Oscillation index, and the mean ratio between upwelling intensity in a season and the number of relaxation events. These results suggest that the decline in gray whale body size may represent a plastic response to changing environmental conditions. Decreasing body size could have cascading effects on the population's demography, ability to adjust to environmental changes, and ecological influence on the structure of their community. This finding adds to the mounting evidence that body size is shrinking in several marine populations in association with climate change and other anthropogenic stressors. Our modeling approach is broadly applicable across multiple systems where morphological data on megafauna are collected using drones.

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.

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.

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.

Scaling of swimming performance in baleen whales.

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.

Temporal and Regional Variability in the Skin Microbiome of Humpback Whales along the Western Antarctic Peninsula.

The skin is the first line of defense between an animal and its environment, and disruptions in skin-associated microorganisms can be linked to an animal's health and nutritional state. To better understand the skin microbiome of large whales, high-throughput sequencing of partial small subunit rRNA genes was used to study the skin-associated bacteria of 89 seemingly healthy humpback whales (Megaptera novaeangliae) sampled along the Western Antarctic Peninsula (WAP) during early (2010) and late (2013) austral summers. Six core groups of bacteria were present in 93% or more of all humpback skin samples. A shift was observed in the average relative abundances of these core bacteria over time, with the emergence of four additional core groups of bacteria that corresponded to a decrease in water temperature, possibly caused by season- or foraging-related changes in skin biochemistry that influenced microbial growth, or other temporal factors. The skin microbiome differed between whales sampled at several regional locations along the WAP, suggesting that environmental factors or population may also influence the whale skin microbiome. Overall, the skin microbiome of humpback whales appears to provide insight into animal- and environment-related factors and may serve as a useful indicator for animal health or ecosystem alterations.IMPORTANCE The microbiomes of wild animals are currently understudied but may provide information about animal health and/or animal-environment interactions. In the largest sampling of any marine mammal microbiome, this study demonstrates conservation in the skin microbiome of 89 seemingly healthy humpback whales sampled in the Western Antarctic Peninsula, with shifts in the bacterial community composition related to temporal and regional variability. This study is important because it suggests that the skin microbiome of humpback whales could provide insight into animal nutritional or seasonal/environment-related factors, which are becoming increasingly important to recognize due to unprecedented rates of climate change and anthropogenic impact on ocean ecosystems.

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.

Shaped by Their Environment: Variation in Blue Whale Morphology across Three Productive Coastal Ecosystems.

Species ecology and life history patterns are often reflected in animal morphology. Blue whales are globally distributed, with distinct populations that feed in different productive coastal regions worldwide. Thus, they provide an opportunity to investigate how regional ecosystem characteristics may drive morphological differences within a species. Here, we compare physical and biological oceanography of three different blue whale foraging grounds: (1) Monterey Bay, California, USA; (2) the South Taranaki Bight (STB), Aotearoa New Zealand; and (3) the Corcovado Gulf, Chile. Additionally, we compare the morphology of blue whales from these regions using unoccupied aircraft imagery. Monterey Bay and the Corcovado Gulf are seasonally productive and support the migratory life history strategy of the Eastern North Pacific (ENP) and Chilean blue whale populations, respectively. In contrast, the New Zealand blue whale population remains in the less productive STB year-round. All three populations were indistinguishable in total body length. However, New Zealand blue whales were in significantly higher body condition despite lower regional productivity, potentially attributable to their non-migratory strategy that facilitates lower risk of spatiotemporal misalignment with more consistently available foraging opportunities. Alternatively, the migratory strategy of the ENP and Chilean populations may be successful when their presence on the foraging grounds temporally aligns with abundant prey availability. We document differences in skull and fluke morphology between populations, which may relate to different feeding behaviors adapted to region-specific prey and habitat characteristics. These morphological features may represent a trade-off between maneuverability for prey capture and efficient long-distance migration. As oceanographic patterns shift relative to long-term means under climate change, these blue whale populations may show different vulnerabilities due to differences in migratory phenology and feeding behavior between regions. Spanish abstract La ecología y patrones de historia de vida de las especies a menudo se reflejan en la morfología animal. Las ballenas azules están distribuidas globalmente, con poblaciones separadas que se alimentan en diferentes regiones costeras productivas de todo el mundo. Por lo tanto, brindan la oportunidad de investigar cómo las características regionales de los ecosistemas pueden impulsar diferencias morfológicas dentro de una especie. Aquí, comparamos la oceanografía física y biológica de tres zonas de alimentación diferentes de la ballena azul: (1) Bahía de Monterey, California, EE. UU., (2) Bahía del sur de Taranaki (BST), Nueva Zelanda, y (3) Golfo de Corcovado, Chile. Adicionalmente, comparamos la morfología de las ballenas azules de estas regiones utilizando imágenes de aeronaves no tripuladas. La Bahía de Monterey y el Golfo de Corcovado son estacionalmente productivos y apoyan la estrategia migratoria de la historia de vida de las poblaciones de ballena azul chilena y del Pacífico Norte Oriental (PNO), respectivamente. Por el contrario, la población de ballena azul de Nueva Zelanda permanece en la menos productiva BST durante todo el año. Las tres poblaciones eran indistinguibles en cuanto a la longitud corporal total. Sin embargo, las ballenas azules de Nueva Zelanda tenían una condición corporal significativamente mayor a pesar de una menor productividad regional, potencialmente atribuible a su estrategia no migratoria que facilita un menor riesgo de desalineación espaciotemporal con oportunidades de alimentación disponibles de manera más consistente. Alternativamente, la estrategia migratoria de las poblaciones de ballenas PNO y chilena puede tener éxito cuando su presencia en las zonas de alimentación se alinea temporalmente con la abundante disponibilidad de presas. Documentamos diferencias en la morfología del cráneo y la aleta caudal entre poblaciones, que pueden estar relacionadas con diferentes comportamientos de alimentación adaptados a las características de hábitat y presas específicas para cada región. Estas características morfológicas pueden representar una compensación entre la maniobrabilidad para la captura de presas y una migración eficiente a larga distancia. A medida que los patrones oceanográficos cambian en términos de mediano a largo plazo debido al cambio climático, estas poblaciones de ballenas azules pueden mostrar diferentes vulnerabilidades debido a diferencias en la fenología migratoria y el comportamiento de alimentación entre regiones.

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.

Baleen whale inhalation variability revealed using animal-borne video tags.

Empirical metabolic rate and oxygen consumption estimates for free-ranging whales have been limited to counting respiratory events at the surface. Because these observations were limited and generally viewed from afar, variability in respiratory properties was unknown and oxygen consumption estimates assumed constant breath-to-breath tidal volume and oxygen uptake. However, evidence suggests that cetaceans in human care vary tidal volume and breathing frequency to meet aerobic demand, which would significantly impact energetic estimates if the findings held in free-ranging species. In this study, we used suction cup-attached video tags positioned posterior to the nares of two humpback whales (Megaptera novaeangliae) and four Antarctic minke whales (Balaenoptera bonaerensis) to measure inhalation duration, relative nares expansion, and maximum nares expansion. Inhalation duration and nares expansion varied between and within initial, middle, and terminal breaths of surface sequences between dives. The initial and middle breaths exhibited the least variability and had the shortest durations and smallest nares expansions. In contrast, terminal breaths were highly variable, with the longest inhalation durations and the largest nares expansions. Our results demonstrate breath-to-breath variability in duration and nares expansion, suggesting differential oxygen exchange in each breath during the surface interval. With future validation, inhalation duration or nares area could be used alongside respiratory frequency to improve oxygen consumption estimates by accounting for breath-to-breath variation in wild whales.

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.

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.