서 론

만경강은 전라북도에 위치한 농업용수 등의 수자원 공 급을 위한 중요한 수체로서 유역면적 1,527.1 km², 유로연 장 77.4 km을 가지고 있고, 고산천, 소양천, 전주천 등의 주 요한 지류를 가지고 있으며, 최종적으로는 새만금으로 유 입되는 하천이다. 만경강 유역은 전북의 주요 농업생산지 대로서, 농경지 및 산림이 유역면적의 대부분을 차지하고 있을 뿐 아니라, 만경강의 중류부와 전주천의 중·상류부 에 도시와 공업단지가 밀집되어 있고, 익산천의 상류역에 대규모 축산단지가 형성되어 있어 비점오염원 및 점 오염 원이 수질특성에 영향을 주는 것으로 보고되고 있다 (Yun et al. 2002).

이 유역내의 오염원인 생활하수, 축산폐수 및 농업배수 등의 유입은 만경강의 수질저하의 주요 원인으로 평가되 어 왔으며, 특히 전주천과 익산천의 합류 이후에 수질이 급격하게 악화되고 있는 것으로 보고되어 왔다 (Yun et al. 2002;Lee et al. 2005). 최근 새만금 유역의 수질개선을 위 해 만경강의 수질에 대한 관심이 높아지고 있으며, 이에 따라 정부에서는 수질을 향상시키기 위해 다양한 노력들 이 진행되어 왔다. 특히, 최근 수질개선 뿐만 아니라 인간 과 자연생태계가 공존할 수 있는 친수공간의 조성을 목적 으로 만경강의 8개 지구에 하천환경정비사업이 진행되고 있어 공사 이전 및 공사단계에서의 지속적인 모니터링이 필요한 실정이다.

본 연구대상 수계에서의 선행연구는 주로, 녹조문제, 오 염부하평가, 수질평가, 생물상 평가 및 군집분석 등에 대한 연구가 수행된 바 있다. 특히, 만경강의 본류 및 지류에서 의 이·화학적 수질변화 (Yun et al. 2002;Oh et al. 2004;Lee et al. 2004;Lee and Park 2008) 및 오염부하량 평가 (Cho 2003;Lee et al. 2005;Kim et al. 2013;Ryu and Lee 2013;Kim et al. 2017) 등이 집중적으로 연구되며, 인 (P), 질소 (N)와 같은 수질오염 및 부영양화의 문제를 주로 다루어 오고 있다. 또한, 본 유역에 서식하고 있는 담수어류 및 저 서무척추동물에 대해서도 생물상 및 군집분석 (Kim and Lee 1998;Hwang 2006;Lee et al. 2008;Park et al. 2009;Jeon 2010;Kim et al. 2015)에 대한 연구가 여러 차례 보고되어 왔다. 그러나 이런 다양한 연구에도 불구하고 대부분의 연 구가 수질 분석 및 군집분석에만 초점을 맞추고 있어, 수 계의 트로픽길드, 생물 내성도 등에 대한 수질 영향 및 수 생태 건강도에 대한 화학적 수질영향에 대한 연구는 극히 미미한 실정이다. Jeon (2010)에 의해 저서동물을 이용한 KSI (Korean saprobic index) 값을 산정하여 수생태계 건강 성이 평가된 바 있지만, 조사기간이 단기간에 그쳐 장기간 생태학적 건강도를 평가하기 위한 장기연구가 절실히 요 구되고 있다.

대부분의 생태계에서 단기간의 연구를 통해 나타나는 결과는 장기간의 연구를 통해 나타나는 결과에 비해 과장 되어 나타나는 경우가 많기 때문에, 단기적인 연구는 생태 계 내에서 일어나는 변화를 합리적으로 파악하기에 불충 분하다고 할 수 있다. 특히, 생물군집과 같이 복잡한 구조 를 가지는 생태계 내에서 일어나는 상호작용은 단기적인 연구를 통해서는 파악할 수 없으므로, 반드시 장기적인 연 구가 진행되어야 한다 (Spellerberg 1991;Lee 2006). 한편, 어류는 다른 지표 생물종에 비해 한 세대가 길어 장기간 의 환경변화를 반영하고 있을 뿐 아니라 이동성을 가지고 있어 다양한 서식조건을 대변한다. 또한, 어류는 하천생태 계 내의 최상위 소비자로서, 먹이사슬을 통한 오염의 최종 적인 영향을 반영하고 있을 뿐 아니라 실제 채집현장에서 비교적 동정이 용이해 생물학적 통합수질 평가를 위한 모 델생물로 널리 인정받아왔다 (Jones et al. 2005;An and Lee 2006). 따라서 이러한 어류를 이용한 생태건강도 평가모 델인 생물통합지수 (Index of Biological Integrity; IBI, Karr 1981)는 Yeom et al. (2000)에 의해 국내실정에 맞게 수 정 및 보완되어 금호강 (Yeom et al. 2000), 평창강 (An et al. 2001a), 갑천 (Bae and An 2006) 등 다양하게 적용되어 왔 으며, 최근 환경부의 국가 수생태계 건강성 평가 프로그램 의 어류평가지수 (Fish Assessment Index; FAI, ME 2016a) 로 보완되어 적벽강의 생태건강도 평가 (Lee et al. 2017) 등에 적용되었다.

