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

Este artículo presenta un intento pionero para rellenar el vacío de dos años en la serie de movimiento polar C01 del IERS correspondiente a 1858.9-1860.9, comparando un modelo astronómico paramétrico con un enfoque basado en Análisis Espectral Singular (SSA) y concluyendo que este último es preferible por basarse en un modelo de movimiento polar más completo.

Lo esencial

¡Claro que sí! Imagina que la Tierra es como un trompo gigante que gira sobre sí mismo. Pero, a diferencia de un trompo de juguete perfecto, la Tierra "bambolea" un poco mientras gira. Este bamboleo se llama movimiento polar.

Los científicos llevan más de 150 años midiendo este bamboleo para entender cómo cambia nuestro planeta (por el clima, los glaciares, el núcleo de la Tierra, etc.). Tienen un registro histórico muy famoso y confiable llamado la serie IERS C01.

Sin embargo, este registro tiene un problema: tiene un "agujero".

El problema: El año perdido

Entre los años 1858 y 1860, faltan dos años completos de datos sobre una de las coordenadas del bamboleo (la llamada $Y_p$).

¿Por qué pasó esto? Imagina que tienes tres amigos (observatorios) en diferentes partes del mundo (Londres, San Petersburgo y Washington) que se turnan para medir el bamboleo. En esos dos años, el amigo de Washington se enfermó y no pudo medir. Los otros dos amigos estaban muy cerca entre sí geográficamente, por lo que, aunque medían, no podían ver bien la dirección del bamboleo que faltaba. Fue como intentar adivinar la dirección del viento solo mirando dos ventanas que están una al lado de la otra; te falta una perspectiva clave.

Sin esos datos, la línea de tiempo tiene un hueco. Para los científicos, esto es como intentar escuchar una canción favorita pero que le falte un verso entero; es difícil analizar la melodía completa si hay un silencio forzado.

La misión: Rellenar el hueco

Los autores de este artículo (tres científicos rusos) se propusieron una misión: inventar los datos que faltan de la manera más inteligente posible, para que la serie quede continua y sin saltos.

Para hacerlo, usaron dos "detectives" diferentes:

1. El Detective de la Fórmula (Modelo Paramétrico)

Este primer método es como un músico que sabe la partitura de memoria.

• La idea: Sabemos que el bamboleo de la Tierra tiene ritmos fijos: uno que dura un año (como las estaciones) y otro que dura unos 14 meses (llamado "bamboleo de Chandler").
• El truco: El detective asume que, en esos años perdidos, la música siguió tocando con la misma canción, pero quizás un poco más fuerte o más suave (cambiando de volumen). Usó matemáticas para "cantar" esa parte faltante basándose en lo que pasó justo antes y justo después del hueco.
• Resultado: Es una buena estimación, pero depende de que la "canción" no haya cambiado de estilo de repente.

2. El Detective de la Inteligencia Artificial (Análisis SSA)

Este segundo método es como un restaurador de arte experto o un chef que sabe el sabor exacto de un plato.

• La idea: En lugar de asumir una fórmula fija, este método (llamado Análisis de Espectro Singular o SSA) mira todo el registro histórico como un todo. Analiza patrones, ritmos y formas ocultas que ni siquiera sabíamos que existían.
• El truco: Imagina que tienes un tapiz antiguo con un trozo roto. Este método no solo mira los hilos de los bordes, sino que entiende la textura, los colores y el diseño general de toda la tela para "tejer" el trozo faltante de forma que encaje perfectamente con el resto. Es más flexible y no asume reglas rígidas.
• Resultado: Este método es más completo porque descubre detalles pequeños que el primero podría ignorar.

¿Quién ganó?

Ambos detectives hicieron un trabajo excelente y sus resultados coincidieron bastante bien (como si dos traductores diferentes tradujeran el mismo libro y salieran casi igual).

Sin embargo, los autores dicen que el Detecte de la Inteligencia Artificial (SSA) es el ganador. ¿Por qué? Porque es más "realista". No asume que el bamboleo de la Tierra sigue una fórmula simple; deja que los propios datos le cuenten la historia. Es como si el primer método fuera un mapa dibujado a mano con reglas, y el segundo fuera una fotografía satelital de alta resolución.

El resultado final

Gracias a este trabajo, ahora tenemos una línea de tiempo continua desde 1846 hasta hoy, sin ese molesto agujero de dos años.

¿Por qué es importante?

1. Ciencia más limpia: Ahora los científicos pueden estudiar cómo ha cambiado el bamboleo de la Tierra durante 180 años sin tener que saltarse esos dos años o usar métodos complicados para ignorar el hueco.
2. Entender el clima: Al tener datos continuos, podemos entender mejor cómo el cambio climático y el movimiento de las masas de agua y hielo afectan la rotación de la Tierra a largo plazo.
3. Precisión: Los datos rellenados son tan precisos que el error es casi el mismo que el de las mediciones originales de la época. Es como si hubiéramos encontrado un fragmento de un rompecabezas que encaja perfectamente.

En resumen: Los científicos encontraron una forma brillante de "reparar" un trozo de la historia de la Tierra que se había perdido hace más de 160 años, permitiéndonos ver el pasado de nuestro planeta con una claridad mucho mayor.

Resumen técnico

Título: Llenado del vacío en la serie de movimiento polar IERS C01 de 1858.9 a 1860.9

1. El Problema

El Servicio Internacional de Rotación Terrestre y Sistemas de Referencia (IERS) proporciona la serie C01 de Parámetros de Orientación Terrestre (EOP), que constituye el registro más largo y fiable de la rotación de la Tierra. Específicamente, la serie de movimiento polar (PM) comienza en 1846 y es fundamental para estudiar variaciones a largo plazo, como el bamboleo de Chandler (CW) y el bamboleo anual (AW).

