대형 지진은 불균형한 전조 사건 뒤에서 발생한다

초록

본 연구는 지진 발생 과정이 자기조직화 임계성(SOC)으로 간주될 때, 사건들의 에너지 분포 불균형을 나타내는 Gini·Kolkata 지수를 이용해 대형 지진의 전조를 탐색한다. 2차원 SOC 모형·열차 모형 시뮬레이션과 전 세계 여러 지역의 실제 지진 카탈로그를 분석한 결과, 큰 지진은 직전 100건의 사건이 높은 불균형(G≈0.87, k≈0.87)을 보인 뒤에 나타나는 경향이 확인되었다. 따라서 불균형 지수를 실시간 모니터링하면 임박한 대형 지진을 예측하고, 시스템이 임계점에 얼마나 근접했는지를 정량화할 수 있다.

상세 분석

이 논문은 지진 현상을 자기조직화 임계성(SOC) 시스템으로 모델링하고, 임계점 근처에서 나타나는 ‘불균형 현상’이 대형 사건의 전조가 될 수 있음을 실증한다. 먼저, 두 가지 전형적인 SOC 기반 모델을 사용한다. 첫 번째는 2차원 격자에서 응력이 일정 속도로 증가하고, 임계값을 초과하면 클러스터가 형성되어 ‘지진’이 발생하는 sandpile‑like 모델이다. 여기서 사건 규모는 클러스터 크기 n에 비례하고, 규모‑지진 진도 변환식 M=3/2 log n을 적용한다. 두 번째는 열차 모델로, 고정된 핀ning 포인트와 탄성 스프링으로 연결된 블록들이 끈질긴 마찰을 극복하며 미끄러지는 과정을 시뮬레이션한다. 두 모델 모두 사건 규모가 파워‑law(굿버그‑리히터) 분포를 보이며, 임계점에 접근할수록 큰 사건이 드물게 발생한다는 점에서 실제 지진과 일맥상통한다.

불균형을 정량화하기 위해 Gini 지수와 Kolkata 지수를 도입한다. 이는 원래 부의 분배 불균형을 측정하는 지표지만, 사건 에너지(또는 규모) 배열에 적용하면 ‘에너지 불균형’ 정도를 나타낸다. 논문에서는 각 사건에 대해 직전 100개의 사건을 슬라이딩 윈도우로 잡고, 그 윈도우 내 에너지 값을 정렬해 Lorenz 곡선을 만든 뒤 Gini(g)와 Kolkata(k)를 계산한다. 이후 계산된 g·k 값을 바로 다음 사건에 할당함으로써, 특정 사건이 발생하기 직전의 불균형 정도를 파악한다.

시뮬레이션 결과, 큰 사건(큰 규모의 avalanche)이 발생하기 전 윈도우의 g·k 값이 0.85~0.90 사이로 급격히 상승한다. 이는 임계점 근처에서 시스템의 응답이 ‘극단적 불균형’ 상태에 놓인다는 기존 이론과 일치한다. 실제 지진 데이터(USGS 카탈로그, 1975‑2025년, 규모 ≥ 4.5)에 대해서도 동일한 분석을 수행했으며, 남부 일본, 동남아시아, 북미 서부, 인도네시아 등 네 지역 모두에서 대형 지진(규모 ≥ 7.5) 직전의 g·k 평균이 약 0.87에 근접함을 확인했다. 특히, 고정 윈도우 방식뿐 아니라 ‘동적 리셋’ 방식(대형 사건이 나타나면 윈도우 시작점을 해당 사건으로 재설정)에서도 높은 불균형이 선행함을 보여, 윈도우 크기에 민감하지 않다는 점을 강조한다.

이러한 결과는 두 가지 중요한 시사점을 제공한다. 첫째, Gini·Kolkata 지수가 임계점 근접성을 정량화하는 보편적 지표로 활용될 수 있다. 둘째, 실시간으로 지진 연속 데이터를 모니터링하면서 불균형 지수를 계산하면, 대형 지진 발생 위험을 사전에 경고할 수 있는 새로운 위험 평가 도구가 된다. 논문은 또한 불균형 지수와 기존의 b‑값 감소, Omori‑Utsu 파라미터 변화 등 다른 전조 현상과의 연계 가능성을 제시하며, 향후 다중 지표 통합 모델링의 필요성을 강조한다.