温泉科学（J. Hot Spring Sci.），63，164-183（2013）

原著

温泉水・温泉付随ガスの地球化学データから見た大分県山香温泉の生成機構と温泉起源流体酒井拓哉 1）#*，大沢信二2），山田誠3），三島壮智2），大上和敏1）（平成 25 年 2 月 12 日受付，平成 25 年 8 月 28 日受理）

A Formation Mechanism and the Origin of Yamaga
Hot Spring in Oita Prefecture Inferred from
Geochemical Data of Hot Spring Waters
and Their Accompanied Gases
Takuya Sakai1）#*, Shinji Ohsawa2）, Makoto Yamada3）
Taketoshi Mishima2） and Kazutoshi Oue1）
Abstract
In order to discuss the formation mechanism and origin of source fluid of Yamaga hot
spring in Kitsuki city, Oita Prefecture, Japan, we sampled hot spring waters and their accompanied gases, and conducted their chemical and isotopic analyses. Geochemical analysesof obtained data revealed that Yamaga hot spring is formed by the mixing between high salinity volcanic fluid of Na-Cl type released from a deep-seated magma and ordinary shallow
ground water of HCO3 type. Moreover, it was suggested that the magma participating in
the formation of Yamaga hot spring is appropriate to be an ancient magma and therefore
the volcanic fluids released from the magma have been reserved for a long period of time
and would be considerably cooled.
Key words : Ashio Mountains, Akagi Volcano, Lowland of the Middle Part of the River Tonegawa, Hot
spring waters, Fluid formation mechanism, Geological and mineralogical interpretations
大分大学大学院教育学研究科教科教育専攻理科専修 〒870-1192 大分県大分市大字旦野原 700．1）Graduate
School of Education, Oita University, 700 Danoharu, Oita, Oita Prefecture 870-1192, Japan. # 現在，日鉄鉱業株式会社資源開発部地熱開発課 〒100-8377 東京都千代田区丸の内二丁目 3 番 2 号（郵船ビル
6F）．# Geothermal Development Section, Resources Development Dept., Nittetsu Mining Co., Ltd.,
Yusen Building, 3-2, Marunouchi 2-Chome, Chiyoda-Ku, Tokyo 100-8377, Japan. *Corresponding
author：E-mail t-sakai@nittetsukou.co.jp, TEL 03-3216-5281, FAX 03-5962-5180.
2）
京都大学大学院理学研究科附属地球熱学研究施設 〒874-0903 大分県別府市野口原．2）Beppu Geothermal
Research Laboratory, Institute for Geothermal Sciences, Graduate School of Science, Kyoto University,
Noguchibaru, Beppu, Oita Prefecture 874-0903, Japan.
3）
奈良女子大学共生科学研究センター 〒630-8506 奈良市北魚屋東町．3）KYOUSEI Science Center, Nara
Women’s University, Kitauoyahigashimachi, Nara, Nara Prefecture 630-8285, Japan.
1）

164

第 63 巻（2013）

温泉水・温泉付随ガスの地球化学データから見た大分県山香温泉の生成機構と温泉起源流体

要旨大分県杵築市山香町のガスを伴って湧出する高塩分の温泉（山香温泉センター，山香温泉風の郷）において温泉水および温泉付随ガスを採取し，それらの化学・同位体分析を行い，分析
結果の地球化学的解析から山香温泉の生成機構と温泉起源流体について考察した．その結果，山香町の温泉は地下深部でマグマから放出された Na と Cl に富む高塩分の火山性流体に起源をもち，浅所の低塩分 HCO3 型地下水の混入を受けて生成していることが示された．また，その生成に関与したマグマは現在のものではなく，古い時代のそれとするのが妥当であり，起源流体と考えられた火山性流体は地下に相当に長期間滞留し，冷却したものであることが示唆された．
キーワード：山香温泉，高塩分水，第三紀，マグマ活動，温泉付随ガス

1． は じ め に
近年，非火山地域に湧出する高塩分の温・鉱泉水を対象とした研究が数多く行われており（例えば，村松ら，2010；大沢ら，2010），大分県でもいくつかの温泉や鉱泉において地球化学的な研究が行われている（大沢，2003，2005；網田ら，2005）．そのような中，杵築市山香町にもガスを伴いながら自噴している高塩分の鉱泉の存在が古くから知られていたが，本格的な研究は行われておらず，山下（1982）による簡単な調査報告があるのみであった．そこで，本研究では山香温泉の生成機構を明らかにするために当該地域の源泉を改めて調査し，泉源の稼動状況など現況把握を行うとともに，温泉水と，今回初めて温泉付随ガスを採取した．温泉水については一般水質，水の安定同位体組成（δD と δ18O），溶存全炭酸（DIC）の濃度とその炭素安定同位体比（δ13C）のデータを，
温泉付随ガスについては二酸化炭素，ヘリウム，水素，アルゴン，窒素，メタン，酸素の化学組成
データとともに希ガス同位体（3He, 4He, 20Ne）濃度，CH4 に対するエタン，プロパンの相対組成とそれら炭化水素ガスの炭素安定同位体比（δ13C）を測定し，得られた地球化学データに基づいて温泉の生成機構について考察したので報告する．

2． 調査泉源と研究対象地域の地質調査の際に利用されているのが確認できた泉源は，山香温泉センターと山香温泉風の郷の 2 泉源のみであり（Fig. 1），いずれも多量のガスを伴いながら自噴していた．この温泉付随ガスは，一般に言われているように，元々地中では温泉水に溶存していたものが井戸管内を伝って上昇してくる
間に温泉水から脱ガス・分離したと考えられる．
山香温泉センターは，かつて中山香駅から約 200 m 下流の立石川の川岸に湧出していた神塩（こうじょう）鉱泉を医療用に活用する目的で昭和 55 年に立石川の道路脇に深度 700 m のボーリングによって得られた温泉である（山下，1982）．一方，山香温泉風の郷は，平成 12 年に 1,500 m の深度ボーリングで新規開発された温泉である．泉源位置周辺の詳しい地質（産業技術総合研究所地質調査総合センター，2011；山下，1982）をみると，山香温泉センターの泉源付近の川底には，新第三紀に多くの鉱床を産んだ熱水溶液による変質をうけた岩石であるプロピライトが露出している．
また，その泉源位置は後期更新世─完新世の堆積岩類の分布域内にあるが，ボーリングで得られた岩石についての研究によれば，地下地質は安山岩らしい原岩が熱水変質を受けており，亀裂が多く，各所に破砕帯があり，粘土鉱物や熱水性鉱物で充填されていることが分かっている（山下，1982）．山香温泉風の郷の地下地質の情報は入手できなかったが，後期中新世─鮮新世の火山岩類の分布域にあることから，山香温泉センターと同様に熱水変質を受けた安山岩からなるのではない165

