Predicting Atmospheric Re-entries Using Segmented M/A Calibration of Public TLE Data

Resumo Técnico: Predição de Reentradas Atmosféricas Usando Calibração Segmentada de M/A de Dados TLE Públicos

Definição do Problema
A predição do tempo preciso de reentradas atmosféricas para objetos de baixa altitude é dificultada pelo ruído intrínseco e pela cadência dos dados de Elementos de Duas Linhas (TLE) publicamente disponíveis. O ajuste de arrasto de arco longo tradicional frequentemente diverge porque acumula o ruído de catálogo de curto período, obscurecendo a verdadeira evolução dinâmica do objeto. Além disso, o estado físico de um objeto em reentrada (atitude, tombamento, condições térmicas e coeficiente balístico) pode variar em escalas de tempo de minutos, tornando a modelagem física detalhada ab-initio obsoleta antes mesmo de sua conclusão. O desafio reside em extrair um sinal dinâmico estável de observações TLE ruidosas e discretas sem exigir conhecimento detalhado da geometria ou do estado aerodinâmico do objeto.

Metodologia: Calibração Segmentada de M/A
O artigo propõe um método empírico e reproduzível que reconstrói a trajetória de decaimento terminal calibrando a razão efetiva massa-área ($M/A$) ao longo de sucessivos arcos de propagação curtos. Os princípios e etapas centrais são:

• Abordagem Segmentada: Em vez de um único ajuste de arco longo, o método utiliza pares de TLEs sucessivos para definir segmentos independentes. O TLE anterior é propagado até a época do TLE posterior.
• Calibração de Parâmetros: A área de arrasto efetiva ($A$) é ajustada para que o decaimento do semieixo maior simulado ($\Delta a_{sim}$) corresponda ao decaimento observado ($\Delta a_{obs}$) entre os dois TLEs. A área ideal é determinada via interpolação linear, resultando em uma razão $M/A$ efetiva para aquele intervalo específico.
• Supressão de Ruído: Para isolar o decaimento secular do ruído de catálogo, impõe-se um espaçamento mínimo de 12 horas entre os TLEs para a calibração padrão. Na fase terminal, onde a queda do semieixo maior torna-se grande o suficiente para dominar as flutuações, os segmentos são aceitos se o semieixo maior diminuir pelo menos 10 km, mesmo com espaçamentos mais próximos.
• Modelo de Força: As propagações são realizadas usando o GMAT da NASA com um modelo de força fixo (gravidade EGM96, densidade NRLMSISE-00, pressão de radiação solar e parâmetros de orientação da Terra). O único parâmetro ajustado é a área efetiva $A$, que atua como um proxy para o comportamento líquido de arrasto, absorvendo variações na geometria, atitude e estado térmico.
• Fonte de Dados: O método baseia-se exclusivamente em dados públicos: TLEs sucessivos e mensagens finais de Previsão de Rastreamento e Impacto (TIP) do Space-Track.

Principais Resultados
O método foi aplicado a reentradas naturais recentes de corpos de foguetes, detritos e satélites não manobráveis, bem como casos específicos da Starlink:

• Precisão: Na fase terminal, o método reproduz consistentemente os tempos TIP publicados com uma precisão de aproximadamente $\pm 1$ hora. O desvio relativo em relação aos TIPs finais frequentemente reduz em um fator de dois a quatro nas últimas 6 a 12 horas de vida orbital em comparação com estimativas anteriores.
• Identificação da Órbita Final: A abordagem identifica de forma confiável a última órbita completa antes da entrada atmosférica, permitindo a exclusão de todas as regiões fora da trajetória final com certeza.
• Estudos de Caso:
  • Decaimento Natural: Para objetos como NORAD 43802 e 61440, o método convergiu suavemente para o TIP final.
  • Variações da Starlink: O método rastreou com sucesso o decaimento natural (STARLINK-2171) e eventos de propulsão (STARLINK-1363). No entanto, para manobras de deorbitagem rápida (STARLINK-34110), a predição em tempo real foi impossível durante a fase de propulsão; os valores de $M/A$ calibrados apresentaram picos (indicando propulsão em vez de arrasto), e a reconstrução precisa só foi possível retrospectivamente após a conclusão da manobra.
• Validação Independente: Várias predições foram confirmadas por observações independentes, incluindo um avistamento da ISS da reentrada do PROGRESS MS-34, relatos de testemunhas no Sri Lanka para o ICEYE-X26 e 91 relatórios da American Meteor Society para o CLARITY-1.

Limitações
O artigo define explicitamente o envelope operacional e as restrições do método:

• Restrições de Catálogo: O principal fator limitante é o ruído intrínseco, a quantização e a cadência de publicação do catálogo TLE, não o modelo atmosférico ou de fluxo solar.
• Manobras: Manobras controladas ou episódios de propulsão invalidam temporariamente as predições de curto prazo porque a evolução do semieixo maior reflete a propulsão em vez do arrasto. As predições recuperam a confiabilidade apenas quando o objeto retorna a um regime puramente balístico.
• Teto de Precisão: Mesmo sob condições favoráveis, o tempo de decaimento final pode tipicamente ser recuperado com uma precisão de cerca de $\pm 1$ hora. Este é um limite fundamental imposto pelo conteúdo de informação dos dados TLE.

Significância e Alegações
O artigo afirma que a calibração segmentada de $M/A$ oferece uma alternativa transparente, leve e operacionalmente robusta ao ajuste de arrasto de arco longo quando apenas dados de rastreamento públicos estão disponíveis. Sua significância reside em:

1. Robustez Empírica: Ao ancorar cada segmento a uma nova medição de catálogo, o método evita a divergência de arco longo e isola o ruído de catálogo como a incerteza dominante.
2. Praticidade: Não requer conhecimento detalhado do estado aerodinâmico ou da geometria do objeto, tornando-o aplicável a uma ampla gama de objetos (detritos, corpos de foguetes, satélites) onde tais dados são desconhecidos.
3. Reprodutibilidade: O método fornece um framework reproduzível para identificação de eventos e reconstrução rápida de decaimento, adequado para verificação independente de reentradas relatadas publicamente.

O autor conclui que, embora o método não possa superar as limitações intrínsecas do catálogo TLE, ele extrai efetivamente o máximo sinal dinâmico possível dos dados públicos, fornecendo precisão na escala de horas na fase terminal crítica de reentrada.