pith. sign in

arxiv: 2601.11943 · v3 · submitted 2026-01-17 · 🌌 astro-ph.SR

A Dynamo Confinement Scenario for the Solar Tachocline and its Implications for Spin-down in the Radiative Spreading Regime

Pith reviewed 2026-05-16 13:22 UTC · model grok-4.3

classification 🌌 astro-ph.SR
keywords solar tachoclinedynamo confinementradiative spreadingMaxwell stressesmagnetic skin effectsolar spin-downdifferential rotation
0
0 comments X

The pith

A dynamo-generated magnetic field confines the solar tachocline against radiative spreading via Maxwell stresses and a skin effect.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The solar tachocline should spread inward over time due to meridional flows burrowing into the radiative zone, yet observations show it remains thin. Global simulations demonstrate that a dynamo operating near the base of the convection zone generates nonaxisymmetric magnetic modes that penetrate downward through a magnetic skin effect. These modes produce Maxwell stresses that hold the shear layer in place, creating a statistically stationary tachocline. As the stable stratification below is made stronger, the dynamo cycles lengthen, the magnetic field reaches deeper, and confinement improves. The same stresses also transmit spin-down from the tachocline region into the deeper radiative interior.

Core claim

In global simulations trending toward solar parameters, the large-scale nonaxisymmetric modes of a dynamo penetrate the tachocline and below through a magnetic skin effect. The resulting Maxwell stresses confine the tachocline against radiative spreading, producing a statistically stationary narrow shear layer. Stronger stratification lengthens dynamo cycle periods, deepens magnetic penetration, and tightens confinement. The burrowing circulation spins down the tachocline, and the stresses transmit this slowing to rigidify the deeper radiative zone.

What carries the argument

The magnetic skin effect of nonaxisymmetric dynamo modes that allows Maxwell stresses to confine the tachocline against radiative spreading.

If this is right

  • The tachocline becomes progressively more confined as stable stratification is strengthened.
  • Maxwell stresses transmit spin-down from the tachocline into the deeper radiative zone, promoting rigid rotation.
  • The dynamo offers a direct pathway linking near-surface spin-down to the deep solar interior.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the skin-effect mechanism operates at solar parameters, it could maintain the observed thin tachocline without requiring additional transport processes.
  • Similar dynamo confinement may operate in other stars that possess tachoclines or strong radial shear layers.
  • Higher-resolution simulations at solar diffusivities could test whether the cycle lengthening and deeper penetration saturate or continue.

Load-bearing premise

That the observed trends with increasing stratification will continue all the way to solar values of magnetic diffusivity, viscosity, and thermal diffusivity.

What would settle it

A simulation performed at solar values of the diffusivities in which the tachocline continues to spread or magnetic penetration fails to deepen would falsify the confinement scenario.

Figures

Figures reproduced from arXiv: 2601.11943 by Loren I. Matilsky, Lydia Korre, Nicholas H. Brummell.

