A Study on the Characteristics of Fishery Resources Distribution in Coastal Waters of Yeongil Bay Using Acoustic Survey

: In this study, the identiﬁcation of dominant ﬁsh species in the East Sea was conducted using the dB-difference method. The survey was conducted using the two frequencies of 38 and 120 kHz in transect 6 of the southern part of the East Sea. Information on ﬁsh species was identiﬁed using ﬁshing gear and e-DNA, and the dominant target ﬁsh species were selected and analyzed as cod, anchovy, common squid, and herring. The dB-difference range for each ﬁsh species was set to − 0.86 dB < ∆ MVBS 38–120 kHz < 0.82 dB for cod and to the range of 2.66 dB < ∆ MVBS 38–120 kHz < 2.84 dB for anchovy. The dB-difference of the common squid was set to − 0.36 dB < ∆ MVBS 38–120 kHz < 1.25 dB and to the range of 0.88 dB < ∆ MVBS 38–120 kHz < 2.28 dB for herring; the ﬁsh species were then identiﬁed in the echograms. When comparing the results of swimming depths by ﬁsh species and previous studies, cod was detected mainly at the bottom of the sea, and anchovy and common squid were detected mainly at a depth of 50 m. Herring was detected to be mainly distributed in water depths from 50 to 150 m.


Introduction
The Korean Peninsula is surrounded by water on three sides with the East Sea, the West Sea, and the South Sea, and the characteristics of the ocean and the distribution of fishery resources differ for each sea area. The East Sea has no islands or bays, and the coastline is relatively straight and featureless. The width of the continental shelf is less than 25 km and is very narrow, with an average of 18 km [1][2][3]. The southern part of the East Sea is affected by the Tsushima Warm Current with high temperature and high salinity and the North Korea Cold Current with low temperature and low salinity. In this sea area, warm-water species and cold-water species show mixed distribution. Thus, various types of water masses are formed in the southern part of the East Sea, and the area is rich in nutrients and food organisms. The main fish species caught include Pacific cods (Gadus macrocephalus), common squids (Todarodes pacificus), herrings (Clupea pallasii), anchovies (Engraulis japonicus), and Japanese amberjacks [4,5]. Currently, the leading countries in the fishery industry have adopted the TAC (total allowable catch) system for the efficient management of fishery resources, managing fish stock for major species. However, accurate information on the fish stocks needs to be obtained for efficient management of fishery resources. The estimation of fish stocks using a hydroacoustic survey uses acoustic data to identify the presence of marine organisms. The hydroacoustic survey enables a survey of large sea areas in a short time as well as acquiring real-time information for all water columns, and it is widely used as a useful tool for estimating fish stocks in major countries in the fishing industry, such as Norway, the United States, and Canada. Among the 2 of 13 techniques of hydroacoustic survey, the dB-difference method is actively utilized as a tool for evaluating the density and biomass of the target organisms, as the method identifies fish species by the difference in frequency characteristics. There are various methods for evaluating fishery resources and the density and biomass of organisms using acoustic survey, such as echo counting, school detect, and time-varied threshold, but in this study, we use dB-difference in order to determine the spatiotemporal distribution of fish species, identifying the target strength of the target fish species and the length distribution of the collected fish species [6][7][8][9].
In this study, acoustic data were collected by means of a hydroacoustic survey in the southern part of the East Sea of South Korea, and the total length-body weight information of fish species was obtained using the bottom trawl survey. Then, seawater was collected, and e-DNA was analyzed to obtain information on fish species. For the efficient management and evaluation of fisheries resources, it is very important to understand the spatiotemporal distribution of fish species. Therefore, this study identifies fish using two frequency dB-difference methods and compares the spatiotemporal distribution of each fish species with that of a prior study.

