일본 INTERMAGNET 관측망을 활용한 임박 지진 예측 모델 연구

초록

본 논문은 일본 INTERMAGNET 관측소(멤반베쓰, 카키오카, 카노야)의 지자기 데이터와 NEIC·NOAA 지진·조석 정보를 결합해 새로운 변수 Schtm(표면 에너지 밀도)을 도입, 과잉결정 방정식으로 지진 규모·심도·진원 거리의 즉시 예측 가능성을 탐색한다.

상세 분석

이 연구는 지진 전조 현상 중 가장 신뢰도가 높다고 주장되는 지자기 변동을 단일 선행 변수로 선택하고, 이를 기반으로 ‘지자기 지진(geomagnetic quake)’ 개념을 정량화한다. 핵심 변수인 Schtm은 지자기 변동의 에너지 밀도를 면적으로 환산한 값으로 정의되며, 관측점마다 시간‑주파수 영역에서 추출된 파워 스펙트럼을 적분해 산출한다. 논문은 세 개의 INTERMAGNET 관측소와 NEIC가 제공하는 지진 카탈로그를 이용해, 특정 시점에 기록된 Schtm 값이 가장 큰 사건을 ‘예측 지진’으로 가정한다.

이 가정을 수학적으로 전개하기 위해 저자는 과잉결정 선형 시스템을 구성한다. 변수는 지진 규모(M), 깊이(D), 관측점으로부터의 거리(R)이며, 각 관측소에서 측정된 Schtm을 함수 f(M, D, R) 형태로 표현한다. 네 개 이상의 관측점이 확보될 경우, 3개의 미지수를 풀기 위한 충분한 방정식이 생성되어 최소 자승법 등으로 해를 구한다는 것이 논문의 핵심 논리다.

방법론적 강점은 다음과 같다. 첫째, INTERMAGNET 네트워크는 고품질, 연속적인 1‑초·1‑분 지자기 데이터를 제공하므로 데이터 신뢰성이 높다. 둘째, Schtm이라는 새로운 정량적 지표를 도입해 기존의 ‘지자기 변동’ 서술을 수치화함으로써 통계적 검증이 가능하도록 만든 점은 의미 있다. 셋째, 과잉결정 시스템을 활용해 다중 관측점의 정보를 동시에 이용하려는 시도는 다변량 예측 모델 설계에 있어 혁신적이다.

하지만 몇 가지 중대한 한계와 의문점도 존재한다. 첫째, Schtm이 실제 물리적 메커니즘(예: 전류 흐름, 스트레스‑전기 변환)과 어떻게 연결되는지에 대한 이론적 근거가 부족하다. 변수 정의가 ‘표면 에너지 밀도’라 함은 직관적이지만, 지자기장 변동이 전도성 매질을 통해 어떻게 에너지 밀도로 전환되는지에 대한 물리적 모델이 제시되지 않는다. 둘째, 논문은 단일 사건(또는 소수의 사건)만을 사례 연구로 제시하고, 통계적 유의성 검증(예: ROC 곡선, 신뢰구간, p‑값)이나 교차 검증이 전혀 이루어지지 않았다. 이는 모델의 일반화 가능성을 판단할 근거가 부족함을 의미한다. 셋째, 과잉결정 방정식의 구체적 형태와 해법(선형/비선형, 조건수, 정규화 방법 등)이 상세히 기술되지 않아 재현 가능성이 낮다. 네 개 이상의 관측점이 필요하다고 주장하지만, 실제로는 일본 내 3개 관측소만을 사용했으며, 700 km 반경 내에 충분히 밀집된 관측소가 존재하지 않을 경우 적용이 제한될 수 있다. 넷째, 조석·조류 데이터와 NOAA의 기타 변수(예: 전리층 전류, 대기 전기장)와의 통합이 언급되지만, 본 연구에서는 전혀 활용되지 않아 ‘다중 전조 변수 네트워크’라는 목표와의 괴리가 있다.

결론적으로, 이 논문은 지자기 데이터를 이용한 임박 지진 예측에 대한 흥미로운 아이디어와 새로운 정량 지표를 제시했지만, 물리적 메커니즘 설명, 통계적 검증, 수학적 구현 상세화가 부족하다. 향후 연구에서는 Schtm과 지진 파라미터 간의 인과관계를 물리‑수학적으로 모델링하고, 대규모 데이터셋을 통한 검증, 그리고 다중 전자기·지구물리학적 전조 변수의 융합을 통해 실용적인 예측 시스템으로 발전시킬 필요가 있다.