発表要旨集
平成29年2月28日（火） 13時30分～16時30分
一橋大学 一橋講堂

気象庁気象研究所

ごあいさつ本日は、気象研究所研究成果発表会にお越し頂き、誠にありがとうございます。
気象研究所は、気象庁の施設等機関として、気象庁が発表する各種情報の改善に役立つ研究や、気象業務の将来を見据えた基盤的な研究など、幅広い研究に取り組み、
気象業務の発展に貢献しております。平成 26 年度から「気象研究所中期研究計画」
を定め、防災情報や地球環境関連情報の一層の精度向上を目指す研究として、３つの柱「台風・集中豪雨等対策の強化に関する研究」「地震・津波・火山対策の強化に関する研究」「気候変動・地球環境対策の強化に関する研究」を立てて、研究を行っています。
こうした私どもの研究活動を広く皆様方にお知らせして、気象研究所をより身近な存在として感じて頂くために、平成 15 年度からこの研究成果発表会を実施しております。今回は、「自然災害に立ち向かう研究の最前線」というテーマで、豪雨・台風の予測、地震発生状況の即時把握や火山活動評価の高度化といった自然災害に関する最先端の研究をご紹介いたします。
最初の発表では、当所で行っている自然災害の軽減に向けた研究につきまして、中期研究計画の研究課題や課題間の関連といった概要をご案内いたします。
２件目は、GNSS（全球測位衛星システム）を使って、海上で水蒸気観測を行うための研究について紹介いたします。周囲を海に囲まれた我が国では、海上からの大量の水蒸気の流入が豪雨発生に重要な役割を果たします。海上での水蒸気観測密度を高めることで、気象予報の精度向上に寄与することが期待されています。
３件目は、気象衛星ひまわり 8 号のデータから算出される風のデータを使用して、台風の強度・構造変化を捉える研究について紹介いたします。2015 年 7 月に運用が開始されたひまわり８号は、大幅な機能の向上が図られており、特に高頻度の観測によって、風の変化が大きな台風領域内でも、時間・空間的に密にデータを得ることができるようになりました。この研究は、将来のより詳細な台風解析の業務への応用が期待されます。
４件目は、伊豆大島での地殻変動観測によるマグマの蓄積を捉える研究について紹
介いたします。伊豆大島火山が現在どのような状態にあるのか、火山活動評価の業務に役立つことが期待されます。
最後の５件目は、地震多発時においても適切に震源を決定する自動処理手法の開発に関する研究について紹介いたします。大規模な地震がひとたび発生すると、それに伴い数多くの地震が発生します。本研究では、いままで、人の手を介して実施していた詳細な震源・マグニチュードの決定を自動的に行う手法について研究しています。
すでに一部は気象庁で業務化されており、本日お話しする昨年 4 月に発生した熊本地震の際も大きな威力を発揮したと高く評価されています。今後、さらに地震活動の推移把握の高度化や、巨大地震の地震像の即時把握にも貢献することが期待されます。
今回の研究成果発表会を通じて、気象研究所の活動についてさらにご理解を賜りますとともに、今後の一層のご支援をお願いいたします。
気象研究所長高野清治

プログラム
國本 美華 氏（気象キャスター）
13:30-13:40 開会の挨拶
気象研究所長13:40-13:50 気象研究所での自然災害の軽減に向けた研究の概要
研究総務官 齊藤 和雄
13:50-14:20 豪雨予測への貢献をめざした海洋での GNSS 水蒸気観測
気象衛星・観測システム研究部 室長 小司 禎教
14:25-14:55 気象衛星ひまわり８号の大気追跡風を用いて捉えた台風の強化
台風研究部 主任研究官 小山 亮

15:00-15:10 休 憩

15:10-15:40 伊豆大島で進行するマグマ蓄積 ～火山性地殻変動は何を明らかにしたか～
火山研究部 室長 山本 哲也
15:45-16:15 熊本地震 ～地震発生状況の即時把握技術の開発～
地震津波研究部 研究官 溜渕 功史
16:20-16:30 閉会の挨拶
気象研究所研究成果発表会実行委員長

気象研究所での自然災害の軽減に向けた研究の概要○齊藤和雄（気象研究所研究総務官）

１．はじめに気象研究所は、気象庁に求められる諸課題に対して、科学

GNSS水蒸気観測と気象衛星ひまわり８号の大気追跡風を用いた台風強度解析について、発表する。

技術に立脚した気象業務への研究開発による貢献を目指し、三つの研究分野（台風・集中豪雨等、地震・火山・津波、及び気候変動・地球環境）の対策強化に関する研究を実施している。研究遂行にあたって、５年以内に業務化のめどをつける
問題解決型の研究・技術開発として「重点研究課題」を、また５年～10年後をめどとした実用化をめざす基盤的な研究・技術開発として「一般研究課題」を設定している。基本方針として研究開発力を向上し、研究連携を促進させ、成果発信を行
い、WMO等国際機関や学術分野への貢献と国民の防災・減第２図：台風・集中豪雨等の研究と気象業務の関連

災・安心安全など国内貢献に努めている（第１図参照）。
研究成果発表会は、気象研究所の成果発信の一つとして、

３．地震・火山・津波研究

一般の方々にその研究成果を知って頂くための場として年一地震・火山・津波研究分野では、B1緊急地震速報の予測回開催している。今年度は、自然災害軽減に向けた研究に関手法の高度化に関する研究、B2地震活動・地殻変動監視のして、研究成果を紹介することにした。

高度化に関する研究、B3津波の予測手法の高度化に関する研究、B4大規模噴火時の火山現象の即時把握及び予測技
術の高度化に関する研究、B5地殻変動観測による火山活動
評価・予測の高度化に関する研究、B6南海トラフ沿いのプレート間固着状態監視と津波地震の発生状況即時把握に関する研究、B7火山ガス観測による火山活動監視・予測に関する研究、の７つの重点研究を行っている。地震、津波及び火山に関する防災情報をより的確なものとし、それらによる災害を防止・軽減するため、発生した現象の推移をより的確に観測・解析する技術を開発するとともに、予知・予測技術の高精度

