지금까지 환경부는 수질오염총량관리제 기간 동안 환경 기초시설확대, 오염원 관리 등 지속적인 유역중심의 예방적 수질관리 정책을 통해 BOD항목에 대한 수질개선을 위해 노력해 왔으며, 그 결과 지역별로 수질이 개선되거나 안정적으로 관리가 이뤄져 왔다. 하지만 인위적 수량통제로 인해 강 합류부나 하구언 등의 관리지점에서는 정체수역이 형성되어 하천 내･외부에서 유입되는 오염물질들에 의한 부영양화로 수량이 적은 저･갈수기에는 여전히 수질이 악화될 가능성이 높다. 일반적으로 하천의 유기물은 외부오염으로부터 기인되는 외부 기원 유기물과 내부의 영양염을 이용하여 성장하는 조류에 의해 기인되는 내부 생성 자생유기물로 나타난다(Kim and Lee, 2001a, 2001b).

낙동강에 합류되는 남강 역시 하류유역은 하상구배가 매우 완만하고, 주변의 도시와 농경지로부터 물 사용량이 많아 저･갈수기에는 유속이 더욱 느려져 이로 인해 수질악화가 가속화되어 수역의 영양단계가 점점 증가하는 부영양화현상이 자주 발생되고 있다(Kim, 2005). Jeong (2010)의 연구결과에 따르면 2006년~2008년까지 남강하류 유역에 Carlson에 의한 TSI(Trophic State Index)로 부영양화 정도를 조사한 결과, 연중 TSI 지수가 61이상으로 부영양상태 (Eutrophic State; TSI > 60)인 것으로 보고하였다. 일반적으로 유속이 빠른 하천에서는 외부에서 유입되는 외부기원 유기물이 수중 유기물의 주요한 공급원이 되며, 정체수역인 호소나 하천 하류 및 하구언 등에서는 수중생물의 1차 생산에 의한 자생 유기물이 주공급원이 된다(Kim et al., 2012; Shin, 2009).

부영양화 된 수역에서 BOD농도는 유역으로부터 발생한 유입 BOD 외에도 흐름이 느린 하천 내에서 미처리된 영양염을 이용한 식물플랑크톤의 생산에 의해 공급되는 자생 BOD의 합으로 나타난다. 지금까지 수중 유기물 농도 중 자생 유기물과 관련한 많은 연구가 이루어져 왔다. Fallon and Thomas (1979)에 따르면 ¹⁴C를 이용하여 조류 대발생(Algal bloom)에서 대발생을 이룬 남조류가 BOD₂₀값의 80%정도를 차지한다고 보고하였다. Kim and Lee (2001a)는 금강에서 식물플랑크톤의 성장이 COD상승에도 기여할 뿐만 아니라, BOD농도에도 기여하고 있음을 보고 하였다. 이는 물의 흐름이 느린 강 하류에서 부유성 식물 플랑크톤의 대량 증식이 일어날 경우, 이들이 내부 기원성 유기물 증가에 크게 기여할 수 있음을 의미한다.

하천수질은 수체 내에서 오염물질의 화학적, 생물학적, 물리적 작용의 결과로 나타나며, 이런 상호작용을 수식화하여 개발한 수질모형을 이용하여 하천의 수질관리를 위한 정책수립과 운영 그리고 평가를 하는데 활용한다. 하천수질 예측에 주로 이용되는 1차원 정상상태 하천수질 모형인 QUAL2E는 전 세계적으로 널리 적용되고 있지만, 조류의 생산 및 사멸에 따른 CBOD의 증가를 고려한 Bottel BOD 개념이 결여되어 있어 정체수역이 많은 우리나라 하천의 특성을 제대로 반영하지 못한다(Choi et al., 2008). 국립환경과학원에서는 이러한 QUAL2E 모형의 한계점을 보완한 QUALKO2 모형을 개발하였다(HRERC (2007)). QUALKO2 모형에서는 Bottle BOD 개념을 추가하여 조류의 생산에 의한 유기물 증가와 실제 정체수역이 많은 하천에서 일어날 수 있는 반응기작을 모의할 수 있도록 보완하였다(Jung, 2009; Kwon, 2009).

본 연구에서는 남강댐 하류유역 본류구간을 대상(Fig. 1)으로 식물플랑크톤 증식에 의한 수중유기물 기여도를 통계적 방법을 통해 분석하고, 아울러 수질예측의 정확도를 통계적으로 평가하는 것과 보정과 검정을 거친 QUALKO2 모형을 이용하여 향후 남강 하류 구간에서 빈번히 발생하는 부영양현상을 남강댐 방류량 증가 시나리오를 통한 수질예측으로 최적의 수질관리 방안을 제시하고자 한다.

남강유역은 남강댐을 경계로 남강 상류유역과 하류유역으로 구분되며 남강댐 하류유역은 남강댐을 기점으로 낙동강 합류 전까지의 남강 본류 77.6 km에 해당하며 유역 인구는 343,000명 정도이고 진주시와 진주상평공단 등 중･소도시 및 공업지역이 입지해 있어 오염물질 발생이 상당히 높다(Fig. 1). 낙동강에 대한 남강 하류유역의 BOD, TP 오염부하기여율은 각각 44.2%, 23.0%로서 매우 높은 편이다 (NRERC, 2009; Yu et al., 2012).

