Vulnerability to climate change of United States marine mammal stocks in the western North Atlantic, Gulf of Mexico, and Caribbean.

Climate change and climate variability are affecting marine mammal species and these impacts are projected to continue in the coming decades. Vulnerability assessments provide a framework for evaluating climate impacts over a broad range of species using currently available information. We conducted a trait-based climate vulnerability assessment using expert elicitation for 108 marine mammal stocks and stock groups in the western North Atlantic, Gulf of Mexico, and Caribbean Sea. Our approach combined the exposure (projected change in environmental conditions) and sensitivity (ability to tolerate and adapt to changing conditions) of marine mammal stocks to estimate vulnerability to climate change, and categorize stocks with a vulnerability index. The climate vulnerability score was very high for 44% (n = 47) of these stocks, high for 29% (n = 31), moderate for 20% (n = 22), and low for 7% (n = 8). The majority of stocks (n = 78; 72%) scored very high exposure, whereas 24% (n = 26) scored high, and 4% (n = 4) scored moderate. The sensitivity score was very high for 33% (n = 36) of these stocks, high for 18% (n = 19), moderate for 34% (n = 37), and low for 15% (n = 16). Vulnerability results were summarized for stocks in five taxonomic groups: pinnipeds (n = 4; 25% high, 75% moderate), mysticetes (n = 7; 29% very high, 57% high, 14% moderate), ziphiids (n = 8; 13% very high, 50% high, 38% moderate), delphinids (n = 84; 52% very high, 23% high, 15% moderate, 10% low), and other odontocetes (n = 5; 60% high, 40% moderate). Factors including temperature, ocean pH, and dissolved oxygen were the primary drivers of high climate exposure, with effects mediated through prey and habitat parameters. We quantified sources of uncertainty by bootstrapping vulnerability scores, conducting leave-one-out analyses of individual attributes and individual scorers, and through scoring data quality for each attribute. These results provide information for researchers, managers, and the public on marine mammal responses to climate change to enhance the development of more effective marine mammal management, restoration, and conservation activities that address current and future environmental variation and biological responses due to climate change.

Quantification by droplet digital PCR and species identification by metabarcoding of environmental (e)DNA from Blainville's beaked whales, with assisted localization from an acoustic array.

Detection and identification of species, subspecies or stocks of whales, dolphins and porpoises at sea remain challenging, particularly for cryptic or elusive species like beaked whales (Family: Ziphiidae). Here we investigated the potential for using an acoustically assisted sampling design to collect environmental (e)DNA from beaked whales on the U.S. Navy's Atlantic Undersea Test and Evaluation Center (AUTEC) in The Bahamas. During 12 days of August 2019, we conducted 9 small-boat surveys and collected 56 samples of seawater (paired subsamples of 1L each, including controls) using both a spatial collection design in the absence of visual confirmation of whales, and a serial collection design in the proximity of whales at the surface. There were 7 sightings of whales, including 11 Blainville's beaked whales (Mesoplodon densirostris). All whales were located initially with the assistance of information from a bottom-mounted acoustic array available on the AUTEC range. Quantification by droplet digital (dd)PCR from the four spatial design collections showed no samples of eDNA above the threshold of detection and none of these 20 samples yielded amplicons for conventional or next-generation sequencing. Quantification of the 31 samples from four serial collections identified 11 likely positive detections. eDNA barcoding by conventional sequencing and eDNA metabarcoding by next-generation sequencing confirmed species identification for 9 samples from three of the four serial collections. We further resolved five intra-specific variants (i.e., haplotypes), two of which showed an exact match to previously published haplotypes and three that have not been reported previously to the international repository, GenBank. A minimum spanning network of the five eDNA haplotypes, with all other published haplotypes of Blainville's beaked whales, suggested the potential for further resolution of differences between oceanic populations.

Using individual-based bioenergetic models to predict the aggregate effects of disturbance on populations: A case study with beaked whales and Navy sonar.

