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.

Ultrasonic antifouling devices negatively impact Cuvier's beaked whales near Guadalupe Island, México.

Widespread use of unregulated acoustic technologies in maritime industries raises concerns about effects on acoustically sensitive marine fauna worldwide. Anthropogenic noise can disrupt behavior and may cause short- to long-term disturbance with possible population-level consequences, particularly for animals with a limited geographic range. Ultrasonic antifouling devices are commercially available, installed globally on a variety of vessel types, and are marketed as an environmentally-friendly method for biofouling control. Here we show that they can be an acoustic disturbance to marine wildlife, as seasonal operation of these hull-mounted systems by tourist vessels in the marine protected area of Guadalupe Island, México resulted in the reduced presence of a potentially resident population of Cuvier's beaked whales (Ziphius cavirostris). Human activities are rapidly altering soundscapes on local and global scales, and these findings highlight the need to identify key noise sources and assess their impacts on marine life to effectively manage oceanic ecosystems.

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.

Recommended snapshot length for acoustic point-transect surveys of intermittently available Cuvier's beaked whales.

Acoustic point-transect distance-sampling surveys have recently been used to estimate the density of beaked whales. Typically, the fraction of short time "snapshots" with detected beaked whales is used in this calculation. Beaked whale echolocation pulses are only intermittently available, which may affect the best choice of snapshot length. The effect of snapshot length on density estimation for Cuvier's beaked whale (Ziphius cavirostris) is investigated by sub-setting continuous recordings from drifting hydrophones deployed off southern and central California. Snapshot lengths from 20 s to 20 min are superimposed on the time series of detected beaked whale echolocation pulses, and the components of the density estimation equation are estimated for each snapshot length. The fraction of snapshots with detections, the effective area surveyed, and the snapshot detection probability all increase with snapshot length. Due to compensatory changes in these three components, density estimates show very little dependence on snapshot length. Within the range we examined, 1-2 min snapshots are recommended to avoid the potential bias caused by animal movement during the snapshot period and to maximize the sample size for estimating the effective area surveyed.

Endangered predators and endangered prey: Seasonal diet of Southern Resident killer whales.

Understanding diet is critical for conservation of endangered predators. Southern Resident killer whales (SRKW) (Orcinus orca) are an endangered population occurring primarily along the outer coast and inland waters of Washington and British Columbia. Insufficient prey has been identified as a factor limiting their recovery, so a clear understanding of their seasonal diet is a high conservation priority. Previous studies have shown that their summer diet in inland waters consists primarily of Chinook salmon (Oncorhynchus tshawytscha), despite that species' rarity compared to some other salmonids. During other times of the year, when occurrence patterns include other portions of their range, their diet remains largely unknown. To address this data gap, we collected feces and prey remains from October to May 2004-2017 in both the Salish Sea and outer coast waters. Using visual and genetic species identification for prey remains and genetic approaches for fecal samples, we characterized the diet of the SRKWs in fall, winter, and spring. Chinook salmon were identified as an important prey item year-round, averaging ~50% of their diet in the fall, increasing to 70-80% in the mid-winter/early spring, and increasing to nearly 100% in the spring. Other salmon species and non-salmonid fishes, also made substantial dietary contributions. The relatively high species diversity in winter suggested a possible lack of Chinook salmon, probably due to seasonally lower densities, based on SRKW's proclivity to selectively consume this species in other seasons. A wide diversity of Chinook salmon stocks were consumed, many of which are also at risk. Although outer coast Chinook samples included 14 stocks, four rivers systems accounted for over 90% of samples, predominantly the Columbia River. Increasing the abundance of Chinook salmon stocks that inhabit the whales' winter range may be an effective conservation strategy for this population.

Diving behaviour of Cuvier's beaked whales exposed to two types of military sonar.

