¿Por qué los elementos orbitales del LRO parecen oscilar constantemente, excepto recientemente?

Tracé algunos de los elementos orbitales del Lunar Reconnaissance Orbiter de la base de datos Horizons del JPL y veo que hay oscilaciones constantes. El período del semieje mayor y los parámetros relacionados con la excentricidad parecen ser de unos 27,25 días (ver último gráfico), lo que coincide con el período orbital de unos 27,32 días. Sin embargo, las oscilaciones de inclinación parecen oscilar con el doble de esa frecuencia.

¿Por qué los elementos orbitales del LRO parecen oscilar constantemente?

Además, ¿los episodios repentinos de excentricidad muy constante son "reales" o simplemente artefactos de empalmar/coser/pegar diferentes simulaciones juntas? Aun cuando la excentricidad parece constante, se aprecian oscilaciones en la inclinación y en menor grado en el semieje mayor.

editar: acabo de agregar información histórica sobre las diversas trayectorias que se han concatenado. El "período de excentricidad tranquila" es de 2016-Oct-21a 2016-Dec-07así que si bien comienza a mediados de 558day_20160907_01.bsp_V0.2, termina en la misma fecha infame en que también termina ese segmento.

Pero recuerda, estoy preguntando tanto sobre los meneos como sobre su ausencia.

ingrese la descripción de la imagen aquí

SPACECRAFT TRAJECTORY: 
  Updated irregularly (on Horizons) or by request. 
  Concatenated historical (reconstructed) trajectories are from PDS. 

  Trajectory name                 Start (TDB)         Stop (TDB)
  ---------------------------  -----------------  -----------------
  Reconstructed trajectory     2009 Jun 18 22:16  2016 Sep 15 00:01
  558day_20160907_01.bsp_V0.2  2016 Sep 15 00:01  2016 Dec 07 00:01 predict
  558day_20161207_01.bsp_V0.2  2016 Dec 07 00:01  2017 Jan 04 00:01 predict
  558day_20170104_01.bsp_V0.1  2017 Jan 04 00:01  2018 Jul 16 00:01 predict
  558day_20170216_01.bsp_V0.1  2017 Feb 16 00:01  2018 Aug 28 00:01 predict

ingrese la descripción de la imagen aquí

ingrese la descripción de la imagen aquí

Abajo: Excentricidad graficada versus tiempo (días) para dos intervalos desplazados por 327 días mostrando una diferencia de 12 oscilaciones. El período extraído es de 27,25, cercano al período orbital de la luna de unos 27,32 días.

ingrese la descripción de la imagen aquí

