招待論文
熱帯降雨観測衛星搭載降雨レーダの研究開発とその観測成果
岡本 謙一† a)

重尚一††

TRMM Precipitation Radar and Its Observation Results
Ken’ichi OKAMOTO†a) and Shoichi SHIGE††

あらまし 熱帯降雨観測衛星（TRMM：Tropical Rainfall Measuring Mission）は，1997 年 11 月 28 日に打ち上げられて以来，現在に至るまで，約 10 年以上も順調に，熱帯・亜熱帯降雨の観測を継続している．TRMM
は，世界で初めてのアクティブフェイズドアレー方式の降雨レーダ（PR）を搭載しており，台風などをはじめとする様々な降雨システムの三次元構造の観測にその威力を発揮してきた．論文では，TRMM 衛星のミッショ
ンの目的，搭載センサの概要，TRMM 降雨レーダシステムの概念設計時の諸検討，TRMM 降雨レーダシステ
ムの概要，TRMM 降雨レーダデータ解析処理アルゴリズムシステム，様々な観測成果（台風，潜熱加熱など），
TRMM を継承する全球降水観測計画（GPM：Global Precipitation Measurement）主衛星搭載 2 周波降水レーダ（DPR）並びに更に次世代の衛星搭載の降雨レーダ技術についての展望などについて解説する．
キーワード

熱帯降雨観測衛星，降雨，降雨レーダ，降雨レーダアルゴリズム，全球降水観測計画衛星

1. ま え が き

を通して，降水量を推定してきたが，赤外放射計は，

降水は，地球的規模の水循環，エネルギー循環，大らの放射輝度温度情報のみから降雨強度を推定して気大循環の駆動源であり，気候変動や異常気象現象とおり，降雨そのものを直接観測することが困難である

深いかかわりがあり，人間活動や社会システムに重要という，原理的な問題を有している．これに対して近な影響を及ぼす．信頼性のある地球的規模の降水観測

年，赤外放射計に比べて降雨からの放射をより直接的

データは，水循環モデルの構築，生態系環境の維持，に観測できる衛星搭載のマイクロ波放射計が全世界で農業生産等の社会基盤にとって必要不可欠なもので，利用されるようになってきた．衛星搭載マイクロ波放人類にとっての昔からの重要な課題であった．地上雨射計の実現以来，降水強度の観測は信頼に値する定量

量計による観測は本質的に点観測であり，全球規模の的なものとなってきた．このことは，熱帯降雨観測衛

降水量の時間空間変動を正確に測定することは困難で星（TRMM：Tropical Rainfall Measuring Mission）

ある．また，地上の降雨観測網は，北半球の先進国地

搭載の世界で初めての降雨レーダの実現によって更に雲の下層部の降雨の有無にかかわらず，雲の上層部か域に偏在している．とりわけ，熱帯地方は，大部分が加速された．以下，TRMM 搭載降雨レーダとその観海洋や熱帯雨林に覆われており，地上からの降雨観測

測成果について紹介する．

が困難である．このため，衛星からのリモートセンシ
ングがほぼ唯一の観測手段となる．既に 1970 年代から，静止気象衛星や NOAA 衛星搭載の赤外の雲画像

2. 熱帯降雨観測衛星（TRMM）衛星の概要
全世界の降雨量の 2/3 以上を占める熱帯降雨は，地†

鳥取環境大学環境情報学部，鳥取市

球規模の大気循環を駆動するエンジンの役割を果たし

Faculty of Environmental and Information Studies, Tottoriている重要な物理量であるにもかかわらず，観測が困

University of Environmental Studies, 1–1–1 Wakabadai-kita,††

Tottori-shi, 689–1111 Japan

難な量である．特に，熱帯太平洋域では広大な海洋が大阪府立大学大学院工学研究科，堺市

広がっており，陸上では熱帯雨林が広がっており，観Graduate School of Engineering, Osaka Prefecture Univer-

測データは非常に少ない領域である．したがって，熱sity, 1–1 Gakuen-cho, Naka-ku, Sakai-shi, 599–8531 Japana) E-mail: kokamoto@kankyo-u.ac.jp

電子情報通信学会論文誌 B Vol. J91–B

帯降雨を観測するためには衛星観測しか方法がないと
c （社）電子情報通信学会 2008
No. 7 pp. 723–733 

723

電子情報通信学会論文誌 2008/7 Vol. J91–B No. 7

いっても過言ではない．このために，降雨観測センサ

計である．観測周波数は，10.65， 19.35， 21.3，37，を搭載した TRMM 衛星が計画された [1]．

85.5 GHz の 5 周波数を有し，21.3 GHz 以外は，垂直
偏波と水平偏波の受信ができる 9 チャネルの全電力タ
イプの受信機から構成される．入射角 49◦ の円錐状機
械的走査を行い，走査幅は 760 km である．水平距離

TRMM は，1997 年 11 月 28 日打ち上げられた日
米共同の人工衛星であり，エルニーニョ現象をはじめとする地球的規模の気候変動に影響を及ぼす熱帯
及び亜熱帯降雨を観測することを目的とする [2], [3]．
TRMM の軌道は，高度 350 km（衛星のミッション寿
命を延ばすために 2001 年 8 月から，大気密度の薄い
高度 402.5 km に上昇），軌道傾斜角 35◦ の円軌道で分解能は，6〜50 km である．海上の降雨観測には，主に低い周波数帯による降雨からの放射の観測データが用いられ，陸上の降雨観測には高い周波数帯による降雨上空の雪氷の散乱の観測データが利用される．VIRSあり，経度 0◦ 〜360◦ ，緯度 ±35◦ の範囲を緯度経度

は，観測波長が，0.63，1.61，3.75， 10.8，12 μm の

5◦ × 5◦ の格子に分けて，各格子ごとの月平均降雨強度を 20%の精度で観測することを目的とした．軌道
傾斜角 35◦ をとったのは，熱帯，亜熱帯降雨（特に5 波長であり，走査幅 720 km のクロストラック（衛星進行方向と直角）走査を行い，水平距離分解能は梅雨前線）の観測頻度を増やすためであったが，これ2 km である．可視から赤外までの五つの波長の放射を受信し，雲画像を取得する．PR は，雨滴によるレーによって同時に，様々な地方時刻で降雨観測を行うこダ送信電波の散乱電力を受信して，降雨強度を推定すとができるようになった．通常の地球観測衛星は，極る．PR は，パルスレーダであり，レンジ方向の分解軌道に近い軌道傾斜角を有する太陽同期軌道をとり，能を有する．これによって降雨の三次元的構造を観測

