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

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.


Introduction
Acoustic methods are widely used for estimating fish abundance , 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 obs [6,11,16], to give the fish density in V obs , ρ a [17], in terms of an echo-integrator equation [3,6,8,10,16].
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 σ bs 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 σ bs and s v was summarized by Korneliussen [12].
Demer and Renfree [14] also used expressions for σ bs 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 xσ bs y (the expected value of σ bs ) 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 σ bs 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.σ bs , 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 σ bs 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 σ bs 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 σ bs and s v given in [11], are available from the literature.That is, no power-budget equations for σ bs 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 σ bs 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 σ bs 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,corr " 10logps a,corr q 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 ´2 a,corr .It follows [22] that the power-budget equations as given in [15] cannot be derived directly from the power-budget equations given in [11].
Fifthly, the derivation given in [23] of the power-budget equations for σ bs and s v being postulated in [15], appears to be incomplete on essential points, such as with respect to the s a,corr parameter.In [23], s a,corr 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 2 a,corr .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 σ bs , 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 σ bs 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, σ bs 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,corr 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.
‚ "Peak voltage echo integration" in σ bs measurement ‚ Full waveform echo integration in s v measurement [22] 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 σ bs 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].

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) same transmit electrical power Π T is used for σ bs 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 pr, θ, φq, 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.

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 pr, θ, φq in the transducer's far field, σ bs , 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]: respectively.σ bs and s v are given in units of m 2 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 " pr p1 `rp2 q{2 is the mid-range of the "ping volume" V p , that is contained within ranges r p1 and r p2 , with thickness dr p " r p2 ´rp1 " 1 2 c 0 τ p [16].In practice, τ p = 1 ms may often be used, corresponding to dr 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 pθ, φq and G 0 (both non-dimensional) are the transducer gain and the axial transducer gain, defined as [11,13,16,29]: respectively, where η is the transducer's electroacoustic conversion efficiency (the non-dimensional ratio of radiated acoustic to transmitted electrical power), and and B i pθ, φq, 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]: ψ 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]: respectively, where "average" refers to averaging over one cycle of the harmonic signal waveform, at the frequency f.V T , V st R , and V v R 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 ptq and integration time τ p are used in single-target (e.g., calibration) and oceanic survey operations, respectively.Z T " R T `iX T is the transducer's input electrical impedance at transmission, when radiating into the fluid, with resistance and reactance R T and X T , respectively.Z R = R R + iX R is the output (internal) electrical impedance of the receiving transducer, and Z E " R E `iX 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, 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.
Dividing σ bs with a reference area (chosen equal to 1 m 2 ), 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, S v , respectively [22].
Equations ( 1) and ( 2) constitute formulation A of the power-budget equations for σ bs and s v .Physical interpretations of Equations ( 1) and ( 2) in terms of power flow are given in [22].

Power Budget Equations
Calculation of the electrical powers Π T , Π st R , and Π v R 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 σ bs 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 where all quantities being a function of the frequency, f, have been indicated by the argument p f q, and is the (non-dimensional) electrical termination factor for the power-budget equations formulated in terms of echo integration [16].F VV accounts for arbitrary electrical termination at signal reception [16].Integration over all frequencies gives Now, assume that the transmitted signal, V T ptq, 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 p f q| 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 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 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 ptq, starting at t = 0. Let τ st " t st2 ´tst1 be the time duration of the voltage echo V st R ptq that is received from a single scattering target, where t st1 and t st2 are the arrival times of the start and end of V st R ptq.Similarly, let τ sph " t sph2 ´tsph1 be the time duration of the voltage echo V sph R ptq that is received from the calibration sphere, where t sph1 and t sph2 are the arrival times of the start and end of V sph R ptq.Let τ g = t g2 ´tg1 be the "gate opening time" used for the voltage waveform V v R ptq that is received from the "gated volume", V g , where t g1 and t g2 are the times of gate opening and closure [16].
τ sph 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 sph1 and t sph2 , determining τ sph .In practice, a slightly different method may often be used.A time window is then set manually, starting just before t sph1 and ending just after t sph2 , 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 st , 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 st1 and t st2 to give τ st .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 ‰ τ st ‰ τ sph ‰ τ 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 giving where σ bs , s v , α, λ, G 0 , G pθ, φq, ψ, and F VV are all evaluated at the carrier frequency, f 0 , and are defined as the "echo integral"-or "time-integral-voltage-squared" [tivs]-values of the transmitted voltage signal (the "pulse"), V T ptq, and three different received voltage waveforms, V st R ptq, V sph R ptq, and V v R ptq, 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,σ bs and s v become essentially independent of pulse duration.By using the same transmit pulse V T ptq 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 V 2 s = W¨ohm¨s = J¨ohm.

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, σ sph bs,theory , calculated at the frequency f 0 in question [or the corresponding target strength TS sph theory " 10logpσ sph bs,theory {r 2  1 q, for a reference area r 2 1 = 1 m 2 ].Consider a calibration situation with the centre of the calibration sphere located at position ´rsph , θ sph , φ sph ¯. (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 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 where C 1 is a constant and known (measured) value for a given calibration measurement.Therefore, by knowing σ 10logpG 0 q " 10logpG 0,old q `1 2 ¨"TS (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,old .σ .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 σ bs 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.

