농업 밀집지역인 해안분지에서 지하수와 하천수의 수리화학적 특성 평가

1. 서 론

최근 전 세계는 극심한 가뭄 및 홍수와 같은 기상이변과 함께 인간 활동으로 인한 수질오염으로 안정적인 수자원 확보에 어려움을 겪고 있다. 이에 많은 나라들은 지속적인 수자원 관리 및 이용하기 위해 물이용 효율성, 지하수 및 하천수의 통합적 관리를 위한 노력을 하고 있다. 이 중 지하수는 많은 나라에서 농업용수 및 생활용수를 공급하는 중요한 수자원 역할을 한다. 현재 전 세계 인구의 약 1/3은 지하수를 음용수로 의존하고 있다(UNEP, 1999). 지하수 자원의 수질 조성 연구는 수자원의 안정적 확보를 위한 기본 자료 제공 및 관리방안의 기본 틀을 마련한다는 측면에서 중요하다.

과거 연구자들은 지하수와 하천수가 각 독립된 수체로 생각하여 연구를 수행하였다(Kalbus et al., 2006; Yang et al., 2014). 그러나 많은 경우, 지하수와 하천수는 수리학적으로 서로 매우 밀접한 상호관계를 가진다. 이를 위해 연구자들은 두 수체의 상호작용을 평가하기 위한 다양한 방법들을 제시하였으며, 지하수와 하천수의 통합관리의 중요성을 보고하였다(Winter et al., 1998; Brodie et al., 2007; Jeon et al., 2015).

지하수 및 하천수의 수리화학적 특성을 비교·분석하는 연구는 지하수의 화학적 성질을 지배하는 지질학적 과정을 이해하는데 도움이 된다. 지하수 및 하천수가 생활/농업용수로서의 수질 적합성을 결정하는데 중요한 척도가 된다. 정기적인 수질평가는 그 지역의 주민과 정책자에게 수질에 대한 많은 정보를 알려주는 역할을 한다. 산소·수소동위원소 조성 분석은 지하수 함양지역 및 시간 추정, 지하수 흐름 경로 및 지역의 함양조건, 증발의 정도, 그리고 두수체간의 상호연결성을 이해하는데 중요한 정보를 제공한다(Kanduč et al., 2014).

양구군 해안분지는 농업밀집 지역으로서 다른 농촌지역에 비해 매년 대량의 화학비료 및 퇴비를 사용하고 있어 지하수 내 나트륨 및 질산성 질소의 오염 가능성이 대두된 바 있다(Kim et al., 2015). 농업지역에 해당되는 낙동강 유역, 제주시, 금강권역 등 여러 지역에서 충적대수층이 질산성 질소로 오염된 연구사례가 있다(Choi et al., 2002; Min et al., 2003; Chae et al., 2009). 이 연구결과들에 따르면 농업활동으로 인해 대부분 지하수 관정의 질산성 질소가 먹는물 수질기준(10 mg/L)을 초과한 것으로 나타났다. 이는 우리나라 농업지역 지하수의 질산성질소에 대한 영향 평가 및 수리적 특성에 대한 이해가 필요함을 시사한다.

양구군 해안분지의 고지대는 농업활동이 활발하게 이루어지는 지역으로, 특히 여름철에 기온역전 현상으로 냉기가 형성되어 대부분 고랭지밭에서 무, 배추, 인삼이 주로 생산된다. 이에 대규모 밭으로 조성되면서 폭우 시 지표유출로 인한 토양유실의 문제가 발생하고 있다. 또한 강우강도에 따라 수일 또는 수주가 탁수현상이 지속되고 있는 상황이다. 이러한 탁수 문제는 연구지역의 수자원 이용, 저수지 수질관리, 수리 생태계에 부정적인 영향을 미친다. 또한 지하수와 하천수의 수질은 농작물의 품질과 생산량 그리고 인간 건강에 상당한 영향을 줄 수 있어 이에 대한 체계적인 분석이 중요하다(Hoek et al., 2001). 현재까지 수화학 측면에서 연구지역의 전반에 걸친 장기 모니터링의 연구는 없었다. 이 연구 목적은 연구지역에서 다년간 수집된 자료를 통해 두 수체의 수리화학적 특성을 다각적으로 평가하고 이를 통해 연구지역의 수리적 특성을 이해하고자 하는데 있다. 이 연구결과가 연구지역의 계절적 변화에 따른 지하수와 하천수의 물수지 상관관계를 평가하는데 있어 기초 자료로 활용될 것으로 기대된다.

2. 연구지역 및 방법

2.1 연구지역

연구지역은 남한과 북한이 대치중인 군사분계선에 인접한 대한민국 강원도 양구군 해안면(38.239~38.329°N, 128.083~128.173°E)에 위치하고 소양강 유역에 해당된다(그림 1a). 2011년 기준 해안분지의 전체면적은 61.75 km²이며 가칠봉(1,242 m), 대우산(807 m), 도솔산(1,147 m) 등 해발고도 1,000 m 내외의 산으로 둘러싸인 폐쇄형 분지이다. 연구지역의 해발고도 범위는 339 m부터 1,321 m 사이를 보인다. 지형적 특성으로 해당 면적에 비해 고도차가 커 전경을 한눈에 볼 수 있고 다른 지역과 격리되어 있어 주변 능선의 경계가 곧 행정구역이다.

Fig. 1.

(a) Location of the Haean basin and water sample points, (b) land use of the study area.

