Rigorous studies of SL for minke whale boing calls have not been published.estimated that the boing call SL was a minimum of 150 dBat 1 m based on received calls from a single whale estimated to be a certain minimum distance from a hydrophone recorder. It is unknown if and how the minke boing call SL varies over time, behavioral state, and environmental conditions. SLs are important, not only for understanding the behavior of an animal, but also for determining the detection range during passive acoustic monitoring for marine mammals and estimating the density or abundance of vocal species. When determining density estimates from acoustic recordings, assumptions are made about the SL of the calls and/or the probability of detecting a call. If SL does change as a function of noise, current methods in acoustic density estimation could be inaccurate. Acoustic density estimation is especially relevant for minke whales since they are a cryptic species and difficult to monitor by traditional visual methods. The purpose of this study was to calculate the SL of minke whale boing calls and compare these SLs over a range of background NLs.

Most marine mammals, including minke whales, have evolved to rely on acoustic calls as their primary form of communication since sound travels efficiently underwater. However, shipping, drilling, military exercises, and other anthropogenic activities have increased the background noise level (NL) in the ocean, which reduces communication range for marine mammals unless they can adapt. Many animals from diverse taxa compensate for noise by increasing the amplitude of their vocalizations (). This phenomenon was first described in humans in 1911 by Etienne Lombard and later termed the Lombard effect (). This increase in amplitude is sometimes also accompanied by increases in frequency, word or call length, and changes in calling rate (). Many experiments have shown evidence for the Lombard effect in birds. For example, the budgerigar or parakeet () increased its call intensity as white background noise within the bandwidth of its calls increased (). Bats change their echolocation behavior in response to noise. Greater horseshoe bats () increased their echolocation amplitude when the noise bandwidth overlapped with their call bandwidth and increased their call frequency in most tested noise conditions (). Some frogs also show evidence for the Lombard effect. Male túngara frogs () increased their call amplitude, calling rate, and call complexity in tank experiments when noise overlapped with their call frequency (). Several species of marine mammals also seem to change their vocalizations in response to noise. Bottlenose dolphins () off Florida increased the apparent output level of their whistles in increased background noise (). Both whistle type and NL seemed to influence the output level of the whistles (). Southern resident killer whales () near the San Juan Islands in Washington increased call duration and the SL of calls as boat noise increased in both short- and long-term studies (). As NLs increased, right whales, especially North Atlantic right whales (), increased the amplitude and frequency of their calls and decreased their calling rate over short time periods and increased the frequency and duration of their calls over long time periods (). Humpback whales () increased the SL of their social calls in increased background noise along their migration route and in a high-latitude feeding area (). Physiological limitations may inhibit animals from fully compensating for background noise, resulting in a reduced communication range. In addition, noise has been linked to increased stress in marine mammals () and decreased reproductive success in songbirds ().

Much is unknown about the minke whale population size, especially the number of individuals that spend the winter in Hawaii, although past studies have used both visual and acoustic methods to estimate their abundance. Since minke whales are small baleen whales, usually seen alone, and produce inconspicuous blows and other surface cues, they are difficult to sight during visual surveys (). No total abundance estimate exists for minke whales in Hawaiian waters since they are rarely seen during summer and fall visual surveys ().estimated minimum density values of 3.64 and 2.77 acoustically active minke whales per hour in a 3780 kmsearch area around the U.S. Navy's Pacific Missile Range Facility (PMRF) before Navy training exercises in February 2011 and 2012, respectively. The proportion of calling animals and the overall calling rate is needed before these minimum density estimates can be converted into animal abundance for the region. In addition, the source level (SL) is necessary to understand the probability of detection in varying noise conditions. These values are difficult to obtain and may change with time, location, season, age, sex, and behavioral state. Visual surveys have been used to estimate seasonal abundance of minke whales in other parts of the eastern North Pacific. Abundance estimates have been calculated off of the west coast of the continental United States (), British Columbia (), the Alaskan Peninsula and Aleutian Islands (), and in the Bering Sea (). Minke whales have also been sighted as far north as the eastern Chukchi Sea in the summer () but sightings were not converted into abundance in those surveys. Most of these abundance estimates are minimum estimates of minke whale abundance since many of the analyses used an assumption that all whales along the trackline were detected (e.g.,). In addition, the probability of visually detecting a minke whale is small so the confidence intervals of these abundance estimates are large.

7. Clarke, J. T., Brower, A. A., Ferguson, M. C., and Willoughby, A. L. (2017a). “ Distribution and relative abundance of marine mammals in the eastern Chukchi and western Beaufort Seas, 2015 and 2016: Final report ,” OCS Study BOEM 2017-019 (National Marine Mammal Laboratory, Alaska Fisheries Science Center, National Marine Fisheries Service, National Oceanic and Seattle, WA).

7. Clarke, J. T., Brower, A. A., Ferguson, M. C., and Willoughby, A. L. (2017a). “ Distribution and relative abundance of marine mammals in the eastern Chukchi and western Beaufort Seas, 2015 and 2016: Final report ,” OCS Study BOEM 2017-019 (National Marine Mammal Laboratory, Alaska Fisheries Science Center, National Marine Fisheries Service, National Oceanic and Seattle, WA).

