2.1. 연구대상지역

본 연구의 대상지역인 영산강은 한반도 남서부 전라 남· 북도에 위치하며, 영산강, 황룡강, 지석천, 고막원천, 함평 천 등 국가하천 5개와 지방하천 163개로 구성되어 있다. 하천에 인접하여 농경지가 발달해 있으며, 강의 중·상류부 에는 광주광역시와 나주시 등의 도심지 및 인구밀집지역이 위치해 있다(^{MLTMA, 2009}). 영산강의 유로연장은 약 129.5 km, 유역면적은 약 3,469 km² 이며, 유역 내에는 국내 최대 규모의 농업용 저수지 장성호, 나주호, 담양호, 광주호가 위 치하고 있다. 4개소 저수지의 유효저수용량의 합이 약 2억 5천만 m³ 정도에 달하며, 이들 저수지의 유역면적이 영산 강 전체 유역면적의 약 9%를 차지한다. 본 연구에서는 영 산강의 주요 제1지류 10개를 중심으로 소유역을 구성하였 으며, 소유역은 상류에서부터 증암천, 풍영정천, 광주천, 황 룡강, 지석천, 만봉천, 고막원천, 함평천, 삼포천, 영암천을 기준으로 구분하였다(Fig. 1).

Fig. 1. The tributaries of Yeongsan river.

2.2. HSPF 모델개요

HSPF는 유역 내에서 장기간에 걸쳐 발생하는 복잡한 수 리·수질에 관한 다양한 현상들을 모의하기 위해 미국 환경 청(U.S. Envioronmental Protection Agency, U.S.EPA)에서 개발한 유역모델로써 투수층(PERLND), 불투수층(IMPLND), 하천(RCHRES)의 3가지 주요 응용모듈로 구성된다. PERLND 모듈은 유역 내 투수 지역에서 발생하는 주요 수문 현상과 수질을 모의한다. 투수지역에서 유역의 수문 현상에 영향을 주는 요소에는 강우, 식물 등에 의한 차단, 지표 저류, 증 발산, 지표수 유출, 중간유출, 지하수 유출 등이 있으며, 이 러한 요소들은 토지피복, 토양, 강우 등에 의해 영향을 받 는다. 수문 현상을 고려한 유역의 물수지 기본 방정식은 아래 식과 같다.

여기서, R은 유출량, P는 강우량, ET는 증발산량, IG는 심층 지하수량, ΔS 저류량 변화를 나타낸다. IMPLND 모 듈은 불투수 지역에서 발생하는 주요 수문 현상과 수질을 모의하며, 투수지역과 달리 침투와 지표하 수문과정을 모의 하지 않아 보다 간단한 구조를 가진다. RCHRES 모듈은 하도구간이나 호소에서의 수문, 유사이동, 수질 기작을 모 의한다. 수체 내 물수지 기본 방정식은 아래 식과 같다.

여기서, VOL은 최종유량, VOLS는 최초유량, IVOL은 유 입량, PRSUPY은 강우량, VOLEV은 증발량, ROVOL은 유 출량을 나타낸다.

유역에서의 토사 침식 현상은 강우, 지표수 유출, 바람 등에 의해 발생하며, 3가지 주요 응용모듈별로 다른 토사 이동 방정식이 적용된다. PERLND 모듈에서 토사 유출량 은 강우에 의해 토양으로부터 분리된 토사 유출량과 유수 에 의한 토양 침식량의 합으로 산정된다. IMPLND 모듈에 서는 불투수층에서 축적되거나 제거되는 토사와 토사이송 능력간의 관계를 통해 유출량이 산정된다. RCHRES 모듈 에서는 수체 내 실트, 점토, 모래의 침식, 이동, 퇴적을 모 의한다. 점착성 퇴적물인 실트와 점토의 침식과 퇴적 작용은 전단응력(shear stress)와 임계전단응력(critical shear stress) 등을 이용하여 모의되며, 비 점착성 퇴적물인 모래의 침식 과 퇴적 여부는 유속에 의해서 정의되는 모래 이송 능력 (sand transport capacity)에 의해서 모의된다.

