レッスン
理論 — アンサンブル天気予報
アンサンブル予報は、開始状態が完全に把握されることは決してないため、わずかに異なる開始状態から同じモデルを何度も実行する。アンサンブル平均はそれらの実行結果の平均であり、最良の単一の推定値である。スプレッドは各メンバーがその平均からどれだけ離れているかを示し、それはその予測がどの程度信頼できるかという、予測自体からの声明である。ここで用いているモデルは Lorenz 96 である。エドワード・ローレンツが 1996 年にカオスのデモンストレーションとしてではなく、トイモデルの大気として書いたものだ。緯度円上の 40 の値が、それぞれ隣接する値からの移流によって押し流され、減衰し、一定の強制力によって駆動される。モデルのひとつの時間単位は五日として読み取られる。これはローレンツ自身の較正によるものだ。そしてひとつのモデル単位は温度異常の一ケルビンに相当する。これは閾値を具体的にするための、このページでの慣例である。
各記号の意味
N- アンサンブルが実行するメンバーの数。50 は運用上の一般的な数字であり、このページでは 2 から 100 までの間で変更できる
σ- スプレッド:自身の平均からの各メンバーの二乗平均平方根距離
E- 誤差:予測と実際に起きたこととの間の二乗平均平方根差。事後にのみ計算可能となる
λ- 小さな誤差の増大率。ここでは一日あたり 0.3154 であり、これは 2.20 日の倍増時間に相当する。実際の大気の振る舞いに近い数値である
P- メンバーから直接読み取られる確率:閾値を超えるメンバーの割合
- 前提
- 完璧なモデル。ここにあるすべてのメンバーは真実を生み出したのと同じ方程式を積分しているため、それらを分かつものは開始地点のみである。解析値(予報士による今日の最良の推定値)は、真実にガウシアンノイズを加えたものだ。そして各メンバーは、その解析値に同規模のノイズをもう一つ加えたものであるため、統計的に見れば真実もまたメンバーのひとりにすぎない。この最後の詳細こそが、確率に意味を持たせている。そうではなく真実そのものを摂動させると、真実が雲の中心に置かれてしまう。その方法で測定すると、15% の頻度で起きた出来事に対してアンサンブルは 24% と予測し、80% の頻度で起きた出来事に対しては 65% と予測する結果になった。
- 成り立たない場合
- 気象センターが頭を悩ませるすべての点。現実のアンサンブルは過小分散になりがちである。すべてのメンバーが同じ不完全な物理法則を抱えているため、スプレッドが誤差よりも小さくなってしまうのだ。これを修復するため、センターはモデルに摂動を加えることに多大な労力を費やしている。また、現実の摂動は無作為に抽出されるわけではない。五十のメンバーが誤差を単にサンプリングするのではなく探索できるように、最も速く成長する方向に沿って選ばれる。そして円上の 40 の数値は決して大気ではない。実際の予測モデルは 10⁹ のオーダーの数値を扱っている。
全プロセスの詳細解説
-
木曜日に氷点下になるか? 7 ステップ
今日は水曜日。来週の木曜日、つまり 8 日先に、場所 19 で氷点下になるかどうかを尋ねられている。アンサンブル平均は +1.44 K を示している。 何を発表するかを決定し、その後に何が起きたかを確認しよう。氷点下になるか?を押して、ツールをこの状態に設定してほしい。
-
まず平均を見てみよう。それが単一の予測から得られるものだからだ。8 日目の場所 19 におけるアンサンブル平均は +1.44 K であり、余裕でゼロを上回っている。モデルを一度しか実行していない予報士であれば、「霜は降りない」と発表して帰宅するだろう。
-
今度は平均する代わりに数えてみよう。50 のメンバーのうち、20 がその場所、その日で氷点下を示している。尋ねられた質問に対するアンサンブルの答えは 40% である。
-
これらふたつの見解は異なる事象について述べており、平均をどれほど慎重に扱っても二番目の答えは得られない。平均値は、その基礎となる数値のうちいくつが基準線を下回っているかについて何も語らない。すべてのメンバーが +1.44 K を示す五十のメンバーと、−4.95 K から +11.25 K まで分布するここでの五十のメンバーは、どちらも平均すれば +1.44 になる。最初のアンサンブルは霜が降りる可能性を完全に否定するが、二番目のアンサンブルは十回に四回の可能性を与える。
-
8 日目の場所 19 で実際に起きたのは −3.86 K だった。平均は答えの符号を間違えていた。しかしアンサンブルは間違っていなかった。アンサンブルはその出来事に十回に四回の可能性があると述べており、十回に四回の出来事が実際に起きることは、予測の失敗ではないからだ。
-
先ほどの文章が妥当だと言えるのは、40%という数字が文字通り40%を意味する場合のみです。そして、たった一つの事例でそれを証明することはできません。実測データを示したのが四番目のパネルです。このリードタイムにおける1,600地点・日を、アンサンブルが提示した確率ごとに分類しました。30%〜40%の範囲にある130件の予測を見ると、実際に起きた割合は29%です。同様に60%〜70%の範囲では66%となっています。アンサンブルの確率をそのまま信頼してよいのは、それがモデルから出力されたものだからではありません。このようにデータと突き合わせて、しっかりと検証されているからです。
-
事例の選択については正直に言っておこう。これは平均の失敗を示すために選ばれた日であり、ほとんどの日において、天候が閾値のどちら側に転ぶかについて平均とアンサンブルは一致する。重要なのは、平均が通常間違っているということではない。平均はそのような形の質問には全く答えられないが、各メンバーなら答えられるということだ。
-
このページの限界について述べておく。現実の霜注意報では、誤った「穏やかな日」のコストと誤った「警告」のコストを天秤にかける。そしてその計算にはふたつの間違いの価値が必要だが、大気はそれを教えてくれない。「40% は本当に 40% か」という疑問については、我々の calibration ページで適切に扱っている。
解答
「霜は降りない」ではなく、40% と発表しよう。平均は +1.44 K であり、50 のメンバー中 20 が氷点下を示し、実際の結果は −3.86 K だった。 平均は、誰も尋ねていない質問に答えていたのである。
-
学習の道すじ
天気予報の仕組み
参考文献 (5)
- The seminar the model comes from, where Lorenz wrote it down as a toy atmosphere rather than as a chaotic curiosity: E. N. Lorenz, "Predictability: A problem partly solved." Proc. ECMWF Seminar on Predictability, Vol. 1, Reading, ECMWF, 1996, pp. 1-18.
- The paper that put the model into general use, and the source of the K = 40, F = 8 settings this page runs: E. N. Lorenz and K. A. Emanuel, "Optimal Sites for Supplementary Weather Observations: Simulation with a Small Model." Journal of the Atmospheric Sciences 55(3), 1998, pp. 399-414. DOI 10.1175/1520-0469(1998)055<0399:OSFSWO>2.0.CO;2
- What an operational ensemble actually is, including why real ones are under-dispersive and why their perturbations are chosen rather than drawn at random: M. Leutbecher and T. N. Palmer, "Ensemble forecasting." Journal of Computational Physics 227(7), 2008, pp. 3515-3539. DOI 10.1016/j.jcp.2007.02.014
- Where the fourth panel comes from: the decomposition that separates "were the probabilities honest" from "did they say anything useful", which is why a reliability curve is read as a curve rather than scored as a number: A. H. Murphy, "A New Vector Partition of the Probability Score." Journal of Applied Meteorology 12(4), 1973, pp. 595-600. DOI 10.1175/1520-0450(1973)012<0595:ANVPOT>2.0.CO;2
- And the 1963 paper the whole subject starts from, which is a different model and a different lesson: E. N. Lorenz, "Deterministic Nonperiodic Flow." Journal of the Atmospheric Sciences 20(2), 1963, pp. 130-141. DOI 10.1175/1520-0469(1963)020<0130:DNF>2.0.CO;2