酒井拓哉，大沢信二，山田誠，三島壮智，大上和敏

温泉科学Fig. 1 Maps showing localities of hot springs investigated in this study (○) with related
geologic information.

かと推測される．
以上の他，温泉関連の地学現象として，本地域には，温泉地域の北側に隣接する辺りに，浅熱水性の金鉱床である馬上（ばじょう）鉱床が存在することが特筆される（星住・森下，1993；Watanabe,2005）．鉱床の分布域には上述のようなプロピライト化作用を受け熱水変質した火山岩類が広い範囲に存在しており，熱水変質と金の鉱化作用は一連の熱水活動によって生じたと考えられている（星住・森下，1993）．

3． 試料採取方法と化学分析
3．1

試料採取

山香温泉センター，山香温泉風の郷の 2 泉源において温泉水試料および温泉付随ガス試料を採取した．試料は後述のように分析目的別に異なる容器に採取し，実験室に持ち帰った．この他に現地では，温泉水の水温と pH を測定した．
全溶存炭酸（DIC＝CO2（aq）＋HCO3－＋CO32－：Dissolved Inorganic Carbon）の化学・同位体分析用の試料水は，容器壁からの CO2 逃散を回避するために，ガスバリア性に優れたニッコー製テクノボトル（BAREXⓇ）に採取した．水の水素・酸素同位体比（δD と δ18O）分析用には，保存中に水の蒸発が起こらないように，ガラスバイアルに採取した．その他の成分の分析用試料水は，一括して汎用的なポリ瓶に採取し，分析まで冷暗所に保存した．
化学分析用のガス試料は，水上置換法によって採取した．気泡状のガスを漏斗を用いて一旦ポリ袋に集めた後，現地で直ちにアルカリ溶液（5M KOH）を含むガラス製浣腸器に移した．その際，ガス試料は繰り返し吸引することによりアルカリ溶液に吸収されないガス成分（残留ガス：以下 R
ガスとする）のボリュームを十分に確保し，実験室に持ち帰った．なお，アルカリ溶液を含むガラス浣腸器の準備は Ohsawa et al.（2000），網田・大沢，（2001）に準じて行った．
ヘリウム同位体比測定用のガス試料も上記の化学分析用のガス試料と同様に水上置換法によって166

第 63 巻（2013）

温泉水・温泉付随ガスの地球化学データから見た大分県山香温泉の生成機構と温泉起源流体採取した．気泡状のガスを漏斗を用いて両端にストップコックの付いた鉛ガラス製採気瓶に採取した．
炭化水素に対するエタン（C2H6），プロパン（C3H8）の相対組成とそれらの炭化水素ガスの炭素安定同位体比（δ13C）分析用試料は，漏斗を用いて気泡状のガスを直接 10 mL または 20 mL のゴム製セプタム付きガラスバイアルに水上置換法によって少量の温泉水と共に採取した．運搬・保存時にはガラス瓶を転倒させ，セプタムとの合わせ部分が温泉水でシールされるようにして外気と遮断し，試料ガスの漏出と空気汚染を防いだ．
3．2

化学分析

試料水の分析では，HCO3－を除く主要化学成分（Na, K, Ca, Mg, Cl, SO4）と Li はイオンクロマ
トグラフィー（DIONEX 社製 DX-120）で，B は分光光度法（アゾメチン H 法，SHIMADZU 社製
UVmini-1240）で行った．DIC 濃度は，試料水を硫酸酸性にして試料中の炭酸成分（CO2（aq）, HCO3－ ,CO32－）をすべて CO2（aq）に変換し，イオンメーター（TOA DKK 社製 MM-60R）に接続した二酸化炭素電極（TOA DKK 社製 CE-2041）を用いて絶対検量線法により測定した．HCO3－ 濃度は DIC
濃度の測定値と水温，pH の値を用いて化学平衡計算から求めた．
水の水素・酸素の同位体比（δD と δ18O）の測定は，質量分析計（GV Instrument 社製 IsoPrime）
を用いてそれぞれ金属亜鉛を用いた還元法，二酸化炭素平衡法により行った．測定精度は δD, δ18Oのそれぞれについて±1.0‰，±0.1‰である．
化学分析用のガス試料は，アルカリ吸収液は H2O2 で酸化処理した後定容とし，分取した検液を硫酸酸性にして，東亜電波製 CE-235 型ガス電極を用いて CO2 を分析した．R ガスは容積を測定後（この時に気温・気圧を測定），一定量を O2 をキャリアーガスとした大倉理研製 AFG-555S 型ガス
クロマトグラフに注入することにより，He, H2, Ar, N2 および CH4 を分析した．O2 は R ガスの体積と試料ガス注入時の気温・気圧，各成分の分析値を基に理想気体の状態方程式を用いて算出した．
DIC および CO2 ガスの炭素安定同位体比（δ13C）は，試料水および CO2 ガスを吸収したアルカリ
溶液に Sr（OH）
2 を加え加熱して沈殿した SrCO3 をろ別し風乾した後，これに濃リン酸を滴下して発生する CO2 を質量分析計（GV Instrument 社製 IsoPrime）に導入し測定した．測定精度は±0.2‰である．
希ガス同位体比（3He/4He, 4He/20Ne）の測定は，磁場型質量分析計（東京大学大学院理学系研究科附属地殻化学実験施設の VG Isotope 製 VG5400）により行った．測定精度は 3He/4He, 4He/20Neのそれぞれについて 1%，10% である．
炭化水素に対するエタン（C2H6），プロパン（C3H8）の相対組成の化学分析は JIS-2301「燃料ガス及び天然ガスの試験分析法」に基づきガスクロマトグラフ（メーカー機器名）を用いて行い，炭素
安定同位体分析はガスクロマトグラフ燃焼前処理装置付安定同位体質量分析計（GV Instrument
社製 IsoPrime-GC）で測定した．測定精度は±0.2‰である．

4． 結
4．1

果

温泉水の溶存化学成分

温泉水の化学・同位体分析結果を，水温，pH の測定値とともに Table 1 に示す．主要溶存化学
成分について分析値をみると，山香温泉センター（以後，温泉センターと省略する），山香温泉風の郷（以後，風の郷と省略する）ともに陽イオンについて卓越する成分はナトリウムイオンであり，陰イオンについては温泉センターが塩化物イオン，風の郷が炭酸水素イオンであった．そして，得167

酒井拓哉，大沢信二，山田誠，三島壮智，大上和敏

温泉科学

Table 1 Chemical and isotopic data of sampled hot spring waters
with information on depth of well and water temperature.
Point Name

