Search

Goose-beaked whales (Ziphius cavirostris) are deep-diving cetaceans known for their elusive nature and specialized foraging behavior. In 2019 and 2020, six telemetry tags were deployed on these whales in Southern California, resulting in 395 h of acoustic and diving data. Foraging dives were manually identified by the presence of echolocation pulses and buzzes, and generalized additive models assessed factors influencing foraging behavior. The median bathymetric depth at foraging sites was 1,419 m (IQR = 359), and the maximum dive depth was highly correlated with bathymetry depth. Individuals started echolocating on descent at a median depth of 410 m (IQR = 74); pulses were not observed shallower than 292 m. Echolocation ceased during the bottom phase for 81.6% of dives, at a median depth of 1,265 m (IQR = 472); pulses were not observed shallower than 587 m on ascent. The median depth of buzzes was 1,215 m (IQR = 479) with 63% occurring during the bottom phase. Deeper dives correlated with longer durations of diving and echolocation, greater echolocation end depths, and wider bottom phase echolocation depth inter-quartile range. The median difference between dive depth and bottom phase median echolocation depth was 98.3 m (IQR = 48.5), suggesting whales in this region forage in a narrow band close to the seafloor. In the San Nicolas Basin, individuals exhibited longer echolocation durations, produced more pulses, and started and ended echolocating at greater depths compared to adjacent regions. These records validate and expand upon previous studies, providing insights into factors influencing foraging behavior in an area with high anthropogenic disturbance.

Seven passive acoustic surveys for marine mammal sounds were conducted by deploying sonobuoys along ship tracks during Antarctic voyages spanning years 2006-2021. These surveys included nearly 330° of longitude throughout Antarctic (south of 60°S) and sub-Antarctic (between 50-60°S) latitudes. Here, we summarise the presence of calls from critically endangered Antarctic blue whales (Balaenoptera musculus intermedia) detected on all seven of these surveys. We describe and compare the spatial distribution of detections of three different types of Antarctic blue whale calls: unit-A, Z-calls, and D-calls. Three sets of voyages partially overlapped spatially but in different years, providing three regions (Indian Sector, Dumont d’Urville Sea, Ross Sea) to investigate differences over time for these three different call types. The proportion of sonobuoys with calls present was significantly higher in the more recent years for seven of the 15 combinations of years, regions, and call type. The proportion of sonobuoys with calls present was significantly lower only for one of the 15 combinations (unit A in the Ross Sea between 2015 vs 2017), and not significantly different for the remaining seven pairwise comparisons. We discuss possible explanations for these observations including: differences in probability of detection, whale behaviour, whale distribution, and abundance. These explanations are not mutually exclusive and cannot yet be resolved without application of complex analytical methods and collection of additional data. Lastly, we discuss future work that could help clarify the contributions of each of these potential drivers of acoustic detection. We propose continued acoustic data collection, application of new analytical methods, and collection of other synergistic data from Antarctic blue whales on their feeding grounds as a basis for future work on this species. This could provide a cost effective and holistic means of monitoring their status after the effects of 20th century industrial whaling, as well as their responses to natural and anthropogenic changes to their main prey, Antarctic krill, and a changing climate.

Abstract 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.

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.

The distribution of wide-ranging cetacean species often cross national or jurisdictional boundaries, which creates challenges for monitoring populations and managing anthropogenic impacts, especially if data are only available for a portion of the species’ range. Many species found off the U.S. West Coast are known to have continuous distributions into Mexican waters, with highly variable abundance within the U.S. portion of their range. This has contributed to annual variability in design-based abundance estimates from systematic shipboard surveys off the U.S. West Coast, particularly for the abundance of warm temperate species such as striped dolphin, Stenella coeruleoalba, which increases off California during warm-water conditions and decreases during cool-water conditions. Species distribution models (SDMs) can accurately describe shifts in cetacean distribution caused by changing environmental conditions, and are increasingly used for marine species management. However, until recently, data from waters off the Baja California peninsula, México, have not been available for modeling species ranges that span from Baja California to the U.S. West Coast. In this study, we combined data from 1992–2018 shipboard surveys to develop SDMs off the Pacific Coast of Baja California for ten taxonomically diverse cetaceans. We used a Generalized Additive Modeling framework to develop SDMs based on line-transect surveys and dynamic habitat variables from the Hybrid Coordinate Ocean Model (HYCOM). Models were developed for ten species: long- and short-beaked common dolphins (Delphinus delphis delphis and D. d. bairdii), Risso’s dolphin (Grampus griseus), Pacific white-sided dolphin (Lagenorhynchus obliquidens), striped dolphin, common bottlenose dolphin (Tursiops truncatus), sperm whale (Physeter macrocephalus), blue whale (Balaenoptera musculus), fin whale (B. physalus), and humpback whale (Megaptera novaeangliae). The SDMs provide the first fine-scale (approximately 9 x 9 km grid) estimates of average species density and abundance, including spatially-explicit measures of uncertainty, for waters off the Baja California peninsula. Results provide novel insights into cetacean ecology in this region as well as quantitative spatial data for the assessment and mitigation of anthropogenic impacts.

