最新の気象予測～その理論と技術～
北海道大学地球惑星科学集中講義
札幌管区気象台 室井ちあし集中講義の内容• １日目– 気象予測の基礎– 気象とコンピュータ– 数値予報モデル• ２日目– 数値予報モデル（続）– データ同化– アンサンブル予報– （セミナー）北海道の気象災害リスクと防災気象情報

データ同化北海道大学地球惑星科学集中講義
札幌管区気象台 室井ちあし気象データ処理の流れ
衛

観星

予報値の検 証航空機

ラジオゾンデ地上観測
入電には遅延がある

測
プロファイラ気象レーダー

海上観測

解析・予報サイクル

観測データ処理デコード

品質管理

GTS lines
Internet

第１推定値

客観解析

数値予報
初期値化

数値予報
GSM,MSM,EPS

情報還元(Feedback)

応用処理
ガイダンス,
各種統計

予報：できる限り早く出したい
解析：できる限り観測を待ちたい速報解析とサイクル（遅延）解析

可視化,
ファクシミリ

外部への情報提供
天気予報

気象情報

データ同化とは5

解析（データ同化）とは• 初期値（現在の天気の状態）が分からなければ予報はできない• 初期値を精度よく再現することが予報精度の向上に（決定的に）重要• 初期値を作成するには圧倒的に観測データが少なく、地域的に偏在• そこで、数値予報を観測データで修正して初期値を作成• 正式には「客観解析」あるいは「データ同化」と言う。 （最近は「データ同化」が一般的）

6データ同化の概要• 観測値＋予報値（第一推定値）→

解析値– 全球解析及びメソ解析では４次元変分法を採用– 局地解析及び毎時大気解析では3次元変分法を採用

観測第一推定値

解析値

データ同化

数値予報の格子とは無関係にまばらに存在

通常は、前の初期時刻の予報値観測を反映させた均一な解析値

データ同化• 観測値– ×：全格子・全要素のデータがそろっていない– △：観測誤差を内包している• ○：が、基本的には現実を反映している• 第一推定値– ○：全格子・全要素データがそろっている– △：予報誤差を内包している• ×：現実に沿っているか定かでない• 双方のいい部分を合わせる– 第一推定値を観測値で修正（同化）する• データ同化

同化される観測データの種類と要素• 直接観測– 地上観測、船、ブイ： 地上気圧
– ゾンデ： 風、気温、相対湿度（←湿数から変換)– 航空機： 風、気温（メソ）• 地上リモートセンシング

