¿Calcular la fecha del equinoccio de primavera, posiblemente usando el sitio web de Horizons de JPL?

Estoy familiarizado con las efemérides de órbita y las formas de usarlas para trazar órbitas y posiciones. Sin embargo, en todos los gráficos que encuentro, parece que no puedo encontrar una manera de determinar cuándo ocurre un equinoccio en otro planeta.

Tomemos como ejemplo los datos que se pueden generar a partir de JPL HORIZONS. Puedo obtener las efemérides y la osculación de la Tierra y Marte. El equinoccio de primavera en la Tierra es fácil, se define como el eje X. Pero dentro de HORIZONTES, parece que no puedo encontrar ningún lugar que determine el "eje X local" para Marte, para determinar el primer día de la primavera.

Ahora, puedo buscar fácilmente en Google "Estaciones de Marte" para determinar la respuesta, pero si quisiera extender este ejercicio para planetas arbitrarios, ¿dónde se almacenan esos datos? ¿Cómo se llama comúnmente ese "desplazamiento" (ángulo entre la referencia X y el X local)?

Para obtener el equinoccio necesitas el ecuador del planeta, lo que significa que necesitas la orientación de su polo. La orientación de los planetas está en Horizons en alguna parte porque puede especificar una posición geográfica en otro planeta y ver otro cuerpo desde él y calcular la altitud/azimut, etc. Pero no estoy seguro de cómo sacar las cosas de allí de tal forma que obtendrías la JD del equinoccio del planeta. Sin embargo, tengo el presentimiento de que alguien podrá responder esto. Gran primera pregunta! +1

Respuestas (1)

Tampoco sé cómo definir un marco inercial en ningún planeta del sistema solar.

Sin embargo, las estaciones también se pueden calcular usando solo el vector polar, que está disponible en el sistema de coordenadas fijas del cuerpo disponible en las efemérides JPL SPICE.

Para hacer eso, se necesita calcular el ángulo entre el vector del polo norte (+Z en el sistema fijo del cuerpo) y el vector del Sol también expresado en el sistema fijo del cuerpo del planeta. Las estaciones se definen entonces como sigue:

  • Ángulo = 90° y decreciente: equinoccio de primavera
  • Ángulo = 90° y creciente: equinoccio de otoño
  • Ángulo en su valor máximo entre dos equinoccios: solsticio de invierno
  • Ángulo en su valor mínimo entre dos equinoccios: solsticio de verano

Hice un pequeño script usando el envoltorio Orekit Python que puede calcular las estaciones en cualquier planeta del sistema solar dentro de un rango de fechas: https://github.com/GorgiAstro/planetary-timemachine/blob/712fe31799faae65b319b2d7ff66685e26ced0f8/seasons-calculation.ipynb The JPL Las efemérides configuradas en este repositorio permiten computar las temporadas desde 1957 hasta 2057.

El guión fue validado contra las estaciones de la Tierra y contra el Calendario de Marte de la Sociedad Planetaria .

Aquí hay una parte de la salida del script para Marte en el rango [2020, 2024] (también se proporciona el ángulo entre el polo norte y el vector solar):

The spring equinoxes are:

    pole_sun_angle_deg
datetime_utc    
2021-02-07  90.097319
2022-12-26  90.089456
2024-11-12  90.084573