資料２－５次世代火山研究・人材育成総合プロジェクト
総合協議会（第７回）
Ｒ２.１２.１５

課題C 火山噴火の予測技術の開発
第7回研究運営委員会

課題成果概要2019年12月～2020年度上期までの成果概要

課題責任者・説明者：北海道大学 中川光弘

課題Cの概要

課題Cでは国内の主要な活火山を対象に噴火履歴の解明と噴火事象の解析を行い、得られた情報を数値シミュレーションで解析することによって噴火の予測技術を開発する。まず個々の火山で中長期予測を行う。そして事象分岐判断基準が伴った噴火事象系統樹を整備するとともに、噴火発生確率の算出に向けた検討を行う。
本課題は、以下の3つのサブテーマの研究が並行して、かつ密に連携しながら実施される。
火山監視・火山防災
噴火対応

課題C: 火山噴火の予測技術の開発

噴火予測・噴火ハザード予測
手法の開発

C-3

噴火確率算定への寄与

噴火事象系統樹の高度化

マグマ移動シミュレーション
噴火ハザードシミュレーション
マグマ過程の再現実験
事象のモデル化⇒パラメータ解析

C-1

C-2

最重点5火山・重点火山
地質学的・物質科学的
地球化学的解析

課題D-3
情報ツール課題D-2
リアルタイム降灰

課題D-1
リアルタイム
無人機観測
観測情報の共有

噴火事象分岐条件の解明

噴火事象系統樹の作成

データ・成果の一元化

データ連携

噴火事象分岐の判断基準の確立

最近の噴火事象物質科学的解析

課題A
一元化DB

課題B
観測データ

中長期噴火予測

長期噴火履歴・噴火推移の解明
マグマ変遷の解明

大学院生受入れ火山研究
人材育成
コンソーシアム

本日の報告〇各サブテーマの成果概要および成果例
サブテーマ２：噴火履歴調査による火山噴火の中長期予測と噴火推移調査に基づく噴火事象系統樹の作成
サブテーマ１：火山噴出物分析による噴火事象分岐予測手法の開発
サブテーマ３：シミュレーションによる噴火ハザード予測手法の開発〇課題C連携研究：玄武岩質噴火（伊豆大島）

東京大学地震研究時HPより引用

サブテーマC-２

噴火履歴調査による火山噴火の中長期予測と噴火推移調査に基づく噴火事象系統樹の作成

1年目

2年目

3年目

4年目

5年目

6年目

7年目

8年目

9年目

10年目

H28

H29

H30

R01

R02

R03

R04

R05R06R07

成果の総括と補填調査検討

最重点火山見直し

地表踏査、トレンチ・ボーリング調査（大学連合・産総研）
物質科学的解析
火山の噴火履歴およびマグマ長期変遷の解明

マグマ変遷解析センター（北大）
立上げ・運用開始

維持・整備

個々の火山の時間―噴出量階段図の作成
地質学的データの一元化データベースへの提供
大規模噴火DB完成

大規模噴火推移DBの構築（産総研）完成

噴火推移DB(仮)※の構築（産総研）※ 採録対象・仕様についてはH32年度末までに決定

手法の確立

中長期噴火予測 手法の開発
中長期噴火
予測・噴火
事象系統樹の作成

成果とりまとめと評価

個々の火山の中長期噴火予測

成果とりまとめと評価
手法の確立
噴火事象系統樹 作成手法の検討 個々の火山の系統樹試作

課題C内・他課題との連携公表

2020年度の活動・成果概要

サブテーマC2〇火山の噴火履歴およびマグマ長期変遷の解明• 噴火活動履歴解明のためのボーリング掘削調査およびトレンチ掘削調査を実施
ex. 有珠山・アトサヌプリ・阿蘇山・雌阿寒岳・秋田焼山・鳥海山・蔵王山・浅間山など• 野外調査に基づく，噴火履歴・噴火活動推移履歴の解明および物質科学的解析による長期マグマ変遷の解明
課題C連携研究ex. 摩周・十勝岳・秋田駒ケ岳・三宅島・霧島山・九重山・白山・伊豆大島・新島など• 活動的火山としての認知度が低い火山について，噴火履歴・噴火推移履歴を解明
ex. ニセコ・鷲羽池火山• 現在起きている噴火活動への対応
ex. 阿蘇山〇大規模噴火データベースの整備（産総研）• 噴火推移データの集約、噴火推移閲覧サイトの試作版構築〇マグマ変遷解析センター（北大）• マグマ変遷解析センターの整備および分析ルーチンの確立• 参加・協力機関のユーザー受入れ〇ボーリングコアの管理システム構築（北大・防災科研）• 管理システム構築のための既存コア試料の集約開始〇その他• 噴火事象系統樹の試作をテーマとした研究集会開催予定• 人材育成コンソーシアムへの貢献（講師としての参画・大学院生の研究テーマ実施）• 地域社会への成果普及活動（トレンチ掘削調査時の説明会など）

