アンサンブル天気予報

ほぼ同じ初期値から五十の予測を走らせ、それらが一致し、やがてバラバラになる様子を観察しよう。このページでは、それぞれの予測の価値、それらの平均の価値、そして予測のばらつきが事前に警告を与えていたかどうかを測定する。

インタラクティブシミュレーションを読み込んでいます...

レッスン

理論 — アンサンブル天気予報

アンサンブル予報は、開始状態が完全に把握されることは決してないため、わずかに異なる開始状態から同じモデルを何度も実行する。アンサンブル平均はそれらの実行結果の平均であり、最良の単一の推定値である。スプレッドは各メンバーがその平均からどれだけ離れているかを示し、それはその予測がどの程度信頼できるかという、予測自体からの声明である。ここで用いているモデルは Lorenz 96 である。エドワード・ローレンツが 1996 年にカオスのデモンストレーションとしてではなく、トイモデルの大気として書いたものだ。緯度円上の 40 の値が、それぞれ隣接する値からの移流によって押し流され、減衰し、一定の強制力によって駆動される。モデルのひとつの時間単位は五日として読み取られる。これはローレンツ自身の較正によるものだ。そしてひとつのモデル単位は温度異常の一ケルビンに相当する。これは閾値を具体的にするための、このページでの慣例である。

各記号の意味

N
アンサンブルが実行するメンバーの数。50 は運用上の一般的な数字であり、このページでは 2 から 100 までの間で変更できる
σ
スプレッド:自身の平均からの各メンバーの二乗平均平方根距離
E
誤差:予測と実際に起きたこととの間の二乗平均平方根差。事後にのみ計算可能となる
λ
小さな誤差の増大率。ここでは一日あたり 0.3154 であり、これは 2.20 日の倍増時間に相当する。実際の大気の振る舞いに近い数値である
P
メンバーから直接読み取られる確率:閾値を超えるメンバーの割合
前提
完璧なモデル。ここにあるすべてのメンバーは真実を生み出したのと同じ方程式を積分しているため、それらを分かつものは開始地点のみである。解析値(予報士による今日の最良の推定値)は、真実にガウシアンノイズを加えたものだ。そして各メンバーは、その解析値に同規模のノイズをもう一つ加えたものであるため、統計的に見れば真実もまたメンバーのひとりにすぎない。この最後の詳細こそが、確率に意味を持たせている。そうではなく真実そのものを摂動させると、真実が雲の中心に置かれてしまう。その方法で測定すると、15% の頻度で起きた出来事に対してアンサンブルは 24% と予測し、80% の頻度で起きた出来事に対しては 65% と予測する結果になった。
成り立たない場合
気象センターが頭を悩ませるすべての点。現実のアンサンブルは過小分散になりがちである。すべてのメンバーが同じ不完全な物理法則を抱えているため、スプレッドが誤差よりも小さくなってしまうのだ。これを修復するため、センターはモデルに摂動を加えることに多大な労力を費やしている。また、現実の摂動は無作為に抽出されるわけではない。五十のメンバーが誤差を単にサンプリングするのではなく探索できるように、最も速く成長する方向に沿って選ばれる。そして円上の 40 の数値は決して大気ではない。実際の予測モデルは 10⁹ のオーダーの数値を扱っている。

五十の予測の平均は、五十の予測すべてに勝る 🖖

ありふれた一週間を押して、16 日目の結果を見てほしい。典型的なメンバーの誤差は 4.35 K、それらの平均の誤差は 2.91 K であり、その平均に劣ったメンバーの数は 50 中 50 である。これは偶然の産物ではない。結果の表がその理由を示している。メンバー自身の平均に対する二乗誤差を展開すると交差項が消滅し、次の恒等式が残る。mean(Eᵢ²) = E(mean)² + (N−1)/N · spread²。メンバーの二乗誤差の平均は、常に、平均値の二乗誤差にメンバー間の不一致を足したものになる。これが、典型的なメンバーが平均を上回ることができない理由であり、単一のメンバーが上回るのが運でしかない理由でもある(穏やかな一週間ではひとつだけそうした例があった)。メンバーの意見が分かれるほど、その差は開く。16 日目には、両辺が 19.1308 = 19.1308 となる。

スプレッドこそが予測である 🖖