Sin embargo, la serie IERS C01 presenta una brecha de dos años en la coordenada de polo $Y_p$ entre los años 1858.9 y 1860.9.

• Causa: Durante este periodo, el Observatorio de Washington no realizó observaciones. Dado que los observatorios de Greenwich y Pulkovo tienen longitudes cercanas (0° y 30°), sus datos solo permiten determinar con precisión la componente $X_p$, pero no $Y_p$.
• Impacto: La falta de datos intermitentes impide el uso de métodos estándar de análisis espectral y estadístico que requieren series temporales uniformemente espaciadas. Esto introduce riesgos de frecuencias espurias y limita la capacidad de analizar la evolución de los modos de rotación terrestre en esa época crítica del siglo XIX.

2. Metodología

Los autores compararon dos enfoques distintos para reconstruir los valores faltantes de $Y_p$ en la brecha, buscando preservar la estructura física de la serie:

A. Modelo Astronómico Paramétrico (Método BCA y BCA2)

• Se desarrolló un modelo basado en la suma de componentes físicas: un sesgo (bias), el bamboleo de Chandler (CW) y el bamboleo anual (AW).
• BCA: Asume amplitudes que cambian linealmente con el tiempo para el CW y el AW.
• BCA2: Una variante que incluye dos componentes de frecuencia para el bamboleo de Chandler (1.181 yr y 1.231 yr) para abordar la posible bifurcación espectral observada.
• Procedimiento: Se ajustaron los parámetros del modelo utilizando el método de mínimos cuadrados (LS) sobre intervalos de datos adyacentes a la brecha (de 12 a 24 años de longitud) para predecir los valores faltantes.

B. Modelo Basado en Datos (Análisis Espectral Singular - SSA)

• Se utilizó el Análisis Espectral Singular (SSA), un método no paramétrico potente para series temporales complejas.
• Enfoque: Se aplicaron tres variantes:
  1. 1D-SSA: Análisis de la serie $Y_p$ por separado.
  2. MSSA (Multivariate SSA): Análisis conjunto de las series $X_p$ e $Y_p$ para capturar sus relaciones físicas.
  3. CSSA (Complex SSA): Tratamiento de la serie como un número complejo ($X_p + iY_p$), lo cual es físicamente más coherente para el movimiento polar y permite una extracción de señal más precisa.
• Optimización: Se probaron diferentes longitudes de ventana ($L$) y órdenes de modelo ($r$) utilizando "brechas artificiales" en los datos existentes para minimizar el error de reconstrucción (RMSE y MAE). Se seleccionó la media de las 7 mejores configuraciones de CSSA para el resultado final.

3. Contribuciones Clave

• Primera solución sistemática: Es el primer intento documentado de llenar específicamente la brecha de $Y_p$ en el IERS C01 de 1858.9–1860.9.
• Comparación de modelos: Se establece una comparación rigurosa entre un modelo físico paramétrico (con suposiciones a priori sobre la linealidad de las amplitudes) y un modelo puramente basado en datos (SSA) que es menos restrictivo.
• Validación de errores: Se demuestra que el error de reconstrucción obtenido es consistente con la incertidumbre inherente de los datos históricos del siglo XIX.
• Descubrimiento de componentes menores: El método SSA permitió identificar componentes de menor amplitud en el espectro del movimiento polar que no eran evidentes en el análisis paramétrico simple.

4. Resultados

• Concordancia: Ambos métodos (BCA y SSA) producen resultados que coinciden dentro de los márgenes de error de la serie IERS C01 (aproximadamente 0.09 segundos de arco).
• Precisión:
  • El error cuadrático medio (RMSE) de la reconstrucción oscila entre 0.08" y 0.09".
  • Este valor es casi idéntico a la incertidumbre declarada en los datos originales de $Y_p$ del IERS C01 para ese periodo, lo que valida la fiabilidad de la reconstrucción.
• Preferencia metodológica: Aunque ambos métodos son válidos, los autores consideran que el enfoque SSA (específicamente CSSA) es preferible.
  • Razón: El modelo paramétrico asume una evolución lineal de las amplitudes en un intervalo corto, mientras que el SSA captura mejor las características a largo plazo y la estructura completa de los datos observados sin depender de suposiciones físicas rígidas.
• Datos generados: Se proporcionan los valores interpolados para los 21 epochas faltantes (Tabla 4 del artículo), permitiendo una serie continua y uniformemente espaciada desde 1846.

5. Significado e Importancia

• Herramienta para la investigación: Proporcionar una serie de datos continua y sin huecos permite a los investigadores utilizar métodos de análisis espectral estándar y realizar análisis detallados "epocha por epocha" sin tener que excluir periodos enteros o usar técnicas de interpolación inadecuadas.
• Estudios de variabilidad a largo plazo: La inclusión de estos datos es crucial para estudiar la variabilidad secular del bamboleo de Chandler (como el ciclo de ~80 años en amplitud y fase) y su relación con cambios climáticos y en la estructura interna de la Tierra.
• Evitar sesgos: Elimina el riesgo de introducir sesgos en los modelos de movimiento polar al ignorar la brecha o tratarla incorrectamente.
• Futuro: Aunque la reconstrucción matemática es exitosa, los autores señalan que la verdadera mejora de la calidad de los datos del siglo XIX requeriría la minería y reprocesamiento de observaciones históricas originales (declinaciones estelares) en archivos de observatorios, un trabajo que excede el alcance de este estudio pero que se identifica como una necesidad futura.

En resumen, el artículo ofrece una solución robusta y validada para completar un eslabón crítico en la historia de la rotación terrestre, facilitando estudios más precisos sobre la dinámica del planeta en el siglo XIX.