"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 σ bs " e 4αpr´r sph q ¨˜r r sph 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 σ bs and s v .

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 dr p " 1 2 c 0 τ p in the observation volume V obs .The volume backscattering from V obs , 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 obs [11,16,17]: representing the backscattering cross section per unit area, within V obs .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]: 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 obs per square nautical mile, is then given as [17] where xσ bs y 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 " where Equation ( 38) is in the form of the traditional echo-integrator equation [3,4,6,8,10,16], where C cal represents the calibration factor used in that equation.E is a "range integrated echo integral" for the observation volume V obs , obtained by integrating the sequence of echo integrals rtivss v R , each associated with a "ping volume" V p .For each V p at range r in V obs , 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 , rtivss v R 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 rtivss v R 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 cal was used solely as an unspecified calibration factor.Here, the functional relationship for C cal is derived, and fully given for small-amplitude conditions in terms of the echo sounder system parameters.

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 σ bs (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.

Power-Budget Equations
Using the approach described above with "peak voltage echo integration" in single-target TS measurements, σ bs and s v given by Equations ( 23) and ( 24) can equivalently be expressed as [22]: where Here [22], rtivss st R,max is a simplified (and erroneous) estimate of the single target echo integral rtivss st R , obtained by multiplying the transmit signal duration τ p with the constant squared voltage value, ˇˇV st R ˇˇ2 max .ˇˇV st R ˇˇmax is the maximum ("peak voltage") value of the magnitude of the single target echo, ˇˇV st R ptq ˇˇ, within the time window [t st1 , t st2 ].G 2 e f f pθ, φq represents an "effective" value of G 2 pθ, φq, where a correction factor has been included to compensate for the error introduced by using rtivss st R,max instead of rtivss st R , to ensure that the resulting σ bs measurement is valid and correct [22].This approach represents a "rectangular echo integration" in single-target TS measurements, i.e., for σ bs .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 p0, 0q.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 "S v transducer gain" (cf.Section 4.3).
The calibration factor G 0,e f f can be expressed by two equivalent expressions [22]: Method 1: Method 2: 10logpG 0,e f f q " 10logpG 0,e f f ,old q `1 2 ¨"TS Here, G 0,e f f ,old is an old (or arbitrary, and possibly erroneous) calibration value used for G 0,e f f initially in the calibration measurement.σ sph bs,old is the corresponding old measured value used initially for the backscattering cross section of the calibration sphere, corresponding to target strength TS sph e f f ,old " 10logpσ sph bs,old {r 2 1 q [22].Similarly, the calibration factor G 0 can be expressed by two equivalent expressions [22]: Method 1: Here, G 0,old is an old (or arbitrary, and possibly erroneous) calibration value used for G 0 initially in the calibration measurement.s sph a,old is the corresponding old measured value used initially for the area backscattering coefficient of the sphere [22].s sph a,theory is a known value of the area backscattering coefficient of the sphere, given from calibration measurements and the known (calculated) value σ sph bs,theory .(A similar expression as Equation ( 51) is given in [11], but without ˇˇB i ´θsph , φ sph ¯ˇˇa nd the factor τ sph {τ p .That expression thus relates to the special case of calibration sphere located on the acoustic axis, θ sph " φ sph " 0, and τ sph " τ p (where the latter condition may not be valid in practice).)Equations ( 41) and (42) for σ bs 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.

Power-Budget Equations
By defining s a,corr " G 0 {G 0,e f f , σ bs and s v given by Equations ( 41) and (42) can equivalently be expressed as [22]: where Equation ( 52) is identical to Equation (41), and [22]: From Equations ( 9) and ( 26) it follows that giving, from Equations ( 43) and ( 54), Here, V st,rms R is the effective (rms) amplitude of the echo received from a single target (e.g., a fish, or a calibration sphere), and Π st R,max is the "maximum average electrical power" of that echo.Equation (57) shows thats 2  a,corr represents the (non-dimensional) ratio of the actual average electrical power, Π st R , in a single-target echo (i.e., from a calibration sphere, or a fish), to the "maximum average electrical power" of that echo, Π st R,max (i.e., the average power obtained by using "peak voltage echo integration", instead of full waveform echo integration).

Formulation D Calibration Factors G 0,e f f and s a,corr
Formulation D involves two calibration factors, G 0,e f f and s a,corr .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,corr 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,corr , that can be used for comparison with the two EK60 calibration factors "Gain" and "S a correction" (cf.Section 4.4).For s a,corr , 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,corr can be expressed by three equivalent expressions [22]: Method 1: S a,corr " 10logpG 0 q ´10logpG 0,e f f q. (62) Equations ( 52) and (53) for σ bs and s v , together with Equations ( 45)-( 47) and ( 58)-(62), respectively, for the calibration factors G 0,e f f and s a,corr , constitute Formulation D of the functional relationship for fish abundance estimation.

