인구의 증가와 각종 산업이 급격히 발달되고 동시에 하천의 오염으로 인한 문제가 심각해지고 있다. 이처럼 산업폐수가 도심하천으로부터 유입이 되고, 농경지와 축산 농가 및 산지로부터 발생된 오염물질이 하천에 큰 영향을 주어 자정능력 범위를 초과하게 되고 수질의 오염이 더욱 증가하게 된다(Rim, 1999). 하천 수질은 유입되는 지류로부터 오염부하, 수문학적 특성, 퇴적물 및 수체내 물질대사, 계절적 요인 등이 상호작용하여 결정되기 때문에 수질 변화의 특성을 파악하기에는 한계가 있으며, 시･공간적으로 변화하는 특성을 지니게 된다(Park et al., 2011). 따라서 효율적인 수질 관리를 위하여 주된 오염인자를 합리적으로 산출하고 신뢰성 있는 수질자료의 해석을 지원하기 위해서는 다변량 분석법을 이용하여 수질관리를 수행할 필요가 있다. 이러한 통계분석 기법은 하천 및 호소의 수질 특성을 평가하고 수질 정책을 위한 방법으로 많이 활용하였으며(Lee and Kwon, 2009; Oh et al., 2004; Seymour et al., 2004; Simeonov et al., 2003; Yu et al., 2003), 농촌유역과 지하수의 수질특성 평가를 위한 방법으로 많은 연구들이 수행 되었다(Kim, Choi et al., 2007; Kim, et al., 2008). 또한 영산강을 비롯하여 한강, 금강, 낙동강 등에서 다변량 해석기법을 이용한 연구가 활발히 진행되어 왔다(Jung et al., 2013; Jung, Lee et al., 2012; Kim, Kim et al., 2007; Kim et al., 1998; Kim, Lee et al., 2007; Kim and Lee, 2011). 그러나 아직까지 섬진강 수계를 대상으로 수행된 연구는 부족한 실정이며, 본 연구에서는 다변량 분석기법의 상관분석(Correlation Analysis)과 주성분분석(Principle Component Analysi s) 및 요인분석(Factor Analysis)을 통해 수질인자간의 상관성과 수질 특성을 파악하였으며, 향후 섬진강 수계의 수질변동에 대해 능동적으로 대처할 수 있는 보다 과학적이고 합리적인 수질관리정책 수립을 위한 기초 자료를 제공하는데 목적이 있다.

섬진강은 남해안 중서부에 위치하며 동경 126° 52' ~ 127° 53', 북위 34° 40' ~ 35° 50' 사이에 위치하며, 유역면적이 4,960 km²로 우리나라에서 4번째로 큰 유역을 가지고 있다. 전반적으로 산지 지형이 우세하며 특히 이 유역의 동측에 위치한 지리산은 높은 고도의 거대한 산체를 이루면서 곧바로 본류와 접하고 있는 것이 특징이다. 유역의 본류 및 보성강을 비롯한 지류들이 비교적 좁게 계곡을 따라 흐르며, 유역의 본류가 지나는 남원과 구례 등 일부 지역을 제외하고는 평야지대의 분포가 매우 협소한 지역이다. 유역 내의 행정구역으로는 전라북도의 정읍시, 남원시, 진안군, 장수군, 임실군, 순창군, 전라남도의 순천시, 광양시, 담양군, 곡성군, 구례군, 화순군, 보성군, 장흥군, 경상남도의 하동군을 포함하여 3개의 도와 4개의 시 그리고 11개 군으로 이루어져 있고, 도별 유역 편입율은 전라남도 47%, 전라북도 44%, 경상남도 9%이다. 주요지류는 섬진강 수계에서 제일 큰 1지류로 국가하천인 보성강이 유입되며 전라북도 남원시 송동면 세전리 지점에서 섬진강 본류의 좌안으로 유입하는 국가하천인 요천이 흐른다(NIER, 2012). 본 연구의 모니터링 지점은 총량관리 단위유역 말단부로서 섬진강수계 내의 본류 지점과 본류로 유입되는 지류 지점이며, 본류 구간의 섬본A(SA), 섬본B(SB), 섬본C(SC), 섬본D(SD), 섬본E(SE)지점과 추령A(CA), 오수A(OA), 요천A(YA), 요천B(YB), 보성A(BA), 동복A(DA), 보성B(BB)의 지류 구간으로서 총 12개 지점이다(Fig. 1).

