Dynamo-action confinement of the solar tachocline: viability in the solar parameter regime and observational tests

The solar tachocline, the thin (~0.05 solar-radius) shear layer at the base of the convection zone, has been an outstanding problem since Spiegel & Zahn (1992) showed that radiative spreading should have thickened it far beyond what helioseismology observes. Matilsky et al. (2022) presented the first 3-D global magnetohydrodynamic simulation in which a self-excited, nonaxisymmetric, quasi-cyclic dynamo confines a statistically steady tachocline, the field penetrating the radiative interior by a magnetic skin effect and local amplification by vigorous horizontal motions — but at three times the solar rotation rate and with diffusivities enhanced by many orders of magnitude. Follow-up work (2024–2026) extended the scenario to the physically relevant radiative-spreading regime and found that confinement improves as stratification trends toward solar values, while also showing that the mechanism appears to require turbulently enhanced magnetic diffusivity (≳3×10^10 cm^2/s) in strong tension with the solar lithium abundance, and that all simulations to date remain 2–3 or more orders of magnitude from solar parameters. The open problem is whether dynamo-action confinement of the tachocline operates in the Sun's actual regime — solar rotation, stratification, and diffusivities consistent with chemical-mixing constraints; the real Sun's multi-timescale dynamo spectrum; dynamically stable field topologies — and what helioseismic or surface-activity signatures could confirm or falsify it.

Background

Helioseismology — the inversion of the Sun's global oscillation frequencies — reveals a thin shear layer of width about $0.05\,R_\odot$ (where $R_\odot$ is the solar radius) near $r\simeq0.7\,R_\odot$, called the tachocline, which separates the differentially rotating convection zone (CZ, the outer, turbulently convecting ~30% by radius) from the radiative zone (RZ, the stably stratified interior rotating nearly as a solid body). Spiegel & Zahn (1992) argued that without an opposing torque, meridional circulation driven by radiative diffusion of the CZ's latitudinal temperature gradient would burrow downward ("radiative spreading") and imprint differential rotation through the RZ over the Sun's lifetime; the observed thinness of the tachocline therefore demands a torque that enforces uniform rotation below it. Three families of confinement scenario dominate the literature. In the fast hydrodynamic scenario, anisotropic turbulence (mostly horizontal, owing to the strong stable stratification) provides the torque on month-to-year timescales. In the slow magnetohydrodynamic scenario, a fossil (primordial) magnetic field threading the RZ does so: Gough & McIntyre (1998) argued that a large-scale primordial field is not merely helpful but "inevitable" if the RZ is to remain uniformly rotating. In the fast magnetic scenario, the cyclic dynamo field of the CZ (22-year magnetic cycle) diffuses downward into a magnetic skin depth — analogous to the skin effect for alternating currents in a conductor — as explored by Forgács-Dajka & Petrovay (2001); Spruit (2002) had additionally proposed that differential rotation in a stably stratified interior can itself sustain a small-scale dynamo whose field feeds back on the shear. A key theoretical constraint is Ferraro's law of isorotation: in a steady axisymmetric state, angular velocity must be constant along poloidal field lines, so any poloidal field penetrating the tachocline exerts Maxwell stresses (magnetic tension and pressure forces) that resist differential rotation.

Matilsky, Hindman, Featherstone, Blume & Toomre (2022) presented the first 3-D, global, anelastic (sound-filtered) magnetohydrodynamic simulation, evolved with the Rayleigh code, that couples a CZ to a stable RZ and produces a statistically stationary tachocline confined by a self-excited dynamo. The confining field is nonaxisymmetric and quasi-cyclic; it enters the RZ partly by diffusive imprinting from above and partly by local inductive amplification driven by vigorous horizontal motions, which the authors identify primarily with equatorial Rossby waves (global-scale planetary-type inertial waves restored by the Coriolis force, with dispersion relation $\omega_{\ell m}=-2\Omega_0 m/[\ell(\ell+1)]$ for azimuthal wavenumber $m$ and spherical-harmonic degree $\ell$). However, the simulation ran at three times the solar Carrington rotation rate ($\Omega_0=3\,\Omega_\odot$, chosen to avoid numerically accessible but unphysical "anti-solar" differential rotation), with magnetic Prandtl number $\mathrm{Pr_m}=4$ (ratio of viscous to magnetic diffusivity) and diffusivities enormously enhanced above solar values; the confined tachocline spread viscously rather than by the physically relevant radiative spreading. The paper itself flags what remains unknown: the definitive origin and parameter dependence of the deep horizontal motions; whether the local amplification operates at solar (microphysical or weakly turbulent) diffusivities; whether the mechanism and its timescales survive solar rotation and diffusivities; whether the real Sun's multi-timescale dynamo spectrum — the primary 22-year cycle plus grand minima and quasi-biennial modulations — delivers comparable cumulative confining torques; and the dynamical stability, long-term evolution, and observational signatures of a radiative-interior dynamo.

Problem Statement

Establish whether the dynamo-confinement mechanism for the solar tachocline — in which a self-excited, nonaxisymmetric, quasi-cyclic dynamo, generated in and above the convective overshoot layer, penetrates the radiative interior through a magnetic skin effect and/or local inductive amplification by horizontal motions, and exerts Maxwell stresses that enforce nearly uniform rotation — operates under the Sun's actual conditions. This one objective comprises four jointly supporting questions, all posed by the source series and left open by it:

(i) Solar parameter regime: does confinement persist at the solar rotation rate $\Omega_0=\Omega_\odot$ and at stable stratification and diffusivities approaching solar values — i.e., molecular diffusivities, or turbulent enhancement no larger than allowed by chemical-mixing constraints such as the lithium-depletion bound $D_T\lesssim500\,\mathrm{cm^2\,s^{-1}}$? In particular, resolve or dissolve the tension that skin-depth confinement appears to require a turbulently enhanced magnetic diffusivity $\eta_T\gtrsim3\times10^{10}\,\mathrm{cm^2\,s^{-1}}$, roughly eight orders of magnitude above the lithium bound.

(ii) Realistic dynamo spectrum: does the Sun's actual dynamo — the primary 22-year magnetic cycle together with its modulations (grand minima, quasi-biennial oscillations) and its nonaxisymmetric components — deliver cumulative confining torques at tachocline depths comparable to those of the idealized simulations?

(iii) Stability and long-term evolution: are the nonaxisymmetric field topologies required for confinement dynamically stable against the shear, baroclinic, and magnetic (e.g., Tayler-type) instabilities expected in the tachocline, and can the scenario simultaneously account for the observed long-term (gigayear) spin-down of the radiative zone?

(iv) Observational signatures: what helioseismic signatures (e.g., quasi-cyclic variations of tachocline depth or of rotation beneath it) or surface-activity signatures (e.g., long-term active-longitude variability) would distinguish dynamo confinement from the fossil-field and hydrodynamic alternatives, and are they present in existing helioseismic and solar-activity data?

An answer may take the form of solar-regime self-consistent simulations, theoretical scaling or incompatibility arguments (including the diffusivity constraints), observational analyses, or a combination; it must decide, or decisively constrain, whether the tachocline is in fact confined by dynamo action in the radiative interior rather than by a fossil field or hydrodynamic turbulence. The four sub-questions (i)–(iv) are jointly supporting rather than separately mandatory; the task is resolved by a decision on, or decisive constraint upon, this central question, drawing on whichever of (i)–(iv) are relevant.

References

1. Loren I. Matilsky, Bradley W. Hindman, Nicholas A. Featherstone, Catherine C. Blume, Juri Toomre (2022), "Confinement of the Solar Tachocline by Dynamo Action in the Radiative Interior", The Astrophysical Journal Letters, DOI: 10.3847/2041-8213/ac93ef, arXiv:2206.12920, https://doi.org/10.3847/2041-8213/ac93ef
2. L. Matilsky, N. Brummell, B. W. Hindman, J. Toomre (2024), "Solar Tachocline Confinement by the Nonaxisymmetric Modes of a Dynamo Magnetic Field", The Astrophysical Journal, DOI: 10.3847/1538-4357/ad18b2, https://doi.org/10.3847/1538-4357/ad18b2
3. Loren I. Matilsky, Lydia Korre, Nicholas H. Brummell (2025), "Dynamo Confinement of a Radiatively Spreading Solar Tachocline Revealed by Self-consistent Global Simulations", The Astrophysical Journal Letters, DOI: 10.3847/2041-8213/adefe3, arXiv:2507.08256, https://doi.org/10.3847/2041-8213/adefe3
4. Loren I. Matilsky, Lydia Korre, Nicholas H. Brummell (2026), "A Dynamo Confinement Scenario for the Solar Tachocline and Its Implications for Spin-down in the Radiative Spreading Regime", The Astrophysical Journal, DOI: 10.3847/1538-4357/ae6705, arXiv:2601.11943, https://doi.org/10.3847/1538-4357/ae6705
5. P. Garaud, D. Gough, L. Matilsky (2025), "Toward a Self-consistent Hydrodynamical Model of the Solar Tachocline", The Astrophysical Journal, DOI: 10.3847/1538-4357/adc72b, https://doi.org/10.3847/1538-4357/adc72b
6. C. Blume, B. W. Hindman, L. I. Matilsky (2024), "Inertial Waves in a Nonlinear Simulation of the Sun's Convection Zone and Radiative Interior", The Astrophysical Journal, DOI: 10.3847/1538-4357/ad27d1, arXiv:2312.14270, https://doi.org/10.3847/1538-4357/ad27d1
7. A. Strugarek, B. Belucz, A. Brun, M. Dikpati, G. Guerrero (2023), "Dynamics of the Tachocline", Space Science Reviews, DOI: 10.1007/s11214-023-01027-0, https://doi.org/10.1007/s11214-023-01027-0
8. D. O. Gough, M. E. McIntyre (1998), "Inevitability of a magnetic field in the Sun's radiative interior", Nature 394, 755, DOI: 10.1038/29472, https://doi.org/10.1038/29472
9. E. A. Spiegel, J.-P. Zahn (1992), "The solar tachocline", Astronomy & Astrophysics 265, 106, bibcode 1992A&A...265..106S, https://articles.adsabs.harvard.edu/pdf/1992A%26A...265..106S
10. H. C. Spruit (2002), "Dynamo action by differential rotation in a stably stratified stellar interior", Astronomy & Astrophysics 381, 923, DOI: 10.1051/0004-6361:20011465, https://doi.org/10.1051/0004-6361:20011465
11. E. Forgács-Dajka, K. Petrovay (2001), "Tachocline confinement by an oscillatory magnetic field", Solar Physics 203, 195, DOI: 10.1023/A:1013389631585, https://doi.org/10.1023/A:1013389631585