Although it is known that minke whales have a wide distribution across the North Pacific, much is still unknown about the minke whale population structure in the eastern North Pacific. Like most mysticetes, minke whales are thought to annually migrate between higher latitude feeding areas and lower latitude wintering areas. In a photo-identification study, almost all minke whales photographed off Vancouver Island and Central British Columbia had scars from tropical and subtropical cookiecutter sharks, and minke whales photographed in multiple years gained new scars in between years (). Minke whales have only been sighted off Vancouver and British Columbia from spring until fall (), and minke whale boing calls have been recorded off Hawaii from fall to spring (). Minke whale boing calls detected in the northeastern Chukchi Sea in fall 2009 fit the criteria of central boings, but those detected in fall 2011 were closer to the description of eastern boings (). Although it is currently unknown where the population of minke whales that winters off Hawaii spends the rest of the year, these indications support the hypothesis that they may travel north to more temperate or polar waters and might feed in the Chukchi Sea in the summer.

North Pacific minke whales () produce boing calls. These calls were first described by, and later by, but it took almost 40 years beforewere able to localize the boing calls to a minke whale source. These calls have been recorded between the southwest coast of North America and Hawaii (), off Hawaii (), in the Southern California Bight (), in the Chukchi Sea (), and near the Mariana Trench (). Boing calls from the central North Pacific region have a peak frequency of about 1.4 kHz () and a duration ranging from 1 to 4.5 s (). The call consists of an initial pulse followed by an amplitude- and frequency-modulated component that decreases in amplitude over its duration ().noticed that boing calls had different characteristics depending on where they were recorded in the region between the southwest coast of North America and Hawaii. Central boings, recorded off Hawaii and other areas west of 135°W, had a pulse repetition rate in the amplitude-modulated component of 114–118 pulses per second, while eastern boings, recorded east of 138°W, had a pulse repetition rate of 91–93 pulses per second (). In addition to the pulse repetition rates being slower, the average duration was longer in the eastern boing than the central boing ().observed that the inter-call intervals of boings seemed to change based on whether another minke whale was producing boings in the area. They reported a median inter-call interval of 6 min for individual boing series and 30 s when two separate boing series were recorded at the same time (). Insufficient evidence exists to conclude the purpose of minke whale boing calls, but long-duration recordings can show patterns over time and environmental conditions, which may help to improve understanding of the function of the call and the subset of the population that is vocally active.

SLs were also analyzed in 5 dB NL bins, and the average SL and variance in each bin were calculated. These bins were used to produce histograms of SLs to examine how the shape and character of the distribution changes as a function of noise. The histograms were fit to the data using nonparametric kernel smoothing distributions evaluated at 100 evenly spaced points covering the range of data for each NL bin.

Masked calls create a detection problem that can potentially bias results. Calls made but not detected due to masking cause the mean call SL in a given noise band to be overestimated. If more calls were masked at higher NLs than at lower levels, then the bias in the mean would increase with NL, artificially inducing or overstating any Lombard effect. The impact of masking was minimized in this study by limiting the localization range to 10 km from the center hydrophone, as demonstrated by the probability of detection and localization calculations. However, it was impossible to ensure all calls were not masked, and therefore the sensitivity to masking was investigated by simulating a range of heavier (left) tailed SL distributions.

A generalized additive model (GAM) with the “mgcv” package in) was used to model the relationship between minke whale SL and ocean NL. A Gaussian distribution for the error terms and an identity link function were used with a smoothing term using cubic regression splines with 5 knots (= 5) capturing nonlinearities in the relationship between the predictor and response variable. The exact number of knots is not critical but was chosen conservatively with the intention of producing biologically meaningful results. To ensure that the number of knots were not over-specified, the effective degrees of freedom were used as a guide ().

SLs were estimated by adding the measured RL of each minke whale boing to the expected TL for the animal's position. For each boing, the RL was measured on the master hydrophone and at least three of the four surrounding hydrophones. Therefore, four or five independent measurements of the RL (and thus SL) were available for each call. The average SL and NL across the contributing hydrophones were recorded for each boing, as well as the standard deviation of the values.

In order to maintain a high probability of detection and localization throughout the study area at observed NLs, the maximum allowable radius ( w 2 ) from the center hydrophone was set to 10 km for tracked calls included in the SL analysis. The minimum radius ( w 1 ) for tracked calls was set to 1 km to avoid uncertainties with directionality of calls produced by minke whales and depth of the calling whale, which would have a greater impact on TL at closer ranges.

With g ( r , θ ) characterized over all likely SLs and NLs, the estimated probability of localization, P ̂ L , within a given area surrounding a hydrophone could be calculated. In order for a call to be localized, the call must be detected on the center hydrophone in the subarray and any three of the four surrounding hydrophones in that subarray. Probabilities of localization functions for each of the center hydrophones were created by multiplying g ( r , θ ) at the center hydrophone with the highest three of four g ( r , θ ) probabilities from the adjacent hydrophones, where r and θ from each adjacent hydrophone differ in order to reference the same position as defined by the center hydrophone. The resulting probability of the localization function, g L ( r , θ ) , at the center hydrophone is inherently lower than the detection function, g ( r , θ ) , due to requirement of the call being detected on the center hydrophone and at least three of the adjacent hydrophones. The probability of localization for each of the four subarrays had similar performance since all the hydrophones were at similar depths and bathymetries.

