シミュレーションとリアルを融合するデータ同化の技術で、これから何が起こるかをより正確に予測する
理研数理は、理化学研究所 計算科学研究センター（R-CCS）データ同化研究チームと連携し、データ同化技術の活用を支援します。
データ同化とは
シミュレーションは現実世界をモデル化して行われるので、その結果と現実世界の間にはどうしてもずれが出てきます。そこで、シミュレーションを実際の観測データとつきあわせ、シミュレーションの軌道を修正して「確からしさ」を高めることが行われています。これが「データ同化」です。理化学研究所R-CCSデータ同化研究チームでは、これまで気象分野で開発してきたデータ同化の手法を高性能スーパーコンピュータで効果的に使えるように高度化するとともに、さらに先端的なデータ同化理論の開発に取り組んでいます。これにより、気象シミュレーションの大規模化・高精度化への貢献が期待される一方、ほかの分野にもデータ同化の適用が広まっていくことが期待されます。
データ同化の手法
データ同化は統計数理学に基づいて行われ、いくつかの手法が知られています。理化学研究所R-CCSデータ同化研究チームでは、 LETKF (Local Ensemble Transform Kalman Filter: 局所アンサンブル変換カルマンフィルタ）法をはじめとする手法を用いて、より精度の高い計算を行うための改善を重ねてきました。
LETKF法
少しばらつきをもった複数の初期値（オレンジの楕円）からシミュレーションを行うと、結果のばらつきは大きくなります（青の楕円）。そこで観測値（赤の円）とつきあわせて結果を絞り込み（緑の楕円）、それをこの時点からのシミュレーションの初期値とします。
データ同化の応用事例
ゲリラ豪雨予測
ゲリラ豪雨を引き起こす積乱雲は、わずか数分の間に急激に発生、発達します。理化学研究所R-CCSデータ同化研究チームでは、3次元的に隙間なくスキャンするフェーズドアレイ気象レーダの特長を生かした立体的な「3D降水ナウキャスト手法」を開発しました。30秒毎に得られるフェーズドアレイ気象レーダの観測データをリアルタイムに処理し、10分後までの予測データを算出することができます。
スーパーコンピュータ「京」を用いたシミュレーションによるゲリラ豪雨の再現
理化学研究所R-CCSデータ同化研究チームは、2014年9月11日に神戸市付近に降った「ゲリラ豪雨」をシミュレーションで再現しました。
「京」と「フェーズドアレイ気象レーダ」の双方から得られる高速かつ膨大なデータを組み合わせることで、解像度100mで30秒ごとに新しい観測データを取り込んで更新する、空間的・時間的に桁違いの天気予報シミュレーションを実現し、実際のゲリラ豪雨の動きを詳細に再現することに成功しました。
詳細は、2016年8月9日プレスリリース「「京」と最新鋭気象レーダを生かしたゲリラ豪雨予測」を参照ください。
台風の予測
気象の予測はスーパーコンピュータを使った天気予報シミュレーションに基づいており、シミュレーションと実測データを融合する「データ同化」が予報の精度を左右します。理化学研究所R-CCSデータ同化研究チームは、ひまわり8号の10分ごとの赤外放射輝度データを、雲域も含めた全天候で数値天気予報に直接利用することに成功し、2015年最強の台風第13号（Soudelor）の急発達の予測の大幅な改善や、2015年9月関東・東北豪雨の雨量予測の改善などを実現しました。
COVID-19
COVID-19の流行は私たちの日々の生活を大きく変えました。様々な分野の研究者たちにとっても、COVID-19とそれがもたらす影響を理解し対抗するため、その知見を駆使することが求められています。理化学研究所R-CCSデータ同化研究チームでは、日々の入院者数、死亡者数、回復者数の観測を用いて、潜在的な無症候性者数および症候性者数を再帰的に分析し、同時に無症候性者罹患率に関するパラメータの推定を行っています。機械学習とデータ同化のアイデアを融合させることで、より信頼性の高い予測につながることが期待されています。
生命医科学における予測の問題に対する深層学習の応用
図は集中治療室のデータから敗血症の発症を予測するための再帰型ニューラルネットワークの基本構造の1例を示しています。最終的な予測モデルはこのユニット6個の合議制でできており、極めて高い予測性能を示しています。しかし深層学習を使っている以上、この強力な予測モデルがどういう仕組みで判断を下しているのか人間には理解ができません。そこでデータ同化など数理的手法を深層学習と融合させることで、説明不可能な予測システムを説明可能に変えていく方法を探求されています。
細胞内のシグナル伝達
細胞性粘菌の走化性シグナル伝達系の鍵となる分子であるPIP3は、走化性誘引物質の濃度勾配があると、その勾配に沿って細胞膜上に極性を形成することが知られています。しかし、走化性誘引物質の刺激がなくても、自発的に非一様分布をランダムな方向に形成し、進行波などのダイナミックな振る舞いを示します。このシグナル伝達系の数理モデルをデータ同化の方法と組み合わせることで、細胞の状態を推定する方法の開発が進んでいます。