A stochastic weather generator for multisite, multivariate climate scenarios with transient climate change

技術要約：マルチサイト・多変数気候シナリオのための確率的気象生成器（非定常な気候変動を伴うもの）

問題提起
気候変動下における水資源の信頼できる評価には、多変数の依存関係、空間的な一貫性、および過渡的な気候シグナルを同時に保持する気候強制力データが必要である。既存の確率的気象生成器は、多くの場合、定常性を仮定しているか、あるいは簡略化された気候変動調整（例：静的なデルタ係数）を適用しており、これらはマルチサイトかつ高次元なアプリケーションにおけるこれらの特性を維持できていない。一般循環モデル（GCM）のアンサンブルは過渡的なシグナルを提供するものの、その規模は限定的であり（ロバストなストレス・テストを阻害する）、空間解像度は粗く、内部的な気候変動を捉えるには不十分であることが多い。一方、定常的な確率的生成器は、大規模なアンサンブルを生成して現実的な依存構造を再現できるが、非定常な温暖化傾向や降水変動の変化を表現することはできない。したがって、水文学的な影響モデリングのために、内部的な気候変動と進化する気候軌道を共に捉えることができる、マルチサイトかつ多変数の大規模アンサンブルを生成する統一されたアプローチが切実に求められている。

手法
著者らは、非パラメトリックなk近傍法（k-NN）再サンプリング手法と過渡的気候摂動モジュールを統合した確率的気象生成フレームワークを開発した。この手法は、週単位の気温および降水量の入力を必要とする、積雪主導型のシステムであるチリ中央部のマイポ川流域に対して適用されている。対象となる257の空間地点におけるデータを用いる。

本手法は、主に以下の2つのコンポーネントで構成される：

1. 定常的マルチサイトk-NN再サンプリング： このコンポーネントは、歴史的な観測値（1979–2019年）から再サンプリングすることで、合成された週単位の時系列を生成する。全257地点および両変数（気温と降水量）にわたって、現在の合成状態と統計的に類似した歴史的な状態を特定するために、完全な地点別状態一致距離指標を採用している。歴史的記録からサンプリングすることにより、このモジュールは、観測された空間依存性、時間構造、および変数間の相関関係を保持し、内部的な気候変動を表現する。
2. 過渡的気候変動摂動： このモジュールは、定常的なk-NNシーケンスを、GCM投影と整合する非定常な軌道へと変換する。バイアス補正およびダウンスケーリングされたCMIP6投影（4つの共通社会経済経路（SSP）にわたる9つのGCM）を用い、固定された1980–2010年のベースラインに対する、時間とともに進化する気候変動シグナルを算出する。
   • 気温シグナルは加法的に適用される。
   • 降水量シグナルは、スケール依存の変動を考慮して乗法的に適用される。
   • 将来の気候特性については、年次変動を平滑化しつつ長期的なトレンドを捉えるために、21年の移動窓が使用される。

最終的なアンサンブルは、各確率的k-NN実現値と特定のGCM–SSP強制経路を組み合わせることによって構築され、その結果、108の過渡的気候シナリオ（3つの確率的実現値 × 36のGCM–SSPの組み合わせ）が得られる。これらのシナリオは、WEAPモデルを強制し、流出量をシミュレートして干魃指数（SPI-12およびSRI-12）を算出するために使用される。

主な貢献

• 統一されたフレームワーク： 本研究は、依存構造のパラメトリックな形式を仮定することなく、内部的な気候変動（k-NN再サンプリングによる）と外部的な過渡的強制力（GCM摂動による）を結合する、スケーラブルな手法を提示している。
• 依存関係の保持： このアプローチは、非定常なトレンドを導入しながらも、257地点にわたる複雑な空間的一貫性と、変数間の相互依存関係（気温–降水量）を維持する。
• 過渡的干魃分析： 本フレームワークは、進化する気候軌道がどのように干魃リスクを変化させるかを評価することを可能にし、「定常的な仮定（固定された気候状態）」と「過渡的なシナリオ（進化するリスクプロファイル）」を区別する。
• 水文学的互換性： 本手法は、さらなる空間的または時間的な解像度変更を必要とせず、プロセスベースの水文学モデル（具体的にはWEAP）に直接適合する週単位のマルチサイトデータを作成する。

結果

• 統計的忠実度： 定常的なk-NN生成器は、気温と降水量の歴史的な季節周期、平均、および変動を密接に再現している。過渡的シナリオは、近未来（2020–2040年）および今世紀末（2080–2100年）の両方のホライゾンにおいて、ダウンスケーリングされたGCMアンサンブルの分布とトレンドを正常に追跡している。
• 空間および変数間相関： 生成されたシナリオは、GCMアンサンブルの空間相関構造を保持している。気温の相関は短距離では1.0に近く、距離が増すにつれてわずかに減少する。降水の相関は距離依存的な減衰を示すが、これは生成器によって適切に再現されている。地形勾配によって引き起こされる気温と降水量の負の相関も、世紀を通じて維持されている。
• 干魃への影響： SSP5-8.5シナリオの下では、過渡的な強制力は定常的な仮定と比較して、より深刻な干魃条件を生み出す。降水量の減少傾向は乾燥のシグナルを示すが、水文学的な応答（SRI-12）は上昇する気温によって著しく増幅される。今世紀末までに、過渡的シナリオと定常的シナリオの間のSRI-12の差は約−5に達し、降水量の減少と気温の上昇という複合的な影響が、降水量の減少単独の場合よりもはるかに深刻な水文学的干魃を引き起こすことを示している。

意義と主張
本論文は、このアプローチが、不確実性下における水文学的な影響モデリングのための、気候投影と確率的気象生成を結びつける計算効率の高いツールを提供すると主張している。内部変動と外部強制力を分離することで、本手法は、あり得る将来の気候の幅広い範囲を系統的に探索することを可能にする。これは、ストレス・テストを行う上で不可欠である。著者らは、インフラ計画において過渡的シナリオが重要であることを強調している。なぜなら、それらは単に離散的な将来の時点における変化の大きさだけでなく、気候の変化のタイミングやシーケンス（例：いつ適応が必要になるのかという具体的な年）を捉えることができるからである。本フレームワークは、深い不確実性に直面している複雑で管理された流域において、堅牢な戦略を評価するために必要な大規模かつ高次元のアンサンブルを生成するための実用的な解決策として提示されている。