Decadal Assessment of Sperm Whale Site-Speciﬁc Abundance Trends in the Northern Gulf of Mexico Using Passive Acoustic Data

: Passive acoustic monitoring has been successfully used to study deep-diving marine mammal populations. To assess regional population trends of sperm whales in the northern Gulf of Mexico (GoM), including impacts of the Deepwater Horizon platform oil spill in 2010, the Littoral Acoustic Demonstration Center-Gulf Ecological Monitoring and Modeling (LADC-GEMM) consortium collected broadband acoustic data in the Mississippi Valley/Canyon area between 2007 and 2017 using bottom-anchored moorings. These data allow the inference of short-term and long-term variations in site-speciﬁc abundances of sperm whales derived from their acoustic activity. A comparison is made between the abundances of sperm whales at speciﬁc sites in different years before and after the oil spill by estimating the regional abundance density. The results show that sperm whales were present in the region throughout the entire monitoring period. A habitat preference shift was observed for sperm whales after the 2010 oil spill with higher activities at sites farther away from the spill site. A comparison of the 2007 and 2015 results shows that the overall regional abundance of sperm whales did not recover to pre-spill levels. The results indicate that long-term spatially distributed acoustic monitoring is critical in characterizing sperm whale population changes and in understanding how environmental stressors impact regional abundances and the habitat use of sperm whales.


Introduction
Sperm whales and other deep-diving marine mammal species produce powerful sounds for echolocation, foraging, and communication [1][2][3][4][5][6]. In the last two decades, several reliable and well-established survey methods have been used to characterize marine mammal behavior and abundance trends such as visual observations, tagging, and passive acoustic monitoring techniques [7][8][9][10][11][12][13][14][15][16][17][18][19][20][21][22]. These techniques are becoming more common in studying the behavior of marine mammals. Although visual observations can identify the presence of marine mammals near the ocean surface and lead to the identification of individuals, they provide no insight into the presence of whales, foraging behavior, or feeding success during deep dives. Visual observations are also limited to daylight times and calm weather conditions. By using acoustic tag data, researchers can obtain information about dive profiles, foraging behavior, and acoustic activity. However, only a small subset of animals gets tagged and the efforts are expensive and time consuming. A few species (e.g., beaked whale species) are very difficult to approach and tag. Passive acoustic monitoring techniques provide a good alternative to simultaneously study many marine mammal species. Passive acoustic systems can record data continuously over long periods and large areas without human involvement. At present, passive acoustic 2 of 16 observations have been successfully used to detect and analyze underwater acoustic signals produced by different marine mammals including Gulf of Mexico (GoM) sperm whales, beaked whales, and dolphins [5][6][7]9,10,[23][24][25].
To assess local abundance trends of marine mammals in the northern Gulf of Mexico, the Littoral Acoustic Demonstration Center-Gulf Ecological Monitoring and Modeling (LADC-GEMM) consortium conducted a series of passive acoustic surveys in the area between 2007 and 2017. One of the research goals was to investigate the impact of the 2010 Deepwater Horizon (DWH) oil spill (20 April 2010-19 September 2010) on the residential population of sperm whales in the northern GoM. Prior to the 2010 oil spill, LADC-GEMM had already collected acoustic data near the oil spill site to study GoM marine mammals using acoustics-based methods [26][27][28][29]. Starting in 2010, the LADC-GEMM consortium redeployed passive acoustic monitoring systems at previous sites to assess the impact of the oil spill on the regional populations of marine mammals [5].
The objective of this paper is to study regional site-specific abundance trends of sperm whales in the past decade and to better understand how abundance variability depends on seasons, locations, and environmental conditions. The remainder of this paper is organized as follows. Section 2 describes the experiment. Section 3 describes the details of the acoustic data processing to detect sperm whale clicks [30] and the density estimation methodology based on the approach developed by Marques et.al [21]. Section 4 presents the results and discussions. The conclusions are drawn and future research directions are outlined in Section 5.

