1.서 론

현재 부산연안은 어항·상업지역·친수공간 등의 다양한 기 능을 제공하고 있으나, 하천 및 하수처리장을 통해 생활하 수 및 산업폐수가 만 내로 유입되고 있어 해역의 오염이 우 려되는 해역이다(Baek, 2012). 그 중 수영만은 대한민국의 대 표적 관광지인 해수욕장이 위치하고 있어 청정해역을 필요 로 한다. 하지만, 2011년 기준 수영만 연안의 COD 농도는 목 표수질인 1 mg/L를 초과하고 있다(MOF, 2012a). 또한, 광안리 및 해운대 해수욕장 인근에서 생산된 수산물은 식중독과 비 브리오 패혈증의 위험이 있다(Kim et al., 2007b). 따라서, 국 민건강과 지역민의 생업 유지를 위해서라도 부산연안 해역 의 수질악화는 저지되어야 한다.

연안해역의 수질 개선을 위해서는 유입되는 부하량을 관 리하는 것이 매우 중요하다. 부하량 감소에 따른 수질의 응 답특성을 평가하는데 생태계 모델은 유용한 수단이 된다 (Hong et al., 2007; Susilowati et al., 2004). 최근 생태계 모델을 이용하여, 부하량 삭감에 따른 수질의 응답특성에 관한 많 은 관심이 모이고 있다(Hong et al., 2007; Gaddis and Voinov, 2010; Zhao et al., 2012). 부산연안 중 수영만을 대상으로 한 수질모델 사례는 Park et al.(1999), Lee and Kwon(1995), Cho et al.(1995) 등이 있으나, 연구범위가 수영만으로 한정적이다. 또한, 부산연안에 관련된 모델연구와 해역 수질관리 정책제 안을 위한 연구는 드물다 . 따라서, 본 연구는 생태계 모델을 이용하여 부산연안으로 유입되는 부하량을 삭감, 이에 따른 해역의 수질개선을 예측하는 것이 목적이다.

2.재료 및 방법

본 연구에서는 생태계 모델을 이용하여 2011년 부산연안 수질을 재현하고, 이를 토대로 2020년 수질 예측 및 부하량 삭감에 따른 2020년 수질개선을 예측하고자 하였다.

2.1.생태계 모델 구성

본 연구에서 사용한 생태계 모델은 EM3(Ecological Model for Marine water quality Management)이다. EM3는 2008년 마산 만 특별관리해역의 제 1차 연안오염총량관리제 도입 및 시 행을 위해 구성되었다. EM3는 마산만·부산연안 오염총량관 리에 사용되고 있다. 또한, Park et al.(2010), Jung et al.(2016)에 의한 연구사례도 있다. 본 연구에서는 EM3의 구성요소중 6 개의 수층 유기태 요소, 4개의 수층 무기태 요소, 5개의 저질 요소 및 2개의 수질 요소를 고려하였다(Fig. 1, Table 1). 생태 계 모델의 격자는 가덕도 동안으로부터 기장 연안사이의 해 역을 대상으로 하였으며, 수심은 해도를 이용하여 산정하였 다. 격자 구성은 동서 및 남북 방향으로 각각 500 m, 수직방 향으로는 5개의 Sigma격자로 구성하였다. 생태계 모델의 계 산간격은 400 sec로 하였다(Table 2) 유동장은 MOF(2012a)와 MOF(2014)에서 수행되었던 해수유동 모델결과를 활용하였 다.

2.2.생태계 모델 입력자료

생태계 모델의 입력자료는 MOF(2012a) 자료를 사용하였 다. 입력값은 2011년 8월부터 2012년 5월까지 분기별로 4회 관측된 자료를 바탕으로 산정하였으며, 경계값은 Neumann boundary condition 방식을 사용하였다(Table 3, Table 4).