Abstract Species distribution models (SDMs) are important management tools for highly mobile marine species because they provide spatially and temporally explicit information on animal distribution. Two prevalent modeling frameworks used to develop SDMs for marine species are generalized additive models (GAMs) and boosted regression trees (BRTs), but comparative studies have rarely been conducted; most rely on presence‐only data; and few have explored how features such as species distribution characteristics affect model performance. Since the majority of marine species BRTs have been used to predict habitat suitability, we first compared BRTs to GAMs that used presence/absence as the response variable. We then compared results from these habitat suitability models to GAMs that predict species density (animals per km2) because density models built with a subset of the data used here have previously received extensive validation. We compared both the explanatory power (i.e., model goodness of fit) and predictive power (i.e., performance on a novel dataset) of the GAMs and BRTs for a taxonomically diverse suite of cetacean species using a robust set of systematic survey data (1991–2014) within the California Current Ecosystem. Both BRTs and GAMs were successful at describing overall distribution patterns throughout the study area for the majority of species considered, but when predicting on novel data, the density GAMs exhibited substantially greater predictive power than both the presence/absence GAMs and BRTs, likely due to both the different response variables and fitting algorithms. Our results provide an improved understanding of some of the strengths and limitations of models developed using these two methods. These results can be used by modelers developing SDMs and resource managers tasked with the spatial management of marine species to determine the best modeling technique for their question of interest.

Advances in mobile autonomous platforms for oceanographic sensing, including gliders and deep-water profiling floats, have provided new opportunities for passive acoustic monitoring (PAM) of cetaceans. However, there are few direct comparisons of these mobile autonomous systems to more traditional methods, such as stationary bottom-moored recorders. Cross-platform comparisons are necessary to enable interpretation of results across historical and contemporary surveys that use different recorder types, and to identify potential biases introduced by the platform. Understanding tradeoffs across recording platforms informs best practices for future cetacean monitoring efforts. This study directly compares the PAM capabilities of a glider (Seaglider) and a deep-water profiling float (QUEphone) to a stationary seafloor system (High-frequency Acoustic Recording Package, or HARP) deployed simultaneously over a 2 week period in the Catalina Basin, California, United States. Two HARPs were deployed 4 km apart while a glider and deep-water float surveyed within 20 km of the HARPs. Acoustic recordings were analyzed for the presence of multiple cetacean species, including beaked whales, delphinids, and minke whales. Variation in acoustic occurrence at 1-min (beaked whales only), hourly, and daily scales were examined. The number of minutes, hours, and days with beaked whale echolocation clicks were variable across recorders, likely due to differences in the noise floor of each recording system, the spatial distribution of the recorders, and the short detection radius of such a high-frequency, directional signal type. Delphinid whistles and clicks were prevalent across all recorders, and at levels that may have masked beaked whale vocalizations. The number and timing of hours and days with minke whale boing sounds were nearly identical across recorder types, as was expected given the relatively long propagation distance of boings. This comparison provides evidence that gliders and deep-water floats record cetaceans at similar detection rates to traditional stationary recorders at a single point. The spatiotemporal scale over which these single hydrophone systems record sounds is highly dependent on acoustic features of the sound source. Additionally, these mobile platforms provide improved spatial coverage which may be critical for species that produce calls that propagate only over short distances such as beaked whales.

The vaquita (Phocoena sinus) is a small porpoise endemic to Mexico. It is listed by IUCN as Critically Endangered because of unsustainable levels of bycatch in gillnets. The population has been monitored with passive acoustic detectors every summer from 2011 to 2018; here we report results for 2017 and 2018. We combine the acoustic trends with an independent estimate of population size from 2015, and visual observations of at least seven animals in 2017 and six in 2018. Despite adoption of an emergency gillnet ban in May 2015, the estimated rate of decline remains extremely high: 48% decline in 2017 (95% Bayesian credible interval (CRI) 78% decline to 9% increase) and 47% in 2018 (95% CRI 80% decline to 13% increase). Estimated total population decline since 2011 is 98.6%, with greater than 99% probability the decline is greater than 33% yr−1. We estimate fewer than 19 vaquitas remained as of summer 2018 (posterior mean 9, median 8, 95% CRI 6–19). From March 2016 to March 2019, 10 dead vaquitas killed in gillnets were found. The ongoing presence of illegal gillnets despite the emergency ban continues to drive the vaquita towards extinction. Immediate management action is required if the species is to be saved.

