Impact of geomagnetic storm intensity on thermospheric density and orbital decay of LEO satellites: a G1–G5 per-event analysis

기술 요약: 지자기 폭풍 강도가 열권 밀도 및 저궤도(LEO) 위성의 궤도 붕괴에 미치는 영향

문제 정의
저지구 궤도(LEO) 위성 군집의 급격한 확장은 더욱 정확한 열권 밀도 특성화에 대한 운영적 필요성을 고조시켰다. 대기 항력은 고도 600km 미만의 위성에 있어 주요한 비중력 섭동이며, 지자기 폭풍은 열권 밀도를 크게 증가시켜 궤도 붕괴를 가속화한다. 2022년 2월 38개의 스타링크 위성 손실과 2024년 10월 스타링크 위성의 조기 재진입과 같은 역사적 사건들은 심각한 운영 위험을 입증한다. 기존의 통계적 방법들은 종종 중첩 에포크 분석(Superposed Epoch Analysis, SEA)에 의존하는데, 이는 여러 사건의 응답을 평균화함으로써 즉각적인 궤도 예측 및 충돌 회피에 필수적인 날카로운 피크 응답을 완화할 가능성이 있다. 또한, 경험적 열권 모델은 강렬한 폭풍 사건 동안 발생하는 밀도 증가를 정확하게 재현하는 데 어려움을 겪는 경우가 많다.

방법론
본 연구는 태양 주기 23, 24, 25(2003–2024)에 걸친 G1에서 G5(최대 $K_p$ 5에서 9) 범위의 22개 지자기 폭풍으로부터 순간적인 피크 응답을 분리하기 위해 "사건별 피크(per-event peak)" 분석 접근법을 채택한다. 데이터셋에는 2024년 5월 개넌(Gannon) G5 폭풍이 포함된다.

• 데이터 소스: 인시츄(in-situ) 중성 질량 밀도 및 정밀 궤도 데이터는 GRACE(2002–2017) 및 GRACE-FO(2018–현재) 쌍둥이 위성의 가속도계 측정값으로부터 도출되었다. 지자기 지수($K_p$, $Dst$, $SYM-H$, $AU$, $AL$) 및 태양풍 매개변수는 NASA/GSFC OMNIWeb에서 확보하였다.
• 사건 선정: 다음 세 가지 기준에 따라 22개의 사건이 선정되었다: (1) 최대 $K_p \ge 5$ (G1 이상); (2) 폭풍 돌발 시작(SSC) 전 48시간부터 폭풍 피크 후 48시간까지의 연속적인 가속도계 데이터 가용성; (3) 중첩되는 선행 폭풍 메인 페이즈 없음. 선정된 데이터에는 G1–G4 범주별로 5개씩의 사건과 사용 가능한 두 개의 G5 사건(2003년 10월 및 2024년 5월)이 포함된다.
• 처리 과정:
  1. 궤도 선정: 각 사건에 대해 "폭풍 궤도"는 폭풍 피크(최소 $Dst$)에 가장 가까운 궤도로 정의되었으며, "정온 참조 궤도(quiet reference orbit)"는 최소 24시간 동안 $K_p \le 3$이 유지된 후 SSC에 가장 가까운 궤도로 정의되었다.
  2. 붕괴율 계산: 원시 10초 고도 데이터를 궤도당 평균한 후, 반일 주기 변동을 제거하기 위해 24시간 이동 평균을 사용하여 평활화하였다. 정온 기간과 폭풍 기간에 선형 적합(linear fit)을 적용하여 고도 붕괴율($R_{quiet}$ 및 $R_{storm}$)을 결정하였다.
  3. 강화 지표: 폭풍 시 궤도 붕괴 강화는 백분율 증가량으로 계산되었다: $\frac{R_{storm} - R_{quiet}}{R_{quiet}} \times 100$.
  4. 통계 모델링: 세 가지 접근 가능한 예측 변수인 정온 시 붕괴율, 위성 운용 고도, 최대 $K_p$를 사용하여 폭풍 시 붕괴율을 예측하는 경험적 다중 선형 회귀(MLR) 모델을 개발하였다.

주요 결과

• 강도에 따른 단조적 스케일링: 평균 궤도 붕괴율 강화는 폭풍 강도에 따라 단조롭게 증가한다. 관찰된 강화 수치는 다음과 같다: 77.2% (G1), 98.5% (G2), 138.1% (G3), 149.6% (G4), 166.5% (G5).
• 임계 조절 인자로서의 고도: 위성 운용 고도는 독립적이고 유의미한 요인으로 작용한다. 저고도 사건(예: ~350 km)은 상대적으로 약한 지자기 강제력 하에서도 강한 항력 응답을 보였으며, 종종 더 높은 고도의 고강도 폭풍의 평균 강화 수치를 초과하였다.
• 사건별 변동성: 강도 범주 내에서 상당한 산포가 존재한다. 특히, 두 개의 G5 사건(2003년 10월 및 2024년 5월)은 동일한 $K_p$ 값과 유사한 고도에도 불구하고 약 2.6배의 붕키 붕괴 강화 차이를 보였다. 저자들은 2024년 개넌 폭풍의 더 높은 강화 수치를 더 낮은 배경 태양 활동(낮은 정온 시 붕괴율 분모)과 짧고 강렬했던 2003년 할로윈 폭풍에 비해 연속적인 CME(>30시간)에 의해 유도된 장기 지속 시간 때문이라고 설명한다.
• MLR 모델 성능: 경험적 모델 $D_s = 1.55D_q - 0.29(h_0 - 450) + 7.39Kp_{max} - 35.60$은 22개 사건 전체에 대해 0.91의 결정 계수($R^2$)와 12.3 m/day의 평균 제곱근 오차(RMSE)를 달s성하였다. 이 모델은 위성 고유의 물리적 매개변수(질량, 단면적, 항력 계수)를 요구하지 않고도 폭풍 강도와 고도에 대한 의존성을 성공적으로 포착한다.

의의 및 주장
본 논문은 SEA 기반 분석에 대한 보완적 벤치마크로서 본 연구의 위치를 설정한다. 평균화 과정을 통해 흐려지는 날카로운 신호를 보존하기 위해 사건별 피크 응답에 집중함으로써, 경험적 열권 모델에 대한 보다 직접적인 검증 참조를 제공한다. 본 연구는 위성 운영자가 표준 추적 데이터(고도, 정온 시 붕괴율)와 우주 기상 지수($K_p$)만을 사용하여 폭풍 시 붕괴율을 추정할 수 있게 하는 단순하고 운영적으로 접근 가능한 MLR 도구와 데이터셋을 제공한다. 저자들은 이러한 결과가 확장되는 저지구 궤도 인구, 특히 현재의 태양 주기 25 맥락에서 개선된 궤도 예측 및 충돌 회피 전략 개발을 지원한다고 주장한다.

저자들은 MLR이 유사한 물리적 특성을 가진 두 개의 특정 위성 플랫폼(GRACE/GRACE-FO)에서 도출되었다는 점을 언급하며 모델의 일반화 가능성에 대해 신중한 범위를 유지한다. 그들은 도구의 적용 범위를 더 다양하고 광범위한 LEO 위성 집단으로 확장하기 위해 위성 질량, 단면적 및 속도를 통합하는 향후 연구가 필요함을 인정한다.