全球：全球解析のみ
メソ：メソ解析のみ– ウィンドプロファイラ： 風– ドップラーレーダー： 風（メソ）– 気象レーダー（解析雨量）： 降水量（メソ）• 衛星– 大気移動風AMV： 風– サウンダ（鉛直探査計）： 放射輝度温度（全球）、気温（メソ）– マイクロ波放射計： 降水強度（メソ）、可降水量（メソ）– マイクロ波散乱計： 海上風• 人工的なデータ– 観測型台風ボーガス： 風 、海面気圧– PAOB(オーストラリア気象局作成のボーガス）： 海面気圧9

客観解析で使えるデータ#1観測値
ブイ観測

地上観測

ウィンドプロファイラ

ドップラーレーダー

高層観測

台風ボーガス

海上観測航空機観測

直接観測

GPS

リモートセンシング
遠隔観測

DMSP

MetOP

NOAA
MTSAT

（準）現業衛星

METEOSAT(MSG)

COSMICGRACE-A
Aqua

TRMM

地球観測衛星GOES

静止軌道衛星

疑似観測

Terra

GRACE

GPSRO衛星

低軌道衛星

観測データの分布例（全球）（１）SYNOP、船、ブイ

ゾンデ、ウィンドプロファイラ

航空機、豪ボーガス

大気移動風AMV

11

観測データの分布例（全球）（２）サウンダ（気温）マイクロ波散乱計

サウンダ（水蒸気）

マイクロ波イメージャ

12
他に静止衛星の晴天輝度温度データがあります。

データ同化とは・・・
最も確からしい推定値を見つけること

情報＝データ ＋信頼度（誤差情報）
観測値

６．０℃？
どこか？

第一推定値

どちらに寄せるか ８．０℃
４．０℃

13

修正量の決定（誤差情報の利用）• 予報値（背景値）xbには誤差Dxがある• 観測値yoにも誤差Dyがある– 統計情報から誤差の期待値は得られる– それらを基に最適な値を探し、得られた値が解析値xa

解析値
予報値
観測値
D値(Departure)
観測誤差
期待値
修正量
（Increment/Innovation）

値
予報誤差
期待値

昔の手法
「最適内挿法」、「修正法」
解析値

第一＋推定値

適切な×（係数観測値－第一推定値 ）※「第一推定値」は解析時刻を予報した前回の予報値解析値が１回の計算で求められる。
ただし、解析値と要素の異なる観測値は使えない。
適切な係数：観測点と格子点の距離のみの関数（修正法）
観測誤差と背景誤差の関係を統計的に考慮（最適内挿法）
15最近の方法• 変分法– 「解析値」を一発で求めることはあきらめ、「解析値の候補」を少しずつ変えて最適なものを求めることにする。
– 「解析値の候補」は「第一推定値」から出発し、繰り返し計算で「評価関数」を小さくする。
– 評価関数が最小となった時点（収束したとみなせる時点）で止める。
評価関数＝（解析値の候補第一 ２－ 推定値 ）＋（解析値の候補から計算した値－

観測値

解析値が１回の計算では求められない。
しかし、解析値と要素の異なる観測値も使える。
16２） × 係数

3次元変分法と4次元変分法• ３次元変分法：時間の変化を考慮しない• 観測値は解析時刻に得られたものと仮定して比較する• ４次元変分法：時間の変化を考慮する• 同化期間、モデルMで時間発展させた値と比較する– これだけの場合：FGAT(First Guess at Appropriate Time)という手法• 差の還元はアジョイントモデル（接線モデルの転置）で行う

４次元変分法の概念図
観測変数～時間方向への拡張～
観測値

第一推定値
数値モデルの時間積分値

予測

観測値
時間

データ同化期間
緑線は第一推定値からの数値モデルの時間積分値、赤線はデータ同化期間内で第一推定値と観測値に最も近い数値モデルの時間積分値で、後者をその期間内の解析値とする。
第一推定値と観測値につけた縦線（エラーバー）は、それぞれの誤差の大きさを示す。
18

3次元変分法と4次元変分法
気温一点同化の例（風はバランス成分の影響）
3D-Var：背景誤差に従い、基本的には等方

4D-Var：場に従って変形し、等温線に沿った方向に伸びている

色：気温のインクリメント／矢印：風のインクリメント

データ同化サイクル（予報解析サイクル）

データ同化サイクル21解析・予報サイクル

解析予報サイクルにより、観測データが無い領域においても合理的に算出することができる。

解析・サイクルによる全域での解析精度向上の説明図
観測情報の同化

観測情報の同化

予報モデルによる同化情報の移流と伝播

同化情報の充填による推定精度の向上
予報と同化の観測データ有機的な結合
次の時刻の観測データ

時間

その次の時刻の観測データ

全域の精度向上
観測データの希薄な地域でも予報・解析の精度が向上

データ同化サイクル• 第一推定値（前時刻の数値予報）を観測データで修正して初期値を作り、• そこからの予報を次の第一推定値とすることをデータ同化（解析・予報）サイクルという。
• この修正量のことをインクリメントという。
• サイクルにより、観測データの情報を伝播させる。
24

速報解析とサイクル解析• 全球データ同化• 速報解析
– 全球予報のためのデータ同化は初期値の時刻の2時間20分後に処理を開始する－それ以降に入電したデータは使えない。→衛星データの入電が十分でない。

• サイクル解析– 十分なデータ入電を待ってサイクルをつなぐための解析で、データ同化サイクルの解析精度を維持するために行う。
25

全球速報解析で処理されるデータの分布• 使用されたデータを種別・衛星ごとにプロット– 黒はQCなどでリジェクトされたもの– 間引きで除かれた衛星データはここにはプロットされていない

SYNOP、船、ブイ

ゾンデ、ウィンドプロファイラ

航空機、豪ボーガス

大気移動風AMV

サウンダ（気温）マイクロ波散乱計

サウンダ（水蒸気）

マイクロ波イメージャ• データカットオフ時間（入電待ち時間）は
2h20m(00/06/12/18UTC)

静止衛星 晴天輝度温度26

全球サイクル解析で処理されるデータの分布• 使用されたデータを種別・衛星ごとにプロット– 黒はQCなどでリジェクトされたもの– 間引きで除かれた衛星データはここにはプロットされていない

SYNOP、船、ブイ

ゾンデ、ウィンドプロファイラ

航空機、豪ボーガス

大気移動風AMV

サウンダ（気温）マイクロ波散乱計

サウンダ（水蒸気）

マイクロ波イメージャ• データカットオフ時間（入電待ち時間）は
11h35m(00/12UTC),
5h35m(06/18UTC)

静止衛星 晴天輝度温度27

現業数値予報における解析システム
全球解析
目的

メソ解析

局地解析

気候解析

各モデルの初期値作成

領域

解析時刻

6時間毎

3時間毎

毎時

6時間毎

手法

4次元変分法

4次元変分法

3次元変分法

3次元変分法解像度（アウター）

約20km/100層

5km/50層

2km/60層

約110km/40層

解像度（インナー）

約55km/100層

15km/40層

5km/50層

約110km/40層

第一推定値

全球3-9時間予報値

メソ0-3時間予報値

メソ予報値（[012]時間予報）+1時間予報x3

全球6時間
予報値

同化ウインドウ

6時間

3時間

1時間x3

6時間

データ打ち切り時間

速報：2h20m
ｻｲｸﾙ:00,12UTC:11h50m
06,18UTC:07h50m

50min

30min

00UTC:36h35m
06UTC:31h00m
12UTC:25h20m
18UTC:19h35m

品質管理

29

同化システムで考慮する観測誤差
原 因

特 徴

対処方法

偶然誤差

現象自体の揺らぎや測定のばらつきなど

正規分布，不可避

データ同化で考慮

系統誤差

測器の調整不良や物理的変換の誤差など

一定方向のず バイアス
補正
れ

ラフエラー

人為的ミス，測器の故障，通信エラーなど

非常に大きな
誤差

気象学的
ノイズ

システムの表現限界と観測 システムにより データ同化での空間代表性のギャップ変わる
考慮

異常値

30

品質管理で修正･棄却

それではこれらにどのように対処するか• バイアス補正• 品質管理（ＱＣ）による棄却– 気候学的にありえない（気候値チェック）– ひとつの通報の中の情報が互いに矛盾（内的整合チェック)– 第一推定値から離れすぎている（外的整合チェック）31