第1図：気象研究所での研究の基本方針

化を進める研究である（第３図）。
本日の講演会では、それらの研究の中から、伊豆大島で２．台風・集中豪雨等の研究
台風・集中豪雨等の研究分野では、A1メソスケール気象予のマグマ蓄積の監視、熊本地震での地震発生状況の即時把握に関する研究について発表する。

測の改善と防災気象情報の高度化に関する研究、A2顕著現象監視予測技術の高度化に関する研究、A3台風の進路予報・強度解析の精度向上に資する研究、A4沿岸海況予測技術の高度化に関する研究、の４つの重点研究課題と一般研
究課題a5大気境界層過程の乱流スキーム高度化に関する研究、を行っている。台風・集中豪雨等に関する予報・警報等の防災気象情報をより的確なものとし、それらによる災害を防止・軽減するため、災害をもたらす現象に対する観測・解析技
術及び予測技術を高度化する研究である（第２図）。
本日の講演会では、それらの研究の中から、海洋での第３図：地震・火山・津波研究と気象業務の関連

豪雨予測への貢献をめざした海洋でのGNSS水蒸気観測○小司禎教（気象衛星・観測システム研究部）

１．はじめにしばしば豪雨や突風など激しい気象現象を引き起こす積乱雲は、地上数百m～1km程度の高度に暖湿な空気が集まることで発生・発達するため、風と水蒸気を、高い時・空間分
解能で観測することが機構解明や予測に重要である。
GPS(Global Positioning System: 全地球測位システム)に代表
さ れ る 人 工 衛 星 を 用 い た 測 位 シ ス テ ム (GNSS: Global
Navigation Satellite System)は、カーナビゲーションやスマー
トフォンの位置情報アプリケーションとして日常生活に浸透している。あまり知られていないが、GNSSは全天候型の高精

第１図： GNSSによる水蒸気量解析の概念を説明する模式

度の水蒸気センサーという側面も有する。国土地理院は図。GNSS解析による可降水量は、地上受信局を頂点とする

1,300 箇 所 以 上 （ 平 均 20km 間 隔 ） か ら な る GNSS 観 測 網

逆円錐空間の代表値である。より多く、ノイズの少ない衛星(GEONET: GNSS Earth Observation Network System)を展開からの電波を受信することが、測位や水蒸気の解析精度向し、日本列島の地殻変動を監視しているが、GEONETから得上をもたらす。

られる水蒸気情報は、気象庁の天気予報にも実利用されている。
周囲を海に囲まれた日本列島では、海上からの湿潤な空気の流入が豪雨発生に重要な役割を果たすため、海上での水蒸気監視能力等の高度化が求められている（国土交通省,
2015）。ここでは、近年のGNSS技術の進展により可能となりつつある、船舶に搭載したGNSSによる海洋上での水蒸気連続観測技術の研究成果を報告する。

とで水蒸気による寄与のみを抽出することができ、さらには水
蒸気量を求めることができる。すなわち、測位解析にとってはノイズとなる水蒸気量を、大気変動のシグナルとして利用することができる。気象庁では各衛星方向の水蒸気量を天頂方向の値に換算した、アンテナを頂点とした逆円錐空間内の代表的な鉛直積算水蒸気量（可降水量、 Precipitable Water
Vapor: PWV）等を天気予報に利用している。この場合、衛星の位置（軌道）と時計情報は、国際GNSSサービス等外部機

２．GNSSを用いた水蒸気解析（第１図）
GNSS による位置計測（測位）では、原子時計を搭載し、高度約 2 万 km を周回する衛星群から発射される電波が、受信機に到達するまでの時間を利用して GNSS 衛星と受信機間の距離を計測する。GNSS 衛星の位置が既知で時計が正確であれば、3 つ以上の衛星からの距離を計測することで受信機の三次元的な位置が求められる。

関の解析結果を利用し、アンテナは解析期間（例えば24時間）中動かないと仮定して解析を行っている。

３．GNSS技術の革新と水蒸気解析
2010年、日本初の衛星測位システムである準天頂衛星(Quasi-Zenith Satellite System: QZSS)初号機（みちびき）が打ち上げられ、2011年にはロシア がGLONASS(Globalnaya

電波は真空中では光速で伝播するが、大気中では速度がNavigatsionnaya Sputnikovaya Sistema)の測位サービスを開遅くなる。遅れの量は経路上の気圧・気温・水蒸気によって決始した。測位衛星数の増加は、固定点に加え、移動体測位のまる。気圧や気温に比べ、水蒸気は時・空間変動が大きく、高精度化・安定化につながっている。また近年のGNSS技術

GNSS測位に誤差をもたらすノイズである。一方、衛星と受信の高度化により、リアルタイムによる高頻度（１秒）測位解析

機の位置を予め求めておけば、衛星から発射された電波が、が行われるようになり、国土地理院による地殻変動観測のた

地上の受信機に到達するまでに経過した時間から、電波のめのGEONETリアルタイム解析システム(REGARD)の開発等

遅れの量（遅延量）がわかる。気圧や気温の寄与を差し引くこに活用されている。
GNSS電波を用いた水蒸気量解析は、衛星からの電波の遅れを利用するため、衛星の位置・時計誤差、GNSS受信機の位置も高精度に解析する必要がある。陸上固定の場合、アンテナは動かないものとしているのに対し、船舶などの移動体に設置した受信機の場合、時間変動する電波の遅れに加え、受信機の位置も同時に推定する必要がある。従来は事後解析による研究が行われてきたが、上述の様な解析技
術の高度化により、移動体におけるリアルタイムでのPWV解析を行える環境が整ってきている。

４．船舶搭載GNSSによる海上でのPWV解析実験
平成27年、4隻の船舶にGNSS受信機を設置し、洋上観測実験を行った（第3図）。移動する受信機位置推定のため、GPS衛星に加え、GLONASS、QZSSを利用し、衛星の軌道情

第２図：準天頂衛星打ち上げ（2010年）以後に実現したGNSS
技術の革新。

