Narrowing the solar surface flux transport parameter space through nonlinear feedbacks

Resumen Técnico: Estrechamiento del espacio de parámetros del transporte de flujo en la superficie solar mediante retroalimentaciones no lineales

Planteamiento del problema
El modelo de Transporte de Flujo Superficial (SFT, por sus siglas en inglés) es el marco principal para describir la evolución del campo magnético fotosférico a gran escala del Sol y predecir las amplitudes de los futuros ciclos solares. Aunque los procesos de transporte lineal (flujo meridional, rotación diferencial y difusión turbulenta) están relativamente bien constreñidos, el espacio de parámetros del modelo sigue siendo amplio e incierto. Estudios de optimización previos (p. ej., Petrovay & Talafha 2019) identificaron "islas" de parámetros admisibles, pero a menudo dependieron de supuestos lineales o descuidaron retroalimentaciones no lineales específicas. Además, el papel de un término de decaimiento de flujo fenomenológico en las formulaciones parametrizadas ha sido objeto de debate: algunas simulaciones de asimilación de datos sugieren que podría ser innecesario, mientras que otras indican que su ausencia conduce a campos dipolares irrealmente persistentes. El impacto específico de los mecanismos de amortiguación (quenching) no lineales —específicamente la amortiguación de la inclinación (TQ) y la amortiguación de la latitud (LQ)— sobre el espacio de parámetros admisible de la velocidad del flujo meridional ($u_0$), la difusividad superficial ($\eta$) y la escala de tiempo de decaimiento ($\tau$) no había sido cuantificado sistemáticamente dentro de un marco de optimización unificado.

Metodología
Los autores extendieron el marco de optimización del espacio de parámetros establecido por Petrovay & Talafha (2019) introduciendo prescripciones analíticas para retroalimentaciones no lineales en el término fuente del SFT.

• Configuración del modelo: La ecuación SFT gobierna la evolución del campo magnético radial $B(\lambda, t)$, incorporando la advección por flujo meridional, la difusión y un término de decaimiento ($-B/\tau$). El término fuente $S(\lambda, t)$ representa la emergencia de regiones magnéticas bipolares (BMR) basadas en un ciclo solar promediado estadísticamente.
• Retroalimentaciones no lineales: Se implementaron dos mecanismos de amortiguación:
  1. Amortiguación de la inclinación (TQ): La separación latitudinal de las BMR (ley de Joy) se reduce para ciclos más fuertes, gobernada por un parámetro de no linealidad $b_{joy} = 0.15$.
  2. Amortiguación de la latitud (LQ): La latitud media de emergencia aumenta con la amplitud del ciclo, gobernada por $b_{lat} = 2.4$.
• Estrategia de optimización: El modelo se resolvió numéricamente sobre una rejilla densa de parámetros ($u_0$, $\eta$, $\tau$) durante 100 ciclos solares para alcanzar un estado cíclico estacionario. Las soluciones admisibles se definieron por su acuerdo con tres restricciones observacionales:
  1. Tiempo de inversión del campo polar ($T_{rev}$).
  2. Relación de amplitud entre el mínimo y el extremo del campo polar ($R_{pol}$).
  3. Latitud del límite del casquete polar ($\lambda_{cap}$).
     Se utilizaron restricciones adicionales basadas en el dipolo ($D_{rev}$, $D_{min/ext}$) con fines diagnósticos.
• Escenarios: El estudio comparó casos lineales (sin amortiguación) frente a casos con solo TQ, solo LQ, y la combinación TQ+LQ, a través de diferentes escalas de tiempo de decaimiento ($\tau = 8$ años y $\tau = \infty$).

Resultos clave