鷲羽池火山：活動的火山の可能性の検討（富山大）

サブテーマC2

中部日本に位置する鷲羽火山の完新世の活動の可能性について検証するために，野外調査を実施した．
【成果】
・広域テフラK-Ah（7300年前）の上位に、鷲羽池火口が給源と推測される水蒸気噴火のテフラ、火砕サージ、噴石を確認（年代測定を予定）・同火口南の硫黄沢ではR2年に少なくとも3回大規模噴気活動が発生（三俣山荘での聞取り調査）

北西から見た硫黄沢の大規模噴気活動

鷲羽池火口が給源と推測されるテフラの等層厚線と、（出典：三俣山荘インスタグラム URL
https://www.instagram.com/explore/tags/三俣山荘） 火砕サージ・噴石の分布範囲（火口東域は未調査）（基図：国土地理院 https://maps.gsi.go.jp/）

大規模噴火データベースの構築（産総研）

サブテーマC2

R01.12～R02.11
プリニー式噴火の噴火推移と大規模噴火の前駆活動と噴火推移の特徴を取りまとめ，噴火推移閲
覧サイトの試作版を構築した．
大規模噴火の前駆活動と噴火推移

プリニー式噴火推移
浅間1783年噴火

7日間90日間

3年間

200年間

有珠1977年噴火

0.5日間

噴火推移閲覧サイトの試作版
3日間

90日間

200年間

Calbuco 2015，Kelud 2014, Puyehue
Cordon-Caulle 2011, Grimsvoton 2011,Chaiten 2008, Pinatubo 1991,
Augustine 1976, 桜島 1914,
Novaraputa-Katmai 1912, Krakatau
1883, Askja 1875, 浅間1783, 富士
1707, Vesuvius 79 (14噴火)

火山灰層7 ＊

縦軸は，VUC (Volcanic Unrest Classification：火山活動強度指数)【成果】 主要プリニー式噴火の噴火推移をとりまとめ，
閲 覧 サ イ ト の 試 作 版 を 作 成 ． 7 月 に GSJ 研 究 資 料 集(no.700)を公開．
プリニー式噴火による降下テフラのGISデータを整備．

【成果】 過去13万年間に発生した46の大規模噴火(VEI>6)の前駆活動の特徴を解析．6月に大規模噴火の前駆活動と噴火推移を取りまとめたGSJ研究資料集(no. 699)を公開．大規模火山噴火のカルデラと降下テフラのGISデータを整備．

サブテーマC-１
火山噴出物分析による噴火事象分岐予測手法の開発
項目

分析・解析プラット
ホームの構築内容(a) 高精度，高効率の分析環境の構築

28年度

29年度

30年度

R1年度

(c) データ保存環境の整
備

機器の公開

(e) 広く研究者や学生に開放するための利用環
境の整備

R5年度

R6年度

R7年度

分析環境のアップデート

データベースの構築
データベース利用環境の整備
データベース利用開始，アップデート
利用規程整備受け入れ体制の構築
受け入れ開始必要に応じての追加

マグマ溜まり環境の実体化

火山噴出物解析

R4年度

標準的な解析プログラム，図化プログラムの整備
データベースとの連携プログラムの整備プログラムのアップデート

試料採取，地質学的検討(d) 対象火山（1１火山＋α）の噴出物解析とカタログ化

R3年度

標準的な定量分析ルーチンの確立
分析の効率化(b) 解析の自動化による作業の効率化・標準化

データベース

R2年度再解析，精密化

マグマの上昇速度，状態の解明 再解析，精密化
混合から噴火に至る時間の解明
カタログ化

再解析，精密化
カタログのアップデート