Power-Budget Equations
By replacing τ g ¨s2 a,corr with τ e f f in the s v equation, σ bs and s v given by Equations ( 52) and (53) can equivalently be expressed as [22]: where Equation ( 63) is identical to Equations ( 41) and (52).From Equations ( 54) and (57) it follows that [22] 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, Π st R , to the "maximum average electrical power", Π st R,max , 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 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 st R ˇˇmax and τ e f f , provided τ st 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: Method 2: Method 3: Here, τ re f is a reference time interval for τ e f f (chosen equal to 1 s).τ e f f ,old 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 σ bs 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.

Results and Discussion
Four equivalent echo integration Formulations B-E of the power-budget equations for σ bs 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.
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 σ bs 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.

Formulation B (Echo Integration, Generic)
Formulation B, Equations ( 23) and (25) for σ bs 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 VV .It also involves several integration time intervals, τ p , τ st , τ sph , 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 σ bs 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 σ bs 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 (σ bs ), 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 (σ bs ) 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 " 10log 10 pG 0,e f f q and S v transducer gain " 10log 10 pG 0 q, and using Equation (51), Equations (47) and (50) become identical to the corresponding expressions postulated for the EK500 calibration factors "TS transducer gain" and "S v 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 σ bs and s v given in [11], the analysis shows that the expressions for σ bs 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 σ bs 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].
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 ´2 a,corr .This indicates that the expressions for σ bs 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,corr in a consistent way.In Formulations C and D, the expressions for σ bs 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,corr .By defining Gain" 10log 10 pG 0,e f f q and S a correction " S a,corr " 10log 10 ps a,corr q, 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 σ bs 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, rtivss T , rtivss v R , and rtivss st R,max (in units of V 2 s = W¨ohm¨s = J¨ohm).Echo integration was thus not accounted for in these expressions.Secondly, the same Π R was used for σ bs 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, Gpθ, φq and G 0 were used instead of G e f f pθ, φq and G 0,e f f in the expressions for σ bs and s v , respectively, which is inconsistent when "peak voltage echo integration" is used for σ bs , and s a,corr 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 Gpθ, φq 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, Gpθ, φq 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 G 0 , 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, G eff (θ, ϕ) and G 0,eff , with a correction factor S a,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,corr parameter; and (c) "harmonize the TS and s A measurements" [19] by employing the same calibration factor G 0,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,corr , 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,corr , 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 "S v transducer gain", and the EK60 etc. calibration factors "Gain" and "S a correction", also appear to have been established and explained.
As the expressions for σ bs and s v derived in [23] correspond exactly to those postulated in [15], the above discussion also applies to [23]  The reasons for introducing the parameter G e f f , s a,corr , 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 (σ bs ) and in survey operation (s v ), only a single calibration factor, G 0 , applies, with no need for G e f f , s a,corr , 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 σ bs , which is compensated for by introducing G 0,e f f instead of G 0 in σ bs , 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 "S v 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 σ bs and s v equations.An error is then introduced in s v , which is compensated for by using s a,corr in s v , referred to as the "S a correction", cf.Equation (53).
In Formulation E, the term τ g s 2 a,corr used in Formulation D (in the s v equation) is denoted τ e f f .As a result, the σ bs 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 Gpθ, φq, G 0 , and τ g .In Formulation E these are Gpθ, φq, 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,corr , 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 rtivss st R,max instead of rtivss st R 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".

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.

Conclusions
On the basis of an average power formulation of power-budget equations for σ bs and s v known from prior literature [13,16] (denoted Formulation A), four new echo integration formulations for σ bs 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,corr , 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,corr , 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 σ bs 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 σ bs 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 rtivss T , rtivss v R , and rtivss st R,max .Similarly, the transducer gains Gpθ, φq and G 0 used by [15,23] are to be replaced by G e f f pθ, φq 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 σ bs 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,corr and τ e f f , which were postulated by [15] and adopted by [23].s a,corr 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 σ bs , 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 σ bs 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 σ bs 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.
sph bs,theory , the relevant calibration value G 0 is given by G 2 0 " C 1 {σ sph bs,theory .Consider a situation where an old (or an arbitrary, and possibly erroneous) calibration value is used initially in the calibration measurement, denoted G 0,old .The corresponding value measured initially for the backscattering coefficient of the calibration sphere is then σ sph bs,old " C 1 {G 2 0,old .It follows that C 1 " .

4. 5 .
On the G e f f , s a,corr , and τ e f f Parameters

Table 1 .
Alternative and equivalent formulations of power budget equations for measurement of σ bs and s v .Calibration factors referred to in the table are defined in Sections 2 and 3.
These quantities are all evaluated at the carrier frequency of the transmitted signal, f 0 .