1. 서 론
2. 재료 및 방법
2.1 대상유역
2.2 SWAT 모형의 개요
2.3 자료수집 및 분석 방법
2.3.1 GIS 및 기상, 수문입력자료
2.3.2 지하수위 실측자료
2.3.3 기후변화 시나리오
2.3.4 연구 방법
3. 결과 및 고찰
3.1 SWAT 모형의 검 ․ 보정
3.2 기후변화에 따른 유출특성 변화
3.3 기후변화에 따른 지하수위 변화
4. 요약 및 결론
1. 서 론
IPCC 4차 보고서(2007)에 따르면, 미래 기후변화로 인한 가장 취약한 부분으로 강수패턴의 시공간 변화로 인한 가용 수자원의 변화를 선정하였으며 특히 IPCC 5차 보고서(2013)에서 아시아지역은 지역별 대응전략 수립, 물 재활용 등 수자원 다양화, 통합형 수자원 관리를 권고하였다. 기후변화 에 따른 영향 중에서 지하수 수문계(Hydrological system)의 변화는 광범위한 부분에서 관측되고 있으며, 수자원의 대부분을 지하수에 의존하는 지역의 경우 특히 민감한 문제가 발생되고 있다(Song et al., 2012).
우리나라의 연간 지하수 사용량은 약 37억 m3으로, 전체 수자원 사용량의 10% 정도에 불과하나, 국내 전체 지하수 이용 가능량의 35%를 상회하고 있다(Lee et al., 2011). 특히 금강유역 내 지하수이용량은 최근 10년(2005~2015)동안 11.6% 증가한 것으로 조사되었다(MOLIT, 2016). 그러나 전국 국가지하수관측망의 과거 10년의 지하수위 자료를 통해 향후 20년간 약 0.58 m의 지하수위 감소를 예측하였으며(Kim et al., 2008), 이러한 원인으로 과다한 지하수 채수 및 기후변화로 인한 함양량 감소임을 정성적으로 추정하였다.
지하수의 변화와 같이 흐름속도가 느리고 지속적인 요소의 경우에는 지표 기후변화의 영향을 쉽게 인식 할 수 없으나 지표변화에 따른 변동이 지하수 환경에서 관측되기 시작하면 그 영향은 지표보다 훨씬 장기적으로 나타남에 따라 지하수의 유출 혹은 유입량을 정량화하려는 시도는 예전부터 지속적으로 이루어져 왔다(Church, 1996). 수문모형을 이용한 지하수 거동평가 국내연구사례를 살펴보면 Kim et al. (2000)은 토양수분추적에 의한 지하수 함양량 추정기법을 적용한 일단위 분포형 모형을 개발하여 보청천 유역에 적용하였으며, Kwon et al. (2005)은 유역의 지형 및 토양특성을 고려한 분포형 일단위 강우유출모형을 이용하여 지하수 개발량을 평가하였다. Kim et al. (2007), Kim et al. (2008), Chung et al. (2010)과 Chung et al. (2011) 등은 SWAT 모형과 MODFLOW 모형을 연계하여 강수로부터 지표유출, 측방유출, 지하수 유출 성분을 계산하고 지하수 함양량을 산정하여 지하수의 정량적인 분석을 수행하였다. Kim et al. (2014)은 제주지역을 대상으로 지하수위 변동 해석모델을 이용한 수리지질 매개변수를 추정하였고 Ahn et al. (2016)은 SWAT 모형의 HRU를 이용하여 지하수위를 검 ․ 보정하여 SWAT 모형 내 지하수위 산정기법을 제시한 바 있다. 국외에서는 Shirmohammadi et al. (2013), Pfannerstill et al. (2014), Izady et al. (2015), Jin et al. (2015), Sun et al. (2016) 등이 SWAT 모형의 모의 결과를 이용하여 지하수 함양량 및 지하수 흐름을 추정한 연구를 수행하였다.
기후변화에 따른 미래 지하수 함양량을 평가한 연구로는 Ahn et al. (2009)이 SLURP 모형을 이용하여 경안천 유역에 대한 지하수 함양량을 산정한 바 있고, Lee et al. (2011)은 상세화된 기후변화 시나리오와 GIS를 이용하여 기후변화 함양량 산정모델을 개발하여 미래 기후변화 시나리오에 따른 전국의 지하수함양량을 추정하였다. Lee et al. (2014)은 해안유역을 대상으로 SWAT 모형을 이용하여 기후변화에 따른 급경사지 유역에서의 지하수 함양량을 전망하였으며, Park et al. (2014)은 기후변화 시나리오와 SWAT 모형을 이용하여 용담댐 유역의 수문요소분석을 통해 지하수유출량을 분석한 바 있다. SWAT을 이용한 지하수분석 관련 국외 연구로는 Dessu et al. (2013)이 기후변화에 따른 유황분석을 및 수문분석을 통해 지하수 흐름을 분석하였고, Awan et al. (2014)은 기후변화에 따른 관개지역의 지하수 함양 분포도를 작성하였다. Fiseha et al. (2014), Verma et al. (2015)은 기후변화에 따른 수문분석을 통해 지하수 함양량을 산정하였고, Goodarzi et al. (2016)은 SWAT 및 MODFLOW 모형을 연계하여 기후변화에 따른 지하수 함양량을 추정하였다. 이처럼 국내 ․ 외에서 수문모형을 통한 지하수 분석 및 기후변화 영향 평가 연구는 활발히 진행되어 왔으나, 기후변화에 따른 지하수위 거동을 평가한 연구는 아직 많이 이뤄지지 않았다.
본 연구의 목적은 물리적 기반의 준분포형 연속강우유출모형인 SWAT 모형을 이용하여 금강유역(9,645.5 km2) 내에 있는 다목적 댐 및 다기능 보의 운영자료 및 지하수위 관측자료를 이용한 검 ․ 보정을 통해 신뢰성 있는 수문순환을 재현하고, 기후변화 시나리오를 적용하여 미래 유출특성 및 지하수위를 추정하는데 있다.
2. 재료 및 방법
2.1 대상유역
본 연구의 대상유역은 금강유역으로 전라북도 장수군에서 발원하여 무주, 금산, 영동, 옥천, 보은, 청주, 대전, 세종, 공주, 청양, 논산, 부여, 서천, 익산을 지나 군산만에서 서해로 유입된다. 유역면적은 9,645.5 km2 이며, 유로연장은 384.8 km이다. 금강유역 내 에는 2개의 다목적 댐인 용담댐(YDD), 대청댐(DCD)과 3개의 다기능 보인 세종보(SJW), 공주보(GJW), 백제보(BJW)가 설치되어 운영 중이다. 유역 전체 면적의 62%가 산림으로 이루어져 있으며, 논이 15%, 밭은 11%를 차지하며 나머지 도시와 초지, 나지 등은 3~5%를 차지한다.
금강유역의 과거 연평균 강수량은 1,323.1 mm, 연평균 기온은 12.2°C이다. 2015년 기준 금강유역의 지하수 이용량은 525.9 백만 m3/년이고, 단위면적당 지하수 이용량은 5개 대권역(한강, 낙동강, 금강, 섬진강, 영산강) 중 2번째로 많은 1,808.3 m3/일/km2로 조사되었다(MOLIT, 2016). Fig. 1은 SWAT 모형에 적용하기 위해 분할된 78개의 소유역과 댐, 보 및 지하수위 검 ․ 보정지점을 나타낸 것이다.
2.2 SWAT 모형의 개요
본 연구에서는 장기간에 걸친 유역의 일 유출을 모의하기 위하여 미농무성 농업연구소 (USDA, Agricultural Research Service, ARS)에서 개발한 SWAT (Arnold et al., 1996)을 선정하였다. SWAT 모형은 물리적 기반의 준분포형 연속강우유출모형으로서, 다양한 종류의 토양과 토지이용 및 토지관리 상태에 따른 유출을 모의한다. 모형은 물수지 방정식에 근거를 두고 강수, 증발산량, 지표유출, 기저유출, 지하수 등에 대한 모의를 각 수문반응단위 (Hydrological Response Unit, HRU) 별로 계산할 수 있다 (Neitsch et al., 2001).
SWAT 모형에서 일별 유출량은 SCS (Soil Conservation Service) CN (Curve Number) 방법으로 산정하고 측방유입량은 Kinematic Storage Model을 이용하며, 침투는 토층을 최대 10개 까지 세분화 하여 선형저수량 추적기법을 사용하여 계산한다. 지하수는 얕은 대수층과 깊은 대수층으로 나뉘며 얕은 대수층은 지표면에서 대수층으로 침루된 유량으로 유역의 하천 회귀수를 공급하고, 얕은 대수층으로 침투된 일부는 깊은 대수층으로 공급되어 유역 밖으로 흘러나간다. 이때 지하수 충전량 및 지하수위는 아래의 식을 이용하여 계산된다.
(1)
(2)
여기서, 
은 매일 하천으로 유입되는 지하수 혹은 기저유출(mm H2O), 
는 대수층의 수리전도도(mm/day), 
는 지하수시스템에 대한 분할경계 및 소유역에서 주 하천까지의 거리(m), 
는 
일의 지하수위(m), 
은 
일의 지하수위, 
는 기저유출 하단계수, 
는 단위시간 (1 일), 
는 
일에 대수층으로 들어가는 침투량, 그리고 
는 얕은 대수층의 비산출율 이다.
2.3 자료수집 및 분석 방법
2.3.1 GIS 및 기상, 수문입력자료
SWAT 모형의 GIS 공간자료로 DEM (Digital Elevation Model)과 토양도는 국가수자원관리종합정보시스템(WAter Management Information System, WAMIS)의 30 m × 30 m DEM, 1:25,000 정밀토양도를 사용하였다(Figs. 2(a) and 2(b)). 토지이용도는 환경부의 2008년 중분류 자료를 이용하였으며, 이를 9개 항목(활엽수림, 침엽수림, 혼효림, 주거지, 논, 밭, 초지, 나지, 수역)으로 분류하여 모형에 입력하였다(Fig. 2(c)). 금강유역 내 토양종류에서 사양토(Sandy loam)와 양토(Loam)가 각각 58%, 24%로 높은 비율을 차지하고 있으며, 토지이용은 산림이 62%, 농지는 15%로 이루어져 있다.
SWAT의 검보정을 위한 기상자료는 9개(청주, 대전, 추풍령, 군산, 보은, 천안, 부여, 금산, 장수)강수량(mm), 최고 및 최저기온(°C), 풍속(m/sec), 상대습도(%), 일사량(MJ/m2) 기상관측소의 11년(2005~2015)동안의 일자료를 사용하였으며, 이 기간의 용담댐과 대청댐의 일별 댐유입량과 저수량 자료, 세종보, 공주보, 백제보는 3년 5개월(2012년 8월~2015년 12월)동안의 일별 보유입량과 저수량 자료를 수집하여 사용하였다(http://www.water.or.kr).
2.3.2 지하수위 실측자료
유역의 수문 및 물수지 평가에 있어 지하수위는 수문학적으로 중요한 부분을 차지하지만 관측자료가 충분하지 않아 경험식에 의한 추정 값을 사용하고 있는 실정이다(Ahn et al., 2016). 따라서 본 연구에서는 모델의 검증을 위해 유역 내에서 측정된 지하수 값을 이용하였다. 지하수위 자료는 국가지하수정보센터(GIMS)의 충적 지하수 관측정에 의해 관측된 지하수위 자료를 수집하였다. 금강유역 내 충적 지하수 관측정 18개소가 운영 중으로 모형의 보정 및 검정기간(2005~2015)에 대한 2년 이상의 결측자료가 없는 5개 지점(장수장수(JSJS), 옥천청성(OCCS), 보은마로(BEMR), 천안수신(CASS), 부여부여(BYBY))을 선정하였다(Fig. 1).
2.3.3 기후변화 시나리오
본 연구에서는 기상청에서 제공하는 HadGEM3-RA RCP (Representative Concentration Pathway) 4.5와 8.5 시나리오 자료(강수량, 최저/최고온도, 풍속, 상대습도)를 이용하였으며, 2020s (2010~2039년), 2050s (2040~2069년), 2080s (2070~2099년)로 구분하여 정리하였다.
Fig. 3과 Table 1은 지상 관측자료(1975~2005, Baseline), RCP 4.5와 8.5 시나리오의 과거자료(1975~2005, Historical), 그리고 2020s, 2050s, 2080s의 월별 강수량과 기온을 비교한 것이다. 2080s RCP 8.5 시나리오에서 미래 기온과 강수량은 Historical에 비해 각각 5.3°C, 11.5% (160.3 mm) 증가하였으며, 봄철에 각각 5.5°C, 63.5% (141.9 mm) 증가, 여름철은 6.1°C 증가, 5.3% (44.3 mm) 감소하는 경향이었다.
Table 1. Future projected seasonal precipitation (mm) and temperature (°C) under the RCP 4.5 and 8.5 scenarios of the HadGEM3-RA  |
2.3.4 연구 방법
본 연구의 방법은 Fig. 4와 같으며, SWAT 모형의 입력자료 및 검 ․ 보정 지점에 대한 관측자료를 수집하였다. 댐 및 보에 대한 제원과 방류량자료를 구축하여 입력한 후 유역의 상류에서 하류로 순차적으로 댐 및 보의 유입량 및 저수량을 검 ․ 보정을 수행했다.
SWAT에서 지하수위는 모의 결과 값으로 출력하지 않음에 따라, 본 연구에서는 지하수위 관측정이 설치된 지점의 HRU의 계산결과를 추출하여 관측자료와 비교했는데 HRU 출력 결과 중 SA_ST (shallow aquifer storage, mm)에서 매개변수인 GWQMN (threshold water level in shallow aquifer for base flow, mm)을 뺀 값을 지하수위 변동량으로 보고, 실측 지하수위 초기 값에 변동량을 더하여 실측자료와 비교하였다. 실측자료와의 비교를 통해 매개변수를 조정했다. 매개변수 추정은 수문요소들을 유기적으로 밸런스 있게 모의하기 위해 시행오차방법(trial and error method)을 적용했다. 기후변화에 따른 지하수위 거동평가를 위해 기상청 RCP 자료를 수집하였으며, 평가기간에 따른 모의 결과를 도출하고 미래 지하수위를 추정하였다.
3. 결과 및 고찰
3.1 SWAT 모형의 검 ․ 보정
본 연구에서는 모형의 검 ․ 보정에 앞서 금강유역에 대한 매개변수 민감도를 수행하였으며, 이를 통해 최적의 매개변수를 선정하였다. Table 2는 모형의 보정 및 검증에 사용된 매개변수를 정리한 것으로 지표유출 관련 매개변수는 CH_N (2), 증발산량 관련 매개변수는 ESCO, 지하수 관련 매개변수는 GW_DELAY, GWQMN, ALPHA_BF, REVAPMN, GW_ REVAP 이다. 저수지 운영매개변수는 7개(RES_ESA, RES_ EVOL, RES_PSA, RES_PVOL, RES_VOL, RES_K, EVRSV)를 선정하고 매개변수의 추정과 검 ․ 보정 방법은 Ahn et al. (2016)이 제안한 방법을 적용하였다.
SWAT 모형의 보정기간은 5년(2005~2009), 검증기간은 6년(2010~2015)으로 설정하였으나, 3개의 보(세종보, 공주보, 백제보)는 2012년 8월부터 운영됨에 따라 2013년을 보정, 2014~2015년을 검증기간으로 하여 분석을 수행하였다. 모형의 적용성 평가는 coefficient of determination (R2), Nash and Sutcliffe (1970) model efficiency (NSE), root-mean square error (RMSE), percent bias (PBIAS)를 사용했다. R2 은 1에 가까울수록 관측 값과 모의 값이 완벽하게 일치함을 나타내고, NSE는 -∞에서 1까지의 범위를 통해 모형의 효율을 평가하는데 0보다 클수록 모형의 적용성이 높다고 설명할 수 있다(Mkhwanazi, 2012). RMSE는 실측치와 모의치의 오차를 나타내는 것으로 0에 가까울수록 오차가 작음을 의미한다. PBIAS는 관측치와 모의치의 오차를 퍼센트로 나타난 것으로 오차값이 작을수록 모형의 효율이 높음을 나타낸다.
Fig. 5는 5개의 댐, 보 지점에 대한 유입량 및 저수량 검 ․ 보정결과를 나타낸 것이며, Table 3은 검 ․ 보정 결과 중 유입량과 저수량에 대한 통계분석 값을 정리한 것이다. 검 ․ 보정 전체기간에 대한 YDD, DCD의 댐 유입량 비(QR)는 실측치보다 +1.4%, +3.4%, 저수량은 실측치보다 2.2%, 3.7% 많게 추정되었다. SJW, GJW 보 유입량 비는 실측치에 비해 +1.4%, +0.4%, 많게 추정되었다. 다만 금강 하류에 위치한 BJW의 경우 -4.3% 적게 모의되었는데 상류인 SJW, GJW의 유입량이 하류인 BJW에 영향을 미친 것으로 판단된다.
댐 유입량 통계분석 결과를 살펴보면, R2는 YDD, DCD, SJW, GJW, BJW 에서 각각 보정기간에 대하여 0.67, 0.75, 0.75, 0.80, 0.81, 검증기간에 대하여 0.70, 0.71, 0.64, 0.61, 0.62로 분석되었다. NSE는 YDD, DCD, SJW, GJW, BJW 에서 각각 보정기간에 대하여 0.55, 0.70, 0.69, 0.76, 0.77, 검증기간에 대하여 0.58, 0.64, 0.64, 0.57, 0.50으로 분석되었다. 검정기간 중 가뭄년인 2014~2015년은 평균 강수량이 941.0 mm로 R2와 NSE가 낮게 나타났다. 전체기간에 대한 RMSE는 YDD, DCD, SJW, GJW, BJW 에서 각각 1.67 mm/day, 1.42 mm/day, 0.94 mm/day, 1.02 mm/day, 1.18 mm/day 였으며, PBIAS는 각각 4.45%, 7.86%, 4.67%, -3.35%, -2.06%의 통계치를 보였고 모든 검 ․ 보정 지점에 대해 통계분석 결과가 통계적으로 유의한 수준인 것으로 분석되었다.
댐 저수량 통계분석 결과 중, R2는 YDD, DCD, SJW, GJW, BJW 에서 각각 보정기간에 대하여 0.68, 0.65, 0.55, 0.55, 0.57, 검증기간에 대하여 0.78, 0.68, 0.70, 0.72, 0.31로 분석되었다. 검정기간 중 가뭄년인 2014~2015년은 평균 강수량이 941.0 mm로 R2가 낮게 나타났다. NSE는 YDD, DCD, SJW, GJW, BJW 에서 각각 보정기간에 대하여 0.79, 0.95, 0.92, 0.98, 0.99, 검증기간에 대하여 0.94, 0.93, 0.98, 0.99, 1.00으로 분석되었다. 전체기간에 대한 RMSE는 YDD, DCD, SJW, GJW, BJW 에서 각각 1.85 mm/day, 2.07 mm/day, 0.13 mm/day, 0.26 mm/day, 0.23 mm/day 였으며, PBIAS는 각각 5.48%, -9.38%, 3.89%, 1.44%, 0.36%의 통계치를 보였고 모든 검 ․ 보정 지점에 대해 통계분석 결과가 통계적으로 유의한 수준인 것으로 분석되었다(Table 3).
지하수위의 변동량 경우, 해당지점의 HRU의 계산결과를 검 ․ 보정에 활용하였다. 수문학적 반응단위인 HRU는 소유역 내 동일한 토지이용과 토양특성의 조합으로 구성되는데, OCCS, BEMR, BYBY는 산림지역(혼효림), JSJS, BYBY는 논 지역으로 분류하여 분석을 수행하였다. 검 ․ 보정 기간 동안의 평균 R2는 JSJS, OCCS, BEMR, CASS, BYBY 각각 0.57, 0.60, 0.53, 0.63, 0.61로 분석되었다(Table 4). 모의값과 실측값의 통계분석 결과로 부터 댐 보 유입량 검 ․ 보정 결과와 동일하게 가뭄을 겪은 해인 2014, 2015년의 5개 지점 R2 평균은 0.40 이었으며, 이는 가뭄 년을 제외한 기간(2005~2013)의 평균인 0.64보다 모형 효율이 낮음을 확인할 수 있었다.
Fig. 6은 지하수위 검 ․ 보정 결과로 논 지역(JSJS, CASS)에 비해 산림지역(OCCS, BEMR, BYBY)에서 강우량에 따른 지하수위 상승 반응이 빠른 것으로 분석되었다. 금강 최상류 지점인 JSJS의 경우 2008~2009년 관측 자료에서 수위가 낮아졌으나 SWAT 모의는 일정하게 유지하는 결과를 나타내었는데, 이는 Table 2에서와 같이 지하수 관련 매개변수인 GW_DELAY를 150 day로 적용하였는데 매개변수의 영향으로 지하수 흐름이 소유역내에서 천천히 흐르도록 계산되었기 때문으로 판단된다. 모의값과 실측값의 통계분석 결과로부터 SWAT은 댐, 보의 유입량 및 저수량 및 지하수위 변통 특성을 잘 반영하고 있는 것으로 판단된다.
3.2 기후변화에 따른 유출특성 변화
기후변화에 따른 지하수위 분석에 앞서 검 ․ 보정기간(2005~2015)의 매개변수를 Historical (1975~2005), 2020s (2010~2039년), 2050s (2040~2069년), 2080s (2070~2099년) 기간에 적용하여 모의하였다. Table 5는 기후변화 시나리오를 사용하여 SWAT 모형을 구동한 결과(강우량, 온도, 총유출, 지하수함양량, 증발산량)를 유역평균값으로 정리한 것이다. 기상인자에 대한 비교를 위해 토지이용도는 2.3.1 장에서 언급한 바와 같이 동일한 자료를 적용하였다.
강수량의 경우 Historical 대비 모든 시나리오 및 평가기간에서 증가하였는데, 특히 RCP 4.5 2080s에서 13.6%로 가장 큰 증가폭을 나타내었다. 강수량의 영향으로 동일한 시나리오의 평가기간에서 총 유출량과 지하수함양량은 Historical 대비 각각 14.6%, 10.2% 증가하였다.
증발산량의 경우 모든 시나리오 및 평가기간에서 증가하였는데, RCP 8.5 2080s에서 17.2%로 가장 크게 증가한 것으로 분석되었다. 이와 더불어 유역 평균 기온의 경우 Historical 대비 5.4°C 상승한 것으로 나타났는데, 이는 기후변화에 따른 기온의 상승으로 증발산량이 증가함에 따라 총 유출 및 지하수함양량에 영향을 준 것으로 판단된다.
Fig. 7은 Historical 대비 지하수함양량의 변화를 표준유역별로 도시한 것으로 대청댐 상류 보청천 유역의 경우 모든 기후변화 시나리오 및 평가기간에서 지하수함양량이 감소하는 것으로 분석되었다. 특히 가장 크게 감소한 RCP 8.5 시나리오의 2020s를 공간적으로 살펴보면 대청댐 상류 뿐만아니라 금강 하류 논산천 하류 주변의 농경지를 포함하는 소유역에서 지하수함양량이 감소하는 것으로 나타났다. 2014년 기준 금강유역 내 지하수이용량 중 49.1%가 농업용수로 사용됨에 따라(MOLIT, 2014) 금강유역의 지하수는 농작물 재배에 크게 기여하고 있다. 이에 따라 기후변화로 인한 지하수함양량의 부족은 농작물 피해 등의 2차 피해를 야기 할 것으로 생각된다.
Table 5. Summary of the future predicted annual hydrological components by climate change scenarios (values in parentheses indicate change based on historical)  |
3.3 기후변화에 따른 지하수위 변화
기후변화에 따른 지하수위 변동량은 5개의 지하수위 검 ․ 보정지점에 대한 동일한 HRU를 추출하여 정리했다. Table 6은 지하수위 변화를 시나리오와 평가기간에 따른 계절별로 정리한 것이고, Fig. 8은 월별 그래프로 나타낸 것이다. 계절은 봄은 3~5월, 여름은 6~8월, 가을은 9~11월, 겨울은 12~2월로 구분했다.
금강 최상류에 위치한 JSJS에서 RCP 8.5 시나리오의 2020s, 2050s, 2080s의 연평균 지하수위는 각각 Historical 대비 4.0 cm, 3.0 cm, 5.0 cm 감소한 것으로 분석되었다. 봄의 지하수위는 상승하였으나, 여름에서의 지하수위가 크게 감소(2020s: 17.0 cm, 2050s: 16.0 cm, 2080s: 24.0 cm)하여 가을과 겨울에 잇달아 영향을 주었기 때문인 것으로 생각된다(Table 5). 대청댐 상류의 지하수관측 지점인 OCCS과 BEMR은 RCP 4.5와 8.5 시나리오 모두 연평균 지하수위가 감소하는 것으로 나타났는데, 특히 5월의 지하수위가 급격히 감소하는 경향을 나타내었다(Fig. 8). OCCS과 BEMR의 HRU 토지이용 속성은 혼효림으로 강수량은 증가했으나, 기온의 상승으로 산림의 수분이용효율이 변화함에 따라 증발산량의 증가로 총 유출 및 지하수함양량이 감소하여 지하수위 역시 감소한 것으로 판단된다.
CASS은 여름을 제외하고 모든 시나리오의 평가기간에서 지하수위가 증가하는데. Fig. 3에서와 같이 미래 여름 기간 중 7월 강수량은 Historical 대비 모두 증가하나 8월은 감소하여 여름의 지하수위가 감소한 반면, 가을은 강수량이 전체적으로 증가됨에 따라 지하수위도 같은 패턴을 나타내었다(Fig. 8). 금강유역 하류에 위치한 BYBY의 모의 결과는 모든 시나리오의 평가기간에서 가을의 지하수위가 모두 증가하였는데, 이는 BYBY가 유역 하류로 경사가 완만한 지역에 위치하고 있어 여름철 강수량으로 인한 지하수가 천천히 흘러 가을까지 영향을 준 것으로 생각된다.
본 연구에서의 5개관측정에 미래 지하수위 변화를 종합하면, 봄(3~5월), 여름(6~8월)의 지하수위는 기온상승에 따른 증발산량의 증가로 인해 감소한 반면 가을(9~11월)은 강수량의 증가로 지하수위가 상승하였고 연평균 지하수위 변동은 RCP 4.5와 8.5 각각 -9.0~8.0 cm, -13.0~4.0 cm 의 범위로 RCP 8.5에서의 지하수위가 전체적으로 낮음을 확인할 수 있었다.
4. 요약 및 결론
본 연구는 물리적 기반의 준분포형 연속강우유출모형인 SWAT 모형을 이용하여 금강 수계(9,645.5 km2) 내에 있는 다목적 댐 및 다기능 보의 운영자료 및 지하수위 관측자료를 이용한 검 ․ 보정을 시도하여 신뢰성 있는 수문순환을 재현하였고 기후변화 시나리오를 적용하여 미래 지하수위를 추정하였다. 본 연구에서 도출된 주요 연구결과는 다음과 같다.
1)금강유역 내 2개의 다목적댐(용담댐, 대청댐) 3개의 다기능 보(세종보, 공주보, 백제보) 지점에 대한 유입량 및 저류량 자료를 이용하여 보정(2005~2009년) 및 검증(2010~ 2015년) 하였다. 모형의 적용성 평가 결과 R2는 0.85, NSE는 0.66, RMSE는 1.24 mm/day, PBIAS는 2.32%의 통계치를 보였고 모든 검 ․ 보정 지점에 대해 통계분석 결과가 통계적으로 유의한 수준인 것으로 분석되었다.
2)기후변화 시나리오의 과거자료(historical)와 과거 기상자료(baseline)과 유사성을 가지도록 편의보정 후, Histo-rical을 기준으로 미래 세 기간 2020s (2010~2039), 2050s (2040~2069), 2080s (2070~2099)에 대해서 분석하였다. RCP 4.5 및 8.5 시나리오별 기후변화를 분석한 결과 Historical 대비 미래 강수량의 경우 RCP 4.5는 4.7~13.6% 증가하였고, RCP 8.5는 1.3~13.0% 증가했다. 미래 기온의 경우 Historical 대비 RCP 4.5는 0.9~2.4°C, RCP 8.5는 1.2~5.3°C 상승했다.
3)기후변화에 따른 유출특성 변화를 분석한 결과 RCP 4.5 2080s에서 총 유출량과 지하수함양량 및 증발산량은 Historical 대비 각각 14.6%, 10.2%, 17.2% 증가한 것으로 분석되었다. 특히 Historical 대비 지하수함양량의 변화를 공간적으로 분석한 결과 대청댐 상류 보청천 유역의 경우 모든 기후변화 시나리오 및 평가기간에서 지하수함양량이 감소하는 것으로 나타났다. 특히 가장 크게 감소한 RCP 8.5 시나리오의 2020s를 공간적으로 살펴보면 대청댐 상류 뿐 만아니라 금강 하류 논산천 하류 주변의 농경지를 포함하는 소유역에서 지하수함양량이 감소하는 것으로 나타났다.
4)최근 10년(2005~2015)동안의 국가지하수조사연보자료를 수집하여 자연로그의 식으로 회귀식을 산정하고 미래 기후변화 평가기간에 해당하는 지하수이용량을 추정한 결과 2015년 대비 각각 4.6%, 7.7%, 9.6% 증가한 것으로 분석되었으나 SWAT 모의 결과는 2020s에 지하수함양량은 RCP 4.5와 RCP 8.5 시나리오 모두 Historical 대비 -0.4%, -6.2% 감소한 것으로 분석되어 추후 지하수이용량의 증가 추세를 고려할 경우 지하수량은 더 감소할 것으로 판단된다.
5)기후변화에 따른 지하수위 변동량을 정리한 결과 봄(3~5월), 여름(6~8월)의 지하수위는 감소한 반면 가을(9~11월)은 강수량의 증가로 상승하였다. 연평균 지하수위 변동은 RCP 4.5와 8.5 각각 -9.0~8.0 cm, -13.0~4.0 cm 의 범위로 분석되었으며, 지하수 이용량 고려 시 지하수위의 변동은 RCP 4.5와 8.5 각각 -9.2 ~7.8 cm, -13.3~3.8 cm으로 나타났다.
지하수위의 변화는 대수층의 지하수량의 증감 파악에 있어 매우 중요한 요소이며, 만일 지하수위가 하강하면 대수층 내의 지하수량이 고갈되고 있는 중이고, 상승하면 대수층 내로 지하수가 유입되어 지하수량이 증가하고 있음을 의미한다(Cho et al., 2012). 본 연구결과 중 지하수 이용량 고려 시 모든 시나리오에서의 지하수위의 변동은 -13.0~8.0 cm로 크게 체감되지 않는 수치이나 가뭄 등의 극심한 물 부족 발생할 경우 기존의 사용량을 초과하여 양수하는 문제가 발생하여 지하수위는 급격히 하강되고, 이는 유역의 수자원관리에 큰 문제를 야기할 수 있다. 또한 Fig. 7과 Table 6에서와 같이 동일한 대권역 내에 위치한 유역일 지라도 해당 소유역의 토지이용상태 혹은 토질에 의해 지하수 함양량 및 지하수위의 변동이 상이 할 수 있다. 따라서 중소유역 별 지하수 관리 지침을 수립하여 향후 기후변화 대응 수자원 확보 측면에서 지하수 함양량 및 지하수위에 대한 연구는 체계적으로 이뤄져야 할 것으로 생각된다.
본 연구의 결과는 지하수위 관측정이 설치된 지점의 HRU의 계산결과를 추출하고 매개변수 값에 변동량을 더해 주어 지하수위를 추정했기 때문에 일부 지하수 하강특성을 완벽히 재현하지 못한 한계점이 있으며 적용한 기후변화 시나리오 및 수문모의 결과의 불확실성을 내포하고 있다. 그러나 기후변화에 따른 지하수 및 지하수위의 거동을 파악하고 예측하는 방법을 유역의 표준단위 유역별, 일/년도/계절별 통계분석결과를 제시함으로써 기후변화 적응을 위한 수자원의 효율적 관리에 적극적으로 활용될 것으로 판단된다.
