RESEARCH ARTICLE

A study on the multi-frequency acoustic target strength of krill using a stochastic distorted-wave born approximation (SDWBA) model

Wuju Son1https://orcid.org/0000-0001-5837-7933, Wooseok Oh2https://orcid.org/0000-0002-1578-6590, Hyoung Sul La1,3,*https://orcid.org/0000-0003-1285-580X, Kyounghoon Lee4,*https://orcid.org/0000-0002-4226-9011
Author Information & Copyright
1Division of Ocean & Atmosphere Sciences, Korea Polar Research Institute, Incheon 21990, Korea
2Institute of Low-Carbon Marine Production Technology, Pukyong National University, Busan 48513, Korea
3Department of Polar Sciences, University of Science and Technology, Daejeon 34113, Korea
4Division of Marine Production System Management, Pukyong National University, Busan 48513, Korea
*Corresponding author: Hyoung Sul La, Division of Ocean & Atmosphere Sciences, Korea Polar Research Institute, Incheon 21990, Korea, Tel: +82-32-760-5810, E-mail:hsla@kopri.re.kr
*Corresponding author: Kyounghoon Lee, Division of Marine Production System Management, Pukyong National University, Busan 48513, Korea, Tel: +82-51-629-5889, Fax: +82-51-629-5886, E-mail:klee71@pknu.ac.kr

Copyright © 2024 The Korean Society of Fisheries and Aquatic Science. This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Oct 14, 2023; Revised: Dec 20, 2023; Accepted: Dec 26, 2023

Published Online: Apr 30, 2024

Abstract

We examined the dB difference in target strength at multiple frequencies (ΔTS) for the identification of Antarctic krill (Euphausia superba) and ice krill (Euphausia crystallorophias) using a stochastic distorted-wave Born approximation model. Our investigation focused on ΔTS patterns at multiple frequencies in relation to size, along with key acoustic properties influencing TS, including density and sound speed contrast, fatness, and orientation. The findings revealed that the orientation and fatness significantly affect the ΔTS patterns. The results provide insight into the importance of the multi-frequency technique for estimating krill biomass and their ecological interactions with environmental features in the Southern Ocean.

Keywords: Antarctic krill; Ice krill; Multi-frequency acoustics; Target strength; Stochastic distorted-wave Born approximation (SDWBA) model

Introduction

Krill, crustacean zooplankton members of the Euphausiacea order, are central to linking primary producers with mid- to upper-trophic levels in Antarctic marine ecosystems (Everson, 2000). Antarctic krill (Euphausia superba) serve as a pivotal indicator species in these ecosystems and are significant fishery resources (Everson, 2000). The Commission for the Conservation of Antarctic Marine Living Resources continuously monitors Antarctic krill biomass for both conservation and resource utilization (Watkins et al., 2004). On the coast of the Southern Ocean, ice krill (Euphausia crystallorophias) predominate, substituting Antarctic krill and supporting neritic marine ecosystems (La et al., 2015a). Compared to Antarctic krill, research on ice krill biomass and spatial distribution is less extensive, largely due to the difficulty of accessing their high-latitude coastal habitats, often obstructed by sea ice. Hence, understanding the spatial distributions of both krill species are crucial for understanding Antarctic marine ecosystems across the open and coastal waters of the Southern Ocean.

Hydroacoustic systems, complemented by net sampling, have been widely used in examining the spatiotemporal distribution of marine organisms across wide areas (Lawson et al., 2008). Because marine organisms exhibit different acoustic scattering properties depending on species shape, size, orientation, and biological characteristics, it is possible to identify acoustic signals from fish and zooplankton using target strength (TS, in dB re 1 m2) variation of the interested organisms (Kang et al., 2002). The stochastic distorted-wave Born approximation (SDWBA) model is a standard model for understanding the TS of krill (SC-CAMLR, 2005), which assumes the krill body shape as a discrete cylinder and stochastically calculates the scattering properties according to the density contrast (g), sound speed contrast (h), fatness, and variations in orientation (Calise & Skaret, 2011). While TS research on krill using the SDWBA model has predominantly focused on Antarctic krill (Calise & Skaret, 2011; Conti & Demer, 2006), the TS characteristics of ice krill remain less explored (La et al., 2015b).

Given the distinct characteristics of ice krill and Antarctic krill in terms of prey preferences, developmental stages, and size distribution (Sala et al., 2002), these species may differ not only in their density and sound speed contrasts (Chu & Wiebe, 2005), but also in their size, fatness, and orientations. Such variations may result in distinct ΔTS patterns, which directly affect for the acoustic identification of these species. Therefore, we analyzed the ΔTS characteristics of both krill species employing the SDWBA model.

Materials and Methods

