４－２

2014 年 1 月房総半島沖スロースリップイベント
The Slow Slip Event off the Boso Peninsula in January, 2014
防災科学技術研究所
National Research Institute for Earth Science and Disaster Prevention2013 年 12 月 31 日より房総半島沖で地震が群発的に発生した．主な地震の発震機構解はフィリ
ピン海プレートと日本島弧との相対運動方向にすべりベクトルを有する低角逆断層型である(第 1図，発震機構解は防災科研高感度地震観測網 Hi-net および AQUA(松村ほか, 2006)1)による)．また，これと同期して房総半島の観測点で傾斜変動が観測された．勝浦東(KT2H)観測点では，北西下がりの傾斜が約 0.4μrad に達した．
房総半島沖ではスロースリップイベントが繰り返し発生することが知られている．房総半島沖ス
ロースリップイベントは群発地震を伴うことが特徴であり，測地データおよび地震活動から，これまでに第 2 図に示した時期の発生が知られている(防災科研, 20032); Hirose et al., 20123))．今回の群発地震は，過去の活動とほぼ同じ場所で発生した(第 1 図)．一方，今回のスロースリップイベントは，過去約 30 年間で最短の 2 年 2 ヶ月で再来した．房総半島沖における応力の蓄積状況をモニタリングするために，スロースリップイベントのすべり履歴を明らかにすることは重要である．そこで，詳細震源分布を決定するとともに，傾斜データより断層モデルの推定を行い，過去の活動と比較した．
第 3 図に示した領域の 30km 以浅で 2005 年 1 月 1 日から 2014 年 2 月 1 日までの期間に発生した地震を対象として，波形相関を用いた Double Difference 法 4)によって詳細震源分布を決定した．初期震源は防災科研 Hi-net による．ただし，一部自動処理結果を含む．再決定結果を第 3 図，時空間分布を第 4 図に示す．今回の活動は，解析期間中に地震が発生した領域(地震発生域とする)の北端周辺で発生した．2014 年 1 月 2 日頃から一宮沖で地震が発生しはじめ，翌日から西に発生域が拡大し，5 日から地震発生域の南端でも地震が少数発生した．
第 5 図に，房総半島の観測点での傾斜記録を日別地震発生数とともに示す．2013 年 12 月 31 日頃から勝浦東 (KT2H) 観測点等で変動が現れはじめ，2014 年 1 月 2 日頃から明瞭になった．群発地
震は傾斜変動と同期して 1 月 2 日から発生数が増加した．2013 年 12 月 31 日から 2014 年 1 月 6 日までの傾斜変化を元に矩形断層モデルを推定したところ，第 6 図に示す解が得られた．スロースリップイベントの規模は Mw6.1 と推定された．
前回の 2011 年および前々回の 2007 年のスロースリップイベントに伴う群発地震の活動を第 7 図に示す．今回の活動と比較すると，以下の特徴が見られる．群発地震がはじめ東の海域で発生し，次いで西に移動する点は 2007 年，2011 年とも共通しており，また 2007 年および 2011 年の期間(I)(はじめの 2 日半)では主に地震発生域の北端周辺で地震が発生した点も今回と類似している．今回推定されたすべり域の位置は，傾斜データから推定された 2007 年および 2011 年の期間(I)のすべり域(防災科研, 20085)の第 5 図; Hirose et al., 20123) の Fig.S1(1)を参照) とほぼ重なる．KT2H の傾
斜変動量は 2007 年(Mw6.4，北西方向に約 0.7μrad)の約 1/2 で，Mw6.1 と規模が小さいことと調和的である．また，2011 年の期間(I)の傾斜変動(西北西方向に約 0.3μradian)と方向・変動量とも近く，この期間の規模(Mw6.2)とほぼ同じである．一方，2011 年は，期間(II)以降に勝浦直下でも－ 108 －地震が多数発生するとともに，KT2H で北北西方向に約 1.0μradian に達する大きな変動が見られ，勝浦沖および勝浦直下でもすべり域が推定されている点(Hirose et al., 20123)の Fig.S1(2),(3)を参照)が今回と異なる．今回のすべり域は，2007 年および 2011 年の期間(I)とほぼ重なるが，全体の規模は小さく，そのために過去約 30 年間で最短の 2 年 2 ヶ月で再来した可能性が考えられる．今後，さらに詳細なすべり過程を明らかにする必要がある．
謝辞：