모델에 입력된 부하량 자료는 2011년 관측값, 2020년 예측 값, 2020년 예측값에서 부하량을 삭감한 값으로 나타내었다. 2011년 부하량 값은 모델의 검증에 필요한 자료로서 관측된 자료이며, 2020년 부하량 값들은 MOF(2012a)에서 예측한 부 하량 값이다. 2020년 예측부하량은 2007년~2011년까지의 자 료를 이용하여 행정구역별, 배출원별 오염원을 전망하였다 고 되어 있으며, 2020년 기준 삭감 부하량은 수질개선계획에 따른 하수처리시설, 분류식 하수관거, 비점오염 저감시설 등 의 확충, 개선, 설치를 통해 예측하였다고 되어있다. 2020년 예측 부하량은 2010년에 비해 약 COD는 +3,488.8 kg/day, T-N 은 +454.9 kg/day, T-P는 -79.1 kg/day로 나타났다. 부하량을 삭 감한 2020년 부하량은 2020년 예측 부하량에 비해 약 COD는 -4,099.4 kg/day, T-N은 -5,165.4 kg/day, T-P는 -359.4 kg/day로 나 타났다. 영역별 부하량의 양을 비교하면 COD는 낙동강, 수 영만 순으로 나타났다. T-N은 수영만, 낙동강 순으로 나타났 다. T-P는 수영만, 낙동강 순으로 나타났다. 유기오염물질은 수영만보다 낙동강이 높게 나타났고, 영양염은 수영만이 낙 동강보다 높게 나타난 것을 볼 수 있다. 수영만과 낙동강을 제외한 동부산과 부산항은 상대적으로 낮은 부하량을 나타 내었다. 2020년 삭감부하량을 보면 COD, T-N, T-P 모두 수영 만이 가장 크게 나타났으며 유입되는 부하량대 삭감비도 가 장 크게 나타났다. 삭감량은 수영만, 낙동강, 부산항, 동부산 순으로 나타났다(Table 5). 본 모델에 사용한 생물학적 계수 (parameter)는 국내 모델연구 및 문헌에서 사용된 자료를 검 토한 후 입력하였다(MOF, 2012a; MOF, 2012b; MOF, 2012c; MOF, 2014; Park et al., 2010).

2.3.생태계 모델을 이용한 수질개선예측

2011년 검증자료로 검증된 생태계 모델을 바탕으로 2020 년 부산연안의 수질개선을 예측하였다. 그 후 2020년 예측부 하량을 반영하여 부산연안의 수질을 예측하고, 부하량을 삭 감했을 때 부산연안의 수질개선을 알아보는 것이 목적이다. 대상수질 항목으로는 유기오염물질인 COD와 영양염인 T-N 과 T-P를 대상으로 하였다.

3.결과 및 고찰

3.1.생태계 모델 검증

모델 검증정점은 Fig. 2에 나타낸 20개 정점을 대상으로 하였다. 검증은 관측치 표층값과 계산치의 상관성 분석 결 과로, 화학적산소요구량(COD)는 0.6440, 용존무기인(DIP)는 0.8620, 용존무기질소(DIN)은 0.7978으로 나타나 비교적 재현 이 잘된 것으로 판단된다(Fig. 3 ~ Fig. 5).Fig. 4

재현된 COD 공간분포를 보면 해운대 해수욕장 왼쪽부터 이기대를 기점으로, 만 내측 방향으로 2.0 mg/L를 초과하는 값을 나타내고 있다. 수영만 외측에서는 1.5 mg/L 이하를 나 타내고 있다. 낙동강하구에서는 2.0 mg/L 부근의 값을 나타 내고 있으며, 동쪽으로 갈수록 농도가 낮아지며, 기장 해역 에서는 1.0 mg/L 이하를 나타내고 있다(Fig. 6). T-N분포는 다 음과 같다. 수영만 하구에서는 1.0 mg/L 이상을 나타내고 상 대적으로 오염된 것을 알 수 있다. 광안리 해수욕장까지 약 0.6 mg/L 이상을 나타내고 있으며, 만 외해로 가면서 농도가 낮아지는 경향을 나타낸다. 낙동강하구에서는 약 0.7 mg/L 이상의 값을 나타내며, 동쪽으로 갈수록 농도가 낮아지는 경향을 나타낸다. 수영만 외측부터 기장해역으로 가면서 0.3 mg/L 이하의 값을 나타내고 있다(Fig. 7). T-P의 분포는 다음 과 같다. 수영만 하구에서는 0.05 mg/L 이상의 값을 나타내고 있으며, 광안리 해수욕장에서는 0.04 mg/L 이상의 값을 나타 내고 있다. 낙동강하구에서는 0.03 mg/L 이상의 값을 나타내 며, 수영만을 제외하면 우측으로 갈수록 낮아지는 경향을 나타낸다(Fig. 8).

