Analysis of Short-term Solar Activity Variability and Estimating the Timings of the Next Enhanced Bursts

Resumen Técnico: Análisis de la Variabilidad de la Actividad Solar a Corto Plazo y Estimación de los Momentos de los Próximos Brotes Intensificados

Planteamiento del Problema
La actividad solar exhibe variabilidad en múltiples escalas temporales. Si bien el ciclo solar decenal (CS) está bien caracterizado y los eventos diarios de fulguraciones/eyecciones de masa coronal (CME) son en gran medida impredecibles a corto plazo, existe una brecha crítica en la previsión de la variabilidad de escala intermedia (semanas a meses). El artículo identifica las "estaciones del clima espacial"—brote de actividad intensificada cuasiperiódicos que duran de 6 a 18 meses—como un factor clave que impulsa el clima espacial peligroso (por ejemplo, fulguraciones mayores y eyecciones de masa coronal). La previsión operativa actual se centra en escalas temporales muy cortas (horas a días) o en la amplitud del ciclo decenal, dejando sin abordar los brotes intermedios. Los autores señalan que, aunque los modelos físicos (por ejemplo, modelos de agua somera magnetohidrodinámica de la tachocla y modelos de dínamo) explican la existencia de estos períodos cuasi-anuales mediante mecanismos como las Oscilaciones No Lineales de la Tachocla (TNO), carecen de la precisión de parámetros para predecir el momento específico, la amplitud y la secuencia hemisférica de los brotes individuales.

Metodología
Los autores proponen una estrategia de previsión híbrida novedosa que integra técnicas numéricas, estadísticas y de aprendizaje automático para predecir estos brotes utilizando datos del número de manchas solares (SSN) suavizado mensualmente desde 1878 hasta 2025. La metodología procede en tres pasos distintos:

1. Detección y Caracterización de Brotes (Ajuste Gaussiano):

   • Los datos de SSN hemisféricos se suavizan utilizando un kernel gaussiano de 4 meses para aislar estructuras subanuales mientras se filtra el ruido a escala mensual.
   • Los brotes intensificados se definen como mejoras de lóbulo único que duran de 6 a 18 meses con un SSN pico $\ge$ 10.
   • Se ajusta un perfil gaussiano con un fondo lineal a los brotes identificados para extraer parámetros clave: tiempo pico, amplitud, Ancho a Media Altura (FWHM) y duración.
   • Se utiliza la Estimación de Densidad de Kernel (KDE) para analizar las distribuciones estadísticas de las amplitudes y duraciones de los brotes en los Hemisferios Norte y Sur y en diferentes fases del CS.

2. Previsión Estadística de Línea Base (SARIMA):

   • Se desarrollan modelos separados de Media Móvil Integrada Autoregresiva Estacional (SARIMA) para los Hemisferios Norte y Sur para capturar las tendencias cíclicas subyacentes y los patrones estacionales.
   • Los parámetros del modelo se optimizan utilizando el Criterio de Información de Akaike (AIC).
   • Hemisferio Norte: El modelo óptimo es SARIMA (3, 0, 5) $\times$ (2, 1, 0) con un período estacional de 9 meses.
   • Hemisferio Sur: El modelo óptimo es SARIMA (3, 0, 1) $\times$ (2, 1, 0) con un período estacional de 10 meses.
   • Para tener en cuenta el declive post-pico del CS25, se aplica un factor de decaimiento multiplicativo (reducción del 20% cada 12 meses) a las previsiones de línea base.

3. Inyección de Brotes mediante Aprendizaje Automático (Random Forest):

   • Se entrena un modelo de regresión de Random Forest (RF) en un catálogo de ~300 brotes históricos para predecir la amplitud y la duración del brote.
   • Las características de entrada incluyen SSN mensuales con retraso (1–5 meses antes del brote) y un indicador de hemisferio.
   • El modelo RF genera amplitudes y duraciones predichas, las cuales se inyectan luego como envolturas gaussianas asimétricas sobre las previsiones de línea base SARIMA.

Resultados Clave
El modelo híbrido realiza con éxito la hindcast y la previsión de eventos de brotes para el actual Ciclo Solar 25 (CS25):

• Hemisferio Norte:
  • Hindcast (nov 2024 – may 2025): El modelo identifica un brote con pico en marzo de 2025 con un SSN de $\sim$70.
  • Previsión: Se predice un brote subsiguiente para diciembre de 2025 con una amplitud ligeramente inferior de $\sim$60 SSN.
• Hemisferio Sur:
  • Comportamiento: Muestra un comportamiento más complejo con múltiples amplitudes. Se identifica un brote principal que comienza a finales de 2024, con un pico en marzo-abril de 2025 y una alta amplitud de 130–140 SSN.
  • Previsión: Se predice un brote terciario para enero de 2026 con una amplitud de $\sim$120 SSN.
• Disco Total:
  • Se prevé que el SSN total combinado alcance un pico de $\sim$205–210 en marzo-abril de 2025, seguido de una meseta de $\sim$170–175 hasta finales de 2025 antes de un declive gradual.
• Rendimiento del Modelo:
  • El modelo RF logra una puntuación $R^2$ de 0.748 para la predicción de amplitud, con un Error Absoluto Medio (MAE) de 12.60 SSN.
  • Los modelos SARIMA muestran asimetría hemisférica en los parámetros óptimos, reflejando las fluctuaciones a corto plazo más regulares del Hemisferio Norte frente a los picos puntuales de alta amplitud del Hemisferio Sur.

Significado y Afirmaciones
El artículo afirma abordar una "brecha crítica" en la previsión operativa del clima espacial al proporcionar plazos de semanas a meses para las "estaciones del clima espacial". Los autores afirman que su enfoque híbrido traduce con éxito el concepto motivado físicamente de las TNO en un flujo de trabajo manejable y basado en datos.

Las contribuciones clave destacadas por los autores incluyen:

• Marco Híbrido: El acoplamiento exitoso de líneas base SARIMA con inyecciones de brotes derivadas del aprendizaje automático permite la resolución de brotes cuasi-estacionales (2–20 meses) que las previsiones decenales tradicionales pasan por alto.
• Asimetría Hemisférica: El estudio modela explícitamente los comportamientos dinámicos distintos de los Hemisferios Norte y Sur, señalando que el Hemisferio Sur domina la actividad total de SSN en el ciclo actual.
• Cuantificación de la Incertidumbre: El marco incorpora intervalos de confianza estadísticos (IC 95%) para la línea base SARIMA y estimaciones heurísticas de incertidumbre para los brotes predichos por RF, proporcionando un producto de previsión probabilística.
• Utilidad Operativa: Al predecir el momento y la amplitud de estos brotes intermedios, el método tiene como objetivo proporcionar un plazo suficiente para que la sociedad se prepare para los impactos de la actividad solar intensificada, distinto de la imprevisibilidad de los eventos diarios de fulguraciones.

Los autores mantienen que, aunque los modelos físicos explican el mecanismo de estos brotes, su enfoque híbrido basado en datos es actualmente el método más efectivo para predecir la ocurrencia y la magnitud específicas de estos eventos.