報としてMADOCA 1)によるリアルタイムプロダクトを適用した解析を試みた。さらに、設置場所に制限のある船舶上で影響の大きい反射波による誤差を含む解析結果を除去する品質管理手法を考案した。
解析の結果、近傍の陸上GNSS点で解析された水蒸気量に匹敵する精度で解析できることが確認できた（第4図）（Shojiet al. 2016）。
このことは、位置が固定しない海上の船舶や浮遊ブイに搭載したGNSSアンテナで、高精度のPWV解析が可能であることを示している。
実用化に向けた課題として、海洋上でのデータ通信がある。
解析する計算機にGNSS観測データ（60MB/日)、衛星軌道情報（150MB/日）を送る必要がある。QZSSは衛星軌道情報を配信する機能を有し、GNSS受信機で受信することができる。

第３図：船舶に設置したGNSSアンテナの一例（関釜フェリー
株式会社の「はまゆう」）。デッキ上部左舷側ポール（赤丸で囲んだ場所）に取り付け、下関と釜山を結ぶ航路上の可降水量を解析した。

4基体制となる2018年には、日本周辺で常時1基以上のQZSSが観測できる。船上で受信したGNSSデータ、衛星軌道情報を用い、船上でPWVが解析できるため、通信データ量の課題
はクリアされると期待される。

参考文献（１） 国土交通省交通政策審議会気象分科会 2015: 提言「新たなステージ」に対応した防災気象情報と観測・予測技術のあり方」，
http://www.jma.go.jp/jma/press/1507/29a/bunkakai150729_2.
html（２）Shoji Yoshinori, Kazutoshi Sato, Masanori Yabuki, and Toshitaka
Tsuda 2016: PWV Retrieval over the Ocean Using Shipborne
GNSS Receivers with MADOCA Real-Time Orbits, SOLA, 12,265-271.
1)

Multi-GNSS Advanced Demonstration Orbit and Clock
Analysis の略称。宇宙航空研究開発機構(JAXA)の開発する
複数 GNSS 対応衛星軌道時刻推定ツール。

第４図：フェリー「はまゆう」で解析された PWV と、近傍（20km以内）の地上 GNSS 点で解析された PWV との比較（赤点：0079（下関）、緑点：0456（上対馬））。

気象衛星ひまわり8号の大気追跡風を用いて捉えた台風の強化〇 小山 亮 （台風研究部）

雲頂付近での風の吹き出し（アウトフロー）１．はじめに台風は、暖かい海面からの水蒸気及び熱をエネルギー源として発達し、その暴風・大雨により甚大な被害をもたらす。
台風の発達、すなわち低気圧性循環（一次循環）の強化には、対流によって水蒸気が凝結することで解放される潜熱と、こ雲頂付近では、中心付近は低気圧性循
環、中心から離れた場所では高気圧性循環が卓越する。

れに伴って強まる鉛直循環（二次循環）（第１図）が重要であることが、理論的、数値的研究によって知られている。一方、地上付近での風の吹き込み（インフロー）

ゾンデ観測をはじめとした現場観測では、時間・空間的に密な観測に基づく台風の発達過程の検証は難しい。このような

青矢印：一次循環、 赤矢印：二次循環、ピンク領域： 暖気核

背景もあり、気象庁他の予報センターの熱帯低気圧強度解

紙面奥から手前に向かう流れ、

析は、主に衛星で観測された雲の形状（眼の明瞭さなど）にも紙面手前から奥に向かう流れ

第１図：台風内部の風の模式図とづいたドボラック法による強度推定に頼っている。
大気追跡風（Atmospheric Motion Vectors: AMV）は、連続する衛星画像から雲及び水蒸気パターンの移動量を算出し、雲の雲頂温度から数値予報モデルの気温プロファイルを参照して高度を求めて得られる風データである。AMVは、現場観測データの少ない海上を含め、広い範囲で定常的に得られる利点があるが、静止気象衛星MTSAT（ひまわり6、7号）、

第２図：2015年7月7日に運用が開始された次世代静止気及び2015年7月7日に運用が開始された次世代静止気象衛

象衛星ひまわり8号（気象庁HP）

星ひまわり8号（第２図）の高頻度観測（2.5～15分間隔）を用いると、風の変化が大きな台風領域内でも、時間・空間的に

27N

接線風

密にデータを得ることができる（小山 2015）。
本研究の目的は、巻雲等の上層雲を追跡して算出される

動径風（アウトフロー）

26N

AMV（第３図）を使用して、台風の雲頂高度付近の一次循環と二次循環を捉え、これらと台風発達との関係を定量的に明らかにすることである。上層AMVを用いた台風への新しいア

25N
100km
24N

プローチによって、台風の発達プロセスに関する新たな知見が得られるとともに、台風解析業務への貢献も期待される。

200km
23N
123E

２．データ・手法

124E

300km
125E

126E

127E

第３図：2015年台風第15号（Goni）の8月23日18UTCのひま

MTSATの北半球領域画像から6時間毎に算出したMTSAT

わり8号ターゲット観測の水蒸気画像（バンド10）から算出した

AMV と 、 ひ ま わ り 8 号 タ ー ゲ ッ ト 観 測 （ 2.5 分 間 隔 、上層AMV（青矢羽、高度は100～300hPa面）。

1000km×1000kmの領域）から10分毎に算出したひまわり8号
AMV（Shimoji 2014）を使用した。AMVの高度は100～300hPa

３．上層AMVが捉えた2013年台風第24号の強化面である。一次循環、二次循環の強さは、半径300km内の各第４図に、MTSAT AMVを使用して行なった、2013年台風第半径で求めた上層AMVの接線風、動径風（アウトフロー）の24号（Danas）の発達期の上層の風の場の解析結果を示す。

最大値として、それぞれ定義した（第３図）。台風の最大風速

Danasは、最大12.8ms-1day-1の発達率（最大風速の前24時間（地上）は、気象庁ベストトラックデータを使用した。

変化）で、最大風速46ms-1 まで発達した。発達期において、

上層最大接線風の時間変化の特徴をより詳しく調べるため、

まり（対流バースト）が2回起こっていたが、これらの対流バー

10分毎のひまわり8号AMVを用いた事例解析を行なった。第

ストの終了頃～6時間後にかけて、接線風と動径風の顕著な