Since the detectability of a call using the GPL detector was dependent on the SL, TL, and NL such that, Monte Carlo simulations were used to characterize the detector performance for the full range of expected TL, SL, and NL. First, the noise was removed from ten minke whale boing signals recorded near the hydrophone with high SNR following the procedure outlined in. Next, 100 h of PMRF noise samples covering the fall, winter, and spring months were collected. The amplitudes of the minke signals were adjusted so that the SLs ranged from 150 to 180 dBat 1 m in 0.5 dB increments. These signals were then reduced in amplitude according to the modeled TL for each location () with TL calculated as described in Sec. II D . The reduced signal was then added to noise, the combined audio was processed with the GPL detector, and the detection was recorded as either detected or missed. Using this technique,was estimated for each hydrophone over all likely combinations of SL and NL. It is important to note that only the amplitudes of the signals were reduced by the expected TL, and any distortions (such as multipath) that may affect the detectability of the call were not simulated. However, for the ranges chosen forand, the signal was minimally distorted by the environment sincewas limited to distances of primarily direct-path propagation.

The potential influence of masking on the detection of minke whale calls was modeled by first estimating the probability of detection at the center hydrophone of each subarray with simulated animal source locations randomly distributed at ranges of 0–20 km from the center hydrophone over all measured SL and NL conditions. The probability of localization was subsequently estimated by combining the detection probabilities over the required number of adjacent hydrophones (at least three) to allow for a valid localization. The estimated probability of detection,, within a given area surrounding a hydrophone was calculated bywhereis the probability density function (PDF) of whale calling locations in the horizontal plane, andis the detection function (). For the purpose of assessing masking, a homogeneous random distribution of animals over the whole area of detection,, was assumed, and so. The accuracy of estimatingrelies on characterizing the range, azimuth, and depth dependent detection function,, in accordance with the detector used. In this paper, the animal was assumed to be near the surface so thatwas taken as a function of range,, and azimuth,, only. The detection function measures the probability of detection from a radial distance from the recording hydrophone () out to a radial distance () over all azimuths. Normally,would be set to zero, with the animal location directly above the hydrophone, but was included here as a variable for reasons explained subsequently. The azimuthal dependence was added to the standard equation to emphasize the complexity caused by bathymetry.

Both the detection and localization processes for any passive acoustic marine mammal monitoring system will be influenced by the acoustic environment, including background noise. Therefore, particular care must be taken when trying to assess the influence of the environment on the marine mammal. Masking, or the addition of natural or artificial sound to the signal of interest that may cover-up or otherwise change the detectability of the sound, is a primary concern when measuring changes in vocal behavior.

Minke whale boing RLs from ranges of 0 to 20 km were measured in this study and used to validate the propagation model. The RL values were binned into 10 m horizontal range increments, and the average values were plotted as a function of range. The TL was calculated for each of the calls using both Peregrine and the geometrical spreading equation, and the average values were also binned in 10 m increments and plotted as a function of range. The TL from horizontal ranges of 0 to 2.5 km was not available from Peregrine due to the inherent limitations of the parabolic equation when modeling high-angle propagation from a source near the surface to a receiver near the bottom in the deep ocean. The TLs calculated from Peregrine and spreading equations were compared over the ranges available.

Both methods of estimating transition loss (the Peregrine model and the geometrical spreading equation) were used to calculate SL by adding TL estimates to the measured RLs of minke whale boing calls. Assuming call SLs are independent of range, correcting the call RLs by the estimated TL should produce similar average SLs over all distances between the animal location and the receiving hydrophones. One caveat to this assumption is that low SL calls at farther ranges could be masked from the detector, limiting detections to those from higher SLs, and therefore causing SLs to trend higher at farther distances (discussed in more detail in Sec. II E ).

The second method for estimating TL was to use the geometrical spreading and attenuation loss equations described by. For slant ranges from the source to the hydrophone greater than the seafloor depth at the source location, the SL was estimated from the RL bywhere SLis the SL (dBat 1 m), RLis the RL (dB),is the transition range in m at which geometrical spreading transitions from spherical to cylindrical,is the attenuation loss coefficient in dB/km, andis the slant range from the whale to the hydrophone in m (). At slant ranges less than, SL was calculated using spherical spreading onlyThe attenuation loss coefficientis primarily influenced by frequency dependent absorption for the relatively short ranges and deep water used in this study. A median frequency of 1400 Hz was used to calculateby means of the method described byand found to be. Since attenuation is minimal compared to geometrical spreading, the “geometrical spreading and attenuation loss equation” is shortened to “geometrical spreading equation” for the remainder of this paper. The transition range,, was hypothesized to be one water depth since the whales are thought to vocalize near the surface, and the hydrophones are raised just above the seafloor. To confirm, a variety of values was tested for, and the binned average SLs were plotted as a function of range. Additionally, TL as a function of range with variousvalues was compared against the Peregrine model.

TL was estimated using two methods. The first method was to use the range dependent acoustic model (RAM) to estimate the TL between each whale call location and each hydrophone location where the call was recorded. The Peregrine software developed by Oasis was utilized for this task. Peregrine is a-language interface to the split-step Pade parabolic equation acoustic propagation code Seahawk (), which is based on RAM (). Peregrine is a general purpose acoustic propagation model suited to the× 2D modeling for this project, where the two-dimensional (2D) model is based on range and depth, but Peregrine adds azimuthal dependence atradials (). The TL over the minke whale boing call bandwidth was calculated by incoherently averaging the TL between the source and receiver in the 1250–1600 Hz band in 5 Hz increments. The TL between the whale location and the hydrophone was approximated by interpolating values derived from TL radials in 60 deg increments for each hydrophone. This azimuthal interpolation was justified because the variability of TL as a function of azimuth was less than 1 dB for ranges and hydrophone locations included in this study. TL was calculated for whale depths between 5 and 100 m, covering the likely depth range at which the animals vocalize. The Peregrine model can import environmental variables that may affect TL. Environmental inputs were interpolated from a variety of four-dimensional [4D; three-dimensional (3D) space plus time] ocean models and bathymetry databases as they were needed in the calculations. Bathymetry information was collected from the National Oceanic and Atmospheric Administration (NOAA) National Geophysical Data Center U.S. Coastal Relief Model () with 3 arc-second resolution. Historical seasonal SSPs were derived from the 2018 World Ocean Atlas (). Within Peregrine, the sediment is treated as an acoustically thick halfspace (implemented as 20 wavelengths at the given frequency, containing an exponential absorptive sponge along the bottom of the sediment layer). Various sediment grain sizes on the Krumbein phi () scale () were chosen as inputs to the model. A sensitivity study was conducted over all likely bottom compositions and SSPs, and TL was calculated for the most likely conditions with bounding values recorded for combinations of SSPs and bottom types that resulted in the highest and lowest TL values.

