종점 기하학이 빙하 흐름 속도 제어

초록

**
그린란드 주요 출구 빙하의 종점 형태가 변하면, 부유 terminus가 짧아질수록 접지선 근처와 그 위쪽 4 km까지 흐름 속도가 급격히 증가한다. 3‑차원 전완전스톡스 모델 실험에서 7 km 길이의 부유 terminus를 단계적으로 제거했을 때, 접지선에서 흐름 속도가 최대 두 배까지 상승한다는 결과가 도출되었다. 이는 종점 기하학이 빙하 가속·감속을 좌우하는 주된 메커니즘이며, 해양 온도 변화에 의한 용융·칼빙이 이러한 기하학 변화를 유발한다는 결론을 뒷받침한다.

**

상세 분석

**
본 논문은 그린란드 출구 빙하의 종점(terminus) 기하학이 흐름 속도에 미치는 영향을 정량적으로 규명하기 위해 3차원 전완전스톡스(full‑Stokes) 모델인 FISMO를 활용하였다. 모델은 비압축성 흐름 방정식과 Glen’s 2‑flow 법칙을 기반으로 하며, 변형률에 따라 점성(η)이 변하는 비선형 시스템을 유한요소법(Q2Q1 등가형 테일러‑후드 요소)으로 해석한다. 경계조건은 베드록과의 접촉면에 무이동(Dirichlet) 조건, 대기와 접촉면에 무응력(Neumann) 조건, 해수와 접촉면에 수압(σ_n = –ρ_water g z) 조건을 적용하였다. 또한, 내부와 외부의 압력 경계조건을 통해 무한 슬래브와 유사한 동작을 재현하였다.

모델 설정에서는 Jakobshavn Isbræ를 대표하는 단순화된 피요르드 형태의 베드록을 사용했으며, 접지선(x = 0)을 기준으로 x > 0은 접지, x < 0은 부유 구역으로 구분하였다. 초기 부유 terminus 길이는 7 km였으며, 표면 및 바닥의 연간 순수량(b)과 해수 아래의 용융률(–50 m a⁻¹)을 적용해 정상 상태에 도달하도록 시뮬레이션을 진행하였다. 이후 부유 terminus를 2 km, 1 km, 0 km 순으로 절단하고, 기하학적 변화를 고정한 채 흐름 속도만 재계산하였다.

결과는 종점 길이가 짧아질수록 접지선 근처와 그 위쪽 4 km까지 평균 흐름 속도가 증가함을 보여준다. 특히 부유 terminus가 완전히 사라졌을 때 접지선 속도가 약 2배로 상승한다. 이는 부유 terminus가 존재할 경우 물압에 의해 발생하는 수평 압축력이 접지선에서의 변형률을 억제하고, 부유 빙상이 없을 경우 물압이 물면 아래에서만 작용해 광범위한 신장(stress) 상태를 초래하기 때문이다. 종점이 존재하면 수평 압축력이 접지선 단면 전체에 고르게 분포하지만, 종점이 사라지면 물압이 물면 이하에서만 작용해 물면 근처에 강한 신장 응력이 집중된다. 이러한 응력 재분배가 속도 증가의 직접 원인으로 작용한다.

또한, 모델은 부유 terminus가 존재할 때와 없을 때의 평균 응력(mean stress)과 편미분 응력(deviatoric stress)의 변화를 정량적으로 제시한다. 부유 terminus가 있을 때는 접지선에서 압축 응력이 우세하지만, 종점이 사라지면 신장 응력이 우세해져 빙하가 더 빠르게 흐른다. 이러한 메커니즘은 기존 연구에서 제시된 “빙상 버팀(buttressing)” 효과가 실제로는 종점 기하학에 의해 주도된다는 점을 강조한다.

논문은 또한 관측된 현상과 모델 결과를 비교한다. Jakobshavn Isbræ와 Helheim Glacier에서 관측된 급격한 속도 상승(칼빙 직후 15분 이내)과 모델에서 재현된 단계적 속도 증가가 일치한다. 다만 관측에서는 접지선 상류 40 km까지 속도 상승이 보이지만, 모델에서는 4 km 이내에 국한된 점은 종점 제거 후 표면 재조정(surface drawdown)을 고려하지 않았기 때문이다. 실제 상황에서는 부유 terminus가 얇아지면서 표면이 내륙으로 내려가고, 이로 인해 응력 재분배가 더 멀리까지 영향을 미친다.

결론적으로, 종점 기하학은 빙하 흐름 속도 변동의 주된 제어인자이며, 해양 온도 상승에 따른 부유 terminus의 용융·칼빙이 이러한 기하학 변화를 유발한다. 따라서 빙하 동역학을 예측하고 해수면 상승을 추정하기 위해서는 종점 형태와 그 변화를 정밀하게 모니터링하는 것이 필수적이다.

**