気象庁のホームページで公開されている気象データを使用させていただきました．記して感謝いたします．
(木村尚紀)
Hisanori Kimura
参考文献
1) 松村稔・伊藤善宏・木村尚紀・小原一成・関口渉次・堀貞喜・笠原敬司: 高精度即時震源パラメータ解析システム(AQUA)の開発, 地震 2, 59, 167-184, 2006.
2) 防災科研, 地殻傾斜の連続観測で捉えた 2002 年 10 月に発生した房総半島東方沖のスロースリップ, 連絡会報, 69, 198-204, 2003.
3)

Hitoshi Hirose, Hisanori Kimura, Bogdan Enescu, and Shin Aoi, 2012, Recurrent slow slip event
likely hastened by the 2011 Tohoku earthquake, Proc. Natl. Acad. Sci. USA, 109, 15157-15161.

4) Waldhauser, F. and Ellsworth, W. L.: A Double-Difference Earthquake Location Algorithm: Methodand Application to the Northern Hayward Fault, California, Bull. Seismol. Soc. Am., 1353-1368,2000.
5) 防災科研, 2007 年 8 月に発生した房総沖スロースリップと群発地震活動, 連絡会報, 79,
123-127, 2008.
6) Tamura, Y., T. Sato, M. Ooe, M. Ishiguro, A procedure for tidal analysis with a Bayesian
information criterion, Geophys. J. Int., 104, 507-516, 1991.

－ 109 －140˚

140.5˚141˚

2. 2014/01/02 22:11

H=22.2km M4.8 H=32km Mw4.9

4. 2014/01/22 20:18

H=24.1km M3.5

35.5˚
Hi-netHi-netAQUA-CMT

1. 2014/01/02 16:48

H=21.3km M4.1

3. 2014/01/07 22:11

H=26.4km M3.9 H=24km Mw4.0

H=30km Mw4.2

35˚
Hi-net

第1図Fig. 1

AQUA-CMT

30 km

Hi-netAQUA-CMT

房総半島沖の房総スロースリップイベントに伴う震央分布．震源要素は 2002 年 9 月 1 日
以降は防災科研 Hi-net，それ以前は防災科研 関東東海観測網による．赤丸は 2013 年 12 月31 日以降の地震を示す(自動処理結果を含む)．主な地震の防災科研 Hi-net およびAQUA 1)による発震機構解をあわせて示す．
The epicentral distribution associated with the Slow Slip Event (SSE) off the Boso Peninsula.
Hypocentral parameters from NIED Hi-net and NIED Kanto-Tokai network are used forearthquakes after and before September, 1, 2002, respectively. Red symbols denote earthquakes
after December 31, 2013 (automatically determined hypocenters are included). Focal
mechanisms determined by NIED Hi-net and AQUA1) are shown for major earthquakes.

－ 110 －6M5

1987年
千葉県東方沖地震

2011年
東北地方太平洋沖地震

4321

1985

1990

1995

2000

2005

2010

6M

543212229 01
2013/1208

15222014/01

第2図29 0108
2014/02

房総半島沖の M-T 図． 第１図の矩形領域について示す．房総半島沖スロースリップイベ
ントを矢印で，スロースリップイベントに同期した地震活動を色付きシンボルで表す．
Fig. 2 The M-T diagram for earthquakes off the Boso Peninsula. Earthquakes in the rectangle in Fig.
1 are shown. Arrow denotes the Boso SSE, and colored symbol shows earthquake associatedwith the Boso SSE.

－ 111 －

35.6˚

431AD35.2˚
CB

34.8˚30 km

205101520253035

Depth(km)20131228 − 20140127 No=108140˚
140.4˚ 140.8˚M
5.0
4.0
3.02.0波形相関データを用いた Double Difference 法 4)による高精度震央分布．再決定処理は 2005
年 1 月 1 日から 2014 年 2 月 1 日までの期間に左図の青線で示した領域の 30km 以浅で発生した地震を対象として行った．初期震源として，防災科研 Hi-net を用いた(一部自動処理結果を含む)．灰丸は解析期間の震央分布を示す．2014 年 1 月房総半島沖スロースリップイベントに伴う群発地震を色付きシンボルで示す．左図の黒四角は解析に使用した観測点を示す．
Fig. 3 The high-precision epicentral distribution determined by Double Difference method4) byincorporating waveform cross corelation data. Earthquakes in the blue rectangle in the left panel
shallower than 30 km from January 1, 2005 to February 1, 2014 are relocated. Hypocentral
parameters from NIED Hi-net are used as initial hypocenters (automatically determined
hypocenters are included). Gray symbols denote earthquakes during the analyzing period.
Colored symbols denote earthquakes associated with the Boso SSE in January, 2014. Solid squarein the left panel shows the seismographic station used for relocation.
第3図－ 112 －Y(km)

