1. 서 론

어류가두리양식장 아래 퇴적물은 잉여 사료, 어류 배설물, 시설물에서 탈락한 부착 생물의 사체 등으로 인하여 상당한 양의 유기물질이 유입되며, 이러한 유기물 축적은 어류가두 리양식장 저층의 퇴적 환경을 악화시키는 주요 요인이 된다 (Gowen and Bradbury, 1987;Karakassis et al., 2000;Gao et al., 2005). 퇴적물 표층의 유기물은 저서동물의 주요한 먹이원이 지만, 과잉 축적되면 저서동물의 생체량과 다양성 감소, 기 회종의 우점 등 저서동물 군집 구조의 변화를 발생시킨다 (Karakassis et al., 2000;Gao et al., 2005;Tomassetti et al., 2009). 또한 한번 교란된 군집은 회복되기까지 여러 단계의 천이과 정을 거쳐야하기 때문에 긴 시간이 걸린다(Pereira et al., 2004, Pearson and Rosenberg, 1978).

퇴적층의 급격한 유기물 증가는 서식지의 지화학적 성질 을 변화시킴으로써 저서동물 군집 구조에 영향을 미치는데, 이때 군집의 종 구성은 민감종의 비율이 감소하고 내성이 강한 종 또는 오염지표종의 비율이 증가하는 것이 일반적이 다(Pearson and Rosenberg, 1978;Hargrave, 2008). 이러한 개념에 기초하여 생태계의 건강도를 측정하고자 하는 다수의 생태 지수들이 개발되어 왔으며, 유럽과 미국 등에서 사용 중인 AMBI(Azti's Marine Biotic Index), M-AMBI(Multivariate Azti's Marine Biotic Index), BENTHIX, BQI(Benthic Quality Index) 등 이 대표적 사례이다(Borja et al., 2000;Muxika et al., 2007;Rosenberg et al., 2004).

생태지수는 출현종의 생태학적 특성을 기반으로 만들어 졌기 때문에, 다양한 환경, 그리고 다양한 교란 반응에 사용 될 수 있다. 생태지수를 활용한 건강도 평가는 단일의 숫자 나 등급으로 건강도를 나타내기 복잡한 과학적 연구 자료의 정보를 보다 직관적인 방식으로 표현할 수 있다(Borja et al., 2008;Pinto et al., 2009). 따라서 모니터링, 건강도 평가 등에 서 정책결정자나 비과학자들이 이해하고 비교하기에 편리 하다는 장점이 있다. 그러나 지수 산출에 앞서 출현종의 생 태학적 특성을 파악하는 과정이 요구되며, 이미 개발된 국 외의 생태지수를 국내에 그대로 적용할 시에는 서식종의 차 이가 있어서 건강도 평가에 활용하기에 적합하지 않다.

우리나라에서는 우리 상황에 맞추어 저서동물지수(BHI)를 개발하였으며(국립수산과학원 고시 2020–15, 어장환경평가 의 방법 및 절차 등에 관한 규정), 이를 어류가두리양식장을 대상으로 하는 어장환경평가에 활용하고 있다.

본 연구에서는 어장환경평가를 통해 조사된 어류가두리 양식장을 대상으로 저서다모류 군집 구조를 파악하고, 저서 동물지수를 활용하여 저서생태계의 건강도를 평가하고자 한다.

2. 재료 및 방법

2.1 시료 채집 및 분석

본 연구에서는 2015년부터 2021년까지 어장환경평가가 이 루어진 어류가두리양식장에 대한 조사 결과 중 일부를 포함 하였다. 각 정점이 위치한 지역 및 수심, 그리고 현장에서 기록한 표층 퇴적물의 특징을 Table 1에 정리하였다. 또한 현장에서 양식장 운영 기간, 양식 어종 등을 기록하여 분석 에 참고하였다.

저서다모류 군집과 총유기탄소(TOC: Total organic carbon) 분석을 위해 0.05 ㎡의 van Veen grab을 사용하여 퇴적물을 채집하였다. 다수의 가두리가 설치된 정점에서는 하나의 가 두리에서 2회씩 채집하여 1회분은 저서다모류 분석을 위한 시료로, 1회분은 표층 총유기탄소 분석을 위한 시료로 사용 하였다. 단일 가두리가 설치된 정점에서는 4회 채집하여 2회 는 저서다모류 분석, 2회는 총유기탄소 분석 시료로 사용하 였다. 저서다모류 분석 시료는 선상에서 1 mm 체에 걸러 퇴 적물을 최대한 제거한 후에 10 %의 중성 포르말린으로 고정 하였으며, 총유기탄소 분석 시료는 채집과 동시에 드라이아 이스에 넣어 냉동하였다.

현장에서 채집한 시료 중 저서동물 시료는 실험실에서 분 류군별로 선별한 이후 다모류를 대상으로 종 수준까지 동정 하였다. 총유기탄소 시료는 실험실에서 동결건조 후 1 N의 HCl을 첨가해 탄산염을 제거하였다. 70℃에서 24시간 동안 수분과 잔류 염소를 날린 이후 105℃에서 2시간 건조과정을 거친 다음 곱게 간 시료를 주석 박막으로 말아 CHN 원소분 석기(Perkin Elmer, model 2400)로 측정하였다.

2.2 자료 분석

정점 별 개체수 자료는 단위면적으로 환산하였으며, 개별 군집의 기본구조 파악을 위한 단변량 변수로 종수, 개체밀도, 종다양도지수(H’)를 산출하였다. 전체 군집의 구조 파악을 위 한 다변량 분석으로 다차원척도법(non-metric multidimensional scaling: MDS)을 적용하였으며, 유사 정점군 간의 그룹화에 임의성을 배제하기 위하여 집괴분석(Cluster analysis)의 결과 를 반영하였다(Clarke and Warwick, 2001;Anderson et al., 2008). 다변량 분석에 앞서 생물자료는 로그 변환하였으며, Bray-Curtis 유사도를 적용하였다. 다변량 분석으로 구분된 각 그룹의 특징종을 찾아내기 위하여 SIMPER 분석을 수행 하였고, 주좌표분석(Principal Coordinates analysis, PCO)을 통해 특징종과 군집구조와의 상관성을 알아보았다. 모든 통계분 석은 PRIMER 패키지 프로그램(ver. 6)을 이용하였다.

어류가두리양식장 건강도 평가는 국립수산과학원에서 개 발한 저서동물지수(Benthic Health Index: BHI)를 사용하였으 며, 계산식은 다음과 같다.

BHI = 25((4×N₁ + 2.68×N₂ + 1.36×N₃ + 0.04×N₄) / N_total)

Ni(i=1−4)는 각 그룹에서 출현한 개체수의 총합을 나타내 며, Ntotal은 전체 출현개체수이다. 유기물 농도와 관련된 각 그룹종의 생태적 특징은 그룹 1이 낮은 유기물 농도에서 높 은 밀도로 출현하거나 출현 빈도와 밀도가 낮은 종(이하 1그 룹종), 그룹 2는 유기물의 농도와 상관없이 고른 분포를 하 는 종(2그룹종), 그룹 3은 비교적 높은 유기물 농도에서 높은 밀도로 출현하는 종(3그룹종), 그룹 4는 계절적으로 무생물 역이 발생하는 해역에서 높은 밀도로 출현하거나 높은 유기 물 농도에서 출현하는 종(4그룹종)을 의미한다. 각 종에 대 한 분류는 국립수산과학원 고시에서 제공하는 종을 참고하 였다. BHI는 0−100 범위에서 산출되며, 0−25는 4등급으로 오염심화 해역(Heavily polluted, H/P), 26−50는 3등급으로 오 염 해역(Moderately polluted, M/P), 51−70은 2등급으로 오염초 기 해역(Slightly polluted, S/P), 71−100는 1등급으로 정상 해 역(Normal, N)을 의미한다(국립수산과학원 고시 2020-15, 어 장환경평가의 방법 및 절차 등에 관한 규정).

3. 결 과

3.1 어류가두리 양식장의 총유기탄소량

국내 어류가두리양식장의 총유기탄소량은 0.80−34.92 mg/g 의 범위였다(Fig. 1a). 총유기탄소량이 10.00 mg/g 이하인 양식 장은 st.01−03, st.32, st.35−39로, 강원도와 경북, 전남 완도 와 신안 그리고 충남 보령의 정점이었다. 어장환경평가에서 유기물 오염이 심각한 퇴적 상태를 나타내는 총유기탄소량 25.01 mg/g 이상인 가두리는 st.12−13, st.20−21, st. 24−26, st.28−29로 경상남도의 정점들이었다. 국내 어류가두리양식 장 표층퇴적물의 총유기탄소량은 남해안의 가두리에서 상대 적으로 농도가 높았고, 동해와 서해 일부 가두리에서 낮았다.

3.2 다모류 군집조성

국내 연안 어류가두리에서 출현한 다모류는 총 224종이 었으며 개체밀도는 평균 8,300 inds./㎡였다. 각 정점당 출현 종수는 6−35 spp./0.05 ㎡ 범위였으며, st.42에서 가장 적었고, st.17에서 가장 많았다(Fig 1b). 출현한 다모류가 10종 미만인 정점은 st.01−03, st.25, st.29, st.34. st.42−43이었다. 개체밀도 는 590−21,735 inds./㎡의 범위였으며, st.01−04, st.08, st.11, st.32, st.36, st.39 정점에서 2,000 inds./㎡ 이하로 출현하여 상 대적으로 밀도가 낮고, st.24, st.26, st.28−29에서 15,000 inds./㎡ 로 밀도가 높았다(Fig 1c). 종다양성지수(H’)는 0.2−3.0 범위 였으며, st.03, st.07, st.25−26, st.28−29, st.31, st.33−34, st.40, st.42−43에서 1.0 이하로 낮은 값을 나타내었다(Fig. 1d).

43개 정점에서 유사한 군집 특성을 가지는 정점군(Group) 을 구분하여 출현한 어류가두리양식장의 다모류 우점종의 경향을 분석하였다. Group의 구분은 개체밀도를 대상으로 집괴분석을 수행하였으며, 분석 결과 크게 5개 그룹(Group A −E)으로 구분되었고, 그룹 E는 다시 2개의 소그룹(Group E1, E2)로 구분되어 총 6개 그룹으로 구분하였다(Fig. 2).

그룹 A에서는 저서동물지수 계산 시 1그룹종인 Spiophanes bombyx(863 inds./㎡, 77.5 %)가 높은 밀도로 출현하였고, 2그 룹종인 Chone teres와 Glycinde gurjanovae가 우점종으로 출현 하였다(Table 2). 그룹 B는 2그룹종인 Magelona japonica(385 inds./㎡, 29.6 %)가 최우점하였으며, 1그룹종인 Mediomastus sp. 와 Heteromastus sp.가 우점하였다. 그룹 C에서는 1그룹종인 Cirrophorus armatus(433 inds./㎡, 19.6 %), Euclymene oerstedii(400 inds./㎡, 18.1 %), Mediomastus sp.(390 inds./㎡, 17.6 %)가 유사한 비 율로 우점하였다. 4번째 우점종은 3그룹종인 Scoletoma longifolia 로 300 inds./㎡(13.7 %)가 출현하였다. Group D, E1, E2에서는 상위 출현종이 유사하였으며, 이들종의 우점순위와 비율에 차이가 있었다. 최우점종은 4그룹종인 Capitella capitata로 동 일하였으며, Group D에서는 627 inds./㎡(39.7 %), Group E1에 서는 1,598 inds./㎡(25.0 %), Group E2에서는 9,074 inds./㎡(78.5 %)가 출현하였다. 차순위 우점종은 Group D에서는 S. longifolia, Sigambra tentacultata가 우점하였고 Group E1에서 는 Schistomeringos rudolphi, S. longifolia, Group E2에서는 S. rudolphi, S. tentacultata였다.

3.3 어류가두리 건강도 평가 (BHI, AMBI 비교)

연구지역의 저서동물지수(BHI)는 2−97 범위로 국내 어류 가두리양시장의 건강도는 정상상태(Normal)에서 오염심화해 역(Heavily polluted) 이었다(Fig. 3a). 어류건강도가 정상상태인 1등급 정점은 st.01−03, st.08, st.11, st.37−39 등 8개였으며, 2 등급인 오염초기해역의 상태(Slightly polluted)인 정점은 st.04, st.16, st.32, st.36로 4개 정점이었다. 3등급인 오염해역 (Moderately polluted)의 건강도를 나타내는 정점은 st.05−06, st.09−10, st.14−15, st.17−18, st.22, st.27, st.30, st.35의 12개 정점이었으며, 나머지 19개 정점은 오염심화해역이었다.

어류양식장에서 BHI를 활용한 건강도의 비교를 위해, AMBI를 산출하여 비교하였다(Fig. 3b). AMBI는 0.9−5.9의 범 위로 정상상태에서 오염심화해역의 건강도를 나타내었다. BHI와 비교하여 23개 정점은 같은 등급의 건강도를 나타내 었고, 13개 정점에서 BHI의 건강도가 높게 나타났고, 7개 정 점에서 AMBI의 건강도가 높게 나타났다. AMBI의 건강도가 낮게 나타난 정점은 st.01−03, st.08, st.11, st.38−39로 모두 BHI는 정상상태를 보였고, AMBI는 약간오염된 상태를 나타 내었다.

3.4 저서다모류 군집 구조

집괴분석으로 구분한 6개 Group을 MDS 분석을 통해 2차 원 평면상에 구현하고 각 그룹의 상대적 위치와 그룹별 특 징을 비교해 보았다. 먼저 Group A로 묶인 정점군은 평면도 의 가장 우측 하단에 위치하였으며, 다른 정점군과 평면상 좌우로 크게 구분되었다. 6개 그룹 중 출현 종수, 개체밀도 그리고 총유기탄소량이 가장 낮았으며, 종다양도 또한 낮았 다(Table 3, Fig. 4a−c). 출현종은 낮은 유기물 농도에서 서식 하는 종(1그룹종, 2그룹종)이 대부분을 차지하였다(Fig. 4d). Group B는 평면상의 좌측 상단 위치하였으며, 다른 그룹에 비해 개체밀도 대비 출현 종수는 많아 다양도가 높은 것이 특징이었으며, 총유기탄소량도 2등급 수준으로 높은 편이었 다. 종조성은 1, 2그룹종이 출현종의 대부분을 차지하였으나 Group A와 비교하여 2그룹종의 비율이 상대적으로 높았다. Group C와 Group D는 평면상의 좌측 중앙에 위치하였다. 두 그룹의 출현 종수와 개체밀도는 6개 그룹 중 중간 정도 수 준으로 상대적으로 높은 다양도를 보였지만, 총유기탄소량 에서는 차이가 있었다. 앞선 두 Group과 비교하였을 때 3, 4 그룹종의 다모류 비율이 높아졌으며, Group D는 3, 4그룹종 의 비율이 50 % 이상을 넘었다. 다른 그룹에 비해 많은 정점 을 포함하는 Group E1과 Group E2는 평면상에서 좌측 하단 에 위치하였다. 두 그룹은 6개 그룹 중 가장 높은 총유기탄 소량을 보이는 공통점이 있었지만 Group E1는 가장 많은 출 현종수를, Group E2는 가장 높은 개체밀도를 보여 다양도에 큰 차이를 보이는 점이 특징적이었다. 출현종은 3, 4종의 비 율이 높았으며, Group E2의 경우에는 4등급종이 86.6 %로 높 은 비율을 차지하였다.

각 Group의 BHI 점수의 범위를 살펴보면, Group A−C는 건강도가 1등급의 정상상태였으며, Group D는 건강도가 2− 3등급이었고, Group E1은 건강도가 3−4등급, Group E2는 4 등급으로, Group E1, E2의 건강도가 악화된 상태였다(Fig. 5). 각 그룹이 대별되는 특징을 종조성 관점에서 알아보기 위해 SIMPER 분석을 수행하였다. 종조성의 차이를 단순 우점종 기준으로 비교할 경우 그룹 간 겹치는 종들의 상대적 중요 도 파악이 어렵기 때문에 그룹 내 그리고 그룹 간의 차이에 있어 상대적 중요도를 개별종 차원에서 제시하는 SIMPER를 적용하였다(Table 4). 먼저, Group A의 경우 대표종은 그룹 내 유사도와 그룹 간 비유사도에 모두 기여한 Spiophanes bombyx 였으며, 특히, 그룹 내 기여율이 90.96으로 매우 높았다. Group B는 그룹 내 기여종인 Prionospio saccifera, Ampharete arctica와 그룹간 비유사도 기여종인 Magelona japonica가 대표종이었다. Group C의 대표종은 그룹 내 기여종인 Mediomastus sp., Prionospio multibranchiata와, 그룹 간 비유사도 기여종인 Cirrophorus armatus, Euclymene oerstedi였다. Group D는 그룹 내 유사도 기여종인 Nephtys oligobranchia와 Heteromastus sp.가 대표종이었다. 그룹 간 비유사도 기여율에서 대별되는 종은 없었으며, 최우점종 인 Capitella capitata도 Group E2에서 밀도가 더 높아 Group D 의 대표종에서 제외하였다. Group E1의 대표종은 그룹 내, 그 룹 간 구분에 모두 기여한 S. longifolia, Schistomeringos rudolphi 와 그룹 간 비유사도에 기여한 Aphelochaeta monilaris이었다. Group E2에서는 그룹 내 기여종인 C. capitata와 그룹간 비유 사도 기여종인 Sigambra tentaculata가 대표종이었다.

각 그룹의 군집 특성을 나타내는 대표종들을 알아보기 위 한 또 하나의 방법으로 주좌표분석(PCO)을 수행하였다(Fig. 6). 그룹의 형성은 다차원척도법과 거의 일치하였으며, 각 그룹의 배치와 각 그룹 형성에 영향을 미친 대표종 또한 SIMPER 분석 결과와 일치하였다.

4. 고 찰

연구 지역의 43개 어류가두리양식장을 대상으로 양식장 의 건강도를 평가한 결과 1등급 양식장은 8개, 2등급 양식장 은 4개, 3등급 양식장은 12개, 4등급 양식장은 19개로 전반적 으로 3, 4등급이 우세하였다. 종조성의 유사도를 기준으로 군집구조를 분석해보면 6개의 그룹으로 구분되었으며, Group A−C의 건강도는 1등급이었고, Group D는 2등급이 우세하였 고, Group E1은 3등급이 우세하였으며, Group E2는 4등급의 건강도였다.

1등급의 건강도를 나타내는 Group A−C의 양식장은 사질 퇴적상(Group A), 사질함량이 높고 분급도가 불량한 퇴적상 (Group C) 그리고 양식 강도가 낮은(Group B) 양식장이 속하 였다.

퇴적물 입도는 종 분포를 결정하는 요인 중 하나이며, 일 부 종은 퇴적물 유형과 상관성이 적으나 일부 종은 특정 퇴 적상을 선호한다(Dauvin et al., 2004). Group A의 최우점종이 자 기여종인 Spiophanes bombyx는 사질 퇴적상을 선호하는 동해안의 대표종이며(Choi and Koh, 1988;Paik et al., 2007), Group C의 우점종인 Cirrophorus armatus, Euclymene oerstedi 역시 국내 펄, 자갈이 혼합된 양식장이나 분급이 불량하고 사질 함량이 높은 퇴적상에서 우점하는 종이다(Koo et al., 2005;An et al., 2006;Jung et al., 2014). 일반적으로 유기물은 세립질 퇴적물과 친화력이 크기 때문에 입도가 세립할수록 유기물 함량이 증가하는 경향이 있다(Horowitz, 1991). 다른 Group과 비교해서 Group A와 C의 낮은 TOC는 사질 특성의 퇴적상의 영향으로 판단되며, 낮은 유기물 농도에서 1그룹 종과 2그룹종이 우점하여 정상해역(Normal)의 건강도를 나타 내었다. Group A와 C를 비교하면, 퇴적물의 이질성이 증가할 수록 생물 다양성이나 개체수가 높다는 연구 결과(Seo, 2003) 와 유사하게 분급도가 불량한 Group C에서 군집의 출현 종 수, 개체밀도, 종다양도가 상대적으로 높았다. Group B는 현 장조사 당시 전년도에 어류를 출하하고 입식을 하지 않은 양식장과 1년간 휴식중인 양식장을 포함하여 다른 그룹과 대별되었다. 이러한 이유로 동일 해역에 위치한 다른 어류 가두리(Group E1, E2)와 비교하여 유기물 농도가 낮고 건강 도가 높았다.

건강도가 2등급인 Group D의 양식장은 육지에서 멀리 이 격된 해역, 양식장 옆으로 수로가 위치한 해역, 만의 입구에 위치한 해역 등 해류의 흐름이 매우 좋은 해역에 위치하였 다. 양식 활동의 영향범위와 강도는 양식장이 설치된 지역 의 수력학적인 조건에 영향을 받는 것으로 알려져 있으며 (Aure and Stigerbrant, 1989;Hevia et al., 1996;Villnäs et al., 2011), 조류나 해류가 강한 해역에서는 양식장 아래에 미치 는 영향이 감소되는 것으로 보고되었다(Frid and Mercer, 1989;Kutti et al., 2008). Group D는 Group E1, E2와 비교하여 유사한 양식강도를 갖지만 어장 주변의 빠른 유속의 영향으 로 유기물이 분산됨으로써 가두리 바로 밑의 유기물 농도는 Group E1, Group E2보다는 많이 낮았다. 그러나 종조성 측면 에서는 유사한 유기물 농도의 Group B보다는 Group E1과 E2 에 가까웠다.

Group E1과 E2는 평면상에서 가까운 거리에 위치하여 유 사도가 높았고, 높은 TOC를 나타내어, 유기물 오염이 진행 된 상태였으나, 건강도는 Group E1은 3등급, E2는 4등급으로 차이가 있었다. Pearson and Rosenberg(1978)은 유기물 농도에 따른 대형저서동물 군집의 천이모델을 제시하였는데, 이 연 구에서는 유기물이 증가함에 따라 종수와 개체밀도가 증가 하는 Transition region, 종수가 감소하고 개체밀도가 급격히 증가하는 Ecotone point, 상대적으로 작은 신체크기를 가지는 소수의 기회주의종이 극우점하여 개체밀도가 최고점인 Peak of opportunists의 상태, 그리고 종수와 개체밀도 모두가 감소 하는 빈약한 군집상태로 천이된다고 제시하였다. Group E1 과 E2는 높은 유기물의 농축으로 Transition region은 이미 넘 어섰으며, Group E1은 Ecotone point 도달하기 이전의 군집 상 태였으며, Group E2는 Peak of opportunists의 상태의 군집으로 판단된다.

Group E와 E2의 최우점종은 Capitella capitata였고, 차순위 우점종은 Schistomeringos rudolphi로 동일하였으나, Group E1 에서는 C. capitata와 S. rudolphi의 밀도 차이가 크지 않았고, Group E2에서 C. capitata가 극우점하는 양상을 나타내었다. C. capitata와 S. rudolphi는 4그룹종이지만, C. capitata는 유기 물에 내성이 있어 높은 유기물 농도에서 고밀도로 출현하지 만 종간 경쟁이 다소 약하며, S. rudolphi는 C. capitata에 비해 오염에 대한 내성은 상대적으로 낮으나 종간 경쟁에서 초기 가입하여 번식하는 차이가 있다(Pearson and Rosenberg, 1978;Dean, 2008). 따라서 두 종의 분포특성에 따르면 C. capitata가 극우점한 Group E2의 유기물 오염이 더 심각한 상태로 판단 된다. 유기물 농도가 일정 농도 이상에서는 TOC만으로 두 Group의 건강도 차이가 구별하기는 어려웠으며, 두 Group의 다모류 군집 조성과 우점종의 경향을 참고하면 Group E2의 유기물 오염이 더 심각한 상태로 판단된다.

BHI와 AMBI의 건강도는 전체 43개의 어류가두리에서 23 개가 같은 등급으로 나타나 53.5 %가 일치함을 보여주었다. 두 지수의 건강도 차이는 지수 계산시 할당된 종그룹이 다 르기 때문이다. AMBI와 비교하여 BHI의 건강도가 높은 정 점은 Group A와 C였는데, 이들 Group의 우점종인 Spiophanes bombyx와 Cirrophorus armatus가 BHI에서는 1등급종에, AMBI 에서는 3등급종에 할당되어 BHI의 건강도가 높게 나타났다. AMBI의 건강도가 높은 정점은 Group E1과 E2의 일부 정점 이었으며 Group E1과 E2의 상위 우점종인 Scoletoma longifolia 가 BHI에서는 3등급종에, AMBI에서는 2등급종에 해당하여 두 지수의 등급차이에 기여하였다. 연구지역에서 두 건강도 지수의 등급과 TOC 등급을 상호비교해보면, 유기물 농도가 낮은 곳에서는 BHI가 TOC와 등급의 차이가 적어 적합한 것 으로 판단되며, 유기물 농도가 높은 곳에서는 AMBI, BHI 와 TOC간의 직접적인 비교는 어려웠다. 다만 BHI가 국내 종을 대상으로 유기물농도와의 상관성에 기초하여 만들어 진 지수이기 때문에, 국내의 양식장에 더 적합한 지수로 판단된다.

국내 어류가두리 양식장의 건강도 평가 결과, 1등급 양 식장은 저층 퇴적상의 영향으로 유기물 축적이 적은 지역과 특정 상황에 의해 양식강도가 낮은 지역이 속하였으며, 2등 급은 해수유통이 원활한 지역이 포함되었다. 그러나 어장환 경평가로 조사된 대부분의 어류가두리양식장의 건강도는 3, 4등급으로 나타나 교란된 상태였다. 국내 연안 및 다른 품종 의 양식장에서 BHI를 활용한 연구 결과와 비교해보면(Table 5), 어류가두리양식장의 건강도는 조간대와 양식장이 없는 해역에 비해 건강도가 낮았으며(Ong et al., 2021;Youn et al., 2021), 패류양식장과 비교하여 유사하거나 낮았다(Park et al. 2017, 2021a). Han et al.(2020)의 연구에서는 국내 환경영향평 가 해역에서 유기물 오염 및 건강도를 평가하였는데, 남해 안이 동해안과 서해안에 비해 상대적으로 건강도가 낮았다. 이는 남해안의 복잡한 해저지형으로 인한 해수유통의 방해 로 상대적으로 유기물 체류가 길어졌기 때문으로 보았다. 본 연구에서도 남해안의 어류양식장 건강도가 상대적으로 낮은 경향이었으나, 양식장의 대부분이 남해안에 위치하여 직접적인 공간 비교는 어려웠으며, 같은 해역에 위치하여도 저층 퇴적상, 양식의 강도, 해수의 수리학적인 요인 등 여러 요인이 어류가두리양식장의 건강도에 영향을 주는 것으로 나타났다.