Anthropogenic activities can lead to changes in animal behavior. Predicting population consequences of these behavioral changes requires integrating short-term individual responses into models that forecast population dynamics across multiple generations. This is especially challenging for long-lived animals, because of the different time scales involved. Beaked whales are a group of deep-diving odontocete whales that respond behaviorally when exposed to military mid-frequency active sonar (MFAS), but the effect of these nonlethal responses on beaked whale populations is unknown. Population consequences of aggregate exposure to MFAS was assessed for two beaked whale populations that are regularly present on U.S. Navy training ranges where MFAS is frequently used. Our approach integrates a wide range of data sources, including telemetry data, information on spatial variation in habitat quality, passive acoustic data on the temporal pattern of sonar use and its relationship to beaked whale foraging activity, into an individual-based model with a dynamic bioenergetic module that governs individual life history. The predicted effect of disturbance from MFAS on population abundance ranged between population extinction to a slight increase in population abundance. These effects were driven by the interaction between the temporal pattern of MFAS use, baseline movement patterns, the spatial distribution of prey, the nature of beaked whale behavioral response to MFAS and the top-down impact of whale foraging on prey abundance. Based on these findings, we provide recommendations for monitoring of marine mammal populations and highlight key uncertainties to help guide future directions for assessing population impacts of nonlethal disturbance for these and other long-lived animals.

Discrete-space continuous-time models of marine mammal exposure to Navy sonar.

Assessing the patterns of wildlife attendance to specific areas is relevant across many fundamental and applied ecological studies, particularly when animals are at risk of being exposed to stressors within or outside the boundaries of those areas. Marine mammals are increasingly being exposed to human activities that may cause behavioral and physiological changes, including military exercises using active sonars. Assessment of the population-level consequences of anthropogenic disturbance requires robust and efficient tools to quantify the levels of aggregate exposure for individuals in a population over biologically relevant time frames. We propose a discrete-space, continuous-time approach to estimate individual transition rates across the boundaries of an area of interest, informed by telemetry data collected with uncertainty. The approach allows inferring the effect of stressors on transition rates, the progressive return to baseline movement patterns, and any difference among individuals. We apply the modeling framework to telemetry data from Blainville's beaked whale (Mesoplodon densirostris) tagged in the Bahamas at the Atlantic Undersea Test and Evaluation Center (AUTEC), an area used by the U.S. Navy for fleet readiness training. We show that transition rates changed as a result of exposure to sonar exercises in the area, reflecting an avoidance response. Our approach supports the assessment of the aggregate exposure of individuals to sonar and the resulting population-level consequences. The approach has potential applications across many applied and fundamental problems where telemetry data are used to characterize animal occurrence within specific areas.

Speciation in the deep: genomics and morphology reveal a new species of beaked whale Mesoplodon eueu.

The deep sea has been described as the last major ecological frontier, as much of its biodiversity is yet to be discovered and described. Beaked whales (ziphiids) are among the most visible inhabitants of the deep sea, due to their large size and worldwide distribution, and their taxonomic diversity and much about their natural history remain poorly understood. We combine genomic and morphometric analyses to reveal a new Southern Hemisphere ziphiid species, Ramari's beaked whale, Mesoplodon eueu, whose name is linked to the Indigenous peoples of the lands from which the species holotype and paratypes were recovered. Mitogenome and ddRAD-derived phylogenies demonstrate reciprocally monophyletic divergence between M. eueu and True's beaked whale (M. mirus) from the North Atlantic, with which it was previously subsumed. Morphometric analyses of skulls also distinguish the two species. A time-calibrated mitogenome phylogeny and analysis of two nuclear genomes indicate divergence began circa 2 million years ago (Ma), with geneflow ceasing 0.35-0.55 Ma. This is an example of how deep sea biodiversity can be unravelled through increasing international collaboration and genome sequencing of archival specimens. Our consultation and involvement with Indigenous peoples offers a model for broadening the cultural scope of the scientific naming process.

Towards a better characterisation of deep-diving whales' distributions by using prey distribution model outputs?

In habitat modelling, environmental variables are assumed to be proxies of lower trophic levels distribution and by extension, of marine top predator distributions. More proximal variables, such as potential prey fields, could refine relationships between top predator distributions and their environment. In situ data on prey distributions are not available over large spatial scales but, a numerical model, the Spatial Ecosystem And POpulation DYnamics Model (SEAPODYM), provides simulations of the biomass and production of zooplankton and six functional groups of micronekton at the global scale. Here, we explored whether generalised additive models fitted to simulated prey distribution data better predicted deep-diver densities (here beaked whales Ziphiidae and sperm whales Physeter macrocephalus) than models fitted to environmental variables. We assessed whether the combination of environmental and prey distribution data would further improve model fit by comparing their explanatory power. For both taxa, results were suggestive of a preference for habitats associated with topographic features and thermal fronts but also for habitats with an extended euphotic zone and with large prey of the lower mesopelagic layer. For beaked whales, no SEAPODYM variable was selected in the best model that combined the two types of variables, possibly because SEAPODYM does not accurately simulate the organisms on which beaked whales feed on. For sperm whales, the increase model performance was only marginal. SEAPODYM outputs were at best weakly correlated with sightings of deep-diving cetaceans, suggesting SEAPODYM may not accurately predict the prey fields of these taxa. This study was a first investigation and mostly highlighted the importance of the physiographic variables to understand mechanisms that influence the distribution of deep-diving cetaceans. A more systematic use of SEAPODYM could allow to better define the limits of its use and a development of the model that would simulate larger prey beyond 1,000 m would probably better characterise the prey of deep-diving cetaceans.

