Finite-Amplitude Power Budget Equations for Acoustic Fish Abundance Estimation

: Finite-amplitude (nonlinear) sound propagation e ﬀ ects in seawater may cause measurement errors in ﬁsh and zooplankton abundance estimation and species identiﬁcation for accessible echo sounder transmit electrical power levels and operating frequencies of about 100 kHz and higher. A su ﬃ ciently validated framework to quantify, control, and compensate for such errors in these applications is not available. The conventional power budget equations in ﬁsheries acoustics are valid for small-amplitude signals only. The study aims to ﬁll this “gap”. The conventional theory is generalized to account for ﬁnite-amplitude incident sound propagation, arbitrary electrical termination, and the range of electrical and acoustical echo sounder parameters. Equations for use in calibration and oceanic surveying are derived in terms of the backscattering cross section, σ bs , and the volume backscattering coe ﬃ cient, s v . The “ﬁnite-amplitude terms” in these expressions can—for relevant transmit electrical power levels of relevant echo sounders—be measured in controlled tank experiments. Alternatively, they can be calculated using numerical models. The resulting equations enable estimation of ﬁnite-amplitude measurement errors in these applications; development of recommended upper limits for echo sounder power levels; controlled reduction of ﬁnite-amplitude errors in calibration and surveying; and development of correction factors for survey data already subjected to such measurement errors.


Conventional Acoustic Fish Abundance Estimation
Acoustic methods for estimating fish stock abundance have been in regular use for several decades , and constitute a key element in national and international regulations of marine resources, such as fish, zooplankton, and krill.For fish aggregated in schools or layers, abundance measurement is based on echo integration [2,4,[6][7][8][9]11], supported by biological sampling.The acoustic methods rely on power budget equations and calibrated echo sounder and sonar systems [3,6,21,27,28].Echo-integrated measurements are used in expressions for the volume backscattering coefficient, v s , and the echo-integrator equation [4,6,7,11], to estimate stock abundance in terms of target (fish) density, a  [10,21,27].Fish abundance is measured using signal frequencies typically in the 10-120 kHz range.Zooplankton measurement employs frequencies above 100 kHz.For Atlantic mackerel (Scomber scombrus L.), the use of 200 kHz is reported to be advantageous in abundance estimation to discriminate against herring (Clupea harengus) [12].Over the recent decades, new developments have been taken into use for zooplankton, krill, and fish species identification (target classification) using multi-frequency measurement data [12][13][14]25,30].The methods are based on interpretation of the relative frequency response of v s , , where v s is measured over a range of frequencies, f, and normalized to its value at, e.g., 38 kHz.Typically, echo sounder frequencies in the range 10-500 kHz or higher are employed.
Recent work has shown that finite-amplitude sound propagation effects may cause errors in abundance and species identification measurements under certain operating conditions in fisheries acoustics [15][16][17][18][19][20].The development of a theory to describe these errors is the subject of this work.
The principal acoustic quantity used in abundance estimation and species identification is v s .
For at-sea calibration of scientific echo sounders prior to oceanic surveys, using a standard target (normally a metal sphere), a related power budget equation for single-target backscattering is utilized, given in terms of the backscattering cross section of the single target, bs  [11,20-28].The same equation is employed in target strength (TS) measurements of individual fish.

bs
 thus applies to the single-target backscattering used in echo sounder calibration and TS measurement of individual fish.v s and a  apply to the multiple-target volume backscattering used in oceanic (field) surveys.The power budget equations for bs  , v s , and a  referred to above, expressed in average power formulations, , will here be referred to as the "conventional generic power budget equations" in fisheries acoustics [21,25] 1 .(The quantities involved in Equations ( 1)-( 3) are defined in Sections 2-3 and summarized in Appendix D.) In this context, "generic" means "instrument independent".Equations ( 1) and (2) give the echo sounder's electrical power transfer function (ratio of received to transmitted average electrical powers) for situations with single-target and multi-target (volume) backscattering, respectively, expressed in terms of bs  and v s .The power "flow" represented by Equations ( 1) and ( 2) is explained in Appendix B of [23] (cf.also Appendix A).Equations ( 1)-(3) constitute the generic fundament for abundance estimation, species identification, and target classification in modern fisheries acoustics, serving as the basis for at-sea echo sounder calibration and survey operation [28][29][30].

