지난 간빙기 해수면 상승 전 지구적 지역적 분석

초록

지난 간빙기(LIG) 시기의 전 지구 평균 해수면(GSL)은 현재보다 69 m 높았을 가능성이 높다. 저자들은 전 세계 300여 개 이상의 지역 해수면 지표와 연대 오류를 포함한 데이터베이스를 구축하고, 가우시안 프로세스 회귀와 마코프 체인 몬테카를로(MCMC) 방법을 결합한 확률론적 모델을 적용했다. 결과는 지역별 해수면 변동성을 동시에 고려했을 때, 전 지구 해수면이 69 m 상승했을 확률이 95 % 이상임을 보여준다. 이는 IPCC가 제시한 4~6 m보다 높은 상승을 의미하며, 온난화 2 °C 수준에서도 심각한 해안 위험이 발생할 수 있음을 시사한다.

상세 분석

본 논문은 지난 간빙기(LIG) 시기의 전 지구 해수면(GSL) 변동을 정량적으로 재구성하기 위해 두 가지 핵심적인 통계적 기법을 도입하였다. 첫째, 저자들은 기존 연구에서 제시된 지역 해수면 지표를 넘어, 전 세계 30개 대륙판과 섬 지역을 포괄하는 350여 개의 지표를 수집·정제하였다. 각 지표는 고도, 퇴적 환경, 동위원소 기록 등 복합적인 지질학적 근거를 바탕으로 선정되었으며, 연대 오류는 방사성 탄소 연대, 우라늄·납 연대, 조석층서 연대 등 다양한 방법으로 추정된 불확실성을 포함한다.

둘째, 이러한 데이터에 가우시안 프로세스(GP) 회귀를 적용하여 공간적·시간적 연속성을 모델링하였다. GP는 관측된 지역 해수면과 잠재적인 전 지구 해수면 사이의 상관 구조를 커널 함수(주로 지오데식 거리와 시간 차이를 결합한 비정상 커널)로 정의함으로써, 관측이 부족한 지역에서도 합리적인 추정을 가능하게 한다. 동시에, 연대 오류를 반영하기 위해 마코프 체인 몬테카를로(MCMC) 샘플링을 수행하였다. MCMC는 각 지표의 연대 분포를 사전 확률로 설정하고, GP 모델과의 결합 가능도(likelihood)를 계산해 사후 분포를 추정한다. 이 과정에서 연대 상호 의존성(예: 동일한 퇴적층 내 다중 표본)과 측정 오류를 동시에 고려함으로써, 기존 연구에서 과소평가된 불확실성을 정량화하였다.

모델 검증 단계에서는 교차 검증(cross‑validation)과 후방 예측 검증(posterior predictive checks)을 수행했으며, 관측된 지역 해수면과 모델이 예측한 값 사이의 평균 절대 오차는 0.4 m 이하로, 기존 선형 보간법 대비 현저히 낮은 편차를 보였다. 또한, 민감도 분석을 통해 커널 파라미터와 연대 사전 분포의 변동이 최종 GSL 추정치에 미치는 영향을 평가했으며, 결과는 주요 결론이 파라미터 선택에 크게 의존하지 않음을 확인하였다.

최종 결과는 전 지구 평균 해수면이 6 m에서 9 m 사이에 머물렀을 확률이 95 % 이상이며, 지역별로는 고위도 대륙붕에서 +12 m, 저위도 열대 해역에서는 +4 m 정도의 차이를 보였다. 이러한 지역 변동성은 중력·탄성·지각 변형 효과를 반영한 지구 물리학적 모델과도 일치한다. 논문은 또한 IPCC 제5차 평가보고서가 제시한 4–6 m 상승 추정이 데이터와 모델링의 한계로 인해 보수적으로 낮게 설정된 것임을 지적한다.

전반적으로, 본 연구는 복합적인 지구 물리·지질 데이터와 최신 베이지안 통계 기법을 결합함으로써, LIG 시기의 해수면 변동을 보다 정밀하게 재구성하고, 향후 기후 변화 시나리오에 대한 해안 위험 평가에 중요한 기준점을 제공한다.