Echolocation click parameters and biosonar behaviour of the dwarf sperm whale (Kogia sima).

Dwarf sperm whales (Kogia sima) are small toothed whales that produce narrow-band high-frequency (NBHF) echolocation clicks. Such NBHF clicks, subject to high levels of acoustic absorption, are usually produced by small, shallow-diving odontocetes, such as porpoises, in keeping with their short-range echolocation and fast click rates. Here, we sought to address the problem of how the little-studied and deep-diving Kogia can hunt with NBHF clicks in the deep sea. Specifically, we tested the hypotheses that Kogia produce NBHF clicks with longer inter-click intervals (ICIs), higher directionality and higher source levels (SLs) compared with other NBHF species. We did this by deploying an autonomous deep-water vertical hydrophone array in the Bahamas, where no other NBHF species are present, and by taking opportunistic recordings of a close-range Kogia sima in a South African harbour. Parameters from on-axis clicks (n=46) in the deep revealed very narrow-band clicks (root mean squared bandwidth, BWRMS, of 3±1 kHz), with SLs of up to 197 dB re. 1 µPa peak-to-peak (µPapp) at 1 m, and a half-power beamwidth of 8.8 deg. Their ICIs (mode of 245 ms) were much longer than those of porpoises (<100 ms), suggesting an inspection range that is longer than detection ranges of single prey, perhaps to facilitate auditory streaming of a complex echo scene. On-axis clicks in the shallow harbour (n=870) had ICIs and SLs in keeping with source parameters of other NBHF cetaceans. Thus, in the deep, dwarf sperm whales use a directional, but short-range echolocation system with moderate SLs, suggesting a reliable mesopelagic prey habitat.

Ambrein: a minor, but common constituent of mammalian faeces?

For nearly 200 years, the only natural source of the alcohol ambrein has been coproliths produced in about 1% of sperm whales and in related jetsam. However, the finding of ambrein in adipocere/faeces of human corpses, led us to hypothesise that ambrein might occur in the faeces of other mammals. Herein, we used a recently developed gas chromatography-mass spectrometry method, with suitable derivatisation of the hindered hydroxy group of ambrein, to screen a number of extracts of mammalian faeces. Minor proportions of ambrein were detected in digested human sewage sludge and in the dung of elephant, domestic cattle, giraffe and buffalo. Whether ambrein formation in the terrestrial species is associated with coprolith formation, is unknown, but solid deposits known as enteroliths and fecaliths occur in humans and some domestic animals.

Oceanographic barriers, divergence, and admixture: Phylogeography and taxonomy of two putative subspecies of short-finned pilot whale.

Genomic phylogeography plays an important role in describing evolutionary processes and their geographic, ecological, or cultural drivers. These drivers are often poorly understood in marine environments, which have fewer obvious barriers to mixing than terrestrial environments. Taxonomic uncertainty of some taxa (e.g., cetaceans), due to the difficulty in obtaining morphological data, can hamper our understanding of these processes. One such taxon, the short-finned pilot whale, is recognized as a single global species but includes at least two distinct morphological forms described from stranding and drive hunting in Japan, the "Naisa" and "Shiho" forms. Using samples (n = 735) collected throughout their global range, we examine phylogeographic patterns of divergence by comparing mitogenomes and nuclear SNP loci. Our results suggest three types within the species: an Atlantic Ocean type, a western/central Pacific and Indian Ocean (Naisa) type, and an eastern Pacific Ocean and northern Japan (Shiho) type. mtDNA control region differentiation indicates these three types form two subspecies, separated by the East Pacific Barrier: Shiho short-finned pilot whale, in the eastern Pacific Ocean and northern Japan, and Naisa short-finned pilot whale, throughout the remainder of the species' distribution. Our data further indicate two diverging populations within the Naisa subspecies, in the Atlantic Ocean and western/central Pacific and Indian Oceans, separated by the Benguela Barrier off South Africa. This study reveals a process of divergence and speciation within a globally-distributed, mobile marine predator, and indicates the importance of the East Pacific Barrier to this evolutionary process.

