Next Article in Journal
Storm Surge Modeling in Large Estuaries: Sensitivity Analyses to Parameters and Physical Processes in the Chesapeake Bay
Next Article in Special Issue
Structure-Borne Wave Radiation by Impact and Vibratory Piling in Offshore Installations: From Sound Prediction to Auditory Damage
Previous Article in Journal
Growth-Prediction Model for Blue Mussels (Mytilus edulis) on Future Optimally Thinned Farm-Ropes in Great Belt (Denmark)
Previous Article in Special Issue
On the Sediment Dynamics in a Tidally Energetic Channel: The Inner Sound, Northern Scotland
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Power-Budget Equations and Calibration Factors for Fish Abundance Estimation Using Scientific Echo Sounder and Sonar Systems

by
Per Lunde
1,3,* and
Rolf J. Korneliussen
2,3
1
Department of Physics and Technology, University of Bergen, P.O. Box 7803, N-5020 Bergen, Norway
2
Institute of Marine Research, P.O. Box 1870, Nordnes, N-5817 Bergen, Norway
3
Christian Michelsen Research AS, P.O. Box 6031, Postterminalen, N-5892 Bergen, Norway
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2016, 4(3), 43; https://doi.org/10.3390/jmse4030043
Submission received: 29 February 2016 / Revised: 25 June 2016 / Accepted: 4 July 2016 / Published: 9 July 2016
(This article belongs to the Special Issue Underwater Acoustics and the Ocean Environment)

Abstract

:
Acoustic methods used in fish abundance estimation constitute a key part of the analytic assessment that makes the basis for abundance estimation of marine resources. The methods rely on power-budget equations and calibrated systems. Different formulations of power-budget equations and calibration factors have been proposed for use in scientific echo sounder and sonar systems. There are unresolved questions and apparent inconsistencies in prior literature related to this field. A generic (instrument independent) and unifying theory is presented that attempts to explain the different power-budget and calibration factor formulations proposed and used in prior literature, and how these are mutually related. Deviations and apparent inconsistencies in this literature appear to be explained and corrected. This also includes different (instrument specific) formulations employed in important modern scientific echo sounder systems, and their relationship to the generic theory of abundance estimation. Prior literature is extended to provide more complete power-budget equations for fish abundance estimation and species identification, by accounting for echo integration, electrical termination, and the full range of electrical and acoustical echo sounder parameters. The expressions provide a consistent theoretical basis for improved understanding of conventional methods and instruments used today, also enabling improved sensitivity and error analyses, and correction possibilities.

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 s v is measured for a sequence of thin spherical shell “ping volumes”, V p , at increasing range. The sequence of s v measurements is integrated over the range of an observation volume V o b s [6,11,16], to give the fish density in V o b s , ρ a [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, σ b s , 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 s v 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 ρ a , 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 σ b s and s v 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 σ b s and s v was summarized by Korneliussen [12].
Demer and Renfree [14] also used expressions for σ b s and s v 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 σ b s (the expected value of σ b s ) and s v , 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 σ b s and s v —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. σ b s , s v , and ρ a 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 σ b s and s v 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 σ b s and s v 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 “ S v 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 σ b s and s v given in [11], are available from the literature. That is, no power-budget equations for σ b s and s v , 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 s A -measurement is based on integration (averaging) of the echo samples”. This difference between σ b s and s v measurements is not reflected in the power-budget equations given in [11]. This situation strongly indicates that the power-budget equations for σ b s and s v 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 s v 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 S a , c o r r 10 log ( s a , c o r r ) introduced by [15] in their logarithmic expression for s v [22]. The ratio of s v as given by [11] and [15] is not equal to 1, but s a , c o r r 2 . 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 σ b s and s v being postulated in [15], appears to be incomplete on essential points, such as with respect to the s a , c o r r parameter. In [23], s a , c o r r is introduced ad hoc by replacing the nominal pulse duration, τ p (in terms of the terminology used here), with an effective pulse duration, defined as τ e f f τ p s a , c o r r 2 . The mathematics and arguments used are not clear, and the expression for s v 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 σ b s , s v , and ρ a , 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 σ b s and s v 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, σ b s and s v are also expressed directly in terms of the basic quantities being measured in calibration. The corresponding expression for ρ a 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 S a , c o r r 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 σ b s and s v 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].

2. Theory

2.1. Assumptions

The analysis is based on the following assumptions for the echo sounder and its environment (cf. [16] and references therein):
(a)
the monostatically operated transducer is passive, reversible, and reciprocal;
(b)
the transmit voltage amplitude is sufficiently small to avoid nonlinear effects in the electroacoustic transducer and electronics (i.e., the transducer and electronics are operated in their linear ranges);
(c)
the fluid medium (seawater) is homogeneous, with constant density and sound velocity;
(d)
the amplitudes of the transmitted sound pressure signals are sufficiently small so that finite-amplitude sound propagation effects in seawater can be neglected;
(e)
targets are in the far field of the transducer;
(f)
possible nonlinear effects in the scattering process at the target itself (involving e.g., fish with gas-filled swim-bladder), can be neglected, so that the linear backscattering theory applies;
(g)
the volume backscattering coefficient can be calculated as a sum of backscattering cross sections (i.e., intensities) per unit volume;
(h)
the scattering objects are uniformly distributed in the observation volume, with
(i)
random phases of the scattered echoes (i.e., random spacing of scattering objects, and movement of objects from one transmission to the next);
(j)
possible multiple-scattering effects and interaction between objects are neglected;
(k)
excess attenuation from power extinction caused by volume scattering is neglected; and
(l)
the same transmit electrical power Π T is used for σ b s and s v measurements, i.e., in calibration and surveying.
These are all common assumptions underlying the traditional theory of fish abundance measurement [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16]. The discussion of their validity is an extensive and complex subject, beyond the scope of the present work, and discussed elsewhere. Assumptions (a) and (b) relate to the transducer and electric components of the echo sounder system, and are normally fulfilled by driving the piezoelectric transducer using low electrical power. In relation to (c): in abundance estimation, the sound velocity is typically taken to be the average value of the sound velocity profile, over the depth range in question [24]. Assumption (d) is addressed by refs. [13,24,25,26], and maximum electrical transmission powers have been suggested [24]. Assumption (f) is discussed e.g., in ref. [27]. The assumptions (g)–(k) are included in the set of assumptions used by Clay and Medwin [1,5] to derive the analogous “in-water” expressions for s v , accounting for acoustic pressures in the sea only. Relatively extensive discussions on the validity of (g)–(k) are given by refs. [1,5,6,28], also summarizing other studies addressing these issues.
A spherical coordinate system is used, with coordinates ( r , θ , ϕ ) , origin at the centre of the transducer front, and with the z axis (i.e., θ = ϕ = 0 ) chosen normal to the transducer’s front surface, and assumed coincident with the transducer’s acoustical beam axis [16]. r is the radial distance, denoted range, θ is the polar angle (rel. to the z axis), and ϕ is the azimuthal angle (rel. to the x axis). Bold-face letters are used to indicate complex-valued quantities.

2.2. Average Power Formulation (“Formulation A”)

