Bellwethers of change: population modelling of North Pacific humpback whales from 2002 through 2021 reveals shift from recovery to climate response

For the 40 years after the end of commercial whaling in 1976, humpback whale populations in the North Pacific Ocean exhibited a prolonged period of recovery. Using mark-recapture methods on the largest individual photo-identification dataset ever assembled for a cetacean, we estimated annual ocean-basin-wide abundance for the species from 2002 through 2021. Trends in annual estimates describe strong post-whaling era population recovery from 16 875 (± 5955) in 2002 to a peak abundance estimate of 33 488 (± 4455) in 2012. An apparent 20% decline from 2012 to 2021, 33 488 (± 4455) to 26 662 (± 4192), suggests the population abruptly reached carrying capacity due to loss of prey resources. This was particularly evident for humpback whales wintering in Hawai'i, where, by 2021, estimated abundance had declined by 34% from a peak in 2013, down to abundance levels previously seen in 2006, and contrasted to an absence of decline in Mainland Mexico breeding humpbacks. The strongest marine heatwave recorded globally to date during the 2014-2016 period appeared to have altered the course of species recovery, with enduring effects. Extending this time series will allow humpback whales to serve as an indicator species for the ecosystem in the face of a changing climate.

Boom-bust cycles in gray whales associated with dynamic and changing Arctic conditions

Climate change is affecting a wide range of global systems, with polar ecosystems experiencing the most rapid change. Although climate impacts affect lower-trophic-level and short-lived species most directly, it is less clear how long-lived and mobile species will respond to rapid polar warming because they may have the short-term ability to accommodate ecological disruptions while adapting to new conditions. We found that the population dynamics of an iconic and highly mobile polar-associated species are tightly coupled to Arctic prey availability and access to feeding areas. When low prey biomass coincided with high ice cover, gray whales experienced major mortality events, each reducing the population by 15 to 25%. This suggests that even mobile, long-lived species are sensitive to dynamic and changing conditions as the Arctic warms.

Migratory destinations and spatial structuring of humpback whales (Megaptera novaeangliae) wintering off Nicaragua

Understanding the migratory patterns of large whales is of conservation importance, especially in identifying threats to specific populations. Migration ecology, including migratory destinations, movements and site fidelity for humpback whales (Megaptera novaeangliae) remain poorly studied in parts of the range of the Central America population, considered endangered under the United States Endangered Species Act. This study aimed to investigate the migratory destinations of humpback whales sighted at two study sites in Nicaragua, which are part of the Central America population. A ten-year photographic database of humpback whales observed off Nicaragua was combined with citizen science contributions and sightings from dedicated research programs. The resulting image collection was compared with available historical photo identifications and databases using an automated image recognition algorithm. This approach yielded 36 years of photographic identification totaling 431 recaptures in Nicaragua (2006-2008 and 2016-2021) and 2539 recaptures (1986-2020) in both feeding and breeding grounds of 176 unique individuals sighted in Nicaragua. Our results showed that photo-identified whales were recaptured between October and April in breeding grounds and year-round in feeding grounds between British Columbia and California, with peak recaptures between June and October. Our study provided first-time evidence on fine-scale site affinity of individual humpback whales within Nicaraguan waters and to other breeding and feeding grounds.

A collaborative and near-comprehensive North Pacific humpback whale photo-ID dataset

We present an ocean-basin-scale dataset that includes tail fluke photographic identification (photo-ID) and encounter data for most living individual humpback whales (Megaptera novaeangliae) in the North Pacific Ocean. The dataset was built through a broad collaboration combining 39 separate curated photo-ID catalogs, supplemented with community science data. Data from throughout the North Pacific were aggregated into 13 regions, including six breeding regions, six feeding regions, and one migratory corridor. All images were compared with minimal pre-processing using a recently developed image recognition algorithm based on machine learning through artificial intelligence; this system is capable of rapidly detecting matches between individuals with an estimated 97-99% accuracy. For the 2001-2021 study period, a total of 27,956 unique individuals were documented in 157,350 encounters. Each individual was encountered, on average, in 5.6 sampling periods (i.e., breeding and feeding seasons), with an annual average of 87% of whales encountered in more than one season. The combined dataset and image recognition tool represents a living and accessible resource for collaborative, basin-wide studies of a keystone marine mammal in a time of rapid ecological change.

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.

Influence of Life-History Parameters on Persistent Organic Pollutant Concentrations in Blubber of Eastern North Pacific Gray Whales (Eschrichtius robustus)

Exposure to persistent organic pollutants (POPs) can significantly impact marine mammal health, reproduction, and fitness. This study addresses a significant 20-year gap in gray whale contaminant monitoring through analysis of POPs in 120 blubber biopsies. The scope of this substantial sample set is noteworthy in its range and diversity with collection between 2003 and 2017 along North America's west coast and across diverse sex, age, and reproductive parameters, including paired mothers and calves. Mean blubber concentrations of polychlorinated biphenyls (∑PCBs), dichlorodiphenyltrichloroethanes (∑DDTs), and chlordanes (∑CHLs) generally decreased since previous reports (1968-1999). This is the first report of polybrominated diphenyl ethers (PBDEs) and select hexachlorocyclohexanes (HCHs) in this species. Statistical modeling of the 19 most frequently detected compounds in this dataset revealed sex-, age-, and reproductive status-related patterns, predominantly attributed to maternal offloading. Mean POP concentrations differed significantly by sex in adults (17 compounds, up to 3-fold higher in males) but not in immatures (all 19 compounds). Mean POP concentrations were significantly greater in adults versus immatures in both males (17 compounds, up to 12-fold) and females (13 compounds, up to 3-fold). POP concentrations were detected with compound-specific patterns in nursing calves, confirming maternal offloading for the first time in this species.