In order to estimate the SL of minke whale boing calls, both the RL and TL between the source and receiver must be estimated, as described bywhere TLis in units of dB (RL at a rangerelative to RL at 1 m from the source), RLis in units of dB, and SLis in units of dBat 1 m. This equation assumes an omni-directional source and receiver.

Monte Carlo simulations were used to estimate the accuracy of the minke whale boing RL measurements over all likely NLs. Thirty high quality audio recordings of minke whales on PMRF were selected. The signals were then reduced in amplitude and added to 100 h of randomly selected ocean noise recorded at PMRF over all likely signal-to-noise ratios (SNRs), defined aswhere NL and RLare both in units of dB. Next, the GPL detector was used to detect the calls inserted into the noise, and the start time and end time were estimated. The results from the two calculation methods were compared against the known RLs and NLs.

Both the noise measurements and minke whale boing call RLs were converted into decibel units, using 10 log 10 ( RL ) ≡ RL dB , where RL is in units of μ Pa 2 , and RL dB is in units of dB re 1 μ Pa . This method calculates the root mean square (RMS) RL, which is the method used for the remainder of this paper.

To estimate the RL of minke whale boings,was again summed from= 1250 to= 1600 Hz over the duration of the signal as determined by the GPL detector. A minke whale boing spectrogram contains both the signal from the whale call and the background noise. Two methods were explored to determine the most accurate way of removing noise from the estimated RL. In the first method, the call templates as described in Sec. II B ofwere used in place of. The call templates were designed to contain only the spectral contributions of the minke whale boings with the background NL removed. For the second method, the full time series, including both the signal and the noise, was fast Fourier transformed and used forto calculate the spectral level. The NL adjacent to the call (as described above) was subtracted separately to estimate the RL.

To estimate the NL, Sp ( f i ) was summed from f 1 = 1250 to f n = 1600 Hz, which was chosen to cover the dominate frequencies of the minke whale boing call and matches the frequencies used for the call templates. The number of time segments nT was chosen to be 60 samples, which corresponds to approximately 2.5 s of noise. The noise sample was constructed using the time period just before and after the call with a 1 s buffer so that any residual minke boing signal not detected would not be included in the noise sample. The noise samples taken before and after the call did not differ significantly, indicating that there was no signal present in the noise measurements.

To estimate the sound pressure spectral level of the received minke whale boing calls and background noise, the spectral density was calculated and integrated over the frequency bandwidth of interest. The spectral level, or mean square received level (RL), iswhereis the sampling frequency,FFT is the number of samples used in each fast Fourier transform (FFT) window, andis the spectral density and summed overfrequency bins. Spectral density is calculated bywhere the factor of 2 accounts for the energy at negative frequencies, andis the number of time segments incoherently averaged to obtain the spectral density estimate. Within the summation, the quantityis the fast Fourier transformed complex value in the frequency bin corresponding to. In the denominator, the sum ofis the sum of the square of all the points in the window applied to each of thetime series segments before Fourier transforming. Dividing by the ratio of the sampling frequency () and the FFT length (FFT) normalizes by the bin width. AnFFT of 2048, an overlap of 87.5%, a sampling frequencyof 6 kHz, and a Hamming window were used in this analysis.

All tracks were validated by an analyst to be minke whale boings by plotting the locations of the calls, the inter-call intervals, and the corresponding spectrograms for a subset of calls along the track.

A semiautomatic tracker, previously described by, was utilized to spatially and temporally associate localized minke whale calls into individual tracks. Localized calls out to 20 km from the center hydrophone were considered for tracks. These calls were recursively examined so that the elapsed time and distance between calls fit reasonable assumptions of minke whale swimming and calling behavior. A minke whale track required a minimum of 12 localized calls and had a maximum of 3 km and 40 min between successive localizations. The number of counted tracks gave a rough estimate of how many individual whale encounters made up the localized calls used in this study. Although the metrics chosen were somewhat subjective, they did not make a difference in the results of this study and are reported for reproducibility only.

To calculate the location of the minke whales producing the calls, minor modifications were made to the methods outlined in Sec. II of. The minke whale boing call primary frequency components ranged from 1250 to 1600 Hz, so these frequencies were used as the template bounds. As with the localization methods used for Bryde's whales (), single minke whale boing call templates were cross correlated to estimate the TDOA of the call between pairs of hydrophones.showed that for single tonal humpback calls, the timing delay errors were on the order of 40 ms, resulting in localization standard deviations of less than 60 m. Unlike the humpback tonal calls, the initial pulse of the minke whale boing vocalization creates a prominent peak when cross correlated, so it is expected that the timing delay errors are no worse than those for humpback tonal calls. Localization accuracy is important for modeling the TL between the whale location and the recording hydrophone, and discussed in Sec. II D

