STORMY : A Real-time Triggering Framework using Yamagawa Solar Spectrograph for Active Solar Emission Observations with the MWA

기술 요약: STORMY – MWA 태양 관측을 위한 실시간 트리거링 프레임워크

문제 정의
태양 라디오 버스트(Solar radio bursts)는 태양 물리학 및 우주 기상에 대한 중요한 통찰력을 제공하지만, 본질적으로 예측 불가능하고 일시적이며 빈도가 낮습니다. SKA의 전신인 Murchison Widefield Array(MWA)와 같은 현대적인 라디오 간섭계는 우수한 이미징 능력을 제공하지만, 관측 시간 할당이 매우 제한적입니다. 블라인드 방식의 정기적인 태양 관측에 의존하는 것은 비효율적이며, 중요한 태양 이벤트를 놓칠 가능성이 높습니다. 아카이브 데이터를 위한 자동 탐지 알고리즘은 존재하지만, 트리거를 생성하기 전에 전체 이벤트 형태가 모두 나타나야 하므로 실시간 트리거링에는 적합하지 않은 경우가 많습니다. 따라서 태양 활동이 감지되는 즉시 이를 포착하여 한정된 관측 시간을 효과적으로 활용할 수 있는 견고한 저지연(low-latency) 프레임워크가 필요합니다.

방법론
저자들은 MWA를 위해 설계된 실시간 트리거링 프레임워크인 STORMY(Solar Triggered Observations of Radio bursts using MWA and Yamagawa)를 개발했습니다. 이 시스템은 다음과 같은 워크플로우로 작동합니다:

1. 데이터 소스 및 지연 관리:

   • 트리거는 70 MHz–9 GHz 대역에서 작동하는 Yamagawa 태양 분광계(일본)의 근실시간 데이터로부터 유도됩니다. 경도가 유사하기 때문에 Yamagawa의 관측 창은 MWA(80–300 MHz)와 완전히 중첩됩니다.
   • 시스템은 2분의 지연 시간을 가진 Yamagawa의 커스텀 데이터 스트림을 활용합니다.
   • 이 지연으로 인한 데이터 손실을 완화하기 위해, MWA는 전압 링 버퍼(최대 160초의 원시 전압 데이터를 저장)를 유지합니다. 트리거가 수신되면 버퍼가 동결되어 디스크로 덤프되며, 이를 통해 트리거 도착 전 발생한 태양 이벤트의 초기 단계를 복구할 수 있습니다.

2. 전처리 파이프라인:

   • 대역 통과 보정(Bandpass Correction): 계통적 스펙트럼 응답을 제거하기 위해 이전 전체 일일 데이터(180초 단위로 분할됨)에서 계산된 중앙값 대역 통과(median bandpass)를 사용하여 원시 동적 스펙트럼을 정규화합니다.
   • RFI 완화: 시스템은 무선 주파수 간섭(RFI)을 제거하기 위해 다음의 특정 알고리즘을 사용합니다:
     • 주기적 RFI: 300초 주기를 갖는 피크 찾기(peak finding)를 통해 식별 및 제거됩니다.
     • 협대역 지속적 RFI: 실제 태양 신호 플래깅을 피하기 위해 MAD(Median Absolute Deviation) 기반의 수정된 z-score를 사용하여 필터링됩니다.
     • 일시적 RFI: 형태학적 연산(morphological operations)을 통해 처리됩니다.
   • 배경 차감(Background Subtraction): 느리게 변하는 배경 노이즈(계통 이득, 시스템 온도, 조용한 태양)를 차감하기 위해 이동 중앙값 필터(윈도우 폭 1800초)가 적용됩니다.

3. 이벤트 탐지:

   • 신호 대 잡음비(SNR) 동적 스펙트럼이 생성됩니다.
   • 이진 동적 스펙트럼은 정밀도와 재현율(precision and recall) 사이의 균형을 최적화하기 위해 결정된 13σ 임계값을 사용하여 생성됩니다.
   • 형태학적 폐쇄(Morphological Closing): 2단계 형태학적 폐쇄 연산(2×2 커널을 사용한 팽창 후 침식)은 RFI 플래깅으로 인해 누락된 픽셀을 채웁니다.
   • 윤곽선 탐지: 닫힌 윤곽선은 OpenCV를 사용하여 식별됩니다. 이벤트는 최소 15 MHz의 주파수 범위와 3초의 시간 범위를 가질 경우 유효한 것으로 간주됩니다.

