Issue 
Acta Acust.
Volume 6, 2022



Article Number  14  
Number of page(s)  12  
Section  Underwater Sound  
DOI  https://doi.org/10.1051/aacus/2022010  
Published online  13 April 2022 
Scientific Article
Finiteamplitude sound propagation effects in volume backscattering measurements for fish abundance estimation
^{1}
University of Bergen, Dept. of Physics and Technology, P.O. Box 7803, 5020 Bergen, Norway
^{2}
Kongsberg Maritime AS, P.O. Box 111, Strandpromenaden 50, 3191 Horten, Norway
^{3}
Institute of Marine Research, P.O. Box 1870, Nordnes, 5817 Bergen, Norway
^{*} Corresponding author: audun.pedersen@uib.no
Received:
6
September
2021
Accepted:
15
March
2022
Scientific advice for governmental management of marine resources relies on acoustical observation methods. Quantification and identification of fish and plankton species are often achieved using multifrequency acoustic data. Accurate measurements of backscattering crosssections and volume backscattering coefficients are essential. Systematic errors from finiteamplitude sound propagation are demonstrated in experimental survey measurements on Atlantic mackerel using 120 kHz and 200 kHz echosounders and high power settings. Finiteamplitude signal distortion causes excess transmission loss that is not accounted for in fisheries acoustics today, other than by fixed limits on the maximum transmitted power. The demonstrated errors are of a magnitude that can seriously bias abundance estimation and species identification. It is shown how the finiteamplitude effects can be modelled and predicted quantitatively, within a framework of electroacoustic power budget equations. A method is provided to calculate related errors in echosounder calibration and oceanic measurement of acoustic volume backscattering. When accounting for finiteamplitude effects in echosounder signal processing, higher transmit powers can be used when needed to improve signaltonoise ratio or extend measurement range. The results indicate that historical survey data can be adjusted for such errors using numerical simulations. The echosounder characteristics relevant to finiteamplitude effects can be determined by laboratory measurements.
Key words: Fisheries acoustics / Finite amplitude sound / Marine acoustics / Nonlinear acoustics / Sonar
© The Author(s), Published by EDP Sciences, 2022
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
A common method to estimate abundance of fish in the oceans is based on measurement of acoustic backscatter from the water column [1, 2]. While 38 kHz is a main measurement frequency for quantification, simultaneous measurements in the range 18 kHz to 400 kHz are typically used to support species identification alongside trawl sampling. Multifrequency echosounder systems can operate at discrete frequencies and lately also across multiple continuous frequency intervals to improve range resolution and frequency response measurements. Some of the higher frequencies are also recommended for abundance estimation, e.g., 200 kHz for Atlantic mackerel [3, 4].
Quantitative measurements of acoustic backscattering are made by calibrated echosounders. Elastic spheres of materials such as copper (Cu) or tungsten carbide (WC) are suspended in the sound beam, and the echosounders’ measurements of their target strengths are adjusted against theoretical values (e.g. [5, 6]). Calibration data thus obtained for single targets on the sound beam axis are extended to the measurement of volume backscatter by applying a factorycalibrated beam pattern for each echosounder transducer [6, 7].
Linear sound propagation is assumed in fisheries acoustics and many other applications, with spherical spreading, sound absorption, and scattering as the only transmission loss mechanisms [8]. This approximation applies to smallamplitude signals. Above an applicationdependent amplitude threshold, nonlinear effects must be considered, and the sound is said to be of finite amplitude [9]. As higher frequencies are employed for species identification and increasingly efficient echosounder systems have been taken into use, the assumption of linear sound propagation is not always valid for transmit power levels accessible in fisheries acoustics. Signal distortion due to finiteamplitude sound propagation has been observed to cause additional transmission loss when the transmit power is increased for measurement at longer distances (several hundred meters). Investigations have demonstrated, first experimentally [10], and later by numerical simulations [10–14] and further experiments [13, 14], that finiteamplitude effects can be significant for echo sounders used for fish abundance estimation and species identification. For available signal transmission powers (up to a few kW) and frequencies, the predominant finiteamplitude effects are convection of the medium and dependence of sound speed upon pressure [9]. Amplitudedependent waveform distortion is observed as increased transmission loss at the measurement frequency [10] (herein referred to as “nonlinear loss”) and generation of harmonics that can interfere with simultaneous measurements at higher frequencies [10]. As multifrequency systems have evolved from operating at discrete frequencies to continuous frequency intervals, it has become more difficult to avoid such interference [15].
Nonlinear loss can lead to error in the evaluation of the backscattering crosssection σ_{bs} and the volume backscattering coefficient s_{v} used in fish abundance estimation and species identification. The waveform distortion caused by finiteamplitude sound propagation accumulates over propagation distance through regions where the sound pressure amplitude is sufficiently high, and depends upon transmit power, signal frequency, transducer directivity, medium parameters, range, and angle off the sound beam axis. The resulting nonlinear loss manifests itself as excess axial transmission loss and flattening of the beam pattern [11–13].
Advice on maximum power settings was given for some echosounder transducers and operating frequencies in 2002 [16] and 2008 [17], to reduce finiteamplitude effects to a presumably insignificant level. However, the recommendations differ and overly simplified extrapolations are used, based on nominal power per transducer surface area [17]. The transmit signal power is adjustable in many of today’s echosounder systems. A recent system also makes it possible to limit the maximum userselectable power setting for each connected transducer. The available recommendations do not indicate how much nonlinear loss may remain when the recommended power settings are used. Further work is needed to manage finiteamplitude effects at a quantitative level, including more accurate calculations specific to the relevant echosounder transducers. A method should be established formally, to quantify finiteamplitude effects so that related measurement errors can be avoided or compensated for. The present work does not replace the existing recommendations directly, but intends to improve the empirical and methodological basis for revisiting such recommendations and their form.
A set of field data, acquired in 2004, shows systematic differences in measured s_{v} in experimental acoustic survey measurements when alternating between “low” and “high” transmit powers. The data presented here are improved in volume and quality compared to what has been reported earlier [13]. Numerical simulations indicate that the observed systematic differences can be ascribed to finiteamplitude sound propagation in seawater. A scheme for modelbased correction for finiteamplitude sound propagation effects was proposed in [13], building on the conventionally used power budget equations (1) and (2). Formal derivations supporting the validity of these equations in acoustics were however not found in scientific literature at the time of the original study. This lack of foundation was an obstacle to publishing the experimental and numerical results and thus led to a series of works [18–21] to establish the required formal derivations. The present paper presents, firstly, a first comparison of experimental volume backscattering field data with numerical simulation results within this newly founded mathematical framework [21]. It is demonstrated that the framework established in refs. [18–21] is functional for describing the measurement of volume backscattering under conditions of significant finiteamplitude sound propagation. Secondly, a method to account for nonlinear loss in echosounder calibration is proposed within the same mathematical framework, to arrive at a complete model for quantifying and correcting for finiteamplitude sound propagation in fisheries research echosounders. A “finiteamplitude correction factor” is introduced to account for nonlinear loss upon both echosounder calibration and measurements of volume backscattering.
The paper is organized as follows: Section 2 describes the theoretical models used. This includes (i) the smallamplitude power budget equations conventionally used in fisheries acoustics, (ii) recent extensions of these power budget equations to finiteamplitude conditions, (iii) a further amendment to account for finiteamplitude effects during echosounder calibration and field operation, and (iv) theory for numerical simulation of finiteamplitude sound propagation in a fluid. The oceanic survey measurements and numerical simulations are described in Sections 3 and 4, respectively. Experimental and numerical results are presented and discussed in Sections 5 and 6, with conclusions given in Section 7. Boldface letters are used to indicate complexnumbered quantities.
2 Theory
2.1 Smallamplitude power budget equations (linear sound propagation)
A spherical coordinate system (r, θ, φ) is used with origin centred at the face of the echosounder transducer and the coordinate axis (r, 0, 0) coincident with the sound beam axis. θ ∈ [0, π] and φ ∈ [0, 2π] are the polar and azimuthal angles, respectively [18].
Under assumption of smallamplitude signal propagation, the electroacoustic power budget equations used conventionally in fish abundance estimation and species identification are [18, 22]
Here, σ_{bs} (units of m^{2}) is the backscattering crosssection of a single target in the far field of the transducer; s_{v} (units of m^{−1}) is the volume backscattering coefficient for a thin farfield spherical shell of thickness , centred at the transducer front; and s_{a} (dimensionless) is the area backscattering coefficient for a farfield interval of observation ranges [r_{1}, r_{2}].
In equation (1), r is the distance from the transducer to the single target. In equation (2), r is the midrange to the thin spherical shell. τ_{p} is the duration of the transmitted voltage signal (τ_{p} = 1 ms may often be used, corresponding to a spherical shell thickness ). α is the acoustic absorption coefficient in seawater, expressed in Np/m. λ = c_{0}/f is the acoustic wavelength at the measurement frequency f, and c_{0} is the smallsignal sound velocity of the fluid propagation medium (seawater). The electrical termination factor F_{Π} accounts for the electrical impedance load at the transducer’s terminals upon signal reception [13].
G(θ, φ) and G_{0} (both nondimensional) are the transducer gain and the axial transducer gain, given as [18]
respectively, where η (nondimensional) is the electroacoustic conversion efficiency of the transducer. B_{i}(θ, φ) = P_{i}(r, θ, φ)/P_{i}(r, 0, 0) is the beam pattern of the incident sound pressure field, P_{i} being the incident sound pressure amplitude assuming smallsignal (linear) sound propagation.
is the axial directivity factor [18, 23]. G(θ, φ) and G_{0} represent the transducer’s oneway electroacoustic conversion efficiency per unit solid angle, in the (θ, φ) and axial directions, respectively, for lossless sound propagation conditions in the fluid [21].
For reciprocal transducers, having the same beam pattern for transmission and reception, ψ (in steradians, sr) is the transducer’s equivalent twoway beam solid angle [18],
where dΩ ≡ sinθdθdφ is a solid angle element. ψ represents an effective beam width of the transducer’s intensity field, expressed in terms of a solid angle, accounting for the combined effect of transmission and reception [2, 13, 18, 21].
and are the average electrical transmit signal powers during transmission of singlefrequency tone bursts for singletarget and volume backscattering measurements, respectively. and are the corresponding average electrical powers received after twoway sound propagation, backscattering, and twoway electroacoustic transduction [21]. “Average” here refers to time averaging over one cycle of the monochromatic signal waveform with frequency f. In conventional operation is often used, and F_{Π} is commonly omitted from equations (1) and (2) [13, 18, 21, 22].
The assumptions underlying equations (1)–(3) are summarized in, e.g., [18–21], with references therein to relevant literature on the subject.
2.2 Finiteamplitude power budget equations (nonlinear sound propagation)
Since the amplitude of the backscattered wave is much weaker than that of the incident wave, finiteamplitude effects are considered most influential in the forwardpropagating (incident) sound beam while the backscattered wave is considered to propagate without nonlinear loss [21]. Amendments to equations (1) and (2) to account for finiteamplitude sound propagation effects in measurement of σ_{bs} and s_{v} were proposed in [13, 24] and more rigorously derived in [21]. Their form as given in [21] is
Three factors in equations (7) and (8) to account for finiteamplitude sound propagation effects. The “axial finite amplitude factor”, defined as [21]
accounts for nonlinear loss on the sound beam axis in the incident complex sound pressure amplitude due to finiteamplitude sound propagation, . The absolute value is in the interval 〈0, 1〉.
The “beam pattern finite amplitude factor”, defined as
accounts for finiteamplitude effects on the beam pattern (beam flattening) in the forwardpropagating (incident) wave. is the beam pattern of the forwardpropagating pressure wave under finiteamplitude sound propagation conditions. for smallamplitude (i.e., linear) sound propagation, and under finiteamplitude conditions. Since and are both unity on the sound beam axis, .
The “beam solid angle finiteamplitude factor” , defined as [21]
accounts for a change in the shape of the radiated sound beam, while the beampattern for reception of the backscattered signal remains unchanged. ψ^{n}(r) is the equivalent twoway beam solid angle under conditions of finiteamplitude propagation of the forwardpropagating (incident) sound.
The finiteamplitude factors , , and can be obtained from hydrophone measurements for relevant echosounder transducers [21].
2.3 Influence of finiteamplitude effects in echosounder calibration
Equations (7) and (8) can be amended to account for possible nonlinear loss when the echosounder is calibrated against a reference target with known backscattering crosssection σ_{bs,ref}. The reference target is positioned on the sound beam axis at a range r = r_{cal} in the far field of the transducer (see, e.g., [6, 22]). Calibration is commonly viewed as determining the axial transducer gain G_{0} [6] so that equation (1) is satisfied with σ_{bs} = σ_{bs,ref} and (r, θ, φ) = (r_{cal}, 0, 0). The value for G_{0} obtained by the calibration measurement is here denoted , given by
Subsequent measurements of singletarget backscattering are carried out using for G(θ, φ) in equation (1) or (7), and volume backscattering measurements are made using for G_{0} in equation (2) or (8).
If finiteamplitude sound propagation effects are significant upon calibration, the propagation loss becomes higher than accounted for by equation (12), leading to error in . The correct value for G_{0} in presence of nonlinear loss is found by applying equation (7) in place of equation (1) to the calibration measurement:
The subscript “cal” indicates that the axial finite amplitude factor is evaluated for the source level and propagation medium upon calibration.
Measurements based on a calibration for which linear sound propagation has been inaccurately assumed, thus inadvertently using a calibration result that is different from the actual transducer gain G_{0}, can be compensated for the aggregated effects of nonlinear loss by substituting equation (13) for G_{0} in equations (7) and (8), yielding
The transducer gain G^{meas} (θ, φ) in equation (14), based on the calibration measurement, and assuming linear sound propagation in both calibration and field operation, is
Equations (14) and (15) account for possible finiteamplitude effects upon both calibration and subsequent field measurements of σ_{bs} and s_{v}, when the calibration has been performed assuming linear sound propagation.
2.4 Finiteamplitude sound propagation modelling
The Khokhlov–Zabolotskaya–Kuznetsov (“KZK”) equation [9, 25, 26] is used here to model nonlinear sound propagation and thus obtain numerical values for , , and for use in equations (14) and (15). The KZK equation is parabolic and suitable for forwardmarching solution by finite differences to the second order in nonlinearity, with validity inside the paraxial region of a sound beam.
A Cartesian coordinate system (x, y, z) is defined with origin and zaxis coinciding with the origin and axis of the previously defined spherical coordinate system. The KZK equation can be written, for the sound pressure p = P − P_{0} [9],
P_{0} and P are the hydrostatic and total pressures, respectively, and t′ = t − z/c_{0} is the retarded time assuming constant velocity c_{0}. β is the coefficient of nonlinearity [9], and ρ_{0} is the ambient density of the fluid propagation medium. Classical absorption of a thermoviscous fluid is accounted for by the diffusivity , where the absorption coefficient α = α_{2} f^{2} and α_{2} is frequencyindependent [23]. α_{2} can be adjusted to specify correct absorption at the operating frequency. The chemical relaxation mechanisms in seawater [27] however make the frequencysquared power law inaccurate for the harmonic frequency components generated through finiteamplitude signal distortion.
Equation (17) is solved in the frequency domain using the “Bergen Code” finitedifference algorithm [29]. The frequency domain approach allows insertion of the Francois–Garrison seawater absorption coefficient [28] for each harmonic component. This yields somewhat improved agreement with the experimental data. 50 harmonics are included in the calculations.
The KZK equation and the Bergen Code have been tested extensively against experimental measurements over the recent 3–4 decades. For an overview of such work in published literature, see, e.g., [9, 13]. Favourable agreement between simulations and experiments has been reported in general, over a wide range of simulation parameters.
3 Measurements
Insitu acoustic measurements on fish schools were performed from the research vessel G.O. Sars in the North Sea [13]. A Simrad EK60 echo sounder system was operated for simultaneous measurements at 38 kHz, 120 kHz, and 200 kHz acoustic frequencies. The three transducers are designed with the same farfield beam width 2θ_{−3dB} = 7° (at main thickness resonance) and mounted close together on a protrudable instrument keel. Effects of air bubbles were minimized by extending the instrument keel 3 m below the hull to 8.5 m depth. The overlap between the 38 kHz sound beam and each of the 120 kHz and 200 kHz sound beams has been estimated to 85% at 30 m depth, assuming that the beams were parallel [17]. Signal transmission was simultaneous between the frequencies so that the ship’s motion did not increase the offset between the ensonified volumes. The echosounder system was calibrated before the survey cruise with reference targets at range r_{cal} = 22 m (about 30 m depth), using electrical power settings of 250 W for the 120 kHz echosounder and 120 W for the 200 kHz echosounder. These are in the following referred to as the “low” power settings, while 1000 W is referred to as the “high” power setting for both operating frequencies. The measurements at 38 kHz were made with constant power setting of 2000 W. Nonlinear loss at this reference frequency is thus constant throughout the measurement series, while factors such as fish density or school size cause similar variability in the measurements at all frequencies.
The vessel crossed repeatedly over an area with several schools while the transmit power for the 120 kHz and 200 kHz transducers were alternated between the “low” and “high” settings for each crossing (Fig. 1). The schools extended predominantly from 20 m to 110 m depth, with more than 90% of the backscattered energy from depths between 40 m and 100 m. The depth distribution was similar throughout the measurement series. There were essentially three schools (purple blobs in Fig. 1) in the area, passed in repeated straight cruise tracks (blue curves in Fig. 1). The schools were remarkably stationary during the 11 h from 13:30 to 00:30 UTC on 6–7 November 2004. The cruise tracks covered 42 nautical miles (nmi) and 70,500 pings over 11 h. The ship followed the schools for another 2.5 h as the schools started to move slowly. Data from this final period were not included in the analysis. A pelagic trawl haul was carried out at depth 70 m–90 m and indicated 95% Atlantic mackerel (Scomber scombrus L.).
Figure 1 Cruise track (blue curves) for the acoustic survey and observed fish schools (purple blobs). The axes indicate longitude and latitude in decimal degrees. The values in the legend are the nautical area scattering coefficient s_{A} [nmi^{−1}] [33]. Note that the cruise tracks do not overlap completely. Each school could be crossed straight over in one leg and on a side in the next, even if the schools did not move. 
Ambient noise was removed from the acoustic data using the method described in [30]. Acoustic data were analysed by means of Bergen Echo Integrator (BEI [31]) during the survey, and in a more thorough analysis immediately after the survey. Schools were encircled, species identified, and the result was stored to a database. The successor of BEI, Large Scale Survey System (LSSS [32]), was used to repeat the analysis and produce the results presented here. A total of 42 passes over schools were extracted from the acoustic survey data. The relative frequency response identified the species as Atlantic mackerel.
The measured s_{v}(r) were integrated with a threshold of S_{v} = 10 log_{10}(s_{v}/1 m^{−1}) = −82 dB re 1 m^{−1} [33] at all frequencies. The removal of ambient noise reduced the average backscatter s_{v} by 0.15%–0.25% depending on frequency. Each of the encircled schools contained 350–700 pings over 0.2 nmi–0.4 nmi path lengths.
The common variability between measurements at 38 kHz, 120 kHz, and 200 kHz cancels partially in the relative frequency response [34],
r_{a}(f) is thus a measurement result with less pointtopoint variability, but with the same dependence upon transmit power at f = 120 kHz and 200 kHz, as the respective s_{a}(f). The experimental results for s_{a}(f) are calculated according to equations (2) and (3).
4 Numerical calculations
The transmitting echosounder transducer is modelled as a circular planar piston, vibrating with uniform axial velocity in a rigid baffle of infinite extent. Within the parabolic approximation, the velocity amplitude U_{s} of the piston can be expressed via a source condition pressure amplitude [9, 29]. The source pressure amplitude and effective piston radius a were calculated from farfield hydrophone measurements with low amplitudes, using echosounder transducers of the same types as for the oceanic survey [13]. The source pressure amplitudes were scaled by measured voltage amplitudes across the electrical terminals of the transducers upon signal transmission [13].
Table 1 gives the parameters used in numerical simulations to compare with the insitu survey measurements. The propagation medium parameters of seawater [28, 35, 36, 37] are calculated based on CTD (conductivity, temperature, and pressure) sampling performed regularly during the survey cruise, averaged over depths from 9 m to 120 m. The measured water temperature was between 9.2 °C and 10.5 °C, and the sound speed was between 1490 m/s and 1492 m/s throughout the water column [35]. The calculated coefficient of nonlinearity [36] was between 3.58 and 3.59.
Parameters used in numerical calculations for comparison with the insitu survey measurements: Nominal power , source condition pressure amplitude P_{s}, effective source radius a, temperature T, salinity S, depth d, density ρ_{0}, pH, ambient sound speed c_{0}, and coefficient of nonlinearity β.
Numerical calculations were made to produce pairs of corresponding simulated sound pressure fields of finiteamplitude sound propagation (corresponding to the “low” and “high” experimental power settings) and where finiteamplitude effects were avoided by setting the source condition pressure amplitude to . From equation (9), simulated values for the axial finiteamplitude factor are calculated as
where the subscript “cal” indicates the sound propagation conditions for the case of echosounder calibration.
The beam pattern finiteamplitude factor is calculated from equation (10) and the definitions of and , yielding
The beam solid angle finiteamplitude factor is calculated by insertion of equation (21) into equation (11). Since the KZK equation is valid under the paraxial approximation and significant finiteamplitude effects are expected only in the main lobe, the integrals in equation (11) are evaluated for θ up to an integration limit θ_{lim} = 20°. Parameter variation shows that the reported results are unchanged with θ_{lim} down to 6°, while the first minimum of the beam pattern is at approximately θ = 8°. For rotational symmetry with respect to φ one has dΩ = 2π sinθ dθ, and the expression for the beam solid angle finiteamplitude factor becomes
Since the nonlinear loss is range dependent, finiteamplitude effects in measurement of r_{a}(f) should be accounted for by using equation (15) to calculate s_{v}(r), before integrating over a range interval [r_{1}, r_{2}] to obtain s_{a} according to equation (3). The experimental measurements of r_{a}(f) were however processed using standard methods as described in Section 3, and therefore no correction for nonlinear loss was made prior to integration.
For clarity in the following, the actual values of s_{v}, s_{a}, and r_{a} are identified with superscript “correct”, while the corresponding (erroneous) measurement results, assuming linear sound propagation, are identified with superscript “meas”. is thus the measurement result obtained using equation (2) (assuming linear sound propagation in the survey measurement), with used in place of (assuming linear sound propagation in echosounder calibration), given as
From equations (15) and (23) one can calculate a “finiteamplitude correction factor” for the total effect of nonlinear loss,
For comparison with the field data on the form of measured under conditions of alternating power settings, the quantity of interest is . By noting that and are both approximately constant at long ranges, for given frequency and source level [21], one can approximate by its average over the full observation range interval, or for a thinner layer, From equations (3) and (24), one finds
for and in the far field. Calibration was performed at one range , using the “low” power setting, so that a single value of applies to each measurement frequency. The average product for a range layer, on the other hand, depends on the power setting during field measurements, so that one can write
Since s_{a} (38 kHz) does not depend on the power settings for the 120 kHz and 200 kHz frequencies, equations (18), (26), and (27) can be combined to
5 Results
Figure 2 shows measurement results using the “low” and “high” power settings for the two measurement frequencies 120 kHz and 200 kHz (relative to 38 kHz), used in the survey measurements on Atlantic mackerel schools, over the observation depth interval [30 m, 110 m]. Results are shown on a logarithmic (dB) scale. The calculated average measurement results for each power setting, 〈r_{a}(f)〉, are shown using dashed and full lines for the “low” and “high” power settings, respectively. The ratios between 〈r_{a}(f)〉 measured using the two power settings, , are 0.75 (−1.2 dB) for the 120 kHz measurements, and 0.48 (−3.2 dB) for the measurements made at 200 kHz.
Figure 2 Measured relative frequency response, r_{a}(f), obtained from oceanic insitu acoustical survey measurements of s_{v} on Atlantic mackerel schools. Data from within an observation depth interval of 30 m to 110 m are included. Each point (x or o) is based on measurements from a whole school. (a) 120 kHz measurement frequency, using 250 W and 1000 W electrical power settings in survey operation (alternating), and 250 W during calibration. (b) 200 kHz measurement frequency, using 120 W and 1000 W electrical power settings in survey operation (alternating), and 120 W during calibration. 
Figure 3a shows numerical calculations of for the “low” and “high” power settings. This represents the effect of nonlinear loss on volume backscattering measurement (Eq. (8)). Results are shown on a logarithmic (dB) scale vs. logarithmic range. For the 120 kHz measurements with the “low” power setting, the nonlinear loss effect grows from about −0.1 dB towards −0.3 dB (factors of 0.98–0.94 in natural units) over the ranges 10 m–300 m. With the “high” power setting the nonlinear loss effect goes from −0.4 dB to −1.0 dB (factors of 0.90–0.78 in natural units) over the same ranges.
Figure 3 (a) Numerically calculated values for , representing the integrated nonlinear loss over the sound beam and thus the nonlinear loss in the s_{v} measurement. The parameters in Table 1 were used for the calculations. (b) Numerical and experimental values for the ratio between volume backscattering measurements using the “high” and “low” power settings. The experimental data are divided in depth bins with 20 m width, indicating mean values and standard deviations for the 11hour measurement series. 
For the 200 kHz operating frequency with the “low” power setting, the nonlinear loss effect goes from −0.2 dB at 10 m range to −0.5 dB at 300 m (factors of 0.95 and 0.89, respectively). The effect increases to −2.1 dB and −3.5 dB (factors of 0.61 and 0.44 in natural units) at the same ranges when the “high” power setting is used.
Figure 3b compares the numerically calculated as a function of depth, to the measurement results separated into 20 m depth bins. The measurement results are shown as mean values and standard deviations at the centre depth of each bin. The numerical and experimental data agree within the experimental uncertainty.
Table 2 compares the measurements and numerical results for the full observation range interval. The first two rows give the average values calculated numerically using equations (19) and (22), for a depth interval [40 m, 100 m], for the “high” and “low” power settings, respectively. The third row in Table 2 gives the simulated from equation (28) and the first two rows of the table. The corresponding insitu measured values , cf. Figure 2, with relative standard deviations, are given in the next two rows. The numbers of datasets exported from LSSS and processed for mean value and standard deviation are given in the last line. The simulated at “high” and “low” power settings agree with the measured to within 0.7 dB.
Summary of simulation and measurement results.
Numerically calculated values for the correction factor (Eq. (24)) are shown as an example in Figure 4. Results for the 1000 W power setting and the parameters in Table 1 are used for both calibration and subsequent measurements of s_{v}(r). Values for five ranges r_{cal} to the calibration reference target have been included to visualize the effect of nonlinear loss during calibration on the total measurement error. Due to the flattening of the main lobe, s_{v}(r) is always overestimated at the calibration range when nonlinear loss is significant but not accounted for. Further finiteamplitude sound propagation beyond the calibration range reduces the total measurement error until it may pass unity and the error from calibration is outweighed.
Figure 4 Numerically calculated example values for the finiteamplitude correction factor , with sound propagation parameters from Table 1 and 1000 W power setting for both calibration at range r_{cal} and s_{v} measurement at range r. The calibration range is indicated with a mark on each curve. Logarithmic scale, with corresponding ratios shown on the righthand ordinate axis. (a) 120 kHz, (b) 200 kHz. 
The same numerical results are used in Figure 5 to illustrate the combined result of axial nonlinear loss, flattening of the main lobe, and echosounder calibration. Dashed curves indicate the twoway intensity beam pattern B_{i}(θ, φ)^{4} for the 200 kHz echosounder without accounting for nonlinear loss, cf. equation (4). The fulldrawn curves show the same twoway beam pattern divided by the ratio
where is calculated using equation (1) and is from equation (14) using the 1000 W power setting. Echosounder calibration is here assumed to take place with r_{cal} = 14 m. The leftmost plot assumes linear sound propagation upon calibration, i.e., . The middle plot illustrates finiteamplitude effects during calibration, where the echosounder is calibrated to measure σ_{bs} correctly on the sound beam axis. The flattening of the main lobe leads to overestimated measurements elsewhere at the same range and thus for s_{v}. The rightmost plot illustrates the effect of nonlinear loss at 300 m range after the same calibration at r_{cal} = 14 m, visualizing how further nonlinear loss between ranges r_{cal} and r > r_{cal} can counteract the overestimate of s_{v} experienced at the calibration range.
Figure 5 Illustration of nonlinear loss in the sound beam, using the numerical results for 200 kHz, 1000 W. Left: Twoway intensity beam pattern at 14 m range, with and without correction for nonlinear loss in the forward propagating wave. Centre: Twoway intensity beam pattern at 14 m range, corrected for nonlinear loss upon both measurement and calibration at the same range. Right: Twoway intensity beam pattern at r = 300 m range, corrected for nonlinear loss upon measurement and upon calibration at r_{cal} = 14 m range. 
6 Discussion
Systematic and significant differences are observed between the insitu r_{a} measurements obtained using “low” and “high” available power settings, for both of the studied operating frequencies. Fair agreement, within 0.5 dB, was found with simulations of finiteamplitude propagation in the forwardpropagating sound beam. This is taken to indicate that the observed amplitude dependency in backscattering measurements is caused predominantly by finiteamplitude sound propagation effects in seawater. This result confirms earlier studies in which numerical data were compared to hydrophone measurements in freshwater [10, 13, 14] and seawater [13], and echosounder measurements on reference targets (calibration spheres) in seawater [13]. The systematic effects are of a magnitude that can cause significant errors in quantitative echosounder measurements at both 120 kHz and 200 kHz measurement frequencies. The 1000 W electrical power setting used here as “high” has been standard in the past, e.g., in the widely used echosounder Simrad EK500 with 120 kHz and 200 kHz transducers with 7° opening angles (2θ_{−3dB}). Figure 4 shows examples of finiteamplitude correction factors that could be used in such cases.
Improved quantitative agreement may be achieved by accounting more accurately for the seawater and operational parameters. Such parameters include the characteristics of the sound source, measurement frequency, source level, range, hydrostatic pressure (depth), and seawater parameters such as the parameter of nonlinearity and the absorption coefficient. The model uncertainty is also important. Factors of relevance in this context are discussed in the following.
Further field data for the same and other relevant operating frequencies would strengthen the experimental basis for quantifying measurement errors across echosounder systems and environmental conditions. Target strength measurements using calibration spheres, comparing results from multiple power settings, have shown fair agreement with simulations on and off the sound beam axis [13]. Absolute sound pressure measurements have been made using hydrophones in freshwater and seawater for distances up to 9 m–12 m [10, 13, 14]. Additional hydrophone measurements of finiteamplitude sound propagation in seawater at longer distances would be useful as a basis for further work if sufficient positioning and acoustic measurement accuracy can be achieved.
With respect to source parameters, scientific echosounders use advanced transducers designed with low sidelobe levels. Their beam patterns are different from that of a uniform circular piston with the same −3 dB beam width, in particular at large angles off the sound beam axis. The paraxial approximation also affects the validity of the numerical calculations at large angles. The importance of this discrepancy is limited since significant nonlinear sound propagation occurs predominantly in the central main lobe in the present application. Some work has been done on the use of alternative source conditions for the simulations [13], and further study can still be useful. Considering the total body of experimental and numerical uncertainties, a possible outcome of such work may however be that the generality of a simple piston source condition outweighs the achievable benefit of improved agreement with specific transducers.
For use of the numerical solution algorithm [29] in the context of the present work, calculations of the factors , , and are based mainly on the pressure amplitude ratio . Their uncertainties are expected to be smaller than for simulations of absolute sound pressure amplitudes, since some systematic model uncertainties cancel each other, partially or completely, in the ratios between the numerical results.
The sensitivity of the numerical results to algorithm configuration parameters, source condition, and propagation medium parameters has been studied in [13]. Further sensitivity study on seawater parameters, including consideration of their spatial and temporal variability, could be a useful contribution. Air bubbles in the upper ocean layer could also be investigated as a possible cause of variability in relation to finiteamplitude effects. Such studies would help to determine how accurately one can estimate the effects of finiteamplitude sound propagation on fishery acoustic survey data, and thus to what extent modelbased compensation would reduce their measurement uncertainty with respect to backscattering cross sections and volume and area backscattering coefficients. Frequencydependent absorption coefficients for seawater were inserted for each harmonic component in place of the frequencysquared power law inherent in the KZK equation being solved. Possible effects of this adhoc adaptation could be investigated, perhaps by comparing to a model derived with relaxation terms.
7 Conclusions
Oceanic survey measurements on Atlantic mackerel have been carried out to investigate possible influence of acoustic signal amplitude on parameters used in fish abundance estimation and species identification. For the relative frequency response, r_{a}, 14% and 47% deviations have been measured at 120 kHz and 200 kHz, respectively, by alternating between “low” and “high” electrical transmit power settings. The “low” power settings were 250 W and 120 W for the 120 kHz and 200 kHz measurement frequencies, respectively. The “high” power settings were 1000 W for both frequencies. Echosounder calibration was done with the “low” power setting for both frequencies.
The simulations strongly indicate that the observed systematic deviations can be ascribed to finiteamplitude sound propagation effects in seawater. Simulations and measurements were compared for the ratio between r_{a}(f) obtained using high and low power settings. The agreement was found to be mainly within the experimental uncertainty. The formal framework established through [18–21] is demonstrated to be functional for describing the measurement of volume backscattering under conditions of significant finiteamplitude sound propagation. A method to account for nonlinear loss in echosounder calibration is proposed within the same framework. The result is a complete method for quantifying and correcting for finiteamplitude sound propagation in fisheries research echosounders or establishing power limits with known and acceptable nonlinear loss for an application at hand. Potentials for further improved agreement are discussed.
If not accounted for, the demonstrated effects can lead to significant bias in biomass abundance estimation and acoustic species identification. Finiteamplitude effects are not accounted for in the conventional methods used today within fishery acoustics, but minimized approximately by recommended limits on maximum transmitted power [16, 17], partially based on overly simplified extrapolations [17]. The effects are however well understood and can, as shown here, be accounted for using available methods.
The results provide an indication of typical errors that can be encountered in fisheries acoustic measurements when high transmit powers are used and finiteamplitude sound propagation effects are not accounted for. The approach demonstrated here may serve as a basis to refine previous recommendations and add quantitative information on nonlinear loss. Nonlinear loss and generation of harmonics can be significant even at relatively low source levels and frequencies. A wellchosen transmitted power thus depends on transducer directivity, efficiency, frequency, and on medium parameters, and should be rooted in a quantitative requirement on measurement uncertainty. Further parametric study could provide a fuller map of measurement uncertainty contributions and potential corrections for finiteamplitude effects in both historical and current acoustic survey data. If finiteamplitude effects are accounted for in echosounder signal processing, higher powers may be used when needed to improve signaltonoise ratio or extend measurement range in echosounder and sonar measurements. The acoustic quantities accounting for finiteamplitude sound propagation effects, , , and , can be calculated for relevant field conditions, and measured for relevant echosounder transducers in a laboratory tank setup.
Conflict of interest
The authors declared no conflict of interest.
Acknowledgments
The work was supported by The Research Council of Norway (MARE 152790) and Norwegian Seafood Research Fund, Norway (FHF 222042). Discussions with M. Vestrheim, University of Bergen, and A.C. Baker, Christian Michelsen Research AS, are acknowledged. The associate editor and two reviewers are thanked for their valuable feedback.
Data availability statement
The research data associated with this article is available in Zenodo (https://doi.org/10.5281/zenodo.6059530 [38]). The data include measurement data exported to ASCII by LSSS and simulation results in a MATLAB .mat file. Accompanying MATLAB scripts produce Figures 2–5 and the results in Table 2 and include guidance to the data formats.
References
 O. Sund: Echo sounding in fishery research. Nature 135, 3423 (1935) 953. [CrossRef] [Google Scholar]
 C.S. Clay, H. Medwin: Acoustical Oceanography: Principles and Applications, John Wiley & Sons. 1977. [Google Scholar]
 ICES: Report of the Planning Group on Aerial and Acoustic Surveys for Mackerel, 57 April 2005, ICES CM 2005/G:13: Bergen, Norway. 2005. [Google Scholar]
 N. Gorska, R.J. Korneliussen, E. Ona: Acoustic backscatter by schools of adult Atlantic mackerel. ICES Journal of Marine Science 64 (2007) 1145–1151. [CrossRef] [Google Scholar]
 K.G. Foote, D.N. MacLennan: Comparison of copper and tungsten carbide calibration spheres. Journal of the Acoustical Society of America 75, 2 (1984) 612–616. [CrossRef] [Google Scholar]
 D.A. Demer, L. Berger, M. Bernasconi, E. Bethke, K. Boswell, D. Chu, R. Domokos, et al.: Calibration of acoustic instruments. ICES Cooperative Research Report No. 326. 2015. [Google Scholar]
 H. Bodholt: The effect of water temperature and salinity on echo sounder measurements, in ICES Symposium on Acoustics in Fisheries and Aquatic Ecology, Montpellier, June, 2002. [Google Scholar]
 J. Simmonds, D.N. MacLennan: Fisheries Acoustics. Theory and Practice, 2nd ed., Blackwell Science Ltd. 2005. [CrossRef] [Google Scholar]
 M.F. Hamilton, D.T. Blackstock (Eds.): Nonlinear Acoustics, Academic Press. 1998. [Google Scholar]
 F.E. Tichy, H. Solli, H. Klaveness: Nonlinear effects in a 200kHz sound beam and consequences for target strength measurement. ICES Journal of Marine Science 60 (2003) 571–574. [CrossRef] [Google Scholar]
 A.C. Baker, P. Lunde: Nonlinear effects in sound propagation from echosounders used in fish abundance estimation. Numerical simulation results. CMR Tech. Note CMRTN02A10008Rev01, Bergen, Norway: Christian Michelsen Research AS. 2002. Declassified revision of CMR Tech. Note CMRTN02F10008, 2002 (confidential). [Google Scholar]
 A.C. Baker, P. Lunde: Nonlinear propagation from circular echosounder transducers. Numerical simulation results. CMR Tech. Note CMRTN01A10010Rev01, Bergen, Norway: Christian Michelsen Research AS. 2011. Declassified revision of CMR Tech. Note CMRTN01F10010, 2001 (confidential).. [Google Scholar]
 A. Pedersen: Effects of nonlinear sound propagation in fisheries research. Ph.D. Thesis, University of Bergen, 2006 [Online]. https://hdl.handle.net/1956/2158. [Google Scholar]
 A. Pedersen, M. Vestrheim, P. Lunde: Nonlinear sound propagation effects in fisheries research echo sounders  Measurements and simulations in fresh water, in Proc. 28th Scandinavian Symposium on Physical Acoustics, Ustaoset, Norway, 2005. [Google Scholar]
 B. Khodabandeloo, E. Ona, G.J. Macaulay, R. Korneliussen: Nonlinear crosstalk in broadband multichannel echosounders. Journal of the Acoustical Society of America 149, 1 (2021) 87–101. [CrossRef] [PubMed] [Google Scholar]
 Kongsberg Maritime: Nonlinear effects. Recommendation for fishery research investigations. [Online]. https://www.kongsberg.com/no/maritime/products/commercialfisheries/ekkoloddsvingere/nonlineareffects/. [Accessed 30 November 2021], 2002. [Google Scholar]
 R.J. Korneliussen, N. Diner, E. Ona, L. Berger, P.G. Fernandes: Proposals for the collection of multifrequency acoustic data. ICES Journal of Marine Science 65, 6 (2008) 982–994. [CrossRef] [Google Scholar]
 P. Lunde, A.O. Pedersen, R.J. Korneliussen, F.E. Tichy, H. Nes: Powerbudget and echointegrator equations for fish abundance estimation. Fisken og havet 10 (2013) https://www.hi.no/hi/nettrapporter/fiskenoghavet/2013/fh_102013_echo_integrator. [Google Scholar]
 P. Lunde, R.J. Korneliussen: A unifying theory explaining different power budget formulations used in modern scientific echosounders for fish abundance estimation. Fisken og havet 7 (2014). https://www.hi.no/hi/nettrapporter/fiskenoghavet/2014/lunde_korneliussen_fh_072014_20150214_mergeddocument. [Google Scholar]
 P. Lunde, R.J. Korneliussen: Powerbudget equations and calibration factors for fish abundance estimation using scientific echo sounder and sonar systems. Journal of Marine Science and Engineering 4 (2016) 43. [CrossRef] [Google Scholar]
 P. Lunde: Finiteamplitude power budget equations for acoustic fish abundance estimation. Journal of Marine Science and Engineering 8 (2020) 98. [CrossRef] [Google Scholar]
 Simrad AS: Operator Manual: SIMRAD EK500 Fishery research echo sounder. Scientific echo sounder: Base version. Doc. no. P2170/Rev. G, Horten. 1997. [Google Scholar]
 L.E. Kinsler, A.R. Frey, A.B. Coppens, J.V. Sanders: Fundamentals of Acoustics. 4th ed., John Wiley & Sons. 2000. [Google Scholar]
 P. Lunde, A.O. Pedersen: Sonar and power budget equations for backscattering of finite amplitude sound waves, with implications in fishery acoustics for abundance estimation of marine resources, in 35th Scand. Symp. Phys. Acoust., Geilo, Norway, Jan 29Feb 1, 2012. [Google Scholar]
 E.A. Zabolotskaya, R.V. Khokhlov: Quasiplane waves in the nonlinear acoustics of confined beams. Soviet Physics Acoustics 15, 1 (1969) 35–40. [Google Scholar]
 V.P. Kuznetsov: Equations of nonlinear acoustics. Soviet Physics Acoustics 16, 4 (1971) 467–470. [Google Scholar]
 R.E. Francois, G.R. Garrison: Sound absorption based on ocean measurements. Part I: Pure water and magnesium sulfate contributions. Journal of the Acoustical Society of America 72, 3 (1982) 896–907. [CrossRef] [Google Scholar]
 R.E. Francois, G.R. Garrison: Sound absorption based on ocean measurements. Part II: Boric acid contribution and equation for total absorption. Journal of the Acoustical Society of America 72, 6 (1982) 1879–1890. [CrossRef] [Google Scholar]
 J. Berntsen, J.N. Tjøtta, S. Tjøtta: Interaction of sound waves. Part IV: Scattering of sound by sound. Journal of the Acoustical Society of America 86, 5 (1989) 1968–1983. [CrossRef] [Google Scholar]
 R.J. Korneliussen: Measurement and removal of echo integration noise. ICES Journal of Marine Science 57 (2000) 1204–1217. [CrossRef] [Google Scholar]
 R.J. Korneliussen: The Bergen Echo Integrator postprocessing system, with focus on recent improvements. Fisheries Research 68 (2004) 159–169. [CrossRef] [Google Scholar]
 R.J. Korneliussen, Y. Heggelund, G.J. Macaulay, D. Patel, E. Johnsen, I.K. Eliassen: Acoustic identification of marine species using a feature library. Methods in Oceanography 17 (2016) 187–205. [CrossRef] [Google Scholar]
 D.N. MacLennan, P.G. Fernandes, J. Dalen: A consistent approach to definitions and symbols in fisheries acoustics. ICES Journal of Marine Science 59 (2002) 365–369. [CrossRef] [Google Scholar]
 R.J. Korneliussen, E. Ona: An operational system for processing and visualizing multifrequency acoustic data. ICES Journal of Marine Science 59 (2002) 293–313. [CrossRef] [Google Scholar]
 C.T. Chen, F.J. Millero: Speed of sound in seawater at high pressures. Journal of the Acoustical Society of America 62, 5 (1977) 1129–1135. [CrossRef] [Google Scholar]
 F.D. Cotaras, C.L. Morfey: Polynomial expressions for the coefficient of nonlinearity β and β/(ρc^{5})^{1⁄2} for fresh water and seawater, Journal of the Acoustical Society of America 94, 1 (1993) 585–588. [CrossRef] [Google Scholar]
 N.P. Fofonoff, R.C. Millard Jr: Algorithms for computation of fundamental properties of seawater. UNESCO Technical Papers in Marine Science 44 (1983). [Google Scholar]
 A.O. Pedersen, P. Lunde, F.E. Tichy, R.J. Korneliussen: Finite amplitude sound propagation effects in volume backscattering measurements for fish abundance estimation (3.0) [Data set]. Zenodo, 2022. https://doi.org/10.5281/zenodo.6059530. [Google Scholar]
Cite this article as: Pedersen AO. Lunde P. Tichy FE. & Korneliussen RJ. 2022. Finiteamplitude sound propagation effects in volume backscattering measurements for fish abundance estimation. Acta Acustica, 6, 14.
All Tables
Parameters used in numerical calculations for comparison with the insitu survey measurements: Nominal power , source condition pressure amplitude P_{s}, effective source radius a, temperature T, salinity S, depth d, density ρ_{0}, pH, ambient sound speed c_{0}, and coefficient of nonlinearity β.
All Figures
Figure 1 Cruise track (blue curves) for the acoustic survey and observed fish schools (purple blobs). The axes indicate longitude and latitude in decimal degrees. The values in the legend are the nautical area scattering coefficient s_{A} [nmi^{−1}] [33]. Note that the cruise tracks do not overlap completely. Each school could be crossed straight over in one leg and on a side in the next, even if the schools did not move. 

In the text 
Figure 2 Measured relative frequency response, r_{a}(f), obtained from oceanic insitu acoustical survey measurements of s_{v} on Atlantic mackerel schools. Data from within an observation depth interval of 30 m to 110 m are included. Each point (x or o) is based on measurements from a whole school. (a) 120 kHz measurement frequency, using 250 W and 1000 W electrical power settings in survey operation (alternating), and 250 W during calibration. (b) 200 kHz measurement frequency, using 120 W and 1000 W electrical power settings in survey operation (alternating), and 120 W during calibration. 

In the text 
Figure 3 (a) Numerically calculated values for , representing the integrated nonlinear loss over the sound beam and thus the nonlinear loss in the s_{v} measurement. The parameters in Table 1 were used for the calculations. (b) Numerical and experimental values for the ratio between volume backscattering measurements using the “high” and “low” power settings. The experimental data are divided in depth bins with 20 m width, indicating mean values and standard deviations for the 11hour measurement series. 

In the text 
Figure 4 Numerically calculated example values for the finiteamplitude correction factor , with sound propagation parameters from Table 1 and 1000 W power setting for both calibration at range r_{cal} and s_{v} measurement at range r. The calibration range is indicated with a mark on each curve. Logarithmic scale, with corresponding ratios shown on the righthand ordinate axis. (a) 120 kHz, (b) 200 kHz. 

In the text 
Figure 5 Illustration of nonlinear loss in the sound beam, using the numerical results for 200 kHz, 1000 W. Left: Twoway intensity beam pattern at 14 m range, with and without correction for nonlinear loss in the forward propagating wave. Centre: Twoway intensity beam pattern at 14 m range, corrected for nonlinear loss upon both measurement and calibration at the same range. Right: Twoway intensity beam pattern at r = 300 m range, corrected for nonlinear loss upon measurement and upon calibration at r_{cal} = 14 m range. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.