The Stochastic distorted-wave Born approximation (SDWBA) model was used to examine ΔTS variations for key variables (length, shape, density and sound speed contrasts, fatness, and orientation distributions) that affect the TS of Antarctic krill and ice krill. Of the key variables used in the SDWBA model, the same values were used as in Table 1 of Calise & Skaret (2011), except for density and sound speed contrasts, fatness, frequencies ranges, and orientation. The density contrast used in the SDWBA model in this study are based on the measured values of 1.0000 and 1.0279 (Chu & Wiebe, 2005) and 1.0357 (Foote et al., 1990) for Antarctic krill and ice krill, respectively, and the sound speed contrasts of 1.0279 (Foote, 1990) and 1.0180 and 1.0480 (Chu & Wiebe, 2005). The fatness was simulated in increments of 0%, 10%, 20%, 30%, and 40%, and for the distribution of krill orientations, the literature values of N (−20°, 28°), N (−8°, 9°), N (0°, 27°), and N (11°, 4°) for Antarctic and ice krill were used (Conti & Demer, 2006; La et al., 2015a; Lawson et al., 2006; SC-CAMLR, 2010). For the SDWBA model simulation, krill lengths were selected in 2 mm intervals ranging from 4 mm to 60 mm, encompassing the length distributions of both Antarctic krill and ice krill. The differentiation between the two krill species was made based on the criteria found in the literature (Sala et al., 2002), categorizing individuals less than 35 mm as ice krill and those greater as Antarctic krill. The variations in ΔTS were analyzed for lengths of 10, 20, 30, 40, 50, and 60 mm. Additionally, the study investigated the variation in multiple frequency differences, including ΔTS120–38 kHz, ΔTS120–70 kHz, ΔTS200–38 kHz, ΔTS200–70 kHz, and ΔTS200–120 kHz, at commonly used frequencies for krill acoustic identification: 38, 70, 120, and 200 kHz.

Results and Discussion

Multi-frequency differences due to changes in density and sound speed contrasts

The effects of density and sound speed contrasts parameters on ΔTS at multiple frequencies were examined. For this purpose, the density and sound speed contrasts for Antarctic krill were set at observed values (g = 1.0279, h = 1.0480; g = 1.0357, h = 1.0279) and for ice krill at (g = 1.0000, h = 1.0180) (Chu & Wiebe, 2005; Foote et al., 1990). Krill fatness was fixed at 40%, and the orientation distribution was established at N (−20°, 28°) (SC-CAMLR, 2010). The multi-frequency difference in density and sound speed contrasts showed a linear decrease with increasing length for Δ TS70–38 kHz, Δ TS120–38 kHz, and Δ TS120–70 kHz (Fig. 1A–1C). For ΔTS200–38 kHz and ΔTS200–70 kHz, the difference declined with lengths of 4–40 mm, escalated with lengths of 40–50 mm, and fell again for lengths > 50 mm (Fig. 1D–1E). At ΔTS200–120 kHz, we observed changes in the multi-frequency difference based on the length of 36 mm (Fig. 1F). Variations in density and sound speed contrasts can significantly influence acoustic scattering characterization in zooplankton through acoustic scattering models, TS calculations, and biomass estimates (Stanton & Chu, 2000). Seasonal changes, such as increased lipid composition in zooplankton during winter for energy storage (Campbell & Dower, 2003), can alter density contrast, and fluctuating water temperatures can modify sound speed contrast (Forman & Warren, 2010). Our study revealed that TS patterns when density and sound speed contrasts were set at 1.0357 and 1.0480, respectively, differed markedly from other cases, underscoring these contrasts as key factors in krill acoustic signal identification. Therefore, accurately determining the density and sound speed contrasts for the target species is critical.

fas-27-4-225-g1
Fig. 1. ΔTS patterns (ΔTS70–38 kHz, ΔTS120–38 kHz, ΔTS120–70 kHz, ΔTS200–38 kHz, ΔTS200–70 kHz, and ΔTS200–120 kHz) versus length considering g and h. The red circle indicates g and h of ice krill measured by Chu & Wiebe (2005), the blue square indicates g and h of Antarctic krill measured by Chu & Wiebe (2005), and the black triangle indicates g and h of Antarctic krill measured by Foote et al. (1990) and Foote (1990). The distribution of orientation is fixed at N (−20°, 28°). Krill measuring less than 35 mm are categorized as ice krill, while those exceeding 35 mm are classified as Antarctic krill, following the criteria established by Sala et al. (2002). TS, target strength.
Download Original Figure
Multi-frequency differences based on fatness