따라서 본 연구는 장기 모니터링을 통해 만경강 본류의 이·화학적 수질 변화 및 담수 어류상의 변화를 파악하고 FAI모델을 적용하여 만경강 본류의 생태건강도를 진단하 는데 목적이 있다.

재료 및 방법

1. 조사시기 및 조사지점 현황

본 연구는 만경강 본류 5개 지점을 선정하여 2009∼ 2016년까지 8년 동안 매년 1회에 걸쳐 어류 현장조사 및 하천 생태건강도 평가를 실시하였다 (Fig. 1). 우리나라의 장마기에 해당하는 7∼8월은 유량증가 및 유속변화 등에 의해 어류군이 이동할 가능성이 있기 때문에, 어류 채집 및 하천 생태건강도 평가는 장마 전 수체가 안정되는 시기 인 4월∼6월에 실시하였다.

본 연구에서 조사지점은 상류에서 최하류역까지 5개 지 점을 선정하였다. 조사지점 중 최상류에 위치한 S1지점은 Strahler (1957)에 의거할 때 3차 하천으로, 대아저수지 수 문 하류에 위치해 있는 지점으로서 하류구간이 인공제방 으로 이루어져 있으나 수변식생이 매우 발달해있어 주변 오염원으로부터 완충작용이 잘 이루어지고 있어 맑고 깨 끗한 상태를 유지하고 있다. S2지점은 4차 하천으로 고산 천이 합류한 후에 위치한 수체로서, 하상이 주로 크고 작 은 돌과 자갈로 이루어져 있으며 자연 하천의 형태가 잘 보존되어 있는 지점이다. 하천 수변부와 하천 내중도에 자 연식생이 잘 발달되어 있으며, 상·하류 구간 모두 여울부 가 잘 발달되어 있다. S3지점은 5차 하천으로 전주시의 생 활하수 및 공단폐수로 인한 수질 악화된 전주천이 합류되 기 전에 위치한 수체로서, 하폭이 넓고 하상이 주로 자갈 과 모래로 이루어져 있으며, 인근에 전미평야가 위치하고 있어, 제내지의 농경지가 오염원으로 작용할 것으로 사료 되는 지점이다. S4지점은 5차 하천으로 전주천이 합류한 후에 위치한 지점이다. 하상이 주로 자갈과 모래로 이루어 져 있으며, 하도 내 하중도가 잘 발달되어 있다. 전주천의 유입으로 인한 수질악화가 예상되며, 특히 전주하수종말 처리장 방류수로 인한 수질악화가 우려된다. S5지점은 조 사지점 중 최하류에 해당하는 지점으로 5차 하천에 해당 하며, 백구제수문 앞에 위치한 지점이다. 이 지점은 왕궁가 축분뇨처리장의 방류수로 인해 수질이 악화된 익산천이 합류한 후에 위치해 있어 가축분뇨처리수 유입으로 인한 수질악화가 우려된다. 세부 조사지점은 다음과 같다.

2. 화학적 수질 자료 분석

