pith. sign in

arxiv: 1907.02762 · v1 · pith:AOV3D6VBnew · submitted 2019-07-05 · 🌌 astro-ph.SR

The need for active region disconnection in 3D kinematic dynamo simulations

T. Whitbread , A. R. Yeates , A. Mu\~noz-Jaramillo This is my paper

Pith reviewed 2026-05-25 02:13 UTC · model grok-4.3

classification 🌌 astro-ph.SR
keywords solar dynamoactive regionssurface flux transportkinematic dynamoconvection zoneturbulent diffusivitymagnetic connectivitysolar cycle
0
0 comments X

The pith

3D kinematic dynamo models of the Sun produce mismatched surface flux evolution because active regions remain connected to the deep toroidal field.

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

The paper demonstrates that 3D simulations differ from calibrated 2D surface flux transport models because magnetic connections from active regions persist down to the base of the convection zone. Raising turbulent diffusivity breaks these connections faster and brings the surface field evolution into closer agreement for individual regions. Yet the same change stops the dynamo from sustaining itself across a full solar cycle. This points to the need for some other process that severs the links without losing net magnetic flux.

Core claim

The discrepancy between 3D kinematic dynamo models and 2D surface flux transport models arises from the connectivity of active regions to the toroidal field at the base of the convection zone. Increasing the turbulent diffusivity profile allows active regions to disconnect more rapidly and improves the match to surface field evolution, but the dynamo cannot be sustained over a full solar cycle under this enhanced diffusivity.

What carries the argument

Persistent magnetic connectivity of emerging active regions to the toroidal field at the base of the convection zone in the 3D model.

If this is right

  • Single active region decay in the 3D model improves when turbulent diffusivity is raised to speed disconnection from the base.
  • Full solar cycle runs with the same raised diffusivity fail to sustain the dynamo.
  • Any successful model must incorporate some alternative disconnection process that conserves total magnetic flux.
  • The 2D surface-only model omits the deep connections that are present in the 3D setup.

Where Pith is reading between the lines

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

  • Other 3D dynamo codes may encounter the same surface mismatch unless they adjust how emerging flux is handled or diffused.
  • Explicit modeling of reconnection events or surface flux cancellation could serve as a flux-preserving disconnection mechanism.
  • The result implies that the real Sun must employ some process that decouples active regions from deep toroidal fields without net flux loss.

Load-bearing premise

The 2D surface flux transport model, calibrated to the real Sun, supplies the correct target behavior that any 3D model must reproduce.

What would settle it

A full-cycle 3D simulation that applies an explicit, flux-conserving disconnection method to active regions and then checks whether surface flux evolution matches the 2D model while the dynamo continues to operate.

Figures

Figures reproduced from arXiv: 1907.02762 by A. Mu\~noz-Jaramillo, A. R. Yeates, T. Whitbread.