본 연구의 12개 모니터링 지점에 대한 수질조사는 2010년 1월부터 2012년 12월까지 3년간 유역 출구 지점에서 평균 8일 간격으로 실측된 수질 자료를 이용하였다. 현장 측정항목은 Water Temperature (WT), pH, Electric Conductivity (EC), Dissolved Oxygen (DO), Biochemical Oxygen Demand (BOD), Chemical Oxygen Demand (COD), Suspended Solid (SS), Total Nitrogen (T-N), Total Phosphorus (T-P), Total Organic Carbon (TOC)의 수질 시료 채취와 분석은 수질오염공정시험기준에 준하여 분석하였다. 수질 특성 분석을 위한 각 측정지점의 통계학적 특성치는 Table 1과 같으며, 수질항목별 결과를 최소값, 최대값, 1사분위수, 3사분위수, 중앙값을 상자그림(Box-Plot)으로 나타내었다(Fig. 2).

본 연구에서의 모든 통계분석은 SPSS(ver. 18.0) 프로그램을 이용하였으며, 상관분석은 변수간의 관련성을 분석하기 위하여 두 변수간의 선형적인 상관관계를 파악하는 것으로 한 변수가 다른 변수와의 관련성이 있는지 여부와 관련성이 있다면 어느 정도의 관련성이 있는지를 알고자 할때 유용한 기법이다. 여기서, 보편적으로 자주 이용되는 척도로는 Pearson 상관계수(Correlation Coefficient, r)가 많이 이용되며, 상관계수는 −1에서 1까지의 값을 취하는데 절대값이 1에 가까울수록 상관성이 강하다는 것을 의미한다.

다변량 분석법 중 주성분 및 요인분석을 실행하여 각 유역의 수질변동 특성을 파악하고자 한다. 주성분분석은 다변량 분석의 여러 가지 기법 중에서도 가장 기본이 되는 분석 방법으로서 변수 간에 상관성을 갖는 여러 종류의 특성치들의 정보를 큰 손실 없이 요약 표현이 가능하며, 상호 유관한 특성치들의 변화 양상으로부터 이들 속에 내재하고 있는 상호 독립적인 고유의 패턴을 도출할 수 있다. 또한 다변량으로 구성되는 데이터로부터 소수의 특징적인 변량을 합성하여 주성분을 산출하고, 차원이 축소되어 데이터 분석이 용이해진 주성분을 사용하여 전체 다변량 데이터의 경향을 분석하게 된다(Park and Rhee, 2012). 즉, 주성분을 구하고 해석하여 원 자료를 축약하는 방법을 주성분분석이라 하고 일단 자료를 간단히 하여 추후에 추가적인 다른 통계적 절차를 적용하기 위한 사전분석으로 볼 수 있다 (Kim, Lee et al., 2007).

요인분석 기법은 특정 변수들 간의 상호 연관성을 보다 간결하게 분석해서 공통적으로 내재된 인자를 추출하여 전체 자료를 대변할 수 있는 변수의 수를 줄이는 기법으로 상관관계가 높은 변수들이 가지고 있는 공통성을 중심으로 하나의 동질적인 요인으로 묶어주는 분석 방법이다(Jung et al., 2013). 한편, 추출된 고유치(Eigen Value)는 요인이 설명할 수 있는 변수들의 분산 크기를 나타내는 것으로 고유치가 1보다 크다는 것은 하나의 요인이 변수 1개 이상의 분석을 설명해 준다는 것을 의미하며, 고유치가 1보다 작다는 것은 요인이 변수 1개의 분산도를 설명해 줄 수 없다는 것을 의미하므로 주요인을 추출하기 위해서는 통상 고유치가 1 이상인 것을 기준으로 결정한다(Park et al., 2013). 그러나 변수들이 어떤 요인들에 대하여 요인 부하량을 나타낼 때 비슷한 부하량을 나타낼 경우 어느 요인에 해당하는지 분류하기가 쉽지 않기 때문에 비슷한 요인 부하량이 어느 한 요인에 높게 나타나도록 하기 위해서는 요인축을 회전시키는 방식을 이용한다. 이러한 회전방식은 사각회전(Oblique Rotation)과 직각회전(Orthogonal Rotation) 방식으로 구분되고 직각회전방식에는 Varimax, Quartimax, Equimax 등이 있으며, 이중에서 보편적으로 많이 이용되는 Varimax 직각회전방식을 사용하였다.