Under the above assumptions, it can be shown from basic acoustic principles that the backscattering cross section of a single scattering target located at position ( r , θ , ϕ ) in the transducer’s far field, σ b s , and the volume backscattering coefficient for a thin spherical shell “ping volume” V p in the far field, s v , are given by the power budget equations [11,13,16]:
σ b s = 16 π 2 r 4 e 4 α r Π R s t G 2 ( θ , ϕ ) λ 2 F Π Π T
s v = 32 π 2 r 2 e 4 α r Π R v G 0 2 ψ λ 2 c 0 τ p F Π Π T ,
respectively. σ b s and s v are given in units of m2 and m−1, respectively. Here, τ p is the duration of the transmitted voltage signal (denoted “pulse duration”). In Equation (1), r is the distance to the single target. In Equation (2), r r p ( r p 1 + r p 2 ) / 2 is the mid-range of the “ping volume” V p , that is contained within ranges r p 1 and r p 2 , with thickness d r p r p 2 r p 1 = 1 2 c 0 τ p [16]. In practice, τ p = 1 ms may often be used, corresponding to d r p = 0.75 m. α is the sound pressure acoustic absorption coefficient of seawater (expressed in Np/m). λ = c 0 / f is the acoustic wavelength, where c 0 is the small-signal sound velocity in seawater.
G ( θ , ϕ ) and G 0 (both non-dimensional) are the transducer gain and the axial transducer gain, defined as [11,13,16,29]:
G ( θ , ϕ ) η D ( θ , ϕ ) = G 0 | B i ( θ , ϕ ) | 2
G 0 G ( 0 , 0 ) = η D 0 ,
respectively, where η is the transducer’s electroacoustic conversion efficiency (the non-dimensional ratio of radiated acoustic to transmitted electrical power), and
D 0 = 4 π 4 π | B i ( θ , ϕ ) | 2 d Ω
D ( θ , ϕ ) = 4 π | B i ( θ , ϕ ) | 2 4 π | B i ( θ , ϕ ) | 2 d Ω = D 0 | B i ( θ , ϕ ) | 2
and B i ( θ , ϕ ) , are the axial directivity factor [30], the directivity factor [29], and the beam pattern (the angular distribution of the sound pressure, normalized to the axial sound pressure) [30], respectively, for the transmitted sound pressure field.
ψ (in steradians, sr) is the equivalent two-way beam solid angle of the transducer, defined by [1,5,10,11,16,29]:
ψ 4 π | B i ( θ , ϕ ) | 4 d Ω = 1 G 0 2 4 π G 2 ( θ , ϕ ) d Ω .
ψ is normally provided by the echo sounder manufacturer, for the frequencies in question.
The average transmitted and received electrical powers at the transducer terminals (W) are given as [16,30]:
Π T = | V T | 2 R T 2 | Z T | 2
Π R s t = | V R s t | 2 R E 2 | Z E | 2
Π R v = | V R v | 2 R E 2 | Z E | 2 ,
respectively, where “average” refers to averaging over one cycle of the harmonic signal waveform, at the frequency f. V T , V R s t , and V R v are the spectral amplitudes of the voltage signals across the transducer’s electrical terminals during transmission and reception, i.e., for the transmitted signal and the single-target and volume backscattered echoes, respectively. Subscripts “T” and “R” indicate “transmit” and “receive”, and superscript “st” and “v” indicate “single target” and “volume backscattering”, respectively. For Equations (1), (2), and (8)-(10) it has been assumed that the same transmit signal V T ( t ) and integration time τ p are used in single-target (e.g., calibration) and oceanic survey operations, respectively.
Z T = R T + i X T is the transducer’s input electrical impedance at transmission, when radiating into the fluid, with resistance and reactance RT and XT, respectively. ZR = RR + iXR is the output (internal) electrical impedance of the receiving transducer, and Z E = R E + i X E is the input electrical impedance of the receiving electric network (the transducer’s “electrical termination impedance”), involving resistances and reactances R R and X R , and R E and X E , respectively. (The roles of these electrical impedances in electrical equivalent circuits of the transmitting and receiving parts of the measurement system are indicated in Figures 3 and 4 of ref. [16].) Z T and Z E should normally be available from the echo sounder manufacturer. The assumption Z R = Z T is often used.
Finally,
F Π 4 R T R E | Z R + Z E | 2
is a (non-dimensional) electrical termination factor for the power-budget equations formulated in terms of average power [13,16]. F Π accounts for arbitrary electrical termination at signal reception.
Equations (1) and (2) are equivalent to the expressions given in the EK500 manual [11], except for the factor F Π , which was neglected there, implicitly implying F Π = 1 [13,16], cf. Section 4.1.
Dividing σ b s with a reference area (chosen equal to 1 m2), multiplying s v by a reference length (chosen equal to 1 m), and applying 10 log on both sides of Equations (1) and (2), leads to logarithmic (dB) versions of Equations (1) and (2), commonly denoted as target strength, TS, and volume backscattering strength, Sv, respectively [22].
Equations (1) and (2) constitute formulation A of the power-budget equations for σ b s and s v . Physical interpretations of Equations (1) and (2) in terms of power flow are given in [22].

2.3. Echo Integration Formulation (“Formulation B”)

2.3.1. Power Budget Equations