We inferred the population densities of blue whales (Balaenoptera musculus) and short-beaked common dolphins (Delphinus delphis) in the Northeast Pacific Ocean as functions of the water-column's physical structure by implementing hierarchical models in a Bayesian framework. This approach allowed us to propagate the uncertainty of the field observations into the inference of species-habitat relationships and to generate spatially explicit population density predictions with reduced effects of sampling heterogeneity. Our hypothesis was that the large-scale spatial distributions of these two cetacean species respond primarily to ecological processes resulting from shoaling and outcropping of the pycnocline in regions of wind-forced upwelling and eddy-like circulation. Physically, these processes affect the thermodynamic balance of the water column, decreasing its volume and thus the height of the absolute dynamic topography (ADT). Biologically, they lead to elevated primary productivity and persistent aggregation of low-trophic-level prey. Unlike other remotely sensed variables, ADT provides information about the structure of the entire water column and it is also routinely measured at high spatial-temporal resolution by satellite altimeters with uniform global coverage. Our models provide spatially explicit population density predictions for both species, even in areas where the pycnocline shoals but does not outcrop (e.g. the Costa Rica Dome and the North Equatorial Countercurrent thermocline ridge). Interannual variations in distribution during El Niño anomalies suggest that the population density of both species decreases dramatically in the Equatorial Cold Tongue and the Costa Rica Dome, and that their distributions retract to particular areas that remain productive, such as the more oceanic waters in the central California Current System, the northern Gulf of California, the North Equatorial Countercurrent thermocline ridge, and the more southern portion of the Humboldt Current System. We posit that such reductions in available foraging habitats during climatic disturbances could incur high energetic costs on these populations, ultimately affecting individual fitness and survival.

For biological populations that form aggregations (or clusters) of individuals, cluster size is an important parameter in line-transect abundance estimation and should be accurately measured. Cluster size in cetaceans has traditionally been represented as the total number of individuals in a group, but group size may be underestimated if group members are spatially diffuse. Groups of false killer whales (Pseudorca crassidens) can comprise numerous subgroups that are dispersed over tens of kilometers, leading to a spatial mismatch between a detected group and the theoretical framework of line-transect analysis. Three stocks of false killer whales are found within the U.S. Exclusive Economic Zone of the Hawaiian Islands (Hawaiian EEZ): an insular main Hawaiian Islands stock, a pelagic stock, and a Northwestern Hawaiian Islands (NWHI) stock. A ship-based line-transect survey of the Hawaiian EEZ was conducted in the summer and fall of 2010, resulting in six systematic-effort visual sightings of pelagic (n = 5) and NWHI (n = 1) false killer whale groups. The maximum number and spatial extent of subgroups per sighting was 18 subgroups and 35 km, respectively. These sightings were combined with data from similar previous surveys and analyzed within the conventional line-transect estimation framework. The detection function, mean cluster size, and encounter rate were estimated separately to appropriately incorporate data collected using different methods. Unlike previous line-transect analyses of cetaceans, subgroups were treated as the analytical cluster instead of groups because subgroups better conform to the specifications of line-transect theory. Bootstrap values (n = 5,000) of the line-transect parameters were randomly combined to estimate the variance of stock-specific abundance estimates. Hawai'i pelagic and NWHI false killer whales were estimated to number 1,552 (CV = 0.66; 95% CI = 479-5,030) and 552 (CV = 1.09; 95% CI = 97-3,123) individuals, respectively. Subgroup structure is an important factor to consider in line-transect analyses of false killer whales and other species with complex grouping patterns.

Species distribution models are now widely used in conservation and management to predict suitable habitat for protected marine species. The primary sources of dynamic habitat data have been in situ and remotely sensed oceanic variables (both are considered “measured data”), but now ocean models can provide historical estimates and forecast predictions of relevant habitat variables such as temperature, salinity, and mixed layer depth. To assess the performance of modeled ocean data in species distribution models, we present a case study for cetaceans that compares models based on output from a data assimilative implementation of the Regional Ocean Modeling System (ROMS) to those based on measured data. Specifically, we used seven years of cetacean line-transect survey data collected between 1991 and 2009 to develop predictive habitat-based models of cetacean density for 11 species in the California Current Ecosystem. Two different generalized additive models were compared: one built with a full suite of ROMS output and another built with a full suite of measured data. Model performance was assessed using the percentage of explained deviance, root mean squared error (RMSE), observed to predicted density ratios, and visual inspection of predicted and observed distributions. Predicted distribution patterns were similar for models using ROMS output and measured data, and showed good concordance between observed sightings and model predictions. Quantitative measures of predictive ability were also similar between model types, and RMSE values were almost identical. The overall demonstrated success of the ROMS-based models opens new opportunities for dynamic species management and biodiversity monitoring because ROMS output is available in near real time and can be forecast.