This paper addresses the problem of automated detection of Z-calls emitted by Antarctic blue whales (B. m. intermedia). The proposed solution is based on a subspace detector of sigmoidal-frequency signals with unknown time-varying amplitude. This detection strategy takes into account frequency variations of blue whale calls as well as the presence of other transient sounds that can interfere with Z-calls (such as airguns or other whale calls). The proposed method has been tested on more than 105 h of acoustic data containing about 2200 Z-calls (as found by an experienced human operator). This method is shown to have a correct-detection rate of up to more than 15% better than the extensible bioacoustic tool package, a spectrogram-based correlation detector commonly used to study blue whales. Because the proposed method relies on subspace detection, it does not suffer from some drawbacks of correlation-based detectors. In particular, it does not require the choice of an a priori fixed and subjective template. The analytic expression of the detection performance is also derived, which provides crucial information for higher level analyses such as animal density estimation from acoustic data. Finally, the detection threshold automatically adapts to the soundscape in order not to violate a user-specified false alarm rate.

1.
A.
Sirovic
,
J. A.
Hildebrand
,
S. M.
Wiggins
, and
D.
Thiele
, “
Blue and fin whale acoustic presence around Antarctica during 2003 and 2004
,”
Mar. Mammal Sci.
25
,
125
136
(
2009
).
2.
A. N.
Gavrilov
and
R. D.
McCauley
, “
Acoustic detection and long-term monitoring of pygmy blue whales over the continental slope in southwest Australia
,”
J. Acoust. Soc. Am.
134
(
3
),
2505
2513
(
2013
).
3.
F.
Samaran
,
K. M.
Stafford
,
T. A.
Branch
,
J.
Gedamke
,
J.-Y.
Royer
,
R. P.
Dziak
, and
C.
Guinet
, “
Seasonal and geographic variation of southern blue whale subspecies in the Indian Ocean
,”
PLoS One
8
(
8
),
1
10
(
2013
).
4.
D.
Ljungblad
,
C. W.
Clark
, and
H.
Shimada
, “
A comparison of sounds attributed to pygmy blue whales Balaenoptera musculus brevicauda recorded south of the Madagascar Plateau and those attributed to true blue whales Balaenoptera musculus recorded off Antarctica
,”
Rep. Int. Whal. Comm.
48
,
431
437
(
1998
).
5.
K. M.
Stafford
,
D. R.
Bohnenstiehl
,
M.
Tolstoy
,
E.
Chapp
,
D. K.
Mellinger
, and
S. E.
Moore
, “
Antarctic type blue whale calls recorded at low latitudes in the Indian and eastern Pacific Oceans
,”
Deep Sea Res., Part I
51
,
1337
1346
(
2004
).
6.
S.
Rankin
,
D.
Ljungblad
,
C.
Clark
, and
H.
Kato
, “
Vocalisations of Antarctic blue whales, Balaenoptera musculus intermedia, recorded during the 2001/2002 and 2002/2003 IWC/SOWER circumpolar cruises, area V, Antarctica
,”
J. Cetacean Res. Manage.
7
(
1
),
13
20
(
2005
).
7.
F.
Samaran
,
O.
Adam
,
J.-F.
Motsch
, and
C.
Guinet
, “
Definition of the Antarctic and pygmy blue whale call templates. Application to fast automatic detection
,”
Can. Acoust.
36
(
1
),
93
103
(
2008
).
8.
A.
Sirovic
,
J. A.
Hildebrand
, and
S. M.
Wiggins
, “
Blue and fin whale call source levels and propagation range in the southern ocean
,”
J. Acoust. Soc. Am.
122
(
2
),
1208
1215
(
2007
).
9.
F.
Samaran
,
C.
Guinet
,
O.
Adam
,
J.-M.
Motsch
, and
Y.
Cansi
, “
Source level estimation of two blue whale subspecies in southwestern Indian Ocean
,”
J. Acoust. Soc. Am.
127
(
6
),
3800
3808
(
2010
).
10.
K. M.
Stafford
, “
Characterization of blue whale calls from the northeast Pacific and development of a matched filter to locate blue whales on U.S. Navy SOSUS (SOund SUrveillance System) arrays
,” Master's thesis,
Oregon State University
,
Corvallis
(
1995
), pp.
1
79
.
11.
D. K.
Mellinger
and
C. W.
Clark
, “
Methods for automatic detection of mysticete sounds
,”
Marine Freshwater Behav. Physiol.
29
(1–4),
163
181
(
1997
).
12.
D. K.
Mellinger
and
C. W.
Clark
, “
Recognizing transient low-frequency whale sounds by spectrogram correlation
,”
J. Acoust. Soc. Am.
107
,
3518
3529
(
2000
).
13.
M. F.
Baumgartner
and
S. E.
Mussoline
, “
A generalized baleen whale call detection and classification system
,”
J. Acoust. Soc. Am.
129
,
2889
2902
(
2011
).
14.
A. M.
Thode
,
K. H.
Kim
,
S. B.
Blackwell
,
C. R.
Greene
,
C. S.
Nation
,
T. L.
McDonald
, and
A. M.
Macrander
, “
Automated detection and localization of bowhead whale sounds in the presence of seismic airgun surveys
,”
J. Acoust. Soc. Am.
131
,
3726
3747
(
2012
).
15.
D.
Harris
,
L.
Matias
,
L.
Thomas
,
J.
Harwood
, and
W. H.
Geissler
, “
Applying distance sampling to fin whale calls recorded by single seismic instruments in the northeast Atlantic
,”
J. Acoust. Soc. Am.
134
(
5
),
3522
3535
(
2013
).
16.
R. P.
Dziak
,
J.-Y.
Royer
,
J. H.
Haxel
,
M.
Delatre
,
D. R.
Bohnenstiehl
,
H.
Matsumoto
,
J.
Goslin
,
C.
Brachet
, and
T.
Lau
, “
Hydroacoustic detection of recent seafloor volcanic activity in the southern Indian Ocean
,” abstract, in
Transactions, American Geophysical Union, Fall Meeting
,
T13 San Francisco, CA
(
2008
).
17.
J.-Y.
Royer
,
R. P.
Dziak
,
M.
Delatre
,
C.
Brachet
,
J. H.
Haxel
,
H.
Matsumoto
,
J.
Goslin
,
V.
Brandon
,
D. R.
Bohnenstiehl
,
C.
Guinet
, and
F.
Samaran
, “
Preliminary results from an hydroacoustic experiment in the Indian Ocean
,” abstract, in
Transactions, American Geophysical Union, Fall Meeting
,
T13, San Francisco, CA
(
2008
).
18.
A. N.
Gavrilov
and
R. D.
McCauley
, “
Steady inter and intra-annual decrease in the vocalization frequency of Antarctic blue whales
,”
J. Acoust. Soc. Am.
131
(
6
),
4476
4480
(
2012
).
19.