• Impacto de la amortiguación: Tanto la TQ como la LQ reducen el dominio de parámetros admisibles, pero la LQ ejerce una influencia significativamente más fuerte. La acción combinada de TQ y LQ produce un efecto de "techo" pronunciado, saturando la amplificación del dipolo axial y evitando campos polares irrealmente grandes en ciclos fuertes.
• Contracción del espacio de parámetros:
  • En el caso lineal (sin amortiguación) con $\tau = 8$ años, las regiones admisibles se concentraron alrededor de $\eta \approx 600\text{--}650$ km$^2$ s$^{-1}$ y $u_0 \approx 18\text{--}19$ m s$^{-1}$, con una región secundaria en valores más bajos.
  • La introducción de solo TQ causó una contracción modesta.
  • La introducción de solo LQ estrechó los rangos significativamente.
  • El caso combinado TQ+LQ con $\tau = 8$ años restringió el dominio admisible a bandas estrechas, específicamente $u_0 \approx 18\text{--}19$ m s$^{-1}$ y $\eta \approx 600\text{--}650$ km$^2$ s$^{-1}$, eliminando efectivamente las islas de soluciones secundarias encontradas en el caso lineal.
• Papel de la escala de tiempo de decaimiento ($\tau$): La inclusión de una escala de tiempo de decaimiento finita ($\tau \approx 8\text{--}10$ años) es crítica. Cuando $\tau = \infty$ (sin decaimiento), el espacio de parámetros admisible se expande considerablemente, incluso con la amortiguación activa. Sin embargo, los autores señalan que en modelos parametrizados sin decaimiento explícito, esto puede conducir a campos dipolares irrealmente persistentes y reversiones retrasadas. El rango óptimo de $\tau \approx 8\text{--}10$ años refleja la escala de tiempo necesaria para reproducir el tiempo de reversión observado en presencia de retroalimentaciones no lineales.
• Equilibrio acoplado: Las correlaciones entre $u_0$, $\eta$ y $\tau$ revelan un equilibrio estrechamente acoplado entre la advección, la difusión y la pérdida de flujo fenomenológica. A medida que $\tau$ disminuye (decaimiento más fuerte), los rangos admisibles para $u_0$ y $\eta$ se desplazan y se estrechan, requiriendo flujos meridionales más lentos o difusividades específicas para mantener la consistencia observacional.

Significado y afirmaciones
El artículo afirma que la incorporación de mecanismos de amortiguación no lineales transforma el paisaje de optimización del SFT de amplias "islas" de soluciones desconectadas en "crestas" estrechas y físicamente consistentes.

• Régimen autolimitado: El estudio concluye que la amortiguación no lineal, combinada con una escala de tiempo de pérdida de flujo fenomenológica, define un régimen de SFT autolimitado. Este régimen restringe la saturación del dínamo al vincular los parámetros de transporte superficial directamente con la retroalimentación no lineal en el término fuente.
• Límites de predictibilidad: Los resultados sugieren que el dínamo solar opera cerca de un umbral de saturación marginalmente estable. Si bien la amplitud media del ciclo está constreñida por la amortiguación, la sensibilidad inherente a los parámetros de transporte superficial implica que los modelos puramente cinemáticos tienen horizontes de predicción limitados debido a la naturaleza estocástica de la emergencia de regiones activas.
• Dependencia del modelo: Los autores declaran explícitamente que sus hallazgos son condicionales a la formulación de la fuente parametrizada adoptada. Reconocen que las simulaciones de asimilación de datos (p. ej., Yeates et al. 2025), que ingieren directamente las regiones activas observadas, pueden no requerir un término de decaimiento explícito, lo que sugiere que la necesidad de $\tau$ depende de la implementación numérica de la fuente.
• Interpretación física: El efecto de "techo" observado se interpreta como un mecanismo físico donde la entrada de flujo adicional ya no se traduce en campos polares más fuertes una vez que la amortiguación limita la eficiencia de la construcción del dipolo. Esto proporciona una imagen unificada que vincula las corrientes de entrada superficiales, la amortiguación y los procesos del dínamo más profundos.

En resumen, este trabajo demuestra que las retroalimentaciones no lineales no son meras correcciones menores, sino fundamentales para estrechar el espacio de parámetros viable del modelo SFT, imponiendo un equilibrio autosistente entre el transporte y la pérdida de flujo que regula la amplitud del ciclo solar.