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

Resumo Técnico: Resposta em Escala de Hora da Altitude de Satélites em Órbita Terrestre Baixa a Tempestades Geomagnéticas Extremas

Definição do Problema
O crescimento rápido da população de satélites em Órbita Terrestre Baixa (LEO) exige previsões precisas e de baixa latência da densidade termosférica para garantir a prevenção de colisões e a sustentabilidade orbital. Embora o decaimento orbital secular seja bem compreendido, os aumentos impulsivos de densidade impulsionados por tempestades geomagnéticas permanecem mal resolvidos em escalas de tempo de ordem de hora. Modelos operacionais convencionais (ex: NRLMSISE-00, DTM-2013, JB2008) dependem de índices geomagnéticos com resolução temporal grosseira (Kp/Ap de 3 horas ou Dst de 1 hora) ou falham em capturar os mecanismos causais específicos do aquecimento durante tempestades. Uma questão crítica ainda não resolvida é qual observável geomagnética serve como o proxy causal mais direto para o aquecimento impulsivo que impulsiona o rápido decaimento orbital, e qual é o timescale de transporte específico do aporte de energia para a expansão atmosférica.

Metodologia
O estudo utiliza telemetria de GNSS de alta cadência de dois satélites de astronomia de raios X operando durante o máximo solar de 2024:

1. NinjaSat: Um CubeSat 6U em órbita heliosíncrona (inclinação de $97^\circ$, altitude de $\sim$500 km) fornecendo dados de posição com cadência de 10 minutos.
2. XRISM: Um satélite da Missão de Imagem e Espectroscopia de Raios X em órbita de baixa inclinação ($31^\circ$, altitude de $\sim$560 km) fornecendo dados com cadência de 30 segundos.

Para isolar as perturbações induzidas por tempestades do fundo (background), os autores construíram um resíduo de altitude ($\Delta a$) definido como a diferença entre a altitude derivada por GNSS ($a_{GNSS}$) e a altitude propagada a partir de Elementos de Duas Linhas (TLEs) do NORAD usando o modelo SGP4 ($a_{SGP4}$). Como o SGP4 captura o decaimento secular, mas não picos impulsivos de densidade, um $\Delta a$ positivo indica aumento de arrasto.

A análise envolveu:

• Remoção de Baseline: Uma regressão linear da taxa de mudança de altitude suavizada ($\dot{g}\Delta a$) contra a irradiância EUV do GOES-18 (140,5 nm) durante períodos geomagneticamente calmos (julho–setembro de 2024) para modelar e subtrair o background impulsionado por EUV.
• Correlação Cruzada: A taxa de mudança de altitude residual do tempo de tempestade foi correlacionada com o índice do Cap Polar (PC) (um proxy para o aporte de energia do vento solar para a ionosfera de alta latitude) e o índice de Distúrbio de Tempo de Tempestade (Dst) (um proxy para a corrente de anel).
• Comparação Latitudinal: Ao comparar a resposta do NinjaSat (cobertura polar) e do XRISM (cobertura de baixa latitude), o estudo avaliou a extensão latitudinal do aquecimento termosférico durante duas grandes tempestades de 2024: o evento "Gannon" de maio e o evento de outubro.

Principais Resultados

• Índice PC como o Principal Driver: A taxa de mudança de altitude durante a tempestade correlaciona-se fortemente com o índice PC ($r \gtrsim 0,8$). Crucialmente, o índice PC antecede a resposta de altitude em aproximadamente 2 a 3 horas (média de lag ponderada pelo inverso da variância de $\tau = -2,2 \pm 2,6$ h; especificamente no evento de outubro $\tau = -3,17 \pm 3,17$ h). Este tempo de antecedência alinha-se com o timescale físico necessário para que o aquecimento Joule depositado a 110–130 km se propague vertical e meridionalmente até os $\sim$500 km da termosfera.
• Índice Dst como uma Consequência, Não uma Causa: Embora a taxa de mudança de altitude também se correlacione fortemente com a magnitude do índice Dst ($|Dst|$, $r \approx 0,82–0,84$), a ordenação temporal é invertida. A queda de altitude precede o mínimo de Dst em $\sim$3 horas. Isso indica que a correlação Dst-altitude surge de um driver comum do vento solar, e não de um acoplamento térmico direto da corrente de anel, que reside centenas de milhares de quilômetros acima da termosfera.
• Dependência Latitudinal:
  • Maio de 2024 (Tempestade Gannon): Tanto o NinjaSat quanto o XRISM exibiram decaimento orbital significativo, indicando que o aquecimento se estendeu para baixas latitudes ($<30^\circ$).
  • Outubro de 2024: Apenas o NinjaSat respondeu significativamente; o XRISM não mostrou desvio de sua tendência secular. Isso sugere que o aquecimento foi confinado ao cap polar e à zona auroral, demonstrando que os impactos de tempestades não são uniformes e podem ser localizados nos polos mesmo durante eventos com mínimos de Dst comparáveis.
• Energética: O poder de aquecimento Joule ($Q_J$) hemisférico derivado do PC durante as tempestades excedeu o baseline de aquecimento EUV concorrente por fatores de 2,5 a 3, estabelecendo o aquecimento Joule geomagnético como a fonte de energia incremental dominante durante esses eventos.

Significância e Alegações
O artigo afirma fornecer a primeira demonstração quantitativa in-orbit de que o aquecimento Joule da ionosfera de alta latitude é o canal de acoplamento decisivo para o decaimento orbital de LEO durante tempestades. Os achados desafiam a utilidade do índice Dst como um preditor oportuno para o arrasto, observando que modelos operacionais que utilizam Dst ou índices Kp/Ap de 3 horas subestimam sistematicamente os aumentos de densidade porque suavizam o aquecimento impulsivo e atrasam a resposta atmosférica real.

O estudo propõe um caminho concreto para melhorar a previsão operacional de arrasto: substituir os proxies Kp, Ap ou Dst pelo índice PC de cadência de 1 minuto em tempo real. Como o índice PC antecede a resposta orbital em 2–3 horas, ele oferece uma janela operacional natural para operadores de satélites recomputarem avaliações de conjunção e planejarem manobras corretivas antes da chegada do pico de densidade. Os autores enfatizam que, embora os resultados sejam robustos para os dois eventos analisados, uma análise de múltiplas tempestades ao longo da fase descendente do Ciclo Solar 25 é necessária para estabelecer plenamente a universalidade dos tempos de lag e das dependências de amplitude.