四国電力，四国総合研究所
研究期報１０３（2015年12月）
９－２１ページ

四国に影響を及ぼす降下火山灰に関するシミュレーション解析㈱四国総合研究所 土木技術部㈱四国総合研究所 土木技術部㈱四国総合研究所 土木技術部
四国電力㈱
土木建築部
四国電力㈱
土木建築部
キーワード：

降下火山灰
シミュレーションTEPHRA2
由布岳
九重山
阿蘇山

Key Words：

辻
智大
池田 倫治
加藤佐代正
西坂 直樹
大西 耕造

air-fall ash
Simulation
TEPHRA2
Yufu Volcano
Kuju Volcano
Aso Volcano

Simulation analysis for air-fall ash influences on Shikoku Island.
Shikoku Research Institute, Inc., Civil Engineering Department
Tomohiro Tsuji, Michiharu Ikeda, Sayomasa Kano
Shikoku Electric Power Co., Civil Department
Naoki Nishizaka, Kozo Onishi
Abstract
If the large volcanic eruption occurred at a volcano on Kyushu, the volcanic ash would be transported to Shikoku
by the westerlies and affect to electric facilities at the downwind region, although there is no volcano in Shikoku.
Especially for nuclear power plants, it is necessary to assess the ash fall hazard to the facilities, considering thetargeted volcanoes. The thickness of volcanic ash around Ikata Power Plant has been evaluated to be almost 0 cm
based on geological survey. However, it is expected that the ash thicker than a few cm will be fallen on the Ikata site,
depend on wind direction. Then we carried out the analytical study using TEPHRA2, simulations of tephra
dispersal.
The tephra fall hazard from a future eruption at Yufu, Kuju and Aso volcanoes are simulated. The most effective
eruption to the site is the VEI5 eruption from Kuju Volcano (Kj-P1, the eruption mass is 2.03 km3). Based on thecalculation using the average wind of each month, the average thickness of the ash fallen on the site is 0.5 cm. The
ash fallen on September was thickest (2.2 cm) in twelve months.
In the case that the wind blows from the vent to the site, 4.5 cm-thick-ash fallen on the site. Given that theeruption mass being 6.2 km3, the ash of 6.9 cm-thick deposited at the site.

－9－四国に影響を及ぼす降下火山灰に関するシミュレーション解析

1.

はじめに

り，堆積物の分布を計算できる（萬年，2013）．

四国には火山が分布しないものの，風上に短時間で簡易的に計算することができるため，位置する九州の火山が大規模な噴火を起こしパラメータスタディおよび確率論的評価に適た場合には，偏西風に乗って四国へ火山灰がしている．解析プログラムは南フロリダ大学

降下し，各種電力設備に影響を及ぼす可能性のがある．特に原子力発電所の安全対策には万（ http://www.cas.usf.edu/~cconnor/vg@usf/t

全を期すことが必須であり，原子力規制委員

ephra.html）からダウンロードすることがで会制定の「原子力発電所の火山影響評価ガイ

きる．計算には，噴火パラメータ，大気パラド」において，原子力発電所に影響を及ぼし

メータ，粒子パラメータ，グリッドパラメー

得る火山を抽出して発電所運用期間中の噴火

タの 4 つが必要となる．計算結果として，単規模に対応する降下火山灰を考慮することが位面積当たりの降灰量（kg/m2）とその粒度組求められている．

成が得られる．

サ

イ

ト

降下火山灰による影響評価においては，まず敷地付近の地質調査により評価を行うが，3. 伊方発電所と四国に影響を及ぼす可能性の

敷地南東の宇和盆地における UT コア（守田

ある活火山の地理的関係

ほか，2014）など連続的な堆積物中にも考慮

伊方発電所の風上に位置する九州中部～

する過去の噴火と対応する火山灰層が認めら

東部には，北東から南西に，鶴見岳・伽藍岳

れないため，敷地付近への火山灰の降下厚さ（敷地からの距離 85 km），由布岳（89 km），

