温泉地学研究所傾斜観測網で観測された地球潮汐について
原田昌武＊、棚田俊收＊、伊東博＊、丹保俊哉＊

Earth tide observed at tiltmeter network of Hot Springs
Research Institute of Kanagawa Prefecture
By＊
Masatake HARADA , Toshikazu TANADA ＊ , Hiroshi ITO ＊ and Toshiya TANBO ＊２. データおよび解析方法

１． はじめに神奈川県温泉地学研究所では、１９８９（平成元）年から箱解析に用いたデータは傾斜観測網の全７観測点の傾斜根カルデラ内の駒ケ岳・湖尻・小塚山の３ヶ所においてデータである。それぞれの観測点では南北、および、東傾斜計による観測を開始し、箱根地域における火山活動のモニタリングを行っている。その後、さらに塔の峰・岩倉・寄・裾野の４観測点を増設し、現在では７観測点による傾斜観測網により箱根火山のみならず、神奈川県西部地域の広域地殻変動をモニタリングしている。傾斜観測網の観測点配置を図 1 に示す。これらの各観測点では、深度約１００ｍの観測井内にボアホール型傾斜計を設置し、南北方向および東西方向の２成分の傾斜観測を西方向の２成分観測を行っているため、全部で１４成分が対象となる。これらのデータはリアルタイムで温泉地学研究所に送信されており、毎分０秒の値を１分値ファイルとして保存している。地球潮汐の解析には、この１分値ファイルから毎時０分の前後３０分間で平均した値を１時間値としてリサンプリングしたものを使用した。解析に用いたデータの期間は１９９９（平成１１）年７月１日から２００３（平成１５）年７月３１日の約４年間のデータである。
通常の傾斜観測のデータに記録される地球潮汐には、行っている（温泉地学研究所、１９９９）。傾斜計は１０-９rad 程度の変動も検出できる非常に高感度なセンサーであるため、地震活動や火山活動に伴う変動を捉えることができる。例えば、２００１（平成１３）年６月から発生した箱根群発
地震活動の際には、これに伴う傾斜変動が検出され、その変動源の推定も行われている（代田ほか、２００２ ; 代田ほ

２つの潮汐変化の合力が観測されている。１つは、固体地球自体が外部天体の引力の作用によって変形が生じさせられる固体地球潮汐である。つまり、引力が地球の固体部分に直接作用していることを示している。２つ目は、地球の流体部分によって引き起こされる海洋潮汐である。
これは、海洋の潮汐変化により海底にかかる海水の荷重か、２００３ａ）。
このように観測された傾斜変化を用いて定量的に地殻活動を議論するためには、傾斜計の感度検定を行う必要がある。すなわち、観測点近傍における変動が傾斜計によって正確に観測されているかどうかを検証しなければが変化し、それによって地殻部分が弾性変形すること（さらに、弾性変形に伴う引力変化や海水自体の引力変化）に起因する潮汐現象である。これら固体地球潮汐と海洋潮汐の影響を理論的に予測するために、地球潮汐・ならない。なぜなら、１０-９rad 程度という極めて微小な量の地殻の傾斜を測定する傾斜計の計器感度を通常の観測35.4°N寄範囲内で直接求めることは困難であり、ほとんどの場合これよりも大きな変動に対する計器の感度から内挿され岩倉たものを通常の観測範囲での計器感度としているためである。これによって得られた変動量の誤差は、伸縮計の小塚山
湖尻場合には±４０％以上になることもある（竹本、１９７１）。そこで本報告では、通常の観測時に記録され、なおかつ、既塔の峰知の周期的な変化をもたらす地球潮汐をシグナルとして利用し、地球潮汐によって引き起こされる傾斜変動の観35.2°N

裾野

駒ケ岳測値とその理論的に予測される値とを比較することによって各観測点での計器特性を明らかにした。また、観測計器の安定性および箱根群発地震活動が地球潮汐に及ぼす影響を検討するために地球潮汐傾斜の振幅の時間変139.0°E化を計算したので報告する。

