Self-consistent 1D modelling of Jupiter's upper atmosphere as an exoplanet analogue

Resumen Técnico: Modelado 1D autoconsistente de la atmósfera superior de Júpiter como un análogo de exoplaneta

Planteamiento del problema
Comprender la estructura térmica y química de la atmósfera superior de Júpiter es crítico no solo para interpretar las observaciones del Sistema Solar, sino también para modelar las atmósferas de exoplanetas gigantes cercanos, como los Júpiter calientes. Si bien Júpiter sirve como el análogo más cercano en masa, radio y composición para el cual existen mediciones fiables, su atmósfera superior presenta interrogantes no resueltos respecto a la energética y la composición. Los modelos existentes enfrentan limitaciones: los modelos fotoquímicos 1D a menudo dependen de perfiles de temperatura prescritos, mientras que los modelos de Equilibrio Radiativo-Convectivo (RCE) no logran contabilizar la química fuera del equilibrio, la conducción térmica y fuentes de calentamiento específicas como el calentamiento por efecto Joule. Por el contrario, los Modelos de Circulación General (GCM) 3D capturan el transporte horizontal pero a menudo carecen de la viabilidad computacional para incluir redes fotoquímicas detalladas y funciones de enfriamiento microfísico por encima de la homopausa. Existe una falta de un modelo radial promedio, autoconsistente, de los perfiles térmico y químico de Júpiter que pueda servir como una base validada para análogos de exoplanetas.

Metodología
Los autores extendieron el código Kompot, un modelo 1D de primeros principios desarrollado originalmente para cuerpos rocosos, para simular atmósferas ricas en hidrógeno de gigantes gaseosos. El modelo resuelve ecuaciones hidrodinámicas, de balance de energía y de cinética fotoquímica acopladas en una rejilla radial 1D que abarca desde la capa de neblina estratosférica (10⁻³ bar) hasta la exosfera.

Los componentes metodológicos clave incluyen:

• Hidrodinámica: El modelo asume un caso hidrostático (velocidad de advección vertical establecida en cero) y calcula los perfiles de densidad integrando desde un límite inferior definido por las mediciones de la sonda Galileo (165 K, 4.36×10¹⁶ cm⁻³). La exosfera se calcula dinámicamente donde el camino libre medio es igual a la escala de altura.
• Química: Una red química exhaustiva incluye reacciones neutras, iónicas y fotoquímicas hasta C₂H₆. Los coeficientes de velocidad se extraen de la Base de Datos Cinética para la Astroquímica (KIDA) y otras redes especializadas (p. ej., VULCAN, Yung et al.). El sistema se resuelve utilizando un integrador semiimplícito de tipo Rosenbrock.
• Balance de Energía: La estructura térmica se determina mediante el equilibrio de las fuentes de calentamiento y enfriamiento.
  • Calentamiento: Incluye radiación solar XUV, calentamiento por fotoelectrones, calentamiento por reacción química y calentamiento por efecto Joule. El calentamiento por efecto Joule se calcula de forma autoconsistente utilizando la conductividad de Pedersen derivada de las frecuencias de ciclotrón iónico y de colisión, escalada a un parámetro de energía de entrada total. Notablemente, el modelo excluye el calentamiento por precipitación de partículas prevalente en las regiones aurorales, centráéndose en un perfil promedio global.
  • Enfriamiento: Dominado por el enfriamiento de la línea de H₃⁺ (por encima de la homopausa) y el enfriamiento radiativo de CH₄ (por debajo de la homopausa). Términos de enfriamiento adicionales incluyen CO₂, NO, O, Lyman-$\alpha$ y H₂O.
• Difusión: El modelo equilibra la difusión molecular y la de remolino (eddy diffusion). Se realizaron cuatro ejecuciones de modelo específicas variando el coeficiente de difusión de remolino ($K_{zz}$): dos utilizando perfiles de Moses et al. (2005) (Modelos A y C) y dos utilizando valores constantes ($10^6$ y $10^7$ cm²/s).
• Transferencia Radiativa: Para generar un espectro de transmisión, los perfiles térmico y químico del Modelo A se introdujeron en el código de transferencia radiativa TauREx, utilizando datos de opacidad de ExoMol, HITRAN y HITEMP.