Onsen center

Kazenosato

Sampring Datedepth of well (m)29/9/2010700

29/9/2010
1500

WT (℃)
pH

32.3 (in the tank)
6.7

39.7
6.9

Li+ (mg/L)
Na+ (mg/L)
K+ (mg/L)
Mg2+ (mg/L)
Ca2+ (mg/L)
Cl－ (mg/L)
SO42－ (mg/L)
HCO3－ (mg/L)
B (mg/L)DIC (mg/L)
δD of H2O (‰)δ18O of H2O (‰)δ13C of DIC (‰)

35.4
9520
716
143
284
12500
0.9
5240
62
4800
－44.4
－1.65
－2.9

17.3
3870
298
210
154
3320
0.7
5290
25
4750
－55.5
－6.74
－2.66

られた分析データを山香町の湧水，河川水，井戸水のデータ（山下，1982）とともに陰イオンの組成図に表すと（Fig. 2a），2 つの温泉水は塩化物イオンと炭酸水素イオンを結ぶ直線上に位置し，酸性鉱山排水を除く山香町の湧水，河川水，井戸水（通常の地下水）は炭酸水素イオンコーナーに偏って分布している．このことから，山香町では温泉センターに湧出するような塩分濃度の高い Na とCl に富む水に塩分濃度の低い HCO3 型水質の地下水が混合して温泉水が生成していることがうかがわれる．
その混合関係は Giggenbach（1988）により提案されている陽イオン組成図（Fig. 2b）にもよく一致する．この解析図は地下の熱水系の貯留状態および貯留層温度を推定する方法として水─岩石
相互作用を仮定した化学平衡により Na-K 地化学温度計と K-Mg 地化学温度計を組み合わせた図であるが，その基礎となる水─岩石相互作用に関与する鉱物は，地熱系に普遍的に存在する曹長石，
カリ長石，緑泥石，カリ雲母，シリカである．これらと熱水の間で以下の平衡関係（1）および（2）の両方が成り立てば完全平衡（Full equilibrium）を示す線上にプロットされる．そして，酸性の水など岩石から熱水へ一方的な溶解を示すものについては未成熟水（immature waters）の領域にプロットされる．
（1）albite＋K+ ⇄ K-feldspar＋Na+（2）0.8K-mica＋0.2chlorite＋5.4silica＋2K+ ⇄ 2.8K-feldspar＋1.6H2O＋Mg2+
また，完全平衡にあった水が低温の地下水の混入を受けるなどして温度低下を経験すると，平衡の反応速度が（1）≪（2）であるため，部分平衡（partially equilibrated waters）領域にプロットされる．なお，この領域内に引かれている線はその水が部分平衡に至る以前に経験した完全平衡状態にあるときの温度を示すもので，プロットされた位置からその温度が読み取れる．
本研究のデータをみると，温泉センターは部分平衡領域にプロットされ，風の郷は未成熟水（immature waters）の領域にプロットされているが，両者を結ぶ直線（図中に破線で示した．）の168

第 63 巻（2013）

温泉水・温泉付随ガスの地球化学データから見た大分県山香温泉の生成機構と温泉起源流体

Fig. 2 Ternary plots for classifying sampled waters based on the relative
proportions of major anions (a) and for recognizing the equilibrium state
and/or a mixing relation on the sample waters (b) with data of othernatural waters in Yamaga reported by Yamashita (1982).

169酒井拓哉，大沢信二，山田誠，三島壮智，大上和敏

温泉科学

一方は，この地域の通常の地下水（湧水，河川水，井戸水）がプロットされる Mg コーナーに達しており，この直線は混合関係を示していると考えられ，主要陰イオン組成から読み取られた内容と整合的である．すなわち，温泉センターのような Na と Cl に富む高塩分水に普通の地下水が混ざって温泉水が生成しているということであり，もし温泉水の起源となる高塩分水が元は上記平衡関係にあって完全平衡状態であったと仮定すると，その温度は 210℃程度と相当に高温であったことになる．
その他の温泉水の化学的特徴として，Li が比較的多量に含まれていることが挙げられる．特に温泉センターの Li 濃度は大分県の温泉の中でも最高値に近い値を示しており，大分県玖珠郡九重町の大岳・八丁原地熱発電所から得られる熱水の Li 濃度（最大 19 mg/L；阿部・茂野，1986）と比べてもどちらの温泉も高い値となっている．また，最近では西尾（2010）により，低温しか経験していない表層水に比べて，高温を経験した深部流体は有意に Li に富んでいることが示されており，これらのことから山香町の温泉水は低い泉温にも関わらず高温の深部熱水としての一面を見せている．
Li は B, Cl と同様に可溶性元素であり，水─岩石相互作用において流体へ移動した後に 2 次的な変化を受けにくく，流体相にとどまるとされているため，温泉水の起源推定などに用いられる．そこで，温泉センター，風の郷のデータを Li-B-Cl の相対組成図上で，火山性深部熱水を起源水にもつ別府の温泉水（NEDO, 1989），山香温泉と同様な高塩分温泉水を賦存させる非火山性の宮崎平野の温泉（大沢ら，2010）や千葉のガス田付随塩水（相川，1991）のデータと比較した（Fig. 3）．山Fig. 3 Relative B-Li-Cl compositions of sampled hot spring waters with data points ofmodern seawater, hot spring waters in Beppu (NEDO, 1989), brine waters accompaniedwith natural gases in Chiba (Aikawa, 1991) and hot spring waters in Miyazaki originatedfrom squeezed pore waters from marine sediments (Ohsawa et al., 2010).

170

第 63 巻（2013）

温泉水・温泉付随ガスの地球化学データから見た大分県山香温泉の生成機構と温泉起源流体

香の温泉水は，海水や古海水を直近の起源とするガス田付随塩水や海底堆積物間隙から排出された海水（大沢ら，2010）とは明確に区別され，別府の熱水性温泉水のデータに近い範囲にプロットされることがわかる．このことから，山香温泉の深部起源水として考えられる Na・Cl に富む高塩分
水は古い時代に地層中に閉じ込められた海水に由来する変質海水ではなく，火山性流体との関連性が示唆される．
4．2

温泉付随ガスの成分