地方時がほぼ同じ時刻において地表を観測する．しかすることが可能であり，これまで観測が困難であった

しながら，降雨の日周変化が顕著な熱帯・亜熱帯（特台風やハリケーンの内部構造を観測することに成功しに陸上）において，特定の地方時のみ降雨を観測すれてきた．また，エルニーニョの時期と通常の時期とでば，観測降雨強度にバイアスを有することになるので，

は，降雨の全世界的な分布が大きく異なることを明確

TRMM 衛星は様々な地方時刻で降雨観測を行う必要に示すことに成功した．また，TRMM 衛星は，このがあった．

ほかに，CERES（Clouds and the Earth’s Radiant（TMI：TRMM Microwave Imager），可視赤外放射

Energy System：雲及び地球放射エネルギー観測装置）と LIS（Lightning Imaging Sensor：雷観測装置）と呼ばれるセンサを搭載している．CERES は，可視か計（VIRS：Visible Infrared Scanner）並びに我が国

ら赤外までの広いスペクトル帯の雲及び地球表面からが世界で初めて開発した衛星搭載用降雨レーダ PRの放射を受信する．LIS は，雷の放電頻度の地球規模TRMM 衛星 の 外 観 を 図 1 に示 す．TRMM 衛星
は ，降雨 観 測 用 の セ ン サ と し て マ イ ク ロ 波放 射 計（Precipitation Radar）を搭載している．TMI は，米での分布を測定する．TRMM 衛星の設計は，3 年で国の軍用気象衛星 DMSP 搭載のマイクロ波放射計

あるが，1997 年 11 月 28 日の打上げ以来，設計寿命

SSM/I の流れをくむ機械的走査型のマイクロ波放射

をはるかに超えた現在に至るまで，約 10 年近く降雨データを収集しつつある．この間，CERES 以外のセンサは，不具合なく動作し，データの収集を続けている．また，衛星の本体であるバス機器も太陽電池パド
ルの駆動系に劣化が見られる以外は順調に動作している．

3. TRMM 降雨レーダシステム
3. 1 TRMM 降雨レーダの概念設計
TRMM 衛星の構想は，1985 年ごろに米国 NASA ゴ
ダード宇宙飛行センターの G. North，T. Wilheit，並びに O. Thiele たちによって提案された．同年，NASA
図 1 TRMM 衛星
Fig. 1 TRMM satellite.

724

ゴダード宇宙飛行センターで開催された TRMM に関する国際科学会合の後，1986 年 6 月に開催された

招待論文／熱帯降雨観測衛星搭載降雨レーダの研究開発とその観測成果
表 1 TRMM 降雨レーダに対するミッション要求条件
Table 1 Mission requirement for the TRMM
precipitation radar.

直下方向を中心にして，衛星の進行方向に直交する平面内を往復走査する方式をとることにした．アンテナ
ビームが走査範囲をすき間なく走査するには，衛星が水平分解能の約 4 km の距離を移動する時間の間（衛星の地表軌跡の移動速度を 7.3 km/秒とするとき，約0.6 秒となる）に ±17◦ の範囲を走査する必要がある．
このような高速でアンテナビームを走査するには，機械的な走査よりも電子的な走査が有利となってくる．

PR の概念設計において，レーダタイプを選択する主
要技術課題として，（1）パルス圧縮方式か，通常の非圧縮パルスレーダか，（2）パッシブアレー方式か，ア
クティブアレー方式か，（3）導波管スロットプラナー
アレーアンテナ方式かオフセットシリンドリカルパラ
ボラアンテナ方式か，の 3 点が比較検討された [6]．結宇宙分野における日米常設幹部連絡会合（SSLG）に論的には，パルス圧縮方式は，地表面クラッタの影響

おいて，日米共同で新規の協力項目として提案，承認を受けないで降雨を観測するためには，レンジサイドされた．同年 12 月より，米国側の強い要請によって，
SSLG のもとで日米両国の科学者を集めて様々な観点

ローブレベルを −60 dB 程度に抑える技術開発の必要から Feasibility Study が約 1 年間実施され，1988 年

径は，アンテナビーム幅，利得，サイドローブ特性な性があることから採用されなかった．アンテナの開口

4 月に報告書が完成した．同 Feasibility Study におい

どによって決まる．アレー素子間隔は，最大走査角にて，当時の郵政省通信総合研究所は降雨レーダを担当

おいてグレーティングローブの生じないための条件よ

し，NASA 側の科学者側によって示された最低限度のり決まる．導波管スロットプラナーアレーアンテナのミッション要求条件に基づいて，様々なタイプの降雨

場合のアンテナ開口径は，約 2.1 m となる．一方，シ

レーダについて検討した [4], [5]．このミッション要求

リンドリカルパラボラアンテナのビーム走査方向の寸条件を表 1 に示す．周波数については，当初は，宇宙法は，ビーム走査を行うため一次ホーンアンテナの配研究用に二次業務として割り当てられていた周波数の列の長さよりもかなり大きく約 3 m となる．このたうちの 2 周波数（13.8 GHz，24.15 GHz）を対象としめに，また，アンテナタイプとしては，搭載性の観点て考えたが，衛星搭載の重量，寸法並びに消費電力のから小さい占有面積をもつ導波管スロットプラナーア

制限から，最終的には，13.8 GHz の 1 周波数のレー

レーアンテナが選ばれた．

ダとなった．降雨減衰の影響を極力減らすためにも，

パッシ ブ ア レ ー 方 式 と は ，送 信 機 及 び 受 信 機 に

水平方向の分解能の約 4 km を達成するためにも，まシ ス テ ム 全 体 で 共 通 の 進 行 波 管 増 幅 器（TWTA：

た H-II ロケットに搭載可能なアンテナサイズ（約 3 m以下）を実現するためにも，13.8 GHz の周波数の選Travelling Wave Tube Amplifier）及び低雑音増幅器（LNA：Low Noise Amplifier）を用い，アレーア

