{"id":"cb092af6-b989-440c-b82f-e5890fc0d1ca","arxiv_id":"2506.01416","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"A solar surface flux transport simulation with observed active regions reproduces cycle 24/25 fields and attributes the anomalous southern leading-polarity poleward flux after 2016 to intermittent active region emergence.","lead":"This paper simulates the Sun's surface magnetic field from 2010 to 2024 by feeding observed active regions into a surface flux transport model. It matches observed dipole and butterfly patterns without extra decay terms, and explains why southern poleward surges were unusually weak after 2016.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'limited impact' claim about radial diffusion and cyclic flow rests on a single hand-tuned parameter set; Section 4 concedes dipole/butterfly agreement cannot constrain source or transport parameters.","rationale":"The reader's weakest_assumption points to transport parameters and ARISE completeness; my concern sharpens the transport-parameter leg. The paper's central 'limited impact' claim is a strong inference from a single simulation with parameters chosen 'after trials' and validated against the same HMI data used to build the source. The paper's own Section 4 explicitly concedes that dipole and butterfly agreement cannot independently constrain source or transport parameters, directly undermining the inference. This is a real soft spot, but it is addressable with sensitivity tests and does not invalidate the otherwise careful anomaly analysis. The public ARISE database and its agreement with independent flux measurements provide genuine support for the source-side attribution, and the S3-AR simulations quantitatively account for the main surge. The CONDITIONAL verdict remains appropriate; no change is needed. I choose 'partial' agreement because I focus on the transport-parameter/radial-diffusion inference rather than the database-completeness branch of the reader's weakest_assumption.","tokens_in":19155,"tokens_out":9065,"duration_ms":95343,"concrete_test":"Using the ARISE database (Zenodo v3.0), run the same 2010-2024 SFT assimilation with (i) a radial-diffusion term B/tau added, tau in 5-15 yr (e.g., 8 yr as in Schrijver et al. 2002), and (ii) a meridional-flow speed modulated by +/-20% around 13 m/s with cycle phase; also scan eta in [250,500] km2/s and u0 in [10,16] m/s. Report the dipole correlation and butterfly-diagram residuals for each case. If a radial-diffusion or cyclic-flow model matches or beats the no-diffusion fit, the 'limited impact' claim is not supported; if the no-diffusion fit is uniquely best and parameters stay interior to the plausible range, the claim is strengthened.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's headline inference—that results 'are achieved without incorporating radial diffusion or cyclic variations in meridional flow speed, suggesting their limited impact'—is load-bearing but not secured by the evidence. Section 2.1 fixes eta=340 km2/s and u0=13 m/s 'after trials' and adopts a single van Ballegooijen flow profile, then validates against HMI dipole and butterfly diagrams that were also used to construct the ARISE source. A good fit from one tuned parameter set shows consistency, not that the omitted terms are unimportant. A radial-diffusion term (B/tau) or cycle-dependent flow variation could be absorbed by the chosen eta and u0. The paper itself states (Section 4) that 'axial dipole strength and visual agreement with butterfly diagrams alone may not be sufficient to independently constrain the source or the transport parameters.' Because that admission applies exactly to the tuned parameters, the 'limited impact' conclusion is an interpretation, not a demonstrated result. If another reasonable parameter set requires radial diffusion to match observations, the headline claim fails.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a continuous surface flux transport (SFT) simulation of the solar surface magnetic field from 2010 to 2024 by assimilating observed active regions (ARs) from the authors' ARISE database into a standard SFT model. The simulation is compared with SDO/HMI synoptic magnetograms, and the authors report good agreement in the axial dipole strength (correlation r=0.98), polar field reversal timing, and the magnetic butterfly diagram, including poleward surges. Based on the success of the simulation without radial diffusion or cyclic meridional flow variations, the paper concludes that these processes have limited impact on multi-cycle surface field evolution. The paper further analyzes anomalous southern-hemisphere poleward flux transport in cycle 24, attributing the post-2016 leading-polarity dominated migrations to intermittent AR emergence and long emergence intervals rather than to anti-Hale or anti-Joy ARs.","tokens_in":19332,"tokens_out":4825,"duration_ms":52832,"significance":"If the claims are fully substantiated, the paper would provide a valuable continuous simulation product spanning two solar cycles, a public database and code, and a plausible explanation for a poorly understood feature of cycle 24. The anomaly analysis, linking active longitudes and emergence timing to the southern flux pattern, is of interest to the SFT and dynamo communities. The open availability of the ARISE database and code is a clear strength. However, the headline inference about the limited impact of radial diffusion and cyclic flow variations is not yet secured, because it rests on a single hand-tuned parameter set and is partly contradicted by the paper's own admission that dipole and butterfly comparisons cannot independently constrain the source or transport parameters.","major_comments":[{"comment":"The conclusion that 'these results are achieved without incorporating radial diffusion or cyclic variations in meridional flow speed, suggesting their limited impact' is not supported by the evidence presented. In §2.1 the diffusivity η=340 km²/s and meridional flow speed u0=13 m/s are adopted 'after trials' against the same HMI observations used for validation. A successful fit with one tuned parameter set shows consistency, not that omitted terms are unimportant; a radial diffusion term B/τ or a cycle-dependent flow speed could be partially absorbed by the chosen η and u0. The paper itself states in §4 that 'axial dipole strength and visual agreement with butterfly diagrams alone may not be sufficient to independently constrain the source or the transport parameters,' which applies directly to this claim. To support the limited-impact claim, the authors should either soften it or perform controlled sensitivity experiments (e.g., adding radial diffusion, varying the flow profile or amplitude, and re-fitting to observations).","section":"Abstract; §2.1; §4"},{"comment":"The source term is constructed from HMI synoptic magnetograms, and the validation is performed against the same HMI observations. The paper admits in §3.1 that the activity belts closely match 'as expected' because the ARs are directly assimilated from those maps. The dipole correlation (r=0.98) is not fully independent evidence either, since the source term directly influences the dipole evolution. The circularity does not invalidate the simulation, but it weakens the claim that the agreement validates the model and the ARISE database. Cross-validation against an independent dataset (e.g., Wilcox Solar Observatory polar fields) or a holdout test (e.g., excluding a portion of the cycle from the source and checking the prediction) would substantially strengthen the reproduction claim.","section":"§2.2; §3.1"},{"comment":"The proposed mechanism for the post-2016 dominance of leading-polarity poleward flux in the southern hemisphere is that long intervals between AR emergences allow leading-polarity flux to spread over a broad latitude range before cancellation. While plausible and consistent with the examples shown in Figures 6 and 8, this mechanism is not directly tested. The paper does not run a controlled experiment, e.g., imposing the northern-hemisphere emergence pattern on the southern hemisphere or a synthetic continuous-emergence scenario, to show that the long intervals are the cause rather than merely a correlated factor. A quantitative test would make the attribution convincing.","section":"§3.2.3"},{"comment":"The conclusion about the anomalous southern flux transport depends on the completeness of the ARISE database for the southern hemisphere during cycle 24. The paper acknowledges known limitations in the database, including low temporal resolution of synoptic maps, absence of far-side ARs, and difficulties detecting new flux within activity complexes (§4). If a significant number of southern ARs after 2016 were missed, the deduced dominance of leading-polarity flux could be an artifact of incomplete source data. The authors should estimate the possible impact of missed ARs on the polar field and the surge patterns, or at least explicitly bound the uncertainty in the key quantities (e.g., the 46% flux fraction and the net dipole contribution).","section":"§4; §2.2"}],"minor_comments":[{"comment":"The header contains a typographical error: 'T able 1' should read 'Table 1'.","section":"Table 1"},{"comment":"The keywords line lacks a space after the colon: 'Keywords:Solar physics' should be 'Keywords: Solar physics'.","section":"Keywords"},{"comment":"The paper references 'Luo et al. (2025, in preparation)' for numerical details of the spectral code. Since the code is central to reproducibility, the authors should provide either a fuller description in the current paper or a public repository link for the code, in addition to the ARISE database.","section":"General"},{"comment":"The correlation coefficient r=0.98 is reported without a confidence interval or significance level; given the serial correlation of the dipole time series, a simple Pearson correlation may be misleading.","section":"Figure 1"},{"comment":"The discussion of the northern polar field reversal discrepancy is qualitative. A quantitative statement (e.g., the difference in reversal dates between simulation and HMI, and a comparison with WSO reversal timings) would clarify the significance of this discrepancy.","section":"§3.1"},{"comment":"The meridional flow formula is presented as an unnumbered equation; numbering it (e.g., as Eq. 2) would make subsequent references to the profile clearer.","section":"§2.1"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern raised by the reader is valid and lands on the manuscript's central claim. The paper's own §4 admission about the non-uniqueness of source and transport parameters directly undermines the headline 'limited impact' conclusion. The manuscript is otherwise solid and the data product is valuable, so I recommend major revision rather than rejection: the authors should soften the overclaim and/or add sensitivity tests, and strengthen the anomaly attribution with a more explicit test of the emergence-interval mechanism."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know up front. This is a competent, honest SFT paper with a genuinely useful deliverable: a continuous 2010-2024 simulation that assimilates HMI-detected ARs from the ARISE database and reproduces the observed axial dipole (r=0.98), reversal timing, and the main butterfly-diagram surges without radial diffusion or cyclic flow-speed variation. The second thing is that the headline 'limited impact' claim about those omitted terms is not actually demonstrated; it is a reasonable interpretation that the authors themselves partially concede in Section 4.\n\nThe genuinely new pieces are the extension through the start of cycle 25, the upgraded ARISE database (now flagging and removing non-fully-emerged ARs), and the account of the southern hemisphere's odd post-2016 behavior: leading-polarity poleward flux becomes dominant because AR emergence is intermittent, with long gaps that let leading-polarity flux spread in latitude before being cancelled. The S3-AR analysis, the active-longitude detection, and the explicit check that anti-Hale/anti-Joy ARs are not responsible are all careful. The paper ships its database and code, which is real and useful.\n\nThe weak spot is the 'limited impact' inference. The transport parameters (eta=340 km^2/s, u0=13 m/s) were chosen 'after trials' against the same HMI observations used for validation, and the source and validation share the same instrument. That is tuning, not a controlled test. The paper's own Section 4 says dipole strength and butterfly agreement alone may not be sufficient to constrain source or transport parameters; that statement applies directly to the parameter choice here. A radial decay term or a different flow profile could plausibly be absorbed into the selected eta and u0. So the claim is plausible but under-supported, and the paper would be improved by sensitivity tests or by softening the abstract.\n\nThe central anomaly analysis does not depend on that claim, though. The intermittent-emergence explanation is well grounded in the ARISE flux and dipole moments, and the comparison with Z.-F. Wang et al. (2020) is fair. This is a paper for SFT and solar-cycle modelers; it deserves a serious referee. The right outcome is likely acceptance after revision to temper the 'limited impact' conclusion, not rejection. I'd bring it to reading group and would cite the database.","headline":"Solid SFT simulation with a useful new anomaly analysis; the 'limited impact' claim about radial diffusion is tuned, not tested.","tokens_in":19918,"tokens_out":2934,"would_cite":true,"duration_ms":29245,"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":"A surface flux transport simulation using observed active regions reproduces the Sun's dipole and polar reversals from 2010 to 2024 without radial diffusion.","keywords":["surface flux transport model","solar cycle 24","axial dipole strength","polar field reversal","active region emergence","poleward flux transport","magnetic butterfly diagram","active longitudes"],"falsifier":"Rerun the same assimilation but replace the observed southern-hemisphere emergence times after 2016 with uniformly spaced synthetic active regions of identical flux, tilt, and longitude, and check whether the leading-polarity poleward migration disappears; if it persists, emergence timing alone is not the mechanism. A second check would be to rerun with a different published meridional flow profile and see whether matching HMI then requires radial diffusion.","tokens_in":18912,"feed_emoji":"☀️","tokens_out":9499,"duration_ms":92111,"temperature":0.7,"pith_summary":"Surface flux transport (SFT) simulations usually need an extra radial diffusion term or cycle-dependent meridional flow changes to match the observed polar field. This paper argues neither is required when the source is a careful assimilation of observed active regions: a continuous 2010-2024 run reproduces the solar axial dipole, polar-field reversal timing, and magnetic butterfly diagram in good agreement with SDO/HMI observations. The paper then uses the same simulation to explain an anomaly of cycle 24: in the southern hemisphere, following-polarity poleward surges were confined to 2011-2016, and after 2016 leading-polarity flux dominated the poleward transport. It attributes that reversal to intermittent active-region emergence, with a burst in September 2013-February 2015 carrying 46% of southern unsigned flux, followed by long gaps that let leading-polarity flux spread over latitude before later active regions cancel it. If correct, the result narrows what solar dynamo and prediction models need to add beyond observed flux emergence.","feed_headline":"A 14-year solar surface simulation matches the Sun's polar reversals","feed_subtitle":"Observed active regions as source reproduce dipole strength and butterfly diagram without an extra decay term.","key_machinery":"The central machinery is the surface flux transport (SFT) equation, which combines advection by differential rotation and a meridional flow with supergranular diffusion, and whose source term is built by assimilating observed active regions from the ARISE database rather than idealized bipolar magnetic regions. The key diagnostic is the axial dipole strength contributed by each active region, computed in initial and final forms ($D_i$ and $D_f$), which lets the authors trace which emergences drive which poleward surges. The argumentative mechanism for the anomaly is temporal intermittency: concentrated emergence produces strong following-polarity surges, while long gaps between emergences give leading-polarity flux time to spread in latitude and reach the poles before cancellation, making it visible in the butterfly diagram.","core_discovery":"Using the ARISE database of active regions detected from SDO/HMI synoptic magnetograms, the authors feed each observed region into the surface flux transport equation at its central-meridian passage time, with constant diffusivity ($\\eta=340\\ \\mathrm{km}^2\\,\\mathrm{s}^{-1}$), constant peak meridional flow speed ($u_0=13\\ \\mathrm{m\\,s}^{-1}$), the van Ballegooijen flow profile, and no radial diffusion term. Starting from the CR 2097 magnetogram, the run through CR 2290 reproduces the HMI axial dipole strength with correlation coefficient $r=0.98$ around 2014-2015, the timing of the cycle 24 and cycle 25 polar reversals, the main poleward surges, and the activity belts in the butterfly diagram. The paper claims this good agreement shows that radial diffusion and cyclic meridional flow variations have limited impact on multi-cycle surface field evolution when the source is well represented. For the southern hemisphere in cycle 24, it finds that the dominant negative (following-polarity) surges occur only during 2011-2016, and the 2014-2016 surge S3 is driven by active regions emerging in Carrington Rotations 2141-2160, which contribute 46% of southern unsigned flux and roughly 1 G of axial dipole, nearly the whole net southern contribution. After 2016, positive (leading-polarity) flux migrations dominate even though most active regions obey Joy's and Hale's laws; the paper attributes this to long intervals between emergences, which allow leading-polarity flux to migrate poleward across a wide latitude range before being canceled by later following-polarity flux.","pith_inferences":["An implication the authors leave implicit is that cycle-prediction schemes built on polar-field proxies could be improved by monitoring the temporal clustering of active-region emergence, since the simulation shows that gaps in emergence, not just total flux or tilt, control which polarity reaches the pole.","The same intermittency mechanism should be testable in earlier cycles: if other cycles show clumped emergence followed by quiet gaps, their butterfly diagrams should show the same leading-polarity dominance, and this could be checked directly in the publicly available database.","The paper's argument that radial diffusion is unnecessary applies to the surface field and axial dipole on a 14-year window; it does not address the deep solar interior or longer-cycle memory, where transport terms beyond surface advection and diffusion may still matter.","Because the paper itself notes that dipole and butterfly agreement cannot uniquely fix source versus transport parameters, a natural next test is to compare simulated and observed full-surface magnetic power spectra or field distributions, which would separate the effect of source completeness from flow choices."],"forward_implications":["Cycle 24 and early cycle 25 surface field evolution, including polar reversals, can be reproduced with constant transport parameters; models do not need radial diffusion or cyclic meridional flow speed to explain the observed axial dipole.","The 2014-2016 southern polar reversal is driven mainly by active regions emerging in Carrington Rotations 2141-2160, which contribute about 1 G to the axial dipole, nearly the entire net contribution from all southern active regions of cycle 24.","The post-2016 dominance of leading-polarity poleward flux in the southern hemisphere is not evidence of anti-Hale or anti-Joy active regions; normal active regions with long emergence gaps can produce it.","The temporal interval between active-region emergences is a first-order factor in poleward flux transport and should be taken into account when interpreting butterfly diagrams and predicting polar fields.","Removing repeated and not-fully-emerged active-region detections yields a source reliable enough for multi-cycle simulation, so source quality matters as much as transport-parameter choices."],"supporting_citations":[{"why":"introduces the ARISE database of active regions detected from magnetograms, the source catalog assimilated in the simulation.","marker":"R. Wang et al. (2023)"},{"why":"supplies the repeat-AR-removal and flux-balancing methods that make the database reliable as a source, and quantifies how BMR approximations misestimate axial dipole contributions.","marker":"R. Wang et al. (2024)"},{"why":"archived version 3.0 of the ARISE database plus AR maps used for the 2010-2024 simulation.","marker":"R. Wang et al. (2025)"},{"why":"provides the meridional flow profile adopted in the model.","marker":"A. A. van Ballegooijen et al. (1998)"},{"why":"introduced the radial diffusion term that this paper argues is unnecessary when the source is assimilated.","marker":"C. J. Schrijver et al. (2002)"},{"why":"a precedent showing polar fields from 1923-1985 reproduced without radial diffusion or cyclic flow variation, supporting the limited-impact claim.","marker":"A. R. Yeates et al. (2025)"},{"why":"identified ARs from CRs 2145-2159 as the dominant sources of the cycle 24 southward surge, the comparison point for the S3 attribution.","marker":"Z.-F. Wang et al. (2020)"},{"why":"gives the dynamo effectivity range used to choose and interpret the transport parameters.","marker":"K. Petrovay et al. (2020)"},{"why":"an earlier HMI-assimilating SFT simulation that reproduced cycle-24 features and showed BMR overestimation, providing the methodological baseline.","marker":"A. R. Yeates (2020)"},{"why":"extends HMI-based SFT to 2010-2023 with idealized BMRs and overestimates cycle-25 dipole, the contrast that motivates direct AR assimilation.","marker":"S. Dash et al. (2024)"}],"fun_headline_variants":["14-year Sun simulation matches polar reversals without diffusion term","Clustered sunspot emergences drove southern pole's anomalous drift","No radial diffusion needed to model Sun's surface field 2010-2024","Anomalous southern flux surge traced to sunspot clustering in cycle 24","Simulation shows active region timing drives poleward flux asymmetry"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The conclusions stand on the assumption that the ARISE database captures essentially all flux emergence that matters and that the chosen constant transport parameters (340 km²/s diffusivity and 13 m/s meridional flow with the van Ballegooijen profile) are representative; if a reasonable alternative flow profile, or missed far-side or complex active regions, would force a radial diffusion term or cycle-dependent flow to match HMI, the paper's central claims about radial diffusion and about emergence timing would no longer follow.","fun_headline_variants_meta":{"raw":{"variants":["14-year Sun simulation matches polar reversals without diffusion term","Clustered sunspot emergences drove southern pole's anomalous drift","No radial diffusion needed to model Sun's surface field 2010-2024","Anomalous southern flux surge traced to sunspot clustering in cycle 24","Simulation shows active region timing drives poleward flux asymmetry"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001054,"raw_usage":{"total_tokens":4554,"prompt_tokens":1201,"completion_tokens":3353,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":817,"completion_tokens_details":{"reasoning_tokens":3261}},"tokens_in":817,"tokens_out":3353,"duration_ms":25131,"temperature":1.0,"reasoning_tokens":3261,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T11:43:40.999304+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Rerun the same assimilation but replace the observed southern-hemisphere emergence times after 2016 with uniformly spaced synthetic active regions of identical flux, tilt, and longitude, and check whether the leading-polarity poleward migration disappears; if it persists, emergence timing alone is not the mechanism. A second check would be to rerun with a different published meridional flow profile and see whether matching HMI then requires radial diffusion.","supporting_citations":[],"review_version":1}