BMR 출현 역추적 두 그룹과 경사각 분석에 미치는 영향

초록

본 연구는 자동 BMR 추적 알고리즘 AutoTAB에 출현 역추적 모듈을 추가하여 12 000개 중 3 080개의 양극성 자기 영역을 초기 출현 시점까지 되돌렸다. 역추적 결과, 플럭스가 크게 증가하는 진정한 출현군과 플럭스 변화가 거의 없는 ‘버려진’ 군으로 나뉘었으며, 후자는 경사각 분포가 무작위이고 Joy’s law와 일치하지 않음을 확인했다. 따라서 BMR 통계 분석 시 버려진 군을 제외해야 물리적 특성을 정확히 파악할 수 있다.

상세 분석

본 논문은 태양 광구에 나타나는 양극성 자기 영역(BMR)의 초기 출현 과정을 정량적으로 파악하기 위해, 기존 자동 검출·추적 시스템인 AutoTAB에 ‘역추적(backtracking)’ 모듈을 구현한 점이 가장 큰 공헌이다. 데이터는 SOHO/MDI와 SDO/HMI의 LOS(선‑오프) 자장 영상을 1996‑2024 기간 동안 사용했으며, HMI 자료는 MDI 해상도에 맞게 재샘플링하고 1.4배 스케일 보정하였다.

검출 단계에서는 각 픽셀을 태양 중심각 코사인으로 나누어 투영 효과를 보정하고, 11 픽셀(≈22″) 박스카 평균을 이용한 적응 임계값을 적용해 강자장을 추출한다. 이후 양극성 균형 비(r = |Φ⁺‑Φ⁻|/(Φ⁺+Φ⁻))가 0 이하인 경우에만 BMR으로 인정함으로써, 초기 단계에서 불균형한 잡음이나 단일극성 구조가 포함되는 것을 방지한다.

추적 단계는 차동 회전 프로파일을 이용해 동쪽에서 서쪽으로 이동하는 동안 예상 위치를 예측하고, 면적 변화 비율(80‑300 %)를 기준으로 전후 프레임을 연결한다. 이 과정에서 데이터 결측이 있으면 연관 기준을 강화해 추적 오류를 최소화한다.

역추적 알고리즘은 첫 검출 시점(T₀)부터 가능한 가장 이른 시점(Tₑ)까지 역방향으로 프레임을 스캔한다. 각 프레임에서 9 픽셀(≈18″) 확장된 사각형 영역을 정의하고, 강자장 픽셀 비율 R(t)와 총 무부호 플럭스 Φ(t)를 계산한다. 두 비율 rₚ = R(t)/R(T₀)와 r_f = Φ(t)/Φ(T₀)를 이용해 성장·소멸 여부를 판단한다. r_f > 1·rₚ > 1이면 비현실적인 성장으로 해당 프레임을 제외하고, r_f ≤ 0.4·rₚ ≤ 0.5이면 급격히 소멸하거나 잡음 수준 이하로 간주해 연속 5프레임 이상이 누락될 경우 역추적을 종료한다. 이러한 임계값은 실험적으로 최적화되었으며, 전체 11 987개의 BMR 중 3 080개(≈25 %)에 대해 성공적인 역추적을 달성했다.

역추적 결과는 두 개의 뚜렷한 군으로 구분된다. ‘출현군’은 Φ와 R이 시간에 따라 지속적으로 증가하며, 초기 단계에서 약 100 G 이상의 강자장 픽셀 비율이 급격히 상승한다. 반면 ‘버려진 군’은 초기 검출 시점과 역추적된 출현 시점 사이에 플럭스 증가가 거의 없으며, 종종 초기 검출 자체가 이미 성숙 단계이거나 잡음에 의해 오인된 경우이다. 통계적으로 버려진 군은 경사각(tilt angle) 분포가 균등하고, 위도와의 상관관계(Joy’s law)가 사라진다. 이는 전통적인 BMR 카탈로그에 포함된 많은 영역이 실제 출현이 아니라 이미 성숙하거나 소멸 단계에 있는 경우가 많음을 시사한다.

이러한 결과는 BMR의 경사각, 플럭스, 성장 속도 등을 이용한 태양 다이너모 이론 검증에 직접적인 영향을 미친다. 특히 Joy’s law를 기반으로 한 플럭스 운반 모델이나 일주기 주기 예측에 버려진 군을 포함하면 통계적 편향이 발생한다. 따라서 향후 BMR 데이터베이스 구축 시 초기 출현 검증 절차를 필수적으로 포함해야 하며, 자동화된 역추적 모듈은 대규모 데이터셋에 적용 가능한 실용적인 도구로 평가된다.

본 연구는 자동 검출·추적 시스템에 초기 출현 검증을 결합함으로써, 기존 카탈로그가 가지고 있던 ‘출현 편향’ 문제를 정량적으로 해결하고, BMR 물리학 연구에 필요한 고품질 샘플을 제공한다는 점에서 의미가 크다. 향후 연구에서는 더 높은 시간 해상도(예: 45 s HMI)와 전면·후면 관측(예: STEREO) 데이터를 결합해 출현 시점을 더욱 정확히 측정하고, 머신러닝 기반의 출현 패턴 인식을 도입함으로써 현재의 임계값 기반 방법을 보완할 여지가 있다.