図１ 傾斜観測点配置＊
神奈川県温泉地学研究所〒２５０-００３１ 神奈川県小田原市入生田５８６
報告，神奈川県温泉地学研究所報告，第３６巻，４７−５２，２００４．

−４７−139.2°E

海 洋 潮 汐 荷 重 計 算 プ ロ グ ラ ム GOTIC（Global Oceanic

駒ケ岳観測点では、東西成分の M2 分潮の位相が９°ず

Tidal Correction；Sato and Hanada, １９８４）の 改 良 版 で あ る

れているが、その他の成分・分潮では３°以内で非常にGOTIC２（Matsumoto et al ., ２００１）を 用 い た。一 方、観 測良く一致している。振幅については東西成分の O1 分潮をデータの潮汐解析には潮汐解析プログラム BAYTAP-G除き１５％程度大きくなる傾向がある。駒ケ岳観測点につ（Bayesian Tidal Analysis Program - Grouping Model；
Ishiguro et al .,１９８１ ; Tamura et al .,１９９１）を用い、地球潮汐によって引き起こされる傾斜変動の振幅と位相を計算した。観測データは固体地球潮汐と海洋潮汐の２つに分離することができないため両者の和である地球潮汐、その中でも特に振幅が大きく他の分潮との分離もよい M2 分潮（月起源の半日周潮）、O1 分潮（月起源の日周潮）について観測値と理論値を比較する。

１８０°ずれていると指摘されており（代田ほか、２００３ｂ）、本解析においても出力極性を反転したデータを用いた結果、観測値と理論値はほとんど調和的な結果が得られた。
湖尻観測点では、東西成分の M2 分潮の位相が１３°程度ずれてはいるが、その他の成分・分潮では３°以内で非常に良く一致し、また、振幅についてもほとんど理論値と同じ値が観測されていることがわかる。
小塚山観測点は、位相は５°以内で良く一致しているが、南北成分の M2 分潮と東西成分の O1 分潮において３. 各観測点での潮汐定数
BAYTAP-G は観測データをトレンド、地球潮汐、気圧または温度等による応答、残差の４つの成分に分解し、その結果を出力する。解析結果の一例として、駒ケ岳観測点の南北成分の傾斜計によって得られたこれらの成分を図２に示す。上から順に観測データ（入力）、トレンド、潮汐成分、残差成分である。なお、各観測点では地中温度の観測も行っているが、計測機器の不調による異常値や欠測が多いため、本報告では温度による応答は計算に含めなかった。
このようにして解析された各観測点での M2 分潮、O1分潮の結果を表１に理論潮汐の計算結果も含めて示す。
また、ｐｈａｓｏｒｄｉ
ａｇｒａｍの形にしたものを図３に示す。原ａｄ）を示し、ま点からの距離が振幅の大きさ（単位は１０-９ｒた、上方０°からの角度が位相（遅れを正）に相当する。
図中の○および□はそれぞれ M2 分潮、O1 分潮の理論値を、また、＋および×はそれぞれ M2 分潮、O1 分潮の観測値の振幅・位相を示している。

いては、傾斜計の観測方位が東西方向・南北方向ともに振幅が１８％程度減少している。
塔の峰観測点は、南北成分のＭ２分潮の位相が７°ずれてはいるものの、その他は良く一致しており、東西成分の位相は若干進む傾向がみられる。また、南北成分の振幅は１５％程度減少する傾向にあるが、東西成分は逆に２０％程度増加しており、何らかの系統的な誤差要因があるかもしれない。その可能性として考えられるものは、（１）観測計器の設置方向の東西・南北方向からのずれ、（２）観測点近傍の局所的な地質構造（厳密には物性パラ
メ ー タ）が、今 回 用 い た１０６６Ａ地 球 モ デ ル（Gilbert andDziewonski, １９７５）とはずれている、（３）実際の地形との相違に起因するモデル計算からのずれ（例えば、Berger andBeaumont, １９７６）、がある。しかし、ここではこのような不均一な地殻構造モデルの影響を含めた地球潮汐の解析を行っていないため、結論付けることはできない。
岩倉観測点では、南北方向の M2 分潮の振幅は１０％減
少し、東西成分の M2 分潮の位相は７°ずれてはいるものの、その他の潮汐定数（振幅および位相）は理論値と非常に調和的である。
寄観測点は、南北成分の M2 分潮の振幅が４０％増加しているが、その他の成分・分潮ではおおよそ調和的であり、位相については２°以内で非常に良く一致している。
裾野観測点については、位相のずれが他の観測点に比べて大きくなる傾向がある。また、南北成分の振幅は増加し、東西成分では減少している。しかし、位相には系統的な変化がみられないため、塔の峰観測点のような単純な問題では解決できないかもしれない。
図４に駒ケ岳観測点、湖尻観測点、小塚山観測点におけるＭ２分潮の振幅の時間変化を示す。これは、地球潮汐を解析する時間ウィンドウを３０日と固定し、１日ずつウィンドウをずらしながら振幅の時間変化を求めたものである。これらの解析結果から、駒ケ岳観測点の南北成分には年周的な変化も見られるが、全成分とも５ nano

図２ ２００３（平成１５）年１月の駒ケ岳観測点（南北成分）の観測データを用いた解析結果の例
上から順に観測データ（入力）、トレンド、潮汐成分、残差成分を示す。

radian 程度の変動に収まっており、安定した観測が継続されていると考えられる。一般的に、地球潮汐は観測点−４８−

表１ 各観測点での地球潮汐定数ධർᚑಽ࿾⃿ầ᳤૏⋧ ᝄ᏷㚤ࠤጪ 0.00 24.99
0.00 25.00
ḓዥዊႦጊ 0.00 25.00
Ⴁߩፄ 0.00 25.00
0.00 25.03
ጤୖ
0.00 25.04
ነ
0.00 24.99
⵿㊁ධർᚑಽ࿾⃿ầ᳤૏⋧ ᝄ᏷㚤ࠤጪ 0.00 7.23
0.00 7.22
ḓዥዊႦጊ 0.00 7.20
Ⴁߩፄ 0.00 7.21
0.00 7.14
ጤୖ
0.00 7.11
ነ
0.00 7.23
⵿㊁᧲⷏ᚑಽ࿾⃿ầ᳤૏⋧ ᝄ᏷㚤ࠤጪ 0.00 43.13
0.00 43.12
ḓዥዊႦጊ 0.00 43.11
Ⴁߩፄ 0.00 43.12
0.00 43.07
ጤୖ
0.00 43.04
ነ
0.00 43.13
⵿㊁
᧲⷏ᚑಽ
࿾⃿ầ᳤
૏⋧ ᝄ᏷
㚤ࠤጪ 0.00 12.70
0.00 12.70
ḓዥዊႦጊ 0.00 12.71
Ⴁߩፄ 0.00 12.71
0.00 12.74
ጤୖ
0.00 12.75
ነ
0.00 12.70
⵿㊁㧹㧞
ᶏᵗầ᳤
૏⋧ ᝄ᏷
-171.49 51.10-173.10 48.72-171.87 49.37-166.99 58.33-163.38 63.33-171.33 45.35-174.40 49.19

ว⸘㧔ℂ⺰୯㧕૏⋧ ᝄ᏷-163.51 26.64-165.93 24.09-163.69 24.87-157.58 34.44-153.06 40.00-160.93 20.94-168.68 24.44

᷹ⷰ୯૏⋧ ᝄ᏷-162.25 30.58-163.33 24.11-163.13 20.35-150.63 29.58-150.50 36.01-161.55 29.16-172.78 31.02

㧻㧝
ᶏᵗầ᳤
ว⸘㧔
ℂ⺰୯㧕
૏⋧ ᝄ᏷ ૏⋧ ᝄ᏷
9.28 23.95
7.13 31.11
8.62 22.55
6.53 29.70
9.08 23.16
6.93 30.29
11.16 28.22
8.90 35.33
12.44 31.46
10.15 38.46
9.01 21.72
6.79 28.76
8.24 22.30
6.23 29.47

᷹ⷰ୯
૏⋧ ᝄ᏷
9.83 34.76
8.78 29.93
8.36 31.75
9.64 29.91
11.60 35.72
7.04 32.76
5.66 34.20

㧹㧞
ᶏᵗầ᳤
ว⸘㧔
ℂ⺰୯㧕
૏⋧ ᝄ᏷ ૏⋧ ᝄ᏷
159.51 36.03
53.37 15.72
166.32 30.37
27.81 15.39
162.69 32.38
38.31 15.54
149.26 48.77
87.23 24.96
157.67 35.19
51.81 17.01
167.27 27.36
20.23 17.43
172.75 27.60
12.46 16.13

᷹ⷰ୯
૏⋧ ᝄ᏷
62.44 18.50
41.14 15.92
32.86 16.84
83.17 29.13
44.68 16.96
22.08 16.69
5.96 14.32

㧻㧝
ᶏᵗầ᳤
ว⸘㧔
ℂ⺰୯㧕
૏⋧ ᝄ᏷ ૏⋧ ᝄ᏷
124.34 22.44
89.89 18.53
126.83 19.02
85.12 15.28
125.66 20.37
87.12 16.57
120.87 29.68
95.66 25.60
124.13 22.32
89.34 18.48
127.78 17.65
82.10 14.09
128.95 17.01
81.39 13.38

᷹ⷰ୯
૏⋧ ᝄ᏷
88.91 19.69
86.49 15.74
90.71 13.86
94.19 31.02
87.30 17.59
82.53 15.71
89.90 13.01

周辺の弾性定数にも影響を受けるため、観測点近傍でのを及ぼすほどの規模ではなかった可能性が考えられる。

地震活動に伴う弾性定数の変化（ダイラタンシーモデル）や火山活動に伴う地殻構造の変化があった場合には、４. まとめ

傾斜計が観測する地球潮汐も変化をすることが予想され傾斜観測網の全観測成分について地球潮汐の解析を行る（観測された事例については例えば、Nishimura, １９５０ ;い、理論潮汐と比較した結果以下のことが確認された。

Latynina and Rizaeva, １９７６ ; Mikumo et al ., １９７７ ;各観測点における潮汐定数の位相は、湖尻観測点にYamauchi, １９８９）。しかし、
２００１（平成１３）年６月から箱根おける東西成分の M2 分潮を除き、観測値と理論値

火山で始まった、現在のテレメータ観測を開始した１９８９（平成元）年４月以降で最大の群発地震活動（棚田ほか、が非常に良く一致している。
振幅についても観測値と理論値は調和的であり、当所の傾斜観測網は、地球潮汐の周期帯でみる限り地２００２）に関連するような潮汐定数の変化は観測されてい殻の傾斜変動を非常に良く捉えていることが明らかない。これは、群発地震活動そのものが地球潮汐に影響となった。

−４９−図３ 地球潮汐定数（Phasor diagrams）
原点からの距離が振幅の大きさ（単位は１０-９rad）を示し、また、上方０°からの角度が位相（遅れを正）に相当する。
図中の○および□はそれぞれ M2 分潮、O1 分潮の理論値を、また、＋および×はそれぞれ M2 分潮、O1 分潮の観測値の振幅・位相を示している。

−５０−

図３ 地球潮汐定数（Phasor diagrams）
原点からの距離が振幅の大きさ（単位は１０-９rad）を示し、また、上方０°からの角度が位相（遅れを正）に相当する。
図中の○および□はそれぞれ M2 分潮、O1 分潮の理論値を、また、＋および×はそれぞれ M2 分潮、O1 分潮の観測値の振幅・位相を示している。

潮汐定数の時間変化は５nano radian 程度の変動に収

代田寧、棚田俊收、伊東博(２００２) ２００１
（平成１３）

まっており、長期間安定した観測が継続している。

年箱根群発地震活動に関連した地殻傾斜変動，温地研

地震・火山活動に伴う変動を捉えるためには、
この ような報告，３４，３５-４４.

高精度の地殻変動観測を今後も継続させる必要がある。

代田寧、丹保俊哉、棚田俊收、伊東博(２００３a)２００１年箱根火山で生じた地殻傾斜変動（その２），日本
参考文献

火山学会２００３年度秋季大会講演予稿集，１
３２.
代田Berger, J. and Beaumont, C.（１９７６）An Analysis of Tidal

寧、棚田俊收、原田昌武、伊東博(２００３b)Strain Observations from the United States of America II.

潮汐変化と遠地地震を用いた温泉地学研究所傾斜計セ

The Inhomogeneous Tide, Bull. Seism. Soc. Am., ６６, １８２１-ンサー方位の検証， 温地研報告， ３５，３３-４０.

１８４６.

−５１−

温泉地学研究所 （１９９９） 温泉地学研究所における「神奈川県西部地震」の取り組み，温地研報告，２９，３４０.
Sato, T. and Hanada, H.（１９８４）A Program for theComputation of Oceanic Tidal Loading Effects 'GOTIC',Publ. Int. Lat. Obs. Mizusawa, １８, ２９-４７.
竹本修三（１９７１）レーザー干渉計による伸縮計・ひずみ地震 計 の Calibration に つ い て , 京 大 防 災 研 究 所 年 報 ,１４A, ７-１３.
Tamura, Y., Sato, T., Ooe, M. and Ishiguro, M.（１９９１）AProcedure for Tidal Analysis with a Bayesian InformationCriterion, Geophys. J. Int.,１０４, ５０７５１６.
棚田俊收、代田寧、伊東博、袴田和夫（２００２）
２００１（平成１３）年箱根火山の群発地震活動について , 温地研
観測だより , ５２, １-４.
Yamauchi, T.（１９８９）Earthquake Occurrence and Temporal
Variations in Earth Tidal Strain Amplitudes, J. Geod. Soc.,
Japan, ３５,１４９-１５７．

図４ Ｍ２分潮の振幅の時間変化(a) 駒ケ岳、(b) 湖尻、(c) 小塚山の解析結果

Gilbert, F. and Dziewonski, A. M.（１９７５）An Application ofNormal Mode Theory to the Retrieval of Structural
Parameters and Source Mechanisms from Seismic Spectra,
Phil. Trans. R. Soc. Lond., A２７８, １８７-２６９.
Ishiguro, M., Akaike, H., Ooe, M. and Nakai, S.（１９８１）A
Bayesian Approach to the Analysis of Earth Tides, Proc.
９th Int. Sympos. Earth Tides, New York, ２８３-２９２.
Latynina, L. A. and Rizaeva, S. D.（１９７６）On tidal-strain
variations before earthquakes, Tectonophysics ３１, １２１-１２７.
Matsumoto, K., Sato, T., Takanezawa, T. and Ooe, M.
（２００１）GOTIC2: A Program for Computation of Oceanic
Tidal Loading Effect, J. Geod. Soc., Japan, ４７, ２４３-２４８.
Mikumo, T., Kato, M., Doi, H., Wada, Y., Tanaka, T., Shichi,
R. and Yamamoto, A.（１９７７）Possibility of Temporal
Variations in Earth Tidal Strain Amplitudes Associatedwith Major Earthquakes, J. Phys. Earth ２５, １
２３-１３６.
Nishimura,E.（１９５０）Onearthtides,Trans.

Am.

Geophys.Union ３
１, ３
５７-３７６.

−５２−