択は妥当なものであったと考えられる．走査角（あるンテナの動作に必要な給電部の位相制御を受動素子

いは走査幅）は，大きいほどよいわけであるが，走査で行う方式のことをいう．移相器としては，5 ビット

角が大きくなると，アンテナビーム幅とパルス幅で決のフェライト移相器が選ばれた．パッシブアレー方

まるレンジビンの鉛直分解能が急激に劣化すること，

式は，構成が単純で軽量であるという利点がある反サイドローブから混入するグランドクラッタが増加す

面，TWTA の故障が直ちにレーダシステム全体の機ること，一つのアングルビン（アンテナビームが指向能低下に結び付くという欠点を有する．パルス TWTAする方向）方向の独立サンプル数を確保する必要があ

は，特殊な高圧の電源回路を必要とし，高圧対策及びることなどの条件から走査角はむやみには大きくす

他の機器に与える EMC 対策も重要な技術課題となることはできない．また，マイクロ波放射計で採用さる．アクティブアレー方式とは，アレー素子のそれぞれているコニカルスキャンは，独立サンプル数を確保

れに移相器と固体増幅器（SSPA：Solid State Powerするためには不利であり，結局アンテナビームは，鉛Amplifier），LNA の能動素子を組み合わせる方式をい725

電子情報通信学会論文誌 2008/7 Vol. J91–B No. 7

う．PR では，アレー素子数は 128 となるので，送信用増幅器として小型軽量で位相特性の安定した GaAs

表 2 TRMM 降雨レーダの主要諸元
Table 2 Main system parameters of PR.

FET の SSPA の使用を考えた．移相器としては，損失が約 5 dB と大きいものの，小型軽量な 5 ビットのPIN ダイオード移相器（PHS）を使用した．アクティ
ブアレーの場合，送受信の移相器の損失はそれぞれ後段に位置する SSPA 及び前段に位置する LNA で補償することができる．アクティブアレー方式は，素子数が多く，構成が複雑で重量が大きくなる反面，少数のアレー素子が故障しても，それがレーダシステム全体の大きな機能低下に結び付かないという利点を有する．
結局のところ，素子数が多く構成が複雑で重量が重くなるという欠点を有するものの，フェイルセーフ機能をもつアクティブアレー方式を採用することにした．

3. 2 TRMM 降雨レーダの概要
PR は，Feasibility Study の後，1988〜1991 年ごろにかけて当時の郵政省通信総合研究所において，導波管スロットアレーアンテナ，SSPA，LNA，PHS な
どの Critical Components の試作が行われ，1/16 のスケールの Bread Board Model（BBM）を製作した．
BBM による機能性能を確認後，Engineering Model，Proto Flight Model の開発のため，プロジェクトは当表 3 四つの観測モードにおけるレーダデータ収集パラ
メータの要約
Table 3 Summary of radar data collection parameters
in four observation modes.

時の宇宙開発事業団に引き継がれた．
宇宙開発事業団で開発した PR の主要諸元を表 2 に示す [7], [8]．PR は，周波数が 13.8 GHz 帯で動作する 128 素子の一次元アクティブフェイズドアレーレー
ダであり，5 ビットの PHS を制御して，電子的にアンテナビームを衛星の進行方向と直角方向に，地上軌跡にすき間が出ないように走査する．走査幅は，215 km
である．アンテナビーム中心の走査角は，±17◦ であ

高度のミラーイメージが 250 m のレンジ分解能で観測

り，アンテナビーム幅 0.71◦ を有するアングルビンがされる．ミラーイメージとは，地表面で散乱された電走 査 角 ±17◦ の 範 囲 に 49 個 存 在 す る．送 受 信 系 は，波が上空の降雨によって地表面方向に散乱され，それ128 素子の SSPA，LNA 及び PHS より構成されており，128 本の導波管スロットアンテナを並べた平面アが地表面によって再度散乱されて受信されるエコーを

レーアンテナと結合されている．アンテナ開口径の寸の位置に見かけ上現れる降雨エコーである．他の観測意味し，地表面（主に海面）をミラーとした鏡面対称法は，2.1 m × 2.1 m である．送信ピーク電力は，約モードには，レンジ分解能の半分の 125 m の間隔で，700 W である．高度 350 km より，雨域の頂きで降雨
強度 0.5 mm/h の降雨を観測することができる．直下点の水平距離分解能は，4.3 km，レンジ分解能は250 m である．四つの観測モードについて，表 3 に示限られた走査角の範囲において表面エコー並びに降雨に，周波数アジリティという手法を用いている．これす [7], [8]．通常，走査角 ±17◦ の範囲について，地表

は，13.796 GHz 及び 13.802 GHz の 6 MHz（パルス

エコーを詳細に観測する（オーバサンプル）モードがある．PR は，各独立サンプル数 64 を確保するため

面から高度 15 km の範囲の降雨並びに地表面エコー

幅 1.67 μs の逆数の約 10 倍）離れた 2 周波数のパル

を 250 m のレンジ分解能で観測する．ミラーイメージ

スのペア 32 組を各アングルビン方向に送信して，独の観測モードでは，直下点についてのみ，0〜5 km の立サンプル数 64 を確保するという方式である．降雨726

招待論文／熱帯降雨観測衛星搭載降雨レーダの研究開発とその観測成果

図 2 TRMM 降雨レーダの構成
Fig. 2 Block diagram of the PR.

エコーの受信強度は，空間内の降雨粒子の配置の変動によって一発一発の各パルスごとに変動するため，何図 3 TRMM 降雨レーダアルゴリズムフロー
Fig. 3 TRMM PR algorithm ﬂow.
表 4 TRMM 降雨レーダ標準アルゴリズム
Table 4 TRMM standard PR algorithms.

発かのパルスに対応する独立した降雨エコーを平均化することによって受信強度の標準偏差を軽減して平均受信強度を求めることができる．降雨の受信強度は指数分布に従うので，64 個の受信エコーの平均を行う
ことによって，受信エコーの標準偏差を約 0.7 dB 程度に抑えることができる．このために，TRMM では，
独立サンプル数 64 を確保している．また，同レーダのブロック図を図 2 に示す．PR は，アンテナ系，送受信系，信号処理系，構造系，熱制御系，計装系の六つのサブシステムから構成される．送信経路では，周波数変換・IF 部（FCIF）からの RF 信号は，分合波
器 1（DIV/COMB1）で 16 分割され，更に分合波器 2（DIV/COMB2）で 8 分割された後，128 分割された
信号が 128 本の導波管スロットアンテナに入力される．