Calculation of the electrical powers Π T , Π R s t , and Π R v involved in the average power formulation of the power-budget equations, Equations (1) and (2), is normally made using time integration of voltage signals, denoted echo integration [3,4,6,7,8]. An echo integration formulation of the power-budget equations for σ b s and s v (denoted “formulation B”) is given in the following, derived from Equations (1) and (2).
From Equations (1), (2), and (8)–(10) one obtains
| V R s t ( f ) | 2 = σ b s ( f ) G 2 ( θ , ϕ , f ) λ 2 ( f ) F V V 2 ( f ) 16 π 2 r 4 e 4 α ( f ) r | V T ( f ) | 2
| V R v ( f ) | 2 = c 0 τ p s v ( f ) G 0 2 ( f ) ψ ( f ) λ 2 ( f ) F V V 2 ( f ) 32 π 2 r 2 e 4 α ( f ) r | V T ( f ) | 2 ,
where all quantities being a function of the frequency, f, have been indicated by the argument ( f ) , and
F V V 2 R T | Z E | | Z R + Z E | | Z T |
is the (non-dimensional) electrical termination factor for the power-budget equations formulated in terms of echo integration [16]. F V V accounts for arbitrary electrical termination at signal reception [16]. Integration over all frequencies gives
| V R s t ( f ) | 2 d f = 1 16 π 2 r 4 σ b s ( f ) G 2 ( θ , ϕ , f ) λ 2 ( f ) F V V 2 ( f ) e 4 α ( f ) r | V T ( f ) | 2 d f
| V R v ( f ) | 2 d f = c 0 τ p 32 π 2 r 2 s v ( f ) G 0 2 ( f ) ψ ( f ) λ 2 ( f ) F V V 2 ( f ) e 4 α ( f ) r | V T ( f ) | 2 d f
Now, assume that the transmitted signal, V T ( t ) , is a narrowband signal with carrier frequency f 0 (e.g., a tone burst of duration » 1 / f 0 ). The magnitude frequency spectrum | V T ( f ) | then has a peak value at f 0 , and a narrow main pass-band around f 0 , with relatively low side band level. Consequently, the dominating contributions to the integrands of Equations (15)–(16) may be assumed to come from the frequency region around f 0 . Under this assumption Equations (15)–(16) are approximately given as
| V R s t ( f ) | 2 d f σ b s ( f 0 ) G 2 ( θ , ϕ , f 0 ) λ 2 ( f 0 ) F V V 2 ( f 0 ) 16 π 2 r 4 e 4 α ( f 0 ) r | V T ( f ) | 2 d f
| V R v ( f ) | 2 d f c 0 τ p s v ( f 0 ) G 0 2 ( f 0 ) ψ ( f 0 ) λ 2 ( f 0 ) F V V 2 ( f 0 ) 32 π 2 r 2 e 4 α ( f 0 ) r | V T ( f ) | 2 d f
From Parseval’s theorem [31] the energy contained in a signal’s frequency spectrum is equal to the energy contained in the signal itself. Consequently, integration over frequency can be expressed as integration over time, giving
| V R s t ( t ) | 2 d t σ b s ( f 0 ) G 2 ( θ , ϕ , f 0 ) λ 2 ( f 0 ) F V V 2 ( f 0 ) 16 π 2 r 4 e 4 α ( f 0 ) r | V T ( t ) | 2 d t
| V R v ( t ) | 2 d t c 0 τ p s v ( f 0 ) G 0 2 ( f 0 ) ψ ( f 0 ) λ 2 ( f 0 ) F V V 2 ( f 0 ) 32 π 2 r 2 e 4 α ( f 0 ) r | V T ( t ) | 2 d t
In practice, the time integrations of Equations (19)-(20) are evaluated over finite time intervals, corresponding to the time gatings used for the transmitted and received signals.
As defined above, τ p is the time duration of the transmitted signal V T ( t ) , starting at t = 0. Let τ s t t s t 2 t s t 1 be the time duration of the voltage echo V R s t ( t ) that is received from a single scattering target, where t s t 1 and t s t 2 are the arrival times of the start and end of V R s t ( t ) . Similarly, let τ s p h t s p h 2 t s p h 1 be the time duration of the voltage echo V R s p h ( t ) that is received from the calibration sphere, where t s p h 1 and t s p h 2 are the arrival times of the start and end of V R s p h ( t ) . Let τ g = t g 2 t g 1 be the “gate opening time” used for the voltage waveform V R v ( t ) that is received from the “gated volume”, V g , where t g 1 and t g 2 are the times of gate opening and closure [16].
τ s p h is always larger than τ p , due to “ringing” (transient decay) caused by limited transducer bandwidth, possible electronic filtering, the finite dimensions of the sphere, possible “ringing” due to elastic sphere vibration modes and circumferential waves at the sphere surface, etc., in the sphere echo. Calibration spheres are normally made to produce strong echoes, and for typical calibration distances in the range 10–25 m, the signal-to-noise ratio (SNR) is very high. Sphere echoes are significantly stronger than fish echoes, and possible false echoes (from fish or other objects) may be excluded from signal processing. If the presence of fish becomes a problem, the calibration exercise is normally paused. A voltage detection threshold can then be used to electronically estimate t s p h 1 and t s p h 2 , determining τ s p h . In practice, a slightly different method may often be used. A time window is then set manually, starting just before t s p h 1 and ending just after t s p h 2 , sufficiently wide to include all calibration sphere echoes from different sphere positions over the transducer’s beam, in the vicinity of the acoustic axis. Which method to use is a question of uncertainty.
With respect to t s t , the time duration of single-target echoes (e.g., individual fish), a fish produces a weaker echo than a calibration sphere, so for the same measurement distance, the SNR is lower. The voltage detection threshold method can still be used in most cases, determining t s t 1 and t s t 2 to give τ s t . One has to ensure that possible false echoes (from other fish or objects) are excluded from processing of the signal.
τ g may be smaller, equal to, or larger than τ p , depending on the echo integration solutions implemented in the echo sounder system [1,5,16]. A choice τ g τ p may often be relevant in practice. The operator of the instrument typically sets the “channel depth” (equal to 1 2 c 0 τ g ) of the “gated volume” ( V g ). 10 m may be common in open sea (corresponding to τ g 133 ms), and 1–2 m close to the sea bottom (corresponding to τ g 13–27 ms). τ g depends on the measured sound velocity, c 0 . In general one thus has τ p τ s t τ s p h τ g .
Consequently, due to different time gatings of the different signals, one cannot compare the signal energies directly, since these depend on the integration times. However, the time averaged energies can be compared, i.e., the effective (rms) values of the time gated signals. Equations (19)-(20) may thus be approximated by
1 τ s t t s t 1 t s t 1 | V R s t ( t ) | 2 d t σ b s ( f 0 ) G 2 ( θ , ϕ , f 0 ) λ 2 ( f 0 ) F V V 2 ( f 0 ) 16 π 2 r 4 e 4 α ( f 0 ) r 1 τ p 0 τ p | V T ( t ) | 2 d t
1 τ g t g 1 t g 2 | V R v ( t ) | 2 d t c 0 τ p s v ( f 0 ) G 0 2 ( f 0 ) ψ ( f 0 ) λ 2 ( f 0 ) F V V 2 ( f 0 ) 32 π 2 r 2 e 4 α ( f 0 ) r 1 τ p 0 τ p | V T ( t ) | 2 d t ,
giving
σ b s = 16 π 2 r 4 e 4 α r τ p [ t i v s ] R s t G 2 ( θ , ϕ ) λ 2 τ s t F V V 2 [ t i v s ] T
s v = 32 π 2 r 2 e 4 α r [ t i v s ] R v G 0 2 ψ λ 2 c 0 τ g F V V 2 [ t i v s ] T ,
where σ b s , s v , α , λ , G 0 , G ( θ , ϕ ) , ψ , and F V V are all evaluated at the carrier frequency, f 0 , and
[ t i v s ] T 0 τ p | V T ( t ) | 2 d t
[ t i v s ] R s t t s t 1 t s t 2 | V R s t ( t ) | 2 d t
[ t i v s ] R s p h t s p h 1 t s p h 2 | V R s p h ( t ) | 2 d t
[ t i v s ] R v t g 1 t g 2 | V R v ( t ) | 2 d t
are defined as the “echo integral”—or “time-integral-voltage-squared” [tivs]—values of the transmitted voltage signal (the “pulse”), V T ( t ) , and three different received voltage waveforms, V R s t ( t ) , V R s p h ( t ) , and V R v ( t ) , for a single target echo (e.g., individual fish), a calibration sphere echo, and the echo received from a “gated volume” in oceanic surveying, respectively [16]. In this method, by using the effective (rms) values of the time gated signals, σ b s and s v become essentially independent of pulse duration. By using the same transmit pulse V T ( t ) in calibration and field operation, s v becomes essentially independent of pulse shape, under the assumptions (b) and (d) of small-amplitude sound propagation.
The definition of the echo integrals corresponds to that given by [3,4,6]. The [tivs] notation used here represents a generalization of the [tips] (“time-integral-pressure-squared”) notation used by Medwin and Clay [5], to account for voltage instead of “in-water” sound pressure signals [16]. The echo integrals ([tivs]) represent energy times electrical impedance, in units of V2s = W·ohm·s = J·ohm.

2.3.2. Formulation B Calibration Factor G 0

Now, consider calibration of the echo sounder or sonar system using a single metal sphere, with a theoretically known value for the backscattering cross section of the sphere, σ b s , t h e o r y s p h , calculated at the frequency f 0 in question [or the corresponding target strength T S t h e o r y s p h 10 log ( σ b s , t h e o r y s p h / r 1 2 ) , for a reference area r 1 2 = 1 m2]. Consider a calibration situation with the centre of the calibration sphere located at position ( r s p h , θ s p h , ϕ s p h ) . (Quantities and parameters that are measured or assumed known at calibration using the metal calibration sphere, are given sub- or superscript “sph”.)
Formulation B involves use of a single calibration fator, G 0 . Two approaches are used to derive G 0 expressions for use in calibration. The first approach (denoted “Method 1”) gives G 0 directly in terms of the basic quantities being measured in calibration. The second approach (denoted “Method 2”) is used to establish an alternative relationship for use in calibration that is similar to those used in [11] for the two calibration factors of the EK500 system, cf. Equations (47) and (50).
The basic quantities being measured in echo sounder calibration are r s p h , θ s p h , ϕ s p h , τ p , τ s p h , [ t i v s ] R s p h , [ t i v s ] T , | B i ( θ s p h , ϕ s p h ) | and c 0 , s p h . F V V , s p h and σ b s , t h e o r y s p h are typically assumed to be known. These quantities are all evaluated at the carrier frequency of the transmitted signal, f 0 .
Method 1:
From Equations (3), (4), and (23), G 0 can be determined from calibration measurements as
G 0 2 = 16 π 2 r s p h 4 e 4 α r s p h τ p [ t i v s ] R s p h λ s p h 2 F V V , s p h 2 τ s p h σ b s , t h e o r y s p h | B i ( θ s p h , ϕ s p h ) | 4 [ t i v s ] T .
Method 2:
An alternative and equivalent method to determine G 0 by echo sounder calibration using the metal sphere, may be derived as follows. From Equations (3), (4), and (23) one has
G 0 2 σ b s s p h = 16 π 2 r s p h 4 e 4 α r s p h τ p [ t i v s ] R s p h λ s p h 2 F V V , s p h 2 τ s p h | B i ( θ s p h , ϕ s p h ) | 4 [ t i v s ] T C 1 ,
where C 1 is a constant and known (measured) value for a given calibration measurement. Therefore, by knowing σ b s , t h e o r y s p h , the relevant calibration value G 0 is given by G 0 2 = C 1 / σ b s , t h e o r y s p h .
Consider a situation where an old (or an arbitrary, and possibly erroneous) calibration value is used initially in the calibration measurement, denoted G 0 , o l d . The corresponding value measured initially for the backscattering coefficient of the calibration sphere is then σ b s , o l d s p h = C 1 / G 0 , o l d 2 . It follows that C 1 = G 0 2 σ b s , t h e o r y s p h = G 0 , o l d 2 σ b s , o l d s p h , giving
G 0 2 = G 0 , o l d 2 σ b s , o l d s p h σ b s , t h e o r y s p h
10 log ( G 0 ) = 10 log ( G 0 , o l d ) + 1 2 [ T S o l d s p h T S t h e o r y s p h ] .
(Since from instrument manuals and relevant literature (cf. e.g., [11,18,19,20,21]) the calibration factors of modern scientific echo sounder and sonar systems are essentially given on logarithmic form, many users of such equipment know these factors in logarithmic form only. For completeness and convenience, both the linear and logarithmic expressions for the calibration factors are thus given here and in Section 3.)
The calibration value G 0 may thus be calculated from the old calibration value, G 0 , o l d . σ b s , o l d s p h corresponds to the target strength T S o l d s p h 10 log ( σ b s , o l d s p h / r 1 2 ) [22].
From these equations it follows that Equations (31) and (32) are equivalent to Equation (29). It thus follows that Methods 1 and 2 for determination of G 0 under Formulation B, are both generally valid approaches, leading to the same calibration value for G 0 .
Equations (23) and (24) for σ b s and s v , together with Equation (29) or Equation (31) for the calibration factor G 0 , constitute formulation B of the functional relationship for fish abundance estimation.