섬진강 수계의 수질특성을 알아보기 위하여 각 지점별 특성치와 Box plot을 이용한 조사결과를 Fig. 2와 Table 1에 나타내었다. 섬진강 수계의 전체 지점에 대한 수질 농도를 살펴보면 BOD₅의 농도 범위는 0.4~6.8 mg/L, 평균 1.5 mg/L이며, COD의 경우 1.4~14.3 mg/L, 평균 3.8 mg/L로 나타나 BOD₅와 COD는 하천 생활환경기준 Іb등급(좋음)에 해당한다. TOC의 농도 범위는 0.2~7.5 mg/L, 평균 2.3 mg/L 로 나타났으며, SS의 경우 0.2~294 mg/L, 평균 6.9 mg/L이며, T-N은 0.4~4.6 mg/L, 평균 1.9 mg/L로 나타났다. T-P의 농도 범위는 0.005~0.386 mg/L, 평균 0.040 mg/L를 나타내어 하천 생활환경기준 Ib등급(좋음)에 해당한다. 섬진강 수계를 전체 지점의 수질 농도로 평가한다면 양호한 수준으로 볼 수 있다. 그러나 섬진강 본류 지점의 SC, SD와 지류 지점의 OA, YB에서 COD는 하천 생활환경기준 II등급(약간 좋음)에 해당하며, T-P의 경우 II등급(약간 좋음)과 III등급(보통)으로 조사되었다. 전반적으로 OA와 YB지점에서 COD, TOC, SS, T-P의 수질 항목에 대한 평균 농도가 타지점과 비교해 볼 때 높은 것으로 조사 되었으며, 이러한 이유는 지점별 특성으로서 OA지점의 경우 농경지 및 축산 농가 등이 유역 인근에 위치하는 특성을 가지고 있으며, YB지점은 도시생활하수가 유입되는 도심하천인 동시에 요천의 말단부에는 축산 농가 등이 하천 인근에 위치하는 특성을 가지고 있기 때문인 것으로 판단된다. 이에 섬진강 수계의 효율적인 수질관리를 위해서는 섬진강의 본류로 유입되는 지방하천에 해당하는 오수천과 국가 하천인 요천에 대한 집중관리가 선행되어야 할 것으로 사료된다.