Experiment and Data Collection
The LADC-GEMM consortium collected several passive acoustic datasets in the northern GoM in the past decade as part of long-term monitoring efforts. Acoustic recordings were collected at three different sites, further referred to as the Northern (N), Southern (S), and Western (W) sites (shown in Figure 1). The sites were about 16, 40 and 80 km away from the oil spill site, respectively. At each site, one or more bottom-anchored deepwater moorings containing an Environmental Acoustic Recording System (EARS) buoy were deployed. The acoustic data were collected in 2007 (active recording period between 6 July and 16 July at the Northern and Southern sites), 2010 (active recording period between 10 September and 23 September at all three sites), 2015 (active recording period between 26 June and 22 October at all three sites), 2016 (active recording period between 10 September and 30 January 2017 at the Southern site) and 2017 (active recording period between 9 June and 9 October at all three sites), respectively. Table 1 summarizes the details of each survey year and site deployment. Table 1.
Geographic locations of EARS buoy moorings and the start/deployment and end/retrieval dates. The EARS buoys were continuously recording acoustic data at a 192 kHz sampling rate. The data were stored on internal hard drives using a 16-bit analogue-to-digital converter. Using EARS and hydrophone calibration parameters, the recorded digitized data were later converted to pressures (in microPascals, µPa). For each deployment, the mooring included two paired single-channel EARS buoys; each buoy was connected to a hydrophone with an external cable providing a 1 m hydrophone separation for future localization studies (not discussed in this paper). Figure 2 shows a schematic of the EARS mooring used in the experiment. The moorings were deployed in different water depths at three sites, as shown in Table 1. The hydrophone depths were chosen to target the . The distance between the Northern and Southern sites is 27 km, the distance between the Northern and Western sites is 55 km, and the distance between the Southern and Western sites is 35 km. Yellow triangle indicates the study area.

Year/Station
The EARS buoys were continuously recording acoustic data at a 192 kHz sampling rate. The data were stored on internal hard drives using a 16-bit analogue-to-digital converter. Using EARS and hydrophone calibration parameters, the recorded digitized data were later converted to pressures (in microPascals, µPa). For each deployment, the mooring included two paired single-channel EARS buoys; each buoy was connected to a hydrophone with an external cable providing a 1 m hydrophone separation for future localization studies (not discussed in this paper). Figure 2 shows a schematic of the EARS mooring used in the experiment. The moorings were deployed in different water depths at three sites, as shown in Table 1. The hydrophone depths were chosen to target the foraging depth of marine mammals (such as sperm and beaked whales) and to maximize the probability of recording their directional clicks.

Methods
Sperm whale clicks are typically emitted in series of tens to hundreds of clicks with a regular or regularly changing repetition rate [2]. These clicks include regular clicks for localization and long-range prey echolocation, buzzes for short-range prey echolocation, and codas for communication [2,3,14,31]. A typical sperm whale click train and associated spectrogram are shown in Figure 3. The interclick interval (ICI) varied between 0.5 and 2

Methods
Sperm whale clicks are typically emitted in series of tens to hundreds of clicks with a regular or regularly changing repetition rate [2]. These clicks include regular clicks for localization and long-range prey echolocation, buzzes for short-range prey echolocation, and codas for communication [2,3,14,31]. A typical sperm whale click train and associated spectrogram are shown in Figure 3. The interclick interval (ICI) varied between 0.5 and 2 s [4]. The main acoustic energy of clicks was in a frequency band between 3 and 20 kHz. The time series and spectrogram showed that a higher amplitude direct-path click was followed by lower amplitude clicks from surface-or bottom-reflected paths.

Population Density Estimation
A method of regional site-specific abundance density estimation from sperm whale clicks detected in the EARS continuous datastream was based on the methodology proposed by Marques et al. [17,21]: where nc is the number of detected cues over a time period, T (a cue is defined here as a foraging click), is the estimated proportion of false positive detections, w is the distance away from the hydrophones beyond which cues are assumed not to be detected (the maximum cue detection range), is the estimated average probability of detecting a cue out to distance w, r is the estimated cue production rate by a single mammal and K is the number of hydrophones used for a common detection range (K = 1 for this study). Formula (1) yields a point estimate of the population density of the animals within the detection region of the hydrophone.