PHS は，分合波器 2 に含まれている．受信経路は，送信経路の逆をたどる．信号処理系は，FCIF 部とシス
テム制御・データ処理部（SCDP）から構成される．

4. TRMM 降雨レーダデータ解析処理ア
ルゴリズム

性的な情報（ブライトバンドの有無，降雨タイプの分類）を与える．13.8 GHz の周波数は，降雨減衰の影
響を受けるため，降雨減衰の補正を行ったレーダ反射

PR データ解析処理用の標準アルゴリズムは，研究

因子及び降雨強度を算出する必要があり，このための公募によって選ばれた日米合同の TRMM PR サイエ

特別のアルゴリズムが整備されている．Level 3 は，ンスチームによって開発された [9]．アルゴリズム及びLevel 2 のプロダクトを時間的，空間的に平均し，緯そのプロダクトは，処理のレベルに応じて，Level 1，

度経度 5◦ × 5◦ の各格子における降雨パラメータの月

2，3 と分類される．Level 1，2 は，各瞬時視野方向の平均値を算出する．代表例は，月平均降雨強度である．

データである．Level 1 は，レーダ固有の工学値（受図 3 に，PR アルゴリズムシステムのフロー並びに相信電力値，ノイズレベル）及び降雨減衰を含んだレー

互のアルゴリズムの入出力関係を示す．また，各アルダ反射因子を与える．Level 2 は，降雨強度，地表面

ゴリズムの名称及びプロダクトを表 4 に示す．

の規格化散乱断面積（NRCS），及び降雨に関する定

1B21 は，レーダエコーのディジタル値を工学値に
727

電子情報通信学会論文誌 2008/7 Vol. J91–B No. 7

変換し，レーダ受信電力値（信号 + 雑音）及び雑音値

Technique）を用いている．衛星からレーダで降雨域（ノイズレベル）を与える．また，受信電力の大きさを観測すると，降雨の背景に地表面の強いエコーが受とあらかじめ定められたしきい値を比較して降雨の信される．同一入射角において降雨域と無降雨域の地有無を判定する．1C21 は，降雨減衰の影響を含んだ

表面のエコーの強さの差は，すべて降雨減衰によるもレーダ反射因子 Zm のレンジプロフィールを与える．

のと仮定して降雨減衰量を求め，制約条件とするのが2A21 は，降雨域におけるパス積分減衰値（PIA）及びその信頼度の値を与える．また，2A21 は無降雨時

表面参照法の考え方である．つまり，降雨によって，においては，地表面の規格化散乱断面積の値を算出す

表面参照法によって推定される降雨減衰量 P IA は，る．PIA は，2A25 において表面参照法（SRT）を用

地表面までの距離を rs とおくと，

地表面の散乱係数は変化しないと仮定するわけである．

いて降雨強度を算出するときに用いられる．無降雨時

の NRCS は，PIA 算出の基準として用いられる [10]．

P IA = exp −0.2 ln 10

 rs



k(s) ds(4)

0

2A23 は，ブライトバンドの有無を判定し，ブライト
バンドがあるときはその高さ，厚さ，レーダ反射因子

であり，これと降雨エコーから計算される地表面までを求める．降雨タイプは，レーダ反射因子 Z の鉛直の降雨減衰量，すなわち分布及び水平分布を考慮して，層状性降雨，対流性降雨，及びその他の降雨に分類される [11]．ブライトバZm (rs )= 1 − 0.2εβ ln 10
Ze (rs )ンドのあるときは，層状性降雨と分類している．また，

 rs

1/β
β

α(s)Zm (s) ds
0(5)

孤立した背の低い雨は，対流性の降雨と分類している．
が一致するように ε を選び，その ε を用いて，式 (3)2A25 は，降雨減衰補正を施したレーダ反射因子 Zeを推定し，その Ze から妥当と考えられる Z-R 関係
を用いて降雨強度 R に変換するという方針をとっている [12]〜 [14]．
降雨レーダから距離 r の点における降雨減衰を含ん

降雨減衰補正を必要とする最も遠い距離 rs での降雨

だレーダ反射因子 Zm は，降雨減衰係数 k を用いて，降雨減衰量そのものが少ない場合は，表面参照から推





 r

Zm (r) = Z(r) exp −0.2 ln 10

k(s) ds

より任意の距離 r における Ze を計算するものである．
減衰量を制約条件としているため，そこに至るまでの任意の距離での減衰量が安定に求められる．ただし，定させる P IA 自体の誤差が大きくてあてにならない(1)0こと，並びにそのように降雨減衰量が少ないときは，と書き表せる．ここで，Z は，降雨減衰を含まない

Hitschfeld-Bordan の解が不安定にならないことなどから，2A25 では実際には，推定された減衰の大きさ

レーダ反射因子の真値である．

に応じて表面参照法で得られた降雨減衰量の制限条件

2A25 では，Z を求める際に，k と Z の関係をの重みを変化させている．すなわち，式 (4) と式 (5)k = εαZ βを等しいとおいたときに得られる ε を ε0 と書くと，(2)のように仮定して，このときの Z = Ze とおいていε = (1 − w) + wε0(6)る．すなわち，式 (2) を式 (1) に代入し，Ze に関するで定義される ε を用いて，降雨減衰が少ない場合に，微分方程式に変換し，表面参照法の誤差に補正が引きずられることを防いでZe (r)=

Zm (r)

 rいる．このようにして補正された Ze に対して，妥当

1/β (3)α(s)Zm (s)β ds

1−0.2εβ ln 100と考えられる Z-R 関係を使って降雨強度を推定している．

を求めている．ここで ε = 1 が Hitschfeld-Bordan の3A25 は，1C21，2A21，2A23，2A25 の緯度経度
5◦ × 5◦ または，0.5◦ × 0.5◦ の領域における月ごとの解に相当する [15]．この解は，降雨減衰が大きいとき