予測指標の検討（新たな分析要素の検討）
今後の展開
R2年度で分析・解析プラットホームの基礎的なハードウエアの導入と分析ルーチン立ち上げ「分析解析プラットホームの構築」 完了．今後はより使いやすいようにソフトウエアやユーザーインターフェイスの整備に注力．
「データベース」 収集したデータの公開にむけて準備．
「機器の公開」 利用規定や操作マニュアルは整備済み．来年度の利用受け入れに向けて準備中．
「火山噴出物解析」 引き続き11火山の噴出物について解析を続行．個々の火山について，マグマ供給系像の精密化をはかる．

2020年度の活動・成果概要

サブテーマC1(a)火山噴出物の分析・解析
・COVID-19による活動制限の影響で，野外調査や分析に非常に大きな影響が出ている．
例）地震研への分析来所数(4-10月)： 2019年は94人・日，2020年は 6人・日に激減．

・個々の参加組織でこれまでに採取した試料を利用して，マグマ供給系の特徴やマグマの上昇速度の把握に努めている（成果例１・補遺）(b)分析・解析プラットホームの構築• サブ課題外からの利用受け入れにむけて，装置利用手引き等を整備．
• COVID-19対応するため，リモート分析を試行→ 依頼者が正確に分析点情報を分析作業者に伝達できるような仕組みを構築
し，試験的に運用中（補遺）• 熱力学解析ツール（MELTS）の運転環境の整備と結果のデータベース化（成果例２）(c)データサーバーとデータベースの整備• 収集画像データ（COMPO像，X線マップ）をWeb経由で試験的に公開中＞＞

画像の所有者（＋指導教員）とオペレーターのみがサーバー上で画像閲覧可能
画像所有者以外も，どのような画像が存在するのかについての情報は閲覧可能

成果例１．石基組織から探る富士山のマグマ供給系
これまでに構築した分析・解析プラットホームの石基画像解析機能を使った成果

サブテーマC1

・この間，全岩化学組成自体は大きくは変化せず

青：斜長石
水色：ガラスSiO2・浅所のわずかな変化（崩壊体積1-2km3 << 山体体積（数百km3))であるにもかかわらず，噴出物の様態が急変し，その影響が長期間に及ぶ．

S-14-3 石基微斑晶が少ない
青：斜長石
水色：ガラス

斜長石量が急減
ガラス量が急増
S-14-1 石基微斑晶が多い

山体崩壊が発生してから数百年にわたって，小規模ながら，石基微斑晶が乏しく，噴出物がガラス質の爆発的噴火が続く．

2900年前，御殿場山体崩壊

富士山での密度成層を考慮したマグマ供給系モデルの再構築

成果例２．MELTS Databaseの構築

サブテーマC1MELTS Databaseとは，課題C-1が対象とする11火山を含む主要な火山について，あらかじめ最近の噴火の噴出物について，様々な温度・圧力条件で結晶分化によるマグマの変化を熱力学ルーツのMELTSを用いて計算しておき，その結果にいつでも簡単に参照できるようにしたもの．
これまでに計算した
組成範囲＜2020年西之島の火山活動＞
石基ガラス
組成
データベースの活用例

液組成変化図

全岩組成

相平衡図
酸素雰囲気：QMF
含水量： 0.5 wt%以下
温度：1075-1090℃
1075 1080

データベースの内容例

1090℃

全岩と石基ガラス組成が判明した段階で，結晶分化でガラス組成を生じさせる条件の絞り込みを行なった．

サブテーマC-３
H28年度

シミュレーションによる噴火ハザード予測手法の開発
H29年度

H30年度

R1年度

R2年度

R3年度

R4年度

R5年度

R6年度

R7年度

①地下におけるマグマ移動シミュ
レーション

火道流・岩脈・結晶化・レオロジーモデル構築

マグマ移動評価システム開発

一元化システムへの統合②噴火ハザードシミュレーションの開発・高度化

降灰・噴煙・溶岩流・噴石評価システム開発

火山ハザード評価システム開発

一元化システムへの統合