The 14 hydrophones used to localize the calling minke whales were divided into 4 subarrays (), each containing five hydrophones as shown in Fig.. The TDOAs were computed between the center hydrophone of each subarray and the nearest four corner hydrophones. A localization was valid if the call was detected on the center hydrophone and any three of the four corner hydrophones. The maximum allowable time delay between the center hydrophone and each adjacent hydrophone in the subarray was limited to the direct path propagation time between them. The subarray configuration was chosen so that a direct path solution on the four hydrophone pairs always existed across the monitored area.

The generalized power-law (GPL) detector () was used to detect minke whale boings. The GPL detector marked the start and end time of each minke boing and used a spectral “templating” procedure that subtracted the underlying noise at each frequency band, leaving only the spectral contents of the signal. These templates were later used in the cross-correlation process to obtain TDOAs.

The process for obtaining whale locations can be divided into three steps: detection and feature extraction, cross correlation of those features to obtain time difference of arrivals (TDOAs) of the call at each hydrophone, and TDOA-based localization. These steps are outlined in detail in other publications using calls from humpback whales () and Bryde's whales (), and therefore are only summarized in this paper.

There are several assumptions used throughout this study. Both the animal source and receiving hydrophones were treated as omni-directional, and therefore the minke whale sound directivity was assumed to be zero. The hydrophones on PMRF were designed and tested to be omni-directional so the omni-directionality of the receivers is supported. The SL and NL measurements presented were limited to the 1250–1600 Hz band, which represents the main components of the minke whale boing. The whales were assumed to be near the surface and NLs were recorded on the bottom hydrophones and are therefore only a proxy for the noise experienced by the whale. The sound propagation models used for this study were checked for sensitivity to variability in sound speed profiles (SSPs), sediment type and thickness, and calling depth. Nevertheless, they cannot fully characterize the complexities of transmission loss (TL) and do not account for sea-surface roughness and internal waves, for example. The validity of these assumptions are discussed in Sec. IV

Five years of data were chosen for this study, spanning August 2012–July 2017. All seasons were included, although minke whale boing calls were only recorded in fall, winter, and spring. Figureshows the recording effort (in h) for each month, along with the number of acoustic localizations of minke whale calls within 1–10 km of the center hydrophones of each array for each month.

The U.S. Navy's PMRF is located off the northwest coast of the island of Kauai in the Hawaiian Islands. Since 2002, an array of time-synchronized hydrophones from the PMRF underwater range have been recorded approximately two days a month, in addition to recordings associated with U.S. Navy mid-frequency sonar training events. Over the years, the number of hydrophones in the array and sampling rate has changed, but from August 2012 to July 2017 the array configuration remained the same, containing 47 broadband hydrophones with a 96 kHz sampling rate. In 2014, opportunistic recordings spanning several weeks were added, recording at a 6 kHz sampling rate. Of the 47 hydrophones, 14 offshore hydrophones, at depths of 3150–4700 m and covering a rectangular-shaped grid approximately 20 km to the east/west and 60 km to the north/south, were selected for localization purposes (Fig.). All data recorded at 96 kHz were down-sampled to 6 kHz before processing for consistency.

SLs grouped in 5 dB NL bins are shown in Fig.and Table I . The median RMS SL of boing calls produced in each NL bin was significantly greater than the SL of calls produced in the NL bin centered 5 dB lower (one-sided Wilcoxon's rank sum test,for all four comparisons). In addition to the median SL increasing for calls produced in each increasing NL bin, the variance of the SL distribution significantly decreased when comparing subsequent NL bins for 75–80 dB, 80–85 dB, and 85–90 dB (one-sided Ansari-Bradley test,for these comparisons). The dashed portions of the histograms in Fig.indicate where masking may artificially suppress the number of calls, which can be seen in the tail of the distributions for the higher NLs.

The lower portion of Fig.shows the slope of the SLs as a function of NLs, which reflects the sensitivity response of the whales to increasing background noise. The GAM model suggests minke whales were most responsive to (i.e., greatest slope at) background NLs greater than or equal to 80 dB RMS. The linear fit suggests minke whales increase their RMS SLs by 0.24 dB per 1 dB increase in background noise in the 1250–1600 Hz band. Masking plays a proportionally larger role above 85 dB RMS, and so the uncertainty of the GAM fit increases. Regardless, it appears the whales have a decreased sensitivity at the highest NLs, suggesting they are unable or unwilling to increase their SLs further. There are also proportionally very few calls above 85 dB RMS, which may suggest the whales reduce calling altogether in these high noise conditions.

Although Figs.andsuggest that masking only occurs well into the tail of the SL distributions, and thus is unlikely to be a serious problem, masking was investigated by simulating a range of heavier (left) tailed SL distributions. Observations in the masked region of the SL distributions were simulated using a functional form, where) is the number of calls simulated in the SL interval () dB RMSat 1 m,is a constant controlling the rate of decay in the tail (= 1 for a triangular distribution), andwas chosen so that the tail distribution generated the detected number of calls just above the masked region and reached a value of zero at SLs of 150 dB RMSat 1 m. Values ofwere used, all of which generated substantially heavier tails than observed [between two (= 4) and five (= 1) times as many observations in the masked region of the SL distribution than observed]. GAMs were fitted to each of the reweighted datasets as described previously. The GAM fits for the original and reweighted datasets can be seen in red in both the upper and lower plots of Fig.. The red line represents the GAM fit if no points were missed due to masking, while the dashed red line represents= 2. The most plausible values lie between these lines, as= 1 produces an unlikely elbow to the distribution.

