茨城県南部における高頻度地震活動のテクトニクス的因果解析:三重プレート干渉・スラブ内部構造・流体移動と強震動応答の定量的評価

  1. Executive Summary
  2. 1. 首都圏直下における複合プレート幾何学と多重沈み込み系
    1. 1.1 北米(陸側)・フィリピン海(PHS)・太平洋(PAC)プレートの3次元接触構造
    2. 1.2 フィリピン海プレート北東端(Subducted PHS Wedge)の空間的限界と舌状伸長
    3. 1.3 関東沈み込み帯における三重会合点近傍の応力場勾配
  3. 2. 震源メカニズム解と震源深さに基づく地震の多層的分節化(Zoning)
    1. 2.1 深度30〜50km:陸側プレート / PHS上面境界型地震
    2. 2.3 深度60〜80km:PAC上面・スラブ内地震
  4. 3. 流体運動と脱水脆化(Dehydration Embrittlement)による地震誘発機構
    1. 3.1 蛇紋岩化かんらん岩および含水鉱物の高圧脱水反応プロセス
    2. 3.2 隙間水圧(Pore Fluid Pressure)の上昇と実効応力の低下メカニズム
    3. 3.3 地震波トモグラフィーによる低速度・高V_p/V_s比異常域
  5. 4. 2011年東北地方太平洋沖地震(M9.0)に伴う余震活動とクーロン破砕応力変化(\Delta CFF)
    1. 4.1 巨視的大震災後の広域応力場再編
    2. 4.2 茨城県南部スラブ境界における \Delta CFF の定量的評価
    3. 4.3 誘発地震の過渡的定常化モデル
  6. 5. 堆積盆地(関東平野)における深部地盤構造と強震動増幅特性
    1. 5.1 千葉・茨城境界域におけるハーフグラーベン構造と基盤岩形状
    2. 5.2 堆積層(S波速度 V_s = 0.3 \sim 1.5\text{ km/s})内の二次波生成と長周期地震動のトランスファー関数
    3. 5.3 震度観測点における震度増幅因子の定量的分析
  7. 6. 実務的インプリケーション:構造物設計・防災・BCPへの適用
    1. 6.1 高密度地震動に耐え力学的マージンを確保するための構造物設計パラメータ
    2. 6.2 確率論的地震動予測(PSHA)に基づく極値リスク評価
    3. 6.3 リアルタイムモニタリングと事業継続計画(BCP)の高度化
  8. 結論と今後の研究課題

Executive Summary

茨城県南部は、日本列島において最も地震活動が活発な領域の一つであり、マグニチュード(M)4〜6クラスの中規模地震が定常的に高頻度発生している。本報告では、当該地域直下における陸のプレート(オホーツクプレート/北米プレート:NA)、フィリピン海プレート(PHS)、ならびに太平洋プレート(PAC)の3次元的多重沈み込み幾何学構造を軸に、震源メカニズムの多層的解剖、流体運動に伴う脱水脆化メカニズム、2011年東北地方太平洋沖地震(M_w 9.0)によるクーロン破砕応力変化(\Delta CFF)、および関東平野深部地盤構造による強震動増幅効果を統合的に解析・評価する。

1. 首都圏直下における複合プレート幾何学と多重沈み込み系

1.1 北米(陸側)・フィリピン海(PHS)・太平洋(PAC)プレートの3次元接触構造

茨城県南部の地下斜面においては、世界的に見ても極めて稀有な3つのプレートが上下方向に重なり合う多重沈み込み構造が形成されている。最上層には陸側のオホーツクプレート(NA)が存在し、その下位に相模トラフから北西方向へと沈み込むフィリピン海プレート(PHS)、さらに最下層には日本海溝から西〜西北西方向へと沈み込む太平洋プレート(PAC)が重層している。

【地下構造の概念的鉛直断面(NW-SE観測)】

標高 0km  ------------------------------------------- [陸側プレート (NA)]
          |  (関 東 平 野 深 部 堆 積 層)
深度20km  -------------------------------------------
          |  陸側プレート / PHS 境界面 (NW-SE圧縮逆断層)
深度40km  ------------------------------------------- [フィリピン海プレート (PHS)]
          |  PHS内部スラブ / PHS-PAC境界 (複合応力場)