Baseband conversion is a classical signal processing method allowing reduction of the sampling frequency of signals whose spectrum is limited to a narrow frequency band. In our case, the spectrum of the signal of interest occupies a 15 Hz bandwidth, centered around 22.5 Hz. The complex baseband signal is obtained by low-pass filtering the received signal downshifted to baseband (i.e., multiplied by e2πj(22.5/250)n, where n denotes the sample index). Its useful spectrum occupies the band [−7.5 7.5] Hz, so that the baseband signal can theoretically be downsampled at a sampling frequency as low as 15 Hz.

20.

To test the noise's stationarity, the methodology described in Ref. 45 that compares the time-frequency structure between the original dataset and an artificially stationarized version of the data was used. With a significance level set to 5%, the test was successful for 88% of the data. The Gaussianity was tested with Lilliefors' test (Ref. 46) with a 5% significance level. 91% of our background noise dataset successfully passed this test.

21.
L. L.
Scharf
,
Statistical Signal Processing: Detection, Estimation, and Time Series Analysis
(
Addison-Wesley
,
Reading, MA
,
1991
), pp.
1
524
.
22.
L. L.
Scharf
and
B.
Friedlander
, “
Matched Subspace Detectors
,”
IEEE Trans. Sign. Process.
42
(
8
),
2146
2157
(
1994
).
23.
S.
Kraut
,
L. L.
Scharf
, and
L. T.
McWhorter
, “
Adaptive subspace detectors
,”
IEEE Trans. Sign. Process.
49
(
1
),
1
16
(
2001
).
24.
Y.
Jin
and
B.
Friedlander
, “
A CFAR adaptive subspace detector for second-order Gaussian signals
,”
IEEE Trans. Sign. Process.
53
(
3
),
871
884
(
2005
).
25.
F.-X.
Socheleau
and
D.
Pastor
, “
Testing the energy of random signals in a known subspace: An optimal invariant approach
,”
IEEE Sign. Process. Lett.
21
(
10
),
1182
1186
(
2014
).
26.
M. A.
McDonald
,
J. A.
Hildebrand
, and
S.
Mesnick
, “
Worldwide decline in tonal frequencies of blue whale songs
,”
Endangered Species Res.
9
(
1
),
13
21
(
2009
).
27.
A. N.
Gavrilov
,
R. D.
McCauley
,
C.
Salgado-Kent
,
J.
Tripovich
, and
C.
Burton
, “
Vocal characteristics of pygmy blue whales and their change over time
,”
J. Acoust. Soc. Am.
130
(
6
),
3651
3660
(
2011
).
28.
B. S.
Miller
,
R.
Leaper
,
S.
Calderan
, and
J.
Gedamke
, “
Red shift, blue shift: Investigating Doppler shifts, blubber thickness, and migration as explanations of seasonal variation in the tonality of Antarctic blue whale song
,”
PloS One
9
(
9
),
e107740
(
2014
).
29.

