용암 돔 분출 지속시간 예측을 위한 베이지안 파레토 모델

이 연구는 용암 돔형 화산 분출의 지속시간을 정량적으로 예측하기 위한 새로운 통계적 방법론을 제시한다. 저자들은 기존에 구축된 DomeHaz 데이터베이스에서 177개의 용암 돔 분출 사건을 선정했으며, 여기에는 완전 종료된 사건 163건과 2015년 기준 진행 중인 사건 14건이 포함된다. ‘분출’은 용암 돔 성장 단계와 그 사이의 비활동(정지) 기간을 모두 포함하는 연속적인 활동으로 정의한다. 이러한 정의는 기존 연구와 달리 장기적인 비활동 기간도 화산의 불안정성을 반영한다는 점에서 의미가 있다. 데이터 분석 결과, 지속시간에 대한 누적 빈도는 로그-로그 플롯에서 거의 직선 형태를 보였으며, 이는 파레토 분포가 데이터에 적합함을 시사한다. χ² 적합 검정에서도 파레토 모델은 p=0.1557로 귀무가설을 기각하지 못했으며, 반면 지수분포는 p=8.25×10⁻²⁷로 완전히 부적합했다. 따라서 저자들은 일반화 파레토 분포(GPD) P(T>t)=(1+t/β)^{-α}를 기본 모델로 채택한다. 여기서 α는 형태 파라미터, β는 규모 파라미터이며, α<1인 경우 평균 지속시간이 무한대가 되어 중위수와 사분위수와 같은 유한 통계량에 초점을 맞춘다. 모델 추정은 베이지안 프레임워크를 사용한다. 사전분포를 최소화한 객관적 베이지안 접근법을 적용해 파라미터 α와 β를 최대우도(MLE)와 헤시안 역행렬을 통해 추정하고, 표준오차를 계산한다. 관측 데이터는 완전 관측(시작·종료 시점 모두 알려진 경우)과 검열 데이터(진행 중인 경우)로 구분되며, 두 종류 모두를 포함하는 전체 가능도 함수를 구축한다. 세 가지 모델 변형을 비교한다. 첫 번째는 모든 사건을 하나의 GPD에 적합한 Aggregate 모델이며, α̂=0.6487±0.0132, β̂=0.7018±0.055 yr를 얻었다. 두 번째는 화산암 성분에 따라 현무암·안산암·유문암/류석 세 그룹으로 나누어 각각 별도 GPD를 적합한 Grouped 모델이다. 현무암군은 α≈0.67, β≈0.37; 안산암군은 α≈1.09, β≈1.21; 유문암/류석군은 지수분포 λ≈0.34가 최적이며, 이는 α가 1에 가까워 파레토 꼬리가 약해짐을 의미한다. 세 번째는 실리카 함량 X(%)에 대한 로그선형 회귀 모델로, α=α₀·e^{γα(X−60)}, β=β₀·e^{γβ(X−60)} 형태를 취한다. 회귀 결과 γα=0.0447, γβ=0.1302가 통계적으로 유의미해, 실리카 함량이 1% 증가할 때 α는 약 4.5%, β는 약 13% 증가한다는 결론을 얻는다. 예측 단계에서는 현재까지 s년 동안 지속된 분출에 대해 남은 기간 T의 조건부 분포가 GPD(α, β+s)임을 이용한다. 이를 통해 진행 중인 14개 화산에 대해 중위수와 90% 신뢰구간을 산출했으며, 특히 산프리레 힐스(중위수 35년, 10–139년)와 시나붕(중위수 4년, 1–17년)의 불확실성이 크게 나타났다. 베이지안 예측은 파라미터 불확실성(인식적)과 지질 과정의 자연 변동성(알레아터리)을 모두 반영한다는 점에서 기존의 평균 기반 예측보다 신뢰성이 높다. 논문의 주요 기여는 다음과 같다. 첫째, 용암 돔 분출 지속시간이 중량 꼬리 분포를 따른다는 실증적 증거를 제공한다. 둘째, 화산암 성분에 따라 파라미터가 달라짐을 정량화하여, 성분별 위험 평가가 가능하도록 한다. 셋째, 객관적 베이지안 접근법을 통해 불확실성을 정량화하고, 진행 중인 화산에 대한 실용적인 남은 기간 예측을 제공한다. 넷째, 이 방법론은 다른 유형의 화산 분출, 지진, 산사태 등 지질재해의 지속시간 예측에도 확장 가능함을 제시한다. 따라서 장기 재난 대비와 위험 관리, 인프라 투자, 주민 대피 계획 등에 직접 활용될 수 있는 과학적 근거를 제공한다.