현재 남강 하류유역의 환경기초처리시설은 하수종말처리장이 6개소(총 시설용량 167,000 m³/day), 소규모 마을하수처리시설이 32개소(총 시설용량 1,792 m³/day)가 있으며, 그 외에도 폐수종말처리시설 5개소 등이 운영 중이다. Table 1에서 보는 바와 같이 2007년~2011년까지 최근 5년간 하수종 말처리장 방류수 농도를 조사한 결과, 진주하수처리장 BOD, TP 방류수 농도가 2007년 평균 6.9 mg/L, 0.931 mg/L에서 2011년 각각 평균 3.1 mg/L, 0.380 mg/L로 대부분의 처리장 방류수 농도는 목표수질 달성을 위해 지속적으로 개선되었음을 알 수 있다.

본 연구에서는 남강 하류유역에서 조사된 2008년~2012년까지 총 5년간의 모니터링 자료를 분석에 사용하였다. 먼저 수질 자료는 국립환경과학원 낙동강물환경연구소 측정자료를 이용하였으며, 분석항목은 BOD, COD_Mn, TOC, TN, TP, Chl-a, Water temperature, pH, DO, Conductivity, Ammonia nitrogen, Nitrate nitrogen, Phosphorate, outflow 등 총 14개 항목을 사용하였고, 유량자료는 국토부(환경부) 정암수위국 (남강D)과 거룡강수위국(남강E) 지점의 각 연도별 유출분석을 통해 산출된 일유량자료를 분석에 사용하였다.

강우량자료는 남강유역에 구축된 기상청 티센망의 티센계수값을 이용한 티센평균 일강우량을 분석에 사용하였고, 남강댐 방류량 자료는 국가수자원관리종합정보 인터넷홈페이지(www.wamis.go.kr)에서 수집된 남강댐 본류 일방류량 자료를 분석하여 사용하였다.

상관성 분석은 수질인자간 연관성이 있는 수질항목을 파악하기 위하여, 인자간의 상관관계를 분석하여 상호 연관성 높은 인자와 항목을 추출하는 방법이다. 상관계수의 크기가 0에 가까우면 두 변수사이에 상관성이 거의 없으며, 크기가 1에 가까우면 두 변수 사이의 상관성이 매우 높다고 해석한다. 다음으로 상관계수의 부호가 양(Positive)이면 두 변수가 같은 쪽으로 변화하고 음(Negative)이면 두 변수는 반대로 변화한다. 요인분석은 주성분분석법과 최소 고유값 기준(고유값 1이상)으로 요인의 수를 결정하였으며, 각 요인의 특성을 알고자 할 때 유용한 베리맥스(Varimax) 회전 방식을 적용하였다.

모델의 정확도를 평가하기 위해서 Table 2에 제시된 통계값을 사용하였다. 결정계수(Coefficient of determination; R²)는 예측값과 실측값의 선형 회귀관계의 적합도를 결정하는 계수로써 0과 1사이의 값을 가지며 1에 가까울수록 적합도 및 모형의 정확도가 높다. 평균제곱근오차(Root mean square error; RMSE)는 예측결과의 평균오차를 나타내는 지표로써 그 값이 0에 가까울수록 즉 값이 작을수록 모형의 정확도가 높고 모형의 성능도 우수하다고 평가할 수 있다 (Choi. et al., 2008).

남강댐의 주요 시설물은 높이 43 m, 길이 1,126 m의 콘크리트 표면차수벽 석괴댐인 본댐과 문비 3문을 갖춘 여수로, 시설용량 14,000 kW의 발전소, 그리고 홍수시 사천만으로 홍수를 방류하기 위한 높이 31 m, 길이 258 m의 콘크리트 중력식 제수문 구조물 등으로 구성된다. 남강댐의 일반 제원으로 유역면적은 2,285.0 km², 저수면적 28.2 km², 유역 연평균유입량 64.4 m³/s, 계획홍수위 46.0 EL.m, 상시만수위 41.0 EL.m, 저수위 32 EL.m이다.

남강댐을 취수원으로 광역상수도 사업을 통하여 용수를 공급받는 곳은 통영시, 사천시, 고성군, 진주시, 거제시, 고성군, 하동군, 남해군 등이다. Table 3과 같이 기본계획공급량에 따르면 남강댐은 총 573.3 백만톤 중에서 생･공업용수 224.4 백만톤, 농업용수 226.8 백만톤, 하천유지용수 122.1 백만톤을 공급하는 것으로 계획되어 있다.

물공급량(Water yield)은 특정 근원으로부터 획득할 수 있는 물의 양을 말한다. 대개 침투나 증발이 없다고 가정하면 하천으로부터 최대 공급량은 연평균 유출량과 같게 된다. 이에 따라 남강댐에서 공급가능한 최대 물공급량은 기본계획상으로는 64.4 m³/s이고, Table 4의 최근 5년간 유입량 분석결과로는 73.5 m³/s (25.2 ~ 106.2 m³/s)까지 공급할 수 있는 것으로 분석되었다.

Table 4에 요약한 남강하류 유역 강우량을 살펴보면 2012년 기준으로 6월 하순부터 8월 중순까지의 여름철 장마기와 9월에 강우가 집중되어 연강우량(1,788.1 mm)의 약 70.7% (1,264.3 mm)가 이 기간 동안 내렸다. 이에 따라 남강댐 방류량도 증가하였고, 그 결과 하류의 유량도 이때 가장 많았다. 한편 2011년 12월과 2012년 1월 ~ 2월까지의 강우량은 매우 적어 가뭄현상을 보였으며 2009년과 2010년 봄철에는 강우 빈도가 산발적으로 많았으나 그 양은 많지 않아 남강댐 방류량과 유량은 적은 양이 방류되고 실측되었다.