2.3.3. “Compact Formulation B” Functional Relationship

A compact version of Formulation B may be derived as follows. Insertion of Equation (29) into Equations (23) and (24) yields
σ b s = e 4 α ( r r s p h ) ( r r s p h ) 4 ( λ s p h λ ) 2 ( F V V , s p h F V V ) 2 τ s p h τ s t [ t i v s ] R s t [ t i v s ] R s p h | B i ( θ s p h , ϕ s p h ) B i ( θ , ϕ ) | 4 σ b s , t h e o r y s p h
s v = e 4 α ( r r s p h ) r 2 r s p h 4 ( λ s p h λ ) 2 ( F V V , s p h F V V ) 2 2 c 0 τ g τ s p h τ p [ t i v s ] R v [ t i v s ] R s p h | B i ( θ s p h , ϕ s p h ) | 4 ψ σ b s , t h e o r y s p h ,
respectively. Equations (33)-(34) represent an alternative and equivalent generic echo integration functional relationship for fish abundance estimation, given in terms of the basic quantities being measured during calibration and oceanic surveying [22].
Equations (33)-(34) reveal which parameters that influence the calibration and oceanic survey measurements. These expressions may be particularly useful for uncertainty or sensitivity studies, such as with respect to effects of sea temperature, or other possible errors or drift in the parameters involved in abundance measurement and species identification. The influence on the measurement accuracy from an uncertainty—or a possible error—in each of the quantities being measured or assumed known in calibration and oceanic survey operation can be investigated by adding an uncertainty term for this (or these) parameter(s), and using Equations (33)-(34) to calculate the resulting changes in σ b s and s v .

2.3.4. Fish Density Estimation

s v as given by Equation (24), or equivalently, Equation (34), represents volume backscattering from the thin spherical shell sub-volume V p of thickness d r p = 1 2 c 0 τ p in the observation volume V o b s . The volume backscattering from V o b s , between ranges r min and r max , is obtained by measuring s v for a continuous sequence of “ping volumes”, V p , and integrating s v over the range of these volumes, giving the area backscattering coefficient of V o b s [11,16,17]:
s a r min r max s v ( r ) d r ,
representing the backscattering cross section per unit area, within V o b s . Due to the integration of s v (in uints of m−1) over distance, s a is a non-dimensional quantity. In echo sounder output, s a is frequently given in terms of the nautical area scattering coefficient (NASC) [6,11,17]:
s A 4 π 1852 2 s a ,
where s a has been multiplied by the surface area of a sphere with radius one nautical mile. The density of targets (fish), expressed as the (non-dimensional) number of specimens in V o b s per square nautical mile, is then given as [17]
ρ a = s A 4 π σ b s ,
where σ b s is the expected value of the backscattering cross section (representing the expected TS) of individual targets (fish) [3,6,8,10]. Insertion of Equation (24) in Equations (35)–(37) yields [16]
ρ a = C c a l ψ σ b s E ,
where
E r min r max r 2 e 4 α r [ t i v s ] R v d r
C c a l 32 π 2 1852 2 G 0 2 λ 2 c 0 τ g F V V 2 [ t i v s ] T .
Equation (38) is in the form of the traditional echo-integrator equation [3,4,6,8,10,16], where C c a l represents the calibration factor used in that equation. E is a “range integrated echo integral” for the observation volume V o b s , obtained by integrating the sequence of echo integrals [ t i v s ] R v , each associated with a “ping volume” V p . For each V p at range r in V o b s , the term r 2 e 4 α r in E is the usual “20log(r) + 2 α ^ r ” TVG (time-varied gain) factor for volume backscattering from V p [3], where α ^ is the absorption coefficient expressed in dB/m. Similarly, for each V p , [ t i v s ] R v is the echo integral for the “gated volume” V g , taken over the time interval τ g , and typically averaged over many transmissions [3,4]. Since [ t i v s ] R v is calculated for each V p , it depends on range, r. In refs. [3,4,6,10] E is referred to as the “echo-integrator output”.
Equations (38)–(40) give the connection between the theory presented here (Formulations A–E), and the traditional echo-integrator equation used e.g., by [3,4,6,8,10]. In that literature, C c a l was used solely as an unspecified calibration factor. Here, the functional relationship for C c a l is derived, and fully given for small-amplitude conditions in terms of the echo sounder system parameters.

3. Application to Echo Sounder and Sonar Systems

For single-target measurements (such as in calibration, or for TS measurement of individual fish), some manufacturers base their σ b s (i.e., TS) measurement on the peak value of the echo samples in the single target echo, and not the echo signal waveform itself, as used for s v measurements in survey operation [11,18,19,20,21]. Cf. Equations (25)-(28). In the following the consequences of such a “peak voltage echo integration” approach in single-target TS measurements are investigated.
For this specific signal processing approach, the generic theory of Section 2 (Formulation B) is used to derive three alternative formulations of the power-budget equations with calibration factors, denoted “Formulation C”, “D”, and “E”, respectively. Since—for these three formulations—a specific signal processing method is employed for integrating the calibration sphere echo, Formulations C—E are less generic than Formulation B, but still equivalent to Formulation B, and fully valid for this approach.
The details of the mathematical derivations are given in [22], and omitted here. Only the resulting expressions for Formulations C–E are summarized, as a basis for the discussion of Section 4. In a few cases equations are duplicated in Section 3. This is done for completeness and convenience in summarizing each of Formulations C–E.

3.1. Echo Integration “Formulation C”

3.1.1. Power-Budget Equations

Using the approach described above with “peak voltage echo integration” in single-target TS measurements, σ b s and s v given by Equations (23) and (24) can equivalently be expressed as [22]:
σ b s = 16 π 2 r 4 e 4 α r [ t i v s ] R , max s t G e f f 2 ( θ , ϕ ) λ 2 F V V 2 [ t i v s ] T
s v = 32 π 2 r 2 e 4 α r [ t i v s ] R v G 0 2 ψ λ 2 c 0 τ g F V V 2 [ t i v s ] T ,
where
[ t i v s ] R , max s t | V R s t | max 2 τ p
G e f f 2 ( θ , ϕ ) G 2 ( θ , ϕ ) τ s t | V R s t | max 2 [ t i v s ] R s t .
Here [22], [ t i v s ] R , max s t is a simplified (and erroneous) estimate of the single target echo integral [ t i v s ] R s t , obtained by multiplying the transmit signal duration τ p with the constant squared voltage value, | V R s t | max 2 . | V R s t | max is the maximum (“peak voltage”) value of the magnitude of the single target echo, | V R s t ( t ) | , within the time window [ t s t 1 , t s t 2 ]. G e f f 2 ( θ , ϕ ) represents an “effective” value of G 2 ( θ , ϕ ) , where a correction factor has been included to compensate for the error introduced by using [ t i v s ] R , max s t instead of [ t i v s ] R s t , to ensure that the resulting σ b s measurement is valid and correct [22].
This approach represents a “rectangular echo integration” in single-target TS measurements, i.e., for σ b s . For s v , full waveform echo integration is used, as in Formulation B, and Equation (42) is identical to Equation (24).