５図に、2016年台風第10号（Lionrock）の解析例を示す。上層

強化が起こっていた。この上層の風の場にみられた特徴的な

最大接線風の時間変化には、雲の形状に基づくドボラック法

変化は、強い対流によって地上付近の大きな角運動量が上で推定した最大風速とベストトラック最大風速からは解析され層まで運ばれたこと、対流に伴う潜熱解放によって、台風中

ない、一日以下の時間スケールの変化が明瞭に解析されて心付近に周囲より気温が高い領域（暖気核）が形成され、構いる。この上層最大接線風強化のタイミングは、極軌道衛星

造が軸対称化したことを示唆している（Oyama et al. 2016）。

NOAA及びMetOpのマイクロ波探査計AMSU-Aで観測された
上層暖気核の発達のタイミングとおよそ一致していることが

27

6

18

4

9

2

0

0

-9

-2

-18

-4

-27
-36

対流
バースト

対流
バースト
00

06 12
Oct. 4

18

00

06 12
Oct. 5

-6
-8
18

00

06 12
Oct. 6

18

00
(UTC)

第４図：2013年台風第24号（Danas）の発達期の上層最大接線風（青線）、上層最大動径風（赤線）、及びMTSAT赤外輝度温度（TB）の半径200km内平均（黒線）の前12時間変化。

分かっており、今後、高分解能非静力学モデルを用いた数値シミュレーションとの比較等を行なって、このメカニズムを明らかにしていきたいと考えている。
2016年台風第10号（Lionrock）

20

50
45

15

40
35

10

30

5

上層AMV最大接線風

25

ベストトラック地上最大風速

20

ドボラック地上最大風速0

地上最大風速（m/s）

8上層AMV最大接線風（m/s）

赤外TB（半径200km内平均）
上層最大接線風
上層最大動径風36

接線風・動径風の前12時間変化（m/s）

赤外TBの前12時間時間変化(K)

MTSATの赤外輝度温度の低下によって示唆される対流の強15
00 12 00 12 00 12 00 12 00 12 00 12 00 12 00 12 00 12
Aug. 22 Aug. 23 Aug. 24 Aug. 25 Aug. 26 Aug. 27 Aug. 28 Aug. 29

４．台風の発達率と上層アウトフローの強さの関係(UTC)

第５図：2016年台風第10号の、ひまわり8号AMVから求め

台風の発達率と上層アウトフローが反映する二次循環のた上層最大接線風（青プロット）、ベストトラック最大風速（黒強さの関係を定量的に明らかにするため、2011～2014年の線）、ドボラック法による最大風速（橙線）の時系列。

44個の台風を対象とした統計調査を、MTSAT AMVを使用して行なった。結果、AMVから求めた上層アウトフローの強さは、６．まとめ発達率が大きな台風ほど強い傾向があることが分かった。ま近年の静止気象衛星観測の高頻度化によって、時間・空た、発達率が大きな台風ほど、台風中心付近に背が高く発達

間的に密に得られるようになったAMVは、台風の発達においした対流雲を伴う傾向もみいだされ、上層アウトフローの強化て重要なプロセスである一次循環及び二次循環の強化の解が主に台風中心付近の対流の活発化によってもたらされて析において、有用な風データであることが示された。また、ひいたことが示唆された。

まわり8号の10分毎のAMVデータを用いることによって、従来以上の結果は、軸対称台風モデルを用いた診断的研究、の観測では捉えることが困難な、一日以下の時間スケール

数値シミュレーションを用いた研究とも整合しており、台風中の台風の強度・構造変化を捉える可能性も示唆された。今後

心付近の活発な対流活動が、台風発達に大きく寄与しているは、これまでの調査結果をもとに、台風解析業務への貢献たことを裏付ける結果である。

めの台風強度診断手法の開発等を進めていく予定である。

５．上層AMVの台風強度診断への活用

参考文献

台風の雲頂高度付近の一次循環の強さは地上付近に比（1）Oyama, R., A. Wada, and M. Sawada, 2016: Intensification of

べて弱いが、大気境界層内（地表～およそ高度1km）での一Typhoon Danas (1324) captured by MTSAT upper tropospheric

次循環の強化、二次循環による角運動量の鉛直輸送を介し

Atmospheric Motion Vectors. SOLA, 12, 135−139.

て強化される。2011～2014年の44個の台風を対象として行な（2）小山亮, 2015: MTSATラピッドスキャン観測データを用いて算出さった統計調査では、地上最大風速とMTSAT AMVから求めたれた台風領域の上層大気追跡風の特性. 天気, 62, 881−894.

上層最大接線風の間には、比較的明瞭な正相関（0.72）がみ（3）Shimoji, K., 2014: Motion tracking and cloud height assignment

られた。この結果は、AMVから求めた上層最大接線風から、

methods for Himawari-8 AMV, Proceedings of 12th International

地上最大風速の強化を診断できる可能性を示している。

Winds Workshop.

伊豆大島で進行するマグマ蓄積～火山性地殻変動は何を明らかにしたか～○山本哲也（火山研究部）
１．はじめに
気象庁では、噴火予報及び警報をはじめとして様々な火山監視のためのデータや国土地理院の電子基準点のデータも用いた。

防災情報を発表している。この業務のためには、観測データに基づいて、火山が現在どのような状態にあるのか、噴火は差し迫っているのかなどの火山活動評価を行う必要がある。
火山の噴火は地下のマグマによって引き起こされる。マグマの蓄積や移動によって火山性地殻変動が生じ、地面のわず
かな移動や傾斜、ひずみ、変形が生じる。これらを地殻変動観測でとらえることは、マグマの蓄積の様子を直接的に捉えることになるため、今日最も有力な火山監視手法のひとつとなっている。我々の研究では、火山における地殻変動のデー
タを活用することで火山活動モニタリングや火山活動評価を高度化することを目指している。

２．伊豆大島の火山性地殻変動
２．１．伊豆大島火山と地殻変動観測
伊豆大島は繰り返し噴火を引き起こしている活火山である。
過去３回の中規模噴火は36～38年の間隔で発生しており、前回の中規模噴火は1986年に発生した。噴火の多くは中央