Study Waters and Period
A training ship named Gaya (1737 G/T) belonging to Pukyong National University was used, and the acoustic survey was performed using the scientific echosounder mounted on the training ship Gaya. Acoustic data from the southern part of the East Sea were collected from 4 to 5 January 2020, and acoustic survey and trawl survey were conducted in parallel with a total of five survey lines from coastal waters of Pohang to coastal waters of Ulsan. Acoustic data were collected while maintaining the speed of the ship at 7-10 knots. During the survey period, fish sampling using the bottom trawl was performed four times, and the survey locations are shown in Table 1 and Figure 1. for all water columns, and it is widely used as a useful tool for estimating fish stocks in major countries in the fishing industry, such as Norway, the United States, and Canada. Among the techniques of hydroacoustic survey, the dB-difference method is actively utilized as a tool for evaluating the density and biomass of the target organisms, as the method identifies fish species by the difference in frequency characteristics. There are various methods for evaluating fishery resources and the density and biomass of organisms using acoustic survey, such as echo counting, school detect, and time-varied threshold, but in this study, we use dB-difference in order to determine the spatiotemporal distribution of fish species, identifying the target strength of the target fish species and the length distribution of the collected fish species [6][7][8][9].
In this study, acoustic data were collected by means of a hydroacoustic survey in the southern part of the East Sea of South Korea, and the total length-body weight information of fish species was obtained using the bottom trawl survey. Then, seawater was collected, and e-DNA was analyzed to obtain information on fish species. For the efficient management and evaluation of fisheries resources, it is very important to understand the spatiotemporal distribution of fish species. Therefore, this study identifies fish using two frequency dB-difference methods and compares the spatiotemporal distribution of each fish species with that of a prior study.

Study Waters and Period
A training ship named Gaya (1737 G/T) belonging to Pukyong National University was used, and the acoustic survey was performed using the scientific echosounder mounted on the training ship Gaya. Acoustic data from the southern part of the East Sea were collected from 4 to 5 January 2020, and acoustic survey and trawl survey were conducted in parallel with a total of five survey lines from coastal waters of Pohang to coastal waters of Ulsan. Acoustic data were collected while maintaining the speed of the ship at 7-10 knots. During the survey period, fish sampling using the bottom trawl was performed four times, and the survey locations are shown in Table 1 and Figure 1.

Acoustic Equipment System and Data Analysis Method
The acoustic system used for the survey was a split-beam scientific echosounder (EK60, Kongsberg Simrad, Norway), and acoustic data were collected for frequencies 38 and 120 kHz. The acoustic data were stored by continuously receiving position information from DGPS (SPR-1400, Samyoung, Seongnam, Korea). The acoustic data collected in situ were subsequently processed in the laboratory using acoustic data analysis software (Echoview V 8.0, Echoview Software Pty Ltd., Hobart, Australia). The flow diagram of data processing for the fish species identification is shown in Figure 2. As shown in the figure, after filtering the noise around the sea surface, sea bottom, and other places, the integration section is set to compose a matrix for each frequency, thus creating a new echogram. In this study, a cell size of 5 ping × 5 m (width × height) was used to examine the dB-difference. When the dB-difference between fish species becomes clear, the range of values is set, the data range bitmap is created, and the mask is created with an echo that matches the size of the cell at 38 kHz. We then extract echoes consistent with the 38 kHz echo, which is removed noise, by dividing them again into ping intervals. Through the above methods, fish species can be identified by separating the echoes of fish suitable for each characteristic.

Acoustic Equipment System and Data Analysis Method
The acoustic system used for the survey was a split-beam scientific echosounder (EK60, Kongsberg Simrad, Norway), and acoustic data were collected for frequencies 38 and 120 kHz. The acoustic data were stored by continuously receiving position information from DGPS (SPR-1400, Samyoung, Korea). The acoustic data collected in situ were subsequently processed in the laboratory using acoustic data analysis software (Echoview V 8.0, Echoview Software Pty Ltd., Australia). The flow diagram of data processing for the fish species identification is shown in Figure 2. As shown in the figure, after filtering the noise around the sea surface, sea bottom, and other places, the integration section is set to compose a matrix for each frequency, thus creating a new echogram. In this study, a cell size of 5 ping × 5 m (width × height) was used to examine the dB-difference. When the dB-difference between fish species becomes clear, the range of values is set, the data range bitmap is created, and the mask is created with an echo that matches the size of the cell at 38 kHz. We then extract echoes consistent with the 38 kHz echo, which is removed noise, by dividing them again into ping intervals. Through the above methods, fish species can be identified by separating the echoes of fish suitable for each characteristic. Using the above method and procedure, each species can be identified by separating the echoes of the respective fish species suitable for the characteristics [10].

Analysis of Sampled Fish Species
In order to identify fish species in the target waters for survey using the acoustic data, biological information of the fish species, such as composition and total length-body weight information of the fish species inhabiting the sea area, is required. Therefore, in Using the above method and procedure, each species can be identified by separating the echoes of the respective fish species suitable for the characteristics [10].