深度60km  ------------------------------------------- [太平洋プレート (PAC)]
          |  PAC上面 / PAC内部スラブ (EW/E-W圧縮逆断層)
深度80km  -------------------------------------------

プレートの相応的な運動ベクトルは以下の通り定義される。\vec{V}_{\text{PAC}} \approx 8.0 \sim 9.0\text{ cm/year}\; (\text{方向: WNW})\vec{V}_{\text{PHS}} \approx 3.0 \sim 5.0\text{ cm/year}\; (\text{方向: NW})

これら運動ベクトルの差異により、相互の境界剪断力およびスラブ間における強い機械的干渉が定常的に発生する。

1.2 フィリピン海プレート北東端(Subducted PHS Wedge)の空間的限界と舌状伸長

地震波トモグラフィー解析および高精度震源再決定(Double-Difference法等)の結果、PHSの北東端は茨城県南部から千葉県北部にかけての地下深度 40〜60 km 付近で先端を形成し、消滅(あるいはPACへ完全収斂)していることが判明している。このPHSの楔形(ウェッジ状)先端部は、上方からの陸側プレートによる押圧と、下方からのPAC上面によるスラブ間力学的せん断を挟み込みの形で受けており、局部的な応力集中(Stress Concentration)の極大領域となっている。

1.3 関東沈み込み帯における三重会合点近傍の応力場勾配

房総沖に位置する相模トラフ・日本海溝の三重会合点(Triple Junction)から北西方向に連なる茨城県南部地域は、巨大地震発生帯の背後地応力場に位置する。主応力軸(\sigma_1)の配向は、浅部から深部へ移行するにつれてNW-SE方向からE-W方向へと連続的に回転する幾何学的勾配を有している。

2. 震源メカニズム解と震源深さに基づく地震の多層的分節化(Zoning)

茨城県南部の地震活動は、単一の断層運動ではなく、震源深さ依存の明確な3層構造に分節化される。 構造区分震源深さ範囲主なテクトニクス境界発震機構(メカニズム解)主応力軸(\sigma_1)方向震源の物理的特徴 第1層(浅部境界)30 km ~ 45 km陸側プレート(NA) / PHS上面逆断層型 (一部横ずれ)NW – SE 圧縮PHS上面のプレート境界すべり。高アスペリティ度。 第2層(中間干渉部)45 km ~ 60 kmPHS内部 / PHS-PAC相互境界逆断層型・複雑型NW-SE ~ WNW-ESE 圧縮PHS先端の挟み込みによる応力集中。 M5\sim6 頻発域。 第3層(深部スラブ)60 km ~ 80 kmPAC上面 / PACスラブ内逆断層型E – W 圧縮太平洋プレート内部屈曲応力および二重深発地震面上面。

2.1 深度30〜50km:陸側プレート / PHS上面境界型地震

本領域で発生する地震(例:2004年7月10日 M4.7 深さ48km、2012年7月16日 M4.8 深さ52km)は、北西ー南東方向に圧力軸を持つ逆断層型が支配的である。陸側プレートの基盤と沈み込むPHS上面のアスペリティ破壊に伴うものであり、プレート運動に伴う剪断応力の蓄積速度が高い。

### 2.2 深度50〜60km:PHSスラブ内およびPHS–PAC境界型地震 最も高頻度かつ M5.5 規模を超える中規模地震が集中する領域である。PHSスラブ下端とPACスラブ上端が物理的に直接接触・せん断作用を及ぼし合っており、PHSスラブの力学的破壊およびPACの押圧による内部変形が重畳して発生する。

2.3 深度60〜80km:PAC上面・スラブ内地震

東西方向に圧力軸を持つ逆断層型地震(例:2004年2月4日 M4.2 深さ65km、2011年5月20日 M4.6 深さ64km)が主である。太平洋プレートスラブが沈み込む過程で被る体積収縮・彎曲応力、ならびにスラブ内部の二次的脆性破壊が原因となる。

3. 流体運動と脱水脆化(Dehydration Embrittlement)による地震誘発機構

一般に地下40kmを超える高圧(\sim 1.5\text{ GPa})・高温環境下では、岩石は延性破壊挙動を示すため脆性破壊(地震)は発生しにくい。しかし茨城県南部では深度50〜70kmに広範な地震活動が存在する。この物理的矛盾を解明するキーが脱水脆化現象(Dehydration Embrittlement)である。