第１図：伊豆大島における1997～2016年のGNSS観測点間の火口丘（三原山）からの噴火であり、山頂部のカルデラ域を影

距離変化（緑：山麓－山麓、青：山頂部－山麓、赤：山頂部－

響範囲とする程度の規模である。ただし、1986年噴火の際に山腹、の観測点間の距離、観測点の位置は第２図を参照）

はカルデラ外の北西山腹にかけて割れ目噴火が発生したため、全住民が島外に避難する措置がとられた。地質学的な調２．２．長期的膨張とマグマの蓄積

査によれば、過去には山腹、山麓における噴火も度々発生し第１図に伊豆大島における過去約20年におよぶGNSS観ており、また約1700年前には全島に影響をおよぼすような火測点間の距離（基線長）の変化を示した。いずれの基線長を砕流を伴う噴火が山頂部で発生したことが知られている。伊みても長期的にはほぼ一定の年率で伸びの傾向がみられ、豆大島は、現在のところ地表では目立った火山活動の兆候

顕著なものではこの間に約25ｃｍ伸びがあった。島の様々なはなく、静穏な状況であるが、前回の噴火から約30年が経過

部分が一様に伸びていることから、伊豆大島は全体として膨し、次の噴火は近づいているとみられる。

張を続けていることがわかる。

伊豆大島では火山監視や研究のために様々な観測が行第２図に、約10年間の各観測点の水平方向の動きを示し

われている。その中で気象研究所は火山性地殻変動の研究た。この間に累積した長期的膨張の様子を表している。地下のためにGNSS（Global Navigation Satellite System）、光波測に球状の圧力源が存在し、それが膨らんだためにこのような距、ひずみなどの観測を行っている。今回はGNSSによる成動きが生じたと考えると、圧力源の位置はカルデラの北部の果を中心に紹介する。GNSSは米国が運用するGPSなどに代深さ5.5ｋｍに推定される。この圧力源はマグマ溜まりに相当表される衛星からの電波を利用した測位システムである。精し、マグマの継続的蓄積によって膨張を続けていると見られ密な解析を行うことで観測点の位置を１cm以内の精度で決ている。推定される10年間で蓄積されたマグマの体積は1.9×定することができ、また、解析結果の長期安定性にも優れて

107m3に及ぶ。1986年の噴火後、この割合で蓄積がすすんでいる。解析では、我々の観測データに加えて、気象庁の火山いたとすると、過去30年間の蓄積量は約6×107m3に達し、前回の噴火で噴出したマグマの量に匹敵もしくは上回る。マグ

発生しているかは、まだ解明されていない。仮説として、マグマの量だけを考えれば次の噴火のためのマグマは既に十分

マから分離した揮発性成分（火山ガス）が地下で蓄積することに準備されているとみられる。

で膨張を引き起こし、それが地殻内に散逸することで収縮するというメカニズムが考えられているが、観測データや数値実

２．３．短期的収縮・膨張とその推移

験などによる検証はこれからである。

第１図の基線長の変化でもう一つ特徴的なのは、周期１～２年程度の縮み・伸びの変化が、各基線長変化に繰り返し見

３．まとめられることである。その縮み・伸びはすべての基線長で同じタ

現在静穏な状況の伊豆大島であるが、火山性地殻変動のイミングで発生しており、伊豆大島全体の短期的な収縮・膨観測データをみると、マグマの継続的な蓄積に伴うとみられる張の変動があることを示している。また、短期的な膨張期に島全体の長期的な膨張が続いている。膨張するマグマ溜まり

は島内の地震活動が活発化することも知られている。

はカルデラ北部の地下に推定されている。前回の噴火以降観測点を増強した2009年以降について、短期的な収縮おのマグマの蓄積量は、すでに前回の噴出量に匹敵するとみらよび膨張の変動それぞれについて圧力源を求めると、その位れる。また、火山性地殻変動として周期１～２年程度の短期

置は長期的膨張と同様にカルデラ北部に推定される。ただし、的収縮・膨張が1990年代から繰り返し起こっている。噴火に圧力源の深さは約４km程度と浅い場所に推定されるので、長むけて変動の振幅が大きくなっている可能性があることから、期的膨張とは別の圧力源である可能性が高い。1回の収縮・今後の推移に注目している。

63膨張で生じる圧力源の体積変動は3×10 m に及ぶ。
この短期的な収縮・膨張は伊豆大島北部山腹に設置された気象庁の体積ひずみ計の観測でもひずみの変動として捉えられている。第３図に、体積ひずみ計の記録から周期３年
以下の帯域の変動を抽出し、面積ひずみに換算したものを示した。同様に伊豆大島北部のGNSS観測点３点の動きから求められた面積ひずみの変動も示した。両者で短期的な変動が発生するタイミングや、その振幅の経年的な変動がよく一致しており、同じ現象を捉えていることが分かる。GNSS観測の開始以前に遡って体積ひずみ計の記録を見ると1990年代前半にはすでに小振幅の短期的変動が発生しており、振幅は2000年代にかけて次第に大きくなっている。また、現用測器ではなく更新前の体積ひずみ計による観測ではあるが、1986年噴火前にも同様の短期的変動があり、しかも振幅がより大きかったことが推定されている。短期的変動の様相は、次の噴火にむけて今後も推移していく可能性があるため、注

第２図：長期的膨張に対応する各観測点の動き（2001年３月

視していく必要がある。

～2011年２月）．赤矢印：観測値、青矢印：球状圧力源による

なお、この短期的な収縮・膨張がどのようなメカニズムで推定値、丸印：圧力源の推定位置（深さは5.5km）

第３図：（左）GNSSおよび体積ひずみ計から求めた面積ひずみの短期的変動、（右）体積ひずみ計および算出に用いたGNSS観測点の位置（緑丸）熊本地震 ～地震発生状況の即時把握技術の開発～○溜渕功史（地震津波研究部）

１．はじめに4月14日21時26分（マグニチュード（M）6.5）、4月16日01時に得られることは、地震活動の特徴を知るうえで重要な手がかりとなる。

