Estimation of Moon Jellyfish Aurelia coerulea Using Hydroacoustic Methods off the Coast of Tongyeong, Korea

통영 인근해역에서 음향기법을 이용한 보름달물해파리(Aurelia coerulea)의 밀도 추정

Abstract

Moon jellyfish Aurelia coerulea are highly abundant off the coast of Tongyeong, Korea. We measured the density of A. coerulea in this area using a scientific echosounder at frequencies of 38 and 120 kHz, and then applied a distorted wave Born approximation (DWBA) model to calculate the target strength of the echosounder at each frequency. Then, we used the frequency difference method to extract jellyfish echo signals and estimate the A. coerulea density. A. coerulea was evenly distributed throughout the water column; the backscattering strength ranged from -75 to -65 dB. In May and August, the A. coerulea densities at survey lines 3 and 4 were estimated at 1.5-1.6 and 0.2-0.9 g/m2, with mean weighted densities of 1.04 and 0.48 g/m2, respectively. In September, the A. coerulea densities estimates in Jaran Bay and Goseong Bay were 0.6-2.1 and 0.1-0.4 g/m2, with mean weighted densities of 1.25 and 0.24 g/m2, respectively.

서 론

최근 지구온난화, 부영양화, 간척지 매립, 수중 건설 등으로젤라틴으로 구성된 해파리가 대량 출현하고 있다. 출현하고 있는 해파리 중 보름달물해파리(Aurelia coerulea)는 연안 해역의 넓은 범위에서 대량 발생하는 종으로 우리나라, 중국, 일본, 대만, 북서 유럽, 북아메리칸, 흑해 등의 전 세계적으로 출현하였다 (Lucas, 2001; Mutlu, 2001; Uye et al., 2003; Dong et al., 2010; Albert, 2011; Aoki et al., 2012). 보름달물해파리는 주로동물플랑크톤, 자치어를 포식하기 때문에 해양 생태계에 교란을 일으키고, 어업 중에 어구에 대량으로 어획되기도 하여 어업 활동을 방해할 뿐만 아니라 발전소에서 해수를 유입하면서 해파리가 취수구에 대량 유입되어 출력을 낮추거나 가동을 중단시키기도 한다 (Toyokawa et al., 2000, 2010; Lucas, 2001; Ishii et al., 2008; Uye, 2011). 이와 같은 피해를 주는 보름달물해파리는 생활사 가운데 폴립 과정 중 횡분열화를 일으켜 많은 에피라를 생성시키고, 무성생식으로 개체수가 대량 증식하기 때문에 대량으로 발생하게 된다. 이러한 피해를 낮추기 위해서는 보름달물해파리의 서식환경, 생활사, 행동 등의 생태학적 특성을 파악하며, 대상생물의 분포 및 현존량을 파악하는 연구가 필요하다. 해파리의 분포 및 현존량의 모니터링은 샘플링 기어를 이용한 채집 조사, 스킨스쿠버를 이용하는 방법과 목시관측이나 항공 관측과 음향 조사가 주를 이루어지고 있다. 샘플링 기어를 이용한 조사는 대상생물의 종과 사이즈, 무게 등을 바로 알 수는 있으나 채집 조사 수행 중에 샘플이 손상되는 경우도 발생하고, 일부 정점과 수심에 대한 샘플 정보를 얻을 수 있어 반 정량적인 방법으로 평가되고 있다. 스킨스쿠버를 이용하는 방법과 목시, 항공 관측의 경우에는 눈으로 바로 확인 할 수는 있으나 관측자의 주관적인 오차가 유입될 수 있으며, 짧은 시간에 인력이 허용되는 수심에서만 조사가 이루어 지는 단점이 있다. 그에 비해서 음향기법을 이용한 조사는 시간에 적게 들고 많은 자료를 얻을 수 있기 때문에 가장 효과적인 방법으로 알려져 있다(Bå mstedt et al., 2003; Brierley et al., 2004; Lee et al., 2007; Colombo et al., 2009; Hirose et al., 2009; Kang et al., 2010; Graham et al., 2010).

