Parameters used in Matlab simulation.

## Abstract

In the past, cross-correlation-based fisheries stock assessment technique utilized the mean and the ratio of standard deviation to the mean of cross-correlation function (CCF) as estimation parameter. However, in this paper, we have utilized only standard deviation of CCF as estimation parameter to estimate the population size. We utilized four acoustic sensors and considered chirp sound which is commonly generated by damselfish (Dascyllus aruanus), humpback whales (Megaptera novaeangliae), dugongs (Dugong dugon), etc., species to accomplish the simulations. We found that a robust estimation can be obtained using standard deviation of CCF as estimation parameter even when the distances between acoustic sensors are small.

### Keywords

- acoustic sensor
- bins
- chirp
- fisheries stock assessment
- standard deviation

## 1. Introduction

Passive acoustic monitoring of fish abundance is an emerging field of research among the conservation researchers and marine ecologists. It has upgraded understanding of the temporal distribution and repertoire of soniferous fish and mammals [1, 2]. Generally, passive acoustic monitoring is used to have an insight about the population size of soniferous fish and mammals, which are problematic to locate using visual sampling techniques [3, 4, 5, 6] in a certain marine area. These types of fishery surveys utilize the advantage of sound production nature of many species of fish and mammals which possess natural acoustic tags. It has the merit of being a non-destructive and non-invasive monitoring technique, unlike the conventional fisheries stock assessment methods, that is, mark recapture techniques, environmental DNA, visual census, echo, minnow traps, etc. [7, 8]. Generally, mechanical instrument-based conventional fishery surveys suffer from poor accuracy, time consuming nature, overly human interaction, costly instruments, etc., which can be overcome by passive acoustic monitoring techniques. Passive monitoring can provide unbiased data on the location and movement of sound producing source in underwater situations [9]. Low-frequency (<10 kHz) acoustic sensors, that is, hydrophones, are used to detect natural sound production by fish and mammals [10]. Usually, fish sound is associated with courtship, feeding or aggressive encounters [10]. Researchers categorized the sound types of fish and mammals by different names, that is, chirps, pops, grunts, whistles, growls, hoots, etc., which are associated with their frequency and temporal characteristics [11].

However, cross-correlation-based fisheries stock assessment technique, a passive acoustic survey technique, was proposed in [12, 13, 14, 15, 16]. In this technique, the sound signals of vocalizing fish and mammals are processed to estimate their population size [17]. This statistics-based technique has the potential to resolve some main drawbacks of conventional techniques like complexity, reliance on human interaction, time consuming nature of estimation, sensitivity, high cost, etc. A simplified block diagram representation of this technique is illustrated in Figure 1.

In the past, the researchers associated with this technique utilized the mean of CCF [12] ratio of standard deviation to the mean of CCF [13, 14, 15, 16, 17, 18, 19, 20] to estimate population size. In this paper, we have introduced standard deviation of CCF as estimation parameter to perform our desired estimation. We considered four acoustic sensors case [21], that is, hydrophones, in this research. For four acoustic sensors case, different types of topologies, that is, acoustic sensors in line, acoustic sensors in a rectangular shape, acoustic sensors in a triangular shape, are possible. Similarly, Acoustic sensors in a triangular shape can be a square shape, a rhombus shape or a trapezoidal shape. In this paper, we considered acoustic sensors in line case (ASL case). The main reason of considering four acoustic sensors is increasing number of cross-correlation function ensures better accuracy in this technique [14]. Likewise, from diverse sound types of fish and mammals, we considered chirp sound which is commonly generated by damselfish (* Dascyllus aruanus*), humpback whales (

*), dugongs (*Megaptera novaeangliae

*), etc., species [11]. We organize this paper as firstly, to state the theoretical procedure of our proposed methodology and finally, the theory will be evaluated by simulation. We used MATLAB simulation environment to accomplish our simulation in this study.*Dugong dugon

## 2. Utilization of the CCF

The formulation of cross-correlation of sound signals of fish and mammals is analogous to the formulation of cross-correlation of Gaussian signal [22], which are the starting materials to estimate the population size. Chirp sound of fish and mammals are received by the acoustic sensor and recorded in the associated computer in which cross-correlation is executed. Transmission and reception of sound signals are performed for a time frame, called “signal length.” Sound (chirp) generating fish and mammals are considered as the sources of sound signals and * N* fish and mammals are distributed over the volume of a large sphere, the center of which lies halfway between the acoustic sensors. A typical scenario of fish and mammals distribution is shown in Figure 2.