Cuvier's beaked whales (Ziphius cavirostris) have stranded in association with mid-frequency active sonar (MFAS) use, and though the causative mechanism linking these events remains unclear, it is believed to be behaviourally mediated. To determine whether MFAS use was associated with behavioural changes in this species, satellite tags were used to record the diving and movements of 16 Cuvier's beaked whales for up to 88 days in a region of frequent MFAS training off the coast of Southern California. Tag data were combined with summarized records of concurrent bouts of high-power, surface-ship and mid-power, helicopter-deployed MFAS use, along with other potential covariates, in generalized additive mixed-effects models. Deep dives, shallow dives and surface intervals tended to become longer during MFAS use, with some variation associated with the total amount of overlapping MFAS during the behaviour. These changes in dives and surface intervals contributed to a longer interval between deep dives, a proxy for foraging disruption in this species. Most responses intensified with proximity and were more pronounced during mid-power than high-power MFAS use at comparable distances within approximately 50 km, despite the significantly lower source level of mid-power MFAS. However, distance-mediated responses to high-power MFAS, and increased deep dive intervals during mid-power MFAS, were evident up to approximately 100 km away.

Estimation of a Killer Whale (Orcinus orca) Population's Diet Using Sequencing Analysis of DNA from Feces.

Estimating diet composition is important for understanding interactions between predators and prey and thus illuminating ecosystem function. The diet of many species, however, is difficult to observe directly. Genetic analysis of fecal material collected in the field is therefore a useful tool for gaining insight into wild animal diets. In this study, we used high-throughput DNA sequencing to quantitatively estimate the diet composition of an endangered population of wild killer whales (Orcinus orca) in their summer range in the Salish Sea. We combined 175 fecal samples collected between May and September from five years between 2006 and 2011 into 13 sample groups. Two known DNA composition control groups were also created. Each group was sequenced at a ~330bp segment of the 16s gene in the mitochondrial genome using an Illumina MiSeq sequencing system. After several quality controls steps, 4,987,107 individual sequences were aligned to a custom sequence database containing 19 potential fish prey species and the most likely species of each fecal-derived sequence was determined. Based on these alignments, salmonids made up >98.6% of the total sequences and thus of the inferred diet. Of the six salmonid species, Chinook salmon made up 79.5% of the sequences, followed by coho salmon (15%). Over all years, a clear pattern emerged with Chinook salmon dominating the estimated diet early in the summer, and coho salmon contributing an average of >40% of the diet in late summer. Sockeye salmon appeared to be occasionally important, at >18% in some sample groups. Non-salmonids were rarely observed. Our results are consistent with earlier results based on surface prey remains, and confirm the importance of Chinook salmon in this population's summer diet.

Persistent Organic Pollutant Determination in Killer Whale Scat Samples: Optimization of a Gas Chromatography/Mass Spectrometry Method and Application to Field Samples.

Biologic sample collection in wild cetacean populations is challenging. Most information on toxicant levels is obtained from blubber biopsy samples; however, sample collection is invasive and strictly regulated under permit, thus limiting sample numbers. Methods are needed to monitor toxicant levels that increase temporal and repeat sampling of individuals for population health and recovery models. The objective of this study was to optimize measuring trace levels (parts per billion) of persistent organic pollutants (POPs), namely polychlorinated-biphenyls (PCBs), polybrominated-diphenyl-ethers (PBDEs), dichlorodiphenyltrichloroethanes (DDTs), and hexachlorocyclobenzene, in killer whale scat (fecal) samples. Archival scat samples, initially collected, lyophilized, and extracted with 70 % ethanol for hormone analyses, were used to analyze POP concentrations. The residual pellet was extracted and analyzed using gas chromatography coupled with mass spectrometry. Method detection limits ranged from 11 to 125 ng/g dry weight. The described method is suitable for p,p'-DDE, PCBs-138, 153, 180, and 187, and PBDEs-47 and 100; other POPs were below the limit of detection. We applied this method to 126 scat samples collected from Southern Resident killer whales. Scat samples from 22 adult whales also had known POP concentrations in blubber and demonstrated significant correlations (p < 0.01) between matrices across target analytes. Overall, the scat toxicant measures matched previously reported patterns from blubber samples of decreased levels in reproductive-age females and a decreased p,p'-DDE/∑PCB ratio in J-pod. Measuring toxicants in scat samples provides an unprecedented opportunity to noninvasively evaluate contaminant levels in wild cetacean populations; these data have the prospect to provide meaningful information for vital management decisions.

