Filling the gap in the IERS C01 polar motion series in 1858.9-1860.9

Este artigo apresenta e compara duas abordagens, um modelo astronômico paramétrico e um modelo baseado em Análise de Espectro Singular (SSA), para preencher a lacuna de dois anos na série de movimento polar C01 do IERS entre 1858,9 e 1860,9, concluindo que o método SSA é preferível por se basear em um modelo de movimento polar mais completo.

O essencial

Imagine que a Terra não é uma bola de bilhar perfeitamente estática, mas sim um pião gigante e um pouco desequilibrado que gira no espaço. Esse pião "balança" enquanto gira. Esse balanço é chamado de movimento polar.

Os cientistas medem esse balanço há mais de 170 anos, criando um diário de bordo chamado série IERS C01. Esse diário é o registro mais longo e confiável que temos de como a Terra gira. É como se fosse a "história da vida" do nosso planeta.

No entanto, ao olhar para esse diário histórico, os autores deste artigo (Malkin, Golyandina e Olenev) notaram um problema chato: faltam duas páginas inteiras.

Entre os anos de 1858 e 1860, a coordenada "Y" (uma das direções do balanço) desapareceu. Por que? Porque naquela época, o observatório em Washington (que era essencial para medir essa direção específica) parou de observar por um tempo. Como as outras estações (Greenwich e Pulkovo) estavam muito próximas em termos de longitude, elas não conseguiam "enxergar" o balanço na direção Y sozinhas.

O resultado? Um buraco de dois anos em um registro que deveria ser contínuo. Para cientistas que querem estudar padrões de longo prazo (como mudanças climáticas ou a estrutura interna da Terra), ter um buraco no meio da linha do tempo é como tentar ouvir uma música clássica inteira, mas com dois compassos de silêncio no meio. Isso atrapalha a análise.

A Missão: Preencher o Buraco

O objetivo deste artigo foi tentar reconstruir essas duas páginas perdidas de forma inteligente, sem inventar dados do nada. Eles usaram duas abordagens diferentes, como se fossem dois detetives com métodos distintos:

1. O Detetive da "Fórmula Matemática" (Modelo Paramétrico)

Imagine que você sabe que o balanço da Terra é feito de três coisas principais:

• Um "viés" (uma inclinação fixa).
• Um balanço anual (como as estações do ano).
• Um balanço de Chandler (uma oscilação que dura cerca de 14 meses).

O primeiro método tentou encaixar uma fórmula matemática nesses padrões conhecidos. Eles assumiram que, naquela época, a força desses balanços estava crescendo ou diminuindo de forma linear (como uma rampa suave). Eles ajustaram essa fórmula aos dados que existiam antes e depois do buraco e usaram a "curva" resultante para preencher o espaço vazio.

O problema: É como tentar adivinhar o que aconteceu em uma festa baseada apenas na música que tocava antes e depois de um intervalo. Você sabe a melodia, mas pode não capturar as nuances exatas daquele momento específico.

2. O Detetive da "Inteligência de Dados" (Análise SSA)

O segundo método foi mais sofisticado. Em vez de forçar uma fórmula pré-definida, eles usaram uma técnica chamada Análise Espectral Singular (SSA).

Pense nisso como um filtro de ruído ultra-poderoso ou um algoritmo de inteligência artificial que olha para todo o histórico de dados (de 1846 até hoje) e aprende os padrões complexos da Terra.

• Ele não assume que o balanço é perfeitamente regular.
• Ele olha para a "assinatura" completa do movimento da Terra, incluindo pequenas variações que a fórmula simples ignoraria.
• Ele usa a relação entre as duas direções do balanço (X e Y) para "adivinhar" o que deveria estar faltando, baseando-se em como elas se comportam juntas no resto da história.

É como se você tivesse um filme de um rio fluindo, mas com um quadro faltando. O método paramétrico tentaria desenhar a água baseada na velocidade média. O método SSA, por outro lado, analisaria a textura da água, as pedras no fundo e o vento, para reconstruir exatamente como a água deveria estar fluindo naquele quadro faltante.

O Veredito

Os dois métodos funcionaram bem e deram resultados muito parecidos. No entanto, os autores concluíram que o método SSA (o "detetive de dados") é o melhor.

Por quê? Porque ele é mais flexível. Ele não fica preso a regras rígidas e consegue capturar a complexidade real do movimento da Terra, que não é perfeitamente previsível. O resultado final é uma série de dados contínua, sem buracos, que os cientistas podem usar para estudar o passado da Terra com mais precisão.

Por que isso importa?

Preencher esse buraco não é apenas um exercício de matemática.

• Precisão Histórica: Permite que os cientistas estudem como o movimento polar mudou ao longo de 180 anos, algo impossível com um buraco de dois anos no meio.
• Entendimento do Clima: O movimento da Terra está ligado à redistribuição de massa no planeta (como gelo derretendo ou oceanos mudando). Ter dados contínuos ajuda a entender essas mudanças climáticas históricas.
• Evitar Erros: Se alguém analisar esses dados sem perceber o buraco, pode encontrar "fantasmas" no espectro (padrões que não existem, apenas porque os dados estão quebrados).

Em resumo: Os autores pegaram um diário histórico da Terra com duas páginas rasgadas, usaram a matemática e a inteligência de dados para reescrever o que estava escrito nelas, e agora temos uma história completa e contínua do nosso planeta girando.

Resumo técnico

1. O Problema

O Serviço Internacional de Rotação da Terra e Sistemas de Referência (IERS) fornece a série de parâmetros de orientação da Terra (EOP) C01, que é o registro mais longo e confiável da rotação da Terra. A série de movimento polar (PM), iniciada em 1846, é fundamental para investigar variações de longo prazo, como o movimento de Chandler (CW) e o movimento anual (AW).