Finite-Amplitude Effects in Fisheries Acoustics
In fisheries acoustics, small-amplitude (linear) sound propagation has been an underlying assumption from the emergence of fishery research echo sounders in the 1930s until recently.Equations ( 1)-( 3) rely on the assumption that no finite-amplitude (nonlinear) effects influence the propagation of sound through seawater.For typical transmit electrical power levels, e.g., some hundred W to about 2 kW, preliminary investigations indicate that this assumption is reasonably good at the lower end of the operational frequency range, such as at 18 and 38 kHz 3 .Finite-amplitude sound propagation effects increase, however, with increasing frequency.For operating frequencies of about 100 kHz and above, such effects may represent a source of measurement error in fisheries acoustics. 2Other small-amplitude power budget equations for and than those given by Equations ( 1) and ( 2) [21] have been proposed [11,20,22,[26][27][28].The various expressions proposed are not necessarily equivalent nor consistent.The expressions given in [20] and [22][23][24] represent further developments of [21] to account for arbitrary electrical termination, and are otherwise consistent with [21], cf.[23,24].The equations postulated by Simmonds and MacLennan [11] [their Equations (3.13) and (3.15)] may be shown [23,24] to differ from those given elsewhere [20][21][22][23][24][25][26][27][28], and power flow balance is not quite preserved (cf.Section 4.1.1).Demer and Renfree [26] used expressions for and that correspond to those of [21], where however the expression for is limited to the acoustic axis [23,24].Ona et al. [27] postulated alternative expressions for and .A derivation of those expressions [27] was presented in [28].The expression for proposed by [27,28] is however not consistent with those given in [20][21][22][23][24], and is not found to be valid [23,24].These issues are all addressed and resolved elsewhere [23,24].
Small-amplitude power budget equations analogous or corresponding to and/or have been addressed in some textbooks, such as [2,9,[31][32][33][34].By Clay and Medwin [2] and Medwin and Clay [9] only "in-water" expressions for and were discussed (i.e., in terms of sound pressures instead of electrical voltages), without accounting for the electroacoustic conversion of the transducer at transmission and reception.Apart from that, the expressions given by Clay and Medwin in [2,9] have been shown to be consistent with Equations ( 1) and ( 2) [22].Similar "in-water" expressions for were given by Garuthers [31], Burdic [32], Lurton [33] and Bjørnø [34].
3 At-sea target strength measurements conducted with the Simrad EK60/38 kHz echosounder system in seawater, using a Cu60 copper sphere reflector at a constant range 15.5 m below the transducer (the transducer being 8 m below the surface), and several transmit electrical power settings in the range 200-2000 W, have shown no sign of significant nonlinear loss.The deviations between target strength measurements at 200 W and 2000 W were less than 0.01 dB.(Pers.comm., R.J.; Korneliussen, Institute of Marine Research, Bergen, Norway, 2011.)Preliminary and tentative simulations using the Bergen Code solution of the KZK equation [38][39][40][41][42]35,16,17,19,20] indicate that-for 2000 W transmit electrical power-nonlinear loss in seawater at 38 kHz operating frequency is less than 0.08 and 0.25 dB, at 10 and 100 m ranges, respectively.(Pers.comm., Pedersen, A., Christian Michelsen Research AS, Bergen, Norway, 2011).Laboratory measurements, oceanic survey measurements, and numerical simulations have demonstrated that finite-amplitude effects in the sea may cause errors in bs  , v s , and thus a  measurements used in fisheries acoustics [15][16][17][18][19][20].During the development of new echo sounders based on composite transducer technology, offering improved dynamic range, efficiency, and bandwidth, signal distortion was experienced at transmit electrical power levels commonly used with earlier types of echo sounders.Investigations have demonstrated, first experimentally [15], and later by simulations [16,17,19,20] and other experiments (in the laboratory and oceanic surveying) [19,20], significant effects caused by finite-amplitude sound propagation in seawater, for operating frequencies in the range of about 100 kHz and above, using echo sounders operated at accessible transmit electrical power levels (some hundred W and higher).The demonstrated effects are at a level that may bias abundance estimation and species identification [20].For example, for a 200 kHz Simrad ES200-7C echo sounder calibrated and operated with the 1 kW power setting, and calibrated with the calibration sphere at 20 m distance, calculations indicate errors in v s (and thus for the abundance estimate) in the range of about 10% to 23% (−0.46 to −1.13 dB), due to finite-amplitude effects, for target ranges r = 20-300 m (cf. Figure 7.4 in [20]).The errors were shown to increase with increasing power and increasing calibration distance.Finite-amplitude sound propagation effects refer to the nonlinear distortion an acoustic pressure signal experiences as it propagates through a fluid [35].As the sound velocity in the fluid depends on the acoustic pressure, the positive pressure peaks of the signal waveform (crests) travel faster than the low-amplitude portions and the negative pressure peaks (rarefactions).This results in a distortion of the signal waveform during propagation in the fluid, and transfer of energy from the fundamental to the higher-frequency components of the signal's frequency spectrum.This effect has been referred to as "nonlinear harmonic distortion" [2,9].The excess pressure loss experienced at the fundamental frequency component due to finite-amplitude effects, in excess of the geometrical spreading and loss caused by fluid absorption and scattering, is referred to here as "nonlinear loss".It increases with increasing sound pressure level, and therefore with the increased electrical power level used for the transducer at transmission.Nonlinear loss and harmonic distortion accumulate with distance from the source until the signal amplitude eventually becomes so small that further finite-amplitude effects are negligible [16,17,19,20,35].
In a transducer sound field, the sound pressure is usually highest along the main lobe.As nonlinear distortion effects are amplitude-dependent, such effects are thus at their strongest along the main lobe, and weaker elsewhere.Nonlinear loss of the transmitted (fundamental) frequency component thus gives rise to a "flattening" and thus "widening" of the main lobe for that frequency component [2,9,16,17,19,20,35].An example can illustrate some typical figures in this respect.In [19], measurements and simulations have been reported for a Simrad EK60 transducer operating at 200 kHz in water.In the far field (about 22 Rayleigh distances), the −3dB beamwidth at the fundamental frequency was measured to increase from about 7.5° to about 9.9° using "small" and "high" amplitudes, respectively (150 and 1500 W electrical power settings), i.e., an increase of 2.4° (32%) 4 . 4These directivity figures apply to a single frequency (the fundamental frequency).Most fish-finding sonars transmit about 10% bandwidth in each pulse.Lurton [33] have discussed the directivity of wideband signals and given expressions for a non-weighted linear array.It was demonstrated that the major lobe is not severely affected by the frequency bandwidth, while the sidelobes are decreased by an amplitude modulating term, and the beampattern throughs are smoothed (filled in) at a degree depending on frequency bandwidth.Although these results for a linear array are not directly applicable to circular echo sounder transducers, it may be expected that similar impacts of frequency bandwidth apply also to circular transducers, at least qualitatively.The figures given in the main text for the increased -3dB beamwidth at the fundamental frequency, using "small" and "high" amplitude, and the corresponding flattening of the major lobe, are thus expected to be representative also for the directivity of a narrowband pulse transmitting e.g., 10% bandwidth.
The conventional expressions for bs  and v s , Equations ( 1) and (2) [21,25], and alternative small-amplitude expressions [11,20,[22][23][24][25][26][27][28], do not account for finite-amplitude sound propagation effects.Thus, if not avoided or corrected for, such effects may cause errors and problems for today's methods of abundance estimation and species identification.Measurements are based on the fundamental frequency component of the received signal.When measurement data subjected to significant finite-amplitude effects are used as input to the conventional expressions for bs  and v s , Equations ( 1) and ( 2), respectively, which are based on small-amplitude (linearized) sound propagation theory, abundance estimates are biased.Similarly, in fish species identification (target classification), possible errors in ) ( f R v at the higher frequencies, due to error in v s , bs  , or both (in survey and calibration, respectively), may lead to erroneous interpretation, when comparing with the ʺreference signature", , for candidate species.In these applications, the control of possible errors caused by finite-amplitude effects is important.
To gain reliable insight into how finite amplitude effects affect fish abundance estimation and species identification, such operations need to be modeled using a full functional relationship for the measurements, where all electrical and acoustical parameters influencing the measurements are accounted for, including the delicate balance between echo sounder calibration measurements and field (survey) operation.The influences of finite-amplitude effects appearing in one or both of the calibration and survey operations of the echo sounders need to be accounted for.
A method to estimate the magnitude of finite-amplitude errors, to avoid or correct such errors, has been proposed and used in numerical calculations by Pedersen [20] (his Equations (7.1)- (7.11)).The conventional small-amplitude power budget equation for v s , Equation ( 2), was generalized to account for finite-amplitude sound propagation.An expression for v s in this form was not explicitly given, but indicated (cf.Equations (7.1)-(7.5) in [20]).On this basis, a finite-amplitude correction factor was proposed and used as a multiplicative factor to the measured v s value (i.e., the erroneous value obtained under the assumption of small-amplitude sound propagation).The analysis accounted for possible finite-amplitude effects during calibration without, however, giving an explicit finite-amplitude expression for bs  .A sufficiently detailed derivation of the correction model was not given, and the expressions constituting the model may possibly not appear highly intuitive.An alternative theory confirming Pedersen's correction model has been presented [36,37], where the mathematical derivations leading to the finite-amplitude power budget expressions for bs  and v s were omitted.
Consequently, the magnitude, importance, and consequences of the error introduced by finite-amplitude effects in abundance estimation and species identification have not yet been sufficiently described and documented.Methods to estimate and compensate for such errors have been presented, accounting for finite-amplitude effects in at-sea echo sounder calibration and survey operation [20,36,37].Recommendations for reduced echo sounder source levels and transmit electrical power levels have been proposed [11,14,18], e.g., based on calculations given elsewhere [16,17,20].However, providing reliable and optimal recommendations for compensation and transmit power levels is difficult without a sufficiently complete and documented analysis of the subject.
The average power model proposed in [20], or the alternative model given in [36], could potentially have served as a candidate method for such analysis and recommendations.However, until mathematical derivations of the proposed expressions are presented, these models cannot be considered sufficiently validated for a reliable analysis.
For control in abundance estimation and species identification, there is thus a need to present reliable and documented power budget expressions for bs  and v s that account for finite-amplitude signal propagation.This need relates to error analysis; the establishment of recommended upper limits for source or transmit electrical power levels; and the development of correction methods for finite-amplitude effects.The derivation of such expressions for bs  and v s is the topic addressed here.

Objectives and Outline
The objective of this article is to derive power budget equations-for single and multiple-target (volume) backscattering, and for the fundamental frequency component of the received signal-that account for finite-amplitude sound propagation in the fluid medium (seawater).
A comment may be needed to clarify the applicability of the analysis.For classical narrowband operation of scientific echo sounders, with narrowband transducer and/or filtering used at transmission and reception (such as using, e.g., Simrad EK500, EK60, and EK80 in "EK60 modus", or similar equipment), the harmonic content of the scattered signal is negligible, and the present analysis is expected to apply.
In modern echo sounder systems (such as Simrad EK80), with possibilities for wideband operation using "chirp" frequency-modulated signals (matched filtering), some of the harmonic contents that may be generated in the case of finite-amplitude signals with carrier frequencies in the range of 100 kHz and above may be received within the wide frequency band employed.If the full wideband waveform signal is used in signal processing, a more complete analysis would be necessary than the one presented here.In that case, the lower harmonic frequency components would have to be accounted for (represented by second moments), as well as frequency dependent absorption, and possibly also wave dispersion.
The handling of possible finite-amplitude effects in wideband operation of fisheries echo sounders is a highly complex matter, and, in practice, precautions to reduce or preferably avoid such effects are preferred.In these applications, the theory presented here may be used, e.g., to establish upper limits for echo sounder power levels, so that finite amplitude effects are reduced or become negligible.
The paper is organized as follows.In Section 2, an electroacoustic power budget equation describing backscattering from a single target is derived (giving bs  ), accounting for finite-amplitude incident sound.From this expression, an electroacoustic power budget equation describing volume backscattering from a multitude of targets is derived in Section 3 (giving v s ), also applicable to finite-amplitude sound propagation conditions.The generic (instrument independent) expressions for bs  and v s are given in terms of average electrical powers, averaged over a single cycle of a monochromatic wave.In Section 4, the results are discussed in relation to prior literature, also summarizing the assumption underlying the analysis.Conclusions are given in Section 5. Interpretations of the derived power budget equations for bs  and v s in terms of average power "flow" are given in Appendix A. Appendices B and C give interpretations of important quantities involved, used for interpretation of the power budget equations and analysis of some prior literature (cf.Appendix A and Section 4.1.1).Symbols and nomenclature are summarized in Appendix D. This article provides a generalization of the conventional expressions given by Equations ( 1) and (2) [21,25] (as well as alternative small-amplitude expressions [11,20,[22][23][24][26][27][28]), to account for finite-amplitude incident sound, for the fundamental frequency component of the received signal.They also represent a generalization of Equations ( 1) and (2) [21,25] (and alternative models [11,[26][27][28]) to account for arbitrary electrical termination at reception.
The description represents a generalization of [22] to account for finite-amplitude sound propagation effects in the fluid medium, and is therefore chosen to relatively closely follow the same sequence of derivation and also some phrasing used in [22], to enable convenient identification of essential similarities, differences, and assumptions in comparison of the two cases of small and finite-amplitude sound propagation in the fluid medium.

Single-Target Backscattering of Finite-Amplitude Incident Sound
A frequency domain description is used, with time harmonic factor t i e  , where is the angular frequency of the monochromatic wave, and t is the time.Bold-face letters are used to indicate complex-numbered quantities, and vectors are represented by underlined characters [22,23].

Acoustic Backscattering from a Single Target in the Farfield
Consider the situation shown in Figure 1.An electrical signal at angular frequency  is fed to an electroacoustic transducer, by which it is converted to an acoustic pressure wave, and radiated into a homogeneous fluid medium, i.e., with constant density and sound velocity [22].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 [14,28].In the fluid, at an arbitrary orientation relative to the transducer (on or off the acoustical axis), consider a single object of unspecified shape and material, or alternatively, a multitude of such objects, of different types, materials, and sizes [22].It is assumed that these are confined to a sufficiently small volume in space, so that the sound backscattered from the object(s) to the transducer appears as if the scattering came from a single target 5 .This object, or small volume of objects, can then be treated as a single target, and is for convenience referred to as "the target", or "the scatterer".A monostatic situation is considered, where the backscattered sound pressure wave is received by the same transducer and converted to an electrical signal.
Figure 1.Sketch of the acoustic system for single-target backscattering, with an electroacoustic transducer operating as a transmitter and receiver of ultrasound, acoustic backscattering from a single scattering object (target) in a homogeneous fluid medium, and the two spherical coordinate systems 1 and 2 used for the transmitted and scattered sound wave fields, respectively.The target center is located at position ( , , ) r   relative to coordinate system no. 1. (Reproduced from [22], with permission from the Institute of Marine Research, 2020.) Two coordinate systems are used to describe this electroacoustic system [22].The origin of coordinate system no. 1, used for the transmitted wave field, is located at the center of the front face of the transducer.The z-axis is chosen normal to the transducer's front surface, and assumed to be co-incident with the transducer's acoustic axis.Coordinate system no.2, used for the scattered wave field, and employing primed coordinates, has its origin located at the center of the target, and the axes parallel to the respective axes of coordinate system no. 1. Figure 1 shows the Cartesian coordinates x, y, and z (x', y', and z') and the spherical coordinates r, θ, and φ (r', θ', and φ') for the two coordinate systems, where r (r') is the radial distance, denoted range, θ (θ') is the polar angle, and φ (φ') is the azimuthal angle.The position vectors in the two coordinate systems are ) , , ( 5 This distinction is important to enable the use of the single-target expression, Equation (50), in the integration over a multitude of scattering objects contained in a volume of finite extent, to derive the volume backscattering coefficient, , given by Equation (56).The same assumption is implicitly underlying the conventional theory [21] for abundance estimation, Equations ( 1)-(3).
Consider a target located in the farfield of the transmitting transducer, where farfield refers here to small-amplitude sound propagation.It is assumed that finite-amplitude sound propagation effects only affect the forward-radiated (transmitted) sound pressure wave (referred to here as the incident sound wave) [20,36,37].At the position of the target, the amplitude of the incident wave is so small that finite-amplitude effects in seawater can be neglected.It is further assumed that possible nonlinear effects in the scattering process at the target itself (involving, e.g., fish with a gas-filled swim-bladder) can be neglected [22][23][24], so linear backscattering theory [2] applies.Consequently, the backscattered wave amplitude is so small that finite-amplitude effects in the scattered field can be neglected, and the scattered field is described by small-amplitude (linearized) theory.In the farfield of the target, the scattered pressure field spreads inversely proportional to range, ' r .Under these assumptions, and motivated, e.g., by numerical simulation results for the finite-amplitude sound pressure field radiated by baffled piston sources [16,17,19,20] using the Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation [35,[38][39][40][41][42] the fundamental frequency components of the finite-amplitude incident pressure wave, n i p , and the scattered pressure wave, s p , are modeled as respectively.Subscripts "i" and "s" are used for incident and scattered waves, respectively, and the superscript "n" for quantities directly subject to finite-amplitude effects.

P
. The functions are the beam patterns of the incident and scattered sound pressure waves, respectively.) 0 , 0 , (r the axial sound pressure amplitude (along the z axis) for the finite-amplitude incident sound field.
is the sound pressure amplitude along the z' axis for the scattered sound field.is the acoustic wavenumber in the fluid medium, 0 c is the small-amplitude (also called small-signal) sound velocity of the fluid [35], and  is the acoustic attenuation coefficient of the fluid for the sound pressure. accounts for acoustic absorption of the fluid [43] and possible excess attenuation due to volume scattering experienced in sound propagation between the transducer and the target [2,9].(In the conventional power budget equations, Equations ( 1)-( 3), and in practical abundance estimates, only acoustic absorption is normally accounted for by  .)For small-amplitude signals, the amplitude of the incident sound pressure is given as [43] ) , ( ) , , ( where are the beam pattern and axial sound pressure amplitude of the incident wave, respectively, under such small-amplitude sound propagation conditions.
The "axial finite-amplitude factor" of the incident pressure wave, defined as represents the deviation from spherical spreading and attenuation (i.e., absorption and scattering) along the acoustical axis (i.e., the deviation from the small-amplitude case), as a measure of axial nonlinear loss.

C
for small-amplitude incident waves, and

C
under finite-amplitude sound propagation conditions.) (r n i

C
depends on range, r, since nonlinear loss increases with increasing r, until the wave amplitude eventually becomes so small that further nonlinear loss is negligible [16,17,19,20,35].Hence, at long ranges

C
becomes approximately independent of r , and thus approximately constant, for a given frequency and source level.
The use of ) 0 , 0 , ' (r s P as the normalization pressure amplitude in Equation ( 9) may need comment, since the z' axis is not necessarily the direction of maximum scattering.This approach has been chosen for convenience and without any loss of generality, since the results derived in the following become independent of the choice of normalization direction for ) ' , ' (   s B [22].
is not.This is due to the finite amplitude of the incident pressure wave, with flattening of the main lobe relative to the corresponding small-amplitude sound wave, which changes with distance from the transducer [16,17,19,20,35].At long ranges, the wave amplitude becomes so small that further flattening is negligible, and ) , , (

B
becomes approximately invariant to r.Under conditions of small-amplitude sound propagation, ) , , (

B
becomes independent of r, and reduces to its range-independent small-amplitude counterpart, ) , From Equations ( 5) and ( 10), the incident sound pressure amplitude can be expressed as where the "beam pattern finite-amplitude factor", defined as represents the finite-amplitude effects on the beam pattern relative to the small-amplitude conditions (i.e., the deviation from the small-amplitude case).

B
depends on range, r, since the flattening of the main lobe due to finite-amplitude effects increases with increasing r until the wave amplitude eventually becomes so small that further flattening is negligible [16,17,19,20].Hence, at large ranges, ) , , (

B
becomes approximately independent of r ,  , and  , and thus approximately constant, for a given frequency and source level.
From Equations ( 5) and ( 12), the incident free-field pressure amplitude at the target position can be written as where is the axial sound pressure amplitude, at the axial reference range 0 r (e.g., 1 m) from the transducer front, under small-amplitude sound propagation conditions, extrapolated spherically from the farfield.


, is given as where is the scattering function 6 [2], and A is the cross-sectional area of the scattering target, viewed from the transducer.Note that s S is independent of range, r, since the ratio From Figure 2, the backscattering direction is given by are the backscattering cross section of the target (m 2 ) and the backscattering function, respectively [2].bs  depends in general on frequency, the direction ) , (   of the incoming (incident) wave, and the shape of the target.Note that in this description-since bs  is a property of the scatterer-terms accounting for attenuation and spherical spreading of the backscattered field (between the scattering object and 0 ' r ) are not to be included in bs  .These are accounted for in 0 , bs I , as seen from Equation (22).
By combining Equations ( 18), (19), and ( 23  From Equation ( 7), the backscattered pressure amplitude at a range r', referred to coordinate system no.2, can be written as By inserting Equations ( 15), (16), and (25) into Equation ( 26), and setting r' equal to the target range, r, the magnitude of the amplitude of the backscattered free-field sound pressure in the fluid, at the center of the transducer front, ) ' , ' , ( Equation ( 27

B C
 accounts for finite-amplitude effects on the incident axial pressure and beam pattern, and represents the deviation from the small-amplitude sound pressure.

Electroacoustic Transmit-Receive Transfer Functions for Single-Target Backscattering
In the following, Equation ( 27) is used to develop electroacoustic transmit-receive transfer functions for backscattering from a single target in the farfield, by accounting for (a) the transmit electrical power, (b) the transducer's electroacoustic conversional efficiency, (c) the transmitting and receiving responses of the transducer, (d) the beam pattern upon reception, (e) farfield spherical reciprocity, (f) the input and output electrical impedances of the transducer, and (g) the input electrical impedance of the echo sounder's receiving electronics.
Assume that the transducer is linear, passive, and reversible, and fulfills the reciprocity relationships [44].The transmit voltage amplitude is thus assumed to be sufficiently small to avoid nonlinear effects in the transducer and the electronics.The transducer's axial transmitting current response, I S , and free-field open-circuit receiving voltage sensitivity, V M , are given as [43,44] T i I ) , ( 0 ) , is the beam pattern of the transducer upon reception, which is equal to the small-amplitude beam pattern upon transmission [44], and thus given by Equation (11).Insertion of Equations ( 28) and ( 29) into ( 27) leads to the magnitude of the transmit-receive current-to-voltage transfer function under open-circuit conditions: The transducer's (one-way) electroacoustic conversion efficiency under conditions of small-amplitude and lossless sound propagation in the fluid is defined as where [43] is the average electrical power delivered to the transducer during transmission, averaged over one vibration cycle of the monochromatic wave (here denoted "average transmit electrical power").T V is the voltage amplitude across the transducer's electrical terminals at transmission, and is the transducer's input electrical impedance when radiating into the fluid with resistance and reactance T R and T X , respectively, cf. Figure 3a.a  is the average acoustic power radiated from the transducer into the fluid medium, averaged over one vibration cycle of the monochromatic wave, under conditions of small-amplitude and lossless sound propagation in the fluid, given in the farfield as [43]  where is the unit solid angle.

D
denote the axial directivity factor [43] and the directivity factor, respectively, for the transmitted (incident) sound field under conditions of small-amplitude and lossless sound propagation.These are defined as the dimensionless ratios of the transducer's axial intensity, respectively the intensity in a given direction     , to the intensity of an omnidirectional (point) source generating the same acoustic power, and can be expressed as [43]  respectively.Transducer gain,     ,

G
, and axial transducer gain (or "peak gain" [21]), 0 G , are defined as [20,21,45]   and axial directions, respectively, for lossless and small-amplitude sound propagation conditions in the fluid (cf.Appendix B, interpretation 3).From Equations ( 5), ( 31), (33), and (37), one obtains Combining Equation (28) with Equations ( 32), (37), and (38) leads to Since 0 , i P on which I S is based (cf.Equation ( 28)) is extrapolated from the transducer's farfield, the spherical-wave reciprocity relationship applies, stating that 7 [44,46] where s J is the spherical-wave reciprocity parameter, giving the relationship between the transmit and receive sensitivities of the echo sounder transducer (under the stated assumption that the transducer is linear, passive, and reversible, fulfilling the reciprocity relationships [44]).f c 0   is the acoustic wavelength in the fluid medium.Equations ( 29), (39), and (40) give Insertion of Equations ( 39) and ( 41) into Equation (30) gives, for the magnitude of the open-circuit transmit-receive transfer function, where 7 Foldy and Primakoff [44] considered a lossless fluid medium, using time dependence exp(it), as here, and the lossless version of Equation ( 40) was given, corresponding to  = 0 [cf.their Equation ( 50)].By repeating their derivation with fluid absorption accounted for, using a complex wavenumber  = k − i in the fluid instead of the real wavenumber which was used by Foldy and Primakoff [44], Equation ( 40) results.
To include the effects of non-ideal electrical termination at the receiver (i.e., a finite electrical termination load), consider the situation indicated in Figure 3b, which can be represented electrically by the Helmholtz-Thevenin equivalent circuit shown in Figure 4. Here, is the output (internal) electrical impedance of the receiving transducer, and is the input electrical impedance of the receiving electric network, involving resistances and reactances R R and R X , and E R and E X , respectively.For the transducer, the assumption R T  Z Z is often used, but for generality in the description, R Z is here distinguished from T Z .R V and R I are the voltage and current amplitudes at the transducer's electrical terminals at reception for single-target backscattering.From Figure 4, R V is given by Insertion of Equation ( 44) and these electrical impedance definitions into Equation (42) yields ) , , ( for the magnitudes of four transmit-receive transfer functions of interest, for a single target located at position ) , , (   r in the farfield.Here, are electrical termination factors for the respective transfer functions in Equations (45).They represent the effect of the finite electrical termination load on the receiving transducer in terms of the electrical impedances of the transducer ( R Z ) and the receiving electronics network ( E Z ) [22].Note that VV F and II F are dimensionless, whereas IV F and VI F are given in units of Ω and Ω -1 , respectively.

Electroacoustic Power Budget Equation for Single-Target Backscattering
The average electrical power delivered by the transducer to the receiving electronics, averaged over one vibration cycle of the monochromatic wave (here denoted "average received electrical power"), for single-target measurements, is given as [43 Insertion of Equations ( 32) and ( 47) into the former of Equations (45a) yields the transmit-receive electrical power transfer function, , 2 ( , , ) where the (dimensionless) electrical termination factor for the electrical power transfer function is defined as [20] 2 Here, and in the following, the symbol T  used in Equation ( 32) has been replaced by st T  , to distinguish between the transmit electrical powers used in single-target and volume backscattering situations (cf.Equation (54) and the accompanying text).st T  is the average transmit electrical power used in single-target measurements (e.g., sphere calibration and fish TS measurement) [22][23][24].
From Equations ( 48) and ( 43), the backscattering cross section of the single target is given as ) Equation ( 48), or equivalently, Equation ( 50), is here denoted the "average power formulation" of the electroacoustic power budget equation for backscattering from a single target located at position ) , , (   r in the transducer's farfield, for the fundamental frequency component of the received signal, under conditions of finite-amplitude incident sound.Appendix A.1 gives an interpretation of Equation (50) in terms of power "flow".It is noted that, under the assumptions used here (cf.Section 4.2), bs  is a property of the scatterer, and is as such invariant to the sound pressure amplitude.The factor   ( ) ( , , ) in Equation ( 50) accounts for finite-amplitude effects in the axial pressure and beam pattern of the incident field, and represents the deviation from the small-amplitude case for bs  .Thus, accounts for the corresponding influence of finite-amplitude effects on the average received electrical power.

Volume Backscattering of Finite-Amplitude Incident Sound
Now, consider backscattering from a spherical shell volume in the farfield, obs V (denoted "observation volume"), between ranges min r and max r , as shown in Figure 5. obs V contains a distribution of scattering objects of different types (e.g., different types of fish, krill, and zooplankton).
In the present section, Equation ( 50) is used to derive an expression for the volume backscattering coefficient, v s , for a spherical shell sub-volume, p V , in obs V .

Electroacoustic Power Budget Equation for Volume Backscattering
In the following, assume that (a) the scattered echoes from different objects in obs V have random phases, (b) multiple scattering effects and interaction between objects can be neglected, and (c) excess attenuation from power extinction [2] caused by volume scattering in obs V (also referred to as the "shadow effect" [11]) can be neglected.Assumption (a) corresponds to random spacing of objects in one "ping", and movement of the objects to the next "ping" [2,11].Assumption (b) means that only echoes backscattered directly from the objects are significant, so that those backscattered via other objects (second-order effects) can be ignored [9,11,47].Assumption (c) may be a reasonable approximation except for strong scatterers at high target densities distributed over an extended volume 8 [9,11,[48][49][50].
For a multitude of small objects in obs V , the echoes from individual objects cannot be resolved, but combine to form a received voltage signal with varying amplitude.The echo intensity is still a measure of the biomass in the volume [1,4,11].Under the above assumptions, the total echo intensity is the incoherent sum of the individual echo intensities [9].The volume backscattering coefficient, v s , is the backscattering cross section per unit volume [2].Consequently, for a multitude of scattering object types, v s can be calculated as a sum over backscattering cross sections of the individual types per unit volume [2,50] so that 8 Excess attenuation from power extinction caused by volume scattering in (also referred to as the "shadow effect") have been discussed e.g., by [9,11,[48][49][50], and in references given therein.The shadow effect of caged fish at relevant fish densities, power levels, and four frequencies in the common operational frequency range used in fisheries acoustics, 38-120 kHz, was investigated by [50] (also reproduced and discussed in [11]), indicating errors in fish density (no. of fish per m 2 ) of typically less than 10%.In a couple of cases errors in the range 14-16% were found, and in one case about 29%.
In comparison with the figures given by Pedersen [20] (cf.Section 1.2), the impact of finite amplitude effects may thus be at the level of the shadow effect due to power extinction, or larger, depending on frequency, power setting, and calibration distance.
Here, N is the number of scattering object types, V  is the unit volume, is the number of scattering objects of type j per unit volume, j m is the number of scattering objects of type j in V  , and j bs,  is the backscattering cross section for an object of type j, j = 1, …, N.
From Equation (51), From Equations ( 51)-( 53), it can be seen that bs d represents the backscattering cross section of a multitude of objects in the (infinitely small) unit volume dV, including objects of different types and objects of the same type with different sizes.
As explained in Section 2, Equations ( 48) and ( 50) apply not only to a single scattering object in the farfield, but also to a multitude of farfield objects of different types, materials, and sizes, confined to a sufficiently small volume in space, so that the backscatter at the transducer appears as if the scattering came from a single target in the farfield.Let the unit volume dV represent such a farfield "effective single target" containing a multitude of scattering objects.For backscattering from dV at range r in obs V , Equations ( 48) and ( 43) thus yield  Assume a uniform distribution of scattering objects in the volume p V , so that bs d as given by Equation ( 53) can be used everywhere in p V .In practice, this means that backscatter is assumed to be the same anywhere in the "sampled volume" n s V , which represents that portion (sub-volume) of p V that is effectively insonified by the acoustic beam upon radiation and reception combined [23], cf. Figure 6 where 2 ) ( is the mid-radius range of p V . The simplifications made in Equation (55) may need a comment.As explained above, the unit volume dV represents a farfield "effective single target" containing a multitude of scattering objects, sufficiently small so that the backscatter at the transducer appears as if the scattering came from a single target.bs d represents the backscattering cross section of this "effective single target".In the derivation of Equation (55), the integration over "effective single targets" represented by bs d in the (thin) spherical shell volume p V has (by substituting Equation ( 53)) been replaced by volume integration using the unit volume   drd r dV 2 , so that the integration over p V is now made in terms of the range r and solid angle Ω of the spherical coordinate system.Since   drd r dV 2 , dV is limited in space by the unit thickness dr and the unit solid angle, . To arrive at the final expression of Equation (55), v s has been (a) assumed to be independent of Ω (i.e., the polar and azimuthal angles) within p V (in practice within the sub-volume n s V ), and (b) approximated by its mid-radius range of p V , p r r  .For this simplification to be valid, it has thus been assumed that the ratio bs d /dV (= v s ) is independent of Ω (i.e., the polar and azimuthal angles) within p V .Since dV contains a multitude of objects (including objects of different types, and objects of the same type with different sizes), this assumption corresponds to assuming a uniform distribution of the scattering objects within p V (in practice within n s V ).
For convenience, the subscript "p" in p r is hereafter understood and omitted (with some exceptions).As should be clear, integration over the range r in the finite volume obs V ( min r to max r ) has not been carried out at this stage.(This range integration is accounted for in the calculation of the fish density contained in obs V , a  , cf.Equation ( 3) [22]).

Volume Backscattering Coefficient
By rearranging Equation (55) to solve for v s , one obtains for volume scattering from the "ping volume" p V , at distance where the definitions and relationships have been used, and the latter of the two expressions given in each of Equations ( 57) and ( 58) follow from Equations ( 36), (37), and (15). and ) (r n  are the equivalent two-way beam solid angles of the transducer, under conditions of small-amplitude and finite-amplitude incident sound, respectively.Whereas  is commonly used, e.g., in the fisheries and marine acoustics literature [2,6,9,11,20-28,31-  34] and is normally provided by the echo sounder manufacturer, ) (r n  has been introduced [20] as a range-dependent generalization of  to account for finite-amplitude effects. (The equivalent two-way beam solid angle,  (also denoted the "integrated beam width [2], or "integrated beam pattern" [9]), represents the effective beam width of the transducer's intensity field in terms of a solid angle for the combined effect of transmission and reception [11].) The "beam solid angle finite-amplitude factor", ) (r n rel  , is introduced here as a measure of the transducer's two-way effective beam width (in terms of a solid angle), at finite-amplitude relative to under small-amplitude conditions.) (r n rel  depends on range, r, since the flattening of the main lobe due to finite-amplitude effects increases with increasing r until the wave amplitude eventually becomes so small that further flattening is negligible [16,17,19,20].Hence, at long ranges ) (r n rel  becomes approximately constant for a given frequency and source level. The latter expression in Equation (55), or equivalently, Equation (56), is denoted the "average power formulation" of the electroacoustic power budget equation for volume backscattering, for the fundamental frequency component of the received signal, under conditions of finite-amplitude incident sound.It applies to the thin spherical shell sub-volume p V (the "ping volume") of thickness

Consistency with Prior Literature
Consistency of Equations ( 50) and ( 56) with prior literature is described in the following for small-amplitude (Section 4.1.1)and finite-amplitude (Section 4.1.2)signals.

Small-Amplitude Signals (Linear Sound Propagation)
For a sufficiently small amplitude of the transmitted field so that the sound propagation is governed by the linearized set of acoustic field equations, ) 0 , 0 , (r       , i.e., the same electrical transmit power is used in calibration and survey operations, Equations ( 60) and (61) become identical to the expressions that were derived under the assumption of small-amplitude (linear) sound propagation [22][23][24].Furthermore, if 1 F   , Equations ( 60) and (61) become identical to Equations ( 1) and ( 2).As discussed in [22,23], this corresponds to specific cases of electrical termination: either (i) E  * R Z Z (i.e., conjugate matched electrical termination, to maximize power transfer from the transducer to the receiving electronic circuit); or (ii) E R  Z Z and 0 T X  Ω (to minimize signal reflections from the receiving electronic circuit in a frequency band close to the series resonance frequency of the transducer vibration mode used).Consistency with the conventional generic power budget equations is discussed first.In the Simrad EK500 manual [21], expressions for the spherical scattering cross section [2,10]  (somewhat incorrectly referred to as volume backscattering strength, which is defined [10] as ) were derived, where 1 r is a reference range (chosen to be equal to 1 m).These expressions correspond to Equations ( 1) and (2).Electrical termination and impedances were not addressed, which implicitly corresponds to setting  F equal to 1 as discussed above.Apart  and v T  , Simrad's expressions [21] are consistent with Equations ( 60) and (61).
Pedersen [20] was the first to account for arbitrary electrical termination at the receiver, and derived average power expressions for bs Regarding comparison with alternative power budget equations, alternative small-amplitude power-budget equations have been proposed in [11,[26][27][28].These theories are, however, either incomplete or not consistent with [20][21][22][23][24], as discussed in the following.
Simmonds and MacLennan [11] postulated power budget equations (their Equations (3.13) and (3.15)) that (when using  60) and (61), as well as with [20][21][22][23][24][25][26][27][28].This is due to a missing factor   4 2 involved in the transducer's "effective receiving area" (cf.Appendix A.2 and C), and a missing factor Demer and Renfree [26], with reference to the EK500 manual [21], gave expressions for the spherical scattering cross section bs  4 and the function (which somewhat misleadingly were referred to as the backscattering cross section and the volume backscattering coefficient, respectively).In their expression for , implying that the equation is limited to a single target (e.g., a calibration sphere, or an individual fish) on the acoustical axis only.Electrical termination was not addressed, implicitly implying 1   F as discussed above.
Ona et al. [27] postulated expressions for bs  and v s in logarithmic (dB) form.By converting their logarithmic expression for v s to normal units, it contains (in terms of the terminology used here) a factor )) corr a s  were referred to as the "nominal pulse duration" and the "integration correction", respectively.The apparent "ad hoc" introduction of corr a S , in the expression for v s was not explained nor motivated, except for a statement that "the sum of p  and corr a S , equals the effective pulse duration".In [23,24], the use of the factor corr a S , has been explained and justified.However, as also shown in [23,24], there seems to be some inconsistencies in the expressions given in [27] with respect to transmit and receive electrical power in relation to echo integration and the use of 0 As the expressions for bs  and v s derived and used in [28] correspond exactly to those postulated in [27], the above discussion also applies to [28], cf.[23,24].The mentioned issues with [11,[26][27][28] have thus been addressed and resolved [23,24].If the inconsistencies in [11,27,28] were corrected, and the limitations in [26] avoided, as derived and explained in [23,24], the expressions would become consistent with Equations (60) and (61).It has been found [23,24] that for small-amplitude sound propagation conditions, Equations (60) and (61) are the correct expressions under assumptions commonly used in fisheries acoustics (cf.Section 4.2).
It follows that Equations ( 50) and (56) represent a consistent generalization of single-target and volume backscattering theory for small-amplitude sound propagation conditions and conditions of finite-amplitude incident sound.

Finite-Amplitude Signals (Nonlinear Sound Propagation)
Power budget expressions for bs  and v s in average power form that account for finite-amplitude effects have been addressed by [20,36,37].
As discussed in Section 1.2, Pedersen [20] proposed a finite-amplitude expression for v s , cf.
Section 7.1.2in [20].A sufficiently detailed derivation of this model was, however, not given, and no explicit finite-amplitude expression was given for bs  , although nonlinear effects in echo sounder calibration were discussed and accounted for, cf.Section 7.1.2in [20].Lunde and Pedersen [36] presented expressions for bs  and v s which in the terminology used here become , where [36]       and the superscripts "T" and "R" denote one-way "transmit" and "receive" transducer gains, respectively.As mentioned above (Section 1.2), the mathematical derivation leading to these expressions [36] has not been published.Equation ( 63) can be shown to be equivalent to the expression for v s proposed by Pedersen [20]   to account for finite-amplitude incident sound.This is seen by introducing his Equations (7.3) and (7.4) into his Equation (7.1), as indicated in [20].The resulting v s expression was, however, not given explicitly by Pedersen.
From Equations (64)-( 68), (15), and ( 36)-( 37), one has ) Equations ( 62) and (63) are thus identical to Equations ( 50) and (56), respectively.It follows that the average power formulations of the power budget equations derived here are equivalent to the power budget expressions given (without sufficient detailed and documented derivations, however) for v s by Pedersen [20], and for bs  and v s by Lunde and Pedersen [36].
The present theory thus constitutes a derivation also of those expressions.

Summary of Assumptions Underlying the Analysis
In Sections 2 and 3, a number of assumptions have been used to derive the expressions for backscattering cross section, bs  , and volume backscattering coefficient, v s , Equations ( 50) and (56), respectively.With some exceptions (cf.(d) below), these are the same as those summarized in [22][23][24] for small-amplitude sound propagation.The assumptions used here include: (a) the monostatically operated transducer is passive, reversible, and reciprocal, fulfilling the reciprocity relationships [44]; (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 is homogeneous, with constant density and sound velocity; (d) finite-amplitude sound propagation effects in seawater are influent only for the forward-radiated (incident) sound pressure wave; (e) targets are stationary during a single transmission, and of any shape; (f) targets are in the farfield of the transducer; (g) possible nonlinear effects in the scattering process at the target itself (involving, e.g., fish with a gas-filled swim-bladder) can be neglected, so linear backscattering theory applies; (h) the volume backscattering coefficient can be calculated as a sum of backscattering cross sections per unit volume; (i) the scattering objects are uniformly distributed in the insonified ("sampled") part of the observation volume, with (j) random phases of the scattered echoes (i.e., random spacing of scattering objects, and movement of objects from one transmission to the next); (k) possible multiple-scattering effects and interaction between objects are neglected; and (l) excess attenuation from power extinction caused by volume scattering is neglected.Except for (d), these are all common assumptions underlying models used in fish abundance measurement [1][2][3][4][5][6][7][8][9][10][11][12][13][14][20][21][22][23][24][25][26][27][28][29][30].Assumptions (a) and (b) relate to the transducer and electric components of the echo sounder system.Assumption (g) is discussed elsewhere [51].Assumptions (h)-(l) are included in the set of assumptions used by Clay and Medwin [2,9] to derive the analogous smallamplitude "in-water" expressions for v s , accounting for acoustic pressures in the sea only.Relatively extensive discussions on the validity of (h)-(l) are given in other studies [2,9,11], which also summarize key literature in this field.
The discussion of the validity of these assumptions is an extensive and complex subject, beyond the scope of this article.One objective here was to clearly point out the assumptions on which the theory relies, and at which step in the derivation each of them is applied, as described in Sections 2 and 3.

Conclusions
The conventional and generic small-amplitude power budget equations for fish abundance estimation and species identification, Equations ( 1) and ( 2) [21,25], and further developments of that theory [20,[22][23][24], have been generalized to account for effects of finite-amplitude sound propagation in the transmitted (incident) sound field for the fundamental frequency component of the received signal.Expressions are derived for the backscattering cross section, bs  , and the volume backscattering coefficient, v s , in terms of power budget equations for single-target and volume backscattering, respectively, averaged over a cycle of a monochromatic wave.These expressions, given by Equations ( 50) and (56), constitute a generic functional relationship for abundance measurement under conditions of finite-amplitude incident sound, accounting for possible finite-amplitude effects in at-sea echo sounder calibration, in survey operation, or both.Arbitrary electrical termination, and the range of electrical and acoustical echo sounder parameters involved in calibration and oceanic surveying, are accounted for.For finite-amplitude signals, the expressions for bs  and v s derived here, Equations ( 50) and (56), are shown to be equivalent to expressions that were previously proposed [20,36,37] without sufficiently detailed and documented derivations.Consequently, the present analysis represents a derivation of, and theoretical fundament, also for those expressions.
The establishment of equations for bs  and v s with such capabilities enable the evaluation of errors caused by finite-amplitude effects in abundance estimation and species identification.Such error analysis can be based on measurements or calculations, or both.For actual transducers or echo sounders, the functions ( ) B   and  routinely provided by echo sounder manufacturers.Alternatively, they can be calculated using numerical models.These capabilities provide a fundament to avoid, or, if necessary, to compensate for such errors.This includes the future establishment of recommended upper limits for echo sounder transmit electrical power levels to obtain a controlled reduction of finite-amplitude errors in calibration and surveying.It also enables the development of correction factors for oceanic survey data (current or historic) already subject to finite-amplitude errors [20,36].The consequences of the theory and results presented here-e.g., in terms of evaluation of measurement errors, power recommendations, and correction factors, accounting for echo sounder calibration prior to survey operation and the survey operation itself-are to be addressed elsewhere.
The finite-amplitude effects of the type addressed here, and the associated errors in fishery research [15][16][17][18][19][20]36,37], originate from the properties of the fluid propagation medium (in this case, seawater) and apply to sonar and echo sounder technology in general, irrespective of manufacturer.
Although this article was motivated by challenges in fisheries acoustics and abundance estimation, and is primarily related to finite-amplitude effects, the equations presented here are quite general and are not limited to this field.Under the assumption of validity of the assumptions stated in Section 4.2, the theory and results also apply to the use of echo sounders and sonar for the measurement of acoustic single-target and volume backscattering more generally, for cases where finite-amplitude effects of the incident field are significant, as well as for conditions of small-amplitude sound propagation.
where 4 sph bs    is the spherical scattering cross-section [2,10] for the equivalent omnidirectional scatterer of the single target.

A.2. Volume backscattering
Similarly, for physical interpretation of the various terms in Equation (56), the power budget equation describing volume backscattering from the spherical shell subvolume p V , is conveniently made using the re-arrangement where have been used.First, an interpretation of the quantity , v n sph  is useful.The equivalent two-way beam solid angle, , represents the transducer's effective beam width (in terms of a solid angle) for transmission and reception combined, including finite-amplitude effects in the beam pattern of the incident wave.From the definition of a solid angle, the portion of the surface area of a sphere with radius at p r that is effectively insonified by the equivalent two-way beam solid angle equals
Three alternative and equivalent interpretations of     ,

G
, as discussed in the following, may provide useful insight, such as for physical interpretation of the various terms in the power budget equations for bs  and v s (cf.Appendix A).

G
may alternatively be interpreted as "the transducer's one-way electroacoustic conversion efficiency per unit solid angle in the     , direction, for lossless and small-amplitude sound propagation conditions in the fluid".This "directional efficiency" interpretation is used in Section 2.2.

G
is (a) a measure of how well the transducer converts input electrical power into acoustic waves headed in the     , direction, or (b) vice versa, by reciprocity, converts acoustic waves arriving from the     , direction into electrical power [22], where both (a) and (b) apply to lossless and small-amplitude sound propagation conditions in the fluid.
( , ) 4 for a pressure wave arriving in the ( , )   direction from a distant scattering object.For a scattered pressure wave arriving at normal incidence to the transducer, it thus follows that    1) and (2).(An alternative possible explanation for the missing factor may be that-by impedance matching, and without stating it-1 F   may possibly have been used for the Simrad EK500 echo sounder over the relevant narrow operational frequency band in question, cf. the discussion in Section 4.1.1 and [22,23].) In [28], the simplified expression is used, stated to be applicable for the special case . From the analysis of cf.Section 4.1.1,this is seen to be correct only for conditions at which 0 T X  Ω, i.e., at (or in the vicinity of) the series resonance frequency of the employed transducer vibration mode.In this case, 1 F   .It may be noted that this simplified expression can be used also for

Appendix D. Symbols and nomenclature
Table D1 summarizes the symbols, nomenclature, and units for the quantities used in the text.Bold-type symbols indicate a complex-valued quantity.

Symbol
wavenumber vectors in the fluid medium, where k e and ' k e are the (position dependent) unit vectors normal to the incident and scattered wavefronts, respectively.

B
for all pressure amplitudes and ranges, r.Away from the axis,

Figure 2 .
Figure 2. Sketch of an arbitrary transducer: target positions giving the relationship between ( , )   and ( ', ')   for backscattering from a single target at arbitrary location ( , , ) r   .The backscattering direction is given as '      , '      .(Reproduced from [22], with permission from the Institute sound pressure amplitude at the reference range ' 0 r from the target. ) gives the sound pressure amplitude in the fluid that is backscattered from a single target located at position ) , , (   r in the farfield, for the fundamental frequency component of the sound field, under conditions of finite-amplitude incident sound.The factor ) respectively.Here, T I is the input electric current amplitude delivered to the transducer during transmission.0 V is the output voltage amplitude across the transducer's electrical terminals at reception under open-circuit conditions.ax V M is the free-field open-circuit receiving voltage sensitivity for pressure waves incident along the acoustical axis (normally incident waves,

Figure 3 .
Figure 3. Sketch of the electrical connections for the electroacoustic transducer and electrical termination operating in (a) transmit and (b) receive modes.(Reproduced from [22], with permission from the Institute of Marine Research, 2020.)Let

G
represent the transducer's one-way electroacoustic conversion efficiency per unit solid angle, in the     ,

Figure 4 .
Figure 4. Helmholtz-Thevenin electric circuit for the electroacoustic transducer upon signal reception.(Reproduced from [22], with permission from the Institute of Marine Research, 2020.)

Figure 5 .
Figure 5. Sketch of the acoustic system under analysis, with an electroacoustic transducer operating as a transmitter and receiver of ultrasound, and acoustic volume backscattering from a multitude of scattering objects in a spherical shell observation volume, obs V .(Reproduced from [22], with permission from the Institute of Marine Research, 2020.) backscattering cross section for scatterers of type j in the unit volume V  .Consequently, backscattering cross section over all scatterer types in the unit volume V  .From Equations (51) and (52), it follows that

. 1 p r and 2 pr , cf. Figure 6 .
electrical power received due to volume backscattering from the unit volume dV .Here, the symbol st T  used in Equation (48) has been replaced by v T  , the average transmit electrical power used in oceanic (field) surveys (volume scattering).In practice, v s t T T    is often used, but for generality in the description, v T  is distinguished here from st T  .Now, assume that the observation volume obs V in the farfield is insonified using a pulse (e.g., a tone burst, denoted here as "ping") of time duration p  and angular carrier frequency  .The spatial extension of the pulse is Within the spherical shell volume obs V , the tone burst then covers a spherical shell sub-volume, p V (the "ping volume"), contained within ranges, say, The arrival times for the start and end of the tone burst are e., 0.75 m for a 1 ms pulse, often used in practice).

Figure 6 .
Figure 6.Sketch of the spherical shell sub-volume p V in the spherical shell observation volume, obs V in the observation volume obs V .An interpretation of Equation (56) in terms of power "flow" is given in Appendix A.2.Under the assumptions used here (cf.Section 4.2), v s is a property of the scatterers in p V , at range r, and is as such invariant to the sound pressure amplitude.The factor   accounts for finite-amplitude effects in the axial pressure and beam pattern of the incident field, and represents the deviation from the small-amplitude case, for v s . of finite-amplitude effects on the average received electrical power.
from the missing factor  F , and the absence of distinguishing between st R  and v s that include the factor  F .Apart from the fact that Pedersen did not distinguish between st R  and v R  , and between st T  and v T  , his Equations (2.26) (slightly rearranged) and (2.35) for small-amplitude sound propagation conditions are identical to Equations (60) and (61), respectively.
is the number of targets (fish) per unit volume.Apart from an apparent sign misprint in the absorption term (the exponent) of their Equation (3.15), these expressions for bs

F
conservation of power upon transmission in lossless media, cf.Appendix A.2. Electrical termination was not addressed, implicitly implying 1   as discussed above.
electrical power levels-be measured in controlled tank or laboratory experiments, in addition to the measurements of i( , )


is the transmitted electrical power, averaged over one cycle of the monochromatic wave, at the angular frequency  in question.Multiplying with     , G gives the acoustic power produced by a point source radiating an (omnidirectional) intensity that is equal to the transducer's radiated intensity in the     , direction at range r under lossless and small-amplitude sound propagation conditions in the fluid, cf.Appendix B (Equation (B.1) and Interpretation 2).Multiplying with yields the transducer's radiated intensity at the target position, conditions, and with absorption accounted for.


axial and beam pattern effects of the incident sound field at the target position.Multiplication with sph  gives the acoustic power scattered by the target, represented here by the target's equivalent omnidirectional scatterer.Multiplying with acoustic power density (i.e., the intensity) of the scattered field at the centre of the transducer front, with absorption accounted for.Multiplication with the "effective receiving area" (cf.Appendix C) of the receiving transducer, received electrical power at the transducer's electrical terminals for the particular electrical termination case 1 at the transducer's electrical terminals for an arbitrary electrical termination load.
scattering cross section of the effective sampled volume n s V in p V , where change in the beam width due to finite-amplitude sound propagation effects is accounted for [52].In Equation (A.2), v T  is the transmitted electrical power, averaged over one cycle of the monochromatic wave, at the angular frequency  in question.By following the reasoning used for interpretation of Equation (A.1), gives the transducer's radiated intensity in the axial direction, at the p V range, p r r  , under small-amplitude conditions, and with absorption accounted for (cf.Appendix B; Equation (B.1) and Interpretation 1).


finite-amplitude effects of the incident beam.Multiplication with the effective spherical scattering cross section of the sampled volume n s V , , v n sph  , gives the acoustic power scattered from p V , including finite-amplitude effects on the incident beam width.Multiplying with yields the free-field acoustic power density (i.e., the intensity) of the scattered field at the center of the transducer front.Multiplication with the "effective receiving area" of the receiving transducer (cf.Appendix C), for normally incident sound to the transducer ( at the transducer's electrical terminals for an arbitrary electrical termination load.Appendix B. Transducer gain,     ,

Interpretation 2 :. 2 . 3 :
radiated by the transducer in the     , direction and at range r, under small-amplitude conditions, and including the effects of absorption.From Equation (B.1),     , G may be interpreted as "the ratio of the intensity produced by the transducer in the  (point) source that is radiating the amount of electrical power T  being supplied to the transducer, both under lossless and small-amplitude sound propagation conditions in the fluid".This "in-fluid" intensity interpretation is used in Appendix A.1.(It may be noted that Interpretation 1 corresponds to the definition of     , G used in [45] for electromagnetic waves.)Alternatively, for the given direction     , , the transducer gain     , G may from Equation (B.1) be interpreted as "the ratio of the acoustic power produced by a point source that is radiating an (omnidirectional) intensity equal to the transducer's radiated intensity in the     , direction and at range r, transmitted electrical power, T  , under lossless and small-amplitude sound propagation conditions in the fluid".This "electroacoustic power conversion" interpretation-more closely related to the transducer's conversional efficiency,  -is used in Appendix AInterpretation Thirdly, since  and 0 D are both independent of     , , it follows from Equations (

A
is thus interpreted as the "effective receiving area" of the receiving acoustic transducer.From Equations (C.1) and (C.2), it follows that in electromagnetics and acoustics (e.g., for radar and sonar) the expressions for the "effective receiving area" of the antenna/transducer are different, by the factor F  .By revisiting the derivation of Equations (1)-(3) given in[21], in which the simplified expression used (valid for specific cases only, cf.below), this result may possibly explain the missing factor F  in Equations ( . Section 4.1.1.Use of the general expression, Equation (C.2), would, however, be more correct and accurate, and comparably simple.From Equations (C.2) and (41), it may be noted that relationship between the transducer's free-field open-circuit voltage receiving sensitivity, V M , and ac eff A .

Nomenclature Unit 0  0 c 0 r 0 rpppPPAIII 0 D 2 bsNmMI 0 VVVI
Ambient density of the fluid medium kg/m 3 Small-amplitude sound velocity of the fluid medium m/s  Sound pressure acoustic attenuation coefficient of the fluid medium Np/m  Sound pressure acoustic attenuation coefficient of the fluid medium dB/m f Frequency of the monochromatic wave = fundamental frequency of the finitePosition vector expressed in the spherical coordinate system of the echo sounder ' ( ', ', ') r r    Position vector expressed in the spherical coordinate system of the target r Range from the center of the transducer front surface m Axial reference range from the center of the transducer front surface (e.g., 1 m) m ' Reference range from a single scattering target (e.g., 1 m) m i Incident pressure wave radiated by the transducer, for small-amplitude sound propagation in the fluid Pa n i Incident pressure wave radiated by the transducer, for finite-amplitude sound propagation in the fluid, for the fundamental frequency of the wavefield Pa s Pressure wave scattered by a single target Pa i Amplitude of i p at axial reference range 0 r , extrapolated from the transducer's far-field Pa bs P Free-field sound pressure in the fluid, backscattered from a single target, in the position of the center of the transducer front Pa ,0 bs P Sound pressure amplitude backscattered from a single target, at reference range ' 0 r from the target, extrapolated from the target's far-field Pa i Amplitude constant for the incident scattered field Paintensity of the incident wave at the center position of a single target, for finite-amplitude sound propagation in the fluid, for the fundamental frequency of the wavefield W/m 2 s Intensity of the wave scattered from a single target, in the far-field of the target WIntensity of the wave scattered from a single target, at reference range ' Intensity of the wave backscattered from a single target, at reference range ' 0 r from the target, extrapolated from the target's far-field Wbeam pattern of the incident sound pressure wave = far-field beam pattern of the transducer in transmit and receive operations, for small-amplitude sound propagation in the fluid − ( , , ) n i r   B Far-field beam pattern of the incident sound pressure wave = far-field beam pattern of the transducer in transmit operation, for finite-amplitude sound propagation in the fluid, for the fundamental frequency of the wavefield − ( ', ') s   B Beam pattern of the sound pressure wave scattered from a single target, in the far-field of the target − ( , , ) n rel r   B Beam pattern finite-amplitude factor − Equivalent two-way beam solid angle of the transducer, for small-amplitude sound propagation in the fluid sr n Equivalent two-way beam solid angle of the transducer, for finite-amplitude sound propagation in the fluid sr n rel Beam solid angle finite-amplitude factor sr bs  Backscattering cross section of a single target (e.g., calibration sphere, or fish) m Expected value of the backscattering cross section of a single target (e.g., fish) Number of scattering objects of type j per unit volume j Number of scattering objects of type j in a unit volume j bs,  Backscattering cross section for a scattering object of type j m 2 bs  Backscattering cross section over all scatterer types, in a unit volume m 2 open-circuit voltage receiving sensitivity of the transducer V/Pa ax V Free-field open-circuit voltage receiving sensitivity of the transducer, for normally incident pressure waves V/Pa T Input electric current amplitude delivered to the transducer during transmission A Output voltage amplitude across the transducer's electrical terminals at reception, under open-circuit conditions V T Voltage amplitude across the transducer's electrical terminals at transmission V R Voltage amplitude across the transducer's electrical terminals upon reception, for single-target backscattering V R Current amplitude at the transducer's electrical terminals upon reception, for single-target backscattering A Input electrical impedance of the transducer in transmit operation, with resistance T R and reactance T X ) electrical impedance of the transducer in receive operation, with resistance R R and reactance R X Input electrical impedance of the receiving electronics network, with resistance E R and reactance E X Average electrical power delivered to the transducer at transmission ("average transmit electrical power"), for single-target ( Average electrical power delivered to the transducer at transmission ("average transmit electrical power"), for single-target backscattering operation (i.e., calibration sphere, individual fish) W v TAverage electrical power delivered to the transducer at transmission ("average transmit electrical power"), for volume scattering (field survey) operation W , st n R  Average electrical power delivered by the transducer to the receiving electronics ("average received electrical power"), for single-target backscattering operation W