북서환해대권 지진 발생의 비정통 통계와 복합성 동역학

초록

본 연구는 2002‑2016년 일본 기상청 지진 카탈로그를 이용해, 북서환해대권(류큐·이즈보닌·혼슈 아크)에서 얕은(지각)과 깊은(아래 지각) 배경 지진의 통계적 특성과 시간‑공간 상관성을 비확장 통계 물리(NESP)를 통해 분석한다. 얕은 배경 지진은 약‑중간 정도의 장기 상관성을 보이며, 대규모 사건 전후에 상관도가 변동한다. 반면 깊은 배경 지진은 거의 포아송적(무상관) 특성을 나타낸다. 이러한 차이는 변환형 변곡대(캘리포니아·알래스카)와는 반대이며, 자유 표면과 고정 경계 조건이 복합성 발달에 미치는 영향을 시사한다.

상세 분석

본 논문은 비확장 통계 물리학(Non‑Extensive Statistical Physics, NESP)의 핵심 파라미터인 Tsallis q‑값을 이용해 지진 발생 간 간격(시간·거리)과 규모‑시간의 다변량 확률분포(F‑M‑T, F‑D)를 정량화한다. 먼저 일본 기상청이 제공한 2002‑2016.5년 구간의 완전하고 균일한 카탈로그를 사용했으며, Zh‑Uang 등(2002)의 확률적 declustering 기법을 적용해 전경(후진·여진) 사건을 제거하고 순수한 배경 지진만을 추출하였다. 이렇게 얻어진 배경 시계열에 대해, (1) 간격 시간 분포(F‑T)의 꼬리 지수와 q‑값을 추정해 장기 상관성의 강도를 평가하고, (2) 사건 간 거리 분포(F‑D)를 통해 공간적 상호작용 범위를 파악하였다. 또한, Gutenberg‑Richter b‑값을 함께 분석해 규모‑공간 구조와의 연관성을 검증하였다.

얕은(지각) 배경 지진에 대해서는 q‑값이 1.1~1.3 사이로, 완전한 포아송( q=1 )보다 약간 높은 값을 보였으며, 이는 약‑중간 정도의 장기 상관성을 의미한다. 특히, 2011년 도호쿠 대지진 전후에 q‑값이 일시적으로 상승(≈1.35)하고, 이후 다시 감소하는 동적 변화를 보였는데, 이는 대규모 응력 누적과 방출 과정이 배경 시계열에 일시적 ‘임계화’를 유발한다는 해석을 가능하게 한다. 시간‑거리 상관 함수는 지표면 근처에서 10‑30 km 정도의 유효 상호작용 반경을 나타냈으며, 이는 단층 네트워크가 작은 세계(small‑world) 구조를 형성한다는 기존 시뮬레이션 결과와 일치한다.

반면, 깊은(아래 지각) 배경 지진은 q‑값이 1.00‑1.05 수준으로 거의 포아송적 특성을 보였다. 특히, Wadati‑Benioff 구역 전역에서 관측된 q‑값은 통계적 유의미성을 갖지 못했으며, 거리 상관 함수도 급격히 감소해 5 km 이하에서만 의미 있는 상호작용을 나타냈다. 이는 깊은 영역이 고정 경계 조건(압축성 높은 매질, 제한된 변형률) 하에 놓여, 자유 표면에서 발생하는 복합성 메커니즘이 억제된다는 가설을 뒷받침한다.

지진 발생 모델 측면에서, 전통적인 포아송 기반 ETAS(후진 시퀀스) 모델은 얕은 배경 지진의 약‑중간 상관성을 충분히 설명하지 못한다는 점을 지적한다. 반면, 자기조직임계성(SOC) 혹은 비평형 복합 네트워크 모델은 관측된 q‑값과 시간‑공간 상관 구조를 자연스럽게 재현한다. 특히, 자유 표면을 갖는 변환형 변곡대(캘리포니아·알래스카)와 달리, 북서환해대권의 수렴형 변곡대는 ‘고정 경계’가 복합성 발달을 억제한다는 결론을 도출한다.

결론적으로, 본 연구는 (1) 얕은 배경 지진은 비포아송적, 약‑중간 상관성을 가지며 대규모 사건 전후에 동적으로 변한다, (2) 깊은 배경 지진은 거의 포아송적이며 시간‑공간 상관성이 약하다, (3) 지구물리적 경계 조건이 복합성 발달에 결정적인 역할을 한다는 점을 실증하였다. 이러한 결과는 지진 위험 평가와 장기 시뮬레이션 모델링에 비확장 통계 물리학을 적용할 필요성을 강조한다.