MADOCA プロダクトを利用した可降水量リアルタイム解析実験 Real Time PWV Analysis Experiment using MADOCA Product 小司禎教*1，佐藤一敏*2 Yoshinori Shoji*1, Kazutoshi Sato*2 *1 気象研究所, 2* 宇宙航空研究開発機構 *1 Meteorological Research Institute, *2 Japan Aerospace Exploration Agency (JAXA) １．まえがき JAXA の複数 GNSS 対応高精度軌道時刻推定ツール MADOCA (Multi-GNSS Advanced Demonstration Orbit and Clock Analysis) による軌道情報のリアルタイムインターネット配信を利用し，国土地理院に協力いただき，30 点の GEONET 点の高頻度可降水量(PWV)リアルタイム解析（1 秒間隔で観測後 1 分以内）システムを気象研究所に構築した．解析された可降水量を，近傍の高層ゾンデ観測と比較したところ，RMS で 3mm 程度と高い一致を示した．一方，幾つかの観測点では鉛直座標が時間と共に数十センチ程度ずれていく現象が見られた．鉛直座標のドリフトが大きい場合，天頂遅延量にもバイアスが生じており，修正する必要がある．サイクルスリップの可能性を考え，閾値を変更して実験を行ったが，効果は得られなかった．図 1 2015 年 1 月の比較結果．赤：00UTC，青：12UTC．２．解析概要と結果 GNSS 解析にはオープンソースの RTKLIB ver. 2.4.2 (http://www.rtklib.com/)のツール rtkrcv を利用し，衛星軌道と時計情報は MADOCA リアルタイムプロダクトを利用した．最低仰角は 5 度，等方マッピング関数は GMF を用いた．国土地理院サーバーからストリーミング配信される 1 秒間隔データを精密単独測位（ppp-static）によりリアルタイム解析し，天頂遅延量，大気遅延勾配，及び視線遅延量を逐次解析した．高層ゾンデ観測点近傍の GNSS 観測点で解析された PWV と高層ゾンデによる PWV との散布図を図 1 に示す． GNSS PWV はゾンデの放球時刻から 30 分間平均した．過去の同様の検証調査では，00 と 12UTC の比較で BIAS に違いが見られたので，00 と 12UTC の比較を独立に行った．比較を行った７箇所全てで RMS は 3mm 未満，バイアスも半数以上で 0.5mm 未満であった．図 2 には実験を行った 30 点の 24 日間の鉛直座標解析値の変動を示す．多くの点で，座標が安定した後，変動は 1cm 程度で推移しているが，5 地点（0225，1062，3016， 3018，3019）では 10～数 10cm の変動が発生し，これらの点では天頂遅延量が後半ほど小さな値で解析されている．表 1 には今回の解析地点のアンテナと受信機を示す． TPSCR.G5 と TPS NETG3 の組み合わせの点で鉛直座標のドリフト発生地点が多い．原因の可能性としてサイクルスリップの検出特性が考えられたので，rtkrcv の slipthres パラメータをデフォルトの 50mm から 10mm 単位で減らして解析してみた．図 3 に 3019 の結果を示す．変動の幅に変化は見られるが，大きな効果は得られない．他の地点でも同様の傾向にあり，今後他の原因を探る必要がある．図 2 解析された鉛直座標変動（2015 年 3 月）．解析開始後 1 時間の鉛直座標からの差を示す．表１アンテナ機種，受信機機種による地点の分類図 3 rtkrcv の slipthres パラメータの値を変化させた場合の鉛直座標解析値の変化．3019 点の例．（デフォルトは 50mm）