Figure 1
Figure 1. Figure 1: Three-dimensional image of an emerged active region in KD3. Magnetic field lines are connected to the toroidal field at the base of the convection zone and the radial magnetic field is shown at the transparent surface. SFT -50 0 50 Latitude (deg) -6 -4 -2 0 2 4 6 KD3 0 2 4 6 8 10 12 Time (years) -50 0 50 Latitude (deg) -6 -4 -2 0 2 4 6 [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Top: Longitude-averaged evolution of Br for a single BMR in a 2D SFT model. Bottom: Surface component of the 3D dynamo model showing the equivalent evolution of the same BMR. of [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: Two initial conditions used in this paper: an active re￾gion connected to the toroidal field (left) and an active region disconnected from the toroidal field (right). Neglecting advection, Equation 1.1 reduces to ∂A ∂t = η (z) ∇2A. (2.2) Importantly, we allow the diffusivity η to be a function of z, so that we can investigate the effect of different diffu￾sivity profiles with depth. The effect of advection… view at source ↗
Figure 3
Figure 3. Figure 3: Top: Comparison of unsigned surface flux from the 2D SFT simulation (blue) and 3D dynamo simulation (orange). Bot￾tom: Comparison of northern (solid and dotted lines) and south￾ern (dashed and dash-dotted) polar flux from the same two sim￾ulations, where polar flux is defined as the flux polewards of 70° latitude. the premise that it is the KD3 model that needs to be mod￾ified. In this paper, we show that … view at source ↗
Figure 5
Figure 5. Figure 5: Normalised multi-step diffusion profiles used in this pa￾per, against a log-scale. The solid orange curve is from KD3, the dashed purple curve is the profile which takes into account dif￾fusivity quenching, and the dotted yellow curve is derived from mixing-length theory. and ∆i respectively. The profiles used in this paper are shown in [PITH_FULL_IMAGE:figures/full_fig_p004_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Snapshots of magnetic field lines from simulations using the KD3 (top row), quenching (second row), MLT (third row) and constant (bottom row) diffusion profiles. For the first three profiles, the left-hand column shows the case where the region is connected to the base, and the middle column shows the simulation with a disconnected initial condition. The black dashed line is the top of the domain, above wh… view at source ↗
Figure 6
Figure 6. Figure 6: The left-hand panel shows magnetic field lines where [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗
Figure 8
Figure 8. Figure 8: Top: Simulation of Cycle 23 from Yeates & Muñoz-Jaramillo (2013). Bottom: Optimal but￾terfly diagram of Cycle 23 from Whitbread et al. (2017). similar time and the shape of the profile is close to that of the SFT model when exponential decay is included. In sum￾mary, these experiments with the decay of a single active region suggest that increasing the diffusivity in the bulk of the convection zone can imp… view at source ↗
Figure 9
Figure 9. Figure 9: Top: Unsigned surface flux from 3D simulations of Cycle 23 using the KD3 diffusion profile (orange), quenching profile (purple) and MLT profile (yellow). Bottom: northern polar flux (solid line) and southern polar flux (dashed line) from the same simulations. We run the simulation for 5000 days, using observed BMRs of Cycle 23 from NSO Kitt Peak as input data (Yeates et al. 2007), as in Yeates & Muñoz-Jara… view at source ↗
Figure 11
Figure 11. Figure 11: Top: Toroidal field at the base of the convection zone from a 3D simulation of Cycle 23 using the quenching profile and a strengthened initial toroidal field. Bottom: Radial magnetic field at the surface from the same simulation. ory, because the enhanced diffusivity acts as a means of disconnecting the emerged regions from the toroidal field, and this profile gave the closest match to the surface-only ev… view at source ↗
Figure 12
Figure 12. Figure 12: Top: Toroidal field at the base of the convection zone from a 3D simulation of Cycle 23 using the MLT profile and a strengthened initial toroidal field. Bottom: Radial magnetic field at the surface from the same simulation. sion in the convection zone is too strong and kills off the majority of rising flux tubes and returning poloidal flux. Even when the initial toroidal field is increased by an order of … view at source ↗
Figure 13
Figure 13. Figure 13: Top: Three-dimensional images of a disconnected ac￾tive region in KD3 with depth 0.65 R⊙ (left) and 0.95 R⊙ (right). The radial magnetic field is shown at the transparent surface and field lines below the surface. Bottom: Unsigned surface flux from 3D simulations of a single disconnected active region with depth 0.65 R⊙ (dashed), 0.8 R⊙ (dash-dotted) and 0.95 R⊙ (solid), using the quenching profile. The e… view at source ↗
read the original abstract

In this paper we address a discrepancy between the surface flux evolution in a 3D kinematic dynamo model and a 2D surface flux transport model that has been closely calibrated to the real Sun. We demonstrate that the difference is due to the connectivity of active regions to the toroidal field at the base of the convection zone, which is not accounted for in the surface-only model. Initially, we consider the decay of a single active region, firstly in a simplified Cartesian 2D model and subsequently the full 3D model. By varying the turbulent diffusivity profile in the convection zone, we find that increasing the diffusivity - so that active regions are more rapidly disconnected from the base of the convection zone - improves the evolution of the surface field. However, if we simulate a full solar cycle, we find that the dynamo is unable to sustain itself under such an enhanced diffusivity. This suggests that in order to accurately model the solar cycle, we must find an alternative way to disconnect emerging active regions, whilst conserving magnetic flux.

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

3 major / 1 minor

Summary. The paper investigates a discrepancy in surface magnetic flux evolution between a 3D kinematic dynamo model and a 2D surface flux transport (SFT) model calibrated to the Sun. It attributes this to the connectivity of active regions to the base of the convection zone in the 3D model. Experiments with a single active region in Cartesian 2D and full 3D setups show that higher turbulent diffusivity disconnects active regions faster, improving surface field match to the 2D model. However, full-cycle simulations with enhanced diffusivity fail to sustain the dynamo, leading to the suggestion that an alternative flux-conserving disconnection mechanism is needed for accurate solar cycle modeling.

Significance. If the central claim holds, the work would identify a structural limitation in 3D kinematic dynamo models arising from active-region connectivity, with implications for how emerging flux must be treated to conserve magnetic flux while matching observed surface evolution. The use of both a simplified Cartesian setup and the full 3D geometry to isolate the connectivity effect is a methodological strength.

major comments (3)
  1. [Abstract] Abstract: the claim that increasing diffusivity 'improves the evolution of the surface field' is stated without any quantitative metric (RMS error, correlation coefficient, or similar) comparing the 3D runs to the 2D SFT benchmark, so the magnitude of the improvement cannot be assessed.
  2. [Abstract] Abstract: the conclusion that an 'alternative way to disconnect emerging active regions' is required depends on the 2D SFT model being the correct target; the text provides no direct comparison of either the standard or high-diffusivity 3D runs against solar observations, leaving the necessity of an alternative mechanism untested.
  3. [Abstract] Abstract: no parameter values, functional form, or numerical values are supplied for the 'enhanced diffusivity' profile, nor is any quantitative criterion given for when the dynamo is 'unable to sustain itself,' preventing evaluation of robustness or reproducibility.
minor comments (1)
  1. The abstract refers to a 'simplified Cartesian 2D model' without stating its relation to the full spherical 3D geometry or the boundary conditions employed.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their thorough review and valuable feedback on our manuscript. We address each of the major comments point by point below.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the claim that increasing diffusivity 'improves the evolution of the surface field' is stated without any quantitative metric (RMS error, correlation coefficient, or similar) comparing the 3D runs to the 2D SFT benchmark, so the magnitude of the improvement cannot be assessed.

    Authors: We agree with this observation. The abstract would benefit from a quantitative metric to support the claim. In the revised version, we will include a specific measure, such as the RMS error or correlation coefficient between the 3D and 2D surface field evolutions, to quantify the improvement achieved with enhanced diffusivity. revision: yes

  2. Referee: [Abstract] Abstract: the conclusion that an 'alternative way to disconnect emerging active regions' is required depends on the 2D SFT model being the correct target; the text provides no direct comparison of either the standard or high-diffusivity 3D runs against solar observations, leaving the necessity of an alternative mechanism untested.

    Authors: The 2D SFT model is used as the target because it has been closely calibrated to solar observations in the literature. Our study identifies a limitation in 3D models relative to this established benchmark. We will revise the abstract to explicitly state that the 2D model is observationally calibrated, thereby strengthening the motivation for seeking an alternative disconnection mechanism. However, performing new direct comparisons to observations is outside the scope of the current work. revision: partial

  3. Referee: [Abstract] Abstract: no parameter values, functional form, or numerical values are supplied for the 'enhanced diffusivity' profile, nor is any quantitative criterion given for when the dynamo is 'unable to sustain itself,' preventing evaluation of robustness or reproducibility.

    Authors: We acknowledge that the abstract lacks these specific details due to its brevity. We will update the abstract to include the functional form and approximate values of the enhanced diffusivity profile, as well as a quantitative criterion for the dynamo's inability to sustain itself, such as the toroidal field amplitude falling below a threshold after one cycle. revision: yes

Circularity Check

0 steps flagged

No significant circularity; 3D-2D comparison uses external benchmark calibrated to observations

full rationale

The paper demonstrates a discrepancy between its 3D kinematic dynamo and a 2D surface flux transport model calibrated to the real Sun, attributes it to active-region connectivity, shows that raising diffusivity improves surface evolution at the cost of dynamo failure, and concludes an alternative disconnection mechanism is needed. This chain rests on simulation outputs compared against an independent external benchmark rather than any self-definition, fitted parameter renamed as prediction, or load-bearing self-citation. No equation or claim reduces to its own inputs by construction, and the 2D calibration is treated as an external reference rather than derived within the paper.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The paper rests on standard kinematic dynamo assumptions (turbulent diffusivity profile, flux conservation) and the domain assumption that the calibrated 2D model is the correct reference; no free parameters or invented entities are explicitly introduced in the abstract.

axioms (1)
  • domain assumption The 2D surface flux transport model is closely calibrated to the real Sun and therefore the correct benchmark for surface field evolution.
    Stated directly in the abstract as the basis for identifying the discrepancy.

pith-pipeline@v0.9.0 · 5717 in / 1211 out tokens · 23876 ms · 2026-05-25T02:13:53.959672+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

41 extracted references · 41 canonical work pages · 1 internal anchor

  1. [1]

    Babcock, H. W. 1961, ApJ, 133, 572

    work page 1961

  2. [2]

    2006, A&A, 446, 307

    Baumann, I., Schmitt, D., & Schüssler, M. 2006, A&A, 446, 307

    work page 2006

  3. [3]

    & Nandy, D

    Bhowmik, P. & Nandy, D. 2018, Nat. Comms., 9, 5209

    work page 2018

  4. [4]

    2010, LRSP, 7, 3

    Charbonneau, P. 2010, LRSP, 7, 3

    work page 2010

  5. [5]

    2014, ARA&A, 52, 251

    Charbonneau, P. 2014, ARA&A, 52, 251

    work page 2014

  6. [6]

    1999, ApJ, 527, 445

    Charbonneau, P., Christensen-Dalsgaard, J., Henning, R., et al. 1999, ApJ, 527, 445

    work page 1999

  7. [7]

    Chatterjee, P., Nandy, D., & Choudhuri, A. R. 2004, A&A, 427, 1019

    work page 2004

  8. [8]

    V., et al

    Christensen-Dalsgaard, J., Dappen, W., Ajukov, S. V., et al . 1996, Science, 272, 1286

    work page 1996

  9. [9]

    A., Arge, C

    Dikpati, M., de Toma, G., Gilman, P. A., Arge, C. N., & White, O . R. 2004, ApJ, 601, 1136

    work page 2004

  10. [10]

    Durney, B. R. 1995, Sol. Phys., 160, 213

    work page 1995

  11. [11]

    Durney, B. R. 1997, ApJ, 486, 1065

    work page 1997

  12. [12]

    H., & McClymont, A

    Fan, Y., Fisher, G. H., & McClymont, A. N. 1994, ApJ, 436, 907

    work page 1994

  13. [13]

    Gilman, P. A. & Rempel, M. 2005, ApJ, 630, 615

    work page 2005

  14. [14]

    & de Gouveia Dal Pino, E

    Guerrero, G. & de Gouveia Dal Pino, E. M. 2008, A&A, 485, 267

    work page 2008

  15. [15]

    Guerrero, G., Dikpati, M., & de Gouveia Dal Pino, E. M. 2009, A pJ, 701, 725

    work page 2009

  16. [16]

    2012, MNRAS, 420, L1

    Guerrero, G., Rheinhardt, M., Brandenburg, A., & Dikpati, M . 2012, MNRAS, 420, L1

    work page 2012

  17. [17]

    R., & Miesch, M

    Hazra, G., Choudhuri, A. R., & Miesch, M. S. 2017, ApJ, 835, 39

    work page 2017

  18. [18]

    & Miesch, M

    Hazra, G. & Miesch, M. S. 2018, ApJ, 864, 110

    work page 2018

  19. [19]

    S., Arlt, R., et al

    Jouve, L., Brun, A. S., Arlt, R., et al. 2008, A&A, 483, 949

    work page 2008

  20. [20]

    Karak, B. B. & Miesch, M. 2017, ApJ, 847, 69

    work page 2017

  21. [21]

    Karak, B. B. & Miesch, M. 2018, ApJL, 860, L26

    work page 2018

  22. [22]

    2019, A&A, 623, A54

    Kumar, R., Jouve, L., & Nandy, D. 2019, A&A, 623, A54

    work page 2019

  23. [23]

    F., & Rouillard, A

    Kumar, R., Jouve, L., Pinto, R. F., & Rouillard, A. P. 2018, Fr ontiers in Astronomy and Space Sciences, 5, 4

    work page 2018

  24. [24]

    Leighton, R. B. 1964, ApJ, 140, 1547

    work page 1964

  25. [25]

    Leighton, R. B. 1969, ApJ, 156, 1

    work page 1969

  26. [26]

    & Charbonneau, P

    Lemerle, A. & Charbonneau, P. 2017, ApJ, 834, 133

    work page 2017

  27. [27]

    2015, Ap J, 810, 78

    Lemerle, A., Charbonneau, P., & Carignan-Dugas, A. 2015, Ap J, 810, 78

    work page 2015

  28. [28]

    Mackay, D. H. & Yeates, A. R. 2012, LRSP, 9, 6

    work page 2012

  29. [29]

    Miesch, M. S. & Dikpati, M. 2014, ApJL, 785, L8

    work page 2014

  30. [30]

    Miesch, M. S. & Teweldebirhan, K. 2016, Advances in Space Res earch, 58, 1571 Muñoz-Jaramillo, A., Dasi-Espuig, M., Balmaceda, L. A., & D eLuca, E. E. 2013, ApJL, 767, L25 Muñoz-Jaramillo, A., Nandy, D., & Martens, P. C. 2008, AGU Sp ring Meeting Abstracts, SP41A Muñoz-Jaramillo, A., Nandy, D., & Martens, P. C. H. 2011, ApJ L, 727, L23 Muñoz-Jaramillo, A....

    work page 2016

  31. [31]

    & Choudhuri, A

    Nandy, D. & Choudhuri, A. R. 2001, ApJ, 551, 576

    work page 2001

  32. [32]

    Roadmap for Reliable Ensemble Forecasting of the Sun-Earth System

    Nita, G., Angryk, R., Aydin, B., et al. 2018, arXiv e-prints [arXiv:1810.08728]

    work page internal anchor Pith review Pith/arXiv arXiv 2018

  33. [33]

    Pesnell, W. D. 2016, Space W eather, 14, 10

    work page 2016

  34. [34]

    1925, Zs

    Prandtl, L. 1925, Zs. angew. Math. Mech., 5, 136

    work page 1925

  35. [35]

    Roberts, P. H. & Soward, A. M. 1975, Astronomische Nachricht en, 296, 49

    work page 1975

  36. [36]

    H., Scherrer, P

    Schatten, K. H., Scherrer, P. H., Svalgaard, L., & Wilcox, J. M. 1978, Geophys. Res. Lett., 5, 411

    work page 1978

  37. [37]

    Schrijver, C. J. & Title, A. M. 1999, Sol. Phys., 188, 331 Schüssler, M. & Rempel, M. 2005, A&A, 441, 337

    work page 1999

  38. [38]

    Tobias, S. M. 1996, ApJ, 467, 870 van Ballegooijen, A. A. 1998, in Astronomical Society of the Pacific Conference Series, Vol. 140, Synoptic Solar Physics, ed. K. S. Bal- asubramaniam, J. Harvey, & D. Rabin, 17 van Ballegooijen, A. A. & Mackay, D. H. 2007, ApJ, 659, 1713 W ang, Y.-M., Nash, A. G., & Sheeley, Jr., N. R. 1989, Science, 245, 712 W ang, Y.-M. ...

    work page 1996

  39. [39]

    R., Muñoz-Jaramillo, A., & Petrie , G

    Whitbread, T., Yeates, A. R., Muñoz-Jaramillo, A., & Petrie , G. J. D. 2017, A&A, 607, A76

    work page 2017

  40. [40]

    R., Mackay, D

    Yeates, A. R., Mackay, D. H., & van Ballegooijen, A. A. 2007, Sol. Phys., 245, 87

    work page 2007

  41. [41]

    Yeates, A. R. & Muñoz-Jaramillo, A. 2013, MNRAS, 436, 3366 Article number, page 10 of 10

    work page 2013