섬진강 수계의 상관분석을 위해 각 수질항목별로 약 1550개의 수질 자료를 이용하여 분석한 결과를 Table 2에 나타내었다. 유기물 간접지표인 BOD₅와 COD, TOC의 상관계수를 살펴보면 0.715, 0.719로서 높은 양의 상관성(p<0.01)을 보이며, 특히 COD와 TOC의 경우 상관계수가 0.854로서 가장 높은 상관성을 보였다. 또한, T-P와 BOD₅, COD, TOC, SS에 대한 상관계수는 각각 0.482, 0.599, 0.608, 0.493으로 보통의 상관성을 보였다. 이는 본 연구대상 지점의 경우 특정 수질항목(BOD₅, COD, TOC, T-P, SS) 중 개별 수질이 높으면 나머지 항목도 높아질 것으로 판단되며, Jung, Lim et al. (2012)이 영산강 수계의 소하천을 대상으로 수행한 수질항목별 상관분석의 연구결과와 유사하다. 그러나 T-N과 BOD₅, COD, TOC, T-P, SS의 경우 상관계수가 −0.023, −0.022, 0.010, 0.318, 0.112로서 상관성이 매우 낮은 것으로 보이며, BOD₅와 COD의 경우 T-N의 수질 농도가 높아지면 반대로 다소 낮아지는 음의 상관성은 보였다. 이러한 이유는 BOD₅와 COD의 경우 평상시 보통의 수질을 유지하다가 강우 시에 농도가 높아지게 되며, 반면에 T-N의 경우에는 기저 유출에 의한 영향을 받아 비강우시에도 높은 농도를 나타내어 상대적으로 강우 시에는 영향을 적게 받는 것으로 사료 된다(Kim et al., 2004). 수온과 DO의 상관계수가 −0.837로서 음의 상관성(p<0.01)을 보였는데, 이는 온도가 높을수록 산소의 용해도가 낮아지는 전형적인 경향이며 계절에 의한 영향으로 볼 수 있다.

3.3. 주성분 및 요인분석결과

연구지점에서 3년간 수집된 10개의 수질항목을 대상으로 해당 자료의 전체 수질(2010~2012)과 계절별 수질(봄 3~5월, 여름 6~8월, 가을 9~11월, 겨울 1~2월, 12월)로 구분하여 주성분 분석과 요인분석을 수행하였다. 주성분 수를 결정하기 위해서 요인을 설명하는 분산의 크기를 나타내는 고유치의 값이 1이상인 성분 축만을 고려하였다. 이는 고유치가 1이상일 경우에 하나의 요인이 변수 1개 이상의 분산을 설명하며, 고유치가 1이하인 경우에는 1개의 요인이 변수 1개의 분산을 설명할 수 없으므로 고유치가 1이상이 되는 요인을 추출하는 방법을 사용하였다(Kim, Lee et al., 2007). 전체 수질의 주성분 분석에 대한 결과를 Table 3에 나타내었으며, 요인분석 결과는 요인과 변량과의 상관계수에 따른 요인 구조를 명확하게 하기 위하여 직각 회전 방식인 Varimax 방식을 적용하여 Table 4에 나타내었다. 그 결과 3개의 주성분이 추출되었고 제 1요인의 고유치 3.788로서 37.884% 기여하고 있으며, 나머지 제 2요인과 제 3요인의 고유치와 기여율은 각각 2.019, 1.091과 20.193%, 10.913%로서 제 1요인부터 제 3요인까지 전체 분산의 68.990%를 설명해 주고 있다. 제 1요인은 BOD₅, COD, TOC, SS, T-P로 분류되었고 제 2요인은 T-N과 DO이며, 제 3요인은 pH로 분류되었다.

계절별 수질의 주성분 분석과 요인분석에 대한 결과를 Fig. 4와 Table 3에 나타내었다. 봄과 여름은 4개의 주성분이 추출되었고 고유치는 1이상에서 봄의 누적 기여율은 81.515%, 여름은 73.550%를 설명해 주고 있다. 가을과 겨울의 경우에는 3개의 주성분이 추출되었으며, 고유치 1이상에서 각각 65.072%와 72.721%로서 수질 변동을 설명해 주는 것으로 나타났다. 계절별 수질의 제 1요인은 가을을 제외한 모든 계절에서 전반적으로 BOD₅, COD, TOC, SS, T-P로서 전체 수질의 결과와 같았으며, 전체 결과를 종합해 보면 섬진강 수계의 대상 지점은 유기물과 영양염류, 부유물질의 영향을 동시에 받는 것으로 설명할 수 있다. Lee et al. (2008)에 의하면 T-P의 경우 토양 입자와 결합되어 평소에는 유출이 일어나지 않지만 강우 시 토양입자의 이탈과 함께 수계로 다량 유출된다고 보고 하였다. 즉 T-P와 SS를 비롯하여 BOD₅, COD, TOC가 상관성이 높고 제 1요인으로 추출된 것으로 볼 때 강우 시에 수질 농도가 함께 상승하는 것으로 볼 수 있으며, 점오염원에 의한 영향 보다는 비점오염원의 영향이 더 큰 것으로도 해석할 수 있다. 이러한 결과는 Park et al. (2013)이 섬진강 수계의 주요 지천(추령A, 오수A)을 대상으로 유량-부하량 관계식을 이용하여 수행한 연구 결과와 유사하였다. 그러나 T-N의 경우 제 1요인에서 제외 되었는데 이는 앞서 설명하였듯이 해당수계에 영향을 주는 것은 분명하나 T-N의 경우에는 기저 유출에 의한 영향을 받아 비강우시에도 높은 농도를 나타내어 상대적으로 강우 시에는 영향을 적게 받으므로 제 1요인과 상관성이 낮아지는 경향으로 볼 수 있다.