はいずれもほぼ 0 cm と評価される．しかし，九重山（108 km），阿蘇山（131 km）などの風向きによっては敷地において数 cm 以上の活火山が分布している（Figure 1）．その中で厚さの火山灰が降下することも想定される．

も阿蘇山および九重山は過去に大規模な噴火

そこで，気象条件の不確かさを考慮した火山をしており，九重山の噴火による火山灰は過灰降下リスクを評価するために，降下火山灰

去に四国にも降下している（熊原・長岡，

シミュレーションを用いた解析的検討を実施

2002）．また，伊方発電所からの距離が比較する．

的近い由布岳は中規模の噴火であるが火山灰を放出している．鶴見岳・伽藍岳は溶岩を放出する噴火を主体としており，遠方での火山2. 研究手法
本研究では，降下火山灰シミュレーション

灰の降下は報告されていない．そこで，本研として実績のある TEPHRA2 を利用する．

究では，由布岳，九重山および阿蘇山の 3 火

TEPHRA2 は移流拡散モデル（Bonadonna et

山を対象火山とし，それぞれの火山におけるal., 2005）に基づいたシミュレーションであ既往最大の噴火を想定した降下火山灰シミュ

り，適当な初期パラメータを与えることによレーションを行った．

－10－

四国に影響を及ぼす降下火山灰に関するシミュレーション解析

Fig. 1 Localities of Yufu, Kuju and Aso volcanoes with Ikata Power Plant.
4. 由布岳，九重山および阿蘇山の噴火による

km3（宮縁ほか，2003））である（Table 1）．
由布岳 2 ka 噴火は，噴火規模を示す爆発

影響の検討

的噴火指数（Volcanic Explosivity Index，

4.1. 検討方針
本検討では，まず，由布岳，九重山および

Newhall and Self, 1982）で表すと VEI = 3

阿蘇山の 3 つの火山の噴火による伊方発電所～4 となる．Newhall and Self（1982）によ周辺での降灰層厚を算定し，それぞれの火山

ると，その規模の噴煙柱高さは 3～25 km との噴火による伊方発電所周辺への影響を比較推定される．同様に，九重第一降下軽石およする．その後，伊方発電所への降灰に最も影び草千里ヶ浜降下軽石の噴火規模は VEI = 5

響を及ぼす火山噴火についてケース検討を実程度であり，噴煙柱高さは 20～35 km と推定

施する．

される．本検討では，由布岳による 2 ka 噴火の噴煙柱高さを 15 km，九重山および阿蘇山の噴火による噴煙柱高さを 25 km として降下

4.2. 噴火パラメータ
それぞれの火山の既往最大噴火は，由布岳

火山灰シミュレーションを行った（Table 1）．

の 2 ka 噴火（噴火年代約 2000 年前，降下火各火山噴火の粒子径（最大粒径，最小粒径，

砕物の噴出量 0.05 km3（須藤ほか，2007）），

平均粒径および粒子分散）については，現時

九重山の九重第一降下軽石（噴火年代約 5.4

点で情報が得られていないため，ここでは，万年前（川辺ほか，2015），噴出量 2.03 km3Bonadonna et al.（2005）に記載されている（須藤ほか，2007），阿蘇山の草千里ヶ浜降

TEPHRA2 推奨値を適用した（Table 2）．

下軽石（噴火年代約 3.1 万年前，噴出量 2.39

－11－

四国に影響を及ぼす降下火山灰に関するシミュレーション解析

Table 1 Eruption parameters used in the simulations.

Volcanoes
Height above
sea level (m)Distance from
Ikata site (km)Eruptive
Mass (km3)
Plume Height(km)
Column Steps

Kuju

Aso

(Nakadake)(Nakadake)

1,583

1,791

1,506

89

108

130

0.05

2.03

2.39

15

25

25

100

100

100

-10

-10

-10

10

10

10

Yufu

Maximum
Grain Size(phi)
Minimum
Grain Size(phi)