25分（M7.3）の地震を含む平成28年（2016年）熊本地震の一連の地震活動では、熊本県で最大震度7を2回観測するなど、大きな被害が生じた。また、この一連の地震活動による震度1以上の揺れを伴う地震回数は半年間で4,000回を超え、全体
としては減衰しつつも、なお活発な地震活動が続いている。
本報告では、熊本地震の一連の地震発生状況と、それを即時に把握するための技術開発について紹介する。

２．熊本地震の地震発生状況
２．１．震源の即時推定
大規模な地震がひとたび発生すると、それに伴い数多くの余震が発生する。通常、信頼性の高い震源の推定は解析者が手動で行っているが、非常に多数の地震が発生している状
況下で、迅速に全ての震源を推定することは難しい。そこで、
溜渕・他（2016） 1)は、地震多発時においても適切に震源を決第1図：熊本地震の震央分布図。赤色は解析者が決定した定する自動処理手法を開発し、気象庁システムへの導入を震源。黒は自動処理結果のうち、精度が良いもの。灰色はそ

進めてきた。この自動処理手法は、通常利用するP波、S波のれ以外の自動処理結果。

時刻に加えて、揺れの大きさも考慮することで、通常時であれば内陸のM1.0以上の地震をほぼ全て自動的かつリアルタ
イムに決定することが可能となった。熊本地震では、この自動処理を活用することにより、4月14日から5月末までの期間に、体に揺れを感じない地震も含めて約70,000個もの地震が発生していることが明らかになった（第1図）。これは、解析者が決定した震源（当時約5,000個）に比べて、約14倍もの数に相当する。
震源がリアルタイムに得られることにより、震源域の広がり

第2図：第1図矩形領域内の地震活動経過図。縦軸はマグや地震活動の推移を即時に把握することが可能となる。震源

ニチュード、横軸は時間（日）を示す。プロットの色付けは第１域の広がりに着目すると、4月14日（M6.5）の地震は日奈久断

図中のプロットと同様。

層帯に沿った地震活動が見られた一方で、4月16日（M7.3）の地震以降の震源分布はさらに北東の布田川断層帯や、大分

２．２．自動震源データを用いた地震活動の統計量推

県の別府－万年山断層帯に沿って広がっていることが分か定る。また、決定されたマグニチュードの下限に着目すると、人

地震活動を多数の地震の集まりとして見た場合、その特徴

手で精査した震源はM2より大きな規模の地震が概ね決定さを表す様々な統計量を求めることができる。地震の発生頻度れているのに対して、自動処理による震源はM0やM1程度のと規模にはGutenberg-Richterの法則（G-R則）と呼ばれる経規模の地震についても多数決定できていることが分かる（第2験則があり、地震の規模を表すマグニチュードが大きくなるに図）。規模の小さい地震も含めて、多数の震源がリアルタイム

つれて、その規模の地震の頻度は指数関数的に減少することが知られている。この減少の程度をb値という指標で表す。b

開発を進める。

値が低いと、相対的に規模の大きな地震の発生頻度が多い
ことに対応する。余震確率の計算では、このb値が重要なパラメータとなっている。また、2011年東北地方太平洋沖地震に先行してb値が低下していたことが報告されており2)、規模が大きな地震の発生を予測するうえで重要な指標と考えられる。
自動処理による震源は、人手を介しての処理に比べて圧倒
的に迅速かつ数多く決定できることから、このような統計量の算出にとって、有用なデータとなる。
そこで、自動震源を使用して、b値の推定を行った（第3図）。
通常、b値を求める際には、解析対象のM下限を検知力以上の値に設定して解析するが、リアルタイムかつ客観的に算出することを目的として、G-R則のフィッティング具合から自動的にM下限を設定する手法を用いた3)。
4月14日（M6.5）の地震の発生前の10年間は、b値は概ね0.8程度の値であったが、4月14日の地震後に0.5～0.6程度に

第３図：b値の空間分布低下した。4月16日にはM7.3の地震が発生し、布田川断層帯の北部ではb値は回復した。一方で、その後の地震活動は、

謝辞

ほとんどの領域で地震発生前のb値の状態に回復しているも北海道大学、弘前大学、東北大学、東京大学、名古屋大学、のの、日奈久断層帯の南部では引き続きb値の低い状態が京都大学、高知大学、九州大学、鹿児島大学、国立研究開

継続していることが明らかになった。このb値の空間分布と合発法人防災科学技術研究所、国立研究開発法人産業技術総成開口レーダ（SAR）干渉解析による余効変動（大きな地震の合研究所、国土地理院、青森県、東京都、静岡県、神奈川

後に長期間に渡ってゆっくりと進行する地殻変動）の分布を

県温泉地学研究所、横浜市、国立研究開発法人海洋研究開比較すると、余効変動が大きな領域ではb値は回復しており、発機構及び気象庁の地震波形と、その波形を気象庁と文部日奈久断層帯南部の余効変動が小さな領域では引き続きb

科学省が協力してデータ処理した一元化震源を使用しま

値が低いといった対応関係が見られた。

した。

このような高い分解能でb値の空間分布や時間変化を解析するためには、数多くの震源が必要なため、これまで解析を参考文献行うことが難しかったが、本研究による自動処理技術の高度（1）溜渕功史・森脇健・上野寛・束田進也 2016: ベイズ推定化により可能となりつつある。

を用いた一元化震源のための自動震源推定手法, 験震時報, 79, 1-13.

３．まとめ(2)Nanjo, K. Z., N. Hirata, K. Obara and K. Kasahara 2012:

自動処理手法の開発により、従来よりも迅速かつ大量に震Decade-scale decrease in b value prior to the M9-class

源を推定することが可能となった。これを熊本地震に対して適2011 Tohoku and 2004 Sumatra quakes, Geophys. Res.

用したところ、従来の約14倍にあたる数の震源を自動的に得

Lett., 39, L20304, doi:10.1029/2012GL052997.

ることができた。また、この自動処理によって得られた震源を（3）Wiemer, S. and M. Wyss 2002: Mapping spatial variability

用いて地震活動の統計量の一つであるb値を算出すると、日of the frequency-magnitude distribution of earthquakes, in

