東海大学海洋研究所研究報告第 33 号（2012）13 − 21 頁

Bull. Inst. Oceanic Res. & Develop., Tokai Univ.（2012）, 33, 13 − 21

平成 23 年東北地方太平洋沖地震（Mw 9.0）の前後における奥多摩地磁気 3 成分データの統計的パラメータ変化
竹 内 昭 洋 1）・長 尾 年 恭 1）

Statistical Parameter Variations in the 3 Components of Geomagnetic Field at Okutama
in Central Japan before and after the 2011 off the Paciﬁc Coast of Tohoku Earthquake
Akihiro Takeuchi 1） and Toshiyasu Nagao 1）

Abstract
We have observed the 3-component geomagnetic ﬁelds with a ﬂuxgate magnetometer at 32 Hz samplingin Okutama in central Japan since December 2003. Okutama is located just outside the area where the seismic
activity increased after the 2011 off the Paciﬁc coast of Tohoku Earthquake. In this study, we conﬁrmed
whether or not there are changes in statistical parameters of the observed data before and after the earthquake.
Because the daytime data included noises mainly from trains driven by DC electric power, we used only thenighttime data from 2 : 30 : 00 to 3 : 04 : 08 in JST after 1-sec averaging. The power spectrum densities, especiallyof pc1−3, had seasonal variations. The median, 1 st quartile, and 3 rd quartile were calculated for the differences
between the power spectrum intensities and the seasonal variations from March 12 to October 31 each year.
As a result, we conﬁrmed that the median and 1 st quartile slightly decreased after the earthquake while the
3 rd quartile increased. One of the possible reasons of the decreases is a decline of the use of electric power
after the earthquake. Moreover, one of the possible reasons of the increase may be a change of the local
groundwater conditions after the earthquake.

緒2011），微震動を伴う火山活動も活発化した（Yukutake

言et al., 2011）．

日 本 時 間 で 2011 年 3 月 11 日 14 時 46 分 に， モ ー メ

さらに，東北地震は大規模な津波の発生も伴った

ントマグニチュード 9.0 の平成 23 年東北地方太平洋（Maeda et al., 2011）．この津波により発生した大気重

沖地震（以下，東北地震と略す）が発生した（Hirose et

力波が電離層まで達し，電離層のプラズマを擾乱させal., 2011）．また，余震のみならず，その周辺域におた（Liu et al., 2011；Tsugawa et al., 2011）．また，地ける地殻ひずみの再配分により誘発されたと推測さ

震発生後の地磁気データには比較的大規模な変化も観れるマグニチュード 6 級の地震も多発し（Okada et al.,

測された．この変化を津波という導体が地球磁場中を

1） 東海大学海洋研究所地震予知研究センター〒 424−8610 静岡県静岡市清水区折戸 3−20−1

Earthquake Prediction Research Center, Institute of Oceanic Research and Development, Tokai University, 3−20−1
Orido, Shimizu−ku, Shizuoka 424−8610, Japan（2011 年 12 月 28 日受付／ 2012 年 2 月 18 日受理）

14

竹内昭洋・長尾年恭

移動することにより発生する誘導電流により 2 次的に

磁気 3 成分を常時観測している．計測器本体にデータ誘導される磁場（津波ダイナモ効果）であるとする報告が保存される．観測データは，ISDN 回線を用いて清もある（Utada et al., 2011）．このように，東北地震は

水のセンター研究室内に設置されたデータ解析サーバ

その発生に伴って電磁気学的な変化を誘導した．こにより毎日自動回収されている．回収された観測デー

れらは過渡的な電磁気現象であると言えよう．それでタは，自動処理プログラムにより 1 秒平均値化され，は，もっと長周期または定常的な地磁気脈動のような観測結果が関係者に毎日メールで配信されている．

電磁気学的パラメータに変化はないのであろうか．

観測データの 1 例として，2011 年 3 月 11 日（日本時東海大学海洋研究所地震予知研究センター（以下，セ間）のプロットを Fig. 2 に示す．3 成分観測値（X，Y，ンターと略す）では，東京都の奥多摩町で 2003 年 12 月

Z）に加え，それらから求められる水平成分（H）と全磁より地磁気 3 成分の常時観測を行っている．奥多摩は地震学的活性化領域の外縁に位置しており，東北地震の前後において地磁気特性に変化の現れていることが期待される．そこで本研究では，約 8 年にわたる時系列デー
タのスペクトルを統計解析することにより検証する．

奥多摩地磁気観測データ
Fig. 1 に示すように，センターの奥多摩地磁気観測
点のセンサーは，JR 青梅線の奥多摩駅から日原街道（都道 204 号線）を約 10 km 北上し，街道から北東方向に約 100 m 逸れた東京都水道局の水源かん養林用地内（35°50 56 N，139°02 19 E）にある．地磁気センサー
はフラックスゲート型（テラテクニカ製 FR−725）であり，沢の斜面を横切る小道沿いの地中約 1 m の深さに埋設してある．この小道は普段は渓流釣り客が稀に通る程度である．計測器本体（データロガー）は日原街道沿いの観測小屋に設置され，アンプおよびケーブルは小道の脇に埋設されている．32 Hz サンプリングで地Fig. 1 Location of the Okutama geomagnetic ﬁeld
observation site. This area is about 70 km
west of Tokyo central part. The sensor（35°50 56 N, 139°02 19 E）is a ﬂuxgate
type buried in the ground at 1 m depth.

Fig. 2 The 1-sec averaged geomagnetic ﬁelds at the Okutama site.（Left columns）The 1-day variations ofeach component on March 11, 2011 JST.（Right columns）The 15-minute variations of each component
before and after Tohoku earthquake. H and F components are calculated from X, Y and Z components.

東北地震前後の奥多摩地磁気パラメータ変化

15

力（F）のプロットも示す．昼間の時間帯，特に Z 成分

「R」の spec.pgram() 関数を用いた．2048 秒間の時系列にはノイズが卓越している．一般的にこのようなノイデータから平均値を引いた後に，高速フーリエ変換にズの主な原因は，直流電源により駆動する電車からのよりペリオドグラムを計算した．得られたパワースペ

漏洩電流であると考えられている．東北地震の地震動

クトル密度を周期に対して 6 つに分けた．これらは国により電車の運行が止まったと同時にノイズレベルが際地球電磁気・超高層物理学協会による連続型地磁気

格段に小さくなっていることが確認できる．また，地脈動の分類に準じており，それぞれ，

震波到達時には約 15 秒周期の減衰振動が地磁気各成

pc1（周期 2−5 秒）分に現れる．過去の地磁気や地電位差の観測においてpc2（周期 5−10 秒）

も，このような地震動に伴う信号の検出がしばしば報

pc3（周期 10−45 秒）

告されており，流動電位との関係（Nagao et al., 2000；

pc4（周期 45−150 秒）

Widarto et al., 2009）や大地という導体が地球磁場中

pc5（周期 150−600 秒）を振動することにより発生する誘導電流により 2 次的pc6（周期 600−2048 秒）に誘導される磁場（地震ダイナモ効果：Honkura et al.,に相当する．以下，地磁気 H 成分の pc1（周期 2−5 秒）

2004；Ujihara et al., 2004）などが示唆されている．今を例に取り上げ，Fig. 4 に示される一連の統計解析結回の減衰振動波形が単なるセンサーの振動によるもの

果を得るための解析手順を説明する．

なのか流動電位などの物理現象によるものなのかは，
現在のところ不明である．

地磁気データの統計解析手順

地磁気 H 成分 pc1 の範囲には，1 日あたり 615 のスペクトル値が存在する．このデータセットから中
央値，第 1 四分位数および第 3 四分位数の 3 パラメー
タを求め，その日の値とする．この計算を全観測期間

解析期間は日本時間で 2004 年 1 月 1 日から 2011 年

にわたって行なう．次に，東北地震発生前までの各年

10 月 31 日までの約 8 年間である．本研究では，経験の同一日のデータにおける 3 パラメータのそれぞれの的に電車由来ノイズのほとんどないと思われる夜中

中央値を求める．つまり，例えば，2004 年から 2011の 1 秒平均値化データのみを用いた．具体的には，日年の 1 月 1 日（計 8 日）における中央値，第 1 四分位数

本時間で 2 時 30 分 00 秒から 3 時 04 分 08 秒までの 2048

および第 3 四分位数のそれぞれの中央値を求め，1 月

秒間（データ数）のデータを採用した．また，以下の

1 日の値とする．3 月 12 日以降の値は，2011 年のデー

データ解析では，基本的に H 成分と Z 成分についてタを使用しないため，計 7 日のデータセットから計算解析を実施した．

されることになる．この計算から得られた時系列デー

日本時間で 2011 年 3 月 11 日 2 時 30 分 00 秒から 3 時

タが Fig. 4a および Fig. 4i である．このグラフだけ

04 分 08 秒までのデータに関するパワースペクトルを

縦軸は対数表示になっている．青線は各日の中央値でFig. 3 に示す．スペクトル解析には，統計解析ソフト

あり，黒線は中央値の日付中央値，2 本の灰色線は第1 四分位数と第 3 四分位数の日付中央値である．中央値の日付中央値（黒線）と第 1 四分位数と第 3 四分位数の日付中央値（灰色線）は，1 年毎（365 日毎）に同じ曲
線を繰り返させている．ちなみに，これら日付中央値（365 日周期）の計算のために，2004 年 2 月 29 日と 2008
年 2 月 29 日のデータは除かれている．
次に，中央値（青線）と中央値の日付中央値（黒線）の差分を取った時系列データが Fig. 4b および Fig. 4j
である．さらに，各年の 3 月 12 日から 10 月 31 日における差分値の中央値，第 1 四分位数および第 3 四分位

Fig. 3 Power spectrum densities of H（a）and Z（b）
components of the geomagnetic ﬁeld at the
Okutama site from 2 : 30 : 00 to 3 : 04 : 08（for2048 sec）on March 11, 2011 JST.

数を求め，プロットしたのが Fig. 4c である．
Fig. 4 には，得られた解析結果を後章で検討するため，地磁気活動度，地震活動度，気象条件および太陽16

竹内昭洋・長尾年恭活動度に関する時系列データも加えてある．Fig. 4d

観測点から最寄りとなる防災科学研究所の奥多摩町

および Fig. 4k は，気象庁の 3 ヶ所の地磁気観測点（柿氷川観測点における有感地震のマグニチュードと震岡，女満別，鹿屋）における日本時間で 0−3 時の K 指度を示す．Fig. 4g および Fig. 4n は，奥多摩地磁気数の平均値を示す．Fig. 4e および Fig. 4l は，京都大観測点から最寄りとなる気象庁の小河内観測点におけ学大学院理学研究科附属地磁気世界資料解析センター

る日本時間で 2 時台の平均気温と降水量を示す．Fig.

およびドイツの GeoForschungsZentrum Potsdam に

4h は，ベルギー王立天文台の Solar Inﬂuences Data

より算出された日本時間で 0−3 時の Dst 指数と ap 指Analysis Center により算出された毎年 3−10 月にお数を示す．Fig. 4f および Fig. 4m は，奥多摩地磁気

ける太陽黒点数の平均値を示す．

Fig. 4 （a and i）Median of pc1（2 : 30 : 00−3 : 04 : 08 JST）of H component everyday（blue）with medians of
the median（black）and 1 st and 3 rd quartiles（gray）for each date calculated from the pc1 data before
Tohoku earthquake.（b and j）The difference of the pc1 median from the seasonal variation, obtained
from the blue line minus the black line in（a）and（i）.（c）The median（blue）with the 1 st and 3 rdquartiles（gray）of the pc1 median difference（Mar. 12−Oct. 31 JST）.（d and k）The average of
K-indexes（0 : 00−3 : 00 JST）at Kakioka, Memanbetsu, and Kanaya in Japan.（e and l）Dst -index andap - index（0 : 00−3 : 00 JST）.（f and m）Magnitude and seismic intensity of the felt earthquakes at
Okutamacho-Hikawa site.（g and n）The average of air temperature and rain fall（2 : 00−3 : 00 JST）at
Ogouchi site.（h）The average sunspot number（Mar. 12−Oct. 31）.

東北地震前後の奥多摩地磁気パラメータ変化

議論

地磁気 H 成分 pc1 の中央値（Fig. 4a の青線）は，1

17

的であるが，1・5 月のピークは説明できない．また，降水による地温の低下とも関係がなさそうである．現在のところ，この季節変動の原因は不明である．

月頃と 5 月頃の小さなピークおよび 8 月頃の大きな

地磁気 H 成分 pc1 の中央値から季節変動を差し引

ピークを持つ季節変動を有していることが確認できいた値に対して，日本時間で各年の 3 月 12 日から 10る．また，時折大きなスパイク状の値を示す．特に大月 31 日までの期間だけから求めた中央値，第 1 四分位

きなピークの一つが 2011 年 3 月 8 日に発生している．

数および第 3 四分位数（Fig. 4c）
は，東北地震前後におこの日の観測データを Fig. 5 に示す．統計解析に使いて統計的には誤差の範囲以内でおおよそ一致してい

われた日本時間で 2 時 30 分 00 秒から 3 時 04 分 08 秒のるように見える．しかし，やはり年ごとに若干の変化

時間帯に，異常な乱れのあることが確認できる．断定があり，中央値差の中央値（Fig. 4c の青線）に関してはできないが，偶然に不定期の夜行貨物列車が運行さ

は，東北地震前に比べて東北地震後には減少しているれていたか，計測上の問題とも考えられる．他の小・ことが確認できる．また，四分位数範囲（2 本の灰色

中スパイクに関しても同様の疑念を持つことができる線に挟まれた値）は，東北地震後に大きくなっているが，原因を判別することは現段階では極めて困難であのも確認できる．ここで注意すべきことがある．季節る．しかし，本研究では平均値ではなく中央値を採用

変動からの差をとることでは，1 年以上の周期を持つすることにより，少数派である異常に大きなデータの

太陽活動変化の影響を取り去ることはできない．よっ

影響を矮小化できていると仮定し，解析を進めた．

て，Fig. 4c に示された結果には，そのような影響も他機関による地磁気パラメータ（Fig. 4d，4e，4k

含まれていることが予想される．しかし，中央値差のおよび 4l）には，1・5・8 月にピークを持つような季節

中央値の変遷（Fig. 4c）と太陽黒点数の変遷（Fig. 4h）

変動が確認できない．よって，奥多摩観測点の地磁気との間に明確な相関は確認できない．

データに現れる季節変動は広域的なものではないと考以上の統計解析を他の分類（pc1−6）や地磁気成分えられる．有感地震の発生（Fig. 4f および 4m）にも季（H および Z）に対しても行なった．その結果をまとめ

節的な偏りは確認できない．センサーがフラックス

たものを Fig. 6 および Fig. 7 に示す．H 成分と Z 成ゲート型であるため，たとえ地中に埋めていたとして分とでは，基本的に同様な統計的傾向を示す．比較的も，温度特性の効果（おおよそ線形の関係）を疑うことに短周期な pc の中央値（左コラム）には季節変動が見ができる．確かに最も気温が高くなる 8 月（Fig. 4g おられるが，長周期になるにつれ不透明になる．統計的よび 4n）にスペクトル強度のピークが現れるのは調和

誤差の範囲以内とはいえ，全ての pc に対して，中央Fig. 5 The 1-sec averaged geomagnetic ﬁelds at the Okutama site.（Left columns）The 1-day variations ofeach component on March 8, 2011 JST.（Right columns）The 45-minute variations of each component
before and after Tohoku earthquake. H and F components are calculated from X, Y and Z components.

18

竹内昭洋・長尾年恭

値差の中央値（右コラムの色線）が東北地震後に減少し含まれていることが十分に考えられる．例えば，伊ていることが確認できる．ただし，その傾向は周期が豆半島における地磁気・地電流の 100 秒周期成分を主

短い pc ほど顕著である．また，四分位数範囲（2 本の成分解析すると，第 2 主成分に人間活動時間帯と呼応

灰色線に挟まれた値）に関しても，東北地震後に大きする増減が確認できる（Uyeda et al., 2002）．震災前はくなっているのが確認できる．

24 時間操業していた工場などが震災後には停止してもし，東北地震後における中央値差の中央値（Fig.

おり，昼間だけではなく夜間の漏洩電流量も減少して

6 および Fig. 7 の右コラムの色線）が真に減少していいたはずである．しかし，もし使用電力量の低下が原るとすると，その原因は何であろうか．まず考えられ因ならば，第 3 四分位数（Fig. 6 および Fig. 7 の右コるのが，福島原子力発電所事故に起因する節電要請に

ラムの上側灰色線）や四分位数範囲（Fig. 6 および Fig.

よる使用電力量の低下である．
「奥多摩地磁気観測デー

7 の右コラムの 2 本の灰色線に挟まれた値）も減少す

タ」でも述べたように，電車からの漏洩電流が地磁気

べきではないかと考えられるが，相反して解析結果はデータに現れるノイズの主な原因である．しかし，そ増加を示している．そもそも，大型工場由来の漏洩電れ以外にも大型工場などからの漏洩電流による影響も流は 50 Hz または 60 Hz の周波数を持った交流であFig. 6 （Left columns）Median of pc 1−6（2 : 30 : 00−3 : 04 : 08 JST）of the H component everyday（colored）
with medians of the median（black）and 1 st and 3 rd quartiles（gray）for each date calculated from the
pc1−6 data before Tohoku earthquake.（Right columns）The median（colored）with the 1 st and 3 rdquartiles（gray）of the pc1−6 median differences（Mar. 12−Oct. 31 JST）.

東北地震前後の奥多摩地磁気パラメータ変化

19

ると予想され，特に長周期の pc には影響が現れにく

多摩地域においても，地下水の存在状態が東北地震後

いと考えられる．以上のことから，東北地震後におけに変化していることが期待できる．地殻内には，太陽る中央値差（Fig. 6 および Fig. 7 の右コラム）の変化は，

風と電離層との相互作用により発生する地磁気に対す

電力使用量の低下という人工的なもののためだけではる誘導電場（電流）が存在する．更に，地殻内の局所的ないと考えられる．

な電磁気的不均一性により，この誘導電場が局所的に次に考えられる原因は地下水である．東北地震の前変化し，2 次的に誘導される局所的な磁場が発生して後で東日本全域の地殻応力・歪み分布は大きく変わっいることも期待される．地下水の存在状態が変化すれている（Takahashi, 2011）．また，一般的に，地震後ば，電磁気的不均一性も変化するため，誘導電場の局に井戸水の水位が一時的または永続的に変化すること所的変化および 2 次誘導磁場も変化するであろう．以（Roeloffs, 1998；King et al., 1999；Huang et al., 2004）

上のことから，地下水の存在状態が地震発生後に時や地殻応力や地表面の歪みと地下水位との間に相関の間変化するとともに，2 次誘導磁場の強度も増減し，あること（Kümpel, 1992；Ji and Herring, 2012）が報告

結果として各 pc の中央値差の四分位数範囲が大きくされている．これらのことを総合的に判断すると，奥なったのではないかと予想される．

Fig. 7 （Left columns）Median of pc1−6（2 : 30 : 00−3 : 04 : 08 JST）of the Z component everyday（colored）
with the medians of the median（black）and 1 st and 3 rd quartiles（gray）for each date calculated fromthe pc1−6 data before Tohoku earthquake.（Right columns）The median（colored）with the 1 st and3 rd quartiles（gray）of the pc1−6 median differences（Mar. 12−Oct. 31 JST）.

20

竹内昭洋・長尾年恭結論と今後の課題

ベルギー王立天文台の Solar Inﬂuences Data Analysis
Center より太陽黒点数データを提供して頂いた．

センターが奥多摩において約 8 年にわたって観測し

引 用 文 献ている地磁気の水平成分（H）および垂直成分（Z）において，東北地震の発生前後における統計パラメータ

Hirose, F., K. Miyaoka, N. Hayashimoto, T. Yamazakiの変化を検証した．結果として，統計的には誤差の範and M. Nakamura（2011）： Outline of the 2011 off

囲内での変化しか確認されなかった．具体的には，パthe Paciﬁc coast of Tohoku Earthquake（Mw 9.0）

ワースペクトル強度（周期：2−2048 秒，季節変動から−Seismicity：foreshocks, mainshock, aftershocks,の差）の各年 3 月 12 日−10 月 31 日期間における中央値

and induced activity−. Earth Planets Space, 63,

が若干減少し，四分位数範囲が若干増加した．中央値

513−518.

減少の原因の 1 つとして，夜間電力使用量の減少に伴

Honkura, Y., H. Satoh and N. Ujihara （2004）：

うノイズ強度の減少が考えられる．また，中央値減少

Seismic dynamo effects associated with the M7. 1

および四分位数範囲増加の原因の 1 つとして，局所的

earthquake of 26 May 2003 off Miyagi Prefecture

な地下水の存在状態の変化が考えられる．

and the M6. 4 earthquake of 26 July 2003 in

本研究では，スペクトルの時系列データに時折見られるスパイク状の異常値の影響は，統計解析において

northern Miyagi Prefecture, NE Japan. Earth
Planets Space, 56, 109−114.

中央値を用いることにより矮小化できていると仮定しHuang, F.−Q., C.−I. Jian, Y. Tang, G.−M. Xu, Z.

ている．しかし，より精度よく解析するためには，や−H. Deng, and G.−C. Chi （2004）： Response

はり全ての観測データを精査し，明らかな計測エラー

changes of some wells in the mainland subsurfaceや夜行列車の影響のある日のデータを解析対象から丁ﬂuid monitoring network of China, due to the

寧に取り除くような作業が必要であろう．また，センSeptember 21, 1999, Ms7. 6 Chi−Chi Earthquake.

ターの奥多摩 3 成分地磁気は 32 Hz サンプリングで観

Techtonophys., 390, 217−234.

測されているため，1−32 Hz の周波数帯ではどのよ

Ji, K. H. and T. A. Herring （2012）：Correlation

うな統計パラメータの変化が現れるのか確認する必

between changes in groundwater levels and

要があろう．もし夜間電力使用量との関係があるならsurface deformation from GPS measurements inば，より顕著な変化が現れると予想される．更に，他

the San Gabriel Valley, California. Geophys. Res.

地点の地磁気データや地電位差データでも同様の解析

Lett., 39, L01301.

を行い，解析結果の地域差から統計パラメータに現れKing, C.−Y., S. Azuma, G. Igarashi, M. Ohno, H. Saito

る季節変動の原因や地下水との具体的な関係が議論でand H. Wakita（1999）： Earthquake−related water

きるようになるだろう．

−level changes at 16 closely clustered wells in謝辞

Tono, central Japan. J. Geophys. Res., 104, 13073−
13082.

本研究の一部は，文部科学省による「地震及び火山

Kümpel, H.−J.（1992）： About the potential of wells

噴火予知のための観測研究計画」（課題番号：2501）おto reﬂect stress variations within inhomogeneous

よび東海大学海洋研究所による「個別プロジェクト研

crust. Techtonophys., 211, 317−336.

究，東海地方における地震・火山・津波災害軽減のLiu, J.−Y., C.−H. Chen, C.−H. Lin, H.−F. Tsai, C.−ための学際的研究」の支援を受けて実施された．また，H. Chen and M. Kamogawa（2011）：Ionospheric

気象庁柿岡地磁気観測所および防災科学技術研究所

disturbances triggered by the 11 March 2011より地震データ，K 指数データおよび気象データを

M9.0 Tohoku earthquake. Geophys. Res. Lett., 116,

提供して頂いた．さらに，京都大学大学院理学研究

A06319.

科附属地磁気世界資料解析センターおよびドイツの

Maeda, T., T. Furumura, S. Sakai and S. Shinohara

GeoForschungsZentrum Potsdam よ り Dst 指 数 デ ー（2011）：Signiﬁcant tsunami observed at ocean−タおよび ap 指数データを提供して頂いた．そして，

bottom pressure gauges during the 2011 off the東北地震前後の奥多摩地磁気パラメータ変化

Paciﬁc coast of Tohoku Earthquake. Earth Planets
Space, 63, 803−808.

21

Ujihara, N., Y. Honkura and Y. Ogawa（2004）：Electric
and magnetic ﬁeld variations arising from the

Nagao, T., Y. Orihara, T. Yamaguchi, I. Takahashi, K.

seismic dynamo effect for aftershocks of the M7.1

Hattori, Y. Noda, K. Sayanagi and S. Uyeda（2000）：

earthquake of 26 May 2003 off Miyagi Prefecture,

Co−seismic geoelectric potential changes observedin Japan. Geophys. Res. Lett., 27, 1535−1538.

NE Japan. Earth Planets Space, 56, 115−123.
Utada, H., H. Shimizu, T. Ogawa, T. Maeda, T.

Okada, T., K. Yoshida, S. Ueki, J. Nakajima, N.

Furumura, T. Yamamoto, N. Yamazaki, Y.

Uchida, T. Matsuzawa, N. Umino, A. Hasegawa

Yoshitake and N. Nagamachi（2011）：Geomagnetic

and Group for the aftershock observations of theﬁeld changes in response to the 2011 off the2011 off the Paciﬁc coast of Tohoku Earthqauke

Paciﬁc Coast of Tohoku Earthquake and Tsunami.

（2011）：Shallow inland earthquakes in NE Japan

Earth Plant. Science Lett., 311, 11−27.

possibly triggered by the 2011 off the Paciﬁc coast

Uyeda, S., M. Hayakawa, T. Nagao, O. Molchanov, K.

of Tohoku Earthquake. Earth Planets Space, 63,

Hattori, Y. Orihara, K. Gotoh, Y. Akinaga and H.

749−754.

Tanaka（2002）：Electric and magnetic phenomena

Roeloffs, E.A.（1998）：Persistent water level changes

observed before the volcano−seismic activity inin a well near Parkﬁeld, California, due to local and2000 in the Izu Island Region, Japan. Proc. Nat.

distant earthquakes. J. Geophys. Res. 103, 869−889.

Acad. Sci., USA, 99, 7352−7355.

Takahashi, H. （2011）：Static strain and stress

Yukutake, Y., R. Honda, M. Harada, T. Aketagawa, H.

changes in eastern Japan due to the 2011 off the

Ito and A. Yoshida（2011）： Remotely−triggered

Paciﬁc coast of Tohoku Earthquake, as derived

seismicity in the Hakone volcano following thefrom GPS data. Earth Planets Space, 63, 741−744.

2011 off the Paciﬁc coast of Tohoku Earthquake.

Tsugawa, T., A. Saito, Y. Otsuka, M. Nishioka, T.

Earth Planets Space, 63, 737−740.

Maruyama, H. Kato, T. Nagatsuma and T. Murata

Widarto, D.S., T. Mogi, Y. Tanaka, T. Nagao,（2011）：Ionospheric disturbances detected by GPS

K. Hattori and S. Uyeda （2009）：Co−seismic

total electron content observation after the 2011

geoelectrical potential changes associated with theoff the Paciﬁc coast of Tohoku Earthquake. Earth

June 4, 2000 s earthquake（Mw7. 9）in Bengkulu,

Planets Space, 63, 875−879.

Indonesia. Phys. Chem. Earth, 34, 373−379.