Global Warming Has Been Accelerating Since At Least 1990

技術的要約：地球温暖化は少なくとも1990年以降加速している

問題提起
近年の文献では、世界平均地上気温（GMST）が加速していることが報告されており、一部の研究ではこの加速は2015年頃に始まったと示唆されている。しかし、加速を検出するための既存の手法には限界がある。標準的な二次式トレンドの仕様（$xt = \alpha + \lambda t + \beta t^2 + \epsilon_t$）は、温暖化率が一定の割合で増加するという硬直的な関数形式を課しており、これは現実を反映していない可能性がある。さらに、二次モデルは、トレンド推定の前にエルニーニョや太陽効果をフィルタリングする二段階の手順を必要とすることが多く、これがパラメータ推定の精度を低下させる可能性がある。対照的に、区分線形トレンドは構造的変化を許容するものの、特定の時期における離散的なシフトを仮定するため、連続的に増加する速度を一連のブレイク（構造的変化）と誤認する可能性がある。本論文は、より柔軟な統計的枠組みを開発することで、超線形（線形よりも大きい）トレンドをテストし、加速の開始時期をより堅牢に特定することを目指している。

手法
著者らは、二次トレンドの柔軟な代替案として**線形対数仕様（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$）の制御変数を、二段階のフィルタリング・アプローチではなく、回帰内に直接含める一段階推定手順を採用している。モデルの推定には、系列相関と不均一分散を考慮するために、不均一分散および自己相関に頑健な（HAC）標準誤差（Newey-West推定量）を用いた最小二乗法（OLS）を用いている。

加速の時期を決定するために、著者らは1970–1990年のサンプルから開始し、終了日を2025年まで漸増的に延長していく拡張ウィンドウ・フレームワークを採用している。このアプローチに内在する多重テストの問題に対処するため、最大絶対値 $|t|$ ワイルド・ブートストラップ法を用いて、クリティカル値の算出とp値の調整を行い、偽陽性に対する堅牢性を確保している。

主な貢献

1. 新規の統計的仕様： 本論文は、加速の堅牢なテストとして線形対数仕様（$\beta \log(t)t$）を導入している。これは二次モデルとは異なり、温暖化率の一定の増加を仮定せず、線形項と二次項の間でしばしば見られる多重共線性問題も回避している。
2. 構造的変化に関する理論的洞察： 著者らは、真の基礎となるトレンドが超線形である場合、標準的な構造的変化テストは線形トレンドモデルの傾きの変化を最終的に検出するという理論的証明（定理 0.3）を提供している。これは、GMSTにおける報告された「構造的変化」（例：1970年代や2015年）が、離散的な温暖化率の変化を示す証拠ではなく、線形モデルに対して超線形トレンドを適合させたことによる統計的なアーティファクト（人工的な現象）である可能性を示唆している。
3. 加速時期の精緻化： ブートストラップ調整を伴う拡張ウィンドウ・フレームワークを適用することで、本研究は加速の開始時期を再評価し、それが単なる2015年以降の現象ではないという概念に異議を唱えている。

結果

• 線形対数仕様： 結果は、5つのデータセットすべてにおいて、線形対数仕様を用いた場合に加速（$\beta > 0$）の統計的に有意な証拠があることを示している。この有意性は、追加の制御変数がある場合とない場合の両方で保持されている。決定的なことに、拡張ウィンドウ分析によれば、加速の証拠は少なくとも1990年以降存在しており、より最近のデータ（2025年まで）が含まれるにつれて統計的有意性が強まっている。
• 二次仕様： 対照的に、二次仕様は、最も長い推定ウィンドウにおいて、かつ主に制御変数が含まれる場合にのみ有意な結果をもたらした。制御変数がない場合、二次係数が有意となったのは5つのデータセットのうち3つのみであった。
• 多項式指数の推定： 一般的な多項式トレンド（$t^\delta$）における指数 $\delta$ を推定する追加の演習は、真の指数が2（二次の場合）よりも小さく、かつ1よりも大きいことを示唆している。これは、観測されたダイナミクスに対して線形対数モデルがより適切な下限として機能することを支持している。
• 構造的変化： 理論的結果は、線形トレンドモデルにおける構造的変化の検出が、基礎となる超線形トレンドと整合していることを示唆しており、これまで報告されてきたGMSTトレンドのチェンジポイントに対する代替的な説明を提供している。

意義と主張
本論文は、地球温暖化の加速は単に2015年頃に始まった最近の現象ではなく、少なくとも1990年以降、数十年間にわたって存在しており、検出可能であることを主張している。著者らは、線形対数仕様が、硬直的な二次モデルと比較して、加速をテストするためのより堅牢で柔軟な枠組みを提供すると論じている。

本研究は、観測された加速が時間の経過とともに温暖化率の上昇を意味しており、それが生態系、人類の健康、および経済に影響を及ぼすことを強調している。著者らは、彼らの知見が、近年における人為的な炭素排出の空気中への割合の増加や、地球のエネルギー不均衡の倍増といった他の物理的な証拠と一致していると述べている。

最後に、本論文は、政策およびモデリングにおいて、超線形トレンドと構造的変化を区別することの重要性を強調している。もし真のトレンドが超線形的であれば、区分線形モデルは、連続的に増加する速度を捉えるのではなく、離散的な傾きの変化のみを捉えることで、GMSTの真のダイナミクスを過小評価する可能性がある。著者らは、加速の存在は緩和策および適応策の緊急性を際立たせるものであると結論づけている。