奈久断層帯北部では4月14日（M6.5）の地震に伴いb値が低“Advances in Geophysics”, ed. by R. Dmowska and B.

下し、4月16日（M7.3）の地震後に回復するといった現象が確

Saltzman, 45, Academic Press, 259-302.

認できた。今後、このような地震活動の統計量を迅速かつ安定して推定することにより、地震活動の推移把握の高度化を目指す。さらに、内陸の地震のみならず、海域の地震活動についての自動決定手法を向上させることにより、南海トラフな
どの海溝沿いの巨大地震の地震像を即時に把握する技術の気象研究所で行っている研究について
気象研究所では、気象庁が発表する気象や地震・火山、気候等に関する各種情報の高度化に向けた実用的な研究を重点的に実施しています。また、将来に向けた我が国の気象業務の健全な発達を図るため、長期的に各種情報の高度化等を見据えた基礎的・基盤的な研究や、世界をリードする研究を実施しています。

１

重点研究

各種防災情報や地球環境情報の高度化など気象業務への貢献を目指し、5 年以内に業務化のめどをつける実用的な研究課題で、台風・集中豪雨等対策、地震・火山・津波対策及び気候変動・地球環境対策の強化に資する研究を重点的に実施しています。

１．１

台風・集中豪雨等対策の強化に関する研究

気象災害を防止・軽減するには、災害をもたらす現象の観測・解析技術及び予測技術を高度化し、予報・警報等の防災情報を避難等防災活動の早期準備や迅速・的確な実施に一層活用可能なものにしていく必要があります。このため以下の研究に取り組んでいます。
A1 メソスケール気象予測の改善と防災気象情報の高度化に関する研究
A2 顕著現象監視予測技術の高度化に関する研究
A3 台風の進路予報・強度解析の精度向上に資する研究
A4 沿岸海況予測技術の高度化に関する研究

１．２

地震・火山・津波対策の強化に関する研究

地震、火山及び津波に関する防災情報をより的確なものとし、それらによる災害を防止・軽減するためには、発生した現象の推移をより的確に観測・解析する技術を開発するとともに、地震の予知技術、火山
噴火及び津波の予測技術の高精度化を進める必要があります。このため、以下の研究に取り組んでいます。
B1 緊急地震速報の予測手法の高度化に関する研究
B2 地震活動・地殻変動監視の高度化に関する研究
B3 津波の予測手法の高度化に関する研究
B4 大規模噴火時の火山現象の即時把握及び予測技術の高度化に関する研究
B5 地殻変動観測による火山活動評価・予測の高度化に関する研究
B6 南海トラフ沿いのプレート間固着状態監視と津波地震の発生状況即時把握に関する研究
B7 火山ガス観測による火山活動監視・予測に関する研究

１．３

気候変動・地球環境対策の強化に関する研究

地球温暖化やオゾン層破壊、国境を越えた大気汚染など地球規模の環境問題が顕在化する中、環境と経済を両立し持続可能な社会の発展を実現するため、気候及び地球環境に関する信頼性の高い情報が求められています。このため、以下の研究に取り組んでいます。
C1 気候モデルの高度化と気候・環境の長期変動に関する研究
C2 季節予報の高度化と異常気象の要因解明に関する研究
C3 地球環境監視・診断・予測技術高度化に関する研究

２ 一般研究気象庁が発表する各種情報の高度化等気象業務への将来の実用化を見据えた基礎的・基盤的な研究や、世界をリードする革新的テーマに挑戦する研究を推進しています。

３ 地方共同研究
気象業務の現場において取り組むべき研究課題について、気象研究所と気象官署が共同し地方共同研究として実施しています。地方共同研究により、気象業務の現場における潜在的なニーズを的確に捉え、気象研究所の研究方針や内容に適宜反映させることによって、気象業務の高度化に貢献しています。さらに、研究活動を通じて気象研究所と気象官署の連携を強化し、気象官署における調査業務の支援を図るとともに、職員の資質向上にも貢献しています。

４ 機動的な研究
竜巻、集中豪雨、地震、火山噴火等の大きな災害を伴う顕著現象が発生した場合は、気象庁本庁・管区
気象台・地方気象台等と連携し、現地調査を含む調査研究等を機動的に実施しています。

気象研究所成果報告会 発表題目（平成18年度以降）
H27
H27
H27
H27
H27
H26
H26
H26
H26
H26
H26
H25
H25
H25
H25
H25
H24
H24
H24
H24
H24
H23
H23
H23
H23
H23

気象研究所での気候変動研究概要
気候変動予測研究の過去・現在・未来
温暖化予測情報をどう使うか
エルニーニョなどの海洋の変化を予測するために成層圏の変動の気候への影響について
御嶽山噴火 －水蒸気噴火のメカニズム解明に向けて－
南海トラフ沿い巨大地震は予知できるか？ －地震発生シミュレーションからの知見－
竜巻等突風の探知 －フェーズドアレイレーダーが切り拓く世界－
集中豪雨の発生メカニズム解明に向けて －平成 26 年 8 月 20 日広島豪雨事例－
温暖化に伴い強雨は増えるのか？ －アメダス観測が示す気温と強雨の関係－
CONTRAIL プロジェクト －大型旅客機による地球規模の温室効果ガス観測－
GPS が拓く突風・大雨の予測・監視への道
局地的豪雨等の予測精度向上に向けた取り組み
津波予測の新展開 －迫り来る津波を捉え予測する－
地球温暖化とともに変わりゆく日本の気候を考える
海洋の CO2 増加と酸性化 －西部北太平洋熱帯・亜熱帯域における実態－
地震の揺れのリアルタイム予測 ： 次世代の緊急地震速報を目指して
ヒートアイランドの現在と未来
グリーンランドから地球温暖化をみる
つくば竜巻 ： 二重偏波レーダーによる実態解明
平成 23 年（2011 年）東北地方太平洋沖地震の地震像と巨大地震の早期規模推定に向けた取り組み
平成 23 年（2011 年）東北地方太平洋沖地震に伴う津波と津波即時予測の研究
2011 年霧島山（新燃岳）噴火対応～噴煙エコーの解析と降灰の量的予測に関する研究～
集中豪雨の統計的・解析的研究 －2010 年 10 月 20 日の奄美豪雨－
福島第一原子力発電所事故に伴う大気・海洋の人工放射能の変動