温泉付随ガスの化学・同位体分析結果を Table 2 に示す．ガス試料中の O2 濃度が，温泉センター
が 0.23%，風の郷が 0.65% とわずかな量であり，採取・輸送時の試料の空気汚染はほとんどなく，
良好なサンプリングが行えたことを示している．
最も卓越する付随ガス成分は温泉センター，風の郷ともに二酸化炭素であり，次いで，メタン，
窒素でこれら 3 成分で全体の 99.9％以上を占めていることがわかる（Table 2）．この主要 3 成分の組成を，地理的に最も山香町に近いところにある活火山（鶴見火山）に位置する別府温泉や九重火山（硫黄山）の噴気ガスのデータ（大沢，2000；網田・大沢，2003）とともに見てみると（Fig. 4），
新しい時代の火山活動にともなうものとはやや違う組成を示し（CO2 濃度を 1/10 で評価していることに留意），相対的に CH4 が比較的多くなっていることがわかった．
主要 3 成分のうちの窒素（N2）は大気を構成する主要成分でもあり，同様な成分として温泉付随
ガス中にはアルゴン（Ar）が含まれている．これら 2 つの成分は化学的に不活性であることから，
Ar の全てが大気起源ならば（Giggenbach and Goguel, 1989），N2/Ar 比が大気のそれ（83.6）より大きい場合に大気に由来しない N2 の混入が予想される．温泉センター，風の郷のそれはそれぞれ
112,96 であり，ともに大気によらない部分がガス試料に含まれていることを示している．そこで，深部起源ガスの混入の検討に使われる N2-He-Ar 相対組成図上（Giggenbach and Goguel, 1989）に温泉センター，風の郷のデータポイントを別府温泉の噴気ガスのデータ（大沢，2000），ならびに別府に噴出するマグマ起源ガスに最も近いとされる噴気ガスのデータ（Kita et al., 1993），沈み込

Table 2 Chemical and isotopic data of sampled gases associated by
hot spring waters.
Point Name

Onsen center

Kazenosato

CO2 (%)
N2 (%)
CH4 (%)
Ar (%)He (ppm)
H2 (ppm)
C2H6 (ppm)
C3H8 (ppm)
δ13C of CO2 (‰ )
δ13C of CH4 (‰ )
δ13C of C2H6 (‰ )
δ13C of C3H8 (‰ )3
He/4He4
He/20Ne

91.1
2.47
6.41
0.022
64
0.9
98
3.6
－7.04
－45.6
－43.3
－36.9
7.13×10－6
698

88.8
5.46
5.6
0.057
41
29
57
1.1
－6.5
－40.7
－46.4
n.a.
7.18×10－6
1180※n.a. : not analysis

171酒井拓哉，大沢信二，山田誠，三島壮智，大上和敏

温泉科学Fig. 4 Relative CO2-N2-CH4 contents of sampled gases associated with
hot spring waters with data points of fumaloric gases from Beppu(Ohsawa, 2000) and from Kuju-Iwoyama (Amita and Ohsawa, 2003).

みプレートあるいは堆積物由来のガス，マントルあるいは地殻起源ガス（Goff and Janik, 2002）とともに表した（Fig. 5）．それをみると，マグマ性である別府地域の噴気ガスの組成範囲内あるいはそれに近いところにプロットされ，大気によらない部分が含まれることが明示されている．また，
温泉センターと風の郷のデータポイントを直線で結ぶと，その一方は空気に至ることから，この直線（図中に破線矢印で示してある．）は大気とマグマ性ガスあるいはマントルや地殻といった地球内部に由来するガスが深部起源ガスとの混合関係を表していると妥当である．そして，もしこの端成分をなす深部ガスが，温泉水からも示唆されたような火山性であるとすると，図中に示されるよ
うにそのマグマ起源ガスは別府のものに比べてやや He に富むという性質をもつことが暗示される．これは，後の付随ガスの希ガス同位体組成の解析によって妥当であることが示される．

5． 考
5．1

察

温泉水の同位体組成

前項（4．1）において温泉起源流体の水質や温泉水形成の際の混合関係が示され，深部温泉起源
流体の候補としてマグマ起源の火山性流体の存在が示唆された．本項では，その妥当性や温泉起源
流体の正体を知るために，まず，温泉水の同位体組成の考察を行う．
172

第 63 巻（2013）

温泉水・温泉付随ガスの地球化学データから見た大分県山香温泉の生成機構と温泉起源流体

Fig. 5 Relative N2-He-Ar contents of sampled gases associated with hot spring waters
with data points of magmatic fumarolic gases in Beppu (Kita et al., 1993).  ASW means
air-saturated water.   End members for slab, sediments, and mantle, crustal are citedfrom Goff and Janik (2002).  For details of magmatic gas in Yamaga, see text.

試料水の安定同位体データ（δD・δ18O）を，後述の参照データとともに δD vs. δ18O 図上に示した（Fig. 6）．参照したデータは，九重硫黄山から放出される噴気水蒸気の同位体データ（Mizutani etal., 1986），薩摩硫黄島火山から放出される噴気水蒸気の同位体データ（Goff and McMurty, 2000），
雲仙火山から放出される噴気水蒸気の同位体データ（Ohba et al., 2008），海水や古海水を起源とす
る宮崎平野の温泉のデータ（大沢ら，2010），千葉のガス田付随塩水のデータ（Muramatsu et al.,2001）の他，天水線（MWL：meteoric water line；δD＝8δ18O＋10），現海水（modern seawater）の値，安山岩質マグマ性水蒸気（andesitic magmatic water）がとる同位体比の範囲（Giggenbach,1992）を図上に表した．山香温泉のデータポイントを通る直線は，安山岩質マグマ性水蒸気のとる範囲へ向かい，天水起源の地下水と安山岩質マグマ性水蒸気の混合線であることを示しており，山香町の高塩分温泉水の起源流体は海水やそれに直接由来する高塩分流体ではなく，マグマから放出される火山性流体であるという見解を支持していると考える．
ところが，δD vs. δ18O 図上に表したデータについて，水素安定同位体比（δD）と塩化物イオン（Cl－）の関係（Fig. 7）を見てみると，液性が中性である山香の温泉水は，安山岩質の九重硫黄山，薩摩硫黄山，雲仙火山から放出される HCl, SO2 といった酸性ガスを含む高温噴気ガスとは異なるトレンドにあり，より Cl－に富むという異なった性質を有することがわかる．前項で山香町の温泉水は深部高塩分起源水に浅所の通常的な地下水が混合していることが示されたが，その高塩分水の塩化物イオン濃度は，温泉水のそれよりもさらに濃いことが考えられ，その場合その濃度は Fig. 7 の山
香温泉の 2 つのデータを結ぶ直線を外挿した先にあるものと予想される．山香の温泉起源流体のも173酒井拓哉，大沢信二，山田誠，三島壮智，大上和敏

温泉科学Fig. 6 δD vs. δ18O plot of sampled hot spring waters with range of andesitic magmaticsteam and data of modern seawater, suqueezed pore waters from marine sediments(Ohsawa et al., 2010), brine waters associated by natural gases in Chiba (Muramatsu
et al., 2001) and fumarolic gases from Kuju-iwoyama, Satsuma-iwojima and Unzan(Mizutani et al., 1986 ; Goff and McMurty, 2000 ; Ohba et al., 2008).