수체 내 수질 기작은 RCHRES 모듈의 RQUAL(Constituents Involved in Biochemical Transformations)항목을 통 해 모의된다. RQUAL 항목에서 모의하는 수질 성분은 DO, BOD, PO₄, NO₂, NO₃, NH₄, 식물성과 동물성 플랑크톤, 저 서 조류(Benthic Algae) 등을 포함하며, 각 모듈의 주요 기 능과 구조에 대한 더욱 구체적인 내용은 ^{Bicknell et al. (2001)}에 상세히 기술되어 있다. 한편, HSPF 모델을 구동 하기 위해 수치고도모델(Digital Elevation Model, DEM), 토지이용도, 토성 등의 공간자료와 기상자료 등의 기본 입 력 자료가 필요하다. 공간자료는 BASINS의 GIS 기능을 이 용하여 구축하며, 강수량, 온도, 풍속, 증발량, 일사량, 이슬 점 온도, 운량 등을 포함하는 기상자료는 USEPA에서 개발 된 자료관리 프로그램 WDM(Watershed Data Management) 을 이용하여 구축한다.

2.3. 입력자료 및 모델구축

모델의 주요 공간 입력파일은 DEM, 하천도, 유역도, 토 지이용도, 토성, 기상 및 유량관측소 위치도 등을 포함한다. DEM은 환경부에서 제공하는 30 m 해상도 자료를 이용하 였으며(Fig. 2(a)), 토지이용은 2009년에 제작된 환경부 중 분류 토지 피복도를 대분류 항목으로 재구성하여 이용하였 다(Fig. 2(b)). 단, 농업지역은 재배환경에 따라 강우유출현 상이 다를 것으로 예상하여 논과 밭으로 세분화하였다. 토 성자료는 국립농업과학원에서 제공하는 개략토양도를 이용 하였다(Fig. 2(c)). GIS 분석 결과에 의하면, 영산강 유역의 평균고도는 약 113 m이며, 영산강 지류유역 중 지석천 유 역의 면적이 664 km²으로 가장 크고, 황룡강, 영암천 유역 순으로 큰 것으로 분석되었다. 영산강 유역의 토지이용 형 태를 분석해 보면, 산림이 46%, 논경지 24%, 밭 14%, 도시 지역이 7%를 차지하는 것으로 분석되었다(Table 1). 주요 지류유역별로 분석해 보면, 증암천, 황룡강, 지석천 유역에 서는 산림이 60% 이상으로 많이 분포하고 있으며, 광주천 유역에서는 도시지역이 43%로 큰 비중을 차지하는 것으로 분석되었다. 영산강 하류에 위치한 삼포천과 영암천 유역에 서는 논경지의 비중 각각 47%, 38%로 크게 차지하는 것으 로 분석되었다. 영산강 유역의 토양은 식양질 내지 사양질 이 41%로 대부분을 차지하는 것으로 분석되었다.

Fig. 2. GIS data for HSPF.

기상자료는 기상청에서 운영 중인 광주지방기상청, 목포 기상대, 해남기상관측소, 장흥기상관측소에서 2011년 1년간 1시간 단위로 관측된 기상 자료를 이용하였다(Fig. 3(a)). 강우자료는 국가수자원관리 종합정보시스템(Water Management Information System, WAMIS)에서 제공하는 16개 지점, 기상청에서 운영 중인 자동기상관측소(Automatic Weather Station, AWS) 15개 지점에서 관측된 2011년 1년간 1시간 단위 자료를 이용하였다(Fig. 3(a)). 점오염원 자료는 환경 부에서 제공하는 2010년 동리별 배출유량, BOD, TN, TP 배출부하량을 소유역별로 재산정하여 모델 내 입력하였다. 오염원 자료는 행정구역별로 조사가 수행되어야 하기 때문 에 분석이 완료되기까지 상당한 시간이 소요되며, 이로 인 하여 본 연구에서 2011년 동리별 배출부하량 자료를 사용 하지 못하였다. 그러나 일반적으로 1년 사이 토지이용 및 인구가 크게 변동되지 않기 때문에 2011년의 점오염원 배 출부하량이 2010년의 점오염원 배출부하량과 큰 차이를 보 이지 않을 것으로 판단하였으며, 이러한 가정 하에 2010년 을 기준으로 분석된 점오염원 자료를 모델 내 입력하였다. 유역 밖에서의 용수 공급으로 인한 유출량은 유역 내에서 점오염원의 형태로 유출되기 때문에 점오염원 자료 입력 시 이에 대한 양도 함께 모델에 반영되었으며, 대상유역에 위치한 담양댐, 광주댐, 나주댐, 장성댐의 2011년 일 단위 하천유지용수량, 방류량 자료를 한국농어촌공사로부터 수집 하여 모델에 반영하였다.

