인용:
윤순봉. 2025. 과거 바다의 해수면을 복원하는 과학자들의 탐구 과정. 윤순봉의서재.
과거 바다의 해수면을 복원하는 과학자들의 탐구 과정 v2
목차
- 1.1. 직접적인 지질학적 기록: 해안의 흔적을 읽다
- 1.2. 간접적인 대리자료: 해양의 동위원소 기억을 풀다
- 1.3. 지구물리학적 보정: 빙하 지각 평형 조정의 필수적 역할
- 1.4. 융빙수 기원 추적의 패러다임 전환: 해수면 지문의 물리학
과거 지구의 해수면 변동을 이해하는 것은 미래 기후 변화에 따른 해수면 상승을 예측하고 대비하는 데 있어 핵심적인 과학적 과제이다. 특히 최후간빙기, 최후빙기절정기, 그리고 그 이후의 급격한 해빙기에 발생했던 해수면의 역동적인 변화는 오늘날 온난화되고 있는 지구의 빙상이 어떻게 반응할지에 대한 중요한 단서를 제공한다.
용어해설
최후간빙기(Last Interglacial, LIG): 지금으로부터 약 13만 년 전(약 129–116 ka), 현재보다 기온이 약간 더 따뜻했던 시기.
최후빙기절정기(Last Glacial Maximum, LGM): 약 2만 6천 년 전(약 26.5–19 ka), 지구 역사상 마지막으로 빙하가 가장 넓게 분포했던 가장 추운 시기.
빙상(Ice Sheet): 그린란드나 남극처럼, 대륙을 넓게 덮고 있는 거대한 얼음 덩어리.
이에 최근 수십 년간 고기후학 및 지구물리학 분야에서 이루어진 핵심적인 연구 성과들을 종합하고 비판적으로 분석하여, 후기 제4기 해수면 변동에 대한 최신 과학적 이해를 정리한다.
과거 해수면을 복원하는 작업은 단편적인 지질학적 증거들을 지구물리학적 모델과 결합하는 복잡하고 반복적인 과정을 통해 이루어진다. 해안가에 남겨진 고대 산호초의 흔적부터 심해 퇴적물에 기록된 미세한 동위원소 비율의 변화에 이르기까지, 다양한 종류의 대리자료들은 각기 다른 시간적, 공간적 해상도와 불확실성을 내포하고 있다. 이러한 자료들을 해석하기 위해서는 빙상과 해수의 질량 재분배에 따른 지각의 변형, 즉 빙하 지각 평형 조정과 같은 지구물리학적 과정을 정밀하게 모델링해야 한다.
용어해설
대리자료(Proxy Data): 온도계가 없던 과거의 기온을 나이테나 빙하 코어로 추정하듯, 과거의 환경을 간접적으로 알려주는 모든 종류의 자연 기록.
빙하 지각 평형 조정(Glacial Isostatic Adjustment, GIA): 무거운 빙상이 땅을 누르면 지각이 가라앉고, 빙상이 녹아 가벼워지면 지각이 다시 솟아오르는 현상. 마치 매트리스 위 볼링공을 치우면 매트리스가 서서히 원래 높이로 돌아오는 것과 같다.
본 메모는 먼저 과거 해수면 복원에 사용되는 핵심적인 방법론들을 체계적으로 검토하고, 각 방법론의 원리와 한계를 명확히 할 것이다. 이를 바탕으로, 미래 온난화 시나리오의 유사 사례로 주목받는 최후간빙기의 해수면 최고점, 빙하기 최성기의 해수면 최저점, 그리고 인류 문명에 지대한 영향을 미쳤을 급격한 해빙기 해수면 상승, 특히 ‘융빙수 펄스 1A’의 미스터리를 심층적으로 탐구한다.
용어해설
융빙수 펄스 1A(Meltwater Pulse 1A, MWP-1A): 약 1만 4,600년 전, 빙하가 갑자기 엄청난 속도로 녹아내려 해수면이 짧은 기간에 폭발적으로 상승한 사건.
이 과정에서 최후빙기절정기 해수면 최저점을 약 -134m로 정량화한 연구(Lambeck et al., 2014), 해빙기 중 발생한 MWP-1A의 규모와 속도를 규명한 연구(Deschamps et al., 2012), 그 융빙수의 주된 기여원이 북반구 빙상(북미·스칸디나비아)임을 보이며 남극 기여는 비지배적(0–35% 범위)일 가능성을 제시한 연구(Lin et al., 2021), 그리고 최후간빙기 해수면 최고점이 기존 통설보다 낮았을 가능성을 제시한 최근 연구들(Dyer et al., 2021; Dumitru et al., 2023)에 이르기까지, 각 시대의 해수면 변동 규모와 원인에 대한 과학적 합의가 어떻게 형성되고 수정되어 왔는지를 추적한다.
1. 과거의 흔적을 찾는 과학의 도구들
과거 해수면의 높이와 변화 속도를 재구성하는 과학은 관측 자료와 지구물리학적 모델링 사이의 복잡하고 상호보완적인 관계에 깊이 의존한다. 어떠한 단일 방법론도 완벽한 해답을 제공하지 못하며, 다양한 접근법을 통해 얻은 증거들을 종합하고 교차 검증함으로써 비로소 신뢰할 수 있는 과거의 그림을 그릴 수 있다. 이 장에서는 고해수면 연구의 근간을 이루는 핵심적인 방법론들을 소개하고, 각 방법의 원리와 내재된 불확실성을 비판적으로 검토한다.
1.1. 직접적인 지질학적 기록: 해안의 흔적을 읽다
과거 해수면을 직접적으로 추정하는 가장 강력한 증거는 과거 해안선 근처에서 형성된 지질학적 특징, 즉 ‘해수면 지시자’에서 나온다 (Rovere et al., 2023).
용어해설
해수면 지시자(Sea-level Indicator): 과거 특정 시점의 해수면 높이가 어디였는지를 알려주는 산호초, 해안 동굴, 해안단구와 같은 모든 종류의 지질학적 증거.
해수면 지시자의 서식 수심(indicative meaning)과 측정·표고 오차는 지역 및 생태 조건에 따라 분산이 크며, 메타데이터 표준화와 불확실도 전파가 중요하다. 최근 WALIS v1.0(The World Atlas of Last Interglacial Shorelines)은 최후간빙기 지시자 메타데이터 체계를 정비하여 자료 선정 및 품질 통제를 강화하고 있다(Rovere et al., 2023).
이 중 가장 중요한 기록 보존소는 화석 산호초이다. 산호는 특정 수심 범위, 특히 얕은 바다에서만 서식하기 때문에, 화석 산호초의 발견 위치와 연대는 과거 특정 시점의 해수면 높이에 대한 강력한 제약을 제공한다 (Lambeck et al., 2014).
특히 중요한 것은 ‘원격장’에 위치한 지점의 기록이다. 원격장이란 과거 거대 빙상이 존재했던 북미나 스칸디나비아로부터 아주 멀리 떨어진 지역을 의미하며, 타히티(Tahiti)가 대표적인 예이다. 이러한 지역은 빙상의 직접적인 압력에 의한 복잡한 지각 변동에서 비교적 자유롭고, 느리고 예측 가능한 속도로 지각이 침강하기 때문에 고품질의 해수면 기록을 얻기에 이상적인 장소로 여겨진다 (Camoin et al., 2007). 국제해양탐사프로그램(IODP)의 310차 탐사(Expedition 310)는 타히티의 해저 산호초를 시추하여 얕은 수심 산호의 우라늄–토륨 연대와 성장 수심 정보를 정밀하게 결합함으로써, MWP‑1A의 기간(14,650–14,310년전), 규모(14–18 m), 속도(≥40 mm/yr)를 고해상도로 제약하는 데 결정적인 역할을 했다(Deschamps et al., 2012).
이러한 지질학적 기록의 신뢰성은 연대 측정의 정밀도에 크게 좌우된다. 가장 널리 사용되는 방법은 우라늄-토륨 연대 측정법이다. 이 방법은 산호가 성장할 때 해수로부터 자신의 골격(탄산칼슘, CaCO₃)에 미량의 우라늄(U)을 포함시키지만 토륨(Th)은 거의 포함시키지 않는다는 점을 이용한다. 시간이 지나면서 산호 골격 내의 우라늄(²³⁸U, ²³⁴U)이 방사성 붕괴를 통해 토륨(²³⁰Th)으로 변환된다. 따라서 시료에 포함된 우라늄과 토륨 동위원소의 비율을 측정하면 산호가 형성된 시점을 정확하게 계산할 수 있다. 이 방법은 마치 모래시계처럼, 시간이 흐르면서 위에 있던 우라늄(모래)이 아래의 토륨(모래)으로 쌓이는 원리와 같다.
이 방법의 가장 중요한 가정은 산호가 형성된 이후 외부 환경과 우라늄이나 토륨을 교환하지 않는 ‘닫힌 계’를 유지했다는 것이다 (Edwards, Gallup and Cheng, 2003).
용어해설
닫힌 계(Closed System) / 열린 계(Open System): 연대 측정에서 ‘닫힌 계’는 외부 물질의 출입 없이 내부의 방사성 붕괴만 일어나는 완벽한 ‘타임캡슐’ 상태를 의미한다. 만약 지하수 등에 의해 우라늄이 씻겨나가거나 외부의 토륨이 들어오는 ‘열린 계’가 되면, 시간 기록이 왜곡되어 정확한 연대를 알 수 없게 된다.
만약 지하수 등에 의해 우라늄이 유실되거나 외부 토륨이 유입되는 ‘열린 계’의 특성을 보이면 연대 측정값에 오차가 발생할 수 있다. 따라서 연구자들은 시료의 결정 구조 분석 등을 통해 가장 보존 상태가 좋은 시료를 엄격하게 선별하는 과정을 거친다 (Dumitru et al., 2023).
1.2. 간접적인 대리자료: 해양의 동위원소 기억을 풀다
산호초와 같은 직접적인 지시자는 특정 시기와 장소에 대한 고해상도 정보를 제공하지만, 기록이 불연속적이라는 한계가 있다. 수십만 년에 걸친 연속적인 해수면 변화를 파악하기 위해 과학자들은 심해 퇴적물 코어(deep-ocean sediment core)라는 또 다른 중요한 기록 보존소를 활용한다 (Spratt and Lisiecki, 2016).
핵심 원리는 퇴적물에 포함된 저서성 유공충이라는 미세 해양 생물의 탄산칼슘 껍데기에 기록된 산소 동위원소 비율(δ¹⁸O)을 분석하는 것이다.
용어해설
저서성 유공충(Benthic Foraminifera): 바다 깊은 곳 바닥에 사는, 껍데기를 가진 아주 작은 단세포 생물. 이들의 껍데기는 과거 바다 환경의 정보를 담고 있는 ‘블랙박스’와 같다.
산소 동위원소 비율(Oxygen Isotope Ratio, δ¹⁸O): 물 분자(H₂O)를 구성하는 산소에는 가벼운 산소(¹⁶O)와 무거운 산소(¹⁸O)가 있다. 가벼운 산소를 포함한 물이 더 쉽게 증발한다. 빙하기에는 증발한 수증기가 육지의 거대한 빙상이 되어 바다로 돌아오지 못한다. 그 결과, 바다에는 상대적으로 무거운 산소(¹⁸O)의 비율이 높아진다. 유공충은 이 비율을 자신의 껍데기에 그대로 기록하기 때문에, 껍데기의 δ¹⁸O 값을 분석하면 당시 빙하의 양(즉, 해수면의 높이)을 역으로 추정할 수 있다.
이 δ¹⁸O 값은 주로 두 가지 요인에 의해 결정된다: 유공충이 살았던 당시의 심해 수온과 해수의 산소 동위원소 구성비이다 (Spratt and Lisiecki, 2016). 해수의 동위원소 구성비는 지구 전체의 빙상 부피와 직접적으로 연관된다. 물 분자(H₂O) 중 가벼운 동위원소인 ¹⁶O를 포함한 물 분자는 무거운 ¹⁸O를 포함한 물 분자보다 더 쉽게 증발한다. 빙하기에는 증발된 수증기가 대륙 빙상으로 쌓여 해양으로 돌아오지 못하므로, 해양에는 상대적으로 무거운 ¹⁸O가 풍부해진다. 반대로 간빙기에는 빙상이 녹아 가벼운 ¹⁶O가 풍부한 융빙수가 해양으로 유입되어 해수의 δ¹⁸O 값이 낮아진다.
따라서 유공충 껍데기의 δ¹⁸O 기록에는 수온 변화와 빙상 부피(즉, 해수면) 변화의 신호가 함께 섞여 있다. 이 두 신호를 분리하는 것이 연구의 핵심 과제이다. 이를 위해 과학자들은 유공충 껍데기의 마그네슘/칼슘(Mg/Ca) 비율과 같은 다른 대리자료를 이용하여 과거 수온을 독립적으로 추정한 뒤, δ¹⁸O 신호에서 수온 효과를 제거하여 순수한 빙상 부피 변화 신호를 추출한다 (Spratt and Lisiecki, 2016).
이러한 방법론적 발전을 집대성한 대표적인 연구가 Spratt & Lisiecki (2016)의 연구이다. 그들은 다양한 연구에서 각기 다른 방식으로 복원된 여러 해수면 기록들을 모아 주성분 분석이라는 통계적 기법을 적용했다.
용어해설
주성분 분석(Principal Component Analysis, PCA): 여러 데이터에 공통적으로 숨어있는 가장 중요하고 강력한 신호(주성분)를 찾아내고, 각 데이터의 특이한 값이나 오류 같은 ‘잡음(noise)’은 걸러내는 통계 기법. 이를 통해 전 지구적인 해수면 변화라는 핵심 신호를 명확하게 분리할 수 있다.
이를 통해 여러 기록에 공통적으로 나타나는 전 지구적 해수면 변화 신호를 강화하고, 각 기록의 지역적 특성이나 측정 오차와 같은 ‘잡음(noise)’을 줄인 표준화된 ‘해수면 스택(sea level stack)’을 생성했다 (Spratt and Lisiecki, 2016). 그들의 분석 결과, 첫 번째 주성분이 전체 데이터 분산의 약 80%를 설명하는 것으로 나타났는데, 이는 전 지구 해수면 변동이라는 강력한 공통 신호가 존재함을 시사한다 (Spratt and Lisiecki, 2016).
1.3. 지구물리학적 보정: 빙하 지각 평형 조정의 필수적 역할
해안가의 화석 산호나 심해 퇴적물 코어에서 얻은 기록은 그 자체로 과거의 ‘전 지구 평균 해수면’을 의미하지 않는다. 이는 지구의 지각이 고정된 것이 아니라, 그 위에 놓인 빙상과 해수의 무게 변화에 따라 변형되기 때문이다. 이 과정을 빙하 지각 평형 조정(GIA)이라고 부른다 (Lambeck et al., 2014).
빙하 지각 평형 조정(GIA)의 물리학적 원리는 다음과 같다. 빙하기에 수백만 톤의 얼음이 대륙을 누르면, 그 아래의 맨틀 물질이 옆으로 밀려나면서 지각이 가라앉는다. 동시에, 빙상 주변부의 지각은 맨틀 물질이 밀려옴에 따라 미세하게 융기하는데, 이를 ‘전방 팽창’이라 한다.
용어해설
전방 팽창(Forebulge): 매트리스 위에 무거운 볼링공을 놓으면 공 주변이 살짝 솟아오르는 것처럼, 거대한 빙상이 지각을 누를 때 그 주변부가 미세하게 융기하는 현상.
해빙기가 되어 빙상이 녹으면 정반대의 현상이 일어난다. 얼음의 무게가 사라진 지역의 지각은 천천히 다시 융기하고(지각 반등), 전방 팽창 지역이었던 곳(예: 북미 동부 해안, 바하마 군도)은 반대로 침강한다 (Lambeck et al., 2014). 따라서 특정 지역에서 관측된 과거의 해수면 높이(상대 해수면, Relative Sea Level, RSL)는 전 지구적인 해수량 변화뿐만 아니라 그 지역의 지각 운동까지 포함한 결과이다. 전 지구 평균 해수면 변화를 알기 위해서는 GIA 모델을 이용해 이러한 국지적인 지각 변동 효과를 반드시 보정해야 한다.
문제는 GIA 모델의 결과가 지구 내부 맨틀의 점성(viscosity)이나 과거 빙상의 역사(두께, 범위, 소멸 시기 등)와 같은 가정에 매우 민감하다는 점이다 (Lambeck et al., 2014). 그런데 과거 빙상의 역사는 우리가 바로 해수면 기록을 통해 알아내고자 하는 것이다. 이처럼 GIA 보정은 원인과 결과가 서로 맞물려 있는 복잡하고 반복적인 계산을 요구하는 도전적인 과제이다. 이 문제를 해결하기 위해 연구자들은 초기 빙상 모델로 GIA를 계산하고, 그 결과를 전 지구의 해수면 관측 자료와 비교한 뒤, 불일치하는 부분을 수정하여 빙상 모델을 개선하고 다시 GIA를 계산하는 과정을 반복한다. Lambeck et al. (2014)의 연구는 이러한 반복적 분석의 대표적인 사례로, 전 지구적 자료와 모델을 통합하여 최적의 해답을 찾아가는 과정을 보여준다.
아울러, GIA 외에 맨틀 유동에 의한 동적 지형(dynamic topography) 변화가 최후간빙기 지시자 고도에 수 m 규모의 체계적 편차를 유발할 수 있음이 보고되었다. 이는 정밀한 최후간빙기 최고해수면 추정에 GIA 보정뿐 아니라 장주기의 비등방(non-isostatic) 지각운동 효과도 함께 고려해야 함을 시사한다(Austermann et al., 2017).
1.4. 융빙수 기원 추적의 패러다임 전환: 해수면 지문의 물리학
GIA 개념에서 더 나아가, 최근 고해수면 연구의 패러다임을 바꾼 것이 ‘해수면 지문’ 이론이다 (Mitrovica et al., 2011). 이 이론은 융빙수의 기원에 따라 해수면 상승이 전 지구적으로 균일하게 일어나지 않는다는 점에 주목한다.
용어해설
해수면 지문(Sea-level Fingerprint): 거대한 빙상은 자체의 질량과 중력으로 주변 바닷물을 끌어당긴다. 만약 그린란드 빙상이 녹으면, 그 중력이 약해져 바닷물이 멀리 퍼져나간다. 그 결과, 역설적으로 그린란드 주변의 해수면은 오히려 하강하고, 지구 반대편의 해수면은 평균보다 더 많이 상승한다. 이처럼 녹아내린 빙상의 위치에 따라 전 세계 해수면이 다르게 변하는 독특한 패턴을 ‘해수면 지문’이라고 부른다.
핵심 물리학적 원리는 만유인력이다. 거대한 빙상은 질량을 가지고 주변의 해수를 끌어당긴다. 만약 그린란드 빙상이 녹아 질량을 잃으면, 그린란드 주변 해수를 끌어당기는 힘이 약해진다. 그 결과, 역설적으로 그린란드 바로 근처의 해수면은 하강하고, 그곳에서 멀리 떨어진 남반구와 같은 원격장의 해수면은 전 지구 평균보다 더 많이 상승하게 된다 (Mitrovica et al., 2011).
이러한 현상은 각각의 융빙수 기원지(예: 그린란드, 서남극, 로렌타이드 빙상)가 마치 사람의 지문처럼 고유한 공간적 해수면 변화 패턴을 만들어낸다는 것을 의미한다. 따라서 전 지구적으로 분포하는 여러 지점의 정밀한 상대 해수면 기록을 확보하고, 이를 각 기원지별로 예측되는 해수면 지문 패턴과 비교하면, 과거 특정 시점에 발생한 해수면 상승의 원인이 어느 빙상이었는지를 역으로 추적할 수 있다 (Lin et al., 2021). 이 기법은 4.2장에서 다룰 MWP-1A의 기원 논쟁을 해결하는 데 결정적인 역할을 했다.
이처럼 고해수면 과학의 발전은 관측 자료와 지구물리학적 모델링 사이의 긴밀하고 반복적인 상호작용을 통해 이루어져 왔다. 타히티에서 시추한 하나의 산호 코어 기록은 그 자체로는 단지 한 지점의 상대적인 수심 변화를 나타내는 숫자에 불과하다. 이 숫자가 전 지구적 의미를 갖기 위해서는 GIA 모델을 통한 보정이 필수적이다. 그러나 GIA 모델의 신뢰성은 그 모델이 타히티뿐만 아니라 전 세계 수백, 수천 개의 다른 지점에서 관측된 다양한 지질학적 증거들을 얼마나 잘 설명할 수 있는지에 달려 있다. 이처럼 개별 데이터의 해석은 전 지구적 모델에 의존하고, 전 지구적 모델의 타당성은 수많은 개별 데이터를 설명하는 능력에 의해 검증되는 상호의존적 관계야말로 이 분야 전체를 관통하는 핵심 원리이다.
2. 따뜻했던 과거: 미래의 바다를 비추는 거울
약 12만 9천 년 전부터 11만 6천 년 전까지 지속된 최후간빙기는 지구의 평균 기온이 산업화 이전보다 약간 높았던 시기로, 오늘날 우리가 직면한 온난화된 미래의 모습을 엿볼 수 있는 중요한 과거 유사 사례(paleo-analogue)로 여겨진다 (Dyer et al., 2021). 이 시기 해수면이 현재보다 얼마나 더 높았는지를 규명하는 것은 미래의 빙상 붕괴와 해수면 상승폭을 예측하는 데 있어 핵심적인 과제이다. 최근 연구들은 기존의 폭넓은 추정치를 넘어서, 혁신적인 방법론과 정밀한 연대 측정을 통해 최후간빙기 해수면 최고점 추정의 상한을 낮추는 방향으로 제약을 강화했으며, 기존의 5–10 m 관점과의 긴장 관계 속에서 추정치가 수렴 중임을 보여준다(IPCC, 2021; Dyer et al., 2021; Dumitru et al., 2023).
2.1. 기존의 통설과 그 한계
오랫동안 과학계의 통설은 최후간빙기 동안 전 지구 평균 해수면(GMSL)이 현재보다 ‘likely 5–10 m’ 더 높았다는 것이고(IPCC, 2021; Dutton et al., 2015), 이는 다수의 종합평가에 반영되어 왔다. 최근에는 아래에서 다룰 바하마의 공간 기울기 제약(Dyer et al., 2021)과 고정밀 U–Th 연대 측정(Dumitru et al., 2023)을 통해 상한값을 낮추는 결과(대략 1–3 m 범위)가 제시되고 있으나, 이는 방법론 및 가정에 따른 견해 차이가 남아 있어 합의가 정리되는 과정에 있다(IPCC, 2021; Dyer et al., 2021; Dumitru et al., 2023).
이처럼 폭넓은 추정 범위(5-10 m)는 상당한 불확실성을 내포하고 있었다. 가장 큰 불확실성의 원인은 앞서 언급된 빙하 지각 평형 조정(GIA) 및 기타 장기적인 지각 변동(예: 동적 지형)을 정확하게 보정하는 데 있었다 (Dutton et al., 2015). 서로 다른 지역에서 발견된 해수면 지시자들의 높이 편차를 설명하기 위해 다양한 GIA 모델이 사용되었지만, 어떤 모델이 가장 현실적인지에 대한 제약이 부족하여 해수면 최고점 추정치 역시 넓은 범위에 머물 수밖에 없었다.
2.2. 바하마의 지렛대: 공간적 기울기를 이용한 GIA 모델 검증
이러한 교착 상태를 돌파한 것은 Dyer et al. (2021)의 혁신적인 접근법이었다. 그들은 바하마 군도가 과거 북미 대륙을 덮었던 거대한 로렌타이드 빙상의 전방 팽창 지역에 가까이 위치했다는 점에 주목했다. GIA 모델에 따르면, 해빙기 동안 전방 팽창이 가라앉으면서 바하마 군도는 북쪽에서 남쪽으로 갈수록 침강의 정도가 달라지는, 즉 예측 가능한 공간적 기울기를 보여야 했다 (Dyer et al., 2021).
그들의 연구 방법은 이 가설을 직접 검증하는 것이었다. 먼저 바하마 군도 전역에 걸쳐 최후간빙기에 형성된 해안 지형의 높이를 고정밀 GPS 장비 등을 이용해 정밀하게 측정했다. 그런 다음, 맨틀 점성이나 빙상 붕괴 시나리오 등을 각기 다르게 가정한 576개의 방대한 GIA 모델 예측 결과와 이 관측 자료를 비교했다. 이 과정에서 실제로 관측된 북-남 간의 해수면 높이 기울기를 가장 잘 재현해내는 모델에 더 높은 ‘가중치’를 부여했다. 이는 지질학적 관측 자료를 이용해 수많은 GIA 모델 중 가장 현실적인 모델들을 선별해내는, 일종의 ‘자연 실험’과 같았다 (Dyer et al., 2021).
이러한 확률론적 접근을 통해 GIA 보정의 불확실성을 크게 줄인 결과, 최후간빙기 해수면 최고점에 대한 추정치는 기존 통설보다 현저히 낮아졌다. 그들의 결론은 최후간빙기 GMSL이 현재보다 최소 1.2 m 이상 높았을 확률이 95%이지만, 5.3 m를 초과했을 확률은 5%에 불과하다는 것이었다 (Dyer et al., 2021). 이는 최후간빙기 해수면 상승에 대한 기존의 인식을 바꾸는 중요한 결과였다. 이 연구는 단순히 전 세계에서 더 많은 데이터를 축적하는 것보다, 특정 지역이 가진 독특한 지구물리학적 이점을 전략적으로 활용하는 것이 과학적 난제를 해결하는 데 얼마나 강력한 도구가 될 수 있는지를 보여주었다. 바하마는 단순한 데이터 지점이 아니라, GIA 모델을 시험하고 검증하는 거대한 자연 실험실의 역할을 한 것이다.
2.3. 고정밀 연대 측정과 수정된 해수면 최고점
Dyer et al. (2021)의 연구가 GIA 모델의 불확실성을 줄이는 데 집중했다면, Dumitru et al. (2023)의 연구는 해수면 지시자의 연대 측정 불확실성을 줄임으로써 최후간빙기 해수면 복원의 정밀도를 한 단계 더 끌어올렸다. 이 연구팀은 바하마의 여러 섬에서 최후간빙기 산호 화석을 새로 채취하여 고정밀 U-Th 연대 측정을 수행했다 (Dumitru et al., 2023).
이렇게 얻어진 더 정확한 연대 자료를 Dyer 등의 연구를 통해 개선된 최신 GIA 모델과 결합하여 보정한 결과, 최후간빙기 GMSL 최고점에 대한 추정치는 더욱 좁은 범위로 수렴했다. 그들의 분석에 따르면, 최후간빙기 GMSL은 현재보다 1 m 이상 높았을 가능성이 높지만, 2.7 m를 초과했을 가능성은 매우 낮다 (Dumitru et al., 2023).
또한 이 연구는 해수면이 최고점에 도달한 기간 동안 비교적 안정적으로 유지되었음을 시사했다. 이는 그린란드와 남극 빙상의 융빙 기여가 동시에 일어난 것이 아니라, 북반구와 남반구의 빙상 붕괴가 서로 다른 시기에 발생했을 가능성(out-of-phase)을 제기한다. 이는 지구의 기후 시스템 내에서 양쪽 극지방이 어떻게 상호작용하는지에 대한 중요한 통찰을 제공한다 (Dumitru et al., 2023).
이러한 최근 연구 결과들은 최후간빙기를 기후 유사 사례로 해석하는 방식에 중요한 변화를 가져온다. 해수면 최고점이 기존에 생각했던 6-9 m(또는 5-10 m)가 아니라 1-3 m 범위에 가까울 수 있다는 제안은, 장기간의 온난화에 대한 남극 빙상의 반응이 이전에 가정했던 것보다 더 안정적이었을 가능성을 시사한다. 다만, 최후간빙기는 궤도 강제력(자연적)에 의한 장기간 평균의 산물인 반면, 현대는 온실가스(인위적)에 의한 빠른 가열이라는 경계 조건이 근본적으로 다르다. 따라서 이 시기의 기록을 미래 해수면 상승의 정량적 상한선으로 직접 외삽하기보다는, 기후-빙상 시스템의 물리적 과정과 민감도를 제약하는 데 활용하는 것이 타당하다(IPCC, 2021).
3. 얼음의 시대: 가장 낮았던 바다의 기준선
지구 기후 역사의 가장 극적인 장면 중 하나인 최후빙기절정기(약 26,000년 전 ~ 19,000년 전)는 북반구의 대부분이 거대한 대륙 빙상으로 뒤덮였던 시기이다. 이 시기 해수면이 얼마나 낮아졌는지를 정확히 아는 것은 당시의 기후와 환경을 이해하고, 이후에 이어진 극적인 해빙 과정의 총 규모를 정량화하는 데 있어 필수적인 기준점이 된다. 이 분야에서 기념비적인 연구로 평가받는 Lambeck et al. (2014)의 연구는 전 지구적 자료 합성을 통해 LGM 해수면 최저점에 대한 가장 신뢰도 높은 추정치를 제시했다.
3.1. 전 지구적 역산: 천 개의 증거를 종합하다
Lambeck et al. (2014) 연구의 핵심은 방대한 데이터의 종합과 정교한 모델링의 결합에 있다. 연구팀은 지난 35,000년간의 해수면 변화를 기록하고 있는 전 세계 원격장 지역의 고해수면 지시자 약 1,000개를 수집하여 광범위한 데이터베이스를 구축했다 (Lambeck et al., 2014).
그들은 이 방대하고 지리적으로 다양한 데이터 세트를 동시에 가장 잘 설명할 수 있는 빙상의 역사를 찾기 위해 반복적인 GIA 모델 역산 기법을 사용했다 (Lambeck et al., 2014). 이 접근법은 특정 지역의 국지적인 특성이나 데이터의 불확실성에 의해 결과가 왜곡되는 것을 방지한다. 개별 기록에 내재된 ‘잡음’은 전 지구적으로 분포된 수많은 데이터를 평균화하는 과정에서 상쇄되고, 그 속에 담긴 공통적인 신호, 즉 진정한 전 지구 평균 해수면 변화의 추세가 드러나게 된다. 이들의 연구는 특정 고품질 기록에 의존하기보다, 대규모 데이터셋의 통계적 힘을 이용하여 개별 기록의 한계를 극복할 수 있음을 보여준 대표적인 사례이다.
3.2. LGM 기준선: 해수면 -134 미터의 세계
이러한 종합적인 분석을 통해 도출된 핵심 결론은, 최후빙기절정기 동안 전 지구 평균 해수면이 약 29,000년 전부터 21,000년 전 사이에 현재보다 약 134 m 낮은 지점에 도달했다는 것이다 (Lambeck et al., 2014).
이 -134 m라는 수치는 기존에 널리 인용되던 -120 m ~ -125 m보다 상당히 낮은 값이다. 이러한 차이가 발생한 주된 이유는 Lambeck 연구팀이 GIA 모델링 과정에서 이전 연구들이 종종 간과했던 중요한 물리적 효과들을 정밀하게 고려했기 때문이다. 특히, 빙상이 녹아 해양으로 유입된 막대한 양의 물이 해저 지각을 눌러 침강시키는 효과, 즉 해양 하중(hydro-isostasy)까지 포함한 GIA 역산을 정밀하게 수행한 것이 결정적이었다(Lambeck et al., 2014).
해수면이 134 m 하강했다는 것은 당시 지구상에 현재보다 약 5,200만 km³에 달하는 막대한 양의 얼음이 육상에 빙상 형태로 존재했음을 의미한다 (Lambeck et al., 2014). 이 수치는 LGM 시기 기후 및 빙상 모델을 검증하고 초기 조건을 설정하는 데 있어 매우 중요한 경계 조건(boundary condition)으로 사용된다. 또한, 이 견고한 LGM 기준선은 이후에 이어질 해빙기 해수면 상승의 총량을 정의하는 출발점이 된다. 예를 들어, 다음 장에서 다룰 MWP-1A 시기에 약 15 m의 해수면 상승이 있었다는 사실은, 이 기준선을 통해 전체 해빙기 상승량(약 134 m)의 10% 이상이 불과 수백 년이라는 짧은 기간에 발생했음을 명확히 보여준다. 이처럼 Lambeck et al. (2014)의 연구는 이후 모든 해빙기 연구의 견고한 토대를 마련했다.
4. 거대한 해빙: 갑작스러운 해수면 상승의 미스터리
최후빙기절정기가 정점에 이른 후, 지구 기후는 기나긴 빙하기를 끝내고 오늘날과 같은 따뜻한 홀로세(Holocene)로 전환되는 극적인 변화를 겪었다. 이 거대한 해빙(Great Thaw) 과정에서 대륙을 뒤덮었던 빙상들이 녹아내리면서 해수면은 급격히 상승했다. 이 과정은 결코 순탄하거나 점진적이지 않았으며, 특히 ‘융빙수 펄스 1A’로 알려진 시기에는 인류 역사상 유례를 찾기 힘든 속도로 해수면이 치솟았다. 이 장에서는 해빙기 해수면 상승의 전반적인 양상과, 그중에서도 가장 격렬했던 사건인 MWP-1A의 규모, 시기, 그리고 그 기원을 둘러싼 치열한 과학적 논쟁과 그 해결 과정을 심층적으로 다룬다.
4.1. 해빙기 해수면 상승의 속도와 규모
Lambeck et al. (2014)에 따르면, 해빙기의 급격한 상승의 주요 국면은 대체로 16,500년 전부터 8,200년 전(16.5–8.2 ka)에 해당하며, 최종적인 수렴은 약 8,000년 전에서 6,000년 전(8–6 ka) 사이에 이루어진 것으로 본다(Lambeck et al., 2014). 이 기간 동안 총 해수면 상승량은 대략 120 m 규모이다. 문헌에서 자주 인용되는 ‘장기 평균(long-term average)’ 속도는 1,000년당 약 12 m(12 m/kyr)이지만(Lambeck et al., 2014), 16.5–8.2 ka 구간을 단순 산술 평균하면 약 14–15 m/kyr로 계산된다. 이처럼 평균 속도는 정의하는 구간에 따라 달라지며, 이러한 평균 속도는 과정의 극적인 비선형성을 가리고 있다.
이 비선형성의 정점을 보여준 것이 바로 Deschamps et al. (2012)의 연구이다. 그들은 타히티에서 시추한 고품질의 산호 화석 기록을 정밀하게 분석하여 MWP-1A의 실체를 명확히 밝혔다. 그들의 고해상도 연대 측정 결과, MWP-1A는 약 14,650년 전에 시작되어 14,310년 전 이전에 종료되었으며, 이는 북반구의 기온이 급격히 상승했던 ‘뵐링 온난기(Bølling warming)’와 정확히 일치하는 시점이다 (Deschamps et al., 2012).
더욱 놀라운 것은 상승의 규모와 속도였다. 이 불과 350년이 채 안 되는 기간 동안 해수면은 14 m에서 18 m 가량 상승했다. 이는 연평균 40 mm가 넘는 경이적인 상승 속도로, 현재 관측되는 해수면 상승률(연간 약 3-4 mm)의 10배를 훌쩍 뛰어넘는 수치이다 (Deschamps et al., 2012). 이러한 급격한 해수면 상승은 당시 해안가에 살던 인류에게는 재앙적인 사건이었을 것이며, 지구 기후 시스템 전체에도 엄청난 충격을 주었을 것이다.
4.2. 융빙수 기원 논쟁: MWP-1A의 급류를 추적하다
이처럼 막대한 양의 융빙수는 과연 어디에서 기원했을까? 이 질문은 고기후학계의 오랜 난제였다. 융빙수의 기원을 밝히는 것은 단순히 지질학적 호기심을 넘어, 지구 기후 시스템의 작동 방식을 이해하는 데 매우 중요하다. 예를 들어, 융빙수가 북대서양으로 유입되었다면 해양 순환에 직접적인 영향을 미쳐 북반구 기후를 급격히 변화시켰을 수 있고, 남극에서 기원했다면 그 영향은 또 다른 양상으로 나타났을 것이기 때문이다 (Deschamps et al., 2012).
4.2.1. 핵심 논쟁: 남극인가, 로렌타이드인가
MWP-1A의 기원을 둘러싼 논쟁은 크게 두 가설로 나뉘었다. 첫 번째는 남극 기원설이다. 남극 주변 해저 퇴적물 기록을 분석한 결과, MWP-1A와 거의 같은 시기에 남극 빙상에서 막대한 양의 빙산이 떨어져 나온 흔적(AID6, Antarctic Iceberg Discharge event 6)이 발견되었다 (Weber et al., 2014). 일부 초기 연구들도 남극의 상당한 기여 가능성을 시사했다 (Weaver et al., 2003; Golledge et al., 2014).
두 번째는 북미의 로렌타이드 빙상 기원설이다. 빙상 모델링 연구들은 당시 기온이 급상승하던 북대서양에 인접한 로렌타이드 빙상이 급격히 후퇴하며 막대한 융빙수를 방출했을 가능성을 제시했다 (Peltier, 2005).
이 논쟁이 오랫동안 해결되지 못한 근본적인 이유는 해수면 지문 분석 기법의 한계 때문이었다. 당시 고해상도 연대 자료를 확보할 수 있었던 지점은 타히티, 바베이도스, 순다 대륙붕 등 단 3곳에 불과했다. 이처럼 적은 수의 원격장 데이터만으로는 여러 가능한 융빙수 시나리오를 명확하게 구분해낼 수 없었고, 결과적으로 여러 해답이 가능한 ‘비유일성 문제’에 부딪혔다 (Lin et al., 2021).
4.2.2. 통합된 해법: 확장된 지리적 데이터셋의 힘
이 교착 상태를 해결한 것은 Lin et al. (2021)의 연구였다. 그들은 기존의 해수면 지문 분석 기법을 사용하되, 분석 대상을 6개의 지리적으로 잘 분포된 지점으로 확장함으로써 결정적인 돌파구를 마련했다 (Lin et al., 2021).
그들의 혁신은 단순히 데이터의 수를 늘린 데 있지 않았다. 가장 중요한 것은 과거 빙상에 가까웠던 ‘근접장’인 스코틀랜드와 같은 지역의 기록을 분석에 포함시킨 것이다. 이 전략적인 선택은 논쟁을 해결하는 ‘결정적 증거(smoking gun)’를 제공했다. 해수면 지문 이론에 따르면, 만약 남극 빙상이 대규모로 붕괴했다면 그 질량 감소로 인한 중력 효과 때문에 북반구에 위치한 스코틀랜드의 해수면은 하강해야 한다. 반대로, 북미의 로렌타이드 빙상이 붕괴했다면 스코틀랜드의 해수면은 (전 지구 평균보다는 적게) 상승해야 한다. 스코틀랜드의 지질학적 기록은 당시 해수면이 명백히 상승했음을 보여주었다. 이 관측은 남극이 MWP‑1A의 ‘지배적’ 기원이라는 시나리오를 배제하는 강한 제약을 제공하며, 남극의 기여가 0–35% 범위의 ‘비지배적’ 수준이었음을 지지한다(Lin et al., 2021).
이처럼 결정적인 제약을 포함한 데이터 기반 역산 모델링을 통해, Lin 연구팀은 MWP-1A의 기원에 대한 통합된 해법을 제시했다. 그들의 결론은 다음과 같다: MWP-1A 동안 총 17.9 m (95% 신뢰구간: 15.7–20.2 m)의 전 지구 평균 해수면 상승이 있었고, 기여도는 북미 빙상 12.0 m (5.6–15.4 m; 약 35–85%), 스칸디나비아 빙상 4.6 m (3.2–6.4 m; 약 10–35%), 그리고 남극 빙상 1.3 m (0–5.9 m; 0–35%)로 추정된다 (Lin et al., 2021). 이 결과는 기존의 해수면 기록뿐만 아니라, 현장 조사를 통해 복원된 각 지역의 빙상 후퇴 역사와도 잘 부합하는 것으로 나타났다 (Lin et al., 2021).
MWP-1A 기원 논쟁의 해결 과정은 과학적 논쟁이 종종 단 하나의 혁명적인 발견이 아니라, 기존 분석 기법을 보다 견고하고 전략적으로 선택된 데이터 기반에 적용함으로써 해결된다는 점을 보여준다. Lin 연구팀은 해수면 지문 분석이라는 기존의 도구를 사용했지만, 데이터의 지리적 범위를 확장함으로써 이전 연구들이 겪었던 비유일성 문제를 극복하고, 상충되던 여러 증거들을 하나로 통합하는 설명을 제시할 수 있었다. 이 결론은 지난 2만 년 중 가장 극적인 융빙 사건이 북반구의 급격한 온난화와 직접적으로 연결되어 있음을 명확히 함으로써, 기후 시스템의 ‘티핑 포인트’에 대한 가장 강력한 과거 유사 사례를 제공한다. 이는 21세기를 살아가는 우리에게 시사하는 바가 매우 크다.
5. 과거가 미래에게 보내는 경고
후기 제4기의 세 가지 핵심적인 시기—최후간빙기, 최후빙기절정기, 그리고 해빙기—에 대한 해수면 변동 연구의 최신 성과들을 종합적으로 검토했다. 이 과학적 여정을 통해 우리는 과거 해수면 복원에 대한 이해가 어떻게 발전해 왔는지, 그리고 과거의 기록이 미래를 예측하는 데 어떤 중요한 통찰을 제공하는지를 확인할 수 있었다.
연구의 흐름을 요약하면 다음과 같다. 첫째, Lambeck et al. (2014)은 전 지구적 데이터 합성과 GIA 모델링을 통해 최후빙기절정기의 해수면이 현재보다 약 134 m 낮았다는 견고한 기준선을 확립했다. 둘째, Deschamps et al. (2012)는 타히티의 산호 기록을 통해 해빙기 동안 발생한 MWP-1A가 연간 40 mm가 넘는 경이적인 속도로 해수면을 상승시킨 극단적인 사건이었음을 정량적으로 밝혔다. 셋째, Lin et al. (2021)은 확장된 지리적 데이터셋과 해수면 지문 분석을 통해 이 급격한 융빙수의 주된 기원이 북반구의 로렌타이드 및 스칸디나비아 빙상이었으며, 남극의 기여는 비지배적이었음을 규명했다. 마지막으로, Dyer et al. (2021)과 Dumitru et al. (2023)는 바하마 군도의 독특한 지구물리학적 특성을 활용하고 고정밀 연대 측정을 결합하여, 최후간빙기의 해수면 최고점이 기존의 5-10 m 범위보다 낮은 약 1-3 m 범위일 가능성을 제시하며 학계의 활발한 논의를 이끌고 있다.
이러한 연구 성과들을 관통하는 핵심적인 주제는 점점 더 정밀해지는 지질학적 데이터와 정교해지는 지구물리학적 모델 간의 반복적인 상호작용이 과학적 이해를 진전시키는 원동력이라는 점이다. 어떤 데이터나 모델도 단독으로는 완전한 해답을 줄 수 없다. 새로운 관측 데이터는 모델의 불확실성을 줄이는 제약 조건을 제공하고, 개선된 모델은 기존 데이터를 더 정확하게 해석할 수 있는 틀을 제공한다. 이러한 선순환 구조를 통해 고해수면 과학은 과거의 불확실성을 극복하고 점차 합의된 그림을 그려나가고 있다.
이러한 과거의 기록들은 오늘날 우리가 직면한 기후 변화와 미래 해수면 상승 문제에 시사점을 던진다.
최후간빙기 기록의 재평가: 해수면 최고점에 대한 추정 범위가 좁혀지는 것은 산업화 이전보다 1-2°C 정도 높은 온도가 장기간 지속될 때 빙상이 어떻게 반응하는지에 대한 중요한 기준점을 제공한다. 그러나 과거(자연적 궤도 변화)와 현재(인위적 온실가스)의 기후 강제력이 다르기 때문에, 이 기록을 미래 예측에 직접 외삽하는 데는 신중해야 하며, 빙상 모델의 민감도를 검증하고 개선하는 자료로 활용되어야 한다(IPCC, 2021).
융빙수 펄스 1A의 경고: 이 사건은 기후 시스템이 특정 임계점을 넘었을 때, 빙상이 점진적이 아닌 급격하고 비선형적인 방식으로 붕괴할 수 있음을 보여주는 가장 명백한 과거 사례이다 (Deschamps et al., 2012). 특히, 그 기원이 북반구의 급격한 온난화와 직접적으로 연결되었다는 점은 현재 진행 중인 북극 온난화가 그린란드 빙상에 어떤 잠재적 위험을 초래할 수 있는지에 대한 강력한 경고이다.
결론적으로, 과거 해수면 기록은 지구 시스템이 기후 변화에 얼마나 민감하고 역동적으로 반응할 수 있는지를 보여주는 자연의 실험 기록이다. 비록 과거의 기후 강제력(예: 지구 궤도 변화)이 현재의 인위적인 온실가스 증가와는 다르지만, 그 속에 담긴 물리적 과정과 피드백 메커니즘은 변하지 않는다. 따라서 과거 기록에 대한 지속적인 연구와 이해의 심화는 불확실성으로 가득한 미래 해수면 상승 시나리오의 범위를 좁히고, 전 세계 해안 지역 사회가 더 현명하고 시의적절한 대응책을 마련하는 데 필수적인 과학적 기반을 제공할 것이다 (IPCC, 2021).
참고문헌
Austermann, J., Mitrovica, J.X., Huybers, P., & Rovere, A. (2017). ‘Detection of a dynamic topography signal in last interglacial sea-level records’. Science Advances, 3(7), e1700457.
Camoin, G. F., et al. (2007). ‘IODP expedition 310 reconstructs the history of sea level, climate, and environment in the South Pacific since the last deglaciation’. Scientific Drilling, 5, pp. 4-12.
Deschamps, P., et al. (2012). ‘Ice-sheet collapse and sea-level rise at the Bølling warming 14,600 years ago’. Nature, 483(7391), pp. 559–564.
Dumitru, O. A., et al. (2023). ‘Last interglacial global mean sea level from high-precision U-series ages of Bahamian fossil coral reefs’. Quaternary Science Reviews, 318, 108287.
Dutton, A., et al. (2015). ‘Sea-level rise due to polar ice-sheet mass loss during past warm periods’. Science, 349(6244), aaa4019.
Dyer, B., et al. (2021). ‘Sea-level trends across The Bahamas constrain peak last interglacial ice melt’. Proceedings of the National Academy of Sciences, 118(33), e2026839118.
Edwards, R.L., Gallup, C.D. and Cheng, H. (2003) ‘Uranium-series dating of marine and lacustrine carbonates’. In Uranium-series Geochemistry, Reviews in Mineralogy and Geochemistry, 52(1), pp. 363–405.
Golledge, N.R., et al. (2014) ‘Antarctic contribution to meltwater pulse 1A from reduced Southern Ocean overturning’, Nature Communications, 5, 5107.
IPCC (2021). Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge University Press.
Lambeck, K., et al. (2014). ‘Sea level and global ice volumes from the Last Glacial Maximum to the Holocene’. Proceedings of the National Academy of Sciences, 111(43), pp. 15296–15303.
Lin, Y., et al. (2021). ‘A reconciled solution of Meltwater Pulse 1A sources using sea-level fingerprinting’. Nature Communications, 12(1), 2015.
Liu, J., et al. (2016) ‘Sea-level constraints on the amplitude and source distribution of Meltwater Pulse 1A’, Nature Geoscience, 9(2), pp.130–134.
Mitrovica, J. X., et al. (2011). ‘On the robustness of predictions of sea level fingerprints’. Geophysical Journal International, 187(2), pp. 729–742.
Peltier, W. R. (2005). ‘On the hemispheric origins of meltwater pulse 1a’. Quaternary Science Reviews, 24(14-15), pp. 1655-1671.
Rovere, A., et al. (2023) ‘The World Atlas of Last Interglacial Shorelines (version 1.0)’, Earth System Science Data, 15, pp. 1–23.
Spratt, R. M. and Lisiecki, L. E. (2016). ‘A Late Pleistocene sea level stack’. Climate of the Past, 12(4), pp. 1079–1092.
Weaver, A.J., et al. (2003) ‘Meltwater Pulse 1A from Antarctica as a trigger of the Bølling–Allerød warm interval’, Science, 299(5613), pp. 1709–1713.
Weber, M. E., et al. (2014). ‘Millennial-scale variability in Antarctic ice-sheet discharge during the last deglaciation’. Nature, 510(7503), pp. 134-138.