References
Japan Meteorological
Agency (2013)

Sudo et al. (2007),
Miyabuchi (2011)
Newhall and Self (1982)

Bonadonna et al. (2005)(Felsic: -10, Mafic:-7)
Bonadonna et al. (2005)(Felsic: 10, Mafic: 7)
Plinian eruption (Mount St

Median Grain
Size (phi)

1.0

4.5

4.5

Helens): 4.5 phi; Subplinian
eruption: 1 phi (TEPHRA2
recommendation)
Plinian eruption (Mount St

Sorting (phi)

1.5

3.0

3.0

Helens):3;

Subplinian

eruption: 1.5 phi (TEPHRA2
recommendation)付近の対流圏上層にジェット気流と呼ばれる

4.3. 大気パラメータ
本検討では，至近 22 年間(1988～2010 年)

風速が大きい偏西風が流れている（Figure 2）．

の毎日午前 9 時の地表～高度約 30 km までの特に冬季には風速が大きく，平均的に 30 m/s

各高度（hPa を高度に換算）の風向・風速観程度になり，中には 70 m/s 程度を超えること

測値を月別に平均化した月別平年値（気象庁もある．なお，夏季は冬季の半分程度以下の「高層気象観測年報」）を使用した（Figure 2）．

風速となる．上述した 3 火山の上空において解析には，3 火山および伊方発電所に最も近も，このジェット気流が流れていると考えらい福岡における高層気象観測データを用いたれる．このため，火山が噴火した場合，火山（Figure 1）．我が国では，上空約 8～15 km

灰がこのジェット気流により東方向に運搬さ－12－四国に影響を及ぼす降下火山灰に関するシミュレーション解析

れ，季節によって降灰領域が変化することが較するため，これらのパラメータを 3 つの火予想される．

山で共通して用いた．

4.4. 粒子パラメータ

4.5. グリッドパラメータ

粒子パラメータの設定には，軽石密度，岩本検討では，国土地理院の数値地図（50

片密度，拡散係数（Diffusion Coefficient），

mDEM）と産総研の数値地図（1 kmDEM）

落下時間の閾値（Fall Time Threshold）およを基に 2 km グリッドメッシュを作成し，こび噴煙モデル等の各データが必要である．これを利用した．なお，計算開始点は，噴火がれらについては萬年（2013）を参考にした

既存の火口から始まると仮定した．火口の座（Table 2）．それぞれの噴火による影響を比標については気象庁（2013）を参考にした．

Fig. 2 Upper-air wind velocity (a) and wind direction (b) averaged every month sampled
every 09:00 local time from 1 March 1989 through 31 December 2010 at Fukuoka (Fig.1a)(Japan Meteorological Agency, 2012).
Table 2 Grain parameters used in the simulations

Parameters

Unit

Value

References

Lithic Density

kg/m3

2600

TEPHRA2 recommendation

Pumice Density

kg/m3

1000

TEPHRA2 recommendation

Diffusion Coefficient

m3/s

200

Mannen (2013)

Eddy Constant

m2/s

0.04

Suzuki (1983)

s

3600

Mannen (2013)

Plume Model0 (uniform)

Mannen (2013)

Plume Ratio

0.1

TEPHRA2 recommendation

Fall Time Threshold－13－

四国に影響を及ぼす降下火山灰に関するシミュレーション解析

5. 由布岳，九重山および阿蘇山の噴火による（Table 3）．伊方発電所における降灰層厚は，影響の比較

いずれの火山噴火でも 9 月に最大となり，由由布岳，九重山，阿蘇山の噴火について，

布岳，九重山，阿蘇山でそれぞれ 0.3 cm，2.2

それぞれ月別平年値の風を用いて解析した結cm，1.0 cm である（Figure 4）．3 つの火山

果の概要を Table 3 に示す．なお，各表にお噴火による層厚が異なる要因としては，各火ける降灰層厚は，降灰の密度を 1000 kg/m3

