지진 전조 현상으로서의 지자기 변동 개선 모델: 스코페와 인터매그넷 데이터 기반 연구

초록

본 논문은 스코페(북마케도니아)와 인터매그넷 네트워크의 지자기 관측 데이터를 이용해 일일 평균 1분 표준편차(σ) 변동이 지진 발생 전 조기 전조(Geomagnetic quake)로 작용함을 검증한다. 2001‑2004년 소피아 1‑성분 듀브나 플럭스게이트와 2004‑현재 스코페 3‑성분 덴마크 플럭스게이트 데이터를 분석했으며, 지구 조석(태양‑달 조석) 극값과 지진 발생 사이의 시간 차이가 평균 ±2.7일임을 확인하였다. 또한 INTERMAGNET 관측소 주변 강진 사례에서도 동일한 패턴이 재현되었다. 논문은 지진 발생 시점·위치·규모 예측을 위한 복합 프로그램과 FORTRAN 기반 분석·시각화 코드를 제시한다.

상세 분석

본 연구는 지자기 전조 현상(Geomagnetic quake, Gq)을 정량화하기 위해 일일 평균 1분 표준편차(σ) 값을 사용한다는 점에서 기존 연구와 차별화된다. 첫 단계에서는 2001‑2004년 소피아 관측소에서 1‑성분 듀브나 플럭스게이트를 이용해 σ의 급격한 상승과 지진 발생 사이의 상관관계를 탐색했으며, 이때 지구 조석의 극값(일양·일음)과 σ 급등이 동시 발생한다는 통계적 패턴을 확인했다. 두 번째 단계에서는 2004년 이후 스코페 지진관측소에 설치된 3‑성분 덴마크 플럭스게이트 데이터를 활용해 다중 성분의 σ 변동을 동시에 분석함으로써 전파 방향성 및 복합 변동성을 고려했다. INTERMAGNET 네트워크의 다중 관측소 데이터를 추가로 검증에 활용한 것은 연구의 외부 타당성을 크게 높였다.

통계적 분석에서는 σ 급등이 관측된 시점과 인접 지역(반경 200 km 이내)에서 발생한 강진(M ≥ 5.0) 사이의 시간 차이를 히스토그램으로 나타내어 평균 ±2.7일, 표준편차 1.1일이라는 구간을 도출했다. 이는 기존 지구 조석‑지진 연계 연구에서 제시된 1‑3일 범위와 일치하지만, σ 기반 전조가 보다 명확한 수치적 기준을 제공한다는 점에서 의미가 크다.

그러나 몇 가지 한계점도 존재한다. 첫째, σ 값은 외부 자기교란(태양풍, 전리층 변동)에도 민감하므로, 논문에서는 K‑인덱스와 같은 전리층 교란 지표를 보정했는지 명확히 기술되지 않았다. 둘째, 샘플 수가 제한적이다. 소피아와 스코페 관측 기간을 합쳐도 약 150건의 강진 사례만을 분석했으며, 이는 통계적 유의성을 확보하기에 부족할 수 있다. 셋째, FORTRAN 코드가 제공된다고는 하나, 오픈소스화 혹은 재현성을 위한 상세 매뉴얼이 부재해 다른 연구자가 동일한 파이프라인을 적용하기 어렵다.

이러한 점들을 보완한다면, σ 기반 지자기 전조는 실시간 지진 위험 평가 시스템에 통합될 가능성이 있다. 특히, 실시간 스트리밍 데이터와 머신러닝 기반 이상 탐지를 결합하면, 전조 탐지의 민감도와 특이도를 동시에 향상시킬 수 있다. 향후 연구에서는 다중 관측소 네트워크를 활용한 공간적 상관관계 분석, 전리층 교란 보정 모델 적용, 그리고 장기적인 데이터베이스 구축이 필요하다.