연구지역의 일반지질은 외륜 및 능선인 선캠브리아기 변성암 복합체와 내부는 쥬라기 흑운모 화강암로 구성되며 그 외에 지천을 중심으로 제4기 충적층으로 형성된다(Won et al., 1987, 1988). 이 분지의 생성원인은 기존 선캠브리아기 변성암 복합체 지역에 흑운모 화강암이 분지 내부로 관입한 후 차별침식작용으로 현재의 펀치볼 모양의 분지가 형성된 것으로 추측된다(Won et al., 1987, 1988; Kwon et al., 1990). 분지 내부는 풍화에 약한 흑운모 화강암이 풍화토(saprolite)로 산출되고 분지 외륜는 변성암복합체며 편암, 편마암, 규암 혹은 석류석으로 분포된다.

분지 내 하천은 하천폭 약 1~2 m의 지천들로 구성된 수지상 또는 방사상의 수계이며 주요 하천은 폭 10~15 m인 만대천, 도솔천 그리고 성황천이다(그림 1a). 여러 지천들이 만대천으로 한곳에 수렴하여 연구지역의 동쪽 당물골로 배출되는 단순한 수리학적 시스템이다. 하천의 경사는 상대적으로 가파른 산지에서 상류(20°)로 시작하여 경사가 완만한 분지 내 하류(5°)까지 점차 감소한다. 이 연구지역의 평균 경사는 11°이다(Lee et al., 2011).

2015년 양구군 통계자료에 따르면 해안분지에 신고 지하수 관정은 310개이지만 미신고된 관정까지 고려하면 약 600개 이상으로 추정된다. 이에 연구지역의 지하수 개발밀도는 4.56개/km²로 추정된다. 대부분 지하수 관정은 농업용수 목적으로 주로 사용되며 장기간 상당량의 관개수가 필요한 논에 밀집되어 있다.

연구지역의 기후는 여름엔 덥고 습하며 겨울에는 춥고 건조한 전형적인 동아시아 몬순기후 영향을 강하게 받는다. 1999년부터 2009년까지 연평균 기온은 약 8.5℃이며 최저기온인 1월 -26.9℃과 최고기온인 8월 33.4℃의 범위를 보이는 등 4계절이 뚜렷하게 나타난다. 연강수량은 1,577 mm이며 이 중 70%가 장마로 인해 우기인 6월에서 8월 사이에 내린다(Kettering et al., 2012). 이 지역의 농업활동은 4월에서 10월 사이에 거의 90%가 이루어진다.

2010년 기준으로 해안분지의 토지 이용률은 해발고도에 따라 그 분포가 다르게 나타난다. 해발고도 650 m 이상은 가장 큰 면적을 차지하는 산림(58%)이며 650 m 이하 고랭지밭(27%), 논(9%), 기타(5%), 과수원(0.5%) 순으로 나타난다(그림 1b). 2002년부터 2014년까지 전체 토지 이용면적(6,150 ha) 변화를 보면 큰 면적을 차지하는 산림은 3,580~3,604 ha, 밭은 1,683~1,725 ha 그리고 논은 529~537 ha로 거의 변화가 없었다.

한편 상류 지역인 양구 해안분지에서 발생한 탁수는 하류에 있는 여러 지역에 탁수 피해를 입히고 소양강댐의 수질 영향을 주고 있다(Jo et al., 2010; Lee et al., 2011). 2006년부터 2012년까지 양구군청은 탁수저감을 시키기 위해 침사지 8개를 건설했으나 폭우 시 큰 저감 효과를 보지 못하였다. 탁수의 장기화는 수리 생태계 파괴 및 하천수 수질 저하에 영향을 끼칠 수 있어 체계적인 관리가 필요하다.

3. 연구결과 및 토의

3.1 일반 수리화학적 특성

표 1은 연구지역의 지하수와 하천수의 수리화학적 특성에 대한 기초통계 분석결과를 보여준다. 두수체의 산소·수소동위원소 조성 및 질산염 농도의 영향인자로 추정되는 논물 13개는 표 1에 함께 나타냈으며 그 내용은 산소·수소동위원소 및 질산염 오염의 평가 파트에서 다루겠다. 각 지하수와 하천수의 칼륨농도의 범위는 0.2~6.9 mg/L, 0~13.3 mg/L이며 평균은 1.5와 1.8 mg/L로 서로 비슷한 값을 보인다. 지하수 내 칼륨은 양이온 교환작용에 의해 지하수는 하천수보다 칼륨농도가 보통 낮게 나타나는 것으로 알려져 있다(Griffioen, 2001).

Table 1.

Descriptive statistics of field parameters and ionic compositions for groundwater, stream water and ponded water in the studied period.

분석결과, 연구지역 지하수와 하천수의 나트륨 농도 범위는 0.7~25.5 mg/L와 2.5~44.4 mg/L이며 평균 9.4와 6.9 mg/L로 보인다. Ca와 Mg는 보통 돌로 마이트와 방해석 같은 탄산염 용해로 생성되며 연구지역에서는 지하수와 하천수의 Ca의 농도범위는 각각 4.1~274 mg/L, 3.4~42.2 mg/L이며 평균은 30.9와 15.8 mg/L이다. Mg의 농도범위는 지하수와 하천수 각각 0.1~20.8 mg/L와 0.5~10.3 mg/L이며 평균은 3.7과 2.7 mg/L이다. 전반적으로 지하수와 하천수의 주요 양이온은 Ca와 Na 그리고 주요 음이온은 HCO₃와 NO₃가 우세하게 나타난다.