統計値（度数，平均値，標準偏差，ヒストグラム，相には，分母が 0 に近づき不安定になり，実用にならな関係数）を計算する．この中で特に重要なものは，緯いことが多い．そこで，2A25 では，安定した降雨減

度経度 5◦ × 5◦ の領域の月平均降雨強度である．3A26

衰補正を行うために，表面参照法（Surface Reference

は，多重しきい値法という統計的手法を用いて，降雨728

招待論文／熱帯降雨観測衛星搭載降雨レーダの研究開発とその観測成果

強度分布を対数正規分布と仮定したときの緯度経度

三次元分布であり，台風の目を含む鉛直面内の降雨強

5◦ × 5◦ の領域の月平均降雨強度を 3A25 とは独立に算出する [16]．

度の鉛直分布を色別表示している．壁雲の降雨頂はスパイラルレインバンドに比べて高く，高度約 13 km

5. TRMM 降雨レーダの観測成果にまで達しているのが観測されている．1998 年の観TRMM 衛星が他の気象衛星と比べてユニークな点とである．PR によって海上の台風の中の降雨強度分のデータベースとして，JAXA/EORC の TRMM のHome Page（http://sharaku.eorc.jaxa.jp/TYP DB/
index j.shtml）で整備され，台風の構造に関する研究

布の水平分布，鉛直分布が初めて得られた．図 4 は，に貢献している．

は，降雨の三次元観測が可能な PR を搭載しているこ

測開始以来，PR が観測してきた多くの台風は，台風

アルゴリズム 2A25 の処理結果であり，海上の台風の

低緯度の大気大循環では，水の相変化に伴う潜熱加一例で，2006 年 8 月 30 日に，台風 12 号の目の上を熱が重要な役割を果たしている．このため，熱帯・亜熱TRMM 衛星がちょうど飛行している．左図は，高度2 km における降雨強度の水平分布が背景の TRMM
衛星搭載 VIRS の赤外雲画像に重ね合わせて表示されている．降雨レーダの走査幅 215 km が白い 2 本の帯域における潜熱加熱の四次元構造（三次元空間 + 時間）の推定が TRMM の主要な科学目的の一つであった [17]．降雨は，大気の激しい上昇下降運動に伴って強く降る対流性降雨と，水平に広がった雲から長時間線で示されている．降雨レーダはグランドクラッタのしとしとと降る層状性降雨に大別されるが，両者で潜影響や鉛直分解能の劣化のために広い範囲が走査でき

熱加熱の鉛直分布（プロファイル）の構造が大きく異ないが，可視・赤外画像では得られない雲の下の降雨なるため，潜熱加熱プロファイルの推定にとって，対域を明りょうにとらえて観測することができる．水平

流性・層状性降雨の分類は極めて重要である．PR は，

分解能は，直下点で約 4.9 km，走査の端で約 5.8 km

層状性降雨の特徴である融解層に伴うブライトバンドである．降雨強度の強さは色別に表示されている．台の検出が地上レーダよりも容易であるため，対流性・風の目を囲む壁雲並びに眼の外側のスパイラルレイ

層状性降雨の分類の精度が高く，潜熱加熱プロファインバンド（らせん状降雨帯）において，30 mm/h 以

ルの推定に適している．Tao らは，平均的な対流性・上の降雨強度が観測されている．右図は，降雨強度の

層状性加熱プロファイルを仮定し，PR から得られる

図 4

TRMM 降雨レーダによる台風の三次元観測例（2006 年 8 月 30 日）．JAXA/NICT/NASA の共同プロジェクトによる（JAXA 提供）．左図：高度 2 km における降雨強度の水平分布，右図：台風の目を含む鉛直面内の降雨強度の鉛直分布．

Fig. 4

An example of the three dimensional observation of the typhoon by TRMM
PR on 30 August 2006. Joint project by JAXA/NICT/NASA (Provided by
JAXA). Left image: Horizontal distribution of rain rate at the altitude of
2 km. Right image: Vertical distribution of rain rate in the vertical plane
including the eye of typhoon.

729

電子情報通信学会論文誌 2008/7 Vol. J91–B No. 7

対流性・層状性降雨比から潜熱加熱プロファイルを求PR と LIS を併用して計算される降雨量と発雷の比める Convective-Stratiform Heating（CSH）アルゴ（Rain-yeild flash：RPF）によって，これまでよく知リズムを開発した [18]．しかしながら，この手法は仮られてきたように降雨特性が大陸と海洋で分類される

定された加熱プロファイルに大きく依存し，地域的・だけでなく，両者の中間的な特徴をもつ遷移域が大陸

季節的に大きく変化する対流性加熱プロファイルの変周辺の海洋上に存在すること，また大陸上であっても化を考慮に入れることができない．

モンスーン期の降水システムは遷移域の特徴をもつこ

Shige らは，Takayabu による PR 降水プロファイルのスペクトル表示 [19] を基礎として，PR の鉛直情報を積極的に活用する Spectral Latent Heating（SLH）とを示した [26]．

アルゴリズムを開発した [20]．このアルゴリズムは，

6. 全球降水観測計画（GPM）
全球降水観測計画（GPM：Global Precipitation

数値モデルに基づくシミュレーションデータから作成

Measurement）は，1 機の主衛星と 8 機程度の副衛星

した参照テーブルに基づいて降水に伴う潜熱加熱の鉛により 3 時間ごとの高頻度で高精度の全球降水観測

直プロファイルを推定する．対流性・層状性降雨量比を行い，気候変動や地球規模の水循環に及ぼす降雨ののみから潜熱加熱プロファイルを求める CSH アルゴ

影響を解明するとともに，天気予報精度の向上，水資

リズムに比べて，源管理（洪水予測，河川管理，ダム貯水量の調節，農1. 降雨タイプや地上降雨強度だけでなく降雨の高

業用水の確保），農業生産性予測などの実利用分野へ

さ情報を用いているために，背の低い対流と背の高いの降雨観測データの利用実証を行うことを目的とし対流の潜熱加熱の違いが推定できる．

た，日米を中心とした国際協力によって計画されてい2. 層状性降雨に対しては融解層の降水強度を用いるミッションである [27]．主衛星の打上げは，2013 年ているために，地上無降雨でも潜熱加熱を推定できる．