Oceanic giants dance to atmospheric rhythms: Ephemeral wind-driven resource tracking by blue whales

Trophic transfer of energy through marine food webs is strongly influenced by prey aggregation and its exploitation by predators. Rapid aggregation of some marine fish and crustacean forage species during wind-driven coastal upwelling has recently been discovered, motivating the hypothesis that predators of these forage species track the upwelling circulation in which prey aggregation occurs. We examine this hypothesis in the central California Current Ecosystem using integrative observations of upwelling dynamics, forage species' aggregation, and blue whale movement. Directional origins of blue whale calls repeatedly tracked upwelling plume circulation when wind-driven upwelling intensified and aggregation of forage species was heightened. Our findings illustrate a resource tracking strategy by which blue whales may maximize energy gain amid ephemeral foraging opportunities. These findings have implications for the ecology and conservation of diverse predators that are sustained by forage populations whose behaviour is responsive to episodic environmental dynamics.

Naturally stressed? Glucocorticoid profiles in blubber of blue and gray whales in response to life history parameters

The goal of the present study was to carry out a thorough methodological validation and describe baseline profiles for glucocorticoid hormones (cortisol and corticosterone) in blubber from blue (n = 77) and gray (n = 103) whales from the eastern North Pacific Ocean. For each species, we modelled cortisol and corticosterone concentrations in response to life history parameters (age, sex, reproductive status) and season or geographic location. In blue whales, cortisol concentrations did not vary significantly by age class, sex, or reproductive status, whereas corticosterone was significantly lower in immature than in adult females (p < .001). In gray whales, cortisol concentrations were significantly higher in lactating whales (p < .05), while corticosterone was significantly different between females and males (p = .001) and elevated in calves (p = .003). In gray whales, corticosterone concentrations were significantly lower in males sampled later in the year (August to November) compared to both sexes sampled between March and August (p = .05), but no seasonal trend occurred in blue whales. Our results indicate that glucocorticoid actions vary between species and sex in large whales. Analysis of multiple hormones improves our understanding of the physiology of maintaining metabolic homeostasis or coping with chronic stressors.

Blue whales increase feeding rates at fine-scale ocean features

Marine predators face the challenge of reliably finding prey that is patchily distributed in space and time. Predators make movement decisions at multiple spatial and temporal scales, yet we have a limited understanding of how habitat selection at multiple scales translates into foraging performance. In the ocean, there is mounting evidence that submesoscale (i.e. less than 100 km) processes drive the formation of dense prey patches that should hypothetically provide feeding hot spots and increase predator foraging success. Here, we integrated environmental remote-sensing with high-resolution animal-borne biologging data to evaluate submesoscale surface current features in relation to the habitat selection and foraging performance of blue whales in the California Current System. Our study revealed a consistent functional relationship in which blue whales disproportionately foraged within dynamic aggregative submesoscale features at both the regional and feeding site scales across seasons, regions and years. Moreover, we found that blue whale feeding rates increased in areas with stronger aggregative features, suggesting that these features indicate areas of higher prey density. The use of fine-scale, dynamic features by foraging blue whales underscores the need to take these features into account when designating critical habitat and may help inform strategies to mitigate the impacts of human activities for the species.

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.

Integrating remote sensing methods during controlled exposure experiments to quantify group responses of dolphins to navy sonar

Human noise can be harmful to sound-centric marine mammals. Significant research has focused on characterizing behavioral responses of protected cetacean species to navy mid-frequency active sonar (MFAS). Controlled exposure experiments (CEE) using animal-borne tags have proved valuable, but smaller dolphins are not amenable to tagging and groups of interacting individuals are more relevant behavioral units for these social species. To fill key data gaps on group responses of social delphinids that are exposed to navy MFAS in large numbers, we describe novel approaches for the coordinated collection and integrated analysis of multiple remotely-sensed datasets during CEEs. This involves real-time coordination of a sonar source, shore-based group tracking, aerial photogrammetry to measure fine-scale movements and passive acoustics to quantify vocal activity. Using an example CEE involving long-beaked common dolphins (Delphinus delphis bairdii), we demonstrate how resultant quantitative metrics can be used to estimate behavioral changes and noise exposure-response relationships.

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.

Epizootiology of a Cryptococcus gattii outbreak in porpoises and dolphins from the Salish Sea