山と敷地との位置関係，噴火規模および風向と仮定して換算した値である（例えば，降灰が挙げられる．由布岳は 3 つの火山の中で最量 1 kg/m2 だと降灰層厚は 0.1 cm となる）．

も伊方発電所までの離隔が小さいものの，噴月別の解析の結果，冬季（11 月～5 月）に

火規模が小さいため， 9 月を含めた各月の平おける伊方発電所における降灰層厚は 0～0.2

年風の風速では敷地に厚く降灰しない．阿蘇cm であり，夏季（6 月～10 月）と比較して，

山は 3 つの火山の中で最も伊方発電所との離格段に薄い（Table 3）．この傾向は 3 つの火隔が大きく，位置が南に位置していることか

山噴火で共通していたため，九重山の噴火にら，西風の卓越する冬季には敷地にほとんどよる結果のみを Figure 3 に示す．由布岳，九降灰しない．また，西南西の風が卓越する 9重山および阿蘇山の火山噴火による敷地での月においても降灰の中心は敷地より南側を向

降灰層厚はそれぞれ月平均で 0.1 cm，0.5 cm

いているため，敷地での降灰は 1.0 cm に留ま

および 0.2 cm であり，九重山の噴火で最大とる（Table 3）．

なり，阿蘇山，由布岳と小さくなっている
Table 3 Calculated thickness of volcanic ash fallen on Ikata Power Plant. Yellow marker
shows the maximum thicknesses.
Volcanoes

Yufu

Kuju

Aso

Jan

0.0

0.0

0.0

Feb

0.0

0.0

0.0

Mar

0.0

0.0

0.0

Apr

0.0

0.1

0.1

May

0.1

0.2

0.1

Thickness

Jun

0.2

0.6

0.2

at Ikata

Jul

0.2

1.5

0.7

Site (cm)

Aug

0.0

1.1

0.7

Sep

0.3

2.2

1.0

Oct

0.2

0.2

0.1

Nov

0.0

0.1

0.0

Dec

0.0

0.0

0.0

Average

0.1

0.5

0.2

－14－

四国に影響を及ぼす降下火山灰に関するシミュレーション解析

Fig. 3 Isopack maps of ash from Kuju computed for each month.

Fig. 4 Isopack maps of ash deposits from Yufu, Kuju and Aso computed for September.

－15－

四国に影響を及ぼす降下火山灰に関するシミュレーション解析

6. 九重山における噴火・大気パラメータの検討

九重第一降下軽石の噴火規模については，敷地への降灰に最も影響が大きい九重山

これまで火山灰の等層厚線図を基に噴火規模の噴火に関し，噴煙柱高度および噴出量等のを見積もった須藤ほか（2007）による 2.03

噴火パラメータ，大気パラメータを変化させ

km3 が示されてきた．しかし近年，長岡・奥たケース検討を行い，各パラメータの影響に野（2014）は給源付近における層厚分布のデついて考察し，敷地周辺における火山灰層厚ータを拡充し，算出方法は示していないものの評価を行う．その際，12 ヶ月のうち，敷地の，九重第一軽石の噴火規模として 6.2 km3

で最も厚い火山灰厚さとなる 9 月の風を用いを提案した．自然現象の評価と将来予測にはた降下火山灰シミュレーションを基本ケース

不確かさが残るため，噴出量を基本ケースのとする（Table 4）．

2.03 km3 対して大きく見積もった場合の噴出量 6.2 km3 を用いて検討を行った（Table 4）．

6.1. 噴煙柱高度のケース検討
九重第一降下軽石を放出した噴火は VEI

6.3. 風向・風速のケース検討= 5 であり，その噴煙柱高さは～25 km と推実際の気象条件を考えると，風速を固定し定される．ただし，同程度の規模の噴火であたまま風向のみを変化させることは現実的でっても，噴煙柱高度はばらつきの大きいパラはない．そこで，高層気象観測年報の 9 月のメータであることから，基本ケースである 25