ごろの予定である．GPM 主衛星は，TRMM を継承という特徴をもっている．

する衛星であり，日米共同で開発を行う予定である．

SLH アルゴリズムをグローバルな PR データに適用

GPM 主衛星には，2 周波降水レーダ（DPR）を搭載するために，アルゴリズムを改良した結果，地上観測

し，同時に搭載されるマイクロ波放射計 GMI を校正

値との良い一致が示された [21]．また，熱帯太平洋上するとともに，GMI のアルゴリズム開発を通して，副の潜熱加熱プロファイルを推定したところ，潜熱加熱

衛星搭載のマイクロ波放射計をも間接的に校正する

プロファイルが対流性・層状性降雨比だけでなく，対ことにより，複数の衛星データを総合し，高精度，高流性加熱プロファイルの変化に伴って大きく変化して頻度の全球降水マップを作成する．マイクロ波放射計

いることが示された．SLH アルゴリズムは，潜熱加熱を搭載した副衛星には，米国 NOAA の極軌道衛星，

量ばかりでなく水蒸気消失量の高度分布も推定するこ中国の FY-3（風雲 3 号），インド，フランス共同のとができ [22]，これらのデータは今後，気候モデルのMegha-Tropiques 衛星，我が国の GCOM-W 衛星な検証に役立つと思われる．

どが候補に挙がっている．

TRMM が打ち上げられるまで，降水の三次元デー
限定されたものしかなく，全球規模で地域的・季節的

GPM 主衛星の衛星高度は，約 407 km，軌道傾斜角は，65◦ である．またミッション期間は，設計寿命3 年，目標 5 年である．GPM 主衛星の軌道は，日周

特徴を調べることは不可能であった．PR によって降変化の大きな降水を様々な時間帯で均等に観測するた

水の三次元構造に関する「降水システム気候学」[23]

めに，非太陽同期軌道が選ばれる．また，できるだけが可能となり，対流性・層状性降雨量比 [24] や降水プ

高緯度までの DPR 観測を実現するとともに，DPR のロファイル [25] などの時空間的変動が詳細に調べら

狭い走査幅でも，できるだけ短期間で全球をカバーすれてきた．また，TRMM の特徴として，PR に加えるという要求条件から上記の軌道傾斜角が選ばれた．

て TMI，VIRS，LIS，CERES といった五つの異なGPM 主衛星搭載の DPR は，TRMM 降雨レーダ

タは地上レーダの特別観測による空間的・時間的にるセンサによって同時に降水システムを観測している

PR を継承するアクティブフェイズドアレーレーダ

点が挙げられる．近年，これら複数のセンサを融合的であり，TRMM PR を高性能化，軽量化した Ku 帯に利用して，PR 単独では得られない情報を引き出す（13.6 GHz）降雨レーダ（KuPR）に加え，高感度化，試みが盛んに行われている．例えば，Takayabu は，高精度化を目指した新しい Ka 帯（35.5 GHz）降雨730

招待論文／熱帯降雨観測衛星搭載降雨レーダの研究開発とその観測成果
表 5 GPM 2 周波降水レーダの主要諸元
Table 5 Main system parameters of GPM DPR.

しかし，TRMM PR は，非常に部品数が多く，重量が非常に重い．また，アンテナは狭帯域，単偏波用であり他の周波数，偏波との共用は不可である．将来の降水観測ミッションにおいては，より詳細な降水雲の構造，降水プロセスを観測するために，降水観測用の

13.6，35.5 GHz に加えて雲観測用の 95 GHz 帯の周波数並びに水平・垂直偏波を用いた同時観測が可能な
マルチパラメータレーダに対する要望が出てくるものと思われる．このために，更に次世代をにらんだ降雨
レーダの研究を継続することが期待される．このためには，シリンドリカルパラボラアンテナを用いた電子
走査型の降雨レーダ [28] や，原点に戻って，Aqua 搭
載の AMSR-E のようなパラポラアンテナを円錐走査するタイプのレーダの可能性も含めて可能なシステムを検討する必要があろう [29]．気象庁が現在行っている数値予報において，TRMM TMI をはじめとするマ
イクロ波放射計からの観測データが数値予報の初期値作成のためのデータ同化に利用されているが，PR のデータは利用されていない．これは，TMI の走査幅（760 km）に比べて PR の走査幅が 215 km と狭いかレーダ（KaPR）から構成される．KaPR は，KuPRらである．PR のアンテナ走査角を広くとると，鉛直では測れない弱い雨や雪の検出に有効であり，両周波分解能が劣化し，様々な高度の降雨を同時に受信し，数帯のレーダの異なるダイナミックレンジによって，

高度方向の降雨の非一様性の影響を受けること，アン

熱帯から高緯度までの降雪を含む降水量を高精度で観テナサイドローブからの地表面クラッタの影響を受け測することができる．また 2 周波のアンテナビームをることなどの問題が発生する．アンテナの走査角を広マッチングさせ，送信パルスタイミングを一致させる

げた場合の影響についての研究は最近行われ始めてことによって，2 周波同時観測が可能となり，2 周波

いる [30], [31]．実際の TRMM は，進行方向に垂直方

アルゴリズムによって雨と雪の識別，雨滴粒径分布の向にアンテナビームを走査しているが，この方向の走推定が可能となり，高精度の降雨強度推定が可能とな

査角は ±17◦ と狭いため，具体的な検討においては，る．DPR の主要諸元を表 5 に示す．

TRMM の進行方向に長いデータを切り出し，データを 90◦ 回転させ，このデータを真値として用いてアン

7. む す び

テナビームを広く走査したときの観測距離の増加によ

TRMM 降雨レーダは，世界初の衛星搭載降雨レーるアンテナビームの広がりやグランドクラッタの混入ダであり，我が国が開発したものである．TRMM 衛

によって，どの程度観測データが真値と比べて劣化す

星は，1997 年 11 月 28 日以降，現在に至るまで約 10

るのかをシミュレーション計算を行って評価している．

年間，何らの不具合もなく観測を続けている．アクアンテナは TRMM と同様な導波管スロットアンテナ

ティブフェイズドアレー方式のレーダであるため，機を仮定している．評価では，真値とシミュレーション