穏やかな一週間を押し、次に荒れ模様の一週間を押してほしい。同じモデル、同じ五十メンバーのアンサンブル、同じ 0.3 K の観測誤差である。穏やかな週には、アンサンブル平均は 15.7 日目 でもまだ発表する価値がある。一方、荒れ模様の週では 6.9 日 でその価値を失う。大気の予測可能性は週によって等しくはなく、重要なのは、正解を知る前に自分がどちらの週にいるのかを見分けられるという点だ。結果的に予測が外れるような日には、早い段階でメンバーの意見が分かれる。これは十分に信頼できる兆候である。第三パネルにプロットされた、六日目で検証済みの六十回の実行を見ると、スプレッドと誤差の相関は r = 0.80 となっている。これこそがアンサンブルの目的であり、「分からない」と述べる予測が、失敗しているのではなくむしろ自らの役割を果たしている理由である。

最良の単一の予測とは、決して起こり得ない天気の図である 🖖

リードタイムが延びるにつれて、ありふれた一週間の表の最後の列に注目してほしい。これは単一のメンバーの分散に対する、アンサンブル平均に残された分散の比率である。0 日目は 99%、8 日目は 57%、16 日目は 21% を示している。意見の異なる五十の図を平均するとそれらは平滑化されるが、平らな図は天気のパターンではない。このモデルに記録された 601 日のうち、最も滑らかな日であっても平均分散の 70% を保持しており、半分を下回る日は一日たりとも存在しない。しかし 16 日目の平均図が保持しているのは 21% である。つまり、アンサンブル平均は得られる最良の単一の推定値であると同時に、決して起こらない現象を描いた一枚の絵でもある。気象センターが平均の天気図ではなく降水確率を発表する理由はここにある。

全プロセスの詳細解説

  1. 木曜日に氷点下になるか? 7 ステップ

    今日は水曜日。来週の木曜日、つまり 8 日先に、場所 19 で氷点下になるかどうかを尋ねられている。アンサンブル平均は +1.44 K を示している。 何を発表するかを決定し、その後に何が起きたかを確認しよう。氷点下になるか?を押して、ツールをこの状態に設定してほしい。

    この問題をツールで開く

    1. まず平均を見てみよう。それが単一の予測から得られるものだからだ。8 日目の場所 19 におけるアンサンブル平均は +1.44 K であり、余裕でゼロを上回っている。モデルを一度しか実行していない予報士であれば、「霜は降りない」と発表して帰宅するだろう。

    2. 今度は平均する代わりに数えてみよう。50 のメンバーのうち、20 がその場所、その日で氷点下を示している。尋ねられた質問に対するアンサンブルの答えは 40% である。

    3. これらふたつの見解は異なる事象について述べており、平均をどれほど慎重に扱っても二番目の答えは得られない。平均値は、その基礎となる数値のうちいくつが基準線を下回っているかについて何も語らない。すべてのメンバーが +1.44 K を示す五十のメンバーと、−4.95 K から +11.25 K まで分布するここでの五十のメンバーは、どちらも平均すれば +1.44 になる。最初のアンサンブルは霜が降りる可能性を完全に否定するが、二番目のアンサンブルは十回に四回の可能性を与える。

    4. 8 日目の場所 19 で実際に起きたのは −3.86 K だった。平均は答えの符号を間違えていた。しかしアンサンブルは間違っていなかった。アンサンブルはその出来事に十回に四回の可能性があると述べており、十回に四回の出来事が実際に起きることは、予測の失敗ではないからだ。

    5. 先ほどの文章が妥当だと言えるのは、40%という数字が文字通り40%を意味する場合のみです。そして、たった一つの事例でそれを証明することはできません。実測データを示したのが四番目のパネルです。このリードタイムにおける1,600地点・日を、アンサンブルが提示した確率ごとに分類しました。30%〜40%の範囲にある130件の予測を見ると、実際に起きた割合は29%です。同様に60%〜70%の範囲では66%となっています。アンサンブルの確率をそのまま信頼してよいのは、それがモデルから出力されたものだからではありません。このようにデータと突き合わせて、しっかりと検証されているからです。

    6. 事例の選択については正直に言っておこう。これは平均の失敗を示すために選ばれた日であり、ほとんどの日において、天候が閾値のどちら側に転ぶかについて平均とアンサンブルは一致する。重要なのは、平均が通常間違っているということではない。平均はそのような形の質問には全く答えられないが、各メンバーなら答えられるということだ。

    7. このページの限界について述べておく。現実の霜注意報では、誤った「穏やかな日」のコストと誤った「警告」のコストを天秤にかける。そしてその計算にはふたつの間違いの価値が必要だが、大気はそれを教えてくれない。「40% は本当に 40% か」という疑問については、我々の calibration ページで適切に扱っている。

    解答

    「霜は降りない」ではなく、40% と発表しよう。平均は +1.44 K であり、50 のメンバー中 20 が氷点下を示し、実際の結果は −3.86 K だった。 平均は、誰も尋ねていない質問に答えていたのである。

学習の道すじ

天気予報の仕組み

参考文献 (5)