In the water medium, a constant propagation speed * S* of sound is considered [23]. Figure 3 shows an example of 3D estimation area under water space with a single fish

_{p}

N

_{1}and four acoustic sensors

H

_{1},

H

_{2},

H

_{3}, and

H

_{4}. We considered the coordinates of

H

_{1},

H

_{2},

H

_{3}, and

H

_{4}are (

*,*x

_{1}

y

_{1},

z

_{1}), (

x

_{2},

y

_{2},

z

_{2}), (

x

_{3},

y

_{3},

z

_{3}), and (

x

_{4},

y

_{4},

z

_{4}) respectively, whereas the coordinate of the fish is (

*,*a

*,*b

*). The distance between the acoustic sensors can be calculated as follows:*c

Here, _{DBS12} = distance between _{1} and _{2}, _{DBS23} = distance between _{2} and _{3}, and _{DBS34} = distance between _{3} and _{4}.

Let us consider, a sound signal coming from _{1} is _{1}(* t*), which is finite in length. The signal received by acoustic sensors

H

_{1},

H

_{2},

H

_{3}, and

H

_{4}are

S

_{r11},

S

_{r12},

S

_{r13}, and

S

_{r14}, respectively:

where _{11}, _{12}, _{13}, and _{14} are the attenuation due to absorption and dispersion in the medium, and _{11}, _{12}, _{13}, and _{14} are the respective time delays for the acoustic signals to reach the acoustic sensors. For four acoustic sensors ASL case, the cross-correlation among the acoustic sensors is taken place for three times, i.e., between sensors _{1} and _{2}, _{2} and _{3}, and _{3} and _{4}. So, the total number of CCF is three.

Therefore, the CCFs are:

To find out the CCFs for * N* number of fish and mammals, we have to take the total sound signals received by the four acoustic sensors.

Thus, the composite signals received by _{1}, _{2}, _{3}, and _{4} are:

Therefore, the total CCFs are:

This is the form of series of delta functions because in cross-correlation procedure one sound signal is the delayed copy of another [22].

## 3. Theoretical estimation from standard deviation of CCF

As we considered chirp generating fish and mammals to estimate their population size, an introduction to chirp signal is an important task in this perspective. Chirps belong to a swept-frequency sound signal, which possess a time varying frequency. From a sound analysis of * Plectroglyphidodon lacrymatus* and

*species of damselfish, It was seen that the produced chirps by them was consisted of trains of 12–42 short pulses of 3–6 cycles [12, 24]. The durations varied from 0.6 to 1.27 ms where the peak frequency varied from 3400 to 4100 Hz [25]. Such a sound signal can be represented as [8, 12, 13]:*Dascyllus aruanus

where _{1} = starting frequency in Hz, _{2} = ending frequency in Hz, * d* = duration in second,

*= starting phase, and*P

*= amplitude.*A

However, the mean of CCF can be expressed by ensemble average of the chirp-signal cross-correlation as [22].

where _{T} represents the acoustic power of the received signals from the sources taken to be constant over time and space, * v* is the creation rate of the sources whose unit is unit time per unit volume,

T

_{r}is the total recording time,

Now, the variance of the CCF can be defined as [22]:

where

and

where * G*(.) = Green’s function. The other parameters signify their usual meanings [22].

Therefore, we can get the standard deviation, * σ* of the CCF as we know that standard deviation is the square root of the variance.

However, to analyze the random signal cross-correlation problem to find the standard deviation in the above way is very hard. Therefore, the problem can be reframed as a binomial probability problem which can make the analysis simpler. Since, cross-correlation function follows the binomial probability distribution in which the parameters are the number of balls, that is, fish and mammals, * N*, and one on the number of bins,

*; therefore, the standard deviation,*b

*of the CCF is defined as bellow [22]:*σ

where * N* is the number of fish and mammals and

*is the number of bins. Here,*b

*can be achieved from the following Eq. [22]:*b

where _{R} is the sampling rate, _{DBS} is the distance between equidistant sensors, and _{p} is the speed of sound propagation.

From Eq. (25), we can write the following formula:

Therefore, if * σ* is available from simulation, the estimated population size of fish and mammals,

*will be found from Eq. (26).*N

Now, for four acoustic sensors ASL case, the final standard deviation will be found from the average of _{1}_{2}, and _{3.}

Thus, from Eq. (26), we can obtain that

Therefore, if * σ* is available from simulation,

*will be found from Eq. (28).*N

## 4. Simulation and discussion

Simulations were executed considering that four acoustic sensors lay on the center of a sphere. We also considered a uniform random distribution of fish and mammals. Thousand iterations were averaged to accomplish the simulated results. To ease the simulation, the power difference among the acoustic pulses transmitted by each fish and mammal was considered negligible. Here, we considered _{DBS12} = _{DBS23} = _{DBS34} = _{DBS}. The parameters used in MATLAB simulation are introduced in Table 1.

Parameters | Value |
---|---|

Dimension of the sphere | 2000 m |

_{DBS} | 0.25, 0.5, 0.75, 1 m |

_{P} | 1500 m/s |

_{R} | 60 kSa/s |

Absorption coefficient, | 1 dBm^{−1} |

dispersion factor, | 0 |

19, 39, 59, 79 |

Figure 4 shows the theoretical and corresponding simulated results for the population estimation of fish and mammals in terms of the estimation parameter * σ* of CCF. The solid lines designate the theoretical results, and the stars, circles, squares, and triangles correspond the simulated results. The variations of

*are as results of varying*b

d

_{DBS}in the four different Figures 4(a)–(d). The other parameters are same for all the figures.

Figure 5 shows the difference between theoretical and simulated population size of fish and mammals for * b* = 79. In this figure, the solid line indicates the theoretical results, and the triangles are corresponding to the simulated results. From Figure 5, it can be seen that the theoretical and simulated results are closely stayed to each other, which signifies the strength of this population estimation method. Similarly, we can see that the number of bins,

*has an impact on the estimation parameter, which is obvious from Eq. (28). We can see that the value of the standard deviation is lower in case of higher*b

*and vice-versa and the simulated results are closer with the theoretical lines also. The figures also illustrate that a very short distance, even to place a single fish between them, between the acoustic sensors can also give a good estimation using this technique.*b

However, our work has some limitations, for example, assuming the delays to be integer, negligence of multipath interference, consideration of negligible amount of power difference among the fish sound pulses during transmitting time.

## 5. Conclusion

Passive acoustic monitoring is a potential tool to survey the population size of fish and mammals in a certain marine area. It can overcome the major drawbacks of conventional techniques. Cross-correlation-based stock assessment technique is also a passive acoustic survey technique dedicated to fish and mammals. An investigation on this technique with different estimation parameters was the cardinal goal of this research. To do that, we performed our desired estimation with standard deviation of CCF as estimation parameter. The small difference between theoretical and simulated results proved that it is highly possible to pursue this passive monitoring technique utilizing standard deviation of CCF as estimation parameter. Here, we considered four acoustic sensors because from the previous research, we found that an increasing number of CCF ensures better accuracy using this technique. In this paper, we considered four different numbers of bins to show its impact on estimation also. It is shown that a robust estimation is possible using standard deviation of CCF as estimation parameter even when the distances between acoustic sensors are small. Therefore, during practical implementation of this technique, these findings will contribute significantly.

## References

