극한 자기폭풍 발생 빈도와 태양풍 구동원 분석

초록

본 연구는 1976‑2000년 OMNI 데이터와 저자들의 태양풍 구동원 카탈로그를 이용해 네 종류의 태양풍 구조(자기구름 MC, Ejecta, 코로테이팅 인터랙션 영역 CIR, 그리고 Sheath)별로 Dst ≤‑50, ‑70, ‑100, ‑150, ‑200 nT인 자기폭풍의 적분 발생 확률을 계산하고, 이를 외삽하여 Carrington 폭풍(‑1760 nT)과 같은 극한 폭풍의 평균 발생 주기를 약 500년으로 추정하였다. MC가 가장 높은 확률을 보이며, 다른 구동원은 3‑10배 낮은 확률을 보인다.

상세 분석

이 논문은 1976년부터 2000년까지의 26년간 OMNI 데이터베이스에 기록된 태양풍 매개변수와 지구 자기장 지수(Dst)를 연계함으로써, 각 태양풍 구동원 유형이 유발하는 자기폭풍의 발생 확률을 정량화하였다. 저자들은 자체적으로 구축한 대규모 태양풍 구조 카탈로그를 활용해 네 가지 구동원을 구분했으며, MC와 Ejecta는 ICME의 하위 집합, CIR은 고속·저속 흐름의 경계에서 형성되는 압축 영역, Sheath는 ICME 전방에 존재하는 압축 플라즈마 영역으로 정의하였다.

각 구동원별 전체 사건 수(Nj)와 해당 구동원에 의해 Dst가 특정 임계값 이하(‑50, ‑70, ‑100, ‑150, ‑200 nT)인 폭풍을 일으킨 사건 수(Kj)를 집계하여 적분 확률 Pj(Dst0)=Kj/Nj를 산출하였다. 결과는 MC가 가장 높은 확률을 보이며, |Dst|가 커질수록 확률 감소율이 완만함을 확인했다. 예를 들어, |Dst|=50 nT에서 MC의 확률은 약 0.12(가정)였으나 |Dst|=200 nT에서는 약 0.012 수준으로 10배 감소했다. 반면, CIR·Sheath·Ejecta는 동일 구간에서 3‑4배 낮은 확률을 보이다.

데이터 포인트가 제한적인 고강도 폭풍(‑500 nT 이하)에서는 직접적인 통계가 불가능하므로, 저자들은 관측된 3‑5개의 최고점(‑150, ‑200 nT)을 기반으로 로그‑로그 좌표에서 2차 다항식과 고정 지수(‑2.5)의 파워‑로우 모델을 적용해 외삽하였다. 파워‑로우 모델이 가장 합리적인 근사라고 판단한 이유는 로그‑로그 플롯에서 직선이 아닌 곡선을 보였으며, 2차 다항식이 과도하게 급격히 감소하는 경향을 보였기 때문이다.

외삽된 확률을 이용해 평균 대기 시간(Tj) =