Journal 
서 론
2012년 완공된 4대강 보는 지속적으로 발생하는 지구 온난화 및 빈번한 가뭄과 홍수 피해 등의 자연재해로 인한 수자원 부족 사태를 해결하기 위해 조성되었으며 (Kang et al. 2012;Ji et al. 2015), 이·치수 효과를 극대화하여 수자원 문제를 해결하고자 하였다. 하지만 대형 인공구조물인 보 의 설치는 하천의 유속을 감소시키고, 퇴적물의 축적량을 증가시켜 시간이 지날수록 부영양화 및 수질오염을 가중 시켜 왔다 (Walther and Whiles 2008;Kang et al. 2012). 보에 의한 정체수역의 형성은 보 상류의 퇴적물 축적량을 증가 시키는 것으로 알려져 있으며 (Kang et al. 2012), 특히, 정체 수역 형성으로 인해 퇴적물에 농축되어 있던 유기물 (인, 질소, Chl-a 등)이 환경적인 변화 (수온변화, 유량 증·감)에 의해 재부유, 확산되어 다시 수중으로 용출될 수 있어 하천 환경에 직·간접적인 영향을 미칠 수 있는 것으로 알려져 있다 (Jung 2011a). 4대강 보 설치 이후 실시한 연구는 대부 분 수질 측정을 이용한 수질분석 연구 (Shin et al. 2015;Lee et al. 2017;Cho et al. 2018)가 수행되었으며, 수중생물 어류, 저서성 대형무척추동물 등을 이용한 연구는 부족한 실정 이다.
저서성 대형무척추동물 중 실지렁이류 (Aquatic Oligochaeta) 는 수생태계의 다양한 서식환경에서 가장 흔하게 발견되고, 풍부한 개체군을 유지하는 분류군으로 (Yoon et al. 2000;Jung 2011b), 국내에 서식하는 실지렁이류는 총 36종이 보고되어 있다 (Park et al. 2013). 실지렁이류는 서 식특성이 굴파는 무리 (burrower)로 유속이 느리고 Silt, Sand, Clay와 같이 작은 크기 입자로 이루어진 하상구조의 서식처를 선호하는 것으로 알려져 있다 (Lazim and Learner 1987;Kang et al. 2017). 또한, 오염 내성범위가 넓어 유기물 오염과 서식처 교란 (인공구조물 설치, 하천공사 등)으로 인해 수환경이 악화되어도 개체수가 증가하는 것으로 알 려져 있어 (Shin et al. 2008), 유기물오염평가를 위한 지표 생물로 이용되고 있다 (Brinkhurst and Kennedy, 1965;Lin and Yo 2008).
국내에서 실시한 실지렁이류에 대한 연구는 분류학적 연구 (Yoon et al. 2000;Jung 2011b;Park et al. 2013;Lee and Jung 2014), 생태독성 실험에 관한 연구 (Kang et al. 2016), 실지렁이를 이용한 수질정화 연구 (Choi 2005;Jun and Park 2005), 서식환경에 따른 생태 연구 (Syrovatka et al. 2009;Han et al. 2013;Kang et al. 2017) 등이 수행되었으나 국내 하천에 서식하고 있는 실지렁이류의 전반적인 생태 연구 자료는 부족한 실정이다. 따라서 본 연구는 4대강 보 설치 이후 변화된 서식처에 따른 실지렁이류의 분포 특성 비교 와 물리적·이화학적 요인에 따른 서식처 선호도를 파악하 여 실지렁이류를 이용한 4대강 서식처 관리 방안마련과 국 내 실지렁이류 생태 연구에 필요한 기초자료를 제공하고 자 수행하였다.
재료 및 방법
1. 조사시기
조사시기에 따른 변화양상을 파악하기 위해 2018년 3월 부터 10월까지 월 1회 조사를 실시하였으며, 수온상승으로 인해 개체수가 대량 발생할 가능성이 상대적으로 높을 것 으로 예상되는 6월과 7월은 월 2회 조사를 수행하였다.
2. 조사지점
조사지역은 한강 이포보 (IW: Ipo weir), 금강 세종보 (SW: Sejong weir), 영산강 죽산보 ( JW: Juksan weir), 낙동강1 강 정고령보 (GGW: Gangjeong-Goryeong weir), 낙동강2 달성 보 (DW: Dalseong weir)를 대상으로 조사를 수행하였다. 또 한, 4대강에 설치되어 있는 5개의 보에서 물리적·이화학적 인 서식환경 차이에 따른 실지렁이류의 분포양상을 파악 하기 위해 4개의 세부지점을 선정하였다. 세부지점 선정은 1) 보와 보의 중간지점을 보와보 사이 (B: Between weir to weir), 2) 기존 실지렁이류 개체수가 다수 출현했던 것으로 알려져 있는 지점을 주요발생지점 (M: Main point of outbreak), 3) 보를 기준으로 상류 500 m 이내에서 실시한 지점 을 보상방 (U: Upstream of weir), 4) 하류 500 m 이내에서 실시한 지점을 보하방 (D: Downstream of weir)으로 총 4개 지점을 선정하여 조사를 수행하였다. 각 조사지역 및 세부 지점의 GPS는 다음과 같다 (Table 1, Fig. 1).
3. 조사방법
1) 서식처의 물리적·이화학적 서식환경 분석
서식처 내 물리적 서식환경 분석은 레이저 거리측정계 (LASER 1200S, Nikon, Japan)를 이용하여 조사지역의 하폭 과 유폭을 측정하였으며, Digital water velocity meter (FP111, Global Water BA1100, USA)를 이용하여 유속과 수심을 측 정하였다. 하상구조는 국가 수생태계 건강성 평가 프로그 램 (ME 2013)에서 사용되는 평가 방법을 적용하여 Boulder, Cobble, Pebble, Gravel, Silt, Sand, Clay의 6단계로 구분하여 상대적인 구성 비율을 측정하였으며, 퇴적량 조사는 Core sampler (지름 80 mm×높이 1,000 mm)를 이용하여 실시하 였다.
서식처 내 이화학적 서식환경 분석은 2018년 3월부터 10 월까지 현장에서 수온, 용존산소, pH, 전기전도도를 각 조 사지점에서 휴대용 다항목측정기 (556MPS, YSI, USA)를 이용하여 측정하였으며, 탁도는 휴대용 탁도계 (HI93703, Hanna instruments, USA)를 이용하여 측정하였다. T-N (Total Nitrogen), T-P (Total Phosphorus), TOC (Total Organic Carbon), Chl-a (Chlorophyll a)는 현장수 2 L를 무균 채 수병에 표층수를 채수하여 실험실로 운반하여 수질오염공 정시험기준 (NIER 2016a)에 따라 분석을 실시하였다.
2) 채집 및 동정
실지렁이류 조사방법은 제4차 전국자연환경조사 지 침 (ME 2012), 국가 수생태계 건강성 평가 프로그램 (ME 2013), 생물측정망 조사 및 평가지침 (NIER 2016b)에 의거 하여 실시하였다. 정량 채집을 위해 Ponar grap (316 stainless steel, WIDCO, USA)을 이용하여 수심 0.5 m와 1.0 m에서 각각 2회 (0.1 m2) 채집 후 Sieve (mesh size 0.5 mm)를 이용 하여 체질 (Screening) 후 포집하였으며, 실지렁이류 개체 수는 1.0 m2 (0.1 m2×10)로 환산하여 분석을 실시하였다. 포집된 시료는 현장에서 90% ethyl alcohol 용액으로 고정 한 후 실험실로 운반하여 동정하였다. 실지렁이류의 동정 은 크기가 큰 개체의 경우 표본의 투명화를 위해 글리세린 에 10∼15분간 담가 투명화 실시 후 동정하였고, 작은 개체 들은 그대로 관찰하여 동정하였다. Timm (2009)과 Brinkhurst (1971)의 검색도설을 이용하여 동정을 실시하였으며, 두 검색도설을 통해 본 연구에서 출현한 모든 실지렁이류 의 종 동정이 가능하였다.
3) 건중량 측정
동정을 마친 표본은 60°C 조건에서 48시간 건조시킨 후 Chemical balance를 이용하여 소수점 넷째자리까지 건 중량을 측정하였으며, 측정이 불가능한 표본은 “0” 처리하 였다.
4) 통계분석
채집된 실지렁이류와 물리적·이화학적인 요인들 간의 상관성을 파악하고자 통계프로그램인 IBM SPSS ver. 21 을 이용하여 다변량 상관관계분석 (Multivariate correlation analysis)을 실시하였다.
결과 및 고찰
1. 물리적·이화학적 서식환경 분석
조사지역의 물리적 서식환경을 분석한 결과는 다음과 같다 (Table 2). 조사지점별 하폭은 250.0±0.0 (죽산보: 보 하방)∼1000.0±0.0 (강정고령보: 보하방) m의 범위로 강 정고령보 보하방이 가장 넓고, 죽산보 보하방이 가장 좁은 것으로 측정되었다. 수폭은 186.0±16.1 (죽산보: 보하방) ∼777.0±17.9 (강정고령보: 보상방) m의 범위로 강정고령 보 보상방이 가장 넓고, 죽산보 보하방이 가장 좁은 것으로 확인되었다. 수심은 47.0±9.0 (강정고령보: 주요발생지점) ∼76.0±22.4 (세종보: 보와보 사이) cm의 범위로 세종보 보 와보 사이가 가장 깊고, 강정고령보 주요발생지점이 가장 얕은 것으로 조사되었다. 실지렁이가 서식하는 퇴적토의 깊이는 10.8±4.6 (이포보: 보하방)∼63.8±8.3 (이포보: 주 요발생지점) cm의 범위로 이포보 주요발생지점이 가장 깊 고, 이포보 보하방이 가장 얕은 것으로 나타났다. 실지렁이 서식과 가장 연관성이 높은 하상구조인 Clay 비율은 55.0± 6.7 (세종보: 주요발생지점)∼10.0±0.0 (강정고령보: 보상 방)의 범위, Silt 비율 10.0±0.0 (이포보: 보상방, 강정고령 보: 보하방)∼32.0±4.0 (세종보: 보하방)의 범위, Sand 비율 11.0±2.0 (이포보: 주요발생지점)∼42.0±9.8 (이포보: 보 상방)의 범위로 대부분 보와보 사이, 주요발생지점, 보상방 에서 높은 비율을 나타냈으며, 특히 세종보에서 Clay의 비 율이 높은 것으로 확인되었다. 일반적으로 하천은 상류에 서 하류로 갈수록 Sand, Silt, Clay 등의 단순한 하상구조로 이루어지며, 다양한 유기물질이 축적되는 것으로 알려져 있다 (Allan 1995;Bea et al. 2003). 특히 횡단구조물 (보)과 같은 인공구조물은 물의 흐름을 방해해 유속감소로 인한 정수역화 및 유기 퇴적물의 축적을 가중시키게 된다 (Walther and Whiles 2008;Kang et al. 2012). 실지렁이류는 일반 적으로 Sand보다는 Clay와 같이 작은 크기의 입자로 이루 어진 하상구조를 선호하며 (Ladle 1971;Lazim and Learner 1987;Kang et al. 2017), 자갈 이상의 하상에서는 종수 및 개 체수가 감소하는 것으로 알려져 있다 (Brinkhurst 1971). 따 라서 퇴적토층의 깊이와 Clay의 비율이 높고, 자갈 이상의 하상재료 비율이 낮은 것으로 확인된 세종보는 실지렁이 류가 가장 선호하는 물리적 서식환경을 유지하고 있는 것 으로 판단된다.
조사지역의 이화학적 서식환경을 분석한 결과는 다음과 같다 (Table 3). 수온은 18.1±4.2 (이포보: 보와보 사이)∼ 19.7±3.7 (세종보: 주요발생지점)의 범위로 나타났다. 전 기전도도는 207.1±35.2 (이포보: 보와보 사이)∼327.7± 86.7 (달성보: 보와보 사이)의 범위로 측정되었다. DO는 8.0±2.2 (죽산보: 보와보 사이)∼9.4±2.8 (세종보: 보와보 사이)의 범위로 확인되었다. pH는 값의 차이는 적은 것으 로 나타났으며, 탁도는 7.0±5.1 (세종보: 보와보 사이)∼ 29.6±49.9 (세종보: 보하방) FTU의 범위로 세종보의 보하 방이 다른 지역과 조사지점들에 비해 매우 탁한 것으로 조 사되었다. 조사지점별 T-N, T-P 농도의 차이는 적은 것으로 확인되었다. TOC 농도는 2.5±0.6 (이포보: 보와보 사이)∼ 4.5±1.4 (달성보: 보상방) mg L-1의 범위로 나타났다. Chl-a 는 5.8±2.3 (이포보: 보하방)∼25.7±26.0 (죽산보: 보와보 사이) mg m-3의 범위로 지점에 따른 변화폭이 큰 것으로 확인되었다.
2. 실지렁이류 서식 분포 및 건중량, TOC 농도와의 상관분석
1) 조사지역별 실지렁이류 분포
실지렁이류 조사 결과 총 18종 35,390 ind. m-2 개체가 출 현하였으며, 각 조사대상지 별로 이포보 (한강) 11종 10,010 ind. m-2, 세종보 (금강) 10종 11,280 ind. m-2, 죽산보 (영산 강) 8종 5,500 ind. m-2, 강정고령보 (낙동강1) 7종 1,890 ind. m-2, 달성보 (낙동강2) 13종 6,710 ind. m-2 개체가 확인되었 다 (Table 4).
각 지점별 실지렁이류 종수 분석결과 달성보를 제외하 고 보와보 사이에서 가장 다양한 실지렁이류가 확인되었 으며, 달성보는 보상방에서 가장 많은 종이 출현하였다. 각 지점별 개체수 분석결과 이포보와 세종보의 경우 주요발 생지점에서 실지렁이류 개체수가 가장 많이 출현하였으며, 죽산보와 강정고령보는 보와보 사이, 달성보는 보상방에서 가장 많은 개체가 확인되었다 (Table 4). 일반적으로 실지렁 이류, 깔따구류 등 내성범위가 넓은 종들은 서식처 교란, 하 상구조의 단순화, 수질악화 등 서식처 수환경이 악화될 수 록 개체수가 증가하는 것으로 알려져 있는데 (Choi 2005;Na 2005;Shin et al. 2008;Lee et al. 2012), 조사지역 중 상대 적으로 실지렁이류의 개체수가 높게 출현한 이포보와 세 종보는 불안정한 수환경을 유지하고 있는 것으로 판단되 며, 대부분의 실지렁이 개체수는 보 설치에 따라 정수화된 지역 (보와보 사이, 주요발생지점, 보상방)에서 많은 개체 가 서식하는 것으로 확인되었다.
조사지역에서 출현한 실지렁이류의 우점종 및 아우점종 을 분석한 결과는 다음과 같다 (Table 4). 조사지역별로 이 포보는 실지렁이 (Limnodrilus hoffmeisteri)가 73.5%로 우점, 아가미지렁이 (Branchiura sowerbyi)가 19.5%로 아우점, 세 종보는 실지렁이가 78.0%로 우점, 긴뿔실지렁이 (Limnodrilus claparedeanus)가 9.1%로 아우점, 죽산보는 실지렁이가 73.3%로 우점, 참실지렁이 (Tubifex tubifex)가 15.6%로 아우 점, 강정고령보는 실지렁이가 82.5%로 우점, 참실지렁이가 7.4%로 아우점, 달성보는 실지렁이가 82.7%로 우점, 몽똑 물지렁이 (Dero obtusa)가 4.0%로 아우점하는 것으로 분석 되었다. 조사지점별로 실지렁이가 73.3∼82.7% 범위의 매 우 높은 비율로 우점하고 있었으며, 죽산보 주요발생지점 에서만 참실지렁이가 우점하고 있는 것으로 나타났다.
2) 실지렁이 건중량 및 TOC 농도와의 상관분석
조사지역에서 출현한 실지렁이류의 건중량 분석결과 지 역별로 이포보 총 16.194 g m-2, 세종보 총 10.227 g m-2, 죽 산보 총 2.369 g m-2, 강정고령보 총 1.497 g m-2, 달성보 총 2.131 g m-2로 측정되어 상대적으로 이포보와 세종보에서 출현한 실지렁이류의 건중량이 높은 것으로 확인되었으 며, 이는 실지렁이류의 개체수 출현율과 비례하는 것으로 나타났다. 각 지점별로는 이포보 1.110 (보와보 사이)∼ 7.127 (보하방) g m-2, 세종보 0.903 (보와보 사이)∼3.586 (보하방) g m-2, 죽산보 0.323 (보와보 사이)∼0.931 (보상 방) g m-2, 강정고령보 0.018 (보하방)∼1.388 (보와보 사이) g m-2, 달성보 0.100 (보하방)∼0.823 (보상방) g m-2의 범위 로 분석되었다 (Table 4). 한편, 달성보의 경우 죽산보 보다 실지렁이류의 개체수가 많이 채집되었으나 건중량에서는 죽산보 보다 낮은 것으로 계측되었는데, 이는 달성보에서 채집된 실지렁이류 대부분이 개체 크기 (total length)가 작 은 반면 죽산보에서는 상대적으로 크기가 크고 비대한 아 가미지렁이가 다수 출현했기 때문인 것으로 판단된다.
조사지역에서 출현한 실지렁이류의 과별 건중량은 실지 렁이과 (Tubificidae) 32.245 g m-2, 낚시지렁이과 (Lumbricidae) 0.156 g m-2, 물지렁이과 (Naididae) 0.016 g m-2, 실지렁 이붙이과 (Lumbriculidae) 0.001 g m-2의 순으로 나타났으며, 종별 건중량은 실지렁이 15.507 g m-2, 아가미지렁이 14.978 g m-2로 총 건중량 32.418 g m-2의 94.0%에 해당하는 건중 량을 차지하고 있는 것으로 나타났다 (Table 5). 실지렁이류 출현 종별 건중량과 TOC 농도와의 상관 분석결과 실지렁 이와 아가미지렁이에서 TOC 농도와 통계적인 음의 상관 성을 갖는 것으로 분석되었다 (Table 6). 일반적으로 실지렁 이류는 유기물과 양의 상관성을 갖는 것으로 알려져 있으 나 (Brinkhurst 1971), 본 연구결과에서는 실지렁이와 아가 미지렁이가 TOC 농도와 음의 상관성을 보여 기존 연구와 는 반비례 관계로 나타났다. 이는 실지렁이와 아가미지렁 이 개체수의 증가는 유기물의 양보다는 물리적인 서식환 경 (하상구조), 생활사 차이, 포식자의 증가 등 다양한 요인 들에 의한 영향이 더 큰 것으로 판단되며, 향후 추가적인 연구수행을 통해 개체수 변화의 명확한 요인 파악이 필요 할 것으로 생각된다.
3. 다변량 상관분석 (Multivariate correlation analysis)
물리적·이화학적 요인을 대상으로 실지렁이류의 종수 및 개체수와의 상관관계를 파악하기 위해 다변량 상관분 석을 실시하였으며, 실지렁이류와 물리적 요인들 간의 상 관분석 결과는 다음과 같다 (Table 7). 실지렁이류의 종수 와 퇴적토의 깊이, Clay, Silt의 비율과는 양의 상관성을 갖 는 것으로 나타났으며, Sand, Cobble, Pebble의 비율과는 음 의 상관성을 갖는 것으로 분석되었다. 실지렁이류의 개체 수는 퇴적토의 깊이, Clay의 비율과 양의 상관성을 갖고, Sand의 비율과는 음의 상관성을 갖는 것으로 분석되었다. 따라서 실지렁이류는 물리적인 요인 중 퇴적토의 깊이가 깊고, Clay, Silt의 비율이 높을수록 증가하며, Sand, Cobble, Pebble의 비율이 높을수록 감소하는 것으로 판단된다. 실 지렁이류의 개체수와 통계적으로 높은 상관성을 갖는 퇴 적토의 깊이와 Clay의 비율을 대상으로 비교·분석한 결과 (Fig. 2) 조사지역별 실지렁이 개체수와 퇴적토의 깊이 및 Clay의 비율은 개체수의 증감에 따라 매유 유사한 양상을 보였다. 실지렁이의 개체수가 상대적으로 많이 출현한 세 종보에서 퇴적토의 깊이와 Clay의 비율이 가장 높았으며, 개체수가 적은 강정고령보에서 Clay의 비율이 가장 낮은 것으로 분석되었다. 실지렁이류는 일반적으로 수환경이 악 화되고, 단순한 하상구조에서 개체수가 증가하는 것으로 알려져 있는데 (Shin et al. 2008), 본 조사결과 다른 조사지 역에 비해 상대적으로 많은 개체수가 출현하고, Clay 비율 및 퇴적토의 깊이가 깊은 세종보에서 수환경이 가장 불안 정하게 유지되고 있는 것으로 판단된다.
실지렁이류와 이화학적 요인들 간의 상관분석 결과는 다음과 같다 (Table 8). 실지렁이류의 종수 및 개체수와 pH, TOC는 모두 음의 상관성을 갖는 것으로 분석되었다. 다변 량 상관분석 결과 하상구조가 단순한 Clay 및 Silt의 비율이 높을수록 실지렁이류의 개체수가 증가한 반면 유기물의 함량이 많을수록 증가하는 TOC는 실지렁이류 개체수와 반비례 관계로 분석되었다. 이는 본 연구에서 채집된 실지 렁이류 대부분이 겨울과 봄철 유성생식을 통해 개체수가 증가하는 것으로 알려져 있는 실지렁이과에 속하는 종 (Brinkhurst 1971)들이 개체수이 대부분을 차지하였으며, 무성생식을 하고 유기물의 양이 많을수록 개체수가 빠르 게 (월 2회 이상 증식) 증가하는 물지렁이과에 속하는 종들 의 개체수 출현이 상대적으로 매우 적었기 때문인 것으로 판단된다.
4. 조사시기별 실지렁이류의 개체수, TOC 농도 비교 분석
조사지역에서 출현한 실지렁이류의 개체수와 TOC 농 도 간의 상관성을 파악하고자 조사시기별 분석을 실시하 였다 (Fig. 3). 3월부터 10월까지 실지렁이류 개체수는 감소 하고 TOC 농도는 증가하는 반비례 경향성을 나타냈다. 조 사시기, 실지렁이류 개체수, TOC 농도의 상관성분석 결과 조사시기는 실지렁이류 개체수와 음의 상관성을 나타내 고, TOC 농도와는 양의 상관성을 갖는 것으로 분석되었다 (Table 9). 일반적으로 유성생식을 하는 실지렁이과는 12월 부터 3월까지 산란하기 때문에 봄과 초여름 (1월~6월)에 개체수가 증가하였다가 다시 감소하는 것으로 알려져 있 는데 (Brinkhurst 1971), 전체 조사시기 중 3월부터 6월 1차 조사까지는 실지렁이류의 개체수가 증가하다 이후 감소하 는 경향성을 보인 본 연구 분석결과와 일치하는 것으로 나 타났다. 이는 물지렁이과에서 850개체 (1.6%)가 채집된 반 면, 실지렁이과에서는 53,320개체 (98.2%)가 출현했기 때 문에 나타난 결과로 판단된다. 일반적으로 실지렁이과 종 들은 대부분 Clay로 이루어진 하상에서 생활사를 유지하 는 반면, 물지렁이과 종들은 하천 수변부의 수생식물에 주 로 산란하고 서식하는 생활사 차이를 보이는데 (Brinkhurst 1971), 본 조사시 하천 수변부 조사가 아닌 Clay 비율이 높 고 퇴적토가 축척되어 있는 지역에서 조사가 이루어졌기 때문으로 판단된다.
적 요
본 연구는 4대강 보 설치에 따른 실지렁이류의 분포 특 성을 파악하고자 연구를 실시하였다. 현지조사는 2018년 3월부터 10월까지 총 10회 실시하였다. 조사결과 총 18종 35,390 ind. m-2 개체가 출현하였다. Limnodrilus hoffmeisteri (77.1%)이 우점종, Branchiura sowerbyi (9.5%)이 아우점종 으로 분석되었다. 건중량 분석결과 총 32.418 g m-2이 계측 되었으며, 각 종별 건중량과 TOC 농도와의 상관 분석 결과 실지렁이, 아가미지렁이와 TOC 농도 간 음의 상관성을 갖 는 것으로 나타났다. 실지렁이류와 물리적·이화학적 요인 들 간의 다변량 상관분석 결과 퇴적토 깊이, Clay 비율, Silt 비율, Cobble 비율, Pebble 비율, pH, TOC 항목과 상관성을 갖는 것으로 분석되었다 (p<0.05, p<0.01). 조사지역별 실 지렁이류 개체수와 퇴적토 깊이 및 Clay의 비율과 비교·분 석한 결과 퇴적토의 깊이와 Clay의 비율이 높을수록 개체 수가 많은 것으로 나타났다. 조사시기에 따른 실지렁이류 의 개체수, TOC 비교·분석결과 3월부터 10월까지 개체수 는 감소하는 경향성을 나타냈으며, TOC는 증가하는 경향 성을 갖는 것으로 분석되었다.
