출처:
dataset v.1: Spratt, R.M. and Lisiecki, L.E. (2016)
참고자료:
[논문 요약]
빙하기 해수면 변화, 여러 조각 맞춰 하나의 큰 그림 완성하기
이 연구는 여러 과학자가 각기 다른 방법으로 추정한 과거의 해수면 높이 기록들을 하나로 합쳐, 훨씬 더 정확하고 믿을 만한 해수면 변화 그래프를 완성했다. 연구진은 통계 기법의 일종인 주성분 분석(Principal Component Analysis, PCA)을 이용해, 개별 기록의 오류는 상쇄하고 공통된 핵심 신호, 즉 지구 전체의 평균 해수면 변화만을 성공적으로 추출했다. 이 통합 기록을 통해, 과거 기후의 핵심 지표인 저서성 유공충 껍질의 산소동위원소비(델타 O-18c) 기록에 숨겨진 두 가지 주요 원인, 즉 빙하 부피(해수면)와 심해 온도의 영향력을 정량적으로 분리해냈다. 분석 결과, 장기적인 기후 변동에서 빙하 부피의 기여도는 약 45%에 그쳐, 심해 온도의 변화가 55%로 더 지배적인 요인이었음이 밝혀졌다. 또한, 해수면과 산소동위원소의 관계가 약 40만 년 전을 기점으로 단순한 선형 관계에서 복잡한 2차 함수 관계로 변화했다는 사실을 발견했다. 이는 과거 기후 시스템의 작동 방식이 시간에 따라 변했음을 의미한다. 따라서 이 연구가 제시한 새로운 해수면 스택은 과거 기후 모델의 정확성을 검증하는 새로운 표준이 되며, 빙상과 해양 온도의 복잡한 상호작용을 이해하는 데 결정적인 단서를 제공한다.
1. 서론: 왜 과거의 바다 높이를 알아야 할까?
과거 수십만 년 동안 지구에는 여러 차례의 빙하기와 따뜻한 간빙기가 반복되었다. 이 거대한 기후 변화의 역사를 이해하기 위해 과학자들은 과거의 바다 높이, 즉 지구 평균 해수면(Eustatic Sea Level)의 변화를 복원하는 데 힘써왔다. 여기서 ‘지구 평균 해수면’이란, 땅이 솟아오르거나 가라앉는 등 지역적인 영향은 모두 빼고, 순수하게 지구 전체의 바닷물 양이 늘거나 줄어서 생긴 평균적인 높이 변화를 뜻한다.
과거 해수면을 알아내는 데에는 주로 두 가지 프록시(Proxy)가 사용된다. ‘프록시’란, 과거의 환경을 직접 볼 수 없기 때문에, 당시 정보를 간직한 간접적인 증거(대리지표)를 말한다.
- 산호초 화석: 특정 시기의 바다 높이를 가장 정확하게 알려주지만, 연속적이지 않고 드문드문 발견된다는 단점이 있다.
- 해저 퇴적물: 바다 밑에 쌓인 흙에서 발견되는 아주 작은 생물(유공충)의 껍데기 화석을 분석하는 방법이다. 연속적인 기록을 얻을 수 있지만, 여러 가지 불확실한 가정이 필요하고 지역적인 바다 환경에 영향을 많이 받는다는 한계가 있다.
이 연구의 목표는 바로 이 ‘신호 대 잡음 비율(signal-to-noise ratio)’ 문제를 해결하는 것이다. 우리가 찾고자 하는 전 지구적 해수면 변화라는 ‘신호(Signal)’가 다양한 국지적 요인과 방법론적 한계라는 ‘잡음(Noise)’에 가려져 있기 때문이다. 각기 다른 종류의 잡음을 가진 여러 기록들을 통계적으로 결합함으로써, 개별 기록에 공통적으로 존재하는 진정한 신호를 분리하고 증폭시켜 신호 대 잡음 비율을 획기적으로 개선하는 것이 이 연구의 핵심이다.
2. 연구 방법: 여러 기록을 하나로 합치는 기술, ‘스택’
연구진은 여러 개의 해수면 기록을 합치기 위해 ‘스택(Stack)’ 과 ‘주성분 분석(PCA)’ 이라는 두 가지 핵심 기술을 사용했다.
스택(Stack)이란, 여러 개의 데이터 기록을 말 그대로 차곡차곡 겹쳐서(stack) 평균을 내는 것을 말한다. 마치 여러 장의 울퉁불퉁한 종이를 겹쳐놓으면 전체적으로는 평평해지는 것처럼, 개별 기록이 가진 자잘한 오류(노이즈)는 서로 상쇄되고 공통된 진짜 신호만 뚜렷하게 남게 된다.
주성분 분석(Principal Component Analysis, PCA)은 이 ‘스택’ 작업을 더 정교하게 해주는 통계 도구이다. 여러 정보가 뒤섞인 복잡한 데이터 속에서 가장 중요하고 핵심적인 패턴(주성분)을 찾아내는 기술이다. 예를 들어, 여러 사람이 동시에 말하는 시끄러운 녹음 파일에서 가장 목소리가 큰 사람의 목소리만 골라내는 것과 비슷하다. 이 연구에서는 여러 기록에 공통적으로 들어있는 ‘지구 평균 해수면 변화’라는 가장 큰 목소리를 찾아낸 것이다.
연구 과정은 다음과 같다.
1. 자료 수집: 기존 연구에서 만들어진 7개의 해수면 변화 기록을 수집했다.
2. 시간 맞추기: 각기 다른 시간표를 가진 기록들을 하나의 표준 시간표에 맞추는 작업이 필요했다. 이때 기준이 된 것이 바로 LR04 저서성 유공충 산소동위원소비 스택이다. 이 어려운 용어는 아래와 같이 풀어서 이해할 수 있다.
- 유공충(Foraminifera): 바다에 사는 아주 작은 단세포 생물로, 죽으면 탄산칼슘 껍데기가 해저에 쌓여 화석이 된다.
- 저서성(Benthic): 이 유공충이 바다 깊은 곳, 바닥 근처에 살았다는 뜻이다.
- 산소동위원소비(델타 O-18): 산소에는 보통의 가벼운 산소(O-16)와 아주 약간 더 무거운 산소(O-18)가 있다. 유공충은 껍데기를 만들 때 주변 바닷물 속 산소를 쓰는데, 이때 두 산소의 비율은 당시의 수온과 지구 전체의 빙하 양에 따라 결정된다. 빙하가 많아져 해수면이 낮아지면, 가벼운 산소(O-16)가 대륙 빙하에 갇히기 때문에 바닷물에는 무거운 산소(O-18) 비율이 높아져 산소동위원소비(델타 O-18) 값이 커진다. 또 수온이 낮아져도 이 값은 커진다. 즉, 이 값은 과거 기후의 ‘성적표’와 같다.
- LR04 스택: 전 세계 57곳의 해저 퇴적물에서 얻은 이 산소동위원소비 기록들을 평균 내어 만든, 지난 530만 년간의 기후 변화를 보여주는 ‘표준 시간표’이다.
3. PCA 분석 및 스택 생성: LR04 시간표에 맞춰진 7개 기록에 PCA를 적용했다. 그 결과, 첫 번째 주성분(PC1)이 전체 변화의 약 80%를 설명하는 가장 강력한 공통 신호로 나타났고, 이를 새로운 해수면 스택으로 삼았다.
4. 실제 높이로 변환: 이 상대적인 변화 그래프를 실제 미터(m) 단위로 바꾸기 위해, 과학적으로 잘 알려진 두 시점의 해수면 높이를 기준으로 사용했다. 하나는 마지막 빙하기 중 가장 추웠던 최종 빙기 절정기(Last Glacial Maximum, LGM)(약 2만 4천 년 전)의 해수면(-130m)이고, 다른 하나는 현재와 거의 같은 5천 년 전의 해수면(0m)이다.
연구진은 통계 처리의 엄밀함을 더하기 위해 추가적인 장치를 마련했다. PCA 분석 전, 각 기록의 평균을 0, 표준편차를 1로 맞추는 ‘정규화(Normalization)’ 과정을 거쳐 특정 기록이 분석에 과도한 영향을 미치는 것을 방지했다. 또한, 최종 결과의 불확실성을 계산할 때, 단순히 오차 범위를 평균 내지 않고 각 기록에 1–2천 년의 무작위 시간 이동을 적용하는 시뮬레이션을 포함했다. 이는 고기후 연구에서 가장 큰 오차 원인 중 하나인 ‘연대 모델의 불확실성’까지 고려한 매우 정교한 접근 방식이다.
그림 1: 분석에 사용된 7개의 해수면 기록
그림 1은 이 연구에서 종합한 7개의 기존 해수면 복원 기록을 하나의 그래프에 모아 보여준다.
- X축 (가로축): 시간 (0 ~ 70만 년 전)
- Y축 (세로축): 해수면 높이 (현대 대비 미터)
- 내용: 각기 다른 색상의 선들이 7개의 개별 기록을 나타낸다. 모든 선들이 큰 틀에서는 함께 오르내리지만(빙하기-간빙기 주기), 세부적인 높이나 시점에서 조금씩 차이를 보이는 것을 확인할 수 있다. 연구의 목표는 이 선들의 평균적인 공통 신호를 찾는 것이다.
표 1: 주성분 분석(PCA) 로딩 값
‘로딩 값’은 각 기록이 주성분(공통 신호)을 만드는 데 얼마나 기여했는지를 보여주는 수치이다. 첫 번째 주성분(PC1)이 전체 변화량의 약 80%를 설명하며, 모든 기록이 PC1에 대해 0.34~0.48 사이의 비슷한 로딩 값을 갖는다. 이는 모든 기록이 공통 신호에 고르게 기여했음을 의미한다.
| 기록 출처 | PC1 (Short, 83%) | PC1 (Long, 77%) |
| Bintanja et al., 2005 | 0.40 | 0.48 |
| Elderfield et al., 2012 | 0.34 | 0.44 |
| Shakun et al., 2015 | 0.37 | 0.45 |
| Rohling et al., 2014 | 0.38 | 0.45 |
| Sosdian and Rosenthal, 2009 | 0.35 | 0.42 |
| Waelbroeck et al., 2002 | 0.40 | – |
| Rohling et al., 2009 | 0.40 | – |
3. 주요 연구 결과
3.1. 해수면 스택의 정확도
이렇게 만들어진 새로운 해수면 변화 그래프는 기존의 가장 확실한 증거인 산호초 기록과 놀라울 정도로 잘 들어맞았다. 특히 과거의 주요 간빙기였던 해양 동위원소 스테이지(Marine Isotope Stage, MIS) 5e와 MIS 11 시기의 최고 해수면 높이가 산호초 기록과 불과 5m 이내의 차이만 보였다. ‘해양 동위원소 스테이지(MIS)’란, 위에서 설명한 산소동위원소비 표준 시간표를 기준으로 지구 기후의 역사를 따뜻했던 시기(간빙기, 홀수 번호)와 추웠던 시기(빙하기, 짝수 번호)로 나눈 구간을 말한다. 예를 들어 MIS 1은 현재 우리가 사는 따뜻한 시대, MIS 2는 바로 이전의 마지막 빙하기를 뜻한다.
그림 2: 완성된 해수면 스택과 정확도 검증
그림 2는 세 개의 패널로 구성되어 이 연구의 최종 결과물을 보여준다.
- (a) 패널: 이 연구에서 만든 해수면 스택(파란 선)과 기후 표준 시간표(LR04, 검은 점선)를 비교한다. 두 기록의 오르내림이 거의 일치하는 것을 보여준다.
- (b) 패널: PCA로 만든 해수면 스택(초록 선)과 단순히 모든 기록을 산술 평균한 값(검은 선)을 비교한다. 두 결과가 매우 유사함을 보여주어 방법의 안정성을 입증한다.
- (c) 패널: 최종 해수면 스택(검은 선)과 통계적 오차 분석을 통해 계산한 신뢰구간(여러 색의 띠)을 함께 보여준다. 최종 결과가 대부분의 신뢰구간 안에 안정적으로 위치함을 보여준다.
표 2 & 3: 해수면 최고점/최저점 비교
아래 표는 주요 빙하기/간빙기 때의 해수면 높이를 각 기록과 비교한 것이다. 개별 기록은 추정치가 제각각이지만, 이들을 평균 내고 이 연구의 스택(표 3의 ‘Scaled PC1’)과 비교하면 산호초 기록(‘GIA-corrected’)과 매우 유사한 값을 보인다. 이는 스택이 개별 기록보다 더 정확함을 보여준다.
| 해양동위원소 스테이지 | 평균(개별기록) | GIA-보정(산호초) | Scaled PC1 (이 연구) |
| 2 (LGM) | -123 m | -130 ~ -134 m | -130 m |
| 5e (간빙기) | 7 m | 6 ~ 9 m | 3 m |
| 11 (간빙기) | 18 m | 6 ~ 13 m | 16 m |
3.2. 빙하와 수온, 누가 기후 변화를 이끌었나?
앞서 산소동위원소비(델타 O-18) 값은 ‘빙하 양’과 ‘심해 수온’ 두 가지 요인이 합쳐져 결정된다고 설명했다. 이 연구는 해수면 스택(빙하 양)을 산소동위원소비 기록과 비교하여, 두 요인의 영향력을 분리해냈다. 그 결과, 지난 수십만 년의 긴 시간 동안 산소동위원소비 변화의 원인 중 빙하 양의 변화가 차지하는 비중은 약 45%였고, 나머지 55%는 심해 수온의 변화 때문이었다. 이는 빙하기-간빙기 순환이 단순히 빙하가 커졌다 작아졌다 하는 현상이 아니라, 그에 못지않게 깊은 바다의 온도가 크게 오르내리는 과정이 함께 작용했음을 의미한다.
그림 3: 스펙트럼 분석 결과
이 그래프는 데이터의 변화가 어떤 주기로 반복되는지를 보여준다.
- X축: 주기 (ka⁻¹, 숫자가 작을수록 긴 주기)
- Y축: 스펙트럼 파워 (변화의 강도)
- 내용: 파란 선(해수면 스택)과 검은 점선(산소동위원소비) 모두 약 10만 년(0.01 지점)과 4만 1천 년(약024 지점) 주기에서 강한 피크를 보인다. 하지만 모든 주기에서 검은 점선의 피크가 파란 선보다 훨씬 높다. 이는 산소동위원소비 기록 안에 해수면 변화 외에 다른 요인(즉, 심해 수온)이 크게 작용하고 있음을 시각적으로 보여준다.
3.3. 약 40만 년 전, 기후 시스템에 무슨 일이?
매우 흥미롭게도, 해수면과 산소동위원소비의 관계가 약 40만 년 전을 기점으로 달라진 것이 발견되었다. 40만 년 전보다 과거에는 두 요소가 비교적 단순한 선형 관계(하나가 커지면 다른 하나도 일정 비율로 커지는 관계)를 보였다. 하지만 40만 년 전 이후부터 현재까지는 더 복잡한 2차 함수 관계로 바뀌었다. 이는 MIS 11이라는 유난히 길고 따뜻했던 간빙기를 겪으면서 지구의 기후 시스템 작동 방식에 근본적인 변화가 생겼을 가능성을 시사한다.
그림 4: 해수면과 산소동위원소비의 관계 변화
이 그림은 해수면과 산소동위원소비의 관계가 시간에 따라 변했음을 보여준다.
- (a) 패널: 실제 해수면 변화(파란 선)를 산소동위원소비 기록만으로 예측한 결과이다. 40만 년 이전(오른쪽)은 선형 모델(검은 선)이 잘 맞지만, 40만 년 이후(왼쪽)는 2차 함수 모델(붉은 선)이 더 잘 맞는다.
- (b), (c) 패널: 각 구간의 데이터를 점으로 찍어 보여주는 산점도이다. (b) 0~40만 년 전 구간은 데이터가 붉은 곡선(2차 함수)을 따라 분포하고, (c) 40만~80만 년 전 구간은 데이터가 파란 직선(선형)을 따라 분포하여 관계가 변했음을 명확히 보여준다.
3.4. 공통 신호 외의 차이점
PCA 분석은 공통된 신호(PC1) 외에, 기록들 간의 차이점(PC2, PC3)도 보여주었다. 두 번째 주성분(PC2)은 대서양과 태평양 깊은 바다의 특성 차이를, 세 번째 주성분(PC3)은 바다의 표층과 심층 사이의 차이를 주로 반영했다. 이는 지구 전체의 해수면 변화라는 큰 흐름 속에서도 각 지역과 수심에 따라 독자적인 변화가 있었음을 알려준다.
그림 5: 기록들 간의 차이 (PC2, PC3)
이 그래프는 PCA 분석의 부가적인 결과물이다.
- (a) 패널: 두 번째 주성분(PC2)의 시간에 따른 변화를 보여준다. 이는 주로 대서양과 태평양 심층수 사이의 차이를 반영한다.
- (b) 패널: 세 번째 주성분(PC3)의 변화를 보여준다. 이는 바다 표층과 심층 사이의 차이를 반영한다.
4. 결론
이 연구는 통계적 통합 기법을 통해, 개별 연구 결과들보다 훨씬 정확하고 신뢰도 높은 후기 플라이토세의 표준 해수면 변화 기록을 만들어냈다. 이 해수면 스택은 과거 기후를 재현하는 컴퓨터 시뮬레이션이 얼마나 정확한지 검증하는 ‘정답지’ 역할을 할 수 있다.
나아가 이 연구는 고기후학 분야의 중요한 방법론적 전환을 보여준다. 과거에는 새로운 퇴적물 코어를 시추하여 개별 기록을 생성하는 데 초점이 맞춰졌다면, 이 연구는 기존에 축적된 방대한 데이터를 정교한 통계적 방법으로 ‘종합’하여 새로운 과학적 통찰을 이끌어내는 방향으로의 전환을 상징한다. 이는 개별적인 목소리를 추가하는 것을 넘어, 기존 기록들의 합창 속에서 공통된 멜로디를 찾아내는 것의 중요성을 일깨워 준다.
이 연구의 중요성과 신뢰성은 2025년에 저자들이 ‘버전 2’ 데이터 세트를 공개함으로써 다시 한번 확인되었다. 이는 새로운 과학적 결론을 담은 후속 연구가 아니라, 기존 2016년 연구의 데이터 세트를 업데이트한 것이다. 이처럼 지속적인 데이터 관리는 이 해수면 스택이 고기후학 연구 커뮤니티에서 중요한 표준 자료로 널리 사용되고 있음을 보여주는 증거이다.
결론적으로, 이 연구는 지구의 빙하 양과 심해 온도가 서로 복잡하게 영향을 주고받았으며, 그 관계가 약 40만 년 전을 기점으로 변화했다는 중요한 사실을 밝혀내어 미래 기후 변화를 이해하는 데 새로운 단서를 제공했다.
[부록] 데이터 세트 버전 1(2016) vs. 버전 2(2025) 비교
핵심적인 변화는 단순한 오류 수정을 넘어, 데이터의 정식화, 투명성, 재사용성을 획기적으로 높인 것이다. 특히 모든 시점의 불확실성 값을 데이터화하여 제공한 것이 가장 큰 개선점이다.
| 항목 | 이전 버전 (2015년 PANGAEA / 2016년 논문) | 버전 2 (2025년 PANGAEA) |
| 자료의 기준 | 토론 논문(Discussion paper) 부록 기반의 비정식 데이터 (2015) | 최종 게재 논문(2016) 기준의 공식 데이터 |
| 불확실성 정보 | 논문 본문에만 평균값(1σ 약 9–12 m)으로 서술됨. 데이터 파일에는 불확실성 값 없음. | 모든 데이터 포인트마다 불확실성 열(column) 추가: • 1σ 표준편차 • 95% 신뢰구간 하한·상한(CI low/up) |
| 데이터 구조 | ‘합성 스택’ 열이 포함되어 있었음. | ‘합성 스택’ 열 삭제. 단기(0-430 ka)와 장기(0-798 ka) 스택을 명확히 분리하여 제공. |
| 메타데이터 및 라이선스 | 라이선스 불명확. 큐레이션 수준 낮음. | CC-BY 3.0 라이선스 명시. 향상된 큐레이션(Enhanced curation) 적용. NOAA 원자료와 상호 링크됨. |
1. 활용성: 불확실성 정보의 전면 데이터화 (가장 큰 변화)
이전에는 연구자가 불확실성을 고려하려면 논문에 서술된 평균치(약 9~12m)를 일괄적으로 적용해야 했다. 하지만 v2에서는 모든 시점(1 ka 간격)마다 표준편차와 95% 신뢰구간 값을 별도의 열로 제공한다. 이로써 다른 연구자들은 특정 시점의 오차를 정확히 알 수 있고, 자신의 모델에 오차 전파(error propagation)를 정량적으로 구현하는 등 고차원적인 분석을 수행할 수 있게 되었다.
2. 신뢰성: 데이터의 공식화 및 표준화
v2는 동료 심사를 거친 최종 논문(2016)을 기준으로 데이터를 재정리한 공식 버전이다. 반면, 이전의 2015년 PANGAEA 데이터는 토론 논문 단계의 예비 자료였다. 또한 v2는 데이터 라이선스(CC-BY 3.0)를 명확히 하고, PANGAEA의 ‘향상된 큐레이션’ 기준에 맞춰 파라미터와 단위를 표준화하여 데이터의 신뢰성과 장기 보존성을 크게 높였다.
3. 투명성: ‘합성 스택’ 열의 제거
이전 버전에 있던 ‘합성 스택’ 열은 사용이 편리했지만, 서로 다른 수의 기록(단기 7개 vs 장기 5개)으로 만들어진 데이터를 합쳐놓아 통계적 특성이 혼재될 여지가 있었다. v2에서는 이 열을 제거하고, 단기 스택과 장기 스택을 각각의 불확실성 정보와 함께 투명하게 분리했다. 이제 사용자는 논문의 방법론에 따라 직접 두 스택을 조합하여 사용해야 하며, 이 과정에서 각 데이터의 근거와 한계를 더 명확히 인지할 수 있게 되었다.