Parece que se activó la estabilización activa mediante RCS.
Mi dinero está en una interrupción de telemetría, lo que haría que la onda sinusoidal limpia en los datos de inclinación fuera un artefacto del filtro Kalman.
@Schlusstein La base de datos JPL Horizons contiene órbitas calculadas , el resultado de simulaciones numéricas de trayectorias de naves espaciales basadas en modelos de gravedad y otras fuerzas. Tenga en cuenta que los datos se extienden hasta 2018: ¡no se utilizaron máquinas del tiempo para recuperar la telemetría del futuro! Ahora tomaré ese dinero, ¿cuánto dijiste? :)
@uhoh Sí, y los datos futuros se parecen a los datos sobre los que se pregunta, que es una gran parte de cómo llegué a mi conclusión.
@Schlusstein no hay filtro de Kalman. No puedo entender tu comentario. ¿Cómo puede concluir que hay un filtro de Kalman? Éstos se utilizan típicamente para datos medidos ruidosos para mejorar la interpretación de las mediciones en presencia de cantidades significativas de ruido experimental. Estas son simulaciones numéricas de precisión, espero que cualquier ruido esté más allá del decimoquinto lugar decimal más o menos.
@uhoh ¿Qué datos de entrada se usaron para realizar estas simulaciones? ¿Dónde encontraste los datos que trazaste?
@uhoh, estoy leyendo sobre HORIZONTES. No es simplemente una simulación numérica. Los datos generados a partir de mediciones reales (presuntamente filtradas), con simulación numérica utilizada para llenar los vacíos en los datos. Por supuesto, esto es mucho más preciso que la simulación numérica, ya que permite corregir el error acumulado y las perturbaciones orbitales no gravitacionales. ssd.jpl.nasa.gov/pub/ssd/Horizons_doc.pdf También lo alentaría a obtener más información sobre el filtrado de Kalman, ya que la integración numérica se puede incorporar fácilmente.
@uhoh Y sí, en realidad miré la pregunta. Sin embargo, en los cuatro días que han pasado, olvidé que se trataba de una herramienta web.
@Schlusstein atengámonos a esta pregunta como se hizo. De acuerdo con el resultado de Horizons, iniciar 2016 Sep 15 00:01estos datos es puramente predictivo si entiendo correctamente. Entonces, todo lo que se muestra en el primer conjunto de gráficos a partir de 2016.7 es una integración estrictamente numérica de (la versión relativistamente correcta de) F = metro a . ¿No es eso correcto?
@uhoh No estoy seguro de que eso sea cierto. No sabrá que las regiones de comportamiento extraño terminan precisamente en las fechas finales de estas líneas: "558day_20160907_01.bsp_V0.2 2016 15 de septiembre 00:01 2016 07 de diciembre 00:01 predecir" Creo que es poco probable que sea una coincidencia. Quizás las fechas finales son donde terminan las fases de predicción. Sin embargo, tampoco veo las secciones planas en los datos sin procesar.
@Schlusstein La discusión con usted es muy útil: ¡estábamos haciendo lo mismo al mismo tiempo! Mientras escribía su comentario, actualicé la pregunta con la misma información. Miré el archivo de salida que Horizons guardó en mi disco: en cada línea, justo después del JD y la fecha del calendario, el primer número de punto flotante es la excentricidad, y repentinamente comienza a flotar en aproximadamente 3.23 o 3.24E-02 durante la duración del plano. Fechas spot (JD 2457682.5a 2457729.5).
¿Qué quieres decir con "Estoy preguntando tanto sobre los meneos como sobre su ausencia"? ¿Se refiere esto al comportamiento periódico de los valores de inclinación durante el tiempo en que la excentricidad es esencialmente fija?
@Schlusstein Estoy tan interesado en, por ejemplo, por qué la inclinación del LRO oscila con una amplitud de 1 grado cada 14 días como en por qué esta oscilación cae repentinamente a 0,1 grados por un tiempo y luego comienza de nuevo. Estoy preguntando tanto sobre estas excursiones como sobre sus repentinas paradas y arranques.

Respuestas (2)

JPL HORIZONS presenta órbitas interpoladas a partir de datos reales en lugar de simulaciones puras. Sospecho que lo que sucedió aquí es que no hay datos orbitales para las tres secciones de datos que no exhiben la oscilación. Sabemos con certeza que no hay datos sin procesar para una de esas secciones porque está en el futuro. Hablé con un experto en controles y dijo que mi sospecha probablemente sea correcta.

Cuando usa la herramienta web, dice que hay varias series de datos y todas menos la primera tienen la palabra predecir al final de la línea. Las fechas al final de esas series de datos corresponden a la oscilación devuelta de los datos. Sospecho que se inicia una nueva serie de datos cuando se reanuda el contacto con la nave espacial.

Creo que el período de oscilación de ~14 días corresponde a la mitad del período orbital de la luna. La perturbación gravitacional de la Tierra sería entonces el principal impulsor de este comportamiento. El hecho de que estas oscilaciones ocurran de manera bastante visible en ausencia de datos reales indica que el modelo está haciendo un trabajo decente al dar cuenta de esta perturbación. En cuanto a por qué la amplitud cae considerablemente, no puedo decirlo con certeza ya que no sé lo suficiente sobre el modelo utilizado, pero puedo suponer que está relacionado con la densidad no uniforme no contabilizada de la luna y las perturbaciones no gravitacionales relacionadas con su orientación relativa al sol, como la presión de la luz y la desgasificación.