Resultos Clave

• Estructura Térmica: El modelo reproduce con éxito la inversión térmica observada y la estructura de temperatura de la atmósfera superior de Júpiter. El calentamiento por efecto Joule se identifica como la fuente de calentamiento dominante en la mayor parte de la atmósfera superior (representando ~200 TW), mientras que el calentamiento químico y XUV se vuelven relevantes cerca de la homopausa. El modelo reproduce la transición donde el enfriamiento por H₃⁺ toma el relevo del enfriamiento por CH₄ como el mecanismo de enfriamiento dominante por encima de la homopausa.
• Perfiles Químicos: El modelo reproduce las razones de mezcla volumétrica observadas para CH₄ y otros hidrocarburos (C₂H₂, C₂H₄, C₂H₆). Los perfiles verticales muestran el comportamiento esperado, como la formación máxima de especies portadoras de C₂ cerca de la homopausa debido a la fotólisis de CH₄ y la subsiguiente recombinación. Los resultados muestran una fuerte concordancia tanto con las mediciones in situ de Galileo como con los datos de telescopios terrestres y modelos fotoquímicos existentes.
• Sensibilidad: El estudio destaca que los coeficientes de difusión de remolino influyen significativamente en las abundancias químicas y los perfiles térmicos. Las variaciones en el ángulo cenital y en la entrada de calentamiento por efecto Joule (dentro de límites razonables) resultan en cambios de temperatura de aproximadamente 100 K, lo cual permanece dentro del rango de la variabilidad espacial observada en Júpiter.
• Espectro de Transmisión: El espectro de transmisión sintético (1–15 µm) está dominado por la absorción de CH₄ en el infrarrojo cercano y medio. Aunque las características de C₂H₂ son prominentes a 13.7 µm, otras características menores de hidrocarburos (p. ej., C₂H₄) son difíciles de detectar de forma independiente debido a la proximidad a las bandas más fuertes de CH₄ o a los niveles de ruido. El modelo señala que no se incluyeron la dispersión de aerosoles y nubes de hielo, lo que limita la replicación de ciertas características en el infrarrojo cercano atribuidas al hielo de H₂O cristalino en las observaciones.

Significancia y Conclusiones
El artículo presenta la primera aplicación 1D autoconsistente del modelo Kompot a la atmósfera superior de Júpiter, tratando al planeta como un análogo de exoplaneta orbitando una estrella similar al Sol. La contribución principal es un marco de referencia validado que acopla hidrodinámica, fotoquímica detallada y balance de energía sin depender de perfiles de temperatura prescritos.

Los autores concluyen que:

1. El modelo reproduce eficazmente los perfiles químicos y térmicos promedio consistentes con las mediciones y los modelos de la literatura.
2. El calentamiento por efecto Joule es el principal motor de la estructura térmica de la atmósfera superior, mientras que el enfriamiento radiativo por CH₄ y H₃⁺ dicta el régimen de enfriamiento.
3. La difusión de remolino sigue siendo un parámetro libre crítico que altera significativamente los resultados del modelo, lo que subraya la necesidad de un manejo cuidadoso en el modelado atmosférico.
4. Este marco proporciona una base físicamente consistente para futuras simulaciones de diversas atmósferas de exoplanetas ricos en hidrógeno, incluyendo Júpiteres calientes y templados, con trabajos futuros planificados para ajustar la red química a temperaturas más altas.

El trabajo reconoce limitaciones, específicamente la exclusión de la precipitación de partículas (calentamiento auroral) y la microfísica de aerosoles, lo que restringe su aplicación a regiones globales promedio, no aurorales, y limita la replicación completa de ciertas características espectrales.