資料２－１柏崎刈羽原子力発電所６号炉及び７号炉火山影響評価について平成27年6月5日東京電力株式会社1
2
1．火山影響評価の概要
> 「原子力発電所の火山影響評価ガイド」に従って，柏崎刈羽原子力発電所の火山影響評価を実施した。 >発電所に影響を及ぼし得る火山を抽出して，32火山を抽出した。 >抽出された火山の火山活動に関する個別評価を実施した結果，設計対応不可能な事象が発電所に影響を及ぼす可能性はなし。 >発電所に影響を及ぼし得る火山事象を抽出した結果，降下火砕物以外に影響評価すべき火山事象はない。考慮すべき降下火砕物の層厚は，文献調査，地質調査及びシミュレーションの結果から，30cmとした。余白
2．1　地理的領域内の第四紀火山>地理的領域内（発電所から半径160kmの範囲）における第四紀火山を文献調査等から81火山抽出した。
地理的領域内の第四紀火山36°N 7
本活火山総覧（第４版）」（気象庁(2013) ）による
>地理的領域内（発電所から半径160kmの範囲）における第四紀火山（81火山）のうち、完新世に活動を行っていない火山(65火山)について，将来の火山活動可能性が否定できない火山を抽出した。 >抽出にあたっては，「日本の火山（第３版）」（中野ほか(2013) ）による活動年代等をもとに，最後の活動終了からの期間が過去の最大休止期間より長いとみなせる場合は，将来の活動可能性がないと判断した。 >その結果, 黒岩山, 苗場山, 志賀高原火山群, 新潟金山, 黒姫山, 志賀, 飯縄山, 子持山, 四阿山, 白馬大池, 男体・女峰火山群,烏帽子火山群,鼻曲山,上廊下,鷲羽・雲ノ平及び環諏訪湖の16火山を,完新世に活動を行っていないが将来の火山活動可能性が否定できない火山として抽出した。
※1：「日本の火山（第３版）」（中野ほか(2013））による．※2：敷地から各火山までの距離．  ※3：「日本の火山（第３版）」（中野ほか(2013））, 「第四紀噴火・貴入活動データーベースVer. 1. 00」（西来ほか（2014））および「日本の主要第四紀火山の積算マグマ噴出量階段図」（山本（2014））等による。  ※4：最後の活動終了からの期間が過去の最大休止期間より長いと判断される場合。将来の活動可能性が無いと判断し、それ以外を将来の火活動可能性が否定できない火山として評価した。
完新世に活動を行っていない火山※1：「日本の火山（第３版）」（中野ほか(2013））による．※2：敷地から各火山までの距離．  ※3：「日本の火山（第３版）」（中野ほか(2013））, 「第四紀噴火・貴入活動データーベースVer. 1.00」（西来ほか（2014））および「日本の主要第四紀火山の積算マグマ噴出量階段図」（山本（2014））等による。  ※4 :最後の活動終了からの期間が過去の最大休止期間より長いと判断される場合。将来の活動可能性が無いと判断し、それ以外を将来の火山活動可能性が否定できない火山として評価した。
完新世に活動を行っていない火山※1：「日本の火山（第３版）」（中野ほか(2013））による．※2：敷地から各火山までの距離．  ※3：「日本の火山（第３版）」（中野ほか(2013））, 「第四紀噴火・貴入活動データーベースVer. 1.00」（西来ほか（2014））および「日本の主要第四紀火山の積算マグマ噴出量階段図」（山本（2014））等による。  ※4 :最後の活動終了からの期間が過去の最大休止期間より長いと判断される場合。将来の活動可能性が無いと判断し、それ以外を将来の火山活動可能性が否定できない火山として評価した。
>地理的領域内（発電所から半径160kmの範囲）における第四紀火山を文献調査等から81火山を抽出した。 >完新世に活動を行った火山として，16火山を抽出した。 >完新世に活動を行っていないが将来の火山活動可能性が否定できない火山として，16火山を抽出した。 ▷以上より，原子力発電所に影響を及ぼし得る火山として32火山を抽出した。
>発電所に影響を及ぼし得る火山（32火山）について，設計対応不可能な火山事象が運用期間中に影響を及ぼす可能性について検討した。 >検討は，設計対応不可能な火山事象の影響範囲と発電所から各火山への距離等に着目して行った。発電所に影響を及ぼし得る火山
注1：噴出中心と原子力発電所との距離が、表中の位置関係に記載の距離より短ければ、火山事象により原子力発電所が影響を受ける可能性があるものとする。注2 :降下火砕物に関しては、原子力発電所の敷地及び敷地什近の調査から求められる単位面積あたりの質量と同等の火山区等が降下するものとする。注3 :新火口の開口については、原子力発電所の運用期間中に、新火口の開口の可能性を検討する。注4：火山活動によるこれらの事象は、原子力発電所との位置関係によらず、個々に検討を行う。 3．1　設計対応不可能な火山事象と発電所の位置関係▷溶岩流及び岩屑なだれ等については，敷地から各火山までの距離が50km以上であることから，評価対象外とした。
>第四紀火山の噴出物分布図によれば，仮にこれらの噴出物が火砕物密度流だと考えても，噴出物の分布が山体周辺に限られることから，火砕物密度流が敷地周辺に到達していないと考えられる。
(中野ほか(2013)に一部加筆 )
深部低周波地震>深部低周波地震は，活動的な火山の周辺に限定的に分布している。 ▷敷地周辺に深部低周波地震の分布は認められない。
地理的領域の深部低周波地震気象庁一元化震源カタログ(1997年10月～2014年7月)
地温勾配>活動的な火山の周辺では，地温勾配が比較的大きい。 ▷敷地周辺の地温勾配は，50k/km程度と小さい。発電所に影響を及ぼし得る火山
地理的領域の地温勾配
地殻熱流量>地殻熱流量は，データが少ないが，草津白根山や焼岳周辺で高い値を示している。 >敷地周辺では,  このような大きな地殻熱流量は得られていない。まとめ>新しい火道の開通については，敷地周辺で深部低周波地震の活動がないこと，地温勾配が小さく，また地殻熱流量が小さいことから，新しい火口の開口が発電所に影響を及ぼす可能性はないと判断される。地理的領域の地殻熱流量
>設計対応不可能な火山事象（火砕物密度流, 溶岩流, 岩屑なだれ他, 新しい火口の開口及び地殻変動）が発電所に影響を及ぼす可能性はない。 >既往最大の噴火を考慮しても発電所に影響を及ぼさないと判断できることから、火山活動のモニタリングは不要と判断し
>発電所の安全性に影響を及ぼす可能性のある火山事象について抽出を行った。
>柏崎刈羽原子力発電所に降下した可能性がある広域火山灰は，町田・新井(2011)によれば，７層の分布が示されている。 ▷これらの火山灰については，カルデラ噴火など大規模噴火に伴って噴出したもので，将来同規模の噴火が発生し，敷地に影響する可能性は小さいと考えられる。
>敷地周辺で確認されている降下火山灰は, 以下の通りである。 >給源が特定出来る降下火砕物については，各火山の活動可能性を評価し，同規模の噴火が発生する可能性は小さいと評価した。 >給源不明の降下火砕物については，敷地周辺での分布状況を整理した。
※：中村正芳・新井房夫（1988）による。 4. 1 降下火砕物の影響可能性▷吉水テフラの分布吉水テフラの分布図(図中の数値は層厚（単位：cm）を表す） 4. 1 降下火砕物の影響可能性▷常楽寺テフラの分布常楽寺テフラの分布図(図中の数値は層厚（単位：cm）を表す）
▷不動滝テフラの分布不動滝テフラの分布図(図中の数値は層厚（単位：cm）を表す) 4. 1 降下火砕物の影響可能性>阿相島テフラの分布阿相島テフラの分布図(図中の数値は層厚（単位：cm）を表す） 4. 1
>給源不明なテフラは，西山層，灰爪層，魚沼層中に挟在しており，これら地層は底生有孔虫等の分析結果等から水成層とされており，これらのテフラも同様に水中で堆積したと考えられる。
>黒川(1990)では，新潟新生代堆積盆における上部中新統から下部更新統中に挟在する水底堆積珪長質テフラについて，その形成機構を検討しており，吉水テフラ,常楽寺テフラ,不動滝テフラおよび阿相島テフラは水底堆積テフラとされている。
>青木・黒川(1996)では，新潟県西頸城地域の鮮新統から下部更新統の火山灰層を調査し，当該地域に分布する谷浜層基底部の大菅パミス質[0op]が,    Hap-2 (阿相島テフラに対比）に対比されるおり，これらの層厚変化について検討されている。 >この結果から，阿相島テフラは長野県北部に噴出源が推定され，堆積盆中を北方向および北東方向に水底を流走したテフラと考えられている。以上より，給源不明なテフラは水中で堆積したテフラであり，その分布状況から堆積過程において前期更新世当時の水系等の影響を受けて堆積したものと推定され, 現在の堆積環境とは大きく異なっていることから,これら給源不明なテフラの分布状況に基づいて，将来の降下火砕物の層厚を想定するのは適切でないと考えられる。 (青木・黒川(1996) )
■給源不明の降下火砕物は，いずれも水成堆積物であり，現在の発電所周辺の堆積環境とは異なること，その噴出時期が前期更新世と古く，分布層厚は敷地周辺で大きくばらつきがあることから，プラント運用期間中に，その様な降下火砕物の堆積が敷地周辺に生じる蓋然性を確認するため，以下の観点で評価を実施した。 ○設置変更許可申請時点の評価項目※ ✓文献を用いた評価・噴火実績のある火山の降下火砕物の火口からの距離に応じた堆積速度による試算・噴火実績のある火山の降下火砕物の火口からの距離に応じた堆積量による試算✓既往解析結果の知見（火口からの距離に応じた堆積量による試算） ※：評価方法は，工事計画認可申請書にのみ記載○設置変更許可申請後，上記評価結果の妥当性を確認するため，追加評価を実施✓評価対象火山に対する降下火山灰シミュレーション▷評価対象火山は，降下火砕物によって発電所に影響を与える可能性がある，発電所周辺全ての火山を対象に抽出を実施。 .設置変更許可申請時点の評価対象火山は，火山のおおよその噴火規模（VEI単位），及び発電所との距離のみを考慮して抽出していたが， •　追加評価では，網羅性を高めるため，噴火時に想定される噴出量（Km³単位），発電所との位置関係（距離・方角）を考慮して抽出し，評価を実施した。余白

( Ο )評価対象火山の抽出及び堆積量評価フロー
抽出観点①②最大噴火規模がVEI5の内，発電所から最も近い妙高山より内側に，VEI4以下の火山が存在しないため，VEI4以下の火山は妙高山の評価に包含されるとした。なお，抽出にあたっては，過去の噴火履歴から，噴火活動が衰退傾向にあるものは， 「プラント運用期間中に想定される噴出量」として噴出量の精査を行った。一方,噴火活動が衰退傾向にないものは，過去最大規模の噴火が生じるとして抽出した。抽出観点③噴火規模，発電所からの位置関係より，他の火山の影響に包含されない火山を代表として抽出した。 4． 1． 1　降下火砕物の影響可能性（堆積量の評価）
※噴火規模、発電所からの距離・方角を考慮し、「北八ヶ岳」は「浅間山」に, 「燧ヶ岳」「榛名山」「男体・女峰火山群」「高原山」「那須岳」「安達太良山」は「沼沢」に包含される
▷評価対象火山の抽出（観点②）にあたり，噴火活動が減衰傾向である，飯縄山，草津白根山，高原山について噴火規模を精査した。
■抽出した６火山（妙高山，沼沢，四阿山，赤城山，浅間山，立山）の特徴について記す。
表評価対象火山の位置及び噴火規模堆積量の評価は以下の通り実施する。 (１)文献を用いた評価・噴火実績のある火山の降下火砕物の火口からの距離に応じた堆積速度による試算・噴火実績のある火山の降下火砕物の火口からの距離に応じた堆積量による試算( 2 )既往解析結果の知見（火口からの距離に応じた堆積量による試算） (３)評価対象火山に対する降下火砕物シミュレーション結果→（１），（２）については，評価方法の性質上，堆積量の結果が火口から発電所の距離に依存することから，発電所から最も近い「妙高山」に限り記載する。
( 1 － 1 )文献を用いた評価（堆積速度からの試算）妙高山と同等の噴火規模（VEI5）の実績を持つ富士山，並びに文献に火口から遠距離における堆積量の記載があるピナツボ火山（VEI6），タンボラ火山（VEI7）の噴火について，継続時間，火口からの距離に応じた堆積量が記載されている文献から，堆積速度を推定した。 ▷妙高山を想定した，火口からの距離（74km）における堆積速度は0.20cm/hとなった。

（１－１）文献を用いた評価（堆積速度からの試算）
発電所敷地内の堆積量を，推定した堆積速度に，富士山（宝永噴火）に関する文献に基づく噴火継続時間（16日間）を乗じ，発電所における風向割合（想定火山より約３割※の風が飛来することを想定）を考慮した結果，約23.1cmとなった。 ※：評価対象火山は発電所に対して，南西から東方向に位置していることを踏まえ，気象庁が行っているラジオゾンデの定期観測（観測地点：輪島）データより，原子力発電所へ吹く風は，最も多い場合でも，約３割であることから設定。発電所に対する風向割合文
(１－２)文献を用いた評価（堆積量からの試算）妙高山と同等の噴火規模（VEI5）の実績を持つ富士山の噴火について，火口からの距離に応じた堆積量が読み取れる文献より，堆積量を推定した。妙高山を想定した，以口からの距離における堆積量は約23cmとなった。富士山1707年「宝永噴火」推積深等高線K

(２)既往解析結果の知見
既往解析結果の知見として，富士山ハザードマップ検討委員会報告書に，富士山の宝永噴火に関する解析結果が示されている。当該の解析は，風向などの気象条件を考慮し，堆積量を推定している。発電所から最も近い火山（妙高山）を想定した，堆積量は約15cmとなった。
①解析コード「Tephra2」の概要▷ 「Tephra2」は移流拡散モデルを用いたシミュレーションプログラムであり、降下火砕物の降灰範囲の予測や既往噴火の降灰状況の復元を目的として利用されている。 >移流拡散モデルは図に示す通り，降下火砕物の挙動を重力による落下，風による移動（移流）及び空中で降下火砕物が自発的に散らばる現象（拡散）を計算するものである。 >風は高度毎に水平な一方向にふくものとされ，拡散も水平方向のみが考慮されている。 >降下火砕物は，火口上に仮定された均質な噴煙柱より放出される。 >その他の解析コードとの比較は以下のとおり。
閩「Tephra2」の移流拡散モデルと支配方程式
②－１評価条件■高度別風速及び風向データ▷気象庁が行っているラジオゾンデの定期観測データ（観測地点：輪島）より，30年間（1981～2010年）のデータを対象とし，高度別風速及び風向の月別平均値を入力データとした。 >風向は西風が卓越しているが，風向分布が2極化している場合，評価対象火山が発電所に対して西側※に位置しているため，東へ向かう風となるよう0° ～180°区間を対象に値を平均する補正を行った。 ※：風向は西側が卓越しているため、発電所より東側の火山の降下火砕物は堆積しないと考えられることから評価は実施しない。ただし，後述の不確かさを考慮した評価は実施する。風速データ風向データ（北方向0°,西方向90°,南方向180°,東方向270°）風速及び風向の平均値
■風速及び風向データ（続き）
②－１評価条件■噴煙柱高度について評価対象火山の想定される噴火規模はVEI5であることから、町田ら(2003)によれば噴煙柱高さは25km以上となる。ただし，風速データのピークは標高10～ 15km程度のため、噴煙高さを25kmで設定した。
年短歴(6: 50 19)
■その他の評価条件

※１：沼沢，赤城山のデータは，後述の不確かさを考慮した評価に使用する※２：ケイ質の噴出物に適応
②－２評価結果■妙高山解析結果より，いずれの月も火山位置から東方に向かって細長く火山灰が堆積し，発電所位置では、堆積量が１cm未満となった。図妙高山における堆積シミュレーション結果（１～４月）
②－２評価結果■妙高山（評価結果続き）図妙高山における堆積シミュレーション結果（５～１２月）
②－２評価結果■立山解析結果より，いずれの月も火山位置から東方に向かって細長く火山灰が堆積し，発電所位置では、堆積量が１cm未満となった。図立山における堆積シミュレーション結果（１～４月）
②－２評価結果■立山（評価結果続き）図立山における堆積シミュレーション結果（５～１２月）
②－２評価結果■浅間山解析結果より，いずれの月も火山位置から東方に向かって細長く火山灰が堆積し，発電所位置では、堆積量が１cm未満となった。
②－２評価結果■浅間山（評価結果続き）浅間山における堆積シミュレーション結果（５～１２月）図
②－２評価結果■四阳山解析結果より，いずれの月も火山位置から東方に向かって細長く火山灰が堆積し，発電所位置では、堆積量が１cm未満となった。
②－２評価結果■四阿山（評価結果続き）四阿山における堆積シミュレーション結果（５～１２月）図
③不確かさを考慮した評価（評価フロー）各パラメータ毎に評価を実施し，堆積量が保守的となるパラメータを確認する。最終的な堆積量は，確認したパラメータを組みあわせて算出する。 ○各パラメータの設定の考え方(a)評価対象火山から発電所に向かう風（±11.25°）を抽出し，評価を実施(b)基本ケースの風速を± σ，２σで増減させ，保守的な条件を設定(c)基本ケースの噴煙柱高さを±5kmで増減させ，保守的な条件を設定(d-1) tephra2推奨値の粒径範囲1/2 6～1/2 6に対して評価し，保守的な条件を設定(d-2)岩片・軽石粒子密度を1000kg/m³，2600kg/m³で評価し，保守的な条件を設定(d-3)噴煙放出下限高度を0～0.4で増減させ，保守的な条件を設定
③不確かさを考慮した評価主要なパラメータの不確かさの要素の考慮は以下のとおり。 (a)風向に関する不確かさ対象火山から発電所に向かう風を抽出し評価する。 (b)風速に関する不確かさ基本ケース（平均風速）に対して± σ（標準偏差）， 2 σを考慮する。 (c)噴煙柱高さに関する不確かさ基本ケース（25km）に対して±5kmを考慮する。

(a)風向に関する不確かさの考慮
風向に関する不確かさを考慮し，評価対象火山から発電所に向かう風（θ ±11.25°）を抽出し、シミュレーションを実施した。図風向抽出イメージ（妙高山）
解析は，風向に関する不確かさを考慮した結果をもとに，更にその他の不確かさを考慮して評価を実施する。
(b)評価結果(風速に関する不確かさの考慮)基本ケースの風速（平均風速）に± σ , ２σを増減させて評価を実施した。なお，±２σの評価で95％信頼区間を得られることから，＋３σ以上については参考に評価を実施。
(c)評価結果(噴煙柱高さに関する不確かさの考慮)基本ケースの噴煙柱高さに±５kmを増減させて評価を実施した結果は以下のとおり。
噴煙柱高さに関する不確かさを考慮したシミュレージョン結果国(d)評価結果(その他不確かさの考慮) ▷その他のパラメータについても，同様に感度解析を実施し，各パラメータについて，評価結果が最大になるパラメータの値を確認した。 (基本ケースから変更したパラメータは赤字，風速については，＋３σ以上を参考値として青字記載） ▷それらのパラメータを組み合わせ，最終的な堆積量を求めた。

(4)堆積量まとめ
発電所敷地周辺で確認された程度の降下火砕物が，プラント運用期間中に生じる蓋然性を確認するため，降下火砕物によって発電所に影響を与える可能性がある，発電所周辺全ての火山を対象に，網羅的に評価対象火山を抽出し，堆積量を評価した。結果としては，同程度の規模の降下火砕物が堆積するという結果は得られなかった。 (１) ～ (３)のうち評価結果が最も大きい， (１－１)の文献を用いた評価結果（約23.1cm）に対し、保守性を考慮した30cmを基準火砕物堆積量と設定した。
表 堆積量の評価結果一覧

(1)降下火砕物の粒径
降下火砕物の粒径は，評価対象火山の噴火規模と同等（VEI5）の噴火実績を持つ富士山及び樽前火山について，火口からの距離と粒径分布が記載された文献から評価した。評価対象火山のうち発電所から最も近い火山（妙高山（74km）），最も遠い火山（立山(131km) ）の距離を考慮し，0.17mm～8.0mmを基準以砕物粒径とした。
(2)降下火砕物の密度降下火砕物の密度は，以下の文献調査の結果から，保守性を考慮した1,5g/cm³を基準火砕物密度とした。 ▷アメリカ地質調査所（USGS）の文献によると、噴火時に想定される降下火山灰の乾燥状態の比重は，0.5～1.3g/cm3とされている。 ▷東京大学出版会の文献（宇井[編]（1997））によると「乾燥した火山灰は、密度が0.4～0.7程度であるが、湿ると1.2を超えることがあること」とされている。
>発電所に影響を及ぼし得る火山（32火山）を対象に，火山性土石流等，火山から発生する飛来物，火山ガス及びその他の火山事象の影響可能性を検討した。 >検討は，各火山事象の影響範囲と発電所から各火山への距離や，発電所周辺の地形等に着目して行った。
発電所に影響を与える可能性のある火山事象及び位置関係
注1：噴出中心と原子力発電所との距離が、表中の位置関係に記載の距離より短ければ、火山事象により原子力発電所が影響を受ける可能性があるものとする。注2 :降下火砕物に関しては、原子力発電所の敷地及び敷地什近の調査から求められる単位面積あたりの質量と同等の火山区等が降下するものとする。注3 :新火口の開口については、原子力発電所の運用期間中に、新火口の開口の可能性を検討する。注4：火山活動によるこれらの事象は、原子力発電所との位置関係によらず、個々に検討を行う。
▷発電所に影響を及ぼし得る火山，主に信濃川水系，関川水系および姫川水系の流域に属し，敷地周辺の河川を流域にもつ火山は存在しない。発電所に影響を及ぼし得る火山と日本海に注ぐ河川
>敷地周辺の河川はその流下方向が敷地へ向いておらず，敷地周辺の河川と敷地の間には地形的な高まりが認められることから，仮にこれらの河川の流域に降下火砕物が堆積しても二次的な泥流が発電所に影響を及ぼす可能性はない。 >敷地内で火山性土石流等の痕跡は認められない。 >以上から，火山性土石流、火山泥流及び洪水が発電所に影響を及ぼす可能性はない。
>火山性土石流等,火山から発生する飛来物,火山ガス及びその他の火山事象が発電所に影響を及ぼす可能性はない。

50.
【原子力発電所に影響を及ぼし得る火山の抽出】 ✓敷地を中心とする半径160kmの範囲には，81の第四紀火山がある。 ✓敷地を中心とする半径160kmの範囲の第四紀火山について，完新世の活動の有無，将来の活動可能性の検討を行い，原子力発電所に影響を及ぼし得る火山として，以下の32火山を抽出した。 (完新世に活動を行った火山)妙高山,新潟焼山,燧ヶ岳,沼沢,草津白根山,日光白根山,榛名山，赤城山，浅間山，高原山，那須岳，立山，磐梯山，吾妻山，北八ヶ岳，安達太良山（計16火山） (将来の火山活動可能性が否定できない火山)黒岩山,苗場山,志賀高原火山群,新潟金山,黒姫山，志賀，飯縄山，子持山，四阿山，白馬大池，男体・女峰火山群，烏帽子火山群，鼻曲山，上廊下鷲羽・雲ノ平及び環諏訪湖（計16火山） 【抽出された火山の火山活動に関する個別評価】 ✓敷地との距離，地形的条件，個別評価等の結果から，設計対応不可能な火山事象（火砕物密度流，溶岩流，岩屑なだれ他，新しい火口の開口及び地殻変動）が発電所に影響を及ぼす可能性はない。 ✓既往最大の噴火を考慮しても発電所に影響を及ぼさないと判断されることから，モニタリングの必要性はない。 【原子力発電所に影響を及ぼし得る火山事象の抽出】 ✓考慮すべき降下火砕物の層厚は，文献調査，地質調査及びシミュレーションの結果から，30cmとした✓火山性土石流,飛来物(噴石) ,火山性ガス及びその他の火山事象のうち影響を評価すべき事象はなしい。 ● 原子力規制委員会（2013）：原子力発電所の火山影響評価ガイド●中野俊・西来邦章・宝田晋治・星住英夫・石塚吉浩・伊藤順一・川辺禎久・及川輝樹・古川竜太・下司信夫・石塚治・山元孝広・岸本清行（2013）：日本の火山（第3版）,産業技術総合研究所地質調査総合センター ● 気象庁（2013）： 「日本活火山総覧（第４版）」 ● 西来邦章, 伊藤順一, 上野龍之, 内藤一樹, 塚本 斉（2014）：第四紀噴火・貫入活動データーベースVer.1.00 ● 山元孝広（2014）：日本の主要第四紀火山の積算マグマ噴出量階段図. 地質調査総合センター研究資料集, no. 613,産総研地質調査総合センター, ●気象庁一元化震源カタログ(1997年10月～2014年7月) ●産業技術総合研究所地質調査総合センター（2004）：日本列島及びその周辺域の地温勾配及び地殻熱流量データーベース● 町田洋・新井房夫（2011）：新編火山灰アトラス［日本列島とその周辺］（第2刷）. 東京大学出版会● 中村正芳・新井房夫（1988）：群馬県中央部で発見された前期更新世の含董青石テフラについて, 地球科学, vol.52, pp.153-157 ● 文献:新潟県（1977） : 20万分の1新潟県地質図・同説明書,  493p. ● 宮下ほか(1972) ：日本油田・ガス田図7. 「魚沼」地質説明書.地質調査所, 36p. ● 天然ガス鉱業会・大陸棚石油開発協会（1982）：日本の石油・天然ガス資源.  455p ●黒川勝己（1990）水底に堆積した珪長質テフラの層相モデルとその形成機構ー新潟地域の例一.地球科学, 44, P361-378. ●岸清・宮脇理一郎(1996)新潟県柏崎平野周辺における鮮新世～更新世の摺曲形成史.地学雑誌, 105, P88-112. ●安井賢・小林巌雄・立石雅昭(1983) :新潟県八石油帯・中央油帯に分布する魚沼累層の層序.地球科学, 37, P22–37. ● 新潟平野団体研究グループ(1970) :新潟県小国町の魚沼層群— 新潟県の第四系・そのXII—.新潟大学教育学部高田分校研究紀要, 15, P267-302. ● 黒川勝巳・丸山悦子・沢栗隆之(1989) :新潟県中央油帯北部における椎谷層・西山層中の水底堆積テフラ.新潟大学教育学部紀要(自然科学編), 30, P39–64. ● 沢栗美香子・黒川勝巳(1986) :水底堆積テフラからみた鮮新・更新世火山活動の変遷— 新潟堆積盆小国町地域の例—.地球科学, 40, P117-191. ● 青木豊樹・黒川勝己(1996) ：新潟県西頸城地域の鮮新統〜下部更新統の火山灰層とその対比.地球科学, 50,  P341〜361. ● 柳沢幸夫 他（2001）：飯山地域の地質● 島津光夫 他(1993) ：苗場山地域の地質● 早津賢二（2008）：妙高火山群一多世代火山のライフヒストリー一● 早津賢二（1994）：新潟焼山火山の活動と年代一歴史時代のマグマ噴火を中心として一● 竹内圭史  他（1994）：20万分の１地質図「高田」地質調査総合センター ● 須藤 茂 他（2007） :わが国の降下火山灰データベース作成● 山元孝広（2012）：福島一栃木地域における過去約30万年間のテフラの再記載と定量化● 五十嵐聡(1984) :新潟,長野県境付近の津南一志賀地域の鮮新一更新世の火山岩類● 山元孝広（1999）：福島一栃木地域に分布する30–10万年前のプリニー式降下火砕物：沼沢・燧ヶ岳・鬼怒沼・砂子原火山を給源とするテフラ群の層序● 早川由起夫 他(1989)  :草津白根火山の噴火史● 産業技術総合研究所：日本の活火山1万年噴火イベントデータ集(https://gbank.gsj.jp/volcano/cgi-bin/volcanic.cgi?id=042) ● 鈴木毅彦 他（1994）：テフラからみた日光火山群の噴火史● 飯塚義之（1996）：子持火山の地質と活動年代● 大石雅之（2009）：四阿火山を起源とする噴出物の岩石記載的特徴とテフラ分布柵山雅則（1980）：白馬大池火山の地質● 山元孝広（2013）：栃木一茨城地域における過去約30万年間のテフラの再記載と定量化● 村本芳英（1992）：日光火山群東方地域に分布する中・後期更新世テフラー日光火山群の噴火史一● 高橋 康（2004）：長野県北東部鳥帽子岳とその周辺の地質と火山形成史● 早川由紀夫（2010）：浅間山の風景に書き込まれた歴史を読み解く● 鈴木毅彦（1992）：那須火山のテフロクロノロジー ● 及川輝樹（2003）：飛驒山脈の隆起と火成活動の時空的関連● 木村純一（1987）：長野県における後期更新世の降下火山砕屑物層序● 及川輝樹 他(2003)  :飛騨山脈中央部, 上廊下～雲ノ平周辺の第四紀火山岩類のＫ−Ａｒ年代● 大石雅之  他（2004）：八ヶ岳火山を起源とする新期テフラ群の層序と噴火史● 山元孝弘 他（2000）：テフラ層序からみた安達太良火山, 最近25万年間の噴火活動● Kuniaki Nishiki et al (2011) : Quaternary volcanism and tectonic history of the Suwa- Yatsugatake Volcanic Province, Sentral Japan ● 田島広一他(1977) :浅間・草津白根山周辺の重力測定● 村上  亮 他（2004）：GPS測量および水準測量結果が示唆する草津白根山の収縮源●金子隆之他（1991）：信越高原地域に分布する第四紀火山のK-Ar年代と形成史● 金子隆之 他(1989)  : K-Ar年代から見た信越高原地域の火山活動他（2009）：塩原カルデラより噴出した大規模火砕流堆積物群の層序と年代● 弦巻賢介●伴雅雄他（1992）：東北本州弧，高原火山噴出物の地球化学●奥野充他(1997) :北関東,高原火山の約6500cal yr BPの噴火● Sigurdsson and Carey (1989) : Plinian and co-ignimbrite tephra fall from the 1815 eruption of Tambora volcano ● Lynn et al (1997) : Tephra Falls of the 1991 Eruptions of Mount Pinatubo ● 宮地直道（1984）：富士火山1707年火砕物の降下に及ぼした風の影響宮地直道 他（2011）：富士火山1707年噴火の推移とその噴出物の物理化学的特性の経時変化宮地直道他（2007）：富士火山1707年噴火（宝永噴火）についての最近の研究成果