남강댐은 남강본류와 사천만방수로를 통하여 방류하고 있다. 남강댐은 농업용수 226.8 백만톤(7.2 m³/s) 및 하천유지용수 122.1 백만톤(3.9m³/s)을 연간 348.9 백만톤(11.1 m³/year)을 공급하는 것으로 계획되어 있다. Table 5와 같이 남강댐이 본류로 방류하는 물의 양은 발전방류량과 여수로방류량의 합이며, 일반적으로 발전방류는 농업용수와 하천유지용수를 공급하며, 여수로방류는 홍수피해저감을 위하여 방류한다. Table 5의 남강댐 본류 연간 방류량을 살펴보면, 발전방류량의 평균은 915.0 백만톤(29.0 m³/s)으로 농업용수와 하천유지용수 연간 계획공급량 348.9 백만톤(11.1 m³/s)을 훨씬 넘는 것으로 나타났다. 하지만 Table 4의 월별 방류량을 살펴보면 위 계획공급량 이하의 월이 1월~ 4월과 10월~ 12월에 집중되어 있는 것을 알 수 있다. 여수로 방류는 주로 비가 많이 내리는 우기(6월~ 9월)에 이루어지기 때문에 그 외 월의 방류량은 모두 발전방류량에 해당하는 것으로 남강하류 수질개선을 위해서는 추가적인 댐 방류량 증가가 필요할 것으로 판단된다.

하천의 수질은 물리적 특성, 즉 수리･수문학적 특성에 의해 많은 영향을 받는다. 오염총량제 2단계 기준유량은 1995년 ~ 2004년까지 통합모형화 기법으로 구축한 SWAT 모형을 이용하여 자연유량 및 댐 방류로 인한 효과와 물수지 요소에 의한 영향 등을 고려하여 10년 평균 저수량(Q₂₇₅)을 기준으로 낙동강 수계 41개 오염총량관리 단위유역별로 기준 유량이 산정되었으며, 275번째 해당하는 유량 즉 365일 중 275일 이상을 유지하는 최소유량으로 초과확률(Exceedance Probability) 75%에 해당하는 유량을 의미한다. 저수량을 기준 유량으로 선정한 것은 저수량보다 유량이 풍부한 경우에는 목표수질을 만족할 수 있다는 의미로 볼 수 있다(GNDI, 2010).

Table 6과 같이 지난 5년간(2008년 ~ 2012년)의 실측유량 자료를 바탕으로 유황분석을 실시한 결과 남강D(남강2A) 지점의 저수량은 Q₂₇₅ = 7.3m³/s이며, 남강E(남강4-1) 지점의 저수량은 Q₂₇₅ = 13.8 m³/s으로 나타났다. 오염원으로부터 배출되는 오염물질의 양이 동일하다고 가정한다면, 유량이 저수량(Q₂₇₅) 이하일 때 측정한 수질은 높게 나타날 가능성이 있다. 실제로 상기 조건에서 저수량 이하에서 측정된 수질을 선별하여 조사한 결과 남강D 지점의 평균 수질농도는 BOD = 3.9 mg/L, 남강E 지점의 평균 수질농도는 BOD = 4.2 mg/L로 목표수질 남강D(BOD = 2.5 mg/L)와 남강E(BOD = 3.1mg/L)를 초과하는 것으로 나타났다. 또한 저수기에 해당하는 평균유속을 선별하여 조사한 결과 남강댐 하류평균유속은 0.19 ㎧ 정도의 저유속으로 나타났다. 또한 본 연구에서는 남강댐 하류로의 하상경사 및 유속분포 그리고 유하시간을 산정하고자 HEC-RAS 모형을 이용하였으며, 오염총량제 2단계 기본계획에 적용된 단면자료를 활용하였다.

경계조건(Boundary Condition)은 앞서 수집･분석된 남강댐 방류량 자료 및 수위국 지점들에서 분석된 일평균유량을 통해 산출된 유황분석 결과값을 상류단 경계값으로 사용하였으며, 하류단 경계조건은 Nomal Depths(=0.000271)을 적용하였다. 각 단위유역별 산출된 최종 유황분석결과는 Table 6과 같다.

HEC-RAS 모형을 이용하여 남강댐 하류의 유속을 모의한 결과 저수기(Q₂₇₅)기준 남강D(남강2A) 부근은 0.10 ㎧, 남강 E(남강4-1) 부근은 0.16 ㎧로 실측한 자료와 매우 근사한 값을 나타내었다. 산정된 유하시간은 5.6 day로 2009년 안동조정지댐으로부터 하구까지 물이 도달하는데 대략 12 day이 소요되는 것으로 산정된 것을 감안할 때 상당히 긴 시간으로 나타났다(Fig. 2 참조)(Choi et al., 2008).

Fig. 3에서 보듯 남강댐 하류로 내려 갈수록 하상 경사가 매우 완만하며 다수의 수중보와 교량 등 인공 구조물과 만곡흐름의 특징으로 저유속으로 인한 강의 유하시간이 증가하는 것으로 분석되었다. 하천에서의 유하시간 증가는 곧 체류시간의 증가를 의미한다. 따라서 남강댐의 방류량이 적은 늦가을부터 다음해 봄까지의 갈수기에 남강 하류 수역의 유속도 매우 느려져 이런 조건은 남강 하류 수역에서 식물플랑크톤의 증식을 촉진 할 것으로 예상된다.