2020年度の活動・成果概要（１）地下におけるマグマ移動シミュレーション• 噴火機構シミュレーション技術開発：低結晶度・高結晶度マグマ物性による火道流モデル
シミュレーションによる噴火様式の再現（1986年伊豆大島噴火）• マグマ移動過程シミュレーション技術開発：境界要素法で地形形状を考慮した地殻変動計
算手法の開発 (例：伊豆大島）• 物性モデルの構築：静置による応力オーバーシュートの測定（例：1986年伊豆大島噴火溶岩）（２）噴火ハザードシミュレーションの開発・高度化• 降灰ハザード予測モデルの開発：JMA-ATM開発（降灰および浮遊火山灰のシミュレーション機能）、降下火砕物の形状・密度の粒径依存性調査、移流拡散系の数理構造を理解するための研究• 噴煙柱ダイナミクスモデルの開発：噴煙ダイナミクスと火山灰輸送をカップリングした
SK-3Dモデルを用いた大規模シミュレーション実施• ハザード評価システムの開発：ハザード評価システムの運用開始・グランドデザイン実施、高濃度粒子流のオープンソース2次元数値コードを用いた溶岩ドーム崩壊型火砕流モデ
ルの開発

②噴火ハザードシミュレーションの開発・高度化
a.降灰ハザード予測モデルの開発

サブテーマC3

降灰ハザード予測モデル開発【令和２年度の目標（新しい移流拡散モデルの開発）】
前年度までに行った「b.噴煙柱ダイナミクスモデル」による噴煙放出を初期条件とする新しい移流拡散モデル（JMA-ATM）の開発を進める．
【令和２年度の進捗（新しい移流拡散モデルの開発）】
平成28～30年度に降灰シミュレーションの予備計算で使用したプロトタイプ（JMA-RATM）を更新するため，令和元年度から新しい移流拡散モデル（JMA-ATM）の開発を開始し，昨年度はモデル地形なしのドライモデルによる降灰シミュレーションを試行した．今年度は，モデル地形ありのフルモデルの開発を継続している．JMA-ATMは，従来のJMA-RATMと大気中の浮遊火山灰を主な予測対象としたJMA-GATMを統一しており，降灰だけでなく浮遊火山灰のシミュレーションにも資する設計にしている．一例として，開発期間中に発生した西之島の連続噴火について，浮遊火山灰と降灰シミュレーションの動作確認を行った結果を図に示す．

火山灰雲の予報官解析

予報官による気象衛星ひまわり８号の画像解析（解析時刻：2020年8月2日08時20分）
西之島の航空路火山灰情報2020/191号．8月3日00時発表より開発中のJMA-ATM（フルモデル）による気柱積算（左図）と降灰シミュレーション（右図）の動作確認結果．初期条件はSuzuki1983による，西之島の連続噴火（噴煙高度：
2020年8月1日09時～2日03時19000 ft, 2日03時～5日03時16000 ft, 5日03時～6日21時
13000 ft）を仮定．大気場は気象庁全球予報値（初期時刻：1日00 UTC，予報時間：5.5日）を使用．動画は6時間ごと．

サブテーマC3c. ハザード評価システムの開発

シミュレーションから統合的ハザード情報の提供

地形データ
計算実行
計算条件等
（スパコン,アップロード docker)

可視化(MF)
脆弱性DB統合GIS(kml)出力→JVDN

シミュレーションDB統計処理（確率計算）（AI機能）

入力ファイル形式の汎用化、異なるシミュレーションのアウトプットの汎用化（メッシュ統合等）→ シミュレーション結果のデータベース化→ 複数のシミュレーションによる確率計算→ 対象メッシュにおけるマルチハザード評価

高濃度粒子流のオープンソース2次元数値コードを用いた雲仙1991年6月3日火砕流シミュレーション
数値計算における火砕流高濃度部の厚さ

火砕流ハザードモデル開発

数値計算における最終堆積物の分布

堆積物分布について，
本数値計算結果が観測データを概ね説明できることを確認

1991年6月3日火砕流堆積物分布•(Modified from Fig. 8 of Yamamoto et al., 1993)

数値コード:
faSavageHutterFOAM (Rauteret al., 2018)

課題C連携研究：玄武岩質噴火（伊豆大島）

課題C連携研究

玄武岩質マグマの多様な噴火様式やその遷移を支配する要因について、火道流モデルで検討し、噴火予測・推移予測について課題全体で取り組む
対象：
安永噴火（Y1）1986年噴火・15世紀（Y4）

噴火様式の変遷

上昇中のマグマの結晶化

火道流モデル（C-3）

石基組織の解析（C-1）
マグマの初期条件・物性（Costa et al., 2007）

マグマ供給系（C-2）

火山地質図（川辺，1998）

圧力・温度・含水量（C-1）

噴火活動履歴（C-2）

R02年度計画：サブ課題ごとで対象噴火イベントについて検討し、研究集会で議論
進捗状況：コロナ禍で計画通りには進んでいないできる範囲で個々の検討は進めており、今後連携し発展させていく予定

