Global Warming Has Been Accelerating Since At Least 1990

기술 요약: 지구 온난화는 적어도 1990년 이후로 가속화되어 왔다

문제 제기
최근 문헌들은 지구 평균 지표 온도(GMST)가 가속화되었다고 보고해 왔으며, 일부 연구들은 이 가속화가 2015년경부터 시작되었다고 제안한다. 그러나 가속도를 탐지하는 기존 방법론들은 한계를 지니고 있다. 표준 이차 함수 추세 명세($xt = \alpha + \lambda t + \beta t^2 + \epsilon_t$)는 온난화율이 일정한 비율로 증가한다고 가정하는 경직된 함수 형태를 부과하며, 이는 현실을 반영하지 못할 수 있다. 또한, 이차 모델은 종종 데이터에서 엘니뇨 및 태양 효과를 필터링한 후 추세를 추정하는 2단계 절차를 요구하는데, 이는 매개변수 추정의 정밀도를 떨어뜨릴 수 있다. 반대로, 분절 선형(piecewise linear) 추세는 구조적 변화(structural breaks)를 허용하지만, 특정 날짜에서의 불연속적인 변화를 가정하므로, 지속적으로 증가하는 비율을 일련의 변화점으로 오인할 가능성이 있다. 본 논문은 보다 유연한 통계적 프레임워크를 개발하여 초선형(supralinear, 선형보다 큰) 추세를 테스트하고, 가속화의 시작점을 더욱 견고하게 결정함으로써 이러한 문제들을 해결하고자 한다.

방법론
저자들은 이차 추세의 유연한 대안으로서 **로그-선형 명세(linearithmic specification)**를 제안한다. 이 모델은 다음과 같이 정의된다:
$$x_t = \alpha + \beta \log(t)t + \epsilon_t$$
여기서 $x_t$는 시간 $t$에서의 GMST 편차이다. 가속화의 존재 여부는 $\beta > 0$이라는 가설을 평가함으로써 테스트된다. 이론적으로 로그-선형 항 $\log(t)t$는 선형 추세보다는 빠르게 성장하지만, 차수가 $\delta > 1$인 임의의 다항식 추세보다는 느리게 성장하므로, 초선형 가속화에 대한 하한선을 제공한다.

본 연구는 1970년 1월부터 2025년 12월까지의 기간을 다루는 5개의 주요 온도 데이터셋(HadCRUT, GISTEMP, NOAA, Berkeley Earth, ERA5)을 활용한다. 분석에는 엘니뇨 효과($N_t$)와 태양 복사 에너지($S_t$)에 대한 제어 변수를 2단계 필터링 방식이 아닌 회귀식 내에 직접 포함하는 단일 단계 추정 절차가 적용된다. 모델은 시계열 상관관계와 이분산성을 고려하기 위해 Newey-West 추정량인 이분산 및 자기상관 일치(HAC) 표준 오차를 사용한 최소제곱법(OLS)으로 추정된다.

가속화의 시점을 결정하기 위해, 저자들은 1970~1990년 샘플로 시작하여 종료 날짜를 2025년까지 점진적으로 확장하는 확장 윈도우(expanding window) 프레임워크를 사용한다. 이 접근 방식에 내재된 다중 테스트 문제를 해결하기 위해, 저자들은 max-|t| 와일드 붓스트랩(wild-bootstrap) 절차를 사용하여 임계값을 계산하고 p-값을 조정함으로써 위양성(false positives)에 대한 견고성을 확보한다.

주요 기여

1. 새로운 통계적 명세: 본 논문은 가속화를 위한 견고한 테스트로서 로그-선형 명세($\beta \log(t)t$)를 도입한다. 이차 모델과 달리, 이 모델은 온난화율의 일정한 증가율을 가정하지 않으며, 선형 항과 이차 항 사이에서 흔히 발견되는 다중공선성 문제를 피한다.
2. 구조적 변화에 대한 이론적 통찰: 저자들은 실제 기저 추세가 초선형일 경우, 표준 구조적 변화 테스트가 선형 추세 모델의 기울기 변화를 결국 감지하게 된다는 것을 입증하는 이론적 증명(Theorem 0.3)을 제공한다. 이는 GMST에서 보고된 "구조적 변화"(예: 1970년대 또는 2015년)가 불연속적인 온난화율의 변화가 아니라, 초선형 추세에 선형 모델을 적합시킬 때 나타나는 통계적 인위적 현상일 수 있음을 시사한다.
3. 가속화 시점의 정밀한 재평가: 붓스트랩 조정을 적용한 확장 윈도우 프레임워크을 통해, 본 연구는 가속화의 시작점을 재평가하며, 이것이 단지 2015년 이후의 현상만이 아니라는 점을 밝혀낸다.

결과

• 로그-선형 명세: 결과는 5개 데이터셋 모두에서 로그-선형 명식을 사용할 때 가속화에 대한 통계적으로 유의미한 증거($\beta > 0$)가 있음을 나타낸다. 이러한 유의성은 추가 제어 변수를 포함하거나 포함하지 않은 경우 모두에서 유지된다. 결정적으로, 확장 윈도우 분석에 따르면 가속화의 증거는 적어도 1990년 이후부터 존재하며, 최근 데이터(2025년까지)가 포함됨에 따라 통계적 유의성이 강화된다.
• 이차 명세: 이와 대조적으로, 이차 명식은 가장 긴 추정 윈도우를 사용하고 주로 추가 제어 변수를 포함했을 때만 유의미한 결과를 나타냈다. 제어 변수가 없는 경우, 이차 계수는 5개 데이터셋 중 3개에서만 유의미했다.
• 다항식 지수 추정: 일반적인 다항식 추세($t^\delta$)에서 지수 $\delta$를 추정하는 추가 실험은 실제 지수가 2(이차의 경우)보다 작지만 1보다는 크다는 것을 시사한다. 이는 관측된 역학에 대해 로그-선형 모델을 사용하는 것이 더 적절한 하한선임을 뒷받αυ한다.
• 구조적 변화: 이론적 결과는 선형 추세 모델에서의 구조적 변화 감지가 초선형 추세와 일치함을 시사하며, 이는 이전에 보고된 GMST 추세의 변화점들에 대한 대안적 설명을 제공한다.

의의 및 주장
본 논문은 지구 온난화 가속화가 단순히 2015년경에 시작된 최근의 현상이 아니라, 적어도 1990년 이후부터 수십 년 동안 존재해 왔으며 탐지 가능하다는 점을 주장한다. 저자들은 로그-선형 명식이 경직된 이차 모델에 비해 가속화를 테스트하는 데 있어 더 견고하고 유연한 프레임워크를 제공한다고 주장한다.

본 연구는 관측된 가속화가 시간 경과에 따른 온난화율의 증가를 의미하며, 이는 생태계, 인류 건강 및 경제에 시사하는 바가 크다는 점을 강조한다. 저자들은 이러한 결과가 인위적 탄소 배출의 대기 잔류 분율 증가 및 최근 수십 년간 지구 에너지 불균형의 두 배 증가와 같은 다른 물리적 증거들과 일치한다고 언급한다.

마지막으로, 본 논문은 정책 및 모델링에 있어 초선형 추세와 구조적 변화를 구분하는 것의 중요성을 강조한다. 만약 실제 추세가 초선형이라면, 분절 선형 모델은 불연속적인 기울기 변화만을 포착함으로써 GMST의 실제 역학을 과소평가할 수 있다. 저자들은 가속화의 존재가 완화 및 적응 전략의 시급성을 강조한다고 결론짓는다.