械的な可動部分がないこと，固体タイプの送受信素子

値の相関係数や RMS 誤差を計算し評価している．予を使ったレーダであり，TWTA のような電子管に比備的な検討では，走査角を 30 度程度に広げた場合でべて寿命が長いことがその理由であろう．TRMM のも相関係数は 0.8 程度の値が得られており走査角を広成功を受けて，GPM 搭載の DPR も TRMM と同タ

げられる可能性があることが報告されている．しかし，

イプのレーダが搭載される．宇宙分野においては，信アンテナ走査角を広げた場合には，降雨域と衛星の距

頼性と実績が重視されるゆえ，うまく動作している

離が増加し，TRMM と同様に 0.5 mm/h 程度の弱いものを他のものに置き換えることは通常敬遠される．

降雨強度を観測するには所用送信ピーク電力が増大し，
731

電子情報通信学会論文誌 2008/7 Vol. J91–B No. 7

その結果消費電力が増加する問題点がある．また，独立サンプル数 64 を確保するには，GPM で用いられ

1988.
[8]

T. Kozu, T. Kawanishi, H. Kuroiwa, M. Kojima, K.
Oikawa, H. Kumagai, K. Okamoto, M. Okumura,ている PRF を走査角ごとに可変にする手法を用いて

H. Nakatsuka, and K. Nishikawa, “Development of

も走査角があまり広げられない問題点がある．更にグ

precipitation radar onboard the Tropical Rainfall

レーティングローブを避けるために導波管の間隔を狭

Measuring Mission (TRMM) satellite,” IEEE Trans.

くする必要があり，導波管の本数が増えて，重量が重Geosci. Remote Sens., vol.39, no.1, pp.102–116, Jan.

くなるような技術的な問題もある．これらの問題点を

2001.
[9]

解決するための今後のシステム検討の研究結果に期待

岡本謙一，石戸喜夫，R. Meneghini，“TRMM 降雨レー
ダアルゴリズム，” 日本リモートセンシング学会誌，vol.18,
no.5, pp.40–51, Dec. 1988.

したい．
衛星搭載降雨レーダは，TRMM，GPM と我が国が[10]

R. Meneghini, T. Iguchi, T. Kozu, L. Liao, K.
Okamoto, J.A. Jones, and J. Kwiatkowski, “Use of

世界をリードして開発を進めてきた．引き続き，次世

the surface reference technique for path attenua-

代の衛星搭載降雨レーダの開発においても，我が国がtion estimates from the TRMM precipitation radar,”

世界をリードして，より先進的な降雨レーダを実現す

J. Appl. Meteor., vol.39, no.12, pp.2053–2070, Dec.
2000.

ることを希望したい．

[11]

J. Awaka, T. Iguchi, and K. Okamoto, “Rain type

TRMM 降雨レーダアルゴリズムについては，やはclassiﬁcation algorithm,” in Measuring Precipitation

り我が国が主導権をもって日米協力のもとでこれまでfrom Space — EURAINSAT and the Future, ed.

開発・改良を続けてきた経緯がある．現在の Version 6

V. Levizzani, P. Bauer, and F.J. Turk, pp.213–224,

において，雨量計と比較したとき，陸上の降雨強度を

Springer, Dordrecht, The Netherlands, 2007.
[12]

20〜50%過小評価している問題点が発生している．現在の PR アルゴリズムは，1 周波数レーダアルゴリズ

Okamoto, “Rain-proﬁling algorithm for the TRMM
precipitation radar,” J. Appl. Meteor., vol.39, no.12,
pp.2038–2052, Dec. 2000.

ムとしては，限界にきていると思われるが，来年の10 月を目指して陸上の過小評価を緩和する Version 7

T. Iguchi, T. Kozu, R. Meneghini, J. Awaka, and K.

[13]

T. Iguchi and R. Meneghini, “Intercomparison of
single-frequency methods for retrieving a vertical rainを開発している最中である．レーダのハードウェア

proﬁle from airborne or spaceborne radar data,” J.

開発とアルゴリズム開発は，不離一体である．ぜひ，Atmos. Oceanic Technol., vol.11, no.6, pp.1507–1516,

GPM の DPR の 2 周波アルゴリズムの開発 [32] においても，世界をリードするアルゴリズムの開発を期待
したい．
文
[1]

[14]

Dec. 1994.
井口俊夫，R. Meneghini，“降雨レーダによる降雨強度推

[15]

定，
” 月刊海洋，vol.31, no.6, pp.324–329, 1999.
W. Hitchfeld and J. Bordan, “Errors inherent in theradar measurement of rainfall at attenuating wave-

献

J. Simpson, ed., Report of the science steering groupfor a tropical rainfall measuring mission (TRMM),

lengths,” J. Merteor., vol.11, pp.58–67, 1954.
[16]

and J. Kwiatkowski, “Statistical methods of estimat-

NASA/Goddard Space Flight Center, 94p., 1988.
[2]

田中 佐，“TRMM 計画の概要，” 日本リモートセンシン
グ学会誌，vol.18, no.5, pp.13–27, Dec. 1988.

[3]

岡本謙一，“熱帯降雨観測衛星，” 信学技報，SANE 2006-[4]

K. Okamoto, ed., “A feasibility study of rain radar for

115, Dec. 2006.

ing average rainfall over large space-timescales using
data from the TRMM precipitation radar,” J. Appl.
Meteor., vol.40, no.3, pp.568–585, March 2001.
[17]

[5]

K. Nakamura, K. Okamoto, T. Ihara, J. Awaka, andsion (TRMM),” Meteor. Atmos. Phys., vol.60, no.1,
pp.19–36, 1996.
[18]

proﬁles of latent heat release:tropical rainfall measuring mission,” Int. J. Satellite
Communications, vol.8, no.3. pp.257–268, May 1990.
岡本謙一，“熱帯降雨観測衛星（Tropical Rainfall Measuring Mission，TRMM）計画における降水測定システ
ムの研究開発，” 天気，vol.41, no.7, pp.361–378, July732

古津年章，黒岩博司，“TRMM 搭載降雨レーダ，” 日本リ
モートセンシング学会誌，vol.18, no.5, pp.28–39, Dec.

Their applications