20 A

3101

06

2014−01
11160−20X(km)−40 B−20C020D

第 4 図 第 3 図 A-B および C-D に投影した時空間分布．
Fig. 4 The space-time plot projected along the A-B and C-D in fig. 3.

－ 113 －2126

CH2H NCH2H EHSNH NHSNH ECBAH N
CBAH EKT2H NKT2H EN, E down
0.1 µ radKTUH NKTUH E

Earthquakes [count/day]

40200100

HA2H atm. pres. [hPa]
15 hPaKatsuura precip. [mm]500

第5図Fig. 51522Dec 2013

2905

12Jan 2014

192013 年 12 月 10 日から 2014 年 1 月 22 日にかけての房総半島の観測点での傾斜時系列．
観測点位置は第 6 図に示した. 上方向への変化が北・東下がりの傾斜変動を表す.
BAYTAP-G 6)による潮汐・気圧成分補正，リニアトレンド除去後の記録を示した. Hi-net
波崎２観測点(HA2H) の気圧観測値を使用した. 灰矢印で示した時刻のステップを補正した．この期間の房総半島沖(第１図矩形領域)の日別地震数, HA2H の気圧, 勝浦での雨量をあわせて表示した.
The time series of tilt records at the stations around the Boso peninsula from December 10,
2013 to .January 22, 2014. Station locations are shown in Fig. 6. `N' and `E' that follow
a four-character station code denote the northward and eastward ground down tilt components,respectively. The displayed tilt records are detided and their atmospheric pressure responses arecorrected with BAYTAP-G6). In analysis, atmospheric pressure change observed at the Hi-net
HA2H station is used. Their linear trends are removed and steps of time-series data accompaniedwith ground motion due to earthquakes are also corrected (gray arrows). The daily numbersof earthquakes around the Boso peninsula (rectangle in Fig. 1), atmospheric pressure change
at HA2H, and daily precipitation at Katsuura are also shown.

－ 114 －lat. 35.21
lon. 140.77
strike 259 °dip 10 °HA2H
0 10 20 30
depth 19 km slip 8.8 cm
leng. 34 km wid. 14 km
CH2HNRTHM0 1.7e+18 Nm MW 6.1
rake 113 °HSNH
CBAH36˚N

35.6˚N35.2˚NKTUH

KT2HKMSH
34.8˚N

SlipTilt
0.3 µ rad10.0 cmObs.
Calc.

34.4˚N139.6˚E

第6図Fig. 6140˚E

140.4˚E

140.8˚E

141.2˚E観測された傾斜変化ベクトル(青矢印)・このデータから推定されたスロースリップの断層モデル(赤矩形・矢印)・モデルから計算される傾斜変化ベクトル(白抜き矢印)．期間中に発生した地震の震央分布を橙色の円で示す．
Observed tilt change vectors (blue arrows), the estimated fault slip (red arrow) and rectangularfault location and geometry (pink rectangle) based on the tilt change vectors, and the calculated
tilt changes due to the fault model (open arrows). Orange circles show epicenters of theearthquake swarm occurred in the same time period.

－ 115 －(A) 20070811 − 20070910 N=3330A

10D

20
30

35.2˚
CB

34.8˚

30 km

M
5.0
4.0
3.0
2.0Y(km)

35.6˚

35.6˚

−20−2006

2011−1111

160−20−40 B

B

30 km

X(km)

CFig. 7

2011−10263101

20 A

35.2˚

第7図D21D34.8˚C020140.8˚A2007−093101
06

−40 B

Y(km)

140.4˚

2007−0821260(B) 20111021 − 20111120 N=123140˚

1620 A

140.8˚X(km)140.4˚

Depth(km)140˚

11−20C020D(I)

(II)

(III)

(A)2007 年 8 月および(B)2011 年 10 月の房総スロースリップイベントに伴う群発地震の震
央分布および時空間分布．図の見方は第 3，4 図と同じ．
The epicentral distribution and the space-time plot for earthquake swarm associated with the Boso
SSE in (A) August, 2007 and (B) October, 2011. Nomenclatures are the same as those inFigs. 3 and 4.

－ 116 －