課題C連携研究

噴火活動履歴解明のための沿岸域調査および
R01.12～R02.11
歴史時代噴火の活動推移の多様性の解明のための地質学的・物質科学的検討
沿岸域調査 産総研

全取得データ

歴史時代大規模噴火の活動推移

東大地震研歴史時代大規模噴火Y1, Y2, Y4, Y6について，層序を再構築し，各噴火の層序毎の噴出物の特徴抽出と分析を進めている。また，1986年噴出物と比較し，大規模噴火と中小規模噴火の違いを生み出す原因についても検討を進めている。
図1

図2

全岩化学組成m

北西沿岸部の海底地形(1mメッシュ)

1986

Y1C

伊豆大島：沿岸海底部の活動度調査

・沿岸部水深400m程度までの範囲で精密海底地形調査を実施（13日間漁船傭船）【成果】これまで知られていない海底部における側火山体や溶岩流の分布を解明

図3・Y1噴火：化学組成や岩石組織の時間変化あり（図１）
・噴火強度が大きいY1C噴出物は，より高密度かつ低マイクロライト量で，1986噴出物などと特徴が大きく異なる(図２・３）

課題C連携研究

伊豆大島1986年噴火や過去の歴史噴火を対象とした多様な噴火様式の火道流数値モデルシミュレーション
1986年伊豆大島噴火“Aマグマ”: 低結晶度 -> 溶岩流“Bマグマ”: 高結晶度 -> サブプリニー式(石橋・種田, 2018)
結晶度

A、Bマグマの条件において、非爆発的、爆発的噴火に対応する火道流を再現

メルト粘性

RhyoliteMELTSによる計算

各マグマの物性を火道流モデルに適用

噴火タイプの存在条件を推定

非爆発的爆発的今後の計画：
安永噴火への適用、時間発展変動解析、地殻変動計算

補遺：
各サブ課題の成果例

サブテーマC2

秋田駒ケ岳：野外調査に基づく
過去1.3万年間の噴火活動履歴の解明（北大，エンバイオ・エンジニアリング）
南部カルデラ２

南部カルデラ１

AK12

南部カルデラ２形成・国見火砕流

AK13

南部カルデラ１形成・生保内火砕流

図１．秋田駒ケ岳の1.3万年間活動の階段図（和知ほか（1997）を改変）

東北地方で最も活動的な火山の一つである秋田駒ケ岳について、約1.3万年間の噴火履歴・噴火様式の再検討を行っている（図１）。本年度では、最も大規模なカルデラ形成噴火について以下のことが明らかになった。
１．時間間隙を挟んで、２回の爆発的噴火により複合カルデラが形成された（図２・３）。
２．それぞれの噴火では火砕流が発生した（図３・４）。
３．それぞれの噴火では異なったマグマ供給系が活動した（図３）。

南岳火砕丘

生保内火砕流図２．秋田駒ケ岳、南部カルデラの地形・構造および生保内火砕流の分布

軽石噴出から縞状軽石・スコリア噴出へ

スコリア噴出図４．南部カルデラ１内部で新たに見出された、南岳火砕丘を覆う国見火砕流（新称）図３．南部カルデラの噴火推移と本質物の時間変化

サブテーマC2
阿蘇山：トレンチ掘削調査による最近数千年間の噴火史と中岳2019～2020年マグマ活動（熊本大）

中岳北東2.7 km地点でのトレンチ調査

中岳2019年7月～2020年6月マグマ噴火

定期的な降灰調査の結果，約11か月間のマグマ噴火
で放出された火山灰の総量は110万トン程度に達することがわかった．

火山灰の構成粒子

ボーリング掘削調査を追加実施予定
阿蘇火山では3600～4000年前頃の時期に断続的な灰噴火のほか，厚い降下スコリアをもたらす準プリニー式噴火，さらに水蒸気噴火などの火山活動が高頻度で発生していることが明らかとなった．

本質物と考えられる新鮮なガラス片の割合は噴火開始から2019年8月中旬頃まで増加
した後は減少傾向にあり，2020年2月以降は1割以下にまで低下した．

0.2 mmサブテーマC1
成果例．富士山のマグマ供給系の特徴：
斑レイ岩捕獲岩から探る富士山のケイ長質マグマだまり
1707年宝永噴火では、富士山で初めてデイサイト質マグマが噴出→富士山での珪長質マグマの発達を示唆• 宝永噴火の斑レイ岩捕獲岩に含まれるガラスから、ケイ長質マグマだまりの状態について調べた。
•1707宝永噴火（現時点で最新噴火）

MELT-B富士山宝永噴火絵図 ｢昼乃景気｣(静岡県立中央図書館）

宝永噴火の斑レイ岩捕獲岩

斑レイ岩の粒間ガラス（大塚他, 2020JpGU）••

斑レイ岩の粒間ガラスの化学組成を分析したところ、玄武岩質安山岩～流紋岩質組成をもつガラスを発見
Alに乏しい流紋岩質メルト（MELT-A）と、Alに富むデイサイト質メルト（MELT-B）の2つのケイ長質端成分メルトが存在し，苦鉄質メルトと混合している

Mafic melt
MELT-A

石英の発見（石橋他, 2020火山学会）••
石英を含むクリスタルマッシュ試料（FjGB18）のBSE像と元素濃度マップ

石英と流紋岩質粒間メルトを含むクリスタルマッシュ状の捕獲岩を発見した（富士山での石英の報告は初）著しく分化したケイ長質メルトが存在する

成果例．リモート分析ルーチンの立ち上げ

サブテーマC1

インターネット経由では，分析機器の操作は直接できないため，分析時には必ずオペレータが介在する．そのため，分析したい点座標や分析条件を機器操作者（オペレーター）に正確に伝達することが鍵となる．
分析依頼者
オペレーター（＠地震研）

申込

受付

薄片送付

画像取得通知分析点登録

薄片全体の画像取得

2019年度に作成した画像DBを利用

定量分析結果通知指定された分析の実行

分析点登録画面
点や矩形指定が行えるツールを新たに作成

指定位置座標はテキストファイルで保存＋画像にマーク

任意のサイズに部分拡大，COMPO像にSEM像を重ねて表示可能

薄片全体のSEM像，COMPO像の取得

個々の画像をWebで見ながら，正確な分析点指定ができる
オペレーターが同じ画面を見ながら作業可能サブテーマC3①地下におけるマグマ移動シミュレーション
b.マグマ移動過程シミュレーション技術開発
境界要素法により地形の形状を考慮した地殻変動計算手法の開発

地殻変動シミュレーション開発

火山周辺の観測点は火山体地形の上にあるため，観測データの解析において地形の影響を考慮することは重要．
境界要素法により地形を考慮した火山性地殻変動を求める手法の開発を進めている．
開発したプログラムを利用しやすくするため，数値計算に必要な山体地形モデルの作成や計算結果を作図用に変換する機能の開発を行った．
開発したプログラムを伊豆大島の気象庁観測点での地殻変動量の計算に適用した．
山体地形・圧力源のメッシュの作成
山体地形： 国土地理院の10mメッシュの数値標高モデルを利用
圧力源の形状：球状・回転楕円体を作成可能

地形を考慮した地殻変動計算の流れ

伊豆大島への適用例

山体地形のメッシュモデル作成
圧力源モデルの作成

節点数: 5224個
メッシュ数: 10436個

地形メッシュモデルの作成例

火山性地殻変動の数値計算
観測点ごとの地殻変動量の抽出
観測点ごとの結果の解析・描画

計算結果例
計算結果の利用
数値計算結果から観測点ごとの傾斜変化量や変位量を抽出することが可能
火山性地殻変動計算プログラムのGUI画面の例．
赤破線内は地形モデル作成のためのパラメータ．
圧力源の位置や大きさといった山体地形のメッシュモデルを作成するためのパラメータおよび計算をGUI上で行えるプログラムを作成した．

気象庁の傾斜計観測点における傾斜変化量を開発したプログラムと標高補正茂木モデルの場合について求めて比較した．
図．球状圧力源の深さを海抜下1000m(左)と海抜上500m(右)とした場合の観測点ごとの傾斜変化量．
二子山北西観測点の傾斜変化量は0.01倍して表示．
圧力源が海抜より上にある場合，数値計算の結果と標高補正茂木モデルによる傾斜変化の違いは大きくなる．

サブテーマC3b. 噴煙柱ダイナミクスモデルの開発

降灰シミュレーションの高精度化
噴出率：4×106 kg/s

噴出率：1×108 kg/s

距離に対する粒子サイズ分布パターンに違い

様々な噴火規模・大気条件の噴煙・火山灰シミュレーションを実施．解析を進めている．