지하수와 하천수의 표준편차는 각각 1,008과 77.1로 EC에서 가장 크게 나타난 반면 산소·수소동위원소에서 각각 4.7과 5.7로 가장 작게 나타난다. 표본표준편차와 표본평균의 비인 지하수와 하천수의 변동계수는 각각 EC (3.17), Cl (1.29)가 가장 크다. 특히 지하수의 Ca와 HCO₃는 하천수에 비해 각각 6.6과 9.4배로 높게 나타났다. 이는 물과 암석 반응으로 용해된 광물이 지하수에 영향을 준 것으로 분석된다.

지하수와 하천수의 현장수질, 주요 양이온, 음이온에 대한 차이와 관계를 평가하기 위해 상관분석을 수행하였다(표 2). 지하수는 Ca와 HCO₃, Mg, Na 그리고 Mg와 HCO₃, Na 그리고 Na와 HCO₃ 순으로 높은 상관관계(r>0.78)를 보인다. 특히 지하수 내 높은 상관성을 보인 Ca, Mg은 몬모릴로나이트, 일라이트 등 점토광물로 기인된 것으로 분석된다(Garrels, 1976). 하천수는 Ca와 Mg, HCO₃ 그리고 EC와 Mg, Ca 그리고 Mg, SO₄ 그리고 Na과 Cl 순으로 높은 상관관계(r>0.70)를 보인다. 하천수에서 EC에 대한 모든 양이온, 음이온의 상관계수가 0.48 이상으로 높게 나타났다. 지하수에 비해 하천수에서 Na와 Cl의 상관계수는 0.90으로 강한 상관성을 보이는데 유속이 빠른 하천수가 농업활동 의한 비료사용 등 인위적인 유입의 영향을 더욱 뚜렷하게 받은 것으로 추정된다.

Table 2.

Pearson correlation coefficients of analyzed chemical parameters of groundwater and stream water.

그림 2는 지하수와 하천수의 수리화학적 자료를 Piper (1944) diagram 을 통해 유형별로 비교하기 쉽게 도시한 결과를 보여준다. 지하수는 114개 시료 중 Ca-HCO₃-Cl (47.0%)와 Ca-HCO₃ (36.5%) 유형에 도시되었다. 그 외에 Ca-Na-HCO₃ (12.2%), Ca-Cl (0.9%), Ca-Na-HCO₃-Cl (2.7%) 그리고 Ca-Na-Cl (0.9%) 유형이 일부 나타난다. 한편 하천수는 주로 Ca-HCO₃-Cl (83.5%) 유형에 도시된다. 그밖에 Ca-Cl (5.2%), Na-Cl, Ca-HCO₃ (2.7%), Ca-Na-HCO₃-Cl (2.7%) 그리고 Ca-Na-Cl (0.9%) 유형 순으로 나타난다. 농업지역에서 Cl은 화학비료 사용으로 인해 질산염과 함께 나타나는 주요 음이온이다. 이에 지하수와 하천수 모두 지표 기원 오염물질인 Cl에 지속적인 영향을 받는다(Ki et al., 2013, Kim et al., 2015). 한편 논물의 유형은 Ca-HCO₃-Cl (61.5%), Ca-Na-HCO₃-Cl (30.8%) 그리고 Ca-Na-Cl (7.3%) 순으로 분포한다.

Fig. 2.

Water types for groundwater and stream water samples in the studied area for 2011~2014.

9회 샘플링 자료를 시기별로 비교·분석한 결과, 지하수와 하천수의 물 유형의 경우 시간에 따라 변화를 보이지 않지만 지하수와 하천수의 도시되는 분포범위는 다르게 나타났다. 지하수는 하천수에 비해 Ca-HCO₃-Cl, Ca-HCO₃ 유형을 걸치면서 넒은 분포를 보이는 반면 하천수는 대부분 Ca-HCO₃-Cl 유형에 집중 분포한다. 이것은 하천수의 유속이 빨라 활발한 수리화학적 반응을 겪을 시간이 짧은 것과 관련이 있을 것으로 추정된다. 지하수의 경우, 상대적으로 넒은 화학조성 분포를 보이는데 이는 지하수는 충분한 화학적 반응과 토지이용별 직접적인 영향을 받기 때문이다.

연구기간에 채취된 지하수와 하천수의 주요 양이온 및 음이온의 구성변화(meq/L)를 Schoeller (1977) diagram을 이용하여 시기별로 도시하였다(그림 3). 전반적으로 지하수와 하천수는 양이온의 경우 Ca>Na>Mg>K, 음이온의 경우 HCO₃>NO₃>Cl>SO₄ 순으로 비슷한 화학조성 비율을 보였으며 이는 두 수체가 수리연결성이 좋다는 것을 시사한다. 선행연구인 Yun et al. (2015)에서는 강우 이벤트 시 기저유출이 단시간 내 발생하여 타 지역에 비해 빠른 수문순환을 보여주는 연구사례가 있다. 지점별 지하수와 하천수의 각 이온들의 변화 정도는 상이하게 나타난다. 지하수의 이온 변화는 상대적으로 불규칙하고 광범위하게 변화를 보이는 반면, 하천수는 좁은 범위 내 일정하게 나타난다. 지하수의 경우 특히, NO₃ 이온의 변화폭이 크게 나타났는데, 이는 작물재배 및 비료사용 등의 인위적인 영향을 받은 것으로 추정된다.

Fig. 3.

Schoeller diagram for groundwater and stream water samples from 2011 to 2014.

지하수와 하천수에 나타나는 강수의 영향, 물-암석 교환, 증발 등의 수리화학적 과정을 이해하기 위해 Gibbs (1970) diagram을 이용하였다. 우리는 TDS에 따른 Na/(Na+Ca)과 Cl/(Cl+HCO₃)의 비율을 각각 계산한 후 (1) 증발 우세 (2) 강수 우세 (3) 암석우세 총 3개 영역으로 도시하였다(그림 4). 지하수와 하천수는 Gibbs diagram에 중간 영역인 암석 우세에 대부분 도시되었으며, 이러한 결과는 전형적 분지형태인 연구지역의 수리지질학적 환경과 관련이 있다. 이 연구지역은 폐쇄형 지하수 함양지역이며 물-암석 교환작용에 의해 지하수 및 하천수의 수화학이 형성된 것으로 보인다. 연구지역에서 지하수는 자유면 대수층 및 피압대수층 상관없이 물-암석 교환작용에 영향을 받는다. 하천수의 일부 시료는 강수 우세를 보이는데 이것은 우기 때 강우영향을 받은 것으로 분석된다.

Fig. 4.

Gibbs diagrams showing the mechanism of groundwater and stream water evolution. (a) TDS versus Na/(Na+Ca), (b) TDS versus Cl/(Cl+ HCO3).

3.2 산소·수소동위원소의 시·공간적 분포 특성

2011년부터 2014년까지 지하수와 하천수의 산소·수소동위원소 범위는 각각 -10.7~-6.5‰, -79.9~-51. 9‰와 -10.6~0.5‰, -80.0~-24.6‰이다. 하천수의 산소·수소동위원소의 변화 범위는 지하수 산소·수소동위원소보다 크다. 추가로 2012년 6월과 2013년 6월에 채취된 13개 논물의 산소·수소동위원소는 각 -8.7~4.8‰와 -67.1~-9.4‰ 범위를 보이며 각 평균 -4.1‰와 -47.0‰이다.

지역순환수선과 함께 지하수 및 하천수를 비교하는 것은 지역적 수리학 연구에서 수원을 결정하는데 유용한 지표로 될 수 있다(Clark and Frjtz, 1997). 이 연구지역의 두 수체는 Park et al. (2006)의 지역순환수선(δD=8.06×δ¹⁸O+12.5, r²=0.95; Local Meteoric Water Line) 그리고 Craig (1961)의 지구순환수선(D=8×δ¹⁸O+10; Global Meteoric Water Line) 부근에 도시된다(그림 5).

Fig. 5.

δ18O versus δD values of groundwater and stream water shown relative to local meteoric water line (LMWL; δD=8.06δ18O+12.5, from Park et al., 2006) and the global meteoric water line (GMWL; δD=8.1δ18O+9.8, from Craig, 1961) in the Haean Basin.

그림 5e, 5f, 5g에서 두 수체가 GMWL 오른쪽 하단에 도시된 것은 강한 증발 영향이 반영된 것이다. 시기별 지하수에 대한 기울기 범위는 5.89~8.54로 보이며 전체 평균 기울기(6.15)는 LMWL (8.06) 기울기보다 작다. 이는 연구지역의 지하수는 주로 강수 기원이며 지하침투 전에 다양한 증발작용을 거쳤음을 의미한다(Friedman et al., 1962; Paul and Wanielista, 2000; Cui and Li, 2014; Wang et al., 2015). Gonfiantini(1986)에 따르면 하천이 강한 증발작용을 겪었다면 하천수의 기울기는 3.9~5.2 범위로 나타난다. 그림 5a, 5f, 5g, 5h에서 하천수의 기울기는 각각 5.04, 3.37, 3.05, 4.92로 강한 증발영향을 받았으며 이것은 일사량이 강했던 계절들과 일치한다.

두 수체의 전반적인 경향을 살펴보면 지하수의 산소·수소동위원소는 시료수 및 시기와 상관없이 하천수에 비해 넓은 분포를 보인다. 이는 증발영향으로 무거워진 논물의 산소·수소동위원소가 지하수와 하천수에 영향을 준 것으로 분석된다(Fernández-Chacón et al., 2010). 일부 시료를 제외한 두 수체의 산소·수소동위원소는 매우 유사하게 나타난다. 이것은 연구지역의 수리학적 수리연결성이 좋으며 빠른 수문순환이 일어나고 있음을 의미한다(Akomako et al., 2011; Lee et al., 2013).

그림 6은 연구지역의 고도별 산소동위원소 조성 변화를 시기별로 보여준다. 이상값을 제외한 고도에 따른 지하수 및 하천수의 산소동위원소 조성변화는 그림 6g, h, i에서 상대적으로 잘 나타난다. 이 중 회귀선이 잘 표현된 산소동위원소는 각 δ¹⁸O=0.009EL-5.46(r²=0.41), δ¹⁸O=0.004EL-7.56 (r²=0.44)로 나타나지만 고도효과는 뚜렷하게 보이진 않는다. 그러나 연구지역의 고도차(170 m)가 작음에도 불구하고 고도가 증가할수록 산소동위원소가 감소하는 경향성을 보여 그 의미가 있다(Lee et al., 2013).

Fig. 6.

Relationships of δ18O with topographic elevations (EL) for groundwater and stream water in the study area. The solid and dotted line represent the regression line for the groundwater and stream water, respectively.

또한 시기별 고도 증가에 대한 지하수의 산소동위원소는 하천수의 산소동위원소보다 수직으로 넓게 분포된 경향성을 보인다. 반면 하천수는 지형학적 특성으로 한 곳에 집중 분포하고 있다. 연구지역 저지대(430~460 m)의 토지이용은 하천을 중심으로 논이 분포하고 있다(그림 1b). 마찬가지로 강한 증발작용 인해 무거워진 논물의 산소·수소동위원소는 지하로 침투하여 지하수내 동위원소에 영향을 미친 것으로 풀이된다(Cui and Li, 2014). 더 나아가 지하수의 수직적 분포는 각 논에 관개수 공급되는 시점 및 증발 시간이 달라 비롯된 것으로 분석된다. 또한 논의 토양 침투속도는 논수위 변화에 따라 다르게 나타난다(Sohn et al., 2004). 이는 충적대수층에 상당한 지하수 함양을 기여한다(Liu et al., 2001; Chen et al., 2002; Huang et al., 2003).

3.3 질산염 오염 가능성 및 관개수의 적합성 평가

국내 먹는물 수질기준으로 질산염 농도는 44.3 mg/L(질산성 질소 농도는 10 mg/L)로 정해져 있다. 연구기간 동안 채취된 지하수내 질산염 농도 범위는 0.1~88.2 mg/L이며 평균 19.2 mg/L로 나타난다(표 1). 연구지역에서 질산성 질소 오염이 다수 발생할 것으로 예상했지만 먹는물 기준에 초과한 지하수 시료는 144개 중 16개(11%)만 나타났다. 하지만 인위적 오염물질의 유입가능성(질산성 질소: 3 mg/L 초과)이 있는 시료는 57개(39.6%)로 높게 나타났다. 한편 2012년 6월과 2013년 6월에 논물의 질산염 농도범위는 0.1~19.3 mg/L이며 평균 6.5 mg/L로 지하수에 비해 낮게 나타났다(표 1).

연구지역 질산염 농도의 시공간적 분포 현황을 파악하기 위해 열지도(heat map) 분석을 수행하였다(DeBoer, 2015). GIS의 열지도를 이용한 지하수 내 질산염 농도의 분포 변화는 그림 7에 시기별로 도시하였다. 분석결과, 연구지역의 질산염 오염 분포 경향의 경우, 시기별로 뚜렷한 차이를 보이지는 않았다. 공간적 분포의 경우, 만대천, 도솔천, 성황천 주변 저지대에서 질산염 농도가 상대적으로 높게 나타났다(그림 7a, 7e, 7g, 7h, 7i). 시료수(n<6)가 적은 시기를 제외한 농한기(n=15) 지하수내 평균 질산염 농도는 6.0 mg/L로 농번기(n>17)의 지하수 질산염 농도보다 약 2~4배가 낮다.

Fig. 7.

Heat map of NO3 concentration in groundwater from 2011 to 2014.

연구지역에서 지하수의 질산염 농도가 낮은 이유는 논의 특성과 관련이 있는 것으로 추정된다. 일반적으로 논토양의 하층은 산소가 부족하여 미생물에 의한 탈질이 강하게 나타나게 된다(Kim et al., 2015). 이는 논에 사용된 질소 비료의 상당량이 탈질작용으로 없어짐을 시사한다. 다만, 이 연구에서 분석한 질산염 데이터의 경우, 시기별 시료 채취된 지점이 상이하고 일부 지점만 갖고 분석된 것으로 전체 연구지역의 질산염 농도 분포를 평가하기엔 다소 한계가 있다.

연구지역에서 지하수와 하천수의 SAR 값은 각 0.1~1.0 meq/L와 0.2~2.4 meq/L의 범위로 매우 낮게 보였다. 연구지역의 SAR에 대한 지하수와 하천수의 세부적 분류는 표 3에 나타냈다. 지하수와 하천수의 시료들은 114개 시료 모두 SAR<10의 값을 보여주었다. 즉 관개수로서의 지하수와 하천수는 나트륨 교환가능 위험이 없는 Excellent quality 등급으로 나타났다. EC 역시 관개수의 염분 위험을 판별할 수 있는 좋은 도구다. EC에 관한 관개수로서 지하수와 하천수의 수질 분류를 살펴보면 지하수는 Excellent에서 88개(77.2%) Good에서 22개(1.9%), Fair에서 3개(2.6%) 그리고 Poor에서 1개(0.9%)로 나타났다(표 3). 하천수는 Excellent에서 101개(88.6%) 그리고 Good에서 13개(11.4%)로 나타났다.

Wilcox diagram에서 EC에 대한 지하수와 하천수의 산출된 초과 나트륨 함량은 그림 8a에 도시하였다. 지하수와 하천수 모두가 Excellent to good 범주에 속하였다. 그 외에 지하수 3개 시료는 Good to permissible 범주에 나타났고 EC가 매우 높은 1개 시료는 Unsuitable로 판명되었다. 이 지점은 식물의 성장에 지장을 줄 수 있는 우려가 있다(Todd, 1980).

Fig. 8.

(a) Plot of percent sodium vs. electrical conductivity (from Wilcox, 1955), (b) diagram for irrigation waters classification (from USSL, 1954).

또한 연구지역에서 지하수와 하천수는 EC와 함께 US Salinity diagram을 이용하여 나타내었다(그림 8b). 대부분 지하수와 하천수의 시료들은 C1-S1과 C2-S1 범주에 속한다. 이것은 낮음에서 중간 염분 그리고 낮은 나트륨 지하수 및 하천수를 나타낸다. 이는 연구지역의 대부분의 지하수 및 하천수는 나트륨 교환 확률이 낮으며, 따라서 대부분의 토양유형에 관개수로 적합함을 의미한다. 그러나 4개의 지하수 시료는 높은 EC (<750 μS/cm)를 보여 교환 가능한 나트륨 위험의 원인이 될 수 있다. 채취된 2지점에서 꾸준히 높은 EC가 나타나 향후 지속적인 모니터링이 필요하다.

Table 3.

Suitability of groundwater and stream water for irrigation based on EC, SAR and Na% classifications.

4. 결 론

농업 밀접지역인 양구 해안분지는 지형학적 및 수리지질학적으로 단순하여 연구가치가 크다. 특히 연구지역은 최근 3년간 기후변화로 인해 농번기에 극심한 가뭄을 겪고 있으며 화학비료 및 퇴비사용으로 수질악화가 우려되고 있다. 이에 본 연구는 두 수체의 수리화학적 특성을 다각적으로 평가하였다.

연구지역의 지하수와 하천수의 일반 수리화학적 특징은 양이온 경우 Ca와 Na 그리고 음이온 경우 HCO₃와 NO₃가 두 수체 모두 우세하게 나타난다. 하천수의 Na와 Cl의 상관계수는 지하수에 비해 강한 상관성을 보였다. 이는 외부에 노출된 하천이 인위적 영향을 받은 것으로 판단된다. 지하수는 Ca-HCO₃-Cl 및 Ca-HCO₃ 유형으로 넓은 분포를 보인 반면 하천수는 Ca-HCO₃-Cl 유형에 집중분포를 하였다. 두 수체는 지표 기원 오염물질인 Cl에 지속적으로 영향을 받은 것으로 분석된다. 또한 지하수와 하천수는 양이온 경우 Ca>Na>Mg>K, 음이온 경우 HCO₃>NO₃>Cl>SO₄ 순으로 유사한 화학적 구성을 보인다. 이를 통해 두 수체간 수리연결성이 상대적으로 좋음을 확인 할 수 있다. 전체적인 연구지역의 지하수 및 하천수의 수화학은 빠른 유동과 함께 물-암석 교환작용을 받아 형성된 것으로 판단된다.

산소·수소동위원소를 통해 연구지역의 지하수는 강수 기원이며 지하 침투 전에 다양한 증발작용을 겪은 것으로 확인되었다. 특히 일부 하천수는 강한 증발의 영향을 받았던 것으로 보인다. 두 수체의 산소·수소동위원소 조성은 유사하게 나타났으며, 이는 수문순환이 빠르고 수리연결성이 좋다는 주장을 뒷받침해준다. 한편 연구지역의 전반적으로 작은 고도차이임에도 불구하고 고도가 증가할수록 산소동위원소가 감소하는 경향을 보인다. 특히 지하수의 광범위한 산소동위원소 분포는 증발작용으로 무거워진 논물의 영향을 받은 것으로 보인다.

시기별 질산염의 농도 분포는 만대천, 도솔천, 성황천 부근에 질산염 농도가 다른 지점에 비해 높게 나타나는 경향을 보였다. 이는 여러 지천의 질산염이 주요 하천으로 모여 높게 나타난 것이다. 한편 지하수의 먹는물 수질 분석 결과 일부 시료만 먹는물 수질 기준치를 초과하는 것으로 나타났으며 이는 탈질의 영향을 받은 것으로 판단된다. SAR, EC 그리고 Na% 지표에 따라 연구지역의 지하수 및 하천수는 관개수로서 매우 적합한 것으로 평가되었다.

해당 연구지역은 농작물의 생산성을 높이기 위해 가축 퇴비뿐만 아니라 다양한 무기질 화학비료와 농약이 쓰이고 있다. 연구지역의 높은 지하수-하천수 간 수리 연결성을 고려하면 이는 해당 지역이 수질오염에 매우 취약할 수 있음을 시사한다. 따라서 향후 지속적인 모니터링과 함께 중금속 및 NH₄, NO₂, PO₄ 이온을 추가 분석하여 종합적 연구수행이 필요하다.

Acknowledgments

이 논문은 2017년도 정부(미래창조과학부)의 재원으로 한국연구재단의 지원을 받아 수행된 연구입니다(No. NRF-2015R1A4A1041105). 또한 2016년 강원대학교 대학회계 학술연구조성비로 연구하였습니다(관리번호-520160395). 아울러 심사과정에서 발전적인 비판과 조언을 해주신 익명의 심사자님들께 깊은 감사를 드립니다.

REFERENCES

• Akomako, D., Osae, S., Akiti, T.T., Faye, S., and Maloszewski, P., (2011), Geochemical and isotopic studies of groundwater conditions in the Densu River Basin of Ghana, Environment Earth Science, 62, p1071-1084.
• Alam, F., (2014), Evaluation of hydrogeochemical parameters of groundwater for suitability of domestic and irrigational purposes: a case study from central Ganga Plain, India, Arabian Journal of Geosciences, 7, p4121-4131. [https://doi.org/10.1007/s12517-013-1055-6]
• Brodie, R., Sundaram, B., Tottenham, R., Hostetler, S., and Ransley, T., (2007), An overview of tools for assessing groundwater-surface water connectivity, Bureau of Rural Sciences, Canberra, Australia, p131.
• Clark, I.D., and Frjtz, P., (1997), Environmental isotopes in hydrogeology, Lewis Publishers, New York.
• Chae, G.T., Yun, S.T., Mayer, B., Choi, B.Y., Kwon, J.S., and Yu, S.Y., (2009), Hydrochemical and stable isotopic assessment of nitrate contamination in an alluvial aquifer underneath a riverside agricultural field, Agricultural Water Management, 96, p1819-1827. [https://doi.org/10.1016/j.agwat.2009.08.001]
• Chen, S.K., Liu, C.W., and Huang, H.C., (2002), Analysis of water movement in paddy rice fields (II) simulation studies, Journal of Hydrology, 268, p259-271. [https://doi.org/10.1016/s0022-1694(02)00180-4]
• Choi, W.J., Han, G.H., Ro, H.M., Yun, S.H., and Lee, S.M., (2002), Evaluation of nitrate contamination sources of unconfined groundwater in the North Han River basin of Korea using nitrogen isotope ratios, Geosciences Journal, 6, p47-55. [https://doi.org/10.1007/bf02911335]
• Craig, H., (1961), Isotopic variations in meteoric waters, Science, 133, p1702-1703. [https://doi.org/10.1126/science.133.3465.1702]
• Cui, B.L., and Li, X.Y., (2014), Characteristics of stable isotope hydrochemistry of the groundwater around Qinghai Lake, NE Qinghai-Tibet Plateau, China, Environment Earth Science, 71, p1159-1167.
• DeBoer, M., (2015), Understanding the heat map, Cartogrphic Perspectives, 80, p39-43.
• Fernández-Chacón, F., Benavente, J., Rubio-Campos, J.C., Kohfahl, C., Jiménez, J., Meyer, H., Hubberten, H., and Pekdeger, A., (2010), Isotopic composition (d¹⁸O and D) of precipitation and groundwater in a semi-arid, mountainous area (Guadiana Menor basin, Southeast Spain), Hydrology Process, 24, p1343-1356.
• Friedman, I., Machta, L., and Soller, R., (1962), Water vapour exchange between a water droplet and its environment, Journal of Geophysical Research, 67, p2761-2766. [https://doi.org/10.1029/jz067i007p02761]
• Garrels, R.M., (1976), A survey of low temperature water mineral relations, In: Interpretation of environmental isotope and hydrogeochemical data in groundwater hydrology, vienna, International Atomic Energy Agency, p65-84.
• Gibbs, R.J., (1970), Mechanism controlling world water chemistry, Science, 170, p1088-1090. [https://doi.org/10.1126/science.170.3962.1088]
• Gonfiantini, R., (1986), Environmental isotopes in lake studies, Elsevier, New York, p133. [https://doi.org/10.1016/b978-0-444-42225-5.50008-5]
• Griffioen, J., (2001), Potassium adsorption ratios as an indicator for the fate of agricultural potassium in groundwater, Journal of Hydrology, 254(1-4), p244-254. [https://doi.org/10.1016/s0022-1694(01)00503-0]
• Hoek, W., Konradsen, F., Ensink, J.H.J., Mudasser, M., and Jensen, P.K., (2001), Irrigation water as a source of drinking water: is safe use possible?, Tropical Medicine & International Health, 6(1), p46-54.
• Huang, H.C., Liu, C.W., Chen, S.K., and Chen, J.S., (2003), Analysis percolation and seepage through paddy bunds, Journal of Hydrology, 284, p13-25.
• Jeon, W.H., Lee, J.Y., Cheong, W.Y., Park, Y.H., Oh, S.H., Eum, D.H., and Park, J.Y., (2015), A multi-method approach revealing the groundwater-stream water interaction in the Inbuk stream, Korea, Geosciences Journal, 19(2), p325-340. [https://doi.org/10.1007/s12303-014-0043-5]
• Jo, K.W., Lee, H.J., Park, J.H., and Owen, J.S., (2010), Effects of Monsoon rainfalls on surface water quality in a mountainous watershed under mixed land use, Korea Journal of Agricultural and Forest Meteorology, 12(3), p197-206. [https://doi.org/10.5532/kjafm.2010.12.3.197]
• Kalbus, E., Reinstorf, F., and Schirmer, M., (2006), Measuring methods for groundwater-surface water interactions: a review, Hydrology and Earth System Sciences, 10, p873-887. [https://doi.org/10.5194/hess-10-873-2006]
• Karanth, K.R., (1987), Groundwater assessment, development and management, Tata McGraw Hill Education, New Delhi, p720.
• Kanduč, T., Grassa, F., McIntosh, J., Stibilj, V., Ulrich-Supovec, M., Supovec, I., and Jamnikar, S., (2014), A Geochemical and stable isotope investigation of groundwater/surface- water interaction in the Velenje Basin, Slovenia, Hydrogeology Journal, 22, p971-984.
• Kettering, J., Park, J.H., Lindner, S., Lee, B., Tenhunen, J., and Kuzyakov, Y., (2012), N fluxes in an agricultural catchment under monsoon climate: a budget approach at different scales, Agriculture, Ecosystem & Environment, 161, p101-111. [https://doi.org/10.1016/j.agee.2012.07.027]
• Ki, M.G., Koh, D.C., Yoon, H., and Kim, H.S., (2013), Characterization of nitrate contamination and hydrogeochemistry of groundwater in an agricultural area of Northeastern Hongseong, Journal of Soil and Groundwater Environment, 18(3), p33-51, (in Korean with English abstract). [https://doi.org/10.7857/jsge.2013.18.3.033]
• Kim, H., Kaown, D., Mayter, B., Lee, J.Y., Hyun, Y., and Lee, K.K., (2015), Identifying the sources of nitrate contamination of groundwater in an agricultural area (Haean basin, Korea) using isotope and microbial community analyses, Science of the Total Environment, 533, p566-575. [https://doi.org/10.1016/j.scitotenv.2015.06.080]
• Kwon, Y.S., Lee, H.H., Han, U., Kim, W.H., Kim, D.J., Kim, D.I., and Youm, S.J., (1990), Terrain analysis of Haean basin in terms of earth science, Journal of the Korea Earth Science Society, 11(3), p236-241, (in Korean with English abstract).
• Lee, J., Eom, J.S., Kim, B.C., Jang, W.S., Ryu, J.C., Kang, H., Kim, K.S., and Lim, K.J., (2011), Water quality prediction at Mandae watershed using SWAT and water quality improvement with vegetated filter strip, Journal of the Korean Society of Agricultural Engineers, 53(1), p37-45, (in Korean with English abstract). [https://doi.org/10.5389/ksae.2011.53.1.037]
• Lee, J.Y., Lee, K.S., Park, Y., Choi, H.M., and Jo, Y.J., (2013), Chemical and isotopic compositions of groundwater and stream water in a heavy agricultural basin of Korea, Journal of Geology Society India, 82, p169-180. [https://doi.org/10.1007/s12594-013-0135-2]
• Liu, C.W., Chen, S.K., Jou, S.W., and Kuo, S.F., (2001), Estimation of infiltration rate of a paddy field in Yun-Lin, Taiwan, Agricultural Systems, 68, p41-54.
• Min, J.H., Yun, S.T., Kim, K., Kim, H.S., and Kim, D.J., (2003), Geologic controls on the chemical behavior of nitrate in riverside alluvial aquifers, Korea, Hydrological Processes, 17, p1197-1211. [https://doi.org/10.1002/hyp.1189]
• Park, Y., Lee, K.S., and Yu, J.Y., (2006), Seasonal variation of dissolved ions and oxygen and hydrogen isotopic compositions of precipitation in Chuncheon, Korea, Journal of the Geology Society of Korea, 42(2), p283-292, (in Korean with English abstract).
• Paul, G., and Wanielista, M., (2000), Effects of evaporation enrichment on the stable isotope hydrology of a central Florida (USA) river, Hydrology Process, 14, p1465-1484.
• Piper, A.M., (1944), A graphic procedure in the geochemical interpretation of water analysis, Transactions American Geophysical Union, 25, p914-928.
• Schoeller, H., (1977), Geochemistry of groundwater, Chap, 15, UNESCO, Paris, p1-18.
• Sohn, S.H., Park, K.J., and Chung, S.O., (2004), Analysis of water balance for ponding depth treatment in paddy fields, Journal of the Korean Society of Agricultural Engineers, 46(4), p13-21, (in Korean with English abstract). [https://doi.org/10.5389/ksae.2004.46.4.013]
• Subramani, T., Elango, L., and Damodarasamy, S.R., (2005), Groundwater quality and its suitability for drinking and agricultural use in Chithar river basin, Tamil Nadu, India, Environmental Geology, 47, p1099-1110. [https://doi.org/10.1007/s00254-005-1243-0]
• Todd, D., (1980), Groundwater hydrology, 2nd edition, Wiley, New York.
• United Nations Environment Program (UNEP), (1999), Global environment outlook 2007, Earthscan, UK.
• United States Salinity Laboratory (USSL), (1954), Diagnosis and improvement of saline and alkali soils, Agriculture handbook 60, US Department of Agriculture (USDA), Washington, p69-81.
• Wang, H., Jiang, X.W., Wan, L., Han, G., and Guo, H., (2015), Hydrogeochemical characterization of groundwater flow systems in the discharge area of a river basin, Journal of Hydrology, 527, p433-441. [https://doi.org/10.1016/j.jhydrol.2015.04.063]
• Wilcox, L.V., (1955), Classification and use of irrigation waters, US Department of Agriculture, New York.
• Winter, T.C., Harvey, J.W., Franke, O.L., and Alley, W.M., (1998), Ground water and surface water a single resource, U.S. Geological Survey Circular, p1139.
• Won, J.K., Lee, M.W., and Woo, J.K., (1987), Geologic and geographic characteristics of Haean basin. Memorial Book of 40th Anniversary of Kangwon Nation University, Kangwon National University, p45-57, (in Korean).
• Won, J.K., Na, K.C., and Lee, M.W., (1988), Geology in northern area of Demilitary Zone (Geology), Cultural Heritage Administration of Korea, p114-120, (in Korean).
• Yang, Z., Zhou, Y., Wenninger, J., and Uhlenbrook, S., (2014), A multi-method approach to quantify groundwater/ surface water-interactions in the semi-arid Hailiutu River basin, northwest China, Hydrogeology Journal, 22, p527-541. [https://doi.org/10.1007/s10040-013-1091-z]
• Yun, S.W., Lee, J.Y., and Lee, H.G., (2015), Variation of stream water quality and baseflow contribution from groundwater during rainfall event in the Haean basin, Journal of Geology Society of Korea, 51(6), p611-621, (in Korean with English abstract). [https://doi.org/10.14770/jgsk.2015.51.6.611]