On the Influence of Fault Geometry on Earthquake Ruptures and Predictability

기술 요약: 단층 기하학이 지진 파쇄에 미치는 영향과 예측 가능성에 대하여

문제 제념
대규모 자연 지진을 예측하는 것은 지질학계의 중대한 과제로 남아 있으며, 이는 주로 단층 네트워크 내의 응력 재분배로 인해 발생하는 자기 유지적 복잡성과 혼돈적 거동 때문입니다. 역사적 기록에 기반한 통계적 접근 방식이 존재하지만, 이는 불완전한 데이터와 과거 사건의 파쇄 매개변수를 추정하기 어렵다는 한계가 있습니다. 또한, 순수하게 시간 의존적인 모델들은 미끄러짐(slip) 규모나 파쇄 시기를 정확하게 예측하는 데 실패하는 경우가 많습니다. 결정적이지만 아직 충분히 탐구되지 않은 요인은 단층의 굴곡(bends) 및 분지(branches)와 같은 단층 기하학이 파쇄 세분화(segmentation)와 그 결과로 나타나는 지진 시퀀스의 복잡성에 미치는 영향입니다. 저자들은 복잡한 기하학적 구조가 더 작고 세분화된 파쇄를 유도하며, 예측 가능성을 저하시키는 혼돈적 상호작용을 일으킨다고 상정하며, 따라서 현실적인 단층 구조를 고려한 물리 기반 모델이 필요하다고 주장합니다.

연구 방법론
본 연구는 복잡하고 비평면적인 단층 네트워크를 위해 설계된 준역학적 지진 시뮬레이터인 Quake-DFN을 사용하여 장기 수치 시뮬레이션을 수행합니다. 시뮬레이션은 2025년에 규모 7.7의 지진을 경험하였으며 지질학적 및 지형학적으로 잘 기록되어 있는 미얀마의 사가잉 단층(Sagaign Fault) 시스템에 맞춰 조정되었습니다.

• 시뮬레이션 설정: 저자들은 다음과 같이 다양한 단층 기하학을 적용하여 11회의 시뮬레이션을 수행했습니다:
  1. 평면적인 단층 기하 구조.
  2. 주 단층은 존재하나 분지가 없는 기하 구조.
  3. 굴곡과 "호스테일(horsetail)" 분지 구조를 포함한 전체 사가잉 단층의 복잡성을 반영한 9개의 시뮬레이션.
• 물리적 프레임워크: 단층 미끄러짐은 속도-상태 마찰(rate-and-state friction) 법칙에 의해 제어됩니다. 모델은 순수 탄성 매질을 가정하며, 장기적인 미끄러짐 속도를 연간 약 20mm로 유지하기 위해 지속적으로 미끄러지는 로딩 단층을 통해 원격 하중(far-field loading)을 가합니다. 시뮬레이션은 서로 다른 파쇄 거동을 탐색하기 위해 다양한 마찰 매개변수($a$, $b$, $D_c$)를 활용합니다.
• 예측 접근 방식: 저자들은 미끄러짐 결손(slip-deficit) 기반 예측 방법을 개발했습니다. 이 접근 방식은 대규모 지진이 미끄러짐 결손(장기적인 지구조적 미끄러짐과 역사적 사건으로부터의 표면 미끄러짐 사이의 차이)의 축적에 따라 재발한다고 가정합니다.
  • 예측 기준: 연속된 최소 10km 구간에서 미끄러짐 결손이 상한값(미끄러짐 결손 분포의 표준편차의 +1.73배로 정의됨)을 초과할 경우 지표 파쇄 사건이 발생할 것으로 예측합니다.
  • 파쇄 범위: 핵 생성 지점(nucleation point)이 식별되면, 예측된 파쇄는 미끄러짐 결손이 표준편차를 초과하는 연결된 구간을 따라 주향(strike) 방향으로 확장됩니다.
  • 평가: 성능은 실제 양성(true positives), 가짜 음성(false negatives), 가짜 양성(false positives)의 균형을 맞추는 F-score 지표와 주기성을 측정하는 사건 간 시간의 **변동 계수(COV)**를 사용하여 평가됩니다.

주요 결과