A GAM model was also used to quantify the relationship between SL and NL (Fig., red line). The GAM model explained only 10% of the variability in minke whale SLs (deviance explained) and the range of predicted values was substantially smaller than the range of observed source values. Since a single explanatory variable was not expected to explain a large proportion of the variance in call SLs, this result was not surprising. Residual analysis plots indicated symmetrically distributed residuals that were approximately normal except for a tendency for the model to overpredict SLs at low levels of ocean noise. There was no discernible evidence of heteroskedasticity or unmodeled relationships between residuals and either observed or fitted values of the dependent variable.

TABLE I. The mean, median, and variance of the RMS SL of minke whale boing calls that occurred during a range of NLs. NLs are broken into 5 dB bins and are in units of dB RMS re 1 μ Pa . The bins include the lower NL limit but not the upper limit. Mean and median SLs are in units of dB RMS re 1 μ Pa at 1 m. All calculations were done in the dB domain. The number of calls ( n ) that were recorded during each of these NL bins is also listed.

Minke whale boing calls were estimated to have a median RMS SL of 166 dBat 1 m measured over the 1250–1600 Hz bandwidth and averaged over all NLs. The 25th and 75th percentiles of the RMS SLs were 163 and 168 dBat 1 m, respectively (Table I ). Since each call was recorded at least four separate times (center hydrophone and three or four surrounding hydrophones), the SL was estimated separately and averaged across each hydrophone and the standard deviation calculated. The mean standard deviation of the SL across hydrophones was 0.66 dB. However, the measured increase in SL across NLs was much greater than the standard deviation across hydrophones. Figure(upper) shows the minke whale SL estimates for all 42 159 localized calls, illustrating that SLs increased as background NLs increased. If a linear fit is used (black line), RMS SL increased on average 0.24 dB per 1 dB increased NL (95% confidence interval, 0.23–0.25 dB/dB). The vast majority of SLs are well above the masking zone (shown in black).

Probability of detection and localization were calculated for all background NLs and minke whale boing call SLs to determine what effect masking had on the observed SL results. Figureillustrates the probability of detection and localization for subarrayfor two example NLs, NL = 70 dB and NL = 85 dB, using a distribution of call SLs with a mean of 164 dB RMSat 1 m and variance of 14 dB. This SL distribution was estimated from measured minke boing RLs limited to 10 km in range during very low NLs of 65–70 dB, and so no masking was expected. The high NL maps therefore represent the worst-case scenario for masking since the SL distribution used assumes the minke whales do not change their SL as NL increases. In this scenario, the probability of localization was= 99.8% for NL = 70 dBand= 74.1% for NL = 85 dB, assuming a random spatial distribution of calls. Values forwere determined for all likely combinations of SL and NL for each subarray with the range from the center hydrophone limited to= 1 km and= 10 km. Thevalues for the four subarrays were averaged and, assuming random spatial distribution of animals on the range, provide the average probability of localization for all SLs and NLs over the observed data on PMRF. The resulting “masking zone” can be seen in Fig., where areas of black background indicate= 0%, and areas of white background represent

Although both models are suitable for estimating TL in the study area, the geometrical spreading equation was chosen as the preferred model for the remainder of the study as it predicted no appreciable slope for SL as a function of range, had a faster computation time, and could be computed at closer ranges.

The average minke whale boing RL values as described in Sec. III A were added to the TL estimates from the two models, allowing the SL to be computed as a function of range and shown in the upper portion of Fig.for Peregrine (magenta) and the geometrical spreading equation (black). The average SL derived from the geometrical spreading equation fluctuated by less than 2.9 dB from 1 to 20 km with no appreciable slope (0.029 dB/km using a linear fit). The average SL as a function of range derived from the Peregrine model varied by 5.9 dB or less with a slightly positive slope (0.15 dB/km using a linear fit). In order to eliminate the possibility that masking could bias the SL estimates upward at farther ranges, the same process was repeated using only calls that occurred in noise backgrounds of 70 dBor less (minimal to no masking expected), which resulted in no appreciable change in Fig.. The few calls that were produced within 1 km horizontal range from the center hydrophone were omitted from the analysis because the localization and depth uncertainty of the whale resulted in proportionally more uncertainty in the SL than those at farther ranges. Furthermore, the maximum horizontal range was limited to 10 km in order to reduce uncertainty in the TL and also to minimize the effects of masking (discussed further in Sec. III C ). The red dotted vertical lines in Fig.illustrate the horizontal range over which SL estimates were confined.