誤った観測データを同化した例(a)正常な観測を同化した

正しい観測（ｂ）誤った観測を同化したと仮定した場合

誤った観測を同化したと仮定(a)と（ｂ）の差★の位置の船舶が、海面気圧として約980hPaの観測値を誤って20hPaほど低く通報し、それを利用したと仮定した場合の例

航路チェック

修正

通報緯度の20度ずれ
33

通報緯度の30度ずれ修正

要素間整合性チェック（例）雨雪判別と相対湿度・気温の関係
通報値の値がお互いに矛盾がないかどうかを判定する

100

霙

相対
湿 70度

雨雪

(%) 50024
地上気温（℃）346

鉛直整合性チェック• ゾンデの例：赤が通報値、緑が修正後の値– 4か所すべてが気温の正負の誤り
STATION 91643( 8.31S, 179.13E ) Funafuti
1996.03.01.23UTC
通報値
300

修正値
気圧(hPa)

400
500
600
700

8509251000-50 -40 -30 -20 -100

気温(℃)

10203040

ゾンデバイアス補正• ゾンデ気温観測のバイアスの原因– 日射がセンサーに熱を与える– センサーから熱が逃げる– 測器システム特有の誤差• 放射が原因の場合、昼夜に明らかな差• 放射の影響を補正しない観測がある• そもそも原因不明なバイアスもある• →前月のD値統計から気温のバイアスを補正する– 第一推定値にバイアスがないことが前提• 実際にはモデルバイアスはある36

