
Groundwater level and flow analysis indicates the existence of freshwater-saltwater interface in the aquifer.
The cation-exchange and precipitation reactions occurred in the aquifer due to seawater intrusion.
During dry seasons, the possibility of seawater inflow through groundwater increases.
석호는 사주, 사취의 발달로 바다와 분리된 호수를 말하며, 전 세계 해안에서 흔하게 관찰되는 연안의 지형이다. 화진포는 우리나라에서 가장 큰 규모의 석호로서 홀로세의 해수면 변동으로 인해 사주가 발달하여 약 5,500년 전부터 기수 환경이었을 것으로 추정된다(Yum
해수침투는 해안선 부근에서 해수가 대수층을 따라 담수 지하수계로 침입하는 현상으로, 특히 생태계와 긴밀한 관계가 있어 종합적인 연구가 필요한 분야이다(Werner
화진포의 경우 7번 국도와 해안도로로 인해 바다와의 연결부분이 좁혀져 있어 하구를 통한 자연적인 해수 유입이 어려운 조건이다. 따라서 호소수의 염분도가 유지되는 기작으로 호수 저면을 통한 해수 유입이 가능할 것으로 추측되며, 이는 이탈리아와 스페인에서 수행된 연구(Colombani and Mastrocicco, 2017; Menció
화진포(38°27′N, 128°28′E)는 강원도 고성군 거진읍에 있는 동해 연안의 석호로 유역 면적은 20.527km2, 호소 면적은 2.031km2, 수변 길이는 12.299km, 연평균 수위는 2.17m, 최대 수심은 4.62m이며 가지형 형태의 호소이다(Fig. 1). 화진포는 좁고 얕은 수로에 의해 내호와 외호가 연결되어 있고, 이러한 지형상의 특징으로 인해 두 호수는 서로 다른 깊이에서 염분약층을 형성한다(내호 수심 200cm, 외호 수심 80cm; Yum
화진포의 퇴적층 구조, 지층 특성, 대수층의 수리특성평가 및 지하수 유동 분석을 위하여 지하수 관측정을 설치하였다. 관측정은 수리·수문학적 상호관계(지표수, 해수, 지하수)를 파악하는데 필요한 지역에 설치되었는데, 화진포 외호의 하류의 사주지역 2개소(BH-1~2)와 상류 1개소(BH-3)이다(Fig. 1, Table 1). 설치된 관측정의 제원은 심도 15m, 굴착구경 100mm, 내부자재 직경 50mm이며, 상부 1~2m 구간을 제외한 전 구간에 스크린(유공관)이 설치되어 있으며, 지하수의 수직적인 현장 수질측정 및 수질분석 시료 채취 등에 활용되었다. 연구지역의 수위 분포 및 지하수 흐름 특성을 알아보기 위해 지하수관측정 3개소와 석호에서 수위를 측정하였으며, 수위 측정값을 이용하여 시기에 따른 지하수 흐름 특성을 파악하였다. 지하수 관측정에서는 기준점(지표)에서 지하수 심도를 측정하는데, 측정 기준점의 고도를 고려하여 해수면 기준 수위로 환산하여 이용하였다.
Table 1 . Details of observation wells at Hwajinpo lagoon
Well No. | TMX | TMY | Elevation (m) | K (cm/s) |
---|---|---|---|---|
BH-1 | 325,572.87 | 654,027.59 | 2.17 | 3.8×10-2 |
BH-2 | 325,422.96 | 654,003.36 | 0.97 | 2.2×10-2 |
BH-3 | 324,838.83 | 653,931.24 | 1.13 | 1.6×10-3 |
화진포 지역의 기반암은 경기 편마암복합체이며 지역적으로 쥐라기 대보화강암의 관입을 받았다. 호수 주변지역에서는 제4기 동안의 해안모래층이나 하천기원 충적층이 같은 기반암 위에 나타난다(Yum
연구지역과 인접한 기상관측지점 자료에 따르면(지점명: 대진, 강원도 고성군 현내면 초도리, 화진포에서 약 1.7km 거리) 화진포 주변 지역의 2020년 연 강수량은 2392.5mm이며, 연평균 기온은 13.28℃이다. 이 지역의 강수는 여름에 집중되어 있으므로(Fig. 2) 분석 결과에 계절성 영향을 최소화하기 위해 수위조사의 경우 3월과 10월, 시료채취 및 분석은 3월과 11월에 시행하였다.
지하수 수질의 수직적 분포 특성을 파악하기 위해 화진포 석호 주변에 설치된 지하수 관측정에서 공내 수질검층을 수행하였다. 공내 수질검층은 다항목 수질측정기(YSI-ProDSS, YSI사)를 이용하여 관정 상단에서 바닥까지 총 15m를 1m 간격으로 수행하였다. 지하수 샘플의 현장 수질검층 항목은 온도, 수소이온농도(pH), 전기전도도(EC), 염도(Salinity), 용존산소(DO), 산화환원전위(ORP)이며, 해수 영향을 고려하여 전기전도도(EC)는 50,000μS/cm 보정용액을 이용하여 보정 후 측정하였다.
지하수, 석호, 해수의 수질현황을 파악하기 위하여 시료채취 및 수질분석을 수행하였다. 지하수의 시료채취 방법은 정량이송펌프(Solinst사 Peristraltic pump, model 410)를 이용하여 상부심도에서 하부심도로 시료채취를 진행하여 수질 교란을 최소화하였으며, 실내 분석시료 채취 전에 현장수질을 측정하였다. 주요 양·음이온은 50ml Conical tube, HCO3는 100ml PE시료병 그리고 TOC 및 기타분석용은 1L 멸균채수병에 현장 여과를 진행한 뒤 채취하였다. 채취된 시료는 수질분석기관인 상지대학교 자연과학연구지원센터에 의뢰하여 분석하였으며, 주요 양이온(Ca2+, Mg2+, Na+, K+)은 유도 결합 플라즈마 광학 방출 분광법(ICP-OES; Thermo, iCAP 6500Dou)을 사용하였고 음이온(Cl-, SO42-, Br-)은 이온 크로마토그래피(IC; Metrohm, Eco IC)를 사용하여 측정되었습니다. 물안정동위원소는 Beta Analytic Inc.에 의뢰하여 분석하였고 공동 광자 감쇠 분광법(CRDS; Picarro, L2140-i)을 사용하여 δ18O와 δ2H를 측정하였으며 정밀도는 각각 0.01~0.23‰, 0.05~0.57‰이다.
수위는 해발고도(Elevation, above mean sea level)로 표시하였고, 측정 결과 전반적인 지하수위는 EL. 0m(해수면)~EL. 0.8m 부근에 분포하는 것으로 나타났다. 석호 상류의 BH-3 관측정 지하수위가 EL. 0.54~0.8m, 사주 지역 BH-1 및 BH-2 관측정 지하수위는 EL. -0.05~0.18m로 나타나 상대적으로 사주 지역이 낮은 경향을 보인다. 시기별 측정 수위를 이용하면 지하수위 분포단면도를 작성할 수 있으며, 이를 이용하여 지하수 분포를 파악할 수 있다(Fig. 3). 석호 상류지역에서는 석호 및 바다 쪽으로 흐름이 항시적으로 유지되는 반면, 사주 지역의 지하수위는 사주에서 석호 및 바다쪽으로 양방향 흐름이 나타나기도 하지만 바다쪽에서 사주 및 석호 쪽으로 흐름이 형성되기도 하였다.
해수침투는 담수에 비하여 상대적으로 밀도가 높은 해수가 대수층 내부로 유입되는 현상이기 때문에 두 유체 사이의 경계면은 해안에서 내륙으로 갈수록 깊어지는 형태를 가진다. Ghyben-Herzberg가 제안한 식(Verruijt, 1968)을 이용하여 화진포호 주변 대수층의 담수-염수 경계면을 분석한 결과, 담수-염수 경계면은 해안에서 BH-3 관측정 방향으로 깊어지는 경향을 보인다(Fig. 4).
2020년 10월에 측정된 수위는 내륙 관측정(BH-3)에서는 높게 나타나나, 사주지역의 BH-1, BH-2 관측정은 해수면 부근이며, 석호 수위는 해수면보다 낮다(Fig. 4b). 이 시기는 상대적으로 3월에 비해 온도는 높고, 일정 규모 이상의 강우가 지속되지 않은 기간이 긴 상황이 이어졌기 때문에 석호 수위가 해수면보다 낮아진 것으로 판단된다. 이러한 수위에 만일 Ghyben-Herzberg 식을 그대로 적용하면 전체가 염수일 수밖에 없는데, 이는 현장의 상황은 정류상태를 가정한 Ghyben-Herzberg 식과 달리 평형이 이루어지지 않은 부정류 상태이기 때문이다. 따라서, 이 시기의 담수-염수경계면은 평형에 가까운 3월의 형태에서 어떤 방향으로 이동했을 것인가를 추정하는 수준에서 가능하다. 석호와 사주 지하수의 수위가 하강하였으므로 경계면은 상부로 상승하였을 것으로 추정되며, 특히 석호 하부에서는 경계면이 석호 바닥보다 높게 상승하였을 가능성이 존재한다. 또한 염분의 확산 작용에 의해 담수-염수 경계면은 점이대를 형성하기 때문에 석호 바닥으로 염수(해수)의 유입이 발생할 수 있을 것으로 추정된다.
전기전도도(Electrical Conductivity, EC)는 석호 주변의 염수 영향을 확인할 수 있는 지표 중 하나이므로 석호를 포함한 EC 수직단면 분포도를 분석하면, 석호 하부의 담수-염수 경계면의 위치를 추정할 수 있다(Fig. 5). 담수-염수 경계면은 EC 값이 해수의 절반 수준인 약 25,000μS/cm의 부근에서 분포할 수 있는데, EC 분포에 의하면 바다에서 석호 하부로 쐐기형태의 담수-염수 경계면이 위치하는 것으로 생각된다. 또한 뚜렷한 경계면(sharp interface)보다는 점이대(transition zone)의 모습을 나타내었다. 전기전도도 변동 구간은 BH-1 관측정의 4~6m 심도와 BH-2 관측정의 5~10m 심도에서 전기전도도의 변화가 크게 일어났는데, 이것은 이론적으로 계산한 담수-염수 경계면과 잘 일치한다고 할 수 있다.
대수층에서 지하수와 해수가 혼합되면 지질매체와 다양한 지구화학적 반응이 일어나며 특히 양이온 교환반응(cation exchange)이 중요한 반응 중 하나이다. 이온델타(ionic delta)는 이 양이온 교환반응을 나타내는 지표로 활용할 수 있는데(Liu
포화지수(Saturation index, SI)는 광물의 침전 여부를 확인하는 지표로 활용할 수 있으며 PHREEQC(Parkhurst and Appelo, 1999)를 사용해 계산하였다. 연구지역의 시료분석 결과를 포화지수로 나타내 보면(Fig. 6c, 6d), 마그네슘을 함유하는 광물인 돌로마이트(Dolomite)의 경우 평형상태인 0보다 큰 값을 주로 보여 침전된 상태임을 알 수 있으며, 칼슘을 함유하는 광물인 방해석(Calcite)의 경우 0과 가깝거나 낮은 값을 보이는 경우가 많아 용해된 상태를 의미한다. 또한 해수 혼합비와 각 이온델타 값을 비교해 보았을 때(Fig. 7), 혼합비가 증가함에 따라 Ca2+은 증가하였고, Mg2+과 Na+은 감소하는 추세를 보여 앞에서 분석한 수리지구 화학적 특성이 해수 혼합비와 관계에서 비례적으로 나타난다고 볼 수 있다. 이러한 결과는 지질매체가 가지고 있던 Ca2+이 해수 기원인 Mg2+과 Na+으로 서서히 대체되는 과정을 보여주는 것으로 이해할 수 있다(Andersen
파이퍼 다이어그램(Piper diagram)으로 연구지역의 지하수의 수질유형을 살펴보면(Fig. 8) 수질분석 결과 대부분 Na-Cl 수질유형이며, 다수의 지하수가 해수와 동일한 유형으로 나타나고 있다. 반면에, BH-3 관측정의 수질 분석결과는 주요 양이온은 해수에서 우세한 Na++K+이 높지만, 음이온은 담수의 특성인 HCO3-의 농도가 가장 높게 나타났다. 이는 담수와 해수의 혼합수 중 해수인 지하수에서 담수가 유입되었을 경우에 발생하는 형태이며, BH-3 관측정의 경우 해수 혹은 호소수로 이루어진 지하수에 배후지로부터 담수가 유입되고 있음을 지시하고 있다.
물의 산소, 수소 안정동위원소는 지표수 또는 지하수 수체의 기원과 혼합을 추적할 때 유용하게 이용할 수 있다. 특히 지하수는 강수 기원이기 때문에 동위원소 분별작용에 의해 상대적으로 가벼운 동위원소 조성(음의 값)을 보이는 특성이 있다. 그리고 해수는 기준값(0‰) 부근의 값을 가지므로 해수의 침투 및 혼합을 평가할 때 유용하다. 또한 상온에서는 암석 또는 광물 속의 물과 반응하지 않는 특성이 있어 물안정동위원소 값이 변하지 않기 때문에 추적자로서 효과적으로 활용할 수 있다(Lee, 2018).
이번 연구의 수질시료 채취 및 분석과정에서도 산소, 수소 안정동위원소 분석을 수행하였다(Fig. 9). 산소·수소 동위원소 조성은 해수와 담수 성분 사이의 혼합된 직선(seawater mixing line)이 형성되어 해수침투에 따른 영향을 받은 것으로 판단된다. 이번 연구에서 대수층에 해수가 얼마나 침투하였는지 파악하기 위해 혼합비를 계산하였고, 추적자로는 염소이온과 산소동위원소를 활용하였다. Fig. 10은 염소이온과 산소동위원소 분석 결과를 활용해 추정한 혼합비를 나타낸 것으로, 혼합비의 평균값은 염소이온의 경우 0.28, 산소동위원소의 경우 0.32이고 전체적으로는 0.3의 평균값을 가진다(Table 2). 각 결과값은 1:1의 관계를 이루어 추적자로 활용한 두 방법 사이에 유사한 경향을 가지는 것으로 이해할 수 있다. 해안에서 가장 멀리 떨어진 BH-3 관측정의 경우 혼합비가 매우 낮아 담수와 가까운 성분을 보이며, 사주 지역의 두 관측정은 낮은 심도에서는 0.2 이하의 낮은 값을 보이나 깊은 심도에서는 0.4~0.6 정도의 혼합비를 보인다. 특히 이론적으로는 바다에 가까운 BH-1 관측정에서 더 높은 혼합비가 나타나야 하는데, BH-2 관측정의 깊은 심도에서 혼합비가 0.6 이상으로 매우 높은 값을 보이고 있다. 이를 이해하기 위해서는 추가적인 연구를 수행하여 BH-2 관측정의 15m 부근의 대수층 특성을 파악할 필요가 있다.
Table 2 . Seawater mixing ratio based on Cl concentrations and δ18O
Sample | Cl | δ18O | |
---|---|---|---|
Mar. 2020 | BH-1 (5m) | 0.09 | 0.13 |
BH-1 (10m) | 0.35 | 0.41 | |
BH-1 (15m) | 0.44 | 0.52 | |
BH-2 (3m) | 0.17 | 0.18 | |
BH-2 (8m) | 0.35 | 0.41 | |
BH-2 (15m) | 0.66 | 0.73 | |
BH-3 (7m) | 0.00 | 0.00 | |
BH-3 (15m) | 0.00 | 0.00 | |
Lagoon | 0.29 | 0.45 | |
Nov. 2020 | BH-1 (5m) | 0.09 | 0.01 |
BH-1 (10m) | 0.60 | 0.61 | |
BH-1 (15m) | 0.57 | 0.61 | |
BH-2 (3m) | 0.03 | 0.00 | |
BH-2 (8m) | 0.36 | 0.37 | |
BH-2 (15m) | 0.87 | 0.87 | |
BH-3 (7m) | 0.00 | 0.00 | |
BH-3 (15m) | 0.00 | 0.00 | |
Lagoon | 0.24 | 0.40 | |
Mean | 0.28 | 0.32 | |
0.30 |
이번 연구에 의하면 연구 대상지역인 화진포의 사주는 두께 15m 이상의 모래층이 분포하고 있다. 이 모래층은 투수성이 매우 양호한 대수층으로 바다에서 석호 방향으로 쐐기형태의 담수-염수 경계면(점이대)이 분포하며, 위치는 석호 및 지하수 수위 변화에 영향을 받는 것으로 확인된다. 강수량이 적은 기간이 지속되면, 석호 및 사주지역에 대수층에서 증발산에 의하여 수위가 해수면 부근 또는 그보다 낮게 분포하는 시기가 존재하는데, 이러한 수위 하강은 담수-염수 경계면이 내륙방향으로 이동(진출)하게 하며, 석호 하부에서는 경계면의 상승을 유발할 수 있다. 특히 석호 수위가 해수면보다 낮은 기간이 지속되면, 담수-염수 경계면이 상승하여 석호 바닥으로 염수(해수)의 유입이 발생할 수 있다.
수리지구 화학적 특징으로는 해수침투가 발생함에 따라 대수층의 지질매체가 가지고 있던 칼슘 이온이 해수 기원인 마그네슘과 나트륨 이온으로 대체되는 현상이 나타났다. 이온델타와 포화지수를 통해 분석한 결과 해수침투로 인한 수리지구 화학적 변화가 일어나 칼슘을 함유하는 방해석은 용해되고 마그네슘을 함유하는 돌로마이트는 침전이 일어났음을 알 수 있다. 또한 파이퍼 다이어그램으로 수질 특성을 분석한 결과 대부분의 지하수와 호소수가 해수와 동일한 유형으로 나타나 해수의 영향을 크게 받았음을 확인하였다. 물안정동위원소 분석 결과 해수 침투에 따른 영향을 받은 것 판단된다. 또한 산소동위원소를 추적자로 활용하여 계산한 해수의 혼합비는 염소이온을 추적자로 활용한 결과와 유사한 경향을 보였다. 해수 혼합비는 해안과의 거리가 가까울수록, 심도가 깊을수록 높은 값을 보이는 경향이 있으나 일부 결과는 그 부근의 대수층 특성으로 인해 예상한 바와 상이하였다.
이번 연구의 결과는 주로 정성적인 측면에서 대수층을 통한 해수침투 기작과 그 영향을 검토하였기 때문에 결과의 타당성 및 보편성을 확보하기 위해서는 추가적인 연구가 필요할 것으로 판단된다. 특히 지속적인 수위 및 수질 모니터링을 수행하여 계절적인 변화를 분석한다면 다양한 조건에서의 해수침투 영향을 파악할 수 있을 것이다. 또한 담수-염수 경계면 계산식은 평형상태를 가정한 것으로, 실제 현상은 부정류 상태이므로 수위변화를 고려한 부정류 수치모사를 통하여 현재의 시나리오를 구현할 필요가 있다.
본 연구는 환경부 “지중환경오염·위해관리기술개발사업; 2020 SEM (Subsurface Environment Management) projects-2020002470001”으로 지원받은 과제임
![]() |
![]() |