Hour-scale response of low-Earth-orbit satellite altitude to extreme geomagnetic storms revealed by the NinjaSat CubeSat during the 2024 solar maximum

기술 요약: 극한 지자기 폭풍 발생 시 저궤도 위성 고도의 시간 단위(Hour-scale) 반응

문제 정의
저궤도(LEO) 위성 인구의 급격한 증가는 충돌 회피 및 궤도 지속 가능성을 보장하기 위해 정밀하고 지연 시간이 짧은 열권 밀도 예측을 필요로 한다. 초세속적(Secular) 궤도 감쇠는 잘 이해되어 있으나, 지자기 폭풍에 의해 유도되는 충격적인 밀도 증가는 시간 단위 규모에서 여전히 미해결 상태로 남아 있다. 기존의 운영 모델(예: NRLMSISE-00, DTM-2013, JB208)은 조악한 시간 해상도(3시간 단위 Kp/Ap 또는 1시간 단위 Dst)를 가진 지자기 지수를 활용하거나, 폭풍 시 발생하는 특유의 가열 메커니즘을 포착하지 못한다. 핵심적인 미해결 과제는 어떤 지자기 관측값이 급격한 궤도 감쇠를 유발하는 충격적 가열의 가장 직접적인 인과적 프록시(Proxy) 역할을 하는지, 그리고 에너지 입력으로부터 대기 팽창까지의 구체적인 수송 시간 규모가 얼마인지이다.

연구 방법론
본 연구는 2024년 태양 극대기 동안 운용된 두 개의 X선 천문 위성으로부터 얻은 고빈도 GNSS 텔레메트리를 활용하였다:

1. NinjaSat: 태양 동기 궤도($97^\circ$ 경사각, 약 500 km 고도)를 도는 6U 큐브위성으로, 10분 간격의 위치 데이터를 제공한다.
2. XRISM: 저경사 궤도($31^\circ$, 약 560 km 고도)를 도는 X선 이미징 및 분광 미션 위성으로, 30초 간격의 데이터를 제공한다.

폭풍에 의한 섭동을 배경 현상으로부터 분리하기 위해, 저자들은 GNSS 유도 고도($a_{GNSS}$)와 SGP4 모델을 사용하여 전파된 NORAD TLE 고도($a_{SGP4}$) 사이의 차이로 정의된 고도 잔차($\Delta a$)를 구축하였다. SGP4는 초세속적 감쇠는 포착하지만 충격적인 밀도 스파이크는 포착하지 못하므로, 양(+)의 $\Delta a$는 강화된 항력을 나타낸다.

분석 과정은 다음과 같다:

• 배경 제거: 지자기적으로 안정된 기간(2024년 7월~9월) 동안 GOES-18 EUV 복사 강도(140.5 nm)에 대한 평활화된 고도 변화율($\dot{g}\Delta a$)의 선형 회귀를 수행하여 EUV 유도 배경을 모델링하고 차감하였다.
• 교차 상관(Cross-Correlation): 잔여 폭풍 시기 고도 변화율을 극지 캡(Polar Cap, PC) 지수(고위도 이온층으로의 태양풍 에너지 입력을 나타내는 프록시) 및 Dst 지수(환 전류의 프록시)와 교차 상관 분석하였다.
• 위도별 비교: NinjaSat(극지 영역 커버리지)과 XRISM(저위도 영역 커버리지)의 반응을 비교함으로써, 두 차례의 주요 2024년 폭풍 동안 발생한 열권 가열의 위도별 범위를 평가하였다.

주요 결과

• 주요 동인으로서의 PC 지수: 폭풍 시 고도 변화율은 PC 지수와 강한 상관관계($r \gtrsim 0.8$)를 보였다. 결정적으로, PC 지수는 고도 반응보다 약 2~3시간 앞서 나타났다(결합 역분산 가중 평균 지연 시간 $\tau = -2.2 \pm 2.6$ 시간; 10월 이벤트의 경우 $\tau = -3.17 \pm 3.17$ 시간). 이 선행 시간은 110–130 km에서 퇴적된 줄 가열(Joule heating)이 약 500 km의 열권까지 수직 및 자오면 방향으로 전파되는 데 필요한 물리적 시간 규모와 일치한다.
• 결과로서의 Dst 지수: 고도 변화율은 Dst 지수의 크기($|Dst|$, $r \approx 0.82–0.84$)와도 강한 상관관계를 보였으나, 시간적 순서가 역전되어 있었다. 고도 하강이 Dst 최솟값보다 약 3시간 앞서 발생했다. 이는 DST-고도 간의 상관관계가 환 전류(Ring current)로부터의 직접적인 열적 결합이 아니라, 공통된 태양풍 구동원으로부터 기인했음을 나타낸다. 환 전류는 열권보다 수십만 킬로미터 위에 존재하기 때문이다.
• 위도 의존성:
  • 2024년 5월 (Gannon Storm): NinjaSat와 XRISM 모두 유의미한 궤도 감쇠를 보였으며, 이는 가열이 저위도($<30^\circ$)까지 확장되었음을 의미한다.
  • 2024년 10월: NinjaSat만 유의미하게 반응하였고, XRISM은 초세속적 추세에서 벗어나지 않았다. 이는 가열이 극지 캡과 오로라 존에 국한될 수 있음을 보여주며, 폭풍의 영향이 동일한 Dst 최솟값을 가진 사건 중에서도 균일하지 않고 극 지역에 국지화될 수 있음을 입증한다.
• 에너지론: 폭풍 기간 동안 PC 기반의 반구 줄 가열 전력($Q_J$)은 동시기 EUV 가열 베이스라인보다 2.5배에서 3배 더 높았으며, 이는 지자기 줄 가열이 해당 사건 동안 지배적인 증분 에너지원임을 확립한다.

의의 및 주장
본 논문은 고위도 이온층의 줄 가열이 폭풍 시 저궤도 궤도 감쇠를 결정하는 핵심 결합 채널임을 최초로 정량적으로 입증하였다고 주장한다. 연구 결과는 Dst 또는 3시간 단위 Kp/Ap 지수를 사용하는 운영 모델이 충격적 가열을 평균화하여 누락시키고 실제 대기 반응보다 뒤처지기 때문에, 항력 예측에 있어 Dst의 유용성에 의문을 제기한다.

본 연구는 운영상의 항력 예측을 개선하기 위한 구체적인 경로를 제안한다. 즉, Kp, Ap, 또는 Dst 프록시를 실시간 1분 단위의 PC 지수로 대체하는 것이다. PC 지수가 궤도 반응보다 2~3시간 앞서기 때문에, 위성 운영자가 밀도 급증이 도래하기 전에 충돌 평가를 재계산하고 수정 기동을 계획할 수 있는 자연스러운 운영 창(Operational window)을 제공한다. 저자들은 결과가 분석된 두 사건에 대해서는 견고하지만, 지연 시간과 진폭 의존성의 보편성을 완전히 확립하기 위해서는 태양 주기 25의 하강 단계 동안 다수의 폭풍에 대한 분석이 필요함을 강조한다.