3.2.수질개선 예측 결과

검증된 수질결과를 바탕으로 본 모델에 2020년 예측 부하 량을 적용하고, 부하량 삭감시의 부산연안의 수질개선을 예 측하였다.

2020년 예측 부하량을 입력한 모델결과를 보면 2011년과 부하량 차이가 크지 않기 때문에, 해역의 수질 또한 2011년 과 유사한 경향을 나타내었다(Fig. 9 ~ Fig. 11).Fig. 10

부하량을 삭감한 COD수질 분포결과는 다음과 같다. 광안 리 하부와 이기대 사이에서는 여전히 2.0 mg/L 이상의 농도 를 나타내지만, 삭감 전과 비교하면 그 범위가 상당히 많이 줄어 약 0.5 mg/L 개선된 것을 알 수 있다. 수영만 외측에서 는 1.5 mg/L 이하의 값을 나타내고 있다. 낙동강 하구에서는 2.0 mg/L 이하의 값을 나타내어 수질이 개선된 것을 알 수 있다. 낙동강 하구 우측부터 이기대 아래쪽 까지는 상대적 으로 미약하지만 수질이 개선된 것을 알 수 있다. 기장해역 부근에서는 크게 개선된 결과를 나타내지는 않았다. T-N분 포는 다음과 같다. 수영만 하구에서는 0.7 mg/L 부근의 값을 나타내고 있어, 삭감 전에 비해 수질이 약 0.2 mg/L 개선된 것으로 나타났다. 낙동강 하구에서는 0.7 mg/L 이상이었던 수질이 0.5 mg/L 이하로 개선된 것을 알 수 있다. 낙동강 하 구 동쪽으로는 크게 개선되지는 않은 것을 알 수 있다. T-P 의 분포는 다음과 같다. 수영강 하구에서는 0.05 mg/L 이상의 값에서 0.04 mg/L 부근의 값으로 약 0.01 mg/L 이상 개선된 것으로 나타났다. 낙동강 하구에서는 0.03 mg/L 이상의 값에 서 0.02 mg/L 부근의 값을 나타내어 수질이 개선된 것으로 나타났다. 연안 동쪽으로 갈수록 약간의 개선은 있지만 그 값이 크지 않은 것으로 나타났다.

2020년을 대상으로 부하삭감량을 적용했을 때 부산연안 의 수질은 광안리 해수욕장이 포함된 수영만에서 가장 큰 개 선효과를 나타내었고, 낙동강 하구에서도 수질 개선이 나타 났다. 수영만을 제외하고 낙동강 하구 동측부터 기장해역까 지는 뚜렷한 수질개선이 나타나지 않았다. COD 수질은 전반 적으로 2.0 mg/L 이하, T-N은 0.5 mg/L 이하, T-P는 0.03 mg/L 이하의 값을 나타내었다(Fig. 12 ~ Fig. 14).Fig. 13