3.1 蛇紋岩化かんらん岩および含水鉱物の高圧脱水反応プロセス

沈み込む海洋プレート(PACおよびPHS)の表層海洋地殻(ガブロ・玄武岩)や上部マントルは、海洋底水熱変成作用により多量の含水鉱物(蛇紋岩 Serpentinite、ローソナイト Lawsonite、緑泥石 Chlorite、角閃石 Amphibole)を含有する。

これらの含水鉱物が沈み込みに伴う温度・圧力上昇曲線(P-T パス)に従って相変態を起こす際、結晶構造内の水(\text{H}_2\text{O})を分離・放出する。\t[span_2](start_span)ext{Serpentine (蛇紋岩)} \xrightarrow{\Delta T, \Delta P} \text{Forsterite (かんらん石)} + \text{Enstatite (輝石)} + \text{H}_2\text{O} (\text{脱水流体})

3.2 隙間水圧(Pore Fluid Pressure)の上昇と実効応力の低下メカニズム

放出された高温高圧の流体(超臨界水)は、プレート境界およびスラブ内部の微細な間隙に充填される。クーロン破壊基準(Coulomb Failure Criterion)は以下のように表される。\tau_[span_3](start_span)f = C + \mu' \sigma_n' = C + \mu ( \sigma_n - \mathcal{P}_f )

ここで: * [span_4](start_span)\tau[span_4](end_span)_f:せん断破壊強度

  • C:粘着力
  • \mu:内面摩擦係数
  • \mu’:実効摩擦係数
  • \sigma_n:全垂直応力
  • \sigma_[span_5](start_span)n’:実効垂直応力
  • \mat[span_5](end_span)hcal{P}_f:間隙流体圧(Pore Fluid Pressure)

脱水反応によって間隙流体圧 \mathcal{P}_f が全垂直応力 \sigma_n に極めて近いレベルまで急上昇(\mathcal{P}_f \to \sigma_n)すると、実効垂直応力 \sigma_n’ が極小化する。結果として、せん断応力 \tau の増加を伴うことなくせん断破壊強度 \tau_f が急激に低下し、微小な差応力下でも脆性滑り(地震発起)が誘発される。

3.3 地震波トモグラフィーによる低速度・高V_p/V_s比異常域

高密度地震観測網(Hi-net)の波形データを用いた3次元地震波速度トモグラフィー解析によれば、茨城県南部地下の深度40〜60km付近に明瞭な「低 V_p・低 V_s・高 V_p/V_s 比(\sim 1.85 以上)」領域が確認されている。高 V_p/V_s 比の存在は、地盤骨格内に高圧間隙流体(脱水流体)が充満している直接的なインシチュ(in-situ)証拠であり、この流体濃縮帯と高頻度震源分布の空間的整合性は極めて高い。

4. 2011年東北地方太平洋沖地震(M9.0)に伴う余震活動とクーロン破砕応力変化(\Delta CFF)

2011年3月11日の東北地方太平洋沖地震(M_w 9.0)の発生は、関東地方北部の応力場構造を根幹から変容させた。

4.1 巨視的大震災後の広域応力場再編

M_w 9.0 のメガトラスト破壊により、日本列島東部は巨大な東向きの変位(茨城県沿岸で最大数メートル)を受けた。これにより、太平洋プレートと陸側プレート間の固着状態および広域のせん断応力場が劇的に再編された。

4.2 茨城県南部スラブ境界における \Delta CFF の定量的評価

巨大地震が周縁断層の破壊限界に与えた影響を定量化するため、クーロン破砕応力変化(\Delta CFF)を以下のように評価する。\Delta[span_6](start_span) CFF = \Delta \tau + \mu' \Delta \sigma_n

ここで:

  • [span_6](end_span)\D[span_7](start_span)elta[span_7](end_span) \tau:受動断層面上のすべり方向のせん断応力変化(促進を正とする)
  • \Delta \sigma_n:断層面に対する法線応力変化(引張を正とする)
  • \mu’:実効摩擦係数(\approx 0.4)

\begin{aligned} \Delta CFF > 0 &\implies \text{断層の破壊を促進(地震誘発)} \\ \Delta CFF < 0 &\implies \text{断層の破壊を抑制(影領域)} \end{aligned}