Note that this relation is rather pessimistic as it assumes statistical independence between overlapping windows. In practice, the actual false alarm will be slightly lower than what is specified.

30.
XBAT
, “
eXtensible BioAcoustic Tool
,” www.birds.cornell.edu/brp/, Cornell Laboratory of Ornithology, Ithaca, NY (Last viewed 6/14/30).
31.
O.
Boisseau
,
D.
Gillespie
,
R.
Leaper
, and
A.
Moscrop
, “
Blue (Balaenoptera musculus) and fin (B. physalus) whale vocalisations measured from northern latitudes of the Atlantic Ocean
,”
J. Cetacean Res. Manage
10
(
1
),
23
30
(
2008
).
32.
F. W.
Shabangu
and
K.
Findlay
, “
Overview of the IWC IDCR/SOWER cruise acoustic survey data
,”
Technical Report No. SC/65b/Forinfo 18
, International Whaling Commission,
2014
.
33.
David K.
Mellinger
, “
A comparison of methods for detecting right whale calls
,”
Can. Acoust.
32
(
2
),
55
65
(
2004
).
34.
P. J.
Huber
and
E. M.
Ronchetti
,
Robust Statistics
, 2nd ed. (
Wiley
,
Hoboken, NJ
,
2009
), pp.
1
354
.
35.
D.
Pastor
and
F.-X.
Socheleau
, “
Robust estimation of noise standard deviation in presence of signals with unknown distributions and occurrences
,”
IEEE Trans. Sign. Process.
60
(
4
),
1545
1555
(
2012
).
36.
F.-X.
Socheleau
,
D.
Pastor
, and
A.
Aissa-El-Bey
, “
Robust statistics based noise variance estimation: Application to wideband interception of non-cooperative communications
,”
IEEE Trans. Aerosp. Electron. Syst.
47
(
1
),
746
755
(
2011
).
37.
A. M.
Zoubir
,
V.
Koivunen
,
Y.
Chakhchoukh
, and
M.
Muma
, “
Robust estimation in signal processing: A tutorial-style treatment of fundamental concepts
,”
IEEE Sign. Process. Mag.
29
(
4
),
61
80
(
2012
).
38.
The approximation becomes an equality when N tends to infinity.
39.
R. M.
Gray
,
Toeplitz and Circulant Matrices: A Review
(
NOW
,
Hanover
,
2006
), pp.
1
90
.
40.
E. L.
Lehmann
and
J. P.
Romano
,
Testing Statistical Hypotheses
, 3rd ed. (
Springer
,
New York
,
2005
), pp.
1
784
.
41.
C.
Walck
, Handbook on Statistical Distribution for Experimentalists, Internal Report No. SUF-PFY/96-01, Stockholm,
2000
, pp.
1
190
.
42.
S.
Broda
and
M. S.
Paolella
, “
Evaluating the density of ratios of noncentral quadratic forms in normal variables
,”
Comput. Stat. Data Anal.
53
(
4
),
1264
1270
(
2009
).
43.
T. Y.
Al-Naffouri
and
B.
Hassibi
, “
On the distribution of indefinite quadratic forms in Gaussian random variables
,” in
Proceedings of the IEEE International Symposium on Information Theory (ISIT'09)
(
2009
), pp.
1744
1748
.
44.
R.
Butler
and
M.
Paolella
, “
Uniform saddlepoint approximations for ratios of quadratic forms
,”
Bernoulli
14
,
140
154
(
2008
).
45.
P.
Borgnat
,
P.
Flandrin
,
P.
Honeine
,
C.
Richard
, and
J.
Xiao
, “
Testing stationarity with surrogates: A time-frequency approach
,”
IEEE Trans. Sign. Process.
58
(
7
),
3459
3470
(
2010
).
46.
H.
Lilliefors
, “
On the Kolmogorov-Smirnov test for normality with mean and variance unknown
,”
J. Am. Stat. Assoc.
62
,
399
402
(
1967
).
You do not currently have access to this content.