3.1.2. Formulation C Calibration Factors G 0 , e f f and G 0

Formulation C involves two calibration factors, G 0 , e f f and G 0 , where G 0 , e f f G e f f ( 0 , 0 ) . Two approaches are used to derive expressions for these. The first approach (denoted “Method 1”) gives G 0 , e f f and G 0 in terms of the basic quantities being measured in calibration. The second approach (denoted “Method 2”) is used to (i) establish alternative relationships for use in calibration; and (ii) establish relationships for the two Formulation D calibration factors, G 0 , e f f and G 0 , that can be used for comparison with the two EK500 calibration factors “TS transducer gain” and “Sv transducer gain” (cf. Section 4.3).
The calibration factor G 0 , e f f can be expressed by two equivalent expressions [22]:
Method 1:
G 0 , e f f 2 16 π 2 r s p h 4 e 4 α r s p h [ t i v s ] R , max s p h λ s p h 2 σ b s , t h e o r y s p h | B i ( θ s p h , ϕ s p h ) | 4 F V V , s p h 2 [ t i v s ] T .
Method 2:
G 0 , e f f 2 = G 0 , e f f , o l d 2 σ b s , o l d s p h σ b s , t h e o r y s p h
10 log ( G 0 , e f f ) = 10 log ( G 0 , e f f , o l d ) + 1 2 [ T S e f f , o l d s p h T S t h e o r y s p h ] .
Here, G 0 , e f f , o l d is an old (or arbitrary, and possibly erroneous) calibration value used for G 0 , e f f initially in the calibration measurement. σ b s , o l d s p h is the corresponding old measured value used initially for the backscattering cross section of the calibration sphere, corresponding to target strength T S e f f , o l d s p h 10 log ( σ b s , o l d s p h / r 1 2 ) [22].
Similarly, the calibration factor G 0 can be expressed by two equivalent expressions [22]:
Method 1:
G 0 2 = 16 π 2 r s p h 4 e 4 α r s p h τ p [ t i v s ] R s p h λ s p h 2 F V V , s p h 2 τ s p h σ b s , t h e o r y s p h | B i ( θ s p h , ϕ s p h ) | 4 [ t i v s ] T .
Method 2:
G 0 2 = G 0 , o l d 2 s a , o l d s p h s a , t h e o r y s p h
10 log ( G 0 ) = 10 log ( G 0 , o l d ) + 1 2 10 log ( s a , o l d s p h s a , t h e o r y s p h ) ,
where
s a , t h e o r y s p h = τ s p h τ p | B i ( θ s p h , ϕ s p h ) | 4 ψ r s p h 2 σ b s , t h e o r y s p h .
Here, G 0 , o l d is an old (or arbitrary, and possibly erroneous) calibration value used for G 0 initially in the calibration measurement. s a , o l d s p h is the corresponding old measured value used initially for the area backscattering coefficient of the sphere [22]. s a , t h e o r y s p h is a known value of the area backscattering coefficient of the sphere, given from calibration measurements and the known (calculated) value σ b s , t h e o r y s p h . (A similar expression as Equation (51) is given in [11], but without | B i ( θ s p h , ϕ s p h ) | and the factor τ s p h / τ p . That expression thus relates to the special case of calibration sphere located on the acoustic axis, θ s p h = ϕ s p h = 0 , and τ s p h = τ p (where the latter condition may not be valid in practice).)
Equations (41) and (42) for σ b s and s v , together with Equations (45)–(51) for the calibration factors G 0 , e f f and G 0 , constitute Formulation C of the functional relationship for fish abundance estimation.

3.2. Echo Integration “Formulation D”

3.2.1. Power-Budget Equations

By defining s a , c o r r G 0 / G 0 , e f f , σ b s and s v given by Equations (41) and (42) can equivalently be expressed as [22]:
σ b s = 16 π 2 r 4 e 4 α r [ t i v s ] R , max s t G e f f 2 ( θ , ϕ ) λ 2 F V V 2 [ t i v s ] T
s v = 32 π 2 r 2 e 4 α r [ t i v s ] R v G 0 , e f f 2 ψ λ 2 c 0 τ g s a , c o r r 2 F V V 2 [ t i v s ] T ,
where Equation (52) is identical to Equation (41), and [22]:
s a , c o r r 2 G 0 G 0 , e f f = [ t i v s ] R s t τ s t | V R s t | max 2 = τ p τ s t [ t i v s ] R s t [ t i v s ] R , max s t .
From Equations (9) and (26) it follows that
Π R s t = R E | Z E | 2 ( V R s t , r m s ) 2 = R E | Z E | 2 1 τ s t t s t 1 t s t 2 | V R s t ( t ) | 2 d t = R E | Z E | 2 [ t i v s ] R s t τ s t
Π R , m a x s t = R E | Z E | 2 | V R s t | max 2 ,
giving, from Equations (43) and (54),
s a , c o r r 2 = Π R s t Π R , max s t .
Here, V R s t , r m s is the effective (rms) amplitude of the echo received from a single target (e.g., a fish, or a calibration sphere), and Π R , m a x s t is the “maximum average electrical power” of that echo. Equation (57) shows that s a , c o r r 2 represents the (non-dimensional) ratio of the actual average electrical power, Π R s t , in a single-target echo (i.e., from a calibration sphere, or a fish), to the “maximum average electrical power” of that echo, Π R , m a x s t (i.e., the average power obtained by using “peak voltage echo integration”, instead of full waveform echo integration).

3.2.2. Formulation D Calibration Factors G 0 , e f f and s a , c o r r

Formulation D involves two calibration factors, G 0 , e f f and s a , c o r r . Similarly to Formulation C, two approaches are used to derive expressions for these. The first approach (denoted “Method 1”) gives G 0 , e f f and s a , c o r r in terms of the basic quantities being measured in calibration. The second approach (denoted “Method 2”) is used to (i) establish alternative relationships for use in calibration; and (ii) establish relationships for the two formulation D calibration factors, G 0 , e f f and s a , c o r r , that can be used for comparison with the two EK60 calibration factors “Gain” and “Sa correction” (cf. Section 4.4). For s a , c o r r , a third approach (denoted “Method 3”) is also used, which may be of interest if an EK60 type echo sounder (or similar) is operated as an EK500 type echo sounder (i.e., using calibration factors G 0 , e f f and G 0 ).
The calibration factor G 0 , e f f can be expressed by two equivalent expressions, given by Equations (45) and (46) [22].
The calibration factor s a , c o r r can be expressed by three equivalent expressions [22]:
Method 1:
s a , c o r r 2 = 16 π 2 r s p h 4 e 4 α r s p h τ p [ t i v s ] R s p h G 0 , e f f 2 λ s p h 2 τ s p h σ b s , t h e o r y s p h | B i ( θ s p h , ϕ s p h ) | 4 F V V , s p h 2 [ t i v s ] T .
Method 2:
s a , c o r r 2 = s a , c o r r , o l d 2 s a , o l d s p h s a , t h e o r y s p h
S a , c o r r = S a , c o r r , o l d + 1 2 10 log ( s a , o l d s p h s a , t h e o r y s p h ) .
Method 3:
s a , c o r r = G 0 G 0 , e f f
S a , c o r r = 10 log ( G 0 ) 10 log ( G 0 , e f f ) .
Equations (52) and (53) for σ b s and s v , together with Equations (45)–(47) and (58)–(62), respectively, for the calibration factors G 0 , e f f and s a , c o r r , constitute Formulation D of the functional relationship for fish abundance estimation.

3.3. Echo Integration “Formulation E”

3.3.1. Power-Budget Equations

