Journal 
서 론
온배수 (Thermal discharges)란 발전소에서 사용한 냉각수, 온천, 폐수처리장 등에서 이용되는 자연 수온보다 높은 온 도의 상태로 주변의 하천, 호소 등에 방출되는 방류수를 의 미한다 (Kim 2000). 수온은 하천생태계에서 물리적 특성을 결정하고 생리기능, 생물의 분포 및 진화에 영향을 미치며 (Leeper and Taylor 1995), 자연 수온보다 높은 온배수는 하 천 생태계의 인위적 온도변화 요인으로 하천, 호소의 수온 을 상승시킨다 (Kim 2000). 온배수 방류로 인한 수온 상승은 겨울철에 하천이 결빙하지 않는 현상, 어류의 집단 폐사, 온 배수 방류지역에 안개가 빈번해지는 현상 등 하천 생태계의 일차생산 및 동·식물들의 종조성과 분포에 영향을 미친다 (Whitehouse 1971; Waite 1981; Marcogliese et al. 1989). 또 한 수온 상승은 하천의 용존산소량이 낮아지게 되어 하천생 물들의 호흡 비율을 증가시키고, 산소 소비를 빠르게 하며, 기생충 및 질병, 독성화학물질에 대한 민감성을 증가시키는 등 부정적인 영향을 미치게 한다 (Kil and Yoo 1999). 방출 되는 온배수가 주위 환경에 순기능의 역할을 하는 경우에는 열이득 (Thermal enrichment)이라고 하며 (Moore 1958), 수 온이 낮은 겨울에 동·식물의 활동을 활발하게 하여 메기, 연어, 굴양식 등에 이용되고 있다 (Shin and Kim 1995). 하지 만 역기능을 유발하는 열오염 (Thermal pollution)은 여름철 에 증식하는 우렁쉥이가 온배수 방류로 인하여 계절에 상관 없이 증식한다는 보고가 있었으며 (Shin and Kim 1995), 수 온이 높은 여름철에 북극 민물송어의 신진대사 증가로 인해 카드뮴과 납의 섭취가 증가하고, 빠른 성장으로 인하여 조기 사망에 이르게 되는 (Theurer et al. 1984; Kock et al. 1996) 등의 사례가 있다. 이처럼 하천 생태계의 수온 변화는 생물 들의 변화를 야기시키고 있으며 (Kwak et al. 2004), 도시하 수 및 비점오염원에서 유출되는 오탁수는 방류되는 지역뿐 만 아니라 인근 생태계에도 높은 영향을 미친다 (Field 1985; Vigon 1985; Hill 1986; Meyers 1986; Higgins et al. 1989).
하천 생태계의 생물군집은 동·식물성 플랑크톤, 연체동 물, 환형동물 등의 하등생물부터 담수어류 등과 같은 고등 동물까지 다양한 생물군으로 구성되어 있으며 (McCafferty 1981), 이러한 생물군 중에서 저서성 대형무척추동물은 하 천 생태계의 주요 변화요인인 물리적 요인, 화학적 요인, 생물학적 요인 등의 환경변화에 민감하게 작용한다 (Kim 2013). 대부분이 고착생활을 하고, 제한된 이동성을 가진 저 서성 대형무척추동물은 환경의 급격한 변화에 능동적으로 대처하는 능력이 낮고 (Im 2015), 비교적 뚜렷한 내성범위를 갖고 있어 환경 변화를 평가하는 지표생물로써 매우 효과적 인 것으로 알려져 있다 (Hynes 1960; Wiederholm 1983; Bea et al. 2003). 수온의 증가와 같은 환경의 변화는 저서성 대형 무척추동물의 개체군과 군집구조에 영향을 미치기 때문에 종조성 및 생물상의 변화는 하천 생태계의 변화를 유추하는 데 유용하게 사용된다 (Gray 1981). 하천의 생물군집은 홍수 와 가뭄 등의 자연적 교란의 영향 (Power et al. 1988; Resh et al. 1988; Kwak et al. 2004) 및 공해, 보 설치, 물의 이용 등에 지장을 받고 있으며, 수중생물에 영향을 주는 수질오 염 (Baek et al. 1996) 등의 인위적 교란에 종 특이성을 나타 내게 된다 (Kwak et al. 2004, 2008; Lee 2014). 또한 수온이 증가함에 따라 저서성 대형무척추동물의 종수가 증가하지 만 (Jacobsen et al. 1997) 민감종들의 분포는 제한하는 것으 로 알려져 있다 (Burgmer et al. 1997). 이러한 문제에도 불구 하고 온배수 방류에 의한 영향은 해양 생태계의 연구에 국 한되어 있으며 (Kim et al. 2000; Kim 2000, 2003; Go 2002; Lee 2016), 담수 생태계의 온배수 영향에 대한 연구는 전무 한 실정이다.
따라서 본 연구는 경기도 포천시에 위치한 온배수 배출시 설에서 방류되는 온배수가 저서성 대형무척추동물의 군집 구조 및 기능군에 미치는 영향을 파악하고, 온배수의 수온이 하천 건강성에 미치는 영향을 평가하여 하천생태계의 복원 및 보전을 위한 기초자료로 제공하고자 수행하였다.
재료 및 방 법
1.조사시기
조사시기는 온배수의 영향이 명확하게 나타날 수 있는 저 수온기를 고려하여 2015년 12월부터 2016년 2월까지 총 3 회에 걸쳐 조사를 실시하였고, 조사시기는 다음과 같다.
2.조사지점
조사지역은 온배수 배출시설에서 방류하는 온배수가 유입 되고 있는 경기도 포천시 영북면에 위치한 부소천을 대상으 로 선정하였다. 온배수 배출시설의 시설용량은 43 m3 일-1이 며, 주로 수영장 및 목욕용수로 이용되고 있었다. 조사지점 은 온배수가 유입되어 합류하는 혼합역 (MIX ZONE)을 기 준으로 상류 (UPPER ZONE)와 하류 (DOWN ZONE)에 각 각 1개 지점을 선정하여 정량 및 정성조사를 실시하였다 (Fig. 1).
-
UPPER ZONE (St. 1) E: 38°04′00.04″, N: 127°18′58.18″
-
MIX ZONE (St. 2) E: 38°04′03.02″, N: 127°18′49.12″
-
DOWN ZONE (St. 3) E: 38°04′06.37″, N: 127°18′40.36″
3.조사방법
1)수질 분석
수질 항목 중에서 수온, 용존산소 (DO), 탁도 및 전기전도 도는 멀티수질측정기인 Horiba U-22XD와 YSI를 이용하여 현장에서 직접 측정하였고, 생물학적 산소요구량 (BOD) 및 화학적 산소요구량 (CODMn), 총질소 (T-N), 총인 (T-P) 등의 항목은 현장수를 채수하여 시료병에 옮긴 후 냉장 보관하여 실험실 내에서 분석하였다.
2)채집 및 분류
저서성 대형무척추동물의 정량채집은 유량 및 물리적 환 경의 특성을 고려하여 계류형 정량채집망 (Surber sampler 30 cm×30 cm, 망목 0.2 mm)을 이용하였으며, riffle, run/ pool에서 지점별 총 3회를 실시하였다. 조사지점별 정확한 저서성 대형무척추동물상을 채집하기 위하여 Bottom kick net, Hand net (지름 18 mm, 망목 0.5 mm)을 사용하여 돌 밑, 수변부 등 가능한 다양한 미소서식처에서 정성채집을 병행 하였다. 채집된 저서성 대형무척추동물은 500 ml Vial에 담 은 후 현장에서 99% Ethanol로 고정하여 실험실로 운반하였 다. 실험실에서 유기물로부터 저서성 대형무척추동물을 골 라내어 Vial (10~25 ml)에 옮겨 80% Ethanol로 보존하였다.
저서성 대형무척추동물은 현미경을 이용하여 분류목록 및 검색표를 참고로 하여 동정하였다. 수서곤충은 McCafferty (1981), Yoon (1988, 1995), Merritt et al. (2008), Won et al. (2005) 등을 참고하였고 연체동물은 Kwon (1990), 거머리류 는 Song (1995)을 참고하였다. 깔따구류 (Chironomidae spp.) 는 정확한 종과 속 수준의 동정이 제한되어 과 수준 (Family level)으로 동정하였다.
3)EPT-group
하루살이목 (Ephemeroptera), 강도래목 (Plecoptera), 날도 래목 (Trichoptera)으로 구성된 EPT-group은 대부분 수환경 이 양호한 상태일수록 높은 출현을 나타내는 것으로 알려져 있으며 (Lenat 1988; Davis et al. 2003; Peitz 2003), 본 연구 에서는 EPT-group을 이용하여 수환경을 평가하였다.
4)군집 분석
조사시기별 정량적으로 채집된 저서성 대형무척추동물을 대상으로 출현한 분류군의 수를 비교하여 출현종수, 출현개 체수, 우점종, 우점도지수 (McNaughton 1967), 다양도지수 (Shannon-Weaver 1949), 균등도지수 (Pielou 1975), 풍부도지 수 (Margalef 1958)를 산출하였다.
5)기능군 분석
저서성 대형무척추동물을 대상으로 기능군 분포를 통하여 온배수의 영향을 분석하고자 출현종을 기준으로 섭식기능 군 (Funtional Feeding Groups, FFGs) 및 서식기능군 (Habitat Orientation Groups, HOGs)을 분석하였다 (Ro and Chun 2004; Merritt et al. 2008).
6)하천 건강성 평가
조사지역의 수생태계의 환경상태를 평가하기 위하여 한국 오수생물지수 (Korean Saprobic Index, KSI)를 이용하여 분 석하였다. KSI는 각 지표생물군의 오탁계급치 및 지표가중 치의 산정을 통하여 산출하였다.
7)통계 분석
통계 분석은 SYSTAT (Version No. 11.0 for windows)을 이용하여 상·하류 지점 및 온배수가 유입되는 혼합역에서 Correspondence analysis를 통해 수온과 저서성 대형무척추 동물의 관련성을 분석하였다. 수온과 개체수, 종수, 군집지 수, EPT-group, 냉수성종과의 통계적인 차이와 다양한 요인 들 간의 상관성을 파악하기 위해 SPSS (Version No. 18.0 for windows)를 이용하여 상관관계 분석 (Pearson Correlations Analysis)을 실시하였으며, 수온과 냉수성종의 관계를 분석 하기 위해 Multivariate analysis를 실시하였다. 또한 Paired samples t-test를 통해 조사지점별 수질항목의 유의성을 검증 하였다.
결과 및 고 찰
1.수질 분석
온배수가 유입되어 합류하는 혼합역 및 상·하류에서의 수질 분석 결과는 다음과 같다 (Table 1). 평균 수온은 상류 지점에서 2.83 (±1.89)으로 가장 낮았고, 혼합역에서 6.31 (±2.83)로 가장 높은 것으로 분석되었다. 상류와 하류지점 의 수온은 유의한 차이가 있는 것으로 나타났으며 (p<0.05), 이는 혼합역에서 유입되는 온배수의 영향에 따른 결과로 판 단된다. 평균 용존산소량 (DO)은 상류에서 14.61 (±1.83) 로 가장 풍부한 것으로 나타났고, 혼합역과 하류지점에서는 큰 차이가 나타나지 않았다 (p>0.05). 평균 생물학적 산소 요구량 (BOD)은 지점별 큰 차이는 없었지만 (p>0.05) 상류 에서 0.67 (±0.06)로 가장 낮은 것으로 분석되었고, 하류에 서 0.87 (±0.15)로 가장 높은 것으로 분석되었다. 평균 화학 적 산소요구량 (COD)은 상류와 혼합역에서는 큰 차이가 나 타나지 않았고 (p>0.05), 하류에서 1.71 (±0.34)로 가장 높 은 것으로 분석되었다. 평균 총 질소 (T-N)는 2.98 (±1.01), 평균 총 인 (T-P)은 0.13 (±0.14)으로 혼합역에서 가장 높은 것으로 나타났으며, 상류와 하류지점은 유의한 차이가 있는 것으로 분석되었다 (T-N, p<0.01; T-P, p<0.05). 질소는 생 물학적 고정, 공업적 질소고정 등에 의해 하천으로 유입될 수 있지만 인은 자연 상태에서 순환하지 못하고 폐수의 유 입 등 인위적인 교란에 의해 하천으로 유입되며 (Choi and Yu 1998), 일상생활 및 사업 활동 등으로 발생한 폐수가 공 공 수역으로 방류될 경우 유기물, 질소 등의 함량이 높은 것 으로 알려져 있다 (Jung et al. 1999). 본 조사지역에서 배출되 는 온배수는 대부분 목욕용수로 사용되고 있으며, 온배수 유 출로 인한 유기물의 유입이 T-N 및 T-P를 증가시키는 원인 인 것으로 판단된다.
2.종조성
조사기간 동안 출현한 저서성 대형무척추동물은 총 3문 4 강 11목 30과 50종 4,015개체가 출현하였다. 분류군별 출현 비율은 하루살이목 (Ephemeroptera)이 6과 16종 (32.0%)으로 가장 높게 나타났으며, 파리목 (Diptera) 5과 10종 (20.0%), 비 곤충류 (Non-Insects) 7과 9종 (18.0%), 날도래목 (Trichoptera) 6과 8종 (16.0%), 잠자리목 (Odonata) 1과 2종 (4.0%), 강도래 목 (Plecoptera), 딱정벌레목 (Coleoptera) 2과 2종 (4.0%), 뱀 잠자리목 (Megaloptera) 1과 1종 (2.0%)이 출현하였다. 분류 군별 개체수 출현비율은 파리목이 1,848개체 (46.03%)로 가 장 높게 나타났으며, 하루살이목 1,572개체 (39.15%), 비곤 충류 270개체 (6.72%), 강도래목 168개체 (4.18%), 날도래목 136개체 (3.39%), 잠자리목, 뱀잠자리목 8개체 (0.20%), 딱 정벌레목 5개체 (0.12%)가 출현하였다. 조사지점별 출현비 율은 하루살이목에서 가장 높은 것으로 나타났으며, 상류에 서 5과 8종, 혼합역에서 6과 8종, 하류에서 5과 13종으로 하 류지점에서 상대적으로 출현종 비율이 가장 높은 것으로 조 사되었다 (Fig. 2a). 날도래목은 유수생태계에서 매우 중요 한 위치를 차지하는 분류군으로 알려져 있다 (Hwang 2006). 본 조사시 상류와 하류지점에 비해 혼합역에서 상대적으로 매우 낮은 비율을 나타냈으며, 이는 유수생태계를 나타내는 상류 및 하류지점에서 보다 다양한 날도래목이 출현한 결 과로 판단된다. 조사지점별 개체수 출현비율은 하루살이목 의 경우, 상류와 하류지점에 비해 혼합역에서 상대적으로 매 우 낮은 비율을 나타내고 있었으며, 반대로 파리목의 경우 에는 혼합역에서 개체수 비율이 증가하는 것으로 분석되었 다 (Fig. 2b). 혼합역에서 개체수가 높게 나타난 파리목은 깔 따구과 (Chironomidae)에서 대부분 차지하고 있었으며, 깔따 구과는 청정수역부터 오염된 하천까지 광범위하게 서식하 는 수서곤충으로 내성범위가 강한 종으로 알려져 있다 (Na 2005). 특히, 깔따구류 (Chironomidae spp.)는 오염도가 높은 도시하천에서도 출현하며 (Yoo et al. 2005; Beak 2008), 낮 은 용존산소량, 넓은 온도범위 등 열악한 환경에서도 생존 이 가능한 특징을 가진다 (Ingersoll and Nelson 1990). 본 조 사시 혼합역에서 다른 분류군에 비해 고온의 방류수에 내성 이 높은 깔따구류가 상대적으로 높게 나타났으며, 이로 인하 여 파리목의 개체수 비율이 증가한 것으로 판단된다. 또한 오염에 내성이 비교적 강하고, 오염된 하천에서 넓은 범위의 수온대를 가지는 비곤충류의 실지렁이 (Limnodrilus gotoi) (Choi 2005; Lee et al. 2012) 역시 상·하류지점에 비해 온 배수가 유입되는 혼합역에서 높은 출현율을 나타내고 있었 다. 조사지점별 종조성은 유수역 및 냉수성을 선호하는 하루 살이목의 부채하루살이 (Epeorus pellucidus)와 민하루살이 (Cincticostella levanidovae) (Lee et al. 2012; Kim 2013)의 경우, 고온의 방류수로 인해 수온이 상대적으로 높은 혼합 역 보다 상류지점에서 높은 출현율을 나타냈으며, 혼합역에 서 하류지점으로 갈수록 하루살이목의 비율이 증가하는 것 으로 조사되었다. 상류지점과 혼합역에서 출현하지 않은 냉 수성을 선호하는 종인 감초하루살이 (Baetis silvaticus), 깜장 하루살이 (Nigrobaetis bacillus), 몽땅하루살이 (Ecdyonurus bajkovae)는 하류지점에서만 출현하였다. 또한 혼합역에 비 해 상·하류지점에서 유수역의 대표적인 분류군인 날도래 목의 꼬마줄날도래 (Cheumatopsyche brevilineata), 광택날 도래 KUa (Glossosoma KUa), 물날도래 KUa (Rhyacophila KUa), 검은머리물날도래 (Rhyacophila nigrocephala) 등의 종이 출현하였다. 조사지점별 개체수 비율이 높게 나타난 파 리목의 경우, 혼합역은 깔따구류가 대부분을 차지하고 있었 으나 상류와 하류지점에서는 유수역에서 서식하는 냉수성 선호종인 먹파리류 (Simulium sp.) (Kim 2013, 2014)가 높은 개체수를 구성하고 있었다. 또한 넓은 범위의 수온대를 가지 는 비곤충류의 실지렁이는 혼합역에 비해 하류지점에서 개 체수가 감소하는 것으로 나타났다. 전체적으로 상류와 하류 지점에서는 콩알하루살이 (Acentrella sibirica), 감초하루살 이, 깜장하루살이, 몽땅하루살이, 두점하루살이 (Ecdyonurus kibunensis), 민강도래 KUa (Nemoura KUa), 먹파리류, 물날 도래 KUa, 검은머리물날도래, 네모집날도래 KUa 등 냉수성 을 선호하는 종의 출현이 높은 것으로 나타났으며, 반대로 온도에 내성이 강한 실지렁이 및 깔따구류와 비교적 내성 의 범위가 넓은 동양하루살이 (Ephemera orientalis), 등딱지 하루살이 KUa (Caenis KUa) 등은 혼합역에서 높게 출현하 였다. 결과적으로 혼합역에서 상대적으로 출현종 및 개체수 가 가장 낮은 것으로 조사되었으며, 오염 및 수온의 내성 범 위가 높은 종들이 대부분을 차지하고 있었다. 따라서 고온의 방류수가 하천에 유입될 경우, 냉수성 종에 대한 서식환경의 교란으로 인해 고수온의 영향권내에 서식하는 저서성 대형 무척추동물의 군집구조에 부정적인 영향을 미치고 있는 것 으로 생각되며, 이로 인해 냉수성 종들의 분포에 제한을 받 는 것으로 판단된다.
3.EPT (Ephemeroptera-Plecoptera-Trichoptera)-group
하천생태계의 곤충강에서 주요 분류군으로 알려진 EPTgroup (Lee et al. 2012; Kim et al. 2013)은 출현비율에 따 라 수질의 청정도를 나타내는 것으로 알려져 있다 (Kwak et al. 2008). 조사기간 동안 출현한 EPT-group의 개체수는 총 1,876개체로 46.72%의 출현율을 나타냈다. 조사지점별 평 균 EPT-group의 비율은 혼합역에서 19.05% (±9.58)로 가장 낮고, 하류지점에서 65.66% (±30.27)로 가장 높은 것으로 분석되었다 (Fig. 3). 조사시기별로는 1차 조사 시 상류에서 218개체 (84.17%), 2차 조사 시 하류에서 574개체 (82.71%), 3차 조사 시 하류에서 440개체 (30.70%)로 혼합역에 비해 상류 및 하류지점에서 높은 비율을 나타내는 것으로 분석되 었다. 일반적으로 하천생태계에서 EPT-group의 비율이 높을 수록 양호한 상태를 나타내는 것으로 알려져 있으며 (Davis et al. 2003; Peitz 2003), 상류와 하류지점에 비해 낮은 EPT 비율을 나타낸 혼합역은 고온의 방류수로 인한 수온 교란으 로 인하여 청정수역을 선호하는 EPT-group에 불리한 서식 환경을 제공하고 있는 것으로 판단된다.
4.군집 분석
조사시기에 따른 조사지점별 우점종 분석 결과 (Table 2), 상류에서는 1차와 2차 조사 시 냉수성 선호종인 민하루살 이, 민강도래 KUa가 우점하고, 3차 조사 시에는 오염에 대 한 내성이 강한 깔따구 sp2.가 우점하여 조사시기별 다소 차 이가 있었다. 하류에서는 1차와 2차 조사 시 등줄하루살이 (Uracanthella rufa), 네점하루살이 (Ecdyonurus levis)가 우 점하고, 3차 조사 시에는 먹파리류가 우점하는 것으로 분석 되었다. 반면에 혼합역에서는 조사시기별 오염에 내성이 강 하고 넓은 수온대를 가지는 실지렁이와 깔따구 sp2.가 우점 하여 상·하류지점과 차이를 나타내었다.
조사지역의 군집 분석 결과 우점도지수는 평균 0.63 (± 0.11), 다양도지수는 평균 1.82 (±0.24)로 상류지점에서 상 대적으로 불안정한 군집을 나타내는 것으로 분석되었다 (Table 3). 이는 조사시기 동안 가장 낮은 수온을 나타낸 상 류에서 냉수성을 선호하는 부채하루살이, 민하루살이 등 하 루살이목의 종들이 높게 나타난 결과이며, 혼합역과 하류에 비해 상대적으로 불안정한 군집구조를 나타낸 것으로 판단 된다. 반면에 하류지점에서는 우점도지수 0.54 (±0.08), 다양 도지수 2.08 (±0.15)로 상대적으로 가장 안정적인 군집구조 를 유지하고 있는 것으로 분석되었다. 우점도지수는 군집의 현존량에서 특정종이 차지하는 비율을 나타내어 특정종으로 인한 생물군집의 구조와 기능을 판단하는 척도이며 (Lee et al. 2009), 저서성 대형무척추동물이 상류에서 하류로 갈수 록 나타나는 종 다양성의 차이는 식생 및 수질환경 등에 따 라 다양하다 (Harper et al. 1997; Jung et al. 2008). 또한 도심 하천을 제외한 대부분의 하천에서는 상류보다 하류에서 종 다양성이 높은 것으로 알려져 있으며 (Kil et al. 2010; Savić et al. 2013), 본 연구에서도 상류지점에 비해 하류지점에서 종다양성이 높게 분석되었다. 균등도지수는 평균 0.64 (± 0.08, 상류)~0.69 (±0.14, 혼합역)로 조사지점별 차이가 적 은 것으로 나타났으며, 풍부도지수는 평균 2.94 (±0.33, 상 류)~3.58 (±0.22, 하류)로 우점도지수 및 다양도지수와 마 찬가지로 상류지점 및 혼합역에 비해 하류지점에서 출현종 이 풍부한 것으로 분석되었다.
5.기능군 분석
수서곤충의 섭식기능군 (FFGs)과 서식기능군 (HOGs)의 구성은 서로 다른 하천유역에서 나타나는 서식처의 환경적 요인과 먹이자원의 특성을 나타낸다 (Cummins 1973; Ward 1992; Williams and Feltmate 1992; Kil et al. 2010). 섭식기 능군은 상류와 혼합역에서 Gathering-Collectors의 비율이 높고 하류에서는 Filtering-Collectors의 비율이 상대적으로 높은 것으로 나타났으며, 상류에서 하류로 갈수록 Shredders 의 비율은 감소하는 것으로 분석되었다 (Fig. 4a). 대형입상 유기물질 (Coarse Particulate Organic Matter: >1.0 mm)이 풍부한 상류에서 Shredders의 출현율이 높고, 하류로 갈수 록 미세입상유기물질 (Fine Particulate Organic Matter: >0.5 μm~1.0 mm)이 증가하는 것으로 알려져 있다 (Allan 1995; Bae et al. 2003). 본 연구에서도 낙엽과 같은 대형입상유기 물질이 높은 상류지점에서 Shredders의 비율이 혼합역과 하 류지점에 비해 높은 것으로 분석되었다. 미세입상유기물질 을 증가시키는 요인으로 공장폐수, 축산폐수 등의 점오염 원과 특정 배출지점이 없는 비점오염원이 있으며 (Charles and William 1992; Lee et al. 1994; Jung et al. 1997), 본 연 구에서는 온배수 방류로 인한 유기물질의 유입으로 미세입 상유기물질이 증가하여 유속이 낮은 혼합역에서 Gathering- Collectors가 높고 유속이 높은 하류로 갈수록 Filtering- Collectors가 증가한 것으로 판단된다. Predators는 상류와 하류지점에 비해 혼합역에서 가장 높은 것으로 분석되었고, Scrapers는 조사지점별 유사한 점유율을 나타내었다. 우리나 라의 하천은 발원지역을 제외하고 대부분의 하천이 연중 햇 빛에 노출되어 있는 상태이며, 1차생산력이 높아 Scrapers의 비율은 비교적 일정하게 나타나는 것으로 알려져 있다 (Shin et al. 2006).
서식기능군 분석 결과 상류지점은 Clingers와 Burrowers 가 대부분을 차지하고 있었으며, 혼합역은 Burrowers가 매 우 높고, 하류지점은 Clingers의 비율이 상대적으로 높아 지는 것으로 분석되었다 (Fig. 4b). 이는 하류지점에서 유속 에 적응력이 높은 네점하루살이, 부채하루살이, 민하루살 이, 검은머리물날도래 등이 출현한 결과로 판단된다. 유속 이 느리면 Clingers는 감소하고, Burrowers와 수초의 표면 을 기어오르는 소수의 Climbers와 하천의 바닥을 기어다 니는 Sprawlers가 증가하는 것으로 알려져 있다 (Yoon et al. 1992). 본 연구에서는 비교적 유속이 느린 혼합역에서 Clingers의 비율이 낮고, Burrowers와 Sprawlers의 비율이 상대적으로 높은 것으로 분석되었다.
6.하천 건강성 평가
조사지역의 환경상태를 평가하기 위해 한국오수생물지수 (Korean Saprobic Index, KSI)를 이용하였다. KSI 분석 결과 조사시기별 혼합역에서 가장 높은 지수값을 나타내는 것으 로 분석되었다 (Fig. 5). 조사지점별 평균 KSI 지수는 혼합역 에서 2.73 (±0.41)으로 C등급의 수환경을 나타내어 가장 불 량한 것으로 분석되었다. 상류지점은 1.12 (±1.25)로 B등급 의 양호한 환경상태를 나타냈으며, 하류지점은 0.71 (±0.38) 로 최적의 환경상태인 A등급을 나타내어 가장 양호한 수환 경으로 분석되었다. 1차 조사 시 상류에서 0.33으로 A등급, 2차 조사 시 상류에서 0.47로 A등급, 3차 조사 시에 하류에 서 1.15로 B등급의 가장 양호한 환경상태를 나타내는 것으 로 분석되었다. 전체적으로 혼합역에서 가장 불량한 환경상 태를 나타내었고, 상류와 하류지점에서는 조사시기별로 다 소 차이가 있는 것으로 분석되었지만 혼합역에 비해 상대적 으로 건강한 환경상태를 나타내고 있는 것으로 분석되었다. 수온의 증가로 인한 수환경의 변화는 저서성 대형무척추동 물의 개체군과 군집구조에 영향을 미치는 것으로 알려져 있 으며 (Gray 1981), 본 연구에서도 고온의 온배수 방류로 인한 수온의 변화는 수질에 민감한 저서성 대형무척추동물의 군 집구조 및 서식처 환경에 교란을 유발하는 요인으로 작용하 여 하천 건강성에 부정적인 영향을 미치는 것으로 판단된다.
7.통계 분석
온배수 유입으로 인한 수온의 변화와 조사지점별 서식하 는 저서성 대형무척추동물과의 대응 분석 결과 (Fig. 6) EPTgroup은 수온이 낮은 상류와 하류지점에서 관련성이 있는 것으로 분석되었다. 잠자리목은 수온이 높은 혼합역과 관련 성이 있는 것으로 나타났고, 비곤충류와 파리목은 수온과의 관련성이 적은 것으로 분석되었다. 수온이 가장 낮은 상류지 점에서 날도래목이 관련성이 있는 것으로 분석되었으며, 애 우묵날도래 KUa (Apatania KUa)가 수온과의 관련성이 가 장 높은 것으로 분석되었다. 하류에서는 강도래목과 하루살 이목이 관련성이 있는 것으로 분석되었고, 강도래목의 녹색 강도래 (Sweltsa nikkoensis)가 가장 관련성이 높은 것으로 나 타났다. 한편, 하루살이목의 등딱지하루살이 KUa는 혼합역 과 관련성이 있는 것으로 분석되었다. 비교적 수온의 영향이 적은 잠자리목과 관련성이 높게 분석된 혼합역은 냉수성 종 분포에 있어 상·하류지점과 차이가 있는 것으로 분석되었 다.
수온과 개체수, 종수, 군집지수, EPT-group, 냉수성 선호 종과의 상관성 분석 결과 (Table 4), 우점도, 개체수, EPT 비 율, 냉수성 선호종은 수온과의 음의 상관관계를 나타내는 것 으로 분석되었다. 반면 다양도, 균등도, 종수는 수온과 양의 상관관계를 나타내었으며, 풍부도는 유의한 양의 상관관계 를 나타내는 것으로 분석되었다. 또한 수온이 증가함에 따라 냉수성 선호종의 분포가 감소하는 것으로 분석되었다 (Fig. 7). 수온의 증가는 저서성 대형무척추동물의 종수를 증가시 키고, 민감종들의 분포를 제한하여 수생생물의 변화를 야기 하는 것으로 알려져 있다 (Jacobsen et al. 1997; Kwak et al. 2004; Burgmer et al. 2007). 본 연구 결과 온배수 유입으로 인한 수온의 변화는 냉수성종들의 다양한 분포를 제한하고, 수생태계에 교란을 야기시키고 있는 것으로 생각되며, 저서 성 대형무척추동물의 군집구조에 부정적인 영향을 미치고 있는 것으로 판단된다.
적 요
본 연구는 2015년 12월부터 2016년 2월까지 고온배수 배 출시설에서 온배수가 유입되는 경기도 포천시 영북면에 위 치한 부소천 일대에서 조사를 실시하였다. 조사지점은 온 배수가 유입되는 혼합역과 혼합역을 기준으로 상·하류 지 점을 선정하여 수질 분석 및 저서성 대형무척추동물의 종 조성, EPT-group 개체 비율, 군집 분석, 기능군 분석, 생물학 적 하천 건강성을 평가하였다. 조사기간 동안 출현한 저서 성 대형무척추동물은 총 3문 4강 11목 30과 50종 4,015개 체이며, 출현종 비율은 하루살이목이 6과 16종 (32.00)%, 개 체수 비율은 파리목이 1,848개체 (46.03%)로 가장 높게 나 타났다. 하루살이-강도래-날도래목 (EPT-group)이 조사기 간 동안 총 17과 32종 1,876개체 (46.72%)가 출현하였으 며, 1차에서 3차 조사기간 동안 혼합역에서 가장 낮은 출현 율을 나타냈다. 군집 분석 결과 혼합역에서 비교적 넓은 수 온대를 가지는 실지렁이 (Limnodrilus gotoi)와 깔따구 sp2. (Chironomidae sp.2)가 우점하는 것으로 분석되었으며, 네 점하루살이 (Ecdyonurus levis), 민하루살이 (Cincticostella levanidovae), 민강도래 KUa (Nemoura KUa) 등이 우점종으 로 분석된 상·하류 지점과 차이를 나타냈다. 섭식기능군 분 석 결과 상류와 혼합역에서 Gathering-Collectors, 하류에서 는 Filtering-Collectors의 비율이 가장 높은 것으로 분석되 었고, 상류에서 하류로 갈수록 Shredders의 비율이 감소하였 다. 서식기능군은 Clingers와 Burrowers가 대부분을 차지하 였으며, 상류에서 하류로 갈수록 Clingers의 비율이 높아지 는 것으로 분석되었다. 하천 건강성 평가 결과 평균 2.73 (± 0.41)으로 혼합역에서 가장 불량한 상태를 나타내는 것으로 분석되었다. 수온과 냉수성종의 상관성 분석 결과 수온이 증 가할수록 냉수성 선호종 및 EPT-group은 감소하는 것으로 분석되었다. 혼합역에 유입되는 고온의 방류수는 하천생태 계에 교란을 일으키는 요인으로써 저서성 대형무척추동물 군집과 냉수성 선호종 및 민감종의 분포에 부정적인 영향을 미치고 있는 것으로 판단된다.