또한, 해파리의 밀도를 평가하기 위한 조사에는 어군탐지기및 음향카메라를 이용한 방법이 여러 연구를 통해 적용되었다(Brierley et al., 2001; Colombo et al., 2003, 2009; Han and Uye, 2009, Makabe et al., 2012). 본 연구에서의 대상으로 하는 보름달물해파리의 경우, 대부분 서로 밀집하여 분포하기 때문에 어군탐지기를 이용하여 계측하기 위해서는 에코 계수법보다는 에코를 적분하는 방법을 이용하는 것이 더 효율적이다(Colombo et al., 2009). 에코 적분법에서의 밀도를 추정하는 방법은 현장에서 계측된 대상생물의 계측된 체적후방산란강도(volume backscattering strength, SV)를 대상생물의 반사강도(target strength, TS)로 나누면 간편하게 계산할 수 있는데, 해파리 TS는 유영행동에 따른 박동과 유영자세각 등에 의해 산란 단면적이 변화하기 때문에 TS 값을 적용하는데 문제점이 따른다. 또한, 음향카메라를 이용한 방법은 자료 분석 시에 수동적인 방법으로 해파리 에코를 카운팅해야 한다. 이 방법은 좋은 자료를 수집하기 위해서는 조사 시에 선속을 천천히 유지하여 조사하여야 하므로 어군탐지기를 사용하는 것보다 조사보다 많이 소요된다.

본 연구에서는 보름달물해파리로 인해 발생하는 피해 방지를 위한 밀도 추정 모니터링을 위해 보름달물해파리가 대량 출현하는 통영 하구역에서 어군탐지기를 이용하여 정선조사를 실시하였으며, 보름달물 해파리가 집중적으로 출현하는 통영 고성만과 자란만에서 수중음향카메라와 과학어군탐지를 이용한 집중 조사를 수행하여 보름달물해파리의 분포와 밀도를 파악하였다.

 

재료 및 방법

본 조사는 통영 욕지도 및 자란만 해역에서 출현하는 보름달물해파리의 밀도를 추정하기 위해 2018년 5월 14일, 8월 9일, 9월 16일 주간에 총 3회 음향조사를 실시하였다. 5월과 8월의 음향 조사 정선은 총 4개의 수평으로 통영과 남해 사이에 정선을 설계하였고, 정선 길이는 120.19 km, 조사 면적은 181 km2이었다(Fig. 1). 9월에는 해파리 출현에 따른 집중조사를 위해 통영의 고성만과 자란만에서 음향조사를 실시하였으며, 대상생물의 종 식별을 위해 수중 음향 카메라(DIDSON, Ocean Marine Industries Inc., Chesapeake, VA, USA) 조사를 동시에 실시하였다. 음향 조사 시스템은 주파수 38과 120 kHz split-beam 방식의 과학어군탐지기(EK60, Simrad, Norway)를 이용하였고, DGPS (differential global positioning system; DSPR-1400, Samyung, Korea)로 연속적으로 위치 정보를 수신하여 과학어군탐지기 시스템에 입력하였다. 노트북에는 위치정보가 입력된 음향자료를 LAN (local area network) 시스템과 RS-232C 인터페이스를 통하여 연속적으로 하드디스크에 수록하였다. 음향 조사는 주파수 38과 120 kHz 트랜스듀서는 선박의 현측에 부착하여 선속을 5-7 knots로 유지하면서 수행되었다. 조사 선박의 반대쪽 현측에는 음향카메라를 설치하였으며, 1.8 MHz의 초음파를 이용하여 수직방향에 대한 21 프레임의 화상 신호를 저장하였고, 빔 폭 0.3-0.6°간격으로 수평각도 약 29°, 수직각도 12°의 영역을 탐지 할 수 있다.

 

Fig. 1. The moon jellyfsh Aurelia coerulea acoustic survey line performed in the (a) Tongyeong costal sea area and the (b) intensive surveyline in Goseung-Jaran bay.

 

조사 해역에 분포하고 있는 보름달물해파리의 사이즈를 파악하기 위하여 음향조사 수행 중 해파리가 탐지된 정점에서 지름60 cm, 그물코 600 μm의 RN80 네트에 유량계를 부착하여 경사로 채집하였다. 채집된 해파리의 체장은 체장판을 이용하여 우산의 직경을 측정하였다.

해파리의 주파수 차이를 파악하기 위하여 DWBA (distorted wave born approximation) 이론모델을 이용하여 주파수 38과 120 kHz에 대한 TS를 계산하였다(Monger et al., 1998). DWBA 모델은 대상 생물의 체형을 잘게 절단한 원통으로 엄밀하게 근사하여 나타낸 것으로 다른 이론 모델에 비하여 복잡한 체형을 묘사할 수 있다. DWBA 모델 식은 다음 식 (1)에서 (5)로 표현할 수 있다(McGehee et al., 1998).

                \(\begin{equation}f_{\mathrm{bs}}=\int_{\mathrm{r}_{\mathrm{ros}}} \frac{k_{2}^{1} a}{4 k_{2}}\left(\gamma_{\mathrm{k}}-\gamma_{\mathrm{p}}\right) e^{2 i k_{2} \mathrm{r}_{\mathrm{pes}}} \frac{J_{1}\left(2 k_{2} \mathrm{acos} \beta_{\mathrm{tilt}}\right.}{\cos \beta_{\mathrm{tit}}}\left|dr_{\mathrm{pos}}\right|\end{equation}\)   ………(1)

                \(\begin{equation}\sigma_{\mathrm{bs}}=\left|f_{\mathrm{bs}}\right|^{2}\end{equation}\)   ……………………… (2)

                \(\begin{equation}\mathrm{TS}=10 \log \sigma_{\mathrm{bs}}\end{equation}\)   …………………… (3)

여기서, r_pos는 산란체의 중심으로부터의 위치이고, k는 파수로서 k는 2π/λ로 나타낼 수 있으며, λ는 파장이다. α는 절단원통의 반경(m)이고, βt_tilt은 입사하는 음파와 원통 축 사이의 각도, J₁은 제 1종 베셀함수이다. r_k, r_ρ는 해파리의 음속비(h), 밀도비(g)를 이용하여 다음 식 (4)와 (5)로 나타낼 수 있다.

                   \(\begin{equation}r_{\mathrm{k}}=\frac{1}{g h^{2}}-1\end{equation}\)    …………………… (4)

                   \(\begin{equation}r_{\rho}=1-\frac{1}{g}\end{equation}\)   …………………… (5)

보름달물해파리의 체형은 공기중 우산의 직경이 10.1 cm인해파리가 수중에서 박동에 따라 우산의 직경이 최소일 때와 최대일 때를 디지타이징 하였다(Fig. 2). 보름달물해파리의 음속비는 1.0004, 밀도비는 1.0031 (Yoon et al., 2015), 평균 TS를 추정하기 위한 유영자세각의 분포는 평균(±표준편차) 8 0.6° (± 1.0° )을 확률밀도함수로 적용하였다(Graham et al., 2010). 과학어군탐지기로 수집한 음향자료는 후처리 소프트웨어(Echoveiw Ver 8.0, Echoview Software Pty Ltd., Australia)를 이용하여 분석하였다. 주파수 38과 120 kHz를 이용하여 음향자료를 수집할 때 배경 잡음(background noise)와 임펄스 잡음(impulse noise)로 인한 신호 손실이 문제가 된다. 본 조사에서 수심이 깊어짐에 따라 발생하는 배경 잡음은 TVT (time varied threshold) 방법을 적용하여 제거하는데, 인위적으로 배경 잡음을 만든 후에 raw data 자료에서 인위적으로 만든 잡음을 다시 가감한다. 임펄스 잡음은 먼저 해수면과 해저면을 데이터 처리 제외선을 설정하여 잡음소거 알고리즘의 데이터 처리 과정에서 불필요한 부분을 제외시켜준다. 그 후 smoothing 과정을 수행하는데 smoothing은 설정한 수직 샘플 범위의 값을 평균하고 그 값이 설정된 역치 값보다 크면 그것을 잡음의 후보로 식별한다. 본 연구에서는 역치 값을 -160 dB로 설정하였다. 잡음의 후보로 판단된 값은 Anderson et al. (2005)의 two-sided comparison method를 이용하여 잡음 또는 생물 에코로 판별되며, two-sided comparison method는 데이터 샘플 a, b, c가 있다고 할 때, 샘플 b를 정의하기 위해서 양 옆 데이터 샘플인 a와 c를 b와 비교하는 것이다.

음향 자료에 포함된 잡음을 제거한 후에 보름달물해파리의 에코를 분리하기 위하여 먼저, 보름달물해파리의 38과 120 kHz에 대한 주파수 특성 및 차이를 파악하여야 한다. 여기서, 주파수 차이는 다중 주파수에서 평균 체적후방산란강도(mean vol-ume backscattering strength, MVBS)의 차이를 나타내는 것으로 ∆MVBS는 양의 값을 두기 위하여 분리하고자 하는 대상으로 하는 종의 주파수별 TS를 비교하여 TS가 큰 값을 나타내는 주파수에서 작은 값을 나타내는 주파수를 가감하면 된다. 따라서, 행렬로 만들어진 38과 120 kHz의 새로운 에코그램에서∆MVBS는 다음 식 (6)으로 표현할 수 있다.

                ΔMVBS=TS_HF-TS_LF=SV_HF-SV_LF …………(6)

본 연구에서 주파수 차이를 적용하여 보름달물해파리의 에코만 추출하기 위한 셀의 크기 가로×세로를 5 ping×0.5 m를 적용하였다. 주파수 차이를 적용한 후에도 강도가 강한 어군이 완전히 제거되지 못하여 어군의 에코를 제거하기 위하여 군 탐지 기능(detect schools) 기능을 이용하여 다음 파라미터에 맞는 어군을 제거 한 후 마스크를 생성하였다. 어군이 제거 된 에코와 잡음이 제거된 120 kHz자료로 보름달물해파리의 에코를 추출하였다.

주파수 120 kHz에 기여한 보름달물해파리의 에코는 탐지거리를 고려하여 표층으로부터 3 m부터 저층까지 5월과 8월에는 1 n.mile, 9월 자란만과 고성만은 0.1 n.mile의 EDSU (elementary distance sampling unit) 간격으로 적분하고, 추출된 면적산란계수(nautical area scattering coeffcient, NASC; m2/n.mile2) 값을 이용하여 보름달물해파리의 밀도를 평가하였다(Kang et al., 2010). 또한, 보름달물해파리의 밀도 파악을 위한 주파수 120 kHz에 대한 TS-사이즈 관계식은 TS=14.72log(BD)-74.63을 이용하였다(Mutlu et al., 1996).
 

Fig. 2. Example of digitized image of moon jellyfsh Aurelia coeru-lea using the DWBA (distorted wave born approximation) model.

 

결 과

보름달물해파리의 우산 직경에 따른 반사강도

5월 네트에 채집된 보름달물해파리의 우산 직경은 0.8-4.0cm, 평균±표준편차는 2.0±0.8 cm이었으며, 8월은 우산 직경 3.0-11.5 cm, 평균±표준편차는 6.6±0.4 cm이었으며, 9월은우산 직경 2.6-10.0 cm, 평균±표준편차는 5.4±0.1 cm이었다(Fig. 3). DWBA 모델에 의한 보름달물해파리의 우산 직경이 최소일 때와 최대일 때의 사이즈에 따른 주파수 38과 120에 대한 TS를 Fig. 4와 같이 나타내었다. 주파수 38과 120 kHz 모두 보름달물해파리의 사이즈가 증가함에 따라 TS도 증가하였고, TS 사이즈에 따른 관계식은 Table 1에 나타내었다. 조사 해역에서 채집된 보름달물해파리의 우산 직경 0.8-4.0 cm를 적용하여 TS 차이 (120-38 kHz)를 파악한 결과 수중에서 박동에 의하여 보름달물해파리의 우산 직경이 최소일 때는 -1.5~2.2 dB, 최대일 때는 -1.4~1.7 dB이었다. 보름달물해파리의 에코만을 추출하기 위하여 주파수 120-38 kHz 차이의 범위는 -1.5~2.2 dB를 적용하였다.

 

Fig. 3. The length frequency analysis of moon jellyfsh Aureliacoerulea has used a net in Tongyeong coastal area on May, Augustand September. N, number of individuals; Avg., average bell diam-eter (cm); SD, standard deviation (cm).

 

Fig. 4. The target strength variations of frequency 38 and 120 kHzusing DWBA model with bell diameter in water of moon jelly fishAurelia coerulea at (a) minimum and (b) maximum.

 

보름달물해파리의 시·공간 분포

주파수 38과 120 kHz의 차이를 이용하여 해파리의 에코를 추출한 예를 Fig. 5와 같이 나타내었다. Fig. 5의 위에는 주파수 38과 120 kHz의 raw 데이터이고, raw 데이터의 셀을 가로로5 ping 세로로 0.5 m로 나누어 120-38 kHz를 하면 raw 데이터아래의 그림과 같이 나타난다. 여기서, 해파리의 주파수 차이인 1.4-4.8 dB 범위에 들어오는 에코를 120 kHz로 데이터와 매치하여 120 kHz에 기여하는 해파리의 에코만을 추출한다. 보름달물해파리는 표층부터 저층까지 고루 분포하였고, 체적산란강도는 -75에서 -65 dB사이인 것을 알 수 있었다. 이와 같이 해파리의 에코를 추출하여 시공간 분포를 파악하였다.

주파수 차이를 이용하여 보름달물해파리의 에코를 추출한 시·공간 분포를 Fig. 6와 같이 나타내었다. 5월 보름달물해파리의 면적산란계수는 정선 2에 8.1 m²/nm²으로 가장 높았고, 정선 4에서 1.1 m²/nm²으로 가장 낮은 값을 나타내었다. 보름달물해파리는 남해쪽 보다는 통영 쪽에서 상대적으로 많이 출현하는 것을 알 수 있었다. 8월 면적산란계수는 정선 모든 정점에서 1.8-4.9 m²/nm²으로 5월에 비해서 낮은 값을 나타내었다. 9월 집중조사에서는 자란만에서 11.1-3.3 m²/nm²으로 가장 높은 값을 보였으며, 고성만에서는 0.7-2.1 m²/nm²으로 나타났다.

 

Fig. 5. Example of echo extraction of moon jellyfsh Aurelia coerulea with frequency difference method at 38 and 120 kHz. SV, volume backscattering strength.

 

Fig. 6. Spatio-temporal distribution of moon jellyfsh Aurelia coerulea extracted by applying the frequency difference method in May, Au-gust and September, 2018.

 

보름달물해파리의 밀도

5월의 보름달물해파리의 정선별 밀도는 정선 3과 4에서는 1.5-1.6 g/m2로 비슷하였고, 정선 4에서 0.2 g/m² 가장 낮았으며, 조사 해역에서 보름달물해파리의 평균 밀도는 1.04 g/m²로평가되었다. 8월의 정선별 밀도는 0.2-0.9 g/m²로 나타났으며, 평균 밀도는 0.48 g/m²로 평가되었다. 9월 자란만의 정선별 밀도는 0.6-2.1 g/m²이었으며, 평균 밀도는 1.25로 가장 높은 값을 보였으며, 고성만의 정선별 밀도는 0.1-0.4 g/m²으로 나타났으며, 평균 밀도는 0.24 g/m²으로 나타났다(Fig. 7).

 

Fig. 7. Density of moon jellyfsh Aurelia coerulea per survey line in May, August, and September 2018.

 

과학어군탐지기와 수중음향카메라의 수직 분포 비교

9월 고성만과 자란만에서 수행된 집중조사에서 과학어군탐지기와 수중 음향카메라의 수직분포 결과를 다음 Fig. 8과 Fig. 9에 나타내었다. 과학어군탐지기의 경우, 수심1 m 부근에서는 보름달물 해파리를 탐지할 수 없었다. 시간대별 수직분포 결과, 고성만에서 11시 10분에 수심 6 m 부근에서 과학어군탐지기의 밀도가 0.9 m/g²의 분포를 보였으며, 수중음향카메라는 11시 08분부터 11분경에 수심 2-3 m 부근에 약 100-75 inds./m³로 가장 많은 분포를 보였다. 자란만에서는 오후 15시 16분경에 과학어군탐지기의 밀도가 1.5 m/g²로 가장 높은 밀도 분포가 나타났으며, 수중음향카메라에서는 15시 1분부터 17분경에 약75-30 inds./m3로 높은 출현량을 나타내었다. 시간대별 수직분포 결과, 보름달물해파리의 탐지영역은 2-4 m로 비슷한 양상을 나타내는 것으로 나타났다.

 

Fig. 8. Comparison of vertical distribution for moon jellyfsh Aurelia coerulea between by scientifc echosounder and underwater a coustic cameras in Goseoung Bay in September 2018.

 

Fig. 9. Comparison of vertical distribution for moon jellyfsh Aurelia coerulea between by scientifc echosounder and underwater a coustic cameras in Jaran Bay in September 2018.

 

고 찰

음향조사를 이용한 대상생물의 밀도 및 현존량 평가는 대부분어류와 크릴에 한정적이고, 해파리를 대상으로 한 연구는 Co-lombo et al. (2009)와 Graham et al. (2010)에 의하여 시도가 되었다. 해파리의 밀도는 플랑크톤 네트와 스쿠버 다이빙으로 조사가 이루어졌으나, 음향으로도 정량적으로 해파리의 분포 및 밀도를 파악할 수 있음을 증명하였다. 본 연구에서도 음향 모델로 보름달물해파리의 TS를 파악하고, 주파수 차이를 적용하여 수중에 다양하게 분포하는 해양생물들로부터 해파리의 에코신호 만을 추출하여 밀도와 현존량을 파악하였다.

해파리의 TS값은 유영자세각과 박동의 변화에 의하여 차이가 많이 나타난다. DWBA 모델을 이용하여 노무라입깃해파리의 유영자세각 0°에서 180° 사이의 변화에 따라 TS는 38 kHz에서는 약 15.0 dB, 120과 200 kHz에서는 약 30.0 dB의 차이를 보였고(Hirose et al., 2007b), 현수법을 이용한 주파수 38 kHz에서는 약 10.0 dB의 차이를 나타내었다(Hirose et al., 2010). 또한, 수중에서 보름달물해파리의 박동에 따라 Mutlu (1996)와 Mukai et al. (2009)의 연구 결과, 200 kHz에서는 10.0-16.0 dB의 차이를 보였으며, 120 kHz에서는 13.0 dB의 차이를 보고한 바 있다. Lee et al. (2007)과 Yoon et al. (2010)의 연구 결과에서는 노무라입깃해파리의 박동에 의한 실측 TS는 38과 120 kHz에서 7.8-11.7 dB의 차이를 보고한 바 있다. Lee and Hwang(2009)의 연구에서는 노무라입깃해파리를 DWBA 모델로 계산한 TS는 38, 70, 120 kHz에서 8.5-15.1 dB의 변화를 보고하였다. 이와 같이 해파리의 TS값은 유영자세각과 박동에 따라 차이가 많이 나타나므로 본 연구에서는 박동에 따른 수중에서 우산의 움직임이 최소일 때와 최대일 때를 고려하여 TS를 계산하였다. 보름달물해파리의 유영자세각은 현장에서 직접 계측은 이루어지지 않아 이전의 연구 결과를 인용하여 적용하였다. Lee et al. (2007)의 연구에서는 광학 카메라를 이용하여 노무라입깃해파리의 유영자세각을 관측한 결과, 유영자세각이 -90°에서 90°사이에서 나타났으나 대부분 수평 방향으로 이동하는 해파리가 가장 많이 관측되어 본 연구에서 적용한 값과 유사하였으나, 향후 수행될 연구에는 광학 카메라도 같이 이용하여 직접적으로 유영자세각의 관측이 필요할 것으로 사료된다.

현장에서 보름달물해파리의 에코신호를 식별하기 위해서는 주파수 차이를 파악하는 것이 필요한데, 본 연구에서는 보름달물해파리의 ∆MVBS _120-38kHz는 -1.5에서 2.2 dB의 차이를 적용하였지만, Mutlu (1996)의 연구에서 Aurelia autrans가 38 kHz에서는 -54.0 dB (8 cm), -51.7 dB (16 cm)이었고, 120kHz에서는 -54.2 dB (8 cm), -50.1 dB (16 cm) 였으며, ∆TS < sub>120-38kHz는-0.2에서 1.6 dB의 차이로 본 연구 결과와 비슷한 범위를 나타나는 것을 확인하였다. 해양에는 다양한 해양생물들이 서식하고 있는데, Miyashita et al. (2004)의 연구에서 부레가 있는 명태의 평균체적후방산란강도 120-38 kHz가 2.0 dB보다 낮은 값을 나타내었고, Kang et al. (2002)은 -1.0~4.0 dB의 범위로 보고되었다. Murase et al. (2009)의 연구에서는 멸치와 까나리의 성어가 -1.6~6.1 dB으로 나타났고, 멸치, 전갱이 등의 자치어가 우점했을 때에는 -9.6~-4.0 dB, 요각류가 13.7-17.3 dB으로 주파수 차이가 크게 나타났다. Safruddin et al. (2013)의 연구에서는 부레가 없는 까나리의 경우에는 -1.8 dB 으로 나타났으며, Kang et al. (2002)의 연구에서는 크릴의 경우, 10 dB 이상으로 주파수 38과 120 kHz를 이용하면 해파리와 동물플랑크톤은 명확히 분리되지만, 어류와는 분리가 되지 않는다. 본 연구에서도 주파수 차이를 적용한 후에는 보름달물해파리와 어군의 에코신호가 존재하여 어군의 에코를 어군 탐지 기능으로 선택하여 제거한 후에 보름달물해파리의 에코신호를 추출하였다. Colombo et al. (2009)의 연구에서 Aurelia는 주간에는 어군과 분리되어 분포하나, 야간에는 해파리와 어류가 혼재하는 특성을 보였다고 보고하였다. 따라서, 낮에는 주파수 차이를 적용하고 어군 에코를 제거하면 어군을 형성하는 에코신호는 제거될 수 있지만, 개체어로 분포하는 어류의 에코신호는 제거되기 어려울 수도 있을 것으로 판단된다. 따라서, 향후 보름달물해파리와 어류를 명확하게 분리할 수 있도록 추가적으로 다중 주파수를 이용하여 조사가 이루어져야 할 것으로 판단된다.

본 연구는 주간에 조사가 수행 되었고, 보름달물해파리는 저층보다는 표층부근에 많이 분포하였다. Malej et al. (2007)은 해파리가 황혼시에는 표층으로 올라가고, 야간에는 깊은 곳으로 하강하는 수직이동을 한다고 보고하였다. 또한, 해파리는 빛에 민감하고 맑은날에는 표층으로 이동하였다(Albert, 2010). 이와 같이 시간대별로 해파리의 이동 수심의 변화를 나타내기 때문에 시간대별 해파리의 분포수심을 파악하는 것도 필요하다. 보름달물해파리의 분포는 태양의 위치(Hamner et al., 1994) 해류, 조류 등의 흐름과 관련이 있는 것으로 알려져 있다(Graham et al., 2003; Rakow and Graham, 2006). 뿐만 아니라 해파리의 분포는 먹이생물과는 상관관계가 있었고, 클로로필, 수온, 염분과는 상관관계가 나타나지 않았다고 보고한바 있다(Han and Uye, 2009). 추후에는 먹이생물과 해양환경 조사를 추가하여 보름달물해파리의 분포를 파악할 필요가 있을 것으로 사료된다. 또한, 음향 기법은 저주파 트랜스듀서의 근거리 음장으로 인하여 수심 3 m 정도까지는 대상생물의 탐지가 어렵지만, 보름달물해파리는 목시 관측 시 표층에서도 많이 분포하므로 음향으로 탐지되지 않은 영역은 조사선박에서 카메라 기법이나 네트 채집을 이용하여 정도 높은 분포 밀도의 추정이 가능할 것으로 판단된다.

사 사

이 논문은 2019년 해양수산부 재원으로 해양수산과학기술진흥원(해양생태계교란 및 유해해양생물 관리기술 개발, 20190518)의 지원을 받아 의해 수행되었으며, 본 논문을 사려깊게 검토하여 주신 심사위원님들과 편집위원님께 감사드립니다.