つ水素安定同位体比（δD）が安山岩質マグマ性水蒸気のとる同位体比と同様であるとし，上記外挿直線から温泉起源流体のもつ塩化物イオン濃度を予想すると約 33,000 mg/L となる．Figure 7 中に示されている外挿直線の先端に示してある星印（☆）はこのようにして求めたものであり，山香の温泉起源流体は海水やそれに直接由来する高塩分流体とは明らかに異なることが確認できる．
5．2

希ガス同位体データ

前項（4．2）でマグマ由来の He が温泉付随ガスに含まれている可能性が示された．そこで，Heの起源についてさらに言及するために，山香温泉の付随ガスのデータを 3He/4He－4He/20Ne ダイア
グラム上に，別府の噴気ガスのデータ（大沢，2000），西南日本で最も高い値を示した別府のマグ
マ性噴気ガスのデータ（北，1995）とともに表した（Fig. 8）．20Ne はマグマ性気体にはほとんど含
まれず，試料中の 4He/20Ne 比が大きいほどマグマ性 He の混入率が高いことを示す．山香の温泉のデータを見ると，別府のマグマ性噴気ガスと同様な値をとっており，温泉付随ガスへのマントル成分に富んだマグマ性 He の混入が明瞭に示される．また，現在のマグマ活動の産物である別府地域の噴気 He に比べ 3He/4He 比がやや低いという特徴が見られ，地殻起源ガスの混入が考えられる．
これは，ヘリウム同位体（3He，4He）のうち，3He は地球生成時に地球内部のマントルに取り込まれた始源的な成分であるが，4He は地殻岩石中のウランやトリウムの放射壊変に伴う α 粒子起源の174

第 63 巻（2013）

温泉水・温泉付随ガスの地球化学データから見た大分県山香温泉の生成機構と温泉起源流体

Fig. 7 Relationships between δD and concentration of Cl－ of sampled hot spring waters withdata of modern seawater, suqueezed pore waters from marine sediments (Ohsawa et al.,2010), brine waters accompanied with natural gases in Chiba (Muramatsu et al., 2001)
and fumarolic gases from Kuju-iwoyama, Satsuma-iwojima and Unzan (Mizutani et al.,1986 ; Goff and McMurty, 2000 ; Ohba et al., 2008).  For source fluid in Yamaga, see text.

成分であり，図中に指示したように，3He/4He 比が低いほど地殻の放射性起源 He が付加されたものと解釈されるからである．
地殻起源の He の付加を視覚的だけでなく数値的に示すために，Shimizu et al.（2005）に準じ，下式を用いて起源ヘリウムの混合率を推算してみた．
（3He/4He）（3He/4He）（3He/4He）（3He/4He）
S＝
M・M＋
C・C＋
A・A

［1］

1/（4He/20Ne）（4He/20Ne）（4He/20Ne）（4He/20Ne）
S＝1/
M・M＋1/
C・C＋1/A・A

［2］

M＋C＋A＝1

［3］

ここで，S, M, C, A はそれぞれ試料，マントル，地殻，大気を表わし，大文字はそれらの混合率を，
添字をもってそれぞれの 3He/4He 比と 4He/20Ne 比を表わす．それぞれのパラメータに，例えば，
Shimizu et al.（2005）に表してある値を代入して，起源ヘリウムの混合率 M, C, A を計算した．その結果，温泉センターはマントル起源ヘリウムが 64%，地殻起源ヘリウムが 36%，大気起源ヘリウ
ムが 0%，一方，風の郷はマントル起源ヘリウムが 65％で，地殻起源ヘリウムが 35%，大気起源ヘリウムは 0% となった．さらに，大沢（2000）から別府の噴気ガスデータを用いて同様な計算を行ったところ，マントル起源ヘリウムは 69％で，地殻起源ヘリウムは 24%，大気起源ヘリウムが 7%となり，温泉センター，風の郷の He に占める地殻起源/マントル起源比はそれぞれ 0.56，0.54 で，別府のそれは 0.35 であり，山香温泉の He はともに別府地域に比べ地殻起源のヘリウムの寄与が相対的に大きいことが客観的に示された．また，このことは，前節の N2-He-Ar 相対組成図（Fig. 5）
でみた別府の噴気ガスに比べて山香の温泉付随ガスがやや He に富む傾向に合点のいく説明を与175酒井拓哉，大沢信二，山田誠，三島壮智，大上和敏

温泉科学Fig. 8 3He/4He vs. 4He/20Ne plot for associated gases of hot spring waters sampled in thisstudy with fumarolic gases and magmatic fumarolic gases in Beppu (Ohsawa, 2000 ; Kita
et al., 1993).   The upper line represents a mixing line between air and upper mantle
components ; the lower one is a mixing line for air and crust components.  For magmaticgases in Yamaga, see text.

え，それは地殻起源の 4He の付加によるものであるとしてよい．
ここまで見てきた温泉水および温泉付随ガスの地球化学的諸情報は，一致して，山香町の温泉の形成にマグマからの流体が関与したことを前向きに支持しているが，山香町には第四紀の火山は存在せず，広範囲にわたって第三紀の火山岩類の分布が見られるのみである（Fig. 1）．このことから，
山香の深部温泉起源流体は古い時代のマグマに由来するとするのが順当であり，そうであればマグ
マ性流体が長い間地中に滞留することで放射壊変起源の 4He が相当分に付加されても別段不思議ではない．次項では，深部温泉起源流体が長時間滞留していたことを示すと考えられる別の証拠を
提示する．
5．3

炭化水素ガスの化学・同位体データ

山香温泉の温泉付随ガス中には無視できない量の CH4 が存在することは前に述べた（Table 2）．
CH4 の起源を推定するために用いられる炭化水素の化学・同位体データを有機地球化学分野で汎用される Bernard 図に表す（Fig. 9）．早稲田・岩野（2007）によると地層を長距離・長時間ガスが移動した場合，通常，炭素同位体組成は一定のまま，エタン以上の炭化水素が選択的に地層に吸着され，C1/（C2＋C3）比が上昇するとされている．本調査地点がプロットされる位置をみると付随ガス
中の CH4 は有機物の熱分解によって生成したもので，生成場からやや浅所へ上昇してきたことが示される．
通常，高温のマグマと共存する熱水流体では CH4 はほとんど生成せず（Giggenbach, 1997），山香の温泉付随ガスは別府の噴気ガスに比べて CH4/CO2 比は 1～2 桁高く，熱分解起源メタンは温泉
176

第 63 巻（2013）

温泉水・温泉付随ガスの地球化学データから見た大分県山香温泉の生成機構と温泉起源流体

Fig. 9 Relation between δ13C of CH4 and CH4/(C2H6+C3H8) ratio of associated gases ofhot spring waters sampled in this study expressed on the Bernard plot, which is often
used for discrimination of origin of CH4 and interpretation of physical and chemical
processes on hydrocarbon gases.

起源流体となったマグマ性流体が貯留層中で低温化した際に付加されたものと考えるのが妥当である．また，Allis and Yusa（1989）や大沢（2000）によると，別府温泉の本源熱水の温度は 250 ─ 300℃であるとされ，また，網田・大沢（2003）で実測されている九重硫黄山の噴気ガスの最高温度は308℃であることなどから，現在活動中の火山の地下の熱水系の温度は 300℃近くに達するのが普通であると思われる．そのような新しい火山の活動の影響を受けている熱水と比べると，山香町の温泉を形成する深部起源流体が完全平衡を経験した場合でも高々210℃程度と低かったのは，永い貯留による低温化のためであると見ることもでき都合が良い．
5．4

炭酸成分の δ13C と ΣC（CO2）/3He の関係

ここでは，温泉の主要成分である炭酸成分（溶存全炭酸＋付随 CO2）の由来について考察する．
炭酸成分の起源を探る場合，炭素安定同位体比（δ13C）を用いることが有効であり，マントルのδ13C 値（－10‰～－4‰），海成炭酸塩の δ13C 値（約 0‰），堆積性有機物の δ13C 値（約－25‰）が指標となるが，地殻物質のうち堆積性有機物起源炭素と海成炭酸塩起源炭素の適当な割合での混合によりマントルの δ13C 値と区別できなくなる．そのため，それを補うための炭素同位体比（δ13C）以外の情報が必要であり，その問題を解決するために提案されたものが 3He を併用した δ13C と ΣC/3Heの関係から炭素の起源を推定する方法（Sano and Marty, 1995）である．第 2 節で述べたように，
温泉付随ガスは，湧出してくる途中で脱ガス・分離したと考えられるため，ガス相の CO2 と水相の溶存全炭酸（DIC）の間の同位体分別を補正する必要がある．そこで，下式を用いて気液分離以前の δ13C 値を求めた．
δ13Ccor＝δ13Cco2×Xco2＋δ13CDIC×XDIC

［1］

Xco2＋XDIC＝1

［2］
177

酒井拓哉，大沢信二，山田誠，三島壮智，大上和敏

温泉科学

ここで δ13Ccor は補正後の δ13C 値，δ13Cco2，δ13CDIC はそれぞれガス相の CO2 および水相の炭酸成分（DIC）の δ13C 値，Xco2，XDIC はそれぞれガス相の CO2 および水相の炭酸成分（DIC）のモル分率を示す．このモル分率はガス/水比とガス相および水相の CO2 濃度，水相の DIC 濃度から求められる．今回は，温泉井および輸送管の構造の制約からガス/水比の測定が困難であったため，長湯温泉の実測値（岩倉ら，2000）を参照し，最小値 0.57，最大値 3.74，中間値 2.23 を想定して計算に用いた．計算結果を Table 3 に，そして，それら補正値をプロットした δ13C-ΣC/3He 図を Fig. 10に示した（Fig. 10 のプロット A はガス/水比を最小値 0.57 で計算，プロット B はガス/水比を中間値 2.23 で計算，プロット C はガス/水比を最大値 3.74 で計算した場合の結果を示している．）この図には，別府の噴気ガスのデータ（大沢，2000）も同時に表してあるが，プレートの沈み込み運動
によりマントル内に持ち込まれた海成炭酸塩（CaCO3）が分解して CO2 となり，火山・地熱活動
を通して再び地表に噴出するとされている別府地域の噴気ガスに近い範囲にプロットされることが示された．
Sano and Marty（1995）の解法は質量および同位体の保存則を前提とした混合関係の解析であるので，起源炭素の δ13C 値を与えられれば，下式によって起源炭素の混合率の推算を行うことができる．
δ13CS＝δ13CM・M＋δ13CC・C＋δ13CO・O

［1］

1/（CO2/ He）（CO2/ He）
（CO2/ He）
（CO2/ He）
S＝1/
M・M＋1/
C・C＋1/O・O

［2］

M＋C＋O＝1

［3］

3

3

3

3

ここに，S, M, C, O はそれぞれ試料，マントル，沈み込み海成炭酸塩，堆積性有機物を表わし，大文字はそれらの混合率を，添字をもってそれぞれの δ13C 値と CO2/3He 比を表わす．それぞれのパラメータに，Sano and Marty（1995），Nishio et al.（1998），Deines（2002）に表してあるマントル
起源炭素，海成炭酸塩，堆積性有機物の値の平均値を代入して，起源炭素の混合率 M, C, O を計算した．その結果を Table 3 に示す．ガス/水比が極端に低い 0.57 という最小値の場合を除き，温泉センター，風の郷いずれもマントル CO2＝10%，堆積性有機物＝20%，沈み込み海成炭酸塩＝70％

Table 3 Summary of corrected δ13C and C/3He values of hot spring waters before
gas separation and estimated contribution ratios of their source carbon.
gas-water ratio
0.57
Onsen center

δ C
13

※1
cor

CO ,

/ He

※2
2 cor

Kazenosato

3δ13Ccor※1

－4.05

/ He
3

MantleOrganic
Carbonate

178

1.1×10
63
5
32

CO ,

※ 1：corrected δ13C
※ 2：corrected CO2－4.4

MantleOrganic
Carbonate

※2
2 cor

2.23

1.7×10
42
7
50

3.74

－5.76
9

－6.17

4.4×10

7.4×109

16
18
65

10
21
69

9

－5.31
9

－5.69

6.6×101.1×1010

11
18
71

6
20
74

9

第 63 巻（2013）

温泉水・温泉付随ガスの地球化学データから見た大分県山香温泉の生成機構と温泉起源流体

Fig. 10 Correlation between ΣC/3He and δ13C of Yamaga hot spring waters before
separating gas phase.  Data points A, B and C for respective hot spring waters varyin gas/water ratios.  For the assumed gas/water ratio, see text and Table 3.  Dataof fumarolic gases in Beppu are cited from Ohsawa, (2000).  Model end members for
mantle (M-a, Sano and Marty, 1995 ; M-b, Denis, 2002 ; M-c, Nishio et al., 1998),
marine carbonate (Sano and Marty, 1995) and sedimentary organic carbon〔meanvalue of Sano and Marty (1995), Denis(2002) and Nishio et al. (1998)〕are shownin this figure.  The tie lines of the model end members are mixing lines.

と沈み込み帯の火山から放出される噴気ガスと同様に沈み込み海成炭酸塩起源の炭素の寄与が相対的に大きいことが数値的にも示された．
鎌田（2002）によると，九州地方を含む西南日本では約 12 Ma にフィリピン海プレートの沈み込みが始まり，現在まで継続していたとされている．山香町も位置する豊肥火山地域では島弧性の火成活動が開始した約 6 Ma から現在に至るまで活動が連続し（現在は九重や別府でみられる），
これらの火山岩が海洋プレートの沈み込みを伴う島弧に典型的な化学組成を示すことがわかっている．この地質学的な情報を加味すると，山香温泉の炭酸成分（溶存全炭酸＋付随 CO2）は約 6 Ma 以
降のフィリピン海プレートの沈み込みによってマントル内に持ち込まれた海成炭酸塩が分解され，マグマ活動の産物として地表付近にもたらされたということになろう．

6． お わ り に大分県杵築市山香町にガスを伴って湧出する高塩分の温泉を対象とし，本格的に温泉水および温泉付随ガス試料のを多項目の化学・同位体分析を行い，それらのデータの地球化学的解析結果に基づいて温泉の生成機構について考察した．その結果，以下のことがらが示された．
・山香町の温泉は Na と Cl に富む高塩分水と低塩分 HCO3 型の浅層地下水の混合によって形成されている．
・高塩分水の温泉起源流体は，地下深部に潜在したマグマから放出される火山性流体であると考えられる．
・その火山性流体は古い時代のマグマに由来するものである可能性が高く，地下の貯留層に永い時間留まることで相当に低温化したようである．
179酒井拓哉，大沢信二，山田誠，三島壮智，大上和敏

温泉科学Fig. 11 An illustration showing source fluid and formation process of Yamaga hot spring.  
The deep source fluid of Yamaga hot spring waters at present would be a cooling
remnants of saline magmatic aqueous fluid.  Hot spring waters are formed by mixingof the remnants with ordinal shallow ground water of meteoric origin.  The magmatic
fluid is probably to be released from magma being active in the Tertiary.  There is apossibility that Bajo gold deposit, which is classified into low sulphidation epithermal
gold deposit, was formed from the magmatic hydrothermal system at the same time.

本研究では以上のことを明らかにすることができたが，山香温泉のような中性の塩化物泉（Na-Cl
型水質）は，Oki and Hirano（1974）によって提案されている箱根火山の成因モデルにもあるように，最も高温の火山中央の深部においてマグマから放出される NaCl に富む高圧蒸気が起源となっ180

第 63 巻（2013）

温泉水・温泉付随ガスの地球化学データから見た大分県山香温泉の生成機構と温泉起源流体て生成されているとされており，一方，酸性高温噴気ガスは HCl, SO2 といったマグマ揮発性成分
が地下の比較的浅所でマグマから脱ガスしたもので，それが地下水に流入することで酸性の塩化物泉（H-Cl-SO4 型水質）を生成するとされている．5．1 項の δD と Cl－の関係（Fig. 7）は，マグマ性
流体のマグマからの脱ガス深度は流体の液性のみならず，塩分濃度にも違いを生ぜしめることを示唆している．マグマ性流体の脱ガス深度と塩分濃度を知る資料はほとんどないが，深度が増すと塩分濃度（Cl 濃度）は増大すると想像される．
また，冒頭にも示したように，本研究地域には形成年代が約 500 万年前の浅熱水性金鉱床（星住・森下，1993；Watanabe, 2005）が存在し，温泉との関連性に興味がもたれる．浅熱水性鉱床の最高形成温度は 200℃程度であり（森下，1994），それが本研究で示された温泉起源流体の推定温度 210℃に近いことは両者に何らかの関連性があることを暗示しているように思える．しかし一方で，浅熱水性鉱床を形成する熱水は一般に天水性であるとされ（森下，1994，2004），山香町の温泉水にはマグマ性流体の混入が明瞭にみられるという違いがある．これらのことを都合よくとらえ
れば，高塩分のマグマ性熱水流体貯留域の周辺部に生じた天水循環熱水系が金鉱床を生んだということが考えられる（Fig. 11）．この点を想像から脱却させるためには，その有効な手段として高塩分温泉水の年代測定があげられる．温泉水の生成年代を探る手段として，ヨウ素 129（129I）法や放射性炭素同位体（14C）を用いた深層地下水の年代推定法，4He 蓄積法による地下水の年代測定法の応用が考えられるが，それぞれの方法の年代適用範囲や試料採取法などをよく検討し，山香町の温泉水の年代測定を行うことを今後の重要な研究課題として提案したい．
謝辞
現地調査と試料採取は，泉源の所有者や管理者の方々のご理解のもとに円滑に行うことができた．また，水素・酸素・炭素の安定同位体分析についてはニュージーランド地質核科学研究所の安定同位体研究室（The Stable Isotope Laboratory, GNS Science Limited）からの，付随ガスの炭化
水素の化学・同位体分析では石油資源開発㈱技術研究所の，付随ガスの希ガス同位体測定に際しては九電産業株式会社環境部からそれぞれ協力・支援をいただいた．ここに感謝の意を表します．この研究の費用の一部は，独立行政法人日本学術振興会科学研究費補助金基盤研究（A），研究課題
番号：22244067（研究代表者：平島崇男）を用いた．関係各位に感謝申し上げる．
引用文献
阿部喜久男，茂野博（1986）：日本のいくつかの地熱地域における熱水・温泉水の Li, Cs 含有量とその地球化学的意味．地質調査所月報，37，12，625-646.
相川嘉正（1991）：千葉県の温泉の化学─温泉およびガス田付随水の化学成分─．温泉科学，42，2，
22-34.
Allis, R.G. and Yusa, Y. (1989) : Fluid flow processes in the Beppu geothermal system, Japan.
Geothermics, 18, 743-759.
網田和宏，大沢信二（2001）：噴気ガス採取用改良注射器法．九大地熱・火山研究報告，10，158167.
網田和宏，大沢信二（2003）：九重硫黄山噴気地域から放出されるマグマ性ガスへの空気及び地下の水の混合過程．日本地熱学会誌，25，4，245-265.
網田和宏，大沢信二，杜建国，山田誠（2005）：大分平野の深部に賦存される有馬型熱水の起源．

温泉科学，55，2，64-77.
Deines, P. (2002) : The carbon isotope geochemistry of mantle xenoliths. Earth-Science Reviews,
181酒井拓哉，大沢信二，山田誠，三島壮智，大上和敏

温泉科学58, 247-278.
Giggenbach, W.F. (1988) : Geothermal solute equilibria. Derivation of Na-K-Mg-Ca geoindicators.
Geochim. Cosmochim. Acta., 52, 2749-2765.
Giggenbach, W.F. (1992) : Isotopic shifts in waters from geothermal and volcanic system along
convergent plate boundaries and their origin. Earth Planet. Sci. Lett., 113, 495-510.
Giggenbach, W.F. (1997) : Relative importance of thermodynamic and kinetic processes in governing
the chemical and isotopic composition of carbon gases in high-heat flow sedimentary basin.
Geochim. Cosmochim. Acta, 61, 17, 3763-3785.
Giggenbach. W.F. and Goguel. R.L. (1989) : Collection and analysis of geothermal and volcanic
water and gas discharges. Fourth Edition. Report CD2401, Chemistry Division, DSIR, 81,
New Zealand.
Goff, F. and Janik, C.J. (2002) : Gas geochemistry of the Valles caldera region, New Mexico andcomparisons with gases at Yellowstone, Long Valley and other geothermal systems. J. Volcanol. Geotherm. Res., 116, 299-323.
Goff F. and McMurty, G.M. (2000) : Tritium and stable isotopes of magmatic waters. J. Volcanol.
Geotherm. Res., 97, 347-396.
星住英夫，森下祐一（1993）：豊岡地域の地質．地域地質研究報告（5 万分の 1 地質図幅），地質調査所，p. 75.
岩倉一敏，大沢信二，高松信樹，大上和敏，野津憲治，由佐悠紀，今橋正征（2000）：長湯温泉（大分県）から放出される二酸化炭素の起源．温泉科学，50，86-93.
鎌田浩毅（2002）：フィリピン海プレート北縁部に見られる広域テクトニクスの転換 . 京都大学防災研究所
北研究集会 13K-7 報告書，10-22.

逸郎（1995）：別府─島原地溝下のプレートの沈み込み状態の怪─そのマグマ性 N2/Ar と 3He/
He 比の意味─．地質ニュース，488，19-23.

4Kita, I., Nitta, K., Nagao, K., Taguchi, S. and Koga, A. (1993) : Difference in N 2/Ar ratio ofmagmatic gases from northeast and southwest Japan : New evidence for different states ofplate subduction. Geology. 21, 391-394.
Mizutani, Y., Hayashi, S., Sugiura, T. (1986) : Chemical and isotopic compositions of fumarolicgases from Kuju-Iwoyama, Kyushu, Japan. Geochem. J., 20, 273-285.
森下祐一（1994）：マグマによる熱水性鉱床の生成─鉱脈型，スカルン型及び斑岩型鉱床から得られた 2，3 の知見─．地質学論集，43，98-108.
森下祐一（2004）：北薩・串木野地域における一つの金鉱床成因研究物語．地質ニュース，599，
59-68.
Muramatsu, Y., Fehn, U., Yoshida, S. (2001) : Recycling of iodine in fore-arc areas : Evidence fromthe iodine brines in Chiba, Japan : Earth Planet. Sci. Lett., 192, 583-593.
村松容一，近藤史也，千葉仁，早稲田周，長島秀行（2010）：関東山地北縁における非火山性温泉の水質および安定同位体比とその地質学的解釈．温泉科学，60，4-21.
NEDO（1989）：平成元年全国地熱資源総合調査（第 3 次）広域熱水流動系調査 鶴見岳地域 流体
地化学調査，92-93.
西尾嘉朗（2010）：リチウム同位体が拓く地殻流体科学．地質ニュース，670，13-21.
Nishio, T., Sasaki, S., Gamo, T., Hiyagon, H., Sano, Y. (1998) : Carbon and helium isotope systematics
of North Fiji Basin basalt glasses : carbon geochemical cycle in the subduction zone. Earth182

第 63 巻（2013）

温泉水・温泉付随ガスの地球化学データから見た大分県山香温泉の生成機構と温泉起源流体

Planet. Sci. Lett., 154, 127-138.
Ohba, T., Hirabayashi, J., Nogami, K., Kusakabe, M., Yoshida, M. (2008) : Magma degassing processduring the eruption of Mt. Unzen, Japan in 1991 to 1995 : Modeling with the chemicalcomposition of volcanic gas. J. Volcanol. Geotherm. Res., 175, 120-132.
大沢信二（2000）：噴気ガスの化学・同位体組成からみた別府温泉の地熱流体の起源及び性状．大分県温泉調査研究会報告，51，19-28．
大沢信二（2003）：塚野鉱泉の水質形成機構．大分県温泉調査研究会報告，59，295-319．
大沢信二（2005）：大分県南部地域の温泉の同位体地球化学調査─臼杵市の六ヶ迫鉱泉─．大分県
温泉調査研究会報告，56，5-9.
Ohsawa, S., Yusa, Y., Oue, K., Amita, K. (2000) : Entrainment of atmospheric air into the volcanic
system during the 1995 phreatic eruption of Kuju Volcano, Japan. J. Volcanol. Geotherm.
Res., 96, 33-43.
大沢信二，網田和宏，山田誠，三島壮智，風早康平（2010）：宮崎平野の大深度温泉井から流出する温泉水の地化学特性と成因─温泉起源流体としての続成脱水流体─．温泉科学，59，4，
295-319.
Oki, Y. and Hirano, T. (1974) : Hydrothermal System and Seismic Activity of Hakone Volcano.
The Utilization of Volcanic Energy, 13-40, Proc. Us-Japan Seminar.
産業技術総合研究所地質調査総合センター（2011）：20 万分の 1 日本シームレス地質図データベース．産業技術総合研究所研究情報公開データベース DB084，産業技術総合研究所地質調査総合センター，産総研著作物管理番号：H17PRO-316.
Sano, Y. and Marty, B. (1995) : Origin of carbon in fumarolic gas from island arcs. Chem. Geol.,
119, 265-274.
Shimizu, A., Sumino, H., Nagao, K., Notsu, K., and Mitropoulos, P. (2005) : Variation in noble gas
isotopic composition of gas samples from the Aegean arc, Greece. J. Volcanol. Geotherm.
Res., 140, 321-339.
早稲田周，岩野裕継（2007）：ガス炭素同位体組成による貯留層評価．石技誌，72，585-593.
Watanabe, Y. (2005) : Late Cenozoic evolution of epithermal gold metallogenic provinces in
Kyushu, Japan. Mineralium Deposit, 40, 307-32.
山下幸三郎（1982）：山香町内の鉱泉・湧水の水質．「山香町誌 第二編 地誌，第四章」，p. 163-p.
171，山香町誌刊行会，大分．

183