Si observa los datos de excentricidad, verá una oscilación con un período de ~14 días además de una oscilación con un período de ~27 días, con algunos armónicos más altos dispersos. También hay una fuerte oscilación del período de ~27 días en los datos de periapsis y apoapsis. Puede ser interesante aplicar una transformada de Fourier a todo esto. Podría observar las magnitudes relativas de la perturbación por período y tal vez tratar de aislar algún comportamiento con un período de un año.

Interpolado no es exactamente la palabra correcta, ajustado podría ser mejor. Por lo general , los únicos datos de alta precisión que obtenemos de las naves espaciales que orbitan otros cuerpos son los de retardo doppler. Una señal especialmente codificada/modulada se transmite desde la Red de Espacio Profundo (DSN) a la nave espacial, que la retransmite a la Tierra a otra frecuencia, pero sincronizada cuidadosamente en fase con la portadora entrante. El desplazamiento doppler y el retraso absoluto dan la velocidad y la distancia relativas con respecto a la estación terrestre, pero estos no son vectores de estado verdaderos con el cuerpo orbitado todavía.
Las simulaciones se utilizan para generar resultados de tipo doppler de retardo simulados y se ajustan iterativamente hasta que pueden reproducir los datos de doppler de retardo medidos. A veces, también se pueden incorporar datos de naves espaciales (cámaras, acelerómetros), pero por lo general estos no ofrecen suficiente precisión para igualar los datos de retardo doppler. Las tres estaciones terrestres de DSN están ocupadas día y noche vigilando las naves espaciales. Aquí hay algunos datos geniales de doppler de retardo de una nave espacial alrededor de la luna, pero en este caso es un radar pasivo.

Basado en algunas respuestas muy útiles de Jon Giorgini en JPL y un poco de lectura adicional, puedo resumir de la siguiente manera:

El cálculo real de las órbitas de las naves espaciales, que combina los efectos de la gravedad y otras fuerzas en el Sistema Solar más las maniobras orbitales de la propia nave espacial, se combinan/reconcilian con el radar y la telemetría disponibles. El resultado es, por supuesto, la tabla de vectores de estado.

Los elementos osculadores se derivan de los vectores de estado únicamente por conveniencia. No están destinados a ser fuentes primarias de datos de órbita, y cada conjunto de puntos solo está destinado a definir la posición de la nave espacial en ese momento .

Durante un breve período de tiempo, mientras se esperaban más datos o una decisión sobre una próxima maniobra, parece que se aplicó un algoritmo alternativo para los elementos osculadores a un breve lapso de fechas, lo que dio como resultado una apariencia más suave de algunos de los gráficos. .

Esto no significa que sean necesariamente menos precisos en su posición, porque los elementos solo están destinados a definir la posición cerca del momento para el que se calculan. Para empezar, la verdadera órbita de la nave espacial no es realmente una sección cónica, especialmente considerando el complicado campo de gravedad del sistema Tierra-Luna.

Siguiendo el cálculo actualizado, la apariencia vuelve a ser la esperada, aunque sigue siendo igualmente válida.

En cuanto a las variaciones periódicas en sí: Las periodicidades mensuales y bimensuales son manifestaciones de la realidad de orbitar el campo gravitatorio grumoso de la luna bajo los efectos gravitacionales de la Tierra. Dado que la órbita de la Luna es significativamente elíptica en lugar de circular, ese fuerte efecto perturbador es periódico. Estos efectos se hacen particularmente visibles en estas gráficas porque, para empezar, no son simples órbitas keplerianas, y tratar de expresarlas como tales resaltará la diferencia.

A continuación se muestra una comparación rápida de los elementos osculadores derivados antes y después de la actualización programada periódicamente. Todo está bien.

ingrese la descripción de la imagen aquí