The effects of fatness parameter on ΔTS at multiple frequencies were confirmed. For this analysis, krill fatness was set at 40%, with density and sound speed contrasts at 1.0357 and 1.0279, respectively, and orientation at N (−20°, 28°) to determine the multi-frequency differences (Fig. 2). We observed that multi-frequency differences in krill fatness reduced with an increase in length for ΔTS70–38 kHz, ΔTS120–38 kHz, and ΔTS120–70 kHz (Fig. 2A–2C). Additionally, the variation in multi-frequency differences amplified with longer lengths and higher fatness. The pattern for ΔTS200–38 kHz and ΔTS200–70 kHz showed a decline with increasing length, an increase between lengths of 40 and 56 mm, and then a subsequent decrease (Fig. 2D–2E). At the same krill length (excluding lengths < 12 mm), the impact of fatness on the multi-frequency difference was least pronounced at 50 mm and 56 mm for ΔTS200-38 kHz and ΔTS200–70 kHz, respectively.

fas-27-4-225-g2
Fig. 2. ΔTS patterns (ΔTS70−38 kHz, ΔTS120−38 kHz, ΔTS120–70 kHz, ΔTS200−38 kHz, ΔTS200–70 kHz, and ΔTS200–120 kHz) versus length considering fatness. The red circle, blue square, black triangle, green inverted triangle, and magenta diamonds represent 0%, 10%, 20%, 30%, 40% fatness, respectively. The distribution of orientation is fixed at N (−20°, 28°). Krill measuring less than 35 mm are categorized as ice krill, while those exceeding 35 mm are classified as Antarctic krill, following the criteria established by Sala et al. (2002). TS, target strength.
Download Original Figure

The influence of fatness was most evident at ΔTS200–120 kHz (Fig. 2F). Here, the pattern of ΔTS at multiple frequencies, typically decreasing with increasing krill length, shifted at fatness levels of 0%, 10%, 20%, 30%, and 40% for krill lengths of 52, 48, 44, 40, and 36 mm, respectively. In essence, higher fatness altered the ΔTS patterns at comparatively shorter krill lengths. TS variations due to fatness were more significant at larger lengths than at smaller ones. Given that the length distribution of summer Antarctic krill is longer than that of ice krill (Sala et al., 2002), this implies that fatness differences play a more vital role in identifying the acoustic signals of Antarctic krill compared to those of ice krill.

Multi-frequency differences base on orientation

The effects of orientation parameter on ΔTS at multiple frequencies were confirmed. For this analysis, the fatness was set to 40%, and the density and sound speed contrasts were set to 1.0357 and 1.0279, respectively, to calculate the multi-frequency difference. Overall, the length-specific ΔTS in krill according to the distribution of orientations were similar for N (−20°, 28°), (−8°, 9°), and N (0°, 27°), except for N (11°, 4°) (Fig. 3). For the same krill length, the ΔTS with orientation were up to 2.8, 3.7, 5.7, 4.8, 4.3, and 5.9 dB at ΔTS70–38 kHz, ΔTS120–38 kHz, ΔTS120–70 kHz, ΔTS200–38 kHz, ΔTS200–70 kHz, and ΔTS200–120 kHz, respectively. In a study on the acoustic scattering properties of Meganyctiphanes novegica, a major krill species in the northern hemisphere, using SDWBA, the distribution of orientation had the greatest impact among the various parameters of the model (Peña & Calise, 2016). Our findings also suggest that orientation distribution substantially affects both the pattern of multi-frequency differences and TS. Thus, accurately observing krill orientation during acoustic surveys is crucial for enhancing the precision of krill acoustic identification from gathered data. However, a standardized protocol for observing krill orientation distributions is currently lacking, underlining the need for a uniform methodology to minimize errors arising from varied observation techniques.

fas-27-4-225-g3
Fig. 3. ΔTS patterns (ΔTS70–38 kHz, ΔTS120–38 kHz, ΔTS120–70 kHz, ΔTS200–38 kHz, ΔTS200–70 kHz, and ΔTS200–120 kHz) versus length considering distribution of orientation. The red circle indicates the distribution of orientation of Antarctic krill recommended by CCAMLR (SC-CAMLR, 2010). The blue square represents the distribution of orientation of ice krill reported by La et al. (2015a). The black triangle presents the distribution of orientation of Antarctic krill measured by Lawson et al. (2006), and the green inverted triangle indicates the distribution of orientation of Antarctic krill reported by Conti & Demer (2006). Krill measuring less than 35 mm are categorized as ice krill, while those exceeding 35 mm are classified as Antarctic krill, following the criteria established by Sala et al. (2002). TS, target strength.
Download Original Figure

Conclusion

In this study, we employed the SDWBA model to ascertain the ΔTS variations at multiple frequencies in Antarctic and ice krill, focusing on parameters such as density and sound speed contrasts, fatness, and orientation. The study revealed that both krill fatness and orientation play significant roles in influencing the variation of ΔTS in krill. These findings are instrumental for the acoustic identification and biomass estimation of Antarctic krill and ice krill, key biological resources in the Southern Ocean. Additionally, this research lays a foundational groundwork for future applications of broadband acoustic systems in this field.

