Resolving the Tachocline using Inversion of Rotational Splitting Derived from Fitting Very Long and Long Time Series

Resumen Técnico: Resolución de la Taquoclina mediante la Inversión del Desdoblamiento Rotacional Derivado del Ajuste de Series Temporales Muy Largas y Largas

Planteamiento del Problema
La taquoclina solar, la capa de cizalladura que separa el interior radiativo de rotación casi rígida de la zona de convección de rotación diferencial, sigue siendo una región de incertidumbre significativa respecto a su geometría precisa, variación latitudinal y estabilidad temporal. Aunque los análisis heliosísmicos de las últimas tres décadas han localizado la taquoclina cerca de $0.7 R_\odot$, persisten discrepancias sustanciales entre estudios sobre su espesor, perfil latitudinal (prolato frente a esférico) y evolución temporal. Estas inconsistencias surgen de las limitaciones en la resolución espacial, las diferencias en la calidad y longitud de los conjuntos de datos, y las parameterizaciones y metodologías de inversión específicas empleadas. Los autores postulan que resolver estas ambigüedades requiere técnicas de inversión capaces de capturar gradientes agudos y características de escala fina, minimizando al mismo tiempo el suavizado que típicamente se requiere para resolver la naturaleza intrínsecamente mal planteada del problema inverso de la rotación.

Metodología
El estudio utiliza desdoblamientos rotacionales derivados de series temporales de observaciones muy largas (25.2 años) y largas (12.6 y 6.3 años) de los instrumentos GONG, MDI y HMI, computados mediante la metodología independiente de Korzennik (2023). El enfoque técnico central involucra dos técnicas de inversión distintas aplicadas a estos desdoblamientos:

1. Mínimos Cuadrados Regularizados (RLS): Un método rápido y no iterativo que utiliza un único parámetro de compromiso para controlar el suavizado mediante la minimización de la segunda derivada de la solución.
2. Técnica de Reconstrucción Algebraica Simultánea (SART): Un método iterativo que permite un suavizado global o local. Un aspecto novedoso de esta implementación es el uso de un "ajuste inicial informado", donde un ajuste sigmoideo agudizado de la solución RLS sirve como punto de partida para la iteración SART.

Una innovación metodológica crítica es la implementación de un mallado radial variable. Para resolver la taquoclina sin incurrir en costos computacionales prohibitivos, los autores probaron cuatro configuraciones de malla:

• Uniforme: Alta resolución en todas partes, lo que conduce a tamaños de matriz grandes.
• Progresiva: El espaciamiento radial varía suavemente (polinómico en $r$), ofreciendo resoluciones de 0.25% y 0.18% en la taquoclina.
• Agresiva: Una malla con resolución extremadamente alta (0.09%) concentrada únicamente alrededor de la taquoclina, volviéndose significativamente más gruesa en las zonas radiativa y de convección.

Los autores también introdujeron el concepto de núcleos de promedio extendidos para calcular el "centro de gravedad" de los núcleos de inversión a través de todo el rango latitudinal ($-\pi/2$ a $\pi/2$). Esto se hizo para evaluar si la ubicación efectiva de la tasa de rotación inferida se desplaza de la ubicación objetivo debido a la asimetría del núcleo, particularmente cerca del ecuador.

Resultados Clave

• Rendimiento de la Malla: Las simulaciones revelaron que, si bien la malla "agresiva" ofrece alta resolución, introduce artefactos numéricos (fenómenos de tipo Gibbs) y magnificación de ruido, particularmente cuando hay ruido presente. Se encontró que las mallas "progresivas" mitigaban estos artefactos manteniendo una alta resolución, aunque a un costo computacional mayor.
• Ubicación y Geometría de la Taquoclina: Las inversiones de los 25.2 años de datos de GONG revelan una discontinuidad distintiva en la ubicación de la taquoclina. En latitudes bajas ($\phi \le 30^\circ$), la taquoclina se localiza más profundamente (por debajo de la base de la zona de convección por $\sim 0.02 R_\odot$), mientras que en latitudes altas, coincide con la base de la zona de convección. Esta bifurcación se observa tanto en las inversiones RLS como en las SART y persiste a través de diferentes resoluciones de malla.
• Ancho y Cizalladura: El ancho estimado de la taquoclina (Ancho Completo a la Mitad del Salto, FWHJ) es altamente sensible a la resolución de la malla y al suavizado. Las inversiones de alta resolución sugieren una taquoclina muy estrecha (FWHJ $< 1\% R_\odot$). La amplitud de la cizalladura radial muestra un cambio de signo entre latitudes bajas y altas, consistente con la discontinuidad de ubicación.
• Variabilidad Temporal: Los análisis de series temporales más cortas (6.3 y 12.6 años) no muestran variaciones temporales definitivas o significativas en la posición, el ancho o la amplitud de la cizalladura de la taquoclina. La dispersión observada está dominada por diferencias sistemáticas entre las metodologías de inversión (RLS vs. SART) y los distintos grados de suavizado requeridos, más que por la variabilidad solar intrínseca.
• Validación: La variación latitudinal inferida de la posición de la taquoclina y la amplitud de la cizalladura se alinea bien con las estimaciones derivadas del modelado directo (Basu & Antia 2003; Basu et al. 2024), aunque el modelado directo utilizando expansiones polinómicas no logra reproducir la discontinuidad observada.

Significancia y Reivindicaciones
El artículo afirma proporcionar una caracterización cuantitativamente restringida de la taquoclina que va más allá de las descripciones cualitativas generales. La principal significancia radica en la identificación de una estructura bifurcada para la taquoclina, donde su profundidad y la dirección de la cizalladura cambian abruptamente entre latitudes bajas y altas. Este hallazgo desafía los modelos previos que asumían una transición latitudinal suave o una geometría puramente prolata.

Los autores señalan modestamente que, si bien sus resultados concuerdan con las estimaciones del modelado directo, la discontinuidad observada no es reproducida por los modelos teóricos actuales o por el modelado directo polinómico. Reconocen que las diferencias sistemáticas entre los métodos de inversión y el ruido en las series temporales más cortas impiden actualmente la detección de variaciones temporales definitivas vinculadas al ciclo solar. En consecuencia, el artículo concluye que, aunque la bifurcación es una característica robusta de sus inversiones, representa un nuevo paradigma para la forma de la taquoclina que requiere confirmación independiente y desarrollo teórico adicional para ser explicado. El estudio enfatiza que resolver tales características de escala fina requiere un equilibrio cuidadoso entre la resolución de la malla, el suavizado y la gestión del ruido, siendo las mallas progresivas un camino viable hacia adelante.