Cryptococcus gattii is a fungal pathogen that primarily affects the respiratory and nervous systems of humans and other animals. C. gattii emerged in temperate North America in 1999 as a multispecies outbreak of cryptococcosis in British Columbia (Canada) and Washington State and Oregon (USA), affecting humans, domestic animals, and wildlife. Here we describe the C. gattii epizootic in odontocetes. Cases of C. gattii were identified in 42 odontocetes in Washington and British Columbia between 1997 and 2016. Species affected included harbor porpoises Phocoena phocoena (n = 26), Dall's porpoises Phocoenoides dalli (n = 14), and Pacific white-sided dolphins Lagenorhynchus obliquidens (n = 2). The probable index case was identified in an adult male Dall's porpoise in 1997, 2 yr prior to the initial terrestrial outbreak. The spatiotemporal extent of the C. gattii epizootic was defined, and cases in odontocetes were found to be clustered around terrestrial C. gattii hotspots. Case-control analyses with stranded, uninfected odontocetes revealed that risk factors for infection were species (Dall's porpoises), age class (adult animals), and season (winter). This study suggests that mycoses are an emerging source of mortality for odontocetes, and that outbreaks may be associated with anthropogenic environmental disturbance.

Blubber endocrine profiles provide insights into reproductive biology in blue whales from the eastern North Pacific Ocean

The goal of the present study was to complement existing data of testosterone and progesterone in blue whale (Balaenoptera musculus) blubber from the eastern North Pacific Ocean to evaluate effects of seasonality and location on these hormones and to better assess reproductive status of individuals. Physiological parameters regarding reproduction are fundamental for describing population dynamics, and hormones can be a valid tool to estimate those for wildlife populations. In this study, blubber tissue was validated for testosterone and progesterone assays. Hormone concentrations were measured in 69 (35 males and 34 females) blubber samples from live (n = 66) and stranded (n = 3) animals collected between 2002 and 2016 from a known winter reproductive ground in the Gulf of California (GoC) and summer feeding areas along the United States West Coast (USWC), specifically off the states of California and Oregon. Results were combined with sighting histories as a tool to determine reproductive status of individual whales. Testosterone concentrations in adult male blue whales were significantly higher (p < 0.05) in blubber biopsies sampled off the USWC between the months of June and November compared to those sampled in the GoC between February and April. Elevated testosterone concentrations likely indicate physiological preparation for reproductive activity while the animals were present off the USWC. Progesterone concentrations were significantly elevated in pregnant females, confirming progesterone as an indicator of pregnancy in blue whales. Probabilities of being pregnant were estimated for adult females with unknown sighting histories based on progesterone concentrations. Testosterone in females was detected and measured only in pregnant whales suggesting its biosynthesis or metabolism is altered during gestation. These results provide updated and new information on the reproductive cycle of blue whales in the eastern North Pacific, posing new milestones to better estimate the timing of the mating season for this endangered population.

Application of endocrine biomarkers to update information on reproductive physiology in gray whale (Eschrichtius robustus)

Most of our knowledge on reproductive biology of gray whales dates back to scientific research conducted during commercial whaling in the late 1950s and 1960s. The goal of the present study was to provide updated insights on reproductive physiology of gray whales, using progesterone and testosterone as biomarkers. We measured hormone concentrations using enzyme immunoassay (EIA) techniques in blubber biopsies collected from 106 individual whales from March to November over a span of 12 years (2004-2016) between California and Alaska. We found testosterone concentrations in males to increase significantly with age (P = 0.03). Adult males showed significantly elevated testosterone concentrations when sampled in the fall compared to the summer (P = 0.01), likely indicating physiological preparation for mating. We measured testosterone concentrations in females of different age classes, but no statistical differences were found. We found significantly higher progesterone concentrations in pregnant females compared to non-pregnant females and adult males (P< 0.001), indicating progesterone is a valid biomarker for pregnancy in gray whales. Both female and male calves had elevated progesterone concentrations, suggesting maternal transfer via lactation. We fit a mixture of two normal distributions to progesterone data from all non-calf females to identify clusters of high and low progesterone and estimated the probability of being pregnant for whales of unknown reproductive status. With this approach we identified likely pregnant and non-pregnant animals. This study represents an important milestone on reproductive profiles in this population, that can be used to estimate more accurate and precise reproductive parameters to be used for better understanding population dynamics of gray whales.

Population structure in a continuously distributed coastal marine species, the harbor porpoise, based on microhaplotypes derived from poor-quality samples

Harbor porpoise in the North Pacific are found in coastal waters from southern California to Japan, but population structure is poorly known outside of a few local areas. We used multiplexed amplicon sequencing of 292 loci and genotyped clusters of single nucleotide polymoirphisms as microhaplotypes (N = 271 samples) in addition to mitochondrial (mtDNA) sequence data (N = 413 samples) to examine the genetic structure from samples collected along the Pacific coast and inland waterways from California to southern British Columbia. We confirmed an overall pattern of strong isolation-by-distance, suggesting that individual dispersal is restricted. We also found evidence of regions where genetic differences are larger than expected based on geographical distance alone, implying current or historical barriers to gene flow. In particular, the southernmost population in California is genetically distinct (F_ST  = 0.02 [microhaplotypes]; 0.31 [mtDNA]), with both reduced genetic variability and high frequency of an otherwise rare mtDNA haplotype. At the northern end of our study range, we found significant genetic differentiation of samples from the Strait of Georgia, previously identified as a potential biogeographical boundary or secondary contact zone between harbor porpoise populations. Association of microhaplotypes with remotely sensed environmental variables indicated potential local adaptation, especially at the southern end of the species' range. These results inform conservation and management for this nearshore species, illustrate the value of genomic methods for detecting patterns of genetic structure within a continuously distributed marine species, and highlight the power of microhaplotype genotyping for detecting genetic structure in harbor porpoises despite reliance on poor-quality samples.

Context-dependent variability in the predicted daily energetic costs of disturbance for blue whales

Assessing the long-term consequences of sub-lethal anthropogenic disturbance on wildlife populations requires integrating data on fine-scale individual behavior and physiology into spatially and temporally broader, population-level inference. A typical behavioral response to disturbance is the cessation of foraging, which can be translated into a common metric of energetic cost. However, this necessitates detailed empirical information on baseline movements, activity budgets, feeding rates and energy intake, as well as the probability of an individual responding to the disturbance-inducing stressor within different exposure contexts. Here, we integrated data from blue whales (Balaenoptera musculus) experimentally exposed to military active sonar signals with fine-scale measurements of baseline behavior over multiple days or weeks obtained from accelerometry loggers, telemetry tracking and prey sampling. Specifically, we developed daily simulations of movement, feeding behavior and exposure to localized sonar events of increasing duration and intensity and predicted the effects of this disturbance source on the daily energy intake of an individual. Activity budgets and movements were highly variable in space and time and among individuals, resulting in large variability in predicted energetic intake and costs. In half of our simulations, an individual's energy intake was unaffected by the simulated source. However, some individuals lost their entire daily energy intake under brief or weak exposure scenarios. Given this large variation, population-level models will have to assess the consequences of the entire distribution of energetic costs, rather than only consider single summary statistics. The shape of the exposure-response functions also strongly influenced predictions, reinforcing the need for contextually explicit experiments and improved mechanistic understanding of the processes driving behavioral and physiological responses to disturbance. This study presents a robust approach for integrating different types of empirical information to assess the effects of disturbance at spatio-temporal and ecological scales that are relevant to management and conservation.

Isotope-based inferences of the seasonal foraging and migratory strategies of blue whales in the eastern Pacific Ocean

Migratory marine megafauna generally move vast distances between productive foraging grounds and environmentally stable breeding grounds, but characterizing how they use these habitats to maintain homeostasis and reproduce is difficult. We used isotope analysis of blue whale skin strata (n = 621) and potential prey (n = 300) to examine their migratory and foraging strategies in the eastern Pacific Ocean. Our results suggest that most whales in the northeast Pacific use a mixed income and capital breeding strategy, and use the California Current Ecosystem as their primary summer-fall foraging ground. A subset of individuals exhibited migratory plasticity and spend most of the year in the Gulf of California or Costa Rica Dome, two regions believed to be their primary winter-spring breeding grounds. Isotope data also revealed that whales in the southern Eastern Tropical Pacific generally do not forage in the northeast Pacific, which suggests a north-south population structure with a boundary near the equator.

Remoras pick where they stick on blue whales

Animal-borne video recordings from blue whales in the open ocean show that remoras preferentially adhere to specific regions on the surface of the whale. Using empirical and computational fluid dynamics analyses, we show that remora attachment was specific to regions of separating flow and wakes caused by surface features on the whale. Adhesion at these locations offers remoras drag reduction of up to 71-84% compared with the freestream. Remoras were observed to move freely along the surface of the whale using skimming and sliding behaviors. Skimming provided drag reduction as high as 50-72% at some locations for some remora sizes, but little to none was available in regions where few to no remoras were observed. Experimental work suggests that the Venturi effect may help remoras stay near the whale while skimming. Understanding the flow environment around a swimming blue whale will inform the placement of biosensor tags to increase attachment time for extended ecological monitoring.

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.

Revision of fin whale Balaenoptera physalus (Linnaeus, 1758) subspecies using genetics

Three subspecies of fin whales (Balaenoptera physalus) are currently recognized, including the northern fin whale (B. p. physalus), the southern fin whale (B. p. quoyi), and the pygmy fin whale (B. p. patachonica). The Northern Hemisphere subspecies encompasses fin whales in both the North Atlantic and North Pacific oceans. A recent analysis of 154 mitogenome sequences of fin whales from these two ocean basins and the Southern Hemisphere suggested that the North Pacific and North Atlantic populations should be treated as different subspecies. Using these mitogenome sequences, in this study, we conduct analyses on a larger mtDNA control region data set, and on 23 single-nucleotide polymorphisms (SNPs) from 144 of the 154 samples in the mitogenome data set. Our results reveal that North Pacific and North Atlantic fin whales can be correctly assigned to their ocean basin with 99% accuracy. Results of the SNP analysis indicate a correct classification rate of 95%, very low rates of gene flow among ocean basins, and that distinct mitogenome matrilines in the North Pacific are interbreeding. These results indicate that North Pacific fin whales should be recognized as a separate subspecies, with the name B. p. velifera Cope in Scammon 1869 as the oldest available name.

Body Flexibility Enhances Maneuverability in the World's Largest Predator

Blue whales are often characterized as highly stable, open-ocean swimmers who sacrifice maneuverability for long-distance cruising performance. However, recent studies have revealed that blue whales actually exhibit surprisingly complex underwater behaviors, yet little is known about the performance and control of these maneuvers. Here, we use multi-sensor biologgers equipped with cameras to quantify the locomotor dynamics and the movement of the control surfaces used by foraging blue whales. Our results revealed that simple maneuvers (rolls, turns, and pitch changes) are performed using distinct combinations of control and power provided by the flippers, the flukes, and bending of the body, while complex trajectories are structured by combining sequences of simple maneuvers. Furthermore, blue whales improve their turning performance by using complex banked turns to take advantage of their substantial dorso-ventral flexibility. These results illustrate the important role body flexibility plays in enhancing control and performance of maneuvers, even in the largest of animals. The use of the body to supplement the performance of the hydrodynamically active surfaces may represent a new mechanism in the control of aquatic locomotion.

Corrigendum to: Heavy with child? Pregnancy status and stable isotope ratios as determined from biopsies of humpback whales

[This corrects the article DOI: 10.1093/conphys/cow050.].

Behavioral responses of individual blue whales (Balaenoptera musculus) to mid-frequency military sonar

This study measured the degree of behavioral responses in blue whales (Balaenoptera musculus) to controlled noise exposure off the southern California coast. High-resolution movement and passive acoustic data were obtained from non-invasive archival tags (n=42) whereas surface positions were obtained with visual focal follows. Controlled exposure experiments (CEEs) were used to obtain direct behavioral measurements before, during and after simulated and operational military mid-frequency active sonar (MFAS), pseudorandom noise (PRN) and controls (no noise exposure). For a subset of deep-feeding animals (n=21), active acoustic measurements of prey were obtained and used as contextual covariates in response analyses. To investigate potential behavioral changes within individuals as a function of controlled noise exposure conditions, two parallel analyses of time-series data for selected behavioral parameters (e.g. diving, horizontal movement and feeding) were conducted. This included expert scoring of responses according to a specified behavioral severity rating paradigm and quantitative change-point analyses using Mahalanobis distance statistics. Both methods identified clear changes in some conditions. More than 50% of blue whales in deep-feeding states responded during CEEs, whereas no changes in behavior were identified in shallow-feeding blue whales. Overall, responses were generally brief, of low to moderate severity, and highly dependent on exposure context such as behavioral state, source-to-whale horizontal range and prey availability. Response probability did not follow a simple exposure-response model based on received exposure level. These results, in combination with additional analytical methods to investigate different aspects of potential responses within and among individuals, provide a comprehensive evaluation of how free-ranging blue whales responded to mid-frequency military sonar.

Context-dependent variability in blue whale acoustic behaviour

Acoustic communication is an important aspect of reproductive, foraging and social behaviours for many marine species. Northeast Pacific blue whales (Balaenoptera musculus) produce three different call types-A, B and D calls. All may be produced as singular calls, but A and B calls also occur in phrases to form songs. To evaluate the behavioural context of singular call and phrase production in blue whales, the acoustic and dive profile data from tags deployed on individuals off southern California were assessed using generalized estimating equations. Only 22% of all deployments contained sounds attributed to the tagged animal. A larger proportion of tagged animals were female (47%) than male (13%), with 40% of unknown sex. Fifty per cent of tags deployed on males contained sounds attributed to the tagged whale, while only a few (5%) deployed on females did. Most calls were produced at shallow depths (less than 30 m). Repetitive phrasing (singing) and production of singular calls were most common during shallow, non-lunging dives, with the latter also common during surface behaviour. Higher sound production rates occurred during autumn than summer and they varied with time-of-day: singular call rates were higher at dawn and dusk, while phrase production rates were highest at dusk and night.

Conception, fetal growth, and calving seasonality of harbor porpoise (Phocoena phocoena) in the Salish Sea waters of Washington, USA, and southern British Columbia, Canada

We evaluated harbor porpoise (Phocoena phocoena (Linnaeus, 1758)) strandings in the Salish Sea to determine calving seasonality (1980–2015). A total of 443 strandings were analyzed, of which 134 were calves and 53 were neonates. Stranded calves were reported every month, but peaked in July, August, and September. Based on fetal size and an estimated fetal growth rate of 80 mm/month, mean (±SD) conception date (and range) was back-calculated to 11 October ± 30 days (16 August – 31 December) and was later than in most other studies. Using mean (±SD) length at birth (80 ± 5.8 cm), gestation was estimated to be approximately 10.8 months. Estimated birthing period was 16 July – 27 November, with a mean (±SD) birth date of 10 September (±30.7 days) and a birth length of 80.0 cm. Estimated pregnancy rate (0.28–0.29) is lower than reported in other areas and is likely an underestimate due to missed early embryos, poor postmortem condition of a large proportion of the stranded adult females, and potential biases related to the animals that strand and are available. This study of harbor porpoise reproduction and calving in the Salish Sea is the first assessment of calving seasonality for this species in the northeast Pacific Ocean.

Using Digital Tags With Integrated Video and Inertial Sensors to Study Moving Morphology and Associated Function in Large Aquatic Vertebrates

The anatomy of large cetaceans has been well documented, mostly through dissection of dead specimens. However, the difficulty of studying the world's largest animals in their natural environment means the functions of anatomical structures must be inferred. Recently, non-invasive tracking devices have been developed that measure body position and orientation, thereby enabling the detailed reconstruction of underwater trajectories. The addition of cameras to the whale-borne tags allows the sensor data to be matched with real-time observations of how whales use their morphological structures, such as flukes, flippers, feeding apparatuses, and blowholes for the physiological functions of locomotion, feeding, and breathing. Here, we describe a new tag design with integrated video and inertial sensors and how it can be used to provide insights to the function of whale anatomy. This technology has the potential to facilitate a wide range of discoveries and comparative studies, but many challenges remain to increase the resolution and applicability of the data. Anat Rec, 300:1935-1941, 2017. © 2017 Wiley Periodicals, Inc.

Risso's dolphins plan foraging dives

Humans remember the past and use that information to plan future actions. Lab experiments that test memory for the location of food show that animals have a similar capability to act in anticipation of future needs, but less work has been done on animals foraging in the wild. We hypothesized that planning abilities are critical and common in breath-hold divers who adjust each dive to forage on prey varying in quality, location and predictability within constraints of limited oxygen availability. We equipped Risso's dolphins with sound-and-motion recording tags to reveal where they focus their attention through their externally observable echolocation and how they fine tune search strategies in response to expected and observed prey distribution. The information from the dolphins was integrated with synoptic prey data obtained from echosounders on an underwater vehicle. At the start of the dives, whales adjusted their echolocation inspection ranges in ways that suggest planning to forage at a particular depth. Once entering a productive prey layer, dolphins reduced their search range comparable to the scale of patches within the layer, suggesting that they were using echolocation to select prey within the patch. On ascent, their search range increased, indicating that they decided to stop foraging within that layer and started searching for prey in shallower layers. Information about prey, learned throughout the dive, was used to plan foraging in the next dive. Our results demonstrate that planning for future dives is modulated by spatial memory derived from multi-modal prey sampling (echoic, visual and capture) during earlier dives.

Resource partitioning facilitates coexistence in sympatric cetaceans in the California Current

Resource partitioning is an important process driving habitat use and foraging strategies in sympatric species that potentially compete. Differences in foraging behavior are hypothesized to contribute to species coexistence by facilitating resource partitioning, but little is known on the multiple mechanisms for partitioning that may occur simultaneously. Studies are further limited in the marine environment, where the spatial and temporal distribution of resources is highly dynamic and subsequently difficult to quantify. We investigated potential pathways by which foraging behavior may facilitate resource partitioning in two of the largest co-occurring and closely related species on Earth, blue (Balaenoptera musculus) and humpback (Megaptera novaeangliae) whales. We integrated multiple long-term datasets (line-transect surveys, whale-watching records, net sampling, stable isotope analysis, and remote-sensing of oceanographic parameters) to compare the diet, phenology, and distribution of the two species during their foraging periods in the highly productive waters of Monterey Bay, California, USA within the California Current Ecosystem. Our long-term study reveals that blue and humpback whales likely facilitate sympatry by partitioning their foraging along three axes: trophic, temporal, and spatial. Blue whales were specialists foraging on krill, predictably targeting a seasonal peak in krill abundance, were present in the bay for an average of 4.7 months, and were spatially restricted at the continental shelf break. In contrast, humpback whales were generalists apparently feeding on a mixed diet of krill and fishes depending on relative abundances, were present in the bay for a more extended period (average of 6.6 months), and had a broader spatial distribution at the shelf break and inshore. Ultimately, competition for common resources can lead to behavioral, morphological, and physiological character displacement between sympatric species. Understanding the mechanisms for species coexistence is both fundamental to maintaining biodiverse ecosystems, and provides insight into the evolutionary drivers of morphological differences in closely related species.

High mortality of blue, humpback and fin whales from modeling of vessel collisions on the U.S. West Coast suggests population impacts and insufficient protection

Mortality from collisions with vessels is one of the main human causes of death for large whales. Ship strikes are rarely witnessed and the distribution of strike risk and estimates of mortality remain uncertain at best. We estimated ship strike mortality for blue humpback and fin whales in U.S. West Coast waters using a novel application of a naval encounter model. Mortality estimates from the model were far higher than current minimum estimates derived from stranding records and are closer to extrapolations adjusted for detection probabilities of dead whales. Our most conservative model estimated mortality to be 7.8x, 2.0x and 2.7x the U.S. recommended limit for blue, humpback and fin whales, respectively, suggesting that death from vessel collisions may be a significant impediment to population growth and recovery. Comparing across the study area, the majority of strike mortality occurs in waters off California, from Bodega Bay south and tends to be concentrated in a band approximately 24 Nm (44.5 km) offshore and in designated shipping lanes leading to and from major ports. While some mortality risk exists across nearly all West Coast waters, 74%, 82% and 65% of blue, humpback and fin whale mortality, respectively, occurs in just 10% of the study area, suggesting conservation efforts can be very effective if focused in these waters. Risk is highest in the shipping lanes off San Francisco and Long Beach, but only a fraction of total estimated mortality occurs in these proportionally small areas, making any conservation efforts exclusively within these areas insufficient to address overall strike mortality. We recommend combining shipping lane modifications and re-locations, ship speed reductions and creation of ‘Areas to be Avoided’ by vessels in ecologically important locations to address this significant source of whale mortality.

Avoidance responses of minke whales to 1-4kHz naval sonar

Minke whales are difficult to study and little information exists regarding their responses to anthropogenic sound. This study pools data from behavioural response studies off California and Norway. Data are derived from four tagged animals, of which one from each location was exposed to naval sonar signals. Statistical analyses were conducted using Mahalanobis distance to compare overall changes in parameters summarising dive behaviour, avoidance behaviour, and potential energetic costs of disturbance. Our quantitative analysis showed that both animals initiated avoidance behaviour, but responses were not associated with unusual dive behaviour. In one exposed animal the avoidance of the sonar source included a 5-fold increase in horizontal speed away from the source, implying a significant increase in metabolic rate. Despite the different environmental settings and exposure contexts, clear changes in behaviour were observed providing the first insights into the nature of responses to human noise for this wide-ranging species.

Erratum: Heavy with child? Pregnancy status and stable isotope ratios as determined from biopsies of humpback whales

[This corrects the article DOI: 10.1093/conphys/cow050.].

Estimating blue whale skin isotopic incorporation rates and baleen growth rates: Implications for assessing diet and movement patterns in mysticetes

Stable isotope analysis in mysticete skin and baleen plates has been repeatedly used to assess diet and movement patterns. Accurate interpretation of isotope data depends on understanding isotopic incorporation rates for metabolically active tissues and growth rates for metabolically inert tissues. The aim of this research was to estimate isotopic incorporation rates in blue whale skin and baleen growth rates by using natural gradients in baseline isotope values between oceanic regions. Nitrogen (δ¹⁵N) and carbon (δ¹³C) isotope values of blue whale skin and potential prey were analyzed from three foraging zones (Gulf of California, California Current System, and Costa Rica Dome) in the northeast Pacific from 1996–2015. We also measured δ¹⁵N and δ¹³C values along the lengths of baleen plates collected from six blue whales stranded in the 1980s and 2000s. Skin was separated into three strata: basale, externum, and sloughed skin. A mean (±SD) skin isotopic incorporation rate of 163±91 days was estimated by fitting a generalized additive model of the seasonal trend in δ¹⁵N values of skin strata collected in the Gulf of California and the California Current System. A mean (±SD) baleen growth rate of 15.5±2.2 cm y^-1 was estimated by using seasonal oscillations in δ¹⁵N values from three whales. These oscillations also showed that individual whales have a high fidelity to distinct foraging zones in the northeast Pacific across years. The absence of oscillations in δ¹⁵N values of baleen sub-samples from three male whales suggests these individuals remained within a specific zone for several years prior to death. δ¹³C values of both whale tissues (skin and baleen) and potential prey were not distinct among foraging zones. Our results highlight the importance of considering tissue isotopic incorporation and growth rates when studying migratory mysticetes and provide new insights into the individual movement strategies of blue whales.

How Baleen Whales Feed: The Biomechanics of Engulfment and Filtration

Baleen whales are gigantic obligate filter feeders that exploit aggregations of small-bodied prey in littoral, epipelagic, and mesopelagic ecosystems. At the extreme of maximum body size observed among mammals, baleen whales exhibit a unique combination of high overall energetic demands and low mass-specific metabolic rates. As a result, most baleen whale species have evolved filter-feeding mechanisms and foraging strategies that take advantage of seasonally abundant yet patchily and ephemerally distributed prey resources. New methodologies consisting of multi-sensor tags, active acoustic prey mapping, and hydrodynamic modeling have revolutionized our ability to study the physiology and ecology of baleen whale feeding mechanisms. Here, we review the current state of the field by exploring several hypotheses that aim to explain how baleen whales feed. Despite significant advances, major questions remain about the processes that underlie these extreme feeding mechanisms, which enabled the evolution of the largest animals of all time.

Determining forward speed from accelerometer jiggle in aquatic environments

How fast animals move is critical to understanding their energetic requirements, locomotor capacity, and foraging performance, yet current methods for measuring speed via animal-attached devices are not universally applicable. Here we present and evaluate a new method that relates forward speed to the stochastic motion of biologging devices since tag jiggle, the amplitude of the tag vibrations as measured by high sample rate accelerometers, increases exponentially with increasing speed. We successfully tested this method in a flow tank using two types of biologging devices and tested the method in situ on wild cetaceans spanning ∼3 to &gt;20 m in length using two types of suction cup-attached and two types of dart-attached tag. This technique provides some advantages over other approaches for determining speed as it is device-orientation independent and relies only on a pressure sensor and a high sample rate accelerometer, sensors that are nearly universal across biologging device types.

Development of an automated method of detecting stereotyped feeding events in multisensor data from tagged rorqual whales

The introduction of animal‐borne, multisensor tags has opened up many opportunities for ecological research, making previously inaccessible species and behaviors observable. The advancement of tag technology and the increasingly widespread use of bio‐logging tags are leading to large volumes of sometimes extremely detailed data. With the increasing quantity and duration of tag deployments, a set of tools needs to be developed to aid in facilitating and standardizing the analysis of movement sensor data. Here, we developed an observation‐based decision tree method to detect feeding events in data from multisensor movement tags attached to fin whales (Balaenoptera physalus). Fin whales exhibit an energetically costly and kinematically complex foraging behavior called lunge feeding, an intermittent ram filtration mechanism. Using this automated system, we identified feeding lunges in 19 fin whales tagged with multisensor tags, during a total of over 100 h of continuously sampled data. Using movement sensor and hydrophone data, the automated lunge detector correctly identified an average of 92.8% of all lunges, with a false‐positive rate of 9.5%. The strong performance of our automated feeding detector demonstrates an effective, straightforward method of activity identification in animal‐borne movement tag data. Our method employs a detection algorithm that utilizes a hierarchy of simple thresholds based on knowledge of observed features of feeding behavior, a technique that is readily modifiable to fit a variety of species and behaviors. Using automated methods to detect behavioral events in tag records will significantly decrease data analysis time and aid in standardizing analysis methods, crucial objectives with the rapidly increasing quantity and variety of on‐animal tag data. Furthermore, our results have implications for next‐generation tag design, especially long‐term tags that can be outfitted with on‐board processing algorithms that automatically detect kinematic events and transmit ethograms via acoustic or satellite telemetry.

Discrimination of fast click-series produced by tagged Risso's dolphins (Grampus griseus) for echolocation or communication

Early studies that categorized odontocete pulsed sounds had few means of discriminating signals used for biosonar-based foraging from those used for communication. This capability to identify the function of sounds is important for understanding and interpreting behavior; it is also essential for monitoring and mitigating potential disturbance from human activities. Archival tags were placed on free-ranging Grampus griseus to quantify and discriminate between pulsed sounds used for echolocation-based foraging and those used for communication. Two types of rapid click-series pulsed sounds, buzzes and burst pulses, were identified as produced by the tagged dolphins and classified using a Gaussian mixture model based on their duration, association with jerk (i.e. rapid change of acceleration) and temporal association with click trains. Buzzes followed regular echolocation clicks and coincided with a strong jerk signal from accelerometers on the tag. They consisted of series averaging 359±210 clicks (mean±s.d.) with an increasing repetition rate and relatively low amplitude. Burst pulses consisted of relatively short click series averaging 45±54 clicks with decreasing repetition rate and longer inter-click interval that were less likely to be associated with regular echolocation and the jerk signal. These results suggest that the longer, relatively lower amplitude, jerk-associated buzzes are used in this species to capture prey, mostly during the bottom phase of foraging dives, as seen in other odontocetes. In contrast, the shorter, isolated burst pulses that are generally emitted by the dolphins while at or near the surface are used outside of a direct, known foraging context.

Prey-mediated behavioral responses of feeding blue whales in controlled sound exposure experiments

Behavioral response studies provide significant insights into the nature, magnitude, and consequences of changes in animal behavior in response to some external stimulus. Controlled exposure experiments (CEEs) to study behavioral response have faced challenges in quantifying the importance of and interaction among individual variability, exposure conditions, and environmental covariates. To investigate these complex parameters relative to blue whale behavior and how it may change as a function of certain sounds, we deployed multi-sensor acoustic tags and conducted CEEs using simulated mid-frequency active sonar (MFAS) and pseudo-random noise (PRN) stimuli, while collecting synoptic, quantitative prey measures. In contrast to previous approaches that lacked such prey data, our integrated approach explained substantially more variance in blue whale dive behavioral responses to mid-frequency sounds (r2 = 0.725 vs. 0.14 previously). Results demonstrate that deep-feeding whales respond more clearly and strongly to CEEs than those in other behavioral states, but this was only evident with the increased explanatory power provided by incorporating prey density and distribution as contextual covariates. Including contextual variables increases the ability to characterize behavioral variability and empirically strengthens previous findings that deep-feeding blue whales respond significantly to mid-frequency sound exposure. However, our results are only based on a single behavioral state with a limited sample size, and this analytical framework should be applied broadly across behavioral states. The increased capability to describe and account for individual response variability by including environmental variables, such as prey, that drive foraging behavior underscores the importance of integrating these and other relevant contextual parameters in experimental designs. Our results suggest the need to measure and account for the ecological dynamics of predator-prey interactions when studying the effects of anthropogenic disturbance in feeding animals.

Humpback whale diets respond to variance in ocean climate and ecosystem conditions in the California Current

Large, migratory predators are often cited as sentinel species for ecosystem processes and climate-related changes, but their utility as indicators is dependent upon an understanding of their response to environmental variability. Documentation of the links between climate variability, ecosystem change and predator dynamics is absent for most top predators. Identifying species that may be useful indicators and elucidating these mechanistic links provides insight into current ecological dynamics and may inform predictions of future ecosystem responses to climatic change. We examine humpback whale response to environmental variability through stable isotope analysis of diet over a dynamic 20-year period (1993-2012) in the California Current System (CCS). Humpback whale diets captured two major shifts in oceanographic and ecological conditions in the CCS. Isotopic signatures reflect a diet dominated by krill during periods characterized by positive phases of the North Pacific Gyre Oscillation (NPGO), cool sea surface temperature (SST), strong upwelling and high krill biomass. In contrast, humpback whale diets are dominated by schooling fish when the NPGO is negative, SST is warmer, seasonal upwelling is delayed and anchovy and sardine populations display increased biomass and range expansion. These findings demonstrate that humpback whales trophically respond to ecosystem shifts, and as a result, their foraging behavior is a synoptic indicator of oceanographic and ecological conditions across the CCS. Multi-decadal examination of these sentinel species thus provides insight into biological consequences of interannual climate fluctuations, fundamental to advancing ecosystem predictions related to global climate change.

Heavy with child? Pregnancy status and stable isotope ratios as determined from biopsies of humpback whales

Understanding reproductive rates of wild animal populations is crucially important for management and conservation. Assessing pregnancy status of free-ranging cetaceans has historically been difficult; however, recent advances in analytical techniques have allowed the diagnosis of pregnancy from small samples of blubber tissue. The primary objectives of this study were as follows: (i) to test the efficacy of blubber progesterone assays as a tool for diagnosing pregnancy in humpback whales (Megaptera novaeangliae); (ii) to estimate the pregnancy rate of humpback whales in Monterey Bay, California; and (iii) to investigate the relationship between stable isotopes and reproductive status of these whales. Progesterone concentrations of female whales fell into two distinct groups, allowing for diagnostic separation of pregnant and non-pregnant individuals. Pregnancy rate varied between years of the study (48.4%% in 2011 and 18.5% in 2012), but fell within the range of other estimates of reproductive success for this population. Stable carbon and nitrogen isotope ratios were examined to investigate the impacts of pregnancy on these values. Neither δ¹⁵N nor δ¹³C varied in a consistent way among animals of different sex or reproductive status. The relationship between δ¹⁵N and δ¹³C was strongly positive for male and non-pregnant female humpbacks; however, no relationship existed for pregnant whales. This difference may be indicative of the effects of pregnancy on δ¹⁵N, resulting from tissue synthesis and reduced excretion of nitrogenous waste, as well as on δ¹³C through increased mobilization of lipid stores to meet the energetic demands of pregnancy. Ultimately, our results support the use of blubber progesterone assays for diagnosing pregnancy in humpback whales and indicate that, when paired with other approaches (e.g. stable isotope analysis), pregnancy status can be an informative tool for addressing questions about animal physiology, ecology and population biology. This information will provide for more effective management and conservation efforts in a rapidly changing world.

Hydrodynamic properties of fin whale flippers predict maximum rolling performance

Maneuverability is one of the most important and least understood aspects of animal locomotion. The hydrofoil-like flippers of cetaceans are thought to function as control surfaces that effect maneuvers, but quantitative tests of this hypothesis have been lacking. Here we construct a simple hydrodynamic model to predict the longitudinal-axis roll performance of fin whales, and we test its predictions against kinematic data recorded by on-board movement sensors from 27 free-swimming fin whales. We found that for a given swimming speed and roll excursion, the roll velocity of fin whales calculated from our field data agrees well with that predicted by our hydrodynamic model. Although fluke and body torsion may further influence performance, our results indicate that lift generated by the flippers is sufficient to drive most of the longitudinal-axis rolls used by fin whales for feeding and maneuvering.