남강댐 직하류인 남강C(남강)의 수질은 2008년~ 2012년까지 거의 변화가 없었고, Fig. 4에 나타난 남강댐 하류 수질측정망 지점에서의 최근 5년간 연평균 분포를 보면 남강하류로 갈수록 BOD, TP, TN, Chl-a의 농도가 높아지는 것을 알 수가 있다.

특히, Chl-a의 경우 단위유역 남강D(남강2A)지점을 지나면서 급격히 증가하는 경향을 보이는 것으로 나타났다. 이러한 현상은 지형적으로 경사가 완만해지는 중･하류부에서 강의 흐름(유속)이 느려지면서 체류시간이 길어져 미처리된 유기오염물질과 질소, 인 등 과량의 영양염류가 집적되면서 수체 내 부유성 식물플랑크톤이 대량 증식하는 현상을 초래한 것으로 추정된다.

남강댐 하류 단위유역 남강D(남강2A)지점의 BOD농도는 2009년 평균 3.6 mg/L에서 2012년 평균 1.9 mg/L로 남강E(남강4-1)지점의 BOD농도는 2009년 평균 4.3 mg/L에서 2012년 평균 2.4 mg/L로 수질이 개선되어진 것으로 조사되었다. 하지만 Fig. 5에 도시한 바와 같이 월평균 BOD값을 살펴보면 두 지점 모두 대체로 1월~ 5월까지 높은 BOD 값을 보여 봄철 유기물 오염도가 높은 것으로 나타났다. 또한 강우가 집중되는 6월 하순에서 8월 중순 풍수기와 9월에서 11월 평수기를 제외하고 12월부터 저･갈수기인 봄철 5월까지 지속적으로 수질이 악화되어 목표수질 남강D(BOD = 2.5 mg/L)와 남강E(BOD = 3.1 mg/L)를 초과하는 것으로 나타났다. 일반적으로 강우량이 많을 때 하천은 세척효과로 정화되지만 선행무강우일수가 길면 하천은 대부분 하･폐수로 이루어지면서 오염물질이 하천에 집적되어 수질오염은 심화될 수밖에 없다.

Fig. 4에서 BOD와 Chl-a의 경우 상류에서 하류로의 수질분포 특성을 보았을 때 남강댐 하류 수역은 농도가 비교적 낮은 구간 (a)상류(남강, 남강1, 남강2), 농도가 증가하는 구간 (b)중류(남강2A, 남강3), 높은 농도가 유지되는 구간 (c)하류(남강4, 남강4-1) 구간으로 구분되어진다. 이를 바탕으로 수질조사 기간 동안 각 구간별 이화학적 수질값의 최대, 최소값 및 평균값을 Box Plot 형태의 통계값으로 표시하여 나타내었다(Fig. 6 참조). 유기물 오염지표인 BOD, COD_Mn, TOC 분포 특성을 살펴보면 BOD의 중앙값은 1.0 ~ 2.6 mg/L, 제1사분위수(25%)는 0.8 ~ 1.7 mg/L, 제3사분위수(75%)는 1.3 ~ 4.2mg/L의 농도 범위를 보이고 있으며, 중 하류로 내려가면서 농도가 증가하고, 농도 변동 폭도 커지는 것으로 나타났다. COD_Mn, TOC의 경우 상류에서 하류까지의 농도 변화는 BOD와 비슷한 경향을 보이며 농도분포가 중앙값에 밀집된 분포를 보이고 있다. TP의 제1사분위수는 0.01 ~ 0.072 mg/L(7.2배), 제3사분위수는 0.03 ~ 0.135 mg/L(4.5배), 중앙값은 0.02 ~ 0.1 mg/L(5배) 증가한 농도를 보이고 있다. TN은 사분위범위가 중앙값을 기준으로 저 농도로 분포되어 있는 밀집형태를 보이고 있다. 식물플랑크톤 현존량 지표인 Chl-a의 최소 10 %percentile 값과 최대 90 %percentile 값은 상류에서 하류로 각각 2.40 ~ 12.87 μg/L, 5.44 ~ 76.83 μg/L, 6.10 ~ 103.6 μg/L로 다량의 영양염 유입과 완만한 하상구배의 느린 흐름으로 Chl-a 농도의 증가 범위가 연중 2.5 ~ 185.8 μg/L로 매우 큰 변동 폭을 보이고 있으며, 중앙값을 기준으로 흩어짐 정도가 매우 크게 나타났다.

남강댐 하류 수역에서의 수질관측자료의 통계분석을 위해 SPSS 12.0 통계 프로그램을 사용하여 항목 간 상관관계를 분석하였다. 통계분석을 실시하기 전에 주요 수질항목이 정규분포를 따르는지를 검토하였으며, 정규분포를 따르지 않는 항목에 대해서는 적절한 자료의 변환을 검토하였다. 하지만 원 자료의 관측치들이 정규분포를 따르는 경우에는 보다 더 의미가 있겠으나, 이러한 정규성의 가정이 반드시 필요한 것은 아니다(Kim et al., 2007). 그러므로 본 연구에서는 자료의 변환 없이 원 자료를 그대로 표준화 시킨 후 분석하였다.