Analysis of Sampled Fish Species
In order to identify fish species in the target waters for survey using the acoustic data, biological information of the fish species, such as composition and total length-body weight information of the fish species inhabiting the sea area, is required. Therefore, in this study, the total length-body weight information of fish was obtained using the bottom trawl survey, and the eDNA information of fish was acquired by collecting surface seawater within the survey area [11]. As a result of identifying fish species using eDNA information, the sampled species were 33.9% of Gadus macrocephalus, 28.3% of Clupea pallasii, 12.7% of Glyptocephalus stelleri, 4.6% of Engraulis japonicus, 4.2% of sea bass (Lates calcarifer), and 3.0% of jack mackerel Trachurus symmetricus. Therefore, the target fish species of this study were selected to be Gadus macrocephalus, Clupea pallasii, and Engraulis japonicus, which appeared as dominant fish species as a result of eDNA analysis. Todarodes pacificus, a major fish species in the East Sea, was additionally selected as the target species (Table 2).

Acoustic Backscattering Strength
The range of total length of Gadus macrocephalus was 36.4-74.0 cm, and that of Engraulis japonicus was 9.5-15.4 cm. In addition, the range of total length of Todarodes pacificus was 16.0-27.0 cm, and that of Clupea pallasii was 12.0-27.0 cm. The total weight of the collected Gadus macrocephalus was 3459 g, and the total weight of Engraulis japonicus was 1910 g. Todarodes pacificus had a collection weight of 1235 g, and Clupea pallasii had a collection weight of 4617 g.
The frequency characteristics of each fish species for the dominant species were identified, and the results of a previous study were applied for the acoustic scattering characteristics for each frequency. Information on the acoustic scattering strength of the dominant species is presented as follows [12][13][14][15]. For the identification of fish species, the range of total length for each species and acoustic scattering strength were used to represent acoustic scattering strength for each frequency, and the frequency range was set using the dB-difference method. There-Appl. Sci. 2021, 11, 6627 5 of 13 fore, the dB-difference range for Gadus macrocephalus was set to be −0.86 dB < ∆MVBS 38-120 kHz < 0.82 dB, and the dB-difference range for Engraulis japonicus was set to be 2.66 dB < ∆MVBS 38-120 kHz < 2.84 dB. In addition, the dB-difference range for Todarodes pacificus was 0.36 dB < ∆MVBS 38-120 kHz < 1.25 dB, and the range for Clupea pallasii was 0.88 dB < ∆MVBS 38-120 kHz < 2.28 dB. These ranges were set for fish species identification.
For the identification of fish species, it is necessary to first identify the frequency characteristics for the fish's frequencies of 38 and 120 kHz and differences between them. In this case, the dB-difference represents the difference in the mean volume backscattering strength (MVBS) at multiple frequencies. To set ∆MVBS as positive values, target strength (TS) for each frequency of the target species for separation is compared, and the frequency with the smaller TS value is subtracted from the frequency with the large TS value. In general, the frequency of fish showed higher values at 38 kHz than at 120 kHz. Therefore, in a new echogram of 38 kHz and 120 kHz generated in the form of a matrix, ∆MVBS can be expressed by the following (Equation (1)).
The method for identifying each fish species is shown in the figure below by applying the sampled data and the acoustic scattering strength for each fish species (Figures 3-6).
. Sci. 2021, 11, x FOR PEER REVIEW 5 of 12 For the identification of fish species, the range of total length for each species and acoustic scattering strength were used to represent acoustic scattering strength for each frequency, and the frequency range was set using the dB-difference method. Therefore, the dB-difference range for Gadus macrocephalus was set to be −0.86 dB < ∆MVBS 38-120 kHz < 0.82 dB, and the dB-difference range for Engraulis japonicus was set to be 2.66 dB < ∆MVBS38-120 kHz < 2.84 dB. In addition, the dB-difference range for Todarodes pacificus was 0.36 dB < ∆MVBS38-120 kHz < 1.25 dB, and the range for Clupea pallasii was 0.88 dB < ∆MVBS 38-∆MVBS TS 38 kHz -TS 120 kHz SV 38 kHz -SV 120 kHz (1) The method for identifying each fish species is shown in the figure below by applying the sampled data and the acoustic scattering strength for each fish species (Figures 3-6).  For the identification of fish species, the range of total length for each species and acoustic scattering strength were used to represent acoustic scattering strength for each frequency, and the frequency range was set using the dB-difference method. Therefore, the dB-difference range for Gadus macrocephalus was set to be −0.86 dB < ∆MVBS 38-120 kHz < 0.82 dB, and the dB-difference range for Engraulis japonicus was set to be 2.66 dB < ∆MVBS38-120 kHz < 2.84 dB. In addition, the dB-difference range for Todarodes pacificus was 0.36 dB < ∆MVBS38-120 kHz < 1.25 dB, and the range for Clupea pallasii was 0.88 dB < ∆MVBS 38-120 kHz < 2.28 dB. These ranges were set for fish species identification.
For the identification of fish species, it is necessary to first identify the frequency characteristics for the fish's frequencies of 38 and 120 kHz and differences between them. In this case, the dB-difference represents the difference in the mean volume backscattering strength (MVBS) at multiple frequencies. To set ∆MVBS as positive values, target strength (TS) for each frequency of the target species for separation is compared, and the frequency with the smaller TS value is subtracted from the frequency with the large TS value. In general, the frequency of fish showed higher values at 38 kHz than at 120 kHz. Therefore, in a new echogram of 38 kHz and 120 kHz generated in the form of a matrix, ∆MVBS can be expressed by the following (Equation (1)).