Fig. 3. The locations of measurement stations.

모델 보검정을 위한 유량자료는 WAMIS에서 제공하는 7 개 지점에서의 2011년 기간에 대한 일단위 관측 자료를 이 용하였다(Fig. 3(a)). 유량 관측지점은 영산강 상류의 본류 에 위치한 삼지 유량 관측지점과 영산강의 지류인 풍영정 천, 광주천, 황룡강, 지석천, 고막원천, 함평천 내 위치한 월곡, 유덕, 선암, 남평, 학교, 함평 6개 유량 관측지점을 대상으로 하였다. 수질자료는 물환경정보시스템에서 제공하 는 환경부 수질측정지점 6개 지점에서 월 1회 제공하는 관 측 자료를 이용하였으며, 영산강의 지류 말단에 위치한 풍 영정천, 광주천2, 황룡강3-1, 지석천4, 고막원천2-1, 함평천2 지점을 대상으로 하였다(Fig. 3(b)).

지류유역별로 1개에서 8개의 유역으로 구분하여 모델구동 을 위한 유역 면적의 차이가 유역 간 크지 않도록 하였으며, 지류와 지류사이의 작은 유역을 포함하여 영산강유역 전체 를 총 45개 소유역으로 구성하였다(Fig. 4(a)). 모의결과를 효 과적으로 제시하기 위해 Fig. 4(b)와 같이 지류유역과 지류와 지류사이 본류구간 유역 각각에 코드를 부여하였다. 코드는 상류에서부터 하류방향으로 1~10까지는 지류유역, 11~21까지 는 영산강 본류 유역을 대상으로 순차적으로 부여하였다. 각각 의 소유역은 강우 관측소 위치, 토지이용, 토성에 따라 세그 먼트를 구분하였으며, 각 세그먼트별 동일한 매개변수를 적용 하여 영산강 전체 소유역에 대해 일관된 기준을 적용하였다.

Fig. 4. Sub-basins and sub-basin codes of the study area.

3.1. 모델 보정 및 검정

영산강 유역에서 관측된 유량과 수질 자료를 이용하여 모 델을 보정하고 검정하였다. 매개변수의 보정은 유량, SS, BOD, TN, TP 농도의 순서로 시행 착오법을 이용하였다. 유 량 보검정의 정량적 평가는 결정계수 R² (Coefficient of determination)와 NSE(Nash-Sutcliffe efficiency coefficient; ^{Nash and Sutcliffe, 1970})를 사용하여 평가하였다. NSE는 -∞ 와 1사이의 값을 갖으며, R² 값은 0 ~ 1사이의 값을 갖는다. SS, BOD, TN, TP 농도는 환경부 수질관측소에서 월 단위로 제공되기 때문에 강우 시 나타나는 피크 값을 재현하는데 한 계가 있다. 따라서 유량과 달리 정량적 평가보다는 전체적인 경향성을 재현하는 것에 초점을 두고 보검정을 수행하였다. 구체적으로는 월 단위로 제공되는 관측 값을 매월 15일에 임 의로 삽임 후 이를 기저농도로 가정하고 모의치의 기저농도 와 비교함으로써 보정을 실시하였다. 모의기간은 2011년으로 설정하였으며 보정은 1월1일부터 7월 31일까지, 검정은 8월 1일부터 12월 31일까지의 기간을 대상으로 하였다.