Characterizing a Foraging Hotspot for Short-Finned Pilot Whales and Blainville's Beaked Whales Located off the West Side of Hawai'i Island by Using Tagging and Oceanographic Data.

Satellite tagging data for short-finned pilot whales (Globicephala macrorhynchus) and Blainville's beaked whales (Mesoplodon densirostris) were used to identify core insular foraging regions off the Kona (west) Coast of Hawai'i Island. Ship-based active acoustic surveys and oceanographic model output were used in generalized additive models (GAMs) and mixed models to characterize the oceanography of these regions and to examine relationships between whale density and the environment. The regions of highest density for pilot whales and Blainville's beaked whales were located between the 1000 and 2500 m isobaths and the 250 and 2000 m isobaths, respectively. Both species were associated with slope waters, but given the topography of the area, the horizontal distribution of beaked whales was narrower and located in shallower waters than that of pilot whales. The key oceanographic parameters characterizing the foraging regions were bathymetry, temperature at depth, and a high density of midwater micronekton scattering at 70 kHz in 400-650 m depths that likely represent the island-associated deep mesopelagic boundary community and serve as prey for the prey of the whales. Thus, our results suggest that off the Kona Coast, and potentially around other main Hawaiian Islands, the deep mesopelagic boundary community is key to a food web that supports insular cetacean populations.

Nuclear and mitochondrial patterns of population structure in North Pacific false killer whales (Pseudorca crassidens).

False killer whales (Pseudorca crassidens) are large delphinids typically found in deep water far offshore. However, in the Hawaiian Archipelago, there are 2 resident island-associated populations of false killer whales, one in the waters around the main Hawaiian Islands (MHI) and one in the waters around the Northwestern Hawaiian Islands (NWHI). We use mitochondrial DNA (mtDNA) control region sequences and genotypes from 16 nuclear DNA (nucDNA) microsatellite loci from 206 individuals to examine levels of differentiation among the 2 island-associated populations and offshore animals from the central and eastern North Pacific. Both mtDNA and nucDNA exhibit highly significant differentiation between populations, confirming limited gene flow in both sexes. The mtDNA haplotypes exhibit a strong pattern of phylogeographic concordance, with island-associated populations sharing 3 closely related haplotypes not found elsewhere in the Pacific. However, nucDNA data suggest that NWHI animals are at least as differentiated from MHI animals as they are from offshore animals. The patterns of differentiation revealed by the 2 marker types suggest that the island-associated false killer whale populations likely share a common colonization history, but have limited contemporary gene flow.

First long-term behavioral records from Cuvier's beaked whales (Ziphius cavirostris) reveal record-breaking dives.

Cuvier's beaked whales (Ziphius cavirostris) are known as extreme divers, though behavioral data from this difficult-to-study species have been limited. They are also the species most often stranded in association with Mid-Frequency Active (MFA) sonar use, a relationship that remains poorly understood. We used satellite-linked tags to record the diving behavior and locations of eight Ziphius off the Southern California coast for periods up to three months. The effort resulted in 3732 hr of dive data with associated regional movements--the first dataset of its kind for any beaked whale--and included dives to 2992 m depth and lasting 137.5 min, both new mammalian dive records. Deep dives had a group mean depth of 1401 m (s.d. = 137.8, n = 1142) and duration of 67.4 min (s.d. = 6.9). The group mean time between deep dives was 102.3 min (s.d. = 30.8, n = 783). While the previously described stereotypic pattern of deep and shallow dives was apparent, there was considerable inter- and intra-individual variability in most parameters. There was significant diel behavioral variation, including increased time near the surface and decreased shallow diving at night. However, maximum depth and the proportion of time spent on deep dives (presumed foraging), varied little from day to night. Surprisingly, tagged whales were present within an MFA sonar training range for 38% of days locations were received, and though comprehensive records of sonar use during tag deployments were not available, we discuss the effects frequent acoustic disturbance may have had on the observed behaviors. These data better characterize the true behavioral range of this species, and suggest caution should be exercised when drawing conclusions about behavior using short-term datasets.

Trackline and point detection probabilities for acoustic surveys of Cuvier's and Blainville's beaked whales.