본 연구에서 COD 유입 부하량은 20,119 kg/day이며, 부하 삭감량은 4,099 kg/day이다(Table 6). 본 연구와 유사하게 부하 량을 삭감하고 COD 수질변화를 예측한 문헌과 비교를 해보 면 다음과 같다. 시화호를 대상으로 한 Kim et al.(2007a)의 연 구에서는 해역 COD 농도를 Ⅲ 등급 기준인 4mg/L 이하로 하 기 위해서는 유입되는 COD 부하량 20,000 kg/day 중 15,000 kg/day를 삭감해야 되는 것으로 나타냈다. 또한, 유입되는 부 하량을 95%를 삭감해도 COD 농도가 2 mg/L 이하로는 어려 울 것으로 판단했다. 진해만을 대상으로 한 Hong et al.(2007) 의 연구에서는 전반적으로 해역 COD 농도가 2 mg/L 이하를 만족시키기 위해서는 유입되는 COD 부하량 18,774 kg/day 중 5,632 kg/day를 삭감해야 되는 것으로 나타났다. 마산만의 모 든 부하량이 덕동하수처리장에서 처리된다는 가정을 한 Oh et al.(2005)의 연구결과에서는, 마산만 전체해역 COD 농도를 2.5 mg/L 이하로 만족시키기 위해서는, 유입되는 COD 부하 량 33,600 kg/day 중 30,200 kg/day를 삭감해야 될 것으로 판단 하였다. 이와 비교했을 때, 부산연안의 COD 부하 삭감량은 다른 해역에 비해 가장 적게 나타났다. 또한 COD 부하 삭감 비율은 부산연안 16.92%, 시화호 75.00 %, 진해만 29.99 %, 마 산만 89.88%로 부산연안이 가장 작게 나타났다. COD 수질 또한 수영만 연안의 특정부분을 제외하면, COD 농도가 2 mg/L 이하로 부산연안이 가장 낮은 농도를 나타낸 것으로 판단된다(Table 6). 이는 부산연안이 남해에 위치한 다른 해 역에 비해 개방형경계를 가지고 있어서, 물질교환이 빠르기 때문으로 판단된다. 따라서, 부산연안에서는 상대적으로 부 하량을 적게 삭감하여도 수질개선효과가 나타났다.

본 연구에 적용한 부하삭감량을 부산연안에 적용한다면 부산연안 해역의 수질이 개선될 것으로 판단된다. 수영만연 안과 낙동강하구 부근에서 뚜렷한 수질개선을 나타냈으며, 부산연안 대부분 해역에서 양호한 수질을 달성 할 수 있을 것으로 판단된다.

4.결 론

생태계 모델링을 이용하여 부산연안의 2020년 수질개선 을 예측하였다. 방법으로는 2011년의 수질재현을 통해 모델 을 검증하고, 2020년 해역수질 예측 및 유입부하량 삭감에 따른 2020년 수질개선을 예측하였다. 부하량 삭감 후의 예측 결과 상대적으로 수영만과 낙동강하구에서 COD 수질개선이 뚜렷하게 나타났다. 수영만을 제외한 만 중부에서 만 동쪽 해역까지는 수질개선이 거의 나타나지 않았다. T-N과 T-P는 COD와 유사한 경향을 나타냈다.

COD 삭감 결과를 유사한 다른 문헌과 비교해보면 다음과 같다. 부산연안은 해역으로 유입되는 COD부하의 삭감량과, 총유입부하량대비 삭감비율이 16.92%로 상대적으로 다른 해역에 비해 작았지만, 해역의 전반적인 수질이 다른 해역 에 비해 수질이 양호하게 나타났다. 이는 부산연안이 다른 만에 비해 개방형 경계이며 물질교환이 빠르기 때문으로 판 단된다.

본 연구에 적용한 부하삭감량을 부산연안에 적용할 경우, 부산연안 해역에서 수질개선효과를 나타낼 것으로 판단된 다. 추가적으로 오염이 가장 우려되는 하계에 수영만에 위 치한 해수욕장으로 수백만 이상의 피서객이 방문하고 있다. 본 연구에서 수질농도가 가장 높게 나타난 수영만 해역에 대해 추후 세부적인 연구가 필요할 것으로 판단된다.