平年風から，敷地方向に吹く風を抽出し，そkm から±5 km 変化させたケース検討を行っ

れらを平均化して仮想的な風を作成して風向た（Table 4）．以下の検討で用いた粒子パラの影響を検討した（Table 4）．尚，風速のばメータは Table 2 と同じである．

らつきの影響を検討するため，高層気象観測
年報の 9 月の平年風から，風速を±1σ変化させた検討も行った（Table 4）．

6.2. 噴出量のケース検討

Table 4 Eruption and atmospheric parameters used for parameter studies. Yellow markers
show the tested parameters in each case.

Cases
Plume
Height (km)Eruption
Mass (km3)Wind
Direction
Wind
Velocity

PrimalPlume

Eruption

Wind Direction

Wind Velocity

Case

Height Study

Mass Study

Study

Study

25

20, 30

25

25

25

2.03

2.03

6.2

2.03

2.03

Ave. Sep.

Ave. Sep.

Ave. Sep.

Blow to Ikata

Ave. Sep.

Ave. Sep.

Ave. Sep.

Ave. Sep.

Ave. Sep.

－16－+/-1 sigma ofAve. Sep

四国に影響を及ぼす降下火山灰に関するシミュレーション解析

7. 結果
7.1. 噴煙柱高度の検討

7.3. 風向・風速の検討

噴煙柱高度を変化させた結果を比較する

風速を変化させて，敷地周辺に堆積する火と，噴煙柱高さ 20 km の場合でも 25 km の山灰層厚がどのように変化するか検討した．

場合でも，伊方発電所における降灰層厚は 2.2

基本ケースの風速-1σおよび+1σの場合，敷cm で変わらない（Figure 5）．噴煙柱高度 20地での降灰層厚はそれぞれ 2.6 cm および 1.4

km の結果の方が噴煙柱高度 30 km 結果より

cm となった（Figure 7）．基本ケースでは，も，厚く降灰する領域（例えば層厚 25 cm 以

層厚分布の軸は敷地よりやや南側を通ってお上の領域）が遠くまで及んでいるものの，全り，軸から離れるほど降灰層厚は薄くなる．

体的な傾向としては，噴煙柱高度 20～30 km

風速を-1σ小さくした場合，基本ケースよりの範囲では，噴煙柱高度の変化が降灰層厚分

も降灰分布の幅が広くなるため，基本ケース

布に及ぼす影響は比較的小さい（Figure 5）．

よりも敷地における降灰層厚が厚くなった．
一方，風速を+1σ大きくした場合，基本ケー
スよりも狭小な降灰分布となるため，軸から7.2. 噴出量の検討
噴出量を 6.2 km3 にした場合，9 月の平年値の風で敷地周辺における層厚は 6.9 cm と離れると層厚が極端に薄くなり，敷地における降灰層厚が薄くなった．

なった（Figure 6）．噴出量 2.03 km3 の場合，風向が敷地向きとなる仮想的な風を想定

敷地での降灰層厚は 2.2 cm であり，噴出量のした場合，敷地周辺における層厚は 4.5 cm と変化とともに層厚も基本ケースのおよそ 3 倍なった（Figure 8）．風向が層厚に大きく影響程度となった．

することを確認できた．

Fig. 5 Isopack maps of ash from Kuju computed for column heights of 20, 25 and 30 km.

－17－四国に影響を及ぼす降下火山灰に関するシミュレーション解析

Fig. 6 Isopack maps of ash from Kuju computed in the case of 2.03 km3 (left) and 6.2 km3
mass volume (right).

Fig. 7 Isopack maps of ash from Kuju computed for the wind of slower (-1σfrom averagewind of September) (left), of September and of faster (+1σfrom that of September ) (right).

Fig. 8 Isopack maps of ash from Kuju computed for the average wind of September (left) andof the day when the wind blow to Ikata Power Plant (right).

－18－

四国に影響を及ぼす降下火山灰に関するシミュレーション解析

ら東北東方向に 108 km 離れた伊方発電所に