Figure 1
Figure 1. Figure 1: (a) Schematic of our three-dimensional spherical-shell simulation domain. The yellow region depicts the CZ with bulk overturning fluid motions represented by circular arrows. The CZ lies atop a thin stably stratified region (green) with significant nonlinear effects from overshoot and a deep RZ (purple) where the motion is essentially linear, represented by curvy “radiation” arrows. (b) Schematic of the re… view at source ↗
Figure 2
Figure 2. Figure 2: Location in parameter space (Bu and σ) of the different simulations and the Sun. Cases 1 and 2 have the same Bu and σ, but different Raf (see [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Spherical and meridional surfaces cut out of the simulation domain showing snapshots of the simulated flows in Case M6 for (a) the fluctuating axial vorticity ω ′ z and (b) the colatitudinal component of the poloidal magnetic field, Bθ. Red/yellow/blue tones indicate positive/zero/negative values. Each snapshot is taken toward the end of the simu￾lation (t = trun). The inner and outer spherical surfaces bo… view at source ↗
Figure 4
Figure 4. Figure 4: Steady-state differential rotation ⟨2Ω − 1⟩ t for the full simulation suite, plotted in the meridional plane. Red/yellow/blue tones indicate positive/zero/negative values. The outer ticks in each color bar mark the the saturation values (which are labeled) and the inner tick marks zero. Contours are equally spaced (for positive and negative values separately) and the zero contour is dashed. The dashed inte… view at source ↗
Figure 5
Figure 5. Figure 5: (a) Differential rotation contrast in the RZ defined in Equation (32a), plotted with respect to case number. (b) Differential rotation contrast in the CZ defined in Equation (32b). (c) The confinement ratio f defined in Equation (34). Note that because of the chosen ordering, f monotonically decreases for the MHD cases within the radiative spreading regime. (d) The RZ’s shear-layer thickness defined in Equ… view at source ↗
Figure 6
Figure 6. Figure 6: Time-averaged meridional circulation profiles |⟨ρ˜um⟩ t |sgn(⟨Ψ⟩ t ) (see Equations 36 and 37) for the full simulation suite, plotted in the meridional plane. Red tones indicate clockwise circulation and blue tones indicate counterclockwise circulation. The solid contours mark ⟨Ψ⟩ t ≡ 0, which are the circulation cell boundaries and also streamlines of um. The dashed semicircles mark r = rbcz. τmc := − ⟨ρ˜… view at source ↗
Figure 7
Figure 7. Figure 7: Cumulatively-time-integrated torque densities from Equation (40) as functions of time, after volume-averaging over the deep RZ, for (a) Case H0, (b) Case M0, (c) Case H6, and (d) Case M6. In the legend, we label the curves by the torque densities, which are understood to be shorthand for the time-integrated and volume-averaged quantities. For instance, R ⟨τrs⟩ is shorthand for R t 0 ⟨τrs⟩drz (t ′ )dt′ . Th… view at source ↗
Figure 8
Figure 8. Figure 8: (a) σdyn(t) for Cases H3, H4, and H6 and (b) Cases M3, M4, and M6. (c) Angular momentum deficit of the RZ, ∆L(t), for Cases H0, H3, H4, and H6 and (d) Cases M0, M3, M4, and M6. (e) and (f): Shear-layer thickness ∆(t) for the same cases as panels (c) and (d), respectively. See Equations (24), (41), and (43). The vertical colored lines correspond to the values of t3/4 for each case. of the CZ is determined b… view at source ↗
Figure 9
Figure 9. Figure 9: Equilibrated-state values of σdyn, de￾noted by (σdyn(∞), for all of our cases, plotted with respect to case number. We approximate σdyn(∞) ≈ [⟨τv⟩drz,t / ⟨τmc⟩drz,t] 1/2 , where the tempo￾ral average is over the equilibrated state. where m ranges across its discrete values −ℓmax + 1 to ℓmax (see Appendix B). Note that we have chosen the decomposition in Equation (45) such that Parseval’s the￾orem takes the… view at source ↗
Figure 10
Figure 10. Figure 10: Mollweide projections and power spectra in azimuthal wavenumber for the poloidal magnetic field Bpol in the upper RZ (r ∗ /R⊙ = 0.692) for (a–c) Case M3, (d–f) Case M4, and (g–i) Case M6. The power spectra |Bm pol| 2 (Equation 45) have been horizontally and temporally averaged over the equilibrated state. problem, we follow P. Garaud (1999) and define rη := rin + R r rin η˜(r ′ ) −1/2dr′ R rskin rin η˜(r … view at source ↗
Figure 11
Figure 11. Figure 11: Time-latitude diagrams over the interval t = (27500, 52500) of Re(B 1 θ ) in the upper RZ (r ∗ /R⊙ = 0.692) and the corresponding truncated horizontally averaged power spectra in frequency |B1ω θ | 2 for (a–b) Case M3, (c–d) Case M4, and (e–f) Case M6. The full time series t = (0, trun) was used to obtain the frequency spectra and a Hanning window function was employed to reduce high-frequency noise. For … view at source ↗
Figure 12
Figure 12. Figure 12: Scatter symbols: Realized amplitude of [PITH_FULL_IMAGE:figures/full_fig_p022_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Difference of the bulk rotation rate of the RZ from the background frame rate, Ωrz − 1/2, as a function of Bu for the MHD cases (see Equation 33). The dotted line separates the radiative spreading cases on the left from the lone viscous spreading case on the right. magnetic diffusivity in the presence of strong Maxwell stresses. Assuming that skin depths can be significantly en￾hanced, whether through tur… view at source ↗
Figure 14
Figure 14. Figure 14: (a) Horizontally averaged differential rotation, [PITH_FULL_IMAGE:figures/full_fig_p023_14.png] view at source ↗
read the original abstract

At the base of the Sun's convective zone, a narrow shear layer called the tachocline separates strong latitudinal differential rotation above from nearly rigid rotation in the radiative zone below. The observed thinness of the tachocline is a long-standing dynamical puzzle because the tachocline should have spread significantly due to inward-burrowing meridional circulation, also called "radiative spreading." We recently presented the first pair of global simulations to reveal a statistically stationary tachocline confined against radiative spreading by the Maxwell stresses from the large-scale nonaxisymmetric modes of a dynamo, which penetrated into and below the tachocline through a novel magnetic skin effect. In the work presented here, we systematically examine how this "dynamo confinement scenario" works against radiative spreading in a suite of simulations as the governing parameters trend in the direction of the true solar regime. We find that as the stable stratification of the radiative zone is made progressively stronger, the dynamo cycles get longer, the magnetic field consequently penetrates deeper due to the skin effect, and the tachocline becomes more confined. Furthermore, these results have interesting consequences for solar spin-down. In all of our radiatively spreading simulations, the tachocline region spins down due to the burrowing circulation. Below the tachocline, the Maxwell stresses transmit this spin-down further to rigidify the deeper radiative zone. We thus speculate that, in addition to confining the tachocline, the dynamo may provide a pathway to communicate spin-down from the near-surface layers to the deep interior.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript presents a suite of global MHD simulations showing that a self-sustained dynamo can confine the solar tachocline against radiative spreading. Maxwell stresses from large-scale nonaxisymmetric modes penetrate the tachocline via a magnetic skin effect and balance the inward burrowing of meridional circulation. As stable stratification is increased, dynamo cycles lengthen, skin-effect penetration deepens, and the tachocline becomes more confined. The same stresses are reported to transmit spin-down from the tachocline region into the deeper radiative interior.

Significance. If the reported trends persist at solar values of magnetic diffusivity, the work supplies a concrete, simulation-based mechanism for the observed thinness of the tachocline and a pathway for angular-momentum transport across the radiative zone. The global, self-consistent treatment of both dynamo action and radiative spreading constitutes a clear advance over earlier analytic or local models.

major comments (2)
  1. [Abstract and simulation-suite description] The parameter study strengthens only the stable stratification (N²) while holding magnetic diffusivity η, viscosity ν, and thermal diffusivity κ fixed at numerically accessible values orders of magnitude larger than solar. Because the skin depth scales as √(η/ω) and the central confinement claim depends on this depth at solar η, the observed improvement in confinement with stratification cannot be extrapolated without either reduced-η runs or an explicit scaling analysis (see abstract and the description of the simulation suite).
  2. [Results on spin-down and confinement] No torque-balance diagnostics, integrated Maxwell-stress profiles, or direct comparison of Maxwell versus Reynolds stresses against the radiative-spreading term are shown. Without these quantities it is impossible to verify that Maxwell stresses dominate the force balance in the reported regime or in the extrapolated solar regime (see results on spin-down transmission).
minor comments (2)
  1. [Abstract] The abstract states that the tachocline is 'statistically stationary' but provides no information on averaging intervals, convergence criteria, or resolution studies; these details should be added to allow assessment of numerical robustness.
  2. [Figures and captions] Figure captions and axis labels should explicitly state the range of N² values explored and the fixed values of η, ν, and κ so that readers can immediately judge how far the suite lies from solar parameters.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed and constructive report. We address each major comment below and will revise the manuscript to strengthen the presentation of our results.

read point-by-point responses
  1. Referee: [Abstract and simulation-suite description] The parameter study strengthens only the stable stratification (N²) while holding magnetic diffusivity η, viscosity ν, and thermal diffusivity κ fixed at numerically accessible values orders of magnitude larger than solar. Because the skin depth scales as √(η/ω) and the central confinement claim depends on this depth at solar η, the observed improvement in confinement with stratification cannot be extrapolated without either reduced-η runs or an explicit scaling analysis (see abstract and the description of the simulation suite).

    Authors: We agree that η, ν, and κ are held fixed while only N² is varied, and that a direct extrapolation to solar η requires care. However, the manuscript already shows that stronger N² produces longer dynamo cycles (smaller ω), which increases the skin depth √(η/ω) at fixed η. In the revision we will add an explicit scaling analysis in a new subsection that quantifies how the observed trend in cycle period with N² combines with the skin-depth formula to support deeper penetration and improved confinement even at solar η. We do not have the computational resources for reduced-η runs in the current study, but the scaling analysis directly addresses the extrapolation concern. revision: yes

  2. Referee: [Results on spin-down and confinement] No torque-balance diagnostics, integrated Maxwell-stress profiles, or direct comparison of Maxwell versus Reynolds stresses against the radiative-spreading term are shown. Without these quantities it is impossible to verify that Maxwell stresses dominate the force balance in the reported regime or in the extrapolated solar regime (see results on spin-down transmission).

    Authors: We acknowledge that the current manuscript lacks explicit torque-balance diagnostics and integrated stress profiles. In the revised version we will add a new figure (and accompanying text) showing radial profiles of the Maxwell and Reynolds stresses, their integrated torques, and a direct comparison against the radiative-spreading term throughout the tachocline and deeper radiative zone. These diagnostics will confirm that Maxwell stresses dominate the force balance and transmit the spin-down, both in the simulated regime and in the extrapolated solar regime. revision: yes

Circularity Check

0 steps flagged

Numerical simulations with independent parameter trends; no derivation reduces to inputs

full rationale

The manuscript reports results from a suite of global MHD simulations that vary stable stratification while holding diffusivities fixed. All load-bearing claims (tachocline confinement via Maxwell stresses, deeper penetration with longer cycles, spin-down transmission) are direct outputs of the time-dependent runs rather than any analytical step that equates a prediction to a fitted input or self-citation by construction. Prior work is cited only for context; no uniqueness theorem, ansatz, or renaming of known results is invoked to close the argument. The acknowledged gap in solar diffusivity values is an extrapolation issue, not a circularity in the reported trends.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

Central claim rests on numerical integration of standard MHD equations in a spherical shell with imposed stable stratification; no new physical entities are introduced.

free parameters (1)
  • stable stratification strength
    Progressively increased to trend toward solar regime; controls cycle length and penetration depth.
axioms (1)
  • standard math Magnetohydrodynamic equations govern the plasma dynamics
    Standard governing equations for global dynamo simulations in stellar interiors.

pith-pipeline@v0.9.0 · 5594 in / 1290 out tokens · 45100 ms · 2026-05-16T13:22:57.192732+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Superrotation and Jet Migration in Simulations of Jupiter's Convective Zone and Weather Layer

    astro-ph.EP 2026-05 unverdicted novelty 5.0

    3D anelastic simulations find that both convective zone Busse columns and weather-layer PV homogenization produce Jupiter-like jets, with the weather layer causing deviations in thermal wind balance and long-term high...

Reference graph

Works this paper leans on

63 extracted references · 63 canonical work pages · cited by 1 Pith paper

  1. [1]

    A., Garaud, P., & Wood, T

    Acevedo-Arreguin, L. A., Garaud, P., & Wood, T. S. 2013, Mon. Not. R. Astron. Soc., 434, 720, doi: 10.1093/mnras/stt1065

  2. [2]

    S., Miesch, M., & Toomre, J

    Augustson, K., Brun, A. S., Miesch, M., & Toomre, J. 2015, Astrophys. J., 809, 149, doi: 10.1088/0004-637x/809/2/149 Barnab´ e, R., Strugarek, A., Charbonneau, P., Brun, A. S., & Zahn, J.-P. 2017, Astron. Astrophys., 601, A47, doi: 10.1051/0004-6361/201630178

  3. [3]

    2016, Liv

    Basu, S. 2016, Liv. Rev. Sol. Phys., 13, doi: 10.1007/s41116-016-0003-4

  4. [4]

    Batchelor, G. K. 1953, Quart. J. R. Meteor. Soc., 79, 224–235, doi: 10.1002/qj.49707934004

  5. [5]

    , keywords =

    Blume, C. C., Hindman, B. W., & Matilsky, L. I. 2024, Astrophys. J., 966, 29, doi: 10.3847/1538-4357/ad27d1

  6. [6]

    Braginskii, S. I. 1958, Soviet Journal of Experimental and Theoretical Physics, 6, 358

  7. [7]

    doi:10.1016/j.physrep.2005.06.005 , url =

    Brandenburg, A., & Subramanian, K. 2005, Phys. Rep., 417, 1, doi: 10.1016/j.physrep.2005.06.005

  8. [8]

    P., Browning, M

    Brown, B. P., Browning, M. K., Brun, A. S., Miesch, M. S., & Toomre, J. 2010, Astrophys. J., 711, 424, doi: 10.1088/0004-637x/711/1/424

  9. [9]

    M., Christensen-Dalsgaard, J., Dziembowski, W

    Brown, T. M., Christensen-Dalsgaard, J., Dziembowski, W. A., et al. 1989, Astrophys. J., 343, 526, doi: 10.1086/167727

  10. [10]

    S., Miesch, M

    Brun, A. S., Miesch, M. S., & Toomre, J. 2011, Astrophys. J., 742, 79, doi: 10.1088/0004-637x/742/2/79

  11. [11]

    S., & Toomre, J

    Brun, A. S., & Toomre, J. 2002, Astrophys. J., 570, 865, doi: 10.1086/339228

  12. [12]

    , year = 2014, month = aug, volume =

    Charbonneau, P. 2014, Ann. Rev. Astron. Astrophys., 52, 251, doi: 10.1146/annurev-astro-081913-040012

  13. [13]

    Charbonneau, P., & MacGregor, K. B. 1997, Astrophys. J., 486, 502, doi: 10.1086/304485

  14. [14]

    Caulfield, C.-c. P. 2022, J. Fluid Mech., 933, A22, doi: 10.1017/jfm.2021.1060

  15. [15]

    V., Anderson , E

    Christensen-Dalsgaard, J., D¨ appen, W., Ajukov, S. V., et al. 1996, Sci., 272, 1286, doi: 10.1126/science.272.5266.1286

  16. [16]

    Clark, A. 1973, J. Fluid Mech., 60, 561, doi: 10.1017/s0022112073000340

  17. [17]

    1999, Par

    Glatzmaier, G. 1999, Par. Comp., 25, 361, doi: 10.1016/s0167-8191(99)00009-5

  18. [18]

    Cope, L., Garaud, P., & Caulfield, C. P. 2020, J. Fluid Mech., 903, A1, doi: 10.1017/jfm.2020.600

  19. [19]

    R., & Gough, D

    Elliott, J. R., & Gough, D. O. 1999, Astrophys. J., 516, 475, doi: 10.1086/307092

  20. [20]

    A., Edelmann, P

    Featherstone, N. A., Edelmann, P. V. F., Gassmoeller, R., et al. 2021, doi: 10.5281/zenodo.1158289

  21. [21]

    A., & Hindman, B

    Featherstone, N. A., & Hindman, B. W. 2016, Astrophys. J., 818, 32, doi: 10.3847/0004-637x/818/1/32

  22. [22]

    Ferraro, V. C. A. 1937, Mon. Not. R. Astron. Soc., 97, 458, doi: 10.1093/mnras/97.6.458

  23. [23]

    J., & Matt, S

    Finley, A. J., & Matt, S. P. 2018, Astrophys. J., 854, 78, doi: 10.3847/1538-4357/aaaab5 Forg´ acs-Dajka, E. 2004, Astron. Astrophys., 413, 1143, doi: 10.1051/0004-6361:20031569 34 Forg´ acs-Dajka, E., & Petrovay, K. 2001, Sol. Phys., 203, 195, doi: 10.1023/a:1013389631585 Forg´ acs-Dajka, E., & Petrovay, K. 2002, Astron. Astrophys., 389, 629, doi: 10.105...

  24. [24]

    T., Done C., Smith D

    Garaud, P. 1999, Mon. Not. R. Astron. Soc., 304, 583, doi: 10.1046/j.1365-8711.1999.02338.x

  25. [25]

    and LeBlanc , F

    Garaud, P. 2002, Mon. Not. R. Astron. Soc., 329, 1, doi: 10.1046/j.1365-8711.2002.04961.x

  26. [26]

    Garaud, P., & Acevedo-Arreguin, L. A. 2009, Astrophys. J., 704, 1, doi: 10.1088/0004-637x/704/1/1

  27. [27]

    Garaud, P., & Brummell, N. H. 2008, Astrophys. J., 674, 498, doi: 10.1086/524837

  28. [28]

    O., & Matilsky, L

    Garaud, P., Gough, D. O., & Matilsky, L. I. 2025, Astrophys. J., 985, 151, doi: 10.3847/1538-4357/adc72b

  29. [29]

    A., & Glatzmaier, G

    Gilman, P. A., & Glatzmaier, G. A. 1981, Astrophys. J. Supp., 45, 335, doi: 10.1086/190714

  30. [30]

    O., & McIntyre, M

    Gough, D. O., & McIntyre, M. E. 1998, Nat., 394, 755, doi: 10.1038/29472

  31. [31]

    K., de Gouveia Dal Pino, E

    Guerrero, G., Smolarkiewicz, P. K., de Gouveia Dal Pino, E. M., Kosovichev, A. G., & Mansour, N. N. 2016, Astrophys. J., 819, 104, doi: 10.3847/0004-637x/819/2/104

  32. [32]

    Hathaway, D. H. 2015, Liv. Rev. Sol. Phys., 7, 1, doi: 10.12942/lrsp-2010-1

  33. [33]

    H., McIntyre, M

    Haynes, P. H., McIntyre, M. E., Shepherd, T. G., Marks, C. J., & Shine, K. P. 1991, J. Atm. Sci., 48, 651, doi: 10.1175/1520-0469(1991)048⟨0651:otcoed⟩2.0.co;2

  34. [34]

    Hazlehurst, J., & Sargent, W. L. W. 1959, Astrophys. J., 130, 276, doi: 10.1086/146715

  35. [35]

    2009, Liv

    Howe, R. 2009, Liv. Rev. Sol. Phys., 6, 1, doi: 10.12942/lrsp-2009-1

  36. [36]

    2007, Astron

    Kim, E.-J., & Leprovost, N. 2007, Astron. Astrophys., 468, 1025, doi: 10.1051/0004-6361:20065971

  37. [37]

    Korre, L., & Featherstone, N. A. 2021, Astrophys. J., 923, 52, doi: 10.3847/1538-4357/ac2dea

  38. [38]

    Korre, L., & Featherstone, N. A. 2024, Astrophys. J., 964, 162, doi: 10.3847/1538-4357/ad2844

  39. [39]

    Kraft, R. P. 1967, Astrophys. J., 150, 551, doi: 10.1086/149359

  40. [40]

    Kumar, P., & Quataert, E. J. 1997, Astrophys. J., 475, L143, doi: 10.1086/310477

  41. [41]

    Matilsky, L. I. 2023, Mon. Not. R. Astron. Soc. Lett., 526, L100–L104, doi: 10.1093/mnrasl/slad121

  42. [42]

    2024, Astrophys

    Toomre, J. 2024, Astrophys. J., 962, 189, doi: 10.3847/1538-4357/ad18b2

  43. [43]

    C., & Toomre, J

    Blume, C. C., & Toomre, J. 2022, Astrophys. J. Lett., 940, L50, doi: 10.3847/2041-8213/ac93ef

  44. [44]

    I., Korre, L., & Brummell, N

    Matilsky, L. I., Korre, L., & Brummell, N. H. 2025, Astrophys. J. Lett., 991, L1, doi: 10.3847/2041-8213/adefe3

  45. [45]

    I., & Toomre, J

    Matilsky, L. I., & Toomre, J. 2020, Astrophys. J., 892, 106, doi: 10.3847/1538-4357/ab791c

  46. [46]

    2016, Geochem., Geophys., Geosys., 17, 1586, doi: 10.1002/2015gc006159

    Matsui, H., Heien, E., Aubert, J., et al. 2016, Geochem., Geophys., Geosys., 17, 1586, doi: 10.1002/2015gc006159

  47. [47]

    Greene, T. P. 2012, Astrophys. J., 754, L26, doi: 10.1088/2041-8205/754/2/l26

  48. [48]

    Ogura, Y., & Phillips, N. A. 1962, J. Atm. Sci., 19, 173, doi: 10.1175/1520-0469(1962)019⟨0173:saodas⟩2.0.co;2

  49. [49]

    Parker, E. N. 1993, Astrophys. J., 408, 707, doi: 10.1086/172631

  50. [50]

    1987, Geophysical Fluid Dynamics (New York: Springer), doi: 10.1007/978-1-4612-4650-3

    Pedlosky, J. 1987, Geophysical Fluid Dynamics (New York: Springer), doi: 10.1007/978-1-4612-4650-3

  51. [51]

    P., Caulfield, C.-c

    Shah, K., Chini, G. P., Caulfield, C.-c. P., & Garaud, P. 2024, J. Fluid Mech., 998, A60, doi: 10.1017/jfm.2024.678

  52. [52]

    A., & Beloborodov, A

    Skoutnev, V. A., & Beloborodov, A. M. 2025, Astrophys. J., 988, 195, doi: 10.3847/1538-4357/ade547

  53. [53]

    Time Scales for CA II Emission Decay, Rotational Braking, and Lithium Depletion

    Skumanich, A. 1972, Astrophys. J., 171, 565, doi: 10.1086/151310

  54. [54]

    Spiegel, E. A. 1972, in NASA Special Publication, Vol. 300, Physics of the Solar System, ed. S. I. Rassoul, 61

  55. [55]

    A., & Zahn, J.-P

    Spiegel, E. A., & Zahn, J.-P. 1992, Astron. Astrophys., 265, 106

  56. [56]

    1962, Physics of Fully Ionized Gases, 2nd Ed

    Spitzer, L. 1962, Physics of Fully Ionized Gases, 2nd Ed. (New York: Interscience Publishers)

  57. [57]

    S., & Zahn, J.-P

    Strugarek, A., Brun, A. S., & Zahn, J.-P. 2011, Astron. Nachr., 332, 891, doi: 10.1002/asna.201111613

  58. [58]

    Vallis, G. K. 2017, Atmospheric and Oceanic Fluid Dynamics (Cambridge University Press), doi: 10.1017/9781107588417

  59. [59]

    M., Lecoanet, D., Augustson, K., et al

    Vasil, G. M., Lecoanet, D., Augustson, K., et al. 2024, Nat., 629, 769–772, doi: 10.1038/s41586-024-07315-1

  60. [60]

    E., van Horn, H

    Wendell, C. E., van Horn, H. M., & Sargent, D. 1987, Astrophys. J., 313, 284, doi: 10.1086/164968

  61. [61]

    Wood, B. E. 2018, Journal of Physics: Conference Series, 1100, 012028, doi: 10.1088/1742-6596/1100/1/012028

  62. [62]

    S., & Brummell, N

    Wood, T. S., & Brummell, N. H. 2012, Astrophys. J., 755, 99, doi: 10.1088/0004-637x/755/2/99

  63. [63]

    Zahn, J. P. 1992, Astron. Astrophys., 265, 115