EKF EKF EKF EKF・Eꜳ Eꜳ Eꜳ Eꜳ の統合 統合 統合 統合 同化 同化 同化 同化実験の開 実験の開 実験の開 実験の開 * 折口 征二 青梨 和正 川畑 拓矢 國井 勝 (気象研究所) 1. 1. 1. 1. はじに はじに はじに はじに 現在、同化の主流になってきてい、 に基づくと変分法の統 合さた同化実験の開てい。こに 1 つの同化でと変 分法の両方の解析すことが可能とな。た、新し い観測の、位置ず補正間引きなど 解析前処理として実装す、どちの同化方法で 扱え利点持つ。EKF は、解析予報として 実装さてい NHM-LEKF (藤田 2011.11;第 1 図) 使用した。同化さてい観測は、従来型観測 のか、GPS 可降水量、動径風、衛星 輝度温度と豊富であ。1 の処理の流として は、「初期場初期摂動の準備→予報→ (観測の取込)→QC→解析→地表面解析→ 保存」とな。この解析分に対応す Eꜳ の 開し NHM-LEKF に組込という方法 でた。従って、実際は NHM-Eꜳ の開ということ にな。 第 1 図 NHM-LEKF 解析予報 2. NHM 2. NHM 2. NHM 2. NHM-Eꜳ Eꜳ Eꜳ Eꜳ の開 の開 の開 の開 基 本 な は 、 E ꜳꜳꝏꜳ Aꜳꝏ(ꜳ 2005; ꜳ ꜳ. 2008; L ꜳ. 2008)に従い、評価関数は以下の式に 表さ。ここで、 は予報平均値(第一 推定値)、 は予報()、 は解析 値、 は背景誤差共分散行列、 は観測誤差共分散 行列、 は数の次元 N 持つ制御変数で あ。摂動で背景誤差共分散行列の 2 乗 根 構成す。 開した NHM-Eꜳ の実際のは以下 青枠で示す。空間の評価関数とその勾 か LBFGS 法用いて行い、 最な解析値計算す。Eꜳ では勾求 時、摂動用いことにな、接線型 随伴必要としない。た、背景誤差共分散 行列は、同化時刻の予報か構成さ た ffꝏ- とな。 空間の評価関数 勾 3. 1 3. 1 3. 1 3. 1 点観測に 点観測に 点観測に 点観測に同化 同化 同化 同化実験 実験 実験 実験と結果 と結果 と結果 と結果 NHM-Eꜳ の動作確認すた、1 点観 測に同化実験行った。2011 年台風第 12 号 (ꜳꜳ)の事例対象として、台風中心の北東側、高度 5 の 場 所 に 、 風 の 擬 似 観 測 (= - 35.54(/), =45.68(/), ꝏ=0.5(/))与えた。風速は第一 推定値の 2 倍とし、観測値取込すくすた観 測誤差は小さく設定した。領域は四国沖で解像度 10、 格子数 180140、同化時刻は 2011 年 9 月 1 日 1200 a 2 / 1 f a - w P x x f ⋅ = [ ] [ ] f f N f f 2 f f 1 f N f 2 f 1 2 / 1 f , , , 1 1 , , , x x x x x x p p p P - - - - = = L L N [ ] T a N a 2 a 1 a w , , w , w L = w a w ] - ) [H( ] - ) [H( 2 1 ] - [ ] - [ 2 1 ) J( a 1 - T a f a 1 - T f a a y x R y x x x P x x x f + = ] - ) [H( ] - ) [H( 2 1 ) ( 2 1 ) J( a 1/2 f f 1 T a 1/2 f f a T a a y w P x R y w P x w w w ⋅ + ⋅ + + = - ] - ) [H( H J a 1/2 f f T a 1 a a y w P x w R w w ⋅ + ∂ ∂ + = ∂ ∂ - ] , , , [ H H H H f N f 2 f 1 a 1/2 f a a a a a p p p x P x w x x w L ⋅ ∂ ∂ = ⋅ ∂ ∂ = ∂ ∂ ⋅ ∂ ∂ = ∂ ∂ ) H( - ) H( H a f i a f i a x p x p x + ≈ ⋅ ∂ ∂ [ ] f f i f i 1 1 x x p - - = N LBFGS LBFGS LBFGS LBFGS a a ~ w w = step f x a x f P f i x R 2 / 1 f P a 2 / 1 f a ~ ~ w P x x f ⋅ + = 0 J ), J( 0 a a a = ∂ ∂ = step w w w 境界値: mfboundary+間摂動 10kmNHM(40) 10kmNHM 初期摂動は間 6hr window (slot=7) 10kmNHMの再積分() 同化してい観測 従来型観測(地上, 高層, 航空機など) 衛星風 可降水量 GPS, TRMM/TMI, SSM/I, (Aqua/AMSR-E) 動径風 輝度温度 解析 解析 解析 解析 解析 解析 解析 0h 6h 12h 18h