Open access peer-reviewed chapter - ONLINE FIRST

# 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

## 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 (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 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 Sp of sound is considered [23]. Figure 3 shows an example of 3D estimation area under water space with a single fish N1and 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:

dDBS12=x1x22+y1y22+z1z22E1
dDBS23=x2x32+y2y32+z2z32E2
dDBS34=x3x42+y4y42+z3z42E3

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:

Sr11t=α11S11tτ11,E4
Sr12t=α12S12tτ12,E5
Sr13t=α13S13tτ13,E6
Sr14t=α14S14tτ14,E7

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:

C1τ=+S11tS12tτ11dτE8
C2τ=+S12tS13tτ12dτE9
C3τ=+S13tS14tτ13dτE10

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 H1, H2, H3, and H4 are:

Srt1=j=1Nαj1Sjtτj1E11
Srt2=j=1Nαj2Sjtτj2E12
Srt3=j=1Nαj3Sjtτj3E13
Srt4=j=1Nαj4Sjtτj4E14

Therefore, the total CCFs are:

C12τ=+Srt1tSrt2tτdτE15
C23τ=+Srt2tSrt3tτdτE16
C34τ=+Srt3tSrt4tτdτE17

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 Dascyllus aruanus 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]:

Xt=Acos2πf2f1t22d+f1t+PE18

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, v is 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]:

VarCt=C2tCt2,E20

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

Ct2=QT2v2Tr2+Tr2dt+dr1+tdτ1Gr1ratτ1×Gr1rbτ+tτ1×Tr2+Tr2dt~+dr3+tdτ3Gr3ratτ3×Gr3rbτ+t~τ3E21

and

C2t=QTv2Tr2+Tr2dt+dr1+tdτ1Gr1ratτ1×Gr1rbτ+tτ1×Tr2+Tr2dt~+dr3+tdτ3Gr3ratτ3×Gr3rbτ+t~τ3+QTv2Tr2+Tr2dt+dt~+dr1+tdτ1Gr1ratτ1×Gr1rat~τ1×+dr1+tdτ1Gr1ratτ1×Gr1rat~τ1+QTv2Tr2+Tr2dt+dt~+dr1+tdτ1Gr1ratτ1×Gr1rbt~τ1×+dr2+tdτ2Gr2rbtτ2×Gr2rbt~τ2,E22

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.

σ=VarCt=C2tCt2E23

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]:

σ=N×1b×11b,E24

where N is the number of fish and mammals and b is the number of bins. Here, b can be achieved from the following Eq. [22]:

b=2×dDBS×SRSP1,E25

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:

N=b2×σ2b1E26

Therefore, if σ is available from simulation, the estimated population size of fish and mammals, N will 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.

σAverage3CCF=σ1+σ2+σ33E27

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

N=b2×σAverage3CCF2b1E28

Therefore, if σ is available from simulation, N will 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.

ParametersValue
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 b are as results of varying dDBS 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, b 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.

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.

## How to cite and reference

### 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 [Online First], IntechOpen, DOI: 10.5772/intechopen.93240. Available from: