1. 서 론
염 수송(Salt transport)은 연안 및 하구에서의 오염물질의 확산, 해양 쓰레기의 이동, 퇴적물의 이동 등과 같은 물질들의 이동을 결정하는 중요한 요인 중 하나이다(Bowen and Geyer, 2003). 이러한 물질 이동의 특성을 파악하기 위해 염수송과 관련된 다양한 연구가 지속적으로 수행되어 왔다. Gong and Shen(2011)은 조석과 담수에 의한 수괴의 수직적 혼합이 중국 더우먼구 하구에서의 염 수송에 중요한 역할을 담당한다고 제시하였다. Wang et al.(2015)은 미국 허드슨 하구에서의 염 수송은 조석과 염분의 대조-소조의 변동성에 의하여 결정된다고 제시하였다. Kim and Park(2012)에 의하면 북부 멕시코만 하구에서 담수 유입에 의한 성층 강도와 바람에 의한 순압성 흐름의 상호 작용이 염 수송에 대한 메커니즘을 결정한다. 이와 같이, 염 수송은 조석, 담수, 바람, 지형적인 특성 등 여러 가지 외력의 상호작용에 의하여 수송 기작이 변화하기 때문에, 연구 지역의 특성에 따라 물질 이동의 양상이 상이하게 나타날 수 있다.
경기만 한강 하구는 영종도 매립, 송도 매립, 인천대교 건설, 신국제여객부두 건설 등 인위적인 요인에 의하여 해수의 유동 및 해저의 지형이 지속적으로 변화하고 있는 해역이다. 또한, 인간 활동에 의한 오염물질, 해양쓰레기, 준설토 등 여러가지 물질들이 해수 내로 유입되어 연안 생태계에 피해를 주고 있는 상황에서 물질 이동의 양상을 파악하는 것은 매우 중요하다. 경기만 한강 하구는 조석, 담수, 지형적인 특성에 의하여 시-공간적으로 염분의 변화가 뚜렷하게 나타나는 지역이기 때문에(Yoon and Woo, 2012), 주기적으로 변화하는 염분에 의해 물질 이동의 특성도 지속적으로 변화할 것이다. 따라서, 경기만 한강 하구에서 물질 이동의 양상을 파악하기 위해 염 수송에 대한 메커니즘을 정밀하게 분석할 필요가 있다.
경기만 한강 하구에서 물질 이동과 관련된 다양한 수리학적 연구가 수행되어 왔다. Lee et al.(2017)은 유속 관측을 통해 시-공간적으로 변화하는 순압력 구배에 의해 염하수로에서 낙조 방향의 물질 이동이 나타나며, 대조기에 강화되고 소조기에 약화되는 대조-소조의 변동성을 가진다고 제시하였다. Choi et al.(2012)은 석모수로에서의 유속 및 염분 관측을 통해 성층이 강하게 나타나는 지역의 저층에서 창조 방향의 물질 이동이 나타난다고 제시하였으며, Lee et al.(2012)에 의하면 약 8000 m3/s 이상의 담수 유입량은 염하수로의 표층에서 낙조 방향, 저층에서 창조 방향인 2층 흐름 구조의 물질이동을 발생시킨다. Yoon and Woo(2012)는 담수의 계절적인 변동에 의한 염분 분포를 분석하였으며, 염수 침입을 결정하는 주요 메커니즘은 조석, 담수 유입량, 하구의 지형적인 형태라고 제시하였다. 기존 연구를 종합하면, 경기만 한강 하구에서의 물질 이동은 주로 조석과 담수의 상호작용에 의해서 결정된다. 또한 담수 유입량의 계절적인 변동성에 의한 수평/수직적인 염분 분포의 변화가 물질 이동에 크게 영향을 미친다. 하지만, 경기만 한강 하구에서 물질 이동에 대해 중요하게 작용되는 염 수송에 관한 연구는 아직 부족한 상황이다.
염 수송에 대한 메커니즘을 분석하기 위해 단면 유속 및 염분에 대한 관측 자료가 필요하다. 하지만, 시-공간적으로 통일된 고해상도의 관측 자료를 수집하는 것은 실질적으로 제한점이 많다. 또한, 북방 한계선 부근에 위치한 강화도 북부수로는 군사적인 문제로 관측 자료의 수집이 제한적이다. 따라서, 본 연구에서는 3차원 유동 수치모델을 적용하여 경기만 한강 하구에서의 염 수송의 특성을 분석하였다. 현장 관측을 통해 유속 및 염분 시계열 자료를 취득하여 수치모델의 검증 자료로 활용하였다. Park et al.(2002)이 제시한 한강 하구에서의 대표적인 6개의 단면(L1~L6 in Figs. 1b and c)에 대하여 염 수송량(FS)을 산출하였다. 염 수송에 대한 발생 기작을 세부적으로 분석하기 위해 FS를 담수 유입에 의한 염의 이류 수송(QfS0), 수평 및 수직적인 유속 차이에 의한 염의 확산 수송(FE), 조석 변동성을 가진 유속과 염분의 위상 차이에 의한 염의 수송(FT)으로 분해하였다. 이를 통해, 2020년 갈수기와 홍수기 기간 동안의 사례 연구를 수행하여 염 수송의 시-공간적인 변동성을 분석하였다.
2. 연구 지역
대한민국 서해 중부의 연안에 위치한 경기만은 30개 이상의 크고 작은 섬들과 넓은 범위의 조간대가 분포하고 있다(Fig. 1). Datum level(DL) 기준, 수심은 외해에서 약 50 m 이상이며, 내륙에 근접한 대부분의 연안 지역에서는 약 10 m 이하이다. 강화도의 우측에는 약 700 m의 폭과 20 km의 길이를 가진 염하수로가 존재한다. 경기만 한강 하구는 한강의 입구 수로(L1), 강화도 북부 수로(L2), 염하수로 북부(L3)로 구성된 A1 지역(Fig. 1b)과 염하수로 남부(L4), 강화도 남부수로(L5), 영종도 동부 수로(L6)로 구성된 A2 지역으로 나누어진다(Fig. 1c).
경기만은 약 8 m 이상의 조차가 발생하는 대조차 환경이며, 주요 5개의 분조(M2, S2, K1, O1, N2)에 의한 조석이 지배적인 해역이다. 반일주조 성분인 M2와 S2의 합성에 의하여 약 14.76일의 주기를 가지는 대조-소조의 변동성이 나타나며, N2에 의하여 월조부등 현상이 발생한다. 또한, 일주조 성분인 K1과 O1에 의하여 일조부등 현상이 나타난다. 염하수로의 상류에 근접할수록 지형적 특성에 의하여 비선형 천해분조인 M4와 MS4가 증가하며, 조석 비대칭 현상이 발생한다(Yoon et al., 2015).
경기만은 담수 유입량의 계절적인 변동성이 뚜렷하게 나타나는 지역이며, 주요 담수 유입원은 한강, 임진강, 예성강이다(Yoon et al., 2015). 한강을 통한 담수 유입의 경우, 2010년부터 2019년까지 10년 동안의 동계(1월, 2월, 12월)와 하계(7월, 8월, 9월)에 해당하는 평균 담수 유입량은 각각 181, 1013 m3/s이다(Fig. 2). 이러한 담수 유입량의 계절적인 변동성은 한강 하구의 염분 분포와 경압성 흐름을 변화시키는 주요 메커니즘으로 작용된다(Yoon and Woo, 2012).
3. 연구 방법
3.1 현장 관측
영종도와 아라뱃길 사이에 위치하는 염하수로의 정점 M-YJ에서 유속 및 염분 관측을 수행하였다(M-YJ in Fig. 1c). 동계의 관측은 2020년 2월 4일부터 3월 17일까지 42일 동안 수행되었으며, 하계의 관측은 2020년 8월 22일부터 9월 21일까지 30일 동안 수행되었다. Acoustic Doppler Current Profiler(ADCP)를 이용하여 수층 별 유속 및 유향 시계열 자료를 수집하였으며, Conductivity-Temperature-Depth(CTD) sensor 2개를 이용하여 표층과 저층의 수온 및 염분 시계열자료를 수집하였다. ADCP는 Trawl-Resistant Bottom Mount(TRBM) 프레임을 사용하여 해저면에 고정된 상태로 상향식 관측을 수행하였다. 저층을 관측하기 위한 CTD는 TRBM에 추가로 설치하였으며, 표층을 관측하기 위한 CTD는 수면하 약 4 m 부근에 설치하였다. ADCP와 CTD의 세부 설정 값은 표에 제시하였다(Table 1). 현장 관측을 통해 수집한 유속 및 염분 시계열 자료를 바탕으로 수치모델의 결과와 비교하였다.
3.2 수치모델 구축
본 연구에서는 미국 메사추세츠 공과대학(University of Massachusetts Dartmouth, UMASSD)과 우즈홀 해양학 연구소(Woods Hole Oceanographic Institution, WHOI)에서 개발한 비정규 격자 해양순환 모델인 Finite Volume Coastal Ocean Model(FVCOM)을 이용하였다. 삼각형 형태의 비정규격자는 복잡한 해안선과 지형에 관한 문제를 해결하기에 적합하다(Chen et al., 2008). FVCOM은 연속 방정식, 운동 방정식, 열-염 보존 방정식을 포함한 지배 방정식으로 구성되어 있다(Chen et al., 2003; Chen et al., 2013).
여기서, u, v, w는 각각 x, y, z 방향의 유속 성분이며, T와 S는 각각 수온과 염분을 의미한다. ρ와 P는 각각 밀도와 압력을 나타내며, f와 g는 각각 Coriolis 계수, 중력 가속도 성분을 의미한다. Km과 Kh는 각각 수직 점성 계수 및 수직 확산 계수를 나타내며, Du, Dv, DT, DS는 수평적인 운동량 및 열-염의 확산항을 의미한다. FVCOM은 적분형태의 지배 방정식을 차분하여 삼각형 격자 내에서의 플럭스 계산을 통해 수치적으로 해결하는 유한 체적법을 사용한다. 또한, Mode-splitting 방식을 사용하여 운동 방정식을 계산한다. 3-D internal mode에서는 second-order Runge-Kutta time-stepping 방식을 사용하여 경압성 흐름을 계산하고, 2-D external mode에서는 fourth-order Runge-Kutta time-stepping 방식을 사용하여 순압성 흐름을 계산한다. FVCOM에 대한 자세한 설명은 FVCOM 매뉴얼에 제시되어 있다(Chen et al., 2013).
3.2.1 격자 구성
담수가 유입되는 한강, 임진강, 예성강과 강화도 서쪽의 석모수로, 동쪽의 염하수로를 포함하는 경기만 한강 하구에 FVCOM 모델을 적용하였다(Fig. 1). 격자의 동-서 방향의 길이는 약 150 km이며, 남-북 방향의 길이는 약 145 km이다. 외해 개방경계는 영종도로부터 서쪽 방향으로 약 85 km 떨어진 지역에서 남-북 방향으로 연장되어 있다. 한강과 임진강의 담수 유입경계는 각각 뚝섬유원지와 비룡대교 인근이며, 예성강의 담수 유입경계는 강의 입구로부터 북쪽 방향으로 약 25 km 떨어진 지역이다. UTM 좌표계로 구성된 수평 방향 격자의 노드와 셀의 개수는 각각 60882개, 115443개이며, 시그마 좌표계로 구성된 수직 방향 격자의 개수는 20개이다. 격자의 해상도는 한강, 임진강, 예성강 및 염하수로와 같은 좁은 수로에서 약 100~300 m이며, 외해 개방경계에 근접할수록 약 4500 m까지 확장된다. 국립해양조사원(Korea Hydrographic and Oceanographic Agency, KHOA)에서 제공하는 최신 수심 자료를 이용하여 평균해면(Mean Sea Level, MSL) 기준으로 수심을 구성하였다(Fig. 3).
3.2.2 초기 조건 및 경계 조건
수위와 유속에 대한 초기 조건은 모든 격자의 노드와 셀에 대해서 ‘0’으로 설정하였다. 조석에 대한 외해 개방경계 조건을 구성하기 위해 오리건주립대학교(Oregon State University, OSU)에서 개발된 Tidal Model Driver(TMD)를 이용하여 TPX08-atlas 모델 결과의 조화상수를 추출하였다. TPX08-atlas 모델 결과에서 추출한 5개의 분조(M2, S2, K1, O1, N2)에 대한 조화상수를 외해 개방경계에 보간 하여 입력하였다. 국립수산과학원(National Institute of Fisheries Science, NIFS)에서 제공하는 정선해양 관측자료의 염분 자료를 수집하여 외해 개방경계의 염분 조건을 구성하였다.
국가수자원관리종합정보시스템(WAter resource Management Information System, WAMIS)으로부터 2020년 1월 20일부터 9월 20일까지 1일 간격의 실시간 한강대교와 비룡대교의 유량 자료를 수집하였다. 수집한 유량 자료를 바탕으로 수치모델의 한강 및 임진강에서의 담수 입력자료를 구성하였다. 북방 한계선을 넘어선 지역에 위치한 예성강의 유량 자료를 확보할 수 없기 때문에, 인근 하천(한강, 임진강)과의 유역비를 산정하여 한강 유량의 15%를 예성강의 유량으로 활용하였다(Park et al., 2002). 염분에 대한 초기장을 구성하기 위해 2010~2019년 기간에 대한 2월의 평균 담수 유입량과 0 psu의 염분을 한강, 임진강, 예성강에 각각 적용하여 1년 동안 모의하였다. 1년 모의 후, 염분의 변화가 거의 일정한 안정화된 모델의 염분 공간장을 염분 초기장으로 입력하였다.
4. 모델 검증
모델의 모의 기간은 2020년 1월 20일부터 9월 20일이며, 동계와 하계 기간의 조위, 유속, 염분에 대한 모델 결과와 관측 자료를 비교 및 검증하였다. 모델과 관측 자료의 비교에 대한 정량적인 평가를 위해 식(2)를 이용하여 평균 오차(mean error, ME), 평균 절대 오차(mean absolute error, MAE) 및 Willmott(1981)에 의해 정의된 predicted skill(Skill)을 계산하였다.
Mn과 On은 각각 모델 결과와 관측 자료에 대한 n번째 값을 의미한다. O ¯
4.1 조위
국립해양조사원으로부터 2020년 1월 20일부터 3월 20일까지 강화대교, 영종대교, 인천, 인천-송도, 안산, 대산, 영흥도, 굴업도를 포함하는 총 8개의 조위 검조소에서의 1시간 간격의 실시간 조위 시계열 자료를 수집하였다(ST1-ST8 in Fig.1). 8개의 정점 중, 대표적인 4개의 정점 ST1, ST3, ST6, ST8에 대한 모델 결과와 관측 자료의 조위 시계열을 비교하였다(Fig. 4). 모델의 조위가 관측 자료에서 나타나는 대조-소조의 변동성, 일조부등, 월조부등 현상을 유사하게 재현하였다. 8개의 조위 검조소에 대하여, ME와 MAE는 각각 1.34~1.76 cm, 17.95~23.52 cm의 범위로 나타났으며, Skill은 모든 조위 검조소에 대하여 0.99의 값이 산출되었다(Table 2). Matlab 프로그램의 T-Tide toolbox를 이용한 조화분해를 통해 5개의 분조(M2, S2, K1, O1, N2)의 조화상수를 추출하여 모델 결과와 관측 자료를 비교하였다(Fig. 5). 모든 조위 검조소에서 5개의 분조의 진폭과 위상에 대한 모델 결과와 관측 자료가 거의 일치하였다(R2 = 0.99).
4.2 유속
염하수로의 정점 M-YJ(Fig. 1c)에서 ADCP를 통해 수집한 유속 자료와 모델 결과의 U-성분, V-성분을 산점도(Scatter plot)를 통하여 비교하였다(Fig. 6). 모델의 모의 기간 동안, 유속의 방향은 대부분 남-북 방향으로 나타나며, V-성분의 유속이 U-성분의 유속에 비하여 약 10배 이상 높다. 염 수송은 수로의 횡단 방향보다 종단 방향에 대한 역학과 관련이 높기 때문에(Lerczak et al., 2006), 주성분 분석(Principal Component Analysis, PCA)을 통해 변환한 수로 방향의 유속(Along-channel velocity)의 시계열에 대한 모델의 재현성을 검토하였다.
4.3 염분
염하수로의 정점 M-YJ(Fig. 1c)에서 CTD를 통해 수집한 염분 자료와 모델 결과의 염분 시계열을 동계(Figs. 7e and f)와 하계(Figs. 8e and f)의 기간에 대해서 비교하였다. ME는 0.01~1.39 psu의 범위이며, MAE는 0.69~1.47 psu의 범위이다(Table 2). 모델의 염분이 관측 자료에서 나타나는 창조-낙조의 변동성을 유사하게 재현하였으며, 동계와 하계 모두 0.87 이상의 높은 Skill이 산출되었다. 담수 유입량의 범위는 동계 검증 기간에서 약 160~490 m3/s이며(Fig. 7b), 하계 검증 기간에서 약 600~4600 m3/s이다(Fig. 8b). 동계에 대한 관측 자료와 모델 결과의 염분은 각각 24~30, 25~30 psu의 범위이며, 하계에서는 각각 4~25, 3~25 psu의 범위이다. 모델의 염분이 관측 자료에서 나타나는 담수 유입량에 의한 염분의 계절적인 변동성을 유사하게 재현했다.
동계의 일부 대조기 시기와 하계의 저층에서 모델 결과가 일부 과대 산정이 되었다. 모델 결과에 대한 염분의 변동성이 관측 자료와 유사하게 재현되기 위하여, 한강 하구 해역의 정확한 수심 자료가 필요하다. 하지만, 북방 한계선에 위치한 강화도 북부 수로에서의 정확한 수심 자료는 취득이 제한되기 때문에, 본 모델 결과의 염분의 재현성에 대한 한계점으로 작용한다.
염분의 공간적인 분포를 검증하기 위해 해양환경공단(Korea Marine Environment Management Corporation, KOEM)에서 제공하는 해양환경측정망으로부터 2010~2019년 기간의 2월, 5월, 8월에 대한 표층 및 저층의 공간적인 염분 자료를 총 8개의 정점에 대해서 수집하였다(W1-W8 in Fig. 1). 해당 염분 자료는 매년 분기마다 1회씩 제공되기 때문에, 각 정점별로 10년 동안의 2월, 5월, 8월 자료에 대한 통계 결과를 Box plot으로 표현하였으며, 각 정점에 해당하는 모델의 염분 결과와 비교하였다(Fig. 9). 2020년 2월, 5월, 8월에 해당하는 모델의 공간적인 염분 결과가 관측 자료의 최대, 최소 및 평균의 범위 내에서 재현되었다.
5. 염 수송의 메커니즘
경기만 한강 하구에서의 염 수송에 대한 메커니즘을 분석하기 위하여, 식(3)과 (4)를 통해 한강 하구의 6개 단면(L1-L6 in Figs. 1b and c)에 대한 체적 수송량(Qf)과 염 수송량(FS)을 계산하였다.
여기서, u는 단면의 법선 방향 유속 성분을 의미하며, S는 염분을 의미한다. 각 단면에서의 u의 양수 값은 한강 상류 방향이며, 음수 값은 외해 방향이다. < >는 48시간 low-pass filter를 적용하여 조석 성분을 제거한 잔차 성분을 의미한다. 염 수송에 대한 메커니즘을 세부적으로 분석하기 위하여, Lerczak et al.(2006)이 제시한 방법으로 FS를 3가지의 성분으로 분해하였다.
여기서, u0, S0, uE, SE는 조석 성분을 제거한 잔차 성분을 의미하며, uT와 ST는 잔차 성분을 제거한 조석 성분을 의미한다. u0와 S0는 단면-평균을 시행한 성분이며, uE, SE, uT, ST는 단면-평균을 시행하지 않은 성분이다.
QfS0, FE, FT는 조석과 담수의 상호작용에 의하여 결정된다(Kim and Park, 2012; Wang et al., 2015). QfS0는 체적 수송량에 의한 염의 이류 수송량을 의미하며, 담수 유입량의 세기와 평균 염분의 농도가 높을수록 Qf와 S0가 강해지기 때문에 수송량이 강화된다. FE는 수평 및 수직적 유속 차이가 야기시키는 전단류에 의한 염의 확산 수송량을 의미하며, 수괴의 성층에 의한 2층 흐름 구조가 강해질수록 수송량이 강화된다. FT는 조석 변동성을 가진 유속과 염분의 위상 차이에 의하여 발생하는 염 수송량을 의미하며, 창조 및 낙조에 대한 유속의 비대칭적인 특성에 의하여 수송량이 결정된다. 또한, 담수 유입량의 세기에 의한 염분의 수평적인 조석 변동성의 변화가 FT에 영향을 미친다.
5.1 체적 및 염 수송량
모델의 결과는 2020년 2월 4일~7월 17일의 갈수기 기간(Dry season; days 35~200)과 2020년 7월 18일~9월 16일의 홍수기 기간(Wet season; days 200~260)으로 구분된다(Fig.10a). 모델의 모의 기간 동안 L1-L6 단면에서의 Qf와 FS의 시계열을 갈수기와 홍수기의 기간에 대해 구분하여 제시하였다(Figs. 10b, c, d and e). 갈수기에 A1과 A2 지역 모두 Qf와 FS는 뚜렷한 대조-소조의 변동성이 나타난다. 반면에, 홍수기의 강한 담수 유입량은 Qf와 FS의 대조-소조의 변동성을 변형시킨다. 따라서, 갈수기에 Qf와 FS의 순 수송량은 대조-소조의 변동 특성에 의하여 결정될 것이며, 홍수기의 순 수송량은 담수 유입량의 세기에 의하여 지배될 것이다.
Qf와 FS의 순 수송량에 대한 메커니즘을 분석하기 위하여, Qf와 FS는 담수 유입량이 가장 낮은 2020년 4월 11일~5월 11일의 기간(Low discharge, LD; days 102~132)과 담수 유입량이 가장 높은 2020년 7월 28일~8월 27일의 기간(High discharge, HD; days 210~240)에 대해 각각 30일 시간-평균되었다(Q f ¯ F S ¯ Q f ¯ Q f ¯ Q f ¯ Q f ¯ Q f ¯ Q f ¯
LD와 HD 기간 모두 A2 지역에서의 체적 교환은 주로 L4와 L6에서 발생한다. LD 기간 동안 L4에서의 Q f ¯ Q f ¯ Q f ¯ Q f ¯ Q f ¯ Q f ¯
LD 기간 동안 외해로부터 한강 하구로 유입되는 염의 공급원은 L2와 L5이다. L2와 L5를 통한 F S ¯ F S ¯ F S ¯
HD 기간 동안 외해로부터 한강 하구로 유입되는 염의 공급원은 L5이며, 한강 및 임진강 상류에 잔류하고 있던 염이 L1을 통해 한강 하구로 다시 유입된다. L1와 L5를 통한 F S ¯ F S ¯
LD와 HD 기간 모두 L5와 L6에서의 F S ¯ F S ¯ F S ¯
5.2 단면을 통한 염 수송의 특성
5.2.1 A1 지역에서의 염 수송
갈수기 동안 A1 지역에서의 FS의 경향은 대조-소조의 변동성을 가진 QfS0와 FT에 의하여 결정된다. L1, L2, L3에서의 QfS0는 각각 -7.5 × 103~1.1 × 103, -11.4 × 103~2.0 × 103, -2.0 × 103~1.3 × 103 kg/s의 범위로 변화하며, 대부분 외해 방향으로 향한다(Figs. 12c, d and e). 주로 L1과 L2를 통해 Qf가 발생하는 A1 지역의 특성에 의하여(Fig. 12f), 외해 방향의 QfS0는 L3보다 L1과 L2에서 강하다. 반면에, FT는 항상 한강 방향으로 발생한다. L1과 L2에서의 FT는 각각 0.1 × 103 ~3.8 × 103, 0.8 × 103~5.1 × 103 kg/s의 범위이며, L3에서의 FT는 0.1×103~2.3 × 103 kg/s의 범위로 변화한다(Figs. 12c, d and e). A1 지역에서의 FT는 대조기에 강화되고 소조기에 약화되는 특성을 나타내며, L3보다 L1과 L2에서 높다. FE에 의한 염 수송은 항상 한강 방향을 향하지만, QfS0와 FT에 비하여 크기가 10배 이상 낮기 때문에 FS에 대한 기여도가 낮다. 따라서, 갈수기 동안 A1 지역에서의 FS는 QfS0와 FT가 균형을 이루어 결정된다.
홍수기 동안 A1 지역에서의 FS는 갈수기와 마찬가지로 외해 방향으로 향하는 QfS0와 한강 방향으로 향하는 FT가 균형을 이루어 결정된다. 하지만, FS의 변동성이 갈수기에 비해 약 20배 이하로 감소한다. L1, L2, L3에서의 QfS0와 FT의 변동범위는 80 kg/s 이하의 낮은 값을 가지며, 대조-소조의 변동성이 사라진다(Figs. 12c, d and e). 2020년 홍수기 기간에는 약 5.0 × 103~2.1 × 104 m3/s 범위의 첫 번째로 높은 담수 유입량이 나타나는 시기(R1; days 217~227)와 약 4.0 × 103 ~ 4.9 × 103 m3/s 범위의 두 번째로 높은 담수 유입량이 나타나는 시기(R2; days 246~252)가 존재한다(R1 and R2 in Fig.12b). R1과 R2시기에서의 약 4.0 × 103 m3/s 이상의 높은 담수 유입량은 외해 방향의 Qf를 갈수기에 비해 약 20배 이상 증가시킨다(Fig. 12f). 하지만, 담수의 낮은 염분에 의하여, S0가 0.1 psu 이하까지 감소하기 때문에(Fig. 12g), R1과 R2시기의 QfS0, FE, FT 및 FS는 갈수기에 비해 급격히 감소한다. 따라서, 담수 유입원과 가까운 A1 지역에서는 강한 담수 유입에 의하여 평균 염분이 0.1 psu 이하로 감소될 수 있으며, 이는 FS를 급격히 감소시킨다.
5.2.2 A2 지역에서의 염 수송
갈수기 동안 A2 지역에서의 FS는 A1 지역과 마찬가지로 대조-소조의 변동성을 가진 QfS0와 FT가 균형을 이루어 결정된다. L4, L5, L6에서의 QfS0는 각각 -3.2 × 103~3.1 × 103, -4.2 × 103~4.8 × 103, -6.4 × 103~6.1 × 103 kg/s의 범위로 변화하며, FT는 각각 0.3 × 103~1.9 × 103, 0.2 × 103~2.3 × 103, 0.1 × 103~0.8 × 103 kg/s의 범위로 변화한다(Figs. 13c, d and e). A2 지역에서의 FT는 A1 지역과 유사하게 대조기에 강화되고 소조기에 약화되는 특성이 나타나며, L6에서 변동의 폭이 가장 낮다. 갈수기에 FE는 거의 변동이 없으며, FS에 거의 영향을 미치지 않는다. 따라서, 갈수기 동안 A2 지역에서의 FS는 QfS0와 FT에 의하여 결정된다.
홍수기 동안 A2 지역에서의 FS는 강한 담수 유입량에 의해 대조-소조의 변동성이 변형된다. R1과 R2시기의 높은 담수 유입량에 의해 L4와 L6에서의 Qf는 각각 최대 -2.5 × 103, -2.2 × 103 m3/s까지 감소하여 뚜렷한 변화를 보이지만, L5에서의 Qf는 지역적인 특성에 의하여 크게 감소하지 않는다(Fig.13f). 담수의 낮은 염분은 L4, L5, L6에서의 S0를 224일 부근에서 각각 1, 16, 21 psu 이하까지 감소시킨다(Fig. 13g). 하지만, A1 지역과 같이 S0가 0.1 psu 이하로 지속되는 현상은 나타나지 않았다. R1과 R2 시기에 증가된 외해 방향의 Qf에 의하여 A2 지역에서의 QfS0는 장 기간 동안 외해 방향으로 유지된다(Figs. 13c, d and e). 또한, 홍수기의 높은 담수 유입량은 염분의 조석 변동성(ST)을 강화시키기 때문에, A2 지역에서의 FT를 증가시킨다. R1의 강한 담수 유입량에 의해 L4, L5, L6에서의 FT는 각각 최대 6.8 × 103, 8.4 × 103, 1.9 × 104 kg/s까지 증가된다(Figs. 13c, d and e). 따라서, 홍수기 동안 A2 지역에서의 QfS0와 FT의 변동은 FS의 변동에 크게 관여한다. 하지만, L6의 경우, FS의 변동에 대해 FE가 높은 기여도를 나타낸다.
홍수기 동안 L6는 6개의 단면 중에서 FE가 가장 높은 지역이며(Fig. 13e), 이는 염분 성층 및 유속의 수직적인 분포와 연관이 있다. 염분의 수직적인 분포는 조석과 담수 유입량의 세기에 의하여 결정되며, 상대적으로 담수 유입량의 세기가 조석의 세기보다 강할수록 표층과 저층의 염분 차이(∆SE)가 높아진다(Choi et al., 2012). L6에서 상대적으로 조석의 힘이 약한 소조기에는 대조기보다 ∆SE가 높으며, 홍수기의 높은 담수 유입량은 ∆SE를 18 psu까지 증가시켜 강한 성층을 형성한다(Fig. 14b). Choi et al.(2012)에 의하면, 하구에서의 강한 염분 성층은 수직적인 2층 흐름 구조를 야기시킨다. R1시기의 강한 성층은 표층과 저층의 유속 차이(∆uE)를 약 0.4 m/s까지 증가시키며, 2층 흐름 구조를 야기시킨다(Fig. 14a). R1시기의 높은 ∆uE에 의해 FE는 최대 3.1 × 104 kg/s까지 증가한다(Fig. 13e). 따라서, 홍수기 동안 L6 단면에서의 FS는 QfS0, FE, FT 모두가 균형을 이루어 결정된다.
5.3 염 수송의 계절적인 변동성
앞서 기술한 한강 하구에서의 염 수송에 대한 계절적인 변동 특성을 요약하기 위하여, 갈수기 및 홍수기에 해당하는 LD와 HD 기간에 대한 L1~L6의 F S ¯ Q f S 0 ¯ F E ¯ F T ¯
갈수기 기간의 F S ¯ F T ¯ Q f S 0 ¯ F E ¯ Q f S 0 ¯ F T ¯ F S ¯ Q f ¯ Q f S 0 ¯ F T ¯
홍수기 기간의 F S ¯ Q f S 0 ¯ Q f S 0 ¯ F T ¯ Q f S 0 ¯ F E ¯ F T ¯ F E ¯ F T ¯
한강 하구에서 F S ¯ Q f S 0 ¯ F T ¯ F T ¯ F S ¯ Q f S 0 ¯ F S ¯ F T ¯ F S ¯ Q f S 0 ¯ F S ¯ F S ¯ F S ¯
6. 결 론
대조차 환경의 경기만 한강 하구에서 염 수송에 대한 시-공간적인 변동성을 분석하기 위해 유동 수치모델을 적용하였다. 유동 및 염분의 변화에 관한 모델의 재현성을 검토하기위해 관측 자료와 모델 결과의 비교를 통하여 모델 검증을 수행하였다. 한강 하구의 주요 6개 단면을 선정하여 염 수송에 대한 메커니즘을 분석하였다.
경기만 한강 하구에서의 염의 순 수송량은 조석과 담수 유입의 상호작용에 의하여 결정된다. 담수 유입의 세기가 약한 갈수기에는 FT에 의해 한강 상류 방향으로 향하는 염 수송이 우세하며, 담수 유입의 세기가 강한 홍수기에는 QfS0에 의해 외해 방향으로 향하는 염 수송이 우세한다. 하지만, 지역적인 특성에 의하여 담수의 영향이 적은 L5는 갈수기와 홍수기의 구분 없이 한강 상류 방향의 염 수송을 발생시키는 FT가 강하기 때문에 외해로부터 지속적으로 염이 유입되며, 평균 염분이 가장 높은 L6는 갈수기와 홍수기의 구분 없이 외해 방향의 염 수송을 발생시키는 QfS0가 강하기 때문에 지속적으로 염이 외해로 유출된다. 결론적으로, 갈수기와 홍수기 모두 외해로부터 한강 하구로 유입되는 염의 대부분은 L5를 통해 수송되며, L6를 통해 외해로 유출된다. L1~L4에서의 염 수송은 갈수기에 한강 상류 방향으로 발생하며, 홍수기에 외해 방향으로 발생한다.
기존 연구 자료에서는 경기만 한강 하구에서의 물질 이동의 양상을 파악하기 위해 잔차류 및 체적 수송량에 대한 분석을 수행해왔다. 하지만, 본 연구에서 염 수송량을 통해 산정된 물질 이동의 방향과 잔차류 및 체적 수송의 방향은 다른 양상을 나타냈다. 따라서, 담수와 해수의 혼합이 활발하게 이루어지는 해역에서 물질 이동의 특성을 정확하게 파악하기 위해 염 수송에 대한 메커니즘을 정밀하게 분석할 필요가 있다.
PDF Links
PubReader
ePub Link
Full text via DOI
Download Citation 
Print 