The TL was calculated using both the Peregrine model and the geometrical spreading equation for 211 184 individually measured minke whale boings from distances that range between 0 and 20 km from the measuring hydrophone. The average TL can be seen as a function of range in the lower portion of Fig.for Peregrine (purple) and the geometrical spreading equation (black), using 10 m bins. TL estimates from Peregrine could not be reliably estimated for ranges less than about 2.5 km due to limitations inherent with the RAM model formulation at very high propagation angles and are therefore not shown in Fig.. For any given range, the modeled TL will vary slightly for each boing because the hydrophone depths vary between 3150 and 4700 m, and the bathymetry along the path is unique to the animal's position. However, plotting the TL average as a function of range helps to illustrate the differences between the two models. The Peregrine TL model used a seafloor sediment grain size of, assumed animal depth of 30 m, and historical SSPs selected from data that most closely matched the date of the call. Within 10 km, changes in the sediment grain size and SSP had negligible effects on the TL. Beyond 10 km, changing the sediment grain size over all likely values (–8) and all likely SSPs showed 3 dB or less variation. A sediment grain size ofwas chosen from TL experiments using sonobuoys in the region, the details of which are not presented because sediment grain size did not affect TL estimates from Peregrine within 10 km. Changing the assumed animal depth between 5 and 100 m also had negligible effects on TL. The only adjustable parameter for the geometrical spreading model is the transition range,, which was set to the water depth. Choosing values ofgreater than or less than the water depth resulted in less agreement between the two models. Navy surface assets on the range with known SLs and similar frequency ranges were also used to verifyas the predicted SLs from measurements closely matched the known SLs. The largest discrepancies between the geometrical spreading equation and Peregrine occurred at the closest ranges (2.6 dB at a range of 2.5 km). These differences are likely attributable to more complex surface-bottom interactions that are not accounted for with the simpler geometrical spreading model. The close agreement between the models does not necessarily make them both correct, but agreement between the two is reassuring.

The RLs and associated NLs for 42 159 boing calls were measured using the GPL templating technique. The 25th, 50th, and 75th percentiles of RL measurements averaged across the hydrophones were 84, 87, and 90 dB, respectively. NLs associated with these calls had 25th, 50th, and 75th percentiles of 75, 79, and 82 dB, respectively. The measured RLs were added to the TL estimates from Sec. III B to calculate the SLs described in Sec. III D

Two methods (GPL templating and spectrogram noise subtraction) for measuring the RLs were compared over a range of SNRs. Figureshows the measurement accuracy for RL estimates as a function of SNR of the detected call. Both techniques produce RL estimates within 0.5 dB of the true RL at SNRs between -1 and 30 dB. However, the GPL detector more accurately measured the RLs, especially at low SNRs when compared with subtracting the adjacent band-limited noise from the band-limited spectrogram. Therefore, the GPL template measurement results were used for estimating minke whale boing call RLs for the rest of this study. In order to ensure accurate RL estimates, calls with measured SNR less than -5 dB were excluded from the analysis. This cutoff value was subsequently carried forward and accounted for in the probability of detection and localization analysis.

Between 2012 and 2017, 42 159 minke whale boing calls were recorded and tracked at ranges between 1 and 10 km from the center hydrophone of each subarray during opportunistic recordings. SL estimates were derived from these calls by adding the measured RLs and the TL estimates. To verify the TL model, a greater search area was used and 211 184 minke whale boing RLs were measured at ranges between 0 and 20 km. The number of localized and tracked calls per month is shown in Fig.. The hydrophones each recorded a total time of approximately 639 days. An example spectrogram and time series of a minke whale boing call is shown in Fig.. These calls formed 1261 tracks through PMRF.

IV. DISCUSSION Section: Choose Top of page ABSTRACT I. INTRODUCTION II. METHODS III. RESULTS IV. DISCUSSION << V. CONCLUSIONS

These results indicate that minke whales off Hawaii exhibit the Lombard effect in that they increased their boing call intensity in increased background noise. In addition to increasing the average boing call SL, minke whales also decreased the variance of the boing call SL in higher background NLs, suggesting that minke whales were more precise with their call intensity in increased NLs.

7 7 8 Studying the Lombard effect in a natural population of whales is inherently difficult as automated detection and localization systems can produce different results depending on changes in background NLs. The deep bathymetry of the study area has both advantages and disadvantages. Deep water means that propagation between the source and receiver is primarily direct-path over the range considered, which minimizes the impact that uncertainties about sediment type, sediment thickness, SSP, and whale calling depth have on TL (and therefore SL) estimates. Unfortunately, in high noise conditions, the deep hydrophone placement also means that low SL calls, even calls emitted directly above the receiver (one water depth), will be masked from detection. Therefore, determining the true distribution of SLs in high noise conditions is difficult. However, the multitude of sensors available on PMRF allows detection and localization of the vast majority of calls within the ranges and in the noise conditions considered. Because the TL was well characterized at these ranges and the detector performance and estimated RL could be characterized as a function of SNR, the probability of detection and localization and the resulting masking zones in Fig.could be determined with high levels of certainty. Additionally, the symmetrical shapes of the distributions of call SLs as shown in Figs.andindicate that very few calls were likely masked.

The close agreement between the SLs estimated from Peregrine and the geometrical spreading equation, combined with no appreciable slope as a function of range, suggests that the models are able to accurately estimate the TL over the ranges in this study. During this study, the Peregrine model had run time durations of many seconds for each TL range-depth slice as compared to instant computation of simple geometrical spreading. For the locations in this study, Peregrine also did not appear to have any advantage for estimating TL compared with the simpler geometrical spreading equation. Further, the geometrical spreading equation allowed for computation of TL even at very short ranges for this deep ocean location, which is not available from Peregrine.