8. 考察
今回実施した九重山の噴火に関する検討

火山灰が厚く堆積する可能性は極めて低い．

ケースの中で，降灰層厚分布に大きく影響すまた，TEPHRA2 では，ジェット気流だけるパラメータとして，風向・風速および噴出でなく，それより低層の風の変化も結果に大量が挙げられる．噴出量を 2.03 km3 から 6.2

きく影響する．これは TEPHRA2 が均質な噴

km3 へ約 3 倍にすると，降灰量は単純に約 3煙柱を想定しており，そこから火山灰粒子が倍になった（Figure 6）．しかし，風向が敷地一様の落下速度で落下しながら，各標高に仮方向へ向かなければ，
たとえ噴出量が 6.2 km3想した水平層における風向・風速の元で拡散であっても降灰層厚はほぼ 0 cm である．そうすることに起因している．すなわち，ジェッ

いった意味では，敷地における降灰層厚に最ト気流の風は噴煙柱の高層から落下した粒子も影響を与えるパラメータは風向・風速であのみに影響を与えるが，低層の風は噴煙柱のる．特に風向による影響が大きいが，風向と低層から落下した粒子にも高層から落下した

風速は相互に関連しているため，両者を独立粒子にも影響を与えるためである．
噴煙柱高度は，降灰量分布にはさほど大きして扱うことはできない．
例えば，冬季は強いジェット気流の影響でく影響しなかった（Figure 5）．これは，高度指向性が強くなり，層厚分布が狭小となるた

20 km 以上では風速が小さいため，20 km 以

め，火山灰が厚く堆積する領域は少なくなる．

上における噴煙柱高度の変化が結果に影響しこれに加えて，西～西北西の風で安定するこにくいためと考えられる．しかし，本来は噴とにより（Figure 2），伊方発電所より南側へ

煙柱高度が高くなるということは噴出率が大降灰する（Figure 3）．それに対して，夏季はきくなり，傘型噴煙を形成することでより広ジェット気流が弱まり，風向がばらつくため，域に火山灰が広がる事が考えられる（例えば，火山灰が堆積する領域は広がる．一方で，風1991 年 Pinatubo 噴火）．当然，それによって速が弱まるため，火山灰が厚く堆積する領域

降灰量分布も変わるはずである． TEPHRA2

は火口近傍（九州付近）に留まる傾向がある

ではこういった傘型噴煙からの粒子の落下を（Figure 3）．9 月の風は，Figure 2 に示す通模擬できていないことが指摘されており，現り，高度 15 km より低い標高の風がほぼ一様

在解析コードの改善が進められているところに 240°（北方から時計回りの角度）より若である（萬年，2013）．噴煙柱高度の設定に干大きい方角から吹いていることから，このついては，今後，野外での噴出物の粒度分布

方位の風下側に位置する伊方発電所の降灰層に基づき推定する（Carey and Sparks, 1986）

厚に大きく影響していると考えられる．

ことが重要である．

本研究では，地表から高度約 30 km まで
を 25 程度の高度に区切って計測された気象9. まとめ

庁の風向・風速データを用いたが，各高度に伊方発電所の火山影響評価においては，原

おいて風向が敷地方向に揃うというのは極め

子力発電所に影響を及ぼし得る火山を抽出して稀な気象条件であり，実際には高度毎の風て，発電所運用期間中の噴火規模に対応する

向・風速は不揃いである．そのため，火口か降下火山灰を考慮する必要がある．そこで，－19－

四国に影響を及ぼす降下火山灰に関するシミュレーション解析

考慮対象である由布岳，九重山および阿蘇山の噴火について，月別平年値の風で降下火山[参考文献]

灰シミュレーションを行った結果，伊方発電1) 守田益宗・須貝俊彦・古澤明・大野裕

所に最も影響のある噴火は九重第一降下軽石記・西坂直樹・辻（噴出量 2.03 km3）であり，敷地における降田

下火山灰厚さは月平均で 0.5 cm，最大である