数値シミュレーション解析によれば、東北地方太平洋沖地震の巨大な主破壊に伴い、茨城県南部の深度 40〜60 km に位置する「NW-SE圧縮逆断層」および「E-W圧縮逆断層」のいずれに対しても、\Delta CFF は +0.1\text{ MPa} \sim +0.5\text{ MPa} という顕著な正の値を記録した(一般に +0.01\text{ MPa} 以上で地震誘発効果が有意とされる)。

4.3 誘発地震の過渡的定常化モデル

大震災直後、茨城県南部の地震発生率は平常時の数倍〜10倍へと跳ね上がった。Rate-and-State 摩擦法則(Dieterichモデル)に基づく地震発生確率の過渡的応答解析式は以下の通りである。R(t) = \fr[span_8](start_span)ac{r}{1 + \left[ \exp\left( - \frac{\Delta CFF}{A \sigma \ln(10)} \right) - 1 \right] \exp\left( - \frac{t}{t_a} \right)}

ここで r は背景地震発生率、A \sigma は constitutive 定数と有効法線応力の積、t_a は緩和時間(Relaxation Time)である。茨城県南部においては、t_a が数十年に及ぶと試算されており、巨大地震から時間が経過した現在もなお、高い発生頻度が維持される定常的誘発状態が継続している。

5. 堆積盆地(関東平野)における深部地盤構造と強震動増幅特性

茨城県南部で地震が「体感されやすい」「震度が高く出やすい」理由は、発生源のメカニズムだけでなく、地表に至る伝播経路における地質構造的増幅効果にも依存している。

【関東平野深部地盤モデル(概念図)】

地表 (震度観測点)
  │  ▲ 沖積・洪積層 (Vs = 100 ~ 400 m/s) ── [表層増幅・1次共振]
  │
深度 500m ~ 2,000m
  │  ■ 新第三紀・第四紀 堆積層 (Vs = 500 ~ 1,500 m/s) ── [長周期地震動のトランスファー]
  │
深度 2,000m ~ 4,000m
  │  ◆ ハーフグラーベン構造(基盤岩の著しい起伏・凹み)
  │
地下基盤 (Vs ≧ 3,000 m/s) ───────────────── [地震波(S波)の入射]

5.1 千葉・茨城境界域におけるハーフグラーベン構造と基盤岩形状

関東平野の地下深部には、およそ1600万年前の日本海拡大に伴って形成されたハーフグラーベン(半落ち状の基盤凹曲構造)が複数伏在している。地震基盤(S波速度 V_s \ge 3.0\text{ km/s})の深さは、茨城県南部〜埼玉県東部にかけて 2,000\text{ m} \sim 4,000\text{ m} に達し、基盤上に未固結〜半固結の未固結堆積物(上総層群・新井宿層等)が厚く積層している。

5.2 堆積層(S波速度 V_s = 0.3 \sim 1.5\text{ km/s})内の二次波生成と長周期地震動のトランスファー関数

深部基盤と堆積層の間の顕著なインピーダンス比(Z = \rho \cdot V_s の不連続)により、下位から入射したS波は堆積層内で捕捉され、多重反射を起こす。さらに、ハーフグラーベンの不均一な傾斜境界において実質的な体積波(S波)から表面波(ラブ波・レイリー波)への波形変換が促進される。

地表と基盤のレスポンス比を示すトランスファー関数 H(\omega) は、一次元重複反射理論において次式で近似される。\vert{}H(\om[span_10](start_span)ega)\vert{} \approx \frac{1}{\sqrt{\cos^2\left(\frac{\omega H}{V_s}\right) + \left(\frac{\rho_1 V_{s1}}{\rho_2 V_{s2}}\right)^2 \sin^2\left(\frac{\omega H}{V_s}\right)}}

このインピーダンスコントラスト(\frac{\rho_2 V_{s2}}{\rho_1 V_{s1}} \gg 1)により、特定の固有周期(T \approx 1 \sim 3\text{ 秒})における地震動振幅が局所的に 3〜5倍以上に増幅 される。

5.3 震度観測点における震度増幅因子の定量的分析

