일본 지진 최대 진도 분포와 미래 예측

초록

본 연구는 일본(1923‑2007) 지진 카탈로그에 극값 이론을 적용해 미래 일정 기간 동안 발생할 최대 진도(M) 분포를 추정한다. 일반화 극값 분포(GEV)와 일반화 파레토 분포(GPD)의 이중성을 이용해 양자 Qq(τ)를 계산하고, 전통적인 절대 최대 진도 Mmax보다 안정적인 위험 지표임을 입증한다.

상세 분석

본 논문은 극값 통계학의 두 핵심 정리, 즉 최대값 정리와 초과값 정리를 기반으로 하여 지진 규모의 꼬리 분포를 정밀하게 추정한다. 기존 연구인 Pisarenko et al. (2009)의 방법을 개선해 GPD와 GEV 사이의 수학적 이중성을 명시적으로 활용함으로써, 관측된 지진 순간순간의 규모가 아닌 일정 기간 τ 내에 발생할 최대 규모의 확률분포를 직접 도출한다. 일본 전국의 1923년부터 2007년까지의 카탈로그를 사용해 GEV의 형태 매개변수 ξ, 위치 매개변수 μ(τ), 척도 매개변수 σ(τ)를 추정했으며, ξ = –0.1901 ± 0.0717이라는 음수값은 분포가 유한한 상한을 갖는 Weibull형임을 시사한다. 이는 지진 규모가 무한히 커질 가능성이 낮으며, 실제 물리적 한계가 존재한다는 지질학적 가설과 일치한다. 위치 매개변수 μ는 τ=200년 기준 6.3387 ± 0.0380으로, 장기 평균 최대 진도가 약 6.3 수준임을 나타낸다. 척도 매개변수 σ는 0.5995 ± 0.0223으로, 분포의 폭이 비교적 좁아 예측 불확실성이 제한적임을 보여준다.

특히 저자들은 Qq(τ)라는 양자를 도입해, 원하는 신뢰수준 q(예: 0.90)와 미래 기간 τ(1~50년)별 최대 진도 기대값을 제공한다. 예를 들어 τ=10년, q=0.90일 때 Q0.90 = 8.34 ± 0.32로, 다음 10년 동안 90% 확률로 최대 진도가 8.34를 초과하지 않을 것이라는 실용적 정보를 제공한다. 반면 전통적인 Mmax 추정값은 9.57 ± 0.86으로, 동일 신뢰수준에서 변동성이 두 배 이상 크다. 이는 Qq(τ)가 Mmax보다 통계적 안정성이 높으며, 재난 대비 정책에 더 적합함을 의미한다.

방법론적 측면에서 저자들은 최대우도 추정(MLE)과 부트스트랩 재표본화를 결합해 매개변수의 불확실성을 정량화하였다. 또한, 관측 데이터의 독립성 가정을 검증하기 위해 시간 간격을 조정하고, 규모-시간 상관성을 최소화하는 절차를 수행했다. 이러한 절차는 극값 모델링에서 흔히 발생하는 과적합 위험을 감소시킨다.

결과적으로, 본 연구는 지진 위험 평가에 있어 “절대 최대 규모” 대신 “특정 기간·신뢰수준 양자”를 활용하는 새로운 프레임워크를 제시한다. 이는 보험·재보험, 인프라 설계, 국가 재난 대비 전략 등 다양한 분야에서 적용 가능하며, 특히 장기적인 위험 관리와 비용‑편익 분석에 유용하다.