Competing interests

No potential conflict of interest relevant to this article was reported.

Funding sources

This research was supported by the Korea Institute of Marine Science & Technology Promotion (KIMST) through funding provided by the Ministry of Oceans and Fisheries (KIMST-20220547, PM24060), and by the Korea Polar Research Institute (KOPRI) through a grant funded by the Ministry of Oceans and Fisheries (KOPRI PE24110).

Acknowledgements

Not applicable.

Availability of data and materials

Upon reasonable request, the datasets of this study can be available from the corresponding author.

Ethics approval and consent to participate

Not applicable.

References

1.

Calise L, Skaret G. Sensitivity investigation of the SDWBA Antarctic krill target strength model to fatness, material contrasts and orientation. CCAMLR Sci. 2011; 18:97-122.

2.

Campbell RW, Dower JF. Role of lipids in the maintenance of neutral buoyancy by zooplankton. Mar Ecol iProg Ser. 2003; 263:93-9

3.

Conti SG, Demer DA. Improved parameterization of the SDWBA for estimating krill target strength. ICES J Mar Sci. 2006; 63:928-35

4.

Chu D, Wiebe PH. Measurements of sound-speed and density contrasts of zooplankton in Antarctic waters. ICES J Mar Sci. 2005; 62:818-31

5.

Everson IK. Biology, ecology, and fisheries. Oxford: Blackwell Science. 2000.

6.

Foote KG. Speed of sound in Euphausia superba. J Acoust Soc Am. 1990; 87:1405-8

7.

Foote KG, Everson I, Walkins JL, Bone DG. Target strengths of Antarctic krill (Euphausia superba) at 38 and 120 kHz. J Acoust Soc Am. 1990; 87:16-24

8.

Forman KA, Warren JD. Variability in the density and sound-speed of coastal zooplankton and nekton. ICES J Mar Sci. 2010; 67:10-8

9.

Kang M, Furusawa M, Miyashita K. Effective and accurate use of difference in mean volume backscattering strength to identify fish and plankton. ICES J Mar Sci. 2002; 59:794-804

10.

La HS, Lee H, Fielding S, Kang D, Ha HK, Atkinson A, et al. High density of ice krill (Euphausia crystallorophias) in the Amundsen sea coastal polynya, Antarctica. Deep Sea Res I Oceanogr Res Pap. 2015a; 95:75-84

11.

La HS, Lee H, Kang D, Lee SH, Shin HC. Ex situ echo sounder target strengths of ice krill Euphausia crystallorophias. Chin J Oceanol Limnol. 2015b; 33:802-8

12.

Lawson GL, Wiebe PH, Ashjian CJ, Chu D, Stanton TK. Improved parametrization of Antarctic krill target strength models. J Acoust Soc Am. 2006; 119:232-42

13.

Lawson GL, Wiebe PH, Ashjian CJ, Stanton TK. Euphausiid distribution along the Western Antarctic peninsula—part B: distribution of euphausiid aggregations and biomass, and associations with environmental features. Deep Sea Res II Top Stud Oceanogr. 2008; 55:432-54

14.

Peña M, Calise L. Use of SDWBA predictions for acoustic volume backscattering and the Self-Organizing Map to discern frequencies identifying Meganyctiphanes norvegica from mesopelagic fish species. Deep Sea Res 1 Oceanogr Res Pap. 2016; 110:50-64

15.

Sala A, Azzali M, Russo A. Krill of the Ross Sea: distribution, abundance and demography of Euphausia superba and Euphausia crystallorophias during the Italian Antarctic expedition (January–February 2000). Sci Mar. 2002; 66:123-33

16.

SC-CAMLR. Report of the fifth meeting of the subgroup on acoustic survey and analysis methods.In Proceedings of the Report of the Twenty-Ninth Meeting of the Scientific Committee (sc-camlrxXIX). 2010 Hobart, Australia. .

17.

SC-CAMLR. Report of the first meeting of the subgroup on acoustic survey and analysis method (SGASAM).In Proceedings of the Report of the Twenty-Fourth Meeting of the Scientific Committee (SC-CAMLRXXIV/BG/3). 2005 Hobart, Australia. .

18.

Stanton TK, Chu D. Review and recommendations for the modelling of acoustic scattering by fluid-like elongated zooplankton: euphausiids and copepods. ICES J Mar Sci. 2000; 57:793-807

19.

Watkins JL, Hewitt R, Naganobu M, Sushin V. The CCAMLR 2000 survey: a multinational, multi-ship biological oceanography survey of the Atlantic sector of the Southern Ocean. Deep Sea Res II Top Stud Oceanogr. 2004; 51:1205-13