Erbe et al., 2016 Communication masking in marine mammals: A review and research strategy ,” Mar. Pollut. Bull. 103(1-2), 15– 38 . https://doi.org/10.1016/j.marpolbul.2015.12.007 14. Erbe, C., Reichmuth, C., Cunningham, K., Lucke, K., and Dooling, R. (2016). “,” Mar. Pollut. Bull.(1-2), 15– Clark et al., 2009 Acoustic masking in marine ecosystems: Intuitions, analysis, and implication ,” Mar. Ecol. Prog. Ser. 395, 201– 222 . https://doi.org/10.3354/meps08402 6. Clark, C., Ellison, W., Southall, B., Hatch, L., Parijs, S. V., Frankel, A., and Ponirakis, D. (2009). “,” Mar. Ecol. Prog. Ser., 201– Erbe et al., 2016 Communication masking in marine mammals: A review and research strategy ,” Mar. Pollut. Bull. 103(1-2), 15– 38 . https://doi.org/10.1016/j.marpolbul.2015.12.007 14. Erbe, C., Reichmuth, C., Cunningham, K., Lucke, K., and Dooling, R. (2016). “,” Mar. Pollut. Bull.(1-2), 15– The Lombard effect demonstrated by the minke whales showed a maximum response of 0.34 dB increase in SL per 1 dB increase in background NL when the noise was 82 dB, and an average response of 0.24 dB per 1 dB increase in background NL for the full range of noise encountered. Minke whales, therefore, did not fully compensate for the increase in background NL, and as such the minke whale communication space decreased with increasing noise. No data are available for the hearing sensitivity of minke whales, but for this hypothetical example an assumption is made that minke whales need approximately 0 dB of SNR in the 1250–1600 Hz band to effectively transmit information to a conspecific through their boing call (assuming that a greater SNR is needed for information transmission than for simple detection). Using this 0 dB SNR assumption, then their calls would have a maximum allowable TL of 96 dB in NLs of 67.5 dB, and 78 dB in NLs of 87.5 dB, assuming that whales call at the median SL for the associated noise bin (164 dB and 168 dB, respectively). Using the geometrical spreading equation, a minke whale calling in 67.5 dB of background noise would be able to communicate to a range of approximately 114 km, while one calling in 87.5 dB of background noise would be able to communicate to 19 km. Although the detection range in low NLs would be highly dependent on bottom type, bathymetry, and SSP over these long ranges, it is still likely that an increase of 20 dB in NL would result in minke whales losing 1 order of magnitude of effective communication range (distance) and 2 orders of magnitude of communication space defined as area (distance-squared). Sound production in marine mammals has been associated with a variety of behaviors, including social interaction, group cohesion, feeding, and mating, all of which may be negatively impacted by reduced communication space (). However, quantifying these effects in marine mammal species has been difficult (). Since the exact function of the minke whale boing is poorly understood, it is difficult to determine the impact reduced communication space has on the species. For example, if the boing is mainly used for interaction of conspecifics within a 10 km distance, then the impacts of increasing noise are less than if the boing is used to locate and coordinate with conspecifics at greater distances.

With bottom-mounted hydrophones, it is impossible to know the exact noise environment that the whale experiences along its transit track. The seafloor sensors used for this study are only a proxy for the noise experienced by the whale. As most natural sources of noise in this study originate at the sea surface interface, it is likely that the NL experienced by the whale is greater than what is measured on the seafloor with noise potentially increasing at a faster rate at the surface than at the bottom as wind and waves increase. Therefore, the observable Lombard effect presented in this study is likely an upper-bound response.

Buckland et al. (2001) 4. Buckland, S., Anderson, D., Burnham, K., Laake, J., and Thomas, L. (2001). Introduction to Distance Sampling: Estimating Abundance of Biological Populations ( Oxford University Press , New York ), pp. 1– 448 . The Lombard effect observed in whales has important implications for passive acoustic animal density studies. Inherent to the acoustic marine mammal density estimation equation described byis the probability of detection of a call. This value is often determined from propagation modeling or distance sampling methods. If an animal is able to completely compensate for changes in background noise, then the probability of detection for a study area would be resilient to changes in background noise. However, if the whales respond to the noise but do not completely compensate for the changes as in this study, determining the probability of detection becomes a complex task. If not properly accounted for, population estimates could be skewed by 1 order of magnitude, considering relatively small changes in background NLs can change the detection range by large values.

In this study, it was assumed that the minke whale source was omni-directional. If minke whales produce directional calls, then the RLs recorded on bottom-mounted hydrophones and estimated SLs may be less than the levels if the sounds were recorded on-axis. Receivers at closer ranges and in the same plane as the whale are necessary to measure the true directionality of the call.

Hage et al., 2013 Ambient noise induces independent shifts in call frequency and amplitude within the Lombard effect in echolocating bats ,” Proc. Natl. Acad. Sci. U.S.A. 110(10), 4063– 4068 . https://doi.org/10.1073/pnas.1211533110 20. Hage, S. R., Jiang, T., Berquist, S. W., Feng, J., and Metzner, W. (2013). “,” Proc. Natl. Acad. Sci. U.S.A.(10), 4063– Halfwerk et al., 2016 Vocal responses to noise reveal the presence of the Lombard effect in a frog ,” Behav. Ecol. 27(2), 669– 676 . https://doi.org/10.1093/beheco/arv204 21. Halfwerk, W., Lea, A. M., Guerra, M. A., Page, R. A., and Ryan, M. J. (2016). “,” Behav. Ecol.(2), 669– Manabe et al., 1998 Control of vocal intensity in budgerigars (Melopsittacus undulatus): Differential reinforcement of vocal intensity and the Lombard effect ,” J. Acoust. Soc. Am. 103(2), 1190– 1198 . https://doi.org/10.1121/1.421227 36. Manabe, K., Sadr, E. I., and Dooling, R. J. (1998). “,” J. Acoust. Soc. Am.(2), 1190– The SL and NL measurements were limited to the bandwidth corresponding with the main components of the minke whale boing. In previous laboratory studies with other taxa, individuals responded the most to noise in the same band as their calls (e.g.,). Examining the RLs of minke whale calls in wider bands revealed negligible differences supporting the bandwidth selection.