모델 보정을 위해 유량 및 SS, BOD, TN, TP 농도에 상 대적으로 민감하게 반응하는 매개변수에 대해 HSPF 사용자 매뉴얼(^{Bicknell et al., 2001})에서 제시하는 값의 범위 내에 서 수정하였다(Table 2). 유량 보정을 위해 첨두유량, 기저 유량, 총유출량에 큰 영향을 주는 LZSN, INFILT, AGWRC, UZSN, DEEPFR, INTFW, IRC 등의 매개변수를 조절하였다. LZSN과 UZSN은 토양층별 함수 저장능, INFILT는 토양침 투능과 관련된 매개변수이며, AGWRC와 DEEPFR는 지하수, INTFW와 IRC는 중간유출 작용의 모의를 위해 이용된다.

SS 농도 보정을 위해 매개변수 KRER, JRER, KSER, JSER, TAUCD, TAUCS 등을 조절 하였다. KRER, JRER, KSER, JSER은 PERLND에서 발생하는 강우에 의한 토양분리 (detachment)와 유실(washoff) 작용을 모의하기 위해 이용되 며, TAUCD와 TAUCS는 각각 RCHRES 모듈에서 발생하 는 점착성유사인 실트와 점토의 퇴적과 침식 작용을 모의 하기 위해 이용된다. BOD 농도 보정을 위해 PERLND 모 듈의 IFLW, GRND 등의 매개변수를 월별로 조절하였으며, RCHRES 모듈의 KBOD20, KODSET 등의 매개변수를 조 절하였다. IFLW, GRND은 각각 중간유출과 지하수에서의 BOD 농도와 관련된 매개변수이며, KBOD20, KODSET은 각각 BOD 분해와 침전작용을 모의하기 위해 이용된다. TN과 TP 농도 보정을 위해 BOD 농도 보정과 마찬가지로 IFLW, GRND 등의 매개변수를 월별로 조절하였으며, RCHRES 모듈의 KTAM20, KNO220, KNO320 등의 매개변수를 조 절하였다. KTAM20, KNO220, KNO320은 각각 NH₃, NO₂ 의 질산화와 NO₃의 탈질 작용과 관련된 매개변수이다.

수집한 유량 자료를 분석한 결과, 삼지, 학교, 함평 유량 관측소에서 관측값이 강우량 대비 20% 이하 또는 100% 이상의 비이상적인 유출률을 보였다. 이와 같이 유역에서 나타나는 일반적인 유출률보다 상대적으로 과소 또는 과대 한 수치를 보이는 경우, 그 원인을 관측값의 부정확성, 유 역 내 저수지 존재, 유역 밖에서 공급되는 용수의 영향이 클 경우 등으로 가정할 수 있다. 삼지 관측소의 경우, 상류 에 담양댐이 존재하지만 댐 유역을 제외하고 유출률 산정 했음에도 불구하고 상대적으로 낮은 유출률을 보였다. 학교 관측소의 경우, 유역 내 대규모 댐이나 환경기초시설이 부 재함에도 불구하고 상대적으로 높은 유출률을 보였으며, 함 평 관측소의 경우, 관측유량이 0인 날의 수가 상대적으로 많아 관측 자료가 부정확하다고 판단하였다. 따라서 이러한 유역의 경우 유량의 절대값을 재현하기보다 유역의 평균적 인 범위의 유출률과 첨두유량의 시간을 합리적으로 재현할 수 있도록 보정을 실시하였다. 모의유량과 관측유량의 R² 값은 보정 0.77 ~ 0.97의 범위를, 검정 0.46 ~ 0.97의 범위를 보여주었으며, NSE 지수는 보정 0.70 ~ 0.96, 검정 0.80 ~ 0.91 의 범위를 보였다(Table 4). HSPF에서 모델적용성의 평가 기준이 되는 R² 값은 Table 3과 같다(^{Donigian, 2002}). 이 기준에 따르면, 검정기간의 R² 값이 0.46을 나타낸 함평천 유역을 제외하고, 모든 유역에서 전반적으로 Good에서 Very Good의 기준을 보여주어 모의유량이 관측유량에 대 한 기저유량 크기와 첨두유량의 크기 및 시간을 합리적으 로 재현하는 것으로 나타났다(Fig. 5).