∆MVBS TS 38 kHz -TS 120 kHz
SV 38 kHz -SV 120 kHz The method for identifying each fish species is shown in the figure below by applying the sampled data and the acoustic scattering strength for each fish species (Figures 3-6).

Vertical Distribution by Species
Using the dB-differences, the vertical distribution of Gadus macrocephalus is shown in the figure below, and the depth of the survey area was found to be in the range from 50 m to 250 m or more. The horizontal axis represents the time, the number at the top indicates the number of the survey line according to the movement of the survey line, the vertical axis represents the depth of water, and the thin solid line represents the water depth at 50 m intervals. The color scale on the right shows the intensity of sound waves received after they were reflected by the fish school. The higher the fish School, the darker brown that represents it.

Vertical Distribution by Species
Using the dB-differences, the vertical distribution of Gadus macrocephalus is shown in the figure below, and the depth of the survey area was found to be in the range from 50 m to 250 m or more. The horizontal axis represents the time, the number at the top indicates the number of the survey line according to the movement of the survey line, the vertical axis represents the depth of water, and the thin solid line represents the water depth at 50 m intervals. The color scale on the right shows the intensity of sound waves received after they were reflected by the fish school. The higher the fish School, the darker brown that represents it.

Vertical Distribution by Species
Using the dB-differences, the vertical distribution of Gadus macrocephalus is shown in the figure below, and the depth of the survey area was found to be in the range from 50 m to 250 m or more. The horizontal axis represents the time, the number at the top indicates the number of the survey line according to the movement of the survey line, the vertical axis represents the depth of water, and the thin solid line represents the water depth at 50 m intervals. The color scale on the right shows the intensity of sound waves received after they were reflected by the fish school. The higher the fish School, the darker brown that represents it.      As a result of examining the vertical distribution of each fish species, Gadus macrocephalus was mostly detected at the sea bottom, and Engraulis japonicus and Todarodes pacificus were mainly detected at a depth of 50 m. Clupea pallasii was mostly detected at a depth between 50 and 150 m.

Discussion
The results of previous studies on the distribution of fish species by water depth showed that Engraulis japonicus shoal was mainly distributed at a depth range of 12-50 m, and the volume scattering strength distribution of Engraulis japonicus shoal was around −50 dB. In addition, in another previous study, the volume scattering strength of the large Engraulis japonicus shoal with a high level of aggregation was in the range of −44.0 to −28.0 dB, and the large Engraulis japonicus shoal was distributed at a depth range of 10-75 m. Furthermore, the previous study on the distribution of Atlantic herring (Clupea harengus) by water depth showed the  As a result of examining the vertical distribution of each fish species, Gadus macrocephalus was mostly detected at the sea bottom, and Engraulis japonicus and Todarodes pacificus were mainly detected at a depth of 50 m. Clupea pallasii was mostly detected at a depth between 50 and 150 m.

Discussion
The results of previous studies on the distribution of fish species by water depth showed that Engraulis japonicus shoal was mainly distributed at a depth range of 12-50 m, and the volume scattering strength distribution of Engraulis japonicus shoal was around −50 dB. In addition, in another previous study, the volume scattering strength of the large Engraulis japonicus shoal with a high level of aggregation was in the range of −44.0 to −28.0 dB, and the large Engraulis japonicus shoal was distributed at a depth range of 10-75 m. Furthermore, the previous study on the distribution of Atlantic herring (Clupea harengus) by water depth showed the As a result of examining the vertical distribution of each fish species, Gadus macrocephalus was mostly detected at the sea bottom, and Engraulis japonicus and Todarodes pacificus were mainly detected at a depth of 50 m. Clupea pallasii was mostly detected at a depth between 50 and 150 m.