By replacing τ g s a , c o r r 2 with τ e f f in the s v equation, σ b s and s v given by Equations (52) and (53) can equivalently be expressed as [22]:
σ b s = 16 π 2 r 4 e 4 α r [ t i v s ] R , max s t G e f f 2 ( θ , ϕ ) λ 2 F V V 2 [ t i v s ] T
s v = 32 π 2 r 2 e 4 α r [ t i v s ] R v G 0 , e f f 2 ψ λ 2 c 0 τ e f f F V V 2 [ t i v s ] T ,
where Equation (63) is identical to Equations (41) and (52). From Equations (54) and (57) it follows that [22]
τ e f f τ g s a , c o r r 2 = τ g τ p τ s t [ t i v s ] R s t [ t i v s ] R , max s t = τ g Π R s t Π R , max s t .
The latter expression in Equation (65) shows that τ e f f represents a scaled gate opening time, where τ g (used for volume backscattering in oceanic surveying) is scaled by the ratio of the actual average electrical power, Π R s t , to the “maximum average electrical power”, Π R , m a x s t , of the single target echo. τ e f f may thus be interpreted as an “effective echo integration time” [22].
From Equations (26), (27), (43), and (65) it also follows that
[ t i v s ] R s t t s t 1 t s t 2 | V R s t ( t ) | 2 d t = τ e f f τ s t τ g | V R s t | max 2
[ t i v s ] R s p h t s p h 1 t s p h 2 | V R s p h ( t ) | 2 d t = τ e f f τ s p h τ g | V R s p h | max 2
for the echo integrals of a single target and a calibration sphere, respectively. Consequently, for single targets, full waveform echo integration can be carried out using | V R s t | max and τ e f f , provided τ s t and τ g are known.

3.3.2. Formulation E Calibration Factors G 0 , e f f and τ e f f

Formulation E involves two calibration factors, G 0 , e f f and τ e f f . Similarly to Formulations C and D, two approaches are used to derive expressions for these. The first approach (denoted “Method 1”) gives G 0 , e f f and τ e f f in terms of the basic quantities being measured in calibration. The second approach (denoted “Method 2”) is used to establish alternative relationships for use in calibration. For τ e f f , a third approach (denoted “Method 3”) is also given, which may be of interest if an EK60 type echo sounder (or similar) is operated as an EK500 type of echo sounder (i.e., using calibration factors G 0 , e f f and G 0 ).
As for Formulations C and D, the calibration factor G 0 , e f f can be expressed by two equivalent expressions, given by Equations (45) and (46) [22].
The calibration factor τ e f f can be expressed by three equivalent expressions [22]:
Method 1:
τ e f f = 16 π 2 r s p h 4 e 4 α r s p h τ p [ t i v s ] R s p h G 0 , e f f 2 λ s p h 2 σ b s , t h e o r y s p h | B i ( θ s p h , ϕ s p h ) | 4 F V V , s p h 2 [ t i v s ] T .
Method 2:
τ e f f = τ e f f , o l d s a , o l d s p h s a , t h e o r y s p h
10 log ( τ e f f τ r e f ) = 10 log ( τ e f f , o l d τ r e f ) + 10 log ( s a , o l d s p h s a , t h e o r y s p h ) .
Method 3:
τ e f f = τ g G 0 2 G 0 , e f f 2
10 log ( τ e f f τ r e f ) = 10 log ( τ g τ r e f ) + 20 log ( G 0 G 0 , e f f ) .
Here, τ r e f is a reference time interval for τ e f f (chosen equal to 1 s). τ e f f , o l d is an old calibration (or arbitrary, and possibly erroneous) value, used initially for τ e f f in the calibration procedure.
Equations (63) and (64) for σ b s and s v , together with Equations (45)–(47) and (68)–(72) for the calibration factors G 0 , e f f and τ e f f , constitute Formulation E of the functional relationship for fish abundance estimation.

4. Results and Discussion

