Open access peer-reviewed chapter

Cross-Correlation-Based Fisheries Stock Assessment Technique: Utilization of Standard Deviation of Cross-Correlation Function as Estimation Parameter with Four Acoustic Sensors

By Shaik Asif Hossain and Monir Hossen

Submitted: February 12th 2020Reviewed: June 23rd 2020Published: July 13th 2020

DOI: 10.5772/intechopen.93240

Downloaded: 88


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.


  • 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.

Figure 1.

Simplified block diagram of cross-correlation fisheries-based stock assessment system.

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 (Megaptera novaeangliae), dugongs (Dugong dugon), 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.


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 Nfish 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.

Figure 2.

Distribution of fish and mammals with four acoustic sensors, that is, four pluses (++++).

In the water medium, a constant propagation speed Spof sound is considered [23]. Figure 3 shows an example of 3D estimation area under water space with a single fish N1 and four acoustic sensors H1, H2, H3, and H4. We considered the coordinates of H1, H2, H3, and H4 are (x1, y1, z1), (x2, y2, z2), (x3, y3, z3), and (x4, y4, z4) respectively, whereas the coordinate of the fish is (a, b, c). The distance between the acoustic sensors can be calculated as follows:

Figure 3.

Diagram of a single fish with four acoustic sensors,H1,H2,H3, andH4.


Here, dDBS12 = distance between H1 and H2, dDBS23 = distance between H2 and H3, and dDBS34 = distance between H3 and H4.

Let us consider, a sound signal coming from N1 is S1(t), which is finite in length. The signal received by acoustic sensors H1, H2, H3, and H4 are Sr11, Sr12, Sr13, and Sr14, 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 H1 and H2, H2 and H3, and H3 and H4. So, the total number of CCF is three.

Therefore, the CCFs are:


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

Thus, the composite signals received by H1, H2, H3, and H4 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 lacrymatusand Dascyllus aruanusspecies 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]:


where f1 = starting frequency in Hz, f2 = ending frequency in Hz, d = duration in second, P = starting phase, and A = amplitude.

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


where QT represents the acoustic power of the received signals from the sources taken to be constant over time and space, vis the creation rate of the sources whose unit is unit time per unit volume, Tr is the total recording time, rsis the path length of sources from the origin, rais the path length of first acoustic sensor from the origin, and rbis the path length of second sensor from the origin.

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


where Ct2and C2tare defined in Eqs. (21) and (22), respectively, as [22]:




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, b; therefore, the standard deviation, σof the CCF is defined as bellow [22]:


where Nis the number of fish and mammals and bis the number of bins. Here, bcan be achieved from the following Eq. [22]:


where SR is the sampling rate, dDBS is the distance between equidistant sensors, and Sp 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, Nwill be found from Eq. (26).

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, Nwill be found from Eq. (28).

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 dDBS12 = dDBS23 = dDBS34 = dDBS. The parameters used in MATLAB simulation are introduced in Table 1.

Dimension of the sphere2000 m
dDBS0.25, 0.5, 0.75, 1 m
SP1500 m/s
SR60 kSa/s
Absorption coefficient, a1 dBm−1
dispersion factor, k0
b19, 39, 59, 79

Table 1.

Parameters used in Matlab simulation.

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 bare as results of varying dDBS in the four different Figures 4(a)–(d). The other parameters are same for all the figures.

Figure 4.

Number of fish and mammals vs.σof CCF (a)b = 19 (dDBS = 0.25 m andSR = 60 kSa/s) (b)b = 39 (dDBS = 0.5 m andSR = 60 kSa/s), (c)b = 59 (dDBS = 0.75 m andSR = 60 kSa/s), and (d)b = 79 (dDBS = 1 m andSR = 60 kSa/s).

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, bhas 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 band 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.

Figure 5.

Exact number of fish and mammals vs. estimated number of fish and mammals forb = 79 (dDBS = 1 m andSR = 60 kSa/s).

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.

© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Shaik Asif Hossain and Monir Hossen (July 13th 2020). Cross-Correlation-Based Fisheries Stock Assessment Technique: Utilization of Standard Deviation of Cross-Correlation Function as Estimation Parameter with Four Acoustic Sensors, Underwater Work, Sérgio António Neves Lousada, IntechOpen, DOI: 10.5772/intechopen.93240. Available from:

chapter statistics

88total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Diving as a Scientist: Training, Recognition, Occupation - The “Science Diver” Project

By Alexandros Tourtas, Kimon Papadimitriou, Elpida Karadimou and Ralph O. Schill

Related Book

First chapter

Introductory Chapter: Introduction to Sea Level Rise and Coastal Infrastructure

By Yuanzhi Zhang, Kapo Wong, Yijun Hou and Xiaomei Yang

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us