한편, 요인분석의 실행결과에 대한 타당성을 알아보기 위하여 Kaiser-Meyer-Olki n(KMO)과 Bartlett's test를 실행하였으며, KMO test의 경우에는 측정할 샘플에 대한 변수의 크기가 적당한 가를 나타내는 척도로서 1에 근접한 값일수록 요인분석의 타당성이 높은 것으로 볼 수 있으며, 0.5이하일 경우에는 분석이 타당성이 없는 것으로 판단할 수 있다. Bartlett's test는 변수 간의 상관행렬이 단위행렬(대각선이 1이고 나머지는 0인 행렬)인지 아닌지를 검정하는 것을 말하며, 각 변수가 서로 상관성이 존재함을 나타낼 때 유의적인 관계라 하고 이때 0에 가까울수록 유의성이 높다고 할 수 있다(Bernad et al., 2004; Kim, Lee et al., 2007).

본 연구의 KMO test에 대한 결과에서 전체 수질을 대상으로 수행한 분석의 검정 값은 0.725이며, 계절별 수질의 분석에 대한 봄, 여름, 가을, 겨울의 검정 값은 0.756, 0.666, 0.684, 0.789로 나타났고, Bartlett's test는 모두 0.000(유의수준<0.05)이며 변수 간의 상관성이 있음을 나타내어 상관행령이 단위행렬임을 기각하여 요인분석이 가능함을 설명해주고 있다.

본 연구에서는 섬진강 수계를 대상으로 다변량 분석기법을 적용하여 수질 항목간 상관분석과 주성분 및 요인분석을 수행하여 수질변동 특성을 알아보고자 하였으며, 그 결과는 다음과 같다. 섬진강 수계는 본류로 유입되는 지류에 의한 영향이 큰 것으로 볼 수 있으며, 그중에서 OA와 YB지점의 COD, TOC, SS, T-P에 대한 평균 농도가 타 지점과 비교해 볼 때 높은 것으로 조사되어 해당 지점은 집중 관리가 필요할 것으로 판단된다. 수질인자간의 상관분석 결과에서 T-N을 제외한 수질항목(BOD₅, COD, TOC, SS, T-P)들은 전반적으로 높은 상관성을 나타내었으며, 연구대상 지점은 상관성이 높은 수질항목 중 특정 수질항목이 높으면 나머지 항목도 높아지는 것으로 해석할 수 있다. 한편 주성분 및 요인분석의 결과에서 전체 수질을 대상으로 3개의 주성분이 추출되었으며, 전체 분산의 68.990% 기여하였고, 계절별 수질의 경우에는 봄과 여름은 각각 4개의 주성분이 추출되었으며, 누적 기여율은 81.515%, 73.550%이다. 가을과 겨울의 경우 각각 3개의 주성분이 추출되었으며 65.072%와 72.721%의 누적 기여율을 나타냈다. 요인분석의 경우 계절적 특이성은 보이지 않았으며, 공통적으로 제 1요인은 BOD₅, COD, TOC, SS, T-P로 분류되었다. 전반적인 결과를 종합해 보면 섬진강 수계의 대상 지점은 유기물과 영양염류, 부유물질의 영향을 동시에 받는 것으로 사료된다.