台風ボーガス• 予報課で作られる台風位置情報（中心位置・中心気圧・強風半径）をもとに作成される人工的な台風の構造• 擬似的な観測データとして与える。優先度は最高とする。
• 台風の位置、形状が初期値に反映される。
– 防災上、重要– 強度はモデルの分解能に配慮して調整
37

実況の台風中心第一推定値

ボーガス埋込後

データの間引き・選択• 最適内挿法では１解析点ごとに使用するデータ数の上限があることにより、間引きが必要。
• 変分法では観測が独立である限り、間引きは必須ではない。しかし、代表性誤差や観測演算子の誤差が相関を持つ場合があることを考慮すると、同種の観測は解析格子の密度程度に間引くのが適当。
• 実際問題として、衛星データはデータ量が多すぎるため、間引かないと使えない。間引きが解析結果を大きく左右する場合がある。
38

観測データの特殊な空間代表性が問題になる場合
気温の解析インクリメント

島の局地的な高温が周囲の海上に広がっている
39

島の気温データは使わないノンリアルタイム品質監視• 国際協力– 浮遊ブイ品質監視業務– 第Ⅱ地区地上観測データモニタリング業務（リードセンター）– 全球データモニタリングレポート（主要数値予報センター間で交換）• ブラックリスト管理
低品質のデータがデータ同化に入るのを防ぐ

40

ノンリアルタイム品質管理
観 測

ルーチン解析システム

リードセンターWMOを通じて働きかける

品質情報の蓄積

CDA
統計処理
品質統計情報ブラックリスト
QCパラメータ適正化
バイアス補正平均誤差
RMSE
ﾘｼﾞｪｸﾄ率
誤差特性情報

解析パラメータの設定

予報誤差
相関
41

問題観測点のリスト

浮遊ブイのモニタリング例

突然バイアスが大きくなった例

常に一定の観測値を報じる例（故障）

常に一定のバイアスがある例

台風中心が通過した例（異常ではない）

観測の改善例

モンゴルの観測点

CBSLCから問題点を指摘し、修正された。
43

最近の衛星データ事情

44

衛星データの重要性• 広域かつ均質な情報を、周期的・定常的に得ることができる– 現業的なシステムにとって重要– 海洋上や砂漠上、大気上層の観測# データ空白域は解析精度が悪い ⇒ 大気3(4)次元的な観測情報が必要• 直接観測では入手が困難な、様々な物理量を得ることができる– 例：全球降水量、全球オゾン分布、大気上端外向き長波放射– モデルの検証にも有効• 近年、高精度な衛星データが急増している。
– 測器・通信技術の発達• 同化手法の高度化（変分法の導入など）– 輝度温度の直接同化などにより、観測情報がより適切に活用– 多様な衛星データに対応可能• 1990年代後半から、数値予報精度は、モデル改良だけでなく、同化システム・衛星データ利用によって大きく向上• ただし雲や雨の観測情報の同化はまだまだ難しい45

