(57)【要約】【課題】  雨水量のレーダデータ同化を終了した後に予測計算を継続しても、モデル領域内での降雨状況を模擬し続けられるようにする。【解決手段】  気象予測モデルに対して雨水量のレーダデータ同化を行って雨水量を算出した後（Ｓ２１）、モデル領域内のある定義点の周囲に同化された雨水量が既定値以上存在する（水蒸気が豊富に存在する）かどうかの水蒸気量の推定計算を行い（Ｓ２２）、周囲にレーダ同化された雨水が存在するか否かを判定する（Ｓ２３）。この判定で、雨水が存在すると判定された場合には、水蒸気ボーガスを作成してモデル内にデータ同化する（Ｓ２４）。これにより、モデル領域内でレーダデータ同化された雨水が直ちに落下しても、空中に漂う水蒸気をモデル領域に残すことができる。降雨地点は周囲も含めて擾乱が激しいので、その水蒸気が凝結して雲や雨を生成することが期待され、モデル領域内で「雨を降り続かせる」ことが可能になる。【選択図】      図３
2006038583.tif 000002
【特許請求の範囲】
【請求項１】
気象予測モデルに気象レーダの雨水量データを同化する雨水量同化手段と、
前記雨水量データ同化の結果、周囲に同化された雨水量が既定値以上存在する地点を判定する判定手段と、
前記雨水量が既定値以上存在する地点について水蒸気ボーガスを作成して前記気象予測モデル内に同化する水蒸気ボーガス同化手段とを具備し、
前記水蒸気ボーガス同化手段は、前記雨水量が既定値以上存在する地点についてデータ同化前のモデル内で表現されている水蒸気量と同点の気温と気圧で求められる飽和水蒸気量とを比較しナッシング処理して飽和水蒸気量に近づけることを特徴とする気象予測システム。
【請求項２】
前記ナッシング処理は、  
qvassim = qvmod * (1-alpha) + qvsat * alpha  
qvassim：同化後の水蒸気量（未知量）  
qvmod  ：同化前の（モデルで計算される）水蒸気量（既知量）  
qvsat  ：定義点の飽和水蒸気量（既知量）  
alpha  ：重み付けパラメータ（０＜alpha＜１）  
により同化後の水蒸気量をqvassim を求めることを特徴とする請求項１記載の気象予測システム。
【請求項３】
前記ナッシング処理を実行するために、水蒸気量同化を行う際に周囲のレーダ同化された雨水の検索範囲、検索する雨水量のしきい値、水蒸気量同化の際の重み付けパラメータを予め決定しておくことを特徴とする請求項２記載の気象予測システム。
【請求項４】
気象予測モデルに気象レーダの雨水量データを同化する雨水量同化ステップと、
前記雨水量データ同化の結果、周囲に同化された雨水量が既定値以上存在する地点を判定する判定ステップと、
前記雨水量が既定値以上存在する地点について水蒸気ボーガスを作成して前記気象予測モデル内に同化する水蒸気ボーガス同化ステップとを具備し、
前記水蒸気ボーガス同化ステップは、前記雨水量が既定値以上存在する地点についてデータ同化前のモデル内で表現されている水蒸気量と同点の気温と気圧で求められる飽和水蒸気量とを比較しナッシング処理して飽和水蒸気量に近づけることを特徴とする気象予測方法。
【請求項５】
前記ナッシング処理は、  
qvassim = qvmod * (1-alpha) + qvsat * alpha  
qvassim：同化後の水蒸気量（未知量）  
qvmod  ：同化前の（モデルで計算される）水蒸気量（既知量）  
qvsat  ：定義点の飽和水蒸気量（既知量）  
alpha  ：重み付けパラメータ（０＜alpha＜１）  
により同化後の水蒸気量をqvassim を求めることを特徴とする請求項４記載の気象予測方法。
【請求項６】
前記ナッシング処理を実行するために、水蒸気量同化を行う際に周囲のレーダ同化された雨水の検索範囲、検索する雨水量のしきい値、水蒸気量同化の際の重み付けパラメータを予め決定しておくことを特徴とする請求項５記載の気象予測方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、気象予測モデルに各種センサ観測データ、気象レーダデータを入力して気象現象をシミュレーションし予測する気象予測システムに関する。
【背景技術】
【０００２】
現在、水平領域５００ｋｍの気象予測モデルでは、水平領域が数千ｋｍの広域気象予測モデルの予測データを初期値として計算を実行している。また、気象レーダにおいて観測される雨水量、風速ベクトルをデータ同化して将来の予測計算を行うオプションが付加されている（例えば特許文献１参照）。
【０００３】
ところで、上記気象レーダの雨水量データを気象予測モデルにデータ同化して予測計算の実験を多岐に渡って行った結果、レーダデータ同化期間内での雨水量予測はアメダス（ＡＭｅＤＡＳ：Automated Meteorological Data Acquisition System地域気象観測システム）雨量データとよい一致をみせることが多い。しかしながら、レーダデータ同化を終了した後にそのまま予測計算を続けると、モデル内で現れていた雨が直ちにやんでしまい、観測との一致が見られなくなる現象が生じている。
【特許文献１】特開２００３−０９０８８８号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
以上述べたように、従来の気象レーダの雨水量データを利用した気象予測モデルによる気象予測システムでは、レーダデータ同化を終了した後にそのまま予測計算を続けると、モデル内で現れていた雨が直ちにやんでしまい、観測との一致が見られなくなる現象が生じている。
【０００５】
本発明は上記の問題を解決し、雨水量のレーダデータ同化を終了した後に予測計算を継続しても、モデル内での降雨状況を模擬し続け、これによって雨水量に関して高精度な予測が可能な気象予測システム及び気象予測方法を提供することを目的とする。
【課題を解決するための手段】
【０００６】
上記問題を解決するために、本発明に係る気象予測システム及び気象予測方法は、気象予測モデルに気象レーダの雨水量データを同化する雨水量同化手段（ステップ）と、前記雨水量データ同化の結果、周囲に同化された雨水量が既定値以上存在する地点を判定する判定手段（ステップ）と、前記雨水量が既定値以上存在する地点について水蒸気ボーガスを作成して前記気象予測モデル内に同化する水蒸気ボーガス同化手段（ステップ）とを具備し、前記水蒸気ボーガス同化手段（ステップ）は、前記雨水量が既定値以上存在する地点についてデータ同化前のモデル内で表現されている水蒸気量と同点の気温と気圧で求められる飽和水蒸気量とを比較しナッシング処理して飽和水蒸気量に近づけることを特徴とする。
【発明の効果】
【０００７】
本発明によれば、雨水量のレーダデータ同化を終了した後に予測計算を継続しても、モデル内での降雨状況を模擬し続け、これによって雨水量に関して高精度な予測が可能な気象予測システム及び気象予測方法を提供することができる。
【発明を実施するための最良の形態】
【０００８】
以下、図面を参照して本発明の実施の形態を詳細に説明する。
【０００９】
図１は本発明が適用される気象予測システムの概略構成を示すブロック図である。図１において、通信処理部１１は、通信インターフェース１２を通じてネットワークＮＴに接続され、当該ネットワーク上の例えば気象庁等の気象データサーバＤＳ０から配信される全国規模の気象観測データ（以下ＧＰＶ（Grid Point Value）データと記す）と、予測対象地域及びその周辺のレーダサイト等のデータサーバＤＳ１，ＤＳ２，…から提供される局所的な気象観測データを入手する。この通信処理部１１で入手された気象観測データは観測データ格納部１３に格納され、気象予測モデル演算部１４からの要求に応じて選択的に演算部１４に送られる。また、この演算部１４で求められた気象予測データは予測データ格納部１５に蓄積される。
【００１０】
上記気象予測モデル演算部１４は、まず観測データ格納部１３からＧＰＶデータを取り込んで初期値とし、気象予測モデルを演算した後、レーダサイト等の局所的な気象観測値を観測値格納部１３から取り込んで気象予測モデルの空間格子点に内挿し、データ同化処理を行って予測値の補正を行う。補正後の予測値（予測データ）は予測データ格納部１５に格納される。次に、予測データ格納部１５から前回の計算された予測データを取り込み、これを初期値として気象予測モデルを演算する。以後、新たなＧＰＶデータが得られるまで、補正後の予測値を初期値として用い、気象予測モデルを更新する。
【００１１】
上記気象予測モデル演算部１４のデータ同化処理の流れを図２に示す。図２において、まず観測値と予測値を取り込み、品質管理処理を行う（Ｓ１１）。ここでは、観測値と予測値を比較することにより、品質の悪いデータを同化から除外する。
【００１２】
次に、観測値と予測値の統計的な誤差特性をもとに別途計算された内挿重み算出値を用いて、品質管理後の観測値を空間的に気象予測モデルの格子点に内挿する（Ｓ１２）。このとき、最初のモデル計算では、初期値にＧＰＶデータを用い、次のモデル計算からは補正後の予測データを用いる。
【００１３】
内挿重み算出処理後、データ同化処理を行う（Ｓ１３）。このデータ同化処理は、内挿した結果を数値モデルに時間的に連続して取り込む。予測値とこの内挿した観測値を比較し、予測値を補正する。
【００１４】
上記システムにおいて、気象レーダの雨水量データをデータ同期して予測計算を行った場合、データ同化期間内での雨水量予測は実際の観測データとほぼ一致した値が得られるが、データ同期を終了した後にそのまま予測計算を続けると、モデル内で現れていた雨が直ちに止んでしまい、観測との一致がみられなくなってしまう。
【００１５】
この要因としていくつか考察されるが、「雨水量のみをデータ同化しており、雨を生成する水蒸気をデータ同化していない」ことがデータ同化終了後に直ちに雨が止んでしまう主要因であると考えられる。すなわち、雨水量のみをデータ同化したとしても、雨水は有意な落下速度を持ち、直ちに地表へ落ちてしまうことにより、落下後はモデル内に殆ど影響を及ぼさないと考えられ、その後もモデル内で「雨を降り続かせる」ためには、雨水を作る物質、つまり水蒸気が欠かせないと考えられる。
【００１６】
そこで、水蒸気量同化を目的として様々な試みを行ったが、水蒸気の生の観測データを手に入れてモデルにデータ同化することは問題が多い。特に、ＧＰＳ（Global Positioning System：衛星測位システム）衛星等による水蒸気データの同化手法は、データ蓄積を行っているが、実運用のためのルーチンデータを入手できないため、現状では非常に難しい。そこで、本発明は「気象レーダデータを用いた水蒸気ボーカス」を提案する。この「水蒸気ボーカス」とは、レーダデータから水蒸気量を推定し、モデル内にデータ同化することを意味する。以下、図３に本発明に係る水蒸気ボーカスを用いた水蒸気量データ同化処理のフローチャートを示してその処理内容を説明する。
【００１７】
まず、気象レーダにおける雨水量のデータ同化を行うことで、同化される雨水量を算出する（Ｓ２１）。雨水量データの同化を行った後、水蒸気量の推定計算を行う（Ｓ２２）。この「水蒸気量の推定計算」は、モデル計算領域内のある定義点の周囲に同化された雨水量が既定値以上存在すれば、その点では周囲の雨水を生成した水蒸気が飽和水蒸気量を上限として豊富に存在するという考えを基本とする。
【００１８】
水蒸気量の推定計算の結果、周囲にレーダ同化された雨水が存在するか否かを判定する（Ｓ２３）。この判定で、存在しない場合にはそのまま雨水量データ同化処理を終了する。存在すると判定された場合には、水蒸気ボーガスを作成してモデル内にデータ同化する（Ｓ２４）。
【００１９】
具体的には、その定義点においてデータ同化前のモデル内で表現されている水蒸気量と同点の気温と気圧で求められる飽和水蒸気量とを比較し、以下の式を用いてナッシングすることで飽和水蒸気量に近づける、つまり雲水や雨水を生成しやすくする推定方法をとる。
【００２０】
qvassim = qvmod * (1-alpha) + qvsat * alpha  
qvassim：同化後の水蒸気量（未知量）  
qvmod  ：同化前の（モデルで計算される）水蒸気量（既知量）  
qvsat  ：定義点の飽和水蒸気量（既知量）  
alpha  ：重み付けパラメータ（０＜alpha＜１）
但し、上記の方法では、以下の３点をパラメータ実験によって予め決定する必要がある。  
・水蒸気量同化を行う際に周囲のレーダ同化された雨水の検索範囲  
・検索する雨水量のしきい値  
・水蒸気量同化の際の「重み付けパラメータ」
上記の処理を実行することにより、レーダデータ同化された雨水は直ちに落下してモデル領域内に影響を残さないが、空中に漂う水蒸気をモデル領域に残すことができる。雨が降っている地点は周囲も含めて擾乱が激しく、水蒸気量を同化させることによりよ、その水蒸気が凝結し、雲や雨を生成することが期待される。これにより、モデル内で「雨を降り続かせる」ことが可能になると考えられる。
【００２１】
したがって、上記構成による気象予測システムでは、雨水量のレーダデータ同化を終了した後に予測計算を継続しても、モデル内での降雨状況を模擬し続けることができ、これによって雨水量に関して高精度な予測が可能となる。
【００２２】
尚、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
【図面の簡単な説明】
【００２３】
【図１】本発明が適用される気象予測システムの概略構成を示すブロック図。
【図２】図１のシステムの気象予測モデル演算部のデータ同化処理の流れを説明するためのフローチャート。
【図３】図１のシステムに適用される、本発明の特徴となる水蒸気量のデータ同化処理過程を説明するためのフローチャート。
【符号の説明】
【００２４】
１１…通信処理部、
１２…通信インターフェース、
１３…観測データ格納部、
１４…気象予測モデル演算部、
１５…予測データ格納部、
ＮＴ…ネットワーク、
ＤＳ０…気象データサーバ、
ＤＳ１，ＤＳ２…レーダサイトデータサーバ。
【図１】
2006038583.tif 000003
【図２】
2006038583.tif 000004
【図３】
2006038583.tif 000005