{"id":"ead510ea-5c1e-4ed6-8eda-364891fc7dc5","arxiv_id":"1908.04525","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"A nonlinear mean-field dynamo model attributes the extended 22-year torsional oscillation pattern to overlapping magnetic cycles plus cyclic variations of meridional circulation caused by quenching of convective heat flux.","lead":"A solar dynamo model that includes magnetic effects on heat and momentum transport reproduces the 22-year equatorward drift of zonal flows known as the extended solar cycle. It identifies overlap of magnetic cycles and magnetic quenching of convective heat transport as the two essential conditions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 22-year mechanism rests on the unvalidated magnetic quenching of eddy heat conductivity in Eq. (A13); if that closure overestimates the magnetic shadow, condition (b) and the extended pattern collapse.","rationale":"The reader's weakest assumption points to the unvalidated mean-field closures, and I agree that the eddy-heat-conductivity quenching is the truly load-bearing piece: model M8 shows that the extended torsional mode can be produced by the magnetic shadow alone, so the H(0,rho) closure is not essential for the central mechanism. The more specific risk is that Eq. (A13) is evaluated with a volume-averaged β while the real solar convection-zone field is intermittent; this could overestimate the shadow effect and the resulting meridional-circulation variations. The ablation evidence for condition (b) is also less clean than presented, since M7 alters both the heat transport and the dynamo cycle itself, and the promised special run suppressing meridional-circulation variations is not documented. These concerns do not overturn the paper's conditional verdict: the mechanism is plausible, internally consistent, and supported by multiple ablation runs, but it needs a quantitative closure validation and a documented control run before the claim can be accepted as a robust explanation. I therefore keep the reader's CONDITIONAL verdict unchanged.","tokens_in":21915,"tokens_out":8883,"duration_ms":92827,"concrete_test":"Recompute model M8 with the heat-flux quenching functions in Eq. (A13) replaced by an effective quenching using the magnetic filling factor, β_eff² = f β², for f=0.1 and f=0.01, keeping all other physics fixed. If the extended 22-year torsional pattern disappears for f≲0.1, the mechanism depends on an unrealistically strong volume-averaged magnetic shadow; if it persists, the closure concern is not load-bearing.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that magnetic quenching of convective heat transport (Eq. A13) be both correctly modeled and demonstrably necessary. The least secure link is the closure itself. The functions φ_χ^(I),(||)(β) are SOCA/mean-field results for a homogeneous mean field, but the model evaluates β with the volume-averaged axisymmetric B, reaching ~5 kG near the tachocline and ~1 kG at 0.9R. Solar dynamo fields are intermittent, so the volume-averaged field can be much weaker than the field that actually suppresses convective transport; the resulting 'magnetic shadow' and the TPB/meridional-circulation perturbations could be overestimated. No independent validation of Eq. (A13) is given. The ablation supporting condition (b) is also not clean: M7 switches off β-quenching in χij, which changes the thermodynamic state and the dynamo cycle itself (period 25 yr and 1.8 yr delay vs 22.9 yr and 1.2 yr in M1); the paper mentions an undocumented special run with meridional-circulation variations suppressed but does not list it in Table 1 or show its zonal-acceleration diagram. Thus the causal isolation of condition (b) is asserted rather than demonstrated, and the physical magnitude of the effect depends on an unvalidated closure.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a nonlinear mean-field model of the solar convection zone that couples a distributed dynamo, angular momentum transport, and heat transport. The model reproduces cyclic magnetic activity, differential rotation, and meridional circulation, and produces zonal flow variations. The central claim is that the extended 22-year pattern of solar torsional oscillations is explained by two necessary conditions: (a) overlap of successive dynamo waves in the time-latitude diagram, and (b) cyclic variations of the meridional circulation caused by magnetic quenching of convective heat transport (the magnetic shadow effect). The claim is supported by ablation runs: M3 (increased anisotropic diffusivity reducing cycle overlap) loses the extended mode, M7 (no magnetic quenching of heat conductivity) loses it, while M8 (only heat-transport quenching active) retains it. The paper also decomposes the zonal forces and shows that different parts of the torsional oscillation pattern are driven by different forces. The extended pattern is an emergent output, not an imposed boundary condition.","tokens_in":22247,"tokens_out":5694,"duration_ms":53388,"significance":"If correct, the paper provides the first self-consistent mean-field explanation of the extended 22-year torsional oscillation mode, tying it to the overlap of dynamo cycles and the magnetic shadow effect in the bulk of the convection zone. The model is ambitious in its coupling of the dynamo, differential rotation, meridional circulation, and heat transport, and the ablation tests are a genuine strength: they demonstrate that the two conditions are individually necessary in the model. The force decomposition into Lorentz force, turbulent-stress quenching, dynamo-induced Lambda effect, meridional-circulation torque, and inertial force is informative and gives a concrete physical narrative. The main weakness is that the two conditions rest on unvalidated mean-field closures and partly undocumented control runs; nevertheless, the central claim is currently defensible and the missing documentation is fixable within the manuscript's scope.","major_comments":[{"comment":"The central condition (b) rests on the magnetic quenching functions phi_chi^(I),(||)(beta) in Eq. (A13), which are analytic SOCA-style closures for a homogeneous mean field. The model evaluates beta using the volume-averaged axisymmetric B, which reaches about 5 kG near the tachocline and about 1 kG at 0.9R, whereas the real convective turbulence is intermittent; if the small-scale field is much stronger than the mean, the 'magnetic shadow' and the associated Taylor-Proudman balance and meridional-circulation perturbations could be overestimated. No independent validation of Eq. (A13) by DNS or observations is given. Please provide a sensitivity study (e.g., varying the quenching amplitude) or direct numerical validation to show that the mechanism is not an artifact of the closure.","section":"Appendix A, Eq. (A13)"},{"comment":"The paper mentions 'a special run, in which we suppressed the variations' of the meridional circulation and states it gives results 'similar to model M7,' but this run is not listed in Table 1 and its zonal-acceleration diagram is not shown. This is load-bearing because M7 itself removes all beta-quenching in chi_ij, which changes the mean thermodynamic state and the dynamo cycle period (25 yr and 1.8 yr delay vs 22.9 yr and 1.2 yr in M1), so the disappearance of the extended mode in M7 could be due to the altered background state or cycle overlap rather than to the suppression of cyclic meridional circulation variations. Please document the special run (parameters, table entry, figure) so that condition (b) is genuinely causally isolated.","section":"Sec. 3.3, Table 1"},{"comment":"The comparison with observations is qualitative: the paper admits that the tachocline shear is about twice too high, the subsurface shear about twice too small, the latitudinal width of the zonal acceleration pattern narrower than observed, and the pattern inclination different. Since the torsional oscillations are driven by the force balance derived from this background, it is not clear that the extended mode would survive with a closer match to the helioseismic rotation profile. A quantitative comparison (e.g., normalized cross-correlation of the observed and modeled time-latitude acceleration maps over a common window) would establish that the model's extended mode is the same phenomenon as the observed one.","section":"Sec. 4, Fig. 1"}],"minor_comments":[{"comment":"The phrase 'a combinations of magnetic field effects' should be 'a combination of magnetic field effects'.","section":"Abstract"},{"comment":"The entries '-/-' in the M2 row are ambiguous; the caption should state explicitly that '-/-' denotes the same full nonlinear treatment as in M1, or should list each effect.","section":"Table 1"},{"comment":"The last two panels are both labeled 'e)'; the panel showing the Taylor-Proudman balance and delta F_c/F_sun should be labeled 'f)'.","section":"Fig. 5 caption"},{"comment":"The term E dot (nabla x B) in the energy equation should have its sign and physical meaning as the Joule-heating-like term clarified, since it is not discussed elsewhere in the text.","section":"Sec. 2.1, Eq. (1)"},{"comment":"The statement that the model agrees well with helioseismology is qualified by factor-of-two mismatches; it would be helpful to quote the rms difference in the convection zone rather than only the maximum difference.","section":"Sec. 3.1"},{"comment":"The reference list contains inconsistencies (e.g., 'Kitchatinov et al. 1994' appears twice with different author orders); please unify the entries.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of the journal and the central thesis is interesting. The two major issues are (i) the undocumented 'special run' that is needed for causal isolation of condition (b), and (ii) the unvalidated quenching closure in Eq. (A13). Both are addressable. If the authors supply the missing run and an explicit sensitivity test, the paper could become acceptable. I do not see a circularity problem: the extended pattern is an emergent output, not used as input."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis paper deserves a serious look. It is, as far as I can tell, the first self-consistent mean-field dynamo model that reproduces the extended 22-year pattern of solar torsional oscillations, and it does not put that pattern in by hand. The two necessary conditions it claims—overlap of subsequent dynamo cycles, and magnetic quenching of convective heat transport that perturbs the meridional circulation and Taylor-Proudman balance—are each supported by ablation runs. That is a real step beyond the earlier 11-year torsional oscillation models.\n\nThe strongest part is the force accounting. The paper shows that no single force (Lorentz, turbulent stress, meridional circulation advection) produces the extended pattern alone; it emerges from the balance, and the meridional circulation feedback, driven by the magnetic shadow, is required. The M1/M2 comparison and the M3 run with anisotropic diffusivity that reduces cycle overlap make the cycle-overlap condition fairly convincing.\n\nNow the soft spots. The first is the closure itself: the magnetic quenching of eddy heat conductivity, Eq. (A13), comes from SOCA mean-field results and is evaluated with the volume-averaged axisymmetric field, which reaches about 5 kG near the tachocline. Real dynamo fields are intermittent, so the volume-averaged field could easily overestimate the actual suppression of convective transport by a large factor. No independent check is offered. If that overestimate is big, the meridional circulation feedback driving condition (b) is exaggerated.\n\nSecond, the ablation for condition (b) is not clean. Model M7 switches off the beta-quenching in the heat conductivity, but it also changes the dynamo cycle period (25 yr vs 22.9) and the cycle overlap delay (1.8 yr vs 1.2). That is not a controlled experiment. The paper mentions a special run with meridional circulation variations suppressed but doesn't show it or list it in Table 1. That should have been a first-class run, not a footnote.\n\nThird, the comparison with observations is qualitative. The zonal acceleration amplitudes are in the right range, but the latitudinal widths and inclinations differ, and the tachocline shear is about twice the observed value, as the authors admit. With this many free parameters, qualitative agreement is suggestive, not demonstrative. I would want a quantitative residual comparison to helioseismic data, and ideally code/data release.\n\nThese concerns are real but they are not fatal to the main claim. The model is complex and the mechanism is internally consistent; the existence of two orthogonal ablation tests, even if imperfect, plus a dedicated run where the heat-transport quenching alone drives the extended pattern, makes the central claim worth taking seriously. The paper is honest about its mismatches, and the self-citations are to prior papers in the same framework, which is normal practice here.\n\nWho is this for? Anyone working on solar torsional oscillations, solar cycle prediction, or mean-field dynamo theory. I would send it to referees: a serious referee could push for the clean meridional-circulation-suppression run, a quantitative comparison, and a sensitivity study of the heat-conductivity quenching magnitude. That is the right use of referee time.","headline":"First self-consistent mean-field model of the 22-year extended torsional oscillation pattern; the mechanism is plausible but hinges on an unvalidated heat-transport closure and needs a cleaner ablation.","tokens_in":22757,"tokens_out":3667,"would_cite":true,"duration_ms":30573,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The extended 22-year pattern of the Sun's torsional oscillations arises from overlapping dynamo waves plus magnetic quenching of convective heat transport.","keywords":["torsional oscillations","extended solar cycle","solar dynamo","mean-field model","meridional circulation","magnetic quenching","convective heat transport","helioseismology"],"falsifier":"A helioseismic measurement that resolved the solar-cycle variations of the meridional circulation in the latitude band 10 to 60 degrees and found no approximately one meter per second converging flow into the activity belts, or found such variations decoupled from the zonal flow pattern, would falsify the proposed mechanism.","tokens_in":1788,"feed_emoji":"🌞","tokens_out":2519,"duration_ms":78176,"temperature":0.7,"pith_summary":"The paper tries to establish that the Sun's extended 22-year torsional oscillation pattern, the zonal flows that migrate from high latitudes to the equator over a full magnetic cycle, can be explained by a single self-consistent mean-field model of convection-zone dynamics and the dynamo. The model reproduces the observed butterfly diagram, differential rotation, and meridional circulation, and for the first time produces the extended 22-year mode. Two conditions are necessary: the subsequent magnetic cycles must overlap in the time-latitude diagram, and the large-scale magnetic field must quench the turbulent convective heat transport, which modulates the meridional circulation and drives the equatorward migration of the zonal flow pattern. If this is right, the extended solar cycle is not a separate phenomenon but a direct consequence of the same turbulent transport processes that drive the ordinary 11-year activity cycle.","feed_headline":"Sun's 22-year torsional wave traced to overlapping magnetic cycles","feed_subtitle":"A magnetized convection-zone model reproduces the observed extended solar cycle pattern for the first time.","key_machinery":"The carrier of the argument is a coupled system of mean-field equations for the large-scale magnetic field, angular momentum, and heat transport in the solar convection zone, with the magnetically induced Lambda-effect component H(0,rho), a non-dissipative angular-momentum flux arising from density stratification, and the magnetic quenching functions in the eddy heat-conductivity tensor. The key mechanism is the magnetic shadow: the large-scale toroidal magnetic field suppresses the convective energy flux, changing the mean entropy gradient; this perturbs the Taylor-Proudman balance and drives the meridional circulation variations that, combined with the overlapping dynamo waves, produce the 22-year extended torsional oscillation pattern.","core_discovery":"The central discovery is that the extended 22-year pattern of the solar torsional oscillations emerges in a nonlinear mean-field dynamo model when two ingredients are present together: (i) a dynamo wave pattern whose successive 11-year cycles overlap in the time-latitude diagram, so that the high-latitude start of the next cycle begins before the previous cycle's equatorial branch has faded, and (ii) magnetic quenching of the eddy heat conductivity, which creates cyclic variations of the convective energy flux, the magnetic shadow, perturbs the Taylor-Proudman balance, and thereby drives cyclic variations of the meridional circulation. The meridional circulation variations act as a coherent transport term that moves the torsional wave equatorward over the full 22-year cycle. The model also shows that different parts of the observed zonal-acceleration pattern are controlled by different forces: the polar branch by the large-scale Lorentz force and meridional flow, mid-latitudes by the Lorentz force and the magnetically induced Lambda-effect, and the equatorial region by the meridional circulation variations.","pith_inferences":["If this mechanism operates on the Sun, it should also leave a trace in other solar-type stars: the duration and visibility of an extended 22-year (or longer) torsional pattern should correlate with how strongly successive magnetic cycles overlap, a quantity that asteroseismic or activity-cycle observations might constrain.","The magnetic shadow effect implies a specific phase relation between photospheric luminosity variations and the toroidal magnetic field strength; high-precision broadband photometry of the Sun could be checked against the model's predicted convective-flux variations.","Long, deep minima such as the Maunder minimum, where cycle overlap is weak or absent, should show a disrupted or missing extended torsional pattern, a prediction testable by reconstructing zonal flows from historical sunspot or coronal data.","The predicted near-surface meridional circulation variations, converging toward the activity belts with an amplitude of order one meter per second, provide a concrete target for local helioseismology with current and future instruments."],"forward_implications":["The extended solar cycle is a natural output of distributed dynamo models, not requiring a flux-transport dynamo with a separate tachocline storage region.","Both identified conditions are necessary: reducing the cycle overlap (as in model M3) or removing the heat-transport quenching (as in model M7) eliminates the extended torsional mode.","Different latitude bands of the zonal acceleration pattern are forced by different physical agents, so helioseismic maps of acceleration can be used to infer which force dominates in each region.","The model reproduces the observed phase relation between flow deceleration and active regions, and the observed zonal acceleration amplitude of about 2 to 4 times ten to the minus eight meters per second squared.","The extended mode appears for both single-cell and double-cell meridional circulation structures, making the result insensitive to that unresolved aspect of the deep flow."],"supporting_citations":[{"why":"Supplies the helioseismic observations of zonal acceleration and its overlay with the magnetic field that the model is designed to reproduce.","marker":"Kosovichev & Pipin 2019"},{"why":"Establishes the concept of the extended dynamo mode and overlapping magnetic cycles in time-latitude diagrams, which the model's cycle-overlap condition builds on.","marker":"Stenflo 1992"},{"why":"Supplies the mean-field expressions for the magnetic quenching of eddy heat conductivity and angular momentum transport used in the model.","marker":"Kitchatinov et al. 1994"},{"why":"Source of the H(0,rho) contribution to the magnetically induced Lambda-effect, a key non-dissipative angular-momentum flux in the force balance.","marker":"Kueker et al. 1996"},{"why":"The previous nonlinear dynamo model that this work extends by adding the tachocline and full feedback on heat transport.","marker":"Pipin 2018"},{"why":"Provides the anisotropic eddy-diffusion formulation used in model M3 to reduce cycle overlap and demonstrate its necessity.","marker":"Pipin & Kosovichev 2014"},{"why":"Gives earlier heuristic arguments for the role of heat-transport quenching in torsional oscillations, which the present model makes quantitative.","marker":"Rempel 2006"},{"why":"An MHD simulation baseline that showed the Lorentz force drives high-latitude branches but did not reproduce the full extended pattern.","marker":"Guerrero et al. 2016a"}],"fun_headline_variants":["Overlapping magnetic cycles drive Sun's 22-year torsional wave","Magnetic heat quenching explains Sun's extended zonal flow pattern","Model ties solar torsional oscillations to overlapping dynamo cycles","Sun's 22-year wave traced to overlapping cycles and heat quenching","New model reproduces extended solar cycle torsional pattern"],"cache_read_input_tokens":24832,"weakest_assumption_plain":"The whole explanation depends on the assumption that the mathematical formulas used for the magnetic suppression of turbulent heat transport and angular-momentum transport in the convection zone match what really happens in the Sun.","fun_headline_variants_meta":{"raw":{"variants":["Overlapping magnetic cycles drive Sun's 22-year torsional wave","Magnetic heat quenching explains Sun's extended zonal flow pattern","Model ties solar torsional oscillations to overlapping dynamo cycles","Sun's 22-year wave traced to overlapping cycles and heat quenching","New model reproduces extended solar cycle torsional pattern"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000191,"raw_usage":{"total_tokens":1364,"prompt_tokens":986,"completion_tokens":378,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":602,"completion_tokens_details":{"reasoning_tokens":293}},"tokens_in":602,"tokens_out":378,"duration_ms":4650,"temperature":1.0,"reasoning_tokens":293,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:40:00.380665+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A helioseismic measurement that resolved the solar-cycle variations of the meridional circulation in the latitude band 10 to 60 degrees and found no approximately one meter per second converging flow into the activity belts, or found such variations decoupled from the zonal flow pattern, would falsify the proposed mechanism.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the concept of the extended dynamo mode and overlapping magnetic cycles in time-latitude diagrams, which the model's cycle-overlap condition builds on."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Source of the H(0,rho) contribution to the magnetically induced Lambda-effect, a key non-dissipative angular-momentum flux in the force balance."}],"review_version":1}