Fig. 5. Calibration and validation result of discharge in Yeongsan river basin.

SS, BOD, TN, TP 농도는 환경부 수질관측소에서 월 단위 로 제공되기 때문에 강우 시 나타나는 농도의 피크 값의 재 현성은 확인하지 못하지만 전반적으로 기저농도로 가정한 관측 값이 모의치의 기저농도와 잘 일치하여, 모의치가 관측 치의 경향을 적절히 재현하는 것으로 나타났다(Fig. 6 ~ Fig. 9).

Fig. 6. Calibration and validation result of SS concentration in Yeongsan river basin.

Fig. 7. Calibration and validation result of BOD concentration in Yeongsan river basin.

Fig. 8. Calibration and validation result of TN concentration in Yeongsan river basin.

Fig. 9. Calibration and validation result of TP concentration in Yeongsan river basin.

3.2. 유출특성 분석

3.2.1. 소유역별 강우 유출특성

모의결과에 따른 강우 유출특성을 분석해 보면, 2011년 을 기준으로 영산강 유역에서의 강우량 대비 증발산량은 46%, 하천 유출량은 58%를 차지하는 것으로 분석되었다. 소유역별로 강우량 대비 증발산량(Evapotranspiration)은 37% ~ 56%의 범위를, 하천유출량(Runoff)은 37% ~ 435%의 범위를 보이는 것으로 나타났다(Table 5). 특히, 14번과 15 번 소유역에서 강우량 대비 100%가 넘는 유량이 유출되는 것을 보여주었는데, 이는 영산강 유역 밖에서 유역 내로 공급하는 용수의 영향 때문인 것으로 판단된다. 광주시와 목포시 일원은 섬진강 유역 내 위치한 주암본댐에서 64만 m³/day의 생활, 공업 및 관개용수를 공급받고 있으며 (^{NIER, 2011}), 주암본댐에서 공급하는 용수의 영향으로 인 해 일부 소유역에서 강우량 대비 100%가 넘는 유출량을 보이는 것으로 판단된다. 강우 유출량이 상대적으로 많은 14번 소유역은 풍영정천 유입 후, 황룡강 유입 전 영산강 본류를 둘러싼 유역이며, 15번 소유역은 황룡강 유입 후, 지석천 유입 전 영산강 본류를 둘러싼 유역이다. 이들 유 역 내에는 영산강 유역에서 가장 큰 처리장인 광주하수처 리장이 위치해 있다. 광주하수처리장은 일평균 방류량이 600 m³/day로 갈수기에는 본류 유량의 절반정도를 차지하 며, 영산강 본류 수질에 상당한 영향을 미치는 것으로 보 고되고 있다(^{Park, Shin et al., 2012}). 강우량 대비 저류량 (storage)은 -374% ~ 10%의 범위를 보이는 것으로 분석되었 다. 저류량은 유역에서의 하천 유출량과 증발산량을 제외한 나머지 유량에 대한 공급 또는 저류 현황을 나타내며, 저 류량은 광역상수도를 통해 유역 밖에서 유역내로 공급되는 용수의 양이 많을 경우 음의 값을, 유역 내 댐/저수지에 의 해 저류되는 양이 많을 경우 양의 값으로 나타난다(Table 5). 대규모 농업용 저수지인 광주댐, 장성댐, 나주댐, 담양댐 등이 위치한 1, 4, 5, 11번 유역에서 저류량은 4% ~ 10%의 범위를 보이는 것으로 분석되었다.

Table 5. Water balance of sub-basins for Yeongsan river (%)

	Total	Sub-Basin codes
		1	2	3	4	5	6	7	8	9	10	11	12	13	14	15	16	17	18	19	20	21
Evapo transpiration	46	52	41	38	37	45	49	49	47	49	47	46	50	41	39	46	47	47	49	54	54	56
Runoff	58	37	62	62	55	51	52	54	55	55	55	44	57	61	435	116	55	56	54	48	49	48
Storage	-3	10	-4	1	8	4	-1	-3	-2	-3	-1	10	-7	-2	-374	-62	-2	-3	-3	-2	-2	-4