남강하류 유역에서 조사된 인자별 상관분석 결과를 Table 7에 (a)상류, (b)중류, (c)하류 순으로 나타내었다. 상류 구간에서 WT(r = 0.421), BOD(r = 0.340)가 Chl-a 농도와 작지만 통계적으로 유의한 양의 상관관계(p < 0.01)를 나타냈다. 중류 구간과 하류구간에서는 pH(r = 0.547), BOD(r = 0.718), COD_Mn(r = 0.766), TOC(r = 0.606) 등이 Chl-a 농도와 높은 양의 상관관계(p < 0.01)를 가짐으로써 조류성장에 따른 자생 유기오염물질과의 관련성을 추정해 볼 수 있다. 또한 댐(본류)방류량과 BOD(r = 0.388)와는 작지만 통계적으로 유의한 음의 상관관계(p = 0.01)로 나타나 방류량에 의한 오염원 제어 효과가 있는 것으로 판단된다. pH와 Chl-a 농도가 높은 상관성을 보인 것은 조류의 성장과 관련한 pH의 변화로 설명할 수 있다. 조류의 경우 광합성이 활발하면 수중의 용존산소를 증가시키는 한편 조류의 증식으로 인한 CO₂의 소모에 따라 탄산염 시스템(Carbonate system)의 역반응으로 H⁺이온이 감소하면서 pH를 상승시키는 결과로 나타난다. 즉, 대기 중 CO₂가 수중으로 녹아들면, 탄산(H₂CO₃)이 형성되며, 이는 다시 수소이온(H⁺)과 중탄산염(HCO₃^-)으로 분해되는데, 이 중 중탄산염은 수중식물(조류)의 광합성 시 CO₂흡수의 주된 형태이므로 식물플랑크톤이 증식하면서 광합성 양이 증가할수록 중탄산염 흡수가 증가되어 수소(H⁺)이온의 양은 계속 감소하게 된다. 이는 결국 수체의 pH 증가를 유발하게 된다. 또한 BOD, COD_Mn, TOC 등이 Chl-a와 높은 인자 부하량을 가짐으로써 조류성장에 따른 자생 유기오염물질의 발생을 나타내고 있다. 그리고 전기전도도와 NH₃-N(r = 0.519), NO₃-N(r = 0.503)가 양의 상관관계(p < 0.01)로 나타나 수중에서 유기물과 함께 전해질로 작용함을 알 수 있다. 한편 NO₃-N와 PO₄-P는 음의 관계를 보여 식물플랑크톤 증식과정에서 우선적으로 소비되는 이들 무기영양염류는 식물플랑크톤의 증식 과정에서 소모되어 음의 상관을 나타낸 것으로 추정된다. 또한 남강댐 하류 구간에서 영양염류 농도와 Chl-a가 뚜렷한 상관성을 보이지 않고 있다. 이는 하류 수역에 영양염류 농도가 식물플랑크톤이 증식하기에 이미 충분한 상태이기 때문인 것으로 판단된다. 조류발생의 필요조건인 인 농도는 미국 EPA 기준으로 0.02 mg/L, OECD 기준으로 0.035 mg/L로(Kim, 2012) 앞서 구간에서 조사된 인 농도는 (a)상류(0.028 ~ 0.107 mg/L), (b)중류(0.073 ~ 0.152 mg/L), (c)하류(0.071 ~ 0.128 mg/L)로 이미 충분한 상태이다. 한편, 상류구간에서 하류구간까지 TN/TP의 질량비(TN:TP mass ratio)를 보면 각 지점에서 조사기간 평균 54.7, 37.4, 38.1 (범위: 6.2 ~ 757.3)의 비를 보였고, 특히 상대적으로 TP 농도가 낮았던 하천구간에서 더 높은 비를 보여 전 지점에서 모두 인이 제한영양염인 것으로 조사되었다(Forsberg and Ryding, 1980). 일반적으로 TN과 TP의 질량비가 16 이하일 경우 총질소가 제한요인으로 작용하고, 29 이상일 경우 TN에 비하여 총인의 부족으로 1차 생산력에 영향을 미칠 수 있는 잠재성을 갖고 있는 것으로 판단되며 세 구간 모두 TN과 TP의 질량비가 평균 29 이상으로 인이 제한(P - limitation) 요인으로 나타났다(Park et al., 2013). 또한, 조사기간 평균 값으로 볼 때 각 지점에서 TP농도는 0.06 mg/L, 0.11 mg/L, 0.1 mg/L로 부영양 기준을 초과하였다.

남강댐 하류유역 총 세 구간에서 인(P)이 제한영향인자인 것으로 조사되었고, 식물플랑크톤 증식에 영향을 미치는 인자인 것으로 판단이 되나, 강우 및 상류 댐의 방류량 증가로 인한 유속의 증가(체류시간 감소) 또한 중요한 영향을 미치는 것으로 추정된다. 조류의 발생원인은 영양염류의 과다 이외에도 수리학적 인자인 유속, 체류시간, 강우량과 이화학적 인자인 수온, 태양광량, 영양염류, 그리고 생물학적 인자들이 복합적으로 작용하여 발생한다. 따라서 저･갈수기시에 남강댐 방류량이 하천유지유량인 발전방류량이 줄어들어 유속이 느려지면, 하류로는 상류에서 흘러내려온 하･폐수가 집적되어 수질이 악화되는 현상을 초래할 것으로 판단된다. 앞서 언급한 수질 특성분석에서 월별 BOD농도를 분석한 결과 봄철 유기물 오염도가 높은 것으로 나타났는데 이는 상기시기와 일치하는 현상으로 볼 수 있다.