1. 기하학이 파쇄 복잡성에 미치는 영향:

   • 평면 단층: 부분적인 파쇄가 거의 없는, 본질적으로 주기적인 규모 7.7의 완전 파쇄 사건들을 생성했습니다.
   • 복잡한 기하 구조: 기하학적 복잡성(굴곡 및 분지)이 증가함에 따라 파쇄 시퀀스는 현저히 더 혼돈적이고 세분화되었습니다. 파쇄는 미끄러짐 속도의 분할로 인해 단층 굴곡과 분지 접합부에서 멈추는 경향을 보였습니다.
   • 현실성: 복잡한 기하 구조를 가진 특정 시뮬레이션은 150~200년의 재발 간격과 2025년 만달레이 규모 7.7 지진과 유사한 미끄러짐 분포를 포함하여, 사가잉 단층의 역사적 기록과 유사한 파쇄 시퀀스를 성공적으로 재현했습니다.

2. 예측 가능성 및 예측 성능:

   • 데이터 기간 의존성: 예측 정확도는 미끄러짐 결손을 계산하는 데 사용되는 역사적 데이터 기간의 길이에 크게 의존합니다. 복잡한 기하 구조의 경우, 데이터 기간이 대규모 사건의 재발 기간(150년)보다 짧으면 정확도가 낮습니다. 데이터 기간이 재발 기간(예: 200500년)을 초과하면 정확도가 크게 향상됩니다.
   • 기하학 vs 예측 가능성: 예측 가능성은 단층의 복잡성과 강한 상관관계가 있습니다.
     • 평면 단층 시뮬레이션에서는 충분한 데이터가 있을 때 F-score가 0.9를 초과했습니다.
     • 현실적이고 복잡한 기하 구조 시뮬레이션에서는 적절한 데이터가 있을 때 F-score가 0.6에서 0.8 사이였습니다.
     • 가장 낮은 예측 가능성(F-score < 0.5)은 세 개의 단층 상호작용이 발생하는 가장 복잡한 "호스테일" 구조에서 나타났습니다.
   • 규모 및 주기성: 단순한 단층 구간은 더 큰 지진(평균 Mw ~7.53)과 낮은 COV(더 주기적임)를 생성한 반면, 복잡한 구역은 더 작은 지진(평균 Mw ~7.12)과 높은 COV(더 군집화되거나 혼돈적임)를 생성했습니다. 미끄러짐 결손 방식은 준주기성에 의존하므로, 매우 복잡한 구역에서는 효과가 떨어집니다.

3. 예측의 한계:

   • 이 방법은 전통적인 의미에서의 엄격한 시간 예측 가능성(time-predictable) 또는 미끄러짐 예측 가능성(slip-predictable)을 갖지 않습니다 (Shimazaki 및 Nakata, 1980).
   • 예측의 가짜 양성은 "너무 이른" 예측으로 식별되었습니다. 즉, 미끄러짐 결손이 임계값을 초과했으나, 실제 파쇄는 이후 30년의 창(window) 내에서 발생했습니다.
   • 시뮬레이션은 계산 효율성을 위해 상대적으로 큰 핵 생성 길이($D_c$)를 사용했기 때문에 주로 큰 사건(M > 6.5)을 예측하며, 이는 현실에서 발생하는 더 작은 지진들을 놓칠 가능성이 있습니다.

의의 및 주장
본 논문은 단층 기하학이 지진 시퀀스의 복잡성을 유도하는 주요 동력이자 지진 예측 가능성을 제한하는 요소라고 주장합니다. 저자들은 미끄러짐 결손 기반 예측이 단순한 단층 구간의 대규모 준주기적 사건에는 효과적일 수 있지만, 상호작용이 혼돈적 파쇄 패턴을 만들어내는 복잡한 단층 구역에서는 그 효용성이 감소함을 입증했습니다.

본 연구는 효과적인 지진 예측이 단순히 통계적이거나 단순한 시간 의존적 모델을 넘어, 복잡한 단층 기하학과 그에 따른 하중 상호작용을 명시적으로 고려하는 포괄적인 물리 기반 접근 방식을 포함해야 한다고 제안합니다. 저자들은 수치 시뮬레이션이 관측된 단층 기하학 및 마찰 매개변수에 맞춰 조정될 때, 예측 가능성을 평가하고 향엔 앙상블 예측(ensemble forecasting)을 위한 유효한 도구로 기능할 수 있다고 제안합니다. 단, 이는 계산 속도가 넓은 매개변수 공간 탐색과 불확실성 포함을 허용할 때 가능하다는 전제가 붙습니다. 그러나 저자들은 실제 적용에 있어 불완전한 고지진 기록, 단층 특성의 불확실성, 그리고 단층 외 플라스틱 변형(off-fault plastic deformation) 및 마찰 불균질성의 배제와 같은 과제들이 존재함을 언급하며 신중한 태도를 유지했습니다.