4. 트리거링 및 관측:

   • 탐지 시, 웹 기반 메커니즘이 MWA를 트리거합니다.
   • 시스템은 버스트의 스펙트럼 범위를 추정하고, 스펙트럼 드리프트 및 약한 방출을 고려하여 가드 밴드(아래로 30 MHz, 위로 10 MHz)를 추가하여 MWA의 30.72 MHz 대역폭을 구성합니다.
   • MWA는 트리거된 관측을 위해 "피켓 펜스(picket fence)" 모드(로그 간격 주파수 청크)로 전환됩니다.

주요 기여

• STORMY 개발: MWA와 Yamagawa 데이터를 사용하여 태양 라디오 버스트를 위한 최초의 근실시간 트리거링 프레임워크를 구현했습니다.
• 알고리즘 혁신: 노이즈가 많은 실시간 데이터 스트림에서 버스트를 탐지하기 위해 중앙값 기반 디노이징, 수정된 z-score RFI 제거, 형태학적 폐쇄를 결리한 견고한 파이프라인을 결합했습니다.
• 지연 보상: 트리거 생성 지연(약 70~80초) 동안 손실된 이벤트 데이터를 복구하기 위해 MWA의 전압 링 버퍼를 사용하는 실용적인 데모를 제시했습니다.
• 운영 배포: STORMY는 정기적인 태양 관측을 위해 MWA에서 운영 가능한 상태로 구축되었습니다.

결과

• 성능 지표: 226개의 실제 이벤트를 포함하는 8,400개의 2분 데이터 청크를 대상으로 한 테스트 결과, 13σ 임계값이 정밀도와 재현율의 균형을 맞추는 최적의 F1-score를 제공함을 보여주었습니다.
• 운영 성공: 2024년 7월 31일부터 2025년 1월 31일 사이에 STORMY는 약 110회의 태양 버스트를 성공적으로 트리거하고 관측했습니다.
• 사례 연구: 논문은 두 가지 성공적인 캡처 사례를 제시합니다:
  1. 2024년 8월 1일 발생한 M-클래스 X선 플레어: 시스템이 126 MHz에서 관측을 트리거하였으며, SDO/AIA 데이터와 일치하는 라디오 이미지를 포착했습니다.
  2. 2024년 11월 4일 발생한 그룹 Type-III 이벤트: 버퍼 데이터가 성공적으로 상관관계 분석되어 LASCO C2 데이터 위에 중첩된 103 MHz에서의 라디오 이미지를 생성했습니다.

의의 및 주장
본 논문은 STORMY가 블라인드 스케줄링의 비효를 해결하여, 간섭계 데이터를 통해 일시적인 태양 활동을 포착할 수 있는 견고하고 효율적인 수단을 제공한다고 주장합니다. 저자들은 이 프레임워크가 다음을 가능하게 한다고 강조합니다:

• MWA, LOFAR, OVRO-LWA와 같은 기존 장비를 사용하여 이벤트가 풍부한 태양 데이터를 효율적으로 포착할 수 있게 합니다.
• SKAO(특히 SKA-Low)를 이용한 트리거 관측을 위한 길을 열어줍니다. 저자들은 SKA-Low의 "태양 서브어레이(solar subarray)" 개념이 고대역폭 관측을 트리거하기 위해 유사한 전압 버퍼링 및 실시간 탐지를 활용할 수 있다고 제안하며, 이를 통해 다른 과학 관측에 미치는 영향을 최소화하면서 예측 불가능한 태양 이벤트를 포착하는 능력을 극대화할 수 있다고 제안합니다.
• 태양 트리거링이 현재 기술로 가능하다는 것을 입증하며, 진행 중인 태양 극대기 동안 미래의 태양 라디오 천문학을 최적화할 수 있는 경로를 제시합니다.