한편 월별 방류량에 따른 BOD와 Chl-a의 상관관계 분석을 실시한 결과 1월(r = 0.840)(r = 0.783)과 2월(r = 0.833)(r = 0.757), 5월(r = 0.751)(r = 0.658)이 가장 높은 음의 상관관계(p < 0.01)로 나타났으며, 3월(r = 0.586) (r = 0.430)과 4월(r = 0.509)(r = 0.472)이 그 다음 순으로 나타났다(Table 8 참조). 강우가 집중되고 유량이 많아 유속이 빠른 7월과 8월은 음의 관계이지만 유의성은 작은 것으로 분석되었으며, 강우 후 수체가 안정화된 9월과 10월은 유의성이 없는 것으로 나타났다. 대체로 하천에 유량이 작고 유기물 오염도가 높은 시기에 상관성이 매우 높은 것으로 나타났다. 위 결과로 볼 때 남강댐 방류량 증가에 따라 오염농도가 감소하는 것은 분명한 것으로 판단된다. 따라서 남강댐 중･하류 구간의 유기물 오염관리를 위해서는 봄철 수중에서 식물플랑크톤의 증식과 더불어 오염부하량의 제어가 매우 중요하며 이를 위해서는 남강댐의 수문조작을 통한 적정한 유량에 의한 유속을 유지해 주는 것이 효과적일 것으로 판단된다. 미국 Patuxent River 하구의 연구에서 겨울철 조류 발생이 질소(N) 제한조건과 관련이 있으며, 유속증가에 따라 조류 발생이 사라짐을 보고하였다(Sohn et al., 2013).

분류된 세 구간은 서로 다른 수질특성을 나타내므로 각 구간의 수질특성을 파악하기 위하여 주성분 분석 및 요인 분석을 실시하여 그 결과를 Table 9와 10에 나타내었다(Bernard et al., 2004; Kim et al., 2007; Shin et al., 1997; Yoo, 2002). 주성분 수를 결정하기 위해서 고유치(eignvalue)가 1.0 이상인 값을 갖는 주성분 축만을 고려하여 상류 4개, 중･하류 3개의 요인을 추출하였다. 주성분 전체 수질 변동에 대하여 (a)상류 78.531%, (b)중류 78.394%, (c)하류 75.669%를 설명하고 있다.

상류 구간에서 제 1요인은 31.161%, 제 2요인은 52.722%, 제 3인은 68.591%, 제 4요인은 78.531%를 기여하고 있고, 제 1요인은 NH₃-N, Cond, TN로 분류되었고, 제 2요인은 수온(WT), 제 3요인은 BOD로 분류되었다. 중류 구간에서 제 1요인은 27.000%, 제 2요인은 53.461%, 제 3인은 78.394%를 기여하고 있고, 제 1요인은 NH₃-N, 제 2요인은 Chl-a, BOD, COD_Mn, pH, 제 3요인은 TP로 분류되었다. 하류 구간에서 제 1요인은 31.835%, 제 2요인은 58.616%, 제 3인은 75.669%를 기여하고 있고, 제 1요인은 NO₃-N, TN, DO, 제 2요인은 Chl-a, CODMn, BOD, TOC 제 3요인은 PO₄-P, TP로 분류되었다.

상류 구간에서는 암모니아성 질소, 총질소와 전기전도도가 제 1요인으로 분류되었고, 수온이 제 2요인으로 유기물 요인인 BOD가 제 3요인으로 분류되어 진주 도시 및 공업지역 오염특성을 잘 나타내는 것으로 보이며, 중류와 하류구간은 암모니아성 질소가 제 1요인으로 공통 분류되었고, BOD, COD_Mn, Chl-a가 제 2요인으로 공통 분류되었다. 이와 같은 결과로 볼 때 남강 중･하류구간의 수질은 조류 증식에 따른 자생 유기오염물질의 발생이 수질 오염의 요인인 것을 확인할 수 있었다. 이는 상류 도시지역의 생활하수가 흘러내려오면서 하류지역에 도달하여 누적됨에 따라 도시특성과 농가의 축산폐수 및 토지이용에 따른 유기오염물질과 영양염류 유출로 오염도가 극대화됨에 수질관리에 있어서 가장 중요한 변수임을 확인할 수 있었다.

본 연구에서는 2단계(2011년 ~ 2015년) 오염총량관리 기본계획수립에서 사용된 바 있는 QUALKO2 모형을 이용하여 남강댐 방류량 변화에 따른 남강댐 하류 지역의 수질개선 효과를 분석하였다. 먼저 신뢰성 있는 수질모형을 구축하기 위해서는 대상하천의 유속 및 수심과 같은 수리학적 특징이 모형에 잘 구현되어야 한다. 따라서 Reach별 수리학적 계수는 낙동강 오염총량관리 2단계 기본계획에 적용된 단면자료를 이용하여 V = aQ^b, H = αQ^β의 관계식으로부터 도출된 수리계수 값을 적용하였다. Reach는 수리분석 결과에 근거하여 하천 내 수리특성이 유사한 구간을 일정하게 구분하여 14개 구간으로 나누었으며, 소구간 Element 길이는 1 km 등간격으로 구분하였다. 남강 본류에 유입되는 하･폐수처리시설과 취･양수장 및 지류하천은 위치를 확인하여 모두 점오염원 형태로 모형에 적용하였으며, 그 외에 비점오염원으로 각 소유역의 유달부하율(유달율)을 고려하여 수체에 유입하는 것으로 모형을 구성하였다(GyeongSangNamDo, 2009; Jung et al., 2013). 유달율은 2단계 오염총량관리 기본계획수립에서 산정된 자료를 이용하여 각소유역의 하단부에서 점오염원으로 유입되는 형태로 모형에 적용하였다.