Four equivalent echo integration Formulations B–E of the power-budget equations for σ b s and s v 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, F Π .
The expressions given in [11] correspond to setting F Π = 1 in Formulation A, which is valid for electrical termination conditions for which Z E = Z R * (conjugate matched electrical termination), or for Z E = Z R when X T = 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 σ b s and s v 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 σ b s and s v , together with Equations (29) or (31) for the calibration factor G 0 , 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, F V V . It also involves several integration time intervals, τ p , τ s t , τ s p h , and τ g , 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, G 0 . 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 σ b s and s v , together with Equations (45)–(51) for the two calibration factors G 0 , e f f and G 0 , 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 σ b s and s v 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 F Π = 1 ), whereas the EK500 implementation involves full waveform echo integration for s v , “peak voltage echo integration” for single-target TS measurements ( σ b s ), and two calibration factors “TS transducer gain” and “ S v 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 s v , and “peak voltage echo integration” for single-target TS measurements ( σ b s ) in Formulation C, expressions are here obtained that involve the calibration factors G 0 , e f f and G 0 . By defining TS transducer gain 10 log 10 ( G 0 , e f f ) and S v transducer gain 10 log 10 ( G 0 ) , 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 σ b s and s v given in [11], the analysis shows that the expressions for σ b s and s v 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 σ b s and s v 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 G 0 , e f f has been introduced here to account for “peak voltage echo integration” in single-target TS measurements. No such parameter as G 0 , e f f 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 σ b s and s v , together with Equations (45)–(47) and (58)–(62), respectively, for the two calibration factors G 0 , e f f and s a , c o r r 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 s v as given by [11] and [15], respectively, is not equal to 1, but s a , c o r r 2 . This indicates that the expressions for σ b s and s v 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 s a , c o r r in a consistent way. In Formulations C and D, the expressions for σ b s and s v are equivalent, and the ratio of the two expressions for s v equals 1.
Secondly, in Formulation D, expressions are obtained that involve the calibration factors G 0 , e f f and s a , c o r r . By defining Gain 10 log 10 ( G 0 , e f f ) and S a correction S a , c o r r 10 log 10 ( s a , c o r r ) , 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 “ S a correction” calibration factors in the manuals of EK60, etc. [18,19,20,21].
The expressions postulated by Ona et al. [15] for σ b s and s v correspond to Formulation D. There are, however, some deviations. Firstly (in terms of the terminology used here), in [15], electrical powers Π T and Π R (in units of W) were used instead of echo integrals, [ t i v s ] T , [ t i v s ] R v , and [ t i v s ] R , max s t (in units of V2s = W·ohm·s = J·ohm). Echo integration was thus not accounted for in these expressions. Secondly, the same Π R was used for σ b s and s v , which is not consistent with use of “peak voltage echo integration” and full waveform echo integration, respectively, for TS and s v (and thus s A ) measurements. Thirdly, G ( θ , ϕ ) and G 0 were used instead of G e f f ( θ , ϕ ) and G 0 , e f f in the expressions for σ b s and s v , respectively, which is inconsistent when “peak voltage echo integration” is used for σ b s , and s a , c o r r is used in s v . Finally, in [15], electrical termination was not accounted for, implicitly implying F Π = 1 [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 G ( θ , ϕ ) and G 0 , referred to in [15] as “transducer gain” and “on-axis transducer gain”, respectively, are not equal to the transducer gain and axial transducer gain, G ( θ , ϕ ) 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 a , c o r r parameter; and (c) “harmonize the TS and s A measurements” [19] by employing the same calibration factor G0,eff in σbs and s v . 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 G 0 , e f f and s a , c o r r , given by Equations (45)–(47) and (58)–(62), respectively.
Through the present analysis, the relationship between the Formulation D calibration factors G 0 , e f f and s a , c o r r , and the Formulation C calibration factors, G 0 , e f f and G 0 , 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 σ b s and s v derived in [23] correspond exactly to those postulated in [15], the above discussion also applies to [23].

4.5. On the G e f f , s a , c o r r , and τ e f f Parameters

The reasons for introducing the parameter G e f f , s a , c o r r , and τ e f f 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 ( σ b s ) and in survey operation ( s v ), only a single calibration factor, G 0 , applies, with no need for G e f f , s a , c o r r , and τ e f f 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 σ b s , which is compensated for by introducing G 0 , e f f instead of G 0 in σ b s , cf. Equations (41), (43), and (44). By still using G 0 in the expression for s v , cf. Equation (42), two calibration factors G 0 , e f f and G 0 thus become involved in Formulation C, referred to as “TS transducer gain” and “Sv transducer gain”, respectively.
In Formulation D, the calibration factor G 0 used in Formulation C (in the s v expression) is replaced by G 0 , e f f , to “harmonize” the σ b s and s v equations. An error is then introduced in s v , which is compensated for by using s a , c o r r in s v , referred to as the “Sa correction”, cf. Equation (53).
In Formulation E, the term τ g s a , c o r r 2 used in Formulation D (in the s v equation) is denoted τ e f f . As a result, the σ b s and s v 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 G ( θ , ϕ ) , G 0 , and τ g . In Formulation E these are G ( θ , ϕ ) , G 0 , e f f , and τ e f f , due to the “peak voltage echo integration” used for single-target TS measurements.
In [15] and [23], τ e f f is referred to as an “effective pulse duration”. This interpretation may be discussed. In the present analysis, s a , c o r r , used to define τ e f f 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 [ t i v s ] R , max s t instead of [ t i v s ] R s t for these targets. Hence, from Equation (65) it appears that τ e f f 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 τ e f f 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 σ b s and s v known from prior literature [13,16] (denoted Formulation A), four new echo integration formulations for σ b s and s v 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 G e f f , s a , c o r r , and τ e f f type of parameters, and only a single calibration factor G 0 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 G e f f , s a , c o r r , or τ e f f 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 σ b s and s v 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 σ b s and s v , respectively. In terms of the terminology used here, the average electrical powers Π T and Π R used by [15,23] are to be replaced by echo integrals [ t i v s ] T , [ t i v s ] R v , and [ t i v s ] R , max s t . Similarly, the transducer gains G ( θ , ϕ ) and G 0 used by [15,23] are to be replaced by G e f f ( θ , ϕ ) and G 0 , e f f . With these modifications applied to the expressions given in [15,23], the expressions given by Formulations B–D for σ b s and s v 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 s a , c o r r and τ e f f , which were postulated by [15] and adopted by [23]. s a , c o r r and τ e f f 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 σ b s , s v , 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 σ b s and s v 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 σ b s and s v 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.

Acknowledgments

The authors wish to acknowledge Geir Pedersen, Christian Michelsen Research AS, Bergen, Norway; Egil Ona and Gavin J. Macaulay, Institute of Marine Research, Bergen, Norway; and Frank Tichy, Kongsberg Maritime, Horten, Norway, for useful discussions related to topics discussed here.

Author Contributions

The theory (Formulations B–E) was developed by Per Lunde, including description of the relationships between the Formulations A–E, prior literature, and expressions given for scientific echo sounders and sonar systems. Rolf J. Korneliussen contributed with long-time experience from fisheries research surveys and abundance estimation using e.g., the Simrad EK500, EK60, and MS70 systems. Both authors contributed to the assembly and writing of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Clay, C.S.; Medwin, H. Acoustical Oceanography: Principles and Applications; John Wiley & Sons: New York, NY, USA, 1977; pp. 180–182, 205, 220, 225–235, 395–404. [Google Scholar]
  2. Nakken, O.; Ulltang, Ø. A Comparison of the Reliability of Acoustic Estimates for Fish Stock Abundance and Estimates Obtained by Other Assessment Methods in Northeast Atlantic; FAO Fish. Rep. No. 300; Food and Agriculture Organization of the United Nations: Rome, Italy, 1983; pp. 249–260. [Google Scholar]
  3. MacLennan, D.N. Acoustical measurements of fish abundance. J. Acoust. Soc. Am. 1990, 87(1), 1–15. [Google Scholar] [CrossRef]
  4. MacLennan, D.N.; Simmonds, E.J. Fisheries Acoustics; Chapman & Hall: London, UK, 1992; pp. 68–69, 147–151. [Google Scholar]
  5. Medwin, H.; Clay, C.S. Fundamentals of Acoustical Oceanography; Academic Press: Boston, MA, USA, 1998; pp. 153–158, 350–361. [Google Scholar]
  6. Simmonds, J.; MacLennan, D.N. Fisheries Acoustics. Theory and Practice, 2nd ed.; Blackwell Science Ltd.: Oxford, UK, 2005; pp. 32, 47, 59–60, 110, 119–121, 187–191, 197. [Google Scholar]
  7. Dragesund, O.; Olsen, S. On the possibility of estimating year-class strength by measuring echo-abundance of 0-group fish. Fiskeridirektoratets Skrifter, Serie Havundersøkelser 1965, 13(8), 48–75. [Google Scholar]
  8. Dalen, J.; Nakken, O. On the Application of the Echo Integration Method; ICES CM Document No. 1983/B:19; International Council for the Exploration of the Sea: Copenhagen, Denmark, 1983. [Google Scholar]
  9. Foote, K.G. Optimizing copper spheres for precision calibration of hydro acoustic equipment. J. Acoust. Soc. Am. 1982, 71, 742–747. [Google Scholar] [CrossRef]
  10. Foote, K.G.; Knudsen, H.P.; Vestnes, G.; MacLennan, D.N.; Simmonds, E.J. Calibration of Acoustic Instruments for Fish Density Estimation: A Practical Guide; ICES Cooperative Research Report No. 144; International Council for the Exploration of the Sea: Copenhagen, Denmark, 1987. [Google Scholar]
  11. Operator Manual: SIMRAD EK500 Fishery Research Echo Sunder. Scientific Echo Sounder: Base Version. Doc. no. P2170 / Rev. G, Simrad AS (now Kongsberg Maritime AS), Horten, Norway. 1997, pp. 165–169, 204, 211–212, 216–217, 223. Available online: https://bora.uib.no/handle/1956/2158?mode=full&submit_simple=Show+full+item+record (accessed on 24 February 2016).
  12. Korneliussen, R.J. Analysis and Presentation of Multi-Frequency Echograms. Ph.D. Thesis, University of Bergen, Bergen, Norway, 2002. [Google Scholar]
  13. Pedersen, A. Effects of Nonlinear Sound Propagation in Fisheries Research. Ph.D. Thesis, University of Bergen, Bergen, Norway, 2007. Available online: https://www.simrad.com/www/01/NOKBG0397.nsf/AllWeb/9F0DB3E3763D8CD3C125718B00332E5A/$file/043762ag_ek500_operator_manual_english_complete.pdf?OpenElement (accessed on 24 February 2016). [Google Scholar]
  14. Demer, D.A.; Renfree, J.S. Variations in echosounder—transducer performance with water temperature. ICES J. Mar. Sci. 2008, 65, 1021–1035. [Google Scholar] [CrossRef]
  15. Ona, E.; Mazauric, V.; Andersen, L.N. Calibration methods for two scientific multibeam systems. ICES J. Mar. Sci. 2009, 66, 1326–1334. [Google Scholar] [CrossRef]
  16. Lunde, P.; Pedersen, A.O.; Korneliussen, R.J.; Tichy, F.E.; Nes, H. Power Budget and Echo-Integrator Equations for Fish Abundance Estimation. Fisken og Havet no. 10/2013, Institute of Marine Research, Bergen, Norway, 2013. Available online: http://www.imr.no/publikasjoner/andre_publikasjoner/fisken_og_havet/nb-no (accessed on 24 February 2016).
  17. MacLennan, D.N.; Fernandes, P.G.; Dalen, J. A consistent approach to definitions and symbols in fisheries acoustics. ICES J. Mar. Sci. 2002, 59, 365–369. [Google Scholar] [CrossRef]
  18. Simrad EK60. Scientific Echo Sounder. Reference Manual. Release 2.4.X. Doc. no. 164692 / Rev. D, Kongsberg Maritime AS, Horten, Norway, 2012. pp. 16, 48–49, 191, 199, 221, (in electronic document). Available online: https://www.simrad.com/www/01/nokbg0397.nsf/AllWeb/62D6EBE0D8EEB97CC125718E004B41C7?OpenDocument (accessed on 10 July 2016).
  19. Simrad ES60. Fish Finding Echo Sounder. Reference Manual. Doc. no. 857-160970 / Rev.H, Kongsberg Maritime AS, Horten, Norway, 2006, p. 80. Available online: http://www.google.no/url?sa=t&rct=j&q=&esrc=s&frm=1&source=web&cd=1&ved=0CCwQFjAA&url=http%3A%2F%2Fwww.simrad.com%2Fwww%2F01%2FNOKBG0397.nsf%2FAllWeb%2F476DD263FD01A618C12570D100630F2B%2F%24file%2F160970ah_es60_reference_manual_english.pdf%3FOpenElement&ei=W4-EU9m1D4SB4gSn8IHoAw&usg=AFQjCNH8YxK7LuYNCspKNZ5g9bSlu8PgJA&sig2=V4kxvlh0jDwMNNDKxQ7pGQ&bvm=bv.67720277,d.bGE (accessed on 24 February 2016).
  20. Simrad ME70. Scientific Multibeam Echo Sounder. Reference Manual. Release 1.2.5. Doc. no. 328588 / Rev. B, Kongsberg Maritime AS, Horten, Norway, 2012, pp. 140, 143, 152, 154, 158, 210, 246, 448 (in electronic document). Available online: http://www.simrad.com/www/01/nokbg0397.nsf/AllWeb/DCB276EA00D5BD46C125719A003AAD45?OpenDocument (accessed on 24 February 2016).
  21. Simrad MS70. Scientific Multibeam Sonar. Operator MANUAL. Doc. no. 301088 / Rev. C, Kongsberg Maritime AS, Horten, Norway, 2007, pp. 164, 171, 179. Available online: http://www.google.no/url?sa=t&rct=j&q=&esrc=s&frm=1&source=web&cd=1&ved=0CC8QFjAA&url=http%3A%2F%2Fwww.imr.no%2Ffilarkiv%2F2006%2F08%2F301088ac_ms70_operator_manual_english.pdf%2Fnn-no&ei=65OEU_OeLoPm4QTB9IHYBQ&usg=AFQjCNEMV_2RNtMRZSqEXuo1QAuZkj3h5A&sig2=WRfF6vgyIjMjYk7yUVthkA&bvm=bv.67720277,d.bGE (accessed on 24 February 2016).
  22. Lunde, P.; Korneliussen, R.J. A Unifying Theory Explaining Different Power-Budget Formulations Used in Modern Scientific Echo Sounders for Fish Abundance Estimation. Fisken og Havet no. 7/2014, Institute of Marine Research, Bergen, Norway, 2014. Available online: http://www.imr.no/publikasjoner/andre_publikasjoner/fisken_og_havet/nb-no (accessed on 24 February 2016).
  23. Calibration of Acoustic Instruments. ICES Cooperative Research Report No. 326, International Council for Exploration of the Sea: Copenhagen, Denmark, May 2015, pp. 12–16, 62. Available online: http://www.google.no/url? sa=t&rct=j&q=&esrc=s&source=web&cd=1&ved=0ahUKEwjdyLWZrJDLAhWh93IKHVU9CpoQFgggMAA&url=http%3A%2F%2Fices.dk%2Fsites%2Fpub%2FPublication%2520Reports%2FCooperative%2520Research%2520Report%2520(CRR)%2Fcrr326%2FCRR326.pdf&usg=AFQjCNHpagCdqvZT8zLLQCVcKB8Jej-XFQ&bvm=bv.114733917,d.bGQ (accessed on 24 February 2016).
  24. Korneliussen, R.J.; Diner, N.; Ona, E.; Berger, L.; Fernandes, P.G. Proposals for the collection of multifrequency acoustic data. ICES J. Mar. Sci. 2008, 65, 982–994. [Google Scholar] [CrossRef]
  25. Tichy, F.E.; Solli, H.; Klaveness, H. Nonlinear effects in a 200-kHz sound beam and consequences for target strength measurement. ICES J. Mar. Sci. 2003, 60, 571–574. [Google Scholar] [CrossRef]
  26. Lunde, P.; Pedersen, A.O. Sonar and Power-Budget Equations for Backscattering of Finite Amplitude Sound Waves, with Implications in Fishery Acoustics for Abundance Estimation of Marine Resources. In Proceeding of 35th Scand. Symp. Phys. Acoust., Geilo, Norway, 29 January–1 Feburary 2012; Available online: http://www.ntnu.edu/sspa/program (accessed on 24 February 2016).
  27. Foote, K.G. Rather-high-frequency sound scattering by swimbladdered fish. J. Acoust. Soc. Am. 1985, 78, 688–700. [Google Scholar] [CrossRef]
  28. Foote, K.G. Linearity of fisheries acoustics, with addition theorems. J. Acoust. Soc. Am. 1983, 73(6), 1932–1940. [Google Scholar] [CrossRef]
  29. Balanis, C.A. Antenna Theory: Analysis and Design, 3rd ed.; John Wiley & Sons: Hoboken, NJ, USA, 2005; pp. 50, 66, 89–95. [Google Scholar]
  30. Kinsler, L.E.; Frey, A.R.; Coppens, A.B.; Sanders, J.V. Fundamentals of Acoustics, 4th ed.; John Wiley & Sons: New York, NY, USA, 2000; pp. 15, 181–184,189. [Google Scholar]
  31. Brigham, E.O. The Fast Fourier Transform; Prentice-Hall: Englewood Cliffs, NJ, USA, 1974; p. 64. [Google Scholar]
  32. Ona, E.; (Institute of Marine Research, Bergen, Norway). Personal communication, 2014.
Table 1. Alternative and equivalent formulations of power budget equations for measurement of σ b s and s v . Calibration factors referred to in the table are defined in Section 2 and Section 3.
Table 1. Alternative and equivalent formulations of power budget equations for measurement of σ b s and s v . Calibration factors referred to in the table are defined in Section 2 and Section 3.
Form No.Type of FormulationGeneric or NotCalibr. FactorsKey CharacteristicsRefs.
AAverage powerInstrument independent (generic)None• Not directly suitable for signal processing (echo integration) in this formulation[11]
[13]
[16]
BEcho integration------ “ ------ G 0 • Full waveform echo integration in σ b s and s v measurements[22]
C------ “ ------Instrument dependent G 0 , e f f
G 0
• “Peak voltage echo integration” in σ b s measurement
• Full waveform echo integration in s v measurement
• Equations related to Simrad EK500 type of echo sounder system [11]
[22]
D------ “ ----------- “ ------ G 0 , e f f
s a , c o r r
• “Peak voltage echo integration” in σ b s measurement
• Full waveform echo integration in s v measurement
• Equations related to Simrad EK60, ES60, ME70, and MS70 type of scientific echo sounder and sonar systems [15,18,19,20,21]
[22]
E------ “ ------------ “ ------ G 0 , e f f
τ e f f
• “Peak voltage echo integration” in σ b s measurement
• Full waveform echo integration in s v measurement
[22]

Share and Cite

MDPI and ACS Style

Lunde, P.; Korneliussen, R.J. Power-Budget Equations and Calibration Factors for Fish Abundance Estimation Using Scientific Echo Sounder and Sonar Systems. J. Mar. Sci. Eng. 2016, 4, 43. https://doi.org/10.3390/jmse4030043

AMA Style

Lunde P, Korneliussen RJ. Power-Budget Equations and Calibration Factors for Fish Abundance Estimation Using Scientific Echo Sounder and Sonar Systems. Journal of Marine Science and Engineering. 2016; 4(3):43. https://doi.org/10.3390/jmse4030043

Chicago/Turabian Style

Lunde, Per, and Rolf J. Korneliussen. 2016. "Power-Budget Equations and Calibration Factors for Fish Abundance Estimation Using Scientific Echo Sounder and Sonar Systems" Journal of Marine Science and Engineering 4, no. 3: 43. https://doi.org/10.3390/jmse4030043

APA Style

Lunde, P., & Korneliussen, R. J. (2016). Power-Budget Equations and Calibration Factors for Fish Abundance Estimation Using Scientific Echo Sounder and Sonar Systems. Journal of Marine Science and Engineering, 4(3), 43. https://doi.org/10.3390/jmse4030043

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop