1. 서 론

어류가두리양식장에서 배출된 어류 배설물과 잉여 사료 등의 양식 부산물은 양식장 바로 아래 퇴적물에 높은 유기 물 농축을 일으킨다(Black, 2001;Gowen and Bradbury, 1987;Hall et al., 1990;Macleod et al., 2006). 양식장 아래에 서식하 는 대형저서동물군집은 환경에 민감성, 정착 습성, 약한 이 동성으로 인해 유기물 농축의 영향을 직접 받으며(Pearson and Rosenberg, 1978), 생체량과 다양성의 감소, 기회종의 우 점, 군집 구조의 변화 등의 결과를 나타낸다(Karakassis et al., 2000;Gao et al., 2005;Tomassetti et al., 2009).

양식 활동으로 인한 부정적인 영향을 줄이고 양식장 아래 퇴적물의 회복을 돕기 위하여 일정기간 동안 양식 활동을 중단하는 방법을 사용하며, 이것은 퇴적물의 점진적인 악화 를 줄이고 장기간의 양식 활동을 가능하게 한다(Black et al., 2008;Beveridge, 1987;Pereira et al., 2004). 온대 해양의 주요 양식어종인 연어의 경우, 전적으로 가두리에서 양식하며 양식 장 아래에 유기물을 대량 방출하였다(Beveridge, 1987;Pereira et al., 2004;Lee et al., 2006;Kutti et al., 2007;Huang et al., 2012). 이에 양식이 종료되면 일정 기간의 휴지기를 주는 양식 방법을 택하고 있다(Aguado-Gimenez et al., 2012;Black et al., 2008;Macleod et al., 2006). 국내 어류가두리양식장은 어장관 리법 제11조의2에 따라 어업면허 또는 허가의 유효기간이 끝나기 1년 전까지 어장환경평가를 실시하며, 어장환경평가 결과는 어장관리법 시행규칙 제3조의3에 따라 환경개선 행 정조치를 받는다. 어장환경평가 결과 3등급과 4등급의 어장 은 면허 연장 시 면허의 유효기간이 감소하며, 어장환경 개 선 및 어장시설물 위치 이동 지시를 받는다. 국내에서는 양 식장 휴지와 관련된 법은 없으나, 어장 시설물 이동을 통해 양식장 바로 아래 퇴적물은 양식장에 휴지를 설정한 상태와 유사한 결과를 얻을 수 있다.

가두리양식장에서 휴지기 설정이나 양식시설물 제거 이 후 시공간적인 변화에 대한 연구는 국외에서 많이 진행되어 왔으며(Liao et al., 2022;Macleod et al., 2006), 양식장 저층 환 경이 회복되는 속도와 영향범위는 양식 품종, 양식장의 위 치, 양식 주기, 양식장의 운영기간 및 강도, 퇴적물 특성 등 에 영향을 받았다(Brooks et al., 2004;Macleod et al., 2006;Norkko et al., 2006;Zhulay et al., 2015). 양식장 휴지기 이후 회복은 생물학적인 회복 속도와 화학적인 회복 속도에 차이 가 있었으며 생물학적인 회복속도가 상대적으로 오래 걸렸 다. 그러나 일단 양식이 중단되면 퇴적물이 회복되는 것은 분명하였다(Brooks et al., 2004;Johannessen et al, 1994;Karakassis et al., 1999;Pereira et al., 2004;Pohle et al., 2001). 국내 가두리 양식장에서는 양식장 휴지와 관련된 연구가 부족하다. 패류양식장에서 휴지기동안 퇴적환경의 회복에 관한 연구가 수행된 바 있으나(Park et al., 2021;Sim et al., 2020) 패류양식장에 비해 교란이 더욱 심한 것으로 알려져 있는 어류양식장의 연구는 없었다.

본 연구에서는 양식장 이동 이후에 저서다모류 군집의 모 니터링을 수행하였으며, 다모류 군집을 이용하여 어장환경 평가 등급을 산정하였다. 이를 통해 양식장 이동 이후 퇴적 환경의 생물학적인 회복을 살펴보았다.

2. 재료 및 방법

2.1 연구 지역

본 연구지역은 경상남도 하동군 금남면 중평항 근처에 위치한 어류가두리양식장으로 2006년에 양식을 시작하였고, 2015년 국립수산과학원에서 수행하는 어장환경평가에서 3등 급을 받은 후 어장환경 개선 및 시설물 위치 이동 지시의 행 정조치를 받았다. 2017년 10월 ~ 11월 사이에 기존 양식장에서 약 150 m 이격된 위치로 모든 시설물을 이동하였다(Fig. 1). 양식장 철거로 휴지기가 주어진 기존 정점에 3개 정점 (Farm1 ~ 3)을 선정하였으며, 새롭게 설치된 양식장과 가까운 정점을 Farm1, 중간 거리 정점을 Farm2, 먼 정점을 Farm3으로 하였다. Farm1과 양식장과의 거리는 약 80 m 정도였으며, Farm3과 양식장과는 약 150 m 이격되어 있다. 기존 정점들 과의 비교를 위하여 양식시설물이 없는 3개 정점(Con1 ~ 3) 을 대조 정점으로 선정하였다. 양식장 이동 이전인 2017년 10월 사전 조사로 첫 조사를 수행하였고, 이후 2018년 10월 까지는 매달, 그 이후에는 2 ~ 3개월 간격으로 현장조사를 수행하였다.

2.2 퇴적물 유기물 분석 및 저서동물 분석

퇴적물 시료는 표면적 0.05 ㎡의 van Veen grab을 이용하여 정점 당 3회씩 반복 채집하였으며, 1회분은 표층퇴적물의 화 학적 특성 분석용으로, 2회분은 저서다모류 군집 분석용으 로 사용하였다. 퇴적물 특성과 관련하여 총유기탄소량(TOC), 총질소량(TN), 산휘발성황화물(AVS), 강열감량(IL), 화학적산 소요구량(COD), 산화환원전위(ORP)를 측정하였으며, 각 변 수들에 대한 분석은 해양환경공정시험기준(MOF, 2013)에 준 하여 측정하였다. 퇴적물과 인접한 저층에서 CTD(SBE19, Seabird, USA)를 이용하여 수온(Temp), 염분(Sal)을 측정하였 고, 다목적 수질측정기(YSI6000, YSI, USA)를 이용하여 용존 산소(DO),수소이온농도(pH)을 측정하였다. 분석된 항목 중 TOC는 저서동물자료와 함께 저서생태계 건강도 평가기법에 적용하였고, TOC를 포함한 모든 측정 항목은 저서다모류 군 집구조와의 상관성 분석에 활용하였다.

2.3 자료 분석

각 해역에서 분류와 종 동정 과정을 통해 얻어진 다모류 자료는 군집 특성을 파악하기 위하여 단변량 분석과 다변량 분석을 수행하였고, 저서생태계 건강도 등급을 산출하였다. 단변량 분석을 위한 생태지수로서 출현종수, 개체밀도, 그리 고 다양도(Shannon and Weaver, 1963)를 활용하였다. 저서다모 류 군집구조 파악을 위해 출현종의 현존량 자료를 로그 변 환 후 Bray-Curtis 유사도 지수를 구하고, 유사도 행렬에 기초 하여 시·공간 변화에 따른 저서동물군집의 구조적 차이를 살펴보았다. 환경변수와 저서다모류 군집구조와의 상관성을 알아보기 위해 주좌표분석(Principal Coordinates analysis, PCO) 을 수행하였으며(Gower, 1966), 모든 통계 분석은 Primer 패키 지(ver. 6)를 이용하였다.

저서생태계 건강도 평가는 국립수산과학원에서 개발한 저서동물지수(Benthic Health Index, BHI)와 총유기탄소량을 활용한 어장환경평기기법(국립수산과학원 고시 제 2020-15 호)을 사용하였다. 저서동물지수의 계산식은 아래와 같다.

여기서, N은 각 그룹의 개체수를 나타내며, N1은 낮은 유 기물 농도에서 높은 밀도로 출현하거나 출현빈도와 밀도가 낮은 종, N2는 유기물 농도와 상관없이 고른 분포를 하는 종, N3은 비교적 높은 유기물농도에서 높은 밀도로 출현하 는 종, N4는 계절적으로 무생물역이 발생하는 해역에서 높 은 밀도로 출현하거나 높은 유기물 농도에서 출현하는 종이 속한다.

어장환경평가 등급의 점수화 과정은 Table 1과 같으며, 1 등급은 유기물 오염이 없는 상태, 2등급은 유기물 오염의 초 기 상태, 3등급은 유기물 오염이 상당히 진행된 상태, 4등급 은 유기물 오염이 심각한 상태를 나타낸다. 어장환경평가와 관련된 구체적인 절차와 방법은 국립수산과학원 고시 제 2020-15의 방법에 준하였다(NIFS, 2020).

3. 결과 및 고찰

3.1 총유기탄소량 변화

조사 기간 동안 총유기탄소량농도(TOC)는 10.80 ~ 28.28 mg․g^-1 dry weight였고, 정점별로 살펴보면 기존 정점은 16.18 ~ 28.28 mg․g^-1 dry weight, 대조 정점은 10.80 ~ 17.80 mg․ g^-1 범위였다(Fig. 2). 전체 조사 기간 동안 대조 정점은 총유 기탄소량의 큰 변화가 없으나, 기존 정점은 양식장 철거 이 후 0.7 mg/month(R²=0.54)로 점차 감소하는 경향을 나타냈다. 양식장 철거 다음 해인 2018년 4월 조사에서 기존 정점의 총 유기탄소농도가 일시적으로 증가하는 경향을 보였으나 다 시 감소하였으며, 어장 철거 약 1년 이후에는 대조정점과 통 계적으로 유의한 차이를 보이지 않았다(p<0.05).

3.2 저서다모류 군집의 생태지수 변화

조사 기간 동안 기존 3개 정점에서 총 76 종, 월평균 8 종 /0.1㎡의 다모류가 출현하였고, 대조 3개 정점에서는 총 137 종, 월평균 24 종/0.1㎡이 출현하여 기존과 대조 정점 간 차 이를 보였다(Fig 3a). 기존 정점의 출현종수는 Farm1은 0 ~ 18 종/0.1㎡, 평균 6 종/0.1㎡, Farm2는 0 ~ 20 종/0.1㎡, 평균 6종 /0.1㎡, Farm3은 0 ~ 22 종/0.1㎡, 평균 11 종/0.1㎡이 출현하여, 양식장과 떨어진 정점에서 출현종수가 많았다. 기존 정점에 서 월평균 출현종수는 2018년 여름까지 증가하였으나, 여름 부터 가을까지 급격히 감소하여 무생물 군집을 나타내었고, 이러한 양상은 2019년, 2020년에도 반복되었다(Fig. 3a). 대조 3개 정점의 출현종수는 Con1은 8 ~ 41 종/0.1㎡, 평균 23 종 /0.1㎡, Con2는 4 ~ 36 종/0.1㎡, 평균 23 종/0.1㎡, Con3은 10 ~ 45 종/0.1㎡, 평균 25 종/0.1㎡ 이었으며, 월평균 출현 종수는 세부적인 증가와 감소의 차는 있지만 봄철에 증가하고 여름 철 이후에 낮아지는 경향은 기존 정점들과 동일하였고 무생 물 출현 시기는 없었다.

기존 정점에서 출현 다모류의 개체밀도는 Farm1은 0 ~ 8,450 개체/㎡, 평균 1,040 개체/㎡, Farm2는 0 ~ 8,350 개체/㎡, 평균 1,081 개체/㎡, Farm3은 0 ~ 2,740 개체/㎡, 평균 752 개체 /㎡로 출현종수와 달리 양식장에서 거리가 먼 정점인 Farm3 에서 개체밀도가 낮았다(Fig. 3b). 개체 밀도의 시간에 따른 변화는 양식시설물과의 거리에 따라 차이를 보였는데, 양식 시설과 상대적으로 가까운 Farm1, Farm2에서는 2017년 10월 무생물 상태 이후 11월 개체밀도가 급격히 증가하여 대조정 점 이상 밀도를 보였으며, 2018년 2월과 3월에는 약 3배 가 까이 증가하였다. Farm3에서는 조사초기 개체수가 증가하였 으나 Farm1, Farm2와 같은 폭발적인 증가는 없었다. Farm1, Farm2에서 개체밀도의 증가는 어장환경평가 시 4그룹에 속 하는 Capitella capitata 단일 종에 의해 이루어 졌으며, 개체 밀도는 폭발적인 증가 이후 불과 한 달 만에 급격히 감소하 였다. 이후 기존 정점의 개체밀도는 점진적으로 감소하였고, 무생물 상태가 나타났다. 2018년 여름철 무생물 시기의 발생 과 유지 기간 역시 차이를 보여 Farm1은 7월부터 10월까지 였으나 Farm3은 9월에만 무생물 상태였다. 기존 3개 정점에 서 무생물 상태 이후에 나타난 개체밀도의 증가는 2017년과 다르게 완만하였으며, C. capitata 개체밀도의 증가는 Farm1 에서만 나타났다. 대조 정점의 개체밀도는 3개 정점 간 차이 가 적었으며, 공통적으로 봄철에 높고 여름철에 낮은 특성 을 보였다.

기존 정점의 정점별 다양도는 양식장에서 멀리 떨어진 Farm3에서 가장 높았으며, Farm1 ≤ Farm2 < Farm3으로 차이 를 보였다(Fig. 3c). 오염종이 급격히 증가한 시기와 무생물 발생 시기에 다양도가 낮았으며, 4월 전·후에 높았다. 대조 정점의 다양도는 가장 외해에 있는 Con3에서 다소 높았고, 봄철에 높고 여름철에 낮은 특성을 보였다.

3.3 저서다모류 우점종의 변화

기존과 대조 정점에서 출현한 다모류의 우점종을 조사하 고 이들 종의 시·공간 분포 특성을 살펴본 결과 크게 4개 유 형(Type1~IV)로 구분되었다. 먼저 Type I은 양식장 철거 직후 인 조사 초기에 밀도가 급격히 증가하고 조사후기에 밀도가 감소한 유형으로 C. capitata, Minuspio pulchra 가 해당하였다 (Fig. 4a, 4b). C. capitata는 어장환경평가 시 4그룹에 속하는 종으로 극심한 오염역이나 가두리 양식장 바로 아래에서 고 밀도로 출현하는 기회종이다(Pearson and Rogenberg 1978;Seo et al., 2015). 조사 초기에 양식장과 가까운 정점인 Farm1, Farm2에서 밀도가 급증하였는데, 2017년 10월 무생물 군집 이후인 2017년 11월부터 출현하기 시작하여 2020년 4월까지 680 ~ 7,980 개체/㎡, 평균 2,859 개체/㎡ 의 높은 밀도를 나타 내었다. 무생물 시기 이후 2020년 12월에 재가입하여 우점 하였으나 조사 초기보다 개체밀도가 낮았다. 양식장 철거 약 2년이 경과한 2020년에는 대조 정점과 유사한 개체밀도 로 출현하였다. 어장환경평가 시 4그룹 종인 M. pulchra는 기존 정점들 중에서 Farm3에서 우점 하였으며, 시간이 경과 하면서 무생물 상태 이후의 개체밀도 증가는 점차 감소하 였다.

Type II는 양식장 철거 이후 개체밀도가 증가하는 유형으 로 Lumbrineris longifolia 와 Pseudopolydora kempi가 해당하였 다. L. logifolia 는 우리나라 연안의 유기물이 풍부한 해역에 서 우점적으로 출현하는 지표종이다(Lim and Hong 1997;Jung et al., 2007). 기존 정점에서는 평균 35 개체/㎡, 대조 정 점에서는 평균 297 개체/㎡가 출현하여 대조 정점에서 높은 개체밀도를 나타내었으나, 대조 정점에서는 시·공간 분포의 뚜렷한 경향 없이 연중 출현하였고, 기존 정점에서는 양식 장 철거 직후 보다 양식장 휴지가 진행 중인 조사 후기, 봄 철에 개체밀도가 증가하였다. C. capitata의 분포와 유사하게 양식장에서 떨어진 거리에 따라 개체밀도의 차이가 있었으 나, 양식장에서 먼 Farm3에서 개체밀도가 높았다(Fig. 4c). P. kempi은 어장환경평가 시 그룹 1에 속하는 종으로, 주로 건 강한 해역에 출현한다. 대조와 기존 정점에서 시간에 따른 경향은 큰 차이가 없었으며, 밀도의 증가 정도는 기존 정점 이 더 뚜렷하였다(Fig. 4d).

Type III에 해당하는 종은 대조 정점에서만 우점하여 출현 하는 유형으로 Mediomastus sp.와 Aquilaspio krusadensis가 해 당하였다(Fig. 4e, 4f). Mediomastus sp.는 어장환경평가 시 1그 룹에 속하는 종으로 기존 정점들에서는 평균 3 개체/㎡가 출 현하였으나, 대조 정점에서는 평균 65 개체/㎡, 최대 880 개 체/㎡ 가 출현하였다. A. krusadensis는 그룹 3에 속하는 종으 로 대조 정점에서 평균 36 개체/㎡, 최대 140 개체/㎡의 밀도 로 연중 출현하였다. 이들 종 이외에 Notomastus latericeus, Praxillella affinis 등과 같이 1, 2 그룹에 속하는 종들이 Type III에 해당하였다.

Type IV는 특정 시기나 정점에 선호도를 보이지 않는 유형 으로, 환경에 대한 내성이 강한 종이 포함되었으며 3그룹에 속하는 Sigambra tentaculata가 이에 해당하였다. S. tentaculata 는 Farm1과 Farm2에서는 봄철에만 출현하였으나 Farm3과 Con 정점에서는 연중 출현하였다(Fig. 4g)

기존정점에서 Type I과 Type II을 대표하는 두 종인 C. capitata와 L. longifolia의 시·공간적 분포 변화를 살펴보면, 양 식장 철거 직후인 조사초기에 우점한 C. capitata의 개체밀도 는 휴지 기간이 지속될수록 감소하고, 조사초기에 낮은 개 체밀도로 출현한 L. longifolia가 증가하는 경향을 보여, 유기 물이 극심히 오염된 상태에서 일반적인 유기물 오염상태로 점진적으로 변화해 가는 과정임을 나타내었다. 또한 이러한 현상이 이전된 양식장에서 멀수록 명확하였다.

3.4 어장환경평가 적용

어장환경평가에서 사용하는 다모류 4개 그룹을 양호한 그 룹에서 서식하는 1, 2 그룹과 유기물 오염이 진행된 환경에 서 서식하는 3, 4 그룹의 종으로 구분하고, 정점별로 시간에 따라 상대적인 밀도가 어떻게 변화하는지 살펴보았다(Fig. 5). 먼저 기존 정점에서 무생물 시기는 Farm1이 총 8회, Farm2 가 총 5회, Farm3이 총 3회로 양식장과 가까운 정점일수록 많았다(Fig. 5a, 5c, 5e). Farm1, Farm2는 양식장 철거 이전인 2017년 10월 조사에서 무생물 군집 이었으며, 다모류가 출현 한 이후 11월부터 12월까지는 거의 3, 4그룹 종이 우세하다 가 2018년 5월 이후에 1, 2 그룹종이 50 %를 넘었다. 2018년 하계의 무생물 시기 이후에는 2019년 2월 조사에서 1, 2그룹 이 50 %를 넘었으며, 2019년에는 하계 무생물 시기 이후에는 2019년 12월에 50 %를 넘었다(Fig. 5a, 5c). Farm3은 조사 초기 에 기존의 두 정점보다 3, 4그룹의 우세가 상대적으로 낮았 으며 무생물 발생 이후 유지되는 기간도 짧고, 1, 2그룹의 상 대비도 상대적으로 우세하였다(Fig. 5e).

기존 정점들은 하계의 무생물시기, 3, 4 그룹에 우세한 시 기, 1, 2그룹이 우세한 시기로 순환되는 양상을 보였다. 다만 시간이 경과함에 따라 1, 2그룹이 우세하는 시기와 1, 2그룹 의 상대적인 비율은 점차 증가하였고, 양식장에 가까운 정 점에서 변화폭이 상대적으로 컸다. 대조 정점은 여름철에 3, 4그룹의 비율이 상승하였으며, 그 외 시기에는 1, 2그룹의 비율이 우세하였다(Fig. 5b, 5d, 5f).

저서동물지수와 총유기탄소의 시·공간 변화를 살펴보고 어장환경평가 등급을 통해 저층 환경의 건강도 변화를 평가 하였다. 먼저 대조 정점은 저서동물지수와 총유기탄소량이 일부시기를 제외하고 전 조사 기간 동안 2등급을 유지하였 으며, 어장환경평가 등급도 2018년 12월 Con2의 1등급을 제 외하면 전 조사 기간 동안 2등급으로 큰 변화가 없었다(Fig. 6, Table 2).

기존 정점은 시기별 정점별 변화가 뚜렷하였는데, 특히 저서동물지수는 여름철 발생한 무생물 시기가 지수에 반영 되어 지수의 변화가 컸다(Fig. 6a). Farm1, Farm2 정점은 조사 초기부터 2018년 4월까지 Farm3은 2월까지 저서동물지수가 4등급이었고, 4등급 내에서도 지수가 조금씩 상승하였다. 이 후 기존 정점 모두 2등급까지 상승하였다가 무생물 시기가 되며 다시 4등급으로 하락하였다. 무생물 시기 이후 다시 저 서동물지수가 상승하였는데, 이때는 조사 초기와 달리 급격 하게 상승하여 기존 3개 정점 모두 2등급까지 상승한 시기 가 2019년 2월이었다. 저서동물지수와 총유기탄소량을 합산 한 어장환경평가 등급은 Farm1과 Farm2는 2018년 4월까지는 3 ~ 4등급을, Farm3은 3등급이었고, 2018년 5월부터는 기존 3개 정점 모두 2 ~ 3등급을 나타내었다.

3.5 저서다모류 군집구조의 시·공간 변화

연구 지역의 6개 정점, 22회 조사에서 수집한 132개의 케 이스를 대상으로 주좌표분석(Principal coordinates analysis)을 수행하여 유사 특성에 따라 소수 그룹으로 구분하고, 각 그 룹이 생물적, 화학적 요인과 어떠한 상관성을 갖는지 알아 보았다(Fig. 7). 먼저, 이차원 평면상에서 대조 정점의 케이스 들은 거의 모여 있었으나 기존 정점의 케이스들은 평면상에 산포하여 분포하였다. 이들을 유사 특성에 따라 분류해보면, Control Group, GroupⅠ, Ⅱ, Ⅲ으로 구분되었다.

Control Group은 2차원 평면상에서 유기물 오염과 관련된 환경 변수와는 역방향, 어장환경평가 시 1 ~ 3그룹에 속하는 종들과는 정방향의 관계를 보였다. 그러나 Control Group에서 L. longifolia과의 상관성이 상대적으로 컸다는 점에서 유기물 함량이 낮은 편은 아니었다. GroupⅠ은 유기물 관련 환경 변 수에 정방향, 용존산소와 염분도에 역방향, 수온에 정방향의 관계를 보였으며, 2차원 평면상에서 가장 넓게 산포하여 교 란 상태에 있음을 암시하였다. 높은 유기물함량과 낮은 용 존산소량으로 인해 생물이 거의 출현하지 않은 하계와 추계 의 기존 정점들이 포함되었다. GroupⅡ는 유기물 변수와 용 존산소는 정방향, 수온은 역방향의 관계를 보여 심각하게 교란되었던 군집이 1차적으로 천이하는 단계였다. 기회주의 적 특성이 강한 C. capitata가 강한 상관성을 보였으며, 늦가 을부터 초봄의 조사시기가 해당하였다. GroupⅢ은 유기물 변수들과 수온은 역방향, 용존산소는 정방향에 위치하며 봄 철이 해당하였다. 어장환경평가 시 그룹 1 ~ 3에 해당하는 종 들과 높은 상관성을 보였으며, 무생물 시기 이후 나타난 기 회종 우점의 1단계 군집 천이 이후의 다음 단계로의 천이 과 정을 포함하였다. 기존 정점은 시간의 경과에 따라 GroupⅠ, GroupⅡ, Group Ⅲ의 특성으로 변화하며 평면상에서 대조정 점과 가까운 거리로 이동하였으나, 다시 GroupⅠ로 돌아간 이후 매년 순환이 반복되었다. 다만, GroupⅠ은 조사가 진행 될수록 중심과 멀리 떨어져 있던 케이스들이 2차원 평면상 의 중심으로 이동하는 경향을 보였다.

3.6 저서동물 분석을 통한 양식장 환경의 회복성 진단

양식장 휴지는 퇴적물의 완전한 회복을 기대하는 것이 아 니라, 지속가능한 양식장 이용을 보장하기 위한 것이다 (Pereira et al., 2004;Black et al., 2008). 퇴적물의 회복은 크게 화학적 회복과 생물학적 회복의 두 가지 유형을 포함하고 있으며(Brooks et al., 2003, 2004), 여러 연구에서 생물학적 회 복이 화학적 회복에 비해 더디게 진행된다고 보고하였다 (Pereira et al., 2004). 또한 어류양식장 유기물 방출에 영향을 받는 연질의 퇴적물에서는 생물학적 회복이 일어나려면 화학 적 회복이 먼저 필요하다고 보고하였다(Aguado-Gimenez et al., 2012;Brooks et al., 2004, 2003;Pearson and Rosenberg, 1978).

양식시설물 제거 이후 퇴적물 회복의 측면으로 모니터링 한 결과, 총유기탄소는 양식장 제거 약 1년 이후에 대조 정 점과 유사한 농도로 감소하여 화학적 회복은 어느 정도 나 타난 것으로 판단된다. 그러나 저서동물 군집의 생태지수, 우점종, 건강도 평가 등을 바탕으로 한 생물학적 요인들은 회복과 교란이 반복되는 결과를 보여주고 있어, 생물학적 회복을 논의하기에는 부족한 상태이다.

본 연구가 진행된 하동 어류양식장의 저층 지형은 동서로 약 250 m, 남북으로 약 800 m의 골짜기가 형성되어 있으며, 기존의 양식장과 이전된 양식장 모두 이 골짜기 내에 위치 하였다. 골짜기 내, 외측의 수심차이는 약 4 ~ 5 m로, 하계에 수온약층이 형성되면 저층은 빈산소층이 발생하고 퇴적층 은 무생물 상태가 된다. 이러한 계절적 현상이 다변량 분석 에서 보여주듯이 저서동물 군집의 회복과 악화의 반복에 영 향을 준 것으로 판단된다(Fig. 7b).

하동 양식장 조사결과에서 긍정적인 결과는 조사 후반기 에는 무생물 발생 기간이 짧아지고, 회복되는 속도가 빨라 졌다는 점이다. 이러한 변화는 양식장에서부터 멀리 떨어져 있을수록 명확하였다. 따라서 연구지역은 회복에는 유리하 지 않은 지역이나, 추가적으로 유기물이 공급되지 않는다면 회복이 지속적으로 진행될 수 있을 것으로 판단된다. 따라 서 지속적인 모니터링으로 기존 정점의 회복상태를 확인할 필요가 있을 것으로 보인다.