衛星・センサの分類• 衛星– 静止衛星と軌道衛星• 静止衛星：MTSAT, GOES, METEOSAT• 軌道衛星：NOAA, MetOp, DMSP, TRMM, Aqua, Terra• センサ– 能動的と受動的• 能動型は自ら電波を発射• 能動型はマイクロ波散乱計やレーダー、ライダー。他は受動型– 観測波長（周波数）：紫外、可視、赤外、サブミリ波、マイクロ波• データ同化で主に用いられるのは可視、赤外、マイクロ波

46

衛星データ特有の事情• 開発が複数の機関にまたがる– 連携が重要• 寿命が短い– 打ち上げ後直ちに利用開始できるよう、準備は周到にしておく必要がある• データ量が多く、中身が特殊– データの入手方法の開拓、「放射」に詳しい人材確保衛星データ（輝度温度）のインパクト
Z500 アノマリ相関• 衛星データの有無の影響

衛星データあり• アノマリ相関– 衛星データあり– 衛星データなし

Better– Z500（500hPa GPH）予報

衛星データなし

北半球– 衛星データの有無で、特に南半球でスコアの差が大Better

衛星データあり

衛星データなし
南半球全球解析で利用している観測データ量とスコアの経年変化

49イメージャとサウンダ• イメージャ（放射計）：–（基本的に）一枚絵の画像データ• （可視赤外）MTSAT/Imager, GOES/Imager, NOAA/AVHRR, EOS/MODIS, …
• （マイクロ波）Aqua/AMSR-E（左下図）, TRMM/TMI, DMSP/SSMIS(I), …

• サウンダ（鉛直探査計）：–

複数の画像の差分等からプロファイルを求める• （赤外）NOAA/HIRS, Aqua/AIRS, MetOp/IASI,…
• （マイクロ波）NOAA/AMSU（右下図）, NOAA/MHS, DMSP/SSMIS(S), …

下層

AMSU-Aの荷重関数
A
14

上層A
3

大気追跡風• AMV: Atmospheric Motion Vector– 連続する複数の画像から、雲や大気のパターン追跡を行い、風を見積もるプロダクト• 静止衛星の主要プロダクト(60S-60N)• 近年（2003年頃から）極軌道衛星で極域の軌道の重なりから同様のプロダクトが算出されるようになった極域衛星風マイクロ波イメージャデータ•

同化利用されているプロダクト– 輝度温度（全球、海上・晴天域のみ）– 降水量リトリーバル（メソ、海上・晴天域のみ）– 可降水量リトリーバル（メソ、海上・晴天域のみ）– 下層水蒸気にインパクト

降水分布
マイクロ波放射計
推定降水強度とRAの合成：降水データなし

可降水量分布[mm]

AMSU-A、AMSU-Bの荷重関数
マイクロ波域の吸収線とチャンネル

水蒸気酸素
http://www.eumetsat.eu/groups/ops/documents/image/img_jpg_amsua_channels_big.jpg

熱帯AMSU-A ch3～14 の荷重関数(Janssen 1993)

AMSU-B の荷重関数
左：熱帯 右：高緯度

JMA全球解析で使用中のch高緯度

GNSS掩蔽観測• GNSS-RO(GPS-Radio Occultation)– MetOp/GRAS/COSMIC– GPS衛星（高度2万km）が発信する電波を、低軌道衛星（LEO,高度数百km）のレシーバで受信– 電波は伝搬経路の大気の状態に応じて屈折• 測定されるのはそれによって生じる遅延量• 遅延量屈折角屈折率気温・水蒸気量– 個々のリトリーブ処理（）ではある仮定条件が入る– 水平分解能は低い（～200km）が、鉛直分解能は高い（0.5-1km）

αLEO
GPS

ひまわり８号• 気象庁の静止気象衛星• 2014年10月7日打上げ• 2015年7月7日正式運用開始• 画像の解像度が４倍、頻度は１２倍に• 世界初、カラーでの撮影が可能に• 北海道・江別に運用副局– 正局は埼玉県は富山町・東京都

55