Click Detection
In order to detect the acoustic signals of sperm whales and other marine mammal clicks, we developed a multiband spectral energy detection algorithm that was described in detail by Li et al. [30] in the context of detecting and classifying beaked whale clicks. The algorithm tracked received acoustic energy in several representative frequency bands. Clicks from different marine mammal species were discriminated by comparing energy against baseline levels in the bands of interest: low band (3-20 kHz), medium band (25-55 kHz), and high band (60-90 kHz). For the potential detection of sperm whale species, energy was required to be above the background level by a user-specified threshold only in the frequency band of 3-20 kHz. This approach allowed the initial discrimination of sperm whale clicks from beaked whale clicks (when energy was above the threshold in the medium In general, a sperm whale click may contain three or more distinct pulses depending on the angle between the acoustic axis of a whale's head and the propagation direction to a hydrophone due to internal reflections occurring in the head. A single click spectrum has multiple peaks and notches in the frequency band that can be used to determine the animal size and separate individuals in a group [32][33][34][35][36][37][38][39][40][41]. The time spacing between the pulses in a single click is referred to as the interpulse interval (IPI); three or more pulses may be found in a click. The IPI is related to the size of the sperm whale. The value of the IPI varies from 2 ms to 10 ms [40] depending on the size/age of the whales. The duration of a click may reach 20 to 30 ms [2]. Caruso et al. [40] used the IPI to estimate the size distribution of sperm whales detected by a deep hydrophone platform to obtain information on the population structure (size and age/sex) and also to understand changes in a sperm whale population over time. The duration of a click depends on the number of received pulses, which depends on the whale's orientation relative to a receiver.

Population Density Estimation
A method of regional site-specific abundance density estimation from sperm whale clicks detected in the EARS continuous datastream was based on the methodology proposed by Marques et al. [17,21]: where n c is the number of detected cues over a time period, T (a cue is defined here as a foraging click),ĉ is the estimated proportion of false positive detections, w is the distance away from the hydrophones beyond which cues are assumed not to be detected (the maximum cue detection range),P d is the estimated average probability of detecting a cue out to distance w,r is the estimated cue production rate by a single mammal and K is the number of hydrophones used for a common detection range (K = 1 for this study). Formula (1) yields a point estimate of the population density of the animals within the detection region of the hydrophone.

Click Detection
In order to detect the acoustic signals of sperm whales and other marine mammal clicks, we developed a multiband spectral energy detection algorithm that was described in detail by Li et al. [30] in the context of detecting and classifying beaked whale clicks. The algorithm tracked received acoustic energy in several representative frequency bands. Clicks from different marine mammal species were discriminated by comparing energy against baseline levels in the bands of interest: low band (3-20 kHz), medium band (25-55 kHz), and high band (60-90 kHz). For the potential detection of sperm whale species, energy was required to be above the background level by a user-specified threshold only in the frequency band of 3-20 kHz. This approach allowed the initial discrimination of sperm whale clicks from beaked whale clicks (when energy was above the threshold in the medium band) and dolphin clicks (when energy was above the thresholds in all three bands).
The details of sperm whale click detection and counting were as follows. First, for a given time series of acoustic data, a short-time Fourier transform was calculated using a 512 point (2.7 ms) Hamming window with no overlap. Choosing a short FFT window without overlap was based on the practical matter of computation time and to minimize file size: the FFT window without overlap compressed the data by a few orders of magnitude. Furthermore, click detection results matched those of a human expert analysis, therefore it was satisfactory for use.
The value of the sum of the squared spectral amplitudes within each of the three bands of interests for each time window was saved, creating what was referred to as the spectral sum time series. Marine mammal clicks often appear as periodic peaks in the spectral sum time series. The mean and standard deviations of spectral sums were calculated over 10 min time intervals in each band of interest and used to set dynamic detection thresholds. A 10 standard deviation value above the 10 min mean for each band in this study were used as detection thresholds. The thresholds could change every 10 min depending on the acoustic background variability and were unique for each band. Figure 4 shows the spectral sum time series for the 3-20 kHz band with potential sperm whale detections marked by red dots. A potential sperm whale click was defined as any point where the spectral sum exceeded a threshold of 10 standard deviations above the mean in the low frequency band but was below the calculated thresholds in the medium and high bands. The 10 standard deviation threshold allowed the detection of clicks (based on the empirical human inspection and experience) that were assumed to be direct acoustic paths while ignoring the smaller amplitude clicks that were multipath acoustic arrivals (surface-and bottom-reflected paths). When searching for sperm whale clicks, a click event must be present only in the low frequency band. A concurrent detection in the medium frequency or high frequency bands disqualified it as a sperm whale click and the click event was discarded.
To lower the impact of ship noise and other manmade signals that might exist in the low frequency band and cause a high false detection rate, we set restrictions for the frequency and bandwidth of each potential click marked at the spectral energy detection step. For each potential click in the spectral sum time series, the algorithm extracted the peak frequency and bandwidth and fed them into the discriminator algorithm: (i) if the peak frequency was not within 3-20 kHz, a click was disqualified as a sperm whale click; (ii) if the bandwidth (measured at a −10 dB level below the peak of the power spectral density level) was not between 2 and 18 kHz, a click was disqualified as a sperm whale click. The criteria were based on prior knowledge of sperm whale click properties published in peer-reviewed literature [42].
The total number of detected clicks over the specified time interval (T) was then passed to the statistical model Equation (1), which was used for density estimation. thresholds. A 10 standard deviation value above the 10 min mean for each band in this study were used as detection thresholds. The thresholds could change every 10 min depending on the acoustic background variability and were unique for each band. Figure 4 shows the spectral sum time series for the 3-20 kHz band with potential sperm whale detections marked by red dots. A potential sperm whale click was defined as any point where the spectral sum exceeded a threshold of 10 standard deviations above the mean in the low frequency band but was below the calculated thresholds in the medium and high bands. The 10 standard deviation threshold allowed the detection of clicks (based on the empirical human inspection and experience) that were assumed to be direct acoustic paths while ignoring the smaller amplitude clicks that were multipath acoustic arrivals (surfaceand bottom-reflected paths). When searching for sperm whale clicks, a click event must be present only in the low frequency band. A concurrent detection in the medium frequency or high frequency bands disqualified it as a sperm whale click and the click event was discarded. To lower the impact of ship noise and other manmade signals that might exist in the low frequency band and cause a high false detection rate, we set restrictions for the frequency and bandwidth of each potential click marked at the spectral energy detection step. For each potential click in the spectral sum time series, the algorithm extracted the peak frequency and bandwidth and fed them into the discriminator algorithm: (i) if the peak frequency was not within 3-20 kHz, a click was disqualified as a sperm whale click; (ii) if the bandwidth (measured at a −10 dB level below the peak of the power spectral density level) was not between 2 and 18 kHz, a click was disqualified as a sperm whale click. The criteria were based on prior knowledge of sperm whale click properties published in peer-reviewed literature [42].

False Positive Rate
The parameter describing the proportion of false positive detections,ĉ, was estimated by a manual examination of a 10 min interval every two hours across the entire dataset duration for each buoy and each survey period. The average false positive rate for each survey year at each location was used for density estimation.

Probability of Detection
The probability of detection was estimated by using the Monte Carlo simulation methodology proposed by Küsel et al. [19]. The signal-to-noise ratio (SNR) of received sperm whale clicks could be expressed by the source level (SL, measured in dB re 1 µPa at 1 m), directivity loss (DL, in dB), transmission loss (TL, in dB), and the ambient noise level (NL, measured in dB re 1 µPa) as [43]: The SL was the sperm whale's on-axis source level. The average source level of 235 ± 0.5 dB re 1 µPa rms [2] was chosen as the sperm whale's on-axis source level.
The directivity loss (DL) could be obtained by calculating the off-axis attenuation of the source level [19]. The off-axis attenuation of the source level as a function of the off-axis angle for a sperm whale could be estimated by using a piston model [19,44] as follows: where P represented pressure, P 0 was a reference source level, J 1 was the Bessel function of the first kind, f was the frequency, a was the piston radius, θ was the off-axis angle, and c was the speed of sound in the water. The piston model was used to calculate the broadband beam pattern, which could be estimated by integrating P(θ, f ) with respect to the frequency. The broadband beam pattern was given by [44]: where W( f ) was the Gaussian weighting function, which was given by: where f 0 and b were the center frequency and rms bandwidth, respectively. Once the broadband beam pattern was obtained, the directivity index (DI) could be estimated by integrating the modeled broadband beam pattern as follows [44]: Table 2 summarizes the parameters used in this study for the sperm whale directivity loss calculations. The transmission loss (TL) was modeled by taking into account spherical spreading and absorption loss [43]: where r was the range between a hypothetical source and receiver in the Monte Carlo simulations and α f was the absorption coefficient expressed in dB/km. The absorption coefficient could be expressed as [45]: The value for α f was 2.01 dB/km (corresponding to a frequency of 13.4 kHz) as used for the transmission loss calculation.
The ambient noise level (NL) was measured by taking data with minimal sperm whale clicks and anthropogenic noise present by manual inspection. The ambient noise levels were estimated by calculating the mean of the ambient noise power of the data between 3 and 20 kHz in the 24 h period on each hydrophone.
Once all parameters in the sonar equation were obtained, the estimated SNR could be used to model the probability of detection used in the Monte Carlo simulations [19].
The SNRs in this study were randomly simulated for 10,000 clicks using hypothetical animal placement and acoustic axis orientation plus 100 iterations to generate detection probability uncertainty. For each experimental year at each location (Northern, Southern, and Western sites), a 10 min data sample having ample sperm whale clicks was selected. Each click in the data sample was manually marked by an analyst to characterize the detector performance at different SNRs. The probability of detection versus the SNR could be obtained by analyzing the detectability as a function of the measured SNR. The simulated SNRs were then used to estimate the probability of detecting sperm whale clicks. The probability of detection as a function of the SNR could be obtained by assigning the probability of detection corresponding with the detector characterization [19]. Finally, the probability of detection as a function of range could be obtained by converting the probability of detection versus the SNR into the probability of detection versus the range. The average probability of detecting a click within w was then obtained by integrating the product of detection function and the distribution of the horizontal object function [17] over the interval (1, w). The distribution of the horizontal object function was given by h(y) = 2y/w 2 , where y referred to the horizontal distance. The probability of detection (P d ) for each survey year at each location is presented in Section 4.

The Maximum Detection Range
The maximum detection range could be obtained from the probability of detection as a function of range. Figure 5 shows the estimated probability of detection curves of sperm whales as a function of range for the 2015 acoustic dataset at the Northern, Southern, and Western sites, respectively.  Figure 5 illustrates that the modeled maximum detection range was about 25 km, so the maximum detection range, w, in Equation (1) was set to 25 km to estimate the population density and to ensure that on-axis clicks were accounted for longer propagation ranges.
The probability of detection curves at the three sites exhibited differences due to local ambient noise levels and local propagation conditions [46]. These differences were taken into account when the probability of detection was calculated for each site each year.

Production Rate
The click production rate, r, gave the average number of clicks produced per unit time per animal. Click production rates varied with time and location of the acoustic survey [47]. Several published studies have proposed methods for estimating the production rates for different types of marine mammals such as sperm whales and beaked whales [7,17,19,24,47-49].  Figure 5 illustrates that the modeled maximum detection range was about 25 km, so the maximum detection range, w, in Equation (1) was set to 25 km to estimate the population density and to ensure that on-axis clicks were accounted for longer propagation ranges.
The probability of detection curves at the three sites exhibited differences due to local ambient noise levels and local propagation conditions [46]. These differences were taken into account when the probability of detection was calculated for each site each year.

Production Rate
The click production rate,r, gave the average number of clicks produced per unit time per animal. Click production rates varied with time and location of the acoustic survey [47]. Several published studies have proposed methods for estimating the production rates for different types of marine mammals such as sperm whales and beaked whales [7,17,19,24,[47][48][49].
As a first step in the absence of tagged animals, click production rate could be estimated by taking the inverse of the mean ICI. However, according to the tagged data recorded over a typical dive of a sperm whale [31], sperm whale clicks are only produced during the vocally active phase of the dive. Thus, using the inverse of the mean ICI as an estimate of click production rate would result in an overestimated parameter. Therefore, click production rate could be adjusted by taking into account the duration of the active clicking part of the dive cycle:r where t dive represented the duration of the dive cycle, t v represented the vocal duration of the dive, t b represented the time spent producing buzzes during a deep dive, and ICI was the mean interclick interval from the detector. The dive duration of a typical sperm whale is about 45 min and the vocal duration is about 37.4 min with a standard deviation of 7.4 min as obtained from tagged sperm whales [31]. The buzz producing time, t b can be obtained from the product of the average number of buzzes, N b , and the average length of each buzz, L b . In this study, N b was estimated to be 17 per dive with a standard deviation of 8.6 as introduced by Watwood et al. [31], and L b was 5 s as suggested by Fais et al. [50]. Thus, we obtained t b = 1.42 min. The mean ICI in this study was estimated by taking the average of the time intervals between two consecutive detected clicks subject to the condition that the separation between the two consecutive clicks was greater than 0.030 s and less than 2 s in order to be used to calculate intervals between successive clicks.
The production rate was estimated for each site for each survey year due to the variability of the ICI between different locations.

Time Period
The processed time duration used in estimating the density population was set to 45 min, which corresponded with a single dive cycle in the foraging phase [31].

Variance of Density Estimates
To investigate the variance of the density estimations, a non-parametric bootstrap method described in Ackleh et al. [5] was used. The interval estimation ofD to obtain the accuracy of estimates was different from the one followed by Marques et al. [17] where the authors used the first order delta (FOD) method to approximate the variance ofD. The approach implemented in [17] would lead to: The data collection method in [17] was different from this study because several animals were simultaneously tagged in their study. Tagged data allowed direct estimations of variances of individual parameters on the right-hand side of Equation (11). Due to the lack of tagged data, a bootstrap technique was proposed for estimating the variance of D for this study. Equation (1) gave only a point estimate of the site density at a specific location of area.D in (1) was an instantaneous density estimate over a period of time, T. Hence, interval estimates of the site density were obtained by employing the bootstrap approach described in detail in [5]. A total of 5000 iterations at each site and each survey period were implemented to generate the uncertainty and average population density with a 95% confidence interval. Table 3 shows the average number of detected sperm whale clicks over a 45 min window (with an 80% overlap between two consecutive windows) using the detection method described above at the Northern, Southern, and Western sites for each survey year. The average was taken over the entire monitoring period. The estimated proportion of false positives and the probability of click detection are also presented in Table 3. Table 3. Average number of detected clicks in 45 min (n c ), false positive rate (ĉ ) and probability of detection (P d ) at three sites for each survey year.

Year
Northern To investigate the daily and site-specific acoustic activity variability related to a particular site usage, a 24 h average number of sperm whale clicks in every minute, referred to as the average daily click rate on each site, was calculated for each survey year. The average daily click rates of sperm whales at three sites (Northern, Southern, and Western sites) in 2015 over a 110-day observational period from July to October are shown in Figure 6. daily click rates of sperm whales at three sites (Northern, Southern, and Western sites) in 2015 over a 110-day observational period from July to October are shown in Figure 6.  Figure 6 shows that all sites were used during the entire monitoring period. The daily site utilization fluctuated throughout the entire monitoring period and none of the sites showed a strong seasonal preference during summer and fall. Similar dynamics were observed for 2017.  Figure 6 shows that all sites were used during the entire monitoring period. The daily site utilization fluctuated throughout the entire monitoring period and none of the sites showed a strong seasonal preference during summer and fall. Similar dynamics were observed for 2017.
To further investigate regional movements of sperm whales and understand the family clan structure, Figure 7 shows the cross-correlation function for 1-hour detected click numbers between two site pairs. J. Mar. Sci. Eng. 2021, 9, x FOR PEER REVIEW 12 Figure 7. Cross-correlation function for 1 hour detected click numbers between Northern and Southern sites, Northern and Western sites, and Southern and Western sites.
One can conclude from Figure 7 that the cross-correlation coefficient (zero lag va for detections was close to zero for each two-site pair. That was indicative that the sites w not used simultaneously by different groups. Thus, each site could be used to estimat gional abundance independently. The cross-correlation functions between the two showed local maxima that suggested sperm whale group movements across different s Table 4 shows the estimated interclick interval (ICI) and click production rates ( the three sites for each survey year. The ICIs and click production rates () varied slig with time and location. The site-specific ICI and ̂ were used for density estimates.  One can conclude from Figure 7 that the cross-correlation coefficient (zero lag value) for detections was close to zero for each two-site pair. That was indicative that the sites were not used simultaneously by different groups. Thus, each site could be used to estimate regional abundance independently. The cross-correlation functions between the two sites showed local maxima that suggested sperm whale group movements across different sites. Table 4 shows the estimated interclick interval (ICI) and click production rates (r) at the three sites for each survey year. The ICIs and click production rates (r) varied slightly with time and location. The site-specific ICI andr were used for density estimates.  The sperm whale site-specific density estimates were related to the site used by sperm whale groups present in the region. In 2017 (seven years after the spill) the sitespecific density estimates were similar across the three sites, which might be indicative of prey field recovery at the Northern site heavily exposed to leaked oil. In 2017, the average density across the three sites was estimated as 0.616 per 1000 km 2 and could be used as a regional population density measure. The limited tag data available for the GoM sperm whales have shown that a whale can swim between sites in less than one day and 100 km ranging is typical for the GoM groups [51], so a uniform site usage was expected unless there were uneven prey resources. Comparing the density estimates prior to (2007)   The sperm whale site-specific density estimates were related to the site used by sperm whale groups present in the region. In 2017 (seven years after the spill) the site-specific density estimates were similar across the three sites, which might be indicative of prey field recovery at the Northern site heavily exposed to leaked oil. In 2017, the average density across the three sites was estimated as 0.616 per 1000 km 2 and could be used as a regional population density measure. The limited tag data available for the GoM sperm whales have shown that a whale can swim between sites in less than one day and 100 km ranging is typical for the GoM groups [51], so a uniform site usage was expected unless there were uneven prey resources. Comparing the density estimates prior to (2007)  The results indicated that sperm whales probably left the Northern site area, which was in close proximity to the spill, and relocated to the other sites. That may explain the increase at the Southern site. A high concentration of pollutants and prey unavailability were quite likely the cause. By comparing the density estimates in 2010 and 2015, average density of sperm whales at the Northern site increased by 136% from 2010 (0.128/1000 km 2 ) to 2015 (0.302/1000 km 2 ). The results indicated that the density of sperm whales at the site closest to the spill showed recovery from 2010 to 2015, so animals increased the site occupancy, probably due to prey field recovery. A comparison of the densities in 2015 with the densities in 2017 across the same monitoring season showed that the average density in 2017 increased by 111.9%, 144.0%, and 3.47% at the Northern, Southern, and Western sites, respectively. The average site-specific density of sperm whales for each survey year and their associated coefficient of variations (CV) are quantitatively summarized in Table 5. It should be noticed that the smallest variations were at the Western site, which was the least impacted by the spill. The average population density of sperm whales across all sites for each survey year is shown in Figure 8b. To avoid any biased estimates due to limited short-term data in 2007 and 2010, the maximum abundance density was estimated by summing the global maximum density across all three sites (as shown in Figure 8c). Taking into account the data limitations for 2007 (short-term), 2010 (short-term), and 2016 (winter season), the comparison of the 2015 and 2017 results was the most reliable and insightful. The regional use considerably improved between 2015 and 2017, which supported the previously stated hypothesis about the regional ecosystem recovery across all trophic layers.

Conclusions
In this study, the regional site-specific sperm whale density dynamics at three northern GoM sites was investigated using passive acoustic data. The study provided insights into the long-term site usage by sperm whales in the past decade to aid in understanding the impacts of the DWH oil spill on the regional sperm whale population.
The method was based on the analysis of the echolocation click rates of sperm whales. The results indicated that sperm whales were present in the study area throughout the entire monitoring period but site use varied with time and locations. The largest decrease in site usage in 2010 was observed for the Northern site, which was in close proximity to the spill site. The smallest decadal variations in site use were at the Western site, which was not exposed to the spilled oil. Site utilizations leveled in 2017, which might indicate a partial recovery for the lower trophic levels in a food chain where sperm whales were the top predators. This study did not directly address how the spill impacted the long-term health, mortality, fertility and survival of the sperm whale population in the Gulf of Mexico. There are multiple ways by which marine mammals were harmed during the spill, including the consumption of polluted prey, breathing in fumes, and swimming through and ingesting oil. The 2010-2015 studies conducted for coastal bottlenose dolphins in the heavy-oiled Barataria Bay showed an increased mortality (35% greater than expected based on studies of other bottlenose dolphin populations), increased failed pregnancies (46% greater than expected), and an increased likelihood of having adverse health effects (37% greater than expected) [52].
Future research efforts in understanding the long-term impacts of the oil spill on sperm whales and the entire GoM ecosystem are critically dependent on the continued collection of long-term acoustic data at the previously used and new sites. Tagging of GoM marine mammals is of pressing importance to determine the region-specific behavior and parameters used in population density estimates and to benchmark processing tools. The development of new processing algorithms, which would allow the investigation of changes in sperm whale population demographics, seasonal long-range migrations, and feeding successes, and their application to collected data are needed to provide data-based input to the population dynamics models, which will inform the mitigation and recovery efforts for current and future acute and chronic stressors.