Physiological, morphological, and ecological tradeoffs influence vertical habitat use of deep-diving toothed-whales in the Bahamas.

Dive capacity among toothed whales (suborder: Odontoceti) has been shown to generally increase with body mass in a relationship closely linked to the allometric scaling of metabolic rates. However, two odontocete species tagged in this study, the Blainville's beaked whale Mesoplodon densirostris and the Cuvier's beaked whale Ziphius cavirostris, confounded expectations of a simple allometric relationship, with exceptionally long (mean: 46.1 min & 65.4 min) and deep dives (mean: 1129 m & 1179 m), and comparatively small body masses (med.: 842.9 kg & 1556.7 kg). These two species also exhibited exceptionally long recovery periods between successive deep dives, or inter-deep-dive intervals (M. densirostris: med. 62 min; Z. cavirostris: med. 68 min). We examined competing hypotheses to explain observed patterns of vertical habitat use based on body mass, oxygen binding protein concentrations, and inter-deep-dive intervals in an assemblage of five sympatric toothed whales species in the Bahamas. Hypotheses were evaluated using dive data from satellite tags attached to the two beaked whales (M. densirostris, n = 12; Z. cavirostris, n = 7), as well as melon-headed whales Peponocephala electra (n = 13), short-finned pilot whales Globicephala macrorhynchus (n = 15), and sperm whales Physeter macrocephalus (n = 27). Body mass and myoglobin concentration together explained only 36% of the variance in maximum dive durations. The inclusion of inter-deep-dive intervals, substantially improved model fits (R2 = 0.92). This finding supported a hypothesis that beaked whales extend foraging dives by exceeding aerobic dive limits, with the extension of inter-deep-dive intervals corresponding to metabolism of accumulated lactic acid. This inference points to intriguing tradeoffs between body size, access to prey in different depth strata, and time allocation within dive cycles. These tradeoffs and resulting differences in habitat use have important implications for spatial distribution patterns, and relative vulnerabilities to anthropogenic impacts.

Beaked whales respond to simulated and actual navy sonar.

Beaked whales have mass stranded during some naval sonar exercises, but the cause is unknown. They are difficult to sight but can reliably be detected by listening for echolocation clicks produced during deep foraging dives. Listening for these clicks, we documented Blainville's beaked whales, Mesoplodon densirostris, in a naval underwater range where sonars are in regular use near Andros Island, Bahamas. An array of bottom-mounted hydrophones can detect beaked whales when they click anywhere within the range. We used two complementary methods to investigate behavioral responses of beaked whales to sonar: an opportunistic approach that monitored whale responses to multi-day naval exercises involving tactical mid-frequency sonars, and an experimental approach using playbacks of simulated sonar and control sounds to whales tagged with a device that records sound, movement, and orientation. Here we show that in both exposure conditions beaked whales stopped echolocating during deep foraging dives and moved away. During actual sonar exercises, beaked whales were primarily detected near the periphery of the range, on average 16 km away from the sonar transmissions. Once the exercise stopped, beaked whales gradually filled in the center of the range over 2-3 days. A satellite tagged whale moved outside the range during an exercise, returning over 2-3 days post-exercise. The experimental approach used tags to measure acoustic exposure and behavioral reactions of beaked whales to one controlled exposure each of simulated military sonar, killer whale calls, and band-limited noise. The beaked whales reacted to these three sound playbacks at sound pressure levels below 142 dB re 1 µPa by stopping echolocation followed by unusually long and slow ascents from their foraging dives. The combined results indicate similar disruption of foraging behavior and avoidance by beaked whales in the two different contexts, at exposures well below those used by regulators to define disturbance.

Field recordings of Gervais' beaked whales Mesoplodon europaeus from the Bahamas.

The first recordings from free-ranging Gervais' beaked whale (Mesoplodon europaeus) are presented. Nine Gervais' beaked whales were observed visually for over 6 h. Clicks were only detected over a 15 min period during the encounter, which coincided with an 88 min period during which no whales were observed at the surface. Click lengths were typically around 200 microS and their dominant energy was in the frequency range 30-50 kHz. While these clicks were broadly similar to those of Cuvier's and Blainville's beaked whales, the Gervais' beaked whale clicks were at a slightly higher frequency than those of the other species.