- 1.
Riera A, Pilkington JF, Ford JK, Stredulinsky EH, Chapman NR. Passive acoustic monitoring off Vancouver Island reveals extensive use by at-risk resident killer whale ( Orcinus orca ) populations. Endangered Species Research. 2019;39 :221-234 - 2.
Palmer KJ, Brookes KL, Davies IM, Edwards E, Rendell L. Habitat use of a coastal delphinid population investigated using passive acoustic monitoring. Aquatic Conservation: Marine and Freshwater Ecosystems. 2019; 29 :254-270 - 3.
Gebbie J, Siderius M, Allen JS III. A two-hydrophone range and bearing localization algorithm with performance analysis. The Journal of the Acoustical Society of America. 2015; 137 (3):1586-1597 - 4.
Mouy X, Cabrera De Leo F, Juanes F, Dosso SE. Acoustic estimation of the biodiversity of fish and invertebrates. The Journal of the Acoustical Society of America. 2018; 144 (3):1921-1921 - 5.
Mouy X, Rountree R, Juanes F, Dosso SE. Cataloging fish sounds in the wild using combined acoustic and video recordings. The Journal of the Acoustical Society of America. 2018; 143 (5):EL333–EL339 - 6.
Mouy X, Rountree RA, Juanes F, Dosso SE. Passive acoustic localization of fish using a compact hydrophone array. The Journal of the Acoustical Society of America. 2017; 141 (5):3863-3863 - 7.
Putland RL, Mackiewicz AG, Mensinger AF. Localizing individual soniferous fish using passive acoustic monitoring. Ecological Informatics. 2018; 48 :60-68 - 8.
Hossain SA, Hossen M. A technical review on fish population estimation techniques: Non acoustic and acoustic approaches. Akustika. 2019; 31 :87-103 - 9.
Locascio JV, Mann DA. Localization and source level estimates of black drum ( Pogonias cromis ) calls. The Journal of the Acoustical Society of America. 2011;130 (4):1868-1879 - 10.
Luczkovich JJ, Mann DA, Rountree RA. Passive acoustics as a tool in fisheries science. Transactions of the American Fisheries Society. 2008; 137 (2):533-541 - 11.
Amorim MCP. Diversity of sound production in fish. Communication in Fishes. 2006; 1 :71-104 - 12.
Hossain SA, Hossen M, Anower S. Estimation of damselfish biomass using an acoustic signal processing technique. Journal of Ocean Technology. 2018; 13 (2):92-109 - 13.
Hossain SA, Hossen M. Statistically processing of different sounds of vocalizing fish and mammals to estimate their population size with two acoustic sensors. Marine Technology Society Journal. 2019b; 53 (4):68-80 - 14.
Hossain SA, Hossen M. Population size estimation of chirp and grunt generating fish and mammals using cross-correlation based technique with three acoustic sensors. Journal of Ocean Engineering and Science. 2019c:183-191 - 15.
Hossain SA, Hossen M. Biomass estimation of a popular aquarium fish using an acoustic signal processing technique with three acoustic sensors. In: 2018 International Conference on Advancement in Electrical and Electronic Engineering (ICAEEE). IEEE; 2018. pp. 1-4 - 16.
Rana MS, Anower MS, Siraj SN, Haque MI. A signal processing approach of fish abundance estimation in the sea. In: 9th International Forum on Strategic Technology (IFOST). IEEE; 2014. pp. 87-90 - 17.
Hossain SA, Hossen M. Selection of the optimum estimation parameter in cross-correlation based fisheries stock assessment technique. Journal of Ocean Technology. 2020; 15 (2):105-119 - 18.
Hossain SA, Hossen M. Impact of underwater bandwidth and SNR on cross-correlation based population estimation technique of fish and mammals. Underwater Technology. 2019d; 36 (2):19-27 - 19.
Hossain SA, Hossen M. Error calculation in cross-correlation based population estimation technique of fish and mammals. Acoustical Science and Technology. 2019e; 40 (6):402-405 - 20.
Hossain SA, Hossen M. Impact of dispersion coefficient on cross-correlation based population estimation technique of fish and mammals. Journal of Ocean Technology. 2019f; 14 (3):79-92 - 21.
Hossain SM. Cross-correlation based acoustic signal processing technique and its implementation on marine ecology [Doctoral dissertation]. Khulna, Bangladesh: Khulna University of Engineering & Technology (KUET); 2018 - 22.
Anower MS. Estimation using cross-correlation in a communications network [PhD dissertation]. Australian Defence Force Academy; 2011 - 23.
Hossain SA, Mallik A, Arefin MA. A signal processing approach to estimate underwater network cardinalities with lower complexity. Journal of Electrical and Computer Engineering Innovations (JECEI). 2017; 5 (2):131-138 - 24.
Hossain SA, Mallik A, Hossen M. An analytical analysis on fish sounds. Akustika. 2019; 33 :15-23 - 25.
Parmentier E, Vandewalle P, Frederich B, Fine ML. Sound production in two species of damselfishes ( Pomacentridae ):Plectroglyphidodon lacrymatus andDascyllus aruanus . Journal of Fish Biology. 2006;69 (2):491-503