1. Introduction
Acoustic methods are widely used for estimating fish abundance [
1,
2,
3,
4,
5,
6,
7,
8,
9,
10,
11,
12,
13,
14,
15,
16,
17,
18,
19,
20,
21,
22,
23,
24,
25,
26,
27,
28], and constitute a key part of the analytic assessment that makes the basis for international regulations of marine resources. For fish aggregated in schools or layers, echo integration [
7,
8] supported by biological sampling, is the most common method used in oceanic surveys [
3]. Fish abundance is measured using narrowband signals typically in the 10–120 kHz range, whereas species identification of fish and quantification of zooplankton are based on multi-frequency signals extended to 500 kHz or higher. The acoustic methods rely on calibrated systems [
9,
10] and power-budget equations.
In oceanic surveys, a power-budget equation for multi-target (volume) backscattering [
6,
11,
12,
13,
14,
15,
16] is typically used to measure volume backscattering from aggregations of fish or zooplankton. In terms of this equation, the volume backscattering coefficient
is measured for a sequence of thin spherical shell “ping volumes”,
, at increasing range. The sequence of
measurements is integrated over the range of an observation volume
[
6,
11,
16], to give the fish density in
,
[
17], in terms of an echo-integrator equation [
3,
6,
8,
10,
16].
Prior to survey operation, a related power-budget equation for single-target backscattering [
6,
11,
12,
13,
14,
15,
16] is used for at-sea calibration of the echo sounder using a standard target [
3,
6,
9,
10]. The same power-budget equation for single-target backscattering is used to measure the backscattering cross section,
, or equivalently, the target strength,
TS, of individual fish.
Different formulations of power-budget equations and calibration factors have been presented for use in scientific echo sounder and sonar systems.
Clay and Medwin [
1], and later Medwin and Clay [
5], derived expressions corresponding to
that account for sound propagation in the sea, in terms of echo integration of the free-field acoustic pressure signals at the centre of the transmitting and receiving transducer’s front surface [
16].
Dalen and Nakken [
8] proposed an empirical echo-integrator equation for estimation of
, involving a calibration factor to be determined by at-sea calibration using a standard target (e.g., a metal sphere). The expression was modified by Foote et al. [
10], cf. also [
3,
4,
6].
Simrad [
11] extended Clay and Medwin’s [
1,
5] analyses by accounting for the transmitting and receiving transducer responses. Power-budget equations for
and
were given on average power form, in addition to expressions postulated for two calibration factors used in the Simrad EK500 echo sounder system [
11]. The formulation applies to certain conditions of electrical termination of the transducer [
16] (cf.
Section 4.1). The Simrad formulation constitutes an important basis for the methods used in fish abundance measurement today. Simrad’s derivation of the expressions for
and
was summarized by Korneliussen [
12].
Demer and Renfree [
14] also used expressions for
and
that correspond to those given in [
11]. Related power-budget equations were postulated by Simmonds and MacLennan [
6] (see their Equations (3.13) and (3.15)). However, when these are re-arranged to yield expressions for
(the expected value of
) and
, the resulting expressions may be shown to differ somewhat from those given in [
11], and power flow balance is not quite preserved.
In a more detailed derivation, Pedersen [
13] extended Simrad’s [
11] analysis by—in the power-budget equations for
and
—including a factor accounting for arbitrary electrical termination of the transducer.
Lunde et al. [
16] further developed the power-budget equations that are given on average power form by [
11,
13], to account for echo integration.
,
, and
were expressed in terms of time integration of the transmitted and received voltage signals (echo integrals). The resulting expressions were shown to represent a consistent generalization of the echo integration expressions derived by Clay and Medwin [
1,
5] for “in-water” pressure signals.
For the more recent generation of Simrad scientific and fisheries echo sounder and sonar systems, such as EK60, ES60, ME70 and MS70, an apparently different set of power-budget equations is used [
15]. Ona et al. [
15] postulated equations for
and
that are used in these systems [
18,
19,
20,
21]. A derivation of similar power-budget equations and calibration factors as used in [
11,
15,
18,
19,
20,
21] was given by Lunde and Korneliussen [
22]. An alternative derivation of the
and
expressions being postulated in [
15] was presented in [
23].
There are unresolved questions and some apparent inconsistencies related to the prior literature in this field.
Firstly, in [
11], expressions have been postulated for two calibration factors, “
TS transducer gain” and “
transducer gain”, being determined by echo sounder calibration. However, neither definition, nor the mathematical relationship between these calibration factors and the power-budget equations for
and
given in [
11], are available from the literature. That is, no power-budget equations for
and
, expressed in terms of the two calibration factors, are given. It can be shown [
22] that the expressions for the two calibration factors postulated in [
11] cannot be readily derived from the average power formulation of the power-budget equations given in [
11].
Secondly, ref. [
11] states that “the TS-measurement is based on the peak value of the echo samples in the sphere echo, whereas the
-measurement is based on integration (averaging) of the echo samples”. This difference between
and
measurements is not reflected in the power-budget equations given in [
11]. This situation strongly indicates that the power-budget equations for
and
that are actually used in EK500, and which are presumed to be expressed in terms of the two EK500 calibration factors, may not be exactly the same as those given in [
11].
Thirdly, the expression for
postulated by [
15] differs from those given by other authors [
6,
11,
13,
14,
16], and is not consistent with the expression given e.g., in [
11]. In particular this concerns a parameter
introduced by [
15] in their logarithmic expression for
[
22]. The ratio of
as given by [
11] and [
15] is not equal to 1, but
. It follows [
22] that the power-budget equations as given in [
15] cannot be derived directly from the power-budget equations given in [
11].
Fourthly, the set of power-budget equations given in [
15] involves two calibration factors, “Axial transducer gain” and “Integration correction”, being determined by echo sounder calibration. In [
18,
19,
20,
21] these are referred to as “Gain” and “
Sa correction”. In the literature available for these instruments [
15,
18,
19,
20,
21], the calibration factors are not defined, nor given by any expressions.
Fifthly, the derivation given in [
23] of the power-budget equations for
and
being postulated in [
15], appears to be incomplete on essential points, such as with respect to the
parameter. In [
23],
is introduced ad hoc by replacing the nominal pulse duration,
(in terms of the terminology used here), with an effective pulse duration, defined as
. The mathematics and arguments used are not clear, and the expression for
given by [
23] deviates from those given by other authors [
6,
11,
13,
14,
16].
The unclear situation connected to (i) the deviating and apparently inconsistent expressions given in the literature in this field; and (ii) the sparse documentation on the power-budget equations and calibration factors actually employed in widely used echo sounder and sonar systems, has caused uncertainty and confusion among users. The equations used for international regulations of marine resources, including their theoretical basis, should preferably be known and documented in available literature. Moreover, the equations should be sufficiently complete to enable accurate and reliable analyses of error and system drift.
The objective of the present paper is to formulate a generic (instrument independent) and unifying theory of the power-budget and echo-integrator equations for , , and , that is capable of explaining the different power-budget and calibration factor formulations used in the literature, and how these different formulations are mutually related. By “unifying” is meant “giving the connection and relationships between different (and apparently diverging) formulations”.
Table 1 gives an overview of the five formulations addressed here. Based on the average power formulation of the
and
equations given in [
13,
16] (denoted “Formulation A”), an echo integrator formulation is derived for narrowband signals (denoted “Formulation B”). Expressions are presented for a calibration factor that can be measured using current methods for at-sea echo sounder calibration employing a metal sphere. As an alternative,
and
are also expressed directly in terms of the basic quantities being measured in calibration. The corresponding expression for
is derived, applicable to both alternatives.
As an application of this generic theory, alternative power-budget and calibration factor formulations are derived for echo sounders based on using “peak voltage echo integration” in single-target
TS measurements, instead of the full echo waveform itself [
11,
18,
19,
20,
21]. Three equivalent formulations are derived, denoted “Formulation C”, “D”, and “E”, respectively. These are all equivalent to Formulation B, and fully valid approaches, for instruments using this signal processing strategy.
The resulting expressions are shown to be closely related to the equations given for two widely used and important families of scientific echo sounder and sonar systems; (i) the Simrad EK500 [
11] and (ii) the more recent Simrad EK60, ES60, ME70 and MS70 [
15,
18,
19,
20,
21] systems. By using the equivalent Formulations C, D, and E, the deviations and apparent inconsistencies between the expressions given in refs. [
15,
23] and those of refs. [
11,
13,
16] seem to be explained and resolved. This includes the
parameter, as well as definitions of, relationships between, and explanation for use of the various calibration factors involved in the two echo sounder systems, and how these are mathematically related to the power-budget equations of the respective instruments.
Formulation B is derived and employed here as the “generic and unifying theory” since it can be used to derive the connections and relationships between Formulation A and Formulations C, D, and E, including explanation of the different calibration factors. It represents a necessary key approach to derive, explain, and understand Formulations C–E. That is, to reveal and explain the formulations and expressions underlying the Simrad EK500, EK60, etc., scientific echo sounder and sonar systems. In Formulation B, full waveform echo integration (here denoted [
tivs], cf.
Section 2.3.1) is used for both
and
measurements (i.e., in calibration and field measurements), since that represents a key approach to reveal and explain the diverse echo integration methods and calibration factors that are in use in modern scientific echo sounder and sonar systems. The objective is thus not to find a “best” or preferred method or formulation of power budget equations and calibration factors for fish abundance estimation. The intention and objective is to address the challenges encountered and inconsistencies identified in the literature on modern scientific echo sounder and sonar systems, and present a theory that explains and resolves these inconsistencies.
The present paper represents a condensed, partly reformulated, and extended version of a technical report [
22], in which further details on the mathematical derivations can be found. Extensions of [
22] relate to placement of the theory and results e.g., in the context of recent literature, and a more comprehensive derivation of the echo integration formulation (
Section 2.3.1).
The main contribution of the present paper is the generic and unifying theory for echo sounders operated in their linear (small-amplitude) range (Formulation B), used to reveal the relationship between, and to correct, some diverging and inconsistent expressions given in prior literature [
11,
15,
18,
19,
20,
21,
23].
4. Results and Discussion
Four equivalent echo integration Formulations B–E of the power-budget equations for
and
with calibration factors have been derived and described. Formulation B is derived to provide a generic echo integration formulation of the functional relationship for abundance estimation. Formulations C–E have been developed from Formulation B, to investigate the consequences of using “peak voltage echo integration” in single-target
TS measurements, as used by some manufacturers [
11,
15,
18,
19,
20,
21].
In the following, Formulations A and B are briefly discussed, and Formulations C–E are compared to the expressions given for two important and widely used families of echo sounder and sonar systems; (i) the Simrad EK500 [
11] and (ii) the Simrad EK60, ES60, ME70, and MS70 [
15,
18,
19,
20,
21] systems.
4.1. Formulation A (Average Power)
Formulation A, Equations (1) and (2), derived by [
13,
16], extends the power-budget equations given in [
11] by accounting for arbitrary electrical termination [
13,
16], represented by an electrical termination factor,
.
The expressions given in [
11] correspond to setting
in Formulation A, which is valid for electrical termination conditions for which
(conjugate matched electrical termination), or for
when
= 0 (i.e., in a frequency band close to the series resonance frequency of the transducer vibration mode used) [
16]. Otherwise, formulation A is equal to the expressions given for
and
in [
11].
As Formulation A, given in terms of average power, is not directly suited for signal processing, it serves here as the basis for deriving the four echo integration formulations, B–E.
4.2. Formulation B (Echo Integration, Generic)
Formulation B, Equations (23) and (25) for
and
, together with Equations (29) or (31) for the calibration factor
, is derived from Formulation A to provide expressions better suited for practical signal processing in echo sounders [
16]. This formulation in terms of echo integrals, for calculation of the electrical powers involved, leads to an alternative electrical termination factor in the power-budget equations,
. It also involves several integration time intervals,
,
,
, and
, to differentiate between echo integration of the transmitted voltage signal; and echo integration of the single-target (fish), calibration sphere, and volume backscattering voltage echoes, respectively. Formulation B involves only a single calibration factor,
. The generic functional relationship Formulation B constitutes the basis for deriving Formulations C–E, applicable to “peak voltage echo integration” in single-target
TE measurements.
4.3. Formulation C (Echo Integration, “EK500 Type”)
Formulation C, given by Equations (41) and (42) for
and
, together with Equations (45)–(51) for the two calibration factors
and
, appears to correspond to, and generalize, the expressions employed in the Simrad EK500 echo sounder [
11]. Formulation C is thus for convenience referred to as an “EK500 type” of power-budget equations. This may be seen as follows.
Firstly, the present theory reveals that the expressions for
and
that are actually used in EK500 cannot be exactly those given in its manual [
11]. The latter expressions correspond to the average power Formulation A (with
), whereas the EK500 implementation involves full waveform echo integration for
, “peak voltage echo integration” for single-target
TS measurements (
), and two calibration factors “
TS transducer gain” and “
transducer gain” [
11], which neither appear in, nor can be obtained directly from, Formulation A. The echo integration actually used by the manufacturer in the EK500 echo sounder system was never documented in the manual [
11] or elsewhere in the available literature. The present theory aims to reveal and explain expressions that appear to be similar to those having been implemented but not documented.
Secondly, by accounting for full waveform echo integration for
, and “peak voltage echo integration” for single-target
TS measurements (
) in Formulation C, expressions are here obtained that involve the calibration factors
and
. By defining
TS transducer gain and
transducer gain , and using Equation (51), Equations (47) and (50) become identical to the corresponding expressions postulated for the EK500 calibration factors “
TS transducer gain” and “
Sv transducer gain” in [
11]. Formulation C is a prerequisite to obtain these expressions.
Consequently, it is shown here that Formulation C leads to expressions for the calibration factors that are identical to those given for the EK500 echo sounder [
11]. As these calibration factors cannot be derived directly from the average power budget expressions for
and
given in [
11], the analysis shows that the expressions for
and
actually implemented and used in the EK500 system, cannot be those given in [
11]. The fact that Formulation C consistently—through the echo integration Formulation B—gives the connection between the average power expressions for
and
given in [
11] and these calibration factors, including description of “peak voltage echo integration” in
TS measurements, strongly indicates that Formulation C may represent a functional relationship applicable to the EK500 system.
The parameter
has been introduced here to account for “peak voltage echo integration” in single-target
TS measurements. No such parameter as
is used in refs. [
11,
15,
18,
19,
20,
21,
23].
4.4. Formulation D (Echo Integration, “EK60 Type”)
Formulation D, given by Equations (52) and (53) for
and
, together with Equations (45)–(47) and (58)–(62), respectively, for the two calibration factors
and
involved, appears to correspond to, and generalize, the expressions used for the Simrad EK60, ES60, ME70 and MS70 systems [
15,
18,
19,
20,
21]. For convenience, this formulation is thus referred to as an “EK60 type” of power-budget equations. This may be seen as follows.
Firstly, as explained in
Section 1, the power-budget equations as these are given in [
11] and [
15] are not consistent. The ratio of the two different expressions for
as given by [
11] and [
15], respectively, is not equal to 1, but
. This indicates that the expressions for
and
actually being used in EK60 and related instruments may be different from those given in [
15]. Formulation D is used here to clarify this issue. By accounting for “peak voltage echo integration” for single-target
TS measurements in Formulation D, expressions are here obtained that are fully consistent both with [
11] and Formulation C, and account for
in a consistent way. In Formulations C and D, the expressions for
and
are equivalent, and the ratio of the two expressions for
equals 1.
Secondly, in Formulation D, expressions are obtained that involve the calibration factors
and
. By defining
Gain and
correction , Equations (47) and (62) become equal to expressions used in EK60 [
32]. Moreover, the expressions derived here appear to be consistent with the available information given on the “Gain” and “
correction” calibration factors in the manuals of EK60, etc. [
18,
19,
20,
21].
The expressions postulated by Ona et al. [
15] for
and
correspond to Formulation D. There are, however, some deviations. Firstly (in terms of the terminology used here), in [
15], electrical powers
and
(in units of W) were used instead of echo integrals,
,
, and
(in units of V
2s = W·ohm·s = J·ohm). Echo integration was thus not accounted for in these expressions. Secondly, the same
was used for
and
, which is not consistent with use of “peak voltage echo integration” and full waveform echo integration, respectively, for
TS and
(and thus
) measurements. Thirdly,
and
were used instead of
and
in the expressions for
and
, respectively, which is inconsistent when “peak voltage echo integration” is used for
, and
is used in
. Finally, in [
15], electrical termination was not accounted for, implicitly implying
[
13,
16], cf.
Section 4.1.
In other words, for power-budget equations expressed in terms of echo integrals, and for which “peak voltage echo integration” is used in
TS measurements instead of full waveform echo integration, the parameters
and
, referred to in [
15] as “transducer gain” and “on-axis transducer gain”, respectively, are not equal to the transducer gain and axial transducer gain,
and
G0, that are involved in the traditional power-budget equation given in [
11], as implicitly stated in [
15]. Instead, they represent “effective” transducer and axial transducer gains, respectively,
Geff(θ, φ) and
G0,eff, with a correction factor
sa,corr involved, to compensate for the use of “peak voltage echo integration” in
TS measurements instead of full waveform echo integration, cf. Equations (52) and (53).
Through the derivation of Formulation D, expressions have been obtained that (a) seem to explain and resolve the deviation between the power-budget equations given in [
15] and [
11]; (b) explain the introduction and use of the
s parameter; and (c) “harmonize the
TS and
measurements” [
19] by employing the same calibration factor
G0,eff in σ
bs and
. The power-budget expressions that are consistent with [
11], and that may replace those given in [
15], appear to be Equations (52) and (53). This includes calibration factors
and
, given by Equations (45)–(47) and (58)–(62), respectively.
Through the present analysis, the relationship between the Formulation D calibration factors and , and the Formulation C calibration factors, and , is established, cf. Equation (61). This also includes their relationships to the generic Formulations A and B types of description. It follows that the relationships between the EK500 calibration factors “TS transducer gain” and “Sv transducer gain”, and the EK60 etc. calibration factors “Gain” and “Sa correction”, also appear to have been established and explained.
As the expressions for
and
derived in [
23] correspond exactly to those postulated in [
15], the above discussion also applies to [
23].
4.5. On the , , and Parameters
The reasons for introducing the parameter , , and in Formulations C–E, and the need for using two calibration factors for each of these formulations, are explained in the following.
As Formulation B is based on full waveform echo integration both for single-target TS measurements () and in survey operation (), only a single calibration factor, , applies, with no need for , , and type of parameters.
In Formulation C, due to the use of “peak voltage echo integration” instead of full waveform echo integration in single-target TS measurements, an error is introduced in , which is compensated for by introducing instead of in , cf. Equations (41), (43), and (44). By still using in the expression for , cf. Equation (42), two calibration factors and thus become involved in Formulation C, referred to as “TS transducer gain” and “Sv transducer gain”, respectively.
In Formulation D, the calibration factor used in Formulation C (in the expression) is replaced by , to “harmonize” the and equations. An error is then introduced in , which is compensated for by using in , referred to as the “Sa correction”, cf. Equation (53).
In Formulation E, the term used in Formulation D (in the equation) is denoted . As a result, the and expressions of Formulation E are on the same form as in Formulation B, except that different parameters apply to the two formulations. In Formulation B these are , , and . In Formulation E these are , , and , due to the “peak voltage echo integration” used for single-target TS measurements.
In [
15] and [
23],
is referred to as an “effective pulse duration”. This interpretation may be discussed. In the present analysis,
, used to define
in Equation (65), is shown to be essentially a correction factor to compensate for the simplified “peak voltage echo integration” of the received echo in single-target
TS measurements, caused by the use of
instead of
for these targets. Hence, from Equation (65) it appears that
represents a
scaled gate opening time, to compensate for erroneous echo integration at
signal reception for single-target
TS measurements. A more representative interpretation of
might thus be “effective echo integration time”.
4.6. Comments in Relation to Conventional Operation
It is emphasized that, from the available literature and practical use in fish abundance estimation, there is no reason to question the expressions actually
implemented and used by the manufacturer in the mentioned instruments [
11,
18,
19,
20,
21]. The issues addressed here are all related to the
explanation and documentation of the expressions given in the available literature.
It is also emphasized that the derivation and presentation of Formulations B–E by no means indicates that any new method for abundance estimation is proposed, as an alternative to the methods used today. The situation is quite the opposite: The expressions presented here, which are all equivalent, are intended to provide a consistent and unifying theory for improved understanding and control in use of the conventional methods, when employing commonly used echo sounder and sonar systems.
In addition, the more general and complete expressions derived here may constitute an improved basis for evaluation and, if necessary, correction of errors in abundance estimation and species determination. One issue in this respect may be possible system drift, due to changing environment or echo sounder parameters from calibration to oceanic surveying [
14]. By knowing the full functional relationship for the abundance measurement, in terms of echo sounder and environmental parameters for the measurement system in question, the possibilities for reliable error evaluation and correction improve.
5. Conclusions
On the basis of an average power formulation of power-budget equations for
and
known from prior literature [
13,
16] (denoted Formulation A), four new echo integration formulations for
and
and associated calibration factors are derived for narrowband signals, denoted Formulations B–E. These are given in terms of the quantities being measured or assumed in calibration and oceanic surveying.
Formulation B gives a generic (instrument independent) functional relationship for fish abundance estimation, and represents a main result of the present work. When waveform integration is used to calculate the echo integrals both in single-target TS measurements and in survey operation, the generic Formulation B serves to be convenient. There is then no need for , , and type of parameters, and only a single calibration factor applies.
Formulations C–E apply to instruments being based on “peak voltage echo integration” in single-target TS measurements (for calibration, or individual fish). In these formulations, the respective parameters , , or are introduced to provide consistent and valid functional relationships, and two calibration factors are needed for each formulation. It is a matter of preference which of Formulations C–E is used. Formulations B–E are all equivalent and fully valid functional relationships.
The theory presented here appears to resolve deviations and apparent inconsistencies between various power-budget equations published in prior literature [
11,
15,
23]. It is shown that for the
and
expressions given in [
15,
23] to be consistent with [
11], some of the terms involved have to be replaced, to reflect the different “peak voltage echo integration” and full waveform echo integration strategies used for
and
, respectively. In terms of the terminology used here, the average electrical powers
and
used by [
15,
23] are to be replaced by echo integrals
,
, and
. Similarly, the transducer gains
and
used by [
15,
23] are to be replaced by
and
. With these modifications applied to the expressions given in [
15,
23], the expressions given by Formulations B–D for
and
are all consistent with those given in [
11,
13,
14,
15,
16,
23]. In terms of this theory, a clear rationale is given for the introduction and use of parameters such as
and
, which were postulated by [
15] and adopted by [
23].
and
are here consistently derived and explained.
By comparison with expressions published in available literature it is shown that Formulations C and D lead to expressions for
,
, and calibration factors that appear to correspond to and explain those given for the Simrad EK500 [
11], and the Simrad EK60, ES60, ME70, and MS70 systems [
15,
18,
19,
20,
21], respectively. The respective calibration factors involved are here precisely defined. Through Formulations B–C, the mathematical relationship between these calibration factors and the power budget equations given in [
11] is established. The power-budget equations for
and
given on echo integration form are expressed directly in terms of the respective calibration factors.
The consistency in the results obtained here indicates that the power-budget equations for
and
employed in the EK500 and the EK60 etc. systems are not exactly those given in [
11,
15], but rather have forms corresponding to Formulations C and D, respectively. However, the lack of documentation of the power-budget equations and calibration factors actually used in the EK500, EK60, ES60, ME70, and MS70 systems unfortunately prevents a complete comparison of the expressions derived here for Formulations C and D with those used in these systems.
Under the assumption of small-amplitude (linear) sound propagation, the expressions derived here are expected to represent a consistent and relatively complete theoretical basis for improved understanding and control in use of conventional methods and instruments for fish abundance measurement and species identification. Prior literature is extended to provide more complete power-budget equations by accounting for echo integration, electrical termination, and the full range of electrical and acoustical echo sounder parameters. The results are expected to constitute an improved theoretical fundament for measurement, error evaluation, possible error compensation, and uncertainty evaluation of fish abundance methods and equipment in use today.