캘리포니아 지역 ETAS 모델의 공간 가변 파라미터 객관적 추정

초록

ETAS 모델의 파라미터를 공간적으로 변하게 추정하기 위해 EM 알고리즘과 Voronoi tessellation을 결합한 효율적인 방법을 제시한다. 합성 카탈로그와 1981‑2015년 캘리포니아 지진 데이터를 적용한 결과, 분기 비율이 지열 흐름과 양의 상관관계를 보이며, 작은 지진이 트리거링에 주도적 역할을 함을 확인했다.

상세 분석

본 연구는 전통적으로 전역적인 파라미터를 가정해 온 ETAS(Epidemic Type Aftershock Sequence) 모델에 공간적 이질성을 도입함으로써 지진 발생 메커니즘을 보다 정밀하게 파악하고자 한다. 핵심 방법론은 두 단계로 구성된다. 첫 번째는 EM(Expectation‑Maximization) 알고리즘을 이용해 각 지진 사건이 배경 발생(background)인지, 혹은 다른 사건에 의해 유발된(aftershock)인지를 확률적으로 할당하는 과정이다. 이때 기대 단계에서는 현재 추정된 파라미터를 바탕으로 사건 간 트리거링 확률을 계산하고, 최대화 단계에서는 이러한 확률을 가중치로 사용해 파라미터를 업데이트한다. 두 번째는 공간적 파라미터 변동을 모델링하기 위해 Voronoi tessellation을 활용한 다중 모델 앙상블을 구축하는 것이다. 연구자는 다양한 Voronoi 셀 수(즉, 공간 해상도)를 가진 모델들을 생성하고, 각 모델의 로그우도와 자유도에 기반한 BIC(Bayesian Information Criterion)를 통해 모델 복잡도와 적합도를 동시에 평가한다. BIC가 최소인 모델들을 선택한 뒤, 선택된 모델들의 파라미터를 위치별 가중 평균하여 최종 앙상블 파라미터 지도를 만든다. 이러한 절차는 과적합을 방지하면서도 지역별 특성을 충분히 반영할 수 있는 장점을 제공한다. 합성 카탈로그 실험에서는 사전에 설정한 공간적 변동을 정확히 복원함으로써 방법론의 검증 가능성을 입증하였다. 실제 캘리포니아 데이터에 적용한 결과, 특히 분기 비율(n)과 배경 발생률(μ)이 지역에 따라 현저히 다름을 확인했다. n은 지열 흐름이 높은 지역에서 크게 나타났으며, 이는 열유동에 의해 유체 압력이 상승하고, 그에 따라 작은 지진이 더 쉽게 트리거링될 가능성을 시사한다. 반면 μ는 지열 흐름과 무관하게 비교적 균일하게 분포했으며, 이는 장거리 응력 전달에 의한 배경 지진 발생이 지열에 크게 좌우되지 않음을 의미한다. 또한, n과 μ는 깊이와도 상관관계가 없었는데, 이는 연구 영역 내에서 지진 발생 메커니즘이 깊이에 따라 변하지 않음을 보여준다. 마지막으로, 파라미터 추정 결과 작은 규모 지진이 전체 트리거링에 기여하는 비중이 크다는 점을 강조한다. 이는 기존의 정적 응력 변화 연구가 중대 규모 지진에만 초점을 맞추는 한계를 지적하고, 소규모 지진이 누적적인 응력 변화를 일으켜 장기적인 지진활동에 영향을 미칠 수 있음을 암시한다. 전반적으로 본 논문은 공간적 파라미터 변동을 정량화하는 새로운 프레임워크를 제공함으로써, 지진 위험 평가와 물리적 메커니즘 이해에 중요한 도구가 될 것으로 기대된다.