Discussion
The results of previous studies on the distribution of fish species by water depth showed that Engraulis japonicus shoal was mainly distributed at a depth range of 12-50 m, and the volume scattering strength distribution of Engraulis japonicus shoal was around −50 dB. In addition, in another previous study, the volume scattering strength of the large Engraulis japonicus shoal with a high level of aggregation was in the range of −44.0 to −28.0 dB, and the large Engraulis japonicus shoal was distributed at a depth range of 10-75 m.
Furthermore, the previous study on the distribution of Atlantic herring (Clupea harengus) by water depth showed the distribution of the species in the range of 10-125 m. A previous study on the distribution of Todarodes pacificus in the East Sea of South Korea showed the distribution of the species mainly at 60 m depth. Therefore, when comparing the results of this study with those of previous studies, it can be concluded that the distribution of species by depth was similar [14][15][16][17][18][19].
The spatiotemporal distribution of each fish species was represented by extracting the Nautical Area-Scattering Coefficient (NASC m 2 /nm 2 ) and is illustrated in the figure below, and the size of NASC was represented by the height of the bar graph.
The results of the spatiotemporal distribution of Gadus macrocephalus showed a higher level of distribution in the south of the area, and a higher level of distribution was identified in coastal waters than in the open ocean. In addition, Gadus macrocephalus was distributed in the open ocean waters of Ulsan ( Figure 11). studies and the location of the distribution of Gadus macrocephalus obtained using hydro acoustic survey were very similar. However, there is little research on Gadus macrocephalu movement, and, thus, it is thought that further studies are needed to obtain data on var ous routes of Gadus macrocephalus by year ( Figure 12).
The results of the spatiotemporal distribution of Engraulis japonicus showed high ag gregation of the species in the coastal waters and open ocean of Ulsan in the survey are and low aggregation in the adjacent waters of Pohang. In a previous study that identifie the distribution of Engraulis japonicus in the southern part of the East Sea, lower aggrega tion of the species was shown in coastal waters of Pohang, and high aggregation was ob served in coastal waters and open ocean of Ulsan. Comparing the two sets of results, it judged that the two results are considerably similar in many aspects [16]. However, th Engraulis japonicus shoal is a migratory species and shows aggregation at various depth and, thus, further research is thought to be necessary (Figures 13 and 14).
The spatiotemporal distribution of Todarodes pacificus showed a high aggregation i the coastal waters and open ocean in Ulsan and a relatively low aggregation in the coast waters of Pohang ( Figure 15).
The spatiotemporal distribution of Clupea pallasii was relatively high in the coast waters and open ocean of Ulsan and the northern waters of Ulsan, whereas the distribu tion was relatively low in the coastal waters of Pohang and adjacent sea ( Figure 16).    The results of the spatiotemporal distribution of Engraulis japonicus showed high aggregation of the species in the coastal waters and open ocean of Ulsan in the survey area and low aggregation in the adjacent waters of Pohang. In a previous study that identified the distribution of Engraulis japonicus in the southern part of the East Sea, lower aggregation of the species was shown in coastal waters of Pohang, and high aggregation was observed in coastal waters and open ocean of Ulsan. Comparing the two sets of results, it is judged that the two results are considerably similar in many aspects [16]. However, the Engraulis japonicus shoal is a migratory species and shows aggregation at various depths, and, thus, further research is thought to be necessary (Figures 13 and 14).
Appl. Sci. 2021, 11, x FOR PEER REVIEW 10 of Figure 12. The migration route of pacific cod [20].      The spatiotemporal distribution of Clupea pallasii was relatively high in the coastal waters and open ocean of Ulsan and the northern waters of Ulsan, whereas the distribution was relatively low in the coastal waters of Pohang and adjacent sea ( Figure 16).
This study used dB-difference to identify the species, space, and time distribution in the southern part of the East Sea waters of the Republic of Korea using various methods and distribution. Fish species were identified using frequencies of 38 and 120 kHz, but the distribution of fish species with similar target strengths and lengths can be mixed. However, since the fish species mainly distributed in the research waters of this study are the target fish species of the study, it is considered that the rate of other fish species is quite low.  This study used dB-difference to identify the species, space, and time distribution in the southern part of the East Sea waters of the Republic of Korea using various methods and distribution. Fish species were identified using frequencies of 38 and 120 kHz, but the distribution of fish species with similar target strengths and lengths can be mixed. However, since the fish species mainly distributed in the research waters of this study are the target fish species of the study, it is considered that the rate of other fish species is quite low.