みた MIS15 以降の植生変遷史，植生史研

9 月の平年値の風で 2.2 cm であった．また，

究，23，3-19．

智大・池田倫治・柳誠，2014，宇和盆地の花粉分析から正確な見積もりが難しいことや最近の研究動

2) Bonadonna, C., C. B. Connor, B. F.

向を踏まえて噴出量 6.2 km3 と噴火規模を大Houghton, L. Connor, M. Byrne, A.

きく見積もった場合，最大である 9 月の平年

Laing,

値の風で敷地における降下火山灰厚さは 6.9

Probabilistic

cm であった．さらに，敷地方向に吹く仮想的

dispersal: Hazard assessment of a

な風を考慮した場合，敷地における降下火山

multiphase

灰厚さは 4.5 cm であった．

Tarawera, New Zealand. Journal of

and

T.

K.

Hincks,modeling
rhyolitic

of

2005,
tephraeruptionatGeophysical Research, 110, B03203,

四国西部には伊方発電所以外にも多数の

doi:10.1029 / 2003JB002896.

電力設備が立地しており，本研究で示した降下火山灰シミュレーション結果は，送電鉄塔，

3) 萬年一剛，2013，降下火山灰シミュレー

変電所，水力発電所等において火山影響評価

ションコード Tephra2 の理論と現状－第を行う上でも有用である．今後，九重第一軽四紀学での利用を視野に，第四紀研究，石について，給源近傍における降下火山灰厚52，173-187．

さ，粒度分布を詳細に把握するとともに，四4) 熊原康博・長岡信治，2002，四国南西部，

国南西部宿毛市において報告された九重第一松田川流域における九重第一テフラの対軽石の産状についてデータを拡充し，フィー比と低位段丘の年代，第四紀研究，41，

ルドデータに基づき降下火山灰シミュレーシ

213-219．
5) 須藤茂・猪股隆行・佐々木寿・向山栄，

ョンを高度化していくことが重要である．

2007，わが国の降下火山灰データベース
作成，地質調査研究報告，58，261 – 321．

[謝辞]
本研究を実施するにあたり，神奈川県温泉6) 川辺禎久・星住英夫・伊藤順一・山崎誠

地学研究所の萬年一剛博士および西日本技術

子，2015，九重火山地質図．火山地質図

開発株式会社の稲倉寛仁博士に TEPHRA2 に

19，産業技術総合研究所地質調査総合セ

ついてご教示頂いた．また，長崎大学故長岡ンター．

信治教授および電力中央研究所の三浦大助博

7) 宮縁育夫・星住英夫・高田英樹・渡辺一

士には降下火山灰に関してご教示をいただい

徳・徐勝，2003，阿蘇火山における過去た．㈱四国総合研究所の鈴木富美子氏には約９万年間の降下軽石堆積物，火山，48，

TEPHRA2 解析および図化に御協力いただい195-214．

た．ここに記して深く感謝いたします．

8) Newhall, C. G. and Self, S., 1982, The－20－

四国に影響を及ぼす降下火山灰に関するシミュレーション解析

volcanic explosivity index (VEI): An

dispersion of tephra. In Shimozuru, D.

estimate of explosive magnitude for

and Yokoyama, I. (Eds) Arc Volcanism:

historical

Physics and Tectonics, 95-116.

Geophysicalvolcanism.
Research,

Journal
87

of
(C2),

12) 長岡信治・奥野充，2014，九重火山のテ

1231-1238.

フラ層序，月刊地球，36，281-296．

9) 気象庁編，2013，日本活火山総覧（第 413) Carey, S. N. and Sparks, R. J., 1986,
Quantitative models of the fallout and版），気象業務支援センター．
10) 気象庁編，2012，高層気象観測年報，気象業務支援センター．
11) Suzuki, T., 1983, A theoretical model for－21－

dispersaloftephraeruptioncolumns.

fromVolcanology, 48, 109-125.

volcanic

Bulletinof