６−８ 愛知県東部の低周波地震に関連した産総研豊橋・豊橋東観測点の歪変化
Crustal strain changes at TYH and TYE station of AIST associated with lowfrequency earthquakes in eastern Aichi Prefecture

産業技術総合研究所
Geological Survey of Japan, AIST
１．はじめに2005 年 7 月 20 日から 23 日にかけて，愛知県で低周波地震活動が活発化した．それと同じ時期に，愛知県東部から静岡県西部の気象庁・防災科研・産総研の観測点（第１図）で，歪・傾斜に変化が観測された．気象庁の解析によれば，プレート境界（深さ約 30km）において，第１図の赤丸で示した地点のいずれかで，M5.8 相当のゆっくりとしたすべりが発生したとすれば，歪計などに見られる変化をほぼ説明できるとのことである．
愛知県豊橋市にある産総研豊橋観測点（1998 年 4 月観測開始）と豊橋東観測点（2004 年 9 月観測開始）では，ボアホール型多成分歪計による地殻歪に変化が認められた（第２図）．その歪変化の特徴をまとめたので報告する．また，過去にも同様の歪変化が発生していたかどうかの調査を行

２．2005 年 7 月の事例の特徴
豊橋観測点では，2005 年 7 月 20 日から 23 日にかけて，歪３(N342E) 成分に約 5 × 10-8 の縮みの変化が観測された（第２図）．
豊橋東観測点では，2005 年 7 月 20 日から 23 日にかけて，N356E 成分の水平歪で約 4 ～ 7 × 10-8の縮み，N131E 成分の水平歪で約 2 × 10-8 の縮み，N221E 成分の水平歪で約 1 × 10-8 の伸びが認められた（第２図）．この地殻歪変化は南北に縮み，東西に伸びの主歪場であった（第３図）．
気象庁のゆっくりすべりの断層モデル 1) に基づいて計算したところ，豊橋付近では南北に縮みの主歪場が推定された（第４図）．この結果は，観測結果とほぼ一致する．
３．過去の歪変化の調査
2005 年 7 月と同様の歪変化が過去に発生していたかどうかの調査を行った．調査は，1999 年 9月以降に愛知県の低周波地震と同期して気象庁蒲郡・佐久間・浜北観測点で歪変化が見られた 18例についてである．対象は愛知県豊橋市の産総研豊橋・豊橋東観測点の歪計の記録である．
2005 年 7 月の事例では，豊橋市付近では南北に縮みの歪場が観測されたので，18 例について，豊橋観測点の歪３(N342E 方向の成分 ) に縮みの変化が見られたかどうかを調べた（第１表）．また，最近の３例については，豊橋東東観測点の歪 (N356E 方向の成分 ) に縮み変化が見られたかどうか調べた（第２表）．
４．調査結果（１）低周波地震が 2005 年 7 月と同じ領域で発生した時には，豊橋観測点の歪の南北方向に縮みが観測された．また，縮みに先行して南北方向に伸びが認められる例もある．

－410－－410－

ったので，その結果を報告する．

（２）低周波地震が 2005 年 7 月の時よりも北東側の領域で発生した時には，豊橋観測点の歪に変化は（基本的には）観測されなかった．
５．まとめ2005 年 7 月の愛知県での低周波地震活動の活発化と同時期に，産総研豊橋観測点と豊橋東観測点での地殻歪に変化が見られることが分かった．この地殻歪変化はプレート境界（深さ約 30km）におけるゆっくりとしたすべりの断層モデルで説明することができる．豊橋観測点と豊橋東観測点での地殻歪変化は，低周波地震活動と同期してプレート境界でゆっくりすべりが発生したことを支持する．また，1999 年 9 月以降にも 10 回以上，低周波地震の活発化と同期した地殻歪変化が発生（つまりプレート境界でゆっくりすべりが発生）したと推測される．
（北川有一・松本則夫・小泉尚嗣）
参 考 文 献1) 小林昭夫・山本剛靖・中村浩二・吉川澄夫・橋本徹夫・竹中潤・木村一洋・青木玲子・露木貴裕・菅沼一成・高野淳，東海地域で歪計により観測された短期的スロースリップ，2005 年測地学会
秋季大会予稿集，2005
1999 年 9 月以降の豊橋観測点での歪変化
タイプA：2005/ 7/20- 7/22 の時の領域に低周波地震が発生
気象庁蒲郡観測点の歪に変化あり（第４図の断層モデル）
タイプB：タイプ A より北東側の領域で低周波地震が発生
気象庁蒲郡観測点の歪に変化なし（第５図の断層モデル）
Table.1 Crustal strain changes at TYH station since September 1999.
Type A:low-frequency earthquakes occurred in the area where low-frequency earthquakes occurredon July 20-22, 2005. There were strain changes at Gamagouri station of JMA. The fault
model of slow slip is shown in Fig.4.
Type B:low-frequency earthquakes occurred in a northeast area of Type A. There were no strainchanges at Gamagouri station of JMA. The fault model of slow slip is shown in Fig.5.

気象庁の歪変化の期間

タイプ

1999/ 9/ 3- 9/ 9
2000/ 2/26- 2/27
2001/ 4/10- 4/13
2001/ 9/ 6- 9/ 9
2002/ 9/22- 9/24
2002/11/26-11/27
2003/ 2/14- 2/16
2003/ 5/30- 6/ 3
2003/ 9/16- 9/17
2003/10/ 9-10/11
2003/11/ 6-11/ 8
2004/ 1/ 1- 1/ 3
2004/ 2/13- 2/17
2004/ 4/ 3- 4/ 7
2004/ 6/26- 7/ 2
2004/11/21-11/26
2004/12/17-12/19
2005/ 7/20- 7/22

A
A
A
B
B
A
B
A
B
B
A
B
A
B
A
B
A
A

豊橋歪３の縮み変化に先行する南北伸びの有無
1999/ 9/ 2- 9/10
あり (1999/ 8/29- 9/ 1)
2000/ 2/27- 2/28
あり (2000/ 2/23- 2/26)
2001/ 4/ 6- 4/13
なし
2001/ 9/10( 大雨の影響？ ) なし
2001/ 9/22- 9/24
あり (2001/ 9/19- 9/21)
2002/11/26-11/27
あり (2002/11/18-11/25)
2003/ 2/13( ステップ状変化 ) なし
2003/ 5/31- 6/ 2なし
なし
なし
2003/11/ 5-11/ 7
あり (2003/10/30-11/ 4)
なし
2004/ 2/15- 2/17
なし
なし
2004/ 6/20- 6/29
なし
なし
2004/12/18-12/22なし
2005/ 7/20- 7/23
なし
豊橋歪３の縮み変化の有無

－411－

第６図
第７図
第８図
第９図
第 10 図
第 11 図第 12 図
第 13 図
第 14 図
第 15 図
第 16 図
第 17 図
第 18 図
第 19 図第 20 図
第 21 図
第 22 図
第２図

－411－

第1表

第 2 表 2004 年 11 月以降の豊橋東観測点での歪変化
Table.2 Crustal strain changes at TYE station since November 2004.

豊橋・豊橋東豊橋東歪 N356E 成分の縮気象庁の歪変化の
豊 橋 東 歪 N356E 成 分タイプ
み変化に先行する南北伸期間の縮み変化の有無
びの有無
2004/11/21-11/26
2004/12/17-12/19
2005/ 7/20- 7/22

BAAなし
2004/12/14-12/19
2005/ 7/20- 7/23

歪（時間値）( 2005/07/01 00:00 - 2005/08/10 00:00 )豊橋1歪１(N102E)

豊橋1歪２(N222E)

豊橋1歪３(N342E)

�������

第 21 図なし
第 22 図あり (2005/ 7/16- 7/19) 第２図

��������������
豊橋東 歪(N356E)

����

－412－

���
先行する伸び？

●●●●豊橋 豊橋東●●����●

豊橋東歪(N131E)

豊橋東歪(N221E)

���●
●●
●
●●

����
���

���

豊橋・豊橋東観測点の位置及び 2005 年 7 月の低周波地震の震央と推定されたすべりの位置（気象庁資料に加筆した）．緑色の丸は産総研の観測点，それ以外は気象庁の観測点．
Fig.1 Location of TYH and TYE stations, epicenter of low-frequency earthquakes inJuly 2005, and estimated area of slow slip. The figure is added with figure byJMA. Green circles show stations of AIST and the others show stations of JMA.

����������������

���

第2図－412－

第2図

産総研豊橋・豊橋東観測点での歪の結果．BAYTAP-G により潮汐応答・気圧応答を除去し，さらに直線トレンドを除去している．
Fig.2 The results of strain at TYH and TYE. These results are derived from the raw databy the removal of tidal and barometric responses estimated by BAYTAP-G and thelinear trends.

5.0E-8/1cm
伸び 縮み

豊橋

第 3 図 豊橋東観測点で観測された歪変化から求めた主歪Fig.3 Principal strain estimated from observed strain changes at TYE.
第5図－413－

第 4 図とは別タイプのゆっくりすべりの断層モデル 1) と断層モデルにより推定される豊橋東観測点での主歪
Fig.5 The fault model of slow slip different from Fig.4 1) and principal strain at TYE
estimated from the fault model.
豊橋（時間値）( 1999/08/21 00:00 - 1999/09/20 00:00 )豊橋1 歪１(N102E)

������
豊橋1 歪２(N222E)

������
豊橋1 歪３(N342E)

����豊橋���
先行する伸び？

����第6図

����

���

����

����

豊橋観測点の歪の結果．BAYTAP-G により潮汐応答・気圧応答を除去した
後，直線トレンドを除去している．
Fig.6 The results of strain at TYH. These results are derived from the raw data by theremoval of tidal and barometric responses estimated by BAYTAP-G and the linear
trends.

第4図2005 年 7 月のゆっくりすべりの断層モデル 1) と断層モデルにより推定される豊橋東観測点での主歪
Fig.4 The fault model of slow slip in July 2005 1) and principal strain at TYE estimatedfrom the fault model.

－413－豊橋（時間値）

豊橋（時間値）( 2000/02/12 00:00 - 2000/03/13 00:00 )

豊橋1

10-9

( 2001/08/26 00:00 - 2001/09/25 00:00 )歪１(N102E)

豊橋1

歪２(N222E)

豊橋1 歪２(N222E)

10-980

50
豊橋1

10-9

歪１(N102E)

1008010-9豊橋110-9歪３(N342E)

豊橋110-9歪３(N342E)

40800
先行する伸び？

2/122/20

2/27

3/5

8/26

3/12

第7図

9/16

9/23

－414－

豊橋観測点の歪の結果．BAYTAP-G により潮汐応答・気圧応答を除去した
後，直線トレンドを除去している．
Fig.9 The results of strain at TYH. These results are derived from the raw data by theremoval of tidal and barometric responses estimated by BAYTAP-G and the linear
trends.

豊橋（時間値）

豊橋（時間値）( 2001/03/26 00:00 - 2001/04/25 00:00 )
10-9

9/9

第9図

豊橋観測点の歪の結果．BAYTAP-G により潮汐応答・気圧応答を除去した
後，直線トレンドを除去している．
Fig.7 The results of strain at TYH. These results are derived from the raw data by theremoval of tidal and barometric responses estimated by BAYTAP-G and the linear
trends.

豊橋19/2( 2002/09/08 00:00 - 2002/10/08 00:00 )歪１(N102E)

豊橋1歪１(N102E)

豊橋1歪２(N222E)

豊橋1歪３(N342E)

10-98050
豊橋1 歪２(N222E)

10-910-9200400
豊橋1 歪３(N342E)

10-910-912080

先行する伸び？

第8図

3/26

4/1

4/8

4/15

4/22

9/8

第 10 図

豊橋観測点の歪の結果．BAYTAP-G により潮汐応答・気圧応答を除去した
後，直線トレンドを除去している．
Fig.8 The results of strain at TYH. These results are derived from the raw data by theremoval of tidal and barometric responses estimated by BAYTAP-G and the lineartrends.

9/15

9/22

9/29

10/6

豊橋観測点の歪の結果．BAYTAP-G により潮汐応答・気圧応答を除去した後，直線トレンドを除去している．
Fig.10 The results of strain at TYH. These results are derived from the raw data by theremoval of tidal and barometric responses estimated by BAYTAP-G and thelinear trends.

－414－豊橋（時間値）

豊橋（時間値）( 2003/05/14 00:00 - 2003/06/15 00:00 )( 2002/11/11 00:00 - 2002/12/11 00:00 )
豊橋1

10-9

歪１(N102E)

豊橋110-950歪１(N102E)

400
豊橋1 歪２(N222E)

10-9

豊橋1 歪２(N222E)

10-9150200
豊橋110-9歪３(N342E)

豊橋110-940歪３(N342E)

80
先行する伸び？

11/11

11/17

11/24

12/1

5/14

12/8

第 11 図

5/25

6/1

6/8

6/15

第 13 図－415－

豊橋観測点の歪の結果．BAYTAP-G により潮汐応答・気圧応答を除去した後，直線トレンドを除去している．
Fig.11 The results of strain at TYH. These results are derived from the raw data by theremoval of tidal and barometric responses estimated by BAYTAP-G and thelinear trends.

豊橋観測点の歪の結果．BAYTAP-G により潮汐応答・気圧応答を除去した後，直線トレンドを除去している．
Fig.13 The results of strain at TYH. These results are derived from the raw data by theremoval of tidal and barometric responses estimated by BAYTAP-G and thelinear trends.

豊橋（時間値）

豊橋（時間値）( 2003/02/01 00:00 - 2003/03/03 00:00 )( 2003/09/01 00:00 - 2003/10/01 00:00 )

豊橋1

10-9

5/18

歪１(N102E)

10-960豊橋1

歪１(N102E)

豊橋1歪２(N222E)

豊橋1

歪３(N342E)

250
豊橋1

10-9

歪２(N222E)

10-910060
豊橋1 歪３(N342E)

10-9

10-980100

2/1

2/92/16

2/23

3/2

9/1

第 12 図第 14 図

豊橋観測点の歪の結果．BAYTAP-G により潮汐応答・気圧応答を除去した後，直線トレンドを除去している．
Fig.12 The results of strain at TYH. These results are derived from the raw data by theremoval of tidal and barometric responses estimated by BAYTAP-G and thelinear trends.

9/7

9/14

9/21

9/28

豊橋観測点の歪の結果．BAYTAP-G により潮汐応答・気圧応答を除去した後，直線トレンドを除去している．
Fig.14 The results of strain at TYH. These results are derived from the raw data by theremoval of tidal and barometric responses estimated by BAYTAP-G and thelinear trends.

－415－豊橋（時間値）

豊橋（時間値）( 2003/09/27 00:00 - 2003/10/27 00:00 )( 2003/12/16 00:00 - 2004/01/15 00:00 )豊橋1

10-9歪１(N102E)

豊橋1 歪２(N222E)

豊橋1 歪２(N222E)

10-912040
豊橋110-9歪３(N342E)

豊橋1

10-9

9/27

10/5

10/12

10/19

12/16

10/26

第 15 図1/4

1/11－416－

豊橋（時間値）( 2003/10/23 00:00 - 2003/11/22 00:00 )( 2004/02/01 00:00 - 2004/03/02 00:00 )歪１(N102E)

10-92500

豊橋1歪１(N102E)

豊橋1歪２(N222E)

豊橋1歪３(N342E)

250
豊橋1 歪２(N222E)

10-950250
豊橋1 歪３(N342E)

10-912/28

豊橋観測点の歪の結果．BAYTAP-G により潮汐応答・気圧応答を除去した後，直線トレンドを除去している．
Fig.17 The results of strain at TYH. These results are derived from the raw data by theremoval of tidal and barometric responses estimated by BAYTAP-G and thelinear trends.

豊橋（時間値）
豊橋1

12/21

第 17 図

豊橋観測点の歪の結果．BAYTAP-G により潮汐応答・気圧応答を除去した後，直線トレンドを除去している．
Fig.15 The results of strain at TYH. These results are derived from the raw data by theremoval of tidal and barometric responses estimated by BAYTAP-G and thelinear trends.

10-9

歪３(N342E)

4060

10-9

歪１(N102E)

1500

120010-9豊橋110-910-980

150
先行する伸び？

10/23 10/26

11/2

11/9

11/16

2/1

第 18 図2/82/15

2/22

2/29

豊橋観測点の歪の結果．BAYTAP-G により潮汐応答・気圧応答を除去した後，直線トレンドを除去している．
Fig.18 The results of strain at TYH. These results are derived from the raw data by theremoval of tidal and barometric responses estimated by BAYTAP-G and thelinear trends.

第 16 図

豊橋観測点の歪の結果．BAYTAP-G により潮汐応答・気圧応答を除去した後，直線トレンドを除去している．
Fig.16 The results of strain at TYH. These results are derived from the raw data by theremoval of tidal and barometric responses estimated by BAYTAP-G and thelinear trends.

－416－豊橋・豊橋東（時間値）

豊橋（時間値）( 2004/11/10 00:00 - 2004/12/10 00:00 )( 2004/03/20 00:00 - 2004/04/19 00:00 )豊橋1 歪１(N102E)

10-9350豊橋1

10-9歪１(N102E)
10-9

豊橋1歪２(N222E)

豊橋1歪３(N342E)

豊橋東

歪(N356E)

豊橋東歪(N131E)

豊橋東歪(N221E)

150

1000

10-9200

豊橋1 歪２(N222E)

10-9

10-9250

500

10-9

豊橋1

10-9200

歪３(N342E)

10-94080

3/20

3/28

4/4

4/11

4/18

－417－

豊橋観測点の歪の結果．BAYTAP-G により潮汐応答・気圧応答を除去した後，直線トレンドを除去している．
Fig.19 The results of strain at TYH. These results are derived from the raw data by theremoval of tidal and barometric responses estimated by BAYTAP-G and thelinear trends.

豊橋・豊橋東（時間値）

豊橋（時間値）( 2004/12/05 00:00 - 2005/01/04 00:00 )( 2004/06/01 00:00 - 2004/07/15 00:00 )
豊橋1

12/5

11/28

豊橋・豊橋東観測点の歪の結果．BAYTAP-G により潮汐応答・気圧応答を除去した．加えて，豊橋観測点については直線トレンドを除去している．
豊橋東観測点については二次関数トレンドを除去している．
Fig.21 The results of strain at TYH. These results are derived from the raw data by theremoval of tidal and barometric responses estimated by BAYTAP-G. In addition,these results are derived by the removal of the linear trends at TYH and thequadratic curve trends at TYE.

第 19 図10-911/21

11/14

11/10

第 21 図豊橋110-9歪１(N102E)250

歪１(N102E)

豊橋1 歪２(N222E)

10-9120

1000

10-9

豊橋110-9

豊橋1歪３(N342E)

豊橋東歪(N356E)

豊橋東歪(N131E)

豊橋東歪(N221E)

200

歪２(N222E)

10-9100200
10-9豊橋110-9歪３(N342E)

200

10-915080

12/5

6/1

6/6

6/13

6/20

6/27

7/4

第 22 図

12/12

12/19

12/26

1/2

豊橋・豊橋東観測点の歪の結果．BAYTAP-G により潮汐応答・気圧応答を除去した．加えて，豊橋観測点については直線トレンドを除去している．
豊橋東観測点については二次関数トレンドを除去している．
Fig.22 The results of strain at TYH. These results are derived from the raw data by theremoval of tidal and barometric responses estimated by BAYTAP-G. In addition,these results are derived by the removal of the linear trends at TYH and thequadratic curve trends at TYE.

7/11

第 20 図

豊橋観測点の歪の結果．BAYTAP-G により潮汐応答・気圧応答を除去した後，直線トレンドを除去している．
Fig.20 The results of strain at TYH. These results are derived from the raw data by theremoval of tidal and barometric responses estimated by BAYTAP-G and thelinear trends.

－417－