태양 활동과 전 세계 대형 지진의 연관성 분석

초록

본 연구는 1996‑2015년 기간 동안 SOHO 위성에서 측정된 양성자 밀도·속도 자료와 ISC‑GEM 지진 카탈로그(M ≥ 5.8)를 비교하여, 양성자 밀도가 낮아진 뒤 1일 이내에 대형 지진 발생 확률이 유의하게 증가함을 밝혀냈다. 통계적 검증은 10⁵개의 부트스트랩 시뮬레이션을 통해 0.001 % 수준의 유의성을 확보했으며, 진도 임계값이 높을수록 상관 강도가 커지는 경향을 보였다. 저자들은 이를 전리층‑지구 전위 차이에 의한 역압전 효과로 설명한다.

상세 분석

이 논문은 전 세계 대형 지진이 포아송 분포를 따르지 않는다는 기존 연구를 출발점으로, 태양풍에 포함된 양성자 밀도 변화가 지진 발생에 미치는 영향을 정량적으로 검증하고자 한다. 데이터는 두 가지 주요 소스로 구성된다. 첫째, 1996년 1월부터 2015년 12월까지의 ISC‑GEM 카탈로그를 사용해 M ≥ 5.8 이상의 지진을 선별했으며, 이는 규모별 완전성을 검증한 뒤 사용하였다. 둘째, SOHO 위성의 L1 지점에서 수집된 양성자 밀도(ρ)와 속도(v) 데이터를 시간당 평균값으로 일일 평균을 산출하였다. 네 가지 파생 변수(밀도, 속도, 플럭스 ρv, 동압 ρv²/2)를 모두 고려했으며, 각 변수에 대해 최소값·최대값을 정규화한 비차원 평균 V_av_ad를 구했다.

통계적 검증 절차는 다음과 같다. 변수 V에 대해 임계값 V_T를 0.01 간격으로 변화시키며, “조건 C”(예: 변수 감소 후 1일) 만족 일수 D_C와 해당 일에 발생한 지진 수 E_C를 계산한다. 상대 발생률 R = (E_C/D_C)/(E/D) 를 구해, R이 1을 크게 벗어나는 경우를 탐색한다. 여기서 핵심은 부트스트랩이다. 실제 지진 간격 분포를 그대로 재조합해 10⁵개의 합성 카탈로그를 생성하고, 각 합성 데이터에 대해 동일한 R 값을 구해 실험값이 합성값을 초과하는 경우를 ‘유의함’으로 정의한다. 이 방법은 지진 발생이 비포아송적이라는 전제 하에, 실제 데이터와 동일한 생존함수를 유지하면서 무작위성을 검증한다.

결과는 양성자 밀도와 플럭스에 대해서만 ‘조건 1Dy bT’(임계값 이하로 떨어진 뒤 1일)에서 R이 유의하게 1보다 크게 나타났으며, 임계값이 높을수록 R이 증가하는 경향을 보였다. 구체적으로 밀도 ρ가 12.7–15.9 cm⁻² 사이로 떨어진 뒤 24시간 이내에 지진 발생률이 평균 대비 약 1.5배~2배 상승했다. 또한, 지진 규모 임계값을 M ≥ 6.0, 6.5, 7.0 등으로 올릴수록 R의 피크가 더욱 뚜렷해져, 대형 지진일수록 상관성이 강함을 시사한다.

논문은 이러한 통계적 상관을 물리적 메커니즘으로 설명하기 위해 ‘역압전 효과’를 제시한다. 양성자 밀도가 높아지면 전리층과 지표면 사이에 전위 차가 형성되고, 이 전위 차가 대형 단층을 통한 전류 흐름을 유도한다. 전류가 암석 내 압전성 광물(주로 석영)을 통과하면 역압전 효과에 의해 순간적인 응력 펄스가 발생하고, 기존의 응력 누적을 초과해 단층 파열을 촉발한다는 가설이다. 저자들은 이와 연관된 ‘지진광(Earthquake Lights)’이나 전파 방출 현상을 부수적인 증거로 제시한다.

비판적으로 보면, 몇 가지 한계가 존재한다. 첫째, 지진 카탈로그를 비감쇠(비decluster) 상태로 그대로 사용했는데, 이는 주된 사건과 여진을 구분하지 않아 상관성 과대평가 위험이 있다. 둘째, 양성자 밀도 자체가 비백색 잡음(white noise) 스펙트럼을 보이며, 27일 주기의 태양 회전 신호와 겹친다는데, 이는 실제 ‘인과관계’보다 공통 주기성을 반영할 가능성을 배제하지 못한다. 셋째, 부트스트랩 시뮬레이션은 실제 지진 발생 메커니즘을 반영하지 못한 채 단순히 인터‑이벤트 간격을 재배열하는데, 이는 복합적인 응력 누적 과정과는 차이가 있다. 마지막으로, 제시된 역압전 모델은 정량적 계산이 전혀 없으며, 전위 차, 전류 밀도, 응력 변화를 실제 지구 물리적 파라미터와 연결시키지 않아 검증 가능성이 낮다.

종합하면, 논문은 통계적으로 눈에 띄는 상관을 제시했지만, 인과관계를 확립하기 위한 물리적·지구화학적 검증이 부족하다. 향후 연구에서는 고해상도 전리층 전위 측정, 단층 전기전도도 조사, 그리고 역압전 효과를 수치 모델링하는 작업이 필요할 것이다.