No entanto, a série C01 apresenta uma lacuna de dois anos na coordenada do polo Yp no intervalo de 1858,9 a 1860,9.

• Causa: Durante este período, o Observatório de Washington (localizado em longitude -77°, essencial para determinar Yp) não realizou observações. As estações restantes (Greenwich e Pulkovo) estão em longitudes próximas (0° e 30°), o que permite calcular Xp com precisão, mas torna impossível estimar Yp com precisão devido à falta de cobertura longitudinal adequada.
• Impacto: A presença de dados faltantes em séries temporais desafia a análise estatística padrão, especialmente análises espectrais que exigem dados espaçados uniformemente e contínuos para evitar a introdução de frequências espúrias.

2. Metodologia

Os autores propuseram e testaram duas abordagens distintas para preencher a lacuna, visando preservar a estrutura da série original:

A. Modelo Astronômico Paramétrico (BCA e BCA2)

• Conceito: Baseado em um modelo físico que assume que o movimento polar é composto por uma tendência (bias), o movimento de Chandler (CW) e o movimento anual (AW).
• Inovação: Reconhecendo que as amplitudes do CW e do AW não são constantes, o modelo permite que essas amplitudes variem linearmente ao longo do tempo.
• Variantes:
  • BCA: Modelo com 9 parâmetros (bias + CW + AW com amplitudes lineares).
  • BCA2: Modelo com 13 parâmetros, que inclui duas componentes distintas para o CW (baseado em picos observados no espectro de 1,181 e 1,231 anos), além do AW.
• Ajuste: Os parâmetros foram ajustados usando o método dos mínimos quadrados (LS) em várias janelas de tempo (de 12 a 24 anos) centradas na lacuna.

B. Modelo Orientado a Dados (SSA - Análise Espectral Singular)

• Conceito: Utiliza a Análise Espectral Singular (SSA), um método não paramétrico robusto para séries temporais complexas.
• Abordagem:
  • 1D-SSA: Aplica-se apenas à série Yp.
  • MSSA (Multivariate SSA): Aplica-se simultaneamente às séries Xp e Yp para explorar a relação física entre elas.
  • CSSA (Complex SSA): Trata a série como um sinal complexo ($X + iY$), o que é fisicamente mais adequado para coordenadas polares, permitindo uma extração de sinal mais precisa.
• Otimização: Os autores testaram diferentes comprimentos de janela ($L$) e ordens de modelo ($r$, número de componentes SVD) criando "lacunas artificiais" na série para validar a precisão do preenchimento. O método CSSA com $L=119$ e $r=9$ (ou combinações similares) mostrou-se o mais preciso.
• Algoritmo: Um processo iterativo de preenchimento de lacunas foi utilizado, onde valores iniciais são inseridos, o sinal é extraído via SSA, e o processo se repete até a convergência.

3. Contribuições Principais

• Primeira Tentativa Sistemática: Este é o primeiro estudo conhecido a abordar especificamente o preenchimento da lacuna de 1858,9–1860,9 na série IERS C01.
• Comparação de Modelos: Oferece uma comparação rigorosa entre modelos paramétricos (baseados em física) e não paramétricos (baseados em dados), demonstrando que o modelo CSSA é superior por capturar melhor as características de longo prazo sem depender de suposições a priori rígidas sobre a evolução das amplitudes.
• Conjunto de Dados Completo: Fornece os valores preenchidos para as 21 épocas faltantes, permitindo que a série IERS C01 seja tratada como uma série temporal uniformemente espaçada do início (1846) até o presente.

4. Resultados

• Concordância: Os resultados obtidos pelos dois métodos (BCA e SSA) concordam dentro dos erros de incerteza da série original (aproximadamente 0,09 segundos de arco).
• Precisão:
  • O erro quadrático médio (RMSE) do preenchimento foi de 0,08" a 0,09".
  • Este valor é extremamente próximo da incerteza declarada para os dados Yp do século XIX na série IERS C01, validando que o preenchimento não introduz ruído significativo além da precisão inerente das observações históricas.
• Escolha Preferencial: Embora ambos os métodos sejam válidos, os autores consideram o resultado do CSSA mais confiável. Isso porque o modelo paramétrico é sensível a suposições sobre a variação linear das amplitudes em janelas curtas, enquanto o CSSA captura a estrutura global e as variações de longo prazo do movimento polar de forma mais adaptativa.
• Dados Finais: A Tabela 4 do artigo lista os valores preenchidos para Yp, com uma média ponderada sugerida para uso geral.

5. Significado e Conclusão

• Viabilidade Científica: O estudo demonstra que é possível reconstruir dados históricos faltantes com uma precisão compatível com as observações originais, utilizando técnicas modernas de análise de séries temporais.
• Utilidade para Pesquisa: A série preenchida permite a aplicação de métodos de análise espectral e de wavelet que exigem dados uniformemente espaçados, facilitando o estudo de variações de longo prazo (como a periodicidade de 80 anos na amplitude do CW) que seriam distorcidas pela lacuna.
• Limitações e Futuro: Os autores alertam que a precisão é limitada pela qualidade dos dados observacionais do século XIX. Eles sugerem que o próximo passo para melhorar o conhecimento sobre a rotação da Terra seria a mineração e reprocessamento de observações históricas brutas (declinações e latitudes) de arquivos de observatórios do século XIX, especialmente para preencher lacunas anteriores a 1846 e refinar a fase e amplitude do CW no início da série.

Em resumo, o artigo fornece uma solução prática e matematicamente robusta para um problema de dados históricos, permitindo que a comunidade científica utilize a série completa de movimento polar da Terra sem a necessidade de ignorar ou tratar inadequadamente o período de 1858–1860.