만경강 본류 수계의 화학적 수질상태를 분석하기 위 하여 2009년에서 2016년까지 측정된 환경부의 수질자 료를 물환경정보시스템 (http://water.nier.go.kr)에서 획 득하여 분석하였다. 현재 환경부에서 측정하고 있는 여 러 일반 수질 항목 중 하천의 수질특성을 가장 잘 반영 할 수 있는 6개 항목을 선정하여 분석에 이용하였다. 선 정한 6개 항목은 생물화학적 산소요구량 (biochemical oxygen demand; BOD), 전기전도도 (conductivity), 총질 소 (total nitrogen; TN), 총인 (total phosphorus; TP), 용존 산소 (dissolved oxygen; DO), 암모니아성-질소 (ammonia nitrogen; NH₄-N)로서 각각의 수질 값을 지점별로 변화양 상을 분석하였다.

3. 어류 현장조사 방법

어류의 현장조사는 미국의 Wading method (Ohio EPA 1989)에 의거하여 한국에 맞게 개발한 An and Lee (2006) 의 기법 (우리나라 환경부 조사기법은 이 기법에 기초함) 에 의거하였다. 어류 채집의 정량화를 위하여 (Catch Per Unit of Effort; CPUE) 채집거리는 200 m, 조사 소요시간 은 50분으로 한정하였고, 여울 (Riffle), 소 (Pool), 흐르는 곳 (Run)을 포함하여 조사를 실시하였다. 채집도구는 조 사지점의 특징에 따라 투망 (망목: 5 × 5 mm)과 족대 (망 목: 4 × 4 mm)를 이용하였다. 각 지점에서 채집된 어류는 현장에서 Kim and Park (2002)에 의거해 동정하였고 일부 현장에서의 분류가 어려운 경우 10% 포르말린 용액에 고 정하여 실험실로 운반 후 동정을 수행하였다. 어류의 과 (family)의 분류체계는 Nelson (1994)에 따라 분류하였으 며, 비정상어종의 외형적 감별은 Sanders et al. (1999)의 방 법에 따라 기형 (Deformity; DE), 지느러미 손상 (Erosion; EF), 피부손상 (Lesions; LE) 및 종양 (Tumors; TU) 등으로 구분하여 분석하였다.

4. 어류 생태지표 특성 분석법

내성도 길드 분석은 US EPA (1993) 및 Barbour et al. (1999)에 기반을 두었으며, 이를 기반으로 한국의 어종에 맞게 길드를 분석한 An et al. (2001a, b)에 의거하여 민감종 (Sensitive species; SS), 중간종 (Intermediate species; IS), 내 성종 (Tolerant species; TS)으로 구분하였다. 민감종은 수질 의 오염에 의해 쉽게 사라지는 어종으로 환경의 질적 변화 에 민감하게 반응하는 어종이다. 그에 반해 내성종은 수질 오염에도 불구하고 종수 및 분포범위가 증가하는 종이며, 중간종은 민감종과 내성종 사이에 포함되지 않는 종으로 구분하였다.

섭식특성 길드 분석은 Ohio EPA (1989)의 기본 분류 체계를 이용하였고, 이를 기반으로 한국의 어종에 맞 게 섭식길드를 분석한 An et al. (2001a, b)에 의거하여 잡 식종 (Omnivores; O), 충식종 (Insectivores; I), 육식종 (Carnivores; C) 및 초식종 (Herbivores; H)으로 구분하였 다. 잡식종은 동·식물 모두를 지속적으로 섭취하며, 충식 종은 주로 수서무척추곤충을 섭취한다. 육식종은 어류 및 다른 척추동물을 섭취하며, 초식종은 주로 수체 내 플랑크 톤과 부착조류를 섭취하는 종이다. 식성의 분류는 다양한 먹이원이 동시에 존재할 때 1차적으로 이용하는 자원에 근거하여 분류하였다.

5. 어류의 다변수 모델을 이용한 생태건강도 평가

본 연구에서는 Karr and Dudlry (1981)의 어류를 이용 한 생물통합지수 (Index of Biological Integrity; IBI)를 기반 으로 개발된 환경부의 국가 수생태계 건강성 평가 프로그 램의 어류평가지수 (Fish Assessment Index; FAI)에 의거하 여 하천 생태건강도를 평가하였다. 본 연구에서 이용한 메 트릭의 속성 및 모델의 세부적 특성은 ME (2016a)에 자세 히 기술되어 있다. 모델 분석에 이용된 메트릭은 M1 국내 종의 총 종수 (Total number of native fish species), M2 여울 성 저서종수 (Number of riffle benthic species), M3 민감종 수 (Number of sensitive species), M4 내성종의 개체수 비율 (Proportion of individuals as tolerant species), M5 잡식종의 개체수 비율 (Proportion of individuals as omnivores), M6 충식종의 개체수 비율 (Proportion of individuals as native insectivores) M7 채집된 국내종의 총 개체수 (Total number of individuals) M8 비정상종의 개체수 비율 (Proportion of individuals with disease, tumors, fin damage and other anomalies)의 총 8개이다. 계급구간은 “12.5”, “6.25”, “0”으 로 구분하였으며, 이에 대한 최대 점수는 8×12.5=100점 으로 산정하였다. 생태건강도 등급은 매우 좋음 (A, 100∼ 80), 좋음 (B, 80∼60), 보통 (C, 60∼40), 나쁨 (D, 40∼ 20), 매우 나쁨 (E, 20∼0)의 5개 등급으로 대별하여 평가 하였다.

결과 및 고찰

1. 화학적 수질특성

만경강 본류 유역에서의 지난 8년간 수질 자료를 분 석한 결과, 생물화학적 산소요구량 (BOD), 전기전도도 (Conductivity), 총질소 (TN), 총인 (TP), 암모니아성 질소 (NH₄-N) 모두 상류에서 하류로 갈수록 뚜렷하게 악화되 는 경향을 보였다. 특히, 전기전도도와 총인, 그리고 암모 니아성 질소의 값은 하류역 (S4, S5)에서 급격하게 증가하 는 경향을 보였다 (Fig. 2). 각 지점별 평균 생물화학적 산 소요구량 (BOD)은 S1에서 0.9 mg L^-1, S2에서 0.8 mg L^-1, S3에서 1.9 mg L^-1, S4에서 3.8 mg L^-1, S5가 5.9 mg L^-1으로 나타나, 환경정책기본법 하천생활환경 기준상 S1과 S2지 점은 Ia등급 (매우 좋음), S3지점은 Ib등급 (좋음), S4지점 은 III등급 (보통), S5지점은 IV 등급 (약간나쁨)상태에 해 당하는 것으로 나타나 수질 구배현상이 지점별로 뚜렷하 게 나타났다. 즉, S1∼S3지점까지는 맑고 깨끗한 수질상태 를 유지하고 있는 반면, S4∼S5지점으로 갈수록 수질이 저 하되고 있는 것으로 나타났다. 이는 각각 S4와 S5지점이 전주천과 익산천 합류 후에 위치한 지점으로, 전주천에서 유입되는 전주시의 공단폐수 및 하수종말처리장 방류수 와 익산천에서 유입되는 가축분뇨처리장방류수의 영향으 로 인해 수질이 급격하게 악화된 것으로 사료된다 (Oh et al. 2004;Kim 2016;Kim et al. 2017).

2. 어류의 종 조성 및 분포특성

어류 종 조성 분석은 2009∼2016년의 8년간 분석되었 고, 이 기간 동안 채집된 어류는 총 14개과 50종 총 3,200 개체로서, 잉어과 (Cyprinidae)가 30종 60%로 가장 많 이 채집되었고, 미꾸리과 (Cobitidae)가 3종, 검정우럭 과 (Centrarchidae), 동사리과 (Odontobutidae), 동자개과 (Bagridae), 메기과 (Siluridae), 퉁가리과 (Amblycipitidae) 가 각각 2종, 망둑어과 (Gobiidae), 숭어과 (Mugilidae), 꺽 지과 (Centropomidae), 연어과 (Salmonidae), 뱀장어과 (Anguilidae), 송사리과 (Adrianichthyidae), 드렁허리과 (Flutidae)가 각각 1종이 나타났다 (Table 1). 일반적으로 하천 주변부에 점오염원이나 비점오염원과 같은 인위적 교란이 없으면, 하천차수가 증가할수록 어종 수와 개체수 가 증가하는 경향이 나타나는데 (An et al. 2001b), 본 조사 지점의 경우, 하천차수가 증가할수록 어종 수와 개체수가 감소하는 경향을 보여 점오염원 및 비점오염원에 의한 영 향을 많이 받는 것으로 사료된다. 특히 전주천 합류 후에 위치한 S4지점에서 가장 적은 종과 개체수가 채집되어 공 단폐수 및 하수처리장방류수와 같은 점 오염원에 의해 수 질악화 및 수생태계 교란이 심화된 것으로 사료된다.

각 조사 지점별 외래종 (Exotic species), 여울성저서종 (Riffle-Benthic species), 비정상종 (Abnomalities) 및 멸종 위기종 (Endangered species)의 상대빈도를 분석한 결과에 따르면, 상류에서 하류로 갈수록 비정상어종의 비율이 증 가하고 여울성저서종의 비율이 감소하는 경향이 뚜렷하 게 나타났다 (Fig. 3). S1과 S2지점은 외래종 및 비정상어종 의 상대빈도가 아주 낮게 나타고 여울성저서종의 상대빈 도가 높게 나타나 수생태계 교란이 아주 적은 안정된 수 생태계에 해당하는 것으로 사료된다. 특히 S2지점의 경우, 외래종의 서식이 확인되지 않고 (Table 1), 여울성저서종 의 상대빈도가 조사 지점 중 가장 높게 나타나 수생태계 건강성이 아주 우수한 것으로 사료된다. 뿐만 아니라, 멸 종위기종에 해당하는 퉁사리 (Liobagrus obesus) 및 감돌고 기 (Pseudopungtungia nigra)의 서식이 확인되어, 향후 이들 의 보존을 위한 대책이 필요할 것으로 사료된다. 반면, S4 지점은 외래종의 상대빈도가 조사 지점 중 가장 높게 나 타나 외래종에 인한 수생태계 교란이 아주 심한 것으로 사료되며, S5지점은 외래종의 상대빈도는 낮지만 조사지 점 중 비정상어종의 상대빈도가 가장 높게 나타나 수생태 계 교란이 아주 심각한 것으로 나타났다. 암모니아성 질소 농도는 pH가 과도하게 높을 때 (특히, 녹조의 과도한 번성 시기) 강한 어독성을 나타내기 때문에, 수층 내에 암모니 아성 질소의 농도가 높아지게 되면 어류는 지느러미와 아 가미가 손상될 수 있고, 병원균에 대한 저항성이 감소하게 되어 집단폐사의 가능성을 가진 것으로 알려져 있다 (ME 2016b). 또한 S5지점은 2011년에 어류의 집단폐사가 나타 났는데, 이는 가축분뇨 처리수의 유입에 의한 것으로 사료 되므로 향후 축산폐수의 유입 방지를 위한 방안 마련이 필 요할 것으로 사료된다.

3. 어류의 오염 내성도 및 트로픽 길드 분석

일반적으로 하천에서 유기물 오염 및 서식지 파괴 등에 의해 서식지의 질이 하강하게 되면 민감종의 종 수 및 개 체수는 감소하고 내성종의 종 수 및 개체수는 증가하게 된 다 (Karr 1981;US EPA 1991). 각 지점별 내성종 (Tolerant species)과 민감종 (Sensitive species)의 상대빈도 분석 결 과, 상류에서 하류로 갈수록 내성종이 증가하고, 민감종이 감소하는 경향이 뚜렷하게 나타났다 (Fig. 4) 특히, S3지점 부터 내성종의 상대빈도가 급증하는 동시에 민감종이 거 의 채집되지 않았을 뿐 아니라, S4와 S5지점에서는 민감종 이 단 한 종도 채집되지 않았다. 이러한 경향은 S4∼S5지 점으로 갈수록 수질이 저하되고 있는 것으로 분석된 이· 화학적 수질 결과와 일치하므로, 공단폐수 및 하수처리수 와 가축분뇨처리수로 인한 수질 악화로 인하여 S4와 S5지 점에서 내성종이 급증하고 민감종이 채집되지 않은 것으 로 사료된다.

또한, 각 지점별 잡식어종 (Omnivore species)과 충식어종 (Insectivore species)의 상대빈도 분석 결과, 상류에서 하류 로 갈수록 잡식종의 비율이 증가하는 경향이 뚜렷하게 나 타났으며, 충식종의 경우는 하류로 갈수록 감소하는 경향이 나타났다 (Fig. 4). 하천의 자연성 및 서식환경이 악화되면 충식종의 상대비율이 감소하기 때문에 충식종의 상대비율 을 이용해 하천의 서식환경 정도를 간접적으로 파악할 수 있다 (Moon et al. 2010). 조사 결과, S4지점에서의 충식종의 상대빈도가 가장 낮게 나타나 S4지점에서 하천의 자연성 및 서식환경에의 교란이 가장 심한 것으로 사료된다.

4. 화학적 수질에 따른 어류분포 특성 분석

이·화학적 수질의 변화는 어류의 분포에 영향을 주는 데, 일반적으로 수질이 악화될수록 내성종과 잡식종의 비 율이 증가하고, 민감종과 충식종의 비율이 감소하게 된 다. 이·화학적 수질변수와 내성·섭식길드 간의 Pearson 상 관도 분석에 따르면, 전기전도도 및 NH₄-N은 내성종 및 잡식종의 상대빈도와 통계적으로 유의한 양의 상관관계 (r>0.30, p<0.05)나타낸 반면, 민감종 및 충식종의 상대 빈도와는 통계적으로 유의한 음의 상관관계 (r<- 0.30, p<0.01)를 나타냈다 (Table 2). 이에 따라 유기물 및 암모 니아성 질소에 의한 수질악화에 의해 만경강 하류유역에 서 민감종 및 충식종의 비율이 감소한 것으로 사료된다. 특히 암모니아성 질소의 경우, 주로 가축분뇨 및 하수처리 폐수에 의해 유입되므로, 향후 전주천 및 익산천에서 유입 되는 하수처리수 및 가축분뇨처리수 관리 대책이 필요할 것으로 사료된다. 하천의 부영양화 정도를 파악할 수 있는 TN과 TP 역시 내성종의 상대빈도와는 통계적으로 유의 한 양의 상관관계를 나타낸 반면, 민감종의 상대빈도와는 음의 상관관계를 나타냈다. 특히, TP의 경우 TN과 달리 충식종의 상대빈도와 통계적으로 유의한 뚜렷한 음의 상 관관계를 나타냈다.

5. 어류 다변수 모델을 이용한 생태건강도 평가

어류 다변수 모델을 이용하여 만경강 본류의 생태건강 도를 평가한 결과, 어류평가지수 (FAI)는 평균 47 (n=40) 로 “보통상태 (C)”에 해당하는 것으로 나타났다 (Fig. 5). 한 편, 2015년의 FAI 점수는 35로 “나쁨상태 (D)”에 해당하는 것으로 나타나 본 연구기간 중 최저치를 나타냈다. 각 지 점별 FAI 값을 비교한 결과, S1과 S2지점이 각각 “좋음상 태 (B)”, “매우좋음상태 (A)”에 해당하는 것으로 나타나 생 태건강도가 우수한 것으로 나타났다. 반면, S3지점부터는 FAI값이 급격히 낮아져 S3, S4 및 S5지점 모두 “나쁨상태 (D)”에 해당하는 것으로 나타나 생태건강도가 심각하게 악화된 것으로 나타났다. 이는 S3∼S4지점 사이에 유입되 는 전주천과 S4∼S5지점 사이에서 유입되는 익산천에 의 한 수질악화 및 유기물 축적에 의한 하상 단순화 등의 서 식지파괴에 의해 생태건강도가 악화된 것으로 사료된다.

적 요

본 연구에서는 만경강 본류에서 2009∼2016년 동안 어 류의 트로픽 길드, 오염도, 내성도 및 어류평가지수 (Fish Assessment Index; FAI)에 기반을 둔 다변수 생태건강도 에 대한 화학적 수질 영향을 분석하였다. 수질 지표로서 전기전도도, 총인 (TP) 및 암모니아성 질소 (NH₄-N) 농도 는 최하류지점에서 가장 악화되는 것으로 나타났다. 어류 현장 조사 결과, 8년간 만경강 본류에서는 총 14과 50종 이 채집되었으며, 우점종은 내성종으로 잘 알려진 피라미 (Zacco platypus)로서 22.9%를 차지하여 트로픽 특성의 악 화현상을 보였고, 최류역에서는 비정상어종의 상대빈도 가 가장 높게 나타나 생태교란이 아주 심각한 것으로 나타 났다. Pearson 상관도 분석에 따르면, 내성어종 및 잡식어 종의 상대빈도는 전기전도도, BOD 및 NH₄-N와 통계적으 로 유의한 상관관계 (r>0.30, p<0.05)를 보였고, 민감어종 및 충식어종의 상대빈도는 이들 변수들과 유의한 역상관 관계 (r<- 0.3, p<0.01)를 보였다. 어류 다변수 모델 (FAI) 을 이용한 생태건강도로서 평가된 FAI 모델 값은 평균 47 (n=40)로서 건강도는 “보통상태 (C)”로서 나타났고, 하류 역 (S3∼S5) 지역에서 모델값은 모두 “나쁜상태 (D)”로 나 타나 수질악화에 의한 영향으로 평가되었다.