地層区分S波速度 (V_s)密度 (\rho)主な地質構成増幅・応答特性 表層地盤(0–30m)100 \sim 300\text{ m/s}1.6 \sim 1.8\text{ g/cm}^3沖積粘土、沖積砂、鬼怒川・利根川氾濫原堆積物短周期(0.1–0.5秒)振幅の著しい増幅、非線形応答。 浅部堆積層(30–500m)300 \sim 800\text{ m/s}1.8 \sim 2.0\text{ g/cm}^3洪積段丘層、関東ローム層、上総層群上部やや短周期~中周期の定常波形成。 深部堆積層(0.5–3km)800 \sim 1,500\text{ m/s}2.0 \sim 2.3\text{ g/cm}^3上総層群下部、三浦層群やや長周期地震動(1–3秒)の生成・長周期共振。 地震基盤\ge 3,000\text{ m/s}\ge 2.6\text{ g/cm}^3先新第三紀三波川帯・足尾帯基盤岩類基準地震波の入力境界。

6. 実務的インプリケーション:構造物設計・防災・BCPへの適用

茨城県南部における高頻度地震環境は、インフラ・建築構造物の耐震設計および企業の事業継続計画(BCP)において独自の設計思想を要求する。

6.1 高密度地震動に耐え力学的マージンを確保するための構造物設計パラメータ

  • 疲労累積評価の義務化: 単一の巨大地震(L2 地震動)に対する一回限りの終局耐力検証のみならず、M5.5\sim6.5 規模の中規模地震の連続飛来(累積塑性ひずみ累積)に対する低サイクル疲労照査が不可欠である。
  • 長周期スペクトル適合: 固有周期 1.0〜3.0 秒に一次固有周期を持つ高層建築・免震構造物・長大橋梁においては、ハーフグラーベン起因の表面波(長周期地震動)増幅を考慮した応答スペクトル(サイト固有スペクトル)を設定しなければならない。

6.2 確率論的地震動予測(PSHA)に基づく極値リスク評価

茨城県南部では、プレート境界型(深度 40-50km)とスラブ内型(深度 60-70km)の異なる震源モデルが同一地点の直下に重なるため、ポアソンプロセスを仮定した発生確率の評価において以下の二項型複合累積確率密度関数モデルを適用する必要がある。P([span_11](start_span)N(t) \ge 1) = 1 - \exp\left[ - \left( \lambda_{\text{inter}} + \lambda_{\text{intra}} \right) t \right]

ここで \lambda_{\text{inter}} はプレート境界地震の平均発生率、\lambda_{\text{intra}} はスラブ内地震の平均発生率である。両項の和が極めて高いため、再現期間 30 年間における M \ge 6.0 の超過確率は首都圏内でも最高水準に位置づけられる。

6.3 リアルタイムモニタリングと事業継続計画(BCP)の高度化

高密度な地震活動域である茨城県南部においては、P波初動検知型(緊急地震速報)の猶予時間が極めて短い(S波到達まで数秒〜10秒程度)。したがって、施設単位でのオンサイトP波検知型自動制御システム(産業用ガス遮断、精密機械の退避機構、重要データの非同期ミラーリング等)の導入が最優先課題となる。

結論と今後の研究課題

茨城県南部で地震が多発する現象は、以下の要因が複合的に作用した結果である:

  1. オホーツク・フィリピン海・太平洋の3プレートによる複合的三重沈み込み構造。
  2. フィリピン海プレート北東端(ウェッジ部)に対する上下からの強力な応力集中。
  3. 含水鉱物の高圧脱水反応に伴う高圧間隙流体の滞留と、実効垂直応力の低下(脱水脆化)。
  4. 2011年東北地方太平洋沖地震による長期的な正のクーロン応力変化(\Delta CFF > 0)の蓄積。
  5. 関東平野の深部ハーフグラーベン構造に伴う強震動およびやや長周期地震動の特異的増幅。

今後の課題としては、超高精度3次元トモグラフィーによる流体挙動の時空間モニタリング、および連続微小地震のメカニズム変化に基づく大地震前兆的応力状態のセンシング手法の確立が求められる。

茨城県南部を震源とするプレート境界・スラブ内地震の発生メカニズムや、気象庁・ウェザーニュースによる揺れ方の解説については、関東で最大震度4 茨城県沖が震源 プレート境界で発生した地震 で専門的な図解とともに詳しく説明されています。

YouTube 動画の再生履歴は YouTube の履歴に保存されます。このデータは、YouTube の 利用規約 に従い、YouTube によって保存、使用されます

タイトルとURLをコピーしました