모형에서 반응계수의 결정은 실측을 통해 산정하여야 하나 이는 불가능할 뿐만 아니라 많은 오차가 포함될 가능성이 크다. 따라서 대부분 수질모형에서 반응계수 추정은 자연현상과 수학적인 해석간의 차이가 최소가 될 때까지 시행착오법을 반복하여 최적의 반응계수를 결정하게 된다. 따라서 모형의 보정 수질은 기준유량(Q₂₇₅)과 가장 유사한 시기인 2008년 3월 환경부 수질 측정망 자료의 평균값을 이용하였으며 BOD, TN, TP, Chl-a 항목에 대해 수질모형인자의 민감도 분석을 통한 보정을 실시하였다. 민감도 분석은 수질모형인자를 10 ~ 50% 비율로 변화시키면서 모형을 실행 마지막 Element에서 BOD, TN, TP, Chl-a 등 대부분의 수질항목에서 영향정도가 크게 나타난 변수를 중심으로 실측값을 근거로 보정하였다. 민감도 분석 결과 K₁(BOD 분해율 계수), K₃(침전에 의한 BOD 제거율), σ₁(조류의 침전율), μ_max(조류 최대 성장률), ρ(조류의 호흡율), σ₂(용존성 P의 하상 용출율, 용존인의 저질 용출율) 등의 수질모형인자에 민감하게 반응하였다(Fig. 7과 Table 11 참조) 이처럼 조류의 성장 및 사멸에 민감하게 반응하는 이유는 QUALKO2 모형이 조류의 성장에 의한 유기물의 내부 재생산을 적절하게 잘 모의하고 있기 때문으로 판단된다.

Fig. 8에서 보듯이 BOD와 Chl-a의 경우 하류로 가면서 농도가 증가하는 경향을 보이는데, 이는 남강 중･하류에서 조류의 성장에 의한 유기물질의 내부 생산이 원인인 것으로 추정된다. 앞서 언급한 Fig. 5와 Fig. 7에 도시한 바와 같은 경향으로 모형 예측에서도 이러한 현상을 잘 모의하고 있다.

모형의 검증은 보정된 수질모형의 적정 사용여부를 판단하는 과정으로, 보정된 모형이 수질예측에 적절한지를 알아보기 위해 보정에 사용한 매개변수를 그대로 고정하고, 2008년 4월 측정한 유량･수질 자료를 이용하여 모형을 검증하였다. Fig. 9에 도시한 바와 같이 BOD, TN, TP, Chl-a 모든 항목에서 수질모형 예측값이 실측값을 안정적으로 잘 모의하고 있다.

모형의 검･보정 결과에 대한 정확도를 평가하기 위해 통계를 이용한 오차 분석을 실시하였다. BOD, TN, TP, Chl-a 항목의 모형 예측 결과의 정확도에 대한 통계분석 결과 모든 항목에서 양호한 결과를 나타내었다. Chl-a의 RMSE 값이 다소 높은 이유는 다른 항목과 달리 단위를 μg/L로 설정하였기 때문이다(Table 12 참조).

최종 보정 및 검정절차를 거친 남강댐 하류유역 QUALKO2 모형을 이용하여 남강댐 방류량 변화에 대한 시나리오를 설정 목표수질 달성여부를 평가하였다. 시나리오 1은 2008년부터 2012년까지의 남강댐 1월 평균 방류량 인 8.6 m³/sec를 Head water 유량 값으로 설정하였고, 시나리오 2, 3, 4, 5는 시나리오 1의 1.5배, 2배, 2.5배, 3배씩 유량을 증가시키면서 모의(Simulation)하였다.

Fig. 10에서 보듯이 방류량이 증가할수록 남강댐 방류수에 의한 희석효과로 남강댐 하류 수역의 BOD농도가 감소하는 것으로 나타났으며, 남강D 단위유역에서는 17.2 m³/sec인 경우에 목표수질인 2.5 mg/L을 만족하는 것으로 나타났고, 남강E 단위유역에서는 21.5 m³/sec인 경우에 목표수질인 3.1 mg/L을 만족하는 것으로 나타났다.

Table 13에서 방류량이 8.6 m³/sec에서 25.8 m³/sec로 증가한 경우 BOD농도가 각각 13%, 13%, 9%, 6% 감소하는 것으로 나타났고, 특히 영양염류 지표인 TP, TN 보다 BOD, Chl-a 같은 유기 오염물이 방류량 증가에 따른 감소정도가 높은 것으로 나타났다. 한편 시나리오별로 유량을 증가했을 때 감소하는 정도가 점차적으로 줄어들었다. 이는 곧 남강댐 하류 수질을 개선하기위해 남강댐 방류량을 무한정으로 늘리는 것은 효과적이지 않은 것으로 생각된다. 하지만 남강댐 방류량이 어느 정도 양의 수준에서는 수질개선 효과에 도움이 되는 것으로 판단된다.

이러한 연구결과로 볼 때 수문조작에 의한 유량 조절은 연중 특정시기에 편중된 강우량으로 유량변동이 매우 큰 우리나라 강우 특성상 한계가 있으며, 남강댐 하류 수역의 영양염류 농도가 식물플랑크톤 증식의 제한인자로 작용하지 않을 만큼 높다는 점과 식물플랑크톤에 의한 내부 유기물 생산이 하류 수역의 유기물 오염의 주요 기원으로 나타난 것을 감안할 때 근본적으로 식물플랑크톤 증식을 제어하기 위한 영양염류 저감 대책을 장기적으로 마련해야 할 것으로 판단된다.