蒲地 政文
尾瀬智昭
高藪 出
山中 吾郎
黒田 友二
横山 博文
弘瀬 冬樹
楠 研一
加藤 輝之
藤部 文昭
松枝 秀和
小司 禎教
瀬古 弘
林 豊
佐々木 秀孝
石井 雅男
干場充之
青柳曉典
青木輝夫
山内 洋
加藤輝之
吉田康宏
対馬弘晃
新堀敏基
津口裕茂
青山道夫

研究総務官
気候研究部
環境・応用気象研究部
海洋・地球化学研究部
気候研究部
火山研究部
地震津波研究部
気象衛星･観測ｼｽﾃﾑ研究部
予報研究部
環境･応用気象研究部
海洋・地球化学研究部
気象衛星・観測ｼｽﾃﾑ研究部
予報研究部
地震火山研究部
環境・応用気象研究部
海洋・地球化学研究部
地震火山研究部
環境・応用気象研究部
物理気象研究部
気象衛星・観測ｼｽﾃﾑ研究部
予報研究部
地震火山研究部
地震火山研究部
地震火山研究部
予報研究部
地球化学研究部

H22
H22
H22
H22
H22

台風に伴う竜巻等突風の発生機構―2006 年台風第 13 号に伴う延岡竜巻の数値シミュレーション―
気候研究のための長期日降水グリッドデータの作成
数値モデルを用いた地球環境解析
衛星データの利用技術に関する研究【記念講演】気象観測鉄塔

益子 渉
上口賢治
眞木貴史
増田一彦
石原正仁

予報研究部
気候研究部
環境・応用気象研究部
気象衛星・観測ｼｽﾃﾑ研究部
気象衛星・観測ｼｽﾃﾑ研究部

H21
H21
H21
H21
H21
H21

東南アジア地域の気象災害軽減国際共同研究
2008 年台風特別観測（T-PARC）結果
伊勢湾台風再現実験プロジェクト（re-analysis/forecast of Typhoon Vera(1959) project: ReVera）
相似地震の発生予測実験
地震波で地下の時間変動を監視する
温暖化による日本付近の詳細な気候変化予測

齊藤和雄
中澤哲夫
別所康太郎
岡田正実(客員)吉田康宏
行本誠史 他

予報研究部
台風研究部
台風研究部
地震火山研究部
地震火山研究部
気候研究部ほか

H20
H20
H20
H20
H20
H20
H20
H20

海洋における炭素循環の変動に関する観測的研究Ⅱ
物質循環モデルの開発改良と地球環境への影響評価に関する研究
接地境界層における乱流の構造に関する研究
地震・地殻変動観測データの高度利用に関する研究
津波の予測精度向上に関する研究
火山観測データの気象補正等による高精度化に関する研究-ノイズの中から火山活動を抽出する新地殻変動観測手法の研究－スロースリップ検知を目指したレーザー式変位計の開発－
３次元数値モデルによる地震発生シミュレーションの研究

緑川 貴
柴田清孝
毛利英明
前田憲二 他
林 豊
山里 平
勝間田明男
弘瀬冬樹

地球化学研究部
環境・応用気象研究部
物理気象研究部
地震火山研究部
地震火山研究部
地震火山研究部
地震火山研究部
地震火山研究部

H19
H19
H19
H19
H19
H19
H19

氷晶発生過程に関する研究
上陸台風の構造変化過程とそれに伴う暴風、豪雨、高潮の発生に関する研究
日本の異常気象の実態及び気候変動との関連に関する研究
気候システムとその変動特性のモデルによる研究
高解像度(渦解像)海洋大循環モデルの開発とそれによる水塊の形成、維持、及び変動機構の解明
海洋データ同化システムの高精度化と海洋現象の季節から経年変動の解析
アジア大陸の影響による大気微量気体・エーロゾル・降水降下塵の化学組成変動に関する研究

村上正隆
中澤哲夫
藤部文昭 他
楠 昌司
本井達夫
蒲地政文
松枝秀和

物理気象研究部
台風研究部
予報研究部ほか
気候研究部
海洋研究部
海洋研究部
地球化学研究部

H18
H18
H18
H18
H18
H18
H18
H18
H18
H18
H18

季節予測システムの構築と経年変動機構・予測可能性の研究
安田珠幾
エーロゾルの地表面放射への影響評価のためのエーロゾル光学特性の推定について
内山明博
大気エーロゾルが積雪面の放射特性に与える影響
青木輝夫
海洋における炭素循環の変動に関する観測的研究
緑川 貴【特別報告】 太平洋をＵターンした津波 －千島列島付近の地震による津波の数値解析結果－
長谷川洋平
ヒートアイランド現象の再現・予測に関する基礎的研究
三上正男
都市ビル群上での乱流特性の研究
栗田 進
衛星搭載多波長赤外サウンダの利用技術に関する研究
真野裕三【特別報告】2006 年 11 月 7 日の佐呂間町竜巻の現地調査報告
鈴木 修【特別報告】佐呂間町に竜巻をもたらしたスーパーセルの再現実験-水平分解能250m の雲解像モデルの結果から- 加藤輝之
非静力学モデル（NHM）の高度化と同化技術の改善に関する研究中村誠臣

社会的に大きな影響があった顕著現象(平成 24 年)～急発達した低気圧・つくば竜巻・九州北部豪雨～

気候研究部
気候研究部
物理気象研究部
地球化学研究部
地震火山研究部
環境・応用気象研究部
環境・応用気象研究部
気象衛星・観測ｼｽﾃﾑ研究部
気象衛星・観測ｼｽﾃﾑ研究部
予報研究部
予報研究部

気象庁気象研究所
連 絡 先：気象庁気象研究所企画室〒305-0052 茨城県つくば市長峰１－１
TEL 029-853-8546
e-mail：ngmn11ts@mri-jma.go.jp
ホームページアドレス：http://www.mri-jma.go.jp