for TRMM,” J. Meteor. Soc. Japan, vol.71, no.6,
pp.685–700, 1993.
[19]

Y.N. Takayabu, “Spectral representation of rain features and diurnal variations observed with TRMM
PR over the equatorial area,” Geophys. Res. Lett.,

1994.
[7]

W.K. Tao, S. Lang, J. Simpson, and R. Adler,“Retrieval algorithms for estimating the vertical

T. Kozu, “Conceptual design of rain radar for the[6]

J. Simpson, C. Kummerow, W.-K. Tao, and R.F.
Adler, “On the Tropical Rainfall Measuring Mis-

the tropical rainfall measuring mission,” J. Commun.
Res. Lab., vol.35, no.145, pp.109–208, July 1988.

R. Meneghini, J.A. Jones, T. Iguchi, K. Okamoto,vol.29, doi:10.1029/2001GL014113, 2002.
[20]

S. Shige, Y.N. Takayabu, W.-K. Tao, and D.E.

招待論文／熱帯降雨観測衛星搭載降雨レーダの研究開発とその観測成果
Johnson, “Spectral retrieval of latent heating pro-[30]

ﬁles from TRMM PR data. Part I: Development ofa model-based algorithm,” J. Appl. Meteor., vol.43,[31]

Dec. 2006.
渡邉圭太，岡本謙一，重 尚一，“次世代の衛星搭載降雨

S. Shige, Y.N. Takayabu, W.-K. Tao, and C.-L. Shie,

レーダのシステム設計，” 第 44 回日本航空宇宙学会関西支“Spectral retrieval of latent heating proﬁles from
[32]

部・中部支部合同秋季大会講演集，pp.57–60, Nov. 2007.
K. Nakamura and T. Iguchi, “Dual-wavelength radar

no.8, pp.1095–1113, Aug. 2004.
[21]

TRMM PR data. Part II: Algorithm improvement[22]

渡邉圭太，岡本謙一，重 尚一，“次世代の衛星搭載降雨レーダのシステム設計，” 信学技報，SANE 2006-115,

and heating estimates over tropical ocean regions,” J.

algorithm,” in Measuring Precipitation from Space

Appl. Meteor. Climatol., vol.46, no.7, pp.1098–1124,— EURAINSAT and the Future, ed. V. Levizzani,

July 2007.

P. Bauer, and F.J. Turk, pp.225–234, Springer, Dor-

S. Shige, Y.N. Takayabu, and W.-K. Tao, “Spec-drecht, The Netherlands, 2007.

tral retrieval of latent heating proﬁles from TRMM

（平成 19 年 11 月 30 日受付，20 年 3 月 11 日再受付）

PR data. Part III: Estimating apparent moisturesink proﬁles over tropical oceans,” J. Appl. Meteor.
Climatol., vol.47, no.2, pp.620–640, Feb. 2008.
[23]

全球降水観測（GPM）計画—科学的・社会的期待，宇宙研究開発機構地球観測研究センター，2007.

[24]

C. Schumacher and R.A. Houze, Jr., “Stratiform rainin the tropics as seen by the TRMM precipitation
radar,” J. Climate, vol.16, no.11, pp.1739–1756, June

[26]

謙一 （正員：フェロー）

昭 48 東大大学院理学系研究科博士課程
了．同年郵政省電波研究所（現，情報通信

M. Hirose and K. Nakamura, “Spatial and seasonal研究機構）入所．地球環境計測部長，標準計測部長，総合研究官を経て，平 12 退所．

variation of rain proﬁles over Asia observed by space-

平 12〜20 大阪府大大学院工学研究科教授．

borne precipitation radar,” J. Clim., vol.15, no.23,
pp.3443–3458, Dec. 2002.

平 20 鳥取環境大学環境情報学部教授．理博．平 14 武田賞など受賞．人工衛星搭載電波リモートセンサ

Y.N. Takayabu, “Rain-yield per ﬂash calculated fromを用いた地球環境計測に関する研究に従事．

2003.
[25]

岡本

TRMM PR and LIS data and its relationship to thecontribution of tall convective rain,” Geophys. Res.
Lett., vol.33, L18705, doi:10.1029/2006GL027531,
2006.
[27]

E.A. Smith, G. Asrar, Y. Furuhama, A. Ginati, A.
Mugnai, K. Nakamura, R.F. Adler, M.-D. Chou, M.
Desbois, J.F. Durning, J.K. Entin, F. Einaudi, R.R.

重尚一

平 13 京都大学大学院理学研究科博士後期課程了．同年宇宙開発事業団（現宇宙航空研究開発機構）招聘研究員，平 16 大阪

Ferraro, R. Guzzi, P.R. Houser, P.H. Hwang, T.

府大大学院工学研究科助手，平 19 同助教．
博士（理学）．気象・気候学，地球観測に関Iguchi, P. Joe, R. Kakar, J.A. Kaye, M. Kojima,する研究に従事．

C. Kummerow, K.-S. Kuo, D.P. Lettenmaier, V.
Levizzani, N. Lu, A.V. Mehta, C. Morales, P. Morel,
T. Nakazawa, S.P. Neeck, K. Okamoto, R. Oki,
G. Raju, J.M. Shepherd, J. Simpson, B.-J. Sohn,
E.F. Stocker, W.-K. Tao, J. Testud, G.J. Tripoli,
E.F. Wood, S. Yang, and W. Zhang, “International
global precipitation measurement (GPM) programand mission: An overview,” in Measuring Precipitation from Space — EURAINSAT and the Future, ed.
V. Levizzani, P. Bauer, and F.J. Turk, pp.611–615,
Springer, Dordrecht, The Netherlands, 2007.
[28]

Y. Rahmat-Sami, J. Huang, B. Lopez, M. Lou, E.
Im, S. Durden, and K. Bahadori, “Advanced precipitation radar antenna: Array-fed oﬀset membrane
cylindrical reﬂector antenna,” IEEE Trans. Antennas
Propag., vol.53, no.8, pp.2503–2515, Aug. 2005.

[29]

大利達也，岡本謙一，重 尚一，真鍋武嗣，“衛星搭載簡” 第 44 回日本航空宇宙
易型降雨レーダのシステム設計，学会関西支部・中部支部合同秋季大会講演集，pp.61–64,
Nov. 2007.

733