Kernel density estimation of conditional distributions to detect responses in satellite tag data

Abstract

Abstract Background As levels of anthropogenic noise in the marine environment rise, it is crucial to quantify potential associated effects on marine mammals. Yet measuring responses is challenging because most species spend the majority of their time submerged. Consequently, much of their sub-surface behavior is difficult or impossible to observe and it can be difficult to determine if—during or following an exposure to sound—an observed dive differs from previously recorded dives. We propose a method for initial assessment of potential behavioral responses observed during controlled exposure experiments (CEEs), in which animals are intentionally exposed to anthropogenic sound sources. To identify possible behavioral responses in dive data collected from satellite-linked time–depth recorders, and to inform the selection and parameters for subsequent individual and population-level response analyses, we propose to use kernel density estimates of conditional distributions for quantitative comparison of pre- and post-exposure behavior. Results We apply the proposed method to nine Cuvier’s beaked whales (Ziphius cavirostris) exposed to a lower-amplitude simulation of Mid-Frequency Active Sonar within the context of a CEE. The exploratory procedure provides evidence that exposure to sound causes animals to change their diving behavior. Nearly all animals tended to dive deep immediately following exposure, potentially indicating avoidance behavior. Following the initial deep dive after exposure, the procedure provides evidence that animals either avoided deep dives entirely or initiated deep dives at unusual times relative to their pre-exposure, baseline behavior patterns. The procedure also provides some evidence that animals exposed as a group may tend to respond as a group. Conclusions The exploratory approach we propose can identify potential behavioral responses across a range of diving parameters observed during CEEs. The method is particularly useful for analyzing data collected from animals for which neither the baseline, unexposed patterns in dive behavior nor the potential types or duration of behavioral responses is well characterized in the literature. The method is able to be applied in settings where little a priori knowledge is known because the statistical analyses employ kernel density estimates of conditional distributions, which are flexible non-parametric techniques. The kernel density estimates allow researchers to initially assess potential behavioral responses without making strong, model-based assumptions about the data.