모의결과에 따른 유출량을 정량적으로 분석해 보면, 2011 년을 기준으로 영산강유역에서 영산호로 약 24억 m³의 유 량이 유출되는 것으로 분석되었다(Table 6). 영산강의 지류 유역 중 유출량이 가장 많은 유역은 지석천과 황룡강 유역 으로, 각각 약 연 4.5억 m³, 4.0억 m³의 유출량을 보이며, 이는 영산강 유역의 총 유출량 중 19%, 17%의 비중에 해 당한다(Table 6, Fig. 10(a)). 단위면적당 유출 유량은 14번 소유역에서 가장 많은 것으로 분석되었다(Fig 10(b)). 이는 하천으로 유입하는 하수처리장 방류량 영향 때문이며, 14번 소유역은 2011년을 기준으로 평균적으로 하천 유량의 약 86% 가 점오염원 배출 유량이 차지하는 것으로 분석되었다.

Table 6. The runoff and pollution loads of sub-basins for Yeongsan river

		Total	Sub-Basin codes
			1	2	3	4	5	6	7	8	9	10	11	12	13	14	15	16	17	18	19	20	21
Runoff	(106m3/yr)	2406	62	55	87	398	450	61	130	118	79	165	137	28	45	228	63	42	61	70	8.3	55	63
	(%)	100	2.6	2.3	3.6	16.6	18.7	2.5	5.4	4.9	3.3	6.9	5.7	1.2	1.9	9.5	2.6	1.7	2.6	2.9	0.3	2.3	2.6
Sediment	(Mg/yr)	152	2.84	3.99	4.13	34.9	25.3	4.02	9.07	9.76	6.87	13.5	6.69	1.99	3.27	2.59	2.75	3.15	3.76	5.25	0.50	3.10	4.96
	(%)	100	1.9	2.6	2.7	22.9	16.6	2.6	6.0	6.4	4.5	8.8	4.4	1.3	2.1	1.7	1.8	2.1	2.5	3.4	0.3	2.0	3.3
BOD	(ton/yr)	15,721	374	400	535	2,129	2,021	388	913	673	550	974	779	219	308	995	308	550	959	756	104	868	915
	(%)	100	2.4	2.5	3.4	13.5	13.5	2.5	5.8	4.3	3.5	6.2	5.0	1.4	2.0	6.3	2.0	3.5	6.1	4.8	0.7	5.5	5.8
TN	(ton/yr)	10,071	249	171	340	911	911	250	518	451	383	613	483	135	164	1,305	393	332	566	413	54	462	504
	(%)	100	2.5	1.7	3.4	9.0	9.0	2.5	5.1	4.5	3.8	6.1	4.8	1.3	1.6	13.0	3.9	3.3	5.6	4.1	0.5	4.6	5.0
TP	(ton/yr)	563	10.7	8.5	28.8	39.1	39.1	10.7	19.5	18.8	15.8	28.1	22.5	6.5	7.9	91.7	17.6	31.7	54.4	26.6	3.3	35.1	35.9
	(%)	100	1.9	1.5	5.1	7.0	7.0	1.9	3.5	3.3	2.8	5.0	4.0	1.2	1.4	16.3	3.1	5.6	9.7	4.7	0.6	6.2	6.4

Fig. 10. The runoff of sub-basins for Yeongsan river basin.

3.2.2. 소유역별 오염부하 유출량 분석

영산강 유역에서 2011년 연간 배출된 유사, BOD, TN, TP 유출량은 각각 152,327 ton, 15,721 ton, 10,071 ton, 563 ton으로 산정되었다(Table 6). 유사와 BOD는 유역면적 이 상대적으로 큰 황룡강 유역과 지석천 유역에서 많은 양 이 유출되며, 두 소유역에서의 유사유출량은 영산강유역 전 체 유출량의 약 40%를, BOD 유출량은 약 26%를 차지하 는 것으로 나타났다(Table 6, Fig. 11(a), (b)). TN은 지석천 유역과 14번 소유역에서 많은 양이 유출되며, 두 소유역에 서의 TN 유출량은 영산강유역 TN 유출량의 약 27%를 차 지하는 것으로 분석되었다. TP는 14번과 17번 소유역에서 많은 양이 유출되며, 두 소유역에서의 TP 유출량은 영산강 유역 TP 유출량의 약 26%를 차지하는 것으로 나타났다 ((Table 6, Fig. 11(c), (d)). 유역면적은 상대적으로 작으나 14번과 17번 소유역에서 TN과 TP 유출량이 많은 것은 영 산강으로 방류되는 광주하수처리장 처리수의 영향 때문인 것으로 판단된다. ^{Park, Shin et al. (2012)}은 EFDC 모델을 이용하여 광주하수처리장 처리수의 유입 유무에 따라 영산 강 본류 수질을 모의한 결과 광주하수처리장이 영산강 본 류 수질에 상당한 영향을 미치는 것으로 보고하였으며, ^{Jung et al. (2012)}은 영산강 본류 10개 지점에서 10년 간 관측된 수질 자료에 대하여 다변량분석법을 이용하여 요인 분석을 수행한 결과 영산강의 주요 오염원인 광주천, 도시 하수 및 광주하수처리장 처리수에서의 질소, 인 농도가 영 산강 수질에 미치는 영향이 매우 크다고 보고하였다. 본 연구에서의 질소와 인의 모의 결과를 분석해 보면, 광주 제1하수처리장이 위치한 14번 소유역에서 2011년을 기준으 로 산정된 NH₄-N 유출량 중 점오염원이 차지하는 비중은 약 79% 이며, PO₄-P 유출량 중 점오염원이 차지하는 비중 은 약 85%인 것으로 분석되었다. 이는 선행연구와 마찬가 지로 광주하수처리장 처리수가 영산강 본류 수질에 상당한 영향을 미치고 있음을 시사한다.

Fig. 11. The pollution loads of sub-basins (ton/yr).

단위면적당 유출되는 유사, BOD, TN, TP 부하량을 산정 한 결과, 유사는 소유역별로 190 ~ 617 kg/ha/yr, BOD는 25 ~ 232 kg/ha/yr, TN은 16 ~ 305 kg/ha/yr, TP는 0.7 ~ 21 kg/ha/yr 의 범위를 보이는 것으로 분석되었다(Fig. 12). 단위면적당 유사는 황룡강 유역과 본류에 위치한 14번, 15번 소유역에 서 상대적으로 많은 양이 유출되는 것으로 분석되었으며, 단위면적당 BOD, TN, TP 부하량은 광주하수처리장의 영 향을 가장 많이 받는 14번 소유역에서 가장 많이 유출되는 것으로 분석되었다.

Fig. 12. The pollution loads of sub-basins per unit area (kg/ha/yr).

3.2.3. 월별 유출특성

2011년 영산강 유역에서 유출되는 유량, 유사, BOD, TN, TP 부하량은 7월에 가장 많은 양이 유출된 것으로 분석되 었다(Fig. 13). 유량은 5월, 7월, 8월에 각각 전체 유출량의 10% 이상이 유출되는 것으로 분석되었으며, 특히, 7월과 8 월에 연간 총 유출량의 44%가 유출되는 것으로 분석되어 영산강 유역에서 영산호로 유입되는 유량이 우기에 집중적 으로 유출되는 것을 보여주었다. 유사는 연간 총 유출량의 63%가 7월과 8월, 강우기에 집중적으로 유출된 것으로 분 석되었으며, BOD, TN, TP는 7월, 8월, 5월, 11월에 연간 총 부하량의 10% ~ 20%의 범위로 유출되는 것으로 분석되 었다. 이는 유사 유출량이 다른 수질 항목에 비해 상대적 으로 강우에 의한 영향을 많이 받을 수 있음을 의미한다. 한편, BOD, TN, TP 유출량은 강우에 의한 유출량이 많지 않은 1월, 12월 등에도 연간 총 유출량의 3% 이상 유출되 는 것으로 분석되었으며, 이는 강우와 무관하게 하천으로 지속적으로 유입하는 점오염원의 영향인 것으로 판단된다.

Fig. 13. Monthly pollutant load in Yeongsan river basin in 2011 (ton/month).