낙동강수계에 유입되는 지천 중에 가장 큰 오염원을 배출하는 남강댐 하류유역 수질개선을 위해 1차적으로 중점 관리 대상지류 선정 연구(Jung et al., 2013)를 통해 수질관리 필요 지류들을 조사하여 선정하였고, 후속 연구로 남강 댐 하류유역 본류구간을 대상으로 식물플랑크톤 증식에 의한 수중유기물 기여도를 통계분석을 통해 평가하고, 아울러 수질예측의 정확도를 통계적으로 평가하는 것과 보정과 검정을 거친 QUALKO2 모형을 이용하여 향후 남강 하류 구간에서 빈번히 발생하는 부영양현상을 남강댐 방류량 증가 시나리오를 통한 수질예측으로 최적의 수질관리 방안을 분석한 결과 다음과 같은 주요한 연구 성과를 얻을 수 있었다.

1) 강우량 분석결과 총강우량의 약 60 ~ 70%가 우기인 6월~ 9월에 집중되고, 그 외 월은 평균강우량 이하의 강우가 내리는 것으로 분석되었다. 총강우량이 많으면 수질이 개선되고, 적으면 수질이 악화되는 양상을 나타내지만, 우기에 집중되는 것보다는 상대적으로 월별로 균등한 양이 일정하게 내리는 것이 수질관리에 유리할 것으로 판단된다. 강우는 자연현상으로 인위적 제어가 불가능한 사상이므로 남강댐 하류유역 수질개선을 위해서는 강우량이 상대적으로 부족한 기간의 댐(본류) 방류량 증가가 효과적일 것으로 판단된다.2) 남강댐은 남강본류와 사천만방수로를 통하여 방류하고 있다. 이중 남강하류 수질개선은 남강댐 본류방류량과 관련이 있다. 남강댐 본류 방류량은 발전방류량과 여수로방류량의 합이다. 이중 발전방류량의 평균은 915.0 백만톤(29.0 m3/s)으로 농업용수와 하천유지용수 연간 계획 공급량 348.9 백만톤(11.1 m3/s)을 훨씬 넘는 것으로 나타났다. 하지만 월별 방류량을 살펴보면 위 계획공급량 이하의 월이 1월~ 4월과 10월~ 12월에 집중되어 있는 것을 알 수 있다. 여수로 방류는 주로 비가 많이 내리는 우기(6월~ 9월)에 이루어지기 때문에 그 외 월의 방류량은 모두 발전방류량에 해당하는 것으로 남강하류 수질개선을 위해서는 발전방류량의 증가가 필요할 것으로 판단된다.3) 남강댐 직하류인 남강C(남강)의 수질은 2008년~ 2012년까지 거의 변화가 없는 것으로 나타났으나, 남강댐 하류 수질측정망 지점에서의 최근 5년간 연평균 분포를 보면 남강 하류로 갈수록 BOD, TP, TN, Chl-a의 농도가 높아지는 것을 알 수가 있다. 특히, Chl-a의 경우 단위유역 남강D(남강2A)지점을 지나면서 급격히 증가하는 경향을 보이는 것으로 나타났다. 또한 남강D(남강2A)와 남강 말단 조사지점인 남강E(남강4-1) 월평균 BOD값을 살펴본 결과 1월 ~ 5월까지 높은 BOD 값을 보여 봄철 유기물 오염도가 높은 것으로 나타났다.4) 구간별로 수질항목간의 상관관계를 분석한 결과 경사가 완만해져 유속이 느려지기 시작하는 전이대 구간인 중류와 하류 구간에서 BOD와 Chl-a 농도간의 상관관계가 높게 나타나 이 구간에서 조류증식에 의한 물질대사 요인인 Chl-a 즉 내부 유기물 생산량이 이 수역의 유기물 오염의 주요기원으로 추정되며, 영양염류는 전 구간에서 부영양화 기준을 초과하여 남강 하류 수역에서 식물플랑크톤 증식에 제한 요인은 아닌 것으로 조사되었다. 이외에도 식물플랑크톤 증식에 주요한 영향을 미치는 인자로는 유속과 직접적으로 관련이 있는 강우사상과 남강댐 방류량 그리고 그에 따른 수리학적 체류시간으로 나타났다. 또한 각 구간의 수질요인을 분석한 결과 중･하류구간에서 BOD, CODMn, Chl-a 농도간의 높은 양의 인자 부하량을 보여 남강 수역의 유기물 오염관리를 위해서는 수중에서 식물플랑크톤의 증식을 제어하는 것이 매우 중요한 것으로 조사되었다.5) QUALKO2 모형을 이용하여 댐 방류량 변화에 따른 남강 하류 지역의 BOD 수질개선을 위한 필요유량 산정 결과 남강D(남강2A) 지점은 남강댐 방류량이 17.2 m3/sec이상의 경우 목표수질인 2.5 mg/L을 만족하는 것으로 나타났고, 남강E(남강4-1) 지점은 21.5 m3/sec이상의 경우 목표수질인 3.1 mg/L을 만족하는 것으로 나타났다. 따라서 남강 하류의 수질개선을 위한 최대 방류량은 21.5 m3/sec이상이 꾸준히 방류되어야 수질개선 효과가 있는 것으로 분석되었다.6) 남강댐에서 공급 가능한 최대 물공급량은 기본계획상으로는 64.4 m3/s이고, 최근 5년간의 유입량 분석결과 73.5m3/s(25.2 ~ 106.2 m3/s)까지 공급할 수 있는 것으로 분석되었다. 이와 같은 결과로 볼 때 추가적인 댐 방류량 확보에 충분한 여력이 있는 것으로 판단되며, 남강댐 하류유역 수질개선을 위해 댐 방류량 증가가 필요할 것으로 판단된다.