Acoustic survey methods can be used to estimate density and abundance using sounds produced by cetaceans and detected using hydrophones if the probability of detection can be estimated. For passive acoustic surveys, probability of detection at zero horizontal distance from a sensor, commonly called g(0), depends on the temporal patterns of vocalizations. Methods to estimate g(0) are developed based on the assumption that a beaked whale will be detected if it is producing regular echolocation clicks directly under or above a hydrophone. Data from acoustic recording tags placed on two species of beaked whales (Cuvier's beaked whale-Ziphius cavirostris and Blainville's beaked whale-Mesoplodon densirostris) are used to directly estimate the percentage of time they produce echolocation clicks. A model of vocal behavior for these species as a function of their diving behavior is applied to other types of dive data (from time-depth recorders and time-depth-transmitting satellite tags) to indirectly determine g(0) in other locations for low ambient noise conditions. Estimates of g(0) for a single instant in time are 0.28 [standard deviation (s.d.) = 0.05] for Cuvier's beaked whale and 0.19 (s.d. = 0.01) for Blainville's beaked whale.

First direct measurements of behavioural responses by Cuvier's beaked whales to mid-frequency active sonar.

Most marine mammal- strandings coincident with naval sonar exercises have involved Cuvier's beaked whales (Ziphius cavirostris). We recorded animal movement and acoustic data on two tagged Ziphius and obtained the first direct measurements of behavioural responses of this species to mid-frequency active (MFA) sonar signals. Each recording included a 30-min playback (one 1.6-s simulated MFA sonar signal repeated every 25 s); one whale was also incidentally exposed to MFA sonar from distant naval exercises. Whales responded strongly to playbacks at low received levels (RLs; 89-127 dB re 1 µPa): after ceasing normal fluking and echolocation, they swam rapidly, silently away, extending both dive duration and subsequent non-foraging interval. Distant sonar exercises (78-106 dB re 1 µPa) did not elicit such responses, suggesting that context may moderate reactions. The observed responses to playback occurred at RLs well below current regulatory thresholds; equivalent responses to operational sonars could elevate stranding risk and reduce foraging efficiency.

Blue whales respond to simulated mid-frequency military sonar.

Mid-frequency military (1-10 kHz) sonars have been associated with lethal mass strandings of deep-diving toothed whales, but the effects on endangered baleen whale species are virtually unknown. Here, we used controlled exposure experiments with simulated military sonar and other mid-frequency sounds to measure behavioural responses of tagged blue whales (Balaenoptera musculus) in feeding areas within the Southern California Bight. Despite using source levels orders of magnitude below some operational military systems, our results demonstrate that mid-frequency sound can significantly affect blue whale behaviour, especially during deep feeding modes. When a response occurred, behavioural changes varied widely from cessation of deep feeding to increased swimming speed and directed travel away from the sound source. The variability of these behavioural responses was largely influenced by a complex interaction of behavioural state, the type of mid-frequency sound and received sound level. Sonar-induced disruption of feeding and displacement from high-quality prey patches could have significant and previously undocumented impacts on baleen whale foraging ecology, individual fitness and population health.

Acoustic and diving behavior of sperm whales (Physeter macrocephalus) during natural and depredation foraging in the Gulf of Alaska.

Sperm whales have depredated black cod (Anoplopoma fimbria) from demersal longlines in the Gulf of Alaska for decades, but the behavior has recently spread in intensity and geographic coverage. Over a three-year period 11 bioacoustic tags were attached to adult sperm whales off Southeast Alaska during both natural and depredation foraging conditions. Measurements of the animals' dive profiles and their acoustic behavior under both behavioral modes were examined for statistically significant differences. Two rough categories of depredation are identified: "deep" and "shallow." "Deep depredating" whales consistently surface within 500 m of a hauling fishing vessel, have maximum dive depths greater than 200 m, and display significantly different acoustic behavior than naturally foraging whales, with shorter inter-click intervals, occasional bouts of high "creak" rates, and fewer dives without creaks. "Shallow depredating" whales conduct dives that are much shorter, shallower, and more acoustically active than both the natural and deep depredating behaviors, with median creak rates three times that of natural levels. These results suggest that depredation efforts might be measured remotely with passive acoustic monitoring at close ranges.

Near-real-time acoustic monitoring of beaked whales and other cetaceans using a Seaglider™.

In most areas, estimating the presence and distribution of cryptic marine mammal species, such as beaked whales, is extremely difficult using traditional observational techniques such as ship-based visual line transect surveys. Because acoustic methods permit detection of animals underwater, at night, and in poor weather conditions, passive acoustic observation has been used increasingly often over the last decade to study marine mammal distribution, abundance, and movements, as well as for mitigation of potentially harmful anthropogenic effects. However, there is demand for new, cost-effective tools that allow scientists to monitor areas of interest autonomously with high temporal and spatial resolution in near-real time. Here we describe an autonomous underwater vehicle--a glider--equipped with an acoustic sensor and onboard data processing capabilities to passively scan an area for marine mammals in near-real time. The glider was tested extensively off the west coast of the Island of Hawai'i, USA. The instrument covered approximately 390 km during three weeks at sea and collected a total of 194 h of acoustic data. Detections of beaked whales were successfully reported to shore in near-real time. Manual analysis of the recorded data revealed a high number of vocalizations of delphinids and sperm whales. Furthermore, the glider collected vocalizations of unknown origin very similar to those made by known species of beaked whales. The instrument developed here can be used to cost-effectively screen areas of interest for marine mammals for several months at a time. The near-real-time detection and reporting capabilities of the glider can help to protect marine mammals during potentially harmful anthropogenic activities such as seismic exploration for sub-sea fossil fuels or naval sonar exercises. Furthermore, the glider is capable of under-ice operation, allowing investigation of otherwise inaccessible polar environments that are critical habitats for many endangered marine mammal species.

Inferred Paternity and Male Reproductive Success in a Killer Whale (Orcinus orca) Population.

We used data from 78 individuals at 26 microsatellite loci to infer parental and sibling relationships within a community of fish-eating ("resident") eastern North Pacific killer whales (Orcinus orca). Paternity analysis involving 15 mother/calf pairs and 8 potential fathers and whole-pedigree analysis of the entire sample produced consistent results. The variance in male reproductive success was greater than expected by chance and similar to that of other aquatic mammals. Although the number of confirmed paternities was small, reproductive success appeared to increase with male age and size. We found no evidence that males from outside this small population sired any of the sampled individuals. In contrast to previous results in a different population, many offspring were the result of matings within the same "pod" (long-term social group). Despite this pattern of breeding within social groups, we found no evidence of offspring produced by matings between close relatives, and the average internal relatedness of individuals was significantly less than expected if mating were random. The population's estimated effective size was <30 or about 1/3 of the current census size. Patterns of allele frequency variation were consistent with a population bottleneck.

Sperm whale population structure in the eastern and central North Pacific inferred by the use of single-nucleotide polymorphisms, microsatellites and mitochondrial DNA.

We use mitochondrial DNA (mtDNA) (400 bp), six microsatellites and 36 single-nucleotide polymorphisms (SNPs), 20 of which were linked, to investigate population structure of sperm whales (Physeter macrocephalus) in the eastern and central North Pacific. SNP markers, reproducible across technologies and laboratories, are ideal for long-term studies of globally distributed species such as sperm whales, a species of conservation concern because of both historical and contemporary impacts. We estimate genetic differentiation among three strata in the temperate to tropical waters where females are found: California Current, Hawai`i and the eastern tropical Pacific. We then consider how males on sub-Arctic foraging grounds assign to these strata. The California Current stratum was differentiated from both the other strata (P < 0.05) for mtDNA, microsatellites and SNPs, suggesting that the region supports a demographically independent population and providing the first indication that males may exhibit reproductive philopatry. Comparisons between the Hawai`i stratum and the eastern tropical Pacific stratum are not conclusive at this time. Comparisons with Alaska males were statistically significant, or nearly so, from all three strata and individuals showed mixed assignment to, and few exclusions from, the three potential source strata, suggesting widespread origin of males on sub-Arctic feeding grounds. We show that SNPs have sufficient power to detect population structure even when genetic differentiation is low. There is a need for better analytical methods for